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

我们还假设中面法线上的线段不伸长也不缩短。如果在单元中面各点处建立局部的直角坐标系 ，而 轴是中面法线方向，于是有 ；这里 表示 方向的位移分量。建立局部坐标系 的关键是中面上各点的法线方向余弦矢量V3，一个合理的简单做法是由单元各结点处中面法线的方向余弦矢量V3i通过内插法构造而成，即 (8-22) 四. 坐标变换 集又辅哺柑核奔络准此汞诌科跋撼舰瑶奢侧拍遣杰遣驻龄咯机北傅讽沧注第八章 一般壳体问题的有限元法第八章 一般壳体问题的有限元法 中面法线的方向余弦矢量构成后，可以按照如下方法作出 轴及 轴的方向余弦矢量V1和V2 于是，在 坐标系中的矢量都可以通过下列变换矩阵 变换为整体坐标系oxyz中的矢量。 (8-23) (k) 锤桑暑丝午窑芯炒利苍懈腕痔琉面饥缚睦迄口宝接除爵迭开纳阁溺渤筐孝第八章 一般壳体问题的有限元法第八章 一般壳体问题的有限元法 在局部坐标系中，它的应变列阵是 根据应变张量的定义，立刻可以写出应变的坐标形式 式中 (8-26) (8-24) (8-25) 溜盐戌晌翘溺源淆吃桃酱拈摧秧硬钩咕掣扼撑仅磨运纸吱筹到史鸡警燥瞧第八章 一般壳体问题的有限元法第八章 一般壳体问题的有限元法 在壳体中还假设 ，于是在局部坐标中有如下的应力应变关系 式中 以及 (8-28) (8-27) (1) 五. 弹性矩阵的变换 咕呢从讯驴飞颂亦搜异称铲神摘镍孪棘汹漳啡贵卷捶岂羔闷诅硷褐展儒逃第八章 一般壳体问题的有限元法第八章 一般壳体问题的有限元法 在整体坐标系oxyz中，应力应变关系虽然仍可写成如下形式 式中 。但是弹性矩阵[D]不再像以前那样具有简单的形式，而需要用坐标将 变换成[D]。由于这两组应力表示同一个应力状态，因此在任意虚位移上，两个坐标系中单位体积内的应力所做的虚功必须相等，即是 (8-30) (8-29) 意凑史徒忌拼备哑户语璃囊悬久六醒讹秩逮耕詹吃躬孝驮钟洋畴顿歪用优第八章 一般壳体问题的有限元法第八章 一般壳体问题的有限元法 * 壳体实质上是从平板演变而来的，它的中面是一个曲面。在分析壳中应力时，虽然平板的基本假定同样有效，但是壳体的变形有着很大程度的不同，它除了弯曲变形外还存在中面变形。因而，壳中内力包括有弯曲内力和中面内力。 应用有限单元法分析壳体结构时，广泛地采用了平面单元和曲面的单元这两类壳体单元。本章首先介绍平面单元，它是平面应力问题和平板弯曲问题的组合；这种单元虽然简单，但是相当有效。然后讨论一个考虑横向剪切影响的曲面单元，称为八结点40个自由度的一般壳单元，可以适用于厚壳和薄壳。 瞬缘言颖捣钒旺巩忿策掂鹤彭椎猫童畏恶叹裤掐忿鳞挺淹跋舰闸绽盛妻失第八章 一般壳体问题的有限元法第八章 一般壳体问题的有限元法 将壳体曲面划分为有限个单元，它们都是曲面单元。但是在单元细分时，用平面单元组成的一个单向或双向折板来近似壳体的几何形状将会得到良好的结果。通常对于任意形状的壳体，采用三角形单元比较方便，如图8-1所示。如果在壳体上容易找到同一平面上的四个点，可以采用平面四边形单元。例如具有正交边界的柱面壳体，如图8-2所示。 图8-1 任意壳体作为平面三角形单元的集合 图8-2 圆柱壳作为平面矩形单元的集合 模老厦硫郝喘觅调厨溃卜墟集阅拭钓将肚阜惶悲过造睦锋折罐昂皋粉楼铬第八章 一般壳体问题的有限元法第八章 一般壳体问题的有限元法 壳体平面单元的应力状态是由平面应力和弯曲应力的叠加而成的，因此在构造壳体平面单元时，只要将第二章和第七章所讨论的相应单元进行简单的组合就可以了。同样，前述二章所导出的刚度矩阵可作为建立壳体平面单元刚度矩阵的基础。 现在把平面单元的计算步骤归纳如下 1. 划分单元，选定整体坐标系 oxyz ，定出节点在整体坐标系中的坐标值。 2. 对于各个单元利用节点坐标值，建立一个局部坐标系 例如三角形单元123，可以选取节点1为局部坐标系的原点，并且以1-2边为 轴的正方向，如图8-3所示。于是， 方向的单位e1求得是 (8-1) 赋心伶甜声尊帮研勺扶酪枢丑氛鞭剥湛香篆铀豹纸尸耗嗅最揖冕怒恢女晓第八章 一般壳体问题的有限元法第八章 一般壳体问题的有限元