移动温度.txt

我做的一些ansys的练习,命令流,对学习ansys很有帮助

/prep7
/units,si 
*set,rb,0.01 
*set,length,0.2 
*set,width,0.12
*set,height,0.02
et,1,solid70
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 
block,0,length,0, height,0, width 
/view,1,1,1,1
LESIZE,1,,,height/rb
LESIZE,4,,,length/rb
LESIZE,9,,,width/rb
MSHKEY,1 
VATT,1,1,1 
VMESH,1 
Save
fini

/solu
antype,trans,new
trnopt,full
toffst,273
tunif,27 
tref,27
V=0.005
TINC=rb/V 
KK=3 
Qmax=1234567 
R=0.004 
MAX_TIME=length/rb+1 
MAX_X= length/rb+1 
MAX_Y=height/rb+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)*rb 
*ENDDO 
*DO,J,1,MAX_Y,1 
FLUX2(0,J,K)=(J-1)*rb 
*ENDDO 
*DO,I,1,MAX_X,1 
*DO,J,1,MAX_Y,1 
XCENTER=V*(K-1)*TINC 
YCENTER=0.01 
DISTANCE=SQRT(((I-1)*rb-XCENTER)**2+((J-1)*rb-YCENTER)**2) 
*IF,DISTANCE,LE,R,THEN 
FLUX2(I,J,K)=Qmax/EXP(KK*DISTANCE**2/(0.004**2)) 
*ELSE 
FLUX2(I,J,K)=0 
*ENDIF 
*ENDDO 
*ENDDO 
FLUX2(0,0,K)=(K-1)*TINC 
nsel,s,loc,z,width
sf,all,hflux,%flux2% 
allsel,all
!nsel,s,loc,x,0
!nsel,a,loc,x,length
!nsel,a,loc,y,0
!nsel,a,loc,y,height
!nsel,a,loc,z,0
sfa,1,,conv,110,25
sfa,3,,conv,110,25
sfa,4,,conv,110,25
sfa,5,,conv,110,25
sfa,6,,conv,110,25

allsel,all
nsubst,10,20,5
time,(k+1)*rb/v 
solve
*enddo