CN111931310B - 一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 - Google Patents

一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 Download PDF

Info

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
Application number
CN202010881786.8A
Other languages
English (en)
Other versions
CN111931310A (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 CN202010881786.8A priority Critical patent/CN111931310B/zh
Publication of CN111931310A publication Critical patent/CN111931310A/zh
Application granted granted Critical
Publication of CN111931310B publication Critical patent/CN111931310B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR 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
    • G06F17/13Differential equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation 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)。对卷绕的首级而言,以硅钢带厚度方向的中线为卷绕路径的长度标准,其圆角半径会在骨架的尺寸基础上增加硅钢片厚度的一半:
Figure GDA0003106104000000011
后面的层级将有所区别,卷绕堆砌后的半径数值将比上一级增加硅钢片厚度的整长,下面给出第二级、第三级、第四级圆角半径的计算公式:
Figure GDA0003106104000000012
由此可以给出不同层级的圆角半径数值的递推关系,即得到卷铁心第i级的圆角半径Ri的表达式:
Figure GDA0003106104000000013
进一步能够得到卷铁心第i级的整个卷绕路径长度li的表达式:
Figure GDA0003106104000000021
式中,i∈{1,2,3,…,m,m+1,m+2,…,2m},2m为卷铁心的总级数;
由于变压器铁心都是正弦型交流励磁,因此应对磁准静态场的麦克斯韦方程进行求解,同时考虑硅钢片的晶粒取向特性和卷铁心各层级的宽度远大于厚度的特征,得到卷铁心层间短路区域内的涡流场复数形式二阶微分方程:
Figure GDA0003106104000000022
式中,j为虚数单位,Hz(x)代表短路区域内部各点的磁场强度,
Figure GDA0003106104000000023
其中:μ为卷铁心材质的磁导率,σ为卷铁心材质的电导率,ω为角频率,它满足:ω=2πf,f为励磁频率;
根据二阶线性齐次常微分方程的基础算法,得到上述方程的通解:
Figure GDA0003106104000000024
式中,C1和C2为待定常数,通过涡流场的边界条件确定。
2)由于卷铁心所采用的硅钢片是冷轧取向型,其在卷绕过程中无论是在心柱、铁轭还是拐角,都与导磁性能最佳的方向保持一致。相对地,在垂直于叠片的方向,硅钢片的磁性能最差,再加之有绝缘涂层的存在,非短路区域的磁力线几乎不会发生跃迁。结合这种“层间并联,层内串联”的分级磁路模型,可以认为各层级的边界磁场强度与整个卷铁心的平均磁场强度的比值等于该层级的磁路长度与整个卷铁心的平均磁路长度的比值,因此可以计算出卷铁心第i级的边界磁场强度:
Figure GDA0003106104000000025
式中,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的计算式:
Figure GDA0003106104000000031
3)根据安培环路定律的微分形式,对磁场强度函数进行求导并取负,得到卷铁心层间短路区域内涡流密度函数Jy(x)的表达式:
Figure GDA0003106104000000032
4)由于卷铁心发生层间短路时,涡流会扩散到整个故障区域,因此涡流损耗的计算域也需要调整为整个短路区域的总厚度D。根据电磁学对有功功耗的定义,得到卷铁心层间短路区域的涡流损耗计算公式:
Figure GDA0003106104000000033
本发明的有益效果是提出了一种对卷铁心层间短路故障区域涡流损耗进行精准评估的方法,为变压器制造阶段的性能评估和缺陷排查提供了理论基础。
附图说明
图1为本发明中卷铁心整体结构正视图。
图2为本发明中卷铁心局部层间短路示意图。
具体实施方式
下面结合附图对本发明的实施流程作进一步的详述。
图1为本发明中所述的卷铁心整体结构正视图,其卷铁心纵截面形态为圆角矩形,卷绕过程会给定以下参数:心柱长度(L)、铁轭长度(Y)、圆角半径(R),硅钢带厚度(d)。对卷绕的首级而言,以硅钢带厚度方向的中线为卷绕路径的长度标准,其圆角半径会在骨架的尺寸基础上增加硅钢片厚度的一半:
Figure GDA0003106104000000041
后面的层级将有所区别,卷绕堆砌后的半径数值将比上一级增加硅钢片厚度的整长,下面给出第二级、第三级、第四级圆角半径的计算公式:
Figure GDA0003106104000000042
由此可以给出不同层级的圆角半径数值的递推关系,即得到卷铁心第i级的圆角半径Ri的表达式:
Figure GDA0003106104000000043
进一步能够得到卷铁心第i级的整个卷绕路径长度li的表达式:
Figure GDA0003106104000000044
式中,i∈{1,2,3,…,m,m+1,m+2,…,2m},2m为卷铁心的总级数;
由于变压器铁心都是正弦型交流励磁,因此应对磁准静态场的麦克斯韦方程进行求解,同时考虑硅钢片的晶粒取向特性和卷铁心各层级的宽度远大于厚度的特征,得到卷铁心层间短路区域内的涡流场复数形式二阶微分方程:
Figure GDA0003106104000000045
式中,j为虚数单位,Hz(x)代表短路区域内部各点的磁场强度,
Figure GDA0003106104000000046
其中:μ为卷铁心材质的磁导率,σ为卷铁心材质的电导率,ω为角频率,它满足:ω=2πf,f为励磁频率;
根据二阶线性齐次常微分方程的基础算法,得到上述方程的通解:
Figure GDA0003106104000000047
式中,C1和C2为待定常数,通过涡流场的边界条件确定。
由于卷铁心所采用的硅钢片是冷轧取向型,其在卷绕过程中无论是在心柱、铁轭还是拐角,都与导磁性能最佳的方向保持一致。相对地,在垂直于叠片的方向,硅钢片的磁性能最差,再加之有绝缘涂层的存在,非短路区域的磁力线几乎不会发生跃迁。结合这种“层间并联,层内串联”的分级磁路模型,可以认为各层级的边界磁场强度与整个卷铁心的平均磁场强度的比值等于该层级的磁路长度与整个卷铁心的平均磁路长度的比值,因此可以计算出卷铁心第i级的边界磁场强度:
Figure GDA0003106104000000051
式中,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的计算式:
Figure GDA0003106104000000052
根据安培环路定律的微分形式,对磁场强度函数进行求导并取负,得到卷铁心层间短路区域内涡流密度函数Jy(x)的表达式:
Figure GDA0003106104000000053
由于卷铁心发生层间短路时,涡流会扩散到整个故障区域,因此涡流损耗的计算域也需要调整为整个短路区域的总厚度D。根据电磁学对有功功耗的定义,得到卷铁心层间短路区域的涡流损耗计算公式:
Figure GDA0003106104000000061

Claims (1)

1.一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法,其特征在于,铁心材质为高导磁冷轧晶粒取向硅钢片,包括以下步骤:
1)得到卷铁心层间短路区域内的涡流场复数形式微分方程:
Figure FDA0003106103990000011
式中,j为虚数单位,Hz(x)代表短路区域内部各点的磁场强度,中间参数
Figure FDA0003106103990000012
其中:μ为卷铁心材质的磁导率,σ为卷铁心材质的电导率,ω为角频率,它满足:ω=2πf,f为励磁频率;
根据二阶线性齐次常微分方程的基础算法,得到(1)的通解:
Figure FDA0003106103990000013
式中,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的计算式:
Figure FDA0003106103990000014
3)根据安培环路定律的微分形式,以及(2)、(4)对磁场强度的计算结果,得到卷铁心层间短路区域内涡流密度函数Jy(x)的表达式:
Figure FDA0003106103990000015
4)根据电磁学对有功功耗的定义,得到卷铁心层间短路区域的涡流损耗:
Figure FDA0003106103990000021
CN202010881786.8A 2020-08-28 2020-08-28 一种考虑相异磁边值的卷铁心层间短路涡流损耗评估方法 Active CN111931310B (zh)

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)

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

* Cited by examiner, † Cited by third party
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 中国西电电气股份有限公司 一种带辅助铁轭的变压器铁心结构
RU2682364C1 (ru) * 2015-04-20 2019-03-19 Ниппон Стил Энд Сумитомо Метал Корпорейшн Электротехнический стальной лист с ориентированной зеренной структурой
CN106257462A (zh) * 2015-06-16 2016-12-28 长沙理工大学 一种变压器内部短路故障非线性仿真的方法
CN106501619B (zh) * 2016-12-31 2018-04-13 西南交通大学 一种铁心片间短路故障下的铁心损耗试验系统及方法
WO2019151399A1 (ja) * 2018-01-31 2019-08-08 Jfeスチール株式会社 方向性電磁鋼板およびこれを用いてなる変圧器の巻鉄心並びに巻鉄心の製造方法
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 齐鲁工业大学 一种确定变压器直流偏磁瞬态磁场和涡流损耗分布的方法

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) 一种片间短路下变压器卷铁心涡流损耗的计算方法
US11975621B2 (en) Double solenoid EV-DWPT system and parameter optimization method thereof
CN114242425B (zh) 一种计及铁心磁路分级的牵引变压器磁滞损耗求解方法
CN105912774A (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) 一种片间短路下卷铁心多点接地故障的建模方法
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
Kurita et al. Demonstration of loss‐reduction effect of the amorphous cores for larger‐capacity transformers by trial manufacture
CN114168887A (zh) 一种计及铁心磁路分级的牵引变压器涡流损耗求解方法
Moradnouri et al. Amorphous metal triangular cores to improve distribution transformers design
CN114400136A (zh) 一种平面变压器优化设计方法及系统
CN108039271B (zh) 一种优化ei变压器材料的方法
CN116127240B (zh) 一种牵引变压器卷铁心过负荷能力的评估方法
Khakroei et al. Designing of multilayer planar spiral air-core inductor for power electronic applications
Gukov et al. Influence of Scattering Flows on the Distribution of Magnetic Induction in the Magnetic Core of a Power Transformer under Load
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
CN203433981U (zh) 一种用于电气化铁路车辆的电压互感器

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