【求助】这个高斯表面热源哪里错了?
不知道是哪里错了,一直加载不上,分析的时候总是全变红
/FILNAME,refenxi,0
/TITLE,lianxi
/PREP7
!*
ET,1,SOLID70
!*
ET,2,SOLID90
!*
MPTEMP,,,,,,,,
MPTEMP,1,298
MPTEMP,2,673
MPTEMP,3,1073
MPTEMP,4,1473
MPTEMP,5,1623
MPTEMP,6,1670
MPTEMP,7,1696
MPTEMP,8,1708
MPTEMP,9,1737
MPTEMP,10,2273
MPTEMP,11,2773
MPDATA,KXX,1,,14.7
MPDATA,KXX,1,,20.7
MPDATA,KXX,1,,26.1
MPDATA,KXX,1,,31.1
MPDATA,KXX,1,,32.9
MPDATA,KXX,1,,33.4
MPDATA,KXX,1,,33.1
MPDATA,KXX,1,,32.5
MPDATA,KXX,1,,32.5
MPDATA,KXX,1,,30.1
MPDATA,KXX,1,,31.4
MPDATA,DENS,1,,7910
MPDATA,DENS,1,,7750
MPDATA,DENS,1,,7580
MPDATA,DENS,1,,7390
MPDATA,DENS,1,,7300
MPDATA,DENS,1,,7260
MPDATA,DENS,1,,7210
MPDATA,DENS,1,,7170
MPDATA,DENS,1,,7010
MPDATA,DENS,1,,6550
MPDATA,DENS,1,,6180
MPDATA,C,1,,472
MPDATA,C,1,,535
MPDATA,C,1,,603
MPDATA,C,1,,679
MPDATA,C,1,,715
MPDATA,C,1,,736
MPDATA,C,1,,743
MPDATA,C,1,,756
MPDATA,C,1,,771
MPDATA,C,1,,796
MPDATA,C,1,,796
MPDATA,HF,1,,1.02
MPDATA,HF,1,,2.51
MPDATA,HF,1,,4.32
MPDATA,HF,1,,6.37
MPDATA,HF,1,,7.23
MPDATA,HF,1,,7.55
MPDATA,HF,1,,7.88
MPDATA,HF,1,,8.14
MPDATA,HF,1,,9.41
MPDATA,HF,1,,12.32
MPDATA,HF,1,,15.07
BLOCK,0,0.004,0,0.002,0,0.0028
BLOCK,0,0.004,0,0.002,0.0028,0.003
BLOCK,0,0.004,0,0.002,0.003,0.00312
FLST,2,3,6,ORDE,2
FITEM,2,1
FITEM,2,-3
VGLUE,P51X
FLST,5,2,4,ORDE,2
FITEM,5,29
FITEM,5,-30
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,5e-5, , , , , , ,1
!*
FLST,5,1,4,ORDE,1
FITEM,5,42
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,2e-005, , , , , , ,1
!*
MSHAPE,0,3D
MSHKEY,1
!*
CM,_Y,VOLU
VSEL, , , , 5
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
FLST,5,12,4,ORDE,2
FITEM,5,1
FITEM,5,-12
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1,2e-4, , , , , , ,1
!*
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
TYPE, 2
MAT, 1
REAL,
ESYS, 0
SECNUM,
!*
MSHAPE,1,3D
MSHKEY,0
!*
CM,_Y,VOLU
VSEL, , , , 4
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!*
TCHG,90,87,0
ALLSEL,ALL
/SOL
NROPT,FULL,,
ANTYPE,4
TRNOPT,FULL
TUNIF,0,
AUTOTS,ON
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'gaosi'
*SET,_FNCCSYS,0
! /INPUT,gaosi.func,,,1
*DIM,%_FNCNAME%,TABLE,6,36,1,,,,%_FNCCSYS%
!
! Begin of equation: 2*0.3*200/(3.1415926*0.0001^2)*exp(-3*(({Y}-0.0005)^2+(0.
! 0025-3000*{TIME}-{X})^2)/0.0001^2)
*SET,%_FNCNAME%(0,0,1), 0.0, -999
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 0.3, 0, 0, -1
*SET,%_FNCNAME%(0,3,1), 0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 200, 0, 0, -3
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 0.0001, 0, 0, 0
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -4, 0, 1, -1, 17, -3
*SET,%_FNCNAME%(0,9,1), 0.0, -1, 0, 3.1415926, 0, 0, -4
*SET,%_FNCNAME%(0,10,1), 0.0, -3, 0, 1, -1, 3, -4
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 1, -2, 4, -3
*SET,%_FNCNAME%(0,12,1), 0.0, -2, 0, 0, 0, 0, 0
*SET,%_FNCNAME%(0,13,1), 0.0, -3, 0, 1, 0, 0, -2
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -2, 2, -3
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 1, -4, 3, -2
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 0.0005, 0, 0, 3
*SET,%_FNCNAME%(0,18,1), 0.0, -4, 0, 1, 3, 2, -2
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,20,1), 0.0, -5, 0, 1, -4, 17, -2
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 3000, 0, 0, 1
*SET,%_FNCNAME%(0,22,1), 0.0, -4, 0, 1, -2, 3, 1
*SET,%_FNCNAME%(0,23,1), 0.0, -2, 0, 0.0025, 0, 0, -4
*SET,%_FNCNAME%(0,24,1), 0.0, -6, 0, 1, -2, 2, -4
*SET,%_FNCNAME%(0,25,1), 0.0, -2, 0, 1, -6, 2, 2
*SET,%_FNCNAME%(0,26,1), 0.0, -4, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,27,1), 0.0, -6, 0, 1, -2, 17, -4
*SET,%_FNCNAME%(0,28,1), 0.0, -2, 0, 1, -5, 1, -6
*SET,%_FNCNAME%(0,29,1), 0.0, -4, 0, 1, -3, 3, -2
*SET,%_FNCNAME%(0,30,1), 0.0, -2, 0, 0.0001, 0, 0, 0
*SET,%_FNCNAME%(0,31,1), 0.0, -3, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,32,1), 0.0, -5, 0, 1, -2, 17, -3
*SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 1, -4, 4, -5
*SET,%_FNCNAME%(0,34,1), 0.0, -2, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,35,1), 0.0, -3, 0, 1, -1, 3, -2
*SET,%_FNCNAME%(0,36,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 2*0.3*200/(3.1415926*0.0001^2)*exp(-3*(({Y}-0.0005)^2+(0.
! 0025-3000*{TIME}-{X})^2)/0.0001^2)
!-->
ESEL,ALL
ESEL,S,CENT,Z,0.00311
SFE,all,6,HFLUX,,%gaosi%,,,
ALLSEL,ALL
TIME,0.0001
NSUBST,50,,,
OUTRES,ALL,ALL,
KBC,0
ALLSEL,ALL
SOLVE 返回小木虫查看更多
这个是GUI下输入生成热源的命令流,所以各项含义难理解,还是自编命令流吧
经试算和检查,问题出在载荷部分。计算结束后,检查单元热载荷,全部为零,即没有热载荷。应该是函数定义或使用部分有问题。
请问楼主这个方程是什么意思,2*0.3*200/(3.1415926*0.0001^2)*exp(-3*(({Y}-0.0005)^2+(0. 0025-3000*{TIME}-{X})^2)/0.0001^2)本人小白,刚开始学习关于高斯热源加载的分析,