三维薄结构热传导问题分层二次元方程的多水平方法

在利用有限元法对三维薄结构进行分析时,为了减少单元数目,常采用六面体薄单元,相应的高阶单元在计算精度、抗畸变程度等方面具有明显优势.但与低阶元相比,高阶单元需要更多的计算机存储空间,离散化线性系统具有更高的计算复杂性,并且系数矩阵是严重病态的,采用通常的求解方法其效率将大大降低.该文针对三维薄结构稳态热传导问题,利用局部块Gauss-Seidel光滑子和基于"距离矩阵"的DAMG法,为其分层二次元离散系统设计了一种具有更好计算效率和鲁棒性(robustness)的多水平方法.由于采用了分层基,程序实现中不再需要建立判定未知数变量指标与所属几何节点类型对应关系的代数判据,网格转换算子的构造也变得非常简单,从而大大提高了运算效率.数值实验结果验证了该方法的有效性和鲁棒性.