格劃分:掃掠、多區(qū)與膨脹層實戰(zhàn)指南)
簡介這份PDF面向結構、流體、電磁等領域的仿真工程師與CAE初學者系統(tǒng)講解ANSYS Meshing高級網(wǎng)格劃分技術幫助解決復雜幾何模型網(wǎng)格質(zhì)量差、劃分效率低等實際問題。內(nèi)容由安世亞太網(wǎng)格工程師程東明主講涵蓋網(wǎng)格劃分流程、三維與二維網(wǎng)格方法、全局與局部網(wǎng)格控制、網(wǎng)格質(zhì)量檢查等模塊并深入剖析四面體Patch Conforming與Patch Independent算法、掃略與多區(qū)網(wǎng)格劃分、六面體為主及自動網(wǎng)格劃分的適用場景。資源包共1個PDF文件大小約1.49MB便于隨時查閱與離線學習。目前已有108人學習。讀者可從中掌握不同物理場下的全局單元尺寸、相關性、關聯(lián)中心、平滑與過渡等參數(shù)設置思路理解力學、CFD、電磁與顯式分析的網(wǎng)格差異并借助局部控制與質(zhì)量檢查方法提升仿真精度與效率適合需要進階網(wǎng)格劃分技能的產(chǎn)品設計與結構優(yōu)化人員參考。1. 從網(wǎng)格畫不動說起ANSYS Meshing 高級技術到底解決什么問題很多人第一次在 ANSYS Workbench 里點開 Mesh畫完四面體就點 Solve結果要么求解器報負體積要么應力集中區(qū)域的結果明顯失真。問題往往不在求解設置而在網(wǎng)格本身。ANSYS Meshing 高級網(wǎng)格劃分技術要解決的就是幾何復雜、局部梯度大、純自動網(wǎng)格搞不定這三類場景裝配體里薄壁與實體混在一起、流場里邊界層需要貼體加密、結構件上圓角與焊縫需要局部控制。它適合已經(jīng)會點自動劃分、但結果精度或收斂性卡住的人。核心思路不是把網(wǎng)格畫得更密而是把網(wǎng)格畫得更懂幾何——用局部尺寸、映射、掃掠、膨脹層這些手段把計算資源堆到真正需要的地方。下面按選型、實現(xiàn)、實戰(zhàn)、排錯的順序把能復現(xiàn)的操作講清楚。2. ANSYS Meshing 的網(wǎng)格方法與全局控制參數(shù)怎么選2.1 五種網(wǎng)格劃分方法的適用邊界ANSYS Meshing 里能選的方法不止一種選錯方法后面怎么調(diào)都別扭。常見的有自動Automatic、四面體Tetrahedrons、六面體主導Hex Dominant、掃掠Sweep、多區(qū)MultiZone。判斷邏輯是看幾何能不能掃出來如果體是沿一條路徑拉伸或旋轉生成的優(yōu)先掃掠如果是多個規(guī)則塊拼起來的用多區(qū)實在不規(guī)則的復雜實體才退回四面體加膨脹層。方法適用幾何單元類型典型場景Automatic任意混合快速試算Tetrahedrons復雜實體四面體鑄造件、復雜支架Hex Dominant較規(guī)則實體六面體為主一般結構件Sweep可掃掠體六面體軸、管、葉片MultiZone可分解塊體六面體規(guī)則裝配體選方法時先問一句這個體能不能用一條源面沿路徑掃到目標面能就別用四面體。四面體單元剛度偏大同樣的節(jié)點數(shù)精度不如六面體而且邊界層不好控制。2.2 全局尺寸與曲率、鄰近度加密全局參數(shù)在 Mesh 的 Details 里最常動的是 Sizing 下的幾個。Element Size 是基準尺寸但真正決定局部密度的是 Curvature 和 Proximity 兩個捕獲函數(shù)。# Workbench 里對應設置腳本化時用 ACT 或 Python 調(diào)用 # Sizing - Size Function Curvature # Curvature Normal Angle 18 # 每 18 度弧長至少一個單元 # Proximity Accuracy 0.5 # 鄰近度捕獲靈敏度 # Num Cells Across Gap 3 # 窄縫里至少鋪 3 層單元Curvature Normal Angle 越小圓角、孔壁上的單元越密默認 18 度對多數(shù)結構夠用做疲勞分析可以壓到 12 度。Num Cells Across Gap 是薄壁和窄縫的關鍵默認 3如果兩個面之間只鋪了 1 層單元求解器很容易在這里出問題調(diào)到 4 到 5 更穩(wěn)。Proximity Accuracy 越大越敏感但太大會讓網(wǎng)格量暴漲一般 0.5 到 0.7 之間。提示全局尺寸改小之前先用局部尺寸把關鍵區(qū)域控住否則整體網(wǎng)格量翻倍求解時間成倍漲收益卻集中在少數(shù)區(qū)域。2.3 局部尺寸與命名選擇的配合全局參數(shù)定基調(diào)局部尺寸定重點。做法是在模型樹里選中面、邊或體右鍵 Insert - Sizing然后指定 Element Size 或 Number of Divisions。更高效的方式是先在幾何里給關鍵面命名選擇Named Selection再在 Meshing 里按命名選擇批量加局部控制。# 通過命名選擇批量設置局部尺寸的典型流程示意 # 1. 幾何中創(chuàng)建 Named Selection: bolt_holes # 2. Meshing 中 Insert - Sizing # 3. Geometry 選擇該 Named Selection # 4. Type Element Size, Element Size 0.5 mm # 5. Behavior Soft, 允許過渡Behavior 選 Soft 時局部尺寸會和周圍網(wǎng)格平滑過渡選 Hard 則強制該區(qū)域嚴格按尺寸走過渡處可能出現(xiàn)尺寸突變。螺栓孔、圓角這類應力集中區(qū)建議 Soft配合全局曲率捕獲而接觸面、對稱面這種需要單元對齊的地方用 Hard 更可控。3. 掃掠、多區(qū)與膨脹層把六面體和邊界層真正畫出來3.1 掃掠網(wǎng)格的源面、目標面與路徑設置掃掠是拿到六面體最直接的路子。在 Mesh 里選中體Method 設為 Sweep然后指定 Source 面和 Target 面。源面和目標面拓撲要一致否則掃不出來。如果體是沿曲線路徑生成的還要指定 Sweep Path讓單元沿路徑分布。# Sweep 方法關鍵參數(shù) # Method Sweep # Src/Trg Selection Manual # Source face_A # Target face_B # Free Face Mesh Type All Quad # 源面全四邊形 # Sweep Num Divs 20 # 沿路徑層數(shù)Free Face Mesh Type 選 All Quad 是拿到純六面體的前提如果源面被迫用三角形掃出來就是楔形單元。Sweep Num Divs 控制沿路徑的層數(shù)對薄壁件來說這個值直接決定厚度方向有幾層單元一般不少于 3 層。掃不出來時先檢查源面和目標面是否真的拓撲對應常見坑是幾何里有小碎面或倒角導致源面被切成多塊。3.2 多區(qū)方法處理不規(guī)則裝配體多區(qū)MultiZone適合那種整體不規(guī)則、但能拆成幾個規(guī)則塊體的幾何。它內(nèi)部會自動做塊分解能分解成功的區(qū)域給六面體分解不了的退回四面體。# MultiZone 關鍵設置 # Method MultiZone # Mapped Mesh On # 盡量映射 # Free Mesh Type Tetra # 無法映射時用四面體 # Src/Trg Selection AutomaticMapped Mesh 打開后能映射的區(qū)域會強制映射網(wǎng)格質(zhì)量更整齊。如果幾何里有大量細小特征導致塊分解失敗可以先在 DM 或 SpaceClaim 里做幾何清理把不影響結果的倒角、小孔抑制掉再回來畫多區(qū)。這一步的收益往往比在 Meshing 里反復調(diào)參數(shù)大得多。3.3 膨脹層參數(shù)與邊界層貼體做流體或傳熱時壁面附近的邊界層必須用膨脹層Inflation貼體加密。在 Mesh 里右鍵 Insert - Inflation選壁面然后設 First Layer Height、Growth Rate 和 Maximum Layers。參數(shù)含義常用取值First Layer Height第一層厚度按 y 目標反算Growth Rate相鄰層厚度比1.1 ~ 1.3Maximum Layers最大層數(shù)5 ~ 15Inflation Option總厚度控制方式Total Thickness / Smooth TransitionFirst Layer Height 不是拍腦袋定的要先估 y。做 k-epsilon 模型時壁面函數(shù)要求 y 在 30 到 300做 SST 或低雷諾數(shù)模型時要壓到 1 附近。Growth Rate 超過 1.3 容易在層間產(chǎn)生大長寬比單元影響收斂。Maximum Layers 要結合幾何厚度別讓膨脹層互相穿透兩個壁面靠近時用 Total Thickness 并設一個小于間隙一半的值。注意膨脹層畫完后一定要看 Skewness 和 Aspect Ratio邊界層單元長寬比大是正常的但 Skewness 超過 0.95 的區(qū)域要單獨處理。4. 網(wǎng)格質(zhì)量診斷與常見報錯的排查路徑4.1 用 Mesh Metrics 定位壞單元畫完網(wǎng)格別急著求解先看 Mesh Metrics。在 Mesh 下 Insert - Mesh Metric常用的是 Element Quality、Skewness、Aspect Ratio、Jacobian Ratio。Element Quality 取值 0 到 1越接近 1 越好低于 0.1 的單元基本要處理。Skewness 越低越好超過 0.95 的單元在求解時容易出問題。# 查看并定位最差單元 # Mesh - Insert - Mesh Metric - Skewness # 在圖形窗口下方滑塊拖動只顯示 Skewness 0.9 的單元 # 右鍵這些單元 - Go To定位到幾何位置定位到壞單元后處理思路分兩種如果壞單元在應力集中區(qū)加局部尺寸或改方法如果在不重要的角落直接抑制該幾何特征。不要為了幾個角落單元把全局尺寸改小那樣代價太大。4.2 負體積與求解器報錯的對應關系求解器報負體積negative volume或高度扭曲單元根子幾乎都在網(wǎng)格。常見對應關系是負體積多出現(xiàn)在大變形區(qū)域說明該處單元太扁或太斜收斂困難但沒報負體積往往是邊界層長寬比過大或接觸面網(wǎng)格不匹配。# 排查順序 # 1. 檢查 Skewness 0.95 的單元分布 # 2. 檢查接觸對兩側網(wǎng)格尺寸是否接近 # 3. 檢查膨脹層是否在窄縫處穿透 # 4. 檢查大變形區(qū)是否有足夠層數(shù)接觸對兩側網(wǎng)格尺寸差太多時接觸搜索會不穩(wěn)定。做法是給接觸面加匹配的局部尺寸讓兩側單元尺度接近。窄縫處膨脹層穿透把 Maximum Layers 降下來或改用 Total Thickness 控制總厚度。4.3 網(wǎng)格無關性驗證的最小做法網(wǎng)格畫完不代表結束還要做網(wǎng)格無關性驗證。做法是同一模型用三套不同密度的網(wǎng)格求解比較關鍵結果最大應力、最大位移、壓降等。如果兩套較密網(wǎng)格的結果差異小于 5%說明當前密度夠用。網(wǎng)格套數(shù)全局尺寸關鍵結果相對偏差粗2.0 mm基準—中1.4 mm接近 5%細1.0 mm接近 2%三套網(wǎng)格不要只改全局尺寸關鍵區(qū)域的局部尺寸要同步縮放否則驗證的是全局密度而不是關鍵區(qū)密度。做疲勞或斷裂分析時關鍵區(qū)網(wǎng)格要單獨加密到結果穩(wěn)定為止。5. 進階技巧用命名選擇與參數(shù)化把網(wǎng)格流程固化下來真正把 ANSYS Meshing 用順的人最后都會把網(wǎng)格設置參數(shù)化。做法是把關鍵尺寸、層數(shù)、增長率提成參數(shù)在 Workbench 的 Parameter Set 里統(tǒng)一管理幾何一改網(wǎng)格跟著更新不用每次重畫。# 參數(shù)化網(wǎng)格的典型參數(shù) # P1: global_size 2.0 # P2: hole_size 0.5 # P3: inflation_layers 8 # P4: growth_rate 1.2 # 在 Meshing 里把這些值綁定到對應輸入框 # 幾何更新后 Rebuild網(wǎng)格自動重劃配合命名選擇可以把哪些面加密、哪些面做膨脹層這套邏輯固定下來。幾何迭代時只要命名選擇還在網(wǎng)格控制就不會丟。另一個實用技巧是保存網(wǎng)格設置模板在 Meshing 里調(diào)好一套參數(shù)后右鍵 Mesh - Save Mesh Settings下次同類模型直接 Load省掉重復勞動。驗證網(wǎng)格是否真的可用最直接的辦法是看求解器前幾步的殘差和監(jiān)測量。如果殘差在前幾步就震蕩或監(jiān)測量跳變先回來看網(wǎng)格而不是去調(diào)松弛因子。把網(wǎng)格質(zhì)量、邊界層、接觸匹配這三件事做扎實后面求解設置才有意義。本文還有配套的精品資源點擊獲取