[工学]第三章 平面问题的有限元法.ppt

第3章 平面问题的有限元法 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.1 平面三角形单元矩阵推导 3.2 利用平面三角形单元进行整体分析 3.2 利用平面三角形单元进行整体分析 3.2 利用平面三角形单元进行整体分析 3.2 利用平面三角形单元进行整体分析 3.2 利用平面三角形单元进行整体分析 3.2 利用平面三角形单元进行整体分析 3.2 利用平面三角形单元进行整体分析 3.2 利用平面三角形单元进行整体分析 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 3.3 平面三角形单元应用举例 例：如图3-4所示是一个薄板，在右上角处受集中载荷作用，底边受到约束 。该平板可以视为弹性力学平面应力问题加以分析。材料参数为E=3.0×105MPa，μ=0.3，厚度t=1。求该平板的应力分布。 根据本章有关内容，有限元求解的具体过程如下： (1) 建立平面直角坐标系oxy，原点为图示节点1处，水平为x轴，垂直为y轴。列写节点坐标码如下： 20 20 20 10 10 10 0 0 0 y坐标值 10 5 0 10 5 0 10 5 0 x坐标值 9 8 7 6 5 4 3 2 1 节点号 8 7 9 8 5 4 6 5 节点3 9 8 6 5 6 5 3 2 节点2 5 4 5 4 2 1 2 1 节点1 8 7 6 5 4 3 2 1 单元号 （2) 计算各单元的单元刚度矩阵并进行扩展。首先计算出所需的系数。根据上节有关公式，对于单元1，三个节点对应的整体编码为, 求得 应变矩阵 弹性力学平面问题的弹性矩阵为 得到单元1的单元刚度矩阵为 单元1的单元刚度矩阵进行扩展，得到一个18×18阶的方阵。只在上述(1, 2, 5)三个节点对应的元素上有值，其它元素上均为0。 全部8 个单元均按上述同样的过程进行计算。 (4) 在考虑位移约束条件的情况下列写结构节点位移列阵。在本示例中，节点1，2，3 处均为全约束，即这三个节点的x, y方向对应的位移分量为0，即 (5) 考虑结构的外载荷，构造结构载荷列阵。在本示例中，只在节点9处作用水平和垂直载荷，因此可以得到 （3）对上述8 个单元的扩展刚度矩阵进行叠加，得到结构的整体刚度矩阵 (6) 根据本节内容引入边界条件，即根据约束情况修正结构有限元方程，特别是消除整体刚度矩阵的奇异性，得到考虑约束条件的、可解的有限元方程。 (7) 利用线性方程组的数值解法，对上述结构的有限元方程进行求解，得到所有各节点的位移向量。最后根据需要求解单元应力。得到的节点位移值和单元应力值可参见如下程序计算出的结果。 Matlab程序结果如下 (1) 节点位移（U={u1, v1, u2, v2, ......, u8, v8}） U = 1.0e-008 * {0 0 0 0 0 0 0.1721 0.1085 0.1611 0.0331 0.1533 -0.0420 0.3904 0.1418 0.3895 0.0561 0.4073 0.0011} (2) 单元应力 单元 1 2 3 4 5 6 7 8 3.2709 -4.1539 3.4894 -1.9224 -4.9644