CN111931310B - 一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 - Google Patents
一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 Download PDFInfo
- Publication number
- CN111931310B CN111931310B CN202010881786.8A CN202010881786A CN111931310B CN 111931310 B CN111931310 B CN 111931310B CN 202010881786 A CN202010881786 A CN 202010881786A CN 111931310 B CN111931310 B CN 111931310B
- Authority
- CN
- China
- Prior art keywords
- eddy current
- short circuit
- wound core
- magnetic field
- layers
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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/15—Correlation function computation including computation of convolution operations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明公开了一种考虑相异磁场边界条件的卷铁心层间短路涡流损耗评估方法:在给定卷铁心结构、尺寸、材料电磁特性参数的条件下,构建横截面涡流场计算的基本方程,同时考虑了不同层级硅钢片边缘磁场强度的非均匀分布,改进了涡流及其损耗求解的边界条件,得到了卷铁心层间短路区域内的涡流密度函数解析式,以及涡流损耗的计算公式。本发明的有益效果是提出了一种对卷铁心层间短路故障区域的涡流损耗进行精准评估的方法,为变压器生产或制造阶段的性能评估和缺陷排查提供了理论基础。
Description
技术领域
本发明属于电气设备电磁分析与数值计算领域,具体涉及一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法。
背景技术
牵引供电系统中的变压器往往具备短时冲击负荷、空载时间长的运行工况,因而铁心能耗的评估与优化研究具有重要的工程价值。而边缘毛刺或绝缘劣化等引发的牵引变压器卷铁心层级短路是变压器最常见的故障之一,这会导致本来独立于各层级硅钢片内部的涡流将会跃迁到不同层级,形成更大的回环并造成可观的损耗。因此,提出一套针对卷铁心存在局部层间短路时的涡流损耗评估公式,是服务于变压器优化设计和故障诊断的重要前置性研究。
目前针对层间短路故障的涡流损耗评估方案,是在公式推导过程中假定各级硅钢片的边界具有一致的磁场分布。然而对于真实变压器的卷铁心,各层级的磁路长度并不相等,因此铁心内磁场并不是严格的均匀分布。在发生层间短路时,短路区域所在的数片硅钢片形成一个整体的涡流回路,这个回路在边界处的磁场强度显然也是不同的。因此应该将短路区域始末两端所在的硅钢片层级的磁场强度作为短路计算的边界条件,否则将会出现明显的误差。
发明内容
本发明的目的是提供一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法,并通过如下技术手段实现:
1)本案卷铁心纵截面形态为圆角矩形,卷绕过程会给定以下参数:心柱长度(L)、铁轭长度(Y)、圆角半径(R),硅钢带厚度(d)。对卷绕的首级而言,以硅钢带厚度方向的中线为卷绕路径的长度标准,其圆角半径会在骨架的尺寸基础上增加硅钢片厚度的一半:
后面的层级将有所区别,卷绕堆砌后的半径数值将比上一级增加硅钢片厚度的整长,下面给出第二级、第三级、第四级圆角半径的计算公式:
由此可以给出不同层级的圆角半径数值的递推关系,即得到卷铁心第i级的圆角半径Ri的表达式:
进一步能够得到卷铁心第i级的整个卷绕路径长度li的表达式:
式中,i∈{1,2,3,…,m,m+1,m+2,…,2m},2m为卷铁心的总级数;
由于变压器铁心都是正弦型交流励磁,因此应对磁准静态场的麦克斯韦方程进行求解,同时考虑硅钢片的晶粒取向特性和卷铁心各层级的宽度远大于厚度的特征,得到卷铁心层间短路区域内的涡流场复数形式二阶微分方程:
根据二阶线性齐次常微分方程的基础算法,得到上述方程的通解:
式中,C1和C2为待定常数,通过涡流场的边界条件确定。
2)由于卷铁心所采用的硅钢片是冷轧取向型,其在卷绕过程中无论是在心柱、铁轭还是拐角,都与导磁性能最佳的方向保持一致。相对地,在垂直于叠片的方向,硅钢片的磁性能最差,再加之有绝缘涂层的存在,非短路区域的磁力线几乎不会发生跃迁。结合这种“层间并联,层内串联”的分级磁路模型,可以认为各层级的边界磁场强度与整个卷铁心的平均磁场强度的比值等于该层级的磁路长度与整个卷铁心的平均磁路长度的比值,因此可以计算出卷铁心第i级的边界磁场强度:
式中,l0为卷铁心横截面的几何中心所在磁路的长度,Havg是卷铁心横截面的磁场强度平均值,它们由下式确定:
l0=2(L+Y)+2π(R+md)
Havg=Urms/(4.44μfNA)
式中,Urms为变压器绕组一次侧电压的有效值,N为励磁绕组匝数,A为卷铁心的横截面积。
当卷铁心发生层间短路时,本来在各个层级内部流动的涡流,会跃迁到其它层级,在整个短路区域出现大环流,但是由于这种故障往往是小范围的,短路片数一般不会太多,因此仍然可以认为卷铁心的横截面宽度远大于短路区域总厚度。此外,根据分级磁路模型的分析可知,卷铁心不同层级具有相异的边界磁场强度,因此当出现局部层间短路故障时,涡流场计算的边界条件应该采用该短路区域首末层各自级所在的边界磁场强度。设卷铁心在第s级发生层间短路,短路片数为k,则求解短路区域内部各点的磁场强度所需的边界条件为:
Hz(0)=Hs,Hz(D)=Hs+k
式中,s,k∈{1,2,3,…,m,m+1,m+2,…,2m},Hs、Hs+k分别代表卷铁心第s、s+k级的边界磁场强度,D为短路区域总厚度,它满足:D=kd,d为硅钢片厚度。
将边界条件带入磁场强度的通解,即可得到待定常数C1和C2的计算式:
3)根据安培环路定律的微分形式,对磁场强度函数进行求导并取负,得到卷铁心层间短路区域内涡流密度函数Jy(x)的表达式:
4)由于卷铁心发生层间短路时,涡流会扩散到整个故障区域,因此涡流损耗的计算域也需要调整为整个短路区域的总厚度D。根据电磁学对有功功耗的定义,得到卷铁心层间短路区域的涡流损耗计算公式:
本发明的有益效果是提出了一种对卷铁心层间短路故障区域涡流损耗进行精准评估的方法,为变压器制造阶段的性能评估和缺陷排查提供了理论基础。
附图说明
图1为本发明中卷铁心整体结构正视图。
图2为本发明中卷铁心局部层间短路示意图。
具体实施方式
下面结合附图对本发明的实施流程作进一步的详述。
图1为本发明中所述的卷铁心整体结构正视图,其卷铁心纵截面形态为圆角矩形,卷绕过程会给定以下参数:心柱长度(L)、铁轭长度(Y)、圆角半径(R),硅钢带厚度(d)。对卷绕的首级而言,以硅钢带厚度方向的中线为卷绕路径的长度标准,其圆角半径会在骨架的尺寸基础上增加硅钢片厚度的一半:
后面的层级将有所区别,卷绕堆砌后的半径数值将比上一级增加硅钢片厚度的整长,下面给出第二级、第三级、第四级圆角半径的计算公式:
由此可以给出不同层级的圆角半径数值的递推关系,即得到卷铁心第i级的圆角半径Ri的表达式:
进一步能够得到卷铁心第i级的整个卷绕路径长度li的表达式:
式中,i∈{1,2,3,…,m,m+1,m+2,…,2m},2m为卷铁心的总级数;
由于变压器铁心都是正弦型交流励磁,因此应对磁准静态场的麦克斯韦方程进行求解,同时考虑硅钢片的晶粒取向特性和卷铁心各层级的宽度远大于厚度的特征,得到卷铁心层间短路区域内的涡流场复数形式二阶微分方程:
根据二阶线性齐次常微分方程的基础算法,得到上述方程的通解:
式中,C1和C2为待定常数,通过涡流场的边界条件确定。
由于卷铁心所采用的硅钢片是冷轧取向型,其在卷绕过程中无论是在心柱、铁轭还是拐角,都与导磁性能最佳的方向保持一致。相对地,在垂直于叠片的方向,硅钢片的磁性能最差,再加之有绝缘涂层的存在,非短路区域的磁力线几乎不会发生跃迁。结合这种“层间并联,层内串联”的分级磁路模型,可以认为各层级的边界磁场强度与整个卷铁心的平均磁场强度的比值等于该层级的磁路长度与整个卷铁心的平均磁路长度的比值,因此可以计算出卷铁心第i级的边界磁场强度:
式中,l0为卷铁心横截面的几何中心所在磁路的长度,Havg是卷铁心横截面的磁场强度平均值,它们由下式确定:
l0=2(L+Y)+2π(R+md)
Havg=Urms/(4.44μfNA)
式中,Urms为变压器绕组一次侧电压的有效值,N为励磁绕组匝数,A为卷铁心的横截面积。
图2为本发明中卷铁心局部层间短路示意图。当卷铁心发生层间短路时,本来在各个层级内部流动的涡流,会跃迁到其它层级,在整个短路区域出现大环流,但是由于这种故障往往是小范围的,短路片数一般不会太多,因此仍然可以认为卷铁心的横截面宽度远大于短路区域总厚度。此外,根据分级磁路模型的分析可知,卷铁心不同层级具有相异的边界磁场强度,因此当出现局部层间短路故障时,涡流场计算的边界条件应该采用该短路区域首末层各自级所在的边界磁场强度。设卷铁心在第s级发生层间短路,短路片数为k,则求解短路区域内部各点的磁场强度所需的边界条件为:
Hz(0)=Hs,Hz(D)=Hs+k
式中,s,k∈{1,2,3,…,m,m+1,m+2,…,2m},Hs、Hs+k分别代表卷铁心第s、s+k级的边界磁场强度,D为短路区域总厚度,它满足:D=kd,d为硅钢片厚度。
将边界条件带入磁场强度的通解,即可得到待定常数C1和C2的计算式:
根据安培环路定律的微分形式,对磁场强度函数进行求导并取负,得到卷铁心层间短路区域内涡流密度函数Jy(x)的表达式:
由于卷铁心发生层间短路时,涡流会扩散到整个故障区域,因此涡流损耗的计算域也需要调整为整个短路区域的总厚度D。根据电磁学对有功功耗的定义,得到卷铁心层间短路区域的涡流损耗计算公式:
Claims (1)
1.一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法,其特征在于,铁心材质为高导磁冷轧晶粒取向硅钢片,包括以下步骤:
1)得到卷铁心层间短路区域内的涡流场复数形式微分方程:
根据二阶线性齐次常微分方程的基础算法,得到(1)的通解:
式中,C1和C2为待定常数,通过涡流场的边界条件确定;
2)设卷铁心在第s级发生层间短路,短路片数为k,则求解短路区域内部各点的磁场强度所需的边界条件为:
Hz(0)=Hs,Hz(D)=Hs+k (3)
式中,s,k∈{1,2,3,…,m,m+1,m+2,…,2m},2m为卷铁心的总层级数,Hs、Hs+k分别代表卷铁心第s、s+k级的边界磁场强度,D为短路区域总厚度,它满足:D=kd,d为硅钢片厚度;
将(3)带入(2)即可得到待定常数C1和C2的计算式:
3)根据安培环路定律的微分形式,以及(2)、(4)对磁场强度的计算结果,得到卷铁心层间短路区域内涡流密度函数Jy(x)的表达式:
4)根据电磁学对有功功耗的定义,得到卷铁心层间短路区域的涡流损耗:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010881786.8A CN111931310B (zh) | 2020-08-28 | 2020-08-28 | 一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010881786.8A CN111931310B (zh) | 2020-08-28 | 2020-08-28 | 一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111931310A CN111931310A (zh) | 2020-11-13 |
CN111931310B true CN111931310B (zh) | 2021-08-13 |
Family
ID=73308741
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010881786.8A Active CN111931310B (zh) | 2020-08-28 | 2020-08-28 | 一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111931310B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112560301B (zh) * | 2020-11-26 | 2021-10-08 | 东南大学 | 一种基于磁感原理的磁性材料涡流损耗计算方法 |
CN114242425B (zh) * | 2021-10-22 | 2023-06-02 | 西南交通大学 | 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 |
CN114168887B (zh) * | 2021-10-25 | 2022-07-26 | 西南交通大学 | 一种计及铁心磁路分级的牵引变压器涡流损耗求解方法 |
CN116127240B (zh) * | 2022-11-22 | 2023-12-05 | 西南交通大学 | 一种牵引变压器卷铁心过负荷能力的评估方法 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4759949A (en) * | 1987-07-23 | 1988-07-26 | Westinghouse Electric Corp. | Method of insulating ferromagnetic amorphous metal continuous strip |
US7506280B2 (en) * | 2004-11-12 | 2009-03-17 | Tabtronics, Inc. | Magnetic winding and method of making same |
JP4922191B2 (ja) * | 2008-01-09 | 2012-04-25 | Tdkラムダ株式会社 | インダクタンス素子 |
CN101572170A (zh) * | 2009-03-13 | 2009-11-04 | 中国西电电气股份有限公司 | 一种带辅助铁轭的变压器铁心结构 |
BR112017020753B1 (pt) * | 2015-04-20 | 2021-08-10 | Nippon Steel Corporation | Chapa de aço elétrica com grão orientado |
CN106257462A (zh) * | 2015-06-16 | 2016-12-28 | 长沙理工大学 | 一种变压器内部短路故障非线性仿真的方法 |
CN106501619B (zh) * | 2016-12-31 | 2018-04-13 | 西南交通大学 | 一种铁心片间短路故障下的铁心损耗试验系统及方法 |
EP3726543A4 (en) * | 2018-01-31 | 2021-03-03 | JFE Steel Corporation | ALIGNED ELECTRIC STEEL SHEET, WRAPPED TRANSFORMER CORE WITH IT AND PROCESS FOR MANUFACTURING A WRAPPED CORE |
CN108647382B (zh) * | 2018-04-02 | 2019-05-03 | 西南交通大学 | 一种卷铁心牵引变压器涡流损耗的测评方法 |
CN110008645A (zh) * | 2019-04-30 | 2019-07-12 | 广东工业大学 | 一种变压器损耗计算方法 |
CN110110491A (zh) * | 2019-05-30 | 2019-08-09 | 哈尔滨理工大学 | 一种铁心片间短路时垂直叠片方向的等效电导率计算法 |
CN110399695B (zh) * | 2019-07-31 | 2020-08-18 | 西南交通大学 | 一种考虑磁通密度不均匀分布的卷铁心涡流损耗评估方法 |
CN110399693B (zh) * | 2019-07-31 | 2022-06-07 | 西南交通大学 | 一种片间短路下变压器卷铁心涡流损耗的计算方法 |
CN110571024A (zh) * | 2019-09-23 | 2019-12-13 | 北京中热信息科技有限公司 | 一种缠绕式壳式变压器及其制造方法 |
CN110988529A (zh) * | 2019-11-26 | 2020-04-10 | 国网山西省电力公司电力科学研究院 | 一种变压器热点温度多物理场计算中热源分布确定方法 |
CN111104743A (zh) * | 2019-12-18 | 2020-05-05 | 齐鲁工业大学 | 一种确定变压器直流偏磁瞬态磁场和涡流损耗分布的方法 |
-
2020
- 2020-08-28 CN CN202010881786.8A patent/CN111931310B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111931310A (zh) | 2020-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111931310B (zh) | 一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 | |
CN110399695B (zh) | 一种考虑磁通密度不均匀分布的卷铁心涡流损耗评估方法 | |
CN108647382B (zh) | 一种卷铁心牵引变压器涡流损耗的测评方法 | |
CN110399693B (zh) | 一种片间短路下变压器卷铁心涡流损耗的计算方法 | |
CN107294081B (zh) | 线损率影响因素的相关性确定方法 | |
CN114242425B (zh) | 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法 | |
US2489977A (en) | Laminated core | |
Qin et al. | Multi-layer non-uniform series self-resonant coil for wireless power transfer | |
Marques et al. | Optimization of transmitter magnetic structures for roadway applications | |
CN110263500B (zh) | 一种片间短路下卷铁心多点接地故障的建模方法 | |
CN105512507B (zh) | 控制变压器降容保有功方法 | |
Xiangxiu et al. | Magnetizer optimization of dd type coils for ev wireless charging system | |
CN110275056B (zh) | 一种考虑硅钢片截面宽度分级的卷铁心多点接地故障电流计算方法 | |
Lope et al. | Minimization of vias in PCB implementations of planar coils with litz-wire structure | |
CN114168887A (zh) | 一种计及铁心磁路分级的牵引变压器涡流损耗求解方法 | |
Moradnouri et al. | Amorphous metal triangular cores to improve distribution transformers design | |
CN114400136A (zh) | 一种平面变压器优化设计方法及系统 | |
CN108039271B (zh) | 一种优化ei变压器材料的方法 | |
Elhaminia et al. | Magnetic Flux Density Determination in 3D Wound Core Transformer Using H-balance Equation | |
CN112379223B (zh) | 一种配电变压器绕组匝间绝缘劣化的等效电路模型 | |
CN116127240B (zh) | 一种牵引变压器卷铁心过负荷能力的评估方法 | |
Gukov et al. | Influence of Scattering Flows on the Distribution of Magnetic Induction in the Magnetic Core of a Power Transformer under Load | |
CN106298215B (zh) | 一种配电变压器设计方法 | |
Wang et al. | Research on radial stability of autotransformer inner winding under the short-circuit condition | |
Hu et al. | A reluctance-based model for the design of inductive power transfer coils |
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 |