声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1917|回复: 5

[综合讨论] 求助:如何在焊接模拟中加入移动热源

[复制链接]
发表于 2009-4-9 20:38 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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 编辑 ]
回复
分享到:

使用道具 举报

 楼主| 发表于 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 | 显示全部楼层
只是从求解开始的,怎么没前处理的程序段呢?
发表于 2011-1-3 17:39 | 显示全部楼层
支持一下,有同样的问题
发表于 2014-6-6 15:01 | 显示全部楼层
有关焊接模拟的资料比较少。回复的人也不多
发表于 2014-6-6 19:54 | 显示全部楼层
Ansys机械工程应用精华60例第四版
第46例  单元生死应用实例—焊接
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 00:39 , Processed in 0.070403 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表