Ansoft瞬态场计算步骤及与RMprt的结果对比

1、Ansoft 瞬态磁场计算（未考虑外电路）总结整理：2008-6-4于德国Kassel大学目录1、说明电机额定运行时的瞬态场分析与计算电机额定运行时的瞬态场分析与计算Setup Boundaries/source （重点考虑 Source Setup）Setup Solution选择Setup SolutionlOptionsSetup Solution选择Setup SolutionlMotion SetupSolve | Nominal ProblemPost ProcessPost Press/Transistant date：计算平均输出功率： Average Output Power

2、计算相电流有效值：计算输入电功率，由此可以计算效率Post ProcessPost Process/field：齿部磁密分布和磁密平均值计算 齿部磁密分布 齿部磁密平均值计算沿齿弧磁密分布和沿齿磁密分布的差异定子轭部磁密分布和磁密平均值定子轭部磁密分布和磁密平均值2.6.4 气隙磁密分布和磁密最大值3、3、电机空载额定转速运行时的瞬态场分析与计算Setup Boundaries/source （重点考虑 Source Setup）Setup Solution选择Setup SolutionlOptionsSetup Solution选择Setup Sol ution

3、|Motion Setup3.4 计算结果4、考虑铁耗的计算结果4.1 铁耗计算设置额定负载时考虑与不考虑铁耗时的比较负载很小时的比较4.3.1 考虑铁耗时4.3.2 不考虑铁耗时5、Maxwell与RmxPrt计算结果比较磁密及额定值比较额定转速时永磁相电势比较气隙磁密分布6、其他心得1、说明以 16 极 36 槽调速永磁同步电动机为例进行分析电机的瞬态场计算，电机由 RmxPrt 开始，并将该模型加到 MAXWELL 11 中。RmxPrt 的项目为 ad_pmsm(Maxwell File 22KB), Maxwell 中的项目为 ad_pmsm_fem.pjt。2、电机额定运行时的瞬态

4、场分析与计算Setup Boundaries/source (重点考虑 Source Setup)A_phase 绕组源的设置：选择A相绕组；使用AssignlSourcelSolid；指定Solid为 VOltage， Name改为 A_Phase；选择Options，将源的形式Constant改为F unction；选择Function； 选择Add，在函数输入框“=”的左边输入U_Phase_A，“ = ”的右边输入350*sqrt(2 / 3)*sin(360T*5(显 然，B 相 应为：350 * sqrt (2 / 3) * sin(360 *50 * T 120)， C 相应为3

5、 5 0s * r t (2/3)* sin- (T 3 6)0指定Done回到 2D Boundary/Source Manager在Flue框内输入 U_Phase_A；选择Strand；（10）选择Winding，出现Winding Setup指定PhA为Positive, PhReA为Negative；在Resistance框内输入相绕组电阻15.42，在Inductance框内输入绕组端部漏感0.00112H （在RMxPrt中有该值的输出结果）；在Total turns as seen from terminal框内输入每相串联匝数 684；在N umber of Parallel

6、 Branches输入并联支路数 1 ；选择OK退回 2D Boundary/Source Manager（11）指定Assign同样给出A_phase、B_phase绕组源的设置（12）选择FilelSave与FilelExit，保存并退出边界条件编辑器Setup Solution选择Setup Solution|Options（1）选择Manual Mesh（进行自定义剖分，略）；在进行完Manual Mesh后必需的一步是Mesh/Line Mofch,选择主、从边界的边，toensure that the meshing points will match at their match

7、ing boundaries. If they don not, youwill receive an error message about a missing transcript file during the. normal solution（2）在Solber Choice选择中，选 DirectFor problem where all of the boundaries are well defined, the direct solver is the best choice3) Transient analysisi、Solution: Start from time zer

8、o在开始时，还没有任何解，因此只能选择 Start from time zero 如果对这个问题已有解，可以选择 Continue Previous Solution 。问题设置可以以任何方 式改变（除了结构变化），求解从前解结果开始进行。比如初解的结果终止计算时间为 0.2s， 则在调整问题设置时终止计算时间变为0.4s，并且选择Start from time zero，则计算从前面 的 0.2s 计算结果开始。ii、Time StepTime Step 的大小可以根据一个电机齿距范围内求解点数来确定。假定电机的转速为375rpm，电机槽数为36,电机一个定子齿距的求解点为10个，则Time

9、 Step的确定如下:375 rpm = 375 rev / 60 s = 1rev / 0.16 s1 个齿距对应的时间为 0.16/36=0.004444s； Time Step=0.004444s/10=0.000444s。最后确定 Time Step 为 0.0004siii、Stop Time该值关系不大，可以从小值开始，如果电机还没有达到稳定，则可以增大 Stop Time 电机的求解可以从前面解的结果开始进行。iv、Model Depth=105mm电机的轴向铁心长度。注意用 2D 瞬态场求解时，没有也无法考虑电机的斜槽。v、Symmetry Multiplier=4整个电机是求

10、解区域的倍数。缩小求解区域可以降低求解时间。Setup Solution选择Setup SolutionlMotion Setup从Object列表中选择 Band,选择Set Band选择Mechanical SetupInitial Angular Velocity: 375 给电机的实际运行转速。如果计算额定转速时的状态,则给额定转速值Moment of Inertia: 0.0012给电机的转动惯量。这一惯量值在RMxPrt中已计算出；实际上电机的稳态性能与转动惯量没有关系,但在 2D 瞬态场计算时,如果这一惯量太大,可能导致计算结果发散。因此该值一般比 RmxPrt 的计算值要小(实

11、际计算值为 0.009262)。Damping:0.0454由于风阻和其他机械损耗所导致的阻尼，很明显是一个经验值，为了比较路的设计计算结果，该值应该与 RmxPrt 中的给定值相同。Load Torque: -19所要仿真的电机负载转矩，当然可以是额定输出转矩。这一转矩值以负值形式给定Solve | Nominal Problem进行求解。进行求解时，可以随时通过 refresh 观察求解结果，主要观察求解是否收敛。求解结果见本窗口中Solutions中的Transient Date，其中有各种曲线。可以调整曲线下方的Settings，只看部分时间段曲线形状。Torque vs Time5a

12、-ICOUtlSm19Nm0 2-11 生0 iNaeEsalPosition vs Time在给宦时间内电机转了 3Position vs Time在给宦时间内电机转了 3转多TorquePower LossTerminal Voltage Back EMFFlux Linkage IMinding Curt已ntPositionSpeedDamping TorqueCor已 Lossns mU下面的两个图形时2D和RmxPrt得到的相反电动势波形。可以看出，两者 差别不大。注意在2D计算中，没有考虑电机的斜槽。另外，2D计算得到的电 势是时间的函数，可以在2D的后处理中将横坐标由时间变为位

13、置。29 May 2008Ansoft Corporation17:05:17Induced Winding Voltages at Rated SpeedRMxprtDesignlPost ProcessPost Press/Transistant date：在其中可以得到 Solve 中的所有曲线结果，但在其中可以对这些结果进行分析和计算， 其中包括前面的将横坐标变由时间便为位置。计算平均输出功率： Average Output Power在进行 2D 计算时，给定输出负载转矩和转速，因此电机的输出功率很容易计算，但是 也可通过转矩曲线和转速曲线进行计算。在 EMpulse 中，电机的功率

14、 Pout 满足下式：P = P -FWout air _ gap其中，FW,表示机械损耗，Pair_gap为气隙功率，由平均转矩(单位Nm)和转速(单位 rad/second)相乘而得。转速为375rpm=39.27rad/s。以下给出平均转矩的计算方法。Choose Plot/Open, Select torque.dat, choose OK;ChooseTools/Calulator to access the Signal Calculator;Select Torque.dat, and chooseCopy to copy the torque plot into top of

15、stack of calculator.Choose Sample, and define the following parameters:丁 Sample: Time丁Specify b尹Size丁 Start:0.2丁 Stop:03/Size:1000Choose OK to accept the values and return to the signal calculator.Enter 39.27 in the Name/Constant field(给定速度)Choose * to multiply by the speed in radians per second;Cho

16、ose the integrate button (计算在一段时间内转矩之和)Enter 0.1 in the Name/Constant/eM(给计算平均值的时间段：Stop time-Start time=0.1s)Choose “/” to calculate the average;Choose Preview. The last number in this plot is the average value(曲线的最后一点就是所求的平均转矩)；12) Choose max,(给出所求的平均转矩)：816.141Nm因此输出功率为：Pout=816141-70=746141；由于给定

17、转矩为19Nm，而不 是要求的 19.1，因此输出功率不是 750W。252 计算相电流有效值：电流有效值的数值计算公式：I = - fi 2( n)N iChoosePlot/Open, Select current.dat, choose OK;ChooseTools/Calulator to access the Signal Calculator;Selectcurrent.dat:A_phase, and chooseCopy to copy it into top of stack of calculator.ChooseSample, and define the followi

18、ng parameters:丁 Sample: Time丁 Specify by:Size丁 Start:0.2丁 Stop:0.3/Size:10005）Choose OK to accept the values and return to the signal calculator.6）Choose Push to duplicate the entry;7）Choose * to multiply the value in top of the stack by itself（计算 i2）;8）Choose the integrate button （计算在一段时间内 i 之和）；9）

19、Enter 0.1 in the Name/Constantfield（给计算平均值的时间段：Stop time-Start time=0.1s）10）Choose “/” to calculate the average;11）Choose 厂（进行开方运算）（这一步可以在最后一步进行，计算结果差别不大）；12）Choose Preview. The last number in this plot is the average value（曲线的最后一点就是所求 的电流有效值） ；13）Choosemax （得所要求的电流有效值）： 1.568相电流的EMpulse计算值为1.568A,

20、RmxPrt的计算结果为1.539A,两者 差别很小同样可以计算外加相电压有效值， 外加相电压是标准的正弦波，因此该值就是350 定义齿部曲线（定义一条包括所有齿的圆弧，以得到电机的最高齿磁密，这一磁密定义为齿磁密）令 Geometry/Creat/Arc 输入弧的原点坐标fr，y） = （0，0），输入弧半径Rad=60，输入弧的张角Ang=90 , 选择 Enter（X； y） . （Rad, Ang） 作为坐标点，给出一 个，另一个自动给出 对应值Snap To: 两 Vertex 雨 Grid I- Other .Arc Axis:窗ZCenter给出Start Point为（60，0

21、），可以调整改起 给出Start Point为（60，0），可以调整改起始点的坐标，选择Enter。Cancel 在Num Points中给出圆弧点的个数2000；不选Clockwise，表示圆弧从前面的起点开始不是顺时针（是逆时针）；在Angle中给出圆弧角度90 ； Name改为Tooth；选择 Enter。齿部磁密分布是沿圆弧的磁密分布，与实际上的齿部磁选择 Enter。齿部磁密分布是沿圆弧的磁密分布，与实际上的齿部磁密分布并不相同，实际的齿部磁密分布应该是沿弦的分布。 绘制tooth中的磁密分布 Data/Calculator，进入 Field Calculator： out0 窗口Q

22、TY/B:定义关于磁密B的后续操作Geom/Line :选择tooth，沿线定义BUnit Vec/2D-Normal：定义发向磁密DotGeom/Line2D Plot可以看出，齿磁密最大值为1.46T，而RmxPrt计算得到的空载时的齿磁密 为1.635T,可以用EMpulse计算电机空载时的磁密。对齿磁密来说，只计算径向磁密即可，切向磁密很低。齿部磁密平均值计算 应首先定义积分路径，然后根据下式计算这一路径处的平均磁密，计算磁密幅值沿路径的积分，再除以路径的长度，就是平均磁密。J B dlavg J (1.0) dl由于积分变量的单位同分母的单位一致，因此在此没有单位的变化问题，同时可以

23、只计算切向磁密或径向磁密的平均值（见后面的例子）。i）定义齿定义齿部路径ab,该路径名字为Linel,如下图所示ii）计算该齿部的平均磁密Data/Calculator 到 field calculato：out0 窗口选择B矢量，Register中结果如下图 Mag,求B矢量的大小，Register中结果如下图 Geo/ZLise,选择B的积分线，积分线为Linel。Register中结果如下图Integral),对B幅值沿Linel进行积分，Register中结果如下图以上计算出上面计算式中的分子，下面计算分母) NwM/Scalar,给出积分函数1.0，Register中结果如下图 Ge

24、o/ZLine,选择积分函数1.0的积分路径，积分线为Linel。Register中结果如下图(.Integral),对1.0沿Linel进行积分。Register中结果如下图* (divide),分子与分母相除。Register中结果如下图 Eval，显示计算结果。Register中结果如下图计算得到平均磁密为 1.3795T。计算 Line1 中法向磁密平均值QM，选择B矢量；Geom/Lineo选择积分线Linel，目的求单位法向矢量;Unit Vec/2D_normal,求 Linel 的单位法向矢量；Doto矢量B与单位法向矢量点乘(标量积)，得到B的法向分量，该值已经是标量,可以再

25、计算其绝对值。 Geom/Line o 选择积分线 Linel(Integral),对B的法向分量沿Linel进行积分 Nm/SCalar,给出积分函数l.0； Geom/Line，选择积分函数l.0的积分路径，积分线为Linel。Integral)，对 l.0 沿 Linel 进行积分。(divide),分子与分母相除。 Eval，显示计算结果;法向磁密平均值为：1.378T。同前面计算的总磁密平均值1.3795T向差无几，说明齿中磁密 主要是法向磁密。实际上计算得到的切向磁密只有 0.02T。沿齿弧磁密分布和沿齿磁密分布的差异图是分析研齿弦磁密分布和沿各齿磁密分布的示意图。沿弧磁密分布的示

26、意图如下图，图中给出的是径向磁密和切向磁密，可以看出切向磁密可以忽略不计，每齿中磁密分布是均匀的。各齿中磁密分布：B1=-1.17T；B2=-1.1TB3=0.775T， B4=1.38T， B5=-0.375T， B6=1.46T。图 沿齿弧的磁密分布（只给出法向磁密合切向磁密）通过计算各齿沿弦的法向平均磁密分别为： B1=-1.17T；B2=-1.1075T，B3=0.7723T，B4=1.379T, B5=-0.372T, B6=1.46T。与沿弧计算的磁密几乎相同。因此为了计算方面，可 以直接定义圆弧，画出沿圆弧的法向磁密分布，找出最高磁密齿，该齿磁密的平均值就是用 瞬态场计算得到的齿

27、部磁密。值得注意的是,用瞬态场计算得到不同转子位置时的定子齿部磁密,其最大值都是相 同的。给定点的坐标，即可以用（X, Y）,又可以用（Rad, Ang）。定子轭部磁密分布和磁密平均值 由于不知道那个轭部磁密最大,因此需要计算多个轭部磁密,如下图所示。各轭部段 的外点直径为155，内点直径为31.92，各段的角度分别为1 （85。）、2 （75。）、3 （65。）、4（55。）、 5 （45。）。II） 2段轭部磁密分布及平均磁密（合成磁密的平均值为：0.686T）IV）4 段轭部磁密分布及平均磁密（合成磁密的平均值为：0.645 T）结论：轭部最高平均磁密为0.686T, RmxPrt计算得

28、到的额部空载磁密为0.743T定子轭部磁密分布和磁密平均值比较下图所示 1、2、3、4 四个转子轭部磁密的分布和平均值。可以先找出轭部1 所对应的两点坐标，2、3、4只要将坐标分别转动22.5。、45。、67.5。即可。I） 1段转子轭部磁密分布及平均磁密（合成磁密的平均值为：0.726 T）ResulloM Normal 7 ongent0.5-05-I0123567II） 2 段转子轭部磁密分布及平均磁密（合成磁密的平均值为： 0.726 T）IV） 4段转子轭部磁密分布及平均磁密（合成磁密的平均值为：0.709T）结论：转子轭部磁密主要是法向磁密，与位置关系不大，计算得到的最高转子轭部平

29、均磁密为 0.726T，RmxPrt 计算得到的空载转子轭部磁密为 0.7561T。气隙磁密分布和磁密最大值计算平均气隙处的磁密分布和磁密最大值。 磁密最大值为 0.9T，RmxPrt 得到的磁密最 大值为 0.837T。Setup Boundaries/source （重点考虑 Source Setup）将前述负载运行时的电压源都赋为0V,各相绕组电阻都变为无穷大即可, 其它不变。Setup Solution选择Setup SolutionlOptions同前面的额定负载运行设置Setup Solution选择Setup SolutionlMotion Setup从Object列表中选择Ba

30、nd,选择Set Band选择Mechanical Setup不选 Consider Mechnical SetupConstant Angular ieioci（y=375（给定运行转速）3.4 计算结果空载电势与前面额定运行时计算的空载电势完全相同(Back VOltage),因为转速相同；齿部磁密计算可以看出空载和负载时齿部磁密最大值相分布是不相同的。3.4.3气隙磁密可以看出空载负载气隙磁密分布有一定的重合性。用空载磁密代替负载磁密问题不大。4、考虑铁耗的计算结果4.1 铁耗计算设置 在前面的额定运行和空载运行分析时，没有考虑电机的铁耗，下面考虑电机的铁耗分析 电机的额定运行，在输出转

31、矩和转速不变时，比较考虑和不考虑铁耗时电机的输入功率和效 率。在涡流场和瞬态场中，可以计算铁耗，在FEM的铁耗根据磁路的方法进行，并用公式 计算，可以根据已知的铁耗系数不同，分成两种方法，即已知磁滞损耗系数和涡流系数时用 Electrical steel ，在已知总的铁耗系数时采用 Power ferrite 。通常给出的是铁耗系数，因此一 般采用 Power ferrite。在其他同额定负载运行设置完全相同的情况下，设置 Setup Executive Parameters/Core loss选择Object中的定子选择 Copute Core Loss on Object/Power Fe

32、rritePower ferrite core loss is based on the following equation (单位重量的铁耗)p = C f xB ym maxWhere:Cm is constant value determined by experiment.f is the frequency. ( Hz)Bmax is the maximum amplitude of the flux density(T)电机设计中给出的铁耗计算公式为：p = k p GFe a he Fe其中，Gf钢的重量；Fek经验系数，把由于钢片加工、磁通密度分布的不均匀、磁通密度随时间不按

33、正弦规律a变化以及旋转磁化与交变磁化之间的损耗差异等而引起的损耗增加都估计在内。PFe相当于前面FEM中的p,但表达式为： TOC o 1-5 h z f1p = p B 2()1-3 = p ()1-3 f1-3 B2he 10 /505010 /50 50max比较 Fem 中单位重量铁耗公式与上式，可以看出1C = p ()1-3, x = 1.3, y = 2m 10 /50 50其中，p10/5o当B=lT、f=50Hz时，钢单位重量内的铁耗，其值可按硅钢片型号，从下表查取。钢的种类P10/50相当于国产钢号低含硅量硅钢片2.80D12中含硅量硅钢片2.2D22高含硅量硅钢片2.0D23经计算，C眉0.0136 赋值如上图所示 Exit/Save4.2 额定负载时考虑与不考虑铁耗时的比较图 不计铁耗时相电流图 计铁耗时相电流001Torque Time001Torque Time不计铁耗时电流大小为1.568A,计铁耗时相电流为：1.568A,电流不变三相功率：不计铁耗时为：918.669W，计铁耗时为918.294W，两者相差不大。负载很小时的比较在负载转矩只有 1Nm 时，其他同额定运行相同，比较考虑和不考虑铁耗时相电流大小和输入功率大小，比较可以看出，考虑铁耗和不考虑铁耗，对电