ansys：用死活單元模擬澆鑄過程中的溫度分布

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

要想對澆鑄過程的溫度場進(jìn)行分析,必須熟悉下面兩個方面的內(nèi)容:1.瞬態(tài)溫度場的分析,2.單元死活的應(yīng)用。

瞬態(tài)溫度場分析:在進(jìn)行瞬態(tài)溫度分析時(shí),我們常遇到的一個問題是溫度結(jié)果明顯不合理,:計(jì)算得到的溫度高于給定的最高溫度或低于給定的最低溫度。造成這種結(jié)果的原因有兩個:
1、單元不合理,網(wǎng)格太大。
2、時(shí)間積分方法不恰當(dāng)。

針對上述兩個原因,我們采取的措施是:
1、細(xì)劃網(wǎng)格。
2、選擇向后積分法(THETA=1)。

單元死活:應(yīng)用單元死活進(jìn)行分析時(shí),需要用到多個載荷步。在不同的載荷步之間單元的死活狀態(tài)不同,ANSYS的通用動畫生成工具生成的動畫無法表現(xiàn)單元狀態(tài)隨時(shí)間的變化,因此必須使用ANSYS的DISPLAY Utility來生成能夠表現(xiàn)單元死活狀態(tài)的動畫。

下面是一個簡單的模擬澆鑄過程的例子,環(huán)境溫度為70,澆鑄溫度為600,計(jì)算結(jié)果應(yīng)是:最低溫度不低于70,最高溫度不高于600。如果我們劃分的網(wǎng)格不適當(dāng),則會出現(xiàn)最低溫度低于70或最高溫度高于600的情況,因此必須細(xì)劃網(wǎng)格并選擇THETA=1。計(jì)算完成后使用ANSYS的DISPLAY Utility來生成動畫文件。

fini
/cle
/prep7

et,1,55
et,2,55
KXX , 1,1.99/3600
DENS,1,0.00073
C ,1,0.100
KXX,2,7.38E-6
DENS,2,9.35E-5
c,2,1

k,1,-.5,-.5
k,4,6,-.5
k,12
k,13,5.5
k,21,-.5,6
k,22,,6
k,23,5.5,6
k,24,6,6
a,12,13,23,22
a,1,12,22,21
a,1,4,13,12
a,4,24,23,13
lsel,s,,,5,6
lsel,a,,,9,11,2
lesize,all,,,8
lsel,all
lesize,all,,,20
asel,s,,,1
aatt,2,2,2
asel,all
amesh,all
fini
/solu
antype,tran
!solc,on
timint,off
lsel,s,,,7,8
lsel,a,,,10
nsll,s
sf,all,conv,100,70
alls
esel,s,type,,2
nropt,full
ekill,all
alls
tunif,100
time,1
outres,all,all
solve

dtop=0
tim=1
delt=.05
timint,on
tintp,,,,1
*do,ii,1,20
ddele,all,temp
esel,s,type,,2
nsle,s
tim=tim+delt
time,tim
dtop=dtop+.3
nsel,r,loc,y,0,dtop
esln,s,1
ealive,all
nsel,r,loc,y,dtop
esln,s
nsle,s,1
d,all,temp,600
esel,s,live
nall
eall
solve
*enddo