焊接残余应力模拟中的问题

kunull

焊接残余应力模拟中为什么计算应力会远大于材料的屈服应力
下面是我做的焊接残余应力模拟（参考网上一个比较流行的例题），为什么最后的焊接残余应力会远大于材料的屈服强度，计算应力为654MPa，屈服强度仅为202MPa，在计算中我设置了eresx,no                !不选择积分点外推
还有，如果Qmax=Qm*kk*0.25    按道理是不用乘以0.25的，但是不乘的话温度场计算不下去。。
高手们给我指教指教啊！小弟先谢过了！  
下面是我的命令流和结果图

/filname,Q235
/title,Q235  WELDING 2-D FEA
/units,si
/prep7
et,1,plane55
et,2,surf151
keyopt,2,8,4
mptemp,1,20,250,500,750,1000,1500
mptemp,7,1700,2500
mpdata,kxx,1,1,50,47,40,27,30,35
mpdata,kxx,1,7,140,142
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350
mpdata,dens,1,7,7300,7090
mpdata,c,1,1,460,480,530,675,670,660
mpdata,c,1,7,780,820
mpdata,prxy,1,1,0.28,0.29,0.31,0.35,0.4,0.45
mpdata,prxy,1,7,0.48,0.5
mpdata,alpx,1,1,1.10e-5,1.22e-5,1.39e-5,1.48e-5,1.34e-5,1.33e-5
mpdata,alpx,1,7,1.32e-5,1.31e-5
mpdata,ex,1,1,2.05e11,1.87e11,1.50e11,7.0e10,2.0e10,1.9e10
mpdata,ex,1,7,1.8e10,1.2e10
tb,biso,1,6,2
tbtemp,20
tbdata,,2.20e8,2.05e10
tbtemp,500
tbdata,,1.3e8,1.50e10
tbtemp,750
tbdata,,4.0e7,7.0e9
tbtemp,1000
tbdata,,2.50e7,2.0e9,
tbtemp,1700
tbdata,,1.0e6,1.8e9
tbtemp,2500
tbdata,,1.0e5,1.2e9,
pi=3.1415926
I=100 !焊接电流
U=23 !焊接电弧电压
effcient=0.85 !焊接热效率η
r=0.005 !电弧有效加热半径
Qm=3*I*U*effcient/(pi*r**2) !最大热流高斯密度
kk=0.85 !能量集中系数
Qmax=Qm*kk*0.25      
ly=0.06
lx=0.3
lsize=0.005 !划分单元网格长度
v=0.004 !焊接速度
tm=lx/v !焊接总时间
tinc=lsize/v !单步时间长度
max_time=1+tm/tinc
max_x=1+lx/lsize
max_y=1+ly/lsize
rectng,0,lx,0,ly   
lesize,all,LSIZE
mshape,0,2d
mshkey,1
amesh,all
nsll,s
type,2
esurf
allsel,all
save
finish
!以下为输入高斯移动热源
*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
distance=sqrt(((i-1)*lsize-xcenter)**2+((j-1)*lsize-ycenter)**2)
    *if,distance,le,r,then
        flux2(i,j,k)=Qmax/exp(3*distance**2/(r**2))
    *else
        flux2(i,j,k)=0
    *endif
    *enddo
  *enddo
flux2(0,0,k)=(k-1)*tinc
*enddo
finish
/solu
antype,trans
tintp,0.005,,,1,,,   !设置向后积分
kbc,1
timint,off
autots,on
ESEL,S,TYPE,,2
SFE,ALL,,CONV,,100
SFE,ALL,,CONV,2,20    !表面效应单元施加对流系数
ALLSEL,ALL
tref,20
time,1e-6
d,all,temp,20
NSUBST,1
solve               !进行稳态求解，得到初始条件
ddele,all,temp
timint,on
outres,all,all
tsres,erase
kbc,0
esel,s,type,,1
nsle,s
sf,all,hflux,%flux2%
allsel,all
time,73.75
autots,off
deltim,tinc,tinc,tinc
solve
sfdele,all,hflux      
time,75            
solve
time,725                !冷却过程到725s
autots,0
deltim,25,25,25
solve
finish
/prep7
etchg,tts                !转换单元，进行结构分析
nsel,s,loc,x,0
nsel,r,loc,y,0
d,all,all
allsel,all
dl,1,,symm
finish
/solu
eresx,no                !不选择积分点外推
tref,20
time,1e-6
LDREAD,TEMP,1,1,, ,'Q235','rth',' '
solve
*do,t,1.25,73.75,1.25
time,t
LDREAD,TEMP,,,t, ,'Q235','rth',' '
solve
*enddo  
*do,i,1,26,1         
time,50+i*25
LDREAD,TEMP,4,i,, ,'Q235','rth',' '  
solve
*enddo
finish

jesus0341

二维焊接温度场能这样做吗？这样做相当于是在无限薄板上加移动热源，不能代表实际的焊接过程。我觉得要简化成二维，模型应该是三维的横向切面做。

kunull

楼上说的有道理，但是没有解决我的问题！
我这样做正是尝试无限大板加移动热源，和实际焊接过程并不符合，因为这是从网上的一个例题改过来的，我就是想试试这样可不可行。

但是结果就是出现了计算应力远大于材料屈服强度的现象，一直很困惑啊！谁能解释一下啊

xujunchen100

，感觉焊接残余应力的模拟好难啊！

无言

我发现越学越难了

暴戾新域

同感，我也一样觉得沉不进去

平淡

是啊！能尘的进去吗？

蓝色风沙

焊接残余应力啊，我服了你啦

XYh151969

我觉得应该是 输入热流量太大的问题  LZ减小电流 电压试试   不要考虑相变等因素  还有就是单位可能出了问题

lu286598468

解决问题了吗楼主 能不能把例题发一下学习学习