关于移动高斯热源的施加问题

hjma123

我正在做焊接热应力模拟，移动热源的过程中，发现热流只是分布在焊缝内部，即沿着焊缝有变化，而在焊缝外面一点热流也没有，温度场一点也没有扩展到母材上，不知道有没有遇到过相同问题的，急切盼望指点！！谢谢

xbzhang

参数设置有问题，或者是模型有问题。

hjma123

非常着急！我用的do循环语句有问题吗？热源就是不动，请高手们指点！谢谢了。下面是命令流，请帮忙看看哪里有问题。



/prep7
/units,si
length=0.06
v=0.01
step=0.002
count=length/step

ET,1,SOLID70


!材料参数
!建立模型 划分网格
!位移约束
!-->!定义高斯函数

*DEL,_FNCNAME   
*DEL,_FNCMTID   
*SET,_FNCNAME,'gaosi3'  
! /INPUT,gaosi3.func
*DIM,%_FNCNAME%,TABLE,6,26,1
!   
! Begin of equation: 2.1486*10^7*exp(-8.33*10^4*({X}^2+{Y}^2+{Z}^2))
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 10, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 7, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -1, 17, -2   
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 2.1486, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 2, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 8.33, 0, 0, -4  
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 10, 0, 0, 0
*SET,%_FNCNAME%(0,12,1), 0.0, -4, 0, 4, 0, 0, -1
*SET,%_FNCNAME%(0,13,1), 0.0, -5, 0, 1, -1, 17, -4  
*SET,%_FNCNAME%(0,14,1), 0.0, -1, 0, 1, -3, 3, -5   
*SET,%_FNCNAME%(0,15,1), 0.0, -3, 0, 2, 0, 0, 2
*SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 1, 2, 17, -3   
*SET,%_FNCNAME%(0,17,1), 0.0, -3, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,18,1), 0.0, -5, 0, 1, 3, 17, -3   
*SET,%_FNCNAME%(0,19,1), 0.0, -3, 0, 1, -4, 1, -5   
*SET,%_FNCNAME%(0,20,1), 0.0, -4, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,21,1), 0.0, -5, 0, 1, 4, 17, -4   
*SET,%_FNCNAME%(0,22,1), 0.0, -4, 0, 1, -3, 1, -5   
*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 3, -4   
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 7, 1, -3, 0, 0
*SET,%_FNCNAME%(0,25,1), 0.0, -3, 0, 1, -2, 3, -1   
*SET,%_FNCNAME%(0,26,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 2.1486*10^7*exp(-8.33*10^4*({X}^2+{Y}^2+{Z}^2))  
!-->


fini
/solu
!
antype,trans！设置分析选项
trnopt,full
nropt,auto, ,
eqslv,
!
solcontrol,on
autots,on
!
kbc,1  !step
outres,all,5,
!
!初始条件!
time,0.001
deltim,0.001,0.001,0.001
timint,off! time integration=off
tunif,25,
sfa,1,,conv,18,25!空气自然流动
sfa,2,,conv,18,25!空气自然流动
sfa,7,,conv,18,25!空气自然流动
sfa,8,,conv,18,25!空气自然流动
sfa,16,,conv,18,25!空气自然流动
sfa,17,,conv,18,25!空气自然流动
sfa,4,,conv,18,25!空气自然流动
sfa,14,,conv,18,25!空气自然流动
solve

!
WPSTYLE,,,,,,,,0

*do,I,0,count  !循环加载热源
disx=I*step
T=(disx+step)/v
Time,T
autots,1
antype,4
kbc,1
outres,basic,last
trnopt,full
nropt,full,,on
nsubst,5,,,1
lnsrch,1
tsres,erase
timint,1,therm

local,11,1,disx,0,0！定义局部坐标系，选节点
nsel,s,loc,y,-0.003,0.003
nsel,r,loc,x,0,0.002
nsel,r,loc,z,0

nplot
local,12,0,disx,0,0!定义局部坐标系，加载荷
sf,all,hflux,%gaosi3%
allsel,all
solve
sfdele,all,hflux
*enddo
time,30  !冷却时间
nsubst,5,,,1

solve
outres,all,5 !写入所有数据

finish

wxpand2004

论坛里有同样问题的你搜一下就出来了

tnanr

没有找到比较好的啊！

zraforgood

*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)*lsize
*enddo
*do,j,1,max_y,1
flux2(0,j,k)=(j-1)*lsize
*enddo
*do,i,1,max_x,1
*do,j,1,max_y,1
xcenter=v*(k-1)*tinc
ycenter=0.02
distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2)
*if,distance,le,r,then
flux2(i,j,k)=Qmax/exp(-3*kk*distance**2/(r**2))
        !高斯移动热源公式
*else
   flux2(i,j,k)=0
*endif
*enddo
*enddo
flux2(0,0,k)=(k-1)*tinc
*enddo