移动热源加载问题

zzgjtx

大家好，最近做了移动热源的问题，加载使用APDL编程，使用TABLE二维数组定义载荷部分程序如下：
finish
/clear
/PREP7  
!*  
ET,1,PLANE55  
!*  
!*  
MPTEMP,1,20,100,200,300,400,500,600,700,765,800,900,1000  
MPDATA,KXX,1,,41.84,40.88,39.58,37.69,35.35,32.15,28.44,24.61,26.77,17.38,22.07,24.88
MPTEMP,1,20,100,200,300,400,500,600,700,765,800,900,1000  
MPDATA,C,1,,467,491,527,565,611,678,787,1050,1423,672,634,649
MPTEMP,,,,,,,,  
MPTEMP,1,20  
MPDATA,DENS,1,,7850
K,1,0,0,0,  
K,2,4.33e-2,0,0,
K,3,0.0433,-1e-2,0,
K,4,0,-0.01,0,  
FLST,2,4,3  
FITEM,2,1   
FITEM,2,2   
FITEM,2,3   
FITEM,2,4   
A,P51X
LPLOT
LESIZE,1,,,100
LESIZE,3,,,100
LESIZE,4,,,20,0.2
LESIZE,2,,,20,5
AMESH,all  
FINISH  
save
*DIM,Grind,TABLE,51,51,,xpath,Time,
*DO,j,1,51,1
*DO,i,1,51,1
Grind(i,0)=(i-1)*0.866e-3
*enddo
*do,k,1,51,1
*if,k,lt,j,then
  *if,j-k,lt,11,then
  Grind(k,j)=40800000-(j-k)*4080000
  *else
  Grind(k,j)=0
  *endif
*endif
*if,k,eq,j,then
  Grind(k,j)=40800000
*endif
*if,k,gt,j,then
Grind(k,j)=0
*ENDIF
*ENDDO
Grind(0,j)=(j-1)*0.289e-1
*ENDDO
finish
/solu
ANTYPE,TRANS
TUNIF,20
KBC,1  
*SET,DELT,1
TIMINT,ON
AUTOTS,OFF  

*DO,TM,1E-6,1.445,0.289e-1
/REP,FAST  
/GO
!*
TIME,TM
AUTOTS,0

KBC,0
!*
TSRES,ERASE
/GO
!*
!*
LSEL,S,LINE,,1
SFL,all,,HFLUX,%Grind%
/STATUS,SOLU
SOLVE
allsel,all
*ENDDO  
FINISH  

出现的问题是：采用此方法三维数组的加载没有问题，已经验证，但是做平面的加不上载荷，请高手指教，谢谢！

[ 本帖最后由 zzgjtx 于 2009-4-16 16:52 编辑 ]

守护蓝天

问题解决了吗？可以说出来解决方法吗？