ansys求解焊接温度场

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 !* !*

!===============================================================================================!材料属性

!=====================================================================================!316L MPTEMP,,,,,,,, MPTEMP,1,0

MPDATA,DENS,1,,7850 MPTEMP,,,,,,,, MPTEMP,1,20 MPTEMP,2,300 MPTEMP,3,900 MPTEMP,4,1400 MPTEMP,5,2000 MPDATA,KXX,1,,18.6 MPDATA,KXX,1,,21.4 MPDATA,KXX,1,,28.4 MPDATA,KXX,1,,33.9 MPDATA,KXX,1,,48 MPTEMP,,,,,,,, MPTEMP,1,20 MPTEMP,2,600 MPTEMP,3,800 MPTEMP,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.131e9 MPDATA,ENTH,1,,7.347e9 MPDATA,ENTH,1,,9.145e9 MPDATA,ENTH,1,,1.03e10 MPDATA,ENTH,1,,1.272e10

!===============================================================================================!定义常量

WidthBase=0.025 !宽度 HeightBase=0.02 !基底高度 Length=0.09 !长度

WidthClad=0.0015 !宽度

HeightDeposition=0.00375 !覆层高度

Layer=15 !层数 HeightClad=HeightDeposition/layer

dt=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.0067

BLOCK,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, 14

BLOCK,0,Length,0,HeightDeposition,0,WidthClad,

VGLUE,ALL NUMCMP,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,2

LSEL,S,LOC,Y,-DT,-0.0066+DT LSEL,R,LOC,X,Length LESIZE,ALL,,,4,0.5

LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DT LSEL,R, LOC, Z, 0, WidthBase-DT, LESIZE,ALL,,,4,2

LSEL,S,LOC,Y,-0.0066-DT,-HeightBase+DT LSEL,R, LOC, Z, WidthBase, LESIZE,ALL,,,4,0.5

LSEL, 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-DT LSEL,R,LOC,Y,-DT,-HeightBase LESIZE,ALL,,,4,2

LSEL,S,LOC,Z,0.0067+DT,WidthBase-DT LSEL,R,LOC,Y,0 LESIZE,ALL,,,4,0.5

VSEL,S,LOC,Y,0,HeightDeposition !网格划分 TYPE,1 MAT,1

MSHAPE,0,3D MSHKEY,1 VMESH,ALL

VSEL,S,LOC,Y,-1,0 !网格划分 TYPE,1 MAT,1

MSHAPE,0,3D MSHKEY,1 VMESH,ALL

ALLSEL,ALL

!===============================================================================================!基底边界条件、初始条件

NSEL, S, LOC, Y, -HeightBase, 0 !基底初始温度 IC,ALL,TEMP,InitTemp ALLSEL,ALL

NSEL, S, LOC, Y, -HeightBase, 0 !基底侧面，换热边界条件 NSEL, R, LOC, Z, WidthBase SF, ALL, CONV, CoffConv, Temp ALLSEL,ALL

NSEL, S, LOC, Y, -HeightBase, 0 !基底左端面，换热边界条件 NSEL, R, LOC, X, 0

SF, ALL, CONV, CoffConv, Temp ALLSEL,ALL

NSEL, S, LOC, Y, -HeightBase, 0 !基底右端面，换热边界条件 NSEL, R, LOC, X, Length

SF, ALL, CONV, CoffConv, Temp ALLSEL,ALL

NSEL, S, LOC, Y, 0 !基底上表面，换热边界条件 NSEL, R, LOC, Z, WidthClad, WidthBase SF, ALL, CONV, CoffConv, Temp ALLSEL,ALL

FINISH

/SOLU

!===============================================================================================!瞬态分析参数设置 ANTYPE,4 !分析类型：瞬态 !*

TRNOPT,FULL !求解方法：完全的N-R方法 LNSRCH, on !*

!DELTIM,0.01,0.001,0.05 !载荷子步(默认子步时间步长、最小、最大)——载荷步为0.333 NSUBST, 4

CNVTOL,HEAT, ,0.01,2,0.000001, !收敛准则：控制热流 OUTRES,NSOL,LAST !结果输出：所有

!===============================================================================================!杀死单元

NSEL, S, LOC, Y, 0, HeightDeposition !杀死熔覆层单元 ESLN, S, 1, ALL EKILL,ALL ALLSEL,ALL

ESEL,S,LIVE EPLOT

ESEL,S,LIVE !激活单元的上表面，指定为对流换热边条 NSLE,S,1

NSEL,R,LOC,Y,0

NSEL,R,LOC,Z,0,Radius

SF,ALL,CONV,CoffConv,Temp ALLSEL,ALL

!===============================================================================================!预热 *DO, i, 1, 2

m=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,LIVE EPLOT

!===============================================================================================!熔覆

*DO, i, 1, Layer, 1

m=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

本文来源：https://www.bwwdw.com/article/7h6f.html