定义了圆锥热源，却没有加载上去，想不明白，求救了！！！！

xpypeipei

/units,si
/prep7
lx=0.03
ly=0.03
lz=0.004
v=0.01
p=2000
ef=0.3
lsize=0.0002
pi=3.1415926
h=0.003
ri=0.0002
re=0.001
n1_fz=h/lsize+1

mptemp
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPTEMP,8,700
MPTEMP,9,730
MPTEMP,10,750  
MPTEMP,11,800
mptemp,12,1000
mptemp,13,1500
mptemp,14,1840
mptemp,15,1890
mptemp,16,2860
MPDATA,KXX,1,,51.9
MPDATA,KXX,1,,49.4   
MPDATA,KXX,1,,49  
MPDATA,KXX,1,,46.1   
MPDATA,KXX,1,,42.7  
MPDATA,KXX,1,,39.4  
MPDATA,KXX,1,,35.6   
MPDATA,KXX,1,,31.8
MPDATA,KXX,1,,30.1   
MPDATA,KXX,1,,28.9   
MPDATA,KXX,1,,26
mpdate,kxx,1,,27.2
mpdata,kxx,1,,29.7
mpdata,kxx,1,,29.7
mpdata,kxx,1,,29.7
mpdata,kxx,1,,29.7
  

MPDATA,C,1,,450   
MPDATA,C,1,,494   
MPDATA,C,1,,526   
MPDATA,C,1,,566   
MPDATA,C,1,,615   
MPDATA,C,1,,684   
MPDATA,C,1,,1139   
MPDATA,C,1,,1384  
MPDATA,C,1,,1191   
MPDATA,C,1,,931   
MPDATA,C,1,,779
MPDATA,C,1,,400
MPDATA,C,1,,847
MPDATA,C,1,,400
MPDATA,C,1,,400

MPDATA,DENS,1,,7872
MPDATA,DENS,1,,7845
MPDATA,DENS,1,,7816
MPDATA,DENS,1,,7740
MPDATA,DENS,1,,7733
MPDATA,DENS,1,,7711
MPDATA,DENS,1,,7669
MPDATA,DENS,1,,7625
MPDATA,DENS,1,,7612
MPDATA,DENS,1,,7602
MPDATA,DENS,1,,7578
MPDATA,DENS,1,,7552
MPDATA,DENS,1,,7268
MPDATA,DENS,1,,6757
MPDATA,DENS,1,,6715
MPDATA,DENS,1,,5902
mptemp
MPTEMP,1,20  
MPTEMP,2,100
MPTEMP,3,300
MPTEMP,4,600
MPTEMP,5,900
MPTEMP,6,1000
MPTEMP,7,1300
MPTEMP,8,1600
MPTEMP,9,2000
MPTEMP,10,2300
MPTEMP,11,2600
mptemp,12,3000
mptemp,13,3100
mptemp,14,3300
mptemp,15,3500

MPDATA,HF,1,,5
MPDATA,HF,1,,12
MPDATA,HF,1,,32
MPDATA,HF,1,,73
MPDATA,HF,1,,123.8
MPDATA,HF,1,,141.2
MPDATA,HF,1,,198.5
MPDATA,HF,1,,290.8
MPDATA,HF,1,,364
MPDATA,HF,1,,446
MPDATA,HF,1,,532.4
MPDATA,HF,1,,724
MPDATA,HF,1,,819
MPDATA,HF,1,,963

et,1,solid70
et,2,surf152
keyopt,2,4,1
keyopt,2,5,0
keyopt,2,8,3
block,0,lx,0,ly,0,-lz
block,0,lx,0,0.004,0,-lz
block,0,lx,0.004,0.0044,0,-lz
block,0,lx,0.0044,0.008,0,-lz
block,0,lx,0.008,0.009,0,-lz

vovlap,1,2,3,4,5
vglue,6,7,8,9,10
esize,0.0002
vmesh,6
esize,0.0005
vmesh,9
esize,0.002
vmesh,7
mshkey,0
mshape,1,3d
CM,_Y,VOLU  
VSEL, , , ,       8
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  

VMESH,_Y1   

CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2
mshkey,0
mshape,1,3d
CM,_Y,VOLU  
VSEL, , , ,       10
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  

VMESH,_Y1   

CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2
type,2
asel,s,,,31,50,1
asel,a,,,4
amesh,all
save

step=0.0002
n=lx/step
*do,kk,0,n,1
xx=kk*step
tm=(xx+step)/v

/SOLU
ANTYPE,4
TRNOPT,FULL
TUNIF,20,
OUTRES,NSOL,last
NROPT,FULL,,ON
time,tm
NSUBST,5, , ,1   !!!
KBC,1
lnsrch,1     !!!
TSRES,ERASE
TINTP,,,,1,,,,,!!!
TIMINT,0,struct
TIMINT,1,THERM
TIMINT,0,MAG
esel,s,type,,2   !第一个载荷步，有对流和热生成率
nsle,s
sf,all,conv,-1,20
allsel

*afun,deg      !!!选择锥形热源内的节点
local,12,1,xx,0,0
esel,s,type,,1
nsle,s
nsel,r,loc,z,-h,-h
nsel,r,loc,x,0,ri
nsel,r,loc,y,0,180
*do,k,1,n1_fz,1
hk=k*lsize
rk=ri+hk*(re-ri)/h
nsel,a,loc,z,-h+hk,-h+hk
nsel,r,loc,x,0,rk
nsel,r,loc,y,0,180
cm,node_f,node
*enddo
local,14,0,xx,0,0

*DEL,_FNCNAME   !!!!!定义锥形热源函数
*DEL,_FNCMTID   
*DEL,_FNC_C1
*DEL,_FNC_C2
*DEL,_FNC_C3
*DEL,_FNC_C4
*DEL,_FNC_C5
*DEL,_FNCCSYS
*del,'yzx'   
*SET,_FNCNAME,'yzx'
*DIM,_FNC_C1,,1
*DIM,_FNC_C2,,1
*DIM,_FNC_C3,,1
*DIM,_FNC_C4,,1
*DIM,_FNC_C5,,1
*SET,_FNC_C1(1),ef
*SET,_FNC_C2(1),p
*SET,_FNC_C3(1),re
*SET,_FNC_C4(1),ri
*SET,_FNC_C5(1),h
*SET,_FNCCSYS,14
! /INPUT,yzx.func,,,1   
*DIM,%_FNCNAME%,TABLE,7,30,1,,,,%_FNCCSYS%  
!   
! Begin of equation: 2*ef*p*exp(-2*({X}^2+{Y}^2)/(re+(re-ri)*{Z}/h)^2)/{PI}*re^
! 2*h   
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), %_FNC_C1(1)%
*SET,%_FNCNAME%(4,0,1), %_FNC_C2(1)%
*SET,%_FNCNAME%(5,0,1), %_FNC_C3(1)%
*SET,%_FNCNAME%(6,0,1), %_FNC_C4(1)%
*SET,%_FNCNAME%(7,0,1), %_FNC_C5(1)%
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 0, 2, 0, 0, 17
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 17
*SET,%_FNCNAME%(0,3,1),   0, -1, 0, 1, -2, 3, 18
*SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, 0, 0, -2
*SET,%_FNCNAME%(0,6,1), 0.0, -4, 0, 1, -2, 2, -3
*SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 2, 0, 0, -4
*SET,%_FNCNAME%(0,8,1), 0.0, -3, 0, 1, -4, 3, -2
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 2, 0, 0, 2  
*SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, 2, 17, -2   
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,12,1), 0.0, -5, 0, 1, 3, 17, -2   
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 1, -4, 1, -5   
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -3, 3, -2   
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, 19, 2, 20   
*SET,%_FNCNAME%(0,16,1), 0.0, -3, 0, 1, -2, 3, 4
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 1, -3, 4, 21   
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, 19, 1, -2   
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,20,1), 0.0, -5, 0, 1, -3, 17, -2  
*SET,%_FNCNAME%(0,21,1), 0.0, -2, 0, 1, -4, 4, -5   
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,23,1), 0.0, -3, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,24,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -3  
*SET,%_FNCNAME%(0,25,1), 0.0, -2, 0, 1, -3, 4, -1   
*SET,%_FNCNAME%(0,26,1), 0.0, -1, 0, 2, 0, 0, 19
*SET,%_FNCNAME%(0,27,1), 0.0, -3, 0, 1, 19, 17, -1  
*SET,%_FNCNAME%(0,28,1), 0.0, -1, 0, 1, -2, 3, -3   
*SET,%_FNCNAME%(0,29,1), 0.0, -2, 0, 1, -1, 3, 21   
*SET,%_FNCNAME%(0,30,1), 0.0, 99, 0, 1, -2, 0, 0
! End of equation: 2*ef*p*exp(-2*({X}^2+{Y}^2)/(re+(re-ri)*{Z}/h)^2)/{PI}*re^2*
! h
!-->

bf,node_f,hgen,%yzx%                !!!!加载圆锥热源
cmdele,node_f
allsel,all
solve
bfdele,all,hgen
*enddo
save

我运行后，热流很小很小，求各位帮忙指正

xpypeipei

为什么自己定义的热源分布函数没有加载上去，请高手帮忙

wonita

我也想做三维椎体的热源模型的，可是不会。

wonita

一点也梅加上。我的是室温

gehuaifeng

1# xpypeipei
楼主，问题解决了没，我也想学锥形体热源怎么编辑的，一直也看不到相关例题，希望不吝赐教，能否把你的源程序发给我一份，不胜感激，guowei0717@gmail.com