第4章+壳体问题的有限单元法.pdfVIP

第三章 壳体问题的有限单元法 用有限单元法分析弹性薄壳，有两种不同的途径。一个途径是直接采用曲面单元，根据 壳体理论推导单元刚度矩阵。另一途径是，用薄板单元组成的折板系统去代替原来的薄壳， 由平面应力状态和板弯曲应力状态加以组合而得到薄壳的应力状态。下面介绍曲面薄壳单 元，然后介绍用折板代替薄壳的计算方法，如图 1 所示，通常用三角形或矩形薄板单元的组 合去代替壳体，其中以三角形单元应用较广、实用价值较大，因为它可以适应壳体的复杂外 形。 图 1 第一节 基于薄壳理论的轴对称壳体单元 轴对称壳体在工程实际中得到广泛的应用。基于薄壳理论的轴对称壳体单元，在厚度方 向引入了壳体理论中的 Kirchhoff 假设。这样一来，轴对称壳体单元本质上是一维单元，从 而使整个分析大为简化。 轴对称壳体单元，最早提出的是在子午线方向为直线的截锥单元。此种单元表达式简单， 用于实际分析一般也可达到合理的精度。但是此种单元在模拟曲率较大的壳体时，不仅需要 较多的单元，而且还可能产生附加弯矩。因此，以后很多研究工作者又提出一系列在子午线 方向为曲线的单元。 为了进一步将轴对称壳元能用于中等厚度壳体的分析，Zienkiewicz等人又先后提出了 从轴对称实体单元蜕化得到的超参数壳体单元，以及挠度和截面转动作为各自独立函数的壳 体单元，此类单元的特点是单元之间只要保持C 连续性。同时由于考虑横向剪切变形的影响， 0 所以还能用于中厚壳的分析。 一、轴对称薄壳理论的基本公式 图1.1 所示轴对称壳体中面上任一点位置由 φ 和θ确定，位移可由其经（子午）向分量 u 、 周向分量 v 和法向分量 w 确定。 在薄壳理论中，壳体内任一点的应变，根据 Kirchhoff 直法线假设，可通过中面的六个广义 图 1.1 应变分量来描述，它们和中面位移的关系是 177 ⎧∂u w ⎫ ⎪∂s +R ⎪ ⎪ s ⎪ ⎪1 ⎛∂υ ⎞ ⎪ +u +w ⎪ ⎜ sinφ cosφ⎟ ⎪ ⎧ ⎫ε ⎪r ⎝∂θ ⎠ ⎪ s ⎪ ⎪ε ⎪∂υ 1 ⎛∂u ⎞ ⎪ θ + −υ φ ⎪ ⎪ ⎪ ⎜ sin ⎟ ⎪ ⎪ ⎪γ ⎪∂s r ⎝∂θ ⎠ ⎪ ⎪ ⎪ sθ {}ε ⎨ ⎬ ⎨ ∂ ∂⎛ w u ⎞