精选有限元考试试题及答案第一组资料.docx

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

2、些区别？ （5 分）答：轴对称单元是三角形或四边形截面的空间的环形单元， 平面单元是三角形或四边形平 面单元； 轴对称单元内任意一点有四个应变分量， 平面单元内任意一点非零独立应变分量有 三个。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则所求得的应力也不会是常应力。3在薄板弯曲理论中做了哪些假设？薄板单元和厚板单元的基本假设有什么不同？ （ 10分）答：四种假设：1） 变形前的中面法线在变形后仍为弹性曲面的法线。2） 变形前后板的厚度不变。

5、3） 板变形时，中面无伸缩。4） 板内各水平层间互不挤压。不同点：薄板单元假设横向纤维无挤压， 板的中面法线变形后仍保持为直线， 该直线垂直于变形后的中面， 但是厚板单元的假设考虑横向变形的影响， 板的中面法线变形后仍基本保持为直线，但该直线不再垂直于变形后的中面，法线绕坐标轴的转角不再 是挠度的导数，而是独立的变量。三、计算题（3道，共计45分）单元的边长及结点编号见图中所示。求(1)形函数矩阵N(2)应变矩阵B和应力矩阵S单元刚度矩阵Ke ( 12分)解：设图1所示的各点坐标为点1(a，0)，点2(a, a)，点3(0，0)于是，可得单元的面积为A*，及(1)形函数矩阵N为1N| 亍(0

6、ax - ay)N -丨1 N1 IN2 I N3 1-IN 1 N2 N3 1aM 1 (0 dx ay) ； a叫 1 (a2 - ax y)a(2)应变矩阵B和应力矩阵S分别为-a00-a03 =耳0-a,S2 E2_0a,S3 = _E200aaa1-la1 a1 a010-a一 22 -2t2 -S = D 丨 B B2 B31-IS-| S2 S3 I(3)2.如图所示的四结点矩形单元，求出节点3的位移。设厚度t = 1m, i= 0, E单元刚度矩阵Ke3-1-1 -10-2-201 1-131K12 K13k K -K 22 K 23 A-11100-140-20200K 32

7、 K 33 J-2000201-1-1001 一K11e tK = B DBtA= K21K 31为常量。(13分)4-3注：对于四节点矩形单元有:NiN2N3N4J1- 14111 41 141十一 14Ni(1 i )(1 i ) (i =123,4)1.2.k11k12k13k14k21k22k23k24k31k32k33k34“41k42k43k44 _D 治 tdxdy 二Akj LqBi j D Bj tdxdy =abt IB ID Bj dAEt812-ja1-1ij21 1 一j-i j1：lijII，a b i j(i,j =1,2,3,4)解：对于四节点矩形单元有:NiN

8、2N3N411- 14J1 1J i 1 4J 1- 14Ni1=丄(1 i )(1 i )4(i =123,4)k e = Bb IB tdxdy 二Ak11 k21 k31 k41k13k23k33k43-.ek14k24k34k44k L H Bij b Bj tdxdy =abt D tj d A _ 一；b- 2_ 1 -jaEt81 3 i j 1_吒 i j1 i j2a i j 1 -b -(i,j ,2,3,4)2 j1 -i -j3-吒i j13ijkeeRe， 代入边界条件 (11= V= M2= V=出=V= 0,将对应的行和列划掉没剩下的方程为:(J3V-P-P1又

9、Ni =丄(1 】)(1 i ) (i =123,4)，且；=1, 3=1 , a=1,b=14所以xm、 “ 1訂4丁一 I31丁一3所以-P-P鉀/曰爲:8P 1解侍丿 = 丿沁：5E 13.有一如图3(a)所示的剪力墙，墙顶作用竖向荷载 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。对于单元

10、：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_400 0-401几何矩阵B=-020 -200824-2 00-4J单元刚度矩阵弹性矩阵D=一2（9-卩）4（1 + ）4（3-2 円2岸1）4岸-1）2（1 -卩）-8卩-404-168卩084 164（卩-1）8（卩-1）4（卩-1）0 0k1 珂BTDBt.： =8（1 _ JEt216（1）29-門24（1 + 卩）4（3-24）2（4-1）4（卩-1）2（1 _ 門3-8-404-16-8卩08卩161卅-1）8-1）4岸-1）0 0 8（1 -k2 二k1然后合成总刚K。P p整体节点力矢量为F -F1x F1y F2x F2y 0 0节点位移矢量为d = 0 0 0 0 u3 v3 u4 vF二Kd，采用缩减矩阵法划去位移为零的行与列，得-2（9 -卩）1-P/2Et04（3_2卩）V30216（1 -卜）-164（卩-1）2（9 -卩）U4P/2一8卩8（卩-1）4（卩 +1）4（3-2卩）一4 一_ 0 一解得-1.88U4单元的应变t =Bd1 二 P 0.285Et-0.470.56 T单元的应力P 08-1.680.84T