1、【我的硕士论文的一部分】求解温度场!上表面上没有对流换热边界条件!单位制:米、秒、摄氏度!/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,0 MPDATA,DENS,1,7850 MPTEMP, MPTEMP,1,20 MPTEMP,2,300MPTEMP
2、,3,900MPTEMP,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,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,130
3、0 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 !覆层
4、高度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=(Leng
5、th+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,
6、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,ALLNUMCMP,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,2LSEL,S,LOC,Y,-DT,-0.
7、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,
8、 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 !网格划分TYPE,1MAT,1MSHAPE,0,3D MSHKEY,1
9、VMESH,ALLVSEL,S,LOC,Y,-1,0 !网格划分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 !基底侧面，换热边界条件NSEL, R, LOC, Z, WidthBaseSF, ALL, CONV, CoffConv, TempALLSEL,ALLNSEL, S, LOC, Y, -HeightBa
10、se, 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/SO
11、LU!=!瞬态分析参数设置ANTYPE,4 !分析类型：瞬态!* TRNOPT,FULL !求解方法：完全的 N-R 方法LNSRCH, on!* !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,LIV
12、EEPLOTESEL,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
13、, 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
14、, 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,
15、 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,
16、 x, LeftX-2*StepDis, RightX-2*StepDisNSEL, R, LOC, Z, O, RadiusSF,ALL,CONV,CoffConv,TempALLSEL,ALLESEL,S,LIVEEPLOT*ENDDO *ENDIF*ENDDO ESEL,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
17、=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
18、!激活单元的上表面，加热流密度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,
19、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,H
20、FLUX !删除热流密度载荷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, Left
21、X-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 !载荷步内载荷随时间分布：常数
22、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, L
23、OC, 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 !激活单元的表面，如果包含右端面，指定为对流换热
24、边条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,ALLSOLVESA
25、VESFEDELE,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,
26、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*ENDDO ESEL,S,LIVEEPLOT!=!冷却!=!100s*DO, k, 1, 2, 1TIME,TotalTime*(Layer+2
27、)+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
28、!载荷步结束时间NSUBST, 5KBC, 1SOLVESAVE*ENDDOFINISH【我的硕士论文的一部分】求解应力场!修改速度、时间子步步长、载荷文件位置!如果修改基底的热膨胀系数，要修改宏文件!单位制:米、秒、摄氏度/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!* !* !=!材
29、料属性!=!316LMPTEMP, MPTEMP,1,0 MPDATA,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 MPDAT
30、A,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,300 TBPT,0,0 TBPT,5.5E-4,9.35E7TBPT,0.0012,1.27E8 TBPT,0.1,1.84E9TBTEMP,600 TBPT,0,0 TBPT,3.2E-4,4.8E7 TBPT,0.0012,7.19E7 TBPT,0.1,1.54E9TBTEMP,900 TBPT,0,0
31、TBPT,2.5E-4,1.25E7 TBPT,0.0012,5.1E7 TBPT,0.1,5.45E8TBTEMP,1300 TBPT,0,0 TBPT,2.5E-4,1E6 TBPT,0.00375,1.13E7 TBPT,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 !基底
32、高度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) !激光光斑面积
33、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
34、,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,ALLNUMCMP,ALL!=!划分网格LSEL, S, LOC, Y, dt, HeightDeposition-dt, !高度方向 LESIZE, ALL, , , Layer,LSEL,S,LOC,Y
35、,-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, Len
36、gth-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.5VSE
37、L,S,LOC,Y,0,HeightDeposition !网格划分TYPE,1MAT,1MSHAPE,0,3D MSHKEY,1VMESH,ALLVSEL,S,LOC,Y,-1,0 !网格划分TYPE,1MAT,1MSHAPE,0,3D MSHKEY,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,FUL
