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

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分）答：不是常应力和常应变。因为应变与位移分量的关系式为 ：，这里除含有微分算符外，还包含了wrz的倒数项1/r ,则即使位移模式为线性的,但由于该项的存在，使得应变与坐标有关,即不会是常应变。应力应变的物理关系为D ,由于应变不是常应变,则所求得的应力也不会是常应力。3 .在薄板弯曲理论

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

6、点编号见图中所示。形函数矩阵N应变矩阵B和应力矩阵(12 分)单元刚度矩阵Ke解：设图1所示的各点坐标为点于是,可得单元的面积为形函数矩阵NiNiN为 4(0 a 工(0 a1 z 2 (a aax0g(ax(2)应变矩阵-a点 3 (0, 0)1 (a, 0),点 2 (a, a),120a ,及2ay)B和应力矩阵-aB212 a0gy)S分别为B31 2 a-a00IN1in2IN3N1N2N3-aB1 B2B3Sia1 a-a2-a ， S200与 0 a , S3 a1 a 02-a与0 a0(3)单元刚度矩阵KeD B1B2SiS2S3K11 K12 K13e _ TK B DBt

7、A K21K22K23K31 K32 K33Et431113111102 020002 12010012000201110012.如图所示的四结点矩形单元，求出节点3的位移。设厚度t = 1m,仙=0,为常量。(13分)注：对于四节点矩形单元有:1N11141N2114 1N3 - 1141N4 - 1141一 Ni -(1 i )(1 i ) (i 1,2,3,4)4kjk11k12k13ek14eTk21k22k23k24,kB D B tdxdyAk31k32k33k34k41k42k43k44T111BiD Bj tdxdy abt11 B D Bj d dA81Et2b a对于四节点

8、矩形单元有:(i,j 1,2,3,4)1N1- 1141N2 - 114- Ni 1(1 i )(1 i ) (i 1,2,3,4)14N3 - 1141N4114ek11k12k13k14e- T .k21k22k23k24B D B tdxdyAk31k32k33k34k41k42k43k44TkjBiD Bj tdxdyAabtD Bj dEt812b ia13ij1"(i,j1,2,3,4)ke ee，代入边界条件国=1= 22= 2>=国=必=0,将对应的行和列划掉没剩下的方程为:-P-P又Ni4(1)(1(i123,4),且1,31 , a=1,b=1所以k33所以

9、P3-P-PP3-P-P解得P38P5E3.有一如图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 。对于单元：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=44 00 0401几何矩阵B=0 2020082 42 004弹性矩阵D=单元刚度矩阵k1EtW2)2(9 4(1 2() )1)4(3 2 )4(k2 k14(2(94(12(8161)8( ) )1)1)481)4(3 24(然后合成总刚8164(1)K。8(整体节点力矢量为F节点位移矢量为d2(1004()1)481)F1x1)1608(12(1004(1)F1yF2xU31608(1F2yV30U4V4F Kd,采用缩减矩阵法划去位移为零的行与列，得2(9)Et04(3 2 )2一