求助有关加载高斯热源的问题！

落叶无情

我是初学者，在做关于加载高斯热源的问题，问题是我想在某段线上加载高斯热源，如何建模和如何加载比较好，麻烦详细一点，谢谢！

wegean

主要是这条线划分网格的时候，离热源中心的位置网格要细一些

落叶无情

我还是加载不上啊，都困扰我一个星期了

落叶无情

wegean 发表于 2014-4-17 20:13
主要是这条线划分网格的时候，离热源中心的位置网格要细一些

我已经试过加载均匀热源和高斯热源，都是加载不上啊，请问有办法吗？谢谢你了

wegean

落叶无情 发表于 2014-4-19 16:24
我已经试过加载均匀热源和高斯热源，都是加载不上啊，请问有办法吗？谢谢你了 ...

是函数没定义好吗？
先定义再读取，然后再施加的

落叶无情

wegean 发表于 2014-4-19 18:00
是函数没定义好吗？
先定义再读取，然后再施加的

我知道步骤啊，我现在初步怀疑是建模的时候出现问题，你能帮我看一下吗？
以下是我的命令流：
/NOPR   ! Suppress printing of UNDO process
/PMACRO ! Echo following commands to log   
FINISH  ! Make sure we are at BEGIN level   
/CLEAR,NOSTART  ! Clear model since no SAVE found   
! WE SUGGEST YOU REMOVE THIS LINE AND THE FOLLOWING STARTUP LINES   
/input,menust,tmp,'',,,,,,,,,,,,,,,,1   
/GRA,POWER  
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1   
/REPLOT,RESIZE  
WPSTYLE,,,,,,,,0
!*  
/NOPR   
KEYW,PR_SET,1   
KEYW,PR_STRUC,0
KEYW,PR_THERM,1
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,0
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,0
KEYW,PR_CFD,0   
/GO
!*  
!*  
/PREP7  
!*  
ET,1,PLANE55
!*  
MPTEMP,,300,800,1400,2000,2600,3253
MPDATA,KXX,1,,40,16,10.1,9.4,8.6,5.4
MPDATA,C,1,,797,973,1113,1329,1540,1660
MPDATA,ENTH,1,,8.56E7,2.46E9,5.89E9,9.825E9,1.427E10,1.97E10
MPDATA,DENS,1,,3980,3980,3980,3980,3980,3980
MPDATA,KXX,2,,395,395,325,295,230,170   
MPDATA,C,2,,380,437,550,737,968,1272
MPDATA,ENTH,2,,8.56E7,2.058E9,5.54E9,1.137E10,2.012E10,3.328E10
MPDATA,DENS,2,,8930,8930,8930,8930,8930,8930
MPDATA,KXX,3,,0.55,0.7,15,45,68,198
MPDATA,C,3,,3400,3600,3800,4000,4240,4380   
MPDATA,ENTH,3,,4.51E8,9.438E9,2.467E10,4.467E10,6.48E10,8.896E10
MPDATA,DENS,3,,1200,1200,1200,1200,1200,1200
!*  
RECTNG,0,0.0005,0,0.00006,  
RECTNG,0,0.0005,0.00006,0.000065,   
RECTNG,0,0.0005,0.000065,0.000195,  
K,13,0.0002,0.000065,,  
K,14,0.0003,6.5e-005,,  
LSTR,      13,      14  
FLST,2,3,5,ORDE,2   
FITEM,2,1   
FITEM,2,-3  
AGLUE,P51X  
NUMCMP,LINE
NUMCMP,AREA
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       3, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       2
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       2, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
CM,_Y,AREA  
ASEL, , , ,       3
CM,_Y1,AREA
CMSEL,S,_Y  
!*  
CMSEL,S,_Y1
AATT,       1, ,   1,       0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1  
!*  
/PNUM,KP,0  
/PNUM,LINE,1
/PNUM,AREA,0
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
!*  
/PNUM,ELEM,0
/REPLOT
!*  
LPLOT   
FLST,5,11,4,ORDE,2  
FITEM,5,1   
FITEM,5,-11
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1,2e-6, , , , , , ,1   
!*  
MSHKEY,0
FLST,5,3,5,ORDE,2   
FITEM,5,1   
FITEM,5,-3  
CM,_Y,AREA  
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
)/GOP    ! Resume printing after UNDO process   
)! We suggest a save at this point