基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題)

2017-03-11  by:CAE仿真在線  來源:互聯(lián)網(wǎng)


如何在ANSYS中模擬鉆孔,隧道開挖,金屬切削這種類似的問題呢?

在金屬切削過程中,表層的金屬離開了表層,然后從模型中消失。對于這種問題,ANSYS提供了幾種支持,其中一種就是生死單元的方法。

使用這種方法,就是在適當?shù)臅r候,當滿足某種準則的時候,殺死那些不需要的單元,然后接著進行仿真。如果有必要的話,也可以在適當?shù)臅r候復活這些單元。

ANSYS中對于生死單元的例子很少。筆者遍歷了其幫助,也只找到一個例子。本文就同朋友們一起欣賞這個例子,看看ANSYS是如何使用生死單元來解決這種工程問題的。

問題描述

一根兩端固定的桿如下圖所示。

基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題)ansys仿真分析圖片1

材料數(shù)據(jù)如下

基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題)ansys仿真分析圖片2

為了闡述如何使用ANSYS的單元生死技術,決定把該桿等分為3個單元,然后通過控制中間單元的生死,進行如下的熱應力仿真

(1)設置所有單元的材料參考溫度是0度,給所有節(jié)點施加100度,并保持所有單元都存活,做1次仿真

(2)設置所有單元的材料參考溫度是0度,給所有節(jié)點施加100度,殺死中間單元,做1次仿真

(3)設置其它單元的材料參考溫度是0度,給所有節(jié)點施加100度,激活中間單元,并設置該單元的材料參考溫度是100度,做1次仿真

(4)設置其它單元的材料參考溫度是0度,給所有節(jié)點施加0度,保持中間單元存活,并設置該單元的材料參考溫度是100度,做1次仿真

通過上述四次仿真,以說明

(1)如何使用單元的生死技術

(2)當單元激活時,會根據(jù)節(jié)點溫度和該單元的材料參考溫度之差來確定它的初始熱應變。

問題分析

1.該例子來自于ANSYS15 APDL的認證算例《VM194 Element Birth/Death in a Fixed Bar》為了更清晰的闡明思路,本文對其進行了較大幅度的調(diào)整。

2.單元生死技術的使用,關鍵是首先要創(chuàng)建出所有的單元,然后在需要殺死改單元時使用EKILL命令,而在需要激活時使用ELIVE命令。

3.使用LINK180來建模桿。

4.創(chuàng)建2種材料。這兩種材料的彈性模量和泊松比一樣,但是參考溫度不一樣。一個參考溫度是0度,一個是100度。

5.先創(chuàng)建4個節(jié)點,然后創(chuàng)建3個單元。

6.固定兩個端節(jié)點,并給所有節(jié)點固定Z方向自由度,借此模擬二維桿件。

7.按照題目要求進行先后四次的計算和后處理,以考察生死單元的使用。

8.本文采用APDL命令進行講解。

求解過程

1. 建模

1.1 創(chuàng)建單元類型

在命令窗口輸入

/PREP7
ET,1,LINK180

上述命令首先進入到前處理器,然后創(chuàng)建桿單元LINK180,該單元用于模擬二力桿。

1.2 設置材料模型

在命令窗口輸入

MP,EX,1,30E6
MP,ALPX,1,.00005

上述命令定義了材料1的彈性模量和線膨脹系數(shù)。這里并沒有定義泊松比,它的默認值是0.3.

在命令窗口輸入

MP,EX,2,30E6
MP,ALPX,2,.00005
MP,REFT,2,100

上述命令定義了材料2,該材料的彈性模量和線膨脹系數(shù)與材料1一致,但是其參考溫度是100。該材料會用在下面的生死單元上面,設定100的參考溫度,是為了顯示單元復活時,它的初始應變是重新計算的。

1.3 設置截面

在命令窗口輸入

SECTYPE,1,LINK
SECDATA,1

上述命令首先定義了桿的截面是LINK,然后定義該LINK的截面積是1.

1.4 創(chuàng)建節(jié)點

在命令窗口輸入

N,1
N,4,10
FILL

上述命令創(chuàng)建了4個節(jié)點。結果如下圖

基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題)ansys仿真分析圖片3

1.5 創(chuàng)建單元

在命令窗口輸入

E,1,2
EGEN,3,1,-1

上述命令生成了三個單元。要注意,這三個單元都是用的第一種材料。換一句話說,它們的材料類型是一致的。創(chuàng)建完畢后結果如下圖:

基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題)ansys仿真分析圖片4

1.6 創(chuàng)建邊界條件

在命令窗口輸入

D,1,ALL,,,4,3
D,ALL,UZ

FINISH

上述命令

首先固定了邊界上兩個節(jié)點,使其成為兩端固定的桿件

然后對所有節(jié)點約束其Z方向的自由度

最后退出前處理器

結果如下圖

基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題)ansys結果圖片5

2. 第一次分析


在命令窗口輸入

/SOLU
ANTYPE,STATIC

NROPT,FULL
OUTPR,BASIC,ALL

上述命令進入到求解器,然后設置靜力學分析,并設置使用完全的牛頓-拉普森算法求解方程組,最后指出要輸出常用的數(shù)據(jù)。

在命令窗口輸入

TREF,0
TUNIF,100
SOLVE

上述命令

首先設置參考溫度是零度

接著給所有節(jié)點施加100度

然后求解。

由于節(jié)點溫度高于參考溫度,所以3個單元均有伸長的趨勢,可惜兩端被固定,于是產(chǎn)生了壓縮應變。

在命令窗口輸入

/post1
PRESOL,EPTH,COMP

上述命令進入到后處理器,并列出了各單元節(jié)點的熱應變。

可以看到,每個單元節(jié)點均處于被壓縮狀態(tài)。這是由于其伸長趨勢受阻礙而產(chǎn)生的。

3. 第二次分析

在命令窗口輸入

/SOLU
EKILL,2
SOLVE

上述命令

進入到求解器

殺死了2號單元,就是中間這個單元

然后求解

由于中間單元殺死,它的變形數(shù)據(jù)會消失,但是兩邊的單元仍然是存活的,它們依然有應變。

在命令窗口輸入
/post1
PRESOL,EPTH,COMP

上述命令進入到后處理器,并列出了各單元節(jié)點的熱應變。

基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題)ansys結果圖片6

可見,由于中間單元被殺死,它的熱應變就消失了。一個單元被殺死,就好像人死了一樣,無論這個世界如何變化,它不會再有任何反應。

要注意的是,此時每個節(jié)點溫度依然是100度,參考溫度仍然是0度。

4. 第三次分析

在命令窗口輸入

/SOLU
EALIVE,2

MPCHG,2,2

SOLVE

上述命令

首先進入處理器

然后讓2號單元復活

接著改變了2號單元的材料,使得它的材料是2號材料

然后開始求解

要注意的是,此時1,3號單元使用的是1號材料,而2號單元使用的是2號材料

2號材料的參考溫度是100度,而1號材料的參考溫度是0度,現(xiàn)在每個節(jié)點溫度依然是100度

這樣,雖然2號單元復活,但是因為它的節(jié)點溫度等于它所屬材料的參考溫度,因此它不會有應變;

而1,3號單元所屬的1號材料參考溫度0度,節(jié)點溫度100度,存在溫度差,所以這兩個單元依舊會有應變。

在命令窗口輸入
/post1
PRESOL,EPTH,COMP

上述命令進入到后處理器,并列出了各單元節(jié)點的熱應變。

基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題)ansys結果圖片7

可見,此時的結果與上一步驟的分析一樣,單元2依然沒有應變。

但是要注意,本次分析與上次分析有重要的不同。

本次分析2號單元是活的,只是因為情況特殊,它對外界的反應暫時是0!而上次分析單元是死的,無論外界是什么情況,它的反應永遠是0!

5. 第四次分析

在命令窗口輸入

/SOLU
TUNIF,0
SOLVE
上述命令

再次進入到求解器

然后設置所有節(jié)點溫度是0度

然后開始求解

此時。由于1,3號單元的節(jié)點溫度就等于參考溫度,所以這兩個單元沒有熱應變;而2號單元的節(jié)點溫度高于參考溫度,所以2號單元會有熱應變

在命令窗口輸入
/post1
PRESOL,EPTH,COMP

上述命令進入到后處理器,并列出了各單元節(jié)點的熱應變。

基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題)ansys結果圖片8

可見只有2號單元有熱應變。注意上述應變的符號與前面應變的符號相反,這是因為該單元相當于從100度降低到0度,是溫度降低要收縮,但是這種趨勢被阻止,所以產(chǎn)生了拉伸應變,因而與前面的壓縮應變符號相反。


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

相關標簽搜索:基于ANSYS APDL的兩端固定桿的單元生死仿真(適用鉆孔、隧道開挖、金屬切削等類似問題) Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術咨詢 ansys基礎知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓 

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

全國服務熱線

1358-032-9919

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