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

【我的硕士论文的一部分】求解温度场!上表面上没有对流换热边界条件!单位制:米、秒、摄氏度!/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,7850 MPTEMP,MPTEMP,1,20 MPTEMP,2,300MPTEMP,3,900MPTEMP,4,1400 MPTEMP,5,2000 MPDATA,KXX,1,18.6MPDATA,KXX,1,21.4MPDATA,KXX,1,28.4MPDATA,KXX,1,33.9MPDATA,KXX,1,48MPTEMP,MPTEMP,1,20 MPTEMP,2,600MPTEMP,3,800MPTEMP,4,1400 MPTEMP,5,2000 MPDATA,C,1,502 MPDATA,C,1,612 MPDATA,C,1,635 MPDATA,C,1,659 MPDATA,C,1,670 MPTEMP,MPTEMP,1,20 MPTEMP,2,1300 MPTEMP,3,1410 MPTEMP,4,1440 MPTEMP,5,1550 MPTEMP,6,2000 MPDATA,ENTH,1,7.88e7 MPDATA,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, 26 V, 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+DT LSEL,R,LOC,X,0 LESIZE,ALL,4,2LSEL,S,LOC,Y,-DT,-0.0066+DT LSEL,R,LOC,X,LengthLESIZE,ALL,4,0.5LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DT LSEL,R, LOC, Z, 0, WidthBase-DT, LESIZE,ALL,4,2LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DT LSEL,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-DT LESIZE,ALL,4,LSEL,S,LOC,Z,0.0067+DT,WidthBase-DTLSEL,R,LOC,Y,-DT,-HeightBase LESIZE,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,3D MSHKEY,1VMESH,ALLVSEL,S,LOC,Y,-1,0 !网格划分TYPE,1MAT,1MSHAPE,0,3D MSHKEY,1VMESH,ALLALLSEL,ALL!=!基底边界条件、初始条件NSEL, S, LOC, Y, -HeightBase, 0 !基底初始温度IC,ALL,TEMP,InitTempALLSEL,ALLNSEL, 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方法LNSRCH, on!*!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, 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 !删除热流密度载荷 ALLSEL, ALL 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 !如果为偶数层，向左扫描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i-1)+StepTime*k !载荷步结束时间 KBC, 1 !载荷步内载荷随时间分布：常数 LeftX=Length-StepDis*k RightX=Length-StepDis*(k-1) 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 !激活单元的上表面，删除载荷 ALLSEL, ALL 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 *ENDIF*ENDDO ESEL,S,LIVEEPLOT!=!熔覆*DO, i, 1, Layer, 1m=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, HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i NSEL, R, LOC, x, LeftX, RightX ESLN, S, 1 EALIVE,ALL ALLSEL, ALL NSEL, S, LOC, Y, HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i 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 ESEL,S,LIVE !激活单元的表面，如果包含左端面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,X,0 NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i SF,ALL,CONV,CoffConv,Temp ALLSEL,ALL ESEL,S,LIVE !激活单元的表面，如果包含右端面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,X,Length NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i SF,ALL,CONV,CoffConv,Temp ALLSEL,ALL ESEL,S,LIVE !激活单元的侧面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,Z,WidthClad NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i SF,ALL,CONV,CoffConv,Temp ALLSEL,ALL SOLVE SAVE SFEDELE,ALL,4,HFLUX !删除热流密度载荷 ALLSEL, ALL ESEL,S,LIVE !激活单元的上表面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,Y,HeightDeposition/Layer*i NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDis SF,ALL,CONV,CoffConv,Temp ALLSEL,ALL NSEL,S,LOC,Y,HeightDeposition/Layer*(i-1) !激活单元的下表面，删除对流换热边条 ESEL, S, LIVE ESLN, R, 0 NSLE, S, 1 NSEL, R, LOC, x, LeftX-2*StepDis,RightX-2*StepDis NSEL, R, LOC, Z, 0, Radius NSEL, U, LOC, Y, HeightDeposition/Layer*i ESLN, S, 1 NSEL, R, LOC, Y,HeightDeposition/Layer*(i-1) SFDELE, ALL, CONV ALLSEL,ALL ESEL,S,LIVE EPLOT *ENDDO *ELSE !如果为偶数层，向左扫描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i+1)+StepTime*k !载荷步结束时间 KBC, 1 !载荷步内载荷随时间分布：常数 LeftX=Length-StepDis*k RightX=Length-StepDis*(k-1) NSEL, S, LOC, Y, HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i NSEL, R, LOC, x, LeftX, RightX ESLN, S, 1 EALIVE,ALL ALLSEL, ALL NSEL, S, LOC, Y, HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i NSEL, R, LOC, x, LeftX, LeftX+2*Radius NSEL, R, LOC, Z, 0, Radius !激活单元的上表面，加热流密度 ESLN, S, 1 SFE, ALL, 4, HFLUX, , LaserPower*Factor/Area, , , ALLSEL, ALL ESEL,S,LIVE !激活单元的表面，如果包含左端面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,X,0 NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i SF,ALL,CONV,CoffConv,Temp ALLSEL,ALL ESEL,S,LIVE !激活单元的表面，如果包含右端面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,X,Length NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i SF,ALL,CONV,CoffConv,Temp ALLSEL,ALL ESEL,S,LIVE !激活单元的侧面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,Z,WidthClad NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*i SF,ALL,CONV,CoffConv,Temp ALLSEL,ALL SOLVE SAVE SFEDELE,ALL,4,HFLUX !激活单元的上表面，删除载荷 ALLSEL, ALL ESEL,S,LIVE !激活单元的上表面，指定为对流换热边条 NSLE,S,1 NSEL,R,LOC,Y,HeightDeposition/Layer*i NSEL, R, LOC, x, LeftX+2*StepDis, RightX+2*StepDis SF,ALL,CONV,CoffConv,Temp ALLSEL,ALL NSEL,S,LOC,Y,HeightDeposition/Layer*(i-1) !激活单元的下表面，删除对流换热边条 ESEL, S, LIVE ESLN, R, 0 NSLE, S, 1 NSEL, R, LOC, x, LeftX+2*StepDis,RightX+2*StepDis NSEL, R, LOC, Z, 0, Radius NSEL, U, LOC, Y, HeightDeposition/Layer*i ESLN, S, 1 NSEL, R, LOC, Y,HeightDeposition/Layer*(i-1) SFDELE, ALL, CONV ALLSEL,ALL ESEL,S,LIVE EPLOT *ENDDO *ENDIF*ENDDO ESEL,S,LIVEEPLOT!=!冷却!=!100s *DO, k, 1, 2, 1 TIME,TotalTime*(Layer+2)+50*k !载荷步结束时间 NSUBST, 5 KBC, 1 SOLVE SAVE *ENDDO!=!1000s *DO, k, 1, 9, 1 TIME,TotalTime*(Layer+2)+100+100*k !载荷步结束时间 NSUBST, 5 KBC, 1 SOLVE SAVE *ENDDO!=!3000s *DO, k, 1, 10, 1 TIME,TotalTime*(Layer+2)+1000+200*k !载荷步结束时间 NSUBST, 5 KBC, 1 SOLVE SAVE *ENDDO!=!10000s *DO, k, 1, 14, 1 TIME,TotalTime*(Layer+2)+3000+500*k !载荷步结束时间 NSUBST, 5 KBC, 1 SOLVE SAVE *ENDDOFINISH【我的硕士论文的一部分】求解应力场!修改速度、时间子步步长、载荷文件位置!如果修改基底的热膨胀系数，要修改宏文件!单位制:米、秒、摄氏度/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,0 MPDATA,ALPX,1,1.75E-005MPTEMP,MPTEMP,1,20 MPTEMP,2,300MPTEMP,3,600MPTEMP,4,900MPTEMP,5,1300MPDATA,EX,1,2.0E+11 MPDATA,EX,1,1.7E+11 MPDATA,EX,1,1.5E+11 MPDATA,EX,1,5.0E+10MPDATA,EX,1,0.4E+10 MPDATA,PRXY,1,0.25 MPDATA,PRXY,1,0.25 MPDATA,PRXY,1,0.25 MPDATA,PRXY,1,0.25 MPDATA,PRXY,1,0.25 TB,KINH,1,5,4,0 TBTEMP,20 TBPT,0,0 TBPT,7E-4,1.4E8TBPT,0.0012,1.83E8 TBPT,0.1,2.16E9TBTEMP,300TBPT,0,0 TBPT,5.5E-4,9.35E7TBPT,0.0012,1.27E8 TBPT,0.1,1.84E9TBTEMP,600TBPT,0,0 TBPT,3.2E-4,4.8E7 TBPT,0.0012,7.19E7 TBPT,0.1,1.54E9TBTEMP,900TBPT,0,0 TBPT,2.5E-4,1.25E7 TBPT,0.0012,5.1E7TBPT,0.1,5.45E8TBTEMP,1300TBPT,0,0 TBPT,2.5E-4,1E6 TBPT,0.00375,1.13E7TBPT,0.1,7.05E7!=!A3ExpandCoeff=1.75E-005!structural-nonlinear-inelastic-rate independent-kinematic hardening plasticity-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, 26 V, 24, 21, 27, 13, 23, 22, 28, 14BLOCK,0,Length,0,HeightDeposition,0,WidthClad,VGLUE,ALLNUMCMP,ALL!=!划分网格LS