Abaqus有限元分析實(shí)例解析——在棒中的應(yīng)力波傳播(三)

2013-08-14  by:有限元軟件應(yīng)用培訓(xùn)中心  來源:仿真在線

Abaqus有限元分析實(shí)例解析——在棒中的應(yīng)力波傳播(三)

 

9.4.3 網(wǎng)格對穩(wěn)定時間增量和CPU時間的影響

在第9.3節(jié)“自動時間增量和穩(wěn)定性”中,我們討論過網(wǎng)格細(xì)劃對穩(wěn)定極限和CPU時間的影響。這里我們將以波的傳播問題來說明這一影響。我們從方形單元的一種合理的精細(xì)網(wǎng)格開始,沿長度方向劃分50個單元,兩個橫向方向各劃分10個單元。為了說明問題,我們現(xiàn)在采用一種25×5×5單元的粗糙網(wǎng)格,并觀察在各種方向上如何細(xì)劃網(wǎng)格改變CPU時間。四種網(wǎng)格如圖9-11所示。

 

表9-1顯示了本問題的CPU時間隨著網(wǎng)格細(xì)劃的改變(以粗糙網(wǎng)格的模型結(jié)果進(jìn)行了單位化)?;谠诒局改现薪榻B的簡單的穩(wěn)定性方程,表格的前一半提供了期望值;表格的后一半給出了在計算機(jī)工作站上由ABAQUS/Explicit運(yùn)行分析得到的結(jié)果。

 

表 9-1 網(wǎng)格細(xì)劃和求解時間

網(wǎng)格

簡化理論

實(shí)際

 

()

單元數(shù)

CPU時間()

最大()

單元數(shù)

CPU時間()

25×5×5

A

B

C

6.06e-6

625

1

50×5×5

A/2

2B

4C

3.14e-6

1250

4

50×10×5

A/2

4B

8C

3.12e-6

2500

8.33

50×10×10

A/2

8B

16C

3.11e-6

5000

16.67

 

對于理論解答,我們選擇25×5×5的最粗糙網(wǎng)格作為基本狀態(tài),并且定義穩(wěn)定的時間增量、單元數(shù)量和CPU時間分別為變量A、B和C。隨著網(wǎng)格的細(xì)劃,產(chǎn)生了兩種結(jié)果:最小的單元尺寸減小了,以及在網(wǎng)格中的單元數(shù)目增加了,這些影響的每一種都會增加CPU時間。在第一次細(xì)劃的50×5×5網(wǎng)格中,最小單元的尺寸減小了一半,并且單元的數(shù)目增加了一倍,因此增加了CPU時間,使其是前一種網(wǎng)格的四倍。然而,進(jìn)一步加倍網(wǎng)格數(shù)目至50×10×5,沒有改變最小單元的尺寸;而僅僅是加倍了單元數(shù)量。因此,CPU時間是50×5×5網(wǎng)格的二倍。進(jìn)一步地細(xì)劃網(wǎng)格為50×10×10,使單元成為均勻的方形,再一次加倍了單元的數(shù)量和CPU時間。

這種簡單的計算非常好地預(yù)測了網(wǎng)格細(xì)劃如何影響穩(wěn)定時間增量和CPU時間的趨勢。然而,為什么我們沒有將預(yù)測值與實(shí)際的穩(wěn)定時間增量值進(jìn)行比較?這里是有原因的。首先回憶一下我們給出的穩(wěn)定時間增量的近似公式為

                                                         Abaqus有限元分析實(shí)例解析——在棒中的應(yīng)力波傳播(三)ansys結(jié)果圖圖片2

然后我們假定單元特征長度是最小的單元尺寸,而ABAQUS/Explicit實(shí)際上是根據(jù)單元的整體尺寸和形狀來確定單元特征長度的。另外一個因素是ABAQUS/Explicit采用了一個整體穩(wěn)定性估算,它允許使用一個更大的穩(wěn)定時間增量。這些因素使得在運(yùn)行分析之前實(shí)際上難以準(zhǔn)確地預(yù)測穩(wěn)定時間增量。然而,由于預(yù)測的趨勢與簡單的理論符合得很好,因此可以直接地預(yù)測穩(wěn)定時間增量如何隨著網(wǎng)格細(xì)劃而發(fā)生變化。

9.4.4 材料對穩(wěn)定時間增量和CPU時間的影響

同樣的波擴(kuò)展分析在不同的材料中進(jìn)行將需要不同的CPU時間,這取決于材料的波速。例如,如果我們改變材料從鋼到鋁,波速將從5.15×103 m/s變成為

                                    Abaqus有限元分析實(shí)例解析——在棒中的應(yīng)力波傳播(三)ansys結(jié)果圖圖片3

因為剛度和密度幾乎改變了相同的量級,所以從鋁到鋼對穩(wěn)定時間增量只有微小的影響。在鉛的情況下,差別則變得非常大,其波速減小為

                                    Abaqus有限元分析實(shí)例解析——在棒中的應(yīng)力波傳播(三)ansys結(jié)果圖圖片4

這個值大約是鋼的波速的五分之一。鉛棒的穩(wěn)定時間增量將是鋼棒的穩(wěn)定時間增量的五倍。

 


開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:Abaqus有限元分析實(shí)例解析——在棒中的應(yīng)力波傳播(三) 廣州有限元課程 Abaqus分析 有限元設(shè)計仿真 Ansys有限元培訓(xùn) Ansys workbench培訓(xùn) ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術(shù)咨詢 ansys基礎(chǔ)知識 ansys代做 

編輯
在線報名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進(jìn)行在線報名,我們會及時給您回復(fù)電話,謝謝!
驗證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com