第五章等参数单元.pptVIP

其中：[J]称为雅可比（Jacobi）矩阵 式（8-17）表示的是由和推导，的变换式，其逆变换式为 等参数单元 返回 第二十九页，共七十四页，2022年，8月28日 等参数单元 其中，[J]-1是[J]的逆阵 返回 第三十页，共七十四页，2022年，8月28日 第三节 平面等参元 等参数单元 　　平面问题的常用等参元有四结点四边形单元、八结点曲边四边形单元和6~8可变结点曲边四边形单元等，本节以八结点曲边四边形等参元为例介绍平面问题分析过程。 一、母单元 　　八结点曲边四边形等参元的母单元是二维二次单元。八个结点分别为正方形的四个角点和四个边中点，母单元采用直角坐标系（ξ,η）。如图5-7所示。单元的位移模式为 其中，ui和vi是结点i的位移。 返回 第三十一页，共七十四页，2022年，8月28日 图 5-8 等 参 元 1 2 3 4 8 7 5 6 2 3 1 5 4 6 7 8 等参数单元 图 5-7 母 单 元 返回 第三十二页，共七十四页，2022年，8月28日 等参数单元 二、等参元 等参元的整体坐标为直角坐标（x, y）。等参元的任意指定的八个结点的整体坐标值分别为（xi, yi），（i=1, 2 ,…, 8），如图5-8所示。 采用坐标变换可使母单元的八个结点 与等参元的八个结点（xi, yi）一一对应。整体坐标 和局部坐标 的变换式为 (8-34) 其中： 是母单元的形函数。 根据等参元的思想，等参元的位移模式仍取为： 返回 第三十三页，共七十四页，2022年，8月28日 等参数单元 (8-35) 　　这样，就确定了平面八结点曲边四边形等参元的几何形状和位移模式。在实际应用中需要注意以下几个问题： 1）在划分单元时，只需确定单元结点的整体坐标值，而不必画出其抛物线形状的边界。因为在计算中实际使用的只有单元八个结点在整体坐标下的位置坐标（xi, yi）（i=1,2,…,8）。 2）在划分单元和布置结点时，单元的各边长度相差不能太大；各边上结点间距应尽量均匀，以减少计算误差。 3）为了计算简单，当求解区域为曲线边界时，只将位于边界的单元取为曲边四边形，而内部单元仍然划分为直边四边形。这样，即能较好地处理曲边边界，又能提高单元内部插值的精度。 第三十四页，共七十四页，2022年，8月28日 等参数单元 三、单元分析 　　将八结点曲边四边形等参元的位移模式代入平面问题的几何方程，便得到单元应变分量的计算式 其中：是单元的结点位移列阵 ( i=1, 2, …, 8) 返回 第三十五页，共七十四页，2022年，8月28日 等参数单元 是单元应变矩阵 ( i=1, 2, …, 8) 由于，形函数是局部坐标的函数。因此，需要进行偏导数的变换 返回 第三十六页，共七十四页，2022年，8月28日 等参数单元 其中，由式（8-20）给出 根据坐标变换式可知 返回 第三十七页，共七十四页，2022年，8月28日 等参数单元 而以上各式中的 和 ，可由式（8-9）分别对偏微分而求得。这样就把 和 转化成了局部坐标的函数，从而求的应变矩阵[B]和单元应变[ε]。 将单元应变代入平面问题的物理方程式，就得到平面八结点等参元的应力列阵 （i=1, 2, …, 8） 式中，[S]为应力矩阵 返回 第三十八页，共七十四页，2022年，8月28日 等参数单元 利用虚功原理可以得到其单刚矩阵 式中，t为单元厚度。 把（8-43）式写成分块矩阵，可分成8×8个子矩阵，每个子矩阵都是2×2阶矩阵，即 返回 第三十九页，共七十四页，2022年，8月28日 应该指出，上式是对ξ和η的重积分，尽管其积分区域十分简单，但其被积函数却比较复杂，需要采用数值积分法求解（通常是采用高斯积分法）。 等参数单元 其中子矩阵 返回 第四十页，共七十四页，2022年，8月28日 等参数单元 四、等效结点载荷 　　整体结构结点载荷列阵是通过将作用在单元上的集中力，表面力和体积力分别等效移置到结点后，经过组集得到 1. 集中力的等效结点载荷 　　设单元任意点c作用有集中载荷 ，则移置到单元各有关结点上的等效结点载荷为 式中 (Ni)c是形函数Ni在集中力作用点c处的取值，可通过以下 返回 第四十一页，共七十四页，2022年，8月28日 等参数单元 步骤计算： 1）根据作用点c的整体坐标 ，得到其局部坐标 。 　　式中： 均为已知数。解此联立方程式就得到c点的局部坐标