有限元考试试题及答案

《有限元考试试题及答案》由会员分享，可在线阅读，更多相关《有限元考试试题及答案（16页珍藏版）》请在装配图网上搜索。

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

2、 60 分，每题 20 分）1. 一杆件如图 3 所示，杆件上方固定后，在下方受垂直向下的集中力作用，已 知：杆件材料的杨氏模量Ei = E2 = 3x107lbf /in2 ,截面积Ai = 525沁,A2二3-75in2，长度Li = L2 = 12 ，集中力P = I00吋，用有限元方法求解B点 和C点位移。备注：（1） 1 lbf （磅力，libra force） = 4.45 N。（2）杨氏模量、 弹性模量、Young氏弹性模量具有相同含义（10分）2. 如图2所示，有一正方形薄板，沿对角承受压力作用，厚度t=1m，载荷F=20KN/m，设泊松比p=0,材料的弹性模量为E,试求它的应

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

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

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

6、质力学领域。在工程分析中的作用已从分析和校核扩展到优 化设计并和计算机辅助设计技术相结合。可以预计，随着现代力学、计算数学和 计算机技术等学科的发展，有限单元法作为一个具有巩固理论基础和广泛应用效 力的数值分析工具，必将在国民经济建设和科学技术发展中发挥更大的作用，其 自身亦将得到进一步的发展和完善。三、计算题1.解：将杆件分解成两个元素,A x E元素 1 的刚度矩阵 K= - i 二 13.125 x 106Ibf /in，1_，a x E元素 2 的刚度矩阵 K2=2l 2 = 9.375 x 106Ibf /in2总刚度矩阵单元刚度矩阵形成后，应将各单元刚度矩阵组装集合成整体刚度矩阵（

7、即总刚矩阵）。如所示为杆系结构两单元节点编号示意图，可得总刚度矩阵为k 11K = k叫0k G12k G + k G)22 :k(2)3222k G)23k(2)(3-11)33图3-3杆系结构两单元节点编号示意图引入边界条件求解节点位移总刚矩阵Ik 组集完成后，即可获得整个结构的平衡方程为r f 11 = E2FI 3丿1厂-T10-1-T1 h+ / /1- 12-T2-1-T12厂2u1 U2u3丿(3-12)整个结构的边界条件为U已知，三个未知量三个方程，因此U=1/l 一一 F 一2112U=E/l +1FL 3JL 11 21-3上式可求得唯一解。解出节点位移:-(F + F )

8、23u1=0 u2=0.762x10-5in u3=0.18295x10-4in2. 解:（1）建立需要计算的力学模型以及划分单元由于该结构几何对称和受载也对称，故可利用其对称性，只需要取薄板的1/4作为计算对象。为了简单起见，我们把它划分成4个三角形单元，单元和节点编号如图（b）所示。由于对称，节点1，2, 4,不可能有水平位移，节点4,5，6不可能有垂直位移，故施加约束如图（b）所示。图两类单元节点编号取总体x, y坐标并确定各节点的坐标值。由图看出，这里只有两类不同的单 元，一类单元是1，2, 4,另一类单元是3。两类单元节点的编排如图所示。单元1,单元节点编排对应于结构的节点编号1,

9、2, 3。三个节点坐标如下:x 0， ix 0，jx 1m, m代入得：yjxm三角形面积:y 1m ， jy -ymi二-1 ;=1m y -yij- x 1；jx -xim二-1;-x 0iA 丄 m 22同样可算出，以及各单元节点坐标以及单元和节点的编号是原始数据，可用手工输入，也可由计 算机完成。对于单元 2，3，4 定出单元节点 的坐标值后， 单元的面积。(2) 计算个单元的刚度矩阵ke及组集成总刚K由于t 1m ,卩=0，所以Et于是由式可求得单元刚度矩阵为_ k 1k1k 1_ k 1k1k 1iiijim111213k 1k1k 1k 1k1k 1jijjjm212223k 1

10、k1k 1k 1k1k 1mimjmm3132330.250；0.250.25 ；00.2500.500.5000.2500.750.250.50.250.250.50.250.7500.2500| - 0.50 |0.500.250! 0.250.25 !00.25=E同理可得单元 2，4 的刚度矩阵分别为k 2k2k 2 k 4k4k 4222425333536ek 2k2k 2, ke k 4k4k 424244454535556k 2k2k 2k 4k4k 4525455636566由于1, 2, 4单元算出的b ,b等值以及三角形面积均相同，故算出2, 4 ij的单元刚度矩阵与单元

11、1 的刚度矩阵数值完全相同。单元3的节点i, j,m相应于总体编号中的2, 5, 3点，其节点坐标为x = 0, x = 1 m, x = 1 mijmy =1m,y =0,yijm=1m由此得：b = -1m, b = 0, b ic = 0, cijm= -1 m, c= 1 m,= 1 m.从而算出单元刚度矩阵为：k 3k 3k3 1k 3k3k3 1iijim222523k e =3k 3k 3jijjk3k33k3jm3k352k3k3k35553k3k3mimj3k3mm3235330.50 |00j - 0.50 _00.25 !0.250-0.25 -0.2500.25 :0.

12、250-0.25 -0.250000.50 -0.15-0.5-0.25 i-0.250i 0.750.250-0.25 丨0.25-0.5! 0.250.75m=E根据各单元刚度矩阵组集成总刚度矩阵k 为k1k1k100111213k 1 + k 2+ k 3k i + k 3k2k2+ k 322 22222323242525k i + k 3 + k 40k3+ k 43333333535k2k24445(K =对 称）k255+ k55+ k 455010k 4360k 456k 466由以上结果求得总刚度矩阵各元素为k =k i =E11 110.250.5k =k i =E12 1

13、20.250.250.5k =k i =E13130.2510.50.750.25+E0.250 -+ E0.50 -二 E 1.50.250.250.7500.500.250.251.5k =k i +k 2 +k 3 =E22 22 22 22240.50 一+ E0.750.25_+ E0.250 一二 E 1.50.25_00.25_0.250.7500.50.251.51111-0.250 一-0.25-0.25 - 0.5-0.25_k =k3+k4=E+ E二 E353535_- 0.25-0.5_ 0-0.5 _- 0.25-1k =k i +k 3 +k 4 =E333333

14、33k =k4=E36360.250k =k2=E44440.750.250.250.75k =k2=E4545- 0.50- 0.25- 0.25k =k2+k3+k4=E555555550 一+ E0.250 一+ E0.750.25二 E1.50.25_0.2500.5_0.250.750.251.50.50-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 =ki +k3=E232323=k2 =E00.25_+ E00_二 E00.25_00_0.2500.250- 0.5024k

15、=k2+k3=E252525=k4=E5656-0.5_ 0-0.25-0.250.50 _k =k4=E66 6600.25把上面计算出的k，k对号入座放到总刚矩阵bl中去，于是得到Ikii66的具体表达式。（3）计算并代入等效节点载荷及相应的位移边界条件，以建立和求解未知 节点位移的平衡方程组。先求出各项等效节点载荷然后叠加，以形成方程组右端载荷项，但本问题只在节点1有一个集中外载荷R = 10kN/m（取F二20kN/m i的一半）。由结构的对称性，可以看出u二u12=v =0。于是需要求的未6知节点位移分量只有 6 个，即 v ， v ， u ， v ， u1 2 3 3 5代入边界条

16、件及外载荷以及支反力后，其方程组为0.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.25- 0.5- 0.250000- 0.25- 0.5000.250.750- 0.25000000.25- 0.5- 0.25- 0

17、.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.25Ev10v2 u2 v300u50u60F1x-10F2 x000F4 xF4y0F5 y0F6y把左端系数矩阵行列倒换，于是可分块求解。第二种办法是把带有支反力的方程去掉，即把系数矩阵中的第 1， 3， 7 个未知位移的方程式：8， 10， 12 行和列划掉，得出带有 60.50.50-0.51.5- 0 . 250- 0 . 251.50-0.50.25E0-0.50.251.

18、500.25- 0 . 25- 0 . 25 000000、(v )1vJ10、0.25002-0.50u3v0- 0 . 250031.50.5u05-0.50.5丿Iu 6 J0丿解此方程组即得到位置节点位移分量由上方程组求得位移分量如下：u3V3u5U )600-10ts= E 010-10.5 0 -0.5 -0.510000 0.532.52/E、-12.52/E-0.88/E-3.72/E1.76/E 1.76/E ,(4) 求单元应力分量求出节点位移分量后，就可以按式计算单元中的应力 我们略去初应变8 0，于是有对于单元 1，2，4：最后可求得各单元的应力为对于单元 3：0010

19、0-100-1-0.500050.5-10k= E 000- 0.5注意至U u = u = u = v = v = v = 0 ,厂00-1x0=)E010y,T、0.50-0.5xyV0-1-0.50、00.5丿124456-0.88 -20.00 kN/m2、4.40 丿v3丿G 1厂00-10x0=)E010-1yVT丿xyV0.50-0.5-0.5xy1000、00.5丿u3(1.76 - 12.52 kN / m 2(00-10x0=)E010-1yT、0.50-0.5-0.5xy丿1000、00.5丿V3丿(-0.88 -3.72 kN / m2,43.08 丿(o 1(00-1

20、0x0=)E010-1y,TJ xy丿,0.50-0.5-0.5xy1000、00.5丿(0 -3.72 kN / m 2，-1，32 丿u6如图所示标出了各个单元的应力值，而且在单元内是不变的，这就说明了是 一近似解。在单元交界处，应力值有突变，这就可以看出，如将单元分得很细,则突变减小，其结果将会改善。计算后的各单元应力3. 解：N1=L1(2L1-1)； N2=L2(2L2-1)； N3=L3(2L3-1)； N4=4L2L2； N5=4L2L3； N6=4L3L1 ； 在三角形142边上L3=0；可得N3=N5=N6=0，因此有：f0f0f厂、0 =vvX5 =v =vfvy丿0V丿；

21、fJ y 丿0J 丿fJ y 丿0k 丿在三角形142边上建立局部坐标系如图:-1 0 11，4，2 节点在局部坐标系下的形函数为：叫=-2 C-映；忙（1一次+次餐所以：f = 11 qt弔 dg = 1 qlt ；1 2 一61 1 2 f = J1 qt dg = qlt ；42_143111f = J1 qt dg = qlt ；其中，2 2一126l 为三角形 142 边的边长；x1x4f cos01l cos 01yf sin0 =qt s6l sin0 y =一 qt3l sin 04qt X2f cos 02yf sin 021=一 qt s6l cos 0l sin 0V丿y 2 - y12 1卜x X21y 2 y121 yX X21y 2 y121 yX X21其中：0为q与水平方向的夹角。