Ansoft瞬态场计算

11111

Ansoft 瞬态磁场计算（未考虑外电路）

11111 总结整理：

2008-6-4于德国Kassel大学

目录

1、说明

2、电机额定运行时的瞬态场分析与计算

2.1 Setup Boundaries/source（重点考虑Source Setup） 2.2 Setup Solution

选择[Setup Solution]|[Options] 2.3 Setup Solution

选择[Setup Solution]|[Motion Setup] 2.4 Solve | Nominal Problem 2.5 Post Process

Post Press/Transistant date：

2.5.1 2.5.2 2.5.3

计算平均输出功率：Average Output Power 计算相电流有效值：

计算输入电功率，由此可以计算效率

2.6 Post Process

Post Process/field：

2.6.1齿部磁密分布和磁密平均值计算 2.6.1.1 齿部磁密分布 2.6.1.2 齿部磁密平均值计算

2.6.1.3 沿齿弧磁密分布和沿齿磁密分布的差异 2.6.2 定子轭部磁密分布和磁密平均值 2.6.3 定子轭部磁密分布和磁密平均值 2.6.4 气隙磁密分布和磁密最大值

3、电机空载额定转速运行时的瞬态场分析与计算

3.1 Setup Boundaries/source（重点考虑Source Setup） 3.2 Setup Solution

选择[Setup Solution]|[Options]

3.3 Setup Solution

选择[Setup Solution]|[Motion Setup]

3.4 计算结果

4、考虑铁耗的计算结果

4.1 铁耗计算设置

4.2 额定负载时考虑与不考虑铁耗时的比较 4.3 负载很小时的比较 4.3.1 考虑铁耗时 4.3.2 不考虑铁耗时

5、Maxwell 与RmxPrt计算结果比较

5.1 磁密及额定值比较 5.2 额定转速时永磁相电势比较

5.3 气隙磁密分布

6、其他心得

1、说明

以16极36槽调速永磁同步电动机为例进行分析电机的瞬态场计算，电机由RmxPrt 开始，并将该模型加到MAXWELL 11中。RmxPrt 的项目为ad_pmsm(Maxwell File 22KB)，Maxwell 中的项目为ad_pmsm_fem.pjt。

2、电机额定运行时的瞬态场分析与计算

2.1 Setup Boundaries/source（重点考虑Source Setup）

A_phase 绕组源的设置： (1) 选择A相绕组；

(2) 使用[Assign]|[Source]|[Solid]；

(3) 指定[Solid]为Voltage，[Name]改为A_Phase； (4) 选择[Options]，将源的形式[Constant]改为[Function]； (5) 选择[Function]；

(6) 选择[Add]，在函数输入框“=”的左边输入U_Phase_A，“=”的右边输入

350*sqrt(2/3)*sin(360*50*T)

（显然，B

s*qrt(2/3)*sin?(T36，0C相应为：350相应为

350s*qrt(2/3)*sin?(T36）0

(7) 指定[Done]回到2D Boundary/Source Manager (8) 在[Value]框内输入U_Phase_A； (9) 选择[Strand]；

(10) 选择[Winding]，出现[Winding Setup]

指定[PhA]为[Positive]，[PhReA]为[Negative]；

在[Resistance]框内输入相绕组电阻15.42?，在[Inductance]框内输入绕组端部漏感0.00112H（在RMxPrt中有该值的输出结果）；

在[Total turns as seen from terminal]框内输入每相串联匝数684； 在[Number of Parallel Branches]输入并联支路数1； 选择[OK]退回2D Boundary/Source Manager

(11) 指定[Assign]------同样给出A_phase、B_phase绕组源的设置 (12) 选择[File]|[Save]与[File]|[Exit]，保存并退出边界条件编辑器

2.2 Setup Solution

选择[Setup Solution]|[Options]

（1） 选择[Manual Mesh…](进行自定义剖分，略)；

在进行完Manual Mesh 后必需的一步是Mesh/Line Match，选择主、从边界的边，to ensure that the meshing points will match at their matching boundaries. If they don not, you will receive an error message about a missing transcript file during the normal solution. （2） 在[Solber Choice] 选择中，选Direct

For problem where all of the boundaries are well defined, the direct solver is the best choice

（3） Transient analysis i、Solution: Start from time zero

在开始时，还没有任何解，因此只能选择Start from time zero

如果对这个问题已有解，可以选择Continue Previous Solution。问题设置可以以任何方式改变（除了结构变化），求解从前解结果开始进行。比如初解的结果终止计算时间为0.2s，则在调整问题设置时终止计算时间变为0.4s，并且选择Start from time zero，则计算从前面的0.2s计算结果开始。 ii、Time Step

Time Step 的大小可以根据一个电机齿距范围内求解点数来确定。假定电机的转速为375rpm，电机槽数为36，电机一个定子齿距的求解点为10个，则Time Step的确定如下：

375rpm?375rev/60s?1rev/0.16s

1个齿距对应的时间为0.16/36=0.004444s； Time Step=0.004444s/10=0.000444s。

最后确定Time Step 为0.0004s iii、Stop Time

该值关系不大，可以从小值开始，如果电机还没有达到稳定，则可以增大Stop Time，电机的求解可以从前面解的结果开始进行。 iv、Model Depth=105mm

电机的轴向铁心长度。注意用2D瞬态场求解时，没有也无法考虑电机的斜槽。 v、Symmetry Multiplier=4

整个电机是求解区域的倍数。缩小求解区域可以降低求解时间。

2.3 Setup Solution

选择[Setup Solution]|[Motion Setup]

(1) 从[Object]列表中选择Band，选择[Set Band] (2) 选择[Mechanical Setup]

i) Initial Angular Velocity: 375

给电机的实际运行转速。如果计算额定转速时的状态，则给额定转速值

ii) Moment of Inertia: 0.0012

给电机的转动惯量。这一惯量值在RMxPrt中已计算出；实际上电机的稳态性能与

转动惯量没有关系，但在2D瞬态场计算时，如果这一惯量太大，可能导致计算结果发

散。因此该值一般比RmxPrt的计算值要小(实际计算值为0.009262)。 iii) Damping:0.0454

由于风阻和其他机械损耗所导致的阻尼，很明显是一个经验值，为了比较路的设计计算结果，该值应该与RmxPrt中的给定值相同。

iv) Load Torque: -19

所要仿真的电机负载转矩，当然可以是额定输出转矩。这一转矩值以负值形式给定。

2.4 Solve | Nominal Problem

进行求解。进行求解时，可以随时通过refresh 观察求解结果，主要观察求解是否收敛。求解结果见本窗口中[Solutions]中的Transient Date，其中有各种曲线。可以调整曲线下方的[Settings]，只看部分时间段曲线形状。

? （Integral），对1.0沿Line1进行积分。Register 中结果如下图

?

（divide），分子与分母相除。Register 中结果如下图

? Eval ，显示计算结果。Register 中结果如下图

计算得到平均磁密为1.3795T。

计算Line1中法向磁密平均值 ? Qty/B，选择B矢量；

? Geom/Line。选择积分线Line1，目的求单位法向矢量； ? Unit Vec/2D_normal，求Line1的单位法向矢量；

? Dot。矢量B与单位法向矢量点乘（标量积），得到B的法向分量，该值已经是标量，

可以再计算其绝对值。 ? Geom/Line。选择积分线Line1 ?

（Integral），对B的法向分量沿Line1进行积分

? Num/Scalar，给出积分函数1.0；

? Geom/Line，选择积分函数1.0的积分路径，积分线为Line1。 ? ?

（Integral），对1.0沿Line1进行积分。 （divide），分子与分母相除。

? Eval ，显示计算结果；

法向磁密平均值为：1.378T。同前面计算的总磁密平均值1.3795T向差无几，说明齿中磁密主要是法向磁密。实际上计算得到的切向磁密只有0.02T。 2.6.1.3 沿齿弧磁密分布和沿齿磁密分布的差异

下图是分析研齿弦磁密分布和沿各齿磁密分布的示意图。

沿弧磁密分布的示意图如下图，图中给出的是径向磁密和切向磁密，可以看出切向磁密可以忽略不计，每齿中磁密分布是均匀的。各齿中磁密分布：B1=-1.17T；B2=-1.1T，B3=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。与沿弧计算的磁密几乎相同。因此为了计算方面，可以直接定义圆弧，画出沿圆弧的法向磁密分布，找出最高磁密齿，该齿磁密的平均值就是用瞬态场计算得到的齿部磁密。

值得注意的是，用瞬态场计算得到不同转子位置时的定子齿部磁密，其最大值都是相同的。

给定点的坐标，即可以用（X，Y），又可以用（Rad，Ang）。 2.6.2

定子轭部磁密分布和磁密平均值

由于不知道那个轭部磁密最大，因此需要计算多个轭部磁密，如下图所示。各轭部段的外点直径为155，内点直径为31.92，各段的角度分别为1（85?）、2（75?）、3（65?）、4（55?）、5（45?）。

I）1段轭部磁密分布及平均磁密（合成磁密的平均值为：0.062T）

II）2段轭部磁密分布及平均磁密（合成磁密的平均值为：0.686T）

III）3段轭部磁密分布及平均磁密（合成磁密的平均值为：0.232 T）

IV）4段轭部磁密分布及平均磁密（合成磁密的平均值为：0.645 T）

v）5段轭部磁密分布及平均磁密（合成磁密的平均值为：0.407T）

结论：轭部最高平均磁密为0.686T，RmxPrt计算得到的额部空载磁密为0.743T

2.6.3 定子轭部磁密分布和磁密平均值

比较下图所示1、2、3、4四个转子轭部磁密的分布和平均值。可以先找出轭部1所对应的两点坐标，2、3、4只要将坐标分别转动22.5?、45?、67.5?即可。

I）1段转子轭部磁密分布及平均磁密（合成磁密的平均值为：0.726 T）

II）2段转子轭部磁密分布及平均磁密（合成磁密的平均值为：0.726 T）

III）3段转子轭部磁密分布及平均磁密（合成磁密的平均值为：0.689 T）

IV）4段转子轭部磁密分布及平均磁密（合成磁密的平均值为：0.709T）

结论：转子轭部磁密主要是法向磁密，与位置关系不大，计算得到的最高转子轭部平均磁密为0.726T，RmxPrt计算得到的空载转子轭部磁密为0.7561T。 2.6.4

气隙磁密分布和磁密最大值

计算平均气隙处的磁密分布和磁密最大值。 磁密最大值为0.9T，RmxPrt得到的磁密最大值为0.837T。

3、电机空载额定转速运行时的瞬态场分析与计算

3.1 Setup Boundaries/source（重点考虑Source Setup）

将前述负载运行时的电压源都赋为0V，各相绕组电阻都变为无穷大即可，其它不变。 3.2 Setup Solution

选择[Setup Solution]|[Options] 同前面的额定负载运行设置 3.3 Setup Solution

选择[Setup Solution]|[Motion Setup]

? 从[Object]列表中选择Band，选择[Set Band] ? 选择[Mechanical Setup] ? 不选Consider Mechnical Setup

? Constant Angular Velocity=375(给定运行转速)

3.4 计算结果

3.4.1 同； 3.4.2

齿部磁密计算

空载电势与前面额定运行时计算的空载电势完全相同（Back Voltage），因为转速相

可以看出空载和负载时齿部磁密的最大值相差不大，但显然空载和负载各齿中磁密的分布是不相同的。 3.4.3

气隙磁密

可以看出空载负载气隙磁密分布有一定的重合性。用空载磁密代替负载磁密问题不大。

4、考虑铁耗的计算结果

4.1 铁耗计算设置

在前面的额定运行和空载运行分析时，没有考虑电机的铁耗，下面考虑电机的铁耗分析电机的额定运行，在输出转矩和转速不变时，比较考虑和不考虑铁耗时电机的输入功率和效率。

在涡流场和瞬态场中，可以计算铁耗，在FEM的铁耗根据磁路的方法进行，并用公式计算，可以根据已知的铁耗系数不同，分成两种方法，即已知磁滞损耗系数和涡流系数时用Electrical steel，在已知总的铁耗系数时采用Power ferrite。通常给出的是铁耗系数，因此一般采用Power ferrite。

在其他同额定负载运行设置完全相同的情况下，设置Setup Executive Parameters/Core loss

? 选择Object 中的定子

? 选择Copute Core Loss on Object/Power Ferrite

Power ferrite core loss is based on the following equation（单位重量的铁耗）

y p?CmfxBmaxWhere:

Cm is constant value determined by experiment. f is the frequency.（Hz）

Bmax is the maximum amplitude of the flux density(T) 《电机设计》中给出的铁耗计算公式为：

pFe?kapheGFe

其中，

GFe——钢的重量；

ka——经验系数，把由于钢片加工、磁通密度分布的不均匀、磁通密度随时间不按正弦规律变化以及旋转磁化与交变磁化之间的损耗差异等而引起的损耗增加都估计在内。 pFe——相当于前面FEM中的p，但表达式为：

phe?p10/50B2(f1.31)?p10/50()1.3f1.3Bmax2 5050比较Fem 中单位重量铁耗公式与上式，可以看出

Cm?p10/50(11.3),x?1.3,y?2 50其中，p10/50——当B=1T、f=50Hz时，钢单位重量内的铁耗，其值可按硅钢片型号，从下表查取。 钢的种类 低含硅量硅钢片 中含硅量硅钢片 P10/50 2.80 2.2 相当于国产钢号 D12 D22

高含硅量硅钢片 经计算，Cm=0.0136 ? 赋值如上图所示 ? Exit/Save

2.0 D23 4.2 额定负载时考虑与不考虑铁耗时的比较

图 不计铁耗时相电流

图 计铁耗时相电流

不计铁耗时电流大小为1.568A，计铁耗时相电流为：1.568A，电流不变 三相功率：不计铁耗时为：918.669W，计铁耗时为918.294W，两者相差不大。

4.3 负载很小时的比较

在负载转矩只有1Nm时，其他同额定运行相同，比较考虑和不考虑铁耗时相电流大小和输入功率大小，比较可以看出，考虑铁耗和不考虑铁耗，对电机的计算没有影响。 4.3.1

考虑铁耗时

图 气隙转矩

图 相电流（有效值为：1.067A，输入功率66*3W）

图 铁耗

4.3.2

不考虑铁耗时

图 气隙转矩

图 相电流（有效值为：1.0686A，输入功率66*3W）

5、Maxwell 与RmxPrt计算结果比较

5.1 磁密及额定值比较

气隙磁密最大值（T） 齿部磁密（T） 定子额部磁密（T） 转子轭部磁密（T） 定子电流（A） 额定运行效率（%） Maxwell11（空载负载差不多，以负载为主） RmxPrt（空载） 0.9 1.46 0.686 0.726 1.568 81.2 0.837 1.635 0.743 0.756 1.539 79.54 5.2 额定转速时永磁相电势比较

图 Maxwell计算结果

图 RmxPrt计算结果

两者计算结果极为相近，不过在RmxPrt计算中，考虑到了电机的斜槽，而Maxwell中无法考虑斜槽。 5.3 气隙磁密分布

本文来源：https://www.bwwdw.com/article/9yz6.html