请教一个高斯热源移动加载的问题
<P>下面是我用移动高斯热源做热分析的命令流,该激光高斯模式为TEM11,在运行该程序时,flux2没有加载上去.如果用的是基模高斯光束,即把FLUX2(I,J,K)=8*qmax*(xcoord/r)**2*(ycoord/r)**2/EXP(KK*DISTANCE**2/(R**2)) 改写为FLUX2(I,J,K)=2*qmax/EXP(KK*DISTANCE**2/(R**2)) ,则分析结果正常.不知道为什么回出现这种"没有道理"的问题,请各位高手指点,不胜感激!</P><P>/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,solid70<BR>mp,dens,1,2200<BR>mp,prxy,1,0.17<BR>mptemp,1,20,500,1000<BR>mpdata,ex,1,1,7.78e10,8.2e10,8.52e10<BR>mptemp,1,100,200,400,600,800,1000<BR>mptemp,7,12000<BR>mpdata,c,1,1,660.4,836,973.9,1061.7,1086.8,1103.5<BR>mpdata,c,1,7,1128.6 <BR>mptemp,1,20,100,200,400,600,800<BR>mptemp,7,1000,1200<BR>mpdata,kxx,1,1,1.38,1.53,1.65,1.86,2.02,2.15<BR>mpdata,kxx,1,7,2.29,2.43<BR>mptemp,1,100,200,300,400,500,600<BR>mptemp,7,700,800,1000,1200<BR>mpdata,alpx,1,1,5.1e-7,5.85e-7, 6.27e-7,5.35e-7,6.12e-7,6.00e-7<BR>mpdata,alpx,1,7,5.71e-7,5.62e-7,5.56e-7,5.42e-7<BR>block,0,length,0, width, 0, height,<BR>/view,1,1,1,1<BR>LESIZE,1,,,width/rb<BR>LESIZE,4,,,length/rb<BR>LESIZE,9,,,heigh/rb<BR>MSHKEY,1 <BR>VATT,1,1,1 <BR>VMESH,1 <BR>Save<BR>fini<BR>/solu<BR>antype,trans,new<BR>timint,on<BR>autos,on<BR>outres,,all<BR>V=0.002 !光斑运动的速度<BR>TINC=rb/V <BR>KK=2 <BR>P=500<BR>R=0.0005 !光斑半径<BR>Qmax=p/(3.141593*r**2)</P>
<P>MAX_TIME=length/rb+1 <BR>MAX_X= length/rb+1 <BR>MAX_Y=width/rb+1 <BR>*DIM,FLUX2,TABLE,MAX_X,MAX_Y,MAX_TIME,X,Y,TIME <BR>*DO,K,1,MAX_TIME,1 <BR>*DO,I,1,MAX_X,1 <BR>FLUX2(I,0,K)=(I-1)*rb <BR>*ENDDO <BR>*DO,J,1,MAX_Y,1 <BR>FLUX2(0,J,K)=(J-1)*rb <BR>*ENDDO <BR>*DO,I,1,MAX_X,1 <BR>*DO,J,1,MAX_Y,1 <BR>XCENTER=V*(K-1)*TINC <BR>YCENTER=0.06<BR>Xcoord=(I-1)*rb-XCENTER<BR>Ycoord=(J-1)*rb-YCENTER<BR>DISTANCE=SQRT(xcoord**2+ycoord**2) <BR>*IF,DISTANCE,LE,2*R,THEN !在半径为2R的圆内加载高斯光束，其他的地方为零<BR>FLUX2(I,J,K)=8*qmax*(xcoord/r)**2*(ycoord/r)**2/EXP(KK*DISTANCE**2/(R**2)) !不同模式激光加载时的热流密度公式*ELSE <BR>FLUX2(I,J,K)=0 <BR>*ENDIF <BR>*ENDDO <BR>*ENDDO <BR>FLUX2(0,0,K)=(K-1)*TINC <BR>nsel,s,loc,z,height<BR>sf,all,hflux, %flux2% <BR>allsel,all<BR>nsel,s,loc,x,0<BR>sf,all,conv,110,25<BR>allsel,all<BR>nsel,s,loc,x,length<BR>sf,all,conv,110,25<BR>allsel,all<BR>nsel,s,loc,y,0<BR>sf,all,conv,110,25<BR>allsel,all<BR>nsel,s,loc,y,width<BR>sf,all,conv,110,25<BR>allsel,all<BR>nsel,s,loc,z,0<BR>sf,all,conv,110,25<BR>allsel,all<BR>time,(k+1)*rb/v<BR>nsubst,10,20,5<BR>solve<BR>*enddo<BR></P>
回复：（jameshust）请教一个高斯热源移动加载的问题
什么地方不对？我运行了一下能运行阿回复：（jameshust）请教一个高斯热源移动加载的问题
<P>你这个程序是有问题</P>回复：（jameshust）请教一个高斯热源移动加载的问题
程序错误的原因在于，你在定义边界条件的时候，高斯热源边界条件没加上去。<BR>另外，每个面上都用对流边界条件，这是错误的设置方法。 除了激光作用面外,其他面均是对流的边界条件,那这个对流边界条件该怎么加? 具体方法请你指教回复：（jameshust）请教一个高斯热源移动加载的问题
啊有gui操作撒？俺是菜鸟！<BR>谢谢
页:
[1]