- 积分
- 0
- 注册时间
- 2003-2-22
- 仿真币
-
- 最后登录
- 1970-1-1
|
根据论坛上的例子,做了些改进,建立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 |
|