/prep7
/UNITS,SI
/PREP7
ET,1,SOLID70
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPTEMP,10,900
MPTEMP,11,1000
MPTEMP,12,1200
MPTEMP,13,1400
MPDATA,KXX,1,,59.4
MPDATA,KXX,1,,55.3
MPDATA,KXX,1,,53.6
MPDATA,KXX,1,,50.2
MPDATA,KXX,1,,44.7
MPDATA,KXX,1,,37.1
MPDATA,KXX,1,,29.5
MPDATA,KXX,1,,21.9
MPDATA,KXX,1,,28.5
MPDATA,KXX,1,,28.8
MPDATA,KXX,1,,29.3
MPDATA,KXX,1,,29.7
MPDATA,KXX,1,,30.1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,800
MPTEMP,10,900
MPTEMP,11,1000
MPTEMP,12,1200
MPTEMP,13,1400
MPDATA,C,1,,469
MPDATA,C,1,,491
MPDATA,C,1,,519
MPDATA,C,1,,575
MPDATA,C,1,,594
MPDATA,C,1,,613
MPDATA,C,1,,632
MPDATA,C,1,,748
MPDATA,C,1,,962
MPDATA,C,1,,887
MPDATA,C,1,,774
MPDATA,C,1,,661
MPDATA,C,1,,548
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,800
MPTEMP,9,900
MPTEMP,10,1000
MPTEMP,11,1200
MPTEMP,12,1400
MPDATA,KXX,2,,51.8
MPDATA,KXX,2,,50.2
MPDATA,KXX,2,,48.6
MPDATA,KXX,2,,47.0
MPDATA,KXX,2,,42.6
MPDATA,KXX,2,,35.4
MPDATA,KXX,2,,25.4
MPDATA,KXX,2,,25.9
MPDATA,KXX,2,,26.9
MPDATA,KXX,2,,27.8
MPDATA,KXX,2,,29.7
MPDATA,KXX,2,,31.6
MPTEMP,,,,,,,,
MPTEMP,1,0
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,800
MPTEMP,9,900
MPTEMP,10,1000
MPTEMP,11,1200
MPTEMP,12,1400
MPDATA,C,2,,469
MPDATA,C,2,,494
MPDATA,C,2,,519
MPDATA,C,2,,544
MPDATA,C,2,,569
MPDATA,C,2,,594
MPDATA,C,2,,667
MPDATA,C,2,,954
MPDATA,C,2,,881
MPDATA,C,2,,801
MPDATA,C,2,,661
MPDATA,C,2,,521
MPDATA,DENS,1,,7860
MP,DENS,2,7840
BLOCK,0,-0.040,0,-0.0014,0,0.1,
BLOCK,0,0.040,0,-0.0014,0,0.1,
VGLUE,1,2
LESIZE,11,,,100
LESIZE,10,,,100
LESIZE,3,,,2
LESIZE,6,,,2
LESIZE,22,,,100
LESIZE,23,,,100
LESIZE,15,,,2
LESIZE,18,,,2
LESIZE,27,,,40,40
LESIZE,28,,,40,40
LESIZE,25,,,40,40
LESIZE,26,,,40,40
LESIZE,5,,,40,0.025
LESIZE,7,,,40,40
LESIZE,4,,,40,0.025
LESIZE,2,,,40,40
LESIZE,1,,,2
LESIZE,8,,,2
LESIZE,9,,,100
LESIZE,12,,,100
TYPE,1
MAT,2
MSHAPE,0,3D
MSHKEY,1
VMESH,3
TYPE,1
MAT,1
MSHAPE,0,3D
MSHKEY,1
VMESH,1
/SOLU
ASEL,S,AREA,,1
ASEL,S,AREA,,12
ASEL,S,AREA,,13
ASEL,S,AREA,,4
ASEL,S,AREA,,16
ASEL,S,AREA,,5
ASEL,S,AREA,,2
ASEL,S,AREA,,14
ASEL,S,AREA,,3
ASEL,S,AREA,,15
ASEL,S,AREA,,ALL
NSLA,S,1
SF,ALL,CONV,65,20
!
ALLSEL,ALL
! 施加载荷
*set,length,0.1
*set,width,0.08
*set,height,0.014
FINISH
/SOLU
antype,trans,new
timint,on
tunif,20
TOFFST,273
TRNOPT,FULL
LUMPM,OFF
NLGEOM,OFF
OUTPR,BASIC,ALL
OUTRES,ALL,ALL
*SET,pie,3.14
*SET,rb,0.1/100
*SET,V,0.04
*SET,TINC,rb/V
*SET,TM,length/V
*SET,P,1200
*SET,R,0.0002
*SET,Absorptivity,0.6
*SET,MAX_TIME,length/rb+1
*SET,MAX_Z, length/rb+1
*SET,MAX_X,R/rb+1
*DIM,FLUX2,TABLE,MAX_X,MAX_Z,MAX_TIME,X,Z,TIME
*DO,K,1,MAX_TIME,1
*DO,I,1,MAX_X,1
*SET,FLUX2(I,0,K),(I-1)*rb
*ENDDO
*DO,J,1,MAX_Z,1
*SET,FLUX2(0,J,K),(J-1)*rb
*ENDDO
*DO,I,1,MAX_X,1
*DO,J,1,MAX_Z,1
*SET,ZCENTER,V*(K-1)*TINC
*SET,XCENTER,0
*SET,Zcoord,(J-1)*rb-ZCENTER
*SET,Xcoord,(I-1)*rb-XCENTER
*SET,DISTANCE,SQRT(zcoord**2+xcoord**2)
*IF,DISTANCE,LE, R,THEN
*SET,Flux2(i,j,k),(2*Absorptivity*P*exp(-2*Distance**2/R**2))/(pie*R**2)
*ELSE
*SET,FLUX2(I,J,K),0
*ENDIF
*ENDDO
*ENDDO
*SET,FLUX2(0,0,K),(K-1)*TINC
*ENDDO
*DO,T,1E-6,TM,TINC
NSEL,S,LOC,Y,0
NSEL,R,LOC,X,-0.0002,0.0002
SF,ALL,HFLUX, %FLUX2%
ALLSEL,ALL
TIME,t
AUTOTS,OFF ! 时间步自动阶跃
NSUBST,1 !指定子步数
KBC,0 !在一个载荷步中指定多个子步数时,需要指定载荷是逐渐递增还是
TSRES,ERASE
OUTRES,ALL,ALL
SOLVE
*ENDDO
楼主能否帮我看一下我这个程序哪里出错误了,计算结果最高温度一万多度,与实际差别太大了 |