(完整word版)有限元考试试题及答案

江西理工大学研究生考试试卷2012—2013学年第一学期课程名称：有限元及数值模拟考试时间：2012年11月3日考试性质（正考、补考或其它）：［正考］考试方式（开卷、闭卷）：［开卷］试卷类别（A、B）:［A］共九大题温馨提示请考生自觉遵守考试纪律，争做文明诚信的大学生。如有违犯考试纪律，将严格按照《江西理工大学学生违纪处分规疋》（试仃）处理。学院专业学号姓名题号-一一-二二-三四五六七八九十十十总分得分一、简答题（共40分，每题10分）论述单元划分应遵循的原则。说明形函数应满足的条件。说明四边形等参数单元中“等参数”的含义，即为什么要引入等参数单元。4•阐述边界元法的主要优缺点。二、计算题（共60分，每题20分）一杆件如图3所示，杆件上方固定后，在下方受垂直向下的集中力作用，已722知：杆件材料的杨氏模量Ei=E2=3.010lbf/in，截面积A=5.25in，・2・A2=3.75in，长度Ll二L2=12in，集中力P=10°lbf，用有限元方法求解B点和C点位移。备注：（1）1lbf（磅力，libraforce）=4.45N。（2）杨氏模量、弹性模量、Young氏弹性模量具有相同含义（10分）A1L1A2L2如图2所示，有一正方形薄板，沿对角承受压力作用，厚度t=1m，载荷F=20KN/m，设泊松比尸0,材料的弹性模量为E,试求它的应力分布。（15分）等效结点荷载图3、简答题答：1）合理安排单元网格的疏密分布2）为突出重要部位的单元二次划分3）划分单元的个数4）单元形状的合理性5）不同材料界面处及荷载突变点、支承点的单元划分6）曲线边界的处理，应尽可能减小几何误差7）充分利用结构及载荷的对称性，以减少计算量答：形函数应满足的三个条件：必须能反映单元的刚体位移，就是位移模式应反映与本单元形变无关的由其它单元形变所引起的位移。能反映单元的常量应变，所谓常量应变，就是与坐标位置无关，单元内所有点都具有相同的应变。当单元尺寸取小时，则单元中各点的应变趋于相等，也就是单元的形变趋于均匀，因而常量应变就成为应变的主要部分。尽可能反映位移连续性；尽可能反映单元之间位移的连续性，即相邻单元位移协调。答：含义：所谓的等参数单元，就是在确定单元形状的插值函数和确定单元位移场的插值函数中采用了完全相同的形函数。意义：构造出一些曲边地高精度单元，以便在给定地精度下，用数目较少地单元，解决工程实际地具体问题。答：有限单元法是基于变分原理的里兹（Ritz）法的另一种形式，从而使里兹法分析的所有理论基础都适用子有限单元法，确认了有限单元法是处理连续介质问题的一种普遍方法.利用变分原理建立有限元方程和经典里兹法的主要区别是有限单元法假设的近似函数不是在全求解域而是在单元上规定的，面且事先不要求满足任何边界条件，因此它可以用来处理很复杂的连续介质问题。有限单元法中所利用的主要是伽辽金（Galerkin）法。它可以用于已经知道问题的微分方程和22边界条件，但变分的泛函尚未找到或者根本不存在的情况，因而进一步扩大了有限单元法的应用领域。三十多年来，有限单元法的应用已由弹性力学平面问题扩展到空间问题、板壳问题，由静力平衡问题扩展到稳定问题、动力问题和波动问题。分析的对象从弹性材料扩展到塑性、粘弹性、粘塑性和复合材料等，从固体力学扩展到流体力学、传热学等连续介质力学领域。在工程分析中的作用已从分析和校核扩展到优化设计并和计算机辅助设计技术相结合。可以预计，随着现代力学、计算数学和计算机技术等学科的发展，有限单元法作为一个具有巩固理论基础和广泛应用效力的数值分析工具，必将在国民经济建设和科学技术发展中发挥更大的作用，其自身亦将得到进一步的发展和完善。计算题1.解：将杆件分解成两个元素,元素1的刚度矩阵K1=A1E1=13.125106lbf/in，L1元素2的刚度矩阵K2=A—=9.375106Ibf/inL2总刚度矩阵单元刚度矩阵形成后，应将各单元刚度矩阵组装集合成整体刚度矩阵(即总刚矩阵)。如所示为杆系结构两单元节点编号示意图，可得总刚度矩阵为下。)1kO)201=k0)1kQ)2+k(2bk(2bI0k(2)2k(233(3-11)1(il图3-3杆系结构两单元节点编号示意图引入边界条件求解节点位移总刚矩阵K1组集完成后，即可获得整个结构的平衡方程为F2二EF3ll1ll丄11—+—11121I01©1p1U3(3-12)12整个结构的边界条件为U1=0，F2F3已知，三个未知量三个方程，因此上式可求得唯一解解出节点位移:片一F2F3U1=0_5U2=0.762105in_4U3=0.1829510in2•解：（1）建立需要计算的力学模型以及划分单元由于该结构几何对称和受载也对称，故可利用其对称性，只需要取薄板的1/4作为计算对象。为了简单起见，我们把它划分成4个三角形单元，单元和节点编号如图（b）所示。由于对称，节点1，2,4,不可能有水平位移，节点4,5,6不可能有垂直位移，故施加约束如图（b）所示。图两类单元节点编号取总体x,y坐标并确定各节点的坐标值。由图看出，这里只有两类不同的单元，一类单元是1,2,4,另一类单元是3。两类单元节点的编排如图所示。单元1,单元节点编排对应于结构的节点编号1,2,3。三个节点坐标如下:XiXjXm=1m,yi=2m,yj=1m,ym=1m代入得:bi二yj_ym-0；bj二ym一yici-Xm-Xj"；Cj=Xj_Xm-_1；cm=Xj-Xj=°三角形面积：•=-m22单元节点坐标以及单元和节点的编号是原始数据，可用手工输入，也可由计算机完成。对于单元2,3,4定出单元节点的坐标值后，同样可算出，以及各单元的面积。计算个单元的刚度矩阵ke及组集成总刚K由于t=1m，J=0，所以Et__E41-」2厶-2于是由式可求得单元刚度矩阵为k17k；k1ijk1.jjkmjk；mk；mk1k21k12k；3k13k；3■0.250:-0.25-0.25'00.25100.5!■』0_-0.5j00-0.250i0.750.25j-0.5-0.25-0.25-0.50.250.75j0-0.2500-0.5—0「I0.50'0.250'-0.25-0.25'00.25一kmimm=E4同理可得单元2,4的刚度矩阵分别为2425k242k254244k：5255k33|k53_k63k35k：5k465由于1,2,4单元算出的bi,bj…等值以及三角形面积均相同，故算出2,4的单元刚度矩阵与单元1的刚度矩阵数值完全相同单元3的节点i,j,m相应于总体编号中的2,5,3点，其节点坐标为xi=0,Xj=1m,xm=1myi"my二0,ym=1m由此得:biCi--1m,bj=0,bm=1m,=0,Cj--1m,cm=1m.从而算出单元刚度矩阵为:k3二飞3k3k3jk3jjk3kmjk；5k；5k；50.50:0:0-0.50100.25!0.250_-0.25-0.2500.25[0.25I0-0.25-0.2500L-0―…0.50-0.15-0.5-0.25|-0.2500.750.250-0.25'-0.25-0.50.250.75一=EK1为k3^mm根据各单元刚度矩阵组集成总刚度矩阵ki1ik；2k；3k23k23k33-k33k330k2k240k2k44（对称）0k25k35k：5k55'k55'k55k2k25k3k35o【0k360k：k；6一由以上结果求得总刚度矩阵各元素为kn=kAE.|0；5005kl2=鳥=—0.25—0.50.25—0.5k22_k22+总+k22_E0.250.25,厂0.250厂0.500.75E|L00.25：1.50.25]—0.251.5k23_k；3+k；3_E0占-0.25…|-0.5-0.250—1一°.25||-0.25-0.25[-0.25-0.5k35_k24_k；4_Ek25_k25+k25_Ek35+k35_k55_k55+k55+k55_-0.25-0.25IL0-0.500]_J00.25一-一J0.250E严01孑50.25〔孑5°]寸0.25〕0.25_]0.250.75_:00.5_]0.251.5一0005E_0.250k33_k33+k33+k33_E皿11-0.250E®5-0.50-0.25-0.5-0.5一11-0.25-0.25-1k44=kk45=244oo_E°.75110.250.250.75k2_匚「0.5k45=E_o-0.25-0.25E°.5_0OI+EP250.25一]00；e0.25025=e150.75|[0.250.251.5-0.25-0.25一00.25K中去，于是得到K】把上面计算出的kii,…，k66对号入座放到总刚矩阵的具体达式计算并代入等效节点载荷及相应的位移边界条件，以建立和求解未知节点位移的平衡方程组。先求出各项等效节点载荷然后叠加，以形成方程组右端载荷项，但本问题只在节点1有一个集中外载荷Ri=10kN/m(取F=20kN/m解此方程组即得到位置节点位移分量由上方程组求得位移分量如下:解此方程组即得到位置节点位移分量由上方程组求得位移分量如下:的一半）由结构的对称性，可以看出Uj=U2=U4=V4=V5=V6=0。于是需要求的未知节点位移分量只有6个，即Vj，V2，U3，V3，U5。代入边界条件及外载荷以及支反力后，其方程组为■「0.250-0.25-0.2500.25000000■01■F1xl00.50-0.500000000V1-10-0.2501.50.25-1-0.25-0.25-0.2500.25000F2x-0.25-0.50.251.5-0.25-0.50-0.50.25000V2000-1-0.251.50.2500-0.5-0.2500.25U200.250-0.25-0.50.251.500-0.25-100V30E00-0.250000.750.25-0.5-0.25000F4x00-0.25-0.5000.250.750-0.25000F4y0000.25-0.5-0.25-0.501.50.25-0.5-0.25U50000.250-0.25-1-0.25-0.250.251.50-0.250F5y00000000-0.500.50U6000000.25000-0.25-0.2500.250J:F6y把左端系数矩阵行列倒换,于是可分块求解第二种办法是把带有支反力的方程去掉，即把系数矩阵中的第1,3,7,个未知位移的方程式：8，10，12行和列划掉，得出带有6‘0.5-0.500-0.51.5-0.25-0.50-0.251.50.25E0-0.50.251.500.25-0.25-0.25000000、SZ-10A0.250V20-0.50U30-0.250V301.5-0.5U50—0.50.5‘lU6J<0丿2gJ32.52/E、V2—12.52/EU3-0.88/EV3-3.72/EU51.76/E帆、、、1.76/E』求单元应力分量求出节点位移分量后，就可以按式计算单元中的应力我们略去初应变，于是有对于单元1，2，4：-■00-10101b】=E010-1000.50-0.5-0.500.5一1对于单元3:--1000101S=E000-100-1-0-0.5-0.500050.5一注意到Ui=U2=U4=V4=V5=V6=0，最后可求得各单兀的应力为(1)=E-1Vi-0.88-20.00kN/m2-0.5-0.50.5丿4.40丿U3CTx(2)'00-101CTy=E010—10T.xy)©.50-0.5-0.50000.5?V200U3r1.76=-12.52kN/m20」a、:3)r00-101=E010-10cxy」050-0.5-0.50000.5;VU50U3-3.72kN/m(43.08/ZU3缶x'(4)广00-100=E010-1<_xy丿卫.50-0.5-0.5100000.5,V3U50U602-3.72kN/mL・32」如图所示标出了各个单元的应力值,而且在单元内是不变的，这就说明了是一近似解。在单元交界处，应力值有突变，这就可以看出，如将单元分得很细,则突变减小，其结果将会改善计算后的各单元应力3•解：Ni=Li（2Li-1）；N2=L2（2L2-1）;N3=L3（2L3-1）；N4=4L2L2；N5=4L2L3；Ne=4L3Li；在三角形142边上L3=0；可得N3=N5=N6=0，因此有：TOC\o"1-5"\h\z[fxj⑹[fxj⑹[01lfyjO；IfyJ©；gO在三角形142边上建立局部坐标系如图：-1010&-B1,4,2节点在局部坐标系下的形函数为:-21-4=仁「；2=~1所以:TOC\o"1-5"\h\z111仏二一qtId二一qlt;12」16112f4=2Nid^q|t；111f2qt2dqlt;2」6其中，l为三角形142边的边长；fx一I1f一cost1Icost1y^y1tH-qUqtt>.fX4,fy4ify11恥"日J6[Isi*J6[x^xj；=2qt：y2「y」3K—x-tcos芈Zqjcos鋼if4sin「3Jsin日「'f2cos^〕1Icos日〕1'厂>=T2_>=—qt」>=—qt；lf2sin廿j6Jsi*j6[其中：B为q与水平方向的夹角