有限元考试试题及答案.doc

有限元考试试题及答案

江西理工大学研究生考试试卷

20__12__—20__13__学年第___一___学期

课程名称：_____有限元及数值模拟________

考试时间:___2012___年__11__月___3___日

考试性质（正考、补考或其它）：[正考]

考试方式(开卷、闭卷)：[开卷]

试卷类别(A、B):[A]共九大题

温馨提示

请考生自觉遵守考试纪律，争做文明诚信的大学生。如有违犯考试纪律，将严格按照《江西理工大学学生违纪处分规定》（试行）处理。

学院专业学号姓名

题号

一

二

三

四

五

六

七

八

九

十

十一

十二

总分

得分

简答题（共40分，每题10分）

论述单元划分应遵循的原则。

说明形函数应满足的条件。

说明四边形等参数单元中“等参数”的含义，即为什么要引入等参数单元。

阐述边界元法的主要优缺点。

计算题（共60分，每题20分）

1.一杆件如图3所示，杆件上方固定后，在下方受垂直向下的集中力作用，已知：杆件材料的杨氏模量，截面积，，长度，集中力，用有限元方法求解B点和C点位移。备注：（1）1lbf（磅力，libraforce）=4.45N。（2）杨氏模量、弹性模量、Young氏弹性模量具有相同含义（10分）

A1

A1

A2

L1

L2

图1

图1

2.如图2所示，有一正方形薄板，沿对角承受压力作用，厚度t=1m，载荷F=20KN/m，设泊松比μ=0，材料的弹性模量为E，试求它的应力分布。（15分）

图2

3.图示结点三角形单元的124边作用有均布侧压力q，单元厚度为t，求单元的等效结点荷载。

图3-3杆系结构两单元节点编号示意图

引入边界条件求解节点位移

总刚矩阵组集完成后，即可获得整个结构的平衡方程为

(3-12)

整个结构的边界条件为，、已知，三个未知量三个方程，因此上式可求得唯一解。

解出节点位移：

u1=0

u2=0.762×10-5

u3=0.18295×10-4

2.解：

(1)建立需要计算的力学模型以及划分单元

由于该结构几何对称和受载也对称，故可利用其对称性，只需要取薄板的1/4作为计算对象。为了简单起见，我们把它划分成4个三角形单元，单元和节点编号如图（b）所示。由于对称，节点1，2，4，不可能有水平位移，节点4，5，6不可能有垂直位移，故施加约束如图（b）所示。

图两类单元节点编号

取总体坐标并确定各节点的坐标值。由图看出，这里只有两类不同的单元，一类单元是1，2，4，另一类单元是3。两类单元节点的编排如图所示。

单元1，单元节点编排对应于结构的节点编号1，2，3。三个节点坐标如下：

，，，

，，

代入得：

；；

；；

三角形面积:

单元节点坐标以及单元和节点的编号是原始数据，可用手工输入，也可由计算机完成。对于单元2，3，4定出单元节点的坐标值后，同样可算出，以及各单元的面积。

（2）计算个单元的刚度矩阵及组集成总刚

由于，，所以

于是由式可求得单元刚度矩阵为

同理可得单元2，4的刚度矩阵分别为

，

由于1，2，4单元算出的…等值以及三角形面积均相同，故算出2，4的单元刚度矩阵与单元1的刚度矩阵数值完全相同。

单元3的节点相应于总体编号中的2，5，3点，其节点坐标为

由此得：

从而算出单元刚度矩阵为：

根据各单元刚度矩阵组集成总刚度矩阵为

由以上结果求得总刚度矩阵各元素为

把上面计算出的，…，对号入座放到总刚矩阵中去，于是得到的具体表达式。

（3）计算并代入等效节点载荷及相应的位移边界条件，以建立和求解未知节点位移的平衡方程组。先求出各项等效节点载荷然后叠加，以形成方程组右端载荷项，但本问题只在节点1有一个集中外载荷（取的一半）。

由结构的对称性，可以看出。于是需要求的未知节点位移分量只有6个，即，，，，。代入边界条件及外载荷以及支反力后，其方程组为

把左端系数矩阵行列倒换，于是可分块求解。第二种办法是把带有支反力的方程去掉，即把系数矩阵中的第1，3，7，8，10，12行和列划掉，得出带有6个未知位移的方程式：

解此方程组即得到位置节点位移分量。由上方程组求得位移分量如下：

(4)求单元应力分量

求出节点位移分量后，就可以按式计算单元中的应力。

我们略去初应变，于是有

对于单元1，2，4：

对于单元3：

注意到，最后可求得各单元的应力为

如图所示标出了各个单