热应力分析的一点问题

jameshust

<P  align=left>移动热源的热应力分析<BR>/prep7<BR>/units,si <BR>*set,rb,0.01 <BR>*set,length,0.2 <BR>*set,width,0.12<BR>*set,height,0.02<BR>et,1,solid5<p></p></P>
<P  align=left>mp,alpx,1,1.06e-5<p></p></P>
<P  align=left>mp,ex,1,1.75e11<p></p></P>
<P  align=left>mp,prxy,1,0.3<BR>mptemp,1,20,250,500,750,1000,1500,1700,2500 <p></p></P>
<P  align=left>mpdata,c,1,1,460,480,530,675,670,660,780,820 <p></p></P>
<P  align=left>mpdata,kxx,1,1,50,47,40,27,30,35,140,142 <p></p></P>
<P  align=left>mpdata,dens,1,1,7820,7700,7610,7550,7490,7350,7300,7090 <p></p></P>
<P  align=left>mpdata,hf,1,1,100,350,520,1000,1500,3000,3100,3500 <p></p></P>
<P  align=left>block,0,length,0, height,0, width  <p></p></P>
<P  align=left>/view,1,1,1,1<p></p></P>
<P  align=left>LESIZE,1,,,height/rb<p></p></P>
<P  align=left>LESIZE,4,,,length/rb<p></p></P>
<P  align=left>LESIZE,9,,,width/rb<p></p></P>
<P  align=left>MSHKEY,1 <p></p></P>
<P  align=left>VATT,1,1,1 <p></p></P>
<P  align=left>VMESH,1 <p></p></P>
<P  align=left>Save<p></p></P>
<P  align=left>fini<p></p></P>
<P  align=left>/solu<p></p></P>
<P  align=left>antype,trans,new<BR>trnopt,full<p></p></P>
<P  align=left>timint,1,struct<p></p></P>
<P  align=left>timint,1,thern<p></p></P>
<P  align=left>timintm,0,mag<p></p></P>
<P  align=left>timint,0,elect<p></p></P>
<P  align=left>tintp,0.005,,,-1,0.5,0.2<p></p></P>
<P  align=left>outres,all<BR>toffst,273<BR>tunif,27  <BR>tref,27<p></p></P>
<P  align=left>V=0.005<p></p></P>
<P  align=left>TINC=rb/V <p></p></P>
<P  align=left>KK=3 <p></p></P>
<P  align=left>Qmax=1234567 <p></p></P>
<P  align=left>R=0.004 <p></p></P>
<P  align=left>MAX_TIME=length/rb+1 <p></p></P>
<P  align=left>MAX_X= length/rb+1 <p></p></P>
<P  align=left>MAX_Y=height/rb+1 <p></p></P>
<P  align=left>*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME <p></p></P>
<P  align=left>*DO,K,1,MAX_TIME,1 <p></p></P>
<P  align=left>*DO,I,1,MAX_X,1 <p></p></P>
<P  align=left>FLUX2(I,0,K)=(I-1)*rb <p></p></P>
<P  align=left>*ENDDO <p></p></P>
<P  align=left>*DO,J,1,MAX_Y,1 <p></p></P>
<P  align=left>FLUX2(0,J,K)=(J-1)*rb <p></p></P>
<P  align=left>*ENDDO <p></p></P>
<P  align=left>*DO,I,1,MAX_X,1 <p></p></P>
<P  align=left>*DO,J,1,MAX_Y,1 <p></p></P>
<P  align=left>XCENTER=V*(K-1)*TINC <p></p></P>
<P  align=left>YCENTER=0.01 <p></p></P>
<P  align=left>DISTANCE=SQRT(((I-1)*rb-XCENTER)**2+((J-1)*rb-YCENTER)**2) <p></p></P>
<P  align=left>*IF,DISTANCE,LE,R,THEN <p></p></P>
<P  align=left>FLUX2(I,J,K)=Qmax/EXP(KK*DISTANCE**2/(0.004**2)) <p></p></P>
<P  align=left>*ELSE <p></p></P>
<P  align=left>FLUX2(I,J,K)=0 <p></p></P>
<P  align=left>*ENDIF <p></p></P>
<P  align=left>*ENDDO <p></p></P>
<P  align=left>*ENDDO <p></p></P>
<P  align=left>FLUX2(0,0,K)=(K-1)*TINC <p></p></P>
<P  align=left>nsel,s,loc,z,width<BR>sf,all,hflux,%flux2% <p></p></P>
<P  align=left>allsel,all<p></p></P>
<P  align=left>nsel,s,loc,x,0<p></p></P>
<P  align=left>nsel,s,loc,x,length<p></p></P>
<P  align=left>nsel,s,loc,y,0<p></p></P>
<P  align=left>nsel,s,loc,y,height<p></p></P>
<P  align=left>nsel,s,loc,z,0<p></p></P>
<P  align=left>sf,all,conv,110,25<p></p></P>
<P  align=left>allsel,all<BR>time,(k+1)*rb/v<p></p></P>
<P  align=left>nsubst,10,20,5<p></p></P>
<P  align=left>solve<BR>*enddo <p></p></P>
<P  align=left>fini<BR><BR>在运行该程序的过程中回出现这样的警告：SF is not a recognized GEGIN command,abbreviation,or macro.<BR>this command will be ingored.<BR>请问问题出在哪里？谢谢 <BR><BR ><BR ><p></p></P>
<P ><p><FONT face="Times New Roman"> </FONT></p></P>
<P ><p><FONT face="Times New Roman"> </FONT></p></P>