平面结构问题的有限单元法.pptxVIP

1 4.1 平面应力问题 第四章 平面结构问题的有限单元法 4.2 平面应变问题 4.3 平面问题的离散化 4.4 平面三结点三角形单元 4.5 ANSYS平面结构计算示例 2 严格地说，任何弹性体都是处于三维受力状态，因而都是空间问题，但是在一定条件下，许多空间问 题都可以简化成平面问题。 平面问题可以分为两类：平面应力问题和平面应变问题。 图4-1 平面问题应力状态 第四章 平面结构问题的有限单元法 3 4.1 平面应力问题 图4-2(a) 平面应力问题 如图所示的深梁结构，其厚度方向的尺寸远比其它两个方向的尺寸小得多，可视为一薄板。它只承受作用在其平面内的载荷，且沿厚度方向不变，计算时以中性面为研究对象。其力学特点是： 平面应力问题的应力应变转换矩阵即弹性矩阵为： 。 第四章 平面结构问题的有限单元法 4 图4-2(b) 平面应变问题 4.2 平面应变问题 图示为一圆形涵洞的横截面。其长度方向上的尺寸远比其它两个方向上的尺寸大得多，同样，载荷作用在xy坐标面内，且沿z轴方向均 匀分布。其力学特点是： 但一般情况下 平面应变问题的弹性矩阵只需将式(4-1)中的E换成 换成 ， 即可。 。 第四章 平面结构问题的有限单元法 5 无论是平面应力问题还是平面应变问题的应力 与 应变 之间的关系均为： ，其中： 为初应变。 式中 4.3 平面问题的离散化 (a) 三结点三角形单元 (b) 四结点正方形单元 (c) 四结点矩形单元 (d) 四结点四边形单元 图4-3 平面问题单元的主要类型 第四章 平面结构问题的有限单元法 6 图4-4(a)表示的是带有椭圆孔的平板，在均匀压力作用下的应力集中问题。图4-5(b)是利用结构的对称性，采用三结点三角形单元而离散后的力学模型，各单元之间以结点相连。 (a) 均匀受力板力学模型 (b) 力学模型离散化 图4-4 平面问题有限单元法的计算力学模型 第四章 平面结构问题的有限单元法 7 4.4 平面三结点三角形单元 4.1.1 位移函数 图4-5 三角形单元 如果把弹性体离散成为有限个单元体，而且单元很小时，就很容易利用其结点的位移，构造出单元的位移插值函数，即位移函数。 位移函数矩阵形式： 第四章 平面结构问题的有限单元法 8 简写为： 由于位移函数适用于单元中的任意一点，所以带入3个结点的坐标后，得出结点处位移函数为 简写为： 第四章 平面结构问题的有限单元法 9 解出 其中， 是三角形单元的面积，当三角形单元结点i、j、m按逆时针次序排列时，则有 4.4.2 形函数矩阵 第四章 平面结构问题的有限单元法 10 其中记号 表示将i、j、m进行轮换后，可得出另外两组带脚标的a、b、c的公式。 单元位移函数为结点位移的插值函数，即 第四章 平面结构问题的有限单元法 (4-9) 11 令 在式(4-10)中表示的 称为形函数，于是位移函数表达式用形函数表示为： (4-10) (4-11) 写成矩阵形式 (4-12) 第四章 平面结构问题的有限单元法 12 由几何方程知 将式(4-9)代入式(4-13)中，并求偏导数，得 (4-13) 4.4.3 单元的应力与应变 第四章 平面结构问题的有限单元法 13 简写为： (4-14) 由于[B]是常量，单元内各点应变分量也都是常量，这是由于采用了线性位移函数的缘故，这种单元称为常应变三角形单元。 (4-15) 第四章 平面结构问题的有限单元法 14 由弹性力学的物理方程可知，其应力与应变有如下关系： (4-16) 将式(4-14)代入式(4-16)，得 (4-17) 式中 (4-18) [S]称为应力转换矩阵，对平面应力问题，其子矩阵为 (4-19) 由式(4-17)看出，应力分量也是一个常量。在一个三角形单元中各点应力相同，一般用形心一点表示。其应变也可同样表示。 第四章 平面结构问题的有限单元法 15 用虚功原理来建立结点力和结点位移间的关系式，从而得