平面有限元法作业

《平面有限元法作业》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第三章作业3-1：试证明平面三角形单元内任一点的形函数之和恒等于1。证明1：设单元发生X方向的刚体位移，则单元内到处应有位移，有若位移函数不满足此要求，则不能反映单元的刚体位移，不能得到正确的结果。#证明2：设P是三角形内任一点，可用面积坐标表示为。由面积坐标的定义和性质知，且三节点三角形的一点的面积坐标即为其形函数，故平面三角形单元内任一点的形函数之和恒等于1。#3-2：试证明三角形单元的任一边上的一点的三个形函数与第三个顶点的坐标无关。证明1：设k是三角形ij边上的任一点，点k面积坐标得#证明2：三角形单元是协调单元，必须在单元边界上保持连续性，所以在单元边界上的点的位移只能由边上两个节点

2、的形函数来贡献，否则就会撕裂和重叠，即（如在ij边上的点）故三角形的三边上的点的形函数只与边上节点的坐标有关，而与第三点无关。#3-3：证明三角形单元是常应变单元。证明：，10#即三角形单元是常应变单元。3-4：已知单元刚度矩阵，试说明分别是什么矩阵，与单元的那些特性有关？若厚度为t的平面三角形常应变单元ijm的单元刚度矩阵记为：说明子块的物理意义，并证明为对称矩阵。解：是应变矩阵又称几何矩阵，与单元节点坐标有关；为弹性矩阵，与材料的弹性常数E、有关。表示当节点j处产生单位位移，其余节点完全被约束时，在节点i处引起的节点力。利用矩阵的运算关系由于是对称矩阵，所以，即为对称矩阵。#3-5：图示平

3、面等腰三角形单元，若，弹性模量为E，厚度为t，求形函数矩阵、应变矩阵及单元刚度矩阵。（补充题意：平面应力情况）xy解：对平面等腰直角三角形建立图示坐标系。10，形函数形函数矩阵：求应变矩阵，则应力矩阵三角形单元的刚度矩阵10等腰直角三角形的单元刚度矩阵与三角形的面积和节点坐标无关，请同学们记住这个结论，解题时会方便。等边三角形的单元刚度矩阵也有此性质（自行推导）。代入已知数据得#3-6：验证矩形单元的位移模式是否满足位移连续性条件。解：矩形单元如图示，矩形单元的位移模式取为：在单元的边界及上，位移是按线性变化的，而在公共边界上有两个节点相连，这两个公共节点有共同的节点位移值，从而保证了两个相邻

4、单元在其公共边界上位移的连续性。故四节点矩形单元满足位移连续性条件。#3-7：求以下受力单元的等效节点载荷。已知：、q、P，厚度t，P点作用在jm中点处，沿x方向，三角形分布载荷垂直于ij边。10解：q的单元，设厚度为t，如图示，等效节点载荷#3-8：如图a,b,c所示的半带宽各是多少？从带宽优化的角度出发，那种节点编号最好？考虑带宽优化即要求半带宽尽可能的小，故a图节点编号最好。#3-9：写出图3-8题a图网格剖分方案用单元矩阵子块组集成的整体刚度矩阵，标出整体刚度矩阵的阶数。子块编号如下图a，解：10结构离散为12个单元，12个节点，故总体刚度矩阵的阶数为24×24。用单元矩阵子块表示的整

5、体刚度矩阵为：#3-10：图示的正方形薄板，在对角线顶点作用有沿厚度均匀分布的载荷，其合力为2kN，板厚为1mm，为简单起见令，1．根据结构特点和受力特点，确定该结构是平面应力问题还是平面应变问题？2.画出平板的有限元计算模型（包括单元类型选择、网格划分、单元、节点编号、载荷和约束的处理等）；3.写出上述有限元计算模型的节点载荷向量和节点位移向量；4.按对称性，画出简化的有限元模型。写出由单元刚度矩阵子块组集而成的整体刚度矩阵，并确定整体刚度矩阵的阶数；5.写出整体刚度方程，求解方板的变形。解：1）根据薄板的结构特点与受力情况，确定该问题属于平面应力问题。2）对平板用4个单元，5个节点进行结构

6、离散，结构离散及约束和载荷的情况见有限元计算模型如右图所示。3）图示有限元模型的节点载荷向量和节点位移向量为：104）按对称性，简化的有限元计算模型如右图所示。由单元刚度矩阵子块组集而成的整体刚度矩阵：5）单元节点编码i，j，m如果按上图，则各个单元刚度矩阵相同，等腰直角三角形的单元刚度矩阵为（题3-5结果）10总体刚度矩阵边界条件：上式利用降阶后解得10本题没有验证，仅供参考。#3-11题意（略）解：单元节点编码i，j，m如果按上图，则各个单元刚度矩阵相同，等腰直角三角形的单元刚度矩阵为（题3-5结果）10总体刚度矩阵边界条件：，本题只计算一遍，仅供参考。#10