一条焊缝的命令流程序，请高手指教！

飘叶和流星

大家好

我编写了一条焊缝的焊接命令流，但是只是加热阶段的程序，冷却阶段的老是处错误，希望有高手指教下，本人将不胜感激！补充一点，我用的是直接法耦合的，下面是我焊接加热阶段的命令流程序，运行没有错误。（注意记得glue一下）

/PreP7

mPtemP,l,20,250,500,750,1000,1500,1700,2500

mpdate,kxx,l,l,50,47,40,27,30,35,140,142

mpdata,dens,l,l,7820,7700,7610,7550,7490,7350,7300,7090

mpdata,c,1,l,460,480,530,675,670,660,780,820

MPDATA,PRXY,l,l,0.28,0.29,0.31,0.35,0.4,0.45,0.48,0.5

mPdata,alPx,l,l,1.10E-05,1.22E-05,1.39E-05,l.48E-05,l.34E-05,l.33E-05,1.32E-05,1.3lE-05

mPdata,ex,l,l,2.05E+11,1.87E+11,1.50E+11,7.00E+10,2.00E+10,1.90E+10,1.80E+10,l.20E+10

MPTEMP,,,,,,,,

MPTEMP,l,20

MPTEMP,2,1533

MPTEMP,3,1590

MPTEMP,4,1670

MPDATA,ENTH,l,,0

MPDATA,ENTH,l,,7.5E+009

MPDATA,ENTH,l,,9.6E+009

MPDATA,ENTH,l,,1.1E+010

MPDATA,MURX,1,,1

MPTEMP,,,,,,,,

MPTEMP,1,20

MPTEMP,2,250

MPTEMP,3,500

MPTEMP,4,750

MPTEMP,5,1000

MPTEMP,6,1500

MPTEMP,7,1700

MPTEMP,8,2500

MPDATA,HF,1,,1

MPDATA,HF,1,,3.5

MPDATA,HF,1,,5.2

MPDATA,HF,1,,10

MPDATA,HF,1,,15

MPDATA,HF,1,,30

MPDATA,HF,1,,31

MPDATA,HF,1,,35

!下面在各温度点下的屈服应力和屈服后的弹性模量

TB,BISO,l,6,2,

TBTEMP,20

TBDATA,,2.20E+08,2.05E+10,,,,

TBTEMP,500

TBDATA,,1.30E+08,l.50E+10,,,,

TBTEMP,750

TBDATA,,4.00E+07,7.00E+09,,,,

TBTEMP,1000

TBDATA,,2.50E+07,2.00E+09,,,,

TBTEMP,1700

TBDATA,,1.00E+06,1.80E+09,,,,

TBTEMP,2500

TBDATA,,1.00E+05,1.20E+09,,,,

SAVE

/prep7

/units,si

Rb=0.002

Length=0.12

width=0.03

height=0.003

et,l,solid5

et,2,solid98

block,0,length,0,width,0,height

block,0,length,0,3*Rb,0,height

block,0,length,0,6*Rb,0,height

block,0,length,0,12*Rb,0,height

VOVLAP,ALL

/view,1,1,1,1

!下面定义焊缝区的单元尺寸

Lesize,4,rb

Lesize,5,rb

Lesize,9,rb

Lesize,14,rb

Lesize,50,rb

Lesize,10,rb

Lesize,52,rb

Lesize,23,rb

Lesize,51,rb

Lesize,9,rb

Lesize,49,rb

Lesize,24,rb

!定义过渡区的单元尺寸

Lesize,31,2*rb

Lesize,26,2*rb

!定义远离焊缝区的单元尺寸

Lesize,43,4*rb

Lesize,38,4*rb

Lesize,2,4*rb

Lesize,7,4*rb

Lesize,53,4*rb

Lesize,54,4*rb

Lesize,55,4*rb

Lesize,56,4*rb

Lesize,11,1.5*rb

Lesize,12,1.5*rb

Lesize,47,1.5*rb

Lesize,48,1.5*rb

Save

记得glue一下

Type,1

Mshkey,1

Mshape,0

Vmesh,5

Vmesh,6

Type,2

Mshkey,0

Mshape,1

Vmesh,7

Vmesh,8

Save

/SOLU

DA,3,

/GO

DL,53,,ux,

DL,61,,UX,

DL,57,,UX,

DL,49,,UX,

DL,4,,uz,

DL,9,,uy,

/Solu

!时间步长设置

Antype,trans,new!

Nlgeom,on!

Timint,1,struct!

Timint,1,therm!

Timint,0,mag!

Timint,0,elect

TintP,0.005,,,-1,0.5,0.2

Tref, 25

Nropt,auto,,

Autos,on

Outres,,all

v=0.01

tinc=rb/v

!移动热源设置

U=25

I=180

ef =0.75

P=25*180*ef

r=0.006

qmax=p/(3.14159*r**2)

Max_time=length/rb+1

max_x=length/rb+1

max_y=width/rb+1

*dim,flux2,table,max_x,max_y,max_time,x,y,time

*do,k,1,max_time,1

*do,i,1,max_x,1

Flux2(i,0,k)=(i-1)*rb

*enddo

*do,j,1,max_y,1

Flux2(0,j,k)= (j-1)*rb

*enddo

*do,i,1,max_x,1

*do,j,1,max_y,1

!设置热源中心位置坐标

xcenter=v*(k-1)*tinc

ycenter=0

distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*rb-ycenter)**2)

*if,distance,le,r,then

Flux2(i,j,k)=qmax/EXP(3*DISTANCE**2/(R**2))

*else

Flux2(i,j,k)=0

*endif

*enddo

*enddo

Flux2(0,0,k)=(k-1)*tinc

Nsel,s,loc,z,height

Sf,all,hflux,%flux2%

Allsel,all

Nsel,s,loc,x,0

Sf,all,conv,10,20

Allsel,all

Nsel,s,loc,x,length

Sf,all,conv,10,20

Allsel,all

Nsel,s,loc,y,0

Sf,all,conv,10,20

Allsel,all

Nsel,s,loc,y,width

Sf,all,conv,10,20

Allsel,all

Nsel,s,loc,z,0

Sf,all,conv,10,20

Allsel,all

Time,k*rb/v

Nsubst,5

solve

*enddo

save

下面冷却阶段的命令流哪位高手能给指教一下，谢谢了！