求助：移动高斯热源加载在圆柱侧面的问题

baiqinghua

根据论坛上的例子，做了些改进，建立1/4圆周的模型，然后施加高斯热源，可是结果只有1/8是正确的，问题始终不能解决，请教高手！！希望能解决啊！

命令流如下：
/prep7  
/units,si   
*afun,deg   
*set,element,0.005  
*set,inradius,0.05  
*set,outradius,0.06
*set,degree,90  
*set,length,0.15
et,1,solid70
mp,dens,1,2640  
mptemp,1,20,100,200,400,600,800   
mptemp,7,1000,1200   
mpdata,c,1,1,660.4,660.4,836,973.9,1061.7,1086.8,1103.5   
mpdata,c,1,7,1128.6
mpdata,enth,1,1,0,7.583e4,2.592e5,4.619e5,6.875e5,9.528e5   
mpdata,enth,1,7,1.278e6
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
csys,1  
k,1,inradius,0,0
k,2,inradius,degree,0   
k,3,inradius,degree,length  
k,4,inradius,0,length   
k,5,outradius,0,0   
k,6,outradius,degree,0  
k,7,outradius,degree,length
k,8,outradius,0,length  
v,1,2,3,4,5,6,7,8   
/view,1,1,1,1   
lesize,5,,,10   
lesize,6,,,90   
lesize,8,,,60   
mshkey,1
vatt,1,1,1  
vmesh,1
save


*set,delt,1
*SET,qmax,10000000  
*SET,r,0.005
*SET,pi,3.14159
*SET,lsize,2
*SET,v,90   
*SET,tinc,lsize/v   
*SET,Y0,0
*dim,flux2,table,44,30,44,y,z,time  
*do,k,1,44,1
  *do,i,1,44,1
  *SET,flux2(i,0,k),(i-1)*lsize*pi*0.06/180   
  *enddo  
  *do,j,1,30,1
  *SET,flux2(0,j,k),(j-1)*0.005   
  *enddo  
  *do,i,1,44,1
     *do,j,1,30,1
     *SET,ycenter,v*(k-1)*tinc   
     *SET,zcenter,0   
     *afun,deg   
     *SET,distance,sqrt((((i-1)*lsize-ycenter)/(180)*pi*0.06)**2+((j-1)*0.005-zcenter)**2)   
         *if,distance,le,2*r,then
            *SET,flux2(i,j,k),qmax/EXP(DISTANCE**2/(r**2))  
         *else   
            *SET,flux2(i,j,k),0
         *endif  
      *enddo  
    *enddo  
    *SET,flux2(0,0,k),k*tinc  !
*enddo


/SOL
antype,trans,new
trnopt,full
timint,on
tunif,20
nropt,auto,,
autos,on
outres,,all
kbc,1   
allsel,all

*do,t,tinc,1,tinc
*SET,Y0,y0+2
time,t
solcontrol,on
timint,on
autots,on
lnsrch,1
nsubst,5
kbc,1
outres,all,5
tsres,erase

Nsel,s,loc,x,0.06
Nsel,r,loc,y,y0-1-1e-6,y0+1+1e-6  
nplot
Sf,all,hflux,%flux2%
Allsel,all  
outres,all,ALL
SOLVE
nsel,s,loc,x,0.06
Nsel,r,loc,y,y0-1-1e-6,y0+1+1e-6  
sfdele,all,hflux,%flux2%
allsel,all
*enddo

fini

baiqinghua

我实在是没办法了啊，呵呵

gus1977

仔细检查一下你的加载条件，是否只加对了1/8，因为只有1/8对，说明你的模型没问题，那问题肯定出在你的边界条件上了

zxw19851015

来学习下,感觉搞模拟好复杂啊

yanqing0680

我最近也在做这方面工作
可以交流一下
我的QQ:278770766

303627564

看看刚开始建立模型时有没有设置成轴对称啊

Vicly2012

新手路过，表示不懂，观摩一下