一种提高平板型壳单元计算精度的改进算法研究

摘要

壳体结构是工程应用中常见的结构形式之一，广泛应用于土木、水利、汽车、航空航天等各个工程领域中。由于壳体是具有薄膜和弯曲两重受力特性的三维特殊结构，其控制方程是一组形式复杂的偏微分方程，直接解析求解是比较困难的。在电子计算机出现以前只能计算一些结构简单规则的问题，随着以计算机为工具的有限元方法出现以后，经过众多学者的努力，不同类型的壳单元也相继出现。其中，平板型壳单元是最早发展起来的壳单元，它是将壳体离散为一系列折板组成的体系，通过平面膜单元和平板弯曲单元叠加组合来模拟壳体的薄膜和弯曲受力状态。平板型壳单元以其表达格式简单、应用广泛、计算效率高、计算结果可靠等优点在不同类型的壳单元中占有重要的地位，在工程实际中得到了广泛的应用。但是，平板型壳单元是由平面应力单元和平板弯曲单元叠加组合而成的壳单元，采用几何近似性，用折板面代替壳体曲面。因此，在计算不规则壳体结构特别是曲壳结构时不能很好地模拟曲壳形状，从而引起较大的误差，需要将网格合理加密才能达到一定的精度，这会引起较大的计算量。
　　为了进一步提高平板型壳单元的计算性能，本文提出了一种新的单元局部坐标系建立方法。该方法在单元的高斯积分点建立多个局部坐标系，并保证每个局部坐标系都位于单元在高斯点处的切平面上，对每个单元刚度转换矩阵进行及时修正，从而可以有效适应曲面壳体形状。为了能够在新的局部坐标系下求得和单元刚度有关的参数，同时又避免频繁的坐标转换造成的计算误差，本文还引入了计算形函数对局部坐标的导数和积分转换的雅可比行列式的高精度算法。这种高精度算法直接基于三维曲面整体坐标系推导得出，具有较高的计算精度。利用本文提出的改进算法计算平板型壳单元的刚度矩阵以及等效节点载荷列阵，都取得了良好的效果。此外，本文还将一种新的线性方程组求解器应用到平板型壳单元当中，进一步提高其计算性能。本文的主要研究内容可以归纳如下:
　　(1)提出了一种计算平板型壳单元在整体坐标系下的单元刚度矩阵的改进算法，并将其应用到几种不同类型的平板型壳单元当中。通过将几种不同类型的平面膜单元和平板弯曲单元叠加组合形成平板型壳单元，对其计算性能进行讨论，并详细给出了利用本文的改进算法计算单元刚度矩阵的理论公式，从原理上对本文方法在提高平板型壳单元计算性能方面进行分析讨论。
　　(2)对本文所提出的改进算法在平板型壳单元中的应用和实现进行程序开发，给出了详细的编程流程，并通过算例对其计算精度进行比较分析。为了对本文提出的改进算法的计算性能进行验证，将改进算法应用到几种典型的平板型壳单元，并对其进行程序实现，通过具体的算例对其正确性和精确性进行验证分析。
　　(3)将本文提出的改进算法应用扩展到平板型壳单元等效节点载荷的计算上。对平板型壳单元在不同载荷下等效节点载荷的计算进行归纳总结，给出了利用本文方法计算其等效节点载荷的详细的理论公式，并利用程序编制对其应用进行数值实现，通过具体的算例对本文方法在计算等效节点载荷方面的性能进行分析验证。
　　(4)在有限元分析中，引入了一种新的线性方程组直接法求解器——同时消元回代法，提出了新的刚度矩阵组集方法，并将其应用到平板型壳单元问题的求解中。给出了新的刚度矩阵的组集方法和详细计算步骤，并结合同时消元回代法求解器对其求解性能进行分析讨论，给出了程序的编制过程，通过具体算例对其性能进行验证。
　　综上所述，本文提出了一种能够用于提高平板型壳单元计算精度的改进算法，为平板型壳单元计算性能的提高提供了一种新的手段。本文所述方法也可以应用到其他类型的四边形平板型壳单元中，具有广阔的发展前景。

目录

声明

摘要

主要符号表

1 绪论

1．1 研究背景与意义

1．2 板壳有限元发展现状

1．3 平板型壳单元发展现状

1．4 本文的主要研究内容

2 四边形平板型壳单元刚度矩阵的改进算法

2．1 引言

2．2 建立在高斯积分点处的单元局部坐标系

2．2．1 单元局部坐标系建立的传统方法

2．2．2 在单元切平面上建立局部坐标系

2．2．3 在高斯积分点处的切平面上建立局部坐标系

2．3 单元刚度积分相关参数的高精度算法

2．3．1 单元刚度积分相关参数的传统算法

2．3．2 单元刚度积分相关参数的高精度算法

2．4 本章小结

3 改进算法在一般四边形平板型壳单元中的实现

3．1 引言

3．2 四节点平板型壳单元

3．2．1 四节点平板型壳单元的构成

3．2．2 四节点平板型壳单元刚度矩阵的改进算法

3．2．3 程序实现

3．2．4 数值方法验证

3．3 八节点平板型壳单元

3．3．1 八节点平板型壳单元的构成

3．3．2 八节点平板型壳单元刚度矩阵的改进算法

3．3．3 程序实现

3．3．4 数值方法验证

3．4 本章小结

4 改进算法在带旋转自由度的平板型壳单元中的实现

4．1 引言

4．2 带旋转自由度的四节点平面膜单元

4．2．1 四节点任意四边形膜单元GQ12

4．2．2 四节点任意四边形膜单元GQ12M

4．3 带旋转自由度的四节点平板型壳单元刚度矩阵的改进算法

4．4 程序实现

4．5 数值方法验证

4．5．1 四节点膜单元GQ12和GQ12M采用改进算法的算例验证

4．5．2 带旋转自由度的平板型壳单元采用改进算法的算例验证

4．6 本章小结

5 基于改进算法的平板型壳单元等效节点载荷计算

5．1 引言

5．2 平板型壳单元等效节点载荷

5．2．1 平面应力单元等效节点载荷

5．2．2 平板弯曲单元等效节点载荷

5．2．3 平板型壳单元等效节点载荷的改进算法

5．3 程序实现

5．4 数值方法验证

5．5 本章小结

6 可用于平板型壳单元问题求解的有限元新型求解器

6．1 引言

6．2 线性方程组的同时消元回代求解法(SEBSM)

6．3 同时消元回代法在有限元中的应用

6．3．1 有限元整体刚度矩阵组集的新方法

6．3．2 性能讨论

6．3．3 数值方法验证

6．4 同时消元回代法在平板型壳单元中的应用

6．5 本章小结

7 结论与展望

7．1 结论

7．2 创新点

7．3 展望

参考文献

附录

攻读博士学位期间科研项目及科研成果

致谢

作者简介