求助：生死单元后，再次求解不合理
本帖最后由 xiebao143 于 2011-11-3 09:15 编辑激光去除加工过程仿真，在第一次加载后，设定超过熔化部分的全去除（杀死或者改变材料属性），自然冷却，在加载，求解，去除，冷却。重启动重复这个过程，但是在冷却时结果不合理，也就是生死单元后，就出错了，无法冷却求解，期待您的帮助。程序如下：
/filename,killelment/title, the simulate of laser cladding process
/prep7
BLOCK,0,100e-6,0,-50e-6,-50e-6,50e-6, ET,1,SOLID70ET,2,SURF152keyopt,2,4,1keyopt,2,5,0keyopt,2,8,1!hfluxtoffst,273 mp,DENS,1,7840mp,kxx,1,34mp,c,1,500
Mshkey,1
Mshape,0,3d
type,1 Mat,1
esize,5e-6,
Vmesh,all
Eplot
asel,all
type,2 Mat,1
amesh,all
esel,s,type,,1 Asel,s,,,1,6Nsla,s,1
sf,all,conv,800,25
allsel,all
/psf,hflux,,2,0
/SOLU
antype,transient,new
TRNOPT,FULL
Kbc,1
tunif,25neqit,25 !nlgeom,on
autots,on
timint,on,therm
TINTP,,,,1
Lnsrch,on nropt,full
outres,all,all
Allsel
nsel,all Aaa=50e-6
local,11,1, Aaa,0,0,0,-90,0
CSYS,11
esel,s,type,,2
nsel,r,loc,x, 0,10e-6
nsel,r,loc,y,0,180
nsel,r,loc,z,0
sf,,hflux,7e10t=0
dt=5e-6
csys,0
t=t+dt
time,t
deltim,1e-6,0.5e-6,1e-6
Allsel
Solve saveFinish/post1Set,lastetable,melty22,temp,ESEL,S,ETAB,MELTY22,400finish/solu*get,elemnum,elem,0,count
*if,elemnum,gt,0,then
ekill,allesel,s,live nsle,s,1 nsel,invert d,all,all,0
*elseallsel
*endifEsel,s,live
!在这之前都没问题，求解后结果就出错了
sfdele,all,hflux！删除载荷，冷却
t=t+2*dt
time,tAllsel
Solvefinish!后面的就更无法运算/solu resume antype,,rest！重复上述过程
页:
[1]