yulanqiuhui 发表于 2006-9-28 16:34:01

求助一个关于高斯热源加载的问题

这是我的命令流:
/PREP7   
!*
ET,1,SOLID70
!*
ET,2,SOLID90
!*
ET,3,SURF152
!*
ET,4,SURF152
!*
KEYOPT,3,1,0
KEYOPT,3,2,0
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,5,0
KEYOPT,3,6,0
KEYOPT,3,7,0
KEYOPT,3,8,2
KEYOPT,3,9,0
!*
KEYOPT,4,1,0
KEYOPT,4,2,0
KEYOPT,4,3,0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,6,0
KEYOPT,4,7,0
KEYOPT,4,8,2
KEYOPT,4,9,0
!*
!*
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,1760   
MPTEMP,8,1893   
MPDATA,KXX,1,,6.8   
MPDATA,KXX,1,,7.4   
MPDATA,KXX,1,,8.7   
MPDATA,KXX,1,,9.8   
MPDATA,KXX,1,,10.3
MPDATA,KXX,1,,11.8
MPDATA,KXX,1,,23.6
MPDATA,KXX,1,,23.6
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,1760   
MPTEMP,8,1893   
MPDE,KXX,1
MPDATA,KXX,1,,6.8   
MPDATA,KXX,1,,7.4   
MPDATA,KXX,1,,8.7   
MPDATA,KXX,1,,9.8   
MPDATA,KXX,1,,10.3
MPDATA,KXX,1,,11.8
MPDATA,KXX,1,,23.6
MPDATA,KXX,1,,25
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,1760   
MPTEMP,8,1893   
MPDATA,C,1,,611
MPDATA,C,1,,624
MPDATA,C,1,,653
MPDATA,C,1,,674
MPDATA,C,1,,691
MPDATA,C,1,,703
MPDATA,C,1,,3731.3
MPDATA,C,1,,3731.3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,4440
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,DENS,1
MPDATA,DENS,1,,4440
SAVE,'practicemodel3','db','D:\LANNYF~1\'   
VPLOT   
/REPLOT
BLOCK,0,0.03,0,0.06,0,0.003,
BLOCK,0,0.003,0,0.06,0,0.003,   
BLOCK,0,0.01,0,0.06,0,0.003,
VPLOT   
FLST,2,3,6,ORDE,2   
FITEM,2,1   
FITEM,2,-3
VOVLAP,P51X
/PNUM,KP,0
/PNUM,LINE,0
/PNUM,AREA,1
/PNUM,VOLU,1
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0   
!*
/PNUM,ELEM,0
/REPLOT
!*
NUMCMP,VOLU
NUMCMP,AREA
NUMCMP,LINE
NUMCMP,ELEM
NUMCMP,NODE
NUMCMP,NODE
VPLOT   
CM,_Y,VOLU
VSEL, , , ,       1
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
ESIZE,0.001,0,
MSHAPE,0,3D
MSHKEY,1
!*
!*
CM,_Y,VOLU
VSEL, , , ,       1
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
ESIZE,0.001,0,
CM,_Y,VOLU
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*
VPLOT   
CM,_Y,VOLU
VSEL, , , ,       2
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
ESIZE,0.002,0,
!*
CM,_Y,VOLU
VSEL, , , ,       2
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*
VPLOT   
CM,_Y,VOLU
VSEL, , , ,       3
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       1, ,   2,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
MSHAPE,1,3D
MSHKEY,0
!*
CM,_Y,VOLU
VSEL, , , ,       3
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*
SAVE
APLOT   
FLST,5,3,5,ORDE,3   
FITEM,5,6   
FITEM,5,10
FITEM,5,14
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT,       1, ,   3,       0,
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
FLST,2,3,5,ORDE,3   
FITEM,2,6   
FITEM,2,10
FITEM,2,14
AESIZE,P51X,0.002,
MSHAPE,0,2D
MSHKEY,1
!*
FLST,5,3,5,ORDE,3   
FITEM,5,6   
FITEM,5,10
FITEM,5,14
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y
!*
AMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*
SAVE
APLOT   
/VIEW,1,,,-1
/ANG,1
/REP,FAST   
FLST,5,3,5,ORDE,3   
FITEM,5,5   
FITEM,5,9   
FITEM,5,13
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
!*
CMSEL,S,_Y1
AATT,       1, ,   4,       0,
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
FLST,2,3,5,ORDE,3   
FITEM,2,5   
FITEM,2,9   
FITEM,2,13
AESIZE,P51X,0.002,
FLST,5,3,5,ORDE,3   
FITEM,5,5   
FITEM,5,9   
FITEM,5,13
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y
!*
AMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*
/VIEW,1,,,1
/ANG,1
/REP,FAST   
VPLOT   
!!!!!!!!!!!!!!!!!
Q=956                !gsosi2
R=0.003
L=0.06
STEP=0.002         
N=5
V=0.004
!方程加载:3*Q*exp(-3*({Xr}^2+{Yr}^2)/R/R)/{PI}/R/R
!包括2个对流方程 加载在表面效应单元上
!关于选择加载的单元和节点利用GUI实现
*DO,I,0,N
YY=I*STEP
T=(YY+STEP)/V
               
/SOLU
ANTYPE,4
TRNOPT,FULL
TUNIF,293,
OUTRES,NSOL,LAST
TIME,T
NROPT,FULL,,ON
AUTOTS,1
NSUBST,5, , ,1   !!!
KBC,1
lnsrch,1   !!!
TSRES,ERASE
TINTP,,,,1,,,,,!!!
TIMINT,0,struct
TIMINT,1,THERM
TIMINT,0,MAG

LOCAL,11,1,0,YY ,0, , , ,1,2,
Csys,11
NSEL,S,LOC,Y,-90,90      
NSEL,R,LOC,X,0,R         
NSEL,R,LOC,Z,0.003,0.003
LOCAL,13,0,0,YY,0, , , ,1,1,
Csys,13
sf, all ,hflux,%gs%             !在选择的节点上面施加定义好的热流密度函数
allsel,all
solve
*enddo
Save
但是我得到的结果图片显示即使是热源加载过的地方也显示的是热源正作用结果
(不好意思不知道怎么表达了),就是热源虽然移动了,但是看不出来是高斯热源,而是一个长条状的加热带,请高手指点一下吧,谢谢谢谢!!!
页: [1]
查看完整版本: 求助一个关于高斯热源加载的问题