有限元考试试题及答案——第一组

1、有限元考试试题及答案一、简答题（ 5 道，共计 25 分）。1. 有限单元位移法求解弹性力学问题的基本步骤有哪些？（ 5 分）答：（1）选择适当的单元类型将弹性体离散化；（2）建立单元体的位移插值函数；（3）推导单元刚度矩阵；（4）将单元刚度矩阵组装成整体刚度矩阵；（5）代入边界条件和求解。2. 在划分网格数相同的情况下，为什么八节点四边形等参数单元精度大于四边形矩形单元？ （5 分）答：在对于曲线边界的边界单元， 其边界为曲边， 八节点四边形等参数单元边上三个节 点所确定的抛物线来代替原来的曲线，显然拟合效果比四边形矩形单元的直边好。3. 轴对称单元与平面单元有哪些区别？ （5 分）答：轴对

2、称单元是三角形或四边形截面的空间的环形单元， 平面单元是三角形或四边形平面单元； 轴对称单元内任意一点有四个应变分量， 平面单元内任意一点非零独立应变分量有 三个。4. 有限元空间问题有哪些特征？ （5 分）答：（1）单元为块体形状。常用单元：四面体单元、长方体单元、直边六面体单元、曲 边六面体单元 、轴对称单元。 （ 2）结点位移 3 个分量。（ 3）基本方程比平面问题多。 3 个 平衡方程， 6 个几何方程， 6 个物理方程。5. 简述四节点四边形等参数单元的平面问题分析过程。 （5）分） 答：（1）通过整体坐标系和局部坐标系的映射关系得到四节点四边形等参单元的母单元， 并选取单元的唯一模

3、式；（ 2） 通过坐标变换和等参元确定平面四节点四边形等参数单元的几何形状和位移模式；（ 3）将四节点四边形等参数单元的位移模式代入平面问题的几何方程，得到单元应变 分量的计算式，再将单元应变代入平面问题的物理方程，得到平面四节点等参 数单元的应力矩阵；（4）用虚功原理求得单元刚度矩阵，最后用高斯积分法计算完成。二、论述题（ 3 道, 共计 30 分）。1. 简述四节点四边形等参数单元的平面问题分析过程。 （ 10 分）答：（1）通过整体坐标系和局部坐标系的映射关系得到四节点四边形等参单元的母单元， 并选取单元的唯一模式；2）通过坐标变换和等参元确定平面四节点四边形等参数单元的几何形状和位移模

4、式；得到单元应变（3）将四节点四边形等参数单元的位移模式代入平面问题的几何方程,分量的计算式，再将单元应变代入平面问题的物理方程，得到平面四节点等参数单元的应力矩阵；（4）用虚功原理求得单元刚度矩阵，最后用高斯积分法计算完成。2. 轴对称问题的简单三角形单元是否是常应力，常应变？为什么？（10 分）答：不是常应力和常应变。因为应变与位移分量的关系式为r zrzu0rru10 ,rru,这里除含有微分算符外，还包含了rw0_ wzzuwzrzr的倒数项1/r ,则即使位移模式为线性的，但由于该项的存在，使得应变与坐标有关，即不会是常应变。应力应变的物理关系为D，由于应变不是常应变，则所求得的应力

5、也不会是常应力。3. 在薄板弯曲理论中做了哪些假设？薄板单元和厚板单元的基本假设有什么不同？ （ 10分）答：四种假设：1 ）变形前的中面法线在变形后仍为弹性曲面的法线。2）变形前后板的厚度不变。3）板变形时，中面无伸缩。4）板内各水平层间互不挤压。不同点： 薄板单元假设横向纤维无挤压， 板的中面法线变形后仍保持为直线， 该直线垂直于变形后的中面， 但是厚板单元的假设考虑横向变形的影响， 板的中面法线变形 后仍 基本保持为直线， 但该直线不再垂直于变形后的中面， 法线绕坐标轴的转角不再 是挠度的导数，而是独立的变量。三、计算题( 3 道，共计 45 分)。1.如图所示等腰直角三角形单元，其厚度

6、为 t ，弹性模量为 E ，泊松比 0； 单元的边长及结点编号见图中所示。求(1) 形函数矩阵 N(2) 应变矩阵 B 和应力矩阵 S(3) 单元刚度矩阵 Ke(12 分)解：设图1所示的各点坐标为点1 (a, 0),点2 (a, a),点3 (0, 0)于是，可得单元的面积为 ，及(1) 形函数矩阵为(2) 应变矩阵和应力矩阵分别为(3) 单元刚度矩阵2.如图所示的四结点矩形单元，求出节点 3的位移。设厚度t = im 厂0,E为常量。（13分）注：对于四节点矩形单元有:1N11114n2111N4N31114N411142 .k11eTk21kBDB tdxdyAk31k4i1丄(1 i

7、)(1 i ) (i1,234)4k12k13ek14k22k23k24，k32k33k34k42k43k44T11TkjBi D Bj tdxdy abt B D Bj d dAEt 2ai1j 1 3 i j1ab i1j 13 i j2ij1i j28 1a,1a,12ij1i j2 ij 1i j1ij 1-bj3 jbj3(i,j 123,4)解：对于四节点矩形单元有:1N,- 1141N24一 Ni !(1 i )(1 i ) (i 1,2,3,4)14N311341N4114k11k12k13eTk21k22k23kB DB tdxdyAk31k32k33k41k42k43T11

8、 TkijBi D Bjtdxdyabt 11 BD 1Bj d dAb1a2-11i j 1Eta3 i jb81 22ij1i jek14k24k34k441i j2ij1i j3 ja1a12ij 1i j1ij 1-b3b34(i,j123,4) eekRe，代入边界条件P 1=V 1=P 2=v 2=p 4=v 4= 0,将对应的行和列划掉没剩下的方程为：,P-3-P33V-P1又Nj丄（1i )(1i )(i 1,2,3,4)，且 31,31 , a=1,b=1所以111 -k33所以-PE 4 1 (J3-P8 1 4 v-P-P8P 15E 13.有一如图3(a)所示的剪力墙，

9、墙顶作用竖向荷载P。将该剪力墙划分为两个三结点三角形常应力单元，单元和结点编号如图3(b)所示，并将荷载P分成两个P/2作用在3、4结点。已知单元厚度为t，弹性模量为E,泊松比卩=1/3。 求结点3和结点4的位移，以及单元的应变和应力。(20分)解：建立直角坐标系(注Y轴向下为正)，单元i,j,m 对应的节点编号为 3,1,4，单元对 应的节点编号为2,4,1。对于单元：i(0,0),j(0,4),m(2,0)bi=yi-ym=4; bj=ym-yi=0 ; bm=yi_yj=_4ci=xm-xj=2; cj=xi-xm=-2 ; cm=xj-xi=0三角形面积A=1/2*2*4=441几何矩阵B= - 000040202008242004弹性矩阵D=2单元刚度矩阵2(9)44(1)4(32 ) 1TEt2(1)4(1) 2(1)1kB DBt216(1 )8404168081634(1)8(1) 4(1)00 8(1)2(9)24(1)4(3 2)22(1)4(1)2(1)3kkc8404168081614(1) 8( 1)4(1)00 8(1)然后合成总刚K。整体节点力矢量为FPF1xF1yF2xF2y0y y 220节点位移矢量为d0 000 u3 v3 u4 v4FKd，采用缩减矩阵法划去位移为零