第五章等参数单元.pptVIP

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