马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
如题。本人正在做一个焊接模拟,因为想要获得比较可靠的数据,所以想施加移动热源。查了一个文献,里面加载的是移动的高斯热源。由于本人水平有限,没弄明白,现将该程序段复制于下方,希望各位前辈老师指点迷津!或者说有其它建议,望不吝赐教!感激不尽!!!
/solu
antype,trans
trnopt,full
nropt,auto,,
eqslv,
solcontrol,on
autots,on
kbc,1
Load step,apply moving heat source
timint,on
*do,t,tinc1,t1,tinc1
time,t
x1=x1+0.002
autots,0
deltim,tinc1/2,tinc1/2,tinc1/2
outres,all,5
WPSTYLE,,,,,,,,1
CSYS,4
KWPAVE,1609
wpoff,x1
CSWPLA,11,1,1,1,
ASEL,S,,,20
Nsla,r,1
Nsla,r,loc,x,0,r/4+remain
Nplot
bf,all,hgen,hgen1
allsel,all
ASEL,S,,,20
Nsla,r,1
Nsla,r,loc,x,r/4+remain,r/2+remain
nplot
bf,all,hgen,hgen2
allsel,all
ASEL,S,,,20
Nsla,r,1
Nsla,r,loc,x,r/2+remain,3*r/4+remain
nplot
bf,all,hgen,hgen3
allsel,all
ASEL,S,,,20
Nsla,r,1
Nsla,r,loc,x,3*r/4+remain,r+remain
nplot
bf,all,hgen,hgen4
allsel,all
solve
ASEL,S,,,20
Nsla,r,1
Nsla,r,loc,x,0,r/4+remain
Nplot
bf,all,hgen,le-10
allsel,all
ASEL,S,,,20
Nsla,r,1
Nsla,r,loc,x,r/4+remain,r/2+remain
nplot
bf,all,hgen,le-10
allsel,all
ASEL,S,,,20
Nsla,r,1
Nsla,r,loc,x,r/2+remain,3*r/4+remain
nplot
bf,all,hgen,le-10
allsel,all
ASEL,S,,,20
Nsla,r,1
Nsla,r,loc,x,3*r/4+remain,r+remain
nplot
bf,all,hgen,le-10
allsel,all
*enddo
/solu
antype,trans,new
nlgeom,on
timint,1,struct
timint,1,thern
timintm,0,mag
timint,0,elect
tintp,0.005,-1,0.5,0.2
tref,25
nropt,auto,,
autos,on
outres,,all
v=0.008
tinc=rb/v
U=25
I=180
ef=0.75
p=25*180*ef
r=0.006
qmax=p/(3.14159*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
distance=sqrt(((i-1)*rb-xcenter)**2+((j-1)*rb-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=qmax/EXP(3*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,10,20
allsel,all
nsel,s,loc,x,length
sf,all,conv,10,20
allsel,all
nsel,s,loc,y,0
sf,all,conv,10,20
allsel,all
nsel,s,loc,y,width
sf,all,conv,10,20
allsel,all
nsel,s,loc,z,0
sf,all,conv,10,20
allsel,all
time,k*rb/v
nsubst,5
solve
*enddo
save
time,20
nsubst,10
solve
time,44
nsbust,24
solve
time,644
nsbust,20
solve
time,7244
nsubst,100
solve
save
[ 本帖最后由 hhb541788 于 2009-4-9 20:42 编辑 ] |