焊接模拟ansys实例PDF.pdf

【我的硕士论文的一部分】求解温度场!上表面上没有对流换热边界条件!单位制:米、秒、摄氏度!/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【补充】其中的宏命令是加比较复杂的位移边条，可以去掉，然后加上你需要的边界条件即可；。

Ansys静力分析实例：　1 问题描述：　如图所示支架简图，支架材料为结构钢，厚度10mm，支架左侧的两个通孔为固定孔，顶面的开槽处受均布载荷，载荷大小为500N/mm。

　2 启动Ansys Workbench，在界面中选择Simulation启动DS模块。

3 导入三维模型，操作步骤按下图进行，单击“Geometry”，选择“From File”。

　从弹出窗口中选择三维模型文件，如果文件格式不符，可以把三维图转换为“.stp”格式文件，即可导入，如下图所示。

　4 选择零件材料：文件导入后界面如下图所示，这时，选择“Geometry”下的“Part”，在左下角的“Details of ‘Part’”中可以调整零件材料属性。

　5 划分网格：如下图，选择“Project”树中的“Mesh”，右键选择“Generate Mesh”即可。

【此时也可以在左下角的“Details of‘Mesh’”对话框中调整划分网格的大小（“Element size”项）】。

　生成网格后的图形如下图所示：　6 添加分析类型：选择上方工具条中的“New Analysis”，添加所需做的分析类型，此例中要做的是静力分析，因此选择“Static Structural”，如下图所示。

　7 添加固定约束：如下图所示，选择“Project”树中的“Static Structural”，右键选择“Insert”中的“Fixed Support”。

　这时左下角的“Details of ‘Fixed Support’”对话框中“Geometry”被选中，提示输入固定支撑面。

本例中固定支撑类型是面支撑，因此要确定图示6位置为“Face”，【此处也可选择“Edge”来选择“边”】然后按住“CTRL”键，连续选择两个孔面为支撑面，按“Apply”确认，如下图所示。

　8 添加载荷：选择“Project”树中的“Static Structural”，右键选择“Insert”中的“Force”，如下图所示。

焊接模拟ansys实例!下面的命令流进行的是一个简单的二维焊接分析, 利用ANSYS单元生死和热-结构耦合分析功能进!行焊接过程仿真, 计算焊接过程中的温度分布和应力分布以及冷却后的焊缝残余应力。

finish/clear/filnam,1-2D element birth and death/title,Weld Analysis by "Element Birth and Death"/prep7/unit,si !采用国际单位制!******************************************************et,1,13,4 !13号二维耦合单元, 同时具有温度和位移自由度et,2,13,4!1号材料是钢!2号材料是铝!3号材料是铜!铝是本次分析中的焊料, 它将钢结构部分和铜结构部分焊接起来!下面是在几个温度点下, 各材料的弹性模量mptemp,1,20,500,1000,1500,2000mpdata,ex,1,1,1.93e11,1.50e11,0.70e11,0.10e11,0.01e11mpdata,ex,2,1,1.02e11,0.50e11,0.08e11,0.001e11,0.0001e11mpdata,ex,3,1,1.17e11,0.90e11,0.30e11,0.05e11,0.005e11!假设各材料都是双线性随动硬化弹塑性本构关系!下面是各材料在各温度点下的屈服应力和屈服后的弹性模量tb,bkin,1,5tbtemp,20,1tbdata,1,1200e6,0.193e11tbtemp,500,2tbdata,1, 933e6,0.150e11tbtemp,1000,3tbdata,1, 435e6,0.070e11tbtemp,1500,4tbdata,1, 70e6,0.010e11tbtemp,2000,5tbdata,1, 7e6,0.001e11!tb,bkin,2,5tbtemp,20,1tbdata,1,800e6,0.102e11tbtemp,500,2tbdata,1,400e6,0.050e11tbtemp,1000,3tbdata,1, 70e6,0.008e11tbdata,1, 1e6,0.0001e11tbtemp,2000,5tbdata,1,0.1e6,0.00001e11!tb,bkin,3,5tbtemp,20,1tbdata,1,900e6,0.117e11tbtemp,500,2tbdata,1,700e6,0.090e11tbtemp,1000,3tbdata,1,230e6,0.030e11tbtemp,1500,4tbdata,1, 40e6,0.005e11tbtemp,2000,5tbdata,1, 4e6,0.0005e11!!材料密度(假设为常值)mp,dens,1,8030mp,dens,2,4850mp,dens,3,8900! 热膨胀系数(假设为常值)mp,alpx,1,1.78e-5mp,alpx,2,9.36e-6mp,alpx,3,1.66e-5!泊松比(假设为常值)mp,nuxy,1,0.29mp,nuxy,2,0.30mp,nuxy,3,0.30!热传导系数(假设为常值)mp,kxx,1,16.3mp,kxx,2,7.44mp,kxx,3,393!比热(假设为常值)mp,c,1,502mp,c,2,544mp,c,3,385!热膨胀系数(假设为常值)!由于该13号单元还有磁自由度, 此处假设一磁特性, 但并不影响我们所关心的结果mp,murx,1,1mp,murx,2,1mp,murx,3,1!假设焊料(铝)焊上去后的初始温度是1500℃mp,reft,2,1500mp,reft,3,20!******************************************************!下面建立几何模型csys,0k,1,0,0,0k,2,0.5,0,0k,3,1,0,0 !长1米k,4,0,0.3,0 !厚度0.3米（二维中叫做宽度）k,5,0.35,0.3,0k,6,0.65,0.3,0k,7,1,0.3,0a,1,2,5,4a,2,6,5a,2,3,7,6!!!!!!!!!!!!!!!!!!!!!!!!!!划分网格!!!!!!!!!!!!!!!!!esize,0.025type,2mat,2amesh,2!esize,0.05 !网格划分出现问题type,1mat,1amesh,1!mat,3amesh,3eplot!/soluantype,4 ! 瞬态分析trnopt,full!!!!!!!!!!!!!!!!!!!!!!!!在模型的左边界加位移约束!!!!!!!!!!!!!!!!!!!!!!!!!!!nsel,all*get,minx,node,,mnloc,xnsel,s,loc,x,minxd,all,ux,0*get,miny,node,,mnloc,ynsel,r,loc,y,minyd,all,uy,0!*****假设模型的左右边界处温度始终保持在20摄氏度左右*****!其他边界条件如对流和辐射等均可施加，此处因为只是示意而已，故只施加恒温边界条件nsel,all*get,minx,node,,mnloc,xnsel,s,loc,x,minxd,all,temp,20nsel,all*get,maxx,node,,mxloc,xnsel,s,loc,x,maxxd,all,temp,20!!!!!!!!!!!!!!!由于第2个面是焊接所在区域，因此首先将该区域的单元“杀死”!!!!!!!!!!!!!!!!nna=2esel,all*get,emax,elem,,num,maxasel,s,area,,nnaesla*get,nse,elem,,count*dim,ne,,nse*dim,nex,,nse*dim,ney,,nse*dim,neorder,,nsemine=0!**********************************************!下面的do循环用于将焊料区的单元按其形心y坐标排序!以便后面模拟焊料由下向上逐步“生长”过程*do,i1,1,nseesel,u,elem,,mine*get,nse1,elem,,countii=0*do,i,1,emax*if,esel(i),eq,1,thenii=ii+1ne(ii)=i*endif*enddo*do,i,1,nse1*get,ney(i),elem,ne(i),cent,y*get,nex(i),elem,ne(i),cent,x*enddominy=1e20minx=1e20*do,i,1,nse1*if,ney(i),lt,miny,thenminy=ney(i)minx=nex(i)mine=ne(i)*else*if,ney(i),eq,miny,then*if,nex(i),lt,minx,thenminy=ney(i)minx=nex(i)mine=ne(i)*endif*endif*endif*enddoneorder(i1)=mine*enddo!************************************************************** max_tem=1500 !按照前面假设，焊料的初始温度为1500℃dt1=1e-3 !用于建立初始条件的一个很小的时间段dt=5 !焊接一个单元所需的时间t=0 !起始时间esel,alleplot/auto,1/replot*do,i,1,nseekill,neorder(i)esel,s,liveeplot*enddoallsel,alloutres,all,allic,all,temp,20kbc,1timint,0,structtimint,1,thermtimint,0,magtintp,0.005,,,1,0.5,0.2!nsub1=2nsub2=40!**************************************************do,i,1,nseealive,neorder(i)esel,s,liveeplotesel,all!******下面的求解用于建立温度的初始条件******t=t+dt1time,tnsubst,1*do,j,1,4d,nelem(neorder(i),j),temp,max_tem*enddosolve!****下面的求解用于保证初始的升温速度为零****t=t+dt1time,tsolve!*********下面的步骤用于求解温度分布***********do,j,1,4ddele,nelem(neorder(i),j),temp*enddot=t+dt-2*dt1time,tnsubst,nsub1solve*enddot=t+50000 !*********下面的步骤用于冷却过程求解***** time,tnsubst,nsub2solvesavefinish!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!后处理过程!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!/post1!**************下面的一系列命令用于生成应力的动画文件******* /seg,dele/cont,1,15,0,1200e6/16,1200e6/dscale,1,1.0avprin,0,0avres,1/seg,multi,stress1,0.1esel,all*do,i,1,nseesel,u,elem,,neorder(i)*enddo*do,i,1,nseesel,a,elem,,neorder(i)set,(i-1)*3+1,1plnsol,s,eqv*do,j,1,nsub1set,(i-1)*3+3,jplnsol,s,eqv*enddo*enddo*do,i,1,nsub2set,(nse-1)*3+4,iplnsol,s,eqv*enddo/seg,off,stress1,0.1/anfile,save,stress1,avi!**********下面的一系列命令用于生成温度的动画文件************ /seg,dele/cont,1,15,0,1500/16,1500/dscale,1,1.0avprin,0,0avres,1/seg,multi,temp1,0.1esel,all*do,i,1,nseesel,u,elem,,neorder(i)*enddo*do,i,1,nseesel,a,elem,,neorder(i)set,(i-1)*3+1,1plnsol,temp*do,j,1,nsub1set,(i-1)*3+3,jplnsol,temp*enddo*enddo*do,i,1,nsub2set,(nse-1)*3+4,iplnsol,temp*enddo/seg,off,temp1,0.1/anfile,save,temp1,avifinish。

Ansys静力分析实例：　1 问题描述：　如图所示支架简图，支架材料为结构钢，厚度10mm，支架左侧的两个通孔为固定孔，顶面的开槽处受均布载荷，载荷大小为500N/mm。

　2 启动Ansys Workbench，在界面中选择Simulation启动DS模块。

3 导入三维模型，操作步骤按下图进行，单击“Geometry”，选择“From File”。

　从弹出窗口中选择三维模型文件，如果文件格式不符，可以把三维图转换为“.stp”格式文件，即可导入，如下图所示。

　4 选择零件材料：文件导入后界面如下图所示，这时，选择“Geometry”下的“Part”，在左下角的“Details of ‘Part’”中可以调整零件材料属性。

　5 划分网格：如下图，选择“Project”树中的“Mesh”，右键选择“Generate Mesh”即可。

【此时也可以在左下角的“Details of‘Mesh’”对话框中调整划分网格的大小（“Element size”项）】。

　生成网格后的图形如下图所示：　6 添加分析类型：选择上方工具条中的“New Analysis”，添加所需做的分析类型，此例中要做的是静力分析，因此选择“Static Structural”，如下图所示。

　7 添加固定约束：如下图所示，选择“Project”树中的“Static Structural”，右键选择“Insert”中的“Fixed Support”。

　这时左下角的“Details of ‘Fixed Support’”对话框中“Geometry”被选中，提示输入固定支撑面。

本例中固定支撑类型是面支撑，因此要确定图示6位置为“Face”，【此处也可选择“Edge”来选择“边”】然后按住“CTRL”键，连续选择两个孔面为支撑面，按“Apply”确认，如下图所示。

　8 添加载荷：选择“Project”树中的“Static Structural”，右键选择“Insert”中的“Force”，如下图所示。

焊接过程模拟与焊接变形、焊接Ansys应力有限元分析1.1 焊接变形与焊接应力焊接时，加热和冷却循环总会导致一定程度的变形，焊接变形对尺寸稳定性以与结构力学性能都有很大的影响，控制焊接变形在焊接加工中是一个关键的任务。

在钢结构焊接中，焊接工艺会使构件温度场产生不均匀变化，从而在构件中产生复杂的剩余应力分布。

剩余应力是一种自相平衡的力系，当构件承受荷载时，如受拉、受压等，荷载引起的应力将与截面剩余应力相叠加，从而使构件某些部位提前达到屈服强度，并发生塑性变形，故会严重降低构件的刚度和稳定性以与结构疲劳强度。

对构件进展焊接，在焊件上产生局部高温的不均匀温度场，焊接中心处温度可达1600℃，高温区的钢材会发生较大程度的膨胀伸长，但受到相邻钢材的约束，从而在焊件内引起较高的温度应力，并在焊接过程中，随时间和温度而不断变化，称其为焊接应力。

焊接应力较高的部位，甚至将达到钢材的屈服强度而发生塑性变形，因而钢材冷却后将有残存于焊件内的应力，称为焊接剩余应力。

并且在冷却过程中，钢材由于不能自由收缩，而受到拉伸，于是焊件中出现了一个与焊件加热方向大致相反的内应力场。

1.2 Ansys有限元焊接分析为通过对焊接过程的三维有限元模拟分析以与焊接后构件变形与剩余应力分布分析，为评估焊接对焊件的影响提供更加合理、有效、可靠的分析数据，并为焊接工艺提供一定的指导，为采用的焊接过程提供一定的分析依据，采用大型有限元计算软件Ansys作为分析工具对焊接过程与焊件的变形与剩余应力进展了分析。

ANSYS有2种方式来考虑热分析与力学分析之间的耦合,即直接耦合和间接耦合。

间接耦合法的处理思路为先进展温度场的模拟,然后将求出的结点温度作为体载荷施加在结构中,计算焊接剩余应力与变形。

即：(1)使用热分析的手段进展热分析，根据需要可采用瞬态分析与稳态分析模型，此处为瞬态分析。

(2)重新进入前处理中,将热分析单元转换为相应的结构分析单元，设置结构分析中材料属性,如弹性模量、泊松比、热膨胀系数等。

ansys模拟焊接/FILNAME,weldingstress,0 /PREP7!定义焊接参数L=1E-1!焊件的长度W=1E-1!焊件的宽度H=6E-3!焊件的高度U=20!焊接电压I=160!焊接电流V=0.01!焊接速度YITA=0.7!焊接热效率R=0.007!电弧有效加热半径Q=U*I*YITA!电弧热功率Qm=3/3.1415/R**2*Q!加热斑点中心最大热流密度!*ET,1,PLANE55ET,2,SOLID70MPTEMP,,,,,,,,MPTEMP,1,20MPTEMP,2,200MPTEMP,3,500MPTEMP,4,750MPTEMP,5,1000MPTEMP,6,1500MPTEMP,7,1700MPTEMP,8,2500MPDATA,KXX,1,,50 MPDATA,KXX,1,,47 MPDATA,KXX,1,,40 MPDATA,KXX,1,,27 MPDATA,KXX,1,,30 MPDATA,KXX,1,,35 MPDATA,KXX,1,,40 MPDATA,KXX,1,,55 MPDATA,DENS,1,,7820 MPDATA,DENS,1,,7700 MPDATA,DENS,1,,7610 MPDATA,DENS,1,,7550 MPDATA,DENS,1,,7490 MPDATA,DENS,1,,7350 MPDATA,DENS,1,,7300 MPDATA,DENS,1,,7090 MPDATA,C,1,,460 MPDATA,C,1,,480 MPDATA,C,1,,530 MPDATA,C,1,,675 MPDATA,C,1,,670 MPDATA,C,1,,660 MPDATA,C,1,,780 MPDATA,C,1,,820 MPDATA,EX,1,,2.05E11 MPDATA,EX,1,,1.87E11 MPDATA,EX,1,,1.5E11 MPDATA,EX,1,,0.7E11 MPDATA,EX,1,,0.2E11 MPDATA,EX,1,,0.19E2 MPDATA,EX,1,,0.18E2 MPDATA,EX,1,,0.12e2 MPDATA,PRXY,1,,0.28 MPDATA,PRXY,1,,0.29 MPDATA,PRXY,1,,0.31 MPDATA,PRXY,1,,0.35 MPDATA,PRXY,1,,0.4 MPDATA,PRXY,1,,0.45 MPDATA,PRXY,1,,0.48MPDATA,PRXY,1,,0.5 UIMP,1,REFT,,,20 MPDATA,ALPX,1,,1.1e-5 MPDATA,ALPX,1,,1.22e-5 MPDATA,ALPX,1,,1.39e-5 MPDATA,ALPX,1,,1.48e-5 MPDATA,ALPX,1,,1.34e-5 MPDATA,ALPX,1,,1.33e-5 MPDATA,ALPX,1,,1.32e-5 MPDATA,ALPX,1,,1.31e-5 TB,BISO,1,6,2, TBTEMP,20 TBDATA,,220e6,0,,,, TBTEMP,250 TBDATA,,175e6,0,,,, TBTEMP,500 TBDATA,,80e6,0,,,, TBTEMP,750 TBDATA,,40E6,0,,,, TBTEMP,1000 TBDATA,,10E6,0,,,, TBTEMP,1500 TBDATA,,1E-5,0,,,,K,1,0,0,0K,2,0,L,0K,3,-W/2*0.15,L,0K,4,-W/2*0.3,L,0K,5,-W/2*0.5,L,0K,6,-W/2,L,0K,7,-W/2,0,0K,8,-W/2*0.5,0,0K,9,-W/2*0.3,0,0K,10,-W/2*0.15,0,0K,11,0,0,HA,1,2,3,10A,10,3,4,9A,9,4,5,8A,8,5,6,7ESIZE,0.0012 AMESH,1ESIZE,0.0025 AMESH,2ESIZE,0.005 AMESH,3ESIZE,0.0065 AMESH,4TYPE,2 EXTOPT,ESIZE,2,0, EXTOPT,ACLEAR,1!*EXTOPT,ATTR,1,0,0 REAL,_Z4ESYS,0!*VOFFST,1,H,, VOFFST,2,H,, VOFFST,3,H,, VOFFST,4,H,,EPLOTNUMMRG,ALL,,,,LOW /SOL!*ANTYPE,4!*TRNOPT,FULL LUMPM,0*DEL,_FNCNAME*DEL,_FNCMTID*DEL,_FNC_C1*DEL,_FNC_C2*DEL,_FNC_C3*DEL,_FNCCSYS*SET,_FNCNAME,'GAOSI'*DIM,_FNC_C1,,1*DIM,_FNC_C2,,1*DIM,_FNC_C3,,1*SET,_FNC_C1(1),QM*SET,_FNC_C2(1),V*SET,_FNC_C3(1),R*SET,_FNCCSYS,0!/INPUT,HANJIE.func,,,1*DIM,%_FNCNAME%,TABLE,6,19,1,,,,%_FNCCSYS%!!Begin of equation:Qm*exp(-3*({X}^2+({Y}-V*{TIME})^2)/R^2) *SET,%_FNCNAME%(0,0,1),0.0,-999*SET,%_FNCNAME%(2,0,1),0.0*SET,%_FNCNAME%(3,0,1),%_FNC_C1(1)%*SET,%_FNCNAME%(4,0,1),%_FNC_C2(1)%*SET,%_FNCNAME%(5,0,1),%_FNC_C3(1)%*SET,%_FNCNAME%(6,0,1),0.0*SET,%_FNCNAME%(0,1,1),1.0,-1,0,0,0,0,0*SET,%_FNCNAME%(0,2,1),0.0,-2,0,1,0,0,-1*SET,%_FNCNAME%(0,3,1),0,-3,0,1,-1,2,-2*SET,%_FNCNAME%(0,4,1),0.0,-1,0,3,0,0,-3*SET,%_FNCNAME%(0,5,1),0.0,-2,0,1,-3,3,-1*SET,%_FNCNAME%(0,6,1),0.0,-1,0,2,0,0,2*SET,%_FNCNAME%(0,7,1),0.0,-3,0,1,2,17,-1*SET,%_FNCNAME%(0,8,1),0.0,-1,0,1,18,3,1*SET,%_FNCNAME%(0,9,1),0.0,-4,0,1,3,2,-1*SET,%_FNCNAME%(0,10,1),0.0,-1,0,2,0,0,-4*SET,%_FNCNAME%(0,11,1),0.0,-5,0,1,-4,17,-1*SET,%_FNCNAME%(0,12,1),0.0,-1,0,1,-3,1,-5*SET,%_FNCNAME%(0,13,1),0.0,-3,0,1,-2,3,-1*SET,%_FNCNAME%(0,14,1),0.0,-1,0,2,0,0,19*SET,%_FNCNAME%(0,15,1),0.0,-2,0,1,19,17,-1*SET,%_FNCNAME%(0,16,1),0.0,-1,0,1,-3,4,-2*SET,%_FNCNAME%(0,17,1),0.0,-1,7,1,-1,0,0*SET,%_FNCNAME%(0,18,1),0.0,-2,0,1,17,3,-1*SET,%_FNCNAME%(0,19,1),0.0,99,0,1,-2,0,0!End of equation:Qm*exp(-3*({X}^2+({Y}-V*{TIME})^2)/R^2) TUNIF,20,!定义初始温度!定义对流换热边界SFA,15,1,CONV,30,20SFA,20,1,CONV,30,20SFA,9,1,CONV,30,20SFA,14,1,CONV,30,20SFA,19,1,CONV,30,20SFA,24,1,CONV,30,20SFA,23,1,CONV,30,20SFA,7,1,CONV,30,20SFA,12,1,CONV,30,20SFA,17,1,CONV,30,20SFA,22,1,CONV,30,20SFA,1,1,CONV,30,20SFA,2,1,CONV,30,20SFA,3,1,CONV,30,20SFA,4,1,CONV,30,20!施加高斯热源SFA,5,1,HFLUX,%GAOSI% SFA,10,1,HFLUX,%GAOSI% OUTRES,ALL,ALL,TIME,L/V!设置求解时间AUTOTS,-1NSUBST,50,50,50KBC,0TSRES,ERASELSWRITE,1,!写入载荷文件为1 !*TIME,20AUTOTS,1NSUBST,20,20,20KBC,0!*TSRES,ERASELSWRITE,2,TIME,50AUTOTS,1NSUBST,30,30,30KBC,0!*TSRES,ERASELSWRITE,3,TIME,1100AUTOTS,1NSUBST,105,105,105 KBC,0!*TSRES,ERASELSWRITE,4,LSSOLVE,1,4,1,!开始求解!进入热应力求解/PREP7ETCHG,TTS!热到结构分析转换/SOLANTYPE,4TRNOPT,FULLLUMPM,0NLGEOM,1NROPT,FULL,,OFF!定义热应力计算参考温度TREF,20,DA,6,SYMMDL,13,,UY,DL,10,,UY,DL,7,,UY,DL,4,,UY,DL,1,,UZ,*DO,I,1,50LDREAD,TEMP,,,0.2*I,,'weldingstress','rth', ''!读入热分析的计算结果OUTRES,ALL,ALL,TIME,0.2*IDELTIM,0.2,0.075,0.2,1SOLVE*ENDDO*DO,I,1,20LDREAD,TEMP,,,10+I*0.5,,'weldingstress','rth&#3 9;,''OUTRES,ALL,ALL,TIME,10+I*0.5DELTIM,0.5,0.5,1,1SOLVE*ENDDO*DO,I,1,30LDREAD,TEMP,,,20+I,,'weldingstress','rth',& #39;'OUTRES,ALL,ALL,TIME,20+IDELTIM,1,1,1,1SOLVE*ENDDO*DO,I,1,105LDREAD,TEMP,,,50+10*I,,'weldingstress','rth&#3 9;,''OUTRES,ALL,ALL,TIME,50+10*IDELTIM,10,10,10,1SOLVE*ENDDO。