用单元死活模拟浇铸过程中的温度分布
??? 要想对浇铸过程的温度场进行分析,必须熟悉下面两个方面的内容:1.瞬态温度场的分析,2.单元死活的应用。瞬态温度场分析:在进行瞬态温度分析时,我们常遇到的一个问题是温度结果明显不合理,:计算得到的温度高于给定的最高温度或低于给定的最低温度。造成这种结果的原因有两个:1、单元不合理,网格太大。2、时间积分方法不恰当。针对上述两个原因,我们采取的措施是:1、细划网格。2、选择向后积分法(THETA=1)。单元死活:应用单元死活进行分析时,需要用到多个载荷步。在不同的载荷步之间单元的死活状态不同,ANSYS的通用动画生成工具生成的动画无法表现单元状态随时间的变化,因此必须使用ANSYS的DISPLAY Utility来生成能够表现单元死活状态的动画。下面是一个简单的模拟浇铸过程的例子,环境温度为70,浇铸温度为600,计算结果应是:最低温度不低于70,最高温度不高于600。如果我们划分的网格不适当,则会出现最低温度低于70或最高温度高于600的情况,因此必须细划网格并选择THETA=1。计算完成后使用ANSYS的DISPLAY Utility来生成动画文件。fini/cle/prep7et,1,55?et,2,55?KXX , 1,1.99/3600DENS,1,0.00073?C ,1,0.100KXX,2,7.38E-6?DENS,2,9.35E-5?c,2,1?k,1,-.5,-.5?k,4,6,-.5?k,12k,13,5.5k,21,-.5,6?k,22,,6?k,23,5.5,6?k,24,6,6a,12,13,23,22?a,1,12,22,21a,1,4,13,12?a,4,24,23,13lsel,s,,,5,6lsel,a,,,9,11,2lesize,all,,,8?lsel,alllesize,all,,,20?asel,s,,,1?aatt,2,2,2?asel,allamesh,all?fini/solu?antype,tran!solc,ontimint,offlsel,s,,,7,8lsel,a,,,10nsll,s?sf,all,conv,100,70allsesel,s,type,,2?nropt,full?ekill,all?allstunif,100time,1outres,all,allsolvedtop=0tim=1delt=.05timint,ontintp,,,,1*do,ii,1,20ddele,all,tempesel,s,type,,2nsle,s?tim=tim+delttime,timdtop=dtop+.3nsel,r,loc,y,0,dtop?esln,s,1ealive,all?nsel,r,loc,y,dtopesln,snsle,s,1d,all,temp,600esel,s,livenalleallsolve*enddofini/show,tmpplt,f33/post1dtop=0*do,ii,2,21set,iidtop=dtop+.3nsel,r,loc,y,0,dtopesln,s,1esel,a,type,,1nall?plns,temp*enddo