FINISH
/CLEAR,START
/units,si
/TITLE, Test for EDGE-Formulation with SOLID117
!
! Parameters
*SET,a,1e-6 ! 1st addition length of magnetic core
! *** Create Model
!
FINISH
/PREP7
ET,1,117 ! Element type #1 is a magnetic edge element, 117.
MP,MURX,1,1 ! Define relative permeability of material
MP,MURX,2,1 ! Define relative permeability of material
MP,MURX,3,500 ! Define relative permeability of material
MP,MURX,4,1 ! Define relative permeability of material
BLOCK,-(480*a),-(460*a),(63*a),(163*a),-(1000*a),(1000*a)
BLOCK,(60*a),(80*a),(63*a),(163*a),-(1000*a),(1000*a)
BLOCK,-(750*a),(750*a),(438*a),(513*a),-(100*a),(100*a)
BLOCK,-(750*a),(750*a),-(90*a),0,-(100*a),(100*a)
BLOCK,-(750*a),-(700*a),0,(438*a),-(100*a),(100*a)
BLOCK,(700*a),(750*a),0,(438*a),-(100*a),(100*a)
BLOCK,-(2000*a),(2000*a),-(1000*a),(2000*a),-(1000*a),(1000*a)
! Overlap volumes
VOVLAP,ALL
NUMCMP,VOLU
! Assign attributes to the volumes
VSEL,S,VOLU,,1
VATT,1 ! conductor
VSEL,S,VOLU,,2
VATT,2
VSEL,S,VOLU,,3,6,1 ! magnetic core
VATT,3
VSEL,S,VOLU,,7! air ambient
VATT,4
! Mesh volumes
点击MESH , 然后点击VOLUMNS , 然后点击FREE , 划分网格
ALLSEL,ALL
FINISH
/SOLU
!*** Apply current density on conductor ***
ESEL,S,MAT,,1 ! Select coil element
BFE,ALL,JS,1,,,-1e9 ! Apply current density
ESEL,S,MAT,,2 ! Select coil element
BFE,ALL,JS,2,,,1e9 ! Apply current density
!*** Flux parallel boundary condition at all extern area of air ambient ***
ASEL,S,EXT
DA,ALL,AZ,0
ALLSEL,ALL
! Solution parameters
ANTYPE,0
NROPT,AUTO,,
EQSLV,SPAR,,0,
PRECISION,0
MSAVE,0
PIVCHECK,1
/STATUS,SOLU
SOLVE
FINISH
/POST1
!*** Vectorplot of magnetic flux (BSUM) for magnetic core's and conductor's
! elements
/VSCALE,1,1,0
!
!*
PLVECT,B, , , ,VECT,ELEM,ON,0
ESEL,U,MAT,,4
/REPLOT
PLVECT,B, , , ,VECT,ELEM,ON,0
ESEL,U,MAT,,2
/REPLOT
/EFACET,1
PLNSOL, B,SUM, 0
好心人能帮忙不??把他们两个综合起来,,第2个图的导体全部用1代替, |