激光熔覆温度场的一个例子
/CLEAR/BATCH
/UNITS,SI
/FILNAME,Thermal,1
/TITLE,Laser Cladding
/PREP7
LENGTH=0.060
WIDTH=0.053
HEIGHT=0.04
Rcoating=0.0015
PI=3.1415926
Rlaser=0.0015
Vlaser=0.006
Plaser=1000
Absor=0.6
ESIZEunit=0.0002
LENGTHdiv=LENGTH/ESIZEunit
RCOATINGdiv=Rcoating/ESIZEunit
XESIZEunit=0.005
YESIZEunit=0.005
XESIZEdiv=(WIDTH/2-Rcoating)/XESIZEunit
YESIZEdiv=HEIGHT/YESIZEunit
XESIZEratio=5
YESIZEratio=5
ET,1,SOLID70
MP,DENS,1,7140
MPTEMP,1,0,100,400,760,900,1100
MP,DENS,1,7140
MPDATA,KXX,1,1,28,29.5,30.7,26.2 ,23.6,22.9
MPDATA,C,1,1,500,520,622,1170,740,840
MP,DENS,2,8410
MPTEMP,1,27,127,327,527,727,927
MPDATA,KXX,2,1,84,74,61,62,66,70
MPDATA,C,2,1,444,485,592,530,562,594
BLOCK,0,Rcoating,-HEIGHT,0,0,LENGTH
BLOCK,Rcoating,WIDTH/2,-HEIGHT,0,0,LENGTH
CYLIND,Rcoating,,0,LENGTH,0,90
VGLUE,1,2,3
LESIZE,11,,,LENGTHdiv
LESIZE,22,,,LENGTHdiv
LESIZE,7,,,RCOATINGdiv
LESIZE,37,,,RCOATINGdiv
LESIZE,8,,,YESIZEdiv,YESIZEratio
LESIZE,1,,,YESIZEdiv,YESIZEratio
LESIZE,18,,,YESIZEdiv,1/YESIZEratio
LESIZE,15,,,YESIZEdiv,1/YESIZEratio
LESIZE,6,,,YESIZEdiv,1/YESIZEratio
LESIZE,3,,,YESIZEdiv,1/YESIZEratio
LESIZE,41,,,XESIZEdiv,XESIZEratio
LESIZE,38,,,XESIZEdiv,XESIZEratio
LESIZE,39,,,XESIZEdiv,XESIZEratio
LESIZE,40,,,XESIZEdiv,XESIZEratio
TYPE,1
MAT,1
MSHAPE,0,3D
MSHKEY,1
VMESH,1
VMESH,5
TYPE,1
MAT,2
MSHAPE,0,3D
MSHKEY,1
VMESH,4
/SOLU
ANTYPE,TRANS,NEW
TRNOPT,FULL
LUMPM,OFF
NLGEOM,OFF
OUTPR,BASIC,ALL
OUTRES,ALL,ALL
TOFFST,273
TUNIF,20
NSEL,S,LOC,X,0
DSYM,SYMM,X,,
ASEL,S,AREA,,2
ASEL,A,AREA,,23
ASEL,A,AREA,,12
ASEL,A,AREA,,25
ASEL,A,AREA,,24
ASEL,A,AREA,,3
ASEL,A,AREA,,22
ASEL,A,AREA,,1
NSLA,S,1
!选择附在面上的节点
SF,ALL,CONV,65,20
!施加对流
TINC=ESIZEunit/Vlaser
TM=LENGTH/Vlaser
MAX_X=Rcoating/ESIZEunit+1
MAX_Z=LENGTH/ESIZEunit+1
MAX_TIME=TM/TINC+1
*DIM, FLUX1,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME
*DO,K,1,MAX_TIME,1
*DO,I,1,MAX_X,1
FLUX1(I,0,K)=(I-1)*ESIZEunit
*ENDDO
*DO,J,1,MAX_Z,1
FLUX1(0,J,K)=(J-1)*ESIZEunit
*ENDDO
*DO,I,1,MAX_X,1
*DO,J,1,MAX_Z,1
Xcenter=0
Zcenter=Vlaser*(K-1)*TINC+ESIZEunit/2
DISTANCE=sqrt(((I-1)*ESIZEunit-Xcenter)**2+((J-1)*ESIZEunit-Zcenter)**2)
*IF,DISTANCE,LE,Rlaser,THEN
FLUX1(I,J,K)=0.6*Plaser/(PI*Rlaser**2)/EXP(3*DISTANCE**2/(Rlaser**2))
*ELSE
FLUX1(I,J,K)=0
*ENDIF
*ENDDO
*ENDDO
FLUX1(0,0,K)=(K-1)*TINC
*ENDDO
/SOLU
ANTYPE,TRANS
TUNIF,20
KBC,1
*SET,DELT,1
TIMINT,ON
AUTOTS,OF
TT=0.8
*DO,TM,1E-6,TT,TINC
/VIEW,1,1,1,1
/ANG,1
TIME,TM
AUTOTS,0
DELTIM, , , ,1
KBC,0
!*
TSRES,ERASE
OUTRES,ALL,ALL,
SOLVE
*ENDDO
注：这个例子有些问题，我现在想让它顺利的加载下去就可以，麻烦大家出出主意，谢谢了93089904 这个没有生死单元吗
页:
[1]