合肥工业大学力学班有限元试题08A答案.docVIP

合肥工业大学

工程力学专业《工程结构数值方法》2008年试卷A答案

2008

(共8题，满分100分)，考试时间：120分钟

适用专业：工程力学05－1（本科），出题人：牛忠荣

（8分）线弹性力学静力学问题有限元法计算列式的理论推导是如何采用弹性力学问题基本方程？

答：弹性力学有限元的基本过程是：

假设单元的位移场模式

代入到几何方程得到

代入到物理方程得到

代入到虚功方程，得到单元刚度方程

叠加到总刚阵，得到结构的平衡方程

引入位移边界条件后，解第5步得到的方程组，可以得到结点位移

（12分）用有限元法计算图示平面刚架时,试求：

(1)如何进行结点编号,使系统总刚度矩阵的带宽最小?

(2)计算在你的结点编号下的系统总刚阵的半带宽。

(3)在结点编号确定后,按此顺序进行自由度编号,则图中结点水平位移对应的主对角元素在总刚阵中的列和行位置是多少?那些单元对该元素的数值有贡献?

题2图答：

题2图

（1）见图；

（2）半带宽＝；

（3）结点水平位移对应的

主对角元素是；

A点周围的4个单元对

有贡献。

（10分）弹性力学有限元中，平面等参数单元中的“等参数”概念是何意思？该单元在跨相邻单元时，位移场连续吗？应力场连续吗？

答：在有限单元法中最普遍采用的是等参变换，即单元几何形状的变换和单元内的场函数采用相同数目的节点参数及相同的插值函数进行变换。采用等参变换的单元称之为等参元。所谓“等参元”是指几何形状插值形函数和单元上的位移插值形函数相同，参数个数相等。

相邻等参元之间，位移场是连续的，应力场不连续。

（10分）回答下列问题：

弹性力学平面问题8节点等参元，其单元自由度是多少个？单元刚阵元素是多少个？

弹性力学空间问题20节点等参元，其单元自由度是多少个？单元刚阵元素是多少个？

一般的平面刚架结构梁单元的自由度是多少个？单元刚阵元素是多少个？

一般的空间刚架结构梁单元的自由度是多少个？单元刚阵元素是多少个？

答：1）自由度16个，刚阵元素16×16＝256；

2）自由度60个，刚阵元素60×60＝3600；

3）自由度6个，刚阵元素6×6＝36；

4）自由度12个，刚阵元素12×12＝144.

（10分）结构振动问题有限元离散的无阻尼自由振动方程为

式中是刚度矩阵，是质量矩阵，是结构固有频率，是振型向量。试问为什么从上式求出的特征对()中，只有前若干低阶频率和相应振型是可靠的，误差较小。

答：在有限单元法中，采用低阶多项式拟合振型。结构的低阶振型曲线与低阶多项式比较通配，结构的高阶振型曲线与低阶多项式曲线有着显著的差异。因而，有限元法中求出的低阶频率和振型是可信的，而所求出的高阶频率和振型误差较大，甚至无效。

（7分）用有限元法求结构的动力响应通常有振型叠加法和逐步积分法，试问该两种方法的优缺点。

答：振型叠加法适应于求解线性振动问题，需要计算结构的固有频率和振型。适合于时间较长求解。

逐步积分法可用于非线性振动问题的求解，不必需要计算结构的固有频率和振型。适合于时间较短的问题求解，当时间长时计算量大。

(25分)图示矩形薄板，厚度。载荷及约束信息如图示，泊松比,，板的重量密度为,沿轴负向。试采用三角形常应变元求：

（1）建立适宜的有限元计算模型；

（2）结构系统的等效结点力列阵；

（3）采用划行划列法引入已知结点位移，

计算出结点3和4的位移；

（4）单元1和2的各应力分量。

（建议单元编号：

ijm

单元1：342

单元2：213)

题

题7图

题7图解：计算模型

题7图

单元1的等效结点力：

结点342

单元2的等效结点力：

结点213

结构整体等效结点力

结点1234

单元1：

单元2：

单元1应力矩阵

结点342

；

单元2应力矩阵

结点213

；

单元1刚度矩阵：

单元1：

，，

，

342

；

213

；

叠加总刚阵：

结点号：1234

已知结点位移：

＝；

采用划行划列，得到

解得结点位移：，。

代入：

得单元应力（常应力）：

，。

(18分)