各位ansys同仁:
小弟在用ansys做簡單生死單元技術(shù)時,遇到了一些問題,還望各位幫忙解答一下。
模型的APDL程序如下:
finish
/clear
/prep7
k,1,0,0
k,2,0.03,0
k,3,0.03,0.03
k,4,0,0.03
a,1,2,3,4
nt=20
tall=60
tsize=tall/nt
ET,1,Plane13,2
MP,Dens,1,920
mp,c,1,2040
mp,kxx,1,1.8
esize,0.0005
amesh,all
finish
/solu
tunif,300
nsel,s,ext
sf,all,conv,10,2000
antype,4
nropt,full
lumpm,0
neqit,100
autots,on
outres,all,all
time,10*tsize
nsel,all
solve
finish
*do,i,2,4
/solu
resume
antype,,rest
time,i*tsize*10
nsel,s,ext
sfdele,all,conv
nropt,full
allsel
solve
save
finish
/post1
set,last
etable,melty,temp
esel,s,etab,melty,340
finish
/solu
ekill,all
esel,s,live
nsle,s
nsel,inve
D,all,all,0
finish
*enddo
不過程序的結(jié)果有問題,這個程序參考的是University of Alberta一個簡單程序(如上或者附件所示),幾何邊界條件為對流。我單獨用了一個tsize步長求解,為了循環(huán)語句中的resume和antype,,rest。同時為了研究生死單元技術(shù),我在循環(huán)中就將載荷刪除施加的載荷,為的是檢測是否生死單元技術(shù)能正確解決問題。
現(xiàn)在主要問題如下:
1. 得到的時間段有問題,理論上求解時間應該為30(最開始的單獨語句-非*do循環(huán)內(nèi))、60、90和120[s],但實際得到的是30、31、32和33[s]。
2. 在循環(huán)語句中,我已經(jīng)將載荷刪除,且對需要殺死單元的節(jié)點進行約束,按照理論上說被殺死的單元或者節(jié)點在后續(xù)求解時間內(nèi)就不應該有溫度變化才對?(是不是我對有限元或者生死單元技術(shù)理解有問題?)。部分圖片如附件所示。
圖片:31秒的所有單元和僅僅為活單元的溫度場
![]()
![]()
32秒的所有單元和僅僅為活單元的溫度場
![]()
![]()
33秒的所有單元和僅僅為活單元的溫度場
![]()
![]()
不知道有人遇到相同問題沒有?又是如何解決的呢?謝謝您的指點。 |