第五章 3D静态磁场分析(标量法)
- 格式:doc
- 大小:358.50 KB
- 文档页数:31
ANSYS电磁场分析指南(共17章)ANSYS电磁场分析指南第一章磁场分析概述:ANSYS电磁场分析指南第二章 2-D静态磁场分析:ANSYS电磁场分析指南第三章2-D谐波(AC)磁场分析:ANSYS电磁场分析指南第四章2-D瞬态磁场分析:ANSYS电磁场分析指南第五章3-D静态磁场分析(标量法):ANSYS电磁场分析指南第六章3-D静态磁场分析(棱边元方法):ANSYS电磁场分析指南第七章3-D谐波磁场分析(棱边单元法):ANSYS电磁场分析指南第八章3-D瞬态磁场分析(棱边单元法):ANSYS电磁场分析指南第九章 3-D静态、谐波和瞬态分析(节点法):ANSYS电磁场分析指南第十章高频电磁场分析:ANSYS电磁场分析指南第十一章磁宏:ANSYS电磁场分析指南第十二章远场单元:ANSYS电磁场分析指南第十三章电场分析:ANSYS电磁场分析指南第十四章静电场分析(h方法):ANSYS电磁场分析指南第十五章静电场分析(P方法):ANSYS电磁场分析指南第十六章电路分析:ANSYS电磁场分析指南第十七章其它分析选项和求解方法:第一章磁场分析概述1.1磁场分析对象利用ANSYS/Emag或ANSYS/Multiphysics模块中的电磁场分析功能,ANSYS可分析计算下列的设备中的电磁场,如:·电力发电机·磁带及磁盘驱动器·变压器·波导·螺线管传动器·谐振腔·电动机·连接器·磁成像系统·天线辐射·图像显示设备传感器·滤波器·回旋加速器在一般电磁场分析中关心的典型的物理量为:·磁通密度·能量损耗·磁场强度·磁漏·磁力及磁矩· S-参数·阻抗·品质因子Q·电感·回波损耗·涡流·本征频率存在电流、永磁体和外加场都会激励起需要分析的磁场。
=====【热力耦合分析单元简介】======SOLID5-三维耦合场实体具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。
本单元由8个节点定义,每个节点有6个自由度。
在静态磁场分析中,可以使用标量势公式(对于简化的RSP,微分的DSP,通用的GSP)。
在结构和压电分析中,具有大变形的应力钢化功能。
与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。
INFIN9-二维无限边界用于模拟一个二维无界问题的开放边界。
具有两个节点,每个节点上带有磁向量势或温度自由度。
所依附的单元类型可以为PLANE13和PLANE53磁单元,或PLANE55和PLANE77和PLANE35热单元。
使用磁自由度(AZ)时,分析可以是线性的也可以是非线性的,静态的或动态的。
使用热自由度时,只能进行线性稳态分析。
PLANE13-二维耦合场实体具有二维磁场、温度场、电场和结构场之间有限耦合的功能。
由4个节点定义,每个节点可达到4个自由度。
具有非线性磁场功能,可用于模拟B-H曲线和永久磁铁去磁曲线。
具有大变形和应力钢化功能。
当用于纯结构分析时,具有大变形功能,相似的耦合场单元有SOLID5、SOLID98和SOLID62。
LINK31-辐射线单元用于模拟空间两点间辐射热流率的单轴单元。
每个节点有一个自由度。
可用于二维(平面或轴对称)或三维的、稳态的或瞬态的热分析问题。
允许形状因子和面积分别乘以温度的经验公式是有效的。
发射率可与温度相关。
如果包含热辐射单元的模型还需要进行结构分析,辐射单元应当被一个等效的或(空)结构单元所代替。
LINK32-二维传导杆用于两节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
可用于二维(平面或轴对称)稳态或瞬态的热分析问题。
如果包含热传导杆单元的模型还需进行结构分析,该单元可被一个等效的结构单元所代替。
LINK33-三维传导杆用于节点间热传导的单轴单元。
该单元每个节点只有一个温度自由度。
ANSYS电磁场分析指南(共17章)ANSYS电磁场分析指南第一章磁场分析概述:ANSYS电磁场分析指南第二章2-D静态磁场分析:ANSYS电磁场分析指南第三章2-D谐波(AC)磁场分析:ANSYS电磁场分析指南第四章2-D瞬态磁场分析:ANSYS电磁场分析指南第五章3-D静态磁场分析(标量法):ANSYS电磁场分析指南第六章3-D静态磁场分析(棱边元方法):ANSYS电磁场分析指南第七章3-D谐波磁场分析(棱边单元法):ANSYS电磁场分析指南第八章3-D瞬态磁场分析(棱边单元法):ANSYS电磁场分析指南第九章3-D静态、谐波和瞬态分析(节点法):ANSYS电磁场分析指南第十章高频电磁场分析:ANSYS电磁场分析指南第十一章磁宏:ANSYS电磁场分析指南第十二章远场单元:ANSYS电磁场分析指南第十三章电场分析:ANSYS电磁场分析指南第十四章静电场分析(h方法):ANSYS电磁场分析指南第十五章静电场分析(P方法):ANSYS电磁场分析指南第十六章电路分析:ANSYS电磁场分析指南第十七章其它分析选项和求解方法:第一章磁场分析概述1.1磁场分析对象利用ANSYS/Emag或ANSYS/Multiphysics模块中的电磁场分析功能,ANSYS可分析计算下列的设备中的电磁场,如:·电力发电机·磁带及磁盘驱动器·变压器·波导·螺线管传动器·谐振腔·电动机·连接器·磁成像系统·天线辐射·图像显示设备传感器·滤波器·回旋加速器在一般电磁场分析中关心的典型的物理量为:·磁通密度·能量损耗·磁场强度·磁漏·磁力及磁矩·S-参数·阻抗·品质因子Q·电感·回波损耗·涡流·本征频率存在电流、永磁体和外加场都会激励起需要分析的磁场。
第39卷第2期原子能科学技术Vol.39,No.2 2005年3月Atomic Energy Science and TechnologyMar.2005100Me V 紧凑型回旋加速器主磁铁的三维磁场有限元分析技术钟俊晴,张天爵,杨建俊,储诚节(中国原子能科学研究院核技术应用研究所,北京 102413)摘要:100MeV 紧凑型回旋加速器主磁铁的几何结构十分复杂,但为了形成加速器束流动力学所要求的磁场分布,本文对初步设计的磁铁进行必要的简化。
综合采用各种适当的三维有限元网格剖分技术,对该磁铁的磁场进行数值分析,计算精度满足加速器物理设计的要求。
关键词:紧凑型回旋加速器;有限元法;剖分;双标量位法中图分类号:TL542 文献标识码:A 文章编号:100026931(2005)022*******3D Finite Element Analysis Method for Main Magnet Design of 100Me V Compact CyclotronZHON G J un 2qing ,ZHAN G Tian 2jue ,YAN G Jian 2jun ,C HU Cheng 2jie(China I nstitute of A tomic Energ y ,P.O.B ox 27523,B ei j ing 102413,China )Abstract : To get t he magnetic field distribution for beam dynamics calculation in t he 100MeV co mpact cyclot ron feat ured wit h a complex geomet ric magnet ,we simplified t he magnet struct ure and used various 3D finite element meshing met hods in t he initial design stage.After t he result s comparison by different meshing met hods ,we analysed numerically t he t hree dimentional magnetic field and found t hat t his 3D finite element calculation p rocess could meet t he requirement of t he p hysics design for 100MeV com 2pact cyclot ron.K ey w ords :compact cyclot ron ;finite element met hod ;meshing ;double scalar potential收稿日期:2004211225基金项目:国家自然科学基金资助项目(10125218)作者简介:钟俊晴(1978—),男,江西瑞金人,实习研究员,加速器专业 串列加速器升级工程主要由以下几部分组成:100MeV 回旋加速器、在线同位素分离器、超导直线增能器以及现有HI 213串列加速器注入器的改造。
第五章3-D静态磁场分析(标量法)5.1 在3-D静态磁场分析(标量法)中要用到的单元表1三维实体单元:单元维数形状或特性自由度SOLID5 3-D 六面体,8个节点每节点6个:位移、电势、磁标量位或温度SOLID96 3-D 六面体,8个节点磁标量位SOLID98 3-D 四面体,10个节点位移、电势、磁标量位、温度表2三维界面单元单元维数形状或特性自由度INTER115 3-D 四边形,4个节点磁标量位,磁矢量位表3三维连接单元单元维数形状或特性自由度3D杆状( Bar)、弧状(Arc)、线圈SOURC363个节点无(Coil)基元表4三维远场单元单元维数形状或特性自由度四边形,4个节点;INFIN47 3-D磁标量位、温度或三边形,3个节点INFIN111 3-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.3 3-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: Main Menu>Preprocessor>Material Props>Electromag Units ·根据用户设定的单位制,自由空间的相对导磁率将自动设定:在MKS单位制中,,或者根据用命令EMUNIT来设定一个值。
5.3.2设置GUI菜单过滤如果你是通过GUI路径来运行ANSYS,当ANSYS被激活后第一件要做的事情是选择菜单路径:Main Menu>Preferences,在对话框出现后,选择Magnetic-Nodal。
因为ANSYS会根据你选择的参数来对单元进行过滤,选择Magnetic-Nodal以确保能够使用用于3D静态磁场分析的单元。
5.3.3定义材料属性分析模型可有一种或多种材料区域:空气、导磁材料、导电区和永磁体。
每种类型的材料区具有所要求的材料性质。
ANSYS材料库自身带有几种磁性材料,可以直接把这些材料性质读入数据库,不用再手动的逐点输入。
如有必要,可对它们进行修改,以便与所分析的课题相匹配。
在ANSYS材料库中定义的磁性材料如下:材料材料性质文件Copper(铜)emag Copper. SI_MPLM3 steel(钢)emag M3. SI_MPLM54 steel(钢)emag M54. SI_MPLSA1010 steel(钢)emag Sa1010. SI_MPLCarpenter steel(硅钢)emag Silicon. SI_MPLIron Cobalt Vanadium steel(铁-钴-钒-钢)emag Vanad. SI_MPL该表中铜的材料性质定义有与温度有关的电阻率和相对导磁率,所有其他材料的性质均定义为B-H曲线。
对于列表中的材料,在ANSYS材料库内定义的都是典型性质,而且已外推到整个高饱和区。
你所需的实际材料值可能与ANSYS材料库提供值有所不同,因此,必要时可修正所用ANSYS材料库文件以满足用户所需。
5.3.3.1访问材料库文件:下面介绍读写材料库文件的基本过程。
详细参见《ANSYS入门指南》和《ANSYS基本过程手册》。
读材料库文件,进行以下操作:1. 如果你还没有定义好单位制,用/UNITS命令定义。
注意:缺省单位制为MKS,GUI列表只列出当前被激活单位制的材料库文件。
2. 定义材料库文件所在的路径。
(你需要知道系统管理员放置材料库文件的路径)命令:/MPLIB,read,pathdataGUI: Main Menu>Preprocessor>Material Props>Material Library>Library Path3. 将材料库文件读入到数据库中。
命令:MPREAD,filename,,,LIBGUI:Main Menu>Preprocessor>Material Props>Material Library>Import Library Main Menu>Preprocessor>Loads>-Load Step Opts-Other>Change MatProps>Material Library>Import Library写材料库文件,进行以下操作:1. 用MP命令或菜单Main Menu>Preprocessor>Material Props>Isotropic编辑材料性质定义,然后将改后的材料特性写回到材料库文件当中去。
2.在前处理器中执行下列命令:命令:MPWRITE,filename,,,LIB,MATGUI:Main Menu>Preprocessor>Material Props>Material Library>Export Library5.3.4定义材料属性和实常数的一般原则下面讲述关于设置物理模型区域的一般原则。
在“2-D谐波(AC)分析”一章中详细描述了2-D模型中需要设定的一些特殊区域。
5.3.4.1空气:说明相对磁导率为1.0。
命令:MP,murxGUI: Main Menu>Preprocessor>Material Props>Material Models > Electromagnetics > Relative Permeability >Constant5.3.4.2自由空间导磁材料区:说明B-H曲线,可以从库中读出,也可以输入自己定义的B-H曲线:命令:MPREAD,filename,…GUI:Main Menu>Preprocessor>Material Props>Material Library>Import Library 命令:TB,TBPTGUI:Main Menu>Preprocessor>Material Props>Material Models>Electromagnetics>BH Curve如果该材料是线性的,那么只需说明相对磁导率(对于均匀或者各向异性介质)。
如果定义B-H曲线,应满足下列准则,以保证准确地模拟材料属性:*输入B-H曲线必须要遵守的规则:1. B与H要一一对应,且应B随H是单调递增,如第2章中图1所示。
B-H曲线缺省通过原点,即0.0点不输入。
用下面的命令验证B-H曲线:命令:TBPLOTGUI: Main Menu>Preprocessor>Material Props>Material Models>Electromagnetics>BH Curve2. ANSYS计算的µ-H曲线应该是光滑连续的,可用TBPLOT命令画出µ-H曲线来进行验证(参见图2(b))。
B-H曲线应覆盖材料的全部工作范围。
如果需要超出B-H曲线的点,程序按固定斜率进行外推处理,固定斜率应等于或者略大于µr,你可以按照如下方式改变X-轴的范围,并用TBPLOT命令画图来观察其外推情况。
命令:/XRANGEGUI: Utility Menu>PlotCtrls>Style>Graphs对于永磁体和非线性各向异性材料,材料特性的输入与2D情况一样。
5.3.5建立模型《ANSYS建模和分网指南》对3-D标标量位方法进行磁场分析的建模有详细描述,但是对于电流源的处理有一些特殊考虑。
5.3.5.1建立电流传导区可以用基元模拟电流传导区域,不需要材料性质。
在3D标量位方法分析中,电流源不是有限元模型的一个组成的部分(在2D矢量位方法分析中是一个组成部分)。
只需用一个有限元哑元单元,SOURC36,来指明电流源的形状和位置。
可以在模型中的任意位置定义线圈、杆状、弧状电流源,电流源的大小和其他电流源数据可以通过哑元单元的实常数定义给出。
图3为用SOURC36原始元表示的一个电流源。
注意:即使采用半对称或四分之一对称模型,也要为整个电流源建模。
线圈和弧单元的内半径不能设置成为0。
因为SOURC36单元并不是一个真正的有限元,只能通过直接生成来定义它们,而不能通过实体建模的方式。