Abaqus焊接模擬分析程序(包括應(yīng)力場和溫度場).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 !層數(shù)HeightClad=HeightDeposition/layerdt=0.0001 !小量Temp=20 !環(huán)境溫度InitTemp=300 !初始溫度CoffConv=30 !對流換熱系數(shù)!=!定義常量Velocity=0.003 !掃描速度StepDis=0.0015 !每個(gè)載荷步位移LaserPower=700 !激光功率Radius=0.0015 !激光光斑半徑Area=3.14159265*(Radius*2) !激光光斑面積Factor=0.3 !吸收因子StepTime=StepDis/Velocity !每個(gè)載荷步時(shí)間TotalTime=(Length+Radius*2)/Velocity !載荷持續(xù)時(shí)間(掃描一層)StepNum=(Length+Radius*2)/StepDis !載荷步數(shù)!=!建模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!=!劃分網(wǎng)格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 !網(wǎng)格劃分TYPE,1MAT,1MSHAPE,0,3D MSHKEY,1VMESH,ALLVSEL,S,LOC,Y,-1,0 !網(wǎng)格劃分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 !基底側(cè)面，換熱邊界條件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!=!瞬態(tài)分析參數(shù)設(shè)置ANTYPE,4 !分析類型：瞬態(tài)!*TRNOPT,FULL !求解方法：完全的N-R方法LNSRCH, on!*!DELTIM,0.01,0.001,0.05 !載荷子步(默認(rèn)子步時(shí)間步長、最小、最大)載荷步為0.333NSUBST, 4CNVTOL,HEAT, ,0.01,2,0.000001, !收斂準(zhǔn)則：控制熱流OUTRES,NSOL,LAST !結(jié)果輸出：所有!=!殺死單元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!=!預(yù)熱*DO, i, 1, 2m=mod(i,2)*IF,m,EQ,1,THEN !如果為奇數(shù)層，向右掃描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i-1)+StepTime*k !載荷步結(jié)束時(shí)間 KBC, 1 !載荷步內(nèi)載荷隨時(shí)間分布：常數(shù) 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 !如果為偶數(shù)層，向左掃描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i-1)+StepTime*k !載荷步結(jié)束時(shí)間 KBC, 1 !載荷步內(nèi)載荷隨時(shí)間分布：常數(shù) 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 !如果為奇數(shù)層，向右掃描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i+1)+StepTime*k !載荷步結(jié)束時(shí)間 KBC, 1 !載荷步內(nèi)載荷隨時(shí)間分布：常數(shù) 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 !激活單元的側(cè)面，指定為對流換熱邊條 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 !如果為偶數(shù)層，向左掃描 *DO, k, 1, StepNum, 1 TIME,TotalTime*(i+1)+StepTime*k !載荷步結(jié)束時(shí)間 KBC, 1 !載荷步內(nèi)載荷隨時(shí)間分布：常數(shù) 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 !激活單元的側(cè)面，指定為對流換熱邊條 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 !載荷步結(jié)束時(shí)間 NSUBST, 5 KBC, 1 SOLVE SAVE *ENDDO!=!1000s *DO, k, 1, 9, 1 TIME,TotalTime*(Layer+2)+100+100*k !載荷步結(jié)束時(shí)間 NSUBST, 5 KBC, 1 SOLVE SAVE *ENDDO!=!3000s *DO, k, 1, 10, 1 TIME,TotalTime*(Layer+2)+1000+200*k !載荷步結(jié)束時(shí)間 NSUBST, 5 KBC, 1 SOLVE SAVE *ENDDO!=!10000s *DO, k, 1, 14, 1 TIME,TotalTime*(Layer+2)+3000+500*k !載荷步結(jié)束時(shí)間 NSUBST, 5 KBC, 1 SOLVE SAVE *ENDDOFINISH【我的碩士論文的一部分】求解應(yīng)力場!修改速度、時(shí)間子步步長、載荷文件位置!如果修改基底的熱膨脹系數(shù)，要修改宏文件!單位制:米、秒、攝氏度/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 !層數(shù)HeightClad=HeightDeposition/layerdt=0.0001 !小量Temp=20 !環(huán)境溫度InitTemp=300 !初始溫度CoffConv=30 !對流換熱系數(shù)!=!定義常量Velocity=0.003 !掃描速度StepDis=0.0015 !每個(gè)載荷步位移LaserPower=700 !激光功率Radius=0.0015 !激光光斑半徑Area=3.14159265*(Radius*2) !激光光斑面積Factor=0.3 !吸收因子StepTime=StepDis/Velocity !每個(gè)載荷步時(shí)間TotalTime=(Length+Radius*2)/Velocity !載荷持續(xù)時(shí)間(掃描一層)StepNum=(Length+Radius*2)/StepDis !載荷步數(shù)!=!建模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!=!劃分網(wǎng)格LS