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

我们还假设中面法线上的线段不伸长也不缩短。如果在单元中面各点处建立局部的直角坐标系 ，而 轴是中面法线方向，于是有 ；这里 表示 方向的位移分量。建立局部坐标系 的关键是中面上各点的法线方向余弦矢量V3，一个合理的简单做法是由单元各结点处中面法线的方向余弦矢量V3i通过内插法构造而成，即 (8-22) 四. 坐标变换 中面法线的方向余弦矢量构成后，可以按照如下方法作出 轴及 轴的方向余弦矢量V1和V2 于是，在 坐标系中的矢量都可以通过下列变换矩阵 变换为整体坐标系oxyz中的矢量。 (8-23) (k) 在局部坐标系中，它的应变列阵是 根据应变张量的定义，立刻可以写出应变的坐标形式 式中 (8-26) (8-24) (8-25) * 壳体实质上是从平板演变而来的，它的中面是一个曲面。在分析壳中应力时，虽然平板的基本假定同样有效，但是壳体的变形有着很大程度的不同，它除了弯曲变形外还存在中面变形。因而，壳中内力包括有弯曲内力和中面内力。 应用有限单元法分析壳体结构时，广泛地采用了平面单元和曲面的单元这两类壳体单元。本章首先介绍平面单元，它是平面应力问题和平板弯曲问题的组合；这种单元虽然简单，但是相当有效。然后讨论一个考虑横向剪切影响的曲面单元，称为八结点40个自由度的一般壳单元，可以适用于厚壳和薄壳。 将壳体曲面划分为有限个单元，它们都是曲面单元。但是在单元细分时，用平面单元组成的一个单向或双向折板来近似壳体的几何形状将会得到良好的结果。通常对于任意形状的壳体，采用三角形单元比较方便，如图8-1所示。如果在壳体上容易找到同一平面上的四个点，可以采用平面四边形单元。例如具有正交边界的柱面壳体，如图8-2所示。 图8-1 任意壳体作为平面三角形单元的集合 图8-2 圆柱壳作为平面矩形单元的集合 壳体平面单元的应力状态是由平面应力和弯曲应力的叠加而成的，因此在构造壳体平面单元时，只要将第二章和第七章所讨论的相应单元进行简单的组合就可以了。同样，前述二章所导出的刚度矩阵可作为建立壳体平面单元刚度矩阵的基础。 现在把平面单元的计算步骤归纳如下 1. 划分单元，选定整体坐标系 oxyz ，定出节点在整体坐标系中的坐标值。 2. 对于各个单元利用节点坐标值，建立一个局部坐标系 例如三角形单元123，可以选取节点1为局部坐标系的原点，并且以1-2边为 轴的正方向，如图8-3所示。于是， 方向的单位e1求得是 (8-1) 式中 是矢量12的长度。取单元的外法线方向作为 轴的正方向，于是它的单位矢量 图8-3 三角形单元局部坐标系 容易看出，矢量12和13的矢性积的模等于三角形面积Δ的一倍，即|12×13|=2Δ。最后，按右手定则可以决定y轴的正方向,它的单位矢量e2是 e2 = e3 ×e1 (8-3) 利用上述方法确定的局部坐标系，三角形单元123是在 平面内，它的三个角点的局部坐标值是很容易确定的。 对于柱面上的矩形单元，局部坐标的原点 选在矩形的形心，通常选 轴和x轴均沿柱面母线方向。如图8-4中所示，由矢量12确定单位矢量e1,再由矢量14确定单位矢量e2，于是e3 = e1 ×e2。 3.对于各个单元，确立在局部坐标系 中的结点载荷列阵 。壳体载荷可以分解成二组：一组作用在平面内，另一组垂直于平面。为此，在计算各个单元的结点载荷列阵 （包括等效结点力）可以直接引用第二章和第七章中所叙述的载荷计算的相应公式。 各个单元的结点载荷列阵 图8-4 矩形单元局部坐标 求得后，建立变换矩阵公式，从而把 转换到整体坐标系中去求出在整体坐标下的单元节点载荷列阵，然后经各单元的简单叠加可以求出结构在整体坐标下的节点载荷列阵。 显然，平面单元在局部坐标系中，结点i有五个广义位移：即 ，其中前两个对应于平面应力问题，后三个对应于平板弯