请教一个热流密度和热生成率的加载问题

xuetx

问题是这样的，有一个长方形块体（1200km*1200km*112km），下表面给定温度1200度，
上表面给定热流密度0.05mw/m^2，块体中间给定热生成率0.0000025w/m^3。加载过后，
进行计算，后处理时，发现上表面的热流密度值根本不是我给定的值，不知道这是为什
么？难道热流密度和热生成率的加载有冲突吗？
源程序如下，希望高人给予指点，谢谢！！！

!The model simulates heat structure
/title, Heat transfer in rectangle structure
/prep7
L=1000
NUM=8      
et,1,solid70   
mp,kxx,1,3
/trlcy,volu,1
/pnum,volu,on
vplot


block,,1200*l,,1200*l,,-112*l
lesize,1,,,256/NUM-1
lesize,2,,,256/NUM-1
lesize,3,,,256/NUM-1
lesize,4,,,256/NUM-1
lesize,5,,,256/NUM-1
lesize,6,,,256/NUM-1
lesize,7,,,256/NUM-1
lesize,8,,,256/NUM-1
lesize,9,,,24/NUM
lesize,10,,,24/NUM
lesize,11,,,24/NUM
lesize,12,,,24/NUM



type,1
mat,1
mshape,0,3d
mshkey,1
vmesh,1


/solu
sfa,2,,hflux,-0.05

nsel,,loc,z,-112*L
d,all,temp,1200

/psf,hflux,2
/pbf,hgen,1

          *do,k,1,24/NUM
              nsel,s,loc,z,-112*L/(24/NUM)*k-0.1,-112*L/(24/NUM)*(k-1)+0.1
              esln,s,1
              esel,r,type,,1  
              bfe,all,hgen,,0.0000025*exp(-10*(k-1)/2)
          *enddo

time,1000

autots,1
nsubst,10,10,10
outres,all,all
allsel
solve
save
finish

200802172

　看不懂　啊