急！！版主求救：关于热源与对流的问题

ljian25

我做的高频焊接温度场模拟
这是焊管的焊接模拟
在它的横截面上加载线热源
在它的内表面加载空气对流
在它的外表面加载对流
可是总是出错
不知道是怎么回事
望各位高手帮忙看看
附上命令流：
/CLEAR,NOSTART
/FILNAM,2211
/TITLE,The Temperature Field Analysis of ERW Pipes  
/PREP7
/UNITS,SI
   
ET,1,SOLID70
et,2,surf152,,,,1,1
KEYOPT,2,4,1
KEYOPT,2,5,1
KEYOPT,2,8,4                             !施加对流，忽略热流密度                             
MPTEMP,1,20,250,500,750,1000,1500,1700
MPDATA,KXX,1,1,50,47,40,27,30,35,40
MPDATA,ENTH,1,1,0,5226.38E5,6271.41E5,7875E5,9320E6,9880E6,1.1E10
MPDATA,C,1,1,460,480,530,675,670,660,780
MPDATA,DENS,1,1,7820,7700,7610,7550,7490,7350,7300
MPDATA,HF,2,1,24,46,73,140,125,224,400,450,500
LENGTH=0.1
WIDTH=0.004
HEIGHT=6E-3
LSIZE=0.001
V=1 !速度
Z=1E-6
F=300000
K=0.5
CR=4E-6
A=9E-6
PX=1.2E- 8
H0= HEIGHT
LA=503/(K*SQRT(F))*(1e-3)
PRE=0.85
I=12
J0=I/(LA*H0)
Q=((PRE*(J0**2)*PX)/CR)*EXP((-2*Z)/LA)
TINC=LSIZE/V !单步时间长度
BLOCK,0, LENGTH,0, HEIGHT,0, WIDTH
ESIZE,0.001
MSHKEY,1
VATT,1,1,1
VMESH,1
/PNUM,DEFA
EPLOT
ASEL,S,LOC,Y,0.006  生成表面效应单元
type,2
ALLSEL
ASEL,S,LOC,Y,0.006  生成表面效应单元
type,2
VMESH,ALL
ALLSEL

FINISH
FINISH
                 
/SOLU
MAX_TIME=LENGTH/LSIZE+1
MAX_X=LENGTH/LSIZE+1
MAX_Y=HEIGHT/LSIZE+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)*LSIZE
*ENDDO
   
*DO,J,1,MAX_Y,1
FLUX2(0,J,K)=(J-1)*LSIZE
*ENDDO
   
*DO,I,1,MAX_X,1
*DO,J,1,MAX_Y,1
XCENTER=V*(K-1)*TINC
YCENTER=(J-1)*LSIZE
*IF,(I-1)*LSIZE,EQ,XCENTER,THEN
*IF,YCENTER,GE,0,THEN
*IF,YCENTER,LE,HEIGHT,THEN
FLUX2(I,J,K)=Q
*ELSE
FLUX2(I,J,K)=0
*ENDIF
*ENDIF
*ENDIF
*ENDDO
*ENDDO
FLUX2(0,0,K)=(K-1)*TINC
*ENDDO
ANTYPE,TRANS选择瞬态分析
TUNIF,20!初始温度20度
KBC,1 选择stepped（介跃式）分析方法
*SET,DELT,1
TIMINT,ON
AUTOTS,OFF
ALLSEL
Esel,s,type,,2
ASEL,S,LOC,y,0                              
sf,all,conv,5000,20
ALLSEL
Esel,s,type,,2
ASEL,S,LOC,y,0.006
sf,all,conv,5000,20                           
ALLSEL
   
TT= LENGTH/ V
*DO,TM,1E-6,TT,TINC
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
   
FLST,2,1,5,ORDE,1
FITEM,2,2
   
/GO
!*
TIME,TM
AUTOTS,0 !打开自动时间步长
DELTIM, , , ,1 !设置最小、最大计算时间步长
KBC,0 !选择rampped（递进式）分析方式
!*
TSRES,ERASE
!FLST,2,1,5,ORDE,1
!FITEM,2,2
/GO
ASEL,S,LOC,Z, 0.004, 0.004
SFA,ALL,1,HFLUX,%FLUX2%
/STATUS,SOLU
SOLVE
ALLSEL,ALL
OUTRES,ALL,ALL,
/PSF,HFLUX,,2     
/REPLOT
      
*ENDDO
OUTRES,ALL,ALL,
FINISH

ljian25

望各位高手和版主帮忙啊
小弟实在是不好解决了

ljian25

帮帮我啊
实在不会弄了

ljian25

兄弟们呢
救救急啊

ljian25

已经解决了
谢谢daji