CN107688679B - 一种硅钢片铁芯宽频电路模型的建模方法 - Google Patents
一种硅钢片铁芯宽频电路模型的建模方法 Download PDFInfo
- Publication number
- CN107688679B CN107688679B CN201610632962.8A CN201610632962A CN107688679B CN 107688679 B CN107688679 B CN 107688679B CN 201610632962 A CN201610632962 A CN 201610632962A CN 107688679 B CN107688679 B CN 107688679B
- Authority
- CN
- China
- Prior art keywords
- silicon steel
- steel sheet
- iron core
- circuit
- layer
- 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
- 229910000976 Electrical steel Inorganic materials 0.000 title claims abstract description 83
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000005291 magnetic effect Effects 0.000 claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 230000035515 penetration Effects 0.000 claims abstract description 4
- 239000010410 layer Substances 0.000 claims description 55
- 230000035699 permeability Effects 0.000 claims description 18
- 230000005415 magnetization Effects 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000004907 flux Effects 0.000 claims description 4
- 238000004804 winding Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 230000006698 induction Effects 0.000 claims description 3
- 239000002356 single layer Substances 0.000 claims description 3
- 239000002344 surface layer Substances 0.000 claims description 3
- 238000009826 distribution Methods 0.000 abstract description 4
- 230000005672 electromagnetic field Effects 0.000 abstract description 4
- 230000005540 biological transmission Effects 0.000 abstract description 3
- 238000005094 computer simulation Methods 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 description 3
- 239000008358 core component Substances 0.000 description 1
- 230000002500 effect on skin Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000003302 ferromagnetic material Substances 0.000 description 1
- 238000009828 non-uniform distribution Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Soft Magnetic Materials (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
本发明涉及电力电子和电力系统的元件仿真建模领域,具体涉及一种硅钢片铁芯宽频电路模型的建模方法。本方法对每层硅钢片进行合理剖分后,基于等效电路法将铁芯等效为一个梯形电路,此电路可按照一定的规则与外电路相连,从而实现整个系统的分析和计算。本发明根据正常工况下电磁波在铁芯中的透入深度,将铁芯中的每片硅钢片分为3层以确定铁芯电路模型的阶数,计算分层后各层硅钢片的等效涡流电阻Rek和非线性电感Ldk,定义一个无源二端口传递网络,用于将铁芯外电路侧的物理量换算到单硅钢片等效网络侧,最终建立与外电路连接的多铁芯模型。本发明一种硅钢片铁芯宽频电路模型的建模方法采用等效电路法,考虑铁芯内部磁场分布的不均匀性,建立了铁芯的等效电路模型,从而使复杂的电磁场计算转换为简单的电路计算。
Description
技术领域
本发明涉及电力电子和电力系统的元件仿真建模领域,具体涉及一种硅钢片铁芯宽频电路模型的建模方法。
背景技术
随着现代电力电子技术的快速发展,硅钢片铁芯在电气设备(变压器、电抗器、电感器等)中得到了广泛应用。进行系统分析时,需要对此类电气设备精确建模,铁芯作为其中电磁转换的核心部件,其宽频模型的准确与否尤为重要。然而,由于铁磁材料的非线性引起铁芯饱和,高频下的集肤效应引起铁芯内磁通分布不均等问题,导致硅钢片铁芯宽频建模变得困难、复杂。
目前,研究硅钢片铁芯动态模型的方法主要有传统解析法、有限元法和等效电路法等。传统解析法需要求解麦克斯韦方程组,计算过程繁琐,目前主要应用于理论分析。有限元法将计算空间剖分为有限个单元,使对电磁场微分方程的求解转化为代数方程的求解,必须对铁芯中每片硅钢片进行网格剖分,使得单元格数量巨大,仿真时间过长。
发明内容
鉴于铁芯内每片硅钢片的磁场分布规律基本相同,对每层硅钢片进行合理剖分后,基于等效电路法可将铁芯等效为一个梯形电路,此电路可按照一定的规则与外电路相连,从而实现整个系统的分析和计算。
基于上述原理,本发明提出如下一种硅钢片铁芯宽频电路模型的建模方法,其特征在于:对每层硅钢片进行合理剖分后,由等效电路法确定此电路模型为n阶等效梯形电路,其中n为自然数,由划分每层硅钢片的层数来确定;每阶电路由等效涡流电阻Re和非线性电感Ld组成,此电路可按照一定的规则与外电路相连,从而实现整个系统的分析和计算,具体过程包括以下步骤:
步骤一:计算正常工况下电磁波在铁芯中的透入深度,由下面的公式计算:
式中,δ为透入深度,σ为铁芯硅钢电导率,ω=2πf为角频率,μ为铁芯硅钢磁导率;
步骤二:对铁芯中的每片硅钢片分为n层以确定铁芯电路模型的阶数;
步骤三:计算分层后各层硅钢片的等效涡流电阻,记Re1,Re2,……Ren分别为各层的等效涡流电阻,可由以下公式计算等效涡流电阻值:
Rek=2ρw/(lcdk)
式中:k=1,2,3……n,ρ为硅钢片的电阻率,w为硅钢片截面宽度,lc近似为铁芯平均磁路长度,dk(k=1,2,3……n)为硅钢片划分后各层区域的厚度;
步骤四:计算分层后各层硅钢片的非线性电感Ldk;
步骤五:建立与外电路连接的多铁芯模型:设电磁装置具有Nc个铁芯,每个铁芯具有NL层硅钢片,铁芯绕组匝数为N,令Nm=N·NL·Nc,
定义一个无源二端口传递网络,用于将铁芯外电路侧的物理量换算到单硅钢片等效网络侧,令铁芯外电路侧为u1,i1,单硅钢片等效电路侧为u2,i2,
则该传递网络的约束方程为:
应用该约束方程,可将铁芯外电路侧的电路元件进行换算,换算后即可建立与外电路相连接的铁芯电路模型。
由等效电路法确定此电路模型为n阶等效梯形电路,每阶电路由等效涡流电阻和非线性电感组成,其元件数值分别由步骤三和步骤四得到,而梯形电路的阶数则由划分硅钢片的层数来确定,划分的层数越多,则模型的精确度越高。
进一步的,对于工作在50Hz--40kHz频段中的铁芯可按以下原则将步骤二中的硅钢片划分3层:由硅钢片表层开始,第1层剖分厚度d1=0.5δ,第2层剖分厚度d2=1.5δ,剩余至中心部分为第3层,剖分厚度d3=0.5d-2δ,其中d为单层硅钢片的厚度。需要明确指出,该原则只是给出了在工作频率下硅钢片剖分层数的最低要求,若对仿真数值精度要求较高,则应相应地增加剖分层数。
进一步的,所述步骤四中计算分层后各层硅钢片的非线性电感Ldk包括以下两个步骤:
首先,进行基于改进郎之万函数的磁化曲线拟合:根据Jiles-Atherton磁滞模型,有:
式中M为整体磁化强度,Ms为饱和磁化强度,H为整体磁场强度,a为函数中的待定系数;
则磁感应强度B可以表达为:
式中μ0为真空磁导率
根据磁导率的定义μ=B/H,由上式得磁导率表达式
继而得磁导率对磁场强度的导数:
其次,计算非线性电感值:根据安培环路定律,有:
Hck·lc+Ha·δa=ik
经推导得:
式中Ack为硅钢片划分各部分的截面积
则非线性电感可表示为:
本发明一种硅钢片铁芯宽频电路模型的建模方法采用等效电路法,考虑铁芯内部磁场分布的不均匀性,建立了铁芯的等效电路模型,从而使复杂的电磁场计算转换为简单的电路计算。
附图说明
图1为单一硅钢片截面划分示意图。
图2为与外电路连接的传递网络。
图3为与外电路连接的多铁芯等效电路。
具体实施方式
下面结合附图对本发明作进一步详细说明。
一种硅钢片铁芯宽频电路模型的建模方法,对每层硅钢片进行合理剖分后,由等效电路法确定此电路模型为n阶等效梯形电路,其中n为自然数,由划分每层硅钢片的层数来确定;每阶电路由等效涡流电阻Re和非线性电感Ld组成,此电路可按照一定的规则与外电路相连,此电路可按照一定的规则与外电路相连,从而实现整个系统的分析和计算,具体过程包括以下步骤:
步骤一:计算正常工况下电磁波在铁芯中的透入深度,由下面的公式计算:
式中,δ为透入深度,σ为铁芯硅钢电导率,ω=2πf为角频率,μ为铁芯硅钢磁导率;
步骤二:对铁芯中的每片硅钢片分为n层以确定铁芯电路模型的阶数。如图1所示,本实施例对于工作在50Hz--40kHz频段中的铁芯可按以下原则将硅钢片划分3层:由硅钢片表层开始,第1层剖分厚度d1=0.5δ,第2层剖分厚度d2=1.5δ,剩余至中心部分为第3层,剖分厚度d3=0.5d-2δ,其中d为单层硅钢片的厚度。需要明确指出,该原则只是给出了在工作频率下硅钢片剖分层数的最低要求,若对仿真数值精度要求较高,则应相应地增加剖分层数。
步骤三:计算分层后各层硅钢片的等效涡流电阻,记Re1,Re2,……Ren分别为各层的等效涡流电阻,可由以下公式计算等效涡流电阻值:
Rek=2ρw/(lcdk)
式中:k=1,2,3……n,ρ为硅钢片的电阻率,w为硅钢片截面宽度,lc近似为铁芯平均磁路长度,dk(k=1,2,3……n)为硅钢片划分后各层区域的厚度;
步骤四:计算分层后各层硅钢片的非线性电感Ldk。该计算过程具体包括以下两个步骤:
首先,进行基于改进郎之万函数的磁化曲线拟合:根据Jiles-Atherton磁滞模型,有:
式中M为整体磁化强度,Ms为饱和磁化强度,H为整体磁场强度,a为函数中的待定系数;
则磁感应强度B可以表达为:
式中μ0为真空磁导率
根据磁导率的定义μ=B/H,由上式得磁导率表达式
继而得磁导率对磁场强度的导数:
其次,计算非线性电感值:根据安培环路定律,有:
Hck·lc+Ha·δa=ik
经推导得:
式中Ack为硅钢片划分各部分的截面积
则非线性电感可表示为:
步骤五:建立与外电路连接的多铁芯模型:设电磁装置具有Nc个铁芯,每个铁芯具有NL层硅钢片,铁芯绕组匝数为N,令Nm=N·NL·Nc,
定义一个无源二端口传递网络,用于将铁芯外电路侧的物理量换算到单硅钢片等效网络侧,令铁芯外电路侧为u1,i1,单硅钢片等效电路侧为u2,i2,如图2所示,
则该传递网络的约束方程为:
应用该约束方程,可将铁芯外电路侧的电路元件进行换算,换算后即可建立与外电路相连接的铁芯电路模型。
例如:电压源us通过串联电阻R0直接作用在具有多个铁芯的电磁装置上,绕组电流为i,采用上述换算公式,即可得与外电路相连的多铁芯等效电路,如图3所示。图中us/Nm、R0/(NmN)、Ni为换算到单硅钢片等效电路侧的参数,假设每层硅钢片划分为三层,计算得Ld1、Ld2、Ld3为各层硅钢片的非线性电感,Re1、Re2为各层硅钢片的等效涡流电阻,则电路模型即为由Ld1、Ld2、Ld3和Re1、Re2组成的三阶等效梯形电路,其中ie1,ie2为各层内产生的涡流。
本发明一种硅钢片铁芯宽频电路模型的建模方法采用等效电路法,考虑铁芯内部磁场分布的不均匀性,建立了铁芯的等效电路模型,从而使复杂的电磁场计算转换为简单的电路计算。
Claims (3)
1.一种硅钢片铁芯宽频电路模型的建模方法,其特征在于:对每层硅钢片进行合理剖分后,由等效电路法确定此电路模型为n阶等效梯形电路,其中n为自然数,由划分每层硅钢片的层数来确定;每阶电路由等效涡流电阻Re和非线性电感Ld组成,此电路可按照一定的规则与外电路相连,从而实现整个系统的分析和计算,具体过程包括以下步骤:
步骤一:计算正常工况下电磁波在铁芯中的透入深度,由下面的公式计算:
式中,δ为透入深度,σ为铁芯硅钢电导率,ω=2πf为角频率,μ为铁芯硅钢磁导率;
步骤二:对铁芯中的每片硅钢片分为n层以确定铁芯电路模型的阶数;
步骤三:计算分层后各层硅钢片的等效涡流电阻,记Re1,Re2,……Ren分别为各层的等效涡流电阻,可由以下公式计算等效涡流电阻值:
Rek=2ρw/(lcdk)
式中:k=1,2,3……n,ρ为硅钢片的电阻率,w为硅钢片截面宽度,lc近似为铁芯平均磁路长度,dk为硅钢片划分后各层区域的厚度,其中k=1,2,3……n;
步骤四:计算分层后各层硅钢片的非线性电感Ldk;
步骤五:建立与外电路连接的多铁芯模型:设电磁装置具有Nc个铁芯,每个铁芯具有NL层硅钢片,铁芯绕组匝数为N,令Nm=N·NL·Nc,
定义一个无源二端口传递网络,用于将铁芯外电路侧的物理量换算到单硅钢片等效网络侧,令铁芯外电路侧为u1,i1,单硅钢片等效电路侧为u2,i2,
则该传递网络的约束方程为:
应用该约束方程,可将铁芯外电路侧的电路元件进行换算,换算后即可建立与外电路相连接的铁芯电路模型。
2.根据权利要求1所述的一种硅钢片铁芯宽频电路模型的建模方法,其特征在于:将步骤二中的每片硅钢片分为3层:由硅钢片表层开始,第1层剖分厚度d1=0.5δ,第2层剖分厚度d2=1.5δ,剩余至中心部分为第3层,剖分厚度d3=0.5d-2δ,其中d为单层硅钢片的厚度。
3.根据权利要求1或2所述的一种硅钢片铁芯宽频电路模型的建模方法,其特征在于:所述步骤四中计算分层后各层硅钢片的非线性电感Ldk包括以下两个步骤:
首先,进行基于改进郎之万函数的磁化曲线拟合:根据Jiles-Atherton磁滞模型,有:
式中M为整体磁化强度,Ms为饱和磁化强度,H为整体磁场强度,a为函数中的待定系数;
则磁感应强度B可以表达为:
式中μ0为真空磁导率
根据磁导率的定义μ=B/H,由上式得磁导率表达式
继而得磁导率对磁场强度的导数:
其次,计算非线性电感值:根据安培环路定律,有:
Hck·lc+Ha·δa=ik
式中,Hck为硅钢片划分各层区域的磁场强度,Ha为气隙中的磁场强度,δa为铁芯的气隙长度,ik为与φk相交链的所有电流之和,φk为各层的磁通量;
经推导得:
式中Ack为硅钢片划分各部分的截面积
则非线性电感可表示为:
电感Ldk是关于dφk/dHck的函数,dφk/dHck的表达式:
将上一步中计算出的磁导率μ(H)和磁导率导数dμ(H)/dH代上式可得dφk/dHck,即可求得非线性电感Ldk的数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610632962.8A CN107688679B (zh) | 2016-08-04 | 2016-08-04 | 一种硅钢片铁芯宽频电路模型的建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610632962.8A CN107688679B (zh) | 2016-08-04 | 2016-08-04 | 一种硅钢片铁芯宽频电路模型的建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107688679A CN107688679A (zh) | 2018-02-13 |
CN107688679B true CN107688679B (zh) | 2020-09-15 |
Family
ID=61151052
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610632962.8A Active CN107688679B (zh) | 2016-08-04 | 2016-08-04 | 一种硅钢片铁芯宽频电路模型的建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107688679B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112084679B (zh) * | 2020-10-14 | 2022-06-17 | 哈尔滨工程大学 | 一种考虑涡流分布不均匀性电磁执行器等效磁路建模方法 |
CN112420227B (zh) * | 2020-10-30 | 2022-08-23 | 中广核检测技术有限公司 | 一种核电站涡流检测探头的制作方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102983576A (zh) * | 2011-09-07 | 2013-03-20 | 上海市电力公司 | 一种用于无功补偿系统的功率模块的单相等效电路 |
CN105302975A (zh) * | 2015-11-09 | 2016-02-03 | 国网福建省电力有限公司 | 一种电磁式电流互感器谐波传变建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1851556B1 (de) * | 2005-02-22 | 2009-07-01 | Trench Switzerland AG | Ringkern-stromwandler mit phasenkompensationsschaltung |
-
2016
- 2016-08-04 CN CN201610632962.8A patent/CN107688679B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102983576A (zh) * | 2011-09-07 | 2013-03-20 | 上海市电力公司 | 一种用于无功补偿系统的功率模块的单相等效电路 |
CN105302975A (zh) * | 2015-11-09 | 2016-02-03 | 国网福建省电力有限公司 | 一种电磁式电流互感器谐波传变建模方法 |
Non-Patent Citations (2)
Title |
---|
A Generalized Dynamic Circuit Model of Magnetic Cores for Low- and High-Frequency Applications-Part 11: Circuit Model Formulation and Implementation;S. Y. R. Hui等;《IEEE TRANSACTIONS ON POWER ELECTRONICS》;19960331;全文 * |
基于Brune综合法的HVDC换流站阀系统宽频建模;刘磊等;《南方电网技术》;20091231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107688679A (zh) | 2018-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lope et al. | AC power losses model for planar windings with rectangular cross-sectional conductors | |
Sato et al. | Model reduction of three-dimensional eddy current problems based on the method of snapshots | |
Salas et al. | Simulation of the saturation and air-gap effects in a pot ferrite core with a 2-D finite element model | |
CN105302975B (zh) | 一种电磁式电流互感器谐波传变建模方法 | |
De Gersem et al. | Three-dimensional–two-dimensional coupled model for eddy currents in laminated iron cores | |
Preis et al. | FEM analysis of eddy current losses in nonlinear laminated iron cores | |
CN107688679B (zh) | 一种硅钢片铁芯宽频电路模型的建模方法 | |
CN112417727B (zh) | 一种考虑端部效应的高频变压器漏电感参数计算方法 | |
Sato et al. | Loss computation of soft magnetic composite inductors based on interpolated scalar magnetic property | |
CN107977485A (zh) | 一种考虑高频涡流效应的绕组损耗半解析计算方法 | |
Jazebi et al. | Duality-synthesized circuit for eddy current effects in transformer windings | |
Bastos et al. | A differential permeability 3-D formulation for anisotropic vector hysteresis analysis | |
Djuric et al. | A compact planar transformer with an improved winding configuration | |
Kovačić et al. | Small-signal calculation of common-mode choke characteristics using finite-element method | |
Hasan et al. | Stabilized reduced-order model of a non-linear eddy current problem by a gappy-POD approach | |
Le-Duc et al. | General integral formulation for the 3D thin shell modeling | |
Salas et al. | Nonlinear saturation modeling of magnetic components with an RM-type core | |
Haerinia et al. | Investigation of receiving pot core effect on magnetic flux density in inductive coupling-based wireless power transfer | |
Ammouri et al. | Design and modeling of planar magnetic inductors for power converters applications | |
Wang | Modeling and design of planar integrated magnetic components | |
Penin et al. | Impact of mechanical deformations of transformer corners on core losses | |
Klis et al. | Fast simulation of wireless power transfer systems with varying coil alignment | |
CN106066846B (zh) | 一种饱和条件下无取向硅钢片铁心磁导率张量表示方法 | |
Wojda | Thermal analytical winding size optimization for different conductor shapes | |
Alvarez et al. | Analysis of impedance matrix in transformer windings through the Finite Element Method (FEM) |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20231117 Address after: Room 203, Building 4, Courtyard 5, Longxiang Plate Making Group, Huilongguan, Changping District, Beijing, 100085 Patentee after: BEIJING HUADIAN ZHICHENG ELECTRICAL EQUIPMENT Co.,Ltd. Address before: 102206 Dewai zhuxinzhuang, Changping District, Beijing Patentee before: NORTH CHINA ELECTRIC POWER University |
|
TR01 | Transfer of patent right |