各位达人,小弟在做三维的粒子转向装置的设计,遇到了一些问题,想请教一下:我想画出x=0平面内的磁场、x=0.001平面内的磁场,可是不知道怎么画;在后处理中,加入电子并观察电子在其中的轨迹,可是老是提示我电子的质量为零(我用的是国际单位制下的电子质量)。
我采用的是用autocad建模,然后将模型输入,附件里就是建模用的sat文件
我的命令流如下:
/input,menust,tmp,'',,,,,,,,,,,,,,,,1
/GRA,POWER
/GST,ON
/PLO,INFO,3
/GRO,CURL,ON
/CPLANE,1
/REPLOT,RESIZE
WPSTYLE,,,,,,,,0
~SATIN,file,sat,,SOLIDS,0
/NOPR
/GO
ET,2,sourc36 ! Current source element
EMUNIT,MKS ! MKS units
! Define convenient parameters:
*SET,I,2 ! Current (amps)
*SET,N,700 ! Turns
*SET,S,0.039 ! Solenoid length
*SET,R,0.176 ! Solenoid radius
*SET,THK,0.112 ! Solenoid thickness
!
R,2,1,N*I,THK,S ! Real constant set 2: coil type, current
! thickness, length,
!first coil
LOCAL,11,1,0,0,0.0195, , , ,1,1, ! Global cylindrical system
N,1001,R,90 ! Nodes for the source element
N,1002,R
N,1003
TYPE,2 ! Attributes
REAL,2
E,1001,1002,1003 ! Element definition
/ESHAPE,1
/VUP,1,Z
/TRIAD, LBOT
/TYPE,1
!second coil
FLST,4,3,1,ORDE,2
FITEM,4,1001
FITEM,4,-1003
NGEN,2,1003,P51X, , ,0,0,0.18,1,
TYPE,2 ! Attributes
REAL,2
E,2004,2005,2006 ! Element definition
/ESHAPE,1
/VUP,1,Z
/TRIAD, LBOT
/TYPE,1
!third coil
FLST,4,3,1,ORDE,2
FITEM,4,1001
FITEM,4,-1003
NGEN,2,2006,P51X, , ,0,0,0.36,1,
TYPE,2 ! Attributes
REAL,2
E,3007,3008,3009 ! Element definition
/ESHAPE,1
/VIEW,1,2,1,.5
/VUP,1,Z
/TRIAD, LBOT
/TYPE,1
!fourth coil
FLST,4,3,1,ORDE,2
FITEM,4,1001
FITEM,4,-1003
NGEN,2,3009,P51X, , ,0,0,0.54,1,
TYPE,2 ! Attributes
REAL,2
E,4010,4011,4012 ! Element definition
/ESHAPE,1
/VUP,1,Z
/TRIAD, LBOT
/TYPE,1
!fifth coil
LOCAL,12,1,0,0.83,0.8495, ,90, ,1,1, ! Global cylindrical system
N,5001,R ! Nodes for the source element
N,5002,R,90
N,5003
TYPE,2 ! Attributes
REAL,2
E,5001,5002,5003 ! Element definition
/ESHAPE,1
/VUP,1,Z
/TRIAD, LBOT
/TYPE,1
!sixth coil
FLST,4,3,1,ORDE,2
FITEM,4,5001
FITEM,4,-5003
NGEN,2,1003,P51X, , ,0,0,0.18,1,
TYPE,2 ! Attributes
REAL,2
E,6004,6005,6006 ! Element definition
/ESHAPE,1
/VUP,1,Z
/TRIAD, LBOT
/TYPE,1
!seventh coil
FLST,4,3,1,ORDE,2
FITEM,4,5001
FITEM,4,-5003
NGEN,2,2006,P51X, , ,0,0,0.36,1,
TYPE,2 ! Attributes
REAL,2
E,7007,7008,7009 ! Element definition
/ESHAPE,1
/VUP,1,Z
/TRIAD, LBOT
/TYPE,1
!eighth coil
FLST,4,3,1,ORDE,2
FITEM,4,5001
FITEM,4,-5003
NGEN,2,3009,P51X, , ,0,0,0.54,1,
TYPE,2 ! Attributes
REAL,2
E,8010,8011,8012 ! Element definition
/ESHAPE,1
/VIEW,1,2,1,.5
/VUP,1,Z
/TRIAD, LBOT
/TYPE,1
EPLOT
!*
/NOPR
/PMETH,OFF,1
KEYW,PR_SET,1
KEYW,PR_STRUC,0
KEYW,PR_THERM,0
KEYW,PR_FLUID,0
KEYW,PR_ELMAG,1
KEYW,MAGNOD,1
KEYW,MAGEDG,0
KEYW,MAGHFE,0
KEYW,MAGELC,0
KEYW,PR_MULTI,0
KEYW,PR_CFD,0
/GO
!*
/COM,
/COM,Preferences for GUI filtering have been set to display:
/COM, Magnetic-Nodal
!*
SAVE
!*
ET,3,SOLID98,10
!*
mp,murx,1,1
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, 2, 3, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
SAVE
MSHKEY,0
MSHAPE,1,3d
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1
MSHKEY,0
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
MSHKEY,0
MSHAPE,1,3d
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
VMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
MSHKEY,0
MSHAPE,1,3d
FLST,2,1,6,ORDE,1
FITEM,2,1
VIMP,P51X,0
FLST,5,4,5,ORDE,3
FITEM,5,2
FITEM,5,-4
FITEM,5,7
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CMSEL,S,_Y
CMDELE,_Y
!*
!*
AREFINE,_Y1, , ,1,0,1,1
CMDELE,_Y1
!*
FLST,2,1,6,ORDE,1
FITEM,2,1
VIMP,P51X,0
eplot
FLST,2,2,5,ORDE,2
FITEM,2,5
FITEM,2,8
DA,P51X,MAG,0,0
FINISH
/SOL
MAGSOLV,2, , ,0.001,25,1
FINISH
/POST1
!*
/VSCALE,1,1,0
!
!*
PLVECT,B, , , ,VECT,ELEM,ON,0
SAVE
/REPLOT,RESIZE
FLST,2,2,8
FITEM,2,0.5,0.25,0
FITEM,2,-3.55,-0.95,0
/ZOOM,1,SCRN,0.853448,0.406483,0.950000,0.319148
FLST,2,1,8
FITEM,2,-3.4,-0.85,0
FLST,2,1,8
FITEM,2,0,0,0
/AUTO,1
/REP,FAST
FLST,2,1,8
FITEM,2,0,0,0
FLST,2,1,8
FITEM,2,0,0,0 |