*** WARNING *** CP = 17.940 TIME= 14:48:45
Material property KXX of material 1 of element 3681 is evaluated at a
temperature of 0, which is below the supplied temperature range.
Temperature range checking terminates.
*** ERROR *** CP = 5984.198 TIME= 15:40:11
Solution not converged at time 1.216 (load step 13 substep 9).
Run terminated.
*** WARNING *** CP = 5984.432 TIME= 15:40:23
The unconverged solution (identified as time 1.3 substep 999999) is
output for analysis debug purposes. Results should not be used for
any other purpose.
/UNITS,SI/CONFIG,NRES,10000/PREP7ET,1,PLANE55 ET,2,SOLID70
!316LMPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7850 MPTEMP,,,,,,,, MPTEMP,1,20 MPTEMP,2,300MPTEMP,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,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
!45MPTEMP,,,,,,,, MPTEMP,1,20,250,500,750,1000,1500,1700,2500MPDATA,C,2,1,460,480,530,675,670,660,780,820MPDATA,KXX,2,1,50,47,40,27,30,35,140,142MPDATA,DENS,2,1,7820,7700,7610,7550,7490,7350,7300,7090MPTEMP,,,,,,,,MPTEMP,1,0,1533,1595,1670MPDATA,ENTH,2,1,0,7.5e,9.6e9,10.5e10
!*****************定义常量**************WidthBase=0.02 !基底宽度HeightBase=0.008 !基底高度
RClad=0.002 !熔覆半径
LENGTH=0.04 !长度LNUMBER=80 !划分的数量LSIZE=LENGTH/LNUMBER !每一小段
P=800 !激光功率R=0.002 !激光光斑半径pi=3.14159265 V=0.005 !扫描速度Factor=0.5 !吸收因子tinc=LSIZE/V !每小段的时间间隔
Temp=25 !环境温度CoffConv=30 !对流换热系数QMAX=6*Factor*P/pi/SQRT(pi)/R/R/R!*****************建立模型,生成网格**************K,,0,0,0K,,RClad,0,0K,,WidthBase,0,0K,,0,-HeightBase,0K,,RClad,-HeightBase,0K,,WidthBase,-HeightBase,0 A,1,2,5,4A,2,3,6,5
/PNUM,KP,1 /PNUM,LINE,1/PNUM,AREA,1/PNUM,VOLU,0/PNUM,NODE,0/PNUM,TABN,0/PNUM,SVAL,0/NUMBER,0 !* /PNUM,ELEM,0/REPLOT !* GPLOT CYL4,0,0,0,0,RClad,90 AGLUE,All
TYPE,1MAT,1MSHKEY,1ESIZE,0.0004,AMESH,4
lesize,4,0.0004,,,0.2 lesize,2,0.0004,,,5 lesize,6,0.0004,,,5 lesize,5,0.0008,,,8 lesize,7,0.0008,,,0.125
TYPE,1MAT,2MSHKEY,1AMESH,1 TYPE,1MAT,2MSHKEY,1AMESH,2 TYPE, 2 EXTOPT,ESIZE,LNUMBER,0, EXTOPT,ACLEAR,1 !* EXTOPT,ATTR,0,0,0 MAT,1 REAL,_Z4ESYS,0 !* !* VOFFST,4,LENGTH, ,
TYPE, 2 EXTOPT,ESIZE,LNUMBER,0, EXTOPT,ACLEAR,1 !* EXTOPT,ATTR,0,0,0 MAT,2 REAL,_Z4ESYS,0 !* !* VOFFST,1,-LENGTH,, VOFFST,2,-LENGTH,, EPLOT
NUMMRG,ALL, , , ,LOWFINISH
/solu!==!基底边界条件、初始条件NSEL, S, LOC, Y, -HeightBase, 0 !基底初始温度IC,ALL,TEMP,25ALLSEL,ALL
NSEL,S,LOC,Y,-HeightBase,0 !基底侧面,换热边界条件NSEL,R,LOC,x,WidthBaseSF,ALL,CONV,CoffConv,TempALLSEL,ALL
NSEL,S,LOC,Y,-HeightBase,0 !基底左端面,换热边界条件NSEL,R,LOC,Z,0SF,ALL,CONV,CoffConv,TempALLSEL,ALL
NSEL,S,LOC,Y,-HeightBase,0 !基底右端面,换热边界条件NSEL,R,LOC,Z,LENGTHSF,ALL,CONV,CoffConv,TempALLSEL,ALL
NSEL,S,LOC,Y,0 !基底上表面,换热边界条件NSEL,R,LOC,X,RClad,WidthBaseSF,ALL,CONV,CoffConv,TempALLSEL,ALLFINISH!*****************载荷设置选项**************/soluANTYPE,TRANSTRNOPT,FULLNROPT,FULL,,ONNLGEOM,OFF
PRED,ONTIMINT,ONTINTP,0.005,,,1,0.5,0.2TREF,25TOFFST,273
!*****************杀死单元**************VSEL,S,MAT,,1ESLV,SII=0*GET,WENUM,ELEM,,COUNT*GET,emax,ELEM,,NUM,MAX*GET,emin,ELEM,,NUM,MIN*DIM,ENUM,,WENUM*DO,I,emin,emax*IF,ESEL(I),EQ,1,THENII=II+1ENUM(II)=I*ENDIF*ENDDOESEL,ALLVSEL,S,MAT,,1ESLV,SEKILL,ALLESEL,S,LIVEEPLOT!************定义数组维数 建立局部坐标系**********
MAX_X=1+LNUMBER*DO,IM,1,MAX_X,1local,12,0,0,RClad,(IM-1)*LSIZE,,,,,,
*DO,I,1,WENUM VSEL,S,MAT,,1 ESLV,S *GET,RX,ELEM,ENUM(I),CENT,X *GET,RY,ELEM,ENUM(I),CENT,Y *GET,RZ,ELEM,ENUM(I),CENT,Z RR=ABS(SQRT(RX*RX+RY*RY+RZ*RZ)) *IF,RR,LE,R,THEN QR=Qmax*exp((-3)*(RX*RX+RY*RY+RZ*RZ)/R/R) EALIVE,ENUM(I) ESEL,S,ELEM,,ENUM(I) BFE,ALL,HGEN,,QR *ENDIF *ENDDO TM=IM*tinc TIME,TM TIMINT,ON AUTOTS,ON KBC,1 DELTIM,0.005,0.001,0.005 ALLSEL,ALL SOLVE BFDELE,ALL,HGEN*ENDDO
|