CN114242425B - 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 - Google Patents
一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 Download PDFInfo
- Publication number
- CN114242425B CN114242425B CN202111231887.1A CN202111231887A CN114242425B CN 114242425 B CN114242425 B CN 114242425B CN 202111231887 A CN202111231887 A CN 202111231887A CN 114242425 B CN114242425 B CN 114242425B
- Authority
- CN
- China
- Prior art keywords
- magnetic
- traction transformer
- magnetic circuit
- core
- iron 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.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01F—MAGNETS; INDUCTANCES; TRANSFORMERS; SELECTION OF MATERIALS FOR THEIR MAGNETIC PROPERTIES
- H01F27/00—Details of transformers or inductances, in general
- H01F27/34—Special means for preventing or reducing unwanted electric or magnetic effects, e.g. no-load losses, reactive currents, harmonics, oscillations, leakage fields
- H01F27/341—Preventing or reducing no-load losses or reactive currents
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Analysis (AREA)
- Power Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Regulation Of General Use Transformers (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法,通过构建基于非线性媒质关系的硅钢片电磁感应方程,并引入非同频物理量功率积分可忽略的特性,针对铁心各级磁路磁效应互异的特点,提出了适用于磁铁心磁路分级的牵引变压器铁心磁滞损耗的简化计算模型。本发明的有益效果在于有助于提出一种更符合材料物性和运行工况的磁滞损耗评估方法,能为牵引变压器生产优化设计和服役性能测评提供必要的数据保障。
Description
技术领域
本发明属于电气设备电磁分析与数值计算领域,具体涉及一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法。
背景技术
牵引变压器作为牵引供电系统中的关键设备,具备短时冲击负荷大、空载时间长的运行特点,对其铁心能耗的评估与优化研究具有重要的工程价值。而磁滞损耗作为铁心能耗的重要组成部分,为对牵引变压器铁心损耗进行进一步研究,提出一种足够精确的牵引变压器磁滞损耗求解方法具有迫切的工程意义。
永久磁铁、铁磁性材料,以及电磁铁中,磁通经过的闭合路径叫做磁路,磁路分析的主要目的是要确定励磁磁通势和它所产生的磁通的关系。由于变压器卷绕工艺,牵引变压器各级铁心几何尺寸不同,将造成铁心各级磁路磁阻不同,进而使得变压器铁心各级磁场强度、磁通密度分布不均。在传统计算公式中,磁滞损耗计算通常将铁心视为均一化整体,数值上正比于平均磁场强度和平均磁通密度,这种公式无法解释由于磁路分级造成的磁场分布不均,计算误差较大且无法精确描述铁心某点的磁滞损耗,不能满足牵引变压器发展对损耗计算精度更高的要求。为此,提出一种考虑牵引变压器磁路分级的磁滞损耗计算公式显得尤为重要。
发明内容
本发明的目的是提供一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法,并通过如下技术手段实现:
1)由于电力电子设备的存在与铁心自身的非线性特征,导致励磁电流存在显著低次谐波分量,由于励磁电流函数满足狄利克雷充分条件,为分析其谐波性质对其进行傅里叶变换,展开式为:
式中In表征为励磁电流基波与各次谐波的幅值,In数值依据傅里叶分解性质求得:
由于电网为平衡的三相系统,在平衡的三相系统中,偶次谐波互相抵消,可近似忽略励磁电流的偶次谐波,将其表达式简化为:
式中k∈{0,1,2,3…},由于谐波幅值与谐波次数呈反比,高次谐波幅值较小,计算仅考虑基波与三次谐波作用,故进一步简化表达式得:
Ih(t)≈I1cosωt+I3cos3ωt
考虑磁滞损耗计算需对磁场强度和磁通密度的乘积进行积分,磁场强度频率等于励磁电流频率,磁通密度频率等于励磁电压频率,由于励磁电压不畸变,恒为标准工频,且不同频率物理量乘积积分等于零,故在磁滞损耗计算中对于决定磁场强度的励磁电流可仅考虑其基波分量,励磁电流表达式进一步简化为:
Ih(t)≈I1cosωt
式中,Ih(t)代表牵引变压器的励磁电流,I1、In分别代表励磁电流经傅里叶分解后的基波分量幅值和n次谐波分量幅值,ω为角频率,它满足:ω=2πf,f 为励磁频率,t为时间;
2)由于牵引变压器铁心所采用的硅钢片是冷轧取向型,其在卷绕过程中无论是在心柱、铁轭还是拐角,都与导磁性能最佳的方向保持一致,可将全电流定律标量化为式中H为磁场强度、N为线圈匝数、L为牵引变压器铁心横截面的几何中心所在磁路的长度。
同时,考虑牵引变压器铁心各级横截面几何中心所在磁路长度不同,对铁心各级磁场强度分别计算,由此将励磁电流表达式代入得到牵引变压器铁心各级磁场强度表达式:
式中,Hi(t)代表牵引变压器铁心第i级的磁场强度,a和b分别代牵引变压器铁心磁路长度和磁路宽度,Ri代表第i级磁路弧段半径;
3)由线圈电磁感应电动势公式有式中E1为变压器一次侧感应电压有效值,φm代表磁通最大值;由于变压器一次侧电压降较低,有E1≈U,式中U为变压器励磁电压有效值。联立以上二式,并将磁感强度与磁通关系代入,得:
式中,Bm代表牵引变压器铁心磁通密度最大值,S代表磁路横截面积,w 和d分别代牵引变压器铁心各级宽度和硅钢片厚度;
4)由于牵引变压器铁心各级磁路的磁动势F=NI相等,冷轧取向型硅钢片具有单向导磁性且铁心各级间存在绝缘层,在磁路层面,可将铁心各级拓扑视为并联。并联磁路,磁通与磁阻Rm成反比,磁阻表达式为:
式中μ为变压器铁心材质磁导率,由于各级磁路材质相同、横截面积相等,铁心各级磁通可视为与其几何中心所在磁路的长度成反比,即有:
进一步,将第一级磁路磁通密度表达式代入各级磁路磁通密度关系式中,得到牵引变压器各级磁通密度表达式:
磁通密度仅有基频分量是因为磁通密度由励磁电压决定,变压器励磁电压一般为工频电压;
5)将上式代入电磁学对磁滞损耗的定义中,得到高湿度环境下计及铁心磁路分级的牵引变压器平均磁滞损耗P计算式:
本发明的有益效果在于提出了一种更符合材料物性和运行工况的考虑磁路分级的牵引变压器磁滞损耗计算方式,能为牵引变压器生产优化设计和服役性能测评提供必要的数据保障。
附图说明
图1为本发明中所述牵引变压器铁心磁路分级示意图。
图2为本发明中所述牵引变压器铁心各级并联拓扑图。
具体实施方式
下面结合附图对本发明的实施流程作进一步的详述。由于电力电子设备的存在与铁心自身的非线性特征,导致励磁电流存在显著低次谐波分量,由于励磁电流函数满足狄利克雷充分条件,为分析其谐波性质对其进行傅里叶变换,展开式为:
式中In表征为励磁电流基波与各次谐波的幅值,In数值依据傅里叶分解性质求得:
由于电网为平衡的三相系统,在平衡的三相系统中,偶次谐波互相抵消,可近似忽略励磁电流的偶次谐波,将其表达式简化为:
式中k∈{0,1,2,3…},由于谐波幅值与谐波次数呈反比,高次谐波幅值较小,计算仅考虑基波与三次谐波作用,故进一步简化表达式得:
Ih(t)≈I1cosωt+I3cos3ωt
考虑磁滞损耗计算需对磁场强度和磁通密度的乘积进行积分,磁场强度频率等于励磁电流频率,磁通密度频率等于励磁电压频率,由于励磁电压不畸变,恒为标准工频,且不同频率物理量乘积积分等于零,故在磁滞损耗计算中对于决定磁场强度的励磁电流可仅考虑其基波分量,励磁电流表达式进一步简化为 Ih(t)≈I1cosωt
式中,Ih(t)代表牵引变压器的励磁电流,I1、In分别代表励磁电流经傅里叶分解后的基波分量幅值和n次谐波分量幅值,ω为角频率,它满足:ω=2πf,f 为励磁频率,t为时间。
由于牵引变压器铁心所采用的硅钢片是冷轧取向型,其在卷绕过程中无论是在心柱、铁轭还是拐角,都与导磁性能最佳的方向保持一致,可将全电流定律标量化为式中H为磁场强度、N为线圈匝数、L为牵引变压器铁心横截面的几何中心所在磁路的长度。
同时,考虑牵引变压器铁心各级横截面几何中心所在磁路长度不同,对铁心各级磁场强度分别计算。
图1为本发明中所述牵引变压器铁心磁路分级示意图,图中以8级为例,由图可见,各级磁路均由四个矩形和四个四分之一圆构成,磁路长度可视为两倍铁心磁路长度、两倍铁心磁路宽度和圆周长之和,由此将励磁电流表达式代入得到牵引变压器铁心各级磁场强度表达式:
式中,Hi(t)代表牵引变压器铁心第i级的磁场强度,a和b分别代牵引变压器铁心磁路长度和磁路宽度,Ri代表第i级磁路弧段半径。
由线圈电磁感应电动势公式有式中E1为变压器一次侧感应电压有效值,φm代表磁通最大值;由于变压器一次侧电压降较低,有 E1≈U,式中U为变压器励磁电压有效值。联立以上二式,并将磁感强度与磁通关系代入,得:
式中,Bm代表牵引变压器铁心磁通密度最大值,S代表磁路横截面积,w 和d分别代牵引变压器铁心各级宽度和硅钢片厚度。
由于牵引变压器铁心各级磁路的磁动势F=NI相等,冷轧取向型硅钢片具有单向导磁性且铁心各级间存在绝缘层,在磁路层面,可将铁心各级拓扑视为并联。图2为本发明中所述牵引变压器铁心各级并联拓扑图,由图可知并联磁路,磁通与磁阻Rm成反比,磁阻表达式为:
式中μ为变压器铁心材质磁导率,由于各级磁路材质相同、横截面积相等,铁心各级磁通可视为与其几何中心所在磁路的长度成反比,即有:
进一步,将第一级磁路磁通密度表达式代入各级磁路磁通密度关系式中,得到牵引变压器各级磁通密度表达式:
磁通密度仅有基频分量是因为磁通密度由励磁电压决定,变压器励磁电压一般为工频电压。
将上式代入电磁学对磁滞损耗的定义中,得到计及铁心磁路分级的牵引变压器平均磁滞损耗P计算式:
本发明的有益效果在于提出了一种更符合材料物性和运行工况的考虑磁路分级的牵引变压器磁滞损耗计算方式,能为牵引变压器生产优化设计和服役性能测评提供必要的数据保障。
Claims (1)
1.一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法,其特征在于,铁心材质为高导磁冷轧晶粒取向硅钢片,包括以下步骤:
1)依据磁滞计算原则,决定磁场强度的励磁电流仅考虑基波分量且非同频物理量周期积分为零,从而得到励磁电流表达式:
式中,Ih(t)为励磁电流,I1、In分别为励磁电流经傅里叶分解后的基波分量幅值和n次分量幅值,ω为角频率,它满足:ω=2πf,f为励磁频率,t为时间;
2)根据全电流定律,考虑磁路分级,牵引变压器铁心各级磁场强度表达为:
式中,Hi(t)代表牵引变压器铁心第i级的磁场强度,N代表励磁绕组线圈总匝数,L为牵引变压器铁心横截面的几何中心所在磁路的长度,a和b分别代牵引变压器铁心磁路长度和磁路宽度,Ri代表第i级磁路弧段半径;
3)根据线圈感应电动势公式及变压器一次侧感应电动势与励磁电压的约等关系,得到牵引变压器铁心内磁通密度最大值的表达式:
式中,Bm代表牵引变压器铁心磁通密度最大值,φm代表磁通最大值,S代表磁路横截面积,U为变压器励磁电压有效值,w和d分别代牵引变压器铁心各级宽度和硅钢片厚度;
4)根据磁路各级并联拓扑及磁阻比例关系,得到铁心各级磁通密度计算式:
5)根据(2)(4)及电磁学对磁滞损耗的定义,得到计及铁心磁路分级的牵引变压器平均磁滞损耗P计算式:
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111231887.1A CN114242425B (zh) | 2021-10-22 | 2021-10-22 | 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111231887.1A CN114242425B (zh) | 2021-10-22 | 2021-10-22 | 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114242425A CN114242425A (zh) | 2022-03-25 |
| CN114242425B true CN114242425B (zh) | 2023-06-02 |
Family
ID=80743199
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111231887.1A Active CN114242425B (zh) | 2021-10-22 | 2021-10-22 | 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114242425B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115656025B (zh) * | 2022-11-22 | 2024-03-12 | 西南交通大学 | 一种海上变压器铁心抗腐蚀能力的评估方法 |
| CN116127240B (zh) * | 2022-11-22 | 2023-12-05 | 西南交通大学 | 一种牵引变压器卷铁心过负荷能力的评估方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108647382A (zh) * | 2018-04-02 | 2018-10-12 | 西南交通大学 | 一种卷铁心牵引变压器涡流损耗的测评方法 |
| CN110399693A (zh) * | 2019-07-31 | 2019-11-01 | 西南交通大学 | 一种片间短路下变压器卷铁心涡流损耗的计算方法 |
| JP2019204269A (ja) * | 2018-05-23 | 2019-11-28 | 日本製鉄株式会社 | 電磁場解析システム、電磁場解析方法、およびプログラム |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6873239B2 (en) * | 2002-11-01 | 2005-03-29 | Metglas Inc. | Bulk laminated amorphous metal inductive device |
| EP2320439A4 (en) * | 2008-09-03 | 2018-05-02 | Hitachi Industrial Equipment Systems Co., Ltd. | Wound iron core for static apparatus, amorphous transformer and coil winding frame for transformer |
| JP2017040002A (ja) * | 2015-08-21 | 2017-02-23 | 日本金属株式会社 | 高周波用無方向性電磁鋼板およびその製造方法 |
| JP2018084940A (ja) * | 2016-11-23 | 2018-05-31 | 吉川工業株式会社 | 鉄心の設計方法、製造方法及び設計プログラム |
| DE102017205004B4 (de) * | 2017-03-24 | 2022-12-08 | Siemens Aktiengesellschaft | Leistungsschalter |
| CN111270064B (zh) * | 2020-01-23 | 2020-10-02 | 西南交通大学 | 一种卷铁心退火工艺的控制方法及退火效果的测评方法 |
| CN112016204B (zh) * | 2020-08-27 | 2022-04-26 | 西南交通大学 | 一种考虑自耦变压器卷铁心材料非线性磁性能的涡流损耗分析方法 |
| CN111931310B (zh) * | 2020-08-28 | 2021-08-13 | 西南交通大学 | 一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 |
-
2021
- 2021-10-22 CN CN202111231887.1A patent/CN114242425B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108647382A (zh) * | 2018-04-02 | 2018-10-12 | 西南交通大学 | 一种卷铁心牵引变压器涡流损耗的测评方法 |
| JP2019204269A (ja) * | 2018-05-23 | 2019-11-28 | 日本製鉄株式会社 | 電磁場解析システム、電磁場解析方法、およびプログラム |
| CN110399693A (zh) * | 2019-07-31 | 2019-11-01 | 西南交通大学 | 一种片间短路下变压器卷铁心涡流损耗的计算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 高洁,杨富尧等.0.23 mm薄规格取向硅钢谐波损耗与磁致伸缩研究.功能材料.2019,12132-12136. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114242425A (zh) | 2022-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Rezaei-Zare | Enhanced transformer model for low-and mid-frequency transients—Part I: Model development | |
| CN110008645A (zh) | 一种变压器损耗计算方法 | |
| CN112632458B (zh) | 一种计及铁心磁饱和的牵引变压器涡流损耗求解方法 | |
| CN114242425B (zh) | 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 | |
| Najafi et al. | Comparison of core loss and magnetic flux distribution in amorphous and silicon steel core transformers | |
| Wang et al. | Estimation model of core loss under DC bias | |
| CN111931310A (zh) | 一种考虑相异磁场边界条件的卷铁心层间短路涡流损耗评估方法 | |
| Li et al. | Analytical loss model for magnetic cores based on vector magnetic circuit theory | |
| Magdaleno-Adame et al. | Electromagnetic finite element analysis of electrical steels combinations in lamination core steps of single-phase distribution transformers | |
| Li et al. | Numerical model for eddy-current loss of wound core in single-phase transformer | |
| Hane et al. | Dynamic hysteresis modeling taking skin effect into account for magnetic circuit analysis and validation for various core materials | |
| Zhao et al. | A dynamic model for prediction of hysteresis and loss property of soft material under multiple harmonic magnetizations | |
| Wang et al. | A higher-order loss-separation model for fast estimation of core loss in high-frequency transformers | |
| Cheng et al. | Engineering-oriented benchmarking of problem 21 family and experimental verification | |
| Dawood et al. | Numerical and experimental comparison of the no-load losses in the different grain-oriented electrical steel materials | |
| Potgieter et al. | Effects observed in the characterization of soft magnetic composite for high frequency, high flux density applications | |
| Qin et al. | Compatibility analysis among vector magnetic circuit theory, electrical circuit theory, and electromagnetic field theory | |
| Wang et al. | Magnetic shunt design and their effects on transformer winding electromagnetic forces | |
| Deželak et al. | Usage of a simplified and Jiles–Atherton model when accounting for the hysteresis losses within a welding transformer | |
| CN114168887B (zh) | 一种计及铁心磁路分级的牵引变压器涡流损耗求解方法 | |
| Li et al. | Core form transformer topological duality based transient model validation for exciting current calculation within DC bias | |
| Carrander | Magnetizing currents in power transformers: Measurements, simulations, and diagnostic methods | |
| Zirka et al. | Modelling losses in electrical steel laminations | |
| Elhaminia et al. | Magnetic flux density determination in 3D wound core transformer using H-balance equation | |
| Moses et al. | Power loss of non oriented electrical steel under square wave excitation |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |







































