附件内容:
/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 |