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

11111Ansoft 瞬态磁场计算（未考虑外电路）11111总结整理： 2008-6-4于德国Kassel大学目录1、说明2、电机额定运行时的瞬态场分析与计算2.1 Setup Boundaries/source（重点考虑Source Setup）2.2 Setup Solution选择Setup Solution|Options2.3 Setup Solution选择Setup Solution|Motion Setup2.4 Solve | Nominal Problem2.5 Post ProcessPost Press/Transistant date：2.5.1 计算平均输出功率：Average Output Power2.5.2 计算相电流有效值：2.5.3 计算输入电功率，由此可以计算效率2.6 Post ProcessPost 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|Options3.3 Setup Solution选择Setup Solution|Motion Setup3.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，“=”的右边输入（显然，B相应为：，C相应为）(7) 指定Done回到2D Boundary/Source Manager(8) 在Value框内输入U_Phase_A；(9) 选择Strand；(10) 选择Winding，出现Winding Setup指定PhA为Positive，PhReA为Negative；在Resistance框内输入相绕组电阻15.42W，在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 选择中，选DirectFor problem where all of the boundaries are well defined, the direct solver is the best choice（3） Transient analysisi、Solution: Start from time zero在开始时，还没有任何解，因此只能选择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 Step的确定如下：1个齿距对应的时间为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整个电机是求解区域的倍数。缩小求解区域可以降低求解时间。2.3 Setup Solution选择Setup Solution|Motion Setup(1) 从Object列表中选择Band，选择Set Band(2) 选择Mechanical Setupi) 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，只看部分时间段曲线形状。下面的两个图形时2D和RmxPrt得到的相反电动势波形。可以看出，两者差别不大。注意在2D计算中，没有考虑电机的斜槽。另外，2D计算得到的电势是时间的函数，可以在2D的后处理中将横坐标由时间变为位置。2.5 Post ProcessPost Press/Transistant date：在其中可以得到Solve 中的所有曲线结果，但在其中可以对这些结果进行分析和计算，其中包括前面的将横坐标变由时间便为位置。2.5.1 计算平均输出功率：Average Output Power在进行2D计算时，给定输出负载转矩和转速，因此电机的输出功率很容易计算，但是也可通过转矩曲线和转速曲线进行计算。在EMpulse中，电机的功率Pout满足下式：其中，FW，表示机械损耗，Pair_gap为气隙功率，由平均转矩（单位Nm）和转速（单位rad/second）相乘而得。转速为375rpm=39.27rad/s。以下给出平均转矩的计算方法。1） Choose Plot/Open, Select torque.dat, choose OK;2) Choose Tools/Calulator to access the Signal Calculator;3) Select Torque.dat, and choose Copy to copy the torque plot into top of stack of calculator.4) Choose Sample, and define the following 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) Enter 39.27 in the Name/Constant field(给定速度)7) Choose * to multiply by the speed in radians per second;8) Choose the integrate button（计算在一段时间内转矩之和）9) Enter 0.1 in the Name/Constant field(给计算平均值的时间段：Stop time-Start time=0.1s)10) Choose “/” to calculate the average;11) Choose Preview. The last number in this plot is the average value(曲线的最后一点就是所求的平均转矩)；12) Choose max,(给出所求的平均转矩)：816.141Nm因此输出功率为：Pout=816.141-70=746.141；由于给定转矩为19Nm，而不是要求的19.1，因此输出功率不是750W。2.5.2 计算相电流有效值：电流有效值的数值计算公式：1） Choose Plot/Open, Select current.dat, choose OK;2) Choose Tools/Calulator to access the Signal Calculator;3) Select current.dat:A_phase, and choose Copy to copy it into top of stack of calculator.4) Choose Sample, and define the following 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（计算在一段时间内i2之和）;9） Enter 0.1 in the Name/Constant field(给计算平均值的时间段：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）Choose max (得所要求的电流有效值)：1.568相电流的EMpulse计算值为1.568A，RmxPrt的计算结果为1.539A，两者差别很小同样可以计算外加相电压有效值，外加相电压是标准的正弦波，因此该值就是V。2.5.3 计算输入电功率，由此可以计算效率 一相输入平均功率的数值计算公式为：1） Choose Plot/Open, Select current.dat、v_src.dat (外加三相项电压) choose OK;2) Choose Tools/Calulator to access the Signal Calculator;3) Select current.dat:A_phase, and choose Copy to copy it into top of stack of calculator.4) Choose Sample, and define the following parameters（将A相相电流离散化）: 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.(此时Calculator Stack中的A相电流表示为：Sample(Current.dat_phase)6） 将Calculator Stack 中的Sample(Current.dat_phase)Load到Loaded signals 中；7） 将Calculator Stack 中的Sample(Current.dat_phase)Clrs(删除)；8） Select v_src.dat:A_phase, and choose Copy to copy it into top of stack of calculator.9) Choose Sample, and define the following parameters（将A相相电压离散化）: Sample: Time Specify by:Size Start:0.2 Stop:0.3 Size:100010） Choose OK to accept the values and return to the signal calculator.11）选中Loaded signals中的Sample(Current.dat_phase)，并将其Copy到Calculator Stack12）Choose “*”，将上述电压和电流相乘；13）Choose （积分），将各成绩量相加；（Choose the integrate button）14）Enter 0.1 in the Name/Constant field(给计算平均值的时间段：Stop time-Start time=0.1s)15）Choose “/” to calculate the average;16）Choose （进行开方运算）（这一步可以在最后一步进行，计算结果差别不大）；17）Choose Preview. The last number in this plot is the average value(曲线的最后一点就是所求的A相输入功率)；18）Choose max (得所要求的A相输入功率)：306.223W同样可得B相、C相输入功率，因此三相输入功率为：918.669W，RmxPrt中计算得到的输入功率为：942.967W，相差不多。显然场的计算值较小，这是由于在场的计算中没有考虑铁耗，如果考虑铁耗（RmxPrt中的铁耗为13.36W），则两者差别更小。有上的结算结果可以得到EMpulse 中计算得到的运行效率为746.141/918.669*100%=81.2%，RmxPrt计算的效率为79.54，相差不大。2.6 Post ProcessPost Process/field：在Post Process Saved Fields 中选择任一pjt，选择Post Process，进入2D Post Processor。2.6.1 齿部磁密分布和磁密平均值计算2.6.1.1 齿部磁密分布 定义齿部曲线（定义一条包括所有齿的圆弧，以得到电机的最高齿磁密，这一磁密定义为齿磁密） Geometry/Creat/Arc 输入弧的原点坐标(x，y) =（0，0），输入弧半径Rad=60，输入弧的张角Ang=90,选择Enter 在Radius中给出弧的半径60，则自动给出Start Point 为（60，0），可以调整改起始点的坐标，选择Enter。 在Num Points中给出圆弧点的个数2000；不选Clockwise，表示圆弧从前面的起点开始不是顺时针（是逆时针）；在Angle中给出圆弧角度90；Name 改为Tooth；选择Enter。 注意：该路径是一圆弧，即得到的齿部磁密分布是沿圆弧的磁密分布，与实际上的齿部磁密分布并不相同，实际的齿部磁密分布应该是沿弦的分布。 绘制tooth中的磁密分布 Data/Calculator，进入Field Calculator：out0窗口 QTY/B：定义关于磁密B的后续操作 Geom/Line：选择tooth，沿线定义B Unit Vec/2D-Normal：定义发向磁密 Dot Geom/Line 2D Plot（横坐标为机械角度、纵坐标为T）可以看出，齿磁密最大值为1.46T，而RmxPrt计算得到的空载时的齿磁密为1.635T，可以用EMpulse 计算电机空载时的磁密。对齿磁密来说，只计算径向磁密即可，切向磁密很低。2.6.1.2 齿部磁密平均值计算应首先定义积分路径，然后根据下式计算这一路径处的平均磁密，计算磁密幅值沿路径的积分，再除以路径的长度，就是平均磁密。由于积分变量的单位同分母的单位一致，因此在此没有单位的变化问题，同时可以只计算切向磁密或径向磁密的平均值（见后面的例子）。i）定义齿 定义齿部路径ab，该路径名字为Line1，如下图所示ii）计算该齿部的平均磁密 Data/Calculator 到field calculato：out0窗口 Qty/B，选择B矢量，Register 中结果如下图 Mag，求B矢量的大小，Register 中结果如下图 Geom/Line，选择B的积分线，积分线为Line1。Register 中结果如下图 （Integral），对B幅值沿Line1进行积分，Register 中结果如下图（以上计算出上面计算式中的分子，下面计算分母） Num/Scalar，给出积分函数1.0，Register 中结果如下图 Geom/Line，选择积分函数1.0的积分路径，积分线为Line1。Register 中结果如下图 （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.743T2.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空载电势与前面额定运行时计算的空载电势完全相同（Back Voltage），因为转速相同；3.4.2齿部磁密计算可以看出空载和负载时齿部磁密的最大值相差不大，但显然空载和负载各齿中磁密的分布是不相同的。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 FerritePower ferrite core loss is based on the following equation（单位重量的铁耗）Where:Cm is constant value determined by experiment.f is the frequency.（Hz）Bmax is the maximum amplitude of the flux density(T)电机设计中给出的铁耗计算公式为：其中，GFe钢的重量；ka经验系数，把由于钢片加工、磁通密度分布的不均匀、磁通密度随时间不按正弦规律变化以及旋转磁化与交变磁化之间的