CN112417727B - 一种考虑端部效应的高频变压器漏电感参数计算方法 - Google Patents

一种考虑端部效应的高频变压器漏电感参数计算方法 Download PDF

Info

Publication number
CN112417727B
CN112417727B CN202011314420.9A CN202011314420A CN112417727B CN 112417727 B CN112417727 B CN 112417727B CN 202011314420 A CN202011314420 A CN 202011314420A CN 112417727 B CN112417727 B CN 112417727B
Authority
CN
China
Prior art keywords
winding
frequency transformer
leakage
energy
thickness
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
CN202011314420.9A
Other languages
English (en)
Other versions
CN112417727A (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.)
China Three Gorges University CTGU
Original Assignee
China Three Gorges University CTGU
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 China Three Gorges University CTGU filed Critical China Three Gorges University CTGU
Priority to CN202011314420.9A priority Critical patent/CN112417727B/zh
Publication of CN112417727A publication Critical patent/CN112417727A/zh
Application granted granted Critical
Publication of CN112417727B publication Critical patent/CN112417727B/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/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Coils Of Transformers For General Uses (AREA)

Abstract

一种考虑端部效应的高频变压器漏电感参数计算方法,包括:确定高频变压器的铁心窗口结构;将高频变压器模型总的漏磁能量划分为铁心内部和铁心外部两个部分;分析单位长度漏磁能量对变压器结构因子的灵敏度,筛选出的决定性影响因子;通过量纲分析方法,对筛选出的决定性影响因子,进行无量纲化处理,得到多个无量纲化参量;建立高频变压器的参数化有限元仿真模型,提取所有不同无量纲变压器组合情况下的漏磁能量;选择基础拟合函数,进行多变量回归分析,确定修正系数。本发明能应用于多层绕组高频变压器漏电感的精确评估,适用于铜箔绕组以及窄箔片、矩形导体在不同填充率和布置方式下的绕组损耗计算,降低了优化设计所需要的计算量和计算时间。

Description

一种考虑端部效应的高频变压器漏电感参数计算方法
技术领域
本发明属于高频变压器设计领域,特别涉及一种考虑端部效应的高频变压器漏电感参数计算方法。
背景技术
为实现大功率双向全桥DC/DC变换器的零电压开关技术(zero voltageswitching,ZVS),降低开关损耗。通常需要利用在双向全桥DC/DC变换器中起电气隔离和电能分配的高频变压器的漏电感作为谐振电路中的电感来实现串联谐振。该电路采用了移相全桥控制方式,实现零电压开关所需的最小电感值取决于变换器的输入输出电压、开关频率、移相角和输出功率。漏电感过大将会降低变换器输出的效率,漏电感过小则不能实现零电压开关。因此,精确计算高频变压器的漏电感对于整个变换器的设计是非常重要的。
在现有技术中的高频变压器漏电感计算中,通常假定磁场是一维变量,忽略了绕组端部效应对应的横向磁场分量的影响。2005年Robert F采用有限元仿真方法分析了绕组与铁轭之间不同垂直绝缘距离下的绕组区域漏磁场分布,指出端部效应会对绕组区域的磁场和电流密度分布产生影响。目前,国内外学者已经相继针对变压器的漏电感参数提出了多种计算方法,如:
(1)、文献一:Dowell P L.Effects of eddy currents in transformerwindings[J].Proceedings of the Institution of Electrical Engineers,1966,113(8):1387-1394.提出了箔片绕组的磁场能量存储在绕组导体中,将绕组阻抗虚部对应的感抗作为绕组漏电感参数。该方法没有考虑层间绝缘区域的漏磁能量。
(2)、文献二:Dauhajre A.Modeling and estimation of leakage phenomena inmagnetic circuits[D].Pasadena,California:California Institute of Technology,1986.假设铁心窗口内的漏磁场呈一维线性分布,进而推导出了漏电感参数的解析公式,但该公式未考虑导体区域的高频涡流效应。
(3)、文献三:Hurley W G,Wilcox D J.Calculation of leakage inductance intransformer windings[J].IEEE Transactions on Power Electronics,1994,9(1):121-126.基于绕组互感和自感的表达式,推导出了考虑漏电感参数频变特性的解析计算式,但该方法只适用于环形铁心。
(4)、文献四:Wilson P R,Wilcock R.Frequency dependent model of leakageinductance for magnetic components[J].Advanced Electromagnetics,2012,1(3):99-106.提出一种高频变压器漏电感的集中参数等效模型,通过电路仿真得到漏电感参数,该方法未能给出漏电感的解析表达式,只适用于电路仿真。
(5)、文献五:律方成,郭云翔,付超,等.基于磁链分区的大功率中频变压器漏电感参数计算方法[J].电工技术学报,2016,31(05):164-169.提出了一种基于磁链分区的大功率中频变压器漏电感参数的解析计算方法,但该模型仅针对宽箔片绕组。
(6)、文献六:陈彬,李琳,赵志斌.一种考虑频变特性的大容量高频变压器漏电感解析计算方法[J].中国电机工程学报,2017,37(13):3928-3937.推导出了一种计及漏电感频变特性的解析计算公式,并考虑了铁心拐角处绕组圆弧区域对漏磁能量的影响,但并未考虑绕组的端部效应。
(7)、文献七:陈彬,李琳,刘海军,等.基于有限元法的高频变压器漏电感和绕组损耗计算与分析[J].电工电能新技术,2018,37(01):8-14.采用有限元方法计算了高频变压器的漏电感参数,并分析了绕组交叉换位对漏电感参数的影响。有限元计算方法精度虽然较高,但该方法只适用于在电磁场仿真软件中进行建模计算,不便于高频变压器的优化设计。
发明内容
针对高频电力变压器漏电感的精确计算,本发明提供一种考虑端部效应的高频变压器漏电感参数计算方法,综合考虑高频涡流效应和端部效应的影响,给出箔片绕组、扁铜线绕组、方形绕组的高频损耗半经验计算公式,节省了漏电感计算时间和计算量,提高了计算精度。
本发明采取的技术方案为:
一种考虑端部效应的高频变压器漏电感参数计算方法,包括以下步骤:
步骤1:确定高频变压器的铁心窗口结构:
铁心窗口结构涉及的参数包括:箔片厚度d、绕组层间绝缘厚度dins、绕组层数m、副边绕组与铁心之间的水平间距dch、绕组与铁心之间的垂直间距dcv、隔离间距diso、铁心窗口高度hw
步骤2:将高频变压器模型总的漏磁能量划分为铁心内部和铁心外部两个部分,而每个部分又划分为绕组区域和隔离区域;借助控制变量法,分析单位长度漏磁能量Wm对变压器结构因子的灵敏度,筛选出决定性影响因子;
步骤3:通过量纲分析方法,对步骤2筛选出的决定性影响因子,进行无量纲化处理,得到多个无量纲化参量,
步骤4:建立高频变压器的参数化有限元仿真模型,提取所有不同无量纲变压器组合情况下的漏磁能量;
步骤5:选择合适形式的基础拟合函数,并进行多变量回归分析,确定修正系数。
所述步骤2中,灵敏度计算表达式如下:
Figure BDA0002790873640000031
式中:
Figure BDA0002790873640000032
为单位长度漏磁能量Wm对变压器结构因子x的变化率;
N为单个结构因子的等间隔点数目;
xi为结构因子x在第i个间隔点处取值;
xi+1为结构因子x在第i+1个间隔点处取值;
WmFEM(xi)为对应于xi的漏磁能量仿真值。
WmFEM(xi+1)为结构因子x在第i+1个间隔点处时的单位长度漏磁能量;
所述步骤3中,确定5个无量纲化参量,如下:
Figure BDA0002790873640000033
X1为与箔片厚度相关的无量纲化参量;△为归一化厚度;d为箔片的厚度;δ为集肤深度。
Figure BDA0002790873640000034
X2为与垂直绝缘间距相关的无量纲化参量;hw为铁心窗口高度;dcv为绕组与铁心之间的垂直间距。
Figure BDA0002790873640000041
X3为与隔离间距相关的无量纲化参量;hw为铁心窗口高度;diso为隔离间距。
X4=m,
X4为与层数相关的无量纲化参量;m为绕组层数。
Figure BDA0002790873640000042
X5为与绕组层间绝缘厚度相关的无量纲化参量;dins为绕组层间绝缘厚度;hw为铁心窗口高度。
所述步骤4中,采用ANSYS/Maxwell电磁场仿真软件高频变压器进行参数化建模,参数化有限元仿真模型的铁心窗口高度设置为hw,在涡流场求解器中选择电流源激励,在短路试验条件下,对箔片厚度d、绕组层间绝缘厚度dins、绕组层数m、原副边绕组间的绝缘间距diso、绕组与铁心之间的垂直间距dcv做参数扫描,在正弦电流激励且原副边绕组安匝数相等的条件下,计算不同设置值时,各个区域单位长度的漏磁能量。
所述步骤5中,基础拟合函数如下:
1)、绕组区域的基础拟合函数:
Figure BDA0002790873640000043
式中:
Figure BDA0002790873640000044
Wr *为绕组区域单位长度漏磁能量,为α、τ、β、ξ的非线性函数;μ为磁导率;α为以无量纲化参量X1-X5为自变量,最高阶次等于2的多项式函数;τ、β、ξ以无量纲化参量X1-X5为自变量,最高阶次等于1的多项式函数;令X0=1,α、τ、β、ξ表达式如下:
Figure BDA0002790873640000051
2)、原副边绕组间隔离区域的基础拟合函数:
Figure BDA0002790873640000052
Figure BDA0002790873640000053
为隔离区域单位长度漏磁能量,为
Figure BDA0002790873640000054
和γ的非线性函数;
Figure BDA0002790873640000055
和γ为以无量纲化参量X1-X5为自变量,最高阶次等于2的多项式函数;令X0=1,
Figure BDA0002790873640000056
和γ表达式如下:
Figure BDA0002790873640000057
当层数m=1时,K=4,待定系数为60个;当层数m≥2时,K=5,待定系数为81个。
该方法还包括步骤6:将步骤5中绕组区域单位长度漏磁能量表达式Wr *=f1(α,τ,β,ξ)、隔离区域单位长度漏磁能量表达式
Figure BDA0002790873640000058
的应用范围扩展至其它形状导体的绕组,并分析匝间距离对半解析公式计算精度的影响。
将无量纲化参量X1进行修正,能够将半经验公式推广至矩形导体绕组:
Figure BDA0002790873640000059
式中:d为矩形导体厚度;w为矩形导体宽度;v为同一层内相邻矩形导体的匝间距;△′为矩形导体绕组的归一化厚度;δ为集肤深度。
本发明一种考虑端部效应的高频变压器漏电感参数计算方法,技术效果如下:
本发明计算方法能应用于多层绕组高频变压器漏电感的精确评估,适用于铜箔绕组以及窄箔片、矩形导体在不同填充率和布置方式下的绕组损耗计算,降低了优化设计所需要的计算量和计算时间,方便快捷,有利于工程应用。
附图说明
图1为本发明计算方法建立的流程图。
图2(a)为铁心窗口内部几何结构因子定义图;
图2(b)为铁心窗口内部几何结构因子定义图。
图3(a)为铁心内部绕组区域漏磁能量拟合曲线图;
图3(b)为铁心内部隔离区域漏磁能量拟合曲线图;
图3(c)为铁心外部绕组区域漏磁能量拟合曲线图;
图3(d)为铁心外部隔离区域漏磁能量拟合曲线图。
图4(a)为铁心内部绕组区域不同匝间距情况下漏磁能量相对偏差对比图;
图4(b)为铁心内部隔离区域不同匝间距情况下漏磁能量相对偏差对比图;
图4(c)为铁心外部绕组区域不同匝间距情况下漏磁能量相对偏差对比图;
图4(d)为铁心外部隔离区域不同匝间距情况下漏磁能量相对偏差对比图。
图5为芯式模型二维结构示意图。
图6为壳式模型二维结构示意图。
图7(a)为芯式高频变压器漏磁场强度分布图;
图7(b)为芯式高频变压器漏磁场强度分布图。
具体实施方式
一种考虑端部效应的高频变压器漏电感参数计算方法,包括以下步骤:
步骤1:确定高频变压器的铁心窗口结构:
铁心窗口结构涉及的参数包括:箔片厚度d、绕组层间绝缘厚度dins、绕组层数m、副边绕组与铁心之间的水平间距dch、绕组与铁心之间的垂直间距dcv、隔离间距diso、铁心窗口高度hw
步骤2:由于变压器铁心窗口内外漏磁场的分布存在一些差异,为了更加精确地计算变压器铁心内外漏磁能量、明确变压器各区域漏磁能量对各个结构参数的灵敏度以及端部效应对各区域漏磁能量的影响关系。将高频变压器模型总的漏磁能量划分为铁心内部和铁心外部两个部分,而每个部分又划分为绕组区域和隔离区域;借助控制变量法,分析单位长度漏磁能量Wm对变压器结构因子的灵敏度,筛选出决定性影响因子。
灵敏度计算表达式如下:
Figure BDA0002790873640000071
式中:
Figure BDA0002790873640000072
为单位长度漏磁能量Wm对变压器结构因子x的变化率;
N为单个结构因子的等间隔点数目;
xi为结构因子x在第i个间隔点处取值;
xi+1为结构因子x在第i+1个间隔点处取值;
WmFEM(xi)为对应于xi的漏磁能量仿真值。
WmFEM(xi+1)为结构因子x在第i+1个间隔点处时的单位长度漏磁能量;
筛选出的决定性影响因子,决定性影响因子包括箔片厚度d、绕组层间绝缘厚度dins、绕组层数m、绕组与铁心之间的垂直间距dcv、隔离间距diso、铁心窗口高度hw
步骤3:通过量纲分析方法,对步骤2筛选出的决定性影响因子,进行无量纲化处理,得到多个无量纲化参量。具体是:
确定5个无量纲化参量,如下:
Figure BDA0002790873640000073
X1为与箔片厚度相关的无量纲化参量;△为归一化厚度;d为箔片的厚度;δ为集肤深度。
Figure BDA0002790873640000074
X2为与垂直绝缘间距相关的无量纲化参量;hw为铁心窗口高度;dcv为绕组与铁心之间的垂直间距。
Figure BDA0002790873640000081
X3为与隔离间距相关的无量纲化参量;hw为铁心窗口高度;diso为隔离间距。
X4=m;
X4为与层数相关的无量纲化参量;m为绕组层数。
Figure BDA0002790873640000082
X5为与绕组层间绝缘厚度相关的无量纲化参量;dins为绕组层间绝缘厚度;hw为铁心窗口高度。
为了保证半经验公式的实用性,上述无量纲化参量的取值范围应该能够满足不同应用背景下高压高频变压器的设计要求。在电力电子变换器优化设计过程中,内部磁性元件的箔片绕组厚度近似等于集肤深度时,可降低集肤效应的影响,d/δ近似等于或略小于1。变换器中非正弦负载电流含有较多的奇次谐波分量,3次、5次等低阶次谐波幅值较大,其余各阶次谐波的幅值较小,可以忽略。因此,将参量X1的下限设置为0.5,上限设置为6,可以满足高频变压器绕组的设计要求。根据变压器设计的绝缘要求,将无量纲化参量(hw-2dcv)/hw、diso/hw和dins/hw的取值范围分别设置为0.4~1、0.02~0.14和0.01~0.04,可以满足高频变压器的绝缘设计。
步骤4:建立高频变压器的参数化有限元仿真模型,提取所有不同无量纲变压器组合情况下的漏磁能量。
采用ANSYS/Maxwell电磁场仿真软件高频变压器进行参数化建模,参数化有限元仿真模型的铁心窗口高度设置为hw=100mm,在涡流场求解器中选择电流源激励,电流频率为5kHz,在短路试验条件下,对箔片厚度d、绕组层间绝缘厚度dins、绕组层数m、原副边绕组间的绝缘间距diso、绕组与铁心之间的垂直间距dcv做参数扫描,在正弦电流激励且原副边绕组安匝数相等的条件下,计算不同设置值时,各个区域单位长度的漏磁能量。
步骤5:选择合适形式的基础拟合函数,并进行多变量回归分析,确定修正系数。
基础拟合函数如下:
1)、绕组区域的基础拟合函数:
Figure BDA0002790873640000091
式中:
Figure BDA0002790873640000092
Wr *为绕组区域单位长度漏磁能量,为α、τ、β、ξ的非线性函数;μ为磁导率;α为以无量纲化参量X1-X5为自变量,最高阶次等于2的多项式函数;τ、β、ξ以无量纲化参量X1-X5为自变量,最高阶次等于1的多项式函数;令X0=1,α、τ、β、ξ表达式如下:
Figure BDA0002790873640000093
2)、原副边绕组间隔离区域的基础拟合函数:
Figure BDA0002790873640000094
Figure BDA0002790873640000095
为隔离区域单位长度漏磁能量,为
Figure BDA0002790873640000096
和γ的非线性函数;
Figure BDA0002790873640000097
和γ为以无量纲化参量X1-X5为自变量,最高阶次等于2的多项式函数;令X0=1,
Figure BDA0002790873640000098
和γ表达式如下:
Figure BDA0002790873640000099
当层数m=1时,K=4,待定系数为60个;当层数m≥2时,K=5,待定系数为81个。
步骤6:将步骤5中绕组区域单位长度漏磁能量表达式Wr *=f1(α,τ,β,ξ)、隔离区域单位长度漏磁能量表达式
Figure BDA0002790873640000101
的应用范围扩展至其它形状导体的绕组,如矩形扁铜线、方形等,并分析匝间距离对半解析公式计算精度的影响。
将无量纲化参量X1进行修正,能够将半经验公式推广至矩形导体绕组:
Figure BDA0002790873640000102
式中:d为矩形导体厚度;w为矩形导体宽度;v为同一层内相邻矩形导体的匝间距;△′为矩形导体绕组的归一化厚度;δ为集肤深度。
实施例:
本发明一种考虑端部效应的高频变压器漏电感参数计算方法,应用范围不能受限于特定类型的绕组结构,例如:矩形扁铜带、方形导体,单层或多层导体,不同填充率等。除此之外,输入变量应该为一组无量纲化参量。本发明提出的考虑端部效应影响的漏电感参数半经验公式的建立流程如图1所示。
在如下假设下,建立半经验公式法:
①、原副边绕组由相互平行且等厚度的直箔片组成;
②、原副边绕组的端部距离上下铁轭的距离相等,并缠绕于高磁导率的铁心上;
③、原副边绕组的安匝数相等;
④、副边绕组位于零磁动势和最大磁动势之间;
⑤、原边为单层导体,如矩形实心导线或箔片。
步骤一、决定性影响因子筛选:
半经验公式的待定系数和多变量回归分析的复杂程度完全依赖于自变量的数目。因此需要***分析铁心窗口几何结构对漏磁能量W的影响规律,筛选出决定性影响因子。高频变压器铁心窗口的原副边绕组均由单匝箔片构成。铁心窗口结构可以由以下结构确定:箔片厚度d;绕组层间绝缘厚度dins;绕组层数m;副边绕组与铁心之间的水平间距dch;绕组与铁心之间的垂直间距dcv;隔离间距diso;铁心窗口高度hw
由于变压器铁心窗口内外漏磁场的分布存在一些差异,为了更加精确地计算变压器铁心内外漏磁能量、明确变压器各区域漏磁能量对各个结构参数的灵敏度以及端部效应对各区域漏磁能量的影响关系,将变压器模型总的漏磁能量划分为铁心内部和铁心外部两个部分,而每个部分又可划分为绕组区域和隔离区域,如图2(a)、图2(b)所示。
箔片绕组漏磁能量取决于频率f和变压器的结构参数。借助控制变量法,分析单位长度漏磁能量Wm对变压器结构因子的灵敏度。值得注意的是,灵敏度分析是在f=5kHz和hw=100mm条件下进行,当频率f和铁心窗口高度hw改变时,漏磁能量对结构因子的灵敏度大小会发生改变,但并不会影响决定性影响因子的筛选结果。采用ANSYS/Maxwell电磁场仿真软件建立图2(a)、图2(b)所示高频变压器二维仿真模型,利用有限元方法计算出单一结构因子在各个等间隔点处的漏磁能量,其余结构因子保持不变,计算出
Figure BDA0002790873640000111
的平均值。灵敏度的计算表达式如下:
Figure BDA0002790873640000112
式中:N为单个结构因子的等间隔点数目;xi为结构因子x在第i个间隔点处取值;WmFEM(xi)为对应于xi的漏磁能量仿真值。
dch是影响最小的几何变量,因此该结构因子可以剔除。
步骤二、无量纲化参量及其有效范围:
各个区域的漏磁能量与变压器结构参数密切相关,为了简化方程形式,消除原始变量量纲的影响。借助量纲分析方法,对步骤一筛选出的决定性影响因子进行无量纲化处理,最终确定5个无量纲化参量,如下:
Figure BDA0002790873640000113
Figure BDA0002790873640000114
利用上述5个无量纲化参量,可以唯一确定图2所示铁心窗口结构。
为了保证半经验公式的实用性,上述无量纲化参量的取值范围应该能够满足不同应用背景下高压高频变压器的设计要求。在电力电子变换器优化设计过程中,内部磁性元件的箔片绕组厚度近似等于集肤深度时,可降低集肤效应的影响,d/δ近似等于或略小于1。变换器中非正弦负载电流含有较多的奇次谐波分量,3次、5次等低阶次谐波幅值较大,其余各阶次谐波的幅值较小,可以忽略。因此,将参量X1的下限设置为0.5,上限设置为6,可以满足高频变压器绕组的设计要求。根据变压器设计的绝缘要求,将无量纲化参量(hw-2dcv)/hw、diso/hw和dins/hw的取值范围分别设置为0.4~1、0.02~0.14和0.01~0.04,可以满足高频变压器的绝缘设计。最终,确定的各个无量纲化参量的有效范围如表1所示。
表1无量纲化参量及其有效范围
Figure BDA0002790873640000121
步骤三、参数化建模及有限元仿真:
采用ANSYS/Maxwell电磁场仿真软件对图2(a)、图2(b)所示高频变压器进行参数化建模。二维有限元模型的铁心窗口高度设置为hw=100mm,在涡流场求解器中选择电流源激励,电流频率为5kHz。在短路试验条件下,对箔片厚度d、绕组层间绝缘厚度dins、绕组层数m、原副边绕组间的绝缘间距diso、绕组与铁心之间的垂直间距dcv做参数扫描,在正弦电流激励且原副边绕组安匝数相等的条件下,计算不同设置值时各个区域单位长度的漏磁能量。由于导线区域存在集肤效应,在集肤效应层应进行加密剖分,集肤效应层以下的网格可以相对稀疏,将透入深度的剖分层数设置为6层,其余区域采用自适应剖分。
有限元模型中的参数化变量对应于步骤一中决定性影响因子(d,dcv,diso,m,dins),各个参量的变化范围如表1所示。例如,对于无量纲化参量X1,其取值范围为0.5~6(共12个取值),可以通过改变箔片厚度d来实现,频率f及集肤深度δ保持不变(f=5kHz)。类似地,保持铁心窗口高度不变(hw=100mm),改变dcv、diso、dins实现X2、X3、X5在各自有效范围内的取值。经过上述参数扫描计算,最终将得到20592组不同绕组结构下单位长度的漏磁能量的有限元仿真值。
步骤四、拟合函数及回归分析:
该步骤需要确定一个基础拟合函数,通过多变量回归分析方法使拟合函数逼近仿真值。由于传统漏磁能量解析表达式考虑了集肤效应和邻近效应的影响,方程本身具有一定物理意义,且涉及了绕组层数m、归一化厚度△、dins/hw以及diso/hw四个无量纲化参量。经过多次回归分析,最终确定的基础拟合函数如下:
1)绕组区域的基础拟合函数:
Figure BDA0002790873640000131
式中
Figure BDA0002790873640000132
2)原副边绕组间隔离区域的基础拟合函数:
Figure BDA0002790873640000133
由式(3)、(4)可知,漏磁能量的半经验公式是以α、τ、β、ξ、
Figure BDA0002790873640000137
和γ为自变量的非线性函数,如下:
Figure BDA0002790873640000134
令α、τ、β、ξ、
Figure BDA0002790873640000135
和γ为因变量,无量纲化参量X1-X5为自变量,此时半经验公就可以计及各个决定性影响因子的影响。考虑到计算精度和拟合函数复杂程度,α、
Figure BDA0002790873640000136
和γ选择了最高阶次等于2的多项式函数,τ、β和ξ选择了最高阶次等于1的多项式函数。当m=1时,K=4,待定系数为60个;当m≥2时,K=5,待定系数为81个。令X0=1,α、τ、β、ξ、
Figure BDA0002790873640000141
和γ的表达式如下:
Figure BDA0002790873640000142
在参数扫描结果的基础上,待定系数Pij,Qi,Ji、Ti、Fij和Gij可以通过最小二乘法拟合得到,此时拟合值与有限元仿真值之差的平方和最小。
Figure BDA0002790873640000143
采用拟合软件对基础拟合函数在无量纲化参量有效范围内进行多次拟合,最终其残差平方和(SSE)均小于1,相关系数(R)均大于99.9%,四个区域漏磁能量拟合曲线如图3(a)~图3(d)所示。
步骤五、非连续导体的计算精度:
前面已经推导了单箔片绕组的半经验公式,通过调整每层的填充率可以将半经验公式推广到其它形状导体的绕组,如矩形、方形等。紧密缠绕导体绕组的漏磁能量不同于稀疏缠绕导体绕组的漏磁能量。因此,将半经验公式应用于多孔导体层时,需要考虑匝间距离的影响。将无量纲化参量X1按照式(8)进行修正,可以将半经验公式可以推广至矩形导体绕组。
Figure BDA0002790873640000144
式中d为矩形导体厚度;w为矩形导体宽度;v为同一层内相邻矩形导体的匝间距。
在均匀缠绕的分层绕组中,匝间距保持不变,匝间距v为矩形导体半厚度的5%~20%,等于2倍的导线自绝缘厚度。因此,将2v/d设置为0.05~0.2可以满足设计要求。其余无量纲化参量的表达式(X1~X5)及其取值范围与箔片导体相同。
为了研究匝间距v对半经验公式计算精度的影响,建立了一台高频变压器仿真模型,原边绕组为单层矩形实心导线,副边绕组为四层矩形实心导线。匝间距的的取值范围为0~0.1d。铁心窗口的其余几何结构如下:d/δ=0.5~6;每层匝数Nt=4;dcv=8mm;dins=2mm;diso=10mm;hw=100mm。以仿真结果为参照,图4(a)~图4(d)给出了各个区域半经验公式的相对偏差,结果表明:
1):由图4(a)和图4(c)可知,匝间距的变化会对绕组区域的漏磁能量产生影响,并且当导体匝间距较大,导线排列较为稀疏时,相对偏差越小;
2):由图4(b)和图4(d)可知,铁心内、外部隔离区域的漏磁能量受匝间距变化的影响较小;
3):由图4(a)、图4(b)、图4(c)、图4(d)可知,推广后的半经验公式对于原副边绕组隔离区域的计算精度高于绕组区域的计算精度,这可能是由于引入的填充率不够精确的原因。
针对一台芯式和一台壳式高频变压器模型,分别采用新的半经验方法、传统解析计算方法、有限元仿真方法和实验测量方法提取两台高频变压器的漏电感参数,用于研究新方法的计算精度。两台变压器模型铁心材料为纳米晶合金,原副边绕组由矩形扁铜线绕制而成,容量为10kW,电压等级为0.54kV/0.54kV,工作频率为5kHz。芯式高频变压器的二维结构图及实物图如图5所示,壳式高频变压器二维结构图及实物图如图6所示,主要参数见表2所示。
表2芯式与壳式高频变压器模型的主要参数
Figure BDA0002790873640000161
利用电磁场分析软件对高频变压器进行时谐电磁场仿真计算。与二维有限元相比,三维仿真可以更好地描述变压器的几何形状,因此三维有限元模型在漏电感仿真中具有更高的精度。此外,三维仿真模型可以准确反映绕组弯曲部分和窗口以外区域绕组端部的漏磁场。芯式高频变压器的漏磁场仿真云图如图7(a)所示,壳式高频变压器的漏磁场仿真云图如图7(b)所示。
采用Agilent 4294A高精度阻抗分析仪对高频变压器试验模型的漏电感参数进行测量。测量频率范围为40Hz~100kHz,测量前进行开短路校准,提高测量结果的准确性。测量过程中施加低电平电流以确保变压器工作在线性区。由于将试验模型副边绕组短路,所测量得到的芯式和壳式变压器电感为归算到原边绕组的等效电感。
表3列出了半经验公式、有限元仿真与实验测量方法得到的频率5kHz下芯式和壳式高频变压器模型的漏电感参数。由表3可知,以测量值为参考,半经验公式和有限元仿真值与测量值的偏差较小。
表3芯式与壳式高频变压器模型漏电感参数表
Figure BDA0002790873640000162
Figure BDA0002790873640000171

Claims (5)

1.一种考虑端部效应的高频变压器漏电感参数计算方法,其特征在于包括以下步骤:
步骤1:确定高频变压器的铁心窗口结构:
铁心窗口结构涉及的参数包括:箔片厚度d、绕组层间绝缘厚度dins、绕组层数m、副边绕组与铁心之间的水平间距dch、绕组与铁心之间的垂直间距dcv、隔离间距diso、铁心窗口高度hw
步骤2:将高频变压器模型总的漏磁能量划分为铁心内部和铁心外部两个部分,而每个部分又划分为绕组区域和隔离区域;借助控制变量法,分析单位长度漏磁能量Wm对变压器结构因子的灵敏度,筛选出决定性影响因子;
步骤3:通过量纲分析方法,对步骤2筛选出的决定性影响因子,进行无量纲化处理,得到多个无量纲化参量;
步骤4:建立高频变压器的参数化有限元仿真模型,提取所有不同无量纲变压器组合情况下的漏磁能量;
步骤5:选择合适形式的基础拟合函数,并进行多变量回归分析,确定修正系数;
步骤5中基础拟合函数如下:
1)、绕组区域的基础拟合函数:
Figure FDA0003551985700000011
式中:
Figure FDA0003551985700000012
Wr *为绕组区域单位长度漏磁能量,为α、τ、β、ξ的非线性函数;μ为磁导率;α为以无量纲化参量X1-X5为自变量,最高阶次等于2的多项式函数;τ、β、ξ以无量纲化参量X1-X5为自变量,最高阶次等于1的多项式函数;令X0=1,α、τ、β、ξ表达式如下:
Figure FDA0003551985700000021
2)、原副边绕组间隔离区域的基础拟合函数:
Figure FDA0003551985700000022
Figure FDA0003551985700000023
为隔离区域单位长度漏磁能量,为
Figure FDA0003551985700000024
和γ的非线性函数;
Figure FDA0003551985700000025
和γ为以无量纲化参量X1-X5为自变量,最高阶次等于2的多项式函数;令X0=1,
Figure FDA0003551985700000026
和γ表达式如下:
Figure FDA0003551985700000027
2.根据权利要求1所述一种考虑端部效应的高频变压器漏电感参数计算方法,其特征在于:所述步骤2中,灵敏度计算表达式如下:
Figure FDA0003551985700000028
式中:
Figure FDA0003551985700000029
为单位长度漏磁能量Wm对变压器结构因子x的变化率;
N为单个结构因子的等间隔点数目;
xi为结构因子x在第i个间隔点处取值;
xi+1为结构因子x在第i+1个间隔点处取值;
WmFEM(xi)为对应于xi的漏磁能量仿真值;
WmFEM(xi+1)为结构因子x在第i+1个间隔点处时的单位长度漏磁能量。
3.根据权利要求1所述一种考虑端部效应的高频变压器漏电感参数计算方法,其特征在于:所述步骤3中,确定5个无量纲化参量,如下:
Figure FDA0003551985700000031
X1为与箔片厚度相关的无量纲化参量;Δ为归一化厚度;d为箔片的厚度;δ为集肤深度;
Figure FDA0003551985700000032
X2为与垂直绝缘间距相关的无量纲化参量;hw为铁心窗口高度;dcv为绕组与铁心之间的垂直间距;
Figure FDA0003551985700000033
X3为与隔离间距相关的无量纲化参量;hw为铁心窗口高度;diso为隔离间距;
X4=m,
X4为与层数相关的无量纲化参量;m为绕组层数;
Figure FDA0003551985700000034
X5为与绕组层间绝缘厚度相关的无量纲化参量;dins为绕组层间绝缘厚度;hw为铁心窗口高度。
4.根据权利要求1所述一种考虑端部效应的高频变压器漏电感参数计算方法,其特征在于:所述步骤4中,采用ANSYS/Maxwell电磁场仿真软件高频变压器进行参数化建模,参数化有限元仿真模型的铁心窗口高度设置为hw,在涡流场求解器中选择电流源激励,在短路试验条件下,对箔片厚度d、绕组层间绝缘厚度dins、绕组层数m、原副边绕组间的绝缘间距diso、绕组与铁心之间的垂直间距dcv做参数扫描,在正弦电流激励且原副边绕组安匝数相等的条件下,计算不同设置值时,各个区域单位长度的漏磁能量。
5.根据权利要求1所述一种考虑端部效应的高频变压器漏电感参数计算方法,其特征在于:该方法还包括步骤6:将步骤5中绕组区域单位长度漏磁能量表达式Wr *=f1(α,τ,β,ξ)、隔离区域单位长度漏磁能量表达式
Figure FDA0003551985700000041
的应用范围扩展至其它形状导体的绕组,并分析匝间距离对半解析公式计算精度的影响;
将无量纲化参量X1进行修正,能够将半经验公式推广至矩形导体绕组:
Figure FDA0003551985700000042
式中:d为矩形导体厚度;w为矩形导体宽度;v为同一层内相邻矩形导体的匝间距;Δ′为矩形导体绕组的归一化厚度;δ为集肤深度。
CN202011314420.9A 2020-11-20 2020-11-20 一种考虑端部效应的高频变压器漏电感参数计算方法 Active CN112417727B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011314420.9A CN112417727B (zh) 2020-11-20 2020-11-20 一种考虑端部效应的高频变压器漏电感参数计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011314420.9A CN112417727B (zh) 2020-11-20 2020-11-20 一种考虑端部效应的高频变压器漏电感参数计算方法

Publications (2)

Publication Number Publication Date
CN112417727A CN112417727A (zh) 2021-02-26
CN112417727B true CN112417727B (zh) 2022-05-06

Family

ID=74777214

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011314420.9A Active CN112417727B (zh) 2020-11-20 2020-11-20 一种考虑端部效应的高频变压器漏电感参数计算方法

Country Status (1)

Country Link
CN (1) CN112417727B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283071B (zh) * 2021-05-20 2023-01-24 郑州轻工业大学 一种用于获取高频变压器绕组温升的方法
CN113871173A (zh) * 2021-09-27 2021-12-31 上海交通大学 一种固体绝缘磁集成高频变压器及其屏蔽环设计方法
CN114742003A (zh) * 2022-04-01 2022-07-12 广东风华高新科技股份有限公司 一种电感器的感量测试方法、装置、设备及存储介质
CN116502472B (zh) * 2023-06-25 2023-11-21 深圳市铭昱达电子有限公司 一种高频变压器漏感计算方法及***

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105095609A (zh) * 2015-09-21 2015-11-25 武汉大学 基于有限元法的变压器电磁振动噪声计算方法
CN105183947A (zh) * 2015-08-12 2015-12-23 华北电力大学 基于寄生参数效应分析的高频变压器传输特性优化方法
CN106257462A (zh) * 2015-06-16 2016-12-28 长沙理工大学 一种变压器内部短路故障非线性仿真的方法
CN107273573A (zh) * 2017-05-16 2017-10-20 电子科技大学 一种基于环电流的漏磁仿真方法
CN107368653A (zh) * 2017-07-20 2017-11-21 云南电网有限责任公司电力科学研究院 一种变压器漏磁场的二维有限元仿真方法和***
CN108920751A (zh) * 2018-05-24 2018-11-30 西安交通大学 基于拓扑优化的电力变压器绕组变形状态逆向求解方法
CN109697335A (zh) * 2019-01-17 2019-04-30 广东电网有限责任公司 用于变压器的自漏感计算方法、装置、设备及存储介质
CN110532647A (zh) * 2019-08-09 2019-12-03 三峡大学 考虑端部效应的高频变压器绕组损耗半经验计算方法
CN111597749A (zh) * 2020-05-28 2020-08-28 南通大学 一种高频变压器空心绕组交流电阻系数的计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101961172B1 (ko) * 2017-01-18 2019-03-22 한국전자통신연구원 코일 설계를 위한 변수 제공 방법 및 이를 위한 장치

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106257462A (zh) * 2015-06-16 2016-12-28 长沙理工大学 一种变压器内部短路故障非线性仿真的方法
CN105183947A (zh) * 2015-08-12 2015-12-23 华北电力大学 基于寄生参数效应分析的高频变压器传输特性优化方法
CN105095609A (zh) * 2015-09-21 2015-11-25 武汉大学 基于有限元法的变压器电磁振动噪声计算方法
CN107273573A (zh) * 2017-05-16 2017-10-20 电子科技大学 一种基于环电流的漏磁仿真方法
CN107368653A (zh) * 2017-07-20 2017-11-21 云南电网有限责任公司电力科学研究院 一种变压器漏磁场的二维有限元仿真方法和***
CN108920751A (zh) * 2018-05-24 2018-11-30 西安交通大学 基于拓扑优化的电力变压器绕组变形状态逆向求解方法
CN109697335A (zh) * 2019-01-17 2019-04-30 广东电网有限责任公司 用于变压器的自漏感计算方法、装置、设备及存储介质
CN110532647A (zh) * 2019-08-09 2019-12-03 三峡大学 考虑端部效应的高频变压器绕组损耗半经验计算方法
CN111597749A (zh) * 2020-05-28 2020-08-28 南通大学 一种高频变压器空心绕组交流电阻系数的计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"一种基于ANSYS Maxwell的变压器漏电感的计算方法研究";刘梦 等;《变压器》;20200725;第57卷(第07期);第20-25页 *
"一种考虑频变特性的大容量高频变压器漏电感解析计算方法";陈彬 等;《中国电机工程学报》;20161107;第37卷(第13期);第3928-3937页 *
"基于有限元法的高频变压器漏电感和绕组损耗计算与分析";陈彬 等;《电工电能新技术》;20180123;第37卷(第1期);第8-14页 *

Also Published As

Publication number Publication date
CN112417727A (zh) 2021-02-26

Similar Documents

Publication Publication Date Title
CN112417727B (zh) 一种考虑端部效应的高频变压器漏电感参数计算方法
Nan et al. Simplified high-accuracy calculation of eddy-current loss in round-wire windings
CN110532647B (zh) 考虑端部效应的高频变压器绕组损耗半经验计算方法
Kaymak et al. Comparison of analytical methods for calculating the AC resistance and leakage inductance of medium-frequency transformers
CN113283073B (zh) 一种三相高频大功率变压器的多目标优化设计方法
Guillod et al. Magnetic equivalent circuit of MF transformers: Modeling and parameter uncertainties
Rasekh et al. Artificial neural network aided loss maps for inductors and transformers
Vaisambhayana et al. State of art survey for design of medium frequency high power transformer
CN117057297A (zh) 一种基于混合铜厚绕组结构的反激平面变压器设计方法
Eslamian et al. An accurate analytical method for leakage inductance calculation of shell-type transformers with rectangular windings
Das et al. Calculation of model based capacitances of a two-winding high-frequency transformer to predict its natural resonance frequencies
Ammouri et al. Design and modeling of planar magnetic inductors for power converters applications
Luo et al. 2-D winding losses calculation for round conductor coil
Guan et al. Analysis and design of planar inductor and transformer for resonant converter
Ammouri et al. PCB-planar transformers equivalent circuit model identification using genetic algorithm
Barqi Optimum design approach of high frequency transformer: including the effects of eddy currents
González-Teodoro et al. Determination of wire resistance caused by skin effect using modified 3D finite element model
Daniele et al. FEM analysis of a HF coreless transformer for automotive applications
CN116502472B (zh) 一种高频变压器漏感计算方法及***
Wan et al. An Improved Calculation Method of Leakage Inductance for High-frequency Transformer Based on Precise Analysis of Magnetic Field
Wu et al. Reduction of Copper Loss of Planar Transformers Based on Optimized Width Winding
CN117709133B (zh) 一种多路输出反激平面电感的设计方法
Kanakri et al. Advances In Planar Transformer Circuit Models: Stray Capacitance Equivalent Circuit for Two and Four Winding Transformers With H-Parameter Network Model
CN113589213B (zh) 一种电网特殊负荷条件下互感器计量特性的试验方法
Chung et al. Effect of switching frequency and number of winding layers on copper loss of an inductor

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210226

Assignee: Yichang Yizhixing Technology Co.,Ltd.

Assignor: CHINA THREE GORGES University

Contract record no.: X2023980034895

Denomination of invention: A Calculation Method of Leakage inductance Parameters of High Frequency Transformer Considering End Effect

Granted publication date: 20220506

License type: Common License

Record date: 20230426

EE01 Entry into force of recordation of patent licensing contract