多道焊中，焊接一道焊缝后的自然冷却问题~请求高人指点。

bloodyforest

多道焊中，焊接一道焊缝后的自然冷却问题~请求高人指点。


我做一个焊件的焊补模拟， 采用多道焊， 我想把第一道焊缝焊完后的自然冷却过程考虑进去，不知道为什么我的命令算出来有问题，请各位高手指点。我初学，命令流比较简单，希望大家有空可以帮我调试一下。

命令流我贴出来， 下面是焊完第一道焊缝：

!******************************第一部分进行温度场分析******************************
/clear
/com, units: SI, m - degree C - W - kg -sec - J
/title,the crack of welding (V fault)
/filname,the v crack of welding
/prep7
et,1,solid90
mptemp,1,20,500,800,1200,1500                           
mpdata,ex,1,,1.068e11,0.95e11,0.87e11,0.54e11,0.40e11     !*****定义弹性模量
mpdata,ex,2,,1.068e11,0.95e11,0.87e11,0.54e11,0.40e11

mpdata,enth,1,,15,23,50,68,90      !********定义焓
mpdata,enth,2,,15,23,50,68,90
tb,bkin,1,5                        !********定义应力计算的参数
tbtemp,20,1
tbdata,1,315e6,0.209e11
tbtemp,500,2
tbdata,1,212e6,0.172e11
tbtemp,800,3
tbdata,1,153e6,0.133e11
tbtemp,1200,4
tbdata,1,93e6,0.084e11
tbtemp,1500,5
tbdata,1,33e6,0.045e11
tb,bkin,2,5               
tbtemp,20,1
tbdata,1,314e6,0.216e11
tbtemp,500,2
tbdata,1,203e6,0.178e11
tbtemp,800,3
tbdata,1,153e6,0.142e11
tbtemp,1200,4
tbdata,1,63e6,0.104e11
tbtemp,1500,5
tbdata,1,13e6,0.068e11
mp,dens,1,8320                   !定义密度
mp,dens,2,8320
mp,alpx,1,20e-6                  !定义线膨胀系数
mp,alpx,2,20e-6
mp,prxy,1,0.337                  !定义泊松比
mp,prxy,2,0.337
mp,kxx,1,80                      !定义热导率
mp,kxx,2,80
mp,c,1,404                       !定义比热
mp,c,2,404
!************************建立模型******************
k,,0.05,0.02,0.075
k,,0.05,0.02,0.125
k,,0.05,0.0125,0.075
k,,0.05,0.0125,0.125
k,,0.035,0.02,0.075
k,,0.035,0.02,0.125
k,,0.0425,0.0125,0.075
k,,0.0425,0.0125,0.125
l,1,2
l,2,4
l,4,3
l,3,1
l,5,6
l,6,8
l,8,7
l,7,5
l,1,5
l,3,7
l,2,6
l,4,8
lglue,all
al,1,2,3,4
al,5,6,7,8
al,1,9,5,11
al,3,10,7,12
al,2,11,6,12
al,4,9,8,10
aglue,all
va,1,2,3,4,5,6
block,0,0.05,0,0.02,0,0.2
vovlap,all
lesize,1,,,20
lesize,3,,,20
lesize,5,,,20
lesize,7,,,20
lesize,2,,,6
lesize,4,,,6
lesize,6,,,6
lesize,8,,,6
lesize,9,,,6
lesize,10,,,6
lesize,11,,,6
lesize,12,,,6
mat,1
type,1
vmesh,1
esize,0.005
mshkey,0
mshape,1
mat,2
type,1
vmesh,3
                                 
eplot
/view,1,1,1,1
/replot
!**********************死亡*******************************
/solu
antype,4                                   
trnopt,full
alls
tunif,25          !设置统一温度
/auto,1
/replot
alls
eplot
esel,s,mat,,1
eplot
ekill,all       !杀死所有被焊区域
alls
esel,s,live
eplot
outres,all,all
!*******************激活与加载*************************
antype,trans
trnopt,full                    !默认的方法
solcontrol,on                  !默认
autots,on                      !自动步长
tintp,0.005,,,1,0.5,0.2        !向后积分
timint,on                      !瞬态效果
keyopt,,1                      !打开对角比热容矩阵
lnsrch,on                      !打开线性搜索
zmin=0.075                     !定义焊缝Z方向被划分的第一组小单元的坐标
zmax=0.0775
ymin=0.0125                    !同上Y方向
ymax=0.0125+0.0075/2
xmin=0.0425                    !同上X方向
xmax=0.05
zd=0.0025                     !坐标增量
yd=0.0075/2
xd=0.0075/2
!*************************第一层焊缝
j=0.5
timint,off                        !稳态效果
esel,s,mat,,1                     !只选择第一组施焊单元
nsle,s,all
nsel,r,loc,z,zmin,zmax
nsel,r,loc,y,ymin,ymax
esln,s,all                        
ealive,all                        !激活第一组焊缝单元
eplot
nsle,all
d,all,temp,25                 !定义焊缝单元上的初始温度
nsel,inve                        
d,all,temp,25                    !定义除了第一组焊缝单元其他所有单元节点上的温度(此处用D的做法是参考书上的,所以没用IC)
asel,all
asel,u,area,,1,6               
sfa,all,,conv,70,20              !施加除了焊缝区的面其他所有面的对流换热系数
time,1e-3                        !0.001S内求解稳态,确立初始温度场
kbc,1
allsel,all
solve
eplot
ddele,all,temp                   !删除所有焊缝和母材上所有节点温度
esel,s,mat,,1                    !再选择第一组焊缝单元
nsle,s,all
nsel,r,loc,z,zmin,zmax
nsel,r,loc,y,ymin,ymax
esln,s,all
ealive,all
eplot
nsle,all                        
d,all,temp,1500                  !加载1500度

timint,on                         !进行瞬态分析
antype,,rest                      !在上一步基础上求解
time,j                            !时间设置成0.5
deltim,0.5-1e-3,0.5-1e-3,0.5-1e-3 !一步求解到位
kbc,1
outres,erase
outres,all,all
asel,all
asel,u,area,,1,6               
sfa,all,,conv,70,20  
allsel,all
solve
eplot
esel,s,mat,,1                     !再次选择第一组焊缝单元
nsle,s,all
nsel,r,loc,z,zmin,zmax
nsel,r,loc,y,ymin,ymax
esln,s,all
eplot
nsle,all
ddele,all,temp                    !删除所加的约束(此处后是否还需要进行冷却分析? 就是把上面的瞬态求解改小,然后删除温度后冷却传热?)
alls
j=j+0.5
zmin=zmin+zd
zmax=zmax+zd
!************************
*do,i,2,20                        !循环19次,到第一道焊缝结束
esel,s,mat,,1                     !选择第2组焊缝单元
nsle,s,all
nsel,r,loc,z,zmin,zmax
nsel,r,loc,y,ymin,ymax
esln,s,all
ealive,all                        !激活
eplot
nsle,all
d,all,temp,1500                   !加载1500度
alls
eplot
antype,,rest                      !在上一步基础上分析
time,j                            !时间设成了1S
deltim,0.5,0.5,0.5                !一步求解到位
alls
solve
esel,s,mat,,1                     !选择第二组焊缝单元
nsle,s,all
nsel,r,loc,z,zmin,zmax
nsel,r,loc,y,ymin,ymax
esln,s,all
eplot
nsle,s,all
ddele,all,temp                    !删除温度约束
j=j+0.5
zmin=zmin+zd
zmax=zmax+zd
*enddo                            !循环结束
alls
eplot

到此，焊完后，我想进行5S钟的自然冷却，也就是没有任何加载，只有第20子步的温度场为载荷。

我本以为用
*do,i,1,10
time,j
deltim,0.5,0.5,0.5
alls
solve
j=j+0.5
*enddo

就可以了，但是算出来后不对，结果是10S之前都对，10.5开始就全是20度。

请教如何实现这一过程？  能否给出具体的命令流~~ 谢谢各位大侠啊！！

[ 本帖最后由 bloodyforest 于 2008-8-29 16:30 编辑 ]

bloodyforest

我也用了ANTYPE,,REST  命令来做, 也用过LDREAD,  不知道怎么回事,都不行.

如果哪位高手会的话,麻烦贴下命令流.

LSWRITE的相关命令不能做， 因为有生死功能.

xbzhang

中间挖的那个槽就是焊接的部位吗

xbzhang

这个好像不是钢，请简单叙述一下焊接工艺方法，以及焊接时焊条和工件的位置关系。
看了温度场，最高温度在焊缝的右下，不知对否？
继续试算中，请楼主介绍则个

bloodyforest

我做的是铜合金的焊补， 材料的各项参数并不精确，那个最后更正就好了。 拟采用MIG焊，中间的是个V型槽的一半，算是个铸件的缺陷，  我对此缺陷进行多道焊补。 铜合金是比较难焊的，导热系数比较大，容易产生裂纹，所以我的课题，最后是要模拟焊接裂纹，温度场只是前提。

因为是MIG焊，所以无所谓焊条，我做的第一道焊缝是由右侧焊到左侧，第二道会焊回去。 中间有个冷却时间间隔~~   最高温度不会在右下啊？？  就在焊接的加温度的那些节点，因为用了生死，所以显示时不能直接用PLOTCRTL来显示。

我留个Q， 359786990，  欢迎交流。我可以把所以命令给你看一下。

bloodyforest

我的问题，其实说简单了，问题主要出在，不知道在应用生死功能的前提下，如何把上一步计算得到的温度场做为下一步的载荷加入。

有人说用LSWRITE... 然后最后用 lssolve   但这个肯定不行，不支持生死功能。

有人说用ANTYPE，，REST 就可以了。  但我试过了，即使全在SOLVE前加ANTYPE，，REST， 我的温度场也是错的。问题还是一改变加载的温度，或者不加载任何温度，计算后，都得到与上一步温度场无关的温度场。

xbzhang

先问下：层间温度是多少？现在好像就是一层。最高温度1500℃，最高点不在表面。
我先试算下，逐行看看。

bloodyforest

层间温度原则上是100以内，我还没计算到那么精确，现在就是想把程序先走通，具体的时间后面可以改。  现在是一层，后面的我没贴。

你加我Q，我全给你，你一起走一遍就知道我说的意思了。

bloodyforest

没人帮忙啊~  没人帮忙，我帖子永不沉底~

reltek

楼主你的问题解决了吗？能否分享一下？

ETsy

谢谢，希望从事焊缝研究的同志们一起讨论！QQ 530814511