一種求解三維Neumann邊界條件線彈性問題二次有限元方程的并行高效預條件子
發(fā)布時間:2019-04-21 07:43
【摘要】:帶純Neumann邊界條件三維線彈性模型是描述固體力學問題和計算材料力學問題的重要模型,有限元法是數(shù)值求解三維線彈性問題最常用的離散方法,并行代數(shù)多層網(wǎng)格(AMG)法是快速求解三維線彈性問題離散系統(tǒng)的最為有效的方法之一。本文針對一種滿足適定性條件的純Neumann邊界的三維線彈性問題的二次有限元離散代數(shù)系統(tǒng),研究其并行快速求解算法及解法器。首先,在文[33]的基礎上,討論了二次有限元離散化代數(shù)系統(tǒng)的適定性,數(shù)值實驗表明二次有限元誤差函數(shù)在L2(Ω)和H1(Ω)范數(shù)下均具有飽和誤差階。接著,針對二次有限元代數(shù)系統(tǒng)給出了兩類求解算法,其中重點研究了兩種基于AMG法和高斯賽德爾迭代法(GS)的組合型預條件子,并研制了相應的預條件GMRES(k)解法器(BAG-GMRES(k)和BvAG-GMRES(k))。數(shù)值實驗驗證了基于AMG的組合型預條件GMRES(k)法的迭代次數(shù)和求解時間均優(yōu)于常用求解方法ILU(0)-GMRES(k)。相比較而言,BvAG-GMRES(k)解法器更穩(wěn)定高效。最后,在上述串行預條件子算法解法器的基礎上,設計了相應的具有極小化數(shù)據(jù)通信的并行AMG法和GS的組合型預條件算法,數(shù)值試驗驗證了該算法具有良好的擴展性。
[Abstract]:The three-dimensional linear elastic model with pure Neumann boundary conditions is an important model for describing the problems of solid mechanics and computational material mechanics. The finite element method is the most commonly used discrete method for numerical solution of three-dimensional linear elastic problems. The parallel algebraic multi-layer grid (AMG) method is one of the most effective methods for fast solving three-dimensional linear elastic problem discrete systems. In this paper, a quadratic finite element discrete algebraic system for three-dimensional linear elastic problem with pure Neumann boundary satisfying the well-posed condition is studied, and its parallel fast algorithm and solver are studied. Firstly, on the basis of [33], the well-posed quadratic finite element discrete algebraic system is discussed. Numerical experiments show that the quadratic finite element error function has saturation error order in L2 (惟) and H1 (惟) norm. Then, two kinds of solving algorithms are given for quadratic finite element algebraic system, in which two kinds of combinatorial preforms based on AMG method and Gao Si Seidel iterative method (GS) are studied emphatically. The corresponding preconditioned GMRES (k) solver (BAG-GMRES (k) and BvAG-GMRES (k).) is developed. Numerical experiments show that the iteration times and solution time of the combined preconditioned GMRES (k) method based on AMG are better than that of ILU (0)-GMRES (k). Method. In contrast, the BvAG-GMRES (k) solver is more stable and efficient. Finally, on the basis of the subalgorithm solver of serial prebar, the parallel AMG method with minimized data communication and the combined preconditioned algorithm of GS are designed. Numerical experiments show that the algorithm has good expansibility.
【學位授予單位】:湘潭大學
【學位級別】:碩士
【學位授予年份】:2017
【分類號】:O343.2;O241.82
本文編號:2461984
[Abstract]:The three-dimensional linear elastic model with pure Neumann boundary conditions is an important model for describing the problems of solid mechanics and computational material mechanics. The finite element method is the most commonly used discrete method for numerical solution of three-dimensional linear elastic problems. The parallel algebraic multi-layer grid (AMG) method is one of the most effective methods for fast solving three-dimensional linear elastic problem discrete systems. In this paper, a quadratic finite element discrete algebraic system for three-dimensional linear elastic problem with pure Neumann boundary satisfying the well-posed condition is studied, and its parallel fast algorithm and solver are studied. Firstly, on the basis of [33], the well-posed quadratic finite element discrete algebraic system is discussed. Numerical experiments show that the quadratic finite element error function has saturation error order in L2 (惟) and H1 (惟) norm. Then, two kinds of solving algorithms are given for quadratic finite element algebraic system, in which two kinds of combinatorial preforms based on AMG method and Gao Si Seidel iterative method (GS) are studied emphatically. The corresponding preconditioned GMRES (k) solver (BAG-GMRES (k) and BvAG-GMRES (k).) is developed. Numerical experiments show that the iteration times and solution time of the combined preconditioned GMRES (k) method based on AMG are better than that of ILU (0)-GMRES (k). Method. In contrast, the BvAG-GMRES (k) solver is more stable and efficient. Finally, on the basis of the subalgorithm solver of serial prebar, the parallel AMG method with minimized data communication and the combined preconditioned algorithm of GS are designed. Numerical experiments show that the algorithm has good expansibility.
【學位授予單位】:湘潭大學
【學位級別】:碩士
【學位授予年份】:2017
【分類號】:O343.2;O241.82
【參考文獻】
相關期刊論文 前3條
1 肖映雄;周志陽;舒適;;幾類典型網(wǎng)格下三維彈性問題的代數(shù)多層網(wǎng)格法[J];工程力學;2011年06期
2 肖映雄;張紅梅;舒適;;三維彈性問題高次有限元方程的代數(shù)多層網(wǎng)格法[J];計算力學學報;2010年06期
3 張平;肖映雄;舒適;;代數(shù)多層網(wǎng)格法及其在固體力學計算中的應用研究[J];湘潭大學自然科學學報;2008年03期
相關博士學位論文 前1條
1 李政;精細油藏數(shù)值模擬中的高效求解器研究[D];昆明理工大學;2017年
相關碩士學位論文 前4條
1 王錦;一種求解三維Neumann邊界條件線彈性問題線性有限元方程的高效預條件子[D];湘潭大學;2016年
2 張楚琴;一種求解線彈性問題的基于Laplace算子的并行DDM預條件子[D];湘潭大學;2015年
3 岳孝強;幾種并行AMG法及其在輻射擴散問題中的應用[D];湘潭大學;2012年
4 梁文濤;一種求解三維彈性問題有限元方程的并行DDM預條件子[D];湘潭大學;2010年
,本文編號:2461984
本文鏈接:http://sikaile.net/kejilunwen/yysx/2461984.html
最近更新
教材專著