ansys电磁场仿真分析教程.ppt
- 文档编号:11476849
- 上传时间:2023-06-01
- 格式:PPT
- 页数:427
- 大小:7.20MB
ansys电磁场仿真分析教程.ppt
《ansys电磁场仿真分析教程.ppt》由会员分享,可在线阅读,更多相关《ansys电磁场仿真分析教程.ppt(427页珍藏版)》请在冰点文库上搜索。
1-1,目录,第一章电磁场仿真简介.第二章二维静态分析第1节.第2节.第3节.第4节.第5节.第三章二维谐波和瞬态分析第1节.第2节.第四章三维电磁场分析第1节.第2节.第3节.第4节.第5节.第五章耦合场分析概况.,1-42.1-12.2-12.3-12.4-12.5-13.1-13.2-14.1-14.2-14.3-14.4-14.5-15-1,1-2,第一章,教程综述,1-3,ANSYS/EMAG能用于模拟工业电磁装置电磁装置当然是3维,但可简化为2维模型。
模拟可考虑为:
稳态交流(谐波)时变瞬态阶跃电压PWM(脉宽调制)(PulseWidthModulation)任意,1-4,利用轴对称衔铁和平面定子设计致动器的一个实例衔铁旋转衔铁气隙可变化完整模型由2个独立部件组成衔铁模块定子模块,执行:
solen3d.avi看动画,1-5,模拟过程概述,利用如下方式观察装置2D与3D平面与轴对称利用轴对称平面简化模型定义物理区域空气,铁,永磁体等等绞线圈,块导体短路,开路为每个物理区定义材料导磁率(常数或非线性)电阻率矫顽磁力,剩余磁感应,衔铁,线圈,锭子,实体模型,1-6,建实体模型给模型赋予属性以模拟物理区赋予边界条件线圈激励外部边界开放边界实体模型划分网格加补充约束条件(如果有必要)周期性边界条件连接不同网格,有限元网格,1-7,进行模拟观察结果某指定时刻整个时间历程后处理磁力线力力矩损耗MMF(磁动势)电感特定需要,1-8,模拟由3个区域组成衔铁区:
导磁材料导磁率为常数(即线性材料)线圈区:
线圈可视为均匀材料.空气区:
自由空间(r=1).,衔铁,线圈,1-9,性质柱体:
r=1000线圈:
r=1匝数:
2000(整个线圈)空气:
r=1,激励线圈励磁为直流电流:
2安培,单位(mm),衔铁,Coil,长度=35,Y,X,模型轴对称,材料号2,材料号3,1-10,建模设置电磁学预选项(过滤器)对各物理区定义单元类型定义材料性质对每个物理区定义实体模型铁芯线圈空气给各物理区赋材料属性加边界条件,1-11,设置预选过滤掉其它应用的菜单Mainmenupreferences,选择OK,1-12,定义所有物理区的单元类型为PLANE53PreprocessorElementtypeAdd/Edit/Delete选择Add选择磁矢量和8节点53号单元选择OK,1-13,模拟模型的轴对称形状选择Options(选项)Elementbehavior(单元行为)选择Axisymmetric(轴对称)选择OK,1-14,定义材料PreprocessorMaterialPropsIsotropic,定义空气为1号材料(MURX=1),选择Apply(自动循环地定义下一个材料号),选择OK,1-15,定义衔铁为2号材料,选择OK,选择Apply(自动循环地选择下一个材料号),1-16,定义线圈为3号材料(自由空间导磁率,MURX=1),选择OK,选择OK(退出材料数据输入菜单),1-17,建立衔铁面PreprocessorCreateRectangleByDimensions,选择Apply(重复显示和输入)建立线圈面,选择Apply,利用TAB键移动输入窗口,1-18,建立空气面,选择OK,到了这步,建立了全部平面,但它们还没有连接起来.,衔铁,线圈,1-19,用Overlap迫使全部平面连接在一起PreprocessorOperateOverlapAreas按PickAll,现在这些平面被连接了,因此当生成单元时,各区域将共享区域边界上节点,这种操作后,原先平面被删除,而新的平面被重新编号,1-20,这些平面要求与物理区和材料联系起来Preprocessor-Attributes-DefinePickedAreas用鼠标点取衔铁平面选择OK(在选取框内)材料号窗口输入2,选择OK,对于没有明确定义属性的面,其属性缺省为1,1-21,这些平面要求与物理区和材料联系起来Preprocessor-Attributes-DefinePickedAreas选取线圈平面(在选择对话框里)点取OK材料号窗口输入3,点OK,1-22,加通量平行边界条件Preprocessorloadsapply-magnetic-boundary-flux-parl选OnLines并选取相应的线选OK,“所选取的线”,“所选取的线”,注:
未划分单元前,加上这种边界条件,1-23,生成有限元网格利用智能尺寸选项来控制网格大小Preprocessor-Meshing-SizeCntrls-smartsize-basic,选择OK,1-24,Preproc-Meshing-Mesh-Areas-Free在选取框内选择ALL选择OK打开绘制单元的材料属性UtilityPlotCtrlsNumbering,选择OK,1-25,力边界条件标志需要单元部件,即一组具有“名称”的单元把衔铁定义为一个单元组件选择衔铁平面Utilityselectentities,用此选项在图形窗口中选择平面,再次选择用APPLY,一旦衔铁已选好,选择OK(在选取框内),1-26,选择与已选平面相对应的单元,选择OK图示衔铁单元Utilityplotelements,衔铁单元,用“面”,1-27,使单元与衔铁组件联系起来UtilitySelectComp/AssemblyCreateComponent,选择OK,1-28,加力边界条件标志PreprocessorLoadsApply-Magnetic-FlagCompForce,选择OK,施加两个标志,用两个不同的方法来计算力Maxwells应力张量虚功,即使只有一种选项,也要鼠标选取,1-29,以毫米单位生成的模型,最好把模型尺寸变换为国际单位制(变换系数=.001)使整个模型激活UtilitySelectEverything缩放平面-不用拷贝Preprocoperatescaleareas,选择OK,1-30,给线圈平面施加电流密度选择线圈平面UtilitySelectEntity,选择OK(实体选择框)选择线圈平面选择OK(选取框内),1-31,激励线圈要求电流密度,故要得到线圈截面积.PreprocessorOperateCalcGeometricItemsOfAreas选择OK要用线圈面积来计算电流密度,将线圈面积赋予参数CAREAUtilityParameterGetScalarData,选择OK,1-32,下面窗口输入面积的参数名,用于后面电流密度输入,去掉面号(如果有的话),这相应于几何面积总和,选择OK,1-33,把电流密度加到平面上PreprocessorLoadsApplyExcitationOnAreas(因为只激活了线圈平面,可在选取框内选择PickAll),选择OK,1-34,进行计算Solu-solve-electromagnetOpt&Solve,选择OK,这些适用于用BH数据来进行的分析,本题将忽略,1-35,生成磁力线圈Postprocplotresults2Dfluxlines选择OK,使用缺省设置,选择OK,(在通常情况下,可这样做),单元边缘围绕的一个红色输廓表示该区域为同类材料号,1-36,计算力PostprocElec&MagCalcComp.Force,选择OK,衔铁上力是在总体坐标系下表示的,此力的方向为使气隙缩小,必须用鼠标选取,1-37,显示总磁通密度值(BSUM)PostprocPlotResultsNodalSolution,选择OK,1-38,第二章第2节,二维静磁学,1-39,EMAG模拟的概念,模型边界条件有:
磁通量垂直磁通量平行周期性对称*偶对称奇对称根据单元方程式施加边界条件矢量(2D或3D)标量(3D)基于单元边(3D)*在第2章来讨论,简单励磁的平面模型,A,A,B,B,线圈(象征性的),铁芯,空气,1-40,在2D静磁场、交流和瞬态分析中采用磁矢量势方法(MVP)此公式称为MVP,磁通量密度(B)等于矢量势(A)的旋度B=Curl(A)对于二维情况,A只有Z方向分量,在ANSYS中表示为“AZ”自由度模型有二种边界条件描述-Dirichlet条件(AZ约束):
磁通量平行于模型边界Neumann条件(自然边界条件):
磁通量垂直于模型边界,1-41,沿A-A通量平行边界条件需满足:
模型中A-A的左边和右边是相同的几何形状相同材料属性相同左边和右边励磁相位差180度(即方向相反)对称平面边界条件沿A-A必须加约束,B,B,(1/2)对称模型,PoleFace,A,A,Preproc.loadsapplyboundaryfluxparllines,1-42,半对称模型与全模型比较:
磁通量密度是相同的线圈上Lorentz力是相同的贮能为1/2极面上力为1/2加载电流密度与全模型相同,线圈(象征性的),简单导磁体的半对称模型,1-43,沿B-B磁通量垂直边条件需满足B-B线上下两边如下参数是相同的几何形状材料性质B-B线上下两边励磁相同对称面(B-B)边界条件2D磁矢量势(MVP)方式,无须处理加载电流与全模型相同,Quartersymmetrymodelofthesimplemagnetizer,B,B,1-44,1/4模型与全模型比较磁通密度分布相同贮能为1/4所示线圈上的Lorentz力1/2作用在极面上力为1/2,励磁体1/4对称模型,B,B,1-45,单元plane13andplane53用于模拟2D磁场Plane13:
4节点四边形耦合场自由度:
温度,结构,磁电源为Z方向B为线性变化适用于:
Plane13,变压器汇流排传感器线性或任意永磁系统,螺线管磁体(致动器)直线或旋转电机负载机械机械力矩,1-46,plane53:
8节点,四边形耦合场自由度:
磁与电路单元耦合电流为Z方向B可为二次非线性变化通常情况下的推荐使用单元适用于精度要求较高的分析场量分析大型机械力矩,中节点,1-47,定义Plane13的单元类型和单元选项Preprocelementtypeadd/edit/delete选择ADD选择Plane13,用单元类型号给平面赋属性,选择OK,1-48,一旦定义单元类型,要选择单元选项单元选项控制:
2D直流模拟为AZ自由度2D模拟型式轴对称平面点取单元选项,1-49,几何体型式,用于直流模拟,选择OK,用于定义平面属性的参考号,因为plane13用于耦合场模拟,故该单元可以具有应力/应变结构选项,1-50,平面与轴对称比较端部效应平面:
不包括轴对称:
自动包括正向电流方向相反,线圈两种情况都是施加正向电流,铁板,磁流密度矢量显示,铁环,轴对称:
+Z电流方向进平面,平面:
+Z电流方向出平面,1-51,磁力线描述平面:
AZ等值线轴对称:
rAZ等值线,电枢,线圈,定子,平面或轴对称?
平面或轴对称?
1-52,力、能量、电感的描述平面:
单位长度轴对称:
整个圆周上的值力:
轴对称:
无有效径向力(相互平衡)单位弧度力不为零(曲度线圈),1-53,定义Plane53的单元类型和选项Preproc.elementtypeadd/edit/delete选择Add选择8节点Plane53,增加单元类型号以给平面赋属性,选择OK,1-54,定义单元类型后选择单元选项单元选项控制:
绞线圈电压加载选项*连接电路单元与有限元区*模拟运动体的自由度,*包括交流分析,plane53单元模拟运动导体示图,选择OK,1-55,考虑速度效应时,要求增加输入信息在实常数中定义。
选择单元选项后,定义实常数是很方便的Preprocrealconstants.,选择,1-56,平面属性要求实常数设置号,速度单位:
米/秒角速度单位:
弧度/秒,选择OK,平面属性要求赋予单元类型号数,1-57,第二章第3节,2-D静磁学,1-58,求解模型的单位制:
SI,分析中使用的单位制为国际单位制:
SI,力(牛顿)能量(焦耳)功率(瓦)长度(米)时间(秒)质量(公斤),磁通密度B(特斯拉)磁场强度H(安培/米)电流(安培)电阻率(欧姆-米)电压V(伏)电感L(亨)磁导率r(亨/米)电容(法拉),1-59,基本关系式:
B=H,其中=r0可为单一值(线性)各相同性或正交各向异性Preprocmaterialpropsisotropic,平面属性要求赋予材料质性号,r相对磁导率,1-60,可为非线性,以模拟饱和状态BH曲线数据能从ANSYS55材料库中获得缺省的BH材料库在ansys55目录下的matlib子目录中:
Preproc.materialpropsmateriallibrarylibrarypath,通过指定路径可在其它位置得到材料数据,1-61,BH数据可用如下方式输入Preprocmaterialpropsmateriallibraryimportlibrary,选择材料,选择材料属性,选择OK,1-62,BH数据生成图形和列表显示,表示在列表显示中的数据点号,材料号,1-63,数据也可列成表格.这种表格也能人工制成Utilitylistpropertiesdatatables,选择OK,1-64,BH曲线输入指南,数据点(0,0)不要输入定义曲线弯曲处的数据点要密(见M54的数据点)BH曲线要避免生成S形通常M钢定义BH数据到8,000A/m数据需要外推这些曲线的值通常需要附加大量的数据以使得值由大逐渐变到最终斜率最终斜率为空气值(0),1-65,BH数据输入,应用实例:
400系列不锈钢输入如下数据,H(A/m)B(T)790.0.771575.1.102365.1.307875.1.5015750.1.5631500.1.6347245.1.6678740.1.70,1-66,首先定义数据表,然后把BH数据输入数据表中Preprocmaterialpropsdatatablesdefine/activate,平面属性要求赋予材料号,选择OK,1-67,利用编辑激活表格输入BH数据Preprocmaterialpropsdatatableseditactive,输入数据后,用鼠标点取FileApply/Quit图示:
Preprocmaterialpropsdatatablesgraph列表:
Preprocmaterialpropsdatatableslist,1-68,实际求解需要用到d/dB2为避免粗劣的v=Yu条件曲线,-B2应该是单调的。
UtilityplotdatatablesgraphNUvs.B*2,1-69,把该曲线数据放置在库内,以备将来使用。
Preproc.materialpropsmateriallibraryexportmaterial,选择文件名,选取生成的BH数据的材料属性,1-70,应用实例:
轴对称直流致动器,课题描述轴对称线圈为直流供电衔铁居中但悬空在定子上方。
分析顺序用axis2d宏建模完成建模后,加边条件求解后处理,力磁动势误差范数电流,磁力线路径图示能量电感,“气隙”(mm),“线圈”部件,“衔铁”部件,材料号5(同衔铁),1-71,励磁直流施加到线圈:
3安培,性质衔铁/定子:
上述BH曲线线圈:
300匝,26线径,r=1空气:
r=1,单位:
毫米(mm),1-72,对于大多数应用,通常指定电压,线圈电流是算出来的.26线规直径(Dw)=0.404mm(在20摄氏度下)铜电阻率()=17.14E-9-m(在20摄氏度下)匝数(N)=300线圈中径为8mm(Rmid)均匀填充圆线圈的电阻为:
R=16000NRmid/Dw2R=4.03对于静态分析,12V电压相应的电流为2.98安,本分析采用3安。
1-73,参数化建模需要:
参数GAP必须定义在命令行输入gap=.5并回车,点取OK,选择分网密度Preprocsizecntrlbasic,1-74,axis2s宏生成模型衔铁单元部件ARMATURE线圈面积参数ACOND线圈单元部件COIL,在ANSYS命令窗口输入axis2s并回车,以建立模型,1-75,材料号1为空气完善边界条件通量平行边界条件Preprocloadsapplyboundary-fluxparl-lines选择模型边界上的所有线,1-76,如下方式定义材料号1(自由空间磁导率)Preprocmaterialpropsisotropic,选择OK,选择OK,1-77,给线圈平面加载线圈电流Preprocloadsapplyexcitation-currentdensity-areas选择线圈平面,选择OK,1-78,给衔铁加力边界条件标志PreprocessorLoadsApply-Magnetic-FlagCompForce,选择OK,用不同的方法计算力,故加载两种标志Maxwell应力张量虚功,选择ARMATURE,1-79,选择所有几何和有限元实体进行模拟Solutionelectromagneticopt&solve选择OK(采用缺省设置进行求解),请确认,1-80,磁力线Postprocplotresults2Dfluxlines注意漏磁位置线圈区定子上角定子与衔铁交界位置,1-81,计算力Postprocelec&magcalccomp.force,轴对称模型只产生垂直方向力定义单元表项FVW_Y虚功Y方向力FMX_YMaxwell应力Y方向力环状模型力总和,选择OK,1-82,用与衔铁邻接的空气单元来计算衔铁力,并显示首先选择空气单元,1)首先选择空气单元-材料属性为1,选择Apply,2)用Num/Pick从中选取邻近衔铁面空气单元,用框选取,1-83,虚功方法计算垂直力并用等值图显示Postprocplotresultselemtable,在气隙中选取空气单元,选择OK,1-84,用路径图示选项(PATH)能获得沿衔铁面的力的分布图必须定义路径Postprocpathoperationsdefinepathbynodes,点取节点2,给一个任意的名字,增加沿路径的数据采样点的数量,点取节点1,选择OK,1-85,路径定义信息如下路径内的结果插值是在总体坐标系下(与柱坐标系或其它局部坐标系相比)路径由直线组成,1-86,单元表项FVW_Y中的力必须插值到路径上Postprocpathoperationsmapontopath,任意名,选择ETAB.FVW-Y,选择OK,1-87,将FVW_Y沿路径显示Postprocpathoperations-plotpathitems-ongeometry,路径图示迭加在几何体上,已定义,将路径显示图缩放到一个较好的程度,选择OK,1-88,节点,作用在衔铁上的垂直方向力的路径图示,1-89,离路径节点节点1的距离,路径上的力(F_Y)也能打印输出Postprocpathoperationslistpathitems,选择OK,1-90,线圈Lorentz力(JxB)选择线圈区域并定义为一个部件。
Utilityselectcomp/assemblyselectcomp/assembly选择线圈为Lorentz力定义单元表Postproelementtabledefinetable,选择,1-91,任意名,作用于整个圆环上的X方向的Lorentz力,选择OK,选择Add,1-92,线圈X方向Lorentz力的等值图Postprocelementtableplotelemtable,选择OK,1-93,作用在线圈单元上的总力Postproelementtablesumofeachitem,该操作作用于全部激活单元上相当于360圆周上的受力力单位为牛顿:
N,1-94,根据节点磁场值差异估计误差,且作为单元表数据贮存Postpromag&eleccalcerroreval,B_ERR单位(T)H_ERR单位(Amps/m)BN_ERR和HN_ERR由最大值归一化,1-95,BN_ERR能用磁力线图进行等值显示图示BN_ERR单元表项Postproelementtableplotelemtable,选择OK激活NOERASE选择Utilityplotcntrlseraseoptions,1-96,图示磁力线Postproplotresults2Dfluxlines选择OK,1-97,线性和非线性材料的共能计算Postproelec&magcalccoenergy选择OK,1-98,也能计算贮能Postprocelec&magcalcenergy,注:
铁的共能大约是贮能的8倍,表示铁的饱和效应所致,1-99,铁单元的磁导率能用等值图显示Postproelementtable选择ADD,plane53单元在线帮助,选择OK这是绝对磁导率,1-100,为了获得相对磁导率,单元表应乘以MUZ系数将自由空间磁导率赋予参数:
MUZ=12.57x10-7Postproelementtableadditems,用已有名字,自由空间磁导率参数,不需要第二个单元表项,选择OK,1-101,用等值图显示相对磁导率MURPostproelementtableplotelemtable,注意饱和区,选择tOK,1-102,沿闭合线计算磁动势MMF确保整个模型都被激活必须定义围绕线圈的路径Postproelec&magcalcdefinepathbynodes,选取如图所示的7个节点,可从任一节点开始,路径的最终节点应与起始节点是同一个,跨越空气隙时,气隙两边的铁边界上各选取一个节点,1-103,完成路径定义,由于铁与空气的界面处H值不连续,故应增加采样点的数目,选择OK,1-104,绕闭合回线计算MMFPostproelec&magcalcMMF选择OK,MMF正负号由右手定则决定,路径的反时针方向与线圈电流的方向相反(对于轴对称模型,正电流方向为进行平面方向),1-105,为了确定铁芯饱和程度,沿定子的中间部分定义一个路径并计算MMF,选取节点1,选取节点2,MMF=-384A-t,1-106,输入的总安匝数为900,铁芯的中间部位有384安匝,也就是空气隙中只有519安匝(忽略其余铁芯中的磁动势)如果384安匝中的大部分都在空气隙中,磁力会有多大?
对于本问题,电磁力至少会增加2倍。
可用另一种方法显示这一点:
将铁芯的磁导率设为1000,进行线性求解,评述,1-107,检查边界条件的正确与否非常重要模型边界磁力线的检查通量平行(用磁力线图检查)通量垂直(用磁力线图检查)电流观察选择线圈组件Postproelec&magcalccurrent选择OK,1-108,对于线性系统:
Postproelec&magcalcTerminalpar对于非线性系统一两种理论选项割线定义增量定义简易割线计算利用共能(C),L=2Ci2,仿照电机计算更精确的方法:
LMATRIX宏Solumagneticinductance见帮助文档中的说明和实例,电感计算,1-109,第二章第4节,二维静磁学,1-110,永磁体,线性永磁体感应曲线为线性可模拟大部分稀土磁体计算需要有“感应曲线”要求两种材料性质相对磁导率r各向同性正交各向异性矫顽磁力Hc矢量值利用单元坐标系定义材料性质缺省:
总体直角坐标系,H(Amp/m),B(T),Br,Hc,固有曲线,感应曲线,第二象限曲线图,稀土磁体典型曲线,1-111,r和Hc可以是随温度变化磁化方向平行/垂直于磁体中心线相
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ansys 电磁场 仿真 分析 教程