引用本文: | 汪芳宗,朱萌瑶,鄢皓文,等.基于L-稳定的显式积分法的电磁暂态仿真计算方法[J].电力系统保护与控制,2021,49(17):10-17.[点击复制] |
WANG Fangzong,ZHU Mengyao,YAN Haowen,et al.Explicit integration method based on L-stable for electromagnetic transient simulation[J].Power System Protection and Control,2021,49(17):10-17[点击复制] |
|
摘要: |
对电力电子化电力系统的电磁暂态仿真计算,除经常面临状态变量的突变现象以外,还存在间断现象或间断问题。临界阻尼调整方法(Critical Damping Adjustment, CDA)可避免状态变量的突变所导致的数值振荡问题,但CDA方法不太适合处理间断问题。事实上,经典的CDA方法通常是“忽略”间断问题。为进一步优化电磁暂态仿真计算,采用泰勒级数的Padé逼近推导出一类显式且L-稳定的数值方法,既可彻底避免数值振荡问题,也能有效处理间断问题。因此,可以用此类显式方法来替换CDA方法中的隐式欧拉法,其中L-稳定的显式方法专门用于处理突变现象以及间断问题,而隐式梯形积分法主要用于处理其他的状态变量。理论分析以及初步的算例结果表明,相对于经典的CDA方法,这种新的组合方法在面临频繁的间断问题时具有更高的计算效率。 |
关键词: 电磁暂态 数值振荡 临界阻尼调整 间断问题 显式积分法 Padé逼近 |
DOI:DOI: 10.19783/j.cnki.pspc.201406 |
投稿时间:2020-11-13修订日期:2020-12-28 |
基金项目:国家电网公司总部科技项目资助(XTB51201802842) |
|
Explicit integration method based on L-stable for electromagnetic transient simulation |
WANG Fangzong,ZHU Mengyao,YAN Haowen,YE Xiaohui |
(1. College of Electrical Engineering & New Energy, China Three Gorges University, Yichang 443002, China;
2. China Electric Power Research Institute, Beijing 100192, China) |
Abstract: |
The electromagnetic transient simulation of power electronics-dominated power systems often face intermittent phenomena or problems in addition to the abrupt changes of state variables. The Critical Damping Adjustment (CDA) method can avoid the numerical oscillations caused by the abrupt changes of state variables, but the CDA method is found to be inadequate for dealing with intermittent problems. In fact, the classical CDA method usually "ignores" such problems. In order to optimize the electromagnetic transient simulation method, the Padé approximation of Taylor series is used to derive a kind of explicit and L-stable numerical method, which can completely avoid the numerical oscillation problems and also effectively deal with the intermittent problems. Therefore, the implicit Euler method of the CDA method can be replaced by such explicit methods, where the L-stable explicit method is specifically used for abrupt changes or intermittent problems, while the implicit trapezoidal integral rule is mainly used for all other state variables. Theoretical analysis and numerical examples preliminarily verify the effectiveness of this new combined method compared to the CDA approach when facing frequent intermittent problems.
This work is supported by the Science and Technology Project of the Headquarter of State Grid Corporation of China (No. XTB51201802842). |
Key words: electromagnetic transient numerical oscillation critical damping adjustment intermittent problem explicit integration methods Padé approximation |