牛-力学班有限元试题07A答案.doc

合肥工业大学

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

课程代码考试时间：2007年6月3日上午7：30-9：30

答案：

―――――――――――――――――――――――――――――――――――――

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

假设单元的位移场模式

代入到几何方程得到

代入到物理方程得到

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

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

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

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

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

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

4．答：

空间轴对称问题3节点等参元，其单元自由度6个，单元刚度元素是36个；

空间问题8节点等参元，其单元自由度24个，单元刚度元素是576个；

空间问题20节点等参元，其单元自由度60个，单元刚度元素是3600个；

平面刚架结构梁单元〔考虑轴向和横向变形〕的自由度6个，单元刚度元素是36个。

5．解：(1)结点编号如图。

〔2〕整体结点位移列向量

〔3〕系统刚度矩阵半带宽

B=2*3+2=8

〔4〕结构的整体等效结点力列向量

〔5〕约束结点自由度的位移信息

6．解：

〔1〕单元=1\*GB3①的等效结点应力

24

单元=2\*GB3②的等效结点应力

213

结构整体的等效结点力

1234

〔2〕单元=1\*GB3①的刚度矩阵

〔取i=2，j=4，α=90°〕

24

[k]=1\*GB3①=0

单元=2\*GB3②的刚度矩阵

〔取i=2，j=1，m=3〕

213

[k]=2\*GB3②=0

结构总刚

[k]=[k]=1\*GB3①+[k]=2\*GB3②

1 234

=E0

采用划行划列法，引入结点位移得：

解得：

〔3〕由于4节点固定，且

所以杆=1\*GB3①内力为

〔拉力〕

7．解：单元的一致质量矩阵

整体一致质量矩阵

8．解：：

引入，在原振动方程两边同时左乘，得到

(*)

其中：，

记

那么〔*〕式化为：

由杜哈米积分有：

其中为待定常数，由于

，且为常数矩阵

所以初始条件可以转化为

从而可以确定出

由此，

那么：