在Abaqus模拟焊接的方法图文操作
ABAQUS模拟焊接的方法
——子程序dflux方法
作者：yanzongbao110
本文主要通过一个非常简单的平板焊接的例子来介绍abaqus使用子程序dflux模拟焊接的一般方法。
这里采用顺次耦合方法，以下对整个过程做简要介绍。
1. part and mesh
新建model-1,建立如下part（为模型的一半），长，宽，高分别是：100，50，5
图表 1
在model菜单中edit attibute，如下：
图表 2
进入mesh选项卡进行划网格，这里进行粗略划分，网格及单元选择如下图：
图表 3
图表 4
2. material
在material选项卡里新建材料，这里假设焊缝与母材的材料一样，所以只建立一种材料。
材料包括以下部分（本例中没有应用latent heat）：
图表 5
图表 6
3. assembly
在assembly选项中，进行装配。
如下：
图表 7
建立不同的集合及表面，以便后面保用。
建立set-allnode，如下图
图表 8
建立surface-rad，如下图。
注意左侧面为对称面，不选择。
图表 9
4.创建分析步
创建两个热传导分析步，分别是焊缝加热和冷却过程。
如下图：
图表 10
定义第一个分析步：
（a） （b）
图表 11
定义第二个分析步：
(a) (b)
图表 12
5. 创建interaction
在interaction中创建对流和辐射边界条件如下：
(a) (b)
图表 13
(a) (b)
图表 14
6. 边界条件建立
首先建立预定义场，定义初始温度为20℃，选择集合set-allnode，如下图：
(a) (b)
图表 15
建立body heat flux载荷，选择整个实体，设置为user-defined，如下图：
(a) (b)
图表 16
7. 建立名称为weld-thermal的job，并输出.inp。
图表 17
8. 编写dflux子程序
其中热源用的是双椭球热源，具体可参考相关文献。
q＝UIη,这里η代表效率。
x0,y0,z0定义热源的起始点，
a,b,c定义热源的三个半轴，
v定义焊接速度。
本例中，取热源中心坐标为：（0.0，0.0，0.0025）
编写结果如下图所示：
图表 18
9. 提交计算
将产生的inp文件与.for文件放到一个文件夹下，本例中放在D:\weld，在abaqus command窗口进行提交，如下所示：
图表 19
10. 温度场计算结果
图表 20 瞬时HFL分布
图表 21 瞬时温度场分布。
Abaqus 焊缝模拟分析实例
BOOKMARKS
SEARCH
3-Dimensional Repair Weld Simulations – Bead on Plate Comparisons
Dr CD Elcoate1, PJ Bouchard2, and Dr MC Smith2
1) Frazer-Nash Consultancy Limited, 1 Trinity Street, College Green, Bristol, BS1 5TE, UK. 2) British Energy Generation Limited, Barnett Way, Barnwood, Gloucester, GL4 3RS, UK. Full 3-dimensional simulation of multi-pass repair welds is now feasible and practical given the development of improved analysis tools (i.e. ABAQUS) and significantly greater computer power. A single weld bead of limited length, deposited on a flat plate, introduces a strongly varying 3dimensional residual stress field that has similar characteristics to a typical repair weld. This simple welded component is therefore a good vehicle to demonstrate new methods and validate thermal and mechanical predictions. This paper presents a series of results that simulate the deposition of a single weld bead on a flat plate. The work exploits newly developed features within ABAQUS (i.e. plastic strain annealing) which have been demonstrated to be robust and efficient and readily extendible to larger and more complex models. Good comparisons are achieved between simultaneous bead laying and moving heat source approaches whilst it is also noted that the simultaneous bead laying methods offer significant modelling advantages in terms of model creation and execution time.
ABAQUS焊接领域解决方案
Q焊接领域解决方案ABAQUS1国内外焊接变形预测方法固有应变焊接预测法23基于热弹塑性理论的焊接预测法Ab4Abaqus针对焊接行业的解决方案自始至终安心托付1国内外焊接变形预测方法固有应变焊接预测法23基于热弹塑性理论的焊接预测法Ab4Abaqus针对焊接行业的解决方案自始至终安心托付焊接变形预测的重要性焊接是船舶制造最主要的加工手段，焊接水平的高低在很大程度上决定了船体的质量和生产效率，而焊接变形又是焊接过程中最中最难控制的一环。
焊接变形的影响：焊接结构形状变异，尺寸精度下降；承载能力降低；船体在工作载荷作用下引起的附加弯矩和应力集中作用下导致结果失效；船舶结构疲劳降低。
船舶结构疲劳降低因此，对焊接变形的预测及控制已成为船舶生产中迫切需要解决的重要课题。
自始至终安心托付国内外焊接变形的预测方法1.经验（试验）法经验（实验）法是通过试验建立经验公式和数据曲线，用经验公式和数据曲线来估计焊缝的收缩量和角变形量和数据曲线来估计焊缝的收缩量和角变形量。
局限性：z一定条件下的试验或生产实际中得到的，一般被限制在特定的变形模式上变形模式上；z实验受到时间和成本的限制。
真实结构的复杂焊接变形是由多个基本变形组合而成的。
每个基本变形不可能通过有限的实验结果来区分。
自始至终安心托付22.解析法解析法（弹性理论方法）是基于经典弹性理论，忽略热弹塑性的方法。
局限性：由于该方法是建立在平截面假定和其它一些假定的基础上局限性：由于该方法是建立在平截面假定和其它些假定的基础上的，故只能适用于一些焊接是通过熔化金属进行连接的工艺过程。
自始至终安心托付3.数值模拟法焊接数值模拟法又叫焊接计算机仿真，实际上就是热传导有限元解析法和非线性有限元应力解析的组合，已成为线性问题及塑性破坏等非线性问题解析非线性有限元应力解析的组合已成为线性问题及塑性破坏等非线性问题解析不可或缺的手段。
焊接数值模拟是以试验为基础，采用一组控制方程来描述一个焊接过程或个焊接过程的某方面，采用分析或数值方法求解以获得该过程的定量认识。
Abaqus焊接模拟分析程序(包括应力场和温度场)
【我的硕士论文的一部分】求解温度场!上表面上没有对流换热边界条件!单位制:米、秒、摄氏度!/CLEAR,START/FILNAME,temp,0/COM,ANSYS RELEASE 10.0 UP20050718 00:09:52 11/26/2007/CONFIG, NRES, 5000/PREP7/VIEW,1,1,2,3/ANG,1/REP,FAST!*!===============================================================================================!指定单元ET,1,SOLID70!*!*!===============================================================================================!材料属性!=====================================================================================!316LMPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,7850MPTEMP,,,,,,,,MPTEMP,1,20MPTEMP,2,300MPTEMP,3,900MPTEMP,4,1400MPTEMP,5,2000MPDATA,KXX,1,,18.6MPDATA,KXX,1,,21.4MPDATA,KXX,1,,28.4MPDATA,KXX,1,,33.9MPDATA,KXX,1,,48MPTEMP,,,,,,,,MPTEMP,1,20MPTEMP,2,600MPTEMP,3,800MPTEMP,4,1400MPTEMP,5,2000MPDATA,C,1,,502MPDATA,C,1,,612MPDATA,C,1,,635MPDATA,C,1,,659MPDATA,C,1,,670MPTEMP,,,,,,,,MPTEMP,1,20MPTEMP,2,1300MPTEMP,3,1410MPTEMP,4,1440MPTEMP,5,1550MPTEMP,6,2000MPDATA,ENTH,1,,7.88e7MPDATA,ENTH,1,,6.131e9MPDATA,ENTH,1,,7.347e9MPDATA,ENTH,1,,9.145e9MPDATA,ENTH,1,,1.03e10MPDATA,ENTH,1,,1.272e10!===============================================================================================!定义常量WidthBase=0.025 !宽度HeightBase=0.02 !基底高度Length=0.09 !长度WidthClad=0.0015 !宽度HeightDeposition=0.00375 !覆层高度Layer=15 !层数HeightClad=HeightDeposition/layerdt=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.0067BLOCK,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, 26V, 24, 21, 27, 13, 23, 22, 28, 14BLOCK,0,Length,0,HeightDeposition,0,WidthClad,VGLUE,ALLNUMCMP,ALL!===============================================================================================!划分网格LSEL, S, LOC, Y, dt, HeightDeposition-dt, !高度方向LESIZE, ALL, , , Layer,LSEL,S,LOC,Y,-DT,-0.0066+DTLSEL,R,LOC,X,0LESIZE,ALL,,,4,2LSEL,S,LOC,Y,-DT,-0.0066+DTLSEL,R,LOC,X,LengthLESIZE,ALL,,,4,0.5LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DTLSEL,R, LOC, Z, 0, WidthBase-DT,LESIZE,ALL,,,4,2LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DTLSEL,R, LOC, Z, WidthBase,LESIZE,ALL,,,4,0.5LSEL, 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-DTLESIZE,ALL,,,4,LSEL,S,LOC,Z,0.0067+DT,WidthBase-DTLSEL,R,LOC,Y,-DT,-HeightBaseLESIZE,ALL,,,4,2LSEL,S,LOC,Z,0.0067+DT,WidthBase-DTLSEL,R,LOC,Y,0LESIZE,ALL,,,4,0.5VSEL,S,LOC,Y,0,HeightDeposition!网格划分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLVSEL,S,LOC,Y,-1,0 !网格划分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLALLSEL,ALL!============================================================================= ==================!基底边界条件、初始条件NSEL, S, LOC, Y, -HeightBase, 0 !基底初始温度IC,ALL,TEMP,InitTempNSEL, S, LOC, Y, -HeightBase, 0 !基底侧面，换热边界条件NSEL, R, LOC, Z, WidthBaseSF, ALL, CONV, CoffConv, TempALLSEL,ALLNSEL, S, LOC, Y, -HeightBase, 0 !基底左端面，换热边界条件NSEL, R, LOC, X, 0SF, ALL, CONV, CoffConv, TempALLSEL,ALLNSEL, S, LOC, Y, -HeightBase, 0 !基底右端面，换热边界条件NSEL, R, LOC, X, LengthSF, ALL, CONV, CoffConv, TempALLSEL,ALLNSEL, S, LOC, Y, 0 !基底上表面，换热边界条件NSEL, R, LOC, Z, WidthClad, WidthBaseSF, ALL, CONV, CoffConv, TempALLSEL,ALLFINISH/SOLU!===============================================================================================!瞬态分析参数设置ANTYPE,4 !分析类型：瞬态!*TRNOPT,FULL !求解方法：完全的N-R方法!*!DELTIM,0.01,0.001,0.05 !载荷子步(默认子步时间步长、最小、最大)——载荷步为0.333NSUBST, 4CNVTOL,HEAT, ,0.01,2,0.000001, !收敛准则：控制热流OUTRES,NSOL,LAST !结果输出：所有!===============================================================================================!杀死单元NSEL, S, LOC, Y, 0, HeightDeposition !杀死熔覆层单元ESLN, S, 1, ALLEKILL,ALLALLSEL,ALLESEL,S,LIVEEPLOTESEL,S,LIVE !激活单元的上表面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,0NSEL,R,LOC,Z,0,RadiusSF,ALL,CONV,CoffConv,TempALLSEL,ALL!===============================================================================================!预热*DO, i, 1, 2m=mod(i,2)*IF,m,EQ,1,THEN !如果为奇数层，向右扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i-1)+StepTime*k !载荷步结束时间KBC, 1 !载荷步内载荷随时间分布：常数LeftX=StepDis*(k-1)RightX=StepDis*kNSEL, S, LOC, Y, 0ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, RightX-2*Radius, RightXNSEL, R, LOC, Z, 0, Radius !激活单元的上表面，加热流密度ESLN, S, 1SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,ALLSEL, ALLSOLVESAVESFEDELE,ALL,4,HFLUX!删除热流密度载荷ALLSEL, ALLESEL,S,LIVE!激活单元的上表面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,0NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDisNSEL, R, LOC, Z, O, RadiusSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ELSE !如果为偶数层，向左扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i-1)+StepTime*k !载荷步结束时间KBC, 1 !载荷步内载荷随时间分布：常数LeftX=Length-StepDis*kRightX=Length-StepDis*(k-1)NSEL, S, LOC, Y, 0ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, RightX-2*Radius, RightXNSEL, R, LOC, Z, 0, Radius !激活单元的上表面，加热流密度ESLN, S, 1SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,ALLSEL, ALLSOLVESAVESFEDELE,ALL,4,HFLUX!激活单元的上表面，删除载荷ALLSEL, ALLESEL,S,LIVE!激活单元的上表面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,0NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDisNSEL, R, LOC, Z, O, RadiusSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ENDIF*ENDDOESEL,S,LIVEEPLOT!===============================================================================================!熔覆*DO, i, 1, Layer, 1m=mod(i,2)*IF,m,EQ,1,THEN !如果为奇数层，向右扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i+1)+StepTime*k !载荷步结束时间KBC, 1 !载荷步内载荷随时间分布：常数LeftX=StepDis*(k-1)RightX=StepDis*kNSEL, S, LOC, Y, HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALLALLSEL, ALLNSEL, S, LOC, Y, HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iNSEL, R, LOC, x, RightX-2*Radius, RightXNSEL, R, LOC, Z, 0, Radius !激活单元的上表面，加热流密度ESLN, S, 1SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,ALLSEL, ALLESEL,S,LIVE!激活单元的表面，如果包含左端面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,0NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活单元的表面，如果包含右端面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,LengthNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活单元的侧面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Z,WidthCladNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLSOLVESAVESFEDELE,ALL,4,HFLUX!删除热流密度载荷ALLSEL, ALLESEL,S,LIVE!激活单元的上表面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDisSF,ALL,CONV,CoffConv,TempALLSEL,ALLNSEL,S,LOC,Y,HeightDeposition/Layer*(i-1) !激活单元的下表面，删除对流换热边条ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, LeftX-2*StepDis,RightX-2*StepDisNSEL, R, LOC, Z, 0, RadiusNSEL, U, LOC, Y, HeightDeposition/Layer*iESLN, S, 1NSEL, R, LOC, Y,HeightDeposition/Layer*(i-1)SFDELE, ALL, CONVALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ELSE !如果为偶数层，向左扫描*DO, k, 1, StepNum, 1TIME,TotalTime*(i+1)+StepTime*k !载荷步结束时间KBC, 1 !载荷步内载荷随时间分布：常数LeftX=Length-StepDis*kRightX=Length-StepDis*(k-1)NSEL, S, LOC, Y, HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALLALLSEL, ALLNSEL, S, LOC, Y, HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX, LeftX+2*RadiusNSEL, R, LOC, Z, 0, Radius !激活单元的上表面，加热流密度ESLN, S, 1SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , ,ALLSEL, ALLESEL,S,LIVE!激活单元的表面，如果包含左端面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,0NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活单元的表面，如果包含右端面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,X,LengthNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活单元的侧面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Z,WidthCladNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLSOLVESAVESFEDELE,ALL,4,HFLUX!激活单元的上表面，删除载荷ALLSEL, ALLESEL,S,LIVE!激活单元的上表面，指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,HeightDeposition/Layer*iNSEL, R, LOC, x, LeftX+2*StepDis, RightX+2*StepDisSF,ALL,CONV,CoffConv,TempALLSEL,ALLNSEL,S,LOC,Y,HeightDeposition/Layer*(i-1) !激活单元的下表面，删除对流换热边条ESEL, S, LIVEESLN, R, 0NSLE, S, 1NSEL, R, LOC, x, LeftX+2*StepDis,RightX+2*StepDisNSEL, R, LOC, Z, 0, RadiusNSEL, U, LOC, Y, HeightDeposition/Layer*iESLN, S, 1NSEL, R, LOC, Y,HeightDeposition/Layer*(i-1)SFDELE, ALL, CONVALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO*ENDIF*ENDDOESEL,S,LIVEEPLOT!===============================================================================================!冷却!==============================================================================!~100s*DO, k, 1, 2, 1TIME,TotalTime*(Layer+2)+50*k!载荷步结束时间NSUBST, 5KBC, 1SOLVESAVE*ENDDO!==============================================================================!~1000s*DO, k, 1, 9, 1TIME,TotalTime*(Layer+2)+100+100*k!载荷步结束时间NSUBST, 5KBC, 1SOLVESAVE*ENDDO!==============================================================================!~3000s*DO, k, 1, 10, 1TIME,TotalTime*(Layer+2)+1000+200*k!载荷步结束时间NSUBST, 5KBC, 1SOLVESAVE*ENDDO!==============================================================================!~10000s*DO, k, 1, 14, 1TIME,TotalTime*(Layer+2)+3000+500*k!载荷步结束时间NSUBST, 5KBC, 1SOLVESAVE*ENDDO FINISH【我的硕士论文的一部分】求解应力场!修改速度、时间子步步长、载荷文件位置!如果修改基底的热膨胀系数，要修改宏文件!单位制:米、秒、摄氏度/CLEAR,START/FILNAME,stress,0/COM,ANSYS RELEASE 10.0 UP20050718 20:15:52 09/10/2007/CONFIG, NRES, 5000/PREP7/PAGE, 1000, , 1000,/VIEW,1,1,2,3/ANG,1/REP,FAST!*!===============================================================================================!指定单元ET,1,45!*!*!===============================================================================================!材料属性!=====================================================================!316LMPTEMP,,,,,,,,MPTEMP,1,0MPDATA,DENS,1,,7850MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,ALPX,1,,1.75E-005MPTEMP,,,,,,,,MPTEMP,1,20MPTEMP,2,300MPTEMP,3,600MPTEMP,4,900MPTEMP,5,1300 MPDATA,EX,1,,2.0E+11MPDATA,EX,1,,1.7E+11 MPDATA,EX,1,,1.5E+11MPDATA,EX,1,,5.0E+10 MPDATA,EX,1,,0.4E+10MPDATA,PRXY,1,,0.25MPDATA,PRXY,1,,0.25MPDATA,PRXY,1,,0.25MPDATA,PRXY,1,,0.25MPDATA,PRXY,1,,0.25TB,KINH,1,5,4,0TBTEMP,20TBPT,,0,0TBPT,,7E-4,1.4E8TBPT,,0.0012,1.83E8TBPT,,0.1,2.16E9TBTEMP,300TBPT,,0,0TBPT,,5.5E-4,9.35E7TBPT,,0.0012,1.27E8TBPT,,0.1,1.84E9TBTEMP,600TBPT,,0,0TBPT,,3.2E-4,4.8E7TBPT,,0.0012,7.19E7TBPT,,0.1,1.54E9TBTEMP,900TBPT,,0,0TBPT,,2.5E-4,1.25E7TBPT,,0.0012,5.1E7TBPT,,0.1,5.45E8TBTEMP,1300TBPT,,0,0TBPT,,2.5E-4,1E6TBPT,,0.00375,1.13E7TBPT,,0.1,7.05E7!=====================================================================!A3ExpandCoeff=1.75E-005!structural->nonlinear->inelastic->rate independent->kinematic hardeningplasticity->mises plasticity->bilinear!===============================================================================================!定义常量WidthBase=0.025 !宽度HeightBase=0.02 !基底高度Length=0.09 !长度WidthClad=0.0015 !宽度HeightDeposition=0.00375 !覆层高度Layer=15 !层数HeightClad=HeightDeposition/layerdt=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.0067BLOCK,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, 26V, 24, 21, 27, 13, 23, 22, 28, 14BLOCK,0,Length,0,HeightDeposition,0,WidthClad,VGLUE,ALLNUMCMP,ALL!===============================================================================================!划分网格LSEL, S, LOC, Y, dt, HeightDeposition-dt, !高度方向LESIZE, ALL, , , Layer,LSEL,S,LOC,Y,-DT,-0.0066+DTLSEL,R,LOC,X,0LESIZE,ALL,,,4,2LSEL,S,LOC,Y,-DT,-0.0066+DTLSEL,R,LOC,X,LengthLESIZE,ALL,,,4,0.5LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DTLSEL,R, LOC, Z, 0, WidthBase-DT,LESIZE,ALL,,,4,2LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DTLSEL,R, LOC, Z, WidthBase,LESIZE,ALL,,,4,0.5LSEL, 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-DTLESIZE,ALL,,,4,LSEL,S,LOC,Z,0.0067+DT,WidthBase-DTLSEL,R,LOC,Y,-DT,-HeightBaseLESIZE,ALL,,,4,2LSEL,S,LOC,Z,0.0067+DT,WidthBase-DTLSEL,R,LOC,Y,0LESIZE,ALL,,,4,0.5VSEL,S,LOC,Y,0,HeightDeposition!网格划分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLVSEL,S,LOC,Y,-1,0 !网格划分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLALLSEL,ALLVSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DT!删除热影响区外的单元VSEL,A,LOC,Z,0.0067+DT,WidthBase-DTVCLEAR,ALLVDELE,ALL, , ,1ALLSEL,ALLFINISH!===============================================================================================!瞬态分析参数设置/SOLANTYPE,4 !分析类型：瞬态TRNOPT,FULL !求解方法：对于材料非线性，这是唯一的方法NLGEOM,on !大变形分析LNSRCH, onNSUBST, 4NEQIT,30CNVTOL,U,,0.05,2,, !收敛准则：控制位移CNVTOL,F,,0.01,2, !收敛准则：控制力OUTRES,NSOL,LAST !结果输出：TREF, Temp!===============================================================================================!杀死单元NSEL, S, LOC, Y, 0, HeightDeposition !杀死熔覆层单元ESLN, S, 1EKILL,ALLALLSEL,ALLESEL,S,LIVEEPLOTNSEL,S,LOC,Z,0 !对称边条（相当于三个约束，UZ,ROTX,ROTY）D,ALL,UZ,0NSEL,S,LOC,Y,-0.0066 !固定中心点（增加两个约束，UX,UY）NSEL,R,LOC,Z,0NSEL,R,LOC,X,length/2D,ALL,ALL,0NSEL,S,LOC,Y,-0.0066,0 !固定中心线（增加一个约束，ROTZ）NSEL,R,LOC,Z,0NSEL,R,LOC,X,length/2D,ALL,UX,0ALLSEL,ALLSAVE!===============================================================================================!熔覆*DO,m,1,Layerk=mod(m,2)*IF,K,EQ,1,THEN*DO,n,1,StepNum,KBC,0TIME,TotalTime*(m-1)+StepTime*nLeftX=StepDis*(n-1)RightX=StepDis*nNSEL, S, LOC, Y, HeightDeposition/Layer*(m-1),HeightDeposition/Layer*m !熔覆层生长NSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALLALLSEL, ALLLDREAD,TEMP,StepNum*(m+1)+n,last,, ,'temp','rth','F:\temp\differentvelocity\3' !读取体载荷NSEL,S,BF,TEMP,1300,3000BF, ALL, TEMP, 1300ESEL,S,LIVE!显示生单元EPLOTALLSEL,ALLMyDBC!宏命令ALLSEL,ALLSOLVESAVE*ENDDO*ELSEIF,K,EQ,0,THEN*DO,n,1,StepNum,KBC,0TIME,TotalTime*(m-1)+StepTime*nLeftX=Length-StepDis*(n-1)RightX=Length-StepDis*nNSEL, S, LOC, Y, HeightDeposition/Layer*(m-1),HeightDeposition/Layer*mNSEL, R, LOC, x, LeftX, RightXESLN, S, 1EALIVE,ALLALLSEL, ALLLDREAD,TEMP,StepNum*(m+1)+n,last,, ,'temp','rth','F:\temp\differentvelocity\3'NSEL,S,BF,TEMP,1300,3000BF, ALL, TEMP, 1300ESEL,S,LIVE!显示生单元EPLOTALLSEL,ALLMyDBC!宏命令ALLSEL,ALLSOLVESAVE*ENDDO*ENDIF*ENDDO!===============================================================================================!冷却*DO, k, 1, 35, 1TIME,TotalTime*Layer+0.5*k!载荷步结束时间NSUBST, 4KBC, 0 !载荷步内载荷随时间分布：常数LDREAD,TEMP,StepNum*(Layer+2)+K,last,, ,'temp','rth','F:\temp\different velocity\3'NSEL,S,BF,TEMP,1300,3000BF, ALL, TEMP, 1300ESEL,S,LIVE!显示生单元EPLOTALLSEL,ALLMyDBC!宏命令ALLSEL,ALLSOLVESAVE*ENDDOESEL,S,LIVEEPLOTFINISH【补充】其中的宏命令是加比较复杂的位移边条，可以去掉，然后加上你需要的边界条件即可；。
Abaqus软件电子书合集
更新时间：2014-7-31
1.利用Fluent、ABAQUS和MpCCI预处理人体心脏的流体和结构相 互作用（En）.pdf
/Wk_index_fileview_id_5958
2.Abaqus使用技巧总结
/Wk_index_fileview_id_9577
/Wk_index_fileview_id_17650
40.利用ABAQUS explicit对准静态过程涉及复杂的接触在超高 压电缆工程发生的模拟
/Wk_index_fileview_id_18924
41f
/Wk_index_fileview_id_5957
6.Abaqus软件和Mpcci软件文档.ppt
/Wk_index_fileview_id_8529
7.Fluent-ABAQUS联合模拟高压成型（En）.pdf
/Wk_index_fileview_id_18133
42.使用ABAQUS模拟氮化硅熔融金属铸件
/Wk_index_fileview_id_18880
43.Abaqus 操作技巧视频-Abaqus 液体滴水效果.zip
27.abaqus传热与热应激问题3-互动版本.doc
/Wk_index_fileview_id_17641
28.Abaqus 接触问题技巧整理
/Wk_index_fileview_id_18960
29.使用ABAQUS standard预测TIG焊接铝合金失真和残余的应力
48.ABAQUS与FORTRAN论坛问题小结.doc
/Wk_index_fileview_id_18087
49.ABAQUS命令参考
/Wk_index_fileview_id_18784
abaqus多道焊子程序解释
Abaqus是一款广泛应用于工程领域的有限元分析软件，它通过建立虚拟模型，模拟真实的物理过程来解决工程问题。
在焊接领域，Abaqus可以用来模拟多道焊接过程。
本文将对Abaqus中多道焊接程序进行解释与分析。
一、多道焊接概述多道焊接是指在焊缝内部进行多次焊接的过程，通常用于大型结构的焊接，比如船舶、桥梁等。
多道焊接能够降低焊接残余应力和变形，提高焊接质量。
二、Abaqus中多道焊接的模拟1. 模型建立在Abaqus中进行多道焊接的模拟，首先需要建立焊接模型。
这包括定义焊接工艺参数、材料性质、焊接接头几何形状等。
2. 材料本构在多道焊接模拟中，材料的本构关系是非常重要的。
Abaqus提供了多种材料本构模型，可以根据具体的焊接材料选择合适的本构模型。
3. 热分析在多道焊接的过程中，焊缝内部会产生大量的热量。
Abaqus可以通过热分析模块，模拟焊接过程中的温度场分布、热应力等。
4. 结构分析除了热分析，结构分析也是多道焊接模拟中的重要步骤。
通过结构分析，可以得到焊接过程中的应力、变形等信息。
三、多道焊接模拟的应用1. 优化焊接工艺通过Abaqus的多道焊接模拟，可以对焊接工艺进行优化。
比如确定合适的焊接顺序、焊接参数，以达到最佳的焊接质量。
2. 预测残余应力和变形多道焊接模拟可以预测焊接过程中产生的残余应力和变形，帮助工程师更好地设计焊接结构，减少后期修复工作。
3. 焊缝设计通过Abaqus的多道焊接模拟，可以对焊缝进行优化设计，提高焊接连接的强度和稳定性。
四、多道焊接模拟的局限性尽管Abaqus可以进行多道焊接模拟，但也存在一些局限性。
比如对于大型结构、复杂工艺的模拟需要计算资源较大，同时需要对模型的建立和参数选择有一定的专业知识。
总结起来，Abaqus多道焊接程序可以帮助工程师模拟多道焊接过程，优化焊接工艺，预测残余应力和变形，提高焊接结构的质量和可靠性。
然而，工程师在使用Abaqus进行多道焊接模拟时需要注意模型的建立和参数选择，以及计算资源的消耗。
Abaqus操作方法 ppt课件
Abaqus操作方法 ppt课件
9
创建新Part之刚性地面2
➢如左图，画一个 100X100的正方形， 来模拟刚性地面。 ➢点击鼠标中键或 点击按钮 ， 完成。
200
2021/3/26
Abaqus操作方法 ppt课件
10
Part管理器
Part管理器的功能完全可以在窗口左侧模型树的右键快捷菜单实现。
现以Planar为例说明如何Create Wire： 1、选择草绘平面 2、为草绘平面定向 3、草绘Wire
2021/3/26
Abaqus操作方法 ppt课件
14
Create Cut
Create Cut的方法依次为： Extrude、Revolve、Sweep、Loft和Circular hole。
Abaqus操作方法 ppt课件
6
创建Part
如同其他CAE软件，Abaqus的建模功能有限，只适合建立简单Part，如跌落分析中的地面等。
2021/3/26
Part模块专有
特征修改、删除等，很少用到 线、面、体分割工具，辅助网格划分 基准点、线、面及坐标系等
小面修复等，辅助网格划分
备注：如果按钮右下方有小黑三角，左键按住
2021/3/26
Abaqus操作方法 ppt课件
31
网格划分
Mesh Part，即整体划分网格 Mesh Region，即局部划分网格 Delete Part Mesh，即删除整体网格 Delete Region Mesh，即删除局部网 格
2021/3/26
Abaqus操作方法 ppt课件
32
选择已划好网格的Part
4、建立装配体 Assembly
2021/3/26
abaqus 焊点
ABAQUS焊点模拟方法Abaqus中焊点模拟，常用的是fastener方式，单元类型为*CONN3D2。
在Hypermesh 中，生成焊点操作步骤为：（1）进入焊点面板：1D>connectors>spot；（2）将焊点类型切换至fastener；（3）从CATIA中导入焊点文件（制作为焊核几何中心的点）；（4）在num layers中根据焊点的层数，选择tatal 2为两层焊点，total 3为三层焊点；（5）选取焊点后，creat焊点焊点生成后，需要设置焊点属性。
焊点属性分为两个部分，一部分用于指定connector 的属性，即定义*CONN3D2单元的类型；另一部分用于设置焊点的尺寸。
（1）指定*CONN3D2属性为焊点：*CONNECTOR SECTION，ELSET=HMprop_HM_C_1JOIN,ALIGN或者，*CONNECTOR SECTION，ELSET=HMprop_HM_C_1BEAM，（2）指定焊点尺寸：*FASTENER PROPERTY, NAME = HM_P_13.0 ,以上设置完成后，即完成焊点模拟。
在Abaqus中，有多种Connector连接方式。
Join，模拟球铰，只约束三个方向的平动自由度，释放三个转动自由度；Align，仅约束三个转动自由度方向，释放三个平动自由。
因此，Join，Align组合在一起，就可以约束两个点之间的三个平动自由度和三个转动自由度。
而Beam，则能够同时约束两个点的三方向平动自由度和转动自由度。
因此，Beam和Join，Align组合产生的约束效果是一样的。
在Abaqus中，模拟焊点的方法有多种。
1.直接通过刚性单元，如KINCOUP，BEAM等单元来模拟。
2.通过connector的方法来模拟焊点是一个高效的方法。
不考虑模型的计算精度因素，单从建模效率来讲。
采用方法1，需要节点之间对齐，这对于拥有4000~6000焊点的白车身而言，工作量太大了。
ABAQUS-Weld(双枪焊接建模流程)
焊接中的温度变化过程是影响焊接质量的重要因素之一，准确的预测和计算温度的变化是控制其冶金和焊接应力应变变化的前提。
焊接是一个复杂的快速加热、冷却、凝固的过程，它涉及电弧物理、传质传热、冶金和力学的复杂过程。
焊接时各个过程相互影响、相互作用。
热量的传递有三种方式，分别是：热传导、热对流和热辐射。
焊接时，电弧热源是作用在一定的面积上把热能传递给工件的，此面积称为加热斑点，因此高斯热源的热流分布是一种比点热源更加实际的热源分布函数。
焊接热源通过加热斑点把热量传递给焊件，加热斑点上热量分布不均匀，中心热量多，边缘热量少。
ABAQUS 是一套功能强大的工程模拟的有限元软件，其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题。
现采用ABAQUS对双枪焊接过程中温度场分布进行分析。
下面是整个流程1.建模，采用国际标准单位制（m，kg，s，℃）根据现场工艺可知，模型结构为几何对称，边界条件对称，兼顾计算精度及计算时间成本，故建立对称的1/2模型，去焊缝中心所在平面为堆成中心面。
板坯尺寸为500mm（长）×500mm（宽）×10mm（厚），焊缝宽度为8mm。
具体模型参见图1：图1 板坯模型2. 材料属性选用材料为Q235，具体材料参数见下表温度（℃）20 250 500 750 1000 1500 1700 2500导热系数（W/m2℃）50 47 40 27 30 35 140 142 密度（kg/m3）7800 7700 7410 7550 7490 7350 7300 7090比热容（J/kg ℃）460 480 530 675 670 660 780 820具体参数设置如图2所示。
图2 参数设置界面3. 装配，将板坯进行定位，为边界条件施加提供装配模型，对称面为YZ平面，焊道与Z轴平行，焊点沿Z轴移动，如图3所示。
图3 板坯组装图4.分析步设定，本为主要分析双枪焊接过程中，温度场变化情况，故选择分析类型为Heat transfer，根据实际工艺，选择Transient类型，分析步总时间根据现场工艺确定。
基于Marc软件的钢板焊接数值模拟
基于Marc软件的钢板焊接数值模拟◎ 闫玉伯1 杨迪2 姜炜珉3 桂劲松31.曹妃甸新天液化天然气有限公司；2.中交天津港湾工程设计院有限公司；3.大连海洋大学海洋与土木工程学院通讯作者：姜炜珉摘 要：随着计算机技术的快速发展，焊接变形的数值模拟在研究和设计领域得到广泛应用。
本文运用有限元分析软件Marc进行钢板焊接的数值模拟分析，并进行了试验验证，确定了Marc软件进行焊接数值模拟分析的准确性。
利用Marc软件进行焊接的数值模拟步骤简单，焊接参数设置齐全。
模拟结论得出，在符合焊接尺寸标准的一定厚度范围内，焊接构件厚度越大，焊后构件上残余应力越大。
关键词：钢板焊接；数值模拟；残余应力1.引言结构数值模拟方法是与计算机技术一同兴起的一种模拟分析方法，这种方法可以弥补实验条件受限、实验成本高以及复杂模型难以进行理论分析的问题，能够比较好地模拟单元节点在不同材料、荷载、边界条件、性能等环境下的受力情况。
通过对应力和应变计算结果的分析，能够对单元节点处的受力状态、变形程度有一个全面的了解和把握，可以针对性地调整受力薄弱区的边界条件、不同焊接坡口形式、焊接顺序等影响因素，相应调整都可以方便地在软件中实现[1]。
焊接数值模拟通过构建目标数学模型，对其施加初始条件和边界条件，求解相应的微分方程组来解决焊接温度场、应力和变形等问题，并将分析得到的结果通过计算机直观地表达出来，研究人员能够通过模拟焊接过程对焊件进行检验，并对工件结构形式以及焊接工艺参数进行优化。
当前，进行焊接数值模拟的主流软件有Marc、Ansys、Abaqus、Jwrian、Sysweld等，其中Marc软件擅长非线性条件下的大变形计算分析，可以利用其便捷的“生死单元法”直接对焊接过程展开热力耦合计算，其计算精度和计算效率均较好，分析钢板焊接残余应力较为有利。
本文旨在探究用Marc软件进行钢板焊接数值模拟的便捷性以及分析数据的准确性，以验证其模拟焊接过程、分析残余应力的实用性。
abaqus 实例
焊接ABAQUS建模实例一、网格的建立或导入在定义材料属性时，往往要根据单元的集合来定义，这些集合首先需要在part中定义。
当采用hypermesh来划分网格时，在hypermesh中定义的SET只存在于assembly中，而在part中没有，因此需要重新在part中定义这些集合，通常情况下可以通过修改inp文件来实现，具体操作步骤如下：1.导入hypermesh所生成的inp文件，具体操作为：File——import——models，出现Import models 对话框，File filter选择*.inp，然后选择所要导入的inp文件。
2. 建立新的job，重新生成inp文件，具体操作如下：在analysis模型树中，右击job——create，输入新的job名称，例如weld，点击continue，然后点击ok。
单击模型树job前的+，会出现刚才生成的weld，右键job，选择write input，就会在文件夹中生成新的inp文件。
3. 选用记事本或者写字板打开刚才生成的weld.inp文件，选择你需定义的材料单元集合，将它们复制到*end part之前，并删除后面的instance=part-1-1，保存。
4.重新用abaqus打开刚才修改后的inp文件，就可以进行后续的建模和分析了。
二、材料的定义材料的定义分为两个过程，首先定义材料性能参数，然后将这些材料性能参数赋予给具体的单元。
1.材料参数的定义。
右击material——create，出现edit material对话框，输入材料的名称，例如weld、base 、haz。
焊接过程需要定义的材料参数主要有density、elastic、plastic、expansion、conductivity、specific heat 、latent heat。
下面以elasticity为例作一个简单的介绍，单击对话框的mechanical——elasticity——elastic，会出现如下的对话框如果材料性能是温度的函数，把use temperature-dependent data 前的框选上。