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

　 ， ３９ ６ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１８年６月１ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３９，Ｎｏ．６，Ｊｕｎ．１，２０１８ 文章编号：１０００⁃０８８７（２０１８）０６⁃０７００⁃１４ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 三维薄结构热传导问题分层 二次元方程的多水平方法∗ 张　 申，　 肖映雄，　 郭瑞奇 （湘潭大学 土木工程与力学学院，湖南 湘潭４１１１０５） 摘要：　 在利用有限元法对三维薄结构进行分析时，为了减少单元数目，常采用六面体薄单元，相 应的高阶单元在计算精度、抗畸变程度等方面具有明显优势．但与低阶元相比，高阶单元需要更多 的计算机存储空间，离散化线性系统具有更高的计算复杂性，并且系数矩阵是严重病态的，采用通 常的求解方法其效率将大大降低．该文针对三维薄结构稳态热传导问题，利用局部块 Ｇａｕｓｓ⁃Ｓｅｉｄｅｌ 光滑子和基于“距离矩阵”的ＤＡＭＧ法，为其分层二次元离散系统设计了一种具有更好计算效率 和鲁棒性（ｒｏｂｕｓｔｎｅｓｓ）的多水平方法．由于采用了分层基，程序实现中不再需要建立判定未知数变 量指标与所属几何节点类型对应关系的代数判据，网格转换算子的构造也变得非常简单，从而大 大提高了运算效率．数值实验结果验证了该方法的有效性和鲁棒性． 关　 键　 词：　 热传导问题；　 薄体结构；　 分层二次元；　 多水平方法；　 代数多重网格法 中图分类号：　 Ｏ３４３．３；ＴＢ１１５　 　 　 文献标志码：　 Ａ ＤＯＩ：１０．２１６５６／ １０００⁃０８８７．３８００３５ 引　 　 言 薄结构由于具有某些方面的特点，已被广泛应用于各类实际问题中，如涂层问题、各种层 合结构问题以及电子器件的散热问题等．利用通常的体单元（如８节点线性元）分析薄结构问 题时，由于单元最大几何尺度受薄层厚度控制，需要使用更多的单元数才能确保数值解的计算 ［１⁃２］ 精度，但由此也将导致计算工作量剧增 ．实际计算时，需要采用薄单元（其一个方向的尺寸 远小于其他两个方向的尺寸）进行网格剖分，这种单元属于典型的各向异性单元．相应的高阶 单元在计算精度、划分数量、抗畸变程度等方面具有非常明显的优势，使得它们在实际计算中 被广泛使用．但与低阶元相比，高阶单元需要更多的计算机存储空间，离散化线性系统具有更 高的计算复杂性．另一方面，这类薄单元对应的线性代数系统的系数矩阵是严重病态的，采用 通常的求解方法其效率将大大降低． 代数多重网格（ａｌｇｅｂｒａｉｃ ｍｕｌｔｉｇｒｉｄ，ＡＭＧ）法是求解有限元离散化线性代数系统最为有效 的计算方法之一，目前已有非常多的研究工作，如文献［３⁃１０］．对高阶元离散系统ＡＭＧ法（包 括多水平方法及其预条件技术）的研究，也取得了一些很好的成果，如文献［１１⁃１５］．但这些方 ∗ 收稿日期：　 ２０１７⁃０２⁃１６；修订日期：　 ２０１７⁃１１⁃３０ 基金项目：　 国家自然科学基金；湖南省教育厅资助科研项目（１５Ａ１８３） 作者简介：　 张申（１９９０—），男，硕士生（Ｅ⁃ｍａｉｌ：ｚｈａｎｇｓｈｅｎｆ＠１２６．ｃｏｍ）； 肖映雄（１９７０—），男，教授，博士（通讯作者． Ｅ⁃ｍａｉｌ：ｘｙｘ６１０ｘｙｘ＠ｘｔｕ．ｅｄｕ．ｃｎ）； 郭瑞奇（１９９３—），男，硕士生（Ｅ⁃ｍａｉｌ：１９１５２２１７７＠ｑｑ．ｃｏｍ）． ７００ 张　 　 申　 　 　 肖　 映　 雄　 　 　 郭　 瑞　 奇 ７０１ 法主要是针对各向同性网格下的高阶有限元离散代数系统而设计的．对各向异性网格情形，文 献［１６⁃１７］研究了几种有代表性的复杂各向异性网格对 ＡＭＧ法效率的影响，并设计了基于 “距