[工程科技]板壳问题的有限元法.pdfVIP

11 板壳问题的有限元法 本章提要 本章叙述了薄板弯曲问题有限元法的矩形单元、三角形单元的计算原理与方法，给出了 矩形薄板单元、三角形薄板单元结合相应平面应力单元计算板壳问题的基本过程。介绍了适 用于薄板弯曲问题、平面应力问题及板壳结构等分析计算的程序。 11.1 引言 从第10 章我们已经知道，有限单元法用于三维连续体，特别是它们的线性分析，已经 发展得相当成熟。从原则上说，我们也可以利用三维实体单元分析板壳结构问题，并可以避 免弹性力学引入的简化，但是这样做在实际分析中遇到了困难。这是由于在用实体单元对结 构进行离散时，如果网格适应结构的几何特点，这时单元的厚度较其它方向的尺寸要小得多， 于是单元在不同方向的刚度系数相差过大，从而导致结构整体方程的病态或奇异，最后将使 解丧失精度或根本失效。反之，为了避免上述问题，保持单元在各个方向尺度相近，将导致 单元总数过分庞大，而使整个分析过程难以进行。 事实上，目前弹性力学板壳理论已较为成熟，且在工程应用中被证明行之有效。但由 于数学求解上的限制，只有少量几何形状、荷载及边界条件较简单的情况才能获得解析的结 果。因此，自有限元产生开始，以弹性力学板壳理论为基础的板壳问题有限元法得到了广泛 的研究与应用，积累了丰富的文献资料和大量实例分析。本章介绍的是以经典板壳理论为基 础的有限元板壳分析，为进一步深入学习打下基础。 11.2 矩形单元薄板分析 在本章板壳问题的有限元法中，板涉及的是薄板，壳则是用薄板弯曲及平面应力问题 的组合分析来完成。因此本章的板壳问题有限元分析，是以经典薄板理论为基础，经典薄板 问题的基本方程需要引用，但均改用矩阵表示，如广义应变 T ⎡ ∂2 w ∂2 w ∂2 w ⎤ {κ} ⎢− 2 − 2 −2 ⎥ （11.1） x y ⎣ ∂x ∂y ∂ ∂ ⎦ {κ} 中各个分量分别代表薄板弯曲后中面在x 方向的曲率、y 方向的曲率以及在x 和y 方向 的扭率。与之相应的薄板的广义内力是 1 {M } [M x M y M xy ]T （11.2 ） 其中，M ，M 分别是垂直于x轴和垂直于y 轴的截面上单位长度的弯矩，M ( M ) 是垂 x y xy yx 直于x （y ）轴截面上单位长度的扭矩。 广义应力应变关系是 {M } [D ]{κ} （11.3 ） 其中弹性矩阵[D]，对于各向同性材料是 ⎡ ⎤ 1 v 0 3 ⎢ ⎥ [D ] Et 2 ⎢v 1 0 ⎥