Abaqus焊接模拟分析程序包括应力场和温度场文档格式.docx
《Abaqus焊接模拟分析程序包括应力场和温度场文档格式.docx》由会员分享,可在线阅读,更多相关《Abaqus焊接模拟分析程序包括应力场和温度场文档格式.docx(26页珍藏版)》请在冰点文库上搜索。
*
材料属性
=====================================================================================!
316L
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7850
MPTEMP,1,20
MPTEMP,2,300
MPTEMP,3,900
MPTEMP,4,1400
MPTEMP,5,2000
MPDATA,KXX,1,,18.6
MPDATA,KXX,1,,21.4
MPDATA,KXX,1,,28.4
MPDATA,KXX,1,,33.9
MPDATA,KXX,1,,48
MPTEMP,2,600
MPTEMP,3,800
MPDATA,C,1,,502
MPDATA,C,1,,612
MPDATA,C,1,,635
MPDATA,C,1,,659
MPDATA,C,1,,670
MPTEMP,2,1300
MPTEMP,3,1410
MPTEMP,4,1440
MPTEMP,5,1550
MPTEMP,6,2000
MPDATA,ENTH,1,,7.88e7
MPDATA,ENTH,1,,6.131e9
MPDATA,ENTH,1,,7.347e9
MPDATA,ENTH,1,,9.145e9
MPDATA,ENTH,1,,1.03e10
MPDATA,ENTH,1,,1.272e10
定义常量
WidthBase=0.025
!
宽度
HeightBase=0.02
基底高度
Length=0.09
长度
WidthClad=0.0015
HeightDeposition=0.00375
覆层高度
Layer=15
层数
HeightClad=HeightDeposition/layer
dt=0.0001
小量
Temp=20
环境温度
InitTemp=300
初始温度
CoffConv=30
对流换热系数
Velocity=0.003
扫描速度
StepDis=0.0015
每个载荷步位移
LaserPower=700
激光功率
Radius=0.0015
激光光斑半径
Area=3.14159265*(Radius**2)
激光光斑面积
Factor=0.3
吸收因子
StepTime=StepDis/Velocity
每个载荷步时间
TotalTime=(Length+Radius*2)/Velocity
载荷持续时间(扫描一层)
StepNum=(Length+Radius*2)/StepDis
载荷步数
建模
BLOCK,0,Length,0,-0.0066,0,WidthClad,
BLOCK,0,Length,0,-0.0066,WidthClad,0.0067
BLOCK,0,Length,-0.0066,-HeightBase,0,WidthClad,
K,,0,0,WidthBase,
K,,Length,0,WidthBase,
K,,0,-HeightBase,WidthBase,
K,,Length,-HeightBase,WidthBase,
V,
16,
13,
27,
25,
15,
14,
28,
26
24,
21,
23,
22,
14
BLOCK,0,Length,0,HeightDeposition,0,WidthClad,
VGLUE,ALL
NUMCMP,ALL
划分网格
LSEL,S,LOC,Y,dt,HeightDeposition-dt,
高度方向
LESIZE,ALL,,,Layer,
LSEL,S,LOC,Y,-DT,-0.0066+DT
LSEL,R,LOC,X,0
LESIZE,ALL,,,4,2
LSEL,R,LOC,X,Length
LESIZE,ALL,,,4,0.5
LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DT
LSEL,R,LOC,Z,0,WidthBase-DT,
LSEL,R,LOC,Z,WidthBase,
LSEL,S,LOC,X,dt,Length-dt,
长度方向
LESIZE,ALL,,,Length/StepDis,
LSEL,S,LOC,Z,dt,WidthClad-dt,
宽度方向
LESIZE,ALL,,,1,
LSEL,S,LOC,Z,WidthClad+DT,0.0067-DT
LESIZE,ALL,,,4,
LSEL,S,LOC,Z,0.0067+DT,WidthBase-DT
LSEL,R,LOC,Y,-DT,-HeightBase
LSEL,R,LOC,Y,0
VSEL,S,LOC,Y,0,HeightDeposition
网格划分
TYPE,1
MAT,1
MSHAPE,0,3D
MSHKEY,1
VMESH,ALL
VSEL,S,LOC,Y,-1,0
ALLSEL,ALL
基底边界条件、初始条件
NSEL,S,LOC,Y,-HeightBase,0
基底初始温度
IC,ALL,TEMP,InitTemp
基底侧面,换热边界条件
NSEL,R,LOC,Z,WidthBase
SF,ALL,CONV,CoffConv,Temp
基底左端面,换热边界条件
NSEL,R,LOC,X,0
基底右端面,换热边界条件
NSEL,R,LOC,X,Length
NSEL,S,LOC,Y,0
基底上表面,换热边界条件
NSEL,R,LOC,Z,WidthClad,WidthBase
FINISH
/SOLU
瞬态分析参数设置
ANTYPE,4
分析类型:
瞬态
TRNOPT,FULL
求解方法:
完全的N-R方法
LNSRCH,on
DELTIM,0.01,0.001,0.05
载荷子步(默认子步时间步长、最小、最大)——载荷步为0.333
NSUBST,4
CNVTOL,HEAT,,0.01,2,0.000001,
收敛准则:
控制热流
OUTRES,NSOL,LAST
结果输出:
所有
杀死单元
NSEL,S,LOC,Y,0,HeightDeposition
杀死熔覆层单元
ESLN,S,1,ALL
EKILL,ALL
ESEL,S,LIVE
EPLOT
ESEL,S,LIVE
激活单元的上表面,指定为对流换热边条
NSLE,S,1
NSEL,R,LOC,Y,0
NSEL,R,LOC,Z,0,Radius
SF,ALL,CONV,CoffConv,Temp
预热
*DO,i,1,2
m=mod(i,2)
*IF,m,EQ,1,THEN
如果为奇数层,向右扫描
*DO,k,1,StepNum,1
TIME,TotalTime*(i-1)+StepTime*k
载荷步结束时间
KBC,1
载荷步内载荷随时间分布:
常数
LeftX=StepDis*(k-1)
RightX=StepDis*k
NSEL,S,LOC,Y,0
ESEL,S,LIVE
ESLN,R,0
NSLE,S,1
NSEL,R,LOC,x,RightX-2*Radius,RightX
NSEL,R,LOC,Z,0,Radius
激活单元的上表面,加热流密度
ESLN,S,1
SFE,ALL,4,HFLUX,,LaserPower*Factor/Area,,,
ALLSEL,ALL
SOLVE
SAVE
SFEDELE,ALL,4,HFLUX
删除热流密度载荷
ESEL,S,LIVE
NSLE,S,1
NSEL,R,LOC,Y,0
NSEL,R,LOC,x,LeftX-2*StepDis,RightX-2*StepDis
NSEL,R,LOC,Z,O,Radius
SF,ALL,CONV,CoffConv,Temp
ALLSEL,ALL
ESEL,S,LIVE
EPLOT
*ENDDO
*ELSE
如果为偶数层,向左扫描
LeftX=Length-StepDis*k
RightX=Length-StepDis*(k-1)
激活单元的上表面,删除载荷
*ENDIF
*ENDDO
熔覆
*DO,i,1,Layer,1
TIME,TotalTime*(i+1)+StepTime*k
NSEL,S,LOC,Y,HeightDeposition/Layer*(i-1),HeightDeposition/Layer*i
NSEL,R,LOC,x,LeftX,RightX
EALIVE,ALL
激活单元的表面,如果包含左端面,指定为对流换热边条
NSEL,R,LOC,X,0
NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1),HeightDeposition/Layer*i
激活单元的表面,如果包含右端面,指定为对流换热边条
NSEL,R,LOC,X,Length
激活单元的侧面,指定为对流换热边条
NSEL,R,LOC,Z,WidthClad