[转贴]求教：一论坛上的移动热源的例子~谢谢！

8145113

如果我这个时候想只在热源中心再施加一个初始的温度载荷，应该在哪更改此命令流：
/prep7
i=250
u=36
v=0.012
lx=0.4
ly=0.2
lz=0.006
x1=-0.002
ET,1,SOLID70
ET,2,SOLID90
BLOCK,0,lx,0,ly,-lz,0,
BLOCK,0,lx,0,0.008,-lz,0,   
BLOCK,0,lx,0.008,0.014,-lz,0,  

FLST,2,3,6,ORDE,2   
FITEM,2,1   
FITEM,2,-3  
VOVLAP,P51X

LESIZE,4,0.001, , , , , , ,1
LESIZE,5,0.001, , , , , , ,1
LESIZE,9,0.002, , , , , , ,1
LESIZE,10,0.002, , , , , , ,1
LESIZE,14,0.001, , , , , , ,1
LESIZE,19,0.001, , , , , , ,1
LESIZE,23,0.002, , , , , , ,1
LESIZE,24,0.002, , , , , , ,1
LESIZE,37,0.001, , , , , , ,1
LESIZE,38,0.001, , , , , , ,1
LESIZE,39,0.001, , , , , , ,1
LESIZE,40,0.001, , , , , , ,1

  
LESIZE,2,0.006, , , , , , ,1
LESIZE,7,0.006, , , , , , ,1
LESIZE,26,0.006, , , , , , ,1
LESIZE,31,0.006, , , , , , ,1
LESIZE,41,0.006, , , , , , ,1
LESIZE,42,0.006, , , , , , ,1
LESIZE,43,0.006, , , , , , ,1
LESIZE,44,0.006, , , , , , ,1

LESIZE,11,0.003, , , , , , ,1
LESIZE,12,0.003, , , , , , ,1
LESIZE,35,0.003, , , , , , ,1
LESIZE,36,0.003, , , , , , ,1   
VMESH,4
VMESH,5
!*  
TYPE,   2   

MSHKEY,0
MSHAPE,1,3d
CM,_Y,VOLU  
VSEL, , , ,       6
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  


r=0.008

remain=r/10

pi=3.1415926

t1=lx/v

tinc1=0.002/v

ef=0.75
Q=ef*i*u      !热源功率
qm=3*Q/(pi*r**2)
q1=qm*(1-0.82902)*(pi*r**2)/3
q2=qm*(0.82902-0.47236)*(pi*r**2)/3
q3=qm*(0.47236-0.18498)*(pi*r**2)/3
q4=qm*(0.18498-0.04978)*(pi*r**2)/3
hgen1=q1/(0.002*pi*(r/4)**2)
hgen2=q2/(0.002*pi*((r/2)**2-(r/4)**2))
hgen3=q3/(0.002*pi*((3*r/4)**2-(r/2)**2))
hgen4=q4/(0.002*pi*(r**2-(3*r/4)**2))

!mp,ex,1,1.93E+5! N/mm^2
!mp,nuxy,1,0.29 ! unitless
mptemp,1,0,20,100,300,600,765,1000
mpdata,kxx,1,1,41.55,42.95,38.23,30.81,38.38,25.39,25.57!W/(m.K)
mpdata,c,1,1,462,481,530,680,1360,718,604!J/kg.k
mpdata,dens,1,1,7810,7800,7700,7610,7550,7490!kg/mm**3

fini
/solu
!
antype,trans
trnopt,full
nropt,auto, ,
eqslv,
!
solcontrol,on
autots,on
!
kbc,1  !step
outres,all,5
!
!load step 1 ,initial condition  25  degrees
!
time,0.001
deltim,0.001,0.001,0.001
timint,off! time integration=off
tunif,25,
sfa,4,,conv,10,25!空气自然流动
sfa,20,,conv,10,25!空气自然流动
sfa,21,,conv,10,25!空气自然流动
sfa,22,,conv,10,25!空气自然流动
sfa,23,,conv,10,25!空气自然流动
sfa,24,,conv,10,25!空气自然流动
sfa,25,,conv,10,25!空气自然流动
sfa,28,,conv,10,25!空气自然流动
sfa,29,,conv,10,25!空气自然流动
sfa,30,,conv,10,25!空气自然流动
sfa,19,,conv,10,25!空气自然流动
sfa,27,,conv,10,25!空气自然流动
solve
!
!load step 2,  apply moving heat source
timint,on! time integration=on
!

*do,t,tinc1,t1,tinc1

time,t
x1=x1+0.002
autots,0
deltim,tinc1/2,tinc1/2,tinc1/2
kbc,1
outres,all,5
tsres,erase


WPSTYLE,,,,,,,,1
CSYS,4  
KWPAVE,      1609  
wpoff,x1
CSWPLA,11,1,1,1,

ASEL,S, , ,20
nsla,r,1
nsel,r,loc,x,0,r/4+remain
nplot
bf,all,hgen,hgen1

allsel,all

ASEL,S, , ,20
nsla,r,1
nsel,r,loc,x,r/4-remain,r/2+remain
nplot
bf,all,hgen,hgen2

allsel,all

ASEL,S, , ,20
nsla,r,1
nsel,r,loc,x,r/2-remain,3*r/4+remain
nplot
bf,all,hgen,hgen3

allsel,all

ASEL,S, , ,20
nsla,r,1
nsel,r,loc,x,3*r/4-remain,r+remain
nplot
bf,all,hgen,hgen4

allsel,all
solve

ASEL,S, , ,20
nsla,r,1
nsel,r,loc,x,0,r/4+remain
nplot
bf,all,hgen,1e-10

allsel,all

ASEL,S, , ,20
nsla,r,1
nsel,r,loc,x,r/4-remain,r/2+remain
nplot
bf,all,hgen,1e-10

allsel,all

ASEL,S, , ,20
nsla,r,1
nsel,r,loc,x,r/2-remain,3*r/4+remain
nplot
bf,all,hgen,1e-10

allsel,all

ASEL,S, , ,20
nsla,r,1
nsel,r,loc,x,3*r/4-remain,r+remain
nplot
bf,all,hgen,1e-10

allsel,all

*enddo
!
time,7200

nsubst,40
BFDELE, all, hgen

solve
OUTRES,all,5

fini

tnanr

不晓得!有会的吗？