ansys高斯热源焊接模拟,焊接瞬态分析部分,请问大神,焊接速度怎么设置呢
j=1 !定义求解时间
timint,on
*do,i,1,20
esel,s,mat,,1 !只选择材料1的单元
nsle,s,1 !选择当前单元的节点
nsel,r,loc,z,(zmin+zmax)/2-dz,(zmin+zmax)/2+dz !在已选节点中重新去选
nsel,r,loc,y,2.5
nsel,r,loc,x,-dx,dx !在大正方形中选择出中心的正方形
nplot
sf,all,hflux,qm !中心高斯热源
r=1 !给个初值
*do,m,1,3
qr=qm*exp(-3*r*r/16) !高斯热源计算公式,16为等离子体云根部半径的平方(有效加热面)。
esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,2.5
nsel,r,loc,x,-(m+1)*dx,(m+1)*dx
nsel,r,loc,z,(zmin+zmax)/2-(m+1)*dz,(zmin+zmax)/2-m*dz
nplot
sf,all,hflux,qr
esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,2.5
nsel,r,loc,x,-(m+1)*dx,(m+1)*dx
nsel,r,loc,z,(zmin+zmax)/2+m*dz,(zmin+zmax)/2+(m+1)*dz
nplot
sf,all,hflux,qr
esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,2.5
nsel,r,loc,z,(zmin+zmax)/2-(m+1)*dz,(zmin+zmax)/2+(m+1)*dz
nsel,r,loc,x,-(m+1)*dx,-m*dx
nplot
sf,all,hflux,qr
esel,s,mat,,1
nsle,s,1
nsel,r,loc,y,2.5
nsel,r,loc,z,(zmin+zmax)/2-(m+1)*dz,(zmin+zmax)/2+(m+1)*dz
nsel,r,loc,x,m*dx,(m+1)*dx
nplot
sf,all,hflux,qr
r=r+1
*enddo
time,j
deltim,j,j,j
alls
solve
esel,s,mat,,1
nsle,s,1
sfdele,all,hflux
j=j+1
zmin=zmin+zd
zmax=zmax+zd
*enddo
alls
eplot
返回小木虫查看更多
京公网安备 11010802022153号
你的命令流只是一部分,无法进行分析。说一点看法吧:
你施加热载荷的命令流部分是:
esel,s,mat,,1 !只选择材料1的单元
nsle,s,1 !选择当前单元的节点
nsel,r,loc,z,(zmin+zmax)/2-dz,(zmin+zmax)/2+dz !在已选节点中重新去选
nsel,r,loc,y,2.5
nsel,r,loc,x,-dx,dx !在大正方形中选择出中心的正方形
nplot
sf,all,hflux,qm !中心高斯热源
其中通过 zmin 和 zmax 来确定受热的节点,改变 zmin 和 zmax,就可以使加载位置移动。
一般情况,焊接速度是已知的,两个载荷步之间的 zmin 和 zmax 的改变,除以两个载荷步的时间间隔,应该等于焊接速度。
就是说,通过移动载荷和改变丘吉尔终止时间,就能够实现所需的焊接速度,
最后一句有误,修改如下:
就是说,通过移动载荷和适当设置相邻载荷步的终止时间,就能够实现所需的焊接速度。
不知道你的算例复杂在哪里,我也没有时间详细研究你的问题。给你一些资料自己看看吧:
谢谢师兄
ly=0.12
lx=0.4
lz=0.012
lsize=0.01
v=0.00185
tm=lx/v
tinc=lsize/v
max_time=1+tm/tinc
max_x=1+lx/lsize
max_y=1+ly/lsize
师兄,tinc表示每个单元载荷的加载时间,这里定义的max_time是指总共需要多少步吗?还有定义的tm是总的热源移动时间,后面也没有用到tm呀,为什么定义了max_time还要定义tm呢?搞不懂
算例是这个文件三维高斯热源apdl命令流.txt