有限元考试试题及答案

1、江西理工大学研究生考试试卷20122013学年第一学期课程名称：有限元及数值模拟考试时间：2012年11月3日考试性质（正考、补考或其它）：正考考试方式（开卷、闭卷）：开卷试卷类别（A、B）:A共九大题温馨提示请考生自觉遵守考试纪律，争做文明诚信的大学生。如有违犯考试纪律，将严格按照江西理工大学学生违纪处分规定（试行）处理。学院专业学号姓名题号一二三四五六七八九十十十总分得分一、简答题（共40分，每题10分）论述单元划分应遵循的原则。说明形函数应满足的条件。说明四边形等参数单元中“等参数”的含义，即为什么要引入等参数单元。阐述边界元法的主要优缺点。二、计算题（共60分，每题20分）一杆件如图3

2、所示，杆件上方固定后，在下方受垂直向下的集中力作用，已知：杆件材料的杨氏模量E1二E2二107吋/in2,截面积A1二5.25in2,A2二75in2,长度L1=L2=12加，集中力P=100lbf，用有限元方法求解B点和C点位移。备注：（1）1lbf（磅力，libraforce）=4.45N。（2）杨氏模量、弹性模量、Young氏弹性模量具有相同含义（10分）如图2所示，有一正方形薄板，沿对角承受压力作用，厚度t=lm,载荷F=20KN/m，设泊松比p=0,材料的弹性模量为E,试求它的应力分布。（15分）等效结点荷载。一、简答题合理安排单元网格的疏密分布为突出重要部位的单元二次划分划分单元的

3、个数单元形状的合理性不同材料界面处及荷载突变点、支承点的单元划分曲线边界的处理，应尽可能减小几何误差充分利用结构及载荷的对称性，以减少计算量1.答：1)2)3一、简答题合理安排单元网格的疏密分布为突出重要部位的单元二次划分划分单元的个数单元形状的合理性不同材料界面处及荷载突变点、支承点的单元划分曲线边界的处理，应尽可能减小几何误差充分利用结构及载荷的对称性，以减少计算量1.答：1)2)3)4)5)6)7)2.答：形函数应满足的三个条件：必须能反映单元的刚体位移，就是位移模式应反映与本单元形变无关的由其它单元形变所引起的位移。能反映单元的常量应变，所谓常量应变，就是与坐标位置无关，单元内所有点都

4、具有相同的应变。当单元尺寸取小时，则单元中各点的应变趋于相等，也就是单元的形变趋于均匀，因而常量应变就成为应变的主要部分。尽可能反映位移连续性；尽可能反映单元之间位移的连续性，即相邻单元位移协调。答：含义：所谓的等参数单元，就是在确定单元形状的插值函数和确定单元位移场的插值函数中采用了完全相同的形函数。意义：构造出一些曲边地高精度单元，以便在给定地精度下，用数目较少地单元，解决工程实际地具体问题。4.答：有限单元法是基于变分原理的里兹（Ritz）法的另一种形式，从而使里兹法分析的所有理论基础都适用子有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法利用变分原理建立有限元方程和经典里兹

5、法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，面且事先不要求满足任何边界条件，因此它可以用来处理很复杂的连续介质问题。有限单元法中所利用的主要是伽辽金（Galerkin）法。它可以用于已经知道问题的微分方程和边界条件，但变分的泛函尚未找到或者根本不存在的情况，因而进一步扩大了有限单元法的应用领域。三十多年来，有限单元法的应用已由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定问题、动力问题和波动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、传热学等连续介质力学领域。在工程分析中的作用已从分析和校核扩展到优

6、化设计并和计算机辅助设计技术相结合。可以预计，随着现代力学、计算数学和计算机技术等学科的发展，有限单元法作为一个具有巩固理论基础和广泛应用效力的数值分析工具，必将在国民经济建设和科学技术发展中发挥更大的作用，其自身亦将得到进一步的发展和完善。三、计算题解：将杆件分解成两个元素,Ae元素1的刚度矩阵0=t1=13.125106Ibf/in，i_AE元素2的刚度矩阵K2=22=9.375x106Ibf/inLz2总刚度矩阵单元刚度矩阵形成后，应将各单元刚度矩阵组装集合成整体刚度矩阵（即总刚矩阵）。如所示为杆系结构两单元节点编号示意图，可得总刚度矩阵为k11K=k0210kG12k11K=k0210

7、kG12kG+kG)22:22kG)23k(2)(3-11)k(2)3233图3-3杆系结构两单元节点编号示意图引入边界条件求解节点位移总刚矩阵Ik组集完成后，即可获得整个结构的平衡方程为rf1FErf1FE2FV3丿lJ1-l10_1-l1X+ll12_l2_1_l12u1u2uJ3丿(3-12)整个结构的边界条件为u已知，三个未知量三个方程，因此上式可求得唯一解。u1l-F2112uEll+1FL3JL1121-3解出节点位移:U=0u2=0.762x10-5inu3=0.18295x10-4in解:（1）建立需要计算的力学模型以及划分单元由于该结构几何对称和受载也对称，故可利用其对称性，

8、只需要取薄板的1/4作为计算对象。为了简单起见，我们把它划分成4个三角形单元，单元和节点编号如图（b）所示。由于对称，节点1，2,4,不可能有水平位移，节点4,5，6不可能有垂直位移，故施加约束如图（b）所示。图两类单元节点编号图两类单元节点编号取总体x,y坐标并确定各节点的坐标值。由图看出，这里只有两类不同的单元，一类单元是1，2,4,另一类单元是3。两类单元节点的编排如图所示。单元1,单元节点编排对应于结构的节点编号1,2,3。三个节点坐标如下:x=0 x=0,ix=0,jy=2m,y=1m,ij代入得：b=y一y=0；bijmjc=x一x=1；cimjj三角形面积:x1m，my1mmy一

9、y二二1；byy=1mimijxx二=1；cxx-0immji单元节点坐标以及单元和节点的编号是原始数据，可用手工输入，也可由计算机完成。对于单元2，3，4定出单元节点的坐标值后，同样可算出，以及各单元的面积。(2)计算个单元的刚度矩阵ke及组集成总刚K由于t=1m,卩=0，所以Et41-2丿,于是由式可求得单元刚度矩阵为=E=Ek1k1k1k1k1k1iiijim111213=k1k1k1=k1k1k1jijjjm212223k1k1k1k1k1k1mimjmm3132330.2500.250.2500.2500.500.5000.2500.750.250.50.250.250.50.250

10、.7500.25000.500.500.2500.250.2500.25同理可得单元2，4的刚度矩阵分别为k2k2k2k4k4k4222425333536e=k2k2k2,ke=k4k4k424244454535556k2k2k2k4k4k4525455636566由于1,2,4单元算出的b,b等值以及三角形面积均相同，故算出2,4ij的单元刚度矩阵与单元1的刚度矩阵数值完全相同。单元3的节点i,j,m相应于总体编号中的2,5,3点，其节点坐标为x0,x1m,x1mTOC o 1-5 h zijmy1m,y0,y1mijm由此得：b-1m,b0,b1m,ijmc0,c-1m,c1m.ijm从而

11、算出单元刚度矩阵为：77q7qk3k3iiijimek3k3k33jik3kjjk3kjmk3kmimjmmk322k325k323k3k3k35235553k3323k3353k3330.5000-0.5000.250.250-0.25-0.2500.250.250-0.25-0.250000.50-0.15-0.5-0.25-0.2500.750.250-0.25-0.25-0.50.250.75E根据各单元刚度矩阵组集成总刚度矩阵b为(Kk1(Kk111k100013k1+k3k2k2+k302323242525k1+k3+k40k3+k4k4333333353536k2k204445k

12、2+k+k4k455555556k4k112k1+k2+k3222222对称）66由以上结果求得总刚度矩阵各元素为k=k1=E11110.25k=k1=E11110.2500.5k=k1=E12120.2500.250.5k=k1=E13130.250.50.750.25,E0.250,+E0.50,二E1.50.25,0.250.7500.500.250.251.5k=k1k2k3=E2222222224k=k2=E44440.750.250.250.75-0.5-0.25_0.50,24k=k2=E44440.750.250.250.75-0.5-0.25_0.50,+E0.250,+E0

13、.750.25,二E1.50.25,00.25_00.5_0.250.75_0.251.5_-0.250-0.5-0.25k=k2=E4545k=k4=E5656k=k2k3k4=E55555555-0.25k=k4=E66660.50.25-0.5-0.25,+E-0.50,二E-1-0.25,0-0.25_-0.25-0.25_-0.25-0.5_-0.25-0.25k=k1k3=E232323=k2=E00.25,+E00,二E00.25,000.2500.250-0.5024k=k2k3=E2525250.50,+E0.750.25,+E0.250,二E1.50.25,00.250.2

14、50.7500.50.251.5k=k1k3k4=E333333331111-0.250,-0.25-0.25,-0.5-0.25,k=k3k4=E+E二E353535-0.25-0.5_0-0.5_-0.25-10.250k=k4=E3636把上面计算出的k11，k对号入座放到总刚矩阵Ik中去，于是得到Ik把上面计算出的k1166的具体表达式。3）计算并代入等效节点载荷及相应的位移边界条件，以建立和求解未知节点位移的平衡方程组。先求出各项等效节点载荷然后叠加，以形成方程组右端载荷项，但本问题只在节点1有一个集中外载荷R=10kN/m（取F二20kN/mi0v0v的一半)。0v0v0v0v由结

15、构的对称性，可以看出uu12u4v4v5v0。于是需要求的未60v0v代入边界条件及外载荷以代入边界条件及外载荷以知节点位移分量只有6个，即v,v,u,v,u12335及支反力后,其方程组为0v0v0v0v0.250-0.25-0.2500.2500000000.50-0.500000000-0.2501.50.25-1-0.25-0.25-0.2500.2500-0.25-0.50.251.5-0.25-0.50-0.50.2500000-1-0.251.50.2500-0.5-0.2500.250.250-0.25-0.50.251.500-0.25-10000-0.250000.750.

16、25-0.5-0.250000-0.25-0.5000.250.750-0.25000000.25-0.5-0.25-0.501.50.25-0.5-0.25000.250-0.25-1-0.25-0.250.251.50-0.2500000000-0.500.5000000.25000-0.25-0.2500.25E把左端系数矩阵行列倒换,于是可分块求解。第二种办法是把带有支反力的10v2u2v300u50u60F1x-10F2x000F4xFy0Fy0F6y0v0v方程去掉,即把系数矩阵中的第1,3,7个未知位移的方程式：8,10,12行和列划掉,得出带有60v0v0v0v0.5-0.50

17、-0.51.5-0.250-0.251.50-0.50.25E0-0.50.251.500.25-0.25-0.25000000、(v)1vJ10、0.25002-0.50u3v0-0.250031.5-0.5u05-0.50.5丿Iu6J0丿0v0v0v0vuu3V3u5u丿6解此方程组即得到位置节点位移分量。由上方程组求得位移分量如下：32.52/E、-12.52/E-0.88/E-3.72/E1.76/E1.76/E,(4)求单元应力分量求出节点位移分量后,就可以按式计算单元中的应力33我们略去初应变，于是有0对于单元1，2，4：00-1000-10ts=E010-10.50-0.5-0

18、.5100000.53333最后可求得各单元的应力为对于单元3：最后可求得各单元的应力为00100-100-1-0.500050.5-10k=E000-0.5注意至Uu=u=u=v=v=v=0,厂00-1x-E010y、0.50-0.5xyV0-10-1-0.50、00.5丿-0.88-20.00kN/m2、4.40丿V3丿3333厂00-10 x仝E010-1y、0.50-0.5-0.5xyV1000、00.5丿v2003yxy0-1010-10-0.5-0.5厂00、0.51000、00.5丿1.76-12.52kN/m2-0.88-3.72kN/m2、43.08丿9,00-10 x0刍E010-1yT0.50-0.5-0.5xy丿0、00.5丿0-3.72kN/m2-132丿6如图所示标出了各个单元的应力值，而且在单元内是不变的,这就说明了是一近似解。在单元交界处，应力值有突变，这就可以看出，如将单元分得很细,则突变减小，其结果将会改善。计算后的各单元应力3.解：N1=L1