CN114242425B - 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 - Google Patents

一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 Download PDF

Info

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
Application number
CN202111231887.1A
Other languages
English (en)
Other versions
CN114242425A (zh
Inventor
周利军
李沃阳
陈家伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202111231887.1A priority Critical patent/CN114242425B/zh
Publication of CN114242425A publication Critical patent/CN114242425A/zh
Application granted granted Critical
Publication of CN114242425B publication Critical patent/CN114242425B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01FMAGNETS; INDUCTANCES; TRANSFORMERS; SELECTION OF MATERIALS FOR THEIR MAGNETIC PROPERTIES
    • H01F27/00Details of transformers or inductances, in general
    • H01F27/34Special means for preventing or reducing unwanted electric or magnetic effects, e.g. no-load losses, reactive currents, harmonics, oscillations, leakage fields
    • H01F27/341Preventing or reducing no-load losses or reactive currents
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex 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)由于电力电子设备的存在与铁心自身的非线性特征,导致励磁电流存在显著低次谐波分量,由于励磁电流函数满足狄利克雷充分条件,为分析其谐波性质对其进行傅里叶变换,展开式为:
Figure SMS_1
式中励磁电流Ih被分解为直流分量
Figure SMS_2
和相互正交的基波与各次谐波ancosnωt和bnsinnωt。由于供电系统电流不存在直流分量,并具有周期性质,在励磁电流一个周期内,可将上式简化为:
Figure SMS_3
式中In表征为励磁电流基波与各次谐波的幅值,In数值依据傅里叶分解性质求得:
Figure SMS_4
由于电网为平衡的三相系统,在平衡的三相系统中,偶次谐波互相抵消,可近似忽略励磁电流的偶次谐波,将其表达式简化为:
Figure SMS_5
式中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)由于牵引变压器铁心所采用的硅钢片是冷轧取向型,其在卷绕过程中无论是在心柱、铁轭还是拐角,都与导磁性能最佳的方向保持一致,可将全电流定律
Figure SMS_6
标量化为
Figure SMS_7
式中H为磁场强度、N为线圈匝数、L为牵引变压器铁心横截面的几何中心所在磁路的长度。
同时,考虑牵引变压器铁心各级横截面几何中心所在磁路长度不同,对铁心各级磁场强度分别计算,由此将励磁电流表达式代入得到牵引变压器铁心各级磁场强度表达式:
Figure SMS_8
式中,Hi(t)代表牵引变压器铁心第i级的磁场强度,a和b分别代牵引变压器铁心磁路长度和磁路宽度,Ri代表第i级磁路弧段半径;
3)由线圈电磁感应电动势公式有
Figure SMS_9
式中E1为变压器一次侧感应电压有效值,φm代表磁通最大值;由于变压器一次侧电压降较低,有E1≈U,式中U为变压器励磁电压有效值。联立以上二式,并将磁感强度与磁通关系
Figure SMS_10
代入,得:
Figure SMS_11
式中,Bm代表牵引变压器铁心磁通密度最大值,S代表磁路横截面积,w 和d分别代牵引变压器铁心各级宽度和硅钢片厚度;
4)由于牵引变压器铁心各级磁路的磁动势F=NI相等,冷轧取向型硅钢片具有单向导磁性且铁心各级间存在绝缘层,在磁路层面,可将铁心各级拓扑视为并联。并联磁路,磁通与磁阻Rm成反比,磁阻表达式为:
Figure SMS_12
式中μ为变压器铁心材质磁导率,由于各级磁路材质相同、横截面积相等,铁心各级磁通可视为与其几何中心所在磁路的长度成反比,即有:
Figure SMS_13
式中Bi(t)为铁心各级的磁通密度,R1代表第1级磁路弧段半径。由于第一级磁路最短、磁通最大,考虑主磁通滞后励磁电流相位
Figure SMS_14
第一级磁路磁通密度可表示为:
Figure SMS_15
进一步,将第一级磁路磁通密度表达式代入各级磁路磁通密度关系式中,得到牵引变压器各级磁通密度表达式:
Figure SMS_16
磁通密度仅有基频分量是因为磁通密度由励磁电压决定,变压器励磁电压一般为工频电压;
5)将上式代入电磁学对磁滞损耗的定义中,得到高湿度环境下计及铁心磁路分级的牵引变压器平均磁滞损耗P计算式:
Figure SMS_17
本发明的有益效果在于提出了一种更符合材料物性和运行工况的考虑磁路分级的牵引变压器磁滞损耗计算方式,能为牵引变压器生产优化设计和服役性能测评提供必要的数据保障。
附图说明
图1为本发明中所述牵引变压器铁心磁路分级示意图。
图2为本发明中所述牵引变压器铁心各级并联拓扑图。
具体实施方式
下面结合附图对本发明的实施流程作进一步的详述。由于电力电子设备的存在与铁心自身的非线性特征,导致励磁电流存在显著低次谐波分量,由于励磁电流函数满足狄利克雷充分条件,为分析其谐波性质对其进行傅里叶变换,展开式为:
Figure SMS_18
式中励磁电流Ih被分解为直流分量
Figure SMS_19
和相互正交的基波与各次谐波 ancosnωt和bnsinnωt。由于供电系统电流不存在直流分量,并具有周期性质,在励磁电流一个周期内,可将上式简化为:
Figure SMS_20
式中In表征为励磁电流基波与各次谐波的幅值,In数值依据傅里叶分解性质求得:
Figure SMS_21
由于电网为平衡的三相系统,在平衡的三相系统中,偶次谐波互相抵消,可近似忽略励磁电流的偶次谐波,将其表达式简化为:
Figure SMS_22
式中k∈{0,1,2,3…},由于谐波幅值与谐波次数呈反比,高次谐波幅值较小,计算仅考虑基波与三次谐波作用,故进一步简化表达式得:
Ih(t)≈I1cosωt+I3cos3ωt
考虑磁滞损耗计算需对磁场强度和磁通密度的乘积进行积分,磁场强度频率等于励磁电流频率,磁通密度频率等于励磁电压频率,由于励磁电压不畸变,恒为标准工频,且不同频率物理量乘积积分等于零,故在磁滞损耗计算中对于决定磁场强度的励磁电流可仅考虑其基波分量,励磁电流表达式进一步简化为 Ih(t)≈I1cosωt
式中,Ih(t)代表牵引变压器的励磁电流,I1、In分别代表励磁电流经傅里叶分解后的基波分量幅值和n次谐波分量幅值,ω为角频率,它满足:ω=2πf,f 为励磁频率,t为时间。
由于牵引变压器铁心所采用的硅钢片是冷轧取向型,其在卷绕过程中无论是在心柱、铁轭还是拐角,都与导磁性能最佳的方向保持一致,可将全电流定律
Figure SMS_23
标量化为
Figure SMS_24
式中H为磁场强度、N为线圈匝数、L为牵引变压器铁心横截面的几何中心所在磁路的长度。
同时,考虑牵引变压器铁心各级横截面几何中心所在磁路长度不同,对铁心各级磁场强度分别计算。
图1为本发明中所述牵引变压器铁心磁路分级示意图,图中以8级为例,由图可见,各级磁路均由四个矩形和四个四分之一圆构成,磁路长度可视为两倍铁心磁路长度、两倍铁心磁路宽度和圆周长之和,由此将励磁电流表达式代入得到牵引变压器铁心各级磁场强度表达式:
Figure SMS_25
式中,Hi(t)代表牵引变压器铁心第i级的磁场强度,a和b分别代牵引变压器铁心磁路长度和磁路宽度,Ri代表第i级磁路弧段半径。
由线圈电磁感应电动势公式有
Figure SMS_26
式中E1为变压器一次侧感应电压有效值,φm代表磁通最大值;由于变压器一次侧电压降较低,有 E1≈U,式中U为变压器励磁电压有效值。联立以上二式,并将磁感强度与磁通关系
Figure SMS_27
代入,得:
Figure SMS_28
式中,Bm代表牵引变压器铁心磁通密度最大值,S代表磁路横截面积,w 和d分别代牵引变压器铁心各级宽度和硅钢片厚度。
由于牵引变压器铁心各级磁路的磁动势F=NI相等,冷轧取向型硅钢片具有单向导磁性且铁心各级间存在绝缘层,在磁路层面,可将铁心各级拓扑视为并联。图2为本发明中所述牵引变压器铁心各级并联拓扑图,由图可知并联磁路,磁通与磁阻Rm成反比,磁阻表达式为:
Figure SMS_29
式中μ为变压器铁心材质磁导率,由于各级磁路材质相同、横截面积相等,铁心各级磁通可视为与其几何中心所在磁路的长度成反比,即有:
Figure SMS_30
式中Bi(t)为铁心各级的磁通密度,R1代表第1级磁路弧段半径。由于第一级磁路最短、磁通最大,考虑主磁通滞后励磁电流相位
Figure SMS_31
第一级磁路磁通密度可表示为:
Figure SMS_32
进一步,将第一级磁路磁通密度表达式代入各级磁路磁通密度关系式中,得到牵引变压器各级磁通密度表达式:
Figure SMS_33
磁通密度仅有基频分量是因为磁通密度由励磁电压决定,变压器励磁电压一般为工频电压。
将上式代入电磁学对磁滞损耗的定义中,得到计及铁心磁路分级的牵引变压器平均磁滞损耗P计算式:
Figure SMS_34
本发明的有益效果在于提出了一种更符合材料物性和运行工况的考虑磁路分级的牵引变压器磁滞损耗计算方式,能为牵引变压器生产优化设计和服役性能测评提供必要的数据保障。

Claims (1)

1.一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法,其特征在于,铁心材质为高导磁冷轧晶粒取向硅钢片,包括以下步骤:
1)依据磁滞计算原则,决定磁场强度的励磁电流仅考虑基波分量且非同频物理量周期积分为零,从而得到励磁电流表达式:
Figure FDA0004119680530000011
式中,Ih(t)为励磁电流,I1、In分别为励磁电流经傅里叶分解后的基波分量幅值和n次分量幅值,ω为角频率,它满足:ω=2πf,f为励磁频率,t为时间;
2)根据全电流定律,考虑磁路分级,牵引变压器铁心各级磁场强度表达为:
Figure FDA0004119680530000012
式中,Hi(t)代表牵引变压器铁心第i级的磁场强度,N代表励磁绕组线圈总匝数,L为牵引变压器铁心横截面的几何中心所在磁路的长度,a和b分别代牵引变压器铁心磁路长度和磁路宽度,Ri代表第i级磁路弧段半径;
3)根据线圈感应电动势公式及变压器一次侧感应电动势与励磁电压的约等关系,得到牵引变压器铁心内磁通密度最大值的表达式:
Figure FDA0004119680530000013
式中,Bm代表牵引变压器铁心磁通密度最大值,φm代表磁通最大值,S代表磁路横截面积,U为变压器励磁电压有效值,w和d分别代牵引变压器铁心各级宽度和硅钢片厚度;
4)根据磁路各级并联拓扑及磁阻比例关系,得到铁心各级磁通密度计算式:
Figure FDA0004119680530000014
式中Bi(t)为铁心各级的磁通密度,R1代表第1级磁路弧段半径,
Figure FDA0004119680530000016
为主磁通滞后励磁电流相位;
5)根据(2)(4)及电磁学对磁滞损耗的定义,得到计及铁心磁路分级的牵引变压器平均磁滞损耗P计算式:
Figure FDA0004119680530000015
CN202111231887.1A 2021-10-22 2021-10-22 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 Active CN114242425B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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 西南交通大学 一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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