在Abaqus模拟焊接的方法图文操作
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类型，分析步总时间根据现场工艺确定。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
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 瞬时温度场分布。