ANSYS六个例子

1Axisymmetricstaticnon-linearmagneticanalysisofprinterheadInthisexampleswewillsimulateamagneticprinterhead.Figure1showsthedevicewhichconsistsofanonferromagneticsubstrate(green),acoil(purple)throughwhichacurrentispassedandferromagneticpolepiece(red).Thepoleisbroughtintoproximityofaironroller(blue)whoseexternaldiameterisathousandtimeslargerthanthecharacteristicdimensionsofthedevice.Itisassumedthatsaturationwilloccurintheroller.ThiswillrequirethedefinitionofaB-Hdatatable.Figure1Schematicdrawingofproblem.Acylindricalsymmetrymaybeassumedasthemagneticfieldproducedbythecoilwillbeconcentratedintheregionbetweenthepoleandtheroller.Axisymmetricelementswillbeusedtomodelthisproblem.InANSYS,theYaxisistheaxisofsymmetryforallaxisymmetricproblems.Thecoilwillbedefinedashavingasingleturn,thecurrentbeingmultipiedbythenumberofrealturnstoobtaintherequiredmagneticpotentialsource.Thefirststepistodefinethevariousparametersfortheanalysis.ThesearevisibleintheANSYSinputlistingbelow.Owingtoaxialsymmetry,onlyhalfofthesensorneedbemodeled.Parametersrinr=90.0E-6!Coilinnerradiusrexr=110.0E-6!Coilouterradiusrthi=10.0e-6!Coilheightpexr=70.0E-6!Poleexteriorradiuspthi=40.0E-6!Hieghtofpolemrad=300.0E-6!Rollerexteriorradiusmthi=200.0E-6!Thicknessofrolleragap=10.0E-6!Airgapwidth2dely=mthi!Depthofsubstrateamps=20E-3!Currentincoil=20milliampsturns=827!NumberofturnsincoilSolidmodelingandmeshingThefiniteelementmodelwilluseaxisymmetricmagneticelementswith8nodes.Wedefinerelativepermeabilitiesfortheair,coilandpoleusingtheMPcommand.TheB-Hdatafortherollerisdefinedusingthetabledatacommands.Finallywedefinethegeometryandmeshthesolidmodel.ThesestepsareperformedinthepreprocessorasshownintheANSYSinputlistingbelow.Theareaofthecoiliscalculatedforlateruseduringcurrentloading.Componentsaredefinedforeachofthematerials./PREP7EMUNIT,MKS!SIelectromagneticunitsET,1,PLANE53,,,1!8nodeaxisymmetricmagneticelementPLANE53/COM,----------------------Materialproperties--------------------------------MP,MURX,1,1!Relativepermeabilityofmaterial1(air)MP,MURX,2,1!Relativepermeabilityofmaterial2(coil)MP,MURX,3,1.0E5!Relativepermeabilityofmaterial3(iron)!TB,BH,4!B-Htableformaterial4(ferromagneticmaterial)!Usearrayparameteranddo-looptodefineB-Hcurve:*DIM,H,,16!16-longHarrayH(1)=110.2,137.8,165.4,185.0,204.7,220.5,240.2,271.7!First8pointsH(9)=322.8,405.5,626.0,1405.5,3181.1,5590,8000,12000!Next8pointsB=0.0!InitializeBvalue*DO,I,1,16!Do-loopforI=1to16(instepsof1)B=B+0.1!Bvaluesincrementedby0.1TBPT,DEFI,H(I),B!B-Htable*ENDDO!Endofdo-loopTBPLOT,BH,4!PlotB-Hdatatable/COM,----------------------BuildAxisymmetricmodel---------------------------RECTNG,rinr,rexr,0,rthi!DefinerectangleforcoilRECTNG,0,pexr,0,pthi!PoleRECTNG,0,mrad,pthi+agap,pthi+agap+mthi!RollerRECTNG,0,mrad,pthi+agap,-dely!AirRECTNG,0,rexr,0,pthi+agapAOVLAP,ALL!Overlapallareas/COM,----------------------Definecomponents----------------------------------ASEL,S,LOC,Y,pthi+agap,pthi+agap+mthi!Selectareasfory(relative)0topthiCM,roller,AREA!Componentname=rollerAATT,4!Material4assignedtoroller3ASEL,S,LOC,X,0,pexr/2!Selectareasforxbetween0andpexrASEL,R,LOC,Y,0,pthi/2!Reselectareasfory(relative)0topthiCM,pole,AREA!Componentname=poleAATT,3!Material3assignedtopoleASEL,S,LOC,X,rinr,rexrASEL,R,LOC,Y,0,rthiCM,coil,AREA!Componentname=coilAATT,2!Material1assignedtocoilASUM!Calculatesareaofcoil(forlateruse)*GET,coilarea,AREA,,AREA!coilarea=areaofcoilCMSEL,S,coil!SelectcoilcomponentsCMSEL,A,pole!AdditionallyselectpolecomponentsCMSEL,A,roller!AdditionalyselectrolercomponentsASEL,INVE,AREA!InverseselectedareasCM,air,AREA!Componentname=airAATT,1!Material1assignedtoairASEL,ALL/PNUM,MAT,1/NUM,1/TITLE,MaterialattributesofareasmakingupthemodelAPLOT/COM,----------------------Meshmodel-----------------------------------------LESIZE,ALL,,,7!Use7divisionsperlineformeshCMSEL,S,POLE!SelectpoleareaAMESH,ALL!MeshalldefinedareasCMSEL,S,COILAMESH,ALLASEL,ALL!SelectallareasAMESH,ALL/TITLE,Elementsinmodel,colorsrepresentpermeabilityattributesEPLOT!PlotelementsSAVE!SavedatabaseFINISH!ExitpreprocessorBoundaryconditions,loadsandsolution2DandaxisymmetricmagneticsinANSYSuseamagneticvectorpotential(A)formulationwhosesoledegreeoffreedomistheZdirection.AconstantAZresultsinparallelfluxlines.ThisishowtheboundaryconditionsaredefinedontheexternallineoftheFEmodel.Intheinputbelow,weusesolidmodelboundaryconditiondefinitions,whichwillbetransferedtotheFEmodel.Theonlyloadingtodefineisthecurrentdensityontheelementsofthecoil.Thisisdoneusingtheelementbodyforceloadcommand(BFE)withparametersdefiningthecoilcurrent,numberofturnsandcross-sectionarea.ANSYSusestheNewton-Raphsonsolvebydefaultfornon-linearmagnetostaticproblems,sonooptionsneedbeissuedbeforesolution./SOLU!EntersolutionprocessorLSEL,S,EXT!AllexteriorlinesDL,ALL,,ASYM!Flux-parallelconditions(AZ=0)atthoselinesLSEL,ALL4CMSEL,S,coilESLA,S!AllcoilelementsBFE,ALL,JS,3,amps*turns/coilarea!Sourcecurrentdensity(JSZ)loadingALLSEL!RestoresfullsetsofallentitiesSOLVE!SolvecurrentloadstepFINISH!QuitsolutionprocessorReviewingresultsTheresultsofthea