有限单元法参考答案

1、有限单元试题参考答案一、问答题 50 分1. 5 分有限单元位移法求解弹性力学问题的根本步骤有哪些？1选择适当的单元类型将弹性体离散化2建立单元体的位移插值函数3推导单元刚度矩阵4将单元刚度矩阵组装成整体刚度矩阵5代入边界条件和求解2. 5 分有限元法在单元划分的时候应注意哪些问题？ 1集中载荷的作用点、分布载荷的突变点和约束的支撑点都应取为结点 2在应力变化剧烈的区域，单元划分得细一些，其它应力平缓的区域划分得粗 一些3为了防止在计算中产生过大的误差，单元的长细比最好不要大于23. 5 分有限元法中建立位移函数一般有广义坐标法和插值函数法，我们经常 用插值函数的哪些性质来直接建立位移函数？

2、1形函数与位移插值函数是相同次数的多项式2形函数Ni在结点i处等于1,在其它结点上的值等于03在单元任意一点，三个形函数之和为 14. 10分在有限元法中,单元刚度矩阵和整体刚度矩阵具有哪些性质？1 单元刚度矩阵每一列元素表示一组平衡力系,对于平面问题每列元素之和为 零2单元刚度矩阵对角线元素总为正3单元刚度矩阵为对称矩阵4单元刚度矩阵为奇异矩阵 整体刚度矩阵前三条性质和单元刚度矩阵一样。另外：1) 整体刚度矩阵为奇异矩阵,排除刚体位移后为正定矩阵 2整体刚度矩阵是带状矩阵5. 5 分什么是等参数单元？它与三角形单元和矩形单元相比有哪些优势？ 1在建立局部坐标系下的形状规那么的标准单元与整体坐

3、标系下形状复杂的实际 单元之间的变换时, 如果坐标变换函数中的形函数及插值结点与描述单元位移函 数的形函数及插值结点完全相同, 那么这种变换我们成为等参数变换, 当中的实际 单元单元称为等参数单元。 其它描述意思一样也可 2三角形单元和矩形单元不能适应复杂的曲线边界,等参数单元可以。6. 10 分平面三角形单元与轴对称问题的三角形截面单元的不同之处在哪里？ 轴对称问题三角形截面单元刚度方程的推导当中, 为了简化计算和消除在对称轴 上 r=0 引起的麻烦,可怎样处理？1 平面三角形单元的 三个应力分量xy和三个应变分量Yy都为常量，是常应变单元也是常应力单元轴对称问题的三角形截面单元有四个应变分

4、量££zYz，四个应力分量rz，其中JYzrz是常量,£0是坐标r、z的函数。不同于平面三角形单元的常应变特性和常应力特性。2把单元中随点变化的r,z,近似的视为常数，用单元截面形心处的r,z,来近似 表示，这样将各个单元近似的当成常应变单元。7. 10分在薄板弯曲理论中做了哪些假设？如何用中面位移确定板内任一点的 位移？假设：1变形前后板的厚度不变，即&Z=0b2变形前的中面法线在变形后仍为弹性曲面的法线，即丫zx=0,Y zy=03薄板中面只发生弯曲变形，没有面内的伸缩变形，即中面水平位移。(u)z=O =(v)z=O =04板内各水平层间互不挤压，即

5、6z=0薄板的全部位移、应力和应变分量都可以用板的挠度3来表示，而薄板小挠度弯 曲被简化为中面的弯曲问题，只要中面挠度3确定，任何点的位移都可确定。薄板内不等于零的应变分量有如下三个:?u?x-z?2 3?x2?v£ = - y ?y丫 xy?u?v= + =?y?x-2z?2 3?x ?y20 分1. 10分对于四节点矩形单元我们通常建立如下位移函数，请分析此位移函数 下单元的完备性和协调性。u12X 3y 4xyv 56X 7y 8xy这个位移函数是x， y的双线性函数，其中a 1、a 2、a 3，和a 5、a 6、a 7反映 了刚体位移和常应变，所以这种单元是完备单元。另外，在

6、相邻单元的公共边界x =±a和丫=± b上，位移函数按线性变化，而相 邻单元的边界的各点有相同的位移， 这就保证了相邻单元的协调性，因此这种单 元也是协调单元。2. 10分有限元法在拼装整体刚度矩阵时可用扩充单元刚度方程法和对号入座 法。整体刚度矩阵中非零元素集中分布在主对角元素两侧，呈带状分布，其集中程度与节点编号有关。如以下图所示平面问题的两种结点编号方式，第一种编号 方式对应的整体刚度矩阵非零元素的分布已在左边矩阵中标出，请将第二种编号方式非零元素的分布在右边的矩阵中中标出可用对号入座法。并分析哪种编号方式更好。注：半带宽B=相关节点编号最大差值+1X每个结点的自由度

7、、第一种编号方式：B= 2+1 X 2=6第一种编号方式：B= : 4+1X 2=10为了节省电脑存储空间和计算时间， 应该使半带宽尽量小，所以第一种编号方式 好。三.计算题301.20分图1所示为一个平面应力状态的直角三角形单元，弹性模量为E,泊松比卩二0,厚度为t，试求：1形函数矩阵N2应变矩阵B3单元刚度矩阵ke注：10ED2101100212. A 2 XjymXmyjyjym xXma*4z(a,0)图11Xj yaibjX c,y(i,J,m)1.平面应力状态下解:1N1x 0y0ax-y02 a0 x0y0a x-y2b0bj0bm010001011B0Ci0CjCm0 2000

8、1012Aa2GbCjbjCmbm0110112000-200110-1-13k eEt0110-1-1400020-2-2-1-10310-1-1-213答题时需写出必要的推导过程2.10分图2所示的四节点矩形单元，求出节点 3的位移。设厚度 =0， E为常量注：对于四节点矩形单元有：1N,- 1141N2111.4 Ni1N3 - 11341N4- 114141 i1 i(i 123,4)eT2. kBAD B tdxdyT1 1kijBi D Bjtdxdy abt 1 1Ab12-ij 1ijEtai j3 i j81 22i jTB D Bj d dk11k12k13ek14k21k22k23k24k31k32k33k34k41k42k43k441a1ij 1 一 i jbj3 ja1i j2-b(i,j 1