请教一个高斯热源移动加载的问题

jameshust · 发表于 2006-3-2 10:17

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

下面是我用移动高斯热源做热分析的命令流,该激光高斯模式为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)) ,则分析结果正常.不知道为什么回出现这种"没有道理"的问题,请各位高手指点,不胜感激!
/prep7 /units,si *set,rb,0.01 !单元的尺寸（长度） *set,length,0.2 *set,width,0.12 *set,height,0.02 et,1,solid70 mp,dens,1,2200 mp,prxy,1,0.17 mptemp,1,20,500,1000 mpdata,ex,1,1,7.78e10,8.2e10,8.52e10 mptemp,1,100,200,400,600,800,1000 mptemp,7,12000 mpdata,c,1,1,660.4,836,973.9,1061.7,1086.8,1103.5 mpdata,c,1,7,1128.6 mptemp,1,20,100,200,400,600,800 mptemp,7,1000,1200 mpdata,kxx,1,1,1.38,1.53,1.65,1.86,2.02,2.15 mpdata,kxx,1,7,2.29,2.43 mptemp,1,100,200,300,400,500,600 mptemp,7,700,800,1000,1200 mpdata,alpx,1,1,5.1e-7,5.85e-7, 6.27e-7,5.35e-7,6.12e-7,6.00e-7 mpdata,alpx,1,7,5.71e-7,5.62e-7,5.56e-7,5.42e-7 block,0,length,0, width, 0, height, /view,1,1,1,1 LESIZE,1,,,width/rb LESIZE,4,,,length/rb LESIZE,9,,,heigh/rb MSHKEY,1 VATT,1,1,1 VMESH,1 Save fini /solu antype,trans,new timint,on autos,on outres,,all V=0.002 !光斑运动的速度 TINC=rb/V KK=2 P=500 R=0.0005 !光斑半径 Qmax=p/(3.141593*r**2)
MAX_TIME=length/rb+1 MAX_X= length/rb+1 MAX_Y=width/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.06 Xcoord=(I-1)*rb-XCENTER Ycoord=(J-1)*rb-YCENTER DISTANCE=SQRT(xcoord**2+ycoord**2) *IF,DISTANCE,LE,2*R,THEN !在半径为2R的圆内加载高斯光束，其他的地方为零 FLUX2(I,J,K)=8*qmax*(xcoord/r)**2*(ycoord/r)**2/EXP(KK*DISTANCE**2/(R**2)) !不同模式激光加载时的热流密度公式*ELSE FLUX2(I,J,K)=0 *ENDIF *ENDDO *ENDDO FLUX2(0,0,K)=(K-1)*TINC nsel,s,loc,z,height sf,all,hflux, %flux2% allsel,all nsel,s,loc,x,0 sf,all,conv,110,25 allsel,all nsel,s,loc,x,length sf,all,conv,110,25 allsel,all nsel,s,loc,y,0 sf,all,conv,110,25 allsel,all nsel,s,loc,y,width sf,all,conv,110,25 allsel,all nsel,s,loc,z,0 sf,all,conv,110,25 allsel,all time,(k+1)*rb/v nsubst,10,20,5 solve *enddo

frogfish · 发表于 2006-3-4 21:13

什么地方不对？我运行了一下能运行阿

yy_lhsx · 发表于 2006-4-11 22:06

你这个程序是有问题

yy_lhsx · 发表于 2006-4-11 22:08

程序错误的原因在于，你在定义边界条件的时候，高斯热源边界条件没加上去。 另外，每个面上都用对流边界条件，这是错误的设置方法。

jameshust · 发表于 2006-4-21 11:35

除了激光作用面外,其他面均是对流的边界条件,那这个对流边界条件该怎么加? 具体方法请你指教

beckshjh · 发表于 2006-4-27 15:36

啊有gui操作撒？俺是菜鸟！ 谢谢

账号		自动登录	找回密码
密码			我要加入

请教一个高斯热源移动加载的问题

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（jameshust）请教一个高斯热源移动加载的问题

回复：（jameshust）请教一个高斯热源移动加载的问题

回复：（jameshust）请教一个高斯热源移动加载的问题

回复：（jameshust）请教一个高斯热源移动加载的问题