/batch,list
/prep7
/title, Actuator model using surface to surface contact
et,1,53 ! plane53
et,2,169 ! Target
et,3,172,7,2,,2 ! Contact, MPC option
keyopt,3,12,5 ! Bonded contact
mp,murx,1,1
tb,bh,2,,40 ! BH curve
tbpt,,355,.7
,,405,.8
,,470,.9
,,555,1.
,,673,1.1
,,836,1.2
,,1065,1.3
,,1220,1.35
,,1420,1.4
,,1720,1.45
,,2130,1.5
,,2670,1.55
,,3480,1.6
,,4500,1.65
,,5950,1.70
,,7650,1.75
,,10100,1.8
,,13000,1.85
,,15900,1.9
,,21100,1.95
,,26300,2.
,,32900,2.05
,,42700,2.1
,,61700,2.15
,,84300,2.2
,,110000,2.25
,,135000,2.3
,,200000,2.41
,,400000,2.69
,,800000,3.22
mp,murx,3,1 ! coil
! create backiron
rectng,0,.06,0,.01
rectng,0,.01,.01,.02
rectng,.05,.06,.01,.02
! create coil
rectng,.02,.04,-.005,0
rectng,.02,.04,.01,.015
aglue,all
! create lower air domain
rectng,-.02,.08,-.01,.021
aovlap,all
numcmp,area
cm,lower,area
! create pole
rectng,0,.06,.022,.032
! create upper air
rectng,-.02,.08,.021,.04
aovlap,7,8
asel,s,area,,7,9
cm,upper,area
asel,s,area,,4,6
aatt,2,1,1 ! backiron
asel,s,area,,7
aatt,2,1,1 ! pole
asel,s,area,,2,3
aatt,3,1,1 ! coil
cmsel,s,lower
esize,.002
amesh,all
cmsel,s,upper
esize,.003
amesh,all
esla,s
nsle,s
nsel,r,loc,y,.021
type,2
esurf ! mesh target elements
cmsel,s,lower
esla,s
nsle,s
nsel,r,loc,y,.021
type,3
esurf ! mesh contact elements
allsel,all
nsel,s,loc,x,-.02
nsel,a,loc,x,.08
nsel,a,loc,y,-.01
nsel,a,loc,y,.04
d,all,az,0 ! flux parallel
asel,s,area,,3
esla,s
bfe,all,js,,,,5e6 ! coil current density
asel,s,area,,2
esla,s
bfe,all,js,,,,-5e6
allsel,all
finish
/solu
magsolv
finish
这个例子,如果继续加大电流密度会怎样呢?5e6 到5e7 甚至是5e8
发现还是不断增大的!实际上衔铁早已经饱和了,为什么啊? |