Abaqus焊接模拟分析程序(包括应力场和温度场)

1、【我的硕士论文的一部分】求解温度场囹登囹！上表面上没有对流换热边界条件!单位制：米、秒、摄氏度!/CLEAR,START/FILNAME,temp,0/COM,ANSYSRELEASE10.0UP2005071800:09:5211/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,14

2、00MPTEMP,5,2000MPDATA,KXX,1,18.6MPDATA,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,20

3、00MPDATA,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!环境

4、温度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!载荷步数!=!建模BLOC

5、K,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,26BLOCK,0,Length,0,HeightDeposition,0,WidthClad,VGLUE,ALLNUMCM

6、P,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,乙0,WidthBase-DT,LESIZE,ALL,4,2LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+D

7、TLSEL,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-

8、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,InitTempALLSEL,ALLNSEL,S,LOC,Y,-HeightBase,0!基底侧面，换热边界条件NSEL,R,LO

9、C,Z,WidthBaseSF,ALL,CONV,CoffConv,TempALLSEL,ALL基底左端面，换热边界条NSEL,S,LOC,Y,-HeightBase,0NSEL,R,LOC,X,0SF,ALL,CONV,CoffConv,TempALLSEL,ALLNSEL,S,LOC,Y,-HeightBase,0基底右端面，换热边界条NSEL,S,LOC,Y,0NSEL,R,LOC,X,LengthSF,ALL,CONV,CoffConv,TempALLSEL,ALL基底上表面，换热边界条件NSEL,R,LOC,Z,WidthClad,WidthBaseSF,ALL,CONV,CoffC

10、onv,TempALLSEL,ALLFINISH/SOLU=!瞬态分析参数设置ANTYPE,4TRNOPT,FULL!DELTIM,0.01,0.001,0.05!分析类型：瞬态!*!求解方法：完全的N-R方法LNSRCH,on!*!载荷子步（默认子步时间步长、最小、最大）一一载荷步为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

11、,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,

12、LIVEESLN,R,0NSLE,S,1NSEL,R,LOC,x,RightX-2DO,k,1,StepNum,1TIME,TotalTime*(i-1)+StepTime*k!载荷步结束时间载荷步内载荷随时间分布:Radius,RightXNSEL,R,LOC,乙0,Radius!激活单元的上表面，加热流密度ESLN,S,1SFE,ALL,4,HFLUX,LaserPower*Factor/Area,ALLSEL,ALLSOLVESAVESFEDELE,ALL,4,HFLUX!删除热流密度载荷ALLSEL,ALLESEL,S,LIVE!激活单元的上表面，指定为对流换热边条NSLE,S,1NS

13、EL,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*ELSE*ENDDO如果为偶数层，向左扫描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,Ra

14、dius!热流密度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,ALL激活单元的上表面，加激活单元的上表面,激活单元的上表面,ESEL,S,LIVEEPLOT*EN

15、DDO*ENDIF*ENDDOESEL,S,LIVEEPLOT!熔覆*IF,m,EQ,1,THEN*DO,i,1,Layer,1m=mod(i,2)如果为奇数层，向右扫描*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,

16、ALLALLSEL,ALLNSEL,S,LOC,Y,HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iNSEL,R,LOC,x,RightX-2*Radius,RightX激活单元的上表面，加NSEL,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

17、*(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

18、,Y,HeightDeposition/Layer*(i-1),HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLSOLVESAVESFEDELE,ALL,4,HFLUX!删除热流密度载荷ESEL,S,LIVEALLSEL,ALL激活单元的上表面,指定为对流换热边条NSLE,S,1NSEL,R,LOC,Y,HeightDeposition/Layer*iNSEL,R,LOC,x,LeftX-2*StepDis,RightX-2*StepDisSF,ALL,CONV,CoffConv,TempALLSEL,ALLNSEL,S,L

19、OC,Y,HeightDeposition/Layer*(i-1)!激活单元的下表面，删除对流换热边条ESEL,S,LIVEESLN,R,0NSLE,S,1NSEL,R,LOC,x,LeftX-2*StepDis,RightX-2*StepDisNSEL,R,LOC,乙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

20、,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

21、*iNSEL,R,LOC,x,LeftX,LeftX+2*RadiusNSEL,R,LOC,乙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!激

22、活单元的表面，如果包含右端面，指定为对流换热边条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,Te

23、mpALLSEL,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,0NSL

24、E,S,1NSEL,R,LOC,x,LeftX+2*StepDis,RightX+2*StepDisNSEL,R,LOC,乙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,TotalTimeENDDO!=!3000s(Layer+2)+50*k!载荷步结

25、束时间NSUBST,5KBC,1SOLVESAVE*ENDDO!=!-1000s*DO,k,1,9,1TIME,TotalTime*(Layer+2)+100+100*k!载荷步结束时间NSUBST,5KBC,1SOLVESAVETIME,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*ENDDOFINISH【我的硕士论文的一部分】求解应力场史

26、江幻!修改速度、时间子步步长、载荷文件位置!如果修改基底的热膨胀系数，要修改宏文件!单位制：米、秒、摄氏度/CLEAR,START/FILNAME,stress,0/COM,ANSYSRELEASE10.0UP2005071820:15:5209/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.75

27、E-005MPTEMP,MPTEMP,1,20MPTEMP,2,300MPTEMP,3,600MPTEMP,4,900MPTEMP,5,1300MPDATA,EX,1,2.0E+11MPDATA,EX,1,1.7E+11MPDATA,EX,1,1.5E+11MPDATA,EX,1,5.0E+10MPDATA,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

28、.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,2.5E-4,1E6TBPT,0.00375,1.13E7TBPT,0.1,7.05

29、E7=!A3ExpandCoeff=1.75E-005!structural->nonlinear->inelastic->rateindependent->kinematichardeningplasticity->misesplasticity->bilinear=!定义常量WidthBase=0.025!宽度HeightBase=0.02!基底高度Length=0.09!长度WidthClad=0.0015!宽度HeightDeposition=0.00375!覆层高度Layer=15!层数HeightClad=HeightDeposition/lay

30、erdt=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

31、*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,Len

32、gth,0,HeightDeposition,0,WidthClad,VGLUE,ALLNUMCMP,ALL!=!划分网格LSEL,S,LOC,Y,dt,HeightDeposition-dt,!高度方向LESIZE,ALL,Layer,LSEL,S,LOC,Y,-DT,-0.0066+DTLESIZE,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,乙0,WidthBase-DT,LESIZE,ALL,4,

33、2LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DTLSEL,R,LOC,乙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,A

34、LL,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-D

35、TVCLEAR,ALLVDELE,ALL,1ALLSEL,ALLFINISH=!瞬态分析参数设置ANTYPE,4/SOL!分析类型：瞬态TRNOPT,FULL!求解方法:的方法:对于材料非线性，这是唯一NLGEOM,on!LNSRCH,on大变形分析CNVTOL,U,0.05,2,NSUBST,4NEQIT,30!收敛准则：控制位移CNVTOL,F,0.01,2,收敛准则：控制力OUTRES,NSOL,LAST结果输出:TREF,Temp!杀死单元杀死熔覆层单元NSEL,S,LOC,Y,0,HeightDepositionESLN,S,1EKILL,ALLALLSEL,ALLESEL,S,LIVEEPLOTNSEL,S,LOC,Z,0对称边条（相当于三个约束，UZ,ROTX,ROTYD,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固定中心线（增加一个约束，ROTZNSEL,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,K