Research Article  Open Access
Flux Reflection Model of the Ferroresonant Circuit
Abstract
The paper presents a linear model of ferroresonant circuit with flux reflection. The proposed model—flux reflection model—derives from observations of typical flux waveforms of nonlinear coil during ferroresonant steady states. Simulation results of the flux reflection model are compared with simulation results of the usual nonlinear model as well as with measurements carried out on the physical model of the ferroresonant circuit. The flux reflection model enables a novel comprehension of the ferroresonant circuit behavior and simplifies the modeling of the nonlinear coil in the ferroresonant circuit.
1. Introduction
A ferroresonant circuit is a circuit which comprises a linear capacitor in series with a nonlinear coil, driven by a sinewave voltage. The ferroresonant circuit can exhibit several bifurcations and a chaotic steady state [1–5]. Common research procedures comprise the measurements of ferroresonant circuit and the computer simulation of the mathematical model based on stateequations of the circuit.
The nonlinearity of the model caused by the nonlinear magnetization characteristic of the coil is the driving force on the route to the chaotic steady state, which is marked by the bifurcations that are initiated by varying values of a circuit parameter. Complex models, that comprehend more nonlinear elements, are needed to match results of measurements and simulation completely [6, 7]. In order to obtain the complex behavior of the circuit, the nonlinear element that comprises the magnetization nonlinearity is only necessary. Therefore, the model comprises the magnetization nonlinearity , implying the core saturation of the coil. The winding resistance of the coil is neglected, and the resistor, that represents the core losses of the coil, is assumed to be linear. State equations of the ferroresonant circuit, Figure 1, are
The magnetization nonlinearity is oddsymmetric and monotonically increasing. It can be presented in various forms, for example, in a polynomial form [8, 9] or piecewise linear form [10, 11]. Thereby, the polynomial form is obtained often by an interpolation of a piecewise linear form. The magnetization nonlinearity is based on measurements carried out on a nonlinear coil regarding the abovementioned assumptions made about properties of the coil. For example, one ought to rely on the results of standard measurements carried out on the nonlinear coil by the manufacturer [11, 12]. Because of a limited thermal endurance of the coil, measurements cannot be carried out by applying the sinewave voltage for the full range of flux and voltage peak values which occur in some operation modes of ferroresonant circuit. However, the importance of this problem is often ignored by a simple extrapolation of magnetization nonlinearity obtained by measurements.
Further problem of modeling based on standard measurements is that the parameters of model depend on the number of measurements, that is, on the number of measured values on UI characteristic. Namely, a piecewise linear form of the magnetization nonlinearity, that is obtained by a small number of measurements, would not be as smooth as necessary; in the case of polynomial form of the magnetization nonlinearity, a small number of measurements could result with a significant interpolation error.
In the paper a novel kind of modeling the magnetization characteristic of the coil in the ferroresonant circuit is presented. The modeling will be based on a characteristic behavior that nonlinear coil exhibits in a ferroresonant circuit. The behavior will be identified from the flux waveforms obtained by measurements and by computer simulation.
The preliminary purpose of obtained novel model of the ferroresonant circuit is to enable a new comprehension of the circuit behavior and to simplify the modeling of nonlinear coil in the ferroresonant circuit by reducing the number of model parameters.
2. Measurements and Nonlinear Model
In order to notice the characteristic behavior of a coil in ferroresonant circuit, the ferroresonant steady states are obtained by measurements carried out on a ferroresonant circuit that is realized in laboratory and by computer simulation carried out on a typical nonlinear model of the ferroresonant circuit.
The ferroresonant circuit realized in the laboratory is composed of the capacitor and the primary winding of the toroidal ironcored twowindings transformer used as a nonlinear coil. The transformer was designed for the nominal apparent power of 200 VA and for the nominal primary voltage of 30 V. The core is stripwounded, made of NiFe alloy (Trafoperm N3). The autotransformer of 10 kVA nominal apparent power is used as a variable voltage source in all experiments.
To simplify calculations, throughout the rest of the paper all the variables and parameters of the ferroresonant circuit are expressed in relation to reference quantities, that is, in a perunit system:
The simulation is carried out using a typical nonlinear model of ferroresonant circuit, with the polynomial form of magnetization nonlinearity:
The magnetization characteristic and the ironcore losses are derived from measured PU and UI characteristic of the nonlinear coil [12].
During measurements and simulation the steady states are obtained by varying the amplitude of source voltage within the range:
The step of the amplitude variation is in simulation, as well as in measurements. Simulation is carried out by using the 4th order RungeKutta method with 10000 integration steps per period, .
Characteristic steady states and bifurcations obtained by simulation and measurements are shown in Tables 1 and 2, respectively. There is no significant disagreement between results of simulation and measurements with the exception of the pitchfork bifurcation and the periodfour steady state, which are not noticed in measurements.


By that means, Figure 3 implies that the peak value of flux is in all polyharmonic steady states limited by the value , called here as a flux saturation value:
From simulation as well as measurements it is observed that during all steady states the slope of flux, immediately before and after the moment each time the flux reaches the peak value , changes the sign only, that is, holding the same absolute value nearly. As it is seen in Figure 2, for instance, at the moment , we assume
(a)
(b)
This property is named here as a flux reflection.
3. Flux Reflection Model
Novel presentation of magnetization nonlinearity is based on observed limitation (2.4) and reflection (2.5) of flux values. Thereby, the magnetization characteristic for values is approximated by a linear inductance . In a case of the polynomial form of magnetization nonlinearity (2.2b) the linear inductance is approximated as shown on Figure 4.
By the reducing of the nonlinear magnetization characteristic to the linear inductance , the stateequations of the ferroresonant circuit, (2.2a) and (2.2b), become linear between two flux reflections, that is, for flux values :
Thereby, the observed properties of a ferroresonant circuit—the limitation of flux values (2.4) and the flux reflection (2.5)—can be realized easily, using the following pseudocode:
The proposed flux reflection model, defined by stateequations (3.1) and pseudocode (3.2), could be comprehended as a piecewise linear model also. However, contrary to a typical piecewise linear model, this model comprises only one set of stateequations, that is, the magnetization characteristic with only one linear segment. As it is shown on Figures 5 and 6 and determined by the pseudocode (3.2), the segment is limited by the saturation value which triggers the alternation of flux slope sign and, at this way, causes the change of flux trend. For instance, the open arrows on Figure 6 indicate the increasing of flux value before the flux reflection and the decreasing of flux values after the flux reflection. Figures 5 and 6 show the flux reflection for the case only. The depiction of the flux reflection for the case can be derived from the mirror symmetry of flux reflection model easily.
(a)
(b)
4. Comparison of the Flux Reflection Model and the Nonlinear Model
The flux reflection model can be considered as a model of the ferroresonant circuit only if the flux reflection model reveals a complex behavior of the ferroresonant circuit. Therefore, the steady states of the flux reflection model are obtained in the same way as for the ferroresonant circuit realized in laboratory and for the nonlinear model in Section 2, that is, by varying the amplitude with the step of the amplitude variation within a range:
Between two flux reflections the flux reflection model is described by linear state equations (3.1). Thus, these equations can be solved analytically. However, as it is shown in the appendix in detail, at the start of transient state, that is, in moment , it is not possible to determine explicitly the ranges of parameter values and values of initial conditions at which a particular steady state would arise; it is only possible to solve (3.1) for a given set of parameter values and values of initial conditions.
In the paper the steady states solutions are obtained in the same way as for the nonlinear model in Section 2, that is, by simulation that is carried out numerically, using the 4th order RungeKutta method with 10000 integration steps per period, .
The route to the chaotic steady state with bifurcations and corresponding steady states is obtained by varying the source voltage amplitude , (3.1), as shown on Table 3. Figure 7 shows the flux waveforms corresponding to each of steady states of the flux reflection model shown in Table 3.

(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
The results of simulation based on the nonlinear model and results of measurements, shown on Tables 1 and 2, differ from results of simulation based on the flux reflection model shown on Table 3. For instance, the periodtwo steady state as well as periodfour steady state and corresponding perioddoubling bifurcations, are not obtained by simulation based on flux reflection model. The reason for this could be the approximation of the maximum value of flux by a constant saturation value . Namely, on Figure 3 during the polyharmonic steady states the maximum value of flux is constant only nearly, that is, it is not equal to the chosen saturation value ().
The main shortcoming of the flux reflection model is the inherent limitation of inductive component of the coil current
which does not exist in the actual ferroresonant circuit.
In spite of that, the proposed model reveals a complex steady state behavior of the ferroresonant circuit. Besides of a new comprehension of the circuit behavior, the flux reflection model enables a simple modeling of the nonlinear coil. It is sufficient to carry out the simple measurements to determine the value of linear inductance L and the saturation value .
5. Additional Simulation
The additional simulation is carried out in order to investigate the sensitivity of flux reflection model on the value of linear inductance L and on the saturation value . The step of the source voltage amplitude variation is . Simulation is carried out by using the 4th order RungeKutta method with 10000 integration steps per period.
Figure 8 shows results of simulation. For easier visualisation the diagrams are constructed using gray scaled squares where different shades of gray are employed to represent values of source voltage amplitude at which the particular bifurcation and steady state, respectively, has occurred. For instance, the bottom left square on Figure 8(c) () indicates that the pitchfork bifurcation for the parameter values , occur at the source voltage amplitude .
The results of simulation reveal that the flux reflection model is more sensitive to the variation of the saturation value , than to the value of linear inductance L. For instance, for the saturation value and linear inductance values , the reverse ferroresonant jump occurs at the source voltage amplitude values () and (), as it is presented by middle row on Figure 8(b). For the linear inductance value and the saturation values the reverse ferroresonant jump occurs for a wider range of source voltage amplitude values, () and (), as it is presented by middle column on Figure 8(b).
Therefore, it is more important to know the saturation value precisely than the value of linear inductance L.
6. Conclusions
The flux reflection model preserves the basic properties of the magnetization characteristic: linearity and saturation. Despite the noticed shortcomings, the flux reflection model reveals a complex steady state behavior characteristic for the ferroresonant circuit. It enables a new comprehension of the behavior of the ferroresonant circuit and simplifies the modeling of the nonlinear coil in the ferroresonant circuit.
Future work will address the simplification of the identification of parameters of a flux reflection model in order to avoid the identification by using the waveforms of flux values during ferroresonant steady states. Namely, in a flux reflection model a coil is defined by a linear inductance value and by a saturation value, which could be determined by standard measurements experimentally or by number of windings and by core parameters of coil analytically.
Furthermore, in order to expand the purpose of the model from the presented new comprehension of the circuit behavior to a wide application of the model, it will be necessary to identify advantages and disadvantages of the flux reflection model precisely by a more detailed comparison with usual models of ferroresonant circuit.
Appendix
On Analytical Conditions for Steady states
Stateequations (3.1) can be expressed as a secondorder linear differential equation:
If , as it is the case for parameter values used in the paper the solution of the equation can be written as:
The constants K and are determined by values of initial conditions, and . Equation, (A.1), as well as its solution, (A.2), is valid for the flux reflection model only between the reflections, that is, for flux values . In the moment t_{1} the flux reaches the saturation value and the flux slope is alternated, (2.5). After the moment of reflection, the solution, (A.2), is valid again, but through the alternation of flux slope , the initial conditions are changed and, consequently, the constants of the solution.
Let us denote the solution before the first reflection in moment as
The constants and are determined by given values of initial conditions, and . Generally, the solution after n reflections in moments and before the reflection in moment is
Constants and are determined using the initial conditions and , which depend on values of solution and its derivative in the moment of reflection :
Each moment of reflection is determined by the condition
which cannot be expressed explicitly regarding the moment , that is, moments of reflection can be determined only approximately.
Characteristic steady state is established if the following conditions are met according to Table 4:

However, at the start of transient state, that is, in moment , it is impossible to determine explicitly the ranges of parameter values and values of initial conditions, and , at which a particular steady state would arise because a solution , (A.4), can be obtained recursively only, starting from first solution , (A.3), and using approximately determined values of reflection moments , (A.6), in each step of recursion.
References
 C. Hayashi, Nonlinear Oscillations in Physical Systems, McGrawHill Electrical and Electronic Engineering Series, McGrawHill, New York, NY, USA, 1964. View at: Zentralblatt MATH  MathSciNet
 C. Kieny, “Application of the bifurcation theory in studying and understanding the global behavior of a ferroresonant electric power circuit,” IEEE Transactions on Power Delivery, vol. 6, no. 2, pp. 866–872, 1991. View at: Publisher Site  Google Scholar
 B. Lee and V. Ajjarapu, “Perioddoubling route to chaos in an electrical power system,” IEE Proceedings: Generation, Transmission and Distribution, vol. 140, no. 6, pp. 490–496, 1993. View at: Google Scholar
 B. A. Mork and D. L. Stuehm, “Application of nonlinear dynamics and chaos to ferroresonance in distribution sysytems,” IEEE Transactions on Power Delivery, vol. 9, no. 2, pp. 1009–1017, 1994. View at: Publisher Site  Google Scholar
 S. Mozaffari, M. Sameti, and A. C. Soudack, “Effect of initial conditions on chaotic ferroresonance in power transformers,” IEE Proceedings: Generation, Transmission and Distribution, vol. 144, no. 5, pp. 456–460, 1997. View at: Publisher Site  Google Scholar
 M. R. Iravani, A. K. S. Chaudhary, W. J. Giesbrecht et al., “Modeling and analysis guidelines for slow transients—part III: the study of ferroresonance,” IEEE Transactions on Power Delivery, vol. 15, no. 1, pp. 255–265, 2000. View at: Publisher Site  Google Scholar
 J. A. MartinezVelasco and B. A. Mork, “Transformer modeling for low frequency transients: the state of the art,” in Proceedings of the IPST International Conference on Power Systems Transients, New Orleans, La, USA, June 2003. View at: Google Scholar
 B. A. T. Al Zahawi, Z. Emin, and Y. K. Tong, “Chaos in ferroresonant wound voltage transformers: effect of core losses and universal circuit behaviour,” IEE Proceedings: Science, Measurement and Technology, vol. 145, no. 1, pp. 39–43, 1998. View at: Publisher Site  Google Scholar
 D. A. N. Jacobson, P. W. Lehn, and R. W. Menzies, “Stability domain calculations of period1 ferroresonance in a nonlinear resonant circuit,” IEEE Transactions on Power Delivery, vol. 17, no. 3, pp. 865–871, 2002. View at: Publisher Site  Google Scholar
 L. O. Chua, M. Hasler, J. Neirynck, and P. Verburgh, “Dynamics of a piecewiselinear resonant circuit,” IEEE Transactions on Circuits and Systems, vol. 29, no. 8, pp. 535–547, 1982. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. L. A. Neves and H. W. Dommel, “On modelling iron core nonlinearities,” IEEE Transactions on Power Systems, vol. 8, no. 2, pp. 417–425, 1993. View at: Publisher Site  Google Scholar
 I. Flegar, D. Fischer, and D. Pelin, “Identification of chaos in a ferroresonant circuit,” in Proceedings of the International Conference on Electric Power Engineering (PowerTech '99), pp. 1–5, Budapest, Hungary, AugustSeptember 1999. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2009 Kruno Miličević et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.