硅钢片损耗计算.docx

任意频率正弦波条件下铁磁材料 损耗的计算崔杨，胡虔生，黄允凯(东南大学电气工程学院，江苏省南京市四牌楼 2 号 210096)Iron Loss Prediction in Ferromagnetic Materials withSinusoidal SupplyCUI Yang， HU Qian-sheng， HUANG Yun-kai(School of Electrical Engineering, Southeast University, Nanjing 210096, China)摘要： 本文首先介绍了铁耗分立计算模型，随后采用标准规定的用爱泼斯坦方圈测硅钢片损耗的方法对铁 磁材料进行损耗实验，对实验结果数据进行回归分析计算出了铁耗分立模型中的未知参数。并分析了参数 的特性，将其应用于铁耗计算中，所得出的结果非常接近于实际值。在此基础上进一步分析了铁耗各分量 随频率、磁密变化的规律。结论对于铁耗分析有非常重要的参考意义。关键字： 铁耗；铁磁材料；回归分析；爱泼斯坦方圈Abstract: The paper presents loss separation model, then the method of iron loss measurement by means of anEpstein frame prescripted in standard is employed to the loss experiment, parameters in the model are calculated through a method called regression, using the experiment result. Parameters are used in predicting iron loss, there is hardly any discrepancy between the computed and the measured results. In the meantime the relationship bitween the loss contribution and frequency, flux density is discussed based on the computed result. Conclution is very valuable for the loss prediction.Keywords: Iron loss; Ferromagnetic material; Regression; Epstein frame1 引言随着电力电子技术的发展， 各种新型电 机在各行各业得到了广泛的应用， 电机铁耗 的准确计算也成为越来越重要的课题， 引起 不少学者的注意。 目前在国内设计电机中是 假设硅钢片内磁场分布均匀， 利用硅钢片供 应商提供的硅钢片在工频正弦波电源下的 损耗曲线和经验公式来近似计算铁耗。 对于 一般电机， 用此方法进行铁耗计算基本可以 满足要求。 但是在各种特种电机特别是高速 电机中，往往是由高频、非正弦电源供电， 如果电机的铁耗计算仍停留在采用工频时 的方法， 主要频率损耗值通过简单缩放比例形式确定，势必会存在较大的误差。基金项目 ：国家自然科学基金项目(50477021)Project Supported by National Natural Science Foundationof China(50477021)在国外， 已经有不少学者提出了铁耗计 算的两种方向， 一种是采用有限元法来分析 硅钢片内磁场的分布， 进而计算损耗； 另一 种是通过研究铁磁材料的磁特性， 提出铁耗 的模型及计算和测量方法。 第一种方法虽然 准确，但计算工作量巨大，且没有通用性。 另一种方法计算方便，其中以Bertottti铁耗 分立计算模型23应4用最为广泛。用这种 模型计算出来的结果与实测数据相差不大。 它的主要问题在于模型中存在未知参数， 且 难于确定， 参数的大小将直接影响到损耗计 算的结果，要求参数的计算必须非常准确。 而国内目前在这方面没有专门研究。在传统电机设计方法中使用的损耗曲 线在低频条件下是按照国家标准GB/T 36552000用爱泼斯坦方圈测量电工钢片(带)磁性能的方法 5 中所规定的方法测量出来 的， 中频条件下的测试则参照GB 10129 88 电工钢片 (带) 中频磁性能测量方法 6。 (以下统称 标准)但在变频高速电机中， 频率是变化的， 直接通过测试的方法得到其 铁耗值是非常困难的， 因此必须要找到一个 准确的铁耗计算公式。本文介绍了Bertotti 铁耗分立计算模型，再以之为研究对象，通 过回归分析准确计算出了其中的未知参数。2 正弦波供电条件下硅钢片损耗的 计算按照 Bertottti 铁耗分立计算模型，不考 虑集肤效应时，铁芯损耗可以用下式表达：P = P + P + Pir h ec e= afB x + bf 2 B 2 + ef 1.5 B 1.5 (1)m m m式中 Pir 为铁耗W/kg； Ph 为磁滞损 耗W/kg； Pec 为经典涡流损耗 W/kg； Pe 为附加损耗W/kg； f 为频率Hz； Bm 为磁 密幅值T； a、 b、 x、 e 为未知系数。a、 b、 x、 e 显然是用式(1)进行铁耗 计算的关键参数。 对于参数的计算，有几种 方法，各有优劣，本文采用直接拟合计算方 法计算未知参数。2.1 损耗的测量测量参照标准中的规定，使用标准 25cm 爱泼斯坦方圈为硅钢片损耗的主要测 试仪器。 测量电路原理图如图 1 所示。 在工 频条件下，电源为电网接调压器； 在非工频正弦波条件下， 电源为数字信号发生器加线形功率发生器为了测量和计算方便， 使用数 字功率分析仪测量各个量。 为了考察不同频 率条件对硅钢片损耗的影响， 选取了牌号为DW31035 的试样，分别在不同频率条件 下测损耗值，所得的数据如图 2 所示。Epstein方圈功率分析仪图 1 硅钢片损耗测试原理图Fig.1 Schematic of measurement of the electrical steel sheet loss2.2 回归分析利用已知的一系列不同频率下的损耗 值，对式(1)作回归分析，得出其中未知 参数，约束条件为使式(2)达到最小：(2)c = (|(P - P*sisPsi)|2si式中 c 为偏差， P 为实测损耗值，Ps *为预测值。采用硅钢片 DW31035 为试样，其在 正弦波供电下损耗曲线如图 2。为了能够比 较数据选取对参数计算结果的影响， 把数据 按频率分成几组， 每一组分别进行回归分析 计算出参数，具体计算结果见表 1。表 1 直接拟合回归计算的参数结果Tab.1 Direct fitting computed result of parametersb5.867e-0054.952e-0055.030e-005注*：表中“所有频率”是指包括 50、 60、 100、 200、 400、 1000、 2000、 5000Hz 在内的所有频率频率(Hz)50、 60、 100、 200、 400400、 1000、 2000、 5000所有频率*e2.735e-0062.625e-0062.945e-006a0.02380.02670.0247x1.6801.7151.682分组2551W/kg50Hz60Hz100Hz200Hz400Hz1000Hz2000Hz5000HzT0.1 1 1.5图 2 各种频率的正弦波供电条件下 DW310 35 的损耗曲线Fig.2 Loss curve of DW310从表 1 中数据可以看出， 各种情况下计 算出来的参数 a、 b、 x、 e 基本上没有差别， 由此可以认为对于某一型号的硅钢片来说 这些参数是与频率无关的常数。 在回归计算 求解参数时， 为了得到较为准确的结果，应 尽量选取较多的频率数和点数。 第组结果 可以被认为是所有计算结果中最准确的， 可 作为 DW310 35 的损耗计算参数参考值 (表 2)。 将参数值代入式 (1)， 利用公式计 算出来的铁耗值与实际值作比较， 误差见表 3 ，可见两者已基本趋于一致。表 2 DW310 35 的 a、 b、 x、 e 参考值Tab.2 Reference parameters of DW310 35参数 a b xe参考值0.02475.030e-0051.6822.945e-006 35 with sinusoiad supply表 3 损耗计算值误差Tab.3 Errors of loss calculationf50Hz60Hz100Hz200Hz400Hz1000Hz2000Hz5000Hz0.008130.003950.002190.002100.002650.005860.002810.009552.3 扩展讨论将上述计算方法用于不同牌号的硅钢 片，可以计算出对应于每个硅钢片的参数a、 b、 x、 e 的参考值，结果如表 4 所示。 表中 的各参数值各不相同，说明这些参数只与磁 性材料本身的特性和厚度有关 。 因此对应于 每一种铁磁材料， 都有一组参数值，只要计 算出了这些参数值，利用式(1)，就能计算任意频率条件下的损耗值了。表 4 不同牌号的硅钢片的 a、 b、 x、 e 参考值Tab.4 Reference parameters of different silicon steel sheet牌号*abxeDW315500.01122.840e-41.6624.420e-7DW465500.02811.919e-41.7417.470e-7DW540500.03191.494e-41.8981.435e-83 计算结果的分析使用式(1)和表 2 的参数值可以对硅 钢片 DW31035 的铁耗进行分析。铁耗由 磁滞、经典涡流和附加损耗三部分组成，一 般附加损耗所占的比率非常小，约为 0.1 0.2%，因此分析时可忽略不计。综合图 3、 4、 5、 6 可得出以下结论：随着频率、磁密的增加， 总损耗也在大； 频率增加时涡流损耗增加的速度比磁 滞损耗快， 且涡流损耗在总损耗中所占 的比率也在增大；磁密增大时涡流损耗在总损耗中所占 的比率也是在增大。30W/kg磁滞损耗涡流损耗2010T 00.5 1 1.5 2 图 3 损耗分量随磁密变化规律(f=400Hz)Fig.3 Loss contributions change vs. frequency(f=400Hz)W/kg10磁滞损耗涡流损耗741Hz200 400 600 800 1000图 4 损耗分量随频率变化规律(Bm=0.5T)Fig.4 Loss contributions change vs. frequency(Bm=0.5T)%800.5T1.0T1.5T6040Hz200 400 600 800 1000图 5 磁滞损耗所占比例随频率变化示意图Fig.5 Percentage of the hysteresis loss vs. frequency70%500.5T301.0T1.5THz 10200 400 600 800 1000图 6 涡流损耗所占比例随频率变化示意图Fig.6 Percentage of the eddy current loss vs. frequency以上结论说明在高频、 高磁密条件下涡 流耗是损耗的主要分量， 损耗主要由铁此材 料中涡流引起。因此当设计的电机将在高 频、 高磁密条件下运行时，应考虑选用涡流 效应较小的铁磁材料， 如选用厚度较薄的硅 钢片等。4 结论通过利用 Bertotti 铁耗分立计算公式对 多种频率下损耗数据的回归分析提出了有 效的预测损耗的方法， 提出了选取数据应该 遵守的原则。 计算出来的参数直接用于正弦 波条件下任意型号硅钢片的损耗预测计算 时结果误差很小， 其对于非正弦条件下损耗 计算也有一定的参考价值。 同时给出了部分 铁磁材料的参数参考值。 对计算结果的分析 得出了铁耗各分量随频率、磁密变化的规 律。 这些结论可以直接应用于对电机的铁耗 分析及高速电机设计中。参考文献1 陈世坤 . 电机设计 M. 第二版 , 北京 , 机械工业出版社 , 2000, 6.2 G. Bertotti. General properties of power losses in softferromagnetic materialsJ. IEEE Trans. Magn., 1988,24(1): 621 630.3 Hyuk Nam, Kyung-Ho Ha, Jeong-Jong Lee, , et al. AStudy on Iron Loss Analysis Method Considering the Harmonics of the Flux Density Waveform Using Iron Loss Curves Tested on Epstein SamplesJ. IEEE Transactions on Magnetics, 2003, 39(3): 1472- 1475.4 R. Kaczmarek, M. Amar, A general formula forprediction of iron losses under nonsinusoidal supply voltage wavefonnJ, IEEE Trans. On Magnetics, 1995, 31(5): 2505-2509.5 国家质量技术监督局 . GB/T 3655 2000. 用爱泼斯坦方圈测量电工钢片(带)磁性能的方法 S, 北京 : 中国标准出版社 , 2000, 106 国 家 质 量 技 术 监 督 局 . GB 10129 88. 电 工 钢 片(带)中频磁性能测量方法 S, 北京 : 中国标准出版社 , 1988, 127 Torres, A.G.; Cardoso Filho, B.J. A generalized Epsteintest method for the computation of core losses ininduction motorsC. IECON 02, 2002, 2: 1150 1155.8 Yicheng Chen, P. Pillay. An improved formula forlamination core loss calculations in machinesoperating with high frequency and high flux densityexcitationC. Industry Applications Conference, 2002.37th IAS Annual Meeting. Conference Record, 2:759-766.9 Aldo Boglietti, Andrea Cavagnino, Mario Lazzari, et al.Predicting Iron Losses in Soft Magnetic MaterialsWith Arbitrary Voltage Supply: An EngineeringApproachJ. IEEE Transactions on Magnetics, 2003,39( 2): 981-989.10 刘国强 , 赵凌志 , 蒋继娅编著 . Ansoft 工程电磁场有限元分析 M. 北京 , 电子工业出版社 , 2005, 8.11 袁海林， 何松波主编， 微特电机设计手册 M， 上海，上海科学技术出版社， 1998