有限元习习题

1、判断题： 1. 象床单那样薄、那样宽的板用梁单元来模型化 2. 对于高压电线的铁塔那样的框架结构的模型化处理使用梁单元3. 一般自由度多的模型分析成本高4. 使用尽可能多种类单元的模型是一个好的模型5. 杆单元是壳单元的一种6. 不能把梁单元、壳单元和实体单元混合在一起作成模型7. 四边形的壳单元尽可能作成接近正方形形状的单元8. 因为实体单元是3维单元，所以即使有严重的扭曲也没关系9. 将作用有垂直载荷的悬臂梁用多个杆单元作成10. 将作用有垂直载荷的两端自由支持的梁用杆单元来模型化11. 三角形单元和四边形单元不能混在一起使用12. 平面应变单元也好，平面应力单元也好，如果以单位厚来作模型

2、化处理的话会得到一样的答案13. 同样形状的话，使用三角形单元和使用四边形单元解是相同的14. 边长为10cm和边长为100cm 的正方形的板，后者的单元数如果是前者的10倍的话，才行 15. 为了校核连续的相同管子剖面内的应力状态，要使用平面应力单元 16. 对热应力问题，1维单元也好2维单元也好，所求的解都搞不清 17. 对于热传导分析必须输入线膨胀系数 18. 热应力随结构的约束状态而变化 19. FEM分析变形越大应力就越高 20. 在线性分析中，即使变形变大，如果将这部分单元划分得多一些的话，也会保证解的适当正确 21. 为了评价应力集中，在网格划分时应该把整个作成一样的单元尺寸 2

3、2. 板厚并不一致的情况下，一定要用到实体单元 23. 单元数相同的话，1阶单元、2阶单元的解都一样 24. 为了忠实地尽可能表现结构的形状，必须严格按装配顺序来做模型化处理 25. 节点的位置依赖于形态，而并不依赖于载荷的位置 26.一般应力变化大的地方单元尺寸要划的小才好 27.仅用TETRA单元的模型与仅用HEXA单元的模型相比，后者的精度要好 28. 相接的单元尺寸大小不要变化太厉害 29. 在进行特征值分析时，必须输入质量 30. 进行热应力分析时，必须输入线膨胀系数 31. 壳单元表面的应力因为与表面内的应力相比精度会降低所以必须注意 32. 象船和火箭那样的结构因为漂浮在水（空）

4、中而没被固定住，所以，FEM分析不可以使用 33. 约束条件用全固定或许加上铰固定就能表现完全 34. 一般在特征值分析中一定是采用节点编号连续来编的方法，所得精度要高 35. 用固有振动分析求应力，应力高的部分必须要加强 36. 屈曲模态并不依赖于约束条件 37. 自由度有位移自由度和转角自由度 38. 一般在FEM中使用的模型称为刚体模型 39. 对比铁更硬的部分所做模型化处理的单元称为刚体单元 40. 刚体单元和梁单元和板单元组合在一起进行分析是不可以的 41. 一般网格划分过度的话，很费分析时间 42. 对啤酒罐的压缩强度要用固有振动分析来评价 43. 表示自由度的坐标系有局部坐标系和

5、整体坐标系 44. 应力集中的部分是多个载荷所加的部位 45. 在加上热载的情况下，即使是同一个模型，根据约束条件，所发生的应力有很大的不同 46. 用有限元法可以对正在动的（移动）物体的结构进行分析 47. 对膜（membran）单元也可用面压载荷 48. 可对膜（membran）单元可以用集中载荷 49. 施加强迫位移的分析要进行静力分析 50. 一般所给出的载荷的总和与反力的总和相一致 51. 即使将不同的局部坐标系下定义好的节点连起来也可定义单元 52. 所谓自由度是直接翻译degrees of freedom 的 53. 所谓实体单元意味着刚体单元的集合 54. 杨氏系数是纵弹性系数

6、（模量） 55. 共鸣现象与固有频率有关 56. 杨氏系数是评价材令的基值 57. 即使是同一种材料，梁单元和板单元也要输入不同的材料性质数值 58. 泊松比是在纵向加压时发生在纵向的应变和横向的应变的比率 59. 用弹性材料可表现塑性化现象 60. 一般线膨胀系数是作为材料常数之一输入 61. 一般用FEM模型化时，大的结构求得的热变形小 62. 约束条件全都没被定义的结构不能分析 63. X、Y、Z全部方向上的位移都是1时称为刚体变形 64. 分析结果是对称的模型，使用对称条件可以用较少的单元来进行分析 65. 所谓铰约束条件是约束位移自由度而让转角自由度自由 66. 强迫位移是一种约束条

7、件 67. 即使所有的自由度都约束也会发生变形 68. 对于设置了约束的自由度即使输入载荷也不发生位移 69. 有限单元分析约束条件尽量少则精度好 70. 所谓约束就是消去自由度 71. 所谓全约束只要将位移自由度约束住 72. 壳单元与实体单元可约束的自由度不同 73. 线性分析将同样大的载荷加在反向产生位移的绝对值不变 74. 由分析所得的最大应力受网格划分的影响 75. 载荷和应力表示同一件东西 76. 主应力并不依赖于基本坐标系 77. 在应力分析中，应力小的部位单元尺寸要小，大的部位单元尺寸要大来进行模型化处理 78. 实特征值分析是一种求最大应力的手段 79. 具有切口附近的应力集

8、中用FEM不能严密地计算 80. 1阶单元是假定单元内的应力都一样的单元 81. 表现材料的弹性界限是所谓的屈服应力 82. 在屈服曲面内材料表现为弹性行为 83. 位移能用6个矢量成分来表示 84. 转角是一种位移 85. 载荷点的位移通常最大 86. 线性应力分析也可以得到极大的变形 87. 与材料无关的相同变形量产生相同的应力 88. 给出同一载荷杨氏系数越大则变形也越大 89. 对于静力分析质量是不可缺少的数据 90. 实特征值分析中必须定义集中载荷或分布载荷 91. 屈曲分析和固有振动分析是类似的特征值问题 92. 使用同一模型时，一般特征值分析要比线弹性分析化时间 93. 一般求特

9、征值分析所求的模态数多也好少也好，分析时间是一样的 94. 在静力分析中，仅施加左右方向的载荷时，不约束上下方向也可以 95. 卡车通过时，玻璃窗会别别地振动，这是与玻璃的固有频率有关 96. FEM也被用在医学上 97. 有限元法、有限体积法、有限差分法、边界元法这中间FEM是有限差分法 98. 有限元法基本的是求解联立方程式 99. FEM理论1950年前开始就有了 100. 考虑阻尼的特征值问题成了复特征值问题 1.如图所示为一平面应力状态的直角三角形单元，设=1/6。试求：（1）形函数矩阵N；（2）几何矩阵B；（3）应力矩阵S；（4）单元刚度矩阵Ke。xym(0,0)i(a,0)j(0

10、,b)ba解:（1） 形函数矩阵N将i,j,m的坐标代入得九个系数: 三角形面积: 所以, （2）几何矩阵B；（3）应力矩阵S；当=1/6时（4）单元刚度矩阵Ke。2试证明单元刚度矩阵任意一行或一列的所有元素之和等于零。3试证明 ，(i，j，m)所定义的形函数Ni、Nj、Nm具有以下性质： 4如图所示某单元节点编号及节点坐标，写出其形函数Ni、Nj、Nm及单元应变矩阵Bxym(5,6)i(2,2)j(6,3)4N/mO5N/m4题5题图)1 5上题中，若单元边上作用水平方向线性分布载荷，如图所示，求等效节点载荷。答案：6若单元jm边上作用水平方向线性分布载荷q(s)，如图所示，求等效节点载荷。

11、xymijOq(s)ls123456791081112123456791081112123111249581076(a)1(b)1(c)16题图)17题图)17从带宽最小这个角度出发，如图所示三种节点编号那种最好。8已知矩形单元边长2a2b，坐标原点取在形心，设双线性位移模式，导出四个形函数。xyObbaaxyObaijmp8题图9题图9已知正方形单元边长a，坐标原点取在形心，设双线性位移模式，取=0.2，单元厚度1，导出平面应变问题的单元刚度矩阵。答案：对称10已知悬臂深梁如图所示，梁右端面作用均布拉力，其合力为P，采用如图网格，设=1/3.厚度t，求节点位移。1m2mp(a)1(0,0)P

12、/2(b)4(0,1)3(2,1)2(2,0)P/2解：单元，i、j、m分别为1、2、3，坐标如图 ， ，=1单元，i、j、m分别为1、3、4，坐标如图 ， ，=1总体刚度矩阵：对称总体刚度方程：K=F其中=u1,v1, u2,v2, u3,v3, u4,v4T ，另外，位移边界条件：u1=v1= u4=v4=0F=U1,V1, U2,V2, U3,V3, U4,V4T，另外，有节点载荷：U 2= U 3=P/2求解化简后总体刚度方程，得节点位移：11已知悬臂深梁如图所示，梁右端面作用均布剪力，其合力为P，采用如图网格，设=1/3.厚度t，求节点位移。1m2mP(a)1(0,0)P/2(b)4

13、(0,1)3(2,1)2(2,0)P/2答案：12已知深梁两端固定如图所示，梁中间面作用集中力P，采用如图网格，根据对称性取一半结构计算，设=1/6.厚度t=1，求2、3节点位移。1m1mP(a)1(0,0)P(b)4(0,1)2(1,0)1m3(1,1)答案：v2= -0.966P/E，v3= -1.368P/E13已知正方形薄板，边长为a，厚度为t，其左上端受集中力P作用，如图所示，设材料弹性模量E，泊松比=0.25，用有限元法求力P作用点位移。P(a)3(0,0)1(0,a)4(a,0)a2(a,a)aP(b)答案：v1=1.36P/Et，u1= 014三角形截面水坝如图所示，受静水压力作用，取厚度t=1，材料弹性模量E，泊松比=1/3，采用如图网格，求各节点位移，及单元应力。1m1m14题图)1yx1m1m32654q0=20kN/m引入约束后的总刚方程：15已知悬臂深梁如图所示，设泊松比=1/3，厚度t，弹性模量为E，求节点温度改变为T时节点位移和单元应力。1m2mP(a)1(0,0)P/2(b)4(0,1)3(2,1)2(2,0)P/2解：单元，i、j、m分别为1、2、3，单元，i、j、m分别为1、3、4，利用前面问题结果，总体刚度矩阵：总体刚度方程：K=F其中=u1,v1, u2,v2, u3,v3, u4,