CN104809322A - 一种建立变压器铁心频变涡流模型的新方法 - Google Patents
一种建立变压器铁心频变涡流模型的新方法 Download PDFInfo
- Publication number
- CN104809322A CN104809322A CN201410031708.3A CN201410031708A CN104809322A CN 104809322 A CN104809322 A CN 104809322A CN 201410031708 A CN201410031708 A CN 201410031708A CN 104809322 A CN104809322 A CN 104809322A
- Authority
- CN
- China
- Prior art keywords
- frequency
- eddy current
- pole
- varying
- core
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000001419 dependent effect Effects 0.000 title claims description 10
- 238000012892 rational function Methods 0.000 claims abstract description 12
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 11
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 7
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 7
- 230000005672 electromagnetic field Effects 0.000 claims abstract description 5
- 238000003475 lamination Methods 0.000 claims description 17
- 230000005284 excitation Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000005290 field theory Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 230000035699 permeability Effects 0.000 claims description 3
- 238000004804 winding Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 4
- 230000001052 transient effect Effects 0.000 abstract description 4
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 239000000306 component Substances 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种建立变压器铁心频变涡流模型的新方法,属于电力系统元件仿真建模领域。首先从电磁场推得的变压器铁心频变涡流的阻抗公式出发,用改进的矢量匹配法对铁心的频域涡流阻抗表达式进行有理函数拟合。然后在此基础上用电路综合的方法实现涡流的频变电路模型。本发明的优势在于,可以在一个很宽的频段内(10Hz-10MHz),只需要6阶电路就可以对涡流的频变阻抗公式进行十分精确的逼近,逼近最大误差不超过0.2%。而采用原有的建立方法逼近误差最高达到10%以上。在有理函数逼近过程中迭代收敛速度快,对拟合初值的要求不高,数值稳定。本发明方法可以提高电力变压器铁心建模的精确性,在电力变压器暂态宽频仿真中有着重要意义。
Description
技术领域
本发明涉及一种建立变压器铁心频变涡流模型的新方法。具体地说是一种用改进的矢量匹配法对涡流的频变阻抗表达式进行有理函数拟合,并在此基础上实现涡流的电路模型,属于电力系统元件仿真建模领域。
背景技术
在电力系统暂态分析中,变压器模型的建立始终是一个重点和难点。其中变压器建模的困难之一是变压器的某些参数不仅是非线性的而且是随激磁电压的频率变化而变化的。铁心作为变压器电磁转换的核心部件,其在变压器建模过程中的重要性是不言而喻的。铁心受到时变电磁场的作用在叠片中产生涡流,涡流随着激磁电压的频率变化而呈现出频关效应(Frequency Dependent Effect)。其外在表现为电路模型中的电阻、电感、电容参数在频域内是随着频率而变化的,这就给时域仿真带来困难。在变压器电磁暂态研究中,由于涡流引起的频关效应是不容忽视的。
以往变压器建模中,对涡流的处理可以分为以下两种方法。
一是将铁心涡流直接等效为一个常数电阻,将该电阻与励磁支路并联。这种处理方法简单易行,往往应用于对涡流仿真要求不高或者频率较低的仿真中,用于宽频暂态仿真时误差较大。
二是从铁心叠片频变涡流阻抗表达式出发,运用电路综合的方法,用集总参数电路实现。这种方法充分考虑了叠片涡流的频关效应,并且可以方便的在感兴趣的频段内进行电路实现。
发明内容
本发明的目的是提出一种建立变压器铁心频变涡流模型的新方法,即用改 进的矢量匹配法对铁心的频域涡流阻抗表达式进行有理函数拟合,在此基础上用电路综合的方法实现涡流的频变电路模型。
本发明提出的采用改进的矢量匹配法建立铁心频变涡流模型的方法,包括以下步骤:
(1)对于变压器铁心涡流的频变效应,首先要得到涡流阻抗的频域表达式,应用电磁场理论可以推得铁心叠片中涡流的频变阻抗表达式:
Z(s)=Rlξtanhξ=R(s)+jωL(s)
其中,w为铁心叠片的宽度;I为磁通路径的长度;d为铁心叠片的厚度;N为绕组的匝数;μ为铁心叠片的磁导率;σ为铁心叠片的电导率;ω为励磁电压角频率;R为涡流阻抗的频变电阻;L为涡流阻抗的频变电感。
(2)采用改进的矢量匹配法对Z(s)进行有理函数逼近:
其中,pk、rk分别为函数的极点和留数,两者可以是实数或者共轭复数对。g和e为实数,根据具体情况进行选择。
(2-1)矢量匹配法(VF)通过求解下式的线性最小二乘问题来得到Z(s)的极点:
σ(s)Z(s)=p(s)
其中:
(2-2){am}是一组设定的初始极点。可以推得f(s)的极点一定和σ(s)的零点是相同的。并且{am}可以通过求解下式矩阵的特征值计算得到:
{qm}=eig(A-b·cT)
式中,A是包含初始极点{am}的对角阵;b是一个单位列向量;cT是含有留数的行向量。通过(2-1)、(2-2)进行迭代,用每次迭代过程中得到新的极点{qm}来替代之前的极点{am},进行下一次的迭代计算。这样在极点重置后通过2-3次迭代通常会收敛。极点求得之后,用得到的极点通过求解相应的最小二乘问题就能够得到所求留数。
改进的矢量匹配法(MVF)通过对原始的矢量匹配法(VF)进行改进,确保了网络函数的端口无源性,使得电路仿真的稳定性得到了保证。并且MVF只需要让初始极点在感兴趣的频段内线性分布或者对数分布就能达到非常精确的拟合结果,降低了对初始极点的要求。此外通过修正,在拟合对象含有噪声的情况下加快了拟合的收敛速度,使得拟合更加精确。
(3)用MVF对涡流的频变阻抗表达式拟合完成之后,得到的部分分式和就可以运用电路综合的方法进行电路实现。
拟合得到的有理函数均为实数极点,每一个极点和留数对应的电路为RL串联的支路形式,然后将每个支路并联就可以得到涡流的电路模型,其中第k个极点pk和留数rk对应的支路元件参数为:
本发明提出的一种建立变压器铁心频变涡流模型的新方法,与原有建立方法相比,有以下优点:
(1)本发明的建立方法,可以在一个很宽的频段内(10Hz-10MHz),只需要6阶电路就可以对涡流的频变阻抗公式进行十分精确的逼近,逼近最大误差不超过0.2%。而采用原有的建立方法逼近误差最高达到10%以上。具体对比可见下表(Cauer为原来的建立方法,MVF为本发明采用的方法):
(2)本发明建立的方法,在有理函数逼近过程中迭代收敛速度快,一般2-3次迭代后就会收敛为十分理想的范围。并且在迭代过程中对拟合初值的要求不高,数值稳定。在函数逼近上,与之前的方法有很大的优势。
(3)本发明建立的方法,在电路综合实现电路的过程中,能保证端口的无源性,就本例来说能保证每个元件都是正值。这就保证了涡流模型在电力系统仿真中计算的稳定性。
附图说明
图1为涡流频变阻抗MVF拟合结果图。(a)为频变电阻,(b)为频变电感。
图2为MVF拟合误差分布图。
图3为涡流6阶频变电路模型图。
具体实施方案
下面结合附图对本发明做进一步说明。
本发明提出的一种建立变压器铁心频变涡流模型的新方法,其具体实施方案如下:
(1)对于变压器铁心涡流的频变效应,首先要得到涡流阻抗的频域表达式,应用电磁场理论可以推得铁心叠片中涡流的频变阻抗表达式:
Z(s)=Rlξtanhξ=R(s)+jωL(s)
其中,w为铁心叠片的宽度;I为磁通路径的长度;d为铁心叠片的厚度;N为绕组的匝数;μ为铁心叠片的磁导率;σ为铁心叠片的电导率;ω为励磁电压角频率;R为涡流阻抗的频变电阻;L为涡流阻抗的频变电感。
(2)采用改进的矢量匹配法对Z(s)进行有理函数逼近:
其中,pk、rk分别为函数的极点和留数,两者可以是实数或者共轭复数对。g和e为实数,根据具体情况进行选择。
(2-1)改进的矢量匹配法(MVF)通过求解下式的线性最小二乘问题来得到Z(s)的极点:
σ(s)Z(s)=p(s)
其中:
为了避免方程组出现零值解,改进的矢量匹配法又针对最小二乘问题补充了一个判别项:
该判别项将σ(s)在给定的频率点上的实数部分之和限制为一个非零常数。当迭代在任一频率点收敛时,取此时σ(s)被强制收敛为1。
(2-2){am}是一组设定的初始极点。可以推得f(s)的极点一定和σ(s)的零点是相同的。并且{am}可以通过求解下式矩阵的特征值计算得到:
{qm}=eig(A-b·cT)
式中,A是包含初始极点{am}的对角阵;b是一个单位列向量;cT是含有留数的行向量。通过(2-1)、(2-2)进行迭代,用每次迭代过程中得到新的极点{qm}来替代之前的极点{am},进行下一次的迭代计算。这样在极点重置后通过2-3次迭代通常会收敛。极点求得之后,用得到的极点通过求解相应的最小二乘问题就能够得到所求留数。
拟合结果参见附图1.
拟合误差分布参见附图2.
(3)用MVF对涡流的频变阻抗表达式拟合完成之后,得到的部分分式和就可以运用电路综合的方法进行电路实现。
拟合得到的有理函数均为实数极点,每一个极点和留数对应的电路为RL串联的支路形式,然后将每个支路并联就可以得到涡流的电路模型,其中第k个极点pk和留数rk对应的支路元件参数为:
至此可以得到涡流频变阻抗的电路模型,参见附图3。
Claims (2)
1.本发明提出一种建立变压器铁心频变涡流模型的新方法,即用改进的矢量匹配法对铁心的频域涡流阻抗表达式进行有理函数拟合,在此基础上用电路综合的方法实现涡流的频变电路模型。
2.如权利要求1所述一种建立变压器铁心频变涡流模型的新方法,其特征在于该方法包括以下步骤:
(1)对于变压器铁心涡流的频变效应,首先要得到涡流阻抗的频域表达式,应用电磁场理论可以推得铁心叠片中涡流的频变阻抗表达式:
Z(s)=Rlξtanhξ=R(s)+jωL(s)
其中,w为铁心叠片的宽度;I为磁通路径的长度;d为铁心叠片的厚度;N为绕组的匝数;μ为铁心叠片的磁导率;σ为铁心叠片的电导率;ω为励磁电压角频率;R为涡流阻抗的频变电阻;L为涡流阻抗的频变电感;
(2)采用改进的矢量匹配法对Z(s)进行有理函数逼近:
其中,pk、rk分别为函数的极点和留数,两者可以是实数或者共轭复数对,g和e为实数,根据具体情况进行选择;
(2-1)矢量匹配法(VF)通过求解下式的线性最小二乘问题来得到Z(s)的极点:
σ(s)Z(s)=p(s)
其中:
(2-2){am}是一组设定的初始极点,可以推得f(s)的极点一定和σ(s)的零点是相同的,并且{am}可以通过求解下式矩阵的特征值计算得到:
{qm}=eig(A-b·cT)
式中,A是包含初始极点{am}的对角阵;b是一个单位列向量;cT是含有留数的行向量;通过(2-1)、(2-2)进行迭代,用每次迭代过程中得到新的极点{qm}来替代之前的极点{am},进行下一次的迭代计算;这样在极点重置后通过2-3次迭代通常会收敛;求得极点之后,用得到的极点通过求解相应的最小二乘问题就能够得到所求留数;
改进的矢量匹配法(MVF)通过对原始的矢量匹配法(VF)进行改进,确保了网络函数的端口无源性,使得电路仿真的稳定性得到了保证;并且MVF在函数逼近过程中,降低了对初始极点的要求;在拟合对象含有噪声的情况下加快了拟合的收敛速度,使得拟合更加精确;
(3)用MVF对涡流的频变阻抗表达式拟合完成之后,得到的部分分式和就可以运用电路综合的方法进行电路实现;
拟合得到的有理函数均为实数极点,每一个极点和留数对应的电路为RL串联的支路形式,然后将每个支路并联就可以得到涡流的电路模型,其中第k个极点pk和留数rk对应的支路元件参数为:
这样就得到一个二端口的集总参数电路模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410031708.3A CN104809322A (zh) | 2014-01-24 | 2014-01-24 | 一种建立变压器铁心频变涡流模型的新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410031708.3A CN104809322A (zh) | 2014-01-24 | 2014-01-24 | 一种建立变压器铁心频变涡流模型的新方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104809322A true CN104809322A (zh) | 2015-07-29 |
Family
ID=53694139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410031708.3A Pending CN104809322A (zh) | 2014-01-24 | 2014-01-24 | 一种建立变压器铁心频变涡流模型的新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104809322A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105954598A (zh) * | 2016-03-11 | 2016-09-21 | 中国石油大学(北京) | 一种铁芯电抗器参数识别方法 |
CN108763802A (zh) * | 2018-06-04 | 2018-11-06 | 天津农学院 | 计及电力传输线的电力电子换流系统电磁干扰估算方法 |
CN110188381A (zh) * | 2019-04-18 | 2019-08-30 | 中国北方车辆研究所 | 一种用于电磁干扰预测的仿真模型的构建方法及系统 |
CN113935172A (zh) * | 2021-10-18 | 2022-01-14 | 西南交通大学 | 一种端接频变负载的传输线系统瞬态响应分析方法 |
-
2014
- 2014-01-24 CN CN201410031708.3A patent/CN104809322A/zh active Pending
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105954598A (zh) * | 2016-03-11 | 2016-09-21 | 中国石油大学(北京) | 一种铁芯电抗器参数识别方法 |
CN105954598B (zh) * | 2016-03-11 | 2018-11-30 | 中国石油大学(北京) | 一种铁芯电抗器参数识别方法 |
CN108763802A (zh) * | 2018-06-04 | 2018-11-06 | 天津农学院 | 计及电力传输线的电力电子换流系统电磁干扰估算方法 |
CN110188381A (zh) * | 2019-04-18 | 2019-08-30 | 中国北方车辆研究所 | 一种用于电磁干扰预测的仿真模型的构建方法及系统 |
CN113935172A (zh) * | 2021-10-18 | 2022-01-14 | 西南交通大学 | 一种端接频变负载的传输线系统瞬态响应分析方法 |
CN113935172B (zh) * | 2021-10-18 | 2024-05-17 | 西南交通大学 | 一种端接频变负载的传输线系统瞬态响应分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Frequency-independent asymmetric double-$ pi $ equivalent circuit for on-chip spiral inductors: Physics-based modeling and parameter extraction | |
CN102609598B (zh) | 一种对大型电力系统进行电磁暂态仿真的方法 | |
CN104809322A (zh) | 一种建立变压器铁心频变涡流模型的新方法 | |
CN109815607B (zh) | 面向大规模交流电网基于时频坐标变换的全电磁暂态仿真方法 | |
Shimotani et al. | Equivalent-circuit generation from finite-element solution using proper orthogonal decomposition | |
CN105631105A (zh) | 一种含外加场耦合的线缆等效电路建模方法 | |
CN109086533B (zh) | 一种三相变压器直流偏磁频域计算方法 | |
CN105445541A (zh) | 一种任意频率下自适应功率计算方法 | |
CN101975892A (zh) | 一种换流阀用饱和电抗器的损耗计算方法 | |
Hu et al. | New approach for precisely measuring the zero-sequence parameters of long-distance double-circuit transmission lines | |
CN105678640B (zh) | 一种考虑变压器饱和影响的交流电网直流电流分布的预测方法 | |
CN113255285B (zh) | 基于ja模型与对偶原理的单相四柱式变压器的建模方法 | |
CN104933238A (zh) | 一种双绕组电力变压器非线性宽频模型及建立方法 | |
CN105467347B (zh) | 一种电压互感器励磁特性曲线求取方法 | |
CN114372433A (zh) | 一种考虑铁芯深度饱和特性的三相三柱变压器改进bctran仿真结构及方法 | |
CN110749799B (zh) | 一种特高压变压器直流偏磁等效试验方法及系统 | |
Sato et al. | Extended Cauer equivalent circuit model of inductors: Representing multi-resonant characteristics due to parasitic capacitance | |
Illia et al. | Finite element method based electromagnetic modeling of three-phase EMI filters | |
Preetha et al. | An equivalent circuit model for star-star-delta_utilized distribution transformer | |
CN104868874A (zh) | 一种-1/2阶rl分抗电路及其设计方法 | |
CN113190949B (zh) | 一种cvt宽频仿真模型的时域构建方法 | |
Mrad et al. | Conducted EMI prediction for integrated Class D audio amplifier | |
CN105957703B (zh) | 电子式互感器的数字积分方法和数字积分器 | |
Mukhtar et al. | Brune's synthesis of linear lossy distributed one-port and symmetric two-port microwave circuits | |
CN204948036U (zh) | 一种-1/2阶rl分抗电路 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150729 |
|
WD01 | Invention patent application deemed withdrawn after publication |