ansys点焊模拟求解释

/prep7
rect,0/1000,3/1000,0/1000,0.01/1000
rect,3/1000,10/1000,0/1000,0.01/1000
rect,0/1000,3/1000,0.01/1000,0.02/1000
rect,20/1000,3/1000,0.01/1000,0.02/1000
rect,0/1000,3/1000,1.52/1000,0.02/1000
rect,20/1000,3/1000,1.52/1000,0.02/1000
rect,0/1000,3/1000,1.52/1000,1.53/1000
rect,20/1000,3/1000,1.52/1000,1.53/1000
rect,6/1000,8/1000,8.94/1000,22.94/1000
k,100,0/1000,8.94/1000
k,101,6.9/1000,4.94/1000
aglue,all
lstr,28,100
lstr,100,33
larc,27,34,101
al,5,7,8,33,27
aglue,all
save
f=3500
r_znzn=1.355e-5
r_zncu=4.967e-6
et,1,plane67
keyopt,1,3,1
mp,rsvx,3,2.64e-8
mp,kxx,3,390
mp,c,3,39.8
mp,dens,3,8900
mptemp,1,20,200,400,600,800,1300
mpdata,rsvx,1,1,0.142e-6,0.262e-6,0.458e-6,0.734e-6,1.081e-6,1.220e-6
mptemp,1,50,150,300,650,700,750
mpdata,c,1,1,481,523,569,858,1139,960
mptemp,1,0,200,600,800,1000,1200
mpdata,kxx,1,1,65.3,54.9,45.2,36.4,28.5,27.6,29.7
mp,dens,1,7800
mp,c,2,390
mp,dens,2,7133
mp,kxx,2,270
mp,rsvx,2,4.7e-8
mp,c,4,390
mp,dens,4,7133
mp,kxx,4,270
mptemp,1,21,204,426
mpdata,rsvx,4,1,r_znzn+4.7e-8,r_znzn*0.95+4.7e-8,r_znzn*0.9+4.7e-8
mp,c,5,390
mp,dens,5,7133
mptemp,1,21,204,426
mp,kxx,5,270
mptemp,1,21,204,426
mpdata,rsvx,5,1,r_zncu+4.7e-8,r_zncu*0.95+4.7e-8,r_zncu*0.9+4.7e-8
asel,s,,,2,9,7
aatt,3
alls
asel,s,,,13
aatt,5
alls
asel,s,,,15
aatt,2
alls
asel,s,,,12,14,2
aatt,1
alls
asel,s,,,1
asel,a,,,10,11
asel,a,,,16
aatt,4
alls
lsel,s,,,2,6,2
lsel,a,,,14
lsel,a,,,39,40
lsel,a,,,43,44
lsel,a,,,48
lesize,all,,,1
alls
lsel,s,,,33,35,2
lesize,all,,,3
alls
lsel,s,,,34,36,2
lesize,all,,,15
alls
lsel,s,,,41,42
lsel,a,,,46
lesize,all,,,15
alls
lesize,5,,,20,10
lesize,8,,,40,10
lsel,s,,,1,3,2
lsel,a,,,7,11,4
lsel,a,,,19,27,8
lesize,all,,,
alls
lsel,s,,,47,49,2
lesize,all,
alls
lsel,s,,,37,38
lsel,a,,,45,50,5
lesize,all,
alls
amesh,all
lsel,s,,,35
nsll,s,1
cp,2,volt,all
/solu
lsel,s,,,1
dl,all,blank,volt,0
alls
lsel,s,,,7,36,29
sfl,all,conv,3800
alls
lsel,s,,,8,34,26
lsel,a,,,49
sfl,all,conv,25
alls
fk,35,amps,10000
antype,4
kbc,1
autots,on
outress,all,all
nlgeom,1
sstif,1
nropt,full, ,on
eqslv, , ,0,
precision,0
msave,0
toffst,0,
physics,write,thermal
physics,clear
finish
/prep7
etchg,tts
mp,ex,1,0.207e9
mp,prxy,1,0.275
mp,dens,1,7800
mp,ex,2,0.1e9
mp,prxy,2,0.27
mp,dens,2,7133
mptemp,3,20,100,150,250,300,400,500,600,700
mpdata,ex,3,1,0.115e9,0.11e12,0.106e12,0.102e12,0.101e12,0.097e12,0.09e12,0.081e12,0.075e12
mptemp,3,20,100,150,250,300,400,500,600,700
mpdata,prxy,3,1,0.35,0.35,0.35,0.36,0.37,0.37,0.38,0.38,0.39
mp,dens,3,8900
mp,ex,4,0.1e12
mp,prxy,4,0.27
mp,dens,4,7133
mp,ex,5,0.115e12
mp,prxy,5,0.35
mp,dens,5,7800
lsel,s,,,35
nsll,s,1
cp,3,uy,all
alls
lsel,s,,,7
nsll,s,1
cp,4,uy,all
alls
/solu
lsel,s,,,4,5
lsel,a,,,40,44,2
lsel,a,,,36
dl,all,,ux
alls
lsel,s,,,1
lsel,a,,,37
dl,all,,uy
alls
fk,35,fy,-3500
antype,4
kbc,1
autots,on
outress,all,all
physics,write,struct
/solu
physics,read,thermal
ic,all,temp,25
time,0.02
deltim,0.005
solve
*do,a,1,19
*get,elemno,node,,count
*dim,elemtemp,,elemno
*do,aa,1,elemno
*get,elemtemp(aa),node,aa,temp
*enddo
physics,read,struct
/solu
ldread,temp,,, , ,' ','rst',' '
time,0.02
deltim,0.001
solve
/solu
*dim,kill,,50
mm=1
*do,i,1,15
*do,j,1,4
*get,nodeno,elem,i,node,j
*get,nodesy,node,nodeno,s,y
*if,nodesy,gt,0,then
kill(mm)=i
*endif
*enddo
*if,kill(mm),ne,0,then
mm=mm+1
*endif
*enddo
physics,read,thermal
*do,aa,1,elemno
ic,aa,temp,elemtemp(aa)
*enddo
/solu
*do,ii,1,15
*if,kill(ii),ne,0,then
ekill,kill(ii)
*endif
*enddo
time,0.02
deltim,0.005
solve
*do,jj,1,15
*if,kill(jj),ne,0,then
ealive,kill(jj)
*endif
*enddo
*enddo
后面的循环的意义如下：
1）假定钢板之间的接触面积等于钢板的表面积，计算热-电耦合，得出温度场
的分布。在一个小的时间步里面完成这一过程。
2）在同样的时间步里，计算热-应力分析结果，上一步得到的温度场分布作为
该步计算的初始条件。为了在两个顺序耦合的分析中方便的调用数据，两个顺序耦
合场采用同一个模型。在将上一步得到的温度场作为初始条件输入时，每个节点对
应的温度分别输入。
3）这一步主要目的为改变接触面积。在建模的时候，添加的过渡层单元用来模
拟钢板之间的接触状态。热-应力分析之后，确定了每个节点的应力分布。当应力为
压应力时，对应的部分认为处于接触状态。反之，若为拉应力，对应的部分认为处
于非接触状态，而对应的单元即利用 ANSYS 里的单元生死功能将其杀死。所谓的
杀死即为不参与下一步的计算。在之后的每一次循环中，在同样的时间步里，均采
用单元生死技术完成单元的杀死或激活。
4）将热-应力分析转回热-电分析，回到第1）步。然而，在这一步，接触面积
不再等于钢板的表面积，而是通过第3）步确定的。这一步也是在和1）和2）里同
样的时间步完成的。
5）激活在第3）步里被杀死的单元，回到热-应力分析，将第4）步计算得到的
温度场分布作为初始条件输入，在同样的时间步里完成新的热-应力分析。
6）完成循环，直到时间结束或者达到计算所需的收敛极限。
哪位大神帮我解释一下循环每一步的意思从do，a，1,19开始