ansys模拟焊接.docx

1、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,SOLID70 MPTEMP, MPTEMP,1,20 MPTEMP,2,200MPTEMP,3,500MPTEMP,4,750MPTEMP,5,1000

2、 MPTEMP,6,1500 MPTEMP,7,1700 MPTEMP,8,2500 MPDATA,KXX,1,50MPDATA,KXX,1,47MPDATA,KXX,1,40MPDATA,KXX,1,27MPDATA,KXX,1,30MPDATA,KXX,1,35MPDATA,KXX,1,40MPDATA,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 MPD

3、ATA,DENS,1,7090 MPDATA,C,1,460 MPDATA,C,1,480 MPDATA,C,1,530MPDATA,C,1,675MPDATA,C,1,670MPDATA,C,1,660MPDATA,C,1,780MPDATA,C,1,820 MPDATA,EX,1,2.05E11MPDATA,EX,1,1.87E11MPDATA,EX,1,1.5E11 MPDATA,EX,1,0.7E11 MPDATA,EX,1,0.2E11 MPDATA,EX,1,0.19E2MPDATA,EX,1,0.18E2MPDATA,EX,1,0.12e2 MPDATA,PRXY,1,0.28

4、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.48 MPDATA,PRXY,1,0.5 UIMP,1,REFT,20MPDATA,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 MP

5、DATA,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,1000TBDATA,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,0

6、K,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.0012AMESH,1ESIZE,0.0025AMESH,2ESIZE,0.005AMESH,3ESIZE,0.0065AMESH,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, , EPLOTNU

7、MMRG,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,%_F

8、NCNAME%,TABLE,6,19,1,%_FNCCSYS% ! ! Begin of equation: Qm*exp(-3*(X2+(Y-V*TIME)2)/R2)*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%(

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

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

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

12、, 0, 0! End of equation: Qm*exp(-3*(X2+(Y-V*TIME)2)/R2) TUNIF,20, !定义初始温度!定义对流换热边界SFA,15,1,CONV,30,20 SFA,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,CO

13、NV,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,-1 NSUBST,50,50,50KBC,0 TSRES,ERASE LSWRITE,1, !写入载荷文件为1!* TIME,20AUTOTS,1NSUBST,20,20,20KBC,0 !* TSRES,ERASE LSWRITE,2,TIME,50AUTOTS,1NSUBST,30

14、,30,30KBC,0 !* TSRES,ERASE LSWRITE,3,TIME,1100AUTOTS,1NSUBST,105,105,105KBC,0 !* TSRES,ERASE LSWRITE,4, LSSOLVE,1,4,1, !开始求解!进入热应力求解/PREP7 ETCHG,TTS !热到结构分析转换/SOLANTYPE,4TRNOPT,FULL LUMPM,0 NLGEOM,1 NROPT,FULL, ,OFF!定义热应力计算参考温度TREF,20,DA,6,SYMMDL,13, ,UY, DL,10, ,UY,DL,7, ,UY, DL,4, ,UY, DL,1, ,UZ,

15、*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, OUTRES,ALL,ALL, TIME,10+I*0.5DELTIM,0.5,0.5,1,1SOLVE*ENDDO*DO,I,1,30LDREAD,TEMP,20+I, ,weldingstress,rth, OUTRES,ALL,ALL, TIME,20+IDELTIM,1,1,1,1SOLVE*ENDDO*DO,I,1,105LDREAD,TEMP,50+10*I, ,weldingstress,rth, OUTRES,ALL,ALL, TIME,50+10*IDELTIM,10,10,10,1SOLVE*ENDDO