绝对经典－载荷覆盖问题－对流和热流密度

mashilong

毕业设计中的热分析程序：
程序介绍：在平板上（面２）加高斯热源（取１／４平板），平板在煤油中加工，所以想加上对流，所以在面２，４，６上施加了对流载荷（表面上生成了surf152单元，经检查确是生成了），结算后表明高斯热源加上了，但是对流载荷总加不上（运算中出了载荷覆盖的警告）。请您帮我看看程序加载上出了什么问题。最好详细点！谢谢
finish
/clear
/title,transient analysis(gauss heat flux)
/filnam,transientgauss
/TRIAD,RBOT
/UNITS,SI
/PREP7
et,1,solid70
et,2,surf152
keyopt,2,8,2　　　!要在surf152上施加对流
X1=0.2　　　　　　　 !输入平板的能量比例
U=25　　　　　　　!电压
I=3　　　　　　　　　!电流
Ton=100      　　　!脉冲宽度
pi=3.1415926                     
Q=3*X1*U*I                  
lisannum=4           ! 控制载荷布的数量   
dividenum=20          !控制网格划分的密度
block_x=0.0005
block_y=0.0005
block_z=-0.00005
BLOCK,0,block_x,0,block_y,block_z,0,
/VIEW,1,1,1,1
save
mptemp,1,20,250,500,750,1000,1500,1700,2500 　　　　!定义材料属性
mpdata,c,1,1,460,480,530,675,670,660,780,820
mpdata,kxx,1,1,50,47,40,27,30,35,140,142
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090
mpdata,hf,1,1,100,350,520,1000,1500,3000,3100,3500
SAVE
FLST,5,2,4,ORDE,2   　　　　!分网部分
FITEM,5,4   
FITEM,5,-5  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,dividenum,5, , , ,1
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,1   
FITEM,5,8   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,dividenum,5, , , ,1
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,1   
FITEM,5,8   
CM,_Y1,LINE
LSEL, , , ,P51X
*GET,_z1,LINE,,COUNT
*SET,_z2,0  
*DO,_z5,1,_z1   
*SET,_z2,LSNEXT(_z2)
*GET,_z3,LINE,_z2,ATTR,NDNX
*GET,_z4,LINE,_z2,ATTR,SPNX
*get,_z6,line,_z2,attr,kynd
*IF,_z3,GT,0,THEN   
*IF,_z4,NE,0,THEN   
LESIZE,_z2,,,_z3,1/_z4,,,,_z6   
*ENDIF  
*ENDIF  
*ENDDO  
CMSEL,S,_Y1
CMDELE,_Y1  
!*  
FLST,5,3,4,ORDE,3   
FITEM,5,9   
FITEM,5,-10
FITEM,5,12  
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1, , ,dividenum/5,5, , , ,1  
!*  
FLST,5,1,4,ORDE,1   
FITEM,5,10  
CM,_Y1,LINE
LSEL, , , ,P51X
*GET,_z1,LINE,,COUNT
*SET,_z2,0  
*DO,_z5,1,_z1   
*SET,_z2,LSNEXT(_z2)
*GET,_z3,LINE,_z2,ATTR,NDNX
*GET,_z4,LINE,_z2,ATTR,SPNX
*get,_z6,line,_z2,attr,kynd
*IF,_z3,GT,0,THEN   
*IF,_z4,NE,0,THEN   
LESIZE,_z2,,,_z3,1/_z4,,,,_z6   
*ENDIF  
*ENDIF  
*ENDDO  
CMSEL,S,_Y1
CMDELE,_Y1  
!*  
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,';VOLU';   
CMSEL,S,_Y  
MSHKEY,1
VATT,1,1,1
VMESH,1
/PNUM,DEFA
EPLOT
FINISH　　　　　　　　　　　　　!分网结束
*dim,R,array,lisannum,1,1　　　!高斯热流密度的计算
tempr=0.11*I*I+0.73*I+8.52
*IF,tempr,GE,Ton,then               
    maxR=2.85*I**0.53*Ton**0.38           
*ELSE
    maxR=2.85*I**0.53*tempr**0.38                  
*ENDIF
maxenum=dividenum+1
*dim,flux2,table,maxEnum,maxEnum,lisannum,x,y,time   
*DO,k,1,lisannum,1
  t=k*Ton/lisannum      
    R(k)=2.85*I**0.53*t**0.38  
      *IF,R(k),GE,maxR,then                    
       R(k)=maxR
      *ENDIF
  *do,j,1,maxEnum,1
   flux2(0,j,k)=(j-1)*block_x/dividenum
  *enddo   
  *do,i,1,maxEnum,1
   flux2(i,0,k)=(i-1)*block_x/dividenum
  *enddo
  *do,i,1,maxEnum,1
    *do,j,1,maxEnum,1
     xcenter=0
     ycenter=0
     distance=((i-1)*block_x/dividenum-xcenter)**2+((j-1)*block_x/dividenum-ycenter)**2
    *if,distance,le,(R(k)*1E-6)**2,then
        flux2(i,j,k)=(Q/(pi*R(k)*R(k)*1E-12)*exp(-3*distance/R(k)*R(k)*1E-12)) !????????
    *else
        flux2(i,j,k)=0
    *endif   
    *enddo
  *enddo
       flux2(0,0,k)=(k-1)*Ton/lisannum
*enddo
/SOLU   
ANTYPE,trans
TRNOPT,FULL
NROPT,AUTO, ,   
EQSLV,   　　　　　　
SOLCONTROL,ON
AUTOTS,on
KBC,1                           
OUTRES,ALL,ALL,                           
TIME,1e-6
DELTIM,1E-6,1E-6,1E-6   
TUNIF,25      
kbc,1
solve                !先施加25度的开始温度计算，只是第一步起作用
type,2
asel,s,,,2,6,2
NSLA,S,1
ESLN,S,0
ESURF 　　　　　　　!在２，４，６面上生成表面效应单元
esel,s,type,,2
sf,all,conv,-1,25 　!施加对流（hf已经在前面材料属性部分定义了）
allsel,all
nsel,all
eplot
totaltime=Ton*1E-6
inctim=Ton*1E-6/lisannum
autots,1
*do,t,2e-6,totaltime,inctim
time,t  
deltim,inctim,inctim,inctim
kbc,0
outres,all,all,
tsres,erase
asel,s,loc,z,0,0
sfa,all,1,hflux,%flux2%  !施加高斯热源
solve
OUTRES,ALL,ALL,
/PSF,HFLUX,,2
/REPLOT
*enddo 　　　　　　
finish
疑问：表面效应单元生成了，但是并没有生成新的节点，是不是实体载荷和表面效应单元载荷都传到了节点上所以覆盖了？