一种考虑动态磁滞效应的高效稳定时域有限元计算方法

doi:10.19595/j.cnki.1000-6753.tces.221180

摘要
图/表
参考文献
相关文章 (12)

全文: PDF (4384 KB) HTML
输出: BibTeX | EndNote (RIS)

摘要在时域有限元分析中,考虑铁心磁滞特性与损耗等效磁场对磁场分布的反馈作用,对精细化模拟电工装备的暂态电磁特性具有重要意义。针对耦合动态磁滞特性有限元分析中,存在计算时间长、收敛不稳定等问题,该文提出一种基于等效磁阻率的改进固定点迭代策略,并应用于耦合Preisach磁滞模型的时域有限元分析中。首先,鉴于逆Preisach模型易于数值实现且具有较高的计算精度,铁心材料的磁滞行为由逆Preisach模型进行描述,并基于解析一阶回转曲线对模型进行了参数辨识;其次,通过分析动态磁场分量对收敛性的影响,提出基于等效磁阻率的固定点迭代策略,并在有限元迭代过程中引入松弛因子以加强算法的收敛速度与稳定性;最后,以爱泼斯坦方圈为例,计算铁心中的磁场分布特性与瞬时损耗分布特性并进行实验验证。结果表明,该文所提动态磁滞有限元计算方法在具有较高精度的同时具有较好的收敛速度及稳定性。

	服务

	把本文推荐给朋友
	加入我的书架
	加入引用管理器
	E-mail Alert
	RSS
	作者相关文章
	魏鹏
	陈龙
	贲彤
	井立兵
	张献

关键词 ：时域有限元, 固定点技术, 磁滞建模, 铁心损耗

Abstract：In the time-domain finite element analysis, considering the hysteresis properties of the iron core and the feedback effect of the loss equivalent magnetic field on the magnetic field distribution is of great significance for finely simulating the transient electromagnetic characteristics of electrical equipment. Aiming at the problems of long calculation time and divergence in the finite element analysis coupled with dynamic hysteresis characteristics, this paper proposes an improved fixed point iteration strategy based on equivalent reluctivity, and applies it to time domain finite element analysis implemented with Preisach model. With the modified reluctivity and the introduced relaxation factor, the finite element algorithm proposed in this paper has good convergence speed and stability with high accuracy.
Firstly, since numerical application of the inverse Preisach model is easy and has high accuracy, the hysteresis behavior of the core material is described by the inverse Preisach model. The Preisach model is established based on the analytical first order reversal curves, of which the parameters are identified by several symmetric hysteresis loops. Secondly, by analyzing the influence of the dynamic magnetic field component on the convergence,a fixed-point iterative strategy based on equivalent reluctivity is proposed, and a relaxation factor is introduced in finite element iterations to enhance the convergence speed and stability of the algorithm. Finally, taking the Epstein frame as an example, experiments and finite element calculations are carried out on the magnetic field distribution and instantaneous loss characteristics of the iron core. There are 9 analytical polynomial parameters of B30P105 electrical steel sheet, identified by Genetic Algorithm. The static hysteresis loops simulated by the Preisach model is of high accuracy. In the 2-D dimensional computation model, the triangle located at the limb-yoke region is selected as an observation spot. The calculated dynamic hysteresis loops, exciting current and instantaneous loss are compared with measured ones, under voltage excitation amplitude of 9 V and 13 V, respectively. The maximum error of the instantaneous loss is less than 6%. After one period computation, the total iron loss is calculated by integrate the area of dynamic hysteresis loop in each element. The calculated iron losses under different flux density are in substantial agreement. In order to investigate the effect of relaxation factor on the convergence property of the program, various values of relaxation factor are chosen and the results show that the average number of iterations and computation time are inversely proportional to the magnitude of the relaxation factor. However, the iteration process doesn't converge with the relaxation factor set as 0.8.
The following conclusions can be drawn from the simulation analysis: (1) The Preisach hysteresis model based on analytical first order reversal curves can accurately simulate the static hysteresis characteristics of electrical steel sheets. (2) By incorporating the dynamic hysteresis effect into the improved fixed-point iteration process, the proposed finite element algorithm can realize precise calculation of low-frequency dynamic hysteresis magnetic field, with good computation speed and stability. (3) The error of the calculated instantaneous loss is small, which meets engineering application requirements.

Key words： Time-domain finite element method fixed-point technique hysteresis modeling iron losses

收稿日期: 2022-06-21

PACS:

TM153

基金资助:基金项目国家自然科学基金（52007102, 52207012, 51977147）和省部共建电工装备可靠性与智能化国家重点实验室（河北工业大学）开放课题基金（EERIKF2021015）资助项目

通讯作者: 陈龙男,1989年生,博士,讲师,研究方向为磁性材料磁特性模拟、全局优化设计。E-mail：chenlong@ctgu.edu.cn

作者简介: 魏鹏男,1996年生,硕士研究生,研究方向为磁性材料磁特性、电磁场数值分析。E-mail：m18805167298@163.com

引用本文:

魏鹏, 陈龙, 贲彤, 井立兵, 张献. 一种考虑动态磁滞效应的高效稳定时域有限元计算方法[J]. 电工技术学报, 2023, 38(21): 5661-5672. Wei Peng, Chen Long, Ben Tong, Jing Libing, Zhang Xian. An Efficient and Stable Time Domain Finite Element Method Considering Dynamic Hysteresis Effect. Transactions of China Electrotechnical Society, 2023, 38(21): 5661-5672.

链接本文:

https://dgjsxb.ces-transaction.com/CN/10.19595/j.cnki.1000-6753.tces.221180 https://dgjsxb.ces-transaction.com/CN/Y2023/V38/I21/5661

[1] 刘刚, 荣世昌, 武卫革, 等. 基于混合有限元法和降阶技术的油浸式变压器绕组2维瞬态流-热耦合场分析[J]. 高电压技术, 2022, 48(5): 1695-1705.
Liu Gang, Rong Shichang, Wu Weige, et al.Two-dimensional transient flow-thermal coupling field analysis of oil-immersed transformer windings based on hybrid finite element method and reduced-order technology[J]. High Voltage Engineering, 2022, 48(5): 1695-1705.
[2] Cui Han, Ngo K D T. Transient core-loss simulation for ferrites with nonuniform field in SPICE[J]. IEEE Transactions on Power Electronics, 2019, 34(1): 659-667.
[3] 李冰, 王泽忠, 刘恪, 等. 特高压变压器直流偏磁对绕组电流的影响[J]. 电工技术学报, 2020, 35(7): 1422-1431.
Li Bing, Wang Zezhong, Liu Ke, et al.Research on winding current of UHV transformer under DC-bias[J]. Transactions of China Electrotechnical Society, 2020, 35(7): 1422-1431.
[4] 陈龙龙, 魏晓光, 焦重庆, 等. 混合式高压直流断路器分断过程电磁瞬态建模和测试[J]. 电工技术学报, 2021, 36(24): 5261-5271.
Chen Longlong, Wei Xiaoguang, Jiao Chongqing, et al.Electromagnetic transient modeling and test of hybrid DC circuit breaker[J]. Transactions of China Electrotechnical Society, 2021, 36(24): 5261-5271.
[5] 李永建, 闫鑫笑, 张长庚, 等. 基于磁-热-流耦合模型的变压器损耗计算和热点预测[J]. 电工技术学报, 2020, 35(21): 4483-4491.
Li Yongjian, Yan Xinxiao, Zhang Changgeng, et al.Numerical prediction of losses and local overheating in transformer windings based on magnetic-thermal-fluid model[J]. Transactions of China Electrotechnical Society, 2020, 35(21): 4483-4491.
[6] Yamazaki K, Fukushima N.Torque and loss calculation of rotating machines considering laminated cores using post 1-D analysis[J]. IEEE Transactions on Magnetics, 2011, 47(5): 994-997.
[7] 迟青光, 张艳丽, 陈吉超, 等. 非晶合金铁心损耗与磁致伸缩特性测量与模拟[J]. 电工技术学报, 2021, 36(18): 3876-3883.
Chi Qingguang, Zhang Yanli, Chen Jichao, et al.Measurement and modeling of lossand magnetostrictive properties for the amorphous alloy core[J]. Transactions of China Electrotechnical Society, 2021, 36(18): 3876-3883.
[8] Lin D, Zhou P, Fu W N, et al.A dynamic core loss model for soft ferromagnetic and power ferrite materials in transient finite element analysis[J]. IEEE Transactions on Magnetics, 2004, 40(2): 1318-1321.
[9] Zhou Ping, Lin Dingsheng, Lu Chuan, et al.A new algorithm to consider the effects of core losses on 3-D transient magnetic fields[J]. IEEE Transactions on Magnetics, 2014, 50(2): 365-368.
[10] Lin D, Zhou P, Chen Q M, et al.The effects of steel lamination core losses on 3D transient magnetic fields[J]. IEEE Transactions on Magnetics, 2010, 46(8): 3539-3542.
[11] 赵志刚, 马习纹, 姬俊安. 考虑频率效应的正弦及直流偏磁条件电工钢磁滞特性模拟及验证[J]. 中国电机工程学报, 2021, 41(23): 8178-8186.
Zhao Zhigang, Ma Xiwen, Ji Junan.Simulation and verification on hysteresis characteristics of electrical steel under sinusoidal and DC bias conditions considering frequency effects[J]. Proceedings of the CSEE, 2021, 41(23): 8178-8186.
[12] Xiao Meng, Zuo Yuxiang, Li Yang, et al.Core loss calculation of anode saturable reactor in damping oscillation state based on J-A theory[J]. IEEE Transactions on Applied Superconductivity, 2021, 31(8): 1-4.
[13] 刘任, 李琳, 王亚琦, 等. 基于随机性与确定性混合优化算法的Jiles-Atherton磁滞模型参数提取[J]. 电工技术学报, 2019, 34(11): 2260-2268.
Liu Ren, Li Lin, Wang Yaqi, et al.Parameter extraction for Jiles-Atherton hysteresis model based on the hybrid technique of stochastic and deterministic optimization algorithm[J]. Transactions of China Electrotechnical Society, 2019, 34(11): 2260-2268.
[14] 陈龙, 易琼洋, 贲彤, 等. 全局优化算法在Preisach磁滞模型参数辨识问题中的应用与性能对比[J]. 电工技术学报, 2021, 36(12): 2585-2593, 2606.
Chen Long, Yi Qiongyang, Ben Tong, et al.Application and performance comparison of global optimization algorithms in the parameter identification problems of the preisach hysteresis model[J]. Transactions of China Electrotechnical Society, 2021, 36(12): 2585-2593, 2606.
[15] 赵小军, 刘小娜, 肖帆, 等. 基于Preisach模型的取向硅钢片直流偏磁磁滞及损耗特性模拟[J]. 电工技术学报, 2020, 35(9): 1849-1857.
Zhao Xiaojun, Liu Xiaona, Xiao Fan, et al.Hysteretic and loss modeling of silicon steel sheet under the DC biased magnetization based on the preisach model[J]. Transactions of China Electrotechnical Society, 2020, 35(9): 1849-1857.
[16] Hussain S, Lowther D A.An efficient implementation of the classical preisach model[J]. IEEE Transactions on Magnetics, 2018, 54(3): 1-4.
[17] 段娜娜, 徐伟杰, 李永建, 等. 基于极限磁滞回线法的软磁复合材料磁特性模拟[J]. 电工技术学报, 2018, 33(20): 4739-4745.
Duan Nana, Xu Weijie, Li Yongjian, et al.Electromagnetic property modeling of the soft magnetic composite material based on the limiting loop method[J]. Transactions of China Electrotechnical Society, 2018, 33(20): 4739-4745.
[18] 张长庚, 李永建, 杨庆新. 一种考虑磁滞和微观涡流效应的电磁模拟方法[J]. 中国电机工程学报, 2016, 36(21): 5966-5974, 6041.
Zhang Changgeng, Li Yongjian, Yang Qingxin.An electromagnetic simulation method considering hysteresis and micro-eddy current effect[J]. Proceedings of the CSEE, 2016, 36(21): 5966-5974, 6041.
[19] Yue Shuaichao, Anderson P I, Li Yongjian, et al.A modified inverse vector hysteresis model for nonoriented electrical steels considering anisotropy for FEA[J]. IEEE Transactions on Energy Conversion, 2021, 36(4): 3251-3260.
[20] Dlala E, Belahcen A, Arkkio A.Efficient magnetodynamic lamination model for two-dimensional field simulation of rotating electrical machines[J]. Journal of Magnetism and Magnetic Materials, 2008, 320(20): 1006-1010.
[21] Dlala E.Efficient algorithms for the inclusion of the preisach hysteresis model in nonlinear finite-element methods[J]. IEEE Transactions on Magnetics, 2011, 47(2): 395-408.
[22] Sadowski N, Batistela N J, Bastos J P A, et al. An inverse Jiles-Atherton model to take into account hysteresis in time-stepping finite-element calculations[J]. IEEE Transactions on Magnetics, 2002, 38(2): 797-800.
[23] Fallah E, Moghani J S.A new approach for finite-element modeling of hysteresis and dynamic effects[J]. IEEE Transactions on Magnetics, 2006, 42(11): 3674-3681.
[24] Dlala E, Belahcen A, Arkkio A.A fast fixed-point method for solving magnetic field problems in media of hysteresis[J]. IEEE Transactions on Magnetics, 2008, 44(6): 1214-1217.
[25] Dlala E, Belahcen A, Arkkio A.Locally convergent fixed-point method for solving time-stepping nonlinear field problems[J]. IEEE Transactions on Magnetics, 2007, 43(11): 3969-3975.
[26] Li Wei, Koh C S.Investigation of the vector Jiles-Atherton model and the fixed point method combined technique for time-periodic magnetic problems[J]. IEEE Transactions on Magnetics, 2015, 51(4): 1-6.
[27] Dlala E, Arkkio A.Analysis of the convergence of the fixed-point method used for solving nonlinear rotational magnetic field problems[J]. IEEE Transactions on Magnetics, 2008, 44(4): 473-478.