为什么只能加载一面热流密度，另一面加载不上？？？用的高斯热源，高手支招！！！

19850913

命令流如下
/clear
/units,si  
/config,nres,10000
/filname,welding,1

/prep7
et,1,solid90
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1700
MPTEMP,8,2750
MPDATA,DENS,1,,7800
MPDATA,DENS,1,,7700
MPDATA,DENS,1,,7610
MPDATA,DENS,1,,7550
MPDATA,DENS,1,,7490
MPDATA,DENS,1,,7350
MPDATA,DENS,1,,7300
MPDATA,DENS,1,,7090
MPTEMP,,,,,,,,,
MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1600
MPTEMP,8,1700
MPTEMP,9,2500
MPDATA,KXX,1,,50
MPDATA,KXX,1,,47  
MPDATA,KXX,1,,40
MPDATA,KXX,1,,27
MPDATA,KXX,1,,30
MPDATA,KXX,1,,35
MPDATA,KXX,1,,200               
MPDATA,KXX,1,,500               
MPDATA,KXX,1,,1500
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,250
MPTEMP,3,500
MPTEMP,4,750
MPTEMP,5,1000
MPTEMP,6,1500
MPTEMP,7,1700
MPTEMP,8,2750  
MPDATA,C,1,,460
MPDATA,C,1,,480
MPDATA,C,1,,530
MPDATA,C,1,,675
MPDATA,C,1,,670
MPDATA,C,1,,660
MPDATA,C,1,,2000
MPDATA,C,1,,2500
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,1200   
MPTEMP,9,2500   
MPDATA,EX,1,,208000
MPDATA,EX,1,,204000
MPDATA,EX,1,,199000
MPDATA,EX,1,,191000
MPDATA,EX,1,,183000
MPDATA,EX,1,,174000
MPDATA,EX,1,,164000
MPDATA,EX,1,,80000  
MPDATA,EX,1,,80000  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,1500   
MPTEMP,9,2500   
UIMP,1,REFT,,,25  
MPDATA,ALPX,1,,1.37e-5  
MPDATA,ALPX,1,,1.37e-5  
MPDATA,ALPX,1,,1.38e-5  
MPDATA,ALPX,1,,1.41e-5  
MPDATA,ALPX,1,,1.44e-5  
MPDATA,ALPX,1,,1.48e-5  
MPDATA,ALPX,1,,1.5e-5   
MPDATA,ALPX,1,,1.5e-5   
MPDATA,ALPX,1,,1.5e-5   
TB,BISO,1,6,2,  
TBTEMP,20   
TBDATA,,370,732,,,,
TBTEMP,250  
TBDATA,,289,782,,,,
TBTEMP,500  
TBDATA,,230,315,,,,
TBTEMP,650  
TBDATA,,30,50,,,,   
TBTEMP,1200
TBDATA,,3,10,,,,
TBTEMP,2500
TBDATA,,3,10,,,,
TBDE,BISO,1,,,  
TB,BISO,1,6,2,  
TBTEMP,20   
TBDATA,,370,732,,,,
TBTEMP,250  
TBDATA,,289,782,,,,
TBTEMP,500  
TBDATA,,230,315,,,,
TBTEMP,650  
TBDATA,,30,50,,,,   
TBTEMP,1200
TBDATA,,3,10,,,,
TBTEMP,2500
TBDATA,,3,10,,,,

k,1,0,0,0
k,2,0,0,0.12
k,3,0,0.006,0.12
k,4,0,0.006,0
k,5,0,0,0.14
k,6,0,0.006,0.14
k,7,0,0.006,0.146
k,8,0,0,0.146
k,9,0,0,0.154
k,10,0,0.006,0.154
k,11,0,0.006,0.16
k,12,0,0,0.16
k,13,0,0,0.18
k,14,0,0.006,0.18
k,15,0,0.006,0.3
k,16,0,0,0.3
k,17,0,0.012,0.154
k,18,0,0.012,0.146
k,19,0,0.032,0.146
k,20,0,0.032,0.154
k,21,0,0.156,0.154
k,22,0,0.156,0.146
a,1,2,3,4
a,2,5,6,3
a,5,6,7,8
a,7,8,9,10
a,9,10,11,12
a,11,12,13,14
a,13,14,15,16
a,6,7,18
a,7,10,17,18
a,10,11,17
a,17,18,19,20
a,19,20,21,22
aglue,all

/view,1,1,1,1
/replot
vext,all,,,0.3,0,0
vglue,all
lesize,70,,,50,,,,,1
lesize,67,,,20,,,,,1
mat,1
type,1
vmesh,7
lesize,40,,,50,,,,,1
lesize,36,,,20,,,,,1
mat,1
type,1
vmesh,1

lesize,64,,,5,,,,,1
mat,1
type,1
vmesh,6
lesize,44,,,5,,,,,1
mat,1
type,1
vmesh,2
lesize,57,,,2,,,,,1
mat,1
type,1
vmesh,5
lesize,54,,,2,,,,,1
mat,1
type,1
vmesh,4
lesize,47,,,2,,,,,1
mat,1
type,1
vmesh,3

vmesh,9
lesize,81,,,5,,,,,1
mat,1
type,1
vmesh,11
lesize,84,,,31,,,,,1
mat,1
type,1
vmesh,12
lesize,73,,,2,,,,,1
vmesh,8
lesize,78,,,2,,,,,1
vmesh,10
nummrg,all
numcmp,all  
/image,save,vmesh,jpeg  
save

rb=0.006
zb=rb/6
length=0.3
width=0.005
/solu
antype,4  
tunif,25
TRNOPT,FULL
OUTRES,NSOL,all
NROPT,FULL,,ON
AUTOTS,1
KBC,1
lnsrch,1     
TSRES,ERASE
TINTP,,,,1,,,,,
TIMINT,0,struct
TIMINT,1,THERM
TIMINT,0,MAG

allsel,all
nsel,s,ext
sf,all,conv,15,25   
v=0.006
tinc=rb/v
tm=length/v
U=15
I=100
ef=0.7
p=U*I*ef
r=0.005
qmax=3*p/(3.14159*r**2)
max_time=length/rb+1
max_x=length/rb+1
max_y=width/zb+1
local,11,0,0,0.009,0.143,,45,
*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)*zb
    *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)*zb-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  
*enddo
t=0
*do,i,1,5,1
        /solu
        t=t+1
time,t
        NSUBST,2,,,1
        allsel,all
asel,s,area,,44
        nsla,s,1
sfdele,all,conv
sf,all,hflux,%flux2%
allsel,all
solve
asel,s,area,,44
        nsla,s,1
sfdele,all,hflux
sf,all,conv,15,25
allsel,all
*enddo

*DO,I,1,20
        /solu
allsel,all
        t=t+1
TIME,t
        NSUBST,2,,,1
SOLVE
*ENDDO

*DO,I,1,20  
       /solu
       allsel,all
       t=t+1*10
       TIME,t
       NSUBST,2,,,1
       SOLVE
*ENDDO
save
/solu
antype,4  
tunif,25
TRNOPT,FULL
OUTRES,NSOL,all
NROPT,FULL,,ON
AUTOTS,1
KBC,1
lnsrch,1     
TSRES,ERASE
TINTP,,,,1,,,,,
TIMINT,0,struct
TIMINT,1,THERM
TIMINT,0,MAG
rb=0.006
zb=rb/6
length=0.3
width=0.005
v=0.006
tinc=rb/v
tm=length/v
U=15
I=100
ef=0.7
p=U*I*ef
r=0.005
qmax=3*p/(3.14159*r**2)
max_time=length/rb+1
max_x=length/rb+1
max_y=width/zb+1
csdele,11,,1
local,11,0,0,0.009,0.157,,135,
*dim,flux3,table,max_x,max_y,max_time,x,y,time
*do,k,1,max_time,1
    *do,i,1,max_x,1
      flux3(i,0,k)=(i-1)*rb
    *enddo
*do,j,1,max_y,1
      flux3(0,j,k)=(j-1)*zb
    *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)*zb-ycenter)**2)
     *if,distance,le,r,then
      flux3(i,j,k)=qmax/exp(3*distance**2/(r**2))
     *else
      flux3(i,j,k)=0
     *endif
    *enddo
   *enddo
  flux3(0,0,k)=(k-1)*tinc  
*enddo

*do,i,1,5,1
        /solu
        t=t+1
time,t
        NSUBST,2,,,1
        allsel,all
asel,s,area,,49
        nsla,s,1
sfdele,all,conv
sf,all,hflux,%flux2%
allsel,all
solve
asel,s,area,,49
        nsla,s,1
sfdele,all,hflux
sf,all,conv,15,25
allsel,all
*enddo

其中，44面可以加载上，我为了调程序就只取到5秒，另一面总是加载在末端的一个点，摸索了很久，不知道为什么？请高手指点


我没能放图片在上面