zzgjtx 发表于 2009-4-16 16:45

移动热源加载问题

大家好,最近做了移动热源的问题,加载使用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 编辑 ]

守护蓝天 发表于 2011-1-3 15:32

问题解决了吗?可以说出来解决方法吗?
页: [1]
查看完整版本: 移动热源加载问题