有限元考试精彩试题及问题详解——第一组

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. 简述四节点四边形等参数单元的平面问题分析过程。（io分）答：（1）通过整体坐标系和局部坐标系的映射关系得到四节点四边形等参单元的母单元，并选取单元的唯一模式；（2）通过坐标变换和等参元确定平面四节点四边形等参数单元的几何形状和位移模式；（

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

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

6、比=0 ；(12 分)2（a, a），点 3（0, 0）(1)形函数矩阵N为12 (0 ax _ ay) a12 (0 lx ay) ； a2 (a2 - ax 0jy) aB和应力矩阵S分别为01NiNi应变矩阵_a1 |B1 = I 0.a01-aa12a_00'a 0jN - IlN1 IN2 I N3 1-IN 1 N2 N3 1'-a0'.000 ； B = lB1 B 2 民-a单元的边长及结点编号见图中所示。求（1）形函数矩阵N（2）应变矩阵B和应力矩阵S单元刚度矩阵Ke解：设图1所示的各点坐标为点1（a，0），点1 2于是,可得单元的面积为人=丄，及2-

7、a00-a03 =耳0-a,S2 E2_0a,S3 = _E200aaa1-la1 a1 a010-a一 22 -2t2 -S = D 丨 BB2B31-IS-|S2S3 I(3)单元刚度矩阵Ke3-1-10-21 1-131-20-1K12K 22K13IK 23Et _ 4-101-2100200-10K 32K 33 -2000201-1-1001 一K11etK = B DBtA= K21K 312. 如图所示的四结点矩形单元，求出节点3的位移。设厚度t = 1m, i= 0, E为常量。(13分)注：对于四节点矩形单元有:1.NiN2N3N4二 1- 1 -4= 71-441 41=

8、- 14Ni1)(1 i )4(i =1,2,34)2.knk12k13k14k21k22k23k24k31k32k33k34J<41k42k43k44k e = BID IB tdxdy =A1 1k L Bi I D B tdxdyA二abt.“81-卩）2P1占十（1_卩卿j12i j T-打j亠吒,(i,j =1,23,4)解：对于四节点矩形单元有:N1N2N3N414厂=*+駅 1+r4(1-行计)4Ni(1：)(1i )4(i =123,4)knk22k23k24k34k44k e = BD tdxdy 二Ak21k41kj L "Bij D 血 tdxdy =ab

9、t j j IB D Bj d：dn11 3A*Eti j 1 -Ia i j 1 二 丿 b I 2 2j 1 七 j(i,j =123,4)kURe, 代入边界条件(1-1= v=忐=v=也=v= 0,将对应的行和列划掉没剩下的方程为:-p1又 n =丄(1)(1 i )4(i =1,2,3,4)，且 3=1，3=1， a=1,b=1所以-13山勻4J所以f13P-P8 J 4i v-P解得丿f138P 15E 13. 有一如图3(a)所示的剪力墙，墙顶作用竖向荷载P。将该剪力墙划分为两个三结点三角形常应力单元，单元和结点编号如图3(b)所示，并将荷载P分成两个P/2作用在3、4结点。已知

10、单元厚度为t，弹性模量为E,泊松比1 =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=4_4000-40 11几何矩阵B= 200_4 一弹性矩阵D=皂1 以0【01-卩Y单元刚度矩阵k1 =BTD Bt.l =2（9，）4（1+ 巴12 -1） 8卩 164（卩-1）2（9）4（1 + »）_2（4 1） -8-16 卅-1）然后合成总刚K。Et16（1 -二2）k2 可k14（3_2i）4（-1）-4-88（二一1）2（1 -）004（二一1）4（3-2）4（ 二 _1）-4-8J8（-1）2（1 7004（-1）48J08（1 -二）160整体节点力矢量为 F二RF1yF 2xF2y节点位移矢量为d = 0000u3V3P2u4 v4 10解得8.421.50 8.