Microsoft Word - 13-2327_s Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3739-3743 3739 www.etasr.com Khlifi: Explaining the d-q Magnetic Couplings Theoretically in Saturated Smooth Air-Gap AC Machines Explaining the d-q Magnetic Couplings Theoretically in Saturated Smooth Air-Gap AC Machines Mohamed Arbi Khlifi Department of Electrical Engineering, University of Hail, Hail, Saudi Arabia and SIME Laboratory, Ecole Nationale Superieure d' Ingenieurs (ENSIT), University of Tunis, Tunisia mohamedarbi.khlifi@issatm.rnu.tn Abstract—The current paper discusses the theoretically detailed physical analysis of the effects of main flux saturation in smooth air gap AC machines. For that purpose, three sections are presented. The first one deals for presenting equivalence saturation-air gap. The second presents the equivalence of inductances between saturation induction machines and a dismounted salient pole synchronous machine. The method analyzes the leakage inductances that properly represent the corresponding leakage fluxes. The third section treats the physical interpretation of the cross-saturation. It is possible to show that the equations and physical picture in a saturated smooth-air-gap machine are similar to the equations and physical picture in a non-saturated salient-pole machine, where mutual coupling exists between the direct and quadrature axes and the self-inductances are also different on the non-salient part of the machine due to the variation of the reluctance with the rotor angle. Keywords-AC machines; physical interpretation; magnetic saturation I. INTRODUCTION In conventional models used in electrical engineering to represent electrical machines, it is quite common to neglect the magnetic saturation phenomenon, allowing the making of the simplifying assumption that inductances depend on the value of the flux imposed by currents, thanks to the non-linear relation: ( ) ( ) m m m m i L iλ = , where ( ) m m iλ is a picture of the non-linear magnetization characteristic B(H) of the employed ferromagnetic material [1-3]. Saturation is a phenomenon that becomes apparent when the current through an inductor exceeds a certain limit. In effect, saturation of an inductor occurs when the core of the inductor can no longer store magnetic energy [4]. An inductor is usually constructed using a ferromagnetic material such as iron or one of its alloys. A ferromagnetic material is composed of small domains whose dipole moments are randomly oriented. When the material is completely demagnetized, the net magnetic moment over the material is zero. As the current passes through the inductor, the domains in the material tend to align in the direction of the magnetic field intensity [5]. As more domains are aligned, the magnetic flux density increases. When all domains have aligned in the direction of the magnetic field intensity, the flux density reaches its saturation density and cannot increase further in the same fashion. As the magnetic field intensity increases beyond the saturation point, relatively small increases of the magnetic flux density occur. When a large number of magnetic domains align this way, the material is said to be in its saturation region [6]. Saturation representation is critical in the correct modeling of synchronous machines. The main effect of saturation in a synchronous generator is the decrease of its mutual inductances depending on the operating level of the generator. Such a decrease may be considerable as the generator is driven higher into saturation. Therefore, it is imperative that the effect of saturation is modeled in the parameter estimation procedure [7-12]. The representation of saturation is rather cumbersome because it is not possible to represent all the nonlinearities in a generator and the interaction between all circuits within the stator and the rotor (including the fictitious windings) with any single model for all the operating conditions. Therefore, a number of assumptions have to be made in order to enable satisfactory saturation modeling. II. PROBLEM FORMULATION The variation of the level of magnetic saturation sheets can be transformed and described by an air-gap with modulated thickness. Let us deal with the essential component of any revolving electric machine, in this occasion, the iron core coil. The usual structure is the one shown in Figure 1. Fig. 1. Torus without air-gap. The application of the law of Ampere gives: f f H l Ni= (1) In addition, the analysis between electric quantities and magnetic carries out (2): f f F R λ = (2) with Corresponding author: Mohamed Arbi Khlifi Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3739-3743 3740 www.etasr.com Khlifi: Explaining the d-q Magnetic Couplings Theoretically in Saturated Smooth Air-Gap AC Machines F Ni= (3) f f fB Sλ = (4) 0f r fB Hµ µ= (5) Finally, the magnetic reluctance of the iron which crosses the flux is: 0 1 f f r f l R Sµ µ = (6) In addition, one can admit that the same magnetizing force Fcan produce the same flow λf in the not homogeneous circuit of Figure 2 or iron is supposed to be perfect. The ring of the permanent magnet material and the air-gap are excited by a wire coil of N turns and current i. Fig. 2. Permanent magnet as a part of a magnetic circuit consisting of iron and an air gap, which is divided into an air-gap flux and leakage flux. The geometry shown in Figure 2 is an excellent canonical problem to illustrate the salient features of permanent magnet circuits. Assume a small cross-section and ignore all leakage flux. The first topic is to show the steps necessary to determine B in this problem, which is the same in the gap and the magnet because of the boundary condition. Under these conditions, the reluctances R must be equal to their precedents Rf described by (6). That is to say: f R R= (7) The new total reluctance is reduced, it is the one of the air- gap Rα: 0 1 f a a e R Sµ = (8) Since the thickness of the air-gap is usually low, the useful section in the air and iron is unchanged: f air S S S= = (9) The equivalence of (6) and (8) gives: ste f f r r l C e µ µ = = (10) Thus, the thickness of the virtual air-gap ef is inversely proportional to the relative permeability of iron. When the magnetic medium is linear, ef corresponds with a minimal width. A rise in the level of saturation, by increase in current I revolves a fall of µr, synonymous with a widening of the air- gap ef, as it is indicated in Figure 3. Fig. 3. Magnetizing curve, λm, ef, µr. So for magnetic circuit containing a real air-gap, the resulting air-gap will be: e+ef. In the electric machines the excitation i and consequent induction vary sinusoidalement according to time (Figure 4). Fig. 4. Relation sheep between λm, ef, µr. The main flux distribution around the air gap periphery is sinusoidal at any rate of saturation. Also the mmfs and voltages are sinusoidally distributed around the periphery. This assumption is equivalent to the supposition that the permeance of the iron is modeled by a factious variation of the air gap clearance (width). Since the real saturation effect caused by the iron core is to be replaced by supposition of a variable air gap clearance of the threshold of saturation the minimal permeability corresponds to maximum of the magnetomotrice force and the air-gap. III. EXPRESSIONS OF INDUCTANCES IN SALIENT POLE SYCHRONOUS MACHINES A. Formulation of Voltage Equations In undamped three phase salient pole synchronous machines there are four mutually coupled coils: the three phase windings and the field. Therefore, the flux linkage equations for these windings result in the following matrix form: Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3739-3743 3741 www.etasr.com Khlifi: Explaining the d-q Magnetic Couplings Theoretically in Saturated Smooth Air-Gap AC Machines [ ] a ab ac af ba b bc bf ca cb c cf fa fb fc f L L L L L L L L L L L L L L L L L       =         (11) where Ljk is a self-inductance when j=k and a mutual inductance when j k≠ .The stator inductances and the coupled inductances between the stator and the rotor are time varying, while the rotor inductances are constants. Expressions for the inductances in (11) are given below. B. Inductances of the Stator Circuits The inductances can be assimilated with those of a diphase salient pole machine. The inductance of a diphase salient pole synchronous machine is expressed as: ' cos cos(2 )sm svaa L L Lβ θ β= − − (12) with max min max min 2 2 2 sm sv L L L L L L π β + =  − =   =  (13) C. Stator Self-Inductances ( ) ( ) cos 2 cos 2 a sm sv b sm sv L L L L L L θ θ = − = + (14) sm L and sv L are constants and sm sv L L� . The double frequency term (2θ) occurs due to the rotor saliency and the fact that self-inductances are the same for both rotor poles. D. Stator Mutual Inductances The stator to stator mutual inductances are a function of rotor position since they are influenced by rotor saliency. ( )' sin 2svaaL L θ= − (15) It is observed that (11) has time-varying terms which will cause complication when their derivatives are taken. Thus, it is convenient to refer all quantities to a rotor frame of reference through a Park’s transformation. IV. INTERPRETATION AND DISCUSSION A. Interpretation Inductances can be comparable with those of a revolving machine diphase with salient poles whose principal inductances are Lmdy, according to the direct axis (weak air-gap, weak reluctance) and Lm following the axis squaring. The standard mathematical model of an induction machine is given in terms of space vectors (underlined variables), in an arbitrary frame of reference ωα with the following set of equations: s s s s a s d v R i jw dt λ λ= + + (16) ( )rr r r a r d v R i j w dt λ ω λ= + + − (17) where v, i and λ denote voltage, current and flux linkage space vectors respectively, while indices s and r stand for stator and rotor, respectively. Equation of mechanical motion is irrelevant here and is hence omitted. ω is the electrical speed of rotation. It is assumed that leakage flux saturation and main flux saturation can be treated independently. As only main flux saturation is discussed, leakage inductances are constants. Stator and rotor flux linkages are: s s s ml iλ λ= + (18) r r r ml iλ λ= + (19) where index m denotes variables and parameters associated with magnetizing flux. The proposed common approach to introduce magnetic saturation in any d-q existing model, for induction machine, relies only on the knowledge of the winding currents saturated model. For that reason, it will be first shown how to obtain it. Leakage inductances in flux equation of induction machine (16)-(17) are assumed to be constant, only the main flux mλ is subject to saturation. Deriving stator and rotor linkage fluxes, in (16)-(18), lead to the time derivative of the magnetizing flux mλ . , , , s r s r m s r d di d l dt dt dt λ λ = + (20) Therefore, md dt λ has to be described by means of winding currents: ( ) ( ) i im m m m d d e d d e dt dt dt dt α αλ λ λ αλ= = + (21) where α is the angular position of mλ with respect to the d- axis. It also characterizes the position of the magnetizing current mi in the air gap, since the hysteresis angle is neglected. Writing mi , α and md dt λ as: 2 2 m dm qmi i i= + (22) 1 tan qm dm i i α −= (23) m m m m mdy m d d di di L dt di dt dt λ λ = = (24) we get: Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3739-3743 3742 www.etasr.com Khlifi: Explaining the d-q Magnetic Couplings Theoretically in Saturated Smooth Air-Gap AC Machines qmdm dm d dq did di L L dt dt dt λ = + (25) qm qmdm dq q d didi L L dt dt dt λ = + (26) m mdy m d L di λ = is called dynamic inductance by analogy with m m m L i λ = static inductance. Both coefficients are given by the conventional nonlinear magnetization curve, determined experimentally. ( )cos sindq mdy mL L L α α= − (27) cotd m dqL L L α= + (28) tanq m dqL L L α= + (29) B. Discussion Between rotor and stator we choose an integration path as shown in Figure 5. This path is defined by the angle θ and we can notice that due to symmetry, the flux density at the two air- gap segments in the path is the same. If we assume the permeability of iron to be infinite, then 0 iron H = .This means that for a given current i in the coil, the flux density in the air- gap varies sinusoidally with angle. Fig. 5. Representation of salient pole synchronous machines. Fig. 6. A sinusoidal winding on the stator. If the induction machine is saturated mdy m L L≺ and the direct axis of the critical two phase machine coincides with the axis with maximum air gap width, that is the quadrature direction in conventional salient pole machines. ( )cos sin ( )sin2 2 mdy m dq mdy m L L L L L α α α − = − = (30) where α is the angle position of the rotating magnetizing flux or current. cot cos = ( ) sin 2 2 sin 1 cos2 ( )( ) 2 ( ) ( )cos2 2 2 d m dq mdy m m m mdy m mdy m mdy m L L L L L L L L L L L L L α α α α α α = + − + + = + − + − = − (31) tan sin ( ) sin 2 2 cos 1 cos2 ( )( ) 2 ( ) ( )cos2 2 2 q m dq mdy m m m mdy m mdy m mdy m L L L L L L L L L L L L L α α α α α α = + − = + − = + − + − = + (32) For the saturated asynchronous machine static inductance is higher than that dynamic. While posing max( )mL L= and min( )mdyL L= , these inductances can be comparable with those of a diphase salient pole machine whose inductances along the direct axis, have weak air-gap and reluctance and inductances along the axis in squaring have strong reluctance and maximum of air-gap. For the saturated machine, permeance also depends in the state of saturation. The effect of saturation in the parts of iron results in the introduction of an equivalent air-gap with the same reluctance as the saturated parts. In the asynchronous machine, the function of permeance due to saturation varies along the armature. In fact, a saturated asynchronous machine can be regarded as equivalent to a diphase unsaturated salient pole synchronous machine. We presented a concise method of determining transient phenomena in a capacitor-excited induction machine, which, of course, requires taking into account the main field saturation. The new feature in this paper is the incorporation of cross- magnetization effects in the analysis. By this, we added a useful application to a number of previous papers dealing with smooth-air-gap machines. The current applied model describes the induction machine quantities in d-q-components. An essential part of the model is the inductance matrix along with the definitions of its elements. It is true that recently some contradictory views on the phenomenon of cross-magnetization have been stated. Therefore it might be useful to repeat the basic considerations under which the present model holds: The magnetizing current produces a sinusoidal mmf distribution which gives rise to a flux-density curve containing a fundamental and saturation harmonics. The winding distribution is such that only the fundamental flux induces a voltage. Expressing the magnetizing current and the Engineering, Technology & Applied Science Research Vol. 9, No. 1, 2019, 3739-3743 3743 www.etasr.com Khlifi: Explaining the d-q Magnetic Couplings Theoretically in Saturated Smooth Air-Gap AC Machines fundamental flux linkage as space-phasors, their magnitudes are related by the saturation curve. Note that there is a formal resemblance to the main inductances of a salient pole synchronous machine in stator frame representation, where the angular position of the pole axis is analog to the position α of the above mentioned magnetizing current space-phasor. In further consequence of this analogy, the (saturation dependent) (Lmdy−Lm) takes the place of the difference between direct and quadrature axis synchronous inductances. It has been shown that when the currents are used as state- variables, it is possible to show that the equations and physical picture in a saturated smooth-air-gap machine are similar to the equations and physical picture in a non-saturated salient-pole machine, where mutual coupling exists between the direct and quadrature axes and the self-inductances are different on the non-salient part of the machine due to the variation of the reluctance with the rotor angle. However, in contrast to this, in the saturated smooth-air-gap machine, the cross-saturation coupling and self-inductances vary with the angle of the magnetizing current. This is one of the reasons why a specific form of the inductance terms has been presented in [l, 2]. This physical fact is also evidenced by (30)-(33), in contrast to the forms before where the small change of the magnetizing inductance is expressed by utilizing the difference between dynamic and static inductances. We utilized the first derivative of the magnetizing inductance with respect to the magnetizing current. The same derivative has been utilized in the paper under discussion, where an alternative formulation of these inductances has been given. If it is the purpose to focus attention to the analogy described above, or to give a clear analytical representation of the modifications in the impedance matrix, then it is useful to have expressions containing the angle of the magnetizing current space vector. V. CONCLUSION A linear analysis of the machine based on the coupled circuit theory and assuming sinusoidal distribution of both mmf and flux density is performed and the saturated model is derived by making some model parameters dependent on suitable saturation laws. In this paper, a comprehensive winding model of both induction motors and synchronous generators has been proposed where the effects of saturation are incorporated in the main flux path and sometimes attributed to the leakage reactance too. Saturation of the magnetizing reactance must be driven by the main flux, whereas saturation of leakages is mainly dependent on the involved currents. In round rotor machines the influence of saturation only depends on the amplitude of the main flux, whereas in salient pole machines the flux angle may also be involved by considering, if known, different saturation characteristics as the flux moves from the direct to the quadrature axis. At steady state operation, only the main flux path is affected by saturation, whereas the leakages are usually maintained at their unsaturated values. According to the selected model, the saturated magnetizing reactance is defined as the slope of the chord corresponding to the operating point on the proper magnetizing characteristic. The effect of saturation at steady state is of main importance for a correct setting of the state variables before a transient analysis is started. For instance, the initial load angle of synchronous generators is closely related to the stability characteristic of the machines. APPENDIX: LIST OF SYMBOLS N Number of spire of magnetic circuit I Current in magnetic circuit ef Length of the air-gap Rf Magnetic Reluctance λf Flux iron Sf Surface of iron Hf Magnetic excitation lf Length of the field lines B Magnetic induction µr Relative permeability µ0 Permeability of vacuum d, q Indices for direct and quadrature axis i Space vector winding current REFERENCES [1] M. A. Khlifi, H. Rehaoulia, “General modeling of saturated AC machines for industrial drives”, COMPEL: The International Journal for Computation and Mathematics in Electrical and Electronic Engineering, Vol. 35, No. 1, pp. 44-63, 2016 [2] M. B. Slimene, M. A. Khlifi, M. B. Fredj, H. Rehaoulia, “Analysis of saturated self-excited dual stator induction generator for wind energy generation”, Journal of Circuits, Systems and Computers, Vol. 24, No. 9, p. 1550129, 2015 [3] M. B. Slimene, M. A. Khlifi, M. B. Ben, H. Rehaoulia, “Modeling of dual stator induction generator with and without magnetic saturation”, Journal of Magnetics, Vol. 22, No. 5, pp. 145-153, 2015 [4] M. A. Khlifi, M. B. Slimene, M. B. Fredj, H. Rehaoulia, “Performance evaluation of self-excited DSIG as a stand-alone distributed energy resources”, Electrical Engineering, Vol. 98, No. 2, pp. 159-167, 2016 [5] M. B. Slimene, M. A. Khlifi, M. B. Fredj, “Sensorless Speed Control for Dual Stator Induction Motor Drive Using IFOC Strategy with Magnetic Saturation”, Applied Computational Electromagnetics Society Journal, Vol. 32, No. 3, pp. 262-267, 2017 [6] I. lglesias, L. Garcia-Tabares, J. Tamarit, “A D-Q model for the self- commutated synchronous machine considering the effects of magnetic saturation”, IEEE Transactions on Energy Conversion, Vol. 7, No. 4, pp. 768-776, 1992 [7] M. A. Khlifi, H. Rehaoulia “Modelling of salient pole synchronous machine with magnetic saturation” International Journal of Physical Sciences, Vol. 6, no. 21, pp. 4928-4932, 2011 [8] M. A. Khlifi, “Influence of Cross-Saturation on the Various Models of Induction Machine”, International Journal of Automation and Control Engineering. Vol. 4, No. 4, pp. 268-272, 2016 [9] M. H. Shin, D. S. Hyun, S. B. Cho, S. Y. Choe, “An improved stator flux estimation for speed sensorless stator flux orientation control of induction motors”, IEEE Transactions on Power Electronics, Vol. 15, No. 2, pp. 312-318, 2000 [10] L. Harnefors, “Design and analysis of general rotor-flux-oriented vector control systems”, IEEE Transactions on Industrial electronics, Vol. 48, No. 2, pp. 383-390, 2001 [11] E. Levi, “A unified approach to main flux saturation modeling in D-Q axis models of induction machines”, IEEE Transactions on Energy Conversion, Vol. 10, No. 3, pp. 455-461, 1995 [12] M. B. Slimene, M. A. Khlifi, “Performance Limits of Three-Phase Self- Excited Induction Generator (SEIG) as a Stand Alone DER”, Journal of Electrical Engineering and Technology, Vol. 12, No. 1, pp. 145-150, 2017