Next Article in Journal
Assessment of the Gross, Technical and Economic Potential of Region’s Solar Energy for Photovoltaic Energetics
Next Article in Special Issue
A Hybrid Analytical Model for the Electromagnetic Analysis of Surface-Mounted Permanent-Magnet Machines Considering Stator Saturation
Previous Article in Journal
Influence of Phosphorus Speciation on Its Chemical Removal from Reject Water from Dewatering of Municipal Sewage Sludge
Previous Article in Special Issue
Design for Reducing Bearing Force Ripple and Torque Ripple of Integrated Magnetic Bearing Motor through Halbach Array
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Analytical Equivalent Circuit for Single-Sided Linear Induction Motors Considering Secondary Leakage Reactance

1
Department of Energy and Automation, Polytechnic School of University of São Paulo, São Paulo 05508-010, Brazil
2
Mauá Institute of Technology, São Caetano do Sul 09580-900, Brazil
*
Author to whom correspondence should be addressed.
Energies 2023, 16(3), 1261; https://doi.org/10.3390/en16031261
Submission received: 14 December 2022 / Revised: 6 January 2023 / Accepted: 19 January 2023 / Published: 24 January 2023

Abstract

:
Single-sided linear induction motors usually appear in magnetic levitation systems of transportation. Since the beginning of such developments, edge effects represent one of the great challenges to overcome in analytical modelling. For almost four decades, in order to simplify the mathematical treatment of border effects, most analytical models have not considered the secondary leakage flux properly. Although concise and accurate in most cases, such approaches have deficiencies in slotted secondaries. This paper presents an analytical equivalent circuit that considers the secondary reactance for both edge effects, i.e., entry and exit sides. The proposed approach uses an analogical RLC circuit which describes the behavior of magnetizing (exit) and demagnetizing (entry) waves, as well as adapted correction factors for transverse effects. By means of an 8 pole/120–175 N prototype, the measured thrust and vertical forces remarkably validate the model for frequencies equal or higher than 60 Hz. The relevance of secondary reactance in such cases can be explained by accounting its influence, specially, in the rise of the demagnetizing entry wave.

1. Introduction and Literature Review

Linear induction motors (LIM) has been object of study for over a hundred years [1], having its applications mainly in machine tools and transportation. The single-sided linear induction motor, a particular case of LIM, appear specially in magnetic levitation systems [2].
Since the well-known publication on LIM’s [3], intense research has followed on normal and traction forces in this type of machine [4]. The external characteristic of the linear induction motor strongly depends on the behavior of the electromagnetic travelling wave in its air gap, which has also moved scientific research to focus efforts on the so-called edge effects, resulting from the inevitable open magnetic circuit of that configuration [5].
Later, the University of Tokyo advanced in the modeling of linear induction motors with double primary, focusing on equating the edge effects [6]. The focus was not to develop a representative circuit, but to propose more accurate expressions for pulling force, output power, line current, efficiency and power factor, taking into account edge effects. In addition, the aforementioned work proposed a compensating winding located at the entry edge, with the aim of mitigating the undesirable effects, specially in the cases of high-speed applications, where such effects increase appreciably.
Until then, the approaches were established in a two-dimensional way. In the late 1970s, quasi-three-dimensional mathematical models began to be used (a two-dimensional analysis with compensating corrections for the third dimension). The work carried out with such a treatment focused on a linear motor with double primary and massive secondary [7], but also applicable to cases of single primary, considering a solid secondary with an magnetic yoke (back iron). The transverse and longitudinal edge effects were considered by applying Ampère’s Circuit Law in two directions, as well as the effects resulting from the stepped winding, that is, with partially filled slots at the borders of the primary, considering two displaced current blades. A simple equivalent circuit was presented, but only for consideration of the secondary power and induced voltages (in this case considered balanced, even recognizing the unbalance resulting from longitudinal edge effects).
At the same time, equivalent circuit models sought to incorporate the effects of open magnetic circuits through correction factors [8,9]. Since then, some advances were made in dynamic parameter models [10], optimization of the primary winding aiming a better traction performance [11] and equivalent circuits, including effects of partially filled primary slots [12].
One of the most important publications for the present paper made a direct analogy with the equivalent circuit of the rotary induction motor [13]. The focus was on corrections factors, resulting from edge effects in relation to the magnetizing branch. The cited paper had considerable weight in later publications and is still used today.
However, all the works cited addressed the issue considering a secondary without slots, either formed entirely of non-ferrous material or made of non-ferrous sheet, generally aluminum, with a back iron. Few publications considered a slotted secondary in the investigation, some in the analysis of normal and thrust forces in the aforementioned configuration [14] and others in the relevance of the contact resistance existing in the welded cage [15]. Particularly, in a publication on design optimization considering a slotted secondary, eventually called just a “ladder”, considered a circuit for massive secondaries [16]. See Figure 1 for different secondaries concepts.
Throughout the present work, the expression “slotted secondary” will refer to the traditional slotted structure in the form of a cage, and not to the secondary called “ladder-slit”, which consists of a grooved conductive plate, which gained attention from the 1990s onwards [17,18]. With a grooved aluminum plate, it was noticed that there was an increase in the resulting tensile strength and a decrease in material cost. Recently, some attention has been given to optimization in the design of such a configuration [19].
Other more recent efforts [20,21] focused their attention on solid sheet type secondaries, eventually with lateral asymmetry in the secondary [22], proposing tests to obtain parameters [23] or an equivalent circuit for solid secondary width optimization. Yet some works have used the conventional Duncan’s approach with modifications, aiming design refinements [24,25,26,27,28]. Some of the advances obtained in the above mentioned works were considered here. Once again, however, it should be noted that such publications did not consider the effect of secondary leakage, when similar in magnitude to that of the primary, on the demagnetizing and magnetizing waves at the entry and exit edges, respectively.
As far as we are aware, there is no equivalent circuit proposal that considers the secondary leakage effect on both edge effects (entry and exit). A fine job has been initiated in [12], but it did not consider the exit-magnetizing component or the vertical forces. Without intending to completely fill the gap, the present work proposes to be a support to the designer that, having at hand the pre-construction data of the linear motor, could determine the parameters of the equivalent circuit and, with basic computational resources, reach satisfactory results regarding the performance of traction and vertical forces. In addition, one would be able to compare such external characteristics, among others, in relation to the secondary configuration alternatives.
Taking into account that the slotted configuration appears in the prominent MAGLEV COBRA project of the Federal University of Rio de Janeiro [29,30], such an advance in modelling is considered relevant.
Finally, it is also intended to establish a experimental validation of the new analytical model for slotted secondaries, by means of a prototype specially manufactured, where thrust and vertical forces can be measured. It will be seen, as well, that the proposed model can be useful even in solid secondaries, be the continuous or the ladder-slit one, when the frequency is high enough.

2. New Equivalent Circuit

In this section, the objective is to highlight the contribution of the present work. In comparison with the most present analytical model in the literature [13], the proposed one adds the effect of secondary flux dispersion, once neglected. In this way, the demagnetizing wave at the entry border will have its shape affected. In addition to that, the magnetizing wave at the exit border will be accounted for in the circuit, which was not the case before.
In solid (continuous or “ladder-slit”) secondaries, the secondary leakage flux is, in fact, substantially lower than the primary. In slotted secondaries, however, such values become similar. Therefore, the proposed model will have greater weight in the latter cases.
According to [13], the main difference between rotating and linear induction motors is that, in the latter, as the primary moves, the secondary continuously presents “new” conductive material at the entry edge, which will oppose the sudden increase in magnetic flux. Similarly, at the exit border, there will be an opposition in the secondary to the sudden decrease in flux. Such phenomena mainly affect the flux distribution in the air gap, which would justify a model with correction factors in the magnetizing branch, according to the source cited. It is worth mentioning that such considerations do not concern the secondary load component.

2.1. Longitudinal Edge Effects

Longitudinal edge effects at entry and exit borders depend on time constants involving secondary and magnetization inductances and secondary ohmic resistance. Such time constants define rise and fall times of the armature reaction at the edges. The spatial distribution of flux along the air gap will also depend on the displacement speed.
At time zero, the secondary reaction current at the entry edge, referred to the primary circuit, will tend to cancel out the primary magnetizing current, almost equaling it in magnitude. Its increase will be quick, since at the input interface the only inductance to delay it will be that of secondary leakage, since the primary and secondary circuits of the new material are not yet linked.
Then, the decay of the demagnetizing wave will no longer be so fast, since now the time constant contemplates mutual flux between the parts, in addition to the secondary leakage.
Figure 2 illustrates the behavior of the secondary magnetomotive force in this demagnetizing effect at a constant speed, also considering the slots fully filled at the ends.
At the right hand of Figure 2, one can see also the magnetizing component of secondary reaction, which takes place at the exit border.
In solid secondaries, the leakage inductance is very low, and the areas described by the time constant exponentials L 2 L 2 R 2 R 2 = T 2 are neglected without significant loss of precision, leaving only the time constant component L 2 + L m L 2 + L m R 2 R 2 = T 2 as significant. According to [13], although the value of the secondary reaction varies along the length of the motor, its reflection to the primary winding can be considered as an average current in phase opposition to the magnetization current. Neglecting the secondary leakage reactance will lead to the circuit of Figure 3. In addition to the demagnetizing effect, there will be an increase in parasitic losses due to the path taken by the secondary reaction currents. Please see [13] for more details.
In the slotted secondary, leakage inductance is significantly higher than in solid and ladder-slit models. The rise of the demagnetizing secondary component current that reacts to the magnetization at the input edge, as it depends on the secondary dispersion, should now be included in the calculation. Additionally, the magnetizing component that opposes demagnetization at the exit edge must be considered.
The edges of entry and exit waves shown in Figure 2 obviously do not correspond to reality, given their discontinuities. Rounded edges are more coherent—see Figure 4. Therefore, to analytically describe the behavior of the waves, the authors propose an analogous (not equivalent!) RLC circuit—Figure 5, whose transient response will correspond to the demagnetizing and magnetizing waves. In the case of [13], L 2 0 for the purposes above mentioned, even though for losses and vertical forces there is partial consideration. See below that the analogous capacitances are functions of the resistances and inductances.
See the secondary leakage influence comparison at Figure 4. It is easy to see that, as the secondary dispersion increases, the effect of the demagnetizing longitudinal edge may be greater. On the other hand, the magnetizing component at the leading edge is no longer negligible.
The adequate treatment can be given by the transient response of an analog RLC circuit as shown in the Figure 5 (please note that it is not a equivalent circuit!).
The solution to the problem comprises two second-order ordinary differential equations:
d 2 i 2 m e n t t d t 2 + R 2 L 2 d i 2 m e n t t d t + 1 L 2 C e n t i 2 m e n t t = 0
d 2 i 2 m e x i t t d t 2 + R 2 L 2 d i 2 m e x i t t d t + 1 L 2 C e x i t i 2 m e x i t t = 0
Regarding the entry border, the solution will unfold in 3 possible damping cases, which will be considered properly. In applications with relatively large air gaps, it is possible to have a magnetizing inductance smaller than 3 L 2 , what would configure a underdamping case. For the exit border, the case of underdamping is evident, regardless of the values, as will be seen.
These are two series RLC circuits, with different capacitances—one for the entry, the other for exit waves analogy. The analysis of damping for entry border results:
R 2 2 L 2 2 > 1 L 2 C e n t = R 2 2 L 2 2 + L 2 L m , if L m > 3 L 2
R 2 2 L 2 2 < 1 L 2 C e n t = R 2 2 L 2 2 + L 2 L m , if L m < 3 L 2
R 2 2 L 2 2 = 1 L 2 C e n t = R 2 2 L 2 2 + L 2 L m , if L m = 3 L 2
Equations (3)–(5) refer to over (OD), under (UD) and critically damped (CD) situations, respectively. For the exit border, UD is always the case, since:
R 2 2 L 2 2 < 1 L 2 C e x i t = R 2 2 L 2 2
Now, applying the general solution for the entry OD, UD and CD transients, considering initial conditions i 2 m e n t 0 = 0 and v C e n t 0 = 0 , that is, d i 2 m 0 d t = R 2 I m L 2 , one can have, respectively,
i 2 m e n t O D t = I m k sinh k T 2 t e t 2 T 2 , k = 0 , 25 T 2 T 2 T 2 T 2 T 2 = L 2 R 2 ;   T 2 = L 2 + L m R 2
i 2 m e n t U D t = I m k sin k T 2 t e t 2 T 2 , k = T 2 T 2 T 2 0 , 25 T 2 0 , 25
i 2 m e n t C D t = I m T 2 t e t 2 T 2
In the above solutions, it is easy to see that in the critical damping interface, that is, L m = 3 L 2 and k = k = 0 , one have lim k 0 I m k sinh k T 2 t e t 2 T 2 = lim k 0 I m k sin k T 2 t e t 2 T 2 = I m T 2 t e t 2 T 2 .
If L 2 0 , entry wave is always OD, T 2 0 and k 0 , 5 , and so
i 2 m e n t t = lim L 2 0 I m k sinh k T 2 t e t 2 T 2 = I m e t T 2
i 2 m e x i t t = lim T 2 0 I m T 2 t T v e t T v T 2 = 0 ,
which are the expressions obtained by [13].
Then, with respect to the exit transient, the initial conditions are i 2 m e x i t T v + 0 and v C e x i t T v = 0 , that is, d i 2 m e x i t 0 d t = R 2 I m L 2 .
Then, one can obtain
i 2 m e x i t t = I m T 2 t T v e t T v T 2
Therefore, it is also clear that for L 2 0 there is no consideration of the magnetizing effect on the exit border, as, in fact, is done in [13].
For the current I 2 m referred to the primary in phase opposition to the magnetization component I m , it is sufficient to subtract the exit component from the entry one. From here the absolute value will be considered, for future inclusion in the magnetizing branch of the new model (see Figure 6):
I 2 m = 1 T v 0 T v i 2 m e n t t d t lim τ T v τ i 2 m e x i t t d t
Considering the dimensionless primary length on normalized time scales Q = T v T 2 and Q = T v T 2 , the current I 2 m results, for entry OD:
I 2 m O D = I 2 m e n t O D I 2 m e x i t = I m 1 Q 1 Q e Q 2 sinh k Q 2 Q k + cosh k Q Q = I m k m
For entry UD:
I 2 m U D = I 2 m e n t U D I 2 m e x i t = I m 1 Q 1 Q e Q 2 sin k Q 2 Q k + cos k Q Q = I m k m
For entry CD:
I 2 m C D = I 2 m e n t C D I 2 m e x i t = I m 3 Q e Q 2 2 + 4 Q = I m k m
One can see the new magnetizing reactance in Figure 6, where the correction factor can be k m , k m or k m , according to the entry transient damping.
Next, eddy losses must be corrected due to the path taken by the secondary reaction currents. First, one have to consider entry and exit waves components:
R 2 1 T v 0 T v i 2 m e n t 2 t d t + lim τ T v τ i 2 m e x i t 2 t d t 2
Secondly, the component corresponding to the stored energy that dissipates at exit border must be taken into account, considering the value of magnetizing current at this edge:
1 2 L m + L 2 I m + i m 2 e n t T v 2 1 T v = = 1 2 L m + L 2 I m + i m 2 e n t T v 2 R 2 L m + L 2 Q = R 2 2 Q I m + i m 2 e n t T v 2
Again, for entry OD, summing up the cited effects:
R 2 I m 2 1 Q + 1 4 Q e Q 1 2 Q + sinh 2 k Q 4 k Q + + 4 e Q 2 sinh k Q 4 k Q sinh 2 k Q 4 k 2 Q = R 2 k p I m 2
For entry UD:
R 2 I m 2 1 Q + 1 4 Q e Q 1 2 Q + sin 2 k Q 4 k Q + + 4 e Q 2 sin k Q 4 k Q sin 2 k Q 4 k 2 Q = R 2 k p I m 2
For entry CD:
R 2 I m 2 1 2 Q + 9 4 Q e Q Q + 2   + + 2 Q + Q 2 Q 2 e Q 2 Q = R 2 k p I m 2
Similarly, in Figure 6 one can calculate the eddy losses with correction factors k p , k p or k p .

2.2. Tranverse Edge Effects

For solid secondaries with a width greater than the primary, there is the simplified Russell and Norsworthy factor for the ν t h harmonic, modified by [8]:
k R N ν = 1 tanh ν π z 2 2 τ p ν π z 2 2 τ p 1 + 1 + 1.3 t o v d d tanh ν π z 2 2 τ p tanh ν π h o v τ p
The main difference from the slotted secondary, as far as transverse effects are concerned, is that there is no lateral extension of a conductive plate, but of the cage bars.
Equation (20) is for solid secondaries. Considering the adaptations of Figure 7, it can be used in the case of slotted ones.

2.3. Half-Filled Slots

The end slots of the primary core, in double-layer windings, inevitably become half-filled, as shown in Figure 8. The distribution of primary magnetomotive forces is affected and, consequently, the magnetizing branch and secondary parameters referred to the primary as well [19].
The primary mmf distribution is approximately halved in the two sections where the slots are partially filled, so the number of slots and the equivalent length can be calculated as
Q 1 = Q 1 2 y 1 m 1 q 1 + y 1 m 1 q 1 = Q 1 y 1 m 1 q 1
D = 2 p y 1 τ p + y 1 τ p = 2 p τ p = Q 1 τ s 1

2.4. Vertical Forces

Applying the new approach to the vertical attractive forces, the RMS value of the effective magnetization current is calculated, keeping in mind the existing proportionality in the main portion in relation to the stored energy:
F v a t t = 3 2 g L m 1 T v 0 T v I m + i 2 m e n t t 2 d t + lim τ T v τ i 2 m e x i t 2 t d t 2
Again, for entry OD:
F v a t t O D = 3 2 g L m I m 2 1 + 1 4 Q 3 2 Q e Q 1 2 Q + sinh 2 k Q 4 k Q + cosh 2 k Q 1 8 k 2 Q e Q 2 sinh k Q k Q + 2 cosh k Q Q = 3 2 g k v L m I m 2
For entry UD:
F v a t t U D = 3 2 g L m I m 2 1 + 1 4 Q 3 2 Q e Q 1 2 Q + sin 2 k Q 4 k Q cos 2 k Q 1 8 k 2 Q e Q 2 sin k Q k Q + 2 cos k Q Q = 3 2 g k v L m I m 2
For entry CD:
F v a t r C D = 3 2 g L m I m 2 1 + 1 4 Q 6 Q + + 4 e Q 2 1 + 2 Q e Q Q + 2 + 2 Q = 3 2 g k v L m I m 2
As for the repulsion component, the expression according to [13] remains:
F v r e p = k r I 2 2 h 12
The length h 12 now corresponds to the distance between the centers of the primary and secondary conductors. The constant k r must be determined from static tests. The results for this constant are shown in the next section.

3. Experimental Tests and Discussion

The manufacturing of a linear motor for tests involves certain difficulties, mainly with regard to the measurement of traction and vertical forces. In the case of the present work, a round geometry was chosen, as will be shown next.

3.1. Prototype

In order to expand the possibilities of experimental data, a double secondary was built: slotted and solid. The short primary was made in such a way as to be able to move laterally, sometimes coupling with the slotted secondary, sometimes with the solid. Part of the design can be viewed in Figure 9. The solid secondary will be treated with more closely attention in a future work. For the purposes of this present article, the efforts are concentrated in the slloted one.
In addition, vertical axes on the sides, guided by lateral bearings, made it possible to measure vertical forces.
The machine was conceived and built in the facilities of Equacional Elétrica e Mecânica company, where the photographs found in Figure 10 and Figure 11 were taken.
The concentricity of the primary and secondary radii occurs with an air gap equal to 2 mm (5 mm in the case of the solid secondary, considering the thickness of the conductive sheet). Any value other than this results in different air gaps across the primary. Thus, mean values at the experimental points were considered.
In Table 1, one can see the rated and constructive data of the prototype. For the purposes of this paper, only the slotted secondary was considered and tested. The comparison with the solid one will be presented in a future work.
The parameters of the equivalent circuit, in a airgap range of 1.86–4.24 mm, are shown in Table 2. The values do not comprise heat and skin effect corrections—these were conveniently taken into account according to frequency and temperature measured in tests.
One can see that the air-gap range demands different treatments for the entry wave (OD, UD, and CD, eventually).

3.2. Curves Compared

The prototype was tested in order to compare different conditions of source frequency and air-gap length. The test bench is shown in Figure 10 and Figure 11.
The thrust force was measured through a torque transducer in the electrodynamic brake, considering its lever arm length and the secondary diameter. The vertical force, on the other hand, was measured by difference of mass in a calibrated scale. The prototype had two vertical shafts through linear bearings, in order to eliminate undesirable frictions.
In Figure 12, one can see the source frequency influence on thrust force, along with the experimental points. The other curves (“No Correction”, that is, the classic induction motor equivalent circuit; and “Duncan”) serve the purpose of analysis for the new model precision and contribution. Here, the air gaps in the prototype are according to entry OD cases.
In sequence, Figure 13 shows, with the same average gaps, the behavior of resulting vertical forces through frequency variation. The dots, once more, are experimental points.
For the determination of constant k r , the locked secondary tests were necessary. Frequencies from 30 to 150 Hz, with average gaps from 1.86 to 4.24 mm, resulted in values presented in Table 3. It has been observed that not just the air gap and the frequency, but the damping case (OD or UD) has influence on this factor.
The different values of average air gaps need some clarification, and this has to do with difficulties faced during experimental tests. First, values are average because the prototype is round, resulting in eccentric radii for most cases. Second, with the variation of measured vertical forces along the tests, changes in air gap have occurred, since the method of measurement involves vertical displacement. In the prototype design, the primary remains supported on the scale by means of vertical axes and linear bearings. In this way, the different vertical forces along the tests provide different bending in the longer pieces, changing a little bit the distance between primary and secondary. This is a weakness of the prototype, to be mitigated in further studies.
Another difficulty in the experiments was the temperature of primary and specially of secondary windings, because of its great effect on the thrust force shape. Some adapted ventilation system was used so that the measurements took place with approximately the same temperatures in such windings. In the extreme cases, corrections were made to a reference temperature. Be means of a thermographic camera, the secondary temperature could be measured.
In the same way of Figure 12 and Figure 13 (entry overdamped case), thrust and vertical forces are also presented for higher values of air gap in Figure 14 and Figure 15—this time, according to entry underdamped cases.

4. Conclusions

The proposed model have shown reasonable accuracy in thrust and vertical forces prevision. The proposed circuit proved to be useful in cases where the secondary leakage inductance is not negligible.
As the frequency increases, the time required to travel the effective length of the primary will decrease, accentuating the demagnetizing (entry) and magnetizing (exit) effects. However, at the same time, increasing the frequency will also accentuate the skin effect, so as to increase the secondary resistance and decrease the secondary leakage. Depending on the case, there will be a range of frequencies where the model will be much more accurate than [13].
The relevance of the secondary leakage reactance, in the authors calculations, was more substantially present in its influence through the rise of the demagnetizing entry wave.
Finally, the work paves the way for some future paths, such as the optimization of the prototype regarding the measurement of vertical force with constant air gap, better development of a method for determine the equivalent circuit from tests and application of the model to other types of linear induction motors, such as double primary and/or secondary without back iron.

Author Contributions

Conceptualization, D.R.G.; methodology, D.R.G.; validation, D.G and I.E.C.; formal analysis, D.R.G. and I.E.C.; investigation, D.R.G.; resources, D.R.G.; data curation, D.R.G.; writing—original draft preparation, D.R.G.; writing—review and editing, D.R.G. and I.E.C.; supervision, I.E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding by institutions with grant number.

Data Availability Statement

Not applicable.

Acknowledgments

The manufacture of the prototype, as well as space and machinery for testing were donated and/or provided by the company Equacional Elétrica e Mecânica.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Laithwaite, E.R.; Barwell, F.T. Application of Linear Induction Motors to High-speed Transport Systems. Proc. Inst. Electr. Eng. 1969, 116, 713–724. [Google Scholar] [CrossRef]
  2. Alscher, H.; Boldea, I.F.; Eastham, A.R.; Iguchi, M. Propelling Passengers Faster Than a Speeding Bullet: As Maglev Urban Railways Get Tracking, Lines in West Germany, Japan, and Canada Are Expected to Follow England’s Commercial Lead. IEE Spectr. 1984, 21, 57–64. [Google Scholar] [CrossRef]
  3. Laithwaite, E.R. Linear Induction Motors. Proc. IEE Part A Power Eng. 1957, 104, 461–470. [Google Scholar] [CrossRef]
  4. Nasar, S.A. Electromagnetic Fields and Forces in a Linear Induction Motor, Taking Into Account Edge Effects. Proc. Inst. Electr. Eng. 1969, 116, 605–608. [Google Scholar] [CrossRef]
  5. Laithwaite, E.R. Some Aspects of Electrical Machines With Open Magnetic Circuits. Proc. Inst. Electr. Eng. 1968, 115, 1275–1283. [Google Scholar] [CrossRef]
  6. Yamamura, S.; Ito, H.; Ishulawa, Y. Theories of the Linear, Induction Motor and Compensated Linear Induction Motor. IEEE Trans. Power Appar. Syst. 1972, PAS-91, 1700–1710. [Google Scholar] [CrossRef]
  7. Lee, C.H.; Chin, C.Y. A Theoretical Analysis of Linear Induction Motors. IEEE Trans. Power Appar. Syst. 1979, PAS-98, 679–688. [Google Scholar] [CrossRef]
  8. Gieras, J.F.; Dawson, G.E.; Eastham, A.R. A New Longitudinal End Effect Factor for Linear Induction Motors. IEEE Trans. Energy Convers. 1987, EC-2, 152–159. [Google Scholar] [CrossRef]
  9. Pai, R.M.; Boldea, I.; Nasar, S.A. A Complete Equivalent Circuit of a Linear Induction Motor With Sheet Secondary. IEEE Trans. Magn. 1988, 24, 639–654. [Google Scholar] [CrossRef]
  10. Zhang, Z.; Eastham, T.R.; Dawson, G.E. LIM Dynamic Performance Assessment from Parameter Identification. In Proceedings of the Conference Record of the 1993 IEEE Industry Applications Conference Twenty-Eighth IAS Annual Meeting, Toronto, ON, Canada, 2–8 October 1993. [Google Scholar]
  11. Hofmann, R.; Binder, A.; Pfeiffer, R. Investigations on a Linear Induction Machine for Railway Applications. In Proceedings of the IEEE International Electric Machines and Drives Conference, Cambridge, MA, USA, 17–20 June 2001. [Google Scholar]
  12. Woronowicz, K.; Safaee, A. A Novel Linear Induction Motor Equivalent-Circuit With Optimized End Effect Model. Can. J. Electr. Comput. Eng. 2014, 37, 34–41. [Google Scholar] [CrossRef]
  13. Duncan, J. Linear Induction Motor–Equivalent Circuit Model. IEE Proc. B–Electr. Power Appl. 1983, 130, 51–57. [Google Scholar] [CrossRef]
  14. Koseki, T.; Sone, S.; Yokoi, T.; Ebihara, D.; Osawa, S. Investigation on Secondary Slot Pitches of a Cage-type Linear Induction Motor. IEEE Trans. Magn. 1993, 29, 2944–2946. [Google Scholar] [CrossRef]
  15. Park, S.C.; Kim, B.T. Effect of Contact Resistance Between Side-bar and Secondary Conductors in a Linear Induction Motor With a Cage-type Secondary. IEEE Trans. Magn. 2003, 39, 1562–1565. [Google Scholar] [CrossRef]
  16. Ravanji, M.H.; Nasiri-Gheidari, Z. Design Optimization of a Ladder Secondary Single-Sided Linear Induction Motor for Improved Performance. IEEE Trans. Energy Convers. 2015, 30, 1595–1603. [Google Scholar] [CrossRef]
  17. Koseki, T.; Sone, S.; Yokoi, T.; Ebihara, D. Analysis of secondary eddy current path of a linear induction motor. In Proceedings of Linear Drive Symposium; LD-92-100; IEEE: Tokyo, Japan, 1992; pp. 21–29. [Google Scholar]
  18. Yamaguchi, T.; Ito, M.; Matsui, K. Improvement of thrust of linear induction motor using modified ladder slits. In Proceedings of the Power Conversion Conference—-PCC’97, Nagaoka, Japan, 6 August 1997; Volume 2, pp. 563–566. [Google Scholar]
  19. Lv, G.; Zhou, T.; Zeng, D. Design of Ladder-Slit Secondaries and Performance Improvement of Linear Induction Motors for Urban Rail Transit. IEEE Trans. Ind. Electron. 2018, 65, 1187–1195. [Google Scholar] [CrossRef]
  20. Xu, W.; Zhu, J.G.; Zhang, Y.; Li, Y.; Wang, Y.; Guo, Y. An Improved Equivalent Circuit Model of a Single-Sided Linear Induction Motor. IEEE Trans. Veh. Technol. 2010, 59, 2277–2289. [Google Scholar] [CrossRef] [Green Version]
  21. Zare-Bazghaleh, A.; Meshkatoddini, M.Z.; Fallah-Choolabi, E. Derivation of Equivalent Circuit Parameters for Single-Sided Linear Induction Motors. Plasma Sci. IEEE Trans. 2015, 44, 3637–3644. [Google Scholar] [CrossRef]
  22. Lv, G.; Zeng, D.; Zhou, T.; Degano, M. A Complete Equivalent Circuit for Linear Induction Motors With Laterally Asymmetric Secondary for Urban Railway Transit. IEEE Trans. Energy Convers. 2021, 36, 1014–1022. [Google Scholar] [CrossRef]
  23. Hirahara, H.; Inoue, M.; Yamamoto, S. A Method for Determining Equivalent Circuit Constant of Linear Induction Motors Using Locked Mover and Standstill Impedance Tests. In Proceedings of the IEEE International Power Electronics Conference, Himeji, Japan, 15–19 May 2022. [Google Scholar]
  24. Heidari, H.; RassÕLKIN, A.; Razzaghi, A.; Vaimann, T.; Kallaste, A.; Andriushchenko, E.; Belahcen, A.; Lukichev, D.V. A Modified Dynamic Model of Single-Sided Linear Induction Motors Considering Longitudinal and Transversal Effects. Electronics 2021, 10, 933. [Google Scholar] [CrossRef]
  25. Lv, G.; Zeng, D.; Zhou, T. An Equivalent Circuit of the Single-Sided Linear Induction Motor Considering the Discontinuous Secondary. IET Electr. Power Appl. 2019, 13, 31–37. [Google Scholar] [CrossRef]
  26. Teriaiev, V.; Dovbyk, A.; Kornienko, V.; Pechenik, M.; Buryan, S. Generalized Mathematical Model of a Linear Induction Motor. In Proceedings of the 2022 IEEE 41st International Conference on Electronics and Nanotechnology (ELNANO), Kyiv, Ukraine, 10–14 October 2022; pp. 741–745. [Google Scholar] [CrossRef]
  27. Di, J.; Fletcher, J.E.; Fan, Y.; Liu, Y.; Sun, Z. Design and Performance Investigation of the Double-Sided Linear Induction Motor With a Ladder-Slot Secondary. IEEE Trans. Energy Convers. 2019, 34, 1603–1612. [Google Scholar] [CrossRef]
  28. Zhang, L.; Obeid, H.; Laghrouche, S.; Cirrincione, M. Second Order Sliding Mode Observer of Linear Induction Motor. IET Electr. Power Appl. 2019, 13, 38–47. [Google Scholar] [CrossRef]
  29. Sotelo, G.G.; de Andrade, R., Jr.; Dias, D.H.N.; Ferreira, A.C.; Costa, F.; Machado, O.J.; de Oliveira, R.A.H.; Santos, M.D.A.; Stephan, R.M. Tests with One Module of the Brazilian Maglev-Cobra Vehicle. IEEE Trans. Appl. Supercond. 2013, 23, 3601204. [Google Scholar] [CrossRef]
  30. Mattos, L.S.; Rodriguez, E.; Costa, F.; Sotelo, G.G.; de Andrade, R., Jr.; Stephan, R.M. MagLev-Cobra Operational Tests. IEEE Trans. Appl. Supercond. 2016, 26, 3600704. [Google Scholar] [CrossRef]
Figure 1. Difference between solid, ladder, and slotted secondaries.
Figure 1. Difference between solid, ladder, and slotted secondaries.
Energies 16 01261 g001
Figure 2. Secondary mmf reaction behavior through time.
Figure 2. Secondary mmf reaction behavior through time.
Energies 16 01261 g002
Figure 3. Resulting equivalent circuit according to [13].
Figure 3. Resulting equivalent circuit according to [13].
Energies 16 01261 g003
Figure 4. Secondary mmf waves compared, as secondary leakage reactance increases.
Figure 4. Secondary mmf waves compared, as secondary leakage reactance increases.
Energies 16 01261 g004
Figure 5. Analogous circuit to describe secondary mmf reaction waves.
Figure 5. Analogous circuit to describe secondary mmf reaction waves.
Energies 16 01261 g005
Figure 6. New proposed equivalent circuit considering all secondary leakage effects.
Figure 6. New proposed equivalent circuit considering all secondary leakage effects.
Energies 16 01261 g006
Figure 7. Section view of linear motor with slotted secondary.
Figure 7. Section view of linear motor with slotted secondary.
Energies 16 01261 g007
Figure 8. Slot number and lengths with half-filled slots, according to [19].
Figure 8. Slot number and lengths with half-filled slots, according to [19].
Energies 16 01261 g008
Figure 9. Prototype design.
Figure 9. Prototype design.
Energies 16 01261 g009
Figure 10. Prototype on test bench.
Figure 10. Prototype on test bench.
Energies 16 01261 g010
Figure 11. Test bench main items.
Figure 11. Test bench main items.
Energies 16 01261 g011
Figure 12. Comparison of thrust force through frequency variation (OD). Average gaps in ascending order [mm]: 1.86 / 2.15 / 2.21 / 2.32 / 2.35 .
Figure 12. Comparison of thrust force through frequency variation (OD). Average gaps in ascending order [mm]: 1.86 / 2.15 / 2.21 / 2.32 / 2.35 .
Energies 16 01261 g012
Figure 13. Comparison of vertical force through frequency variation (OD). Average gaps in ascending order [mm]: 1.86 / 2.15 / 2.21 / 2.32 / 2.35 .
Figure 13. Comparison of vertical force through frequency variation (OD). Average gaps in ascending order [mm]: 1.86 / 2.15 / 2.21 / 2.32 / 2.35 .
Energies 16 01261 g013
Figure 14. Comparison of thrust force through frequency variation (UD). Average gaps in ascending order [mm]: 3.89 / 3.99 / 4.11 / 4.23 / 4.24 .
Figure 14. Comparison of thrust force through frequency variation (UD). Average gaps in ascending order [mm]: 3.89 / 3.99 / 4.11 / 4.23 / 4.24 .
Energies 16 01261 g014
Figure 15. Comparison of vertical force through frequency variation (UD). Average gaps in ascending order [mm]: 3.89 / 3.99 / 4.11 / 4.23 / 4.24 .
Figure 15. Comparison of vertical force through frequency variation (UD). Average gaps in ascending order [mm]: 3.89 / 3.99 / 4.11 / 4.23 / 4.24 .
Energies 16 01261 g015
Table 1. Rated values of the prototype.
Table 1. Rated values of the prototype.
DataValueDataValue
v [m/s]2.6–13.5 P t h r [kW]0.3–2.4
F t h r [N]120–175 V 1 l i n e [V]87.5–437.5
f [Hz]30–150 2 p 8
m 1 3 N 1 256
l z 1 [mm]80 l z 2 [mm]90
Q 1 / Q 1 26 / 24 Q 2 / Q 2 38 / 35
Table 2. Calculated parameters for the prototype.
Table 2. Calculated parameters for the prototype.
ParameterValueParameterValue
l g [mm]1.86–4.24 k R N 0.844
R 1 [ Ω /ph] 0.646 R 2 [ Ω /ph] 0.524
L 1 [mH/ph]7.58–7.36 L 2 [mH/ph]4.31–4.20
L m [mH/ph]18.1–8.54 L m / L 2 4.2–2.0
Table 3. Experimental k r values for various frequencies and average air gaps.
Table 3. Experimental k r values for various frequencies and average air gaps.
mm/Hz 1.86 / 30 2.15 / 60 2.21 / 90 2.32 / 120 2.35 / 150
k r 0.00637 0.00487 0.00337 0.00330 0.00330
mm/Hz 3.89 / 30 3.99 / 60 4.11 / 90 4.23 / 120 4.24 / 150
k r 0.0160 0.0155 0.0155 0.0120 0.0130
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gomes, D.R.; Chabu, I.E. A Novel Analytical Equivalent Circuit for Single-Sided Linear Induction Motors Considering Secondary Leakage Reactance. Energies 2023, 16, 1261. https://doi.org/10.3390/en16031261

AMA Style

Gomes DR, Chabu IE. A Novel Analytical Equivalent Circuit for Single-Sided Linear Induction Motors Considering Secondary Leakage Reactance. Energies. 2023; 16(3):1261. https://doi.org/10.3390/en16031261

Chicago/Turabian Style

Gomes, Daniel R., and Ivan E. Chabu. 2023. "A Novel Analytical Equivalent Circuit for Single-Sided Linear Induction Motors Considering Secondary Leakage Reactance" Energies 16, no. 3: 1261. https://doi.org/10.3390/en16031261

APA Style

Gomes, D. R., & Chabu, I. E. (2023). A Novel Analytical Equivalent Circuit for Single-Sided Linear Induction Motors Considering Secondary Leakage Reactance. Energies, 16(3), 1261. https://doi.org/10.3390/en16031261

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop