用单元死活模拟浇铸过程中的温度分布 要想对浇铸过程的温度场进行分析,必须熟悉下面两个方面的内容:1 瞬态温度场的分析,2.单元死活的应用 瞬态温度场分析:在进行瞬态温度分析时,我们常遇到的一个问题是温度结果明 显不合理,:计算得到的温度高于给定的最高温度或低于给定的最低温度。造成 这种结果的原因有两个:1、单元不合理,网格太大。2、时间积分方法不恰当 针对上述两个原因,我们采取的措施是:1、细划网格。2、选择向后积分法 (THETA=1) 单元死活:应用单元死活进行分析时,需要用到多个载荷步。在不同的载荷步之 间单元的死活状态不同, ANSYS的通用动画生成工具生成的动画无法表现单元状 态随时间的变化,因此必须使用 ANSYS的 DISPLAY Utility来生成能够表现单元 死活状态的动画 下面是一个简单的模拟浇铸过程的例子,环境温度为70,浇铸温度为600,计算 结果应是:最低温度不低于70,最高温度不高于600。如果我们划分的网格不适 当,则会出现最低温度低于70或最高温度高于600的情况,因此必须细划网格 并选择 THETA=1。计算完成后使用 ANSYS的 DISPLAY Utility来生成动画文件。 /cle /prep t,1,55 KXX,1,1.99/3600 DENS,1,0.00073 C,1,0.100 KXX,2,7.38E-6 DENS,2,9.35E-5 k,1,-.5,-.5 k,4,6,-.5 k,13,5.5 k,21,-.5,6 k,22,,6 k,23,5.5,6 k,24,6,6 12,13,23,22 a,1,12,22,21 a,1,4,13,12 a,4,24,23,13
用单元死活模拟浇铸过程中的温度分布 要想对浇铸过程的温度场进行分析,必须熟悉下面两个方面的内容：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 /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 les l1,,,8 aaaa t,2,2,2 fi imint, off lln el a f, all, conv. 100.70 kill. all all so⊥ve dtop=o delt=0. imint. on tint adele,all, te timstim+delt nsel, r, loc, y, 0, dtop esIn. s. 1 alive. all nsel, r, loc, y, dtop
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
d, all, temp, 600 esel,s. live solve kendo fi /show, tmpplt, f3 dtop=O 2,21 dtop=dtop. 3 nsel, r, loc, y, 0, dtop esIn. s, 1 pIns, temp enddo
esln,s nsle,s,1 d,all,temp,600 esel,s,live nall eall solve *enddo fini /show,tmpplt,f33 /post1 dtop=0 *do,ii,2,21 set,ii dtop=dtop+.3 nsel,r,loc,y,0,dtop esln,s,1 esel,a,type,,1 nall plns,temp *enddo