有限元法 - ok.ppt

应变的矩阵： 应变矩阵的分块矩阵： 由于等参单元的形函数式局部坐标（ξ，η）的函数，因此应变距阵[B]也是局部坐标（ξ，η）的函数，形成邓粲单元的刚度矩阵需要在整体坐标系中对局部坐标的函数进行积分，包括以下步骤 ： 由于局部坐标与整体坐标之间存在坐标转换关系，因此形函数Ni 是局部坐标的函数，同时也可以看作是整体坐标的函数。由复合函数求导法则可得： 令： （1）计算形函数对整体坐标x，y的偏导数 矩阵[J]称为雅可比矩阵（Jacobian Matrix），单元的整体坐标可以形函数来表示，因此用坐标变换公式可以计算雅可比矩阵。 ， （2）将整体坐标系中的面积积分转换为在局部坐标系中的面积积分 在整体坐标系中，面积微元为x方向和y方向微矢量的叉乘的模量， 单元刚度矩阵在局部坐标系中的积分公式为： 单元刚度矩阵中的任意一个分块矩阵的积分公式为， （3）用数值积分计算出单元刚度矩阵中的元素 等参单元刚度矩阵的每个元素都是局部坐标的函数，在有限元程序中不用解析的办法来计算局部坐标系中的积分，而采用数值积分方法。通常采用高斯积分方法计算单元刚度矩阵中的元素。 高斯积分 高斯积分方法预先定义了积分点和相应的加权系数，求出被积分的函数在指定积分点上的数值，加权后求和，就得到了该函数的积分。这种方法具有比较高的计算精度。已经证明，采用n个