三維六面體網格自適應生成算法研究及其應用.pdf

1、隨著計算機技術和數(shù)值方法的發(fā)展，有限元法、有限體積法和有限差分法等數(shù)值方法已成為計算科學與工程問題的重要方法，其數(shù)學基礎是變分原理，基本思想是離散化和分片插值，即用網格去逼近所分析物體的空間。數(shù)值分析或模擬的主要步驟包括幾何建模、模型離散、求解、后處理等過程，其中首要任務是模型離散，其依賴于有限元網格生成技術的發(fā)展。據(jù)調查統(tǒng)計，模型離散(即網格劃分)通常占整個項目完成時間和費用的80％左右，而且網格的質量直接影響數(shù)值分析結果的精度和效率

2、。在三維分析問題中，有四面體和六面體及二者混合單元，目前四面體網格的劃分方法和有關商業(yè)化軟件已漸成熟，而六面體網格由于其復雜的拓撲關系和與物體表面邊界的一致性等方面的限制，現(xiàn)有六面體網格劃分方法中都存在各自的優(yōu)缺點。但六面體網格在計算精度、劃分數(shù)量、抗畸變程度及重劃分次數(shù)等方面比四面體網格具有明顯優(yōu)勢，因此根據(jù)三維實體模型幾何形狀的復雜性，研究在任意三維空間內自動生成過渡均勻、協(xié)調、自適應的六面體網格，建立可靠、高效的六面體網格離散算法

3、，可以為金屬成形、巖土工程、流體力學、生物醫(yī)學工程等領域分析大型復雜科學與工程問題建立三維模型離散平臺，直接推動工程分析軟件的商品化進程。 本文以基于柵格法六面體網格生成方法作為基礎方法，開展了基于實體模型幾何形狀的六面體網格自適應生成方法的研究。建立了基于STL文件進行三維實體模型幾何特征自動識別的關鍵技術；提出了基于實體模型幾何形狀的六面體網格單元密度信息場的生成方法、加密準則和加密模板，將基于實體模型幾何形狀的網格自適應生

4、成技術直接應用于柵格法六面體初始網格的生成；在六面體初始網格自適應生成方法研究的基礎上，研究了六面體網格再生成過程的關鍵技術；建立了三維六面體網格質量評價準則，研究給出了提高網格單元質量的優(yōu)化方法，主要包括網格質量評價模型的建立、網格單元拓撲關系優(yōu)化、網格單元形狀質量優(yōu)化等；研究了三維網格顯示消隱技術、三維剖切面算法等關鍵技術，應用計算機圖形學知識開發(fā)了三維六面體網格自適應生成軟件AUTOMESH-3D。 首先，從力學角度和計算

5、角度并以金屬體積成形為例論述了三維有限元分析中六面體單元的優(yōu)越性。對典型的六面體網格生成方法進行了總結對比，得出基于柵格法六面體網格生成方法具有高度自動化和易于網格局部加密的優(yōu)點，適合于六面體網格自適應生成和再生成。在此基礎上，本文選擇基于柵格法的六面體網格自適應生成方法作為深入研究的方法。采用CAD造型軟件生成的STL文件傳遞實體模型的表面幾何形狀信息，并建立了新的具有拓撲關系的STL數(shù)據(jù)文件，通過計算實體模型表面三角形面片的曲率，實

6、現(xiàn)了幾何特征的精確識別，為獲得能夠全面反映實體模型幾何特征的六面體網格奠定了基礎；提出了基于實體模型的幾何特征建立單元密度信息場技術，給出了基于實體模型表面曲率和局部厚度建立加密源點信息場和單元密度信息場的具體方法和步驟；為保證生成過渡均勻的協(xié)調六面體網格，提出了一套六面體網格加密過渡模板。與Schneiders等人提出的加密模板相比較，本文改進了面加密模板和點加密模板，減少了一次加密后網格單元和節(jié)點的數(shù)目，并增加了兩個拐角模板，有效地

7、控制了加密區(qū)域的蔓延現(xiàn)象。 將幾何自適應技術與基于柵格法六面體網格生成方法相結合，提出了基于實體模型幾何特征的柵格法六面體網格自適應生成方法。主要研究了從外向內柵格法六面體網格生成方法和從內向外柵格法六面體網格生成方法的基本算法和關鍵技術。鑒于目前常用的“最近距離反復搜索法”不能滿足復雜網格模型的邊界擬合，本文提出了穿插法邊界擬合方法，并建立了相應的擬合規(guī)則。研究給出了基于表面四邊形單元的六面體單元空隙填充法，提出了基于雅可比矩

8、陣的最優(yōu)邊界擬合方法，基于網格邊界單元與實體模型特征邊界的相對位置關系，實現(xiàn)了從內向外柵格法六面體網格與實體模型邊界的精確擬合，達到了因地制宜的擬合效果。 網格再劃分是金屬體積成形等大變形模擬中不可缺少的過程，本文在研究六面體初始網格自適應生成的基礎上，深入研究了六面體網格再劃分技術。主要包括六面體網格再劃分的兩個主要標準：工件與模具幾何離散造成的干涉和工件網格單元的變形，并給出了兩個標準的建議值；提出并建立了新舊網格之間包容測

9、試和方程迭代求解物理量傳遞的具體算法，經汽車轉向節(jié)鍛壓成形過程等工程實例計算表明，該算法具有計算速度快、物理量傳遞準確的特點。 網格的質量直接影響著數(shù)值模擬結果的精度和效率，本文針對三維六面體網格質量評價準則和優(yōu)化技術進行了研究。建立了基于六面體單元雅可比矩陣行列式的值和條件數(shù)的網格質量評價準則。單元的雅可比矩陣能夠全面地反映網格的質量，包括角度、面積、長度等網格質量指標。評價準則量化了六面體網格質量，使網格質量的檢測更加簡便和

10、精確，確保獲得符合有限元數(shù)值分析要求的六面體網格單元；研究分析了基于柵格法生成六面體網格的拓撲關系，提出了插入新單元技術和單元退化技術，顯著改善了網格單元間的拓撲關系，為網格形狀質量優(yōu)化奠定了基礎。提出了六種插入新單元模式和三種單元退化模式，分別適合于六面體網格凸特征邊界和凹特征邊界單元拓撲關系的調整；研究提出了基于曲率拉普拉斯邊界優(yōu)化方法，該方法不僅能夠實現(xiàn)網格邊界質量的優(yōu)化，而且保證了優(yōu)化后網格邊界單元很好地描述實體模型的邊界幾何特

11、征；提出了基于單位化雅可比矩陣條件數(shù)為目標函數(shù)的六面體網格表面和內部節(jié)點優(yōu)化方法，建立了基于四邊形單元的雅可比矩陣條件數(shù)的目標函數(shù)，實現(xiàn)了網格表面節(jié)點質量的優(yōu)化，優(yōu)化后網格保留了原網格的表面幾何特征，確保獲得的六面體網格符合有限元數(shù)值模擬計算的要求；建立了通用的三維六面體網格質量評價準則，并將其與網格質量優(yōu)化技術相結合，采用C++語言編寫了三維六面體網格質量優(yōu)化程序，將其應用于基于柵格法自適應生成的六面體網格實例中，驗證了本文提出算法的

12、有效性和開發(fā)的六面體網格優(yōu)化程序的準確性。 深入研究了計算機圖形學基礎及OpenGL技術和可視化技術的基本原理，并給出了自主研發(fā)的六面體網格自適應生成軟件AUTOMESH-3D的結構及其功能特點。主要通過設計合理的輸入輸出數(shù)據(jù)接口，實現(xiàn)了與三維CAD造型軟件和DEFORM-3D等數(shù)值模擬軟件的數(shù)據(jù)銜接；針對基于三維六面體網格的消隱技術、三維剖切算法進行了研究，實現(xiàn)了任意剖切面上物理量的可視化，方便了用戶的操作；通過網格質量檢測模