- 积分
- 0
- 注册时间
- 2008-8-25
- 仿真币
-
- 最后登录
- 1970-1-1
|
多道焊中,焊接一道焊缝后的自然冷却问题~请求高人指点。
我做一个焊件的焊补模拟, 采用多道焊, 我想把第一道焊缝焊完后的自然冷却过程考虑进去,不知道为什么我的命令算出来有问题,请各位高手指点。我初学,命令流比较简单,希望大家有空可以帮我调试一下。
命令流我贴出来, 下面是焊完第一道焊缝:
!******************************第一部分进行温度场分析******************************
/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 编辑 ] |
|