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

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

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

因为应变与位移分量的关系式为u—0rrru10{ }zr wr0zzrzuwzrzr的倒数项1/r ,则即使位移模式为线性的, 即不会是常应变应力应变的物理关系为 则所求得的应力也不会是常应力—，这里除含有微分算符外，还包含了 rw但由于该项的存在，使得应变与坐标有关，D ，由于应变不是常应变，3•在薄板弯曲理论中做了哪些假设？薄板单元和厚板单元的基本假设有什么不同？ （ 10分）答：四种假设：1） 变形前的中面法线在变形后仍为弹性曲面的法线2） 变形前后板的厚度不变3） 板变形时，中面无伸缩4） 板各水平层间互不挤压不同点：薄板单元假设横向纤维无挤压， 板的中面法线变形后仍保持为直线， 该直线垂直于变形后的中面， 但是厚板单元的假设考虑横向变形的影响， 板的中面法线变形后仍基本保持为直线，但该直线不再垂直于变形后的中面，法线绕坐标轴的转角不再 是挠度的导数，而是独立的变量三、计算题（3道，共计45分）1•如图所示等腰直角三角形单元，其厚度为t，弹性模量为E，泊松比 0；单元的边长及结点编号见图中所示求(1) 形函数矩阵N(2) 应变矩阵B和应力矩阵S(12 分)(3) 单元刚度矩阵Kea解：设图1所示的各点坐标为点1(a，0)，点2(a, a)，点3(0，0)于是，可得单元的面积为(1)形函数矩阵N为N12(0axay)aN12(0a0xay)N12 (a2aax0|y)NI N1in2IN3N1N2N3(2)应变矩阵B和应力矩阵S分别为Aa00A0-aA0B1 ~20-a ， B20a ，B3-2 00aaa-aaa00-aB B1 B2 B3a000-a0S1 -Er 0-a ， S20a ， S300-lalala00a2222SD B1 B2 B3SiS2 S3(3)单元刚度矩阵KK11K12K13Ke BTDBtAK21K22K23K31K32K33311021131201Et11100140202002000201110012. 如图所示的四结点矩形单元，求出节点3的位移。

设厚度t = 1m, [i= 0, E为常量13分)團2注：对于四节点矩【形单元有•1.1N1—1141N2411—N1N3—1141N4—1142 .k11eTk21kBDB tdxdyAk31k411丄(1 i )(1 i ) (i 1,234)4k12k13ek14k22k23k24 ，k32k33k34k42k43k44kj Bi D Bj tdxdy abt〔〔 B D Bj d dAEt 2ai1j 1 3 i j1ab i1j 1 3 i j2 ij 1 i j28 1^a,1a,12 ij 1i j2— ij 1 i j1ij 1 -bj 3 jbj 3(i,j 123,4)解：对于四节点矩形单元有:1Ni - 1 1424f Ni 1(1i )(1i ) (i 1,2,3,4)N311414N41141N2ek14k24k34k444k11k12k13e Tk21k22k23k B DB tdxdyAk31k32k33k41k42k43T11 Tkj Bi D Bjtdxdyabt11 BD 1Bj d dAb1a2-11i j 1Et a3 i jb81 22 ij 1i j(i,jee1i j2 ij 1i j3 j1 ja1a12 ij 1 i j1—ij 1 -b3b31,2,3,4)R e，代入边界条件皿=v =国=v=慨=v= 0,将对应的行和列划掉没剩下的方程为:-P-P1又Ni丄（14)(1 i ) (i 1,2,3,4)，且 3 1， 3 1，a=1,b=1所以k33所以-P-P-P E 4 1 悶-> -P 8 1 4 v8P 15E 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对于单元①：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]= - 000042020824200弹性矩阵[D]=单元刚度矩阵2(94(1))1TEt2(1)[k][B] [D][B]t一 2、16(1)8164(1)2(9)4(1)4(32)2“J 2(1) 4(1)2(1[k][k]84016804(1) 8(1)4(1)然后合成总刚[K]) 348 16 10 0 8(1 )4( 1)2(1)140480816 38( 1)4( 1)00 8(1 )24(3 2 )P p整体节点力矢量为{F} [F1x F1y F2x F2y 0 0]节点位移矢量为{d} 0 0 0 0 u3 v3 u4 v4{F} [K]{d}，采用缩减矩阵法划去位移为零的行与列，得。