Abaqus焊接模擬分析程序

1、【我的碩士論文的一部分】求解溫度場(chǎng)!上表面上沒(méi)有對(duì)流換熱邊界條件!單位制:米、秒、攝氏度!/CLEAR,START/FILNAME,temp,0/COM,ANSYS RELEASE 10.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,1

2、400MPTEMP,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,2

3、000MPDATA,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!長(zhǎng)度WidthClad=0.0015!寬度HeightDeposition=0.00375!覆層高度Layer=15!層數(shù)dt=0.0001Temp=20InitTemp=300CoffConv=30Velocity=0.003St

4、epDis=0.0015LaserPower=700Radius=0.0015HeightClad=HeightDeposition/layer!小量!環(huán)境溫度!初始溫度!對(duì)流換熱系數(shù)!定義常量!掃描速度!每個(gè)載荷步位移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

5、,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,HeightDepos

6、ition,0,WidthClad,VGLUE,ALLNUMCMP,ALL=劃分網(wǎng)格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,LE

7、SIZE,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,!長(zhǎng)度方向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-D

8、TLSEL,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!網(wǎng)格劃分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLVSEL,S,LOC,Y,-1,0!網(wǎng)格劃分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLALLSEL,ALL=基底邊界條件、初始條件NSEL, S, LOC, Y, -HeightBase, 0!基底

9、初始溫度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,

10、 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!基底側(cè)面，!基底左端面，!基底右端面，!分析類(lèi)型：瞬態(tài)TRNOPT,FULL!求解方法：完全的N-R方法!*LNSRCH, on!*!DELTI虬0.01,0.001,0.05!載荷子步（默認(rèn)子步時(shí)間步長(zhǎng)、最小、最大）一一載荷步為0.333NSUBST, 4CNV

11、TOL,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!激活單元的上表面，指定為對(duì)流換熱邊條NSLE,S,1NSEL,R,LOC,Y,0NSEL,R,LOC,Z,0,RadiusSF,ALL,CONV,CoffConv,Temp1ALLSEL,ALL!預(yù)熱=*DO, i, 1, 2m=mod(i,2)!如*IF

12、,m,EQ,1,THEN果為奇數(shù)層，向右掃描*DO, k, 1, StepNum, 1TIME,TotalTime*(i-1)+StepTime*k!載荷步結(jié)束時(shí)間KBC, 1!載荷步內(nèi)載荷隨時(shí)間分布：常數(shù)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,

13、, LaserPower*Factor/Area,ALLSEL, ALLSOLVESAVESFEDELE,ALL,4,HFLUX!刪除熱流密度載荷 ALLSEL, ALLESEL,S,LIVE!激活單元的上表面，指定為對(duì)流換熱邊條NSLE,S,1NSEL,R,LOC,Y,0NSEL, R, LOC, x, LeftX-2*StepDis, RightX-2*StepDisNSEL, R, LOC, Z, O, RadiusSF,ALL,CONV,CoffConv,TempESEL,S,LIVEEPLOT*ENDDO*ELSE!如果為偶數(shù)層，向左掃描*DO, k, 1, StepNum, 1TI

14、ME,TotalTime*(i-1)+StepTime*k!載荷步結(jié)束時(shí)間KBC, 1!載荷步內(nèi)載荷隨時(shí)間分布：常數(shù)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

15、, ALLSOLVESAVESFEDELE,ALL,4,HFLUX!激活單元的上表面，刪除載荷ESEL,S,LIVE!激活單元的上表面，指定為對(duì)流換熱邊條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

16、,1,THEN!如果為奇數(shù)層，向右掃描*DO, k, 1, StepNum, 1TIME,TotalTime*(i+1)+StepTime*k!載荷步結(jié)束時(shí)間KBC, 1!載荷步內(nèi)載荷隨時(shí)間分布：常數(shù)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

17、/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!激活單元的表面，如果包含左端面，指定為對(duì)流換熱邊條NSLE,S,1NSEL,R,LOC,X,0NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposit

18、ion/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活單元的表面，如果包含右端面，指定為對(duì)流換熱邊條NSLE,S,1NSEL,R,LOC,X,LengthNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempESEL,S,LIVE!激活單元的側(cè)面，指定為對(duì)流換熱邊條NSLE,S,1NSEL,R,LOC,Z,WidthCladNSEL,R,LOC,Y,HeightDeposition/Layer*(i-

19、1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLSOLVESAVESFEDELE,ALL,4,HFLUX!刪除熱流密度載荷ALLSEL, ALLESEL,S,LIVE!激活單元的上表面，指定為對(duì)流換熱邊條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

20、/Layer*(i-1)!激活單元的下表面，刪除對(duì)流換熱邊條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, CONVESEL,S,LIVEEPLOT*ENDDO*ELSE!如果為偶數(shù)層，向左掃描*DO, k, 1,

21、StepNum, 1TIME,TotalTime*(i+1)+StepTime*k!載荷步結(jié)束時(shí)間KBC, 1!載荷步內(nèi)載荷隨時(shí)間分布：常數(shù)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), He

22、ightDeposition/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!激活單元的表面，如果包含左端面，指定為對(duì)流換熱邊條NSLE,S,1NSEL,R,LOC,X,0NSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,

23、CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活單元的表面，如果包含右端面，指定為對(duì)流換熱邊條NSLE,S,1NSEL,R,LOC,X,LengthNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), HeightDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVE!激活單元的側(cè)面，指定為對(duì)流換熱邊條NSLE,S,1NSEL,R,LOC,Z,WidthCladNSEL,R,LOC,Y,HeightDeposition/Layer*(i-1), Heig

24、htDeposition/Layer*iSF,ALL,CONV,CoffConv,TempALLSEL,ALLSOLVESAVESFEDELE,ALL,4,HFLUX!激活單元的上表面，刪除載荷ALLSEL, ALLESEL,S,LIVE!激活單元的上表面，指定為對(duì)流換熱邊條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/La

25、yer*(i-1)!激活單元的下表面，刪除對(duì)流換熱邊條 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,L

26、IVEEPLOT=冷卻!=!100s*DO, k, 1, 2, 1TIME,TotalTime*(Layer+2)+50*k!載荷步結(jié)束時(shí)間NSUBST, 5KBC, 1SOLVESAVE*ENDDO=!1000s*DO, k, 1, 9, 1TIME,TotalTime*(Layer+2)+100+100*k!載荷步結(jié)束時(shí)間NSUBST, 5KBC, 1SOLVESAVE*ENDDO!=!3000s*DO, k, 1, 10, 1TIME,TotalTime*(Layer+2)+1000+200*k!載荷步結(jié)束時(shí)間NSUBST, 5KBC, 1SOLVESAVE*ENDDO!=!10000s

27、*DO, k, 1, 14, 1TIME,TotalTime*(Layer+2)+3000+500*k!載荷步結(jié)束時(shí)間NSUBST, 5KBC, 1SOLVESAVE*ENDDOFINISH【我的碩士論文的一部分】求解應(yīng)力場(chǎng)!修改速度、時(shí)間子步步長(zhǎng)、載荷文件位置!如果修改基底的熱膨脹系數(shù)，要修改宏文件!單位制:米、秒、攝氏度/CLEAR,START/FILNAME,stress,0/COM,ANSYS RELEASE 10.0UP2005071820:15:5209/10/2007/CONFIG, NRES, 5000/PREP7/PAGE, 1000, , 1000,/VIEW,1,1,2,

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

29、.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,9

30、00TBPT,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!基底

31、高度Length=0.09!長(zhǎng)度WidthClad=0.0015!寬度HeightDeposition=0.00375!覆層高度Layer=15!層數(shù)HeightClad=HeightDeposition/layerdt=0.0001!小量Temp=20!環(huán)境溫度InitTemp=300!初始溫度CoffConv=30I!對(duì)流換熱系數(shù)!=!定義常量=Velocity=0.003!掃描速度StepDis=0.0015!每個(gè)載荷步位移LaserPower=700Radius=0.0015Area=3.14159265*(Radius*2)Factor=0.3!激光功率!激光光斑半徑!激光光斑面積!

32、吸收因子StepTime=StepDis/VelocityTotalTime=(Length+Radius*2)/VelocityStepNum=(Length+Radius*2)/StepDis!每個(gè)載荷步時(shí)間!載荷持續(xù)時(shí)間（掃描一層）!載荷步數(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,

33、-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!=!劃分網(wǎng)格LSEL, S, LOC, Y, dt, HeightDeposition-dt,!高度方向LESIZE, ALL, , , Layer,LSEL,S,LOC,Y,-DT,-0.0066+DTLSEL,R,LOC,X,0LESIZE,ALL,4,

34、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,!長(zhǎng)度方向LESIZE, ALL, , , Length/StepDis,LSEL,

35、 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,0VSEL,S,LOC,Y,0,HeightDeposition!網(wǎng)格劃分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,A

36、LLVSEL,S,LOC,Y,-1,0網(wǎng)格劃分TYPE,1MAT,1MSHAPE,0,3DMSHKEY,1VMESH,ALLALLSEL,ALLVSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DT!刪除熱影響區(qū)外的單元VSEL,A,LOC,Z,0.0067+DT,WidthBase-DTVCLEAR,ALLVDELE,ALL, , ,1ALLSEL,ALLFINISH=瞬態(tài)分析參數(shù)設(shè)置/SOLANTYPE,4!分析類(lèi)型：瞬態(tài)TRNOPT,FULL!求解方法：對(duì)于材料NLGEOM,on非線(xiàn)性，這是唯一的方法!大變形分析LNSRCH, onCNVTOL,U,0.05,2,C

37、NVTOL,F,0.01,2,OUTRES,NSOL,LASTNSUBST, 4NEQIT,30TREF, Temp!收斂準(zhǔn)則：控制位移!收斂準(zhǔn)則：控制力!結(jié)果輸出：!=NSEL, S, LOC, Y,!殺夕匕單兀0, HeightDepositionESLN, S, 1 EKILL,ALL ALLSEL,ALL!殺死熔覆層單元ESEL,S,LIVEEPLOTNSEL,S,LOC,Z,0束，UZ,ROTX,ROTY)D,ALL,UZ,0!對(duì)稱(chēng)邊條（相當(dāng)于三個(gè)約NSEL,S,LOC,Y,-0.0066(增加兩個(gè)約束，UX,UY)NSEL,R,LOC,Z,0NSEL,R,LOC,X,length/2D,ALL,ALL,0!固定中心點(diǎn)NSEL,S,LOC,Y,-0.0066,0束，ROTZ)NSEL,R,LOC,Z,0NSEL,R,LOC,X,leng