焊接移动载荷不移动-附命令流

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

iomega

why not define the heat generation as:

2.1486*10^7*exp(-8.33*10^4*(( {X} - (X0+V*{TIME}) )^2+({Y}-Y0)^2+({Z}-Z0)^2))
X0, Y0 and Z0 are initial locations and V is the moving speed along the x direction.

Then you just need to apply the load once.

hjma123

谢谢指点，不过我把方程改了再加载，怎么还是不动呢？不得其解啊

woodcopper

貌似你的gaosi热源的中心没有随着时间的变化而变化啊！我看错了？

hjma123

问题可能就在这里！但是我不知道怎么才能让热源中心随时间变化。请多指点些。具体该加哪些参数进去？

woodcopper

将你的步增与gaosi热源的中心坐标绑定在一起，比如，原来坐标为x，时间为t，速度为v，则下一时刻的位置为x－v*t（或者step）那么中心就移动了，所以你在定义gaosi热源的时候应该把变量定义进去
不知道江明白没

ecjtuljsh

高手
都是厉害人

ljian25

受教了
谢谢

wuhong_0804

选节点和加载使用同一个局部坐标系吧，我最近也在做焊接移动热源，我的QQ：46971933，欢迎大家交流

wxpand2004

明白多了！谢谢

ihziuhoay

围观同学习

181809216

都要变化吧