再次求助！管子对焊问题

liangzaisimwe

再次求助！管子对焊问题我在论坛上搜到了一个apdl的例子，但是那个例子也是没有最终做出来，我改了改，结果还是不行
焊缝热场也是轴向延伸的，没有体现出沿着圆周移动的效果，请高手指点啊！
下面是我参照一些例子改过的apdl（管子对接焊接），哪位高手帮着看看是哪的问题啊，万分感谢！
/prep7  
/units,si   
*afun,deg   
*set,inradius,0.02  
*set,outradius,0.03
*set,degree,90  
*set,length,0.1
et,1,solid70
mp,dens,1,2640  
mptemp,1,100,200,400,600,800,1000   
mptemp,7,1200   
mpdata,c,1,1,660.4,836,973.9,1061.7,1086.8,1103.5   
mpdata,c,1,7,1128.6
mpdata,enth,1,1,0,7.583e4,2.592e5,4.619e5,6.875e5,9.528e5   
mpdata,enth,1,7,1.278e6
mptemp,1,20,100,200,400,600,800
mptemp,7,1000,1200  
mpdata,kxx,1,1,1.38,1.53,1.65,1.86,2.02,2.15
mpdata,kxx,1,7,2.29,2.43
csys,1  
k,1,inradius,0,0
k,2,inradius,degree,0   
k,3,inradius,degree,length  
k,4,inradius,0,length   
k,5,outradius,0,0   
k,6,outradius,degree,0  
k,7,outradius,degree,length
k,8,outradius,0,length  
v,1,2,3,4,5,6,7,8   
/view,1,1,1,1   
lesize,5,,,5   
lesize,6,,,25   
lesize,8,,,50   
mshkey,1
vatt,1,1,1  
vmesh,1
save

*set,delt,1
*SET,qmax,10000000  
*SET,r,0.004
*SET,pi,3.14159
*SET,lsize,0.002
*SET,v,0.003  
*SET,tinc,lsize/v   
*dim,flux2,table,25,50,25,y,z,time  
*do,k,1,25,1
*do,i,1,25,1
*SET,flux2(i,0,k),(i-1)*lsize/(0.03*pi)*180   
*enddo  
*do,j,1,50,1
*SET,flux2(0,j,k),(j-1)*0.002   
*enddo  
*do,i,1,25,1
  *do,j,1,50,1
  *afun,deg
  *SET,ycenter,v*(k-1)*tinc/(0.03*pi)*180
  *SET,zcenter,length/2
    a=(i-1)*lsize-ycenter
    b=(j-1)*0.002-zcenter
   *SET,distance,sqrt(a**2+b**2)
   *if,distance,le,r,then
   *SET,flux2(i,j,k),qmax/EXP(DISTANCE**2/(r**2))  
   *else   
   *SET,flux2(i,j,k),0
   *endif  
  *enddo  
*enddo   
*SET,flux2(0,0,k),(k-1)*tinc
*enddo
/solu
*afun,deg
antype,trans
nroPt,auto,,!在完全瞬态分析中,指定Newton-raphson选项
*do,t,1e-6,15.71,tinc
time,t
autos,0
deltim,,,,1
kbc,1
allsel,all
asel,s,,,6
nsel,r,loc,y,0,90
nsel,r,loc,z,0,0.1
sf,all,hflux,%flux2%

timint,on                                        ！打开时间积分
tintp,,,,1                                        ！将THETA设定为1
outres,all,all                                        ！输入每个子步的结果
solve
*enddo
finish

liangzaisimwe

没有人理我啊？
大家来帮帮我啊