第八章一般壳体问题的有限元法.pptVIP

式中 而 (j) (8-19) 第29页，共48页，编辑于2022年，星期一 因此，为了计算应变需要算出六个偏导数 、 、 和 、 前三个偏导数可以按照公式（5-54）进行计算，后三个则按定义分别是[J]-1的第三列的三个元素，即 由（5-58）式得到 (8-20) 第30页，共48页，编辑于2022年，星期一 出现在所提到的各式中的S、T和V分别按公式（5-55）所定义。将（8-11）式代入（5-55）式得 (8-21) 第31页，共48页，编辑于2022年，星期一 我们还假设中面法线上的线段不伸长也不缩短。如果在单元中面各点处建立局部的直角坐标系 ，而 轴是中面法线方向，于是有 ；这里 表示 方向的位移分量。建立局部坐标系 的关键是中面上各点的法线方向余弦矢量V3，一个合理的简单做法是由单元各结点处中面法线的方向余弦矢量V3i通过内插法构造而成，即 (8-22) 四. 坐标变换 第32页，共48页，编辑于2022年，星期一 中面法线的方向余弦矢量构成后，可以按照如下方法作出 轴及 轴的方向余弦矢量V1和V2 于是，在 坐标系中的矢量都可以通过下列变换矩阵 变换为整体坐标系oxyz中的矢量。 (8-23) (k) 第33页，共48页，编辑于2022年，星期一 在局部坐标系中，它的应变列阵是 根据应变张量的定义，立刻可以写出应变的坐标形式 式中 (8-26) (8-24) (8-25) 第34页，共48页，编辑于2022年，星期一 第1页，共48页，编辑于2022年，星期一 将壳体曲面划分为有限个单元，它们都是曲面单元。但是在单元细分时，用平面单元组成的一个单向或双向折板来近似壳体的几何形状将会得到良好的结果。通常对于任意形状的壳体，采用三角形单元比较方便，如图8-1所示。如果在壳体上容易找到同一平面上的四个点，可以采用平面四边形单元。例如具有正交边界的柱面壳体，如图8-2所示。 图8-1 任意壳体作为平面三角形单元的集合 图8-2 圆柱壳作为平面矩形单元的集合 §8-1 平面壳体单元 第2页，共48页，编辑于2022年，星期一 壳体平面单元的应力状态是由平面应力和弯曲应力的叠加而成的，因此在构造壳体平面单元时，只要将第二章和第七章所讨论的相应单元进行简单的组合就可以了。同样，前述二章所导出的刚度矩阵可作为建立壳体平面单元刚度矩阵的基础。 现在把平面单元的计算步骤归纳如下 1. 划分单元，选定整体坐标系 oxyz ，定出节点在整体坐标系中的坐标值。 2. 对于各个单元利用节点坐标值，建立一个局部坐标系 例如三角形单元123，可以选取节点1为局部坐标系的原点，并且以1-2边为 轴的正方向，如图8-3所示。于是， 方向的单位e1求得是 (8-1) 第3页，共48页，编辑于2022年，星期一 式中 是矢量12的长度。取单元的外法线方向作为 轴的正方向，于是它的单位矢量 图8-3 三角形单元局部坐标系 第4页，共48页，编辑于2022年，星期一 容易看出，矢量12和13的矢性积的模等于三角形面积Δ的一倍，即|12×13|=2Δ。最后，按右手定则可以决定y轴的正方向,它的单位矢量e2是 e2 = e3 ×e1 (8-3) 利用上述方法确定的局部坐标系，三角形单元123是在 平面内，它的三个角点的局部坐标值是很容易确定的。 对于柱面上的矩形单元，局部坐标的原点 选在矩形的形心，通常选 轴和x轴均沿柱面母线方向。如图8-4中所示，由矢量12确定单位矢量e1,再由矢量14确定单位矢量e2，于是e3 = e1 ×e2。 第5页，共48页，编辑于2022年，星期一 3.对于各个单元，确立在局部坐标系 中的结点载荷列阵 。壳体载荷可以分解成二组：一组作用在平面内，另一组垂直于平面。为此，在计算各个单元的结点载荷列阵 （包括等效结点力）可以直接引用第二章和第七章中所叙述的载荷计算的相应公式。 各个单元的结点载荷列阵 图8-4 矩形单元局部坐标