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

An Efficient and Stable Time Domain Finite Element Method Considering Dynamic Hysteresis Effect

作　　者：魏鹏 陈龙 贲彤 井立兵[2] 张献[3] Wei Peng;Chen Long;Ben Tong;Jing Libing;Zhang Xian(College of Electrical Engineering and New Energy China Three Gorges University,Yichang 443002 China;Hubei Provincial Engineering Technology Research Center for Microgrid China Three Gorges University,Yichang 443002 China;State Key Laboratory of Reliability and Intelligence of Electrical Equipment Hebei University of Technology,Tianjin 300130 China)

机构地区：[1]三峡大学电气与新能源学院,宜昌443002 [2]湖北省微电网工程技术研究中心(三峡大学),宜昌443002 [3]省部共建电工装备可靠性与智能化国家重点实验室(河北工业大学),天津300130

出　　处：《电工技术学报》2023年第21期5661-5672,共12页Transactions of China Electrotechnical Society

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

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

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

分 类 号：TM153[电气工程—电工理论与新技术]