Next Article in Journal
An ALO Optimized Adaline Based Controller for an Isolated Wind Power Harnessing Unit
Previous Article in Journal
Neonatal Jaundice Detection Using a Computer Vision System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

End Effects and Geometric Compensation in Linear Permanent Magnet Synchronous Generators with Different Topologies †

Department of Electrical Engineering, Division of Electricity, Uppsala University, Box 65, 75103 Uppsala, Sweden
*
Author to whom correspondence should be addressed.
This paper is an extended version of our paper published in 2020 24th International Conference on Electrical Machines (ICEM), Gothenburg, Sweden, 23–26 August 2020.
Submission received: 15 September 2021 / Revised: 27 September 2021 / Accepted: 28 September 2021 / Published: 9 October 2021
(This article belongs to the Section Electrical Engineering Design)

Abstract

:
Electricity production from ocean waves with different solutions is a topic of major research interest. Many of such designs are based on linear generators that inherently introduce end forces. In this paper, detent force using Maxwell Stress Tensor and induced voltage is initially investigated for two different winding patterns for a generator topology with buried magnets in a finite element software. Two ways of overcoming the end forces are further examined: the first method reduces the magnetic flux variations of the translator between stator and air. The second method aims at countering the end forces at both ends for full active stator area. A comparison is then made between buried and surface-mounted topologies for the second end effect compensation method. Both no-load and load conditions are investigated in the comparison. The end effect compensation shows promising results for both topologies. Some clear similarities of the extended stator used to counter the end forces are also apparent, where the stator extensions completely cover the outer poles of both topologies. The results also indicate a longer full active stator area for the buried topology for the same pole-pitch and stroke length, resulting in a higher average voltage for partial stator overlap.

1. Introduction

Linear electrical machines, though much less common than rotating electrical machines, can be utilized in several applications such as wave-powered electricity generation. The basic working principles for linear electrical machines and their rotating counterpart are in several aspects the same. However, in a linear electrical machine, there are open ends that cause an uneven flux distribution both in regards to the transition between air and stator iron and, as for a buried topology, the reduced flux for the outer poles [1]. The open ends also mean that the overlap between the translator, the moving part of the linear electrical machine, and the stator can vary. Further, the open ends lead to some special traits of the linear electrical machines, for example, that the stator experiences a varied number of poles for a partial stator overlap. As shown in [2], the flux distribution varies depending on if an even or odd number of poles are used, which means that the flux distribution will change when there is partial stator overlap.
Similar to the rotating permanent magnet (PM) machine topologies, there is a broad range in PM placement and utilization in linear electrical machines as well. The linear electrical PM machines can be constructed with surface-mounted PMs [3,4,5], buried PMs (also called sandwich-mounted PMs) [6], or a Halbach array [7]. The PMs are often placed on the moving translator but can also be placed on the stator to utilize the magnetic gearing effect in a linear Vernier hybrid machine for both surface-mounted [8] and buried (in this case V-shaped) [9] PMs.
When using buried PMs in a linear generator, the open ends for the outer poles of the translator result in different end forces compared to surface-mounted PMs. The flux concentration topology for buried PMs also results in different end forces in regards to the stator ends. The end forces combined with the cogging forces result in detent forces. The detent force can be approximated as a sum of sine-functions with periodicities depending on the different force components [4,10]. In the generator presented in [11], the three phases were separated into different stator packages. An advantage of this is that the detent forces can be reduced as the force peaks occur at different times. To further protect the brittle surfaces of the ferrite magnets, rubber films were placed between the magnets and pole shoes, which also could dampen forces in the longitudinal direction [12,13].
End forces of linear machines and ways of reducing these forces have been investigated in earlier studies. One study approximates the flux distribution to a reluctance model through a magnetic field similarity method and investigates end-teeth variations to counter the end forces under no-load, derived through virtual work [3]. One study separates the stator block into sections with two phases in each block to counter the end forces [14]. Another study adds assistant teeth and notes the different magnitudes in phase voltages for windings closer to the stator ends [15]. For a linear Vernier hybrid machine, the force ripple and cogging have been minimized by optimizing the consequent poles, translator poles, and PM width [16]. The study most similar to this paper studied surface-mounted PMs where the length of the end teeth are changed to reduce the end forces, where the end forces are calculated using inductance formulations derived from a frozen permeability technique [17]. The study by Yao et al. [17], however, has a greater focus on the skewing effects of the linear machine. Most end force optimization studies on linear electrical machines are done for surface-mounted PMs [3,14,16,17].
The translator can, during operation, fully or partially overlap the stator. The partial overlap of the stator can be significant, unless the relative length of the stator is short compared to the translator. The stator can also be made much longer than the translator to have a constant partial overlap during operation [18]. Partial overlap has previously been studied both experimentally [4,6,19] and through simulations [4,5], though there is still more to learn. As for the studies on end-force minimization, most of the studies for partial overlap are also done for surface-mounted PMs [4,5,18,19]. According to [6], the induced voltage should decrease linearly with decreasing stator overlap and thereby the power generation should decrease quadratically. In [4], it was observed that the cogging force is larger when the stator is fully overlapped.
The longitudinal ends of a linear electrical machine influence the flux distribution significantly [2], can cause non-negligible end forces [14], and can have a dominating effect on the thrust ripple [17], making it a topic worth investigating further.
In this paper, the detent forces and induced voltages are initially investigated in a linear synchronous generator with buried PMs through finite element simulations. Alterations of the stator ends to reduce the detent forces are also proposed and investigated through simulations. A preliminary force comparison between Maxwell Stress Tensor (MST) and Virtual Work (VW) is done to show that the results of the two methods match. The linear synchronous generator with PMs (LPMSG) has a stator that is shorter than the translator and can during operation have partial stator overlap. The influence of partial overlap can thus be studied for detent forces, attraction force, and induced voltage. The buried topology is further compared to a topology with surface-mounted PMs. Currents in the stator windings are also introduced to see if similar behavior in end force reduction can be achieved under load. The magnitude of the power associated with the electric current and the detent forces can ultimately be compared to further verify the validity of the results.

2. Method

2.1. Force Calculations

The electromagnetic forces in an electrical machine can be derived in several ways, whereof one common way is to use the principle of Virtual Work (VW), which is used by, e.g., [3,5,7]. When utilizing the principle of VW, the force in the y-direction is calculated as the gradient of the magnetic co-energy: F = W ϕ y . The magnetic co-energy is expressed as
W ϕ = V 0 H B d H d V .
For a linear case, such as air, the inner integration can be simplified to 0 H B d H = μ 2 | H | 2 [20]. The force calculation using virtual work is more often used in rotating machines when the relevant air is enclosed between the rotor and stator. For linear generators, the translator does not enclose the stator, or vice versa, leading to non-negligible flux passing between the stator and the translator on the upper and lower sides of the stator. The flux passing through the translator ends is greater for a linear generator with buried PMs than it is for surfance-mounted PMs.
Another way to calculate the force is the line integration of Maxwell’s stress tensor (MST) [21] (ch.1). This is used by, e.g., [22,23]. The tangential magnetic component of the stress tensor, or in this case the y-direction, is given by
σ y = 1 μ 0 B y B x .
The total tangential force for a rotating electrical machine is expressed as the line integral enclosing the rotor and further averaged over the air gap, resulting in Arkkio’s method [21,24] (p. 56). For a linear electrical machine, the line integral is instead a line enclosing the translator. Ergo, the line integral runs through the air gap and splits the translator from the stator. The line on the opposite side of the translator is arbitrary. Having it at the border of the translator gives, as one might expect, negligible force. The line integration of the tangential stress tensor is, however, sensitive to the mesh [1]. Therefore, the same principle as Arkkio’s method is used, where the line integration is averaged over the air gap by integrating the tangential stress tensor over the air gap surface and dividing by the length of the air gap. Rather than torque, the force F is expressed as
F = l μ 0 l a g S B x B y d S .
where l is the stack length of the machine, μ 0 is the permeability of the vacuum, l a g is the length of the air gap, S is the area of the air gap, and B is the magnetic flux density. The force F is calculated instead of torque, since there ideally should be no torque in a linear electrical machine.
When using Maxwell stress tensor, one does not have to look at all domains with a significant flux density: only a boundary that encloses the translator (or stator). This is illustrated in Figure 1 for the slot per pole and phase, q = 1 , which shows that the detent force from Maxwell stress tensor (MST) is equal to the total detent force from virtual work for all domains (VW All) for an incremental displacement of one pole-pitch τ p . The end position of the displacement is when the generator reaches full stator overlap at t = 0.6 s. It also shows two components of the virtual work outside of the air-gap domain that are not negligible: the domains above and below the stator (ABS) and the domains within the stator slots (VW S s l o t s ).
These two force contributions give a clear representation of the detent forces in the linear generator for q = 1 : cogging force of three times the frequency of the pole pitch and the end force with a period the same as the pole pitch. In the figure, the same two contributions are plotted for the starting position of 50% stator overlap for half a pole pitch. From this, one can see that the cogging force increases with stator overlap, while the magnitude of the end force remains the same with less than 100% stator overlap. As the Maxwell Stress Tensor is based on the values given at the closed surface around the translator, it is more difficult to determine the separate contributions of end effects and cogging effects in comparison to the use of virtual work.
Similar to the tangential forces, one can calculate the attraction forces (normal forces) from the stress tensor elements given in (4). The force is acquired from the surface integral of the same boundaries as in (3).
σ n = 1 2 μ 0 ( H x 2 H y 2 ) ,
where H x and H y are the normal and tangential components of the magnetic field strength, respectively.

2.2. Finite Element Simulations

The software used in the simulations is the finite element software COMSOL MultiphysicsTM (COMSOL Multiphysics® is a Trademark of COMSOL AB). In the simulations for the buried topology, a translator with pole shoes and buried permanent magnets of alternating polarity are used. The PMs are rectangular, while the pole shoes are slightly extruded towards the air gap and have a rounded surface towards the air gap. The PM sizes of the flux concentrating topology are dimensioned to reach an adequate flux density level in the stator. The top and bottom pole shoes were split into halves. An example of buried PMs can be seen in Figure 2, together with the different winding patterns and the coordinate system used in this paper.
The numbers of slots per pole and per phase q used in this paper are 1 and 6/5, respectively. The translator is kept the same, while the number of slots are set to match the two winding topologies and the slot-pitch is determined by q and τ p . The windings closest to the stator yoke is shifted by one for q = 1 , reducing the thirdrd harmonic, with a trade-off of an increased fifth harmonic [1]. The cables are rectangular and take up 33% of the slot-pitch with 1 mm space on either side of the cables, representing the non-magnetic insulation. In the simulation model, this space is filled with air. The length of the stator tooth is then expressed as the remainder of the slot-pitch.
The full active stator area covers 10 pole-pitches, but to allow for partial translator overlap with a maintained full stator overlap, 6 more PMs are added. From Figure 2, one can see that each pole-pitch overlaps exactly three stator slots for q = 1 and slightly more than three slots for q = 6 5 . As all PMs for q = 1 enter a new slot at the same time, the cogging force is much higher for integer q numbers. Translator and stator parameters used in the simulation for the buried topology can be seen in Table 1 and Table 2, respectively.
The formulation for the PMs are based on the relative permeability μ r and the remanent magnetic flux density B r and given as
B = μ 0 μ r H + B r .
Assuming a linear BH-curve in the second quadrant until the point of maximum energy product, the relative permeability μ r is given by
μ r = B r 2 4 μ 0 ( B H ) max
where ( B H ) max is the maximum energy product. In the simulations for the buried topology, ferrites of grade Y40 are used. Values for B r and ( B H ) max are taken from online (https://e-magnetsuk.com/ferrite_magnets/ferrite_grades.aspx, accessed on 3 December 2020).
The stator iron and translator pole shoes are of Silicon Steel NGO 50PN270, taken from COMSOL material library. The non-linear BH-curve can be seen in Figure 3. Induced eddy currents and hysteresis losses are not considered in the study. The translator and its adjacent moving boundary has a constant velocity set to 0.5 m/s. The detent force and induced voltage are calculated for a stack length of 1 m. It can thus directly correspond to the force and voltage per unit of stack length.
The simulations are conducted in a two-dimensional finite element environment with the formulation for the Magnetic vector potential. The translator motion is acquired by using a moving mesh. A continuous boundary between the two bodies enforces a continuous magnetic vector potential at the boundary. The air gap is separated into layers, which, together with a quadratic discretization, increases the accuracy of the model. The coil groups are assumed homogeneous, with a uniform current density for each phase. The mesh is finer in the translator and the stator. The mesh is extra-fine in the air gap, since the boundary between the moving domain and the stationary domain is located in the air gap, and many of the data presented in the results are retrieved from the air gap. The complete mesh consists of a total of 171,056 elements in the simulation files with buried topology and of totally 151,931 elements in the simulations files with surface-mounted topology. The configuration of the finite element environment is also detailed in [1].
A preliminary comparison between experimental results and the simulated linear generator of a buried topology can be seen in [1]. Experimental measurements show similar behaviour for partial stator overlap and the ripple associated with the detent forces.

2.3. End Effect Compensation

Two methods are introduced in this study to reduce the impact of end effects for the buried topology.
The first method is to add soft magnetic material on the outer ends of the generator to keep the poles fully magnetized the entire stroke length. The magnetic material would, of course, have to be composed of thin lamination to minimize the eddy currents similar to the stator steel. An advantage to this is that there is little change in magnetization of the PMs during the whole stroke length and thus less end forces. The PMs are also thereby less subjected to irreversible demagnetization that otherwise is the case when they are surrounded by air. The added pieces, here denoted end plates (EP), are of the same width as the stator yoke.
The second method is to adjust the end regions of the stator with extended steel to coincide with the peak force of the poles leaving and entering the stator. This is, in some aspects, similar to [17], but here it is done for buried PMs in a finite element environment and a different method for the force calculation. This method is, in the paper, denoted extended stator (ES). In [4], the detent force is defined as
F c ( x ) = a 0 + b 1 c o s ( 2 π 1 w p x + θ 1 ) + b 2 c o s ( 2 π 3 w p x + θ 2 ) ,
where a 0 and b 1 are constant amplitudes while b 2 represents the cogging forces that increase in amplitude with active stator area. The introduction of extended stator ends would allow the splitting of θ 1 in (7) into two components and separating those to 180 to cancel out the force peaks. This is only possible during the full stator overlap.
The two end-force compensation methods are illustrated in Figure 4.

2.4. Comparison of Force Reduction between Buried and Surface Topology under No-Load and Load Conditionz

To further investigate how the alterations of the stator ends affect the end forces, a comparison is made between the buried and surface topologies. The hypothesis is that there is a difference in behavior because of the forced flux path through the surface-mounted PMs to the fluctuating proximity of the stator teeth, whereas the buried topology allows for a higher flux concentration at pole face regions closer to the stator teeth.
The dimensions of the surface-mounted translator are set in order to reach similar rms values of the no-load voltage at a time interval within the full active stator area. The same time interval is also used to investigate the impact of introduced currents in the stator windings. In this study, the introduced currents are in phase with the induced voltage. In practice, this kind of phase angle compensation can be accomplished by constant-torque-angle control [25]. As similar behavior of the end force reduction was apparent for both winding patterns, the winding patterns used for the comparison between the two topologies is q = 6 5 . As a way of further evaluating the force computations, the electrical power without resistive losses can be compared to the power associated with the magnetic forces. By having the currents in-phase with the induced voltage E, the load angle ( δ ) equals the phase angle ( θ ). The terminal voltage, expressed as V = ( E R I ) c o s ( δ ) is then used to determine the electrical power
P el = V I c o s ( θ ) = E I R I 2 .
The translator and stator dimensions can be seen in Table 1 and Table 2, respectively. The fill factor for the stator slots are 51%.
For the buried topology, the relative permeability is derived from (6). For the NdFeB magnets used in the surface-mounted topology, the relative pearmeability is set similar to other literature [21]. The difference in topology between surface-mounted and buried results in different flux paths in the translator. Due to the difference in flux path, the initial displacement of the surface-mounted PMs have to be adjusted to meet the same voltage phases. In other words, by adjusting the displacement, one can get the same phasors for the quadrature voltage component. Since the full active stator area is defined by the adjoining and leaving outer poles of the buried topology, the surface-mounted topology has a shorter full active stator area equal to half a pole-pitch on either side of the translator. This is visualized in Figure 5 and later evident especially in the results for the attraction forces. Despite covering the stator at t = 0.54 s for surface-mounted PMs, it is not until t = 0.6 s that it reaches a distribution that, ideally, covers the whole stator. For the two topologies, the stator overlap is therefore the same when the translator enters the stator and different when leaving the stator. The figure also reveals the varying flux in the translator yoke for the surface-mounted topology. The shorter coverage is the cause for a shorter interval of full active stator.

3. Results and Discussion

3.1. Buried Topology without End Region Alterations

Without any alterations to the end regions, the detent force is shown in Figure 6.
From the figure, it is evident that the peaks of the detent forces appear at the same instants. The end forces arise when the poles of the translator enter and leave the stator. When a pole enters the stator, the flux increases remarkably. The peaks of the average B y component for the outer stator tooth face perfectly align with the peak detent forces. This is in accordance with (2), where the time instant of maximum B y coincides with the time when the detent forces reach their maximum. Here, the maximum peak force is when the boundary between pole and PM passes the end of the stator iron.
The induced voltages for each phase can be seen in Figure 7. The induced voltage appears to be proportional to the stator overlap, as predicted in [6]. The induced phase voltages show slightly different magnitudes for a finite stator length. This is because of the different contributions along the stator windings for every phase. A resembling imbalance between phases was mentioned in [15].
This difference in magnitude between phases is more prominent when the ratio of the machine length over stroke length is low, as it is for this machine. Normally, the machine is longer where only the outer windings of the stator are subjected to a reduced flux linkage. The induced voltage still appears proportional to the active stator area. The end forces will also be of less significance compared to the cogging forces for a longer linear machine. In this paper, the end forces do, however, appear significant, and looking at ways to minimize the end forces can be of interest even for longer linear machines with partial overlap.

3.2. Force Reduction for Buried Topology with Different Winding Patterns

The flux distribution of a generator with added soft magnetic material on the stator is shown in Figure 4a, where the flux through each PM is of similar magnitude throughout the translator. The detent forces for q = 1 and q = 6 5 are shown in Figure 8 and Figure 9 with and without the added pieces, respectively.
Since the magnetic flux remains of the same magnitude for the whole stroke length, the end effects are efficiently reduced. What is then left of the detent forces is the cogging force that increases linearly with increasing stator overlap. It also shows that the cogging force of q = 1 has a magnitude close to 50% at 50% stator overlap with an increased amplitude towards 100% stator overlap. This is in accordance with the increased amplitude of b 2 in [4] and is also indicated in Figure 1 for the slot component of the cogging force. A clear disadvantage is the increased attraction forces, as can be seen in Figure 10. The increased attraction force can increase the friction in the generator, which is a drawback.
The added end plates are of the same dimensions as the stator yoke. Still, one can observe from the attraction force that there is a slightly higher flux in the end plates than the stator yoke. This is because of the nonlinear BH curve of the iron and a flux density closer to saturation in the stator teeth.
The second way, having extended stator, is illustrated in Figure 4b and the detent force in Figure 11 for q = 1 and in Figure 12 for q = 6 5 .
By adjusting the end lengths of the stator, one can change θ 1 of (7) to counter the end forces from both ends. In the figures, it is shown that there is a large reduction in the detent force when the linear generator has a full active stator area. It is, however, slightly higher than its counterpart in Figure 8 and Figure 9. The new, but much lower, force peaks arise when one stator end is next to the mid pole and the other stator end is next to the mid PM. The detent force is slightly distorted for q = 6 5 due to the new force peaks. The reduction is still visible, especially for q = 6 5 . Even with less saturation of the outer teeth, there is still a difference in magnitude between the phase voltages, as illustrated in Figure 13.
As earlier mentioned, this is because the winding patterns for the phases have different contributions in the end regions. The extended stator ends give little increase of the attraction forces, as shown in Figure 10.

3.3. Comparison in Force Reduction between Buried and Surface Topology and under Load Condition

The no-load voltage of the surface-mounted machine can be seen in comparison with the buried topology in Figure 14.
The dimensions of the surface-mounted topology is adjusted to match the rms values of the buried topology at full active stator area. The time interval starts from roughly t = 0.7 s and is set when one of the phases intersects the 0 voltage line and continues for one electric period. A matlab script is used to calculate the interpolated value for the start. The difference in mean rms between the two topologies is less than 0.2 promille for the time interval. By not splitting the outer magnets of the surface-mounted PM, the voltage is slightly higher when there is no full active stator area compared to the buried topology. The lower phase voltages for the surface mounted on the right side in Figure 14 are due to the lower active stator area for the same time instant.
The values used for the extended stator of both topologies are derived from parametric sweeps at the full active stator area where the detent force is calculated. To get a more precise value for the stator end extensions, algorithmic optimization could be used with an increased number of allowed decimals. The parametric sweeps gave, however, optimal stator extensions of 15 mm and 20 mm for buried and surface topology, respectively. This equals half of the pole width for both topologies.
The detent force for both topologies with and without the extended stator (ES) is shown in Figure 15.
The range of full active stator area is, as mentioned earlier, different for the surface-mounted topology. This can be realized by the extra negative force peak of the buried topology at t = 1.287 s. The figure shows similar behavior for both topologies when introducing the extended stator; a time shift of the force peak outside the full active stator area and a clear reduction in the forces within the full active stator area. The pair-wise connection of the surface-mounted PMs also introduces a difference in end forces when an even versus odd number of PMs overlap the stator. This phenomena has been known for quite some time [2]. One can also see in the figure that the detent force is slightly higher for the surface-mounted topology.
In Figure 16, the surface-mounted topology can be seen with and without extended stator at t = 0.756 s. This is the time instant where a force peak arises without the extended stator.
Comparing the two topologies, one can see some clear similarities; the extended stator for the buried topology covers the outer edges of the poles, and the extended stator for the surface-mounted topology covers the outer edges of the permanent magnets; i.e., for either topology of this study, a significant reduction in end forces is achieved when the poles closest to the air gap is completely covered.
As for the attraction forces in Figure 17, one can see similar behavior in both topologies, with a smaller change in ripple for the surface topology.
Lastly, the behavior of the detent force is investigated when the machine is subjected to electrical loading in the stator windings. Each phase has a current density of 2 A / mm 2 RMS, taking the fill factor into account and is ramped up for the three phases before they reach their designated phase current. The electromagnetic (EM) power associated with the detent force and constant translator velocity, P EM = F den v and the resistive-free electrical power with constant-torque-angle P el = E I can be seen in Figure 18 for one electric period (or two pole-pitches). From the figure, one can see that the electrical power is of similar magnitude as the EM power. One can also see that the ripple is higher for both topologies without the extended stator. To further analyze the ripple, the detent force ripple for the second pole pitch can be seen in Figure 19. To give a fair comparison of the force ripple, they are all offset by their average number (i.e., the averaged value for each case is zero). The fraction of force ripple with and without extended stator f Δ , t o p , J = r i p p l e E S r i p p l e for each topology (top) and current density (J) is; f Δ , S , 0 = 0.22 , f Δ , S , 2 = 0.45 , f Δ , B , 0 = 0.23 , f Δ , B , 2 = 0.36 . Since the velocity remains constant in the simulation, the same fraction is valid for both mechanical power ripple and force ripple. From these results, one can see that the extended stator reduces the force ripple of both topologies by similar factors at no-load. Under a load, however, the force ripple reduction due to the extended stator is less than under no load, with a better maintained reduction for the buried topology.

4. Conclusions

The electromagnetic forces were investigated for linear machines of buried and surface-mounted topologies, and means of reducing the end forces were included in a simulated FE environment. Initial investigation of the buried topology shows attraction forces and cogging forces that increase with active stator area. Meanwhile, the end forces did not increase in magnitude with an increased active stator area.
It was found that adding soft magnetic material for the whole stroke length reduced the detent force for the whole operating range but almost doubled the attraction forces, making it only viable for a linear machine where stroke length-to-machine length ratio is very low. If, instead, extending the stator ended slightly, the detent force was reduced for full active stator area by countering the end forces on opposite sides of the generator. This end effect compensation gave little increase to the attraction forces and was significant for both q = 1 and q = 6 5 , with a higher reduction in detent force when q was a fractional number.
The detent force of the surface-mounted topology varied in magnitude, depending on whether the number of active poles was odd or even. At full active stator area, the extended stator reduced the detent force in a similar manner as the buried topology. Although the two topologies showed somewhat different behavior in how the detent force acts with the number of active poles, they showed some clear similarities with regards to end effect compensation: the optimal stator extension to reduce the end forces was achieved for both topologies when the poles closest to the air gap overlapped the extended stator on both ends.
Furthermore, during the load condition, it was shown that the extended stator reduced the force (or power) ripple for both topologies. The reduction was, however, lower under a load than under no load and with a slightly better result for the buried topology.
The results also indicate that the buried topology better utilizes the stator, as the surface-mounted topology has a lower full active stator area for the same stroke length, effectively reducing the average induced voltage for a stroke length with partial stator area. A partial translator overlap for the surface-mounted topology inherently introduces magnetic flux variations in the translator, which can be a source for induced eddy currents.
In a broader context, one can see that the end effects that are inherent for the longitudinal ends of linear machines depend on the soft magnetic material of the stator ends. By modifying the stator ends, one can influence the force peaks. The modifications presented in the paper can be done regardless of the winding pattern.

Author Contributions

This was a collaborative work between J.S., A.E.F., M.L. and S.E. J.S. designed and performed the simulations and wrote most of the paper. A.E.F., M.L. and S.E. contributed significantly to the analysis of the results and the revision of the text. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by Uppsala University and the Swedish strategic research program StandUp for Energy.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
ABSAbove and below the stator
PMPermanent Magnet
LPMSGLinear Permanent Magnet Synchronous Generator
BHMagnetic flux density-Magnetic field strength dependency
MSTMaxwell’s Stress Tensor
VWVirtual Work
BBuried (translator topology)
Ssurface-mounted (translator topology)
EPEnd plates
ESExtended stator

References

  1. Sjölund, J.; Frost, A.E.; Leijon, M.; Eriksson, S. End Effects and Geometric Compensation in a Linear Permanent Magnet Synchronous Generator with Buried Topology. In Proceedings of the 2020 International Conference on Electrical Machines (ICEM), Gothenburg, Sweden, 23–26 August 2020; Volume 1, pp. 455–461. [Google Scholar] [CrossRef]
  2. Danielsson, O.; Leijon, M. Flux Distribution in Linear Permanent-Magnet Synchronous Machines including Longitudinal End Effects. IEEE Trans. Magn. 2007, 43, 3197–3201. [Google Scholar] [CrossRef]
  3. Liu, C.; Gao, H.; Xiong, Y.; Zhou, S.; Fu, W. Detent Force Reduction in Permanent Magnet Linear Synchronous Motor Base on Magnetic Field Similarity Method. IEEE Access 2019, 7, 57341–57348. [Google Scholar] [CrossRef]
  4. Lejerskog, E.; Strömstedt, E.; Savin, A.; Boström, C.; Leijon, M. Study of the Operation Characteristics of a Point Absorbing Direct Driven Permanent Magnet Linear Generator deployed in the Baltic Sea. IET Renew. Power Gener. 2016, 10, 1204–1210. [Google Scholar] [CrossRef] [Green Version]
  5. Ekergård, B.; Leijon, M. Longitudinal End Effects in a Linear Wave Power Generator. Energies 2020, 13, 327. [Google Scholar] [CrossRef] [Green Version]
  6. Frost, A.E.; Ulvgård, L.; Sjökvist, L.; Eriksson, S.; Leijon, M. Partial Stator Overlap in a Linear Generator for Wave Power: An Experimental Study. J. Mar. Sci. Eng. 2017, 5, 53. [Google Scholar] [CrossRef] [Green Version]
  7. Yu, H.; Liu, C.; Yuan, B.; Hu, M.; Huang, L.; Zhou, S. A Permanent Magnet Tubular Linear Generator for Wave Energy Conversion. J. Appl. Phys. 2012, 111, 07A741. [Google Scholar] [CrossRef]
  8. Almoraya, A.A.; Baker, N.J.; Smith, K.J.; Raihan, M.A.H. Development of a Double-sided Consequent Pole Linear Vernier Hybrid Permanent-Magnet Machine for Wave Energy Converters. In Proceedings of the 2017 IEEE International Electric Machines and Drives Conference (IEMDC), Miami, FL, USA, 21–24 May 2017; pp. 1–7. [Google Scholar] [CrossRef] [Green Version]
  9. Almoraya, A.A.; Baker, N.J.; Smith, K.J.; Raihan, M.A.H. A New Configuration of a Consequent Pole Linear Vernier Hybrid Machine with V-shape Magnets. In Proceedings of the 2018 XIII International Conference on Electrical Machines (ICEM), Alexandroupoli, Greece, 3–6 September 2018; pp. 2002–2008. [Google Scholar] [CrossRef]
  10. Huang, X.; Ji, T.; Li, L.; Zhou, B.; Zhang, Z.; Gerada, D.; Gerada, C. Detent Force, Thrust, and Normal Force of the Short-Primary Double-Sided Permanent Magnet Linear Synchronous Motor With Slot-Shift Structure. IEEE Trans. Energy Convers. 2019, 34, 1411–1421. [Google Scholar] [CrossRef]
  11. Hultman, E.; Ekergård, B.; Kamf, T.; Salar, D.; Leijon, M. Preparing the Uppsala University Wave Energy Converter Generator for Large-scale Production. In Proceedings of the 5th International Conference on Ocean Energy, Halifax, NS, Canada, 4–6 November 2014. [Google Scholar]
  12. Kamf, T. Automated Production Technologies and Measurement Systems for Ferrite Magnetized Linear Generators. Ph.D. Thesis, Acta Universitatis Upsaliensis, Uppsala, Weden, 2017. [Google Scholar]
  13. Ekergård, B.; Leijon, M. Ideal Analytical Expression of the Magnetic Circuit in a Permanent Magnet Linear Machine. In Proceedings of the 2013 Eighth International Conference and Exhibition on Ecological Vehicles and Renewable Energies (EVER), Monte Carlo, Monaco, 27–30 March 2013; pp. 1–6. [Google Scholar] [CrossRef]
  14. Lu, Q.; Huang, L.; Ye, Y.; Huang, X.y.; Fang, Y. Design of a Novel Permanent Magnet Linear Synchronous Motor with Segmented Armature Core for Ropeless Lifter. COMPEL-Int. J. Comput. Math. Electr. Electron. Eng. 2016, 35, 556–571. [Google Scholar] [CrossRef]
  15. Xia, T.; Yu, H.; Shi, Z.; Guo, R. Comparative Analysis and Experimental Verification of a Linear Tubular Generator for Wave Energy Conversion. Energies 2018, 11, 1707. [Google Scholar] [CrossRef] [Green Version]
  16. Botha, C.D.; Kamper, M.J.; Wang, R.J.; Sorgdrager, A.J. Force Ripple and Cogging Force Minimisation Criteria of Single-Sided Consequent-Pole Linear Vernier Hybrid Machines. In Proceedings of the 2020 International Conference on Electrical Machines (ICEM), Gothenburg, Sweden, 23–26 August 2020; Volume 1, pp. 469–475. [Google Scholar] [CrossRef]
  17. Yao, Y.; Chen, Y.; Lu, Q.; Huang, X.Y.; Ye, Y. Analysis of Thrust Ripple Harmonics on different Electric Loadings in PM Linear Synchronous Machines with Skewed PMs. COMPEL-Int. J. Comput. Math. Electr. Electron. Eng. 2016, 35, 655–669. [Google Scholar] [CrossRef]
  18. Polinder, H.; Damen, M.E.; Gardner, F. Linear PM Generator System for Wave Energy Conversion in the AWS. IEEE Trans. Energy Convers. 2004, 19, 583–589. [Google Scholar] [CrossRef] [Green Version]
  19. Ulvgård, L.; Sjökvist, L.; Göteman, M.; Leijon, M. Line Force and Damping at Full and Partial Stator Overlap in a Linear Generator for Wave Power. J. Mar. Sci. Eng. 2016, 4, 81. [Google Scholar] [CrossRef] [Green Version]
  20. Pile, R.; Devillers, E.; Le Besnerais, J. Comparison of Main Magnetic Force Computation Methods for Noise and Vibration Assessment in Electrical Machines. IEEE Trans. Magn. 2018, 54, 1–13. [Google Scholar] [CrossRef]
  21. Pyrhönen, J.; Jokinen, T.; Hrabovcova, V. Design of Rotating Electrical Machines. In Design of Rotating Electrical Machines, 2nd ed.; John Wiley & Sons, Ltd.: Chichester, UK, 2014; Chapter 1; pp. 1–45. [Google Scholar]
  22. Meessen, K.J.; Paulides, J.J.H.; Lomonova, E.A. Force Calculations in 3-D Cylindrical Structures Using Fourier Analysis and the Maxwell Stress Tensor. IEEE Trans. Magn. 2013, 49, 536–545. [Google Scholar] [CrossRef]
  23. Shin, K.; Kim, K.; Hong, K.; Choi, J. Detent Force Minimization of Permanent Magnet Linear Synchronous Machines Using Subdomain Analytical Method Considering Auxiliary Teeth Configuration. IEEE Trans. Magn. 2017, 53, 1–4. [Google Scholar] [CrossRef]
  24. Arkkio, A. Analysis of Induction Motors based on the Numerical Solution of the Magnetic Field and Circuit Equations. Ph.D. Thesis, Helsinki University of Technology, Helsinki, Finland, 1987. [Google Scholar]
  25. Eriksson, S. Design of Permanent-Magnet Linear Generators with Constant-Torque-Angle Control for Wave Power. Energies 2019, 12, 1312. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The total detent force calculated with Maxwell stress tensor (MST) and virtual work (VW) for one pole pitch τ p , starting at t = 0.48 s and ending at t = 0.6 s. The contribution of VW from all domains (VW All), above and below stator (ABS), and in the stator slots ( S s l o t s ) are plotted.
Figure 1. The total detent force calculated with Maxwell stress tensor (MST) and virtual work (VW) for one pole pitch τ p , starting at t = 0.48 s and ending at t = 0.6 s. The contribution of VW from all domains (VW All), above and below stator (ABS), and in the stator slots ( S s l o t s ) are plotted.
Designs 05 00064 g001
Figure 2. Illustrations for the buried topology. (a) shows the initial position with 50% stator overlap and the coordinate system. The translator moves in the negative y-direction in the simulations. The arrows indicate PMs of alternating direction, and the inclined lines indicate the domains for silicon steel. The winding patterns are shown for q = 1 and q = 6 / 5 in (b) and (c), respectively.
Figure 2. Illustrations for the buried topology. (a) shows the initial position with 50% stator overlap and the coordinate system. The translator moves in the negative y-direction in the simulations. The arrows indicate PMs of alternating direction, and the inclined lines indicate the domains for silicon steel. The winding patterns are shown for q = 1 and q = 6 / 5 in (b) and (c), respectively.
Designs 05 00064 g002
Figure 3. B-H curve for Silicon Steel NGO 50PN270 used in the translator and stator iron. The zoomed-in version shows the hysteresis-free consideration in this study.
Figure 3. B-H curve for Silicon Steel NGO 50PN270 used in the translator and stator iron. The zoomed-in version shows the hysteresis-free consideration in this study.
Designs 05 00064 g003
Figure 4. (a) End plate (EP) above and below the stator to keep the PMs magnetized and thereby reducing the end forces, and (b) extended stator (ES) iron to counter the end forces for full active stator area. The modifications are done to both top and bottom stator ends.
Figure 4. (a) End plate (EP) above and below the stator to keep the PMs magnetized and thereby reducing the end forces, and (b) extended stator (ES) iron to counter the end forces for full active stator area. The modifications are done to both top and bottom stator ends.
Designs 05 00064 g004
Figure 5. Surface-mounted (S) and buried (B) topologies at different time instants.
Figure 5. Surface-mounted (S) and buried (B) topologies at different time instants.
Designs 05 00064 g005
Figure 6. Detent forces of the buried topology for q = 1 and q = 6 5 . The stroke length allows for 50% active stator area. Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Figure 6. Detent forces of the buried topology for q = 1 and q = 6 5 . The stroke length allows for 50% active stator area. Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Designs 05 00064 g006
Figure 7. Phase voltages of the buried topology for q = 1 and q = 6 5 . The stroke length allows for 50% active stator area. Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines. The different phases are indicated by different colors.
Figure 7. Phase voltages of the buried topology for q = 1 and q = 6 5 . The stroke length allows for 50% active stator area. Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines. The different phases are indicated by different colors.
Designs 05 00064 g007
Figure 8. Detent forces for q = 1 with end plates (EP). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Figure 8. Detent forces for q = 1 with end plates (EP). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Designs 05 00064 g008
Figure 9. Detent forces for q = 6 5 with end plates (EP). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Figure 9. Detent forces for q = 6 5 with end plates (EP). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Designs 05 00064 g009
Figure 10. Attraction forces of the buried topology for q = 1 and q = 6 5 , and for both modifications to reduce the end forces. Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Figure 10. Attraction forces of the buried topology for q = 1 and q = 6 5 , and for both modifications to reduce the end forces. Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Designs 05 00064 g010
Figure 11. Detent forces for q = 1 with extended stator (ES). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Figure 11. Detent forces for q = 1 with extended stator (ES). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Designs 05 00064 g011
Figure 12. Detent forces for q = 6 5 with extended stator (ES). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Figure 12. Detent forces for q = 6 5 with extended stator (ES). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines.
Designs 05 00064 g012
Figure 13. Phase voltages of the buried topology with and without the extended stator (ES). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines. The different phases are indicated by different colors.
Figure 13. Phase voltages of the buried topology with and without the extended stator (ES). Full active stator area is from t = 0.6 s to t = 1.32 s, indicated by the dashed vertical lines. The different phases are indicated by different colors.
Designs 05 00064 g013
Figure 14. Phase voltages for q = 6 5 for both topologies. The stroke length allows for 50% active stator area. Full stator overlap for buried (B) and surface (S) topologies is from t = 0.6 s to t = 1.32 s and t = 1.2 s, respectively, indicated by the dashed vertical lines. The different phases are indicated by different colors.
Figure 14. Phase voltages for q = 6 5 for both topologies. The stroke length allows for 50% active stator area. Full stator overlap for buried (B) and surface (S) topologies is from t = 0.6 s to t = 1.32 s and t = 1.2 s, respectively, indicated by the dashed vertical lines. The different phases are indicated by different colors.
Designs 05 00064 g014
Figure 15. Detent force for buried (B) and surface-mounted (S) PMs with and without extended stator (ES). Full stator overlap for buried (B) and surface (S) topologies is from t = 0.6 s to t = 1.32 s and t = 1.2 s, respectively, indicated by the dashed vertical lines.
Figure 15. Detent force for buried (B) and surface-mounted (S) PMs with and without extended stator (ES). Full stator overlap for buried (B) and surface (S) topologies is from t = 0.6 s to t = 1.32 s and t = 1.2 s, respectively, indicated by the dashed vertical lines.
Designs 05 00064 g015
Figure 16. Magnetic flux density for the surface-mounted topology without (a) and with (b) extended stator (ES) at t = 0.756 s.
Figure 16. Magnetic flux density for the surface-mounted topology without (a) and with (b) extended stator (ES) at t = 0.756 s.
Designs 05 00064 g016
Figure 17. Attraction force for buried (B) and surface-mounted (S) PMs with and without extended stator (ES). Full stator overlap is from t = 0.6 s to t = 1.32 s and t = 1.2 s for buried and surface topologies, respectively, indicated by the dashed vertical lines. The far right vertical line indicate 50% active stator area for the surface-mounted topology.
Figure 17. Attraction force for buried (B) and surface-mounted (S) PMs with and without extended stator (ES). Full stator overlap is from t = 0.6 s to t = 1.32 s and t = 1.2 s for buried and surface topologies, respectively, indicated by the dashed vertical lines. The far right vertical line indicate 50% active stator area for the surface-mounted topology.
Designs 05 00064 g017
Figure 18. The lossless electrical power without the extended stator (ES) and the electromagnetic power (associated with the detent force) with and without ES at current density of 2 A / mm 2 RMS.
Figure 18. The lossless electrical power without the extended stator (ES) and the electromagnetic power (associated with the detent force) with and without ES at current density of 2 A / mm 2 RMS.
Designs 05 00064 g018
Figure 19. Detent force ripple for surface (solid) and buried (dashed), offset by the average value for the interval. The extended stator is indicated by orange and purple colors for both topologies. Blue and yellow colors indicate the stator without extension.
Figure 19. Detent force ripple for surface (solid) and buried (dashed), offset by the average value for the interval. The extended stator is indicated by orange and purple colors for both topologies. Blue and yellow colors indicate the stator without extension.
Designs 05 00064 g019
Table 1. Translator parameters.
Table 1. Translator parameters.
BuriedSurface
B r 0.451.2
μ r 1.0161.05
( B H ) max 39.67 kJ / m 3 -
pole-pitch τ p 60 mm60 mm
PM x 90 mm8.6 mm
PM y 30 mm40 mm
Yoke-25 mm
Translator ironSilicon Steel NGO 50PN270
Table 2. Stator parameters.
Table 2. Stator parameters.
Parameter Values
Slot depth81.875 mm
Slot face1.875 mm × 3.75 mm
Copper33%
Cable insulation (air)1 mm
Yoke24.8 mm
Stator ironSilicon Steel NGO 50PN270
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sjölund, J.; Frost, A.E.; Leijon, M.; Eriksson, S. End Effects and Geometric Compensation in Linear Permanent Magnet Synchronous Generators with Different Topologies. Designs 2021, 5, 64. https://0-doi-org.brum.beds.ac.uk/10.3390/designs5040064

AMA Style

Sjölund J, Frost AE, Leijon M, Eriksson S. End Effects and Geometric Compensation in Linear Permanent Magnet Synchronous Generators with Different Topologies. Designs. 2021; 5(4):64. https://0-doi-org.brum.beds.ac.uk/10.3390/designs5040064

Chicago/Turabian Style

Sjölund, Jonathan, Anna E. Frost, Mats Leijon, and Sandra Eriksson. 2021. "End Effects and Geometric Compensation in Linear Permanent Magnet Synchronous Generators with Different Topologies" Designs 5, no. 4: 64. https://0-doi-org.brum.beds.ac.uk/10.3390/designs5040064

Article Metrics

Back to TopTop