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

江西理工大学研究生考试试卷
20__12__—20__13__ 学年 第___一___学期
课程名称：_____有限元及数值模拟________
考试时间:___2012___ 年__11__月___3___日
考试性质（正考、补考或其它）：[ 正考 ]
考试方式(开卷、闭卷)：[ 开卷 ]
试卷类别(A、B):[ A ] 共 九 大题
温 馨 提 示
请考生自觉遵守考试纪律，争做文明诚信的大学生。如有违犯考试纪律，将严格按照《江西理工大学学生违纪处分规定》（试行）处理。
学院 专业 学号 姓名
题号
一
二
三
四
五
六
七
八
九
十
十一
十二
总 分
得分

一、 简答题（共40分，每题10分）
1. 论述单元划分应遵循的原则。
2. 说明形函数应满足的条件。
3. 说明四边形等参数单元中“等参数”的含义，即为什么要引入等参数单元。
4. 阐述边界元法的主要优缺点。
二、 计算题（共60分，每题20分）
1. 一杆件如图3所示，杆件上方固定后，在下方受垂直向下的集中力作用，已知：杆件材料的杨氏模量，截面积，，长度，集中力，用有限元方法求解B点和C点位移。备注：（1）1 lbf（磅力，libra force） = 4.45 N。（2）杨氏模量、弹性模量、Young氏弹性模量具有相同含义（10分）
A1
A2
L1
L2

图1

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

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

图3
一、简答题
1. 答：
1）合理安排单元网格的疏密分布
2）为突出重要部位的单元二次划分
3）划分单元的个数
4）单元形状的合理性
5）不同材料界面处及荷载突变点、支承点的单元划分
6）曲线边界的处理，应尽可能减小几何误差
7）充分利用结构及载荷的对称性，以减少计算量
2. 答：
形函数应满足的三个条件：
a. 必须能反映单元的刚体位移，就是位移模式应反映与本单元形变无关的由其它单元形变所引起的位移。
b. 能反映单元的常量应变，所谓常量应变，就是与坐标位置无关，单元内所有点都具有相同的应变。当单元尺寸取小时，则单元中各点的应变趋于相等，也就是单元的形变趋于均匀，因而常量应变就成为应变的主要部分。
c. 尽可能反映位移连续性；尽可能反映单元之间位移的连续性，即相邻单元位移协调。
3. 答：
含义：所谓的等参数单元，就是在确定单元形状的插值函数和确定单元位移场的插值函数中采用了完全相同的形函数。
意义：构造出一些曲边地高精度单元，以便在给定地精度下，用数目较少地单元，解决工程实际地具体问题。
4. 答：
有限单元法是基于变分原理的里兹（Ritz）法的另一种形式，从而使里兹法分析的所有理论基础都适用子有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法．利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，面且事先不要求满足任何边界条件，因此它可以用来处理很复杂的连续介质问题。有限单元法中所利用的主要是伽辽金（Galerkin）法。它可以用于已经知道问题的微分方程和边界条件，但变分的泛函尚未找到或者根本不存在的情况，因而进一步扩大了有限单元法的应用领域。
三十多年来，有限单元法的应用已由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定问题、动力问题和波动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、传热学等连续介质力学领域。在工程分析中的作用已从分析和校核扩展到优化设计并和计算机辅助设计技术相结合。可以预计，随着现代力学、计算数学和计算机技术等学科的发展，有限单元法作为一个具有巩固理论基础和广泛应用效力的数值分析工具，必将在国民经济建设和科学技术发展中发挥更大的作用，其自身亦将得到进一步的发展和完善。

三、 计算题
1. 解：将杆件分解成两个元素，
元素1的刚度矩阵K1=，
元素2的刚度矩阵K2=
总刚度矩阵
单元刚度矩阵形成后，应将各单元刚度矩阵组装集合成整体刚度矩阵（即总刚矩阵）。如所示为杆系结构两单元节点编号示意图，可得总刚度矩阵为
(3-11)

图3-3 杆系结构两单元节点编号示意图
引入边界条件求解节点位移
总刚矩阵组集完成后，即可获得整个结构的平衡方程为
(3-12)
整个结构的边界条件为，、已知，三个未知量三个方程，因此上式可求得唯一解。

解出节点位移：
u1=0
u2=0.762×10-5in
u3=0.18295×10-4in
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：

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

如图所示标出了各个单元的应力值，而且在单元内是不变的，这就说明了是一近似解。在单元交界处，应力值有突变，这就可以看出，如将单元分得很细，则突变减小，其结果将会改善。

计算后的各单元应力
3. 解：
N1=L1(2L1-1)；N2=L2(2L2-1)；N3=L3(2L3-1)；N4=4L2L2；N5=4L2L3；N6=4L3L1；
在三角形142边上L3=0；可得N3=N5=N6=0，因此有：
； ；
在三角形142边上建立局部坐标系如图：
1 4 2
1
0
-1

1，4，2节点在局部坐标系下的形函数为：
；；
所以：
；
；
；
其中，l为三角形142边的边长；
；
；
；
其中：θ为q与水平方向的夹角。