《有限元考试精彩试题及问题详解——第一组.doc》由会员分享，可在线阅读，更多相关《有限元考试精彩试题及问题详解——第一组.doc（12页珍藏版）》请在三一文库上搜索。
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 分）答：不是常应力和常应变。因为应变与位移分量的关系式为 ：.u.r1.:ruzrz:wcz：u : w:z : r,这里除含有微分算符外，还包含了rW J大全的倒数项1/r,则即使位移模式为线性的，但由于该项的存在，使得应变与坐标有关，即不会是常应变。应力应变的物理关系为匸；-bH，由于应变不是常应
5、变，则所求得的应力也不会是常应力。3. 在薄板弯曲理论中做了哪些假设？薄板单元和厚板单元的基本假设有什么不同？ （ 10分）答：四种假设：1 ）变形前的中面法线在变形后仍为弹性曲面的法线。2）变形前后板的厚度不变。3）板变形时，中面无伸缩。4）板内各水平层间互不挤压。不同点：薄板单元假设横向纤维无挤压，板的中面法线变形后仍保持为直线，该直线垂直于变形后的中面， 但是厚板单元的假设考虑横向变形的影响，板的中面法线变形后仍 基本保持为直线，但该直线不再垂直于变形后的中面，法线绕坐标轴的转角不再是挠度的导数，而是独立的变量。三、计算题(3道，共计45分)1.如图所示等腰直角三角形单元，其厚度为t，弹
6、性模量为E，泊松比=0 ;解：设图1所示的各点坐标为点1( a，0)，点2( a，a),点3( 0，0)1 2于是，可得单元的面积为A=a2,及2(1)形函数矩阵N为1N|2 (0 ax - ay)a1叫 2 (0 lx ay) ；aN - IlN1 IN2 I N3 1-In 1 N2 N3 1M $ (a2 - ax 出y)a应变矩阵B和应力矩阵S分别为_a1 |B1 = I 0.a0-aa12a00'a0la ，B 30'-a.000-aB =B1B 2单元的边长及结点编号见图中所示。求(1)形函数矩阵N(2)应变矩阵B和应力矩阵S（12 分）单元刚度矩阵Ke-a00-a
7、03 =耳0-a,S2 'E20a,S3 = -E200aaa1-la1 a1 a010-a一 22 -2t2 -S = D 丨 BB2B31-IS-|S2S3 Ie(3)单元刚度矩阵"KuK e = BT DBtA 二 K21K 31K12K 22K 32K13K 23K 33=Et_ 43-11-101-2-131-20-1-11100-10-21 1-20-100-1200020001 一2.如图所示的四结点矩形单元，求出节点3的位移。设厚度t = 1m 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 二Ak21k41标准文案kj L "Bij D 血 tdxdy =
9、abt j j 辰D Bj d：dn A*Et811 3b(i,j =1,2,3,4)2j 1仝i j2a i j 11 j31大全keuae， 代入边界条件卩1=V 1= 2 =V 2=卩4 =v 4= 0，将对应的行和列划掉没剩下的方程为:-P(i =1,2,3,4)，且 3=1，3=1， a=1,b=1又 n =丄(1)(1 i )4所以12131313山勻4】所以解得丿3'-P-P8 L1 4V-P8P 15E 13.有一如图3(a)所示的剪力墙，墙顶作用竖向荷载P。将该剪力墙划分为两个三结点三角形常应力单元，单元和结点编号如图3(b)所示，并将荷载P分成两个P/2作用在3、4
10、结点。已知单元厚度为t，弹性模量为E,泊松比卩=1/3。求结点3和结点4的位移，以及单元的应变和应力。(20分)解：建立直角坐标系(注Y轴向下为正)，单元i,j,m 对应的节点编号为3,1,4，单元对应的节点编号为2,4,1。;bj=ym-yi=O ; bm=yi_yj=_4;cj=xi-xm=-2 ; cm=xj-xi=O对于单元：i(0,0),j(0,4),m(2,0)bi=yi-ym=4A=1/2*2*4=4ci=xm-xj=2三角形面积4 0 1几何矩阵B= - 0820 -2-2 0-4-4标准文案弹性矩阵D=皂1 以0【01-卩Y单元刚度矩阵k1 =BTD Bt.l =2（9，）4
11、（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.99-1.88TU4EtF =Kd，采用缩减矩阵法划去位移为零的行与列，得2（9）U3P/2Et04（3_2卩）V3016（1 厂）-164（卩-1）2（9 -卩）U4P/2一8卩8（卩-1）4（卩 +1）4（3-2叽V 一0 一单元的应变H 询d1 唁 0.285 心7 0.56单元的应力J 冷DBKd10.58 -1.68 0.84T