ANSYS电磁场分析指南第五章3-D静态磁场分析(标量法

第五章３-Ｄ静态磁场分析（标量法）5.1在3-D静态磁场分析（标量法）中要用到的单元表1三维实体单元：单元维数形状或特性自由度SOLID53-D六面体，8个节点每节点6个：位移、电势、磁标量位或温度SOLID963-D六面体，8个节点磁标量位SOLID983-D四面体，10个节点位移、电势、磁标量位、温度表2三维界面单元单元维数形状或特性自由度INTER1153-D四边形，4个节点磁标量位，磁矢量位表3三维连接单元单元维数形状或特性自由度SOURC363D杆状(Bar)、弧状(Arc)、线圈(Coil)基元3个节点无表4三维远场单元单元维数形状或特性自由度INFIN473-D四边形，4个节点；或三边形，3个节点磁标量位、温度INFIN1113-D六面体，8个或20个节点磁矢量位、磁标量位、电势、温度SOLID96和SOLID97是磁场分析专用单元，SOLID62、SOLID5和SOLID98更适合于耦合场求解。5.2磁标量位（MSP）法介绍在磁标量位方法中，可使用三种不同的分析方法：简化标势法(RSP)、差分标势法(DSP)和通用标势(GSP)法。·若模型中不包含铁区，或有铁区但无电流源时，用RSP法。若模型中既有铁区又有电流源时，就不能用这种方法。·若不适用RSP法，就选择DSP法或GSP法。DSP法适用于单连通铁区，GSP法适用于多连通铁区。5.2.1单连通区与多连通区单连通铁区是指不能为电流源所产生的磁通量提供闭合回路的铁区，而多连通铁区则可以构成闭合回路。参见图1(a)、(b)“连通域”。数学上，通过安培定律来判断单连通区或是多连通区，即磁场强度沿闭合回路的积分等于包围的电流（或是电动势降MMF）。因为铁的磁导率非常大，所以在单连通区域中的MMF降接近于零，几乎全部的MMF降都发生在空气隙中。但在多连通区域中，无论铁的磁导率如何，所有的MMF降都发生在铁芯中。5.33-D静态磁标势分析的步骤该分析类型与2－D静态分析的步骤基本一样：1.建立物理环境2.建模、给模型区域赋属性和分网格3.加边界条件和载荷（激励）4.用RSP、DSP或GSP方法求解5.观察结果5.3.1创建物理环境首先设置分析参数为“Magnetic-Nodal”，并给出分析题目。然后用ANSYS前处理器定义物理环境包含的项目。即单元类型、KEYOPT选项、材料特性等。3D分析的大部分过程与2D分析一致，本章下面部分介绍3D分析中要特殊注意的事项。·SOLID96单元可为模型所有的内部区域建模，包括：饱和区、永磁区和空气区（自由空间）。对于电流传导区，需用SOURC36单元来表示，关于电流传导区建模，后面有详细讲述。·对于空气单元的外层区域，推荐使用INFIN47单元（4节点边界单元）或INFIN111单元（8节点或20节点边界单元）。INFIN47单元和INFIN111单元可很好地描述磁场的远场衰减，通常比使用磁力线垂直或磁力线平行条件得到的结果更准确。二种单元中，INFIN111更精确一些。·缺省单位制使用MKS单位制(米-千克-秒国际单位制)，可用下列方式改变成其他单位制。一旦选定，所有输入数据都应该使用该单位制。为了方便建模，可以先在其他单位制系统下面建模（如毫米或英寸），然后进行缩放。用下列方式定义单位制：命令：EMUNITGUI:MainMenuPreprocessorMaterialPropsElectromagUnits·根据用户设定的单位制，自由空间的相对导磁率将自动设定：在MKS单位制中，，或者根据用命令EMUNIT来设定一个值。5.3.2设置GUI菜单过滤如果你是通过GUI路径来运行ANSYS，当ANSYS被激活后第一件要做的事情是选择菜单路径：MainMenuPreferences，在对话框出现后，选择Magnetic-Nodal。因为ANSYS会根据你选择的参数来对单元进行过滤，选择Magnetic-Nodal以确保能够使用用于3D静态磁场分析的单元。5.3.3定义材料属性分析模型可有一种或多种材料区域：空气、导磁材料、导电区和永磁体。每种类型的材料区具有所要求的材料性质。ANSYS材料库自身带有几种磁性材料，可以直接把这些材料性质读入数据库，不用再手动的逐点输入。如有必要，可对它们进行修改，以便与所分析的课题相匹配。在ANSYS材料库中定义的磁性材料如下：材料材料性质文件Copper（铜）emagCopper.SI_MPLM3steel（钢）emagM3.SI_MPLM54steel（钢）emagM54.SI_MPLSA1010steel（钢）emagSa1010.SI_MPLCarpentersteel（硅钢）emagSilicon.SI_MPLIronCobaltVanadiumsteel（铁－钴－钒－钢）emagVanad.SI_MPL该表中铜的材料性质定义有与温度有关的电阻率和相对导磁率，所有其他材料的性质均定义为B－H曲线。对于列表中的材料，在ANSYS材料库内定义的都是典型性质，而且已外推到整个高饱和区。你所需的实际材料值可能与ANSYS材料库提供值有所不同，因此，必要时可修正所用ANSYS材料库文件以满足用户所需。5.3.3.1访问材料库文件：下面介绍读写材料库文件的基本过程。详细参见《ANSYS入门指南》和《ANSYS基本过程手册》。读材料库文件，进行以下操作：1.如果你还没有定义好单位制，用/UNITS命令定义。注意：缺省单位制为MKS,GUI列表只列出当前被激活单位制的材料库文件。2.定义材料库文件所在的路径。（你需要知道系统管理员放置材料库文件的路径）命令：/MPLIB,read,pathdataGUI:MainMenuPreprocessorMaterialPropsMaterialLibraryLibraryPath3.将材料库文件读入到数据库中。命令：MPREAD,filename,,,LIBGUI:MainMenuPreprocessorMaterialPropsMaterialLibraryImportLibraryMainMenuPreprocessorLoads-LoadStepOpts-OtherChangeMatPropsMaterialLibraryImportLibrary写材料库文件，进行以下操作：1.用MP命令或菜单MainMenuPreprocessorMaterialPropsIsotropic编辑材料性质定义，然后将改后的材料特性写回到材料库文件当中去。2.在前处理器中执行下列命令：命令：MPWRITE,filename,,,LIB,MATGUI:MainMenuPreprocessorMaterialPropsMaterialLibraryExportLibrary5.3.4定义材料属性和实常数的一般原则下面讲述关于设置物理模型区域的一般原则。在“2-D谐波（AC）分析”一章中详细描述了2-D模型中需要设定的一些特殊区域。5.3.4.1空气：说明相对磁导率为1.0。命令：MP,murxGUI:MainMenuPreprocessorMaterialPropsMaterialModelsElectromagneticsRelativePermeabilityConstant5.3.4.2自由空间导磁材料区：说明B-H曲线，可以从库中读出，也可以输入自己定义的B-H曲线：命令：MPREAD,filename,…GUI:MainMenuPreprocessorMaterialPropsMaterialLibraryImportLibrary命令：TB,TBPTGUI:MainMenuPreprocessorMaterialPropsMaterialModelsElectromagneticsBHCurve如果该材料是线性的，那么只需说明相对磁导率（对于均匀或者各向异性介质）。如果定义B-H曲线，应满足下列准则，以保证准确地模拟材料属性：*输入B-H曲线必须要遵守的规则：1.B与H要一一对应，且应B随H是单调递增，如第2章中图1所示。B-H曲线缺省通过原点，即0.0点不输入。用下面的命令验证B-H曲线：命令：TBPLOTGUI:MainMenuPreprocessorMaterialPropsMaterialModelsElectromagneticsBHCurve2.ANSYS计算的µ-H曲线应该是光滑连续的，可用TBPLOT命令画出µ-H曲线来进行验证（参见图2（b））。B-H曲线应覆盖材料的全部工作范围。如果需要超出B-H曲线的点，程序按固定斜率进行外推处理，固定斜率应等于或者略大于µr，你可以按照如下方式改变X-轴的范围，并用TBPLOT命令画图来观察其外推情况。命令：/XRANGEGUI:UtilityMenuPlotCtrlsStyleGraphs对于永磁体和非线性各向异性材料，材料特性的输入与2D情况一样。5.3.5建立模型《ANSYS建模和分网指南》对3-D标标量位方法进行磁场分析的建模有详细描述，但是对于电流源的处理有一些特殊考虑。5.3.5.1建立电流传导区可以用基元模拟电流传导区域，不需要材料性质。在3D标量位方法分析中，电流源不是有限元模型的一个组成的部分（在2D矢量位方法分析中是一个组成部分）。只需用一个有限元哑元单元，SOURC36，来指明电流源的形状和位置。可以在模型中的任意位置定义线圈、杆状、弧状电流源，电流源的大小和其他电流源数据可以通过哑元单元的实常数定义给出。图3为用SOURC36原始元表示的一个电流源。注意：即使采用半对称或四分之一对称模型，也要为整个电流源建模。线圈和弧单元的内半径不能设置成为0。因为SOURC36单元并不是一个真正的有限元，只能通过直接生成来定义它们，而不能通过实体建模的方式。命令：NGUI:MainMenuPreprocessorCreateNodesInActiveCSMainMenuPreprocessorCreateNodesOnWorkingPlane命令：EGUI:MainmenuPreprocessorCreateElementsThruNodes命令：EGENGUI:MainMenuPreprocessorCopyAutoNumbered直接生成的电流源单元在屏幕上是不显示的，可通过以下命令显示：命令：/ESHAPE,EPLOTGUI:UtilityMenuPlotCtrlsStyleSizeandShapeUtilityMenuPlotElements下面是一个定义电流源的命令流实例：/PREP7ET,2,36!CurrentsourceelementEMUNIT,MKS!MKSunits!Defineconvenientparameters:I=0.025!Current(amps)N=300!TurnsS=0.04!SolenoidlengthR=0.01!SolenoidradiusTHK=0.002!Solenoidthickness!R=2,1,N*1,THK,S!Realconstantset2:coiltype,current!thickness,length,CSYS,1!GlobalcylindricalsystemN,1001,R!NodesforthesourceelementN,1002,R,90N,1003TYPE,2!AttributesREAL,2E,1001,1002,1003!Elementdefinition/ESHAPE,1/VIEW,1,2,1,.5/VUP,1,Z/TRIAD,LBOT/TYPE,1,HIDPEPLOT关于下列命令的详细描述，参见《ANSYS命令手册》：ET,EMUNIT,R,CSYS,N,TYPE,REAL,E,/ESHAPE,/VIEW,/VUP,/TRIAD,和/TYPE5.3.6创建3D“