hhb541788 发表于 2009-4-9 20:38

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

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

回复 楼主 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

守护蓝天 发表于 2011-1-3 15:30

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

sulkar 发表于 2011-1-3 17:39

支持一下,有同样的问题

lingyl2004 发表于 2014-6-6 15:01

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

jxxansys 发表于 2014-6-6 19:54

Ansys机械工程应用精华60例第四版
第46例单元生死应用实例—焊接
页: [1]
查看完整版本: 求助:如何在焊接模拟中加入移动热源