網(wǎng)格收斂的實例教程
摘要:針對車載導(dǎo)航的模態(tài)分析,以其中單個零件的網(wǎng)格尺寸作用變量,而其余零件的網(wǎng)格尺寸不變,然后進行模態(tài)分析,通過對計算結(jié)果的比較,確定模型最終的網(wǎng)格收斂尺寸,解決計算精度與計算時間成本這一矛盾。
關(guān)鍵字:模態(tài)分析;網(wǎng)格尺寸;收斂;固有頻率;車載導(dǎo)航
一、前言
由于公司產(chǎn)品做的最多的是動力學(xué)分析,動力學(xué)分析需要基于模態(tài)特征值上進行計算的,所以我們以最有代表性的模態(tài)分析作為網(wǎng)格收斂尺寸研究課題的分析類型。
相對于靜力分析而言,動力學(xué)分析的模態(tài)求解過程所需要的CPU計算時間更長,我們一般會通過控制網(wǎng)格的尺寸來調(diào)整計算時間成本,即粗化網(wǎng)格。但是如果網(wǎng)格太粗,則網(wǎng)格模型必會忽略或簡化結(jié)構(gòu)上的細(xì)節(jié),從而影響計算的精度。
面對計算精度與計算時間成本這一矛盾,則需要對模型的網(wǎng)格收斂尺寸、計算耗時等數(shù)據(jù)進行研究。
二、計算模型
計算模型采用某一車載導(dǎo)航屏模型,零件包括觸摸屏玻璃、前后面板、SECC支架,另外TFT顯示屏簡化為只包含質(zhì)量的點。
三、網(wǎng)格尺寸控制
導(dǎo)航產(chǎn)品的前后面板為PC+ABS注塑件,結(jié)構(gòu)特征多且比較復(fù)雜,為了盡可能詳盡地模擬結(jié)構(gòu)的細(xì)節(jié),如圓角、肋骨、凹凸臺、小圓孔等),縮小網(wǎng)格模型與原幾何模型的差距,因此網(wǎng)格不能過于粗糙。此外,根據(jù)模型的整體尺寸,最終確定模型的網(wǎng)格研究尺寸為1mm、2mm、3mm、4mm。這四種網(wǎng)格尺寸中,以所有零件的網(wǎng)格尺寸都為3mm的計算結(jié)果作為不變量,而其中單一個零件的網(wǎng)格尺寸作為變量,變量值即為1mm、2mm、3mm、4mm。
部分網(wǎng)格模型圖如下:
四、模型建立及求解的軟件
建立三維有限元模型,用Lanzos方法求解特征值。
展開 一 網(wǎng)格收斂的重要性
由于應(yīng)力集中(區(qū)別于應(yīng)力奇異)的存在,在結(jié)構(gòu)不連續(xù)處存在應(yīng)力較大,而且隨著網(wǎng)格質(zhì)量數(shù)量的增加,應(yīng)力值趨于收斂,據(jù)說收斂與否的應(yīng)力差值可以很大,所以說重要細(xì)節(jié)結(jié)構(gòu)的網(wǎng)格收斂十分重要。
二 WorkBench中網(wǎng)格收斂的實現(xiàn)
WorkBench中在solution選項中設(shè)置網(wǎng)格循環(huán)次數(shù),關(guān)鍵點網(wǎng)格優(yōu)化系數(shù),在求解結(jié)果選項下插入convergence,定義deformation或者stress的收斂系數(shù)。
雖然六面體網(wǎng)格變成四面體網(wǎng)格進行細(xì)分,但是初始網(wǎng)格劃分的尺寸,對結(jié)果仍然有一定影響。而且優(yōu)化的方式和APDL中也有一定差異,結(jié)果與APDL結(jié)果相同。
初始網(wǎng)格
網(wǎng)格收斂后
三 APDL中網(wǎng)格收斂的實現(xiàn)
(1) 建模,注意不要劃分網(wǎng)格,而且3D模型只能使用4面體單元網(wǎng)格;
(2) 加載邊界條件,由于沒有網(wǎng)格,邊界條件只好由面或者線確定;
(3) 啟動ADAPT宏命令,指定能量收斂誤差,最大循環(huán)次數(shù),網(wǎng)格優(yōu)化系數(shù);看起來很厲害的樣子,但是使用方法和命令一樣,只是輸入命令框的不提示有此命令存在;.
(4) 后處理查看結(jié)果。
四 問題描述
!APDL命令:
finish
/clear
/prep7
et,1,solid187 !
展開 這種精度對網(wǎng)格質(zhì)量的不敏感性支持 Mueller 的立場,即細(xì)胞質(zhì)量差是一個穩(wěn)定性問題。因此,STAR-CCM+ 的方法是保守的——選擇穩(wěn)健性而不是準(zhǔn)確性。具體來說,他們正在尋找將導(dǎo)致求解器中被零除的指標(biāo)。影響擴散通量和線性化的偏度就是這樣的一個例子。
Mesher 的觀點
John Steinbrenner 博士和 Nick Wyman 博士采用違反直覺的方法分享了 Pointwise 對與解決方案無關(guān)的質(zhì)量指標(biāo)的看法。您會認(rèn)為網(wǎng)格生成開發(fā)人員會提升先驗指標(biāo)的功效。但 CFD 解中的誤差包括幾何誤差、離散化誤差和建模誤差。幾何錯誤類似于 Dannenhoffer 和 Mueller 關(guān)于正確表示形狀的觀點。建模誤差來自湍流、化學(xué)和熱物理特性。離散化涉及求解器數(shù)值的退化。離散化誤差是由網(wǎng)格和求解器的數(shù)值算法之間的耦合驅(qū)動的。
圖 3:此表總結(jié)了 Pointwise 中可用的網(wǎng)格質(zhì)量指標(biāo)。從參考 1h。
因此,盡管 Fidelity Pointwise 可以計算和顯示許多指標(biāo),但需要注意的是,其中許多指標(biāo)與求解器的數(shù)值沒有直接關(guān)系,因此它們只是解決方案準(zhǔn)確性的松散指標(biāo)。另一方面,這些指標(biāo)計算方便,可以解決 Dannenhoffer 的網(wǎng)格有效性問題,并提供啟動網(wǎng)格改進技術(shù)的機制。它們還構(gòu)成了用戶開發(fā)領(lǐng)域?qū)I(yè)知識的能力的基礎(chǔ)——與他們的特定應(yīng)用領(lǐng)域相關(guān)的指標(biāo)。
結(jié)論
CFD 求解器開發(fā)人員認(rèn)為網(wǎng)格質(zhì)量對收斂的影響遠(yuǎn)大于精度。因此,由于收斂不良或不完全收斂而導(dǎo)致的求解誤差不容忽視。
一位研究人員能夠證明網(wǎng)格質(zhì)量與求解精度之間完全沒有相關(guān)性。為其他求解器和流動條件重現(xiàn)此結(jié)果將很有價值。
使用盡可能多的網(wǎng)格點 (Dannenhoffer, McDaniel)。在許多情況下,分辨率勝過質(zhì)量。
展開 <h1><strong>一、收斂的定義和重要性</strong></h1><h2><strong style="color: rgb(51, 51, 51);">1.收斂的多種含義</strong></h2><p>在有限元分析中,收斂具有多重意義。它包括網(wǎng)格收斂、時間積分精度和非線性程序收斂。</p><p><strong>l 網(wǎng)格收斂</strong>是指增加模型單元數(shù)量會使仿真解趨于解析解。對于線性和非線性問題都適用,AbaqUS 中使用 H 網(wǎng)格自適應(yīng)技術(shù)來輔助實現(xiàn)網(wǎng)格收斂。當(dāng)進一步加密網(wǎng)格時,結(jié)果變化很小或不變時,可認(rèn)為網(wǎng)格達(dá)到收斂。但也存在一些例外情況,如網(wǎng)格奇異解或材料損傷累積在模型特定區(qū)域的局部問題。</p><p><strong>l 時間積分精度</strong>則是針對具有物理時間尺度的瞬態(tài)問題,AbaqUS 提供用戶定義參數(shù),以控制對相關(guān)方程的積分精度。</p><p><strong>l 非線性程序收斂</strong>是本文重點討論的內(nèi)容,要獲得精確解需要滿足網(wǎng)格收斂、瞬態(tài)問題的精確時間積分以及非線性求解過程收斂等條件。</p><h2><strong style="color: rgb(51, 51, 51);">2.收斂對分析結(jié)果的影響</strong></h2><p>收斂性直接關(guān)系到分析結(jié)果的準(zhǔn)確性。如果模型不收斂,得到的結(jié)果可能毫無意義,甚至?xí)`導(dǎo)工程決策和學(xué)術(shù)研究。因此,理解和掌握 ABAQUS 中的收斂問題是正確使用該軟件進行有效分析的基礎(chǔ)。
展開 HLPW-4 /GMGW-3 中探索的高升力 CFD 分析
研討會結(jié)果表明,固定網(wǎng)格 RANS 通常無法實現(xiàn)足夠的網(wǎng)格收斂。在攤位附近,無法得出關(guān)于網(wǎng)格適用性的結(jié)論,即使是在最精細(xì)的網(wǎng)格上也是如此。對于高階方法,最佳實踐網(wǎng)格劃分指南與 RANS 完全不同。高度各向異性元素的三維網(wǎng)格彎曲是一項艱巨的任務(wù)。隱式求解器對于解決高階離散化至關(guān)重要。總的來說,將高階有限元和有限體積方法應(yīng)用于 CRM-HL 配置得到了證明,但非常具有挑戰(zhàn)性。
對于 HRLES,導(dǎo)致 RANS 網(wǎng)格收斂的網(wǎng)格設(shè)計不會導(dǎo)致 HRLES 網(wǎng)格收斂,并且網(wǎng)格分辨率被認(rèn)為對預(yù)測的分離模式有顯著影響。典型的 HRLES 運行成本大約是 RANS 的 10 到 15 倍。在一些風(fēng)洞試驗中,很難匹配所測隧道壁邊界層的形狀和厚度,這會影響求解。
對于壁模型大渦模擬和格子玻爾茲曼 (WMLESLB),在失速處和失速處發(fā)現(xiàn)對網(wǎng)格分辨率和網(wǎng)格特性(例如單元各向異性)的高度敏感。與 HRLES 一樣,網(wǎng)格分辨率對預(yù)測的分離模式有重大影響。典型的 WMLES 運行成本大約是 RANS 的 5 到 10 倍。與試圖代表風(fēng)洞試驗中的實際跳閘機制相比,數(shù)值跳閘更具成本效益。就平均時間而言,高攻角需要比低攻角更長的時間來確保實現(xiàn)平穩(wěn)性。
圖 2. 攻角 α = 7.05°(頂部)和 α = 19.57°(底部)時自由空氣 CL 預(yù)測的統(tǒng)計分析
結(jié)論
HLPW-4 和 GMGW-3 之間的聯(lián)合研討會模型使人們更加關(guān)注網(wǎng)格劃分,這繼續(xù)對 CFD 解決方案產(chǎn)生重大影響。高升力流的幾何準(zhǔn)備和固定網(wǎng)格劃分仍然很困難。為不同的方法、代碼和提升曲線區(qū)域確定固定網(wǎng)格指南也很困難。
展開