等参数单元.pptVIP

第6章 等参数单元 6.1 等参元的基本概念 6.1 等参元的基本概念 6.1 等参元的基本概念 6.1 等参元的基本概念 6.1 等参元的基本概念 6.1 等参元的基本概念 6.1 等参元的基本概念 6.1 等参元的基本概念 6.1 等参元的基本概念 6.1 等参元的基本概念 6.2 等参元的单元分析 6.2 等参元的单元分析 6.2 等参元的单元分析 6.2 等参元的单元分析 6.2 等参元的单元分析 6.3 平面三角形单元的等参元刚度矩阵 6.3 平面三角形单元的等参元刚度矩阵 6.3 平面三角形单元的等参元刚度矩阵 6.3 平面三角形单元的等参元刚度矩阵 6.4 矩形等参元单元的有限元分析 6.4 矩形等参元单元的有限元分析 6.4 矩形等参元的有限元分析 6.4 矩形等参元的有限元分析 6.4 矩形等参单元的有限元分析 6.4 矩形等参单元的有限元分析 6.4 矩形等参单元的有限元分析 6.4 矩形等参单元的有限元分析 6.4 矩形等参单元的有限元分析 6.5 20节点三维空间等参元的单元刚度矩阵分析 6.5 20节点三维空间等参元的单元刚度矩阵分析 6.5 20节点三维空间等参元的单元刚度矩阵分析 6.5 20节点三维空间等参元的单元刚度矩阵分析 6.5 20节点三维空间等参元的单元刚度矩阵分析 6.5 20节点三维空间等参元的单元刚度矩阵分析 6.5 20节点三维空间等参元的单元刚度矩阵分析 6.6 高斯积分法简介 6.6 高斯积分法简介 6.7 应用举例 6.7 应用举例 6.7 应用举例 6.7 应用举例 6.7 应用举例 6.7 应用举例 6.7 应用举例 6.7 应用举例 6.7 应用举例 根据弹性力学几何方程，可以得到单元应变列阵为 (6.52) (6.51) 坐标变换关系可表示为 (6.53) 其中，B是单元的应变矩阵，其分块形式如下 上式中的形函数Ni是局部坐标的函数。在对整体坐标求导时，类似于平面问题，根据复合函数求导数的规则，有以下用雅可比矩阵表达的关系式 (6.54) 其中，三维雅可比矩阵J为 (6.55) 式中每个元素可以分别按如下公式求得 (6.56) 将单元应变代入空间问题的物理方程式，得到单元的应力 (6.57) 式中， S为应力矩阵 （i=1, 2, …, 20） (6.58) 进而利用虚功原理，推导得到单元刚度矩阵为 (6.59) (6.60) 其中 在等参元的单元刚度积分过程中，矩阵中的每个元素都很复杂，必须采用数值积分方法加以完成，如高斯积分法。在这里对高斯积分法的基本原理进行介绍。 高斯积分法是计算复杂的定积分时通常采用的一种数值方法。对于一维定积分问题 (6.61) 可以近似的化为加权为题。所谓数值积分就是把定积分问题近似地化为加权求和问题，在积分区间选定某些点(称为积分点)，求出积分点处的函数值，然后再乘上与这些积分点相对应的求积系数（又称加权系数），再求和，所得的结果认为是被积函数的近似积分值。这种求积方法表达如下 (6.62) 高斯积分法采用以上这种格式，其中积分点坐标 及其对应的加权 可以构造出二维和 系数Hi如 表6-1所示。逐次利用一维高斯求积公式 三维高斯求积公式 高斯积分的阶数通常根据等参元的维数和节点数来选取。例如，平面4节点等参元可取2阶，平面8节点等参元可取3阶，空间8节点等参元可取2阶，而空间20节点等参元可取3阶。 表6-1高斯积分法中的积分点坐标和加权系数 积分点坐标 5688889 0.2369269 0.4786287 0.0000000 ±0.9061798 ±0.5384693 5 3478548 0.6521452 ±0.8611363 ±0.3399810 4 8888889 0.5555556 0000000 ±0.7745967 3 1.0000000 ±0.5773503 2 加权系数Hi 积分点数n (6.63) (6.64) 现有一个受均匀分布载荷作用的薄板结构，几何尺寸及受力情况如图6-6(a) 所示，本节将采用4节点矩形等参元对其求解，读者应掌握具体的求解过程。设该结构的弹性模量E=210Gpa，泊松比 板厚度t=0.025m，均布载荷 试对该结构进行静力学分析。 (b) 图6-6 薄板结构力学模型及有限元模型 为了方便说明，我们将平板仅仅分解为两个单元，6个节点，如图6-6(b)所示。分布载荷的总作用力平均分给节点3和节点6。每个节点力为 求解过程如下: (1)对结构进行离散化 KN 节点坐标见表6-2 表