ANSYS FEA 热仿真，程序不收敛，高手们来帮忙看看吧
FINISH/CLEAR,START
/FILNAME,Drilling,1
/TITLE,LaserPower Drilling
/config,nres,3000000
/PREP7
/UNIT,SI
/UIS,MSGPOP,3
/RGB,INDEX,100,100,100, 0
/RGB,INDEX, 80, 80, 80,13
/RGB,INDEX, 60, 60, 60,14
/RGB,INDEX, 0, 0, 0,15
/REPLOT
/VIEW,1,607,-748,267
/VUP,1,Z
/REP,FAST
*AFUN,DEG
!-------------------------------------------------------------------------------
ET,1,solid70
!-------------------------------------------------------------------------------
MP,DENS,1,2000
MP,KXX,1,4.4
MPTEMP,1,0,7,1000,1001,3000,3001
MPDATA,ENTH,1,1,0,5.25E3*2000,7.5E5*2000,2.55E6*2000,4.05E6*2000,17.7E6*2000
MPTEMP,7,4000,600000
MPDATA,ENTH,1,7,3E7*2000,1E10*2000
!-------------------------------------------------------------------------------
INTemp=7
TOFFSET,293
*ASK,RBase,'The Radius of the Model(mm)',0.02
*ASK,HBase,'The Hight of the Model(mm)',0.01
*ASK,MIS,'The Model Size Parameter,the Greater the Number,the Finer the Model',6
*ASK,LaserPower,'The Power ofThe Laser(W/m2)',4E7
*ASK,LengthR,'The Radius ofThe Laser',0.016
KK=4
R=LengthR
Qmax=LaserPower*2
RPower=R
Factor=0.6
*ASK,dt1,'The Time for Ten Steps',1E-3
*ASK,dt2,'The Time for Working',0.1
!-------------------------------------------------------------------------------
K,1,0,0,0
K,2,RBase,0,0
K,3,RBase,RBase,0
K,4,0,RBase,0
K,5,0,0,HBase
K,6,RBase,0,HBase
K,7,RBase,RBase,HBase
K,8,0,RBase,HBase
L,1,2
L,1,4
L,4,3
L,2,3
L,5,6
L,5,8
L,8,7
L,6,7
L,1,5
L,2,6
L,3,7
L,4,8
A,1,2,3,4
A,5,6,7,8
A,1,2,6,5
A,1,4,8,5
A,4,3,7,8
A,2,3,7,6
VA,1,3,6,5,4,2
LSEL,S,,,1,8,,
CM,L1,LINE
LESIZE,L1,,,MIS*7
CMDELE,L1
LSEL,S,,,9,12
CM,L2,LINE
LESIZE,L2,,,MIS*7*HBase/Rbase,0.8
CMDELE,L2
MSHKEY,1
VMESH,1
ALLSEL,ALL
ACLEAR,ALL
NUMMRG,ALL,,,,LOW
NUMCMP,ALL
DA,5,TEMP,INTEMP
DA,6,TEMP,INTEMP
!-------------------------------------------------------------------------------
emax=0
ESEL,ALL
*GET,emax,ELEM,,NUM,MAX
*GET,Nmax,NODE,,NUM,MAX
*DIM,ne1,,emax
*DIM,ne1y,,emax
*DIM,ne1x,,emax
*DIM,ne1z,,emax
*DIM,ne1r,,emax
*DIM,ne1flag,,emax
*DIM,ne3,,emax
*DIM,ne3r,,emax
*DIM,ne3t,,emax
*DO,i,1,emax
ne1(i)=-1
ne1x(i)=-1
ne1y(i)=-1
ne1z(i)=-1
ne1r(i)=-1
ne1flag(i)=-1
ne3(i)=-1
ne3r(i)=-1
ne3t(i)=-1
*ENDDO
*DO,i,1,emax
ne1(i)=i
*GET,ne1y(i),ELEM,ne1(i),CENT,Y
*GET,ne1x(i),ELEM,ne1(i),CENT,X
*GET,ne1z(i),ELEM,ne1(i),CENT,Z
ne1r(i)=SQRT(ne1x(i)**2+ne1y(i)**2)
*ENDDO
ESEL,ALL
minnode=NODE(0,0,0)
minelem=ENEARN(minnode)
noene=minelem
*GET,noenex,ELEM,noene,CENT,X
*GET,noeney,ELEM,noene,CENT,Y
*GET,noenez,ELEM,noene,CENT,Z
*DIM,Middle,,3
i2=1
i3=1
i4=1
Middle(1)=-1
Middle(2)=-1
Middle(3)=-1
First=-1
Firstx=-1
Firsty=-1
Firstz=-1
*DO,i,1,emax
First=NELEM(ne1(i),1)
*GET,Firstx,NODE,NELEM(ne1(i),1),LOC,X
*GET,Firsty,NODE,NELEM(ne1(i),1),LOC,Y
*GET,Firstz,NODE,NELEM(ne1(i),1),LOC,Z
*IF,Firstx,GT,Middle(1),THEN
Middle(1)=Firstx
*ENDIF
*IF,Firsty,GT,Middle(2),THEN
Middle(2)=Firsty
*ENDIF
*IF,Firstz,GT,Middle(3),THEN
Middle(3)=Firstz
*ENDIF
*IF,ne1r(i),LT,R,THEN
ne3(i3)=ne1(i)
ne3r(i3)=ne1r(i)
Cne3r=(ne3r(i3)/R)*(ne3r(i3)/R)
ne3t(i3)=QMAX/(exp(KK*Cne3r))
*IF,Firstz,EQ,0,THEN
i4=i4+1
*ENDIF
i3=i3+1
*ENDIF
*ENDDO
Lemax=i3-1
LFemax=i4-1
*DIM,ne3x,,Lemax
*DIM,ne3y,,Lemax
*DIM,ne3z,,Lemax
*DO,i,1,Lemax
*GET,ne3x(i),ELEM,ne3(i),CENT,X
*GET,ne3y(i),ELEM,ne3(i),CENT,Y
*GET,ne3z(i),ELEM,ne3(i),CENT,Z
*ENDDO
*DIM,ne5,,LFemax
*DIM,ne5x,,LFemax
*DIM,ne5y,,LFemax
*DIM,ne5z,,LFemax
*DIM,ne5r,,LFemax
*DIM,ne5t,,LFemax
*DIM,checknumber,,10
i5=1
*DO,i,1,Lemax
*IF,ne3z(i),GE,Middle(3)*(1-1E-5),THEN
ne5(i5)=ne3(i)
ne5x(i5)=ne3x(i)
ne5y(i5)=ne3y(i)
ne5r(i5)=SQRT(ne5x(i5)**2+ne5y(5)**2)
ne5z(i5)=ne3z(i)
ne5t(i5)=ne3t(i)
i5=i5+1
*ENDIF
*ENDDO
i5=i5-1
!checknumber(1)=Lemax
!checknumber(2)=LFemax
!checknumber(3)=Lemax/Lfemax
!checknumber(4)=Lemax/i5
!checknumber(5)=MIS*6*HBase/Rbase
!checknumber(6)=checknumber(4)-checknumber(3)
!checknumber(7)=middle(3)
!*status,checknumber
LFemax=i5
FINISH
/SOLU
ESEL,NONE
ALLSEL,ALL
IC,ALL,TEMP,INTemp
*DO,i5,1,LFemax
SFE,ne5(i5),6,HFLUX,,ne5t(i5),,,
*ENDDO
ANTYPE,TRANS
!EQSLV,PCG,1.0e-8,2
TRNOPT,FULL
LUMPM,OFF
KBC,1
OUTRES,,ALL
OUTPR,NSOL,LAST
AUTOTS,ON
TIMINT,ON
TINTP,0.005,,,1,0.5,0
LNSRCH,1
NEQIT,25
NCNV,2
TIME,dt2
DELTIM,0.1*dt1
SOLVE
ANSYS热分析，模型参数：密度2000Kg/m3，比热0.75J/g*K,高斯表面热源4000W/cm2，加热半径0.16cm，模型半径0.2cm，环境温度为7摄氏度。熔点1000，融化热1.8KJ/g，汽化点3000，汽化热13.7KJ/g
页:
[1]