求助：如何在焊接模拟中加入移动热源

hhb541788

如题。本人正在做一个焊接模拟，因为想要获得比较可靠的数据，所以想施加移动热源。查了一个文献，里面加载的是移动的高斯热源。由于本人水平有限，没弄明白，现将该程序段复制于下方，希望各位前辈老师指点迷津！或者说有其它建议，望不吝赐教！感激不尽！！！
/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 编辑 ]

hhb541788

上面包含两种方法，前面部分是第一种
/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

守护蓝天

只是从求解开始的，怎么没前处理的程序段呢?

sulkar

支持一下，有同样的问题

lingyl2004

有关焊接模拟的资料比较少。回复的人也不多

jxxansys

Ansys机械工程应用精华60例第四版
第46例  单元生死应用实例—焊接