{"id":1518,"date":"2023-05-21T21:32:51","date_gmt":"2023-05-21T18:32:51","guid":{"rendered":"https:\/\/fti.dp.ua\/conf\/?p=1518"},"modified":"2023-06-12T08:47:24","modified_gmt":"2023-06-12T05:47:24","slug":"05216-1930","status":"publish","type":"post","link":"https:\/\/fti.dp.ua\/conf\/2023\/05216-1930\/","title":{"rendered":"On Typical Flaws of the Transformer Models for Inrush Current Evaluation"},"content":{"rendered":"\n<h1 class=\"wp-block-heading citation_title\">On Typical Flaws of the Transformer Models for Inrush Current Evaluation<\/h1>\n\n\n\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<h5 class=\"wp-block-heading citation_author\">Sergey Zirka<\/h5>\n\n\n\n<p class=\"citation_author_url\"><em>ORCID: <a href=\"https:\/\/orcid.org\/0000-0001-7607-1436\" target=\"_blank\" rel=\"noopener\" title=\"\">https:\/\/orcid.org\/0000-0001-7607-1436<\/a><\/em><\/p>\n\n\n\n<p><em>Oles Honchar Dnipro National University, Dnipro<\/em><\/p>\n<\/div><\/div>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<div class=\"wp-block-group\"><div class=\"wp-block-group__inner-container is-layout-constrained wp-block-group-is-layout-constrained\">\n<h5 class=\"wp-block-heading citation_author\">Yuriy Moroz<\/h5>\n\n\n\n<p class=\"citation_author_url\"><em>ORCID:&nbsp;<a href=\"https:\/\/orcid.org\/0000-0001-7696-1324\" target=\"_blank\" rel=\"noopener\" title=\"\">https:\/\/orcid.org\/0000-0001-7696-1324<\/a><\/em><\/p>\n\n\n\n<p><em><em>Oles Honchar Dnipro National University, Dnipro<\/em><\/em><\/p>\n<\/div><\/div>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center\">1 Introduction<\/h3>\n\n\n\n<p>For many decades, the main, if not the only, equivalent circuit of a single-phase transformer was its T-model, Fig. 1. To avoid ambiguity and take into account the topology of the core and windings, we will use the terms <em>inner<\/em> and <em>outer<\/em> windings, instead of the primary and secondary windings. So, inductances <em>L<\/em><sub>S1<\/sub> and <em>L<\/em>\u2032<sub>S2<\/sub> are the so-called &#8220;leakage inductances&#8221; of the inner and outer windings, respectively, <em>r<\/em><sub>1<\/sub> and <em>r<\/em>\u2032<sub>2<\/sub> are their resistances. In order to reproduce hysteretic (quasi-static) properties of the core material and to account for the dynamic core losses and saturation, the magnetizing branch <em>L<\/em><sub>m<\/sub> in Fig. 1 is represented by an ATP-inductor DHM [1], which implements, starting with 2019 version, a Dynamic Hysteresis Model (DHM) [2].<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-1-\u2013-a-conventional-equivalent-circuit-of-a-single-phase-transformer.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-1-\u2013-a-conventional-equivalent-circuit-of-a-single-phase-transformer.png\" alt=\"\" class=\"wp-image-1520\" width=\"339\" height=\"106\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-1-\u2013-a-conventional-equivalent-circuit-of-a-single-phase-transformer.png 452w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-1-\u2013-a-conventional-equivalent-circuit-of-a-single-phase-transformer-300x94.png 300w\" sizes=\"auto, (max-width: 339px) 100vw, 339px\" \/><\/a><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 1 \u2013 A conventional equivalent circuit of a single-phase transformer<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Although the <em>T<\/em>-model does not &#8220;disclose the actual distribution of magnetic fluxes&#8221; within the transformer [3], it is widely used to describe its operation under the <em>normal<\/em> conditions and is a dominant transformer model in textbooks and papers.<\/p>\n\n\n\n<p>Shortcomings of the <em>T<\/em>-model are revealed when the transformer core approaches saturation. Examples include inrush current and ferroresonance events, transformer operations under geomagnetically induced currents and those occurring in generator transformers during an out-of-phase synchronization. Each of these operations has its own characteristics and is considered separately. In this paper, we confine ourselves to considering inrush currents, which remain a subject of ongoing publications. Unfortunately, a lot of these papers focus on secondary details, mostly on hysteresis modeling [4], [5] and miss the fact that using the <em>T<\/em>-model in the inrush current evaluation is fraught with serious flaws.<\/p>\n\n\n\n<p>The same can be said about the models of three-phase transformers [6], [7], [8], which neglect the recommendation in [9] that each wound limb must be shunted by an air branch (the latter divides the iron branches into legs and yokes).<\/p>\n\n\n\n<p>The first disadvantage of the <em>T<\/em>-model is the artificial subdivision of the leakage transformer inductance <em>L<\/em><sub>S<\/sub> (=&nbsp;<em>L<\/em><sub>S1<\/sub>&nbsp;+&nbsp;<em>L<\/em>\u2032<sub>S2<\/sub>) into individual (physically non-existing) &#8220;leakage inductances&#8221; of individual windings. As will be shown in Section 2, such a subdivision can lead to underestimated inrush current peaks calculated for the inner winding and to overestimated current peaks for the outer winding.<\/p>\n\n\n\n<p>In the shell and core type transformers (their legs and yokes are saturated differently), these errors are masked by another inherent drawback of the <em>T<\/em>-model, namely by the fact it contains only a <em>single<\/em> magnetization branch. To<em> circumvent this masking effect<\/em> and to better explain the errors due to the leakage subdivision, a toroidal core transformer is first considered in Section 2. Section 3 explains the necessity of introducing an additional magnetizing branch, which converts the <em>T<\/em>-model into a \u03c0-type structure devoid of the mentioned disadvantages.<\/p>\n\n\n\n<p>Although this paper is devoted to generalized transformer models, simple analytical estimates of inrush current peaks are provided in Section 3 to facilitate subsequent explanations.<\/p>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center\">2 Toroidal transformer model<\/h3>\n\n\n\n<p>Our intention in this section is to show the errors in inrush current predictions due to the &#8220;leakage subdivision&#8221;. To avoid the problem of different magnetization levels of the leg and yokes, consider a transformer with a <em>single<\/em> magnetization branch, that is, a long toroidal core with two infinitely thin windings wound tightly over the entire core. The absence of the yokes and a large length <em>l<\/em> of the core allow one to straighten it out mentally into a ferromagnetic saturable rod, Fig.&nbsp;2, whose ends are joined by a magnetic closure having an infinitely large magnetic permeability&nbsp;(\u00b5&nbsp;= \u221e).<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-2-\u2013-toroidal-transformer-unwrapped-into-a-straight-wound-rod.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-2-\u2013-toroidal-transformer-unwrapped-into-a-straight-wound-rod.png\" alt=\"\" class=\"wp-image-1521\" width=\"368\" height=\"306\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-2-\u2013-toroidal-transformer-unwrapped-into-a-straight-wound-rod.png 491w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-2-\u2013-toroidal-transformer-unwrapped-into-a-straight-wound-rod-300x249.png 300w\" sizes=\"auto, (max-width: 368px) 100vw, 368px\" \/><\/a><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 2 \u2013 Toroidal transformer &#8220;unwrapped&#8221; into a straight wound rod<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Neglecting the winding resistance and assuming an idealized voltage source, , switched to the inner or outer winding (each having <em>N<\/em> turns), the inrush current peak <em>I<\/em><sub>m<\/sub> in the device of Fig. 2 can be calculated fairly accurately using inductance <em>L<\/em><sub>sat<\/sub> of the completely saturated wound core [10], [11], [12]:<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"alignright size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/formula-1-2.png\" alt=\"\" class=\"wp-image-1522\" width=\"481\" height=\"70\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/formula-1-2.png 641w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/formula-1-2-300x44.png 300w\" sizes=\"auto, (max-width: 481px) 100vw, 481px\" \/><\/figure><\/div>\n\n\n<p>where<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"alignright size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/formula-2-1.png\" alt=\"\" class=\"wp-image-1523\" width=\"434\" height=\"52\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/formula-2-1.png 578w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/formula-2-1-300x36.png 300w\" sizes=\"auto, (max-width: 434px) 100vw, 434px\" \/><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Here, <em>A<\/em> is the cross section of the core steel, <em>B<\/em><sub>0<\/sub> is its initial flux density; <em>B<\/em><sub>S<\/sub> (= 2.0 \u2012 2.03 T) is the saturation flux density of grain-oriented steel; <em>h<\/em> and <em>d<\/em><sub>eqv<\/sub> are the height and mean diameter of the excited winding.<\/p>\n\n\n\n<p>It should be recalled that the leakage inductance <em>L<\/em><sub>S<\/sub> is determined by the width <em>d<\/em><sub>12<\/sub> of the leakage channel between the inner winding 1 and the outer winding 2.<\/p>\n\n\n\n<p>When the excited coil is the <em>inner<\/em> winding 1 (when <em>d<\/em><sub>eqv<\/sub> = <em>D<\/em><sub>1<\/sub>), the open-circuit outer winding 2 &#8220;does not exist&#8221; [13] and the leakage channel does not participate in the magnetization process. This means that inductance <em>L<\/em><sub>S<\/sub> is not involved in the inrush current evaluation and <em>no<\/em> part of <em>L<\/em><sub>S<\/sub> can be included in the current path between terminals 1 and 0 of the inner winding in Fig. 3. However, a necessary element of the magnetization branch is inductance <em>L<\/em><sub>01<\/sub> which represents the insulation clearance <em>d<\/em><sub>01<\/sub> between the inner winding 1 and the grounded core.<\/p>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-3-\u2013-an-equivalent-circuit-of-a-toroidal-transformer.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-3-\u2013-an-equivalent-circuit-of-a-toroidal-transformer.png\" alt=\"\" class=\"wp-image-1525\" width=\"287\" height=\"168\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-3-\u2013-an-equivalent-circuit-of-a-toroidal-transformer.png 382w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-3-\u2013-an-equivalent-circuit-of-a-toroidal-transformer-300x176.png 300w\" sizes=\"auto, (max-width: 287px) 100vw, 287px\" \/><\/a><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 3 \u2013 An equivalent circuit of a toroidal transformer<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>When the excited coil is the outer winding 2, i.e. <em>d<\/em><sub>eqv<\/sub> = <em>D<\/em><sub>2<\/sub>, the space between winding 2 and the core includes the leakage channel (of width <em>d<\/em><sub>12<\/sub>) and the innermost clearance (of width <em>d<\/em><sub>01<\/sub>). Therefore, the inrush current path between terminals 2 and 0 in Fig. 3 should include the <em>whole<\/em> inductance <em>L<\/em><sub>S<\/sub> plus the just mentioned inductance <em>L<\/em><sub>01<\/sub>.<\/p>\n\n\n\n<p>It should be clear that moving any part of inductance <em>L<\/em><sub>S<\/sub> into the loop of the inner winding 1, that is, returning to the model in Fig. 1 would result in underestimation of the true current peak in the inner winding 1 and, respectively, in overestimation of the current peak in the outer winding 2.<\/p>\n\n\n\n<p>For the <em>toroidal<\/em> transformer considered in this section, the <em>\u0393<\/em>-type model in Fig. 3 is equally suitable for inrush current evaluation and modeling the rated operation. This is because the core length <em>l<\/em> used in the core loss evaluation is equal to the winding height <em>h<\/em> employed in (2). Quite different situation takes place in transformers with legs and yokes considered in the next section.<\/p>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center\">3 The Core and Shell Type Transformers<\/h3>\n\n\n\n<p>It is crucial to note that the core length <em>l<\/em> employed at the transformer design stage (it is used for the core loss evaluation) is several times greater than the height <em>h<\/em> of the winding used in the inrush current prediction. This makes inductance <em>L<\/em><sub>m<\/sub> in the general purpose models of Figs. 1 and 3 appreciably smaller than inductance <em>L<\/em><sub>sat<\/sub> in (2). Therefore, these models would always overestimate the inrush current significantly regardless of which winding is energized. The reason can be explained by different saturation depths of the wound leg and yokes.<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-4-\u2013-main-magnetic-flux-paths-in-a-shell-type-transformer-with-two-concentric-windings.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-4-\u2013-main-magnetic-flux-paths-in-a-shell-type-transformer-with-two-concentric-windings.png\" alt=\"\" class=\"wp-image-1526\" width=\"357\" height=\"310\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-4-\u2013-main-magnetic-flux-paths-in-a-shell-type-transformer-with-two-concentric-windings.png 476w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-4-\u2013-main-magnetic-flux-paths-in-a-shell-type-transformer-with-two-concentric-windings-300x260.png 300w\" sizes=\"auto, (max-width: 357px) 100vw, 357px\" \/><\/a><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 4 \u2013 Main magnetic flux paths in a shell-type transformer with two concentric windings<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>When a transformer is energized, the core leg can reach very deep saturation. In this case, the magnetic flux density in the wound leg and the adjacent air channel(s) can reach values of the order of 3 T and higher. This is due to the fact that most of the leg is in a &#8220;tube&#8221; formed by the magnetizing winding. Regardless of the degree of saturation of the leg, no or a small part of its flux can escape from the tube. <\/p>\n\n\n\n<p>The saturation of the side branches (SBs) of the core (they compose of the lateral limbs and the related yokes) remain moderate in the sense that the flux density in the SBs does not usually exceed 2 T. This is due to the fact that, as the SBs saturate, the magnetic permeability of their material approaches the vacuum permeability (\u00b5<sub>0<\/sub>), resulting in arising shunting magnetic fluxes, designated as <em>\u0424<\/em><sub>0<\/sub>\u00a0in Fig. 4. The lengths of these air paths are shorter than that of the saturated SBs. Thus, the shunting fluxes <em>\u0424<\/em><sub>0<\/sub>\u00a0keeps the SBs from further saturation. Regarding the space occupied by <em>\u0424<\/em><sub>0<\/sub>, it depends on the winding which is energized. When the exited winding is the inner one, the flux in the leakage channel flows downward and merges with the flux flowing beyond the windings. When exiting the outer winding, the flux in the leakage channel flows upward and merges with flux <em>\u03a6<\/em><sub>01<\/sub>.<\/p>\n\n\n\n<p>These flux paths are reflected in the magnetic network in Fig. 5 [12] where <em>F<\/em><sub>1<\/sub> and <em>F<\/em><sub>2<\/sub> are sources of magneto-motive force of the inner and outer windings, respectively. The shaded elements represent hysteretic reluctances of the leg (\u211c<sub>leg<\/sub>) and side branches (\u211c<sub>SB<\/sub>), whereas unshaded elements designate linear reluctances of the paths in air (the subscripts of the \u211c<sub>01<\/sub>, \u211c<sub>12<\/sub> and \u211c<sub>0<\/sub> are the same as those of the magnetic fluxes).<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-5-\u2013-magnetic-equivalent-networks-of-a-two-winding-trans-former.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-5-\u2013-magnetic-equivalent-networks-of-a-two-winding-trans-former.png\" alt=\"\" class=\"wp-image-1527\" width=\"261\" height=\"235\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-5-\u2013-magnetic-equivalent-networks-of-a-two-winding-trans-former.png 348w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-5-\u2013-magnetic-equivalent-networks-of-a-two-winding-trans-former-300x270.png 300w\" sizes=\"auto, (max-width: 261px) 100vw, 261px\" \/><\/a><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 5 \u2013 Magnetic equivalent networks of a two-winding transformer<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Using the duality transformation [14], [15], the magnetic network in Fig. 5 is converted into its electrical equivalent shown between points 1 and 2 in Fig. 6.<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-6-\u2013-p-shaped-electric-equivalent-circuit-of-a-two-winding-transformer.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-6-\u2013-p-shaped-electric-equivalent-circuit-of-a-two-winding-transformer.png\" alt=\"\" class=\"wp-image-1528\" width=\"314\" height=\"170\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-6-\u2013-p-shaped-electric-equivalent-circuit-of-a-two-winding-transformer.png 419w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-6-\u2013-p-shaped-electric-equivalent-circuit-of-a-two-winding-transformer-300x163.png 300w\" sizes=\"auto, (max-width: 314px) 100vw, 314px\" \/><\/a><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 6 \u2013 \u03c0-shaped electric equivalent circuit of a two-winding transformer<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>In deep saturation, the final slope <em>dB<\/em>\/<em>dH<\/em> of <em>B<\/em>&#8211;<em>H<\/em> curve of any material from the DHM library is equal to the magnetic constant \u00b5<sub>0<\/sub> = 4\u03c0\u00d710\u22127 H\/m. It is important that the segment with slope \u00b5<sub>0<\/sub> lies at the level of the technical saturation of transformer steel (near 2 T), but not at 1.7-1.8 T sometimes taken carelessly. Therefore, the model behavior in saturation is determined by inductances <em>L<\/em><sub>01<\/sub> and <em>L<\/em><sub>0<\/sub>. Their proper choice makes the model <em>reversible<\/em>, that is, capable of replicating inrush currents measured on terminals of <em>both<\/em> the windings.<\/p>\n\n\n\n<p>At this point one should be cautious about uncritical use of \u03c0-models with equal magnetizing inductances of the leg and yoke branches. (This model version is suggested for the case of insufficient information about transformer design [16]). This simplification, as well as the absence of inductances <em>L<\/em><sub>01<\/sub> and <em>L<\/em><sub>0<\/sub>, makes the \u03c0-model symmetrical that contradicts different magnetization curves for the leg and yoke in saturation first observed in [17] and then reprinted in Fig. 6.16 of [18].<\/p>\n\n\n\n<p>The difference between two magnetizing branches at large flux linkages results in different magnetization curves and inrush current values determined for terminals of the whole transformer. Such curves calculated for the transformer in [12] are shown in Fig. 7 where flux linkages \u03a8<sub>1<\/sub> and \u03a8<sub>2<\/sub> are integrals of the voltages at points 1 and 2 of the model in Fig. 6.<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-7-\u2013-magnetization-curves-calculated-from-terminals-of-inner-and-outer-windings.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-7-\u2013-magnetization-curves-calculated-from-terminals-of-inner-and-outer-windings.png\" alt=\"\" class=\"wp-image-1529\" width=\"471\" height=\"320\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-7-\u2013-magnetization-curves-calculated-from-terminals-of-inner-and-outer-windings.png 628w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-7-\u2013-magnetization-curves-calculated-from-terminals-of-inner-and-outer-windings-300x204.png 300w\" sizes=\"auto, (max-width: 471px) 100vw, 471px\" \/><\/a><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 7 \u2013 Magnetization curves calculated from terminals of inner and outer windings<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>In their saturation regions, the \u03a8\u2013<em>i<\/em> curves in Fig. 7 can be characterized by differential inductances, <em>L<\/em><sub>sat1<\/sub> = <em>d<\/em>\u03a8<sub>1<\/sub>\/<em>di<\/em> and <em>L<\/em><sub>sat2<\/sub> = <em>d<\/em>\u03a8<sub>2<\/sub>\/<em>di<\/em>. The methods of their measurements are proposed in [19], [20] and [21] where inequality <em>L<\/em><sub>sat2<\/sub> &gt; <em>L<\/em><sub>sat1<\/sub> was clearly demonstrated.<\/p>\n\n\n\n<p>To account for the winding thicknesses, the model in Fig. 6 is provided by two groups of three inductances as shown in Fig.&nbsp;8. This technique was proposed in [22] and then used in [23], [24] where inductances <em>L<\/em><sub>p<\/sub> are shown to be proportional to one sixth of the thickness of the corresponding winding.<\/p>\n\n\n\n<p>When introducing inductances <em>L<\/em><sub>p<\/sub>, the values of <em>L<\/em><sub>01<\/sub> and <em>L<\/em><sub>0<\/sub> of the model are recalculated so as to maintain the reversibility of the model. This keeps unchanged the magnetization curves in Fig. 7 and thus inrush currents.<\/p>\n\n\n\n<p>In concluding the paper, we can evaluate the influence of the hysteresis properties of the core material on the transformer behavior during inrush current events. For definiteness, the transformer is demagnetized thoroughly before each energization. Fig. 9 shows inrush currents calculated with the models in Fig. 6 and 8 excited from the inner and outer windings (assuming opposite winding open-circuited).<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-8-\u2013-p-equivalent-circuit-of-a-two-winding-transformer-which-accounts-for-the-winding-thicknesses.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-8-\u2013-p-equivalent-circuit-of-a-two-winding-transformer-which-accounts-for-the-winding-thicknesses.png\" alt=\"\" class=\"wp-image-1530\" width=\"413\" height=\"229\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-8-\u2013-p-equivalent-circuit-of-a-two-winding-transformer-which-accounts-for-the-winding-thicknesses.png 550w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-8-\u2013-p-equivalent-circuit-of-a-two-winding-transformer-which-accounts-for-the-winding-thicknesses-300x166.png 300w\" sizes=\"auto, (max-width: 413px) 100vw, 413px\" \/><\/a><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 8 \u2013 \u03c0-equivalent circuit of a two-winding transformer, which accounts for the winding thicknesses<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n<div class=\"wp-block-image\">\n<figure class=\"aligncenter size-full is-resized\"><a href=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-9-\u2013-inrush-currents-on-terminal-of-the-inner-and-outer-windings.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-9-\u2013-inrush-currents-on-terminal-of-the-inner-and-outer-windings.png\" alt=\"\" class=\"wp-image-1531\" width=\"472\" height=\"244\" srcset=\"https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-9-\u2013-inrush-currents-on-terminal-of-the-inner-and-outer-windings.png 629w, https:\/\/fti.dp.ua\/conf\/wp-content\/uploads\/sites\/14\/2023\/05\/figure-9-\u2013-inrush-currents-on-terminal-of-the-inner-and-outer-windings-300x155.png 300w\" sizes=\"auto, (max-width: 472px) 100vw, 472px\" \/><\/a><\/figure><\/div>\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>Figure 9 \u2013 Inrush currents on terminal of the inner and outer windings<\/p>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<p>It was found with certainty that both the current peaks in Fig. 9 and the rate of their decay practically do not depend on the hysteresis properties of the core and on the core loss at all. This is consistent with conclusions made in the recent modeling of ferroresonance processes in small voltage transformers [25] and in studying effects of geomagnetically induced currents in large 400 MVA units [26]. At the same time, accounting for material saturation is found necessary in all of these cases. For this reason, the saturation curve is included in all materials contained in the DHM library.<\/p>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center\">4 Conclusion<\/h3>\n\n\n\n<p>We have attracted attention to typical mistakes made in modeling transformer inrush currents. It is explained that the main mistake is the use of the convenient <em>T<\/em>-model and their three-phase derivatives. Their disadvantages are rooted in using separate leakage models of primary and secondary circuits, as well as in inability to reproduce different magnetization levels of the legs and yoke. We show the benefits of the \u03c0-model, but at the same time, caution against its oversimplification. Although we use an advanced hysteresis model to describe processes in the legs and yokes, it is shown that accounting hysteresis and core losses is completely optional in modeling inrush currents. These model properties have no visible impact neither on the current peaks nor their decay in time.<\/p>\n\n\n\n<h3 class=\"wp-block-heading has-text-align-center citation_references\">References<\/h3>\n\n\n\n<ol class=\"wp-block-list\" type=\"1\">\n<li>H.Kr. H\u00f8idalen, L. Prikler, F. Pe\u00f1aloza, ATPDraw, version 7.0 for Windows, User&#8217;s Manual, 2019.<\/li>\n\n\n\n<li>S. E. Zirka, Y. I. Moroz, N. Chiesa, R. G. Harrison, and H. Kr. Hoidalen, &#8220;Implementation of inverse hysteresis model into EMTP \u2013 Part II: Dynamic model,&#8221; <em>IEEE Trans. Power Del<\/em>., vol. 30, no. 5, pp. 2233-3241, Oct. 2015. <strong>doi:&nbsp;<\/strong>10.1109\/TPWRD.2015.2416199.<\/li>\n\n\n\n<li>A. Boyajian, Impedance characteristics of transformers. Ch. IV in Transformer Engineering L.F.&nbsp;Blume, G. Camilli, A. Boyajian, V.M. Montsinger, N.Y. John Wiley &amp; Sons, 1938.<\/li>\n\n\n\n<li>H. Altun, S. S\u00fcnter,\u00b7\u00d6. Aydo\u011fmu\u015f, &#8220;Modeling and analysis of a single-phase core-type transformer under inrush current and nonlinear load conditions,&#8221; <em>Electrical Engineering<\/em>.<em> doi:<\/em><em> <\/em>10.1007\/s00202-021-01283-9.<\/li>\n\n\n\n<li>A. Toki\u0107, I. Ugle\u0161i\u0107, and G. \u0160tumberger, &#8220;Simulations of transformer inrush current by using BDF-based numerical methods,&#8221; Hindawi Publishing Corp., Math. Problems in Eng., vol. 2013, Article ID 215647. http:\/\/dx.doi.org\/10.1155\/2013\/215647.<\/li>\n\n\n\n<li>J. Pedra, L. Sainz, F. C\u00f3rcoles, R. Lopez, and M. Salichs, &#8220;PSPICE computer model of a nonlinear three-phase three-legged transformer,&#8221; <em>IEEE Trans. Power Del<\/em>., vol. 19, no. 1, Jan. 2004. pp. 200-207. doi: 10.1109\/TPWRD.2003.820224.<\/li>\n\n\n\n<li>A. D. Theocharis, J. Milias-Argitis, and T. Zacharias, &#8220;Three-phase transformer model including magnetic hysteresis and eddy currents effects,&#8221; <em>IEEE Trans. Power Del<\/em>., vol. 24, no. 3, pp. 1284\u20131294, Jul. 2009. doi<strong>:&nbsp;<\/strong>10.1109\/TPWRD.2009.2022671.<\/li>\n\n\n\n<li>P.&nbsp;S. Moses, M.&nbsp;A.&nbsp;S. Masoum, and H.&nbsp;A. Toliyat, &#8220;Dynamic modeling of three-phase asymmetric power transformers with magnetic hysteresis: No-load and inrush conditions,&#8221; <em>IEEE Trans. Energy Convers.<\/em>, vol. 25, no. 4, pp. 1040\u20131047, Dec. 2010.<strong> doi: <\/strong>10.1109\/TEC.2010.2065231.<\/li>\n\n\n\n<li>C.&nbsp;M. Arturi, &#8220;Transient simulation of a three-phase five-limb step-up transformer following an out-of-phase synchronization,&#8221; <em>IEEE Trans. Power Del<\/em>., vol. 6, no. 1, pp. 196\u2013207, Jan. 1991. doi: 10.1109\/61.103738.<\/li>\n\n\n\n<li>N. Chiesa, B.&nbsp;A. Mork, H.&nbsp;K. H\u00f8idalen, &#8220;Transformer model for inrush current calculations: Simulations, measurements and sensitivity analysis,&#8221; <em>IEEE Trans. Power Del<\/em>., vol. 25, no. 4, 2010, pp. 2599\u20132608. doi: 10.1109\/TPWRD.2010.2045518.<\/li>\n\n\n\n<li>S.&nbsp;V. Kulkarni and S.&nbsp;A. Khaparde, <em>Transformer Engineering: Design and Practice<\/em>. New York: Marcel Dekker, 2004.<\/li>\n\n\n\n<li>S. E. Zirka, Y. I. Moroz, C. M. Arturi, N. Chiesa, and H. K. H\u00f8idalen, \u201cTopology-correct reversible transformer model,\u201d <em>IEEE<\/em> <em>Trans. Power Del<\/em>., vol. 27, no. 4, pp. 2037\u20132045, Oct. 2012. doi: 10.1109\/TPWRD.2012.2205275.<\/li>\n\n\n\n<li>S. Jazebi, F. de Leon, N. Wu. &#8220;Enhanced analytical method for the calculation of the maximum inrush currents of single-phase power transformers,&#8221; <em>IEEE<\/em> <em>Trans. Power Del<\/em>., vol. 30, no. 6, pp. 2590\u20132599, Dec. 2015. doi: 10.1109\/TPWRD.2015.2443560.<\/li>\n\n\n\n<li>E.&nbsp;C. Cherry, &#8220;The duality between interlinked electric and magnetic circuits and the formation of transformer equivalent circuits,&#8221; Proc. Phys. Soc. B 62. pp. 101\u2013111, 1949. doi: 10.1088\/0370-1301\/62\/2\/303.<\/li>\n\n\n\n<li>S. Jazebi, S.&nbsp;E. Zirka, M.&nbsp;Lambert, A. Rezaei-Zare, <em>et al<\/em>., &#8220;Duality derived transformer models for low-frequency electromagnetic transients\u2014Part I: Topological models,&#8221; <em>IEEE Trans. Power Del<\/em>., vol. 31, no. 5, pp. 2410-2419, 2016. doi: 10.1109\/TPWRD.2016.2517327.<\/li>\n\n\n\n<li>F. de Le\u00f3n, P. G\u00f3mez, J.&nbsp;A. Martinez-Velasco, and M. Rioual, \u201cTransformers,\u201d Chapter 4 in <em>Power System Transients<\/em>. Parameter Determination, J.&nbsp;A. Martinez-Velasco (ed.), CRC Press, 2009.<\/li>\n\n\n\n<li>M.&nbsp;Kh.&nbsp;Zikherman, &#8220;Magnetization characteristics of large transformers,&#8221; <em>Elektrichestvo<\/em>, no. 3, pp. 79-82, 1972. [in Russian].<\/li>\n\n\n\n<li>H.W. Dommel, <em>Electromagnetic transients program reference manual<\/em> (<em>EMTP Theory Book<\/em>); Bonneville Power Administration, Portland, OR, USA, 1986.<\/li>\n\n\n\n<li>S. Calabro, F. Coppadoro and S. Crepaz, &#8220;The measurement of the magnetization characteristics of large power transformers and reactors through DC excitation&#8221;, <em>IEEE Trans. Power Del.<\/em>, vol. 1, no. 4, pp. 224-232, Oct. 1986. doi: 10.1109\/TPWRD.1986.4308052.<\/li>\n\n\n\n<li>S. Jazebi<em>, <\/em>F. de Le\u00f3n<em>, <\/em>A. Farazmand, and D. Deswal, &#8220;Dual reversible transformer model for the calculation of low-frequency transients,&#8221; <em>IEEE<\/em> <em>Trans. Power Del<\/em>., vol. 28, no. 4, pp. 2509-2517, Oct. 2013. doi: 10.1109\/TPWRD.2013.2268857.<\/li>\n\n\n\n<li>W. Sima, B. Zou, M. Yang, F. de Le\u00f3n, &#8220;New method to measure deep-saturated magnetizing inductances for dual reversible models of single-phase two-winding transformers,&#8221; <em>IEEE<\/em> <em>Trans. Power Del<\/em>., vol. 36, no. 1, pp. 488-491, Feb. 2021. doi: 10.1109\/TEC.2020.3023596.<\/li>\n\n\n\n<li>H.&nbsp;Edelmann, &#8220;Descriptive determination of&nbsp; transformer equivalent circuits,&#8221; Archiv f\u00fcr elektrische \u00dcbertragung, vol. 13, no. 6, pp. 253-261, 1959. [in German].<\/li>\n\n\n\n<li>S.&nbsp;E. Zirka, Y.&nbsp;I. Moroz, C.&nbsp;M. Arturi, &#8220;Accounting for the influence of the tank walls in the zero-sequence topological model of a three-phase, three-limb transformer,&#8221; <em>IEEE Trans Power Del<\/em>., vol. 29, no. 5, pp. 2172\u20132179, Oct. 2014. doi: 10.1109\/TPWRD.2014.2307117.<\/li>\n\n\n\n<li>J.&nbsp;Zhao, S.&nbsp;E. Zirka, Y.&nbsp;I. Moroz, C.&nbsp;M. Arturi, &#8220;Structure and properties of the hybrid and topological transformer models,&#8221; <em>Int. J. Electr. Power &amp; Energy Systems<\/em>, vol. 118, 2020, 105785. doi: 10.1016\/j.ijepes.2019.105785.<\/li>\n\n\n\n<li>S.&nbsp;E. Zirka, Y.&nbsp;I. Moroz, A.&nbsp;V. Zhuykov, D.&nbsp;A. Matveev, <em>et al<\/em>., &#8220;Eliminating VT uncertainties in modeling ferroresonance phenomena caused by single phase-to-ground faults in isolated neutral network,&#8221; <em>Int. J. Electr. Power &amp; Energy Systems<\/em>, 2021, vol. 133. Paper 107275. doi.org\/10.1016\/j.ijepes.2021.107275.<\/li>\n\n\n\n<li>S.&nbsp;E. Zirka, Y.&nbsp;I. Moroz, J. Elovaara, M. Lahtinen, <em>et al<\/em>., &#8220;Simplified models of three-phase, five-limb transformer for studying GIC effects,&#8221; <em>Int. J.<\/em> <em>Electr. Power &amp; Energy Systems<\/em>, vol. 103, pp. 168-175, 2018. https:\/\/doi.org\/10.1016\/j.ijepes.2018.05.035.<\/li>\n<\/ol>\n\n\n\n<div style=\"height:3em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity is-style-default\"\/>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<div class=\"wp-block-group is-vertical is-content-justification-right is-layout-flex wp-container-core-group-is-layout-b6c475e2 wp-block-group-is-layout-flex\">\n<div class=\"wp-block-group is-content-justification-right is-nowrap is-layout-flex wp-container-core-group-is-layout-fd526d70 wp-block-group-is-layout-flex\"><div class=\"taxonomy-post_tag wp-block-post-terms\"><a href=\"https:\/\/fti.dp.ua\/conf\/tag\/cims-2023\/\" rel=\"tag\">CIMS 2023<\/a><\/div>\n\n<div class=\"wp-block-post-date\"><time datetime=\"2023-05-21T21:32:51+03:00\">May 21, 2023<\/time><\/div><\/div>\n\n\n<div class=\"taxonomy-category wp-block-post-terms\"><a href=\"https:\/\/fti.dp.ua\/conf\/session\/electronics\/\" rel=\"tag\">Electronics and telecommunications<\/a><\/div>\n\n\n<p><em>URI: <a href=\"https:\/\/fti.dp.ua\/conf\/2023\/05216-1930\/\">https:\/\/fti.dp.ua\/conf\/2023\/05216-1930\/<\/a><\/em><\/p>\n<\/div>\n\n\n\n<div style=\"height:1em\" aria-hidden=\"true\" class=\"wp-block-spacer\"><\/div>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity is-style-default\"\/>\n","protected":false},"excerpt":{"rendered":"<p>On Typical Flaws of the Transformer Models for Inrush Current Evaluation Sergey Zirka ORCID: https:\/\/orcid.org\/0000-0001-7607-1436 Oles Honchar Dnipro National University, Dnipro Yuriy Moroz ORCID:&nbsp;https:\/\/orcid.org\/0000-0001-7696-1324 Oles Honchar Dnipro National University, Dnipro 1 Introduction For many decades, the main, if not the only, equivalent circuit of a single-phase transformer was its T-model, Fig. 1. To avoid ambiguity and take into account the topology of the core and windings, we will use the terms inner and outer windings, instead of the primary and &hellip; <\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[20],"tags":[23],"class_list":["post-1518","post","type-post","status-publish","format-standard","hentry","category-electronics","tag-cims-2023"],"_links":{"self":[{"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/posts\/1518","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/comments?post=1518"}],"version-history":[{"count":2,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/posts\/1518\/revisions"}],"predecessor-version":[{"id":2323,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/posts\/1518\/revisions\/2323"}],"wp:attachment":[{"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/media?parent=1518"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/categories?post=1518"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/fti.dp.ua\/conf\/wp-json\/wp\/v2\/tags?post=1518"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}