基于ANSYS經(jīng)典界面的帶孔薄板的自適應(yīng)網(wǎng)格劃分

我們知道,有限元技術(shù)的基本思想是用分片插值來逼近真實(shí)的函數(shù)。在一般情況下,真實(shí)函數(shù)是什么樣的這很難知道,但是通過縮小單元尺寸可以越來越逼近真實(shí)解。上述思想在有限元分析中就體現(xiàn)在加密網(wǎng)格以得到精確解。

為了進(jìn)行網(wǎng)格加密,一般有限元軟件提供了具備網(wǎng)格加密的方式,有些也提供了自適應(yīng)網(wǎng)格劃分方法來幫助用戶迅速找到收斂解。

本文說明ANSYS是如何使用自適應(yīng)網(wǎng)格劃分技術(shù)來自動得到收斂解的。

【問題描述】

一塊帶孔方板,一端固定,另外一端施加分布壓力,要求其中的應(yīng)力分布。

已知方板長200mm,寬100mm,在正中間鉆一通孔,半徑為10mm。

【問題分析】

• 靜力學(xué)問題,平面應(yīng)力,最簡單的線彈性問題。

• 為了得到問題的真實(shí)解,需要細(xì)分網(wǎng)格。

• 如果對整體細(xì)分網(wǎng)格,則會面臨一個問題:在左邊固定端的上下角點(diǎn)處,由于應(yīng)力集中,此處的應(yīng)力會隨著網(wǎng)格劃分細(xì)密而無限增大。真正應(yīng)該細(xì)分的是中間空周圍。

• 如果是手工劃分方式來細(xì)分網(wǎng)格,也是可行的。不過這要手工細(xì)分多次,這里使用自適應(yīng)網(wǎng)格劃分方法來自動細(xì)分網(wǎng)格。

• 為了只對中間關(guān)注區(qū)域進(jìn)行網(wǎng)格細(xì)分,這里把整個面分為三部分,然后選中中間一塊面,對它進(jìn)行自適應(yīng)網(wǎng)格劃分。

• 計(jì)算的結(jié)果就是收斂的結(jié)果,不需要再次細(xì)分網(wǎng)格。

1. 建模

1.1創(chuàng)建單元類型,設(shè)置材料模型

/PREP7
ET,1,PLANE42
MP,EX,1,2e11
MP,PRXY,1,0.3

上述命令確定用PLANE42單元來建模,并給出了彈性模量和泊松比。

1.2繪制帶孔方板

RECTNG,0,0.03,0,0.1
RECTNG,0.03,0.17,0,0.1
RECTNG,0.17,0.2,0,0.1
AGLUE,ALL
CYL4,0.1,0.05,0,,0.01
ASBA,4,2

上述命令首先以此創(chuàng)建三個相連的面,然后把它們粘接在一起,接著在中間創(chuàng)建一個圓面,然后減去該圓面,從而得到帶孔的方板模型。

這里把整體劃分為三個部分,是為了只對中間部分進(jìn)行網(wǎng)格細(xì)分,因?yàn)榛诹W(xué)理論我們知道,中間圓形周圍才是最危險處。如果我們對整體細(xì)分網(wǎng)格,會發(fā)現(xiàn)固定端角落點(diǎn)的應(yīng)力會無限上升,這將使得自適應(yīng)網(wǎng)格劃分失去意義。

1.3劃分網(wǎng)格

SMRT,OFF
AMESH,ALL

上述命令對整體劃分網(wǎng)格。網(wǎng)格劃分結(jié)果如下圖

1.4設(shè)置邊界條件

DL,4,,ALL
SFL,10,PRES,10E6
FINISH

上述命令固定了左端面,然后對右邊施加10MPa的均布壓力。結(jié)果如下圖

2. 求解

/SOLU
ASEL,S,,,3
ADAPT,10,5,,0.25,2

上述命令選擇方板中間區(qū)域進(jìn)行自適應(yīng)網(wǎng)格劃分,并啟動靜力學(xué)分析。

在分析過程中,宏ADAPT會根據(jù)指定的誤差5%,來進(jìn)行迭代計(jì)算。當(dāng)前后兩次的能量誤差小于該值時,網(wǎng)格細(xì)分停止,此時認(rèn)為結(jié)果已經(jīng)收斂。

3. 后處理

/POST1
ALLS
PLNSOL,S,EQV,0,1

上述命令查看等效應(yīng)力云圖。

可見,中間孔上下邊緣最危險,應(yīng)力達(dá)到29.5Mpa,該值是收斂值。

如果我們是為公司做項(xiàng)目分析,可以出具該值作為最終結(jié)果。

【結(jié)論】

對于每一個有限元分析來說,都需要確保計(jì)算收斂。很多有限元分析案例,只是進(jìn)行了一次網(wǎng)格劃分,然后給出了結(jié)果,就認(rèn)為結(jié)果是正確的。對于這種做法是高度懷疑的。從仿真實(shí)踐中我們發(fā)現(xiàn),當(dāng)網(wǎng)格細(xì)分時,有時候結(jié)果會發(fā)生很大的改變,在沒有明確的理論根據(jù)的前提下,就把一次網(wǎng)格劃分的結(jié)果認(rèn)定為最終結(jié)果,這是很難令人信服的。

ANSYS的自適應(yīng)網(wǎng)格劃分解決了這個難題。不過該技術(shù)還存在諸多限制。例如只能用于線性靜力學(xué)結(jié)構(gòu)分析和線性穩(wěn)態(tài)熱分析等,不過這種限制,隨著ANSYS版本的更新,在逐漸減少。希望隨著ANSYS的發(fā)展,最終能夠?qū)τ谌我獾姆治龆寄軌蜃龅竭@一點(diǎn),這對于用戶來說無疑是相當(dāng)重要的,我們翹首企盼好了。