Next Article in Journal
Ab Initio Computations of O and AO as well as ReO2, WO2 and BO2-Terminated ReO3, WO3, BaTiO3, SrTiO3 and BaZrO3 (001) Surfaces
Previous Article in Journal
Crossover in Extended Newtonian Gravity Emerging from Thermodynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Research on Detent Force Characteristics of Winding Segmented Permanent Magnet Linear Synchronous Motor Based on Analytical Model

1
Department of Mechanical Engineering, Xi’an Jiaotong University, Xi’an 710049, China
2
Department of Mechanical and Electrical Engineering, Xi’an Polytechnic University, Xi’an 710048, China
*
Author to whom correspondence should be addressed.
Submission received: 21 April 2022 / Revised: 8 May 2022 / Accepted: 15 May 2022 / Published: 20 May 2022

Abstract

:
The winding segmented permanent magnet linear synchronous motor (WS-PMLSM) has not only a general slotting effect but also a unique two-end effect caused by the segmented stator and the finite length passive mover, which makes the detent force of WS-PMLSM possess the qualities of all types of linear motors. This paper adopts a physics-based multiple analytical method to study the detent force properties. First, an accurate sub-domain model is used to obtain the magnetic field distribution of the tooth slot, and the results of the analytical solution and the graph of the tooth slot effect are obtained using the scalar magnetic potential. Subsequently, the effect of the end unenergized core on the magnetic field is analyzed using Maxwell’s theory for the stator end effect, revealing that the end effect causes pulsating and reversing magnetic fields in the motor, which in turn creates mutual inductance with unequal phases, and generates thrust fluctuations. Finally, a complex relative air-gap permeability is introduced on the basis of the conformal mapping to equate the end effect generated by the permanent magnet. The accuracy of the proposed multiple analytical models was verified by the results of the finite element method (FEM) and experimental model.

1. Introduction

An important application of winding segmented permanent magnet linear synchronous motors (WS-PMLSM) in modern industry is material handling and processing [1]. The main purpose is to reduce labor. A common form is made into a ring structure so that it can be recycled, and the transmission moving parts are called movers, which are passive and do not require any energy or information transfer. The track consists of segmented stators, equivalent to multiple motors connected in series, each stator being separately powered and controlled [2]. Figure 1 shows an example of an application in a battery-powered line, where multiple movers operate highly independently on the same track, with four stations (P1–P4) distributed along the track at the desired locations [3].
Obviously, the detent force of WS-PMLSM is a problem in the processing area, and it also tends to cause workpiece loosening and running noise, which seriously affects the smoothness of the transfer system and is even more impermissible in the case of transferring liquids. This detent force is mainly caused by the slot effect and the longitudinal end effect [4]. The slot effect is more common and the corresponding literature is more abundant. This literature can be summarized into two methods; one is to improve the motor structure [5,6,7,8], and the other is to use special control methods [9,10,11].
However, compared to conventional linear motors, the WS-PMLSM generates two types of end effects, defined as secondary end effects and primary end effects, due to the limited length of the secondary module of the permanent magnet and the primary module of the winding, which increases the complexity of implementing thrust fluctuation suppression for such motors.
The secondary end effect is an important factor for the thrust fluctuation generated by WS-PMLSM, which is also a unique phenomenon of WS-PMLSM [12]. It is caused by the magnetic field of the permanent magnets at the secondary end, mainly in the form of the magnetic field of the secondary pole affecting the magnetic field distribution inside the motor through the magnetic field association of the adjacent poles.
The primary end effect is related to the discontinuous distribution of the winding, similar to the commonly analyzed end effect in linear induction motors [13], but it is reflected in the change of electromagnetic parameters only when there is a current excitation in the winding, due to the abnormal magnetic field distribution at the end [14]. Due to the presence of closely connected cores outside the energized modules and the variation of magnetic fields, less of the literature has focused on this aspect.
To effectively predict the detent forces due to end effects and slotting effects, different methods have been developed, which can be divided into four main categories: analytical, semi-analytical, numerical, and hybrid methods (a combination of all three) [15]. Numerical methods are very flexible for a wide range of geometries, including nonlinear and non-homogeneous materials, and can be combined with other physics. One of the more representative numerical methods is the FEM [16,17]; however, the high level of complexity in FEM leads to very time-consuming simulations that slow down the design process. Therefore, it seems necessary to develop a simple analytical model for the initial design that is very fast in computation and provides more physical insight. Several analytical methods have been proposed, such as the magnetic circuit method [18], the hierarchical model method [19], the integral transformation method [20], the field–circuit coupled method [21], the air gap magnetic permeability method, the winding function method [22] and the Maxwell mathematical physics method [23]. Among all analytical techniques, a semi-analytical method called the subdomain model (SDM) has been developed that can calculate magnetic fields very accurately and much faster than numerical methods (e.g., FEM), making it a good method to study and deploy.
The SDM approach is to simplify the motor model by dividing it into several regular domains in which the problem is linearized. Liu et al. [24] first introduced an accurate subdomain modeling approach for slot models with pole pair conversion. Guo and Hu [25] proposed a new analytical model for long dynamic magnetic linear motors, which uses a method that scales up the diameter of the rotating motor to a certain degree so as to approach the linear motor model, the results obtained from this analytical model depend on the choice of radius. An effective analytical model for PMLSM is presented in [26,27], which takes into account the end effects by adding virtual slots and slotting effects by using relative air-gap permeability. These methods try to solve all magnetic resistance problems of linear motors; however, the nonlinear magnetic field distribution at the end of WS-PMLSM requires the division of more subdomains. Therefore, in this paper, only the slotting effect is studied by using the mature SDM method, and the end effect is analyzed by other methods.
In view of the limitations of SDM in the application of WS-PMLSM, it is necessary to integrate various magnetic field solution methods to analyze and study the end effect of WS-PMLSM magnetically. The Conformal Map Method (CM) [28] is a powerful analytical tool for determining two-dimensional magnetic fields. Using the CM method, a complex magnetic field structure can be transformed into a simpler linear magnetic field structure so that it is easier to solve the differential equations. For example, a special CM, called the Schwarz–Christoffel (SC) mapping method [29], is used to determine end effects, but all of the literature is on linear-induced motors, and there is very little literature on linear motors where the secondary is a permanent magnet.
For the end effect of WS-PMLSM with open core and segmented power supply, there is even less literature, generally using the direct table look-up method of current-flux linkage [30,31], FEM [32], etc. The WS-PMLSM has the process of moving through the stator armature segmentation [10] and also the pulsating air gap magnetic field caused by the inductance imbalance caused by the adjacent segmented cores [33], which further complicates the end effect.
In this paper, multi-analysis models of WS-PMLSM considering the effect of both ends, and the effect of the slot are proposed. In Section 2, the effect of the slot is analyzed with an exact subdomain model, and the SD model of the SW-PMLSM is developed to derive the analytical solution of the air-gap density and obtain the detent force of the slot. In Section 3, the secondary end effect is further investigated based on the conformal map and the complex relative air-gap permeability, which leads to the derivation of the effect law of the permanent magnet end on the detent force. In Section 4, based on the classical magnetic potential theory, the magnetic potential of the single-phase winding is derived by considering the energized section of the stator and the non-energized stator at its two ends as a whole, which leads to end effects in the windings. In Section 5, the relevant conclusions are verified by finite elements and experiments.

2. Analysis of the Slot Effect

The effect of the slots is analyzed using an accurate subdomain model in which a scalar magnetic potential is used so that the effect of the slots can be simulated using the magnetic field distribution in each subdomain.

2.1. Basic Assumptions

The following basic assumptions are made for WS-PMLSM, the structure is as Figure 2.
  • The permeability of the core and support of the stator is infinite.
  • The permeability of the permanent magnet is equal to the permeability of the air, i.e., μ r = μ 0 .
  • Permanent magnet conductivity and eddy currents are ignored.
  • The motor is symmetric and the effect in the z-direction is negligible.
  • Permanent magnets can only be magnetized in the x- or y-direction.
The exact subdomain model of WS-PMLSM is shown in Figure 3. Where l m is the pole length, τ is the pole pitch, h m is the permanent magnet height, δ is the air gap length, h s is the slot height, b s is the slot width, b t is the tooth width, t s is the tooth pitch, and the total number of slots is z 1 .
The relationship between the magnetic flux density B and the magnetic field strength H in the region of air and permanent magnets is as follows.
In the PM subdomain
B 1 = μ 0 μ r H 1 + μ 0 M
In the air gap subdomain
B = μ 0 H
where μ r is the relative permeability of PM, μ 0 is the permeability of air gap, M is the magnetization of PM. In this way, the vector equation of PM can be written in the following form:
M = M x i + M y j
where M x and M y denote the components of M in the x and y directions, respectively, and i and j are the unit vectors in the x and y directions, and in this paper, M x = 0.
Currently, there are two main methods for the numerical solution of the permanent magnetic field. One is to use the vector magnetic potential A as the solution variable of the whole solution field, and this algorithm has many degrees of freedom, takes up a large amount of memory, takes a long time to solve, and requires high computer resources; the other is to use the scalar magnetic potential ϕ as the solution variable of the solution field, and this algorithm has the advantages of fewer variables, less memory, and fast computation. Therefore, the scalar magnetic potential ϕ is used to describe the magnetic field distribution, the magnetic field strength can be expressed in terms of the magnetic scalar potential as follows:
H x = ϕ x H y = ϕ y
The number of slots under a pair of poles is:
z p = z 1 p
where p is the number of polar pairs. Taking the center of the z p th slot as the initial position and fixing the two-dimensional Cartesian coordinates on the stator, the initial position of the ith slot is defined as:
x i = b s 2 + t s ( i 1 ) i = 1 , 2 , , z p

2.2. Scalar Magnetic Potential Flux Solution for Each Subdomain

According to the relationship between the flux density and the magnetic field strength, the magnetic field of the motor can be divided into three regions: permanent magnet subdomain-1, air gap subdomain-2 and slot subdomain-3i i = 1 , 2 , , z p .
Here, Cartesian coordinates are used to calculate the magnetic field, and when the Laplace equation is used for each of the above three subdomains with the scalar magnetic potential as the variable, the control equation for each subdomain is as follows.

2.2.1. The PM Subdomain

According to Equation (1), the scalar magnetic potential ϕ of region 1 can be rewritten as:
· B 1 = μ 0 μ r 2 ϕ 1 + · μ 0 M = 0
2 ϕ 1 = · H 1 μ r
Since the magnet pole is magnetized only in the Y direction, M x = 0 , and for the convenience of description, let M ( x ) = M y ( x ) , the magnetization intensity distribution function is shown in Figure 4.
The permanent magnet is equated to the magnetization intensity distribution function M ( x ) , for which a Fourier series expansion is performed to obtain [34]:
M ( x ) = n = 1 , 3 M n sin K n x x 0
where
M n = ( 1 ) n + 3 2 4 n π M 0 sin n π 2 α p
where α p is the permanent magnet pole arc coefficient, x 0 is the initial position of the PM mover, α p = l m / τ , M 0 is the permanent magnet magnetization strength, M 0 = B r / μ 0 , B r is the permanent magnet remanence, μ 0 is the air permeability, and k n = n π / τ . Let a = h m , b = h m + δ , c = h m + δ + h s , yield: b = a + δ , c = b + h s = a + δ + h s . Therefore, the M ( x ) can also be expanded to the following equation:
M ( x ) = n = 1 , 2 , 3 M x n s sin K n x M x n c cos K n x
M x n s = M n cos K n x 0 M x n c = M n sin K n x 0
· M = M x x + M y y = 0
According to Equation (8), the Poisson equation for the PM region can be written as:
2 ϕ 1 x 2 + 2 ϕ 1 y 2 = 0

2.2.2. The Air Gap and Slot Regions

The scalar magnetic potential of air gap subdomain 2:
2 ϕ 2 x 2 + 2 ϕ 2 y 2 = 0
The ith ( i = 1 , 2 , , z p ) slot subdomain 3 i :
2 ϕ 3 i x 2 + 2 ϕ 3 i y 2 = 0
The permanent magnet subdomain 1 scalar magnetic potential flux solution [35] is
ϕ 1 ( x , y ) = n = 1 , 2 , A 1 ( n ) e k n y + B 1 ( n ) e k n y cos k n x + n = 1 , 2 , C 1 ( n ) e k n y + D 1 ( n ) e k n y sin k n x
The air gap subdomain 2 scalar magnetic potential flux solution is
ϕ 2 ( x , y ) = n = 1 , 2 , A 2 ( n ) e k n y + B 2 ( n ) e k n y cos k n x + n = 1 , 2 , C 2 ( n ) e k n y + D 2 ( n ) e k n y sin k n x
For the slot region, assuming that the depth of the slot is infinite, the boundary conditions for the scalar magnetic potential are as follows:
ϕ 3 i = 0 , x = x i 0 , x = x i + b s ϕ s ( x , y ) , y = b 0 , y = c
where ϕ s ( x , y ) is the magnetic potential of the toothed surface.
ϕ 3 i = m C 3 i ( m ) e F m y + D 3 i ( m ) e F m y · sin F m x x i
where F m = ( m π / b s ) . Due to ϕ 3 i y = + = 0 , so C 3 = 0 , Equation (20) can also be written as:
ϕ 3 i = m D 3 i ( m ) e F m y sin F m x x i

2.3. Boundary Conditions

As shown in Figure 3, according to the relationship between the boundaries of each sub-domain, the equations of the coefficients of each boundary condition are established. The boundary conditions at the junction of each sub-domain are as follows:
H 1 x y = 0 = 0 H 1 x y = h m = H 2 x y = h m B 1 y y = h m = B 2 y y = h m
Furthermore, on the interface between the stator and the air gap:
H 3 i x y = h m + δ = H 2 x y = h m + δ x x i , x i + b s 0
B 3 i y y = h m + δ = B 2 y y = h m + δ
By applying the above boundary conditions and interface conditions, the unknown coefficients from A 1 to D 1 and from A 2 to D 2 can be determined, and the derivation process can be referred to in Appendix A.

2.4. Analysis of Magnetic Fields in Each Subdomain

Substituting A 1 ( n ) , B 1 ( n ) , C 1 ( n ) , D 1 ( n ) , A 2 ( n ) , B 2 ( n ) , C 2 ( n ) , D 2 ( n ) into ϕ 1 ( x , y ) , ϕ 2 ( x , y ) and ϕ 3 ( x , y ) gives the magnetic potential of each subdomain. According to the relationship between scalar magnetic potential and magnetic density:
B 1 y ( x , y ) = μ 0 μ r n = 1 , 2 k n A 1 ( n ) e k n y B 1 ( n ) e k n y cos k n x + n = 1 , 2 k n C 1 ( n ) e k n y D 1 ( n ) e k n y sin k n x + μ 0 M y ( x )
B 1 x ( x , y ) = μ 0 μ r n = 1 , 2 k n A 1 ( n ) e k n y + B 1 ( n ) e k n y sin k n x + n = 1 , 2 k n C 1 ( n ) e k n y + D 1 ( n ) e k n y cos k n x
B 2 y ( x , y ) = μ 0 n = 1 , 2 k n A 2 ( n ) e k n y B 2 ( n ) e k n y cos k n x + n = 1 , 2 k n C 2 ( n ) e k n y D 2 ( n ) e k n y sin k n x
B 2 x ( x , y ) = μ 0 n = 1 , 2 k n A 2 ( n ) e k n y + B 2 ( n ) e k n y sin k n x + n = 1 , 2 k n C 2 ( n ) e k n y + D 2 ( n ) e k n y cos k n x
B 3 i y ( x , y ) = μ 0 m = 1 , 2 C 3 i ( m ) F m e F m y 1 + e 2 F m ( c y ) sin F m x x i
B 3 i x ( x , y ) = μ 0 m = 1 , 2 C 3 i ( m ) F m e F m y 1 + e 2 F m ( c y ) cos F m x x i

3. Analysis of Windings End Effect

Figure 5 shows a WS-PMLSM model with a three-phase full-pitch winding and an effective pole number of 2p. To simplify the calculation process, Where a , b , c are the magnetic flux paths of the three regions, τ is the pole pitch, H c is the stator height, δ is the air gap length. the following assumptions are made:
  • Neglecting the gearing effect of the winding, the magnetic field generated by each phase winding in the air gap can be seen as a result of the traveling wave current layer generated;
  • The magnetic field at x < l or x > 2 l is zero, the studied motor section is divided into three regions: ①, ②, ③, where the energized region is ①, and the non-energized regions are ①, ②
  • The air gap contains only the y-direction magnetic field, and the gap between segments is ignored.

Magnetic Field Equation

In zone ①, the equivalent piece line current density:
j 1 = J sm sin ω t α x
When considering only region ①:
B y 1 ( t ) = B m cos ω t α x + c B m = μ 0 τ π g s J sm
where α = π / τ , c is the end effect of ②, ③ on ①, J sm is the equivalent piece line current density amplitude, τ is the motor pole pitch, ω is the power supply angular frequency, applying Ampere’s law of cycles along the circuit shown yields:
B y 0 ( t ) g s μ 0 B y 1 ( t ) g s μ 0 = 0 x J sm cos ω t α x d x
where B y 1 ( t ) —① air gap magnetic density, B y 0 ( t ) —air gap magnetic density at x = 0 , μ 0 is the vacuum magnetic permeability, g s is the equivalent air gap length after taking into account the slot effect and the static lateral side end effect. From Equation (33), the following equation can be obtained by the quadratic differential of x:
d 2 B y 1 ( t ) d x 2 = μ 0 α J sm 2 g s sin ω t α x
Similarly, the magnetic field equations of the air gap in ② and ③ can be obtained:
d 2 B y 2 ( t ) d x 2 = 0 , d 2 B y 3 ( t ) d x 2 = 0
where B y 2 ( t ) and B y 3 ( t ) are the magnetic density of the air gap in ② and ③, respectively. The boundary conditions satisfied by region① ② ③ are:
B y 2 ( t ) x x = l = B y 3 ( t ) x x = 2 l = 0 B y 1 ( t ) B y 2 ( t ) x = 0 = 0 B y 1 ( t ) B y 3 ( t ) x = l = 0 l 0 B y 1 ( t ) d x + 0 l B y 2 ( t ) d x + l 2 l B y 3 ( t ) d x = 0
Applying Ampere’s law of loops at x = 0 , l
d B y 1 ( t ) d x + d B y 2 ( t ) d x = μ 0 J sm 2 g s cos ω t α x d B y 1 ( t ) d x + d B y 3 ( t ) d x = μ 0 J sm 2 g s cos ω t α x
From the above differential equations and the corresponding boundary conditions, we can obtain the air-gap magnetic field in each region in WS-PMLSM without considering the half-full slot and odd poles.
B y 1 ( t ) = B m cos ω t α x ( 1 ) ( p ) k s B m c h β x cos ω t + ( 1 ) ( p ) k h B m s h β x cos ω t B y 2 ( t ) = B y 3 ( t ) = ( 1 ) ( p ) k s B m cos ω t β = 2 μ 0 g s H c μ fe k s = s h β l s h β p τ
When the iron core is not saturated β 0 ,   k s = 2 ,   k h = 0
B y 1 ( t ) = B m cos ω t α x ( 1 ) ( p ) 2 B m cos ω t
From the above equation, it can be seen that due to the presence of unenergized segmented cores on both sides of the supply section, the magnetic field generated by each phase winding has a DC component in the energized segmented interval, resulting in a DC pulsating component in the synthesized magnetic field, forming a pulsating air-gap magnetic field. The effect of the length of the end core on the ratio of the pulsating magnetic field and the amplitude of the current traveling wave field formed by the static end effect can be expressed as:
k m 1 = l cos w t p τ cos w t θ 0
From Equation (40), it can be seen that the longer the core length, the larger the amplitude of the pulsating magnetic field; increasing the number of pole pairs can suppress the effect of the end effect, and increasing the length of the single-section primary and decreasing the length ratio of the secondary and primary can weaken the effect of the end effect. For WS-PMSM, the phase imbalance characteristic is an important characterization of the longitudinal end effect of the primary module, which has a significant impact on the control characteristics of the motor.

4. End Effects Caused by Secondary Finite Length

The secondary longitudinal end effect can also be called the secondary end effect. It is mainly manifested as the magnetic field of the secondary pole affecting the magnetic field distribution through the discontinuity of two adjacent stator sections, which in turn causes the thrust fluctuation of the motor.

Modeling of Secondary Longitudinal End Effects

To analyze the end magnetic field of the WS-PMLSM, the Schwarz–Christoffel transformation is used, and the specific variation process is shown in Figure 6a–c. In the derivation process, assuming that the scalar potential of the secondary is ω = 1 and the scalar potential of the primary armature core is ω = 1 , then there is a line of symmetry with a potential of ω = 0 at the center of the air gap, as shown in Figure 6a, and the problem can be reduced to a polygonal solution containing a right-angle coordinate system, such that by the Schwarz–Christoffel transformation and logarithmic transformation, the magnetic field in the z plane can be turned into a semicircular region in the ω plane, and finally, into a region of parallel line segments in the t plane to obtain a regular magnetic field pattern.
The transformation between the z plane, w plane, and t plane coordinate systems is performed according to Equation (41). The g in the equation represents the air gap.
z w : z = g 2 π 2 ( w + 1 ) 1 2 + ln ( w + 1 ) 1 2 1 ( w + 1 ) 1 2 + 1 w t : w = e π t or t = 1 π ln w z t : z = g π 2 e π t + 1 + ln e π t + 1 1 e π t + 1 + 1
In the t plane, the magnetic field is regular and the magnetic induction at point b in the z plane is considered to be maximum and uniform, and the magnetic potential set value A = 1 yields B m a x = 2 μ 0 / g (in the z-plane, when passing through the point B, across the trip g / 2 ), then the magnetic density at any point within the end field is:
B = μ 0 H = μ 0 d t d z = μ 0 d t d w d w d z = B max ( w + 1 ) 1 2
The magnetic density distribution in this region varies with position in an exponential-like pattern, and the expression is obtained from the curve fitting.
B e n d = B max e x g e x 0 B max e x + ( 2 p 1 ) τ + L m g e x ( 2 p 1 ) τ + L m
where g e is the length of the electromagnetic air gap after the introduction of the cardinality factor. Then, the equivalent air gap length at the secondary end is:
g e n d = g e e x g e x 0 g e e x [ ( 2 p 1 ) ) + L m g e x ( 2 p 1 ) τ + L m
The change of the air-gap ratio permeability function with the position in this area can be obtained by curve fitting as follows [36]:
B e n d B max = 0.6635 e 1.597 x + 0.2012 e 0.2335 x x > 0 0.1782 e 4.385 x + 1.009 e 0.0007946 x τ < x 0
The influence of the end effect caused by the finite length of the permanent magnet on the air-gap magnetic field of WS-PMLSM is shown in Figure 7.
Figure 7 shows magnetic flux density normalized curve B B max = f x g in the longitudinal end region. When x = g , the magnetic inductance B m a x , near the transverse edge of the armature core, and the magnetic inductance of the air gap gradually decays, and when x = −0.25 g, B = 0.95 B m a x , the magnetic inductance decays to 95% of the maximum value. When x = 0 , corresponding to the transverse edge of the armature core, i.e., point A in the z-plane, the magnetic inductance is B = 0.833 B m a x . That is, the original magnetic field is weakened, which corresponds to the generation of an additional magnetic field that weakens the main air gap field and causes a change in its distribution. When the WS-PMLSM is at no load, the additional magnetic field generated by the end effect is a constant. In terms of the magnetic circuit, the end effect changes the flux path of the end poles, the air-gap flux density under the surface of the end poles decreases, and the N-pole motor forms an (N + 1) magnetic circuit, as shown in Figure 8. Where B , B 1 , B 2 is B is the magnetic flux density test point, h m is the permanent magnet height, g is the air gap length, h s is the slot height, h c is the height of the stator support.
It can be seen that a total of three regions are generated: the end magnetic field region II ( x < 0 ) and the effective magnetic field region I ( 0 < x < τ ) , which forms an entire magnetic circuit system with area I and area III ( τ < x < 2 τ ) .

5. Validation of the Analytical Methods by the FEM and Experiments

5.1. Finite Element Flux Simulation and Experimental Setup

In order to evaluate the merits of the proposed analytical model, the FEM method and experimental validation were used. The main parameters of the prototype are listed in Table 1. In the FEM analysis method, the material of the primary and secondary backing iron was steel-1010, and the JMAG software was used for modeling and simulation. A passive mover with a width of 150 mm and a stator with a length of 1000 mm were manufactured. The stator was composed of five sets of three-phase winding modules. The passive mover was supported by a linear bearing consisting of LM tracks. The passive mover at the origin of the permanent magnet was measured by a non-contact linear encoder.
Figure 9 shows the experimental prototype of the proposed system.
Figure 10 shows the magnetic flux simulation. We chose full coverage and found that the two end effects and the tooth slot effect are very obvious. For further verification and comparison, the analytical solution and the finite element diagram are compared below.

5.2. Slot Flux Density Distribution in Each Subdomain

The initial positions are chosen at three different locations, namely, the surface of the permanent magnet, the surface at the junction of the air gap and the tooth slot, and the center of the air gap. Their normal and tangential flux densities, as determined by the analytical results and the FEM simulations, are compared in Figure 11, Figure 12 and Figure 13.
It can be seen that the flux density of each subdomain obtained by the accurate subdomain analysis model basically overlaps with the finite element results, which proves the accuracy of the analysis method. From the above figure, it can be seen that the winding slotting causes the change of air gap permeability and the change of air-gap flux density integrates the effect of all slotting effects, and the closer to the slotting surface, the more obvious the slotting effect is.

5.3. End Effects Caused by PM

From Figure 14, it can be seen that the secondary end effect reduces the magnetic density of the air-gap field of the end magnets, resulting in an uneven distribution of the air-gap field under the surface of the pole plate. Due to the structural symmetry of the “outgoing” and “incoming” ends of the motor, the magnitude of the magnetic induction losses at the two ends is equal. At the same time, the secondary end effect of the air-gap magnetic field of the WS-PMLSM causes the flux path at the ends to change, while the magnetic induction strength of the air-gap under the pole surface at the ends decreases.

5.4. End Effects Caused by Winding Segmentation

As shown in Figure 15 below, when three symmetrical voltages are applied to the motor, the impedance of each phase of the motor is asymmetrical due to the effect of winding segmentation, which will lead to the asymmetrical three-phase current of the motor. The air-gap magnetic field of the motor is an asymmetrical magnetic field, and the positive sequence current component of the motor will generate positive direction electromagnetic thrust. The negative sequence current component of the motor will generate a reverse direction electromagnetic thrust. At the same time, it will produce the asymmetry of the motor inductance matrix and impedance matrix, which will cause the pulsation of the air gap magnetic field and thus the thrust fluctuation.

5.5. Detent Force Calculation

On basis of the aforementioned analytical model, the x- and y-direction components of flux density can be accurately obtained. Using the Maxwell stress theory, the cogging force along the x-direction and the normal magnetic force along the y-direction caused by the effect of the slot and PM mover can be expressed as [37]
F x = L z μ 0 x 1 x 2 B 2 x B 2 y d x F y = L z 2 μ 0 x 1 x 2 B 2 y 2 B 2 x 2 d x
where L z is the longitudinal length of the linear motor along the z-direction, and B x and B y are the x- and y-direction components of flux density in region 2, respectively. The lower and upper limits of integration are x 1 = L 2 and x 2 = L 2 , respectively. For the proposed analytical method, the distance between x 1 and x 2 includes the adjacent unit windings not energized at both sides of the mover, as shown in Figure 8.
Finally, the three effects were verified in practice on the experimental platform. We used a thrust test to verify the fluctuations caused by three conditions. The characteristics of the cogging force in the x direction at different positions are shown in Figure 16. In the starting section, the deviation of the force is mainly caused by the end effect, then the amplitude decays as a function in the second transition section, and finally, it is completely generated by the slot in the third section. Fluctuations caused by the cogging effect are always present, which is in accordance with the experimental data.
The above results show that the cogging effect is caused by the structure and can be solved from the motor body design and control method, while the fluctuations generated by the end effect with high harmonics can be handled by the Taylor expansion method, which is not described here.

6. Conclusions

In this paper, the characteristic law of defense force in motor design and operation is revealed by the method of analytics with WS-PMLSM as the research object. The magnetic field analysis operation is carried out for two types of end effects and slot effects in the motor, which provides a theoretical basis for studying the fluctuation characteristics of electromagnetic force in this type of motor. The following research work is mainly accomplished in this paper.
(1) An accurate subdomain analysis model is proposed instead of the hierarchical model for solving the magnetic field distribution and thrust fluctuations of the slotted teeth of the WS-PMLSM motor. Using this analysis model, the magnetic field of the magnet, including the magnetized PM, and the effect of different slot widths on the thrust of the motor can be solved.
(2) For the end effect caused by the finite length of WS-PMLSM, the secondary is investigated, the longitudinal end effect formed by the permanent magnet secondary is equated with the edge effect at the end, and the Schwarz–Christoffel transform is used to solve such problem. This method is conventional but effective.
(3) In WS-PMLSM, the end effect caused by winding discontinuity is less studied. In this paper, Maxwell’s equations are used to study it, and some conclusions are obtained; namely, under the action of three symmetrical voltages, the motor generates three-phase asymmetric currents, and negative sequence currents generate reverse electromagnetic thrust. The impedance matrix of the motor is not a cyclic symmetric array at this time, so the sequence impedance cannot be transformed into a diagonal array.
In the next step, the authors will compare the finite element analysis with experimental results for the conclusions obtained in this paper. Thanks to the expression of the WS-PMLSM cogging torque waveform, the rotor eddy current loss due to the cogging effect is completed in the next study. Since the local values of flux density or vector potential are not easy to measure, this experiment is not sufficient, but it is possible to compare the integral magnitude, such as counter-electromotive force or torque, which will be completed in the next study.

Author Contributions

Conceptualization, T.Z.; methodology, T.Z.; validation, T.Z. and X.M.; formal analysis, T.Z.; investigation, X.M.; writing—original draft preparation, T.Z.; writing—review and editing, X.M.; supervision, X.M.; funding acquisition, X.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural ScienceFoundation of China under the grant 51735010.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The dataset can be accessed upon request to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

From the boundary condition in Equation (22), the following equations can be given
A 1 ( n ) + B 1 ( n ) = 0 C 1 ( n ) + D 1 ( n ) = 0
According to the boundary conditions in Equation (23), the following equations can be obtained:
A 1 ( n ) e k n h m + B 1 ( n ) e k n h m = A 2 ( n ) e k n h m + B 2 ( n ) e k n h m C 1 ( n ) e k n h m + D 1 ( n ) e k n h m = C 2 ( n ) e k n h m + D 2 ( n ) e k n h m
From the boundary condition Equations (1) and (23), the following equations can be obtained:
μ r k n A 1 ( n ) e k n h m B 1 ( n ) e k n h m M xc n = k n A 2 ( n ) e k n h m B 2 ( n ) e k n h m μ r k n C 1 ( n ) e k n h m D 1 ( n ) e k n h m M xs n = k n C 2 ( n ) e k n h m D 2 ( n ) e k n h m
From the boundary condition Equation (24), the magnetic potential at the slot surface and air gap boundary is continuous
ϕ 2 ( x , y ) y = m + δ = ϕ s ( x , y )
ϕ s ( x , y ) = ϕ 3 i ( x , y ) y = h m + δ = ϕ 3 i ( x , y ) y = h m + δ x x i , x i + b s 0 other
According to Equation (18)
ϕ 2 ( x , y ) y = h m + δ = n = 1 , 2 , A 2 ( n ) e k n h m + δ + B 2 ( n ) e k n h m + δ cos k n x + n = 1 , 2 , C 2 ( n ) e k n h m + δ + D 2 ( n ) e k n h m + δ sin k n x
The magnetic level on the surface of the slot is
ϕ s ( x , y ) = ϕ 3 i ( x , y ) y = h m + δ = m = 1 , 2 , C 3 i ( m ) χ m sin F m x x i = m = 1 , 2 , C e i ( m ) sin F m x x i x x i , x i + b s
where
C e i ( m ) = C 3 i ( m ) χ m χ m = e F m h m + δ 1 e 2 F m h s
The magnetic position of the tooth surface is
ϕ s ( x , y ) = 0 x x i , x i + b s
The Fourier series expansion of ϕ s ( x , y )
ϕ s ( x , y ) = n = 1 A s ( n ) cos k n x + B s ( n ) sin k n x
Although the period of the slot is not the same as the harmonic period of the PM, at boundary condition Equation (23), the harmonics of ϕ 2 ( x , y ) originate from the PM, so the period of the Fourier expansion is the period of the PM.
A s ( n ) = 1 τ x i x i + b s m = 1 , 2 , C 3 i ( m ) χ m sin F m x x i cos k n x d x = i = 1 z p m = 1 , 2 , C 3 i ( m ) χ m η si ( m , n )
B s ( n ) = 1 τ x i x i + b s m = 1 , 2 , C 3 i ( m ) χ m sin F m x x i sin k n x d x = i = 1 z p m = 1 , 2 , C 3 i ( m ) χ m ξ si ( m , n )
η s i ( m , n ) = η s ( m , n ) cos k n x i ξ s ( m , n ) sin k n x i
ξ si ( m , n ) = η s ( m , n ) sin k n x i + ξ s ( m , n ) cos k n x i
in the formula
η s ( m , n ) = m b s τ 1 cos ( m π ) cos k n b s π m 2 τ 2 n 2 b s 2 m τ n b s b s 2 t s sin m π b s x i m τ = n b s
ξ s ( m , n ) = m b s τ cos ( m π ) sin k n b s π m 2 τ 2 n 2 b s 2 m τ n b s b s 2 t s cos m π b s x i m τ = n b s
According to the boundary conditions Equation (24), the following equation set can be obtained:
A 2 ( n ) e k n h m + δ + B 2 ( n ) e k m h m + δ = A s ( n ) = i = 1 z p m = 1 , 2 , C 3 i ( m ) χ m η si ( m , n ) C 2 ( n ) e k n h m + δ + D 2 ( n ) e k n h m + δ = B s ( n ) = i = 1 z p m = 1 , 2 , C 3 i ( m ) χ m ξ si ( m , n )
From the boundary condition Equation (24), the normal magnetic flux density at the surface boundary of the air gap subdomain and the tooth slot subdomain is continuous.
B 2 y y = h m + δ = μ 0 ϕ 2 y y = h m + δ = μ 0 n = 1 , 2 , k n A 2 ( n ) e k n h m + δ B 2 ( n ) e k n h m + δ cos k n x μ 0 n = 1 , 2 , k n C 2 ( n ) e k n h m + δ D 2 ( n ) e k n h m + δ sin k n x
B 3 i y y = h m + δ = μ 0 ϕ 3 y y = h m + δ = μ 0 m = 1 , 2 C 3 i ( m ) δ m sin F m x x i
in the formula
δ m = F m e F m h m + δ 1 + e 2 F m h s
Since B 2 y is in the form of a general solution, a Fourier series expansion is made at the i-th slotted region x i x x i + b s
B 2 y y = h m + δ = μ 0 m = 1 , 2 , C s i ( m ) sin F m x x i x x i , x i + b s
where
C si ( m ) = 2 b s x i x i + b s n = 1 , 2 , k n A 2 ( n ) e k n h m + δ B 2 ( n ) e k n h m + δ . cos k n x sin F m x x i d x + 2 b s x i x i + b s n = 1 , 2 , k n C 2 ( n ) e k n h m + δ D 2 ( n ) e k n h m + δ . sin k n x sin F m x x i d x
Simplifying C si ( m ) yields
C si ( m ) = n = 1 , 2 , k n A 2 ( n ) e k n h m + δ B 2 ( n ) e k n h m + δ σ si ( m , n ) + n = 1 , 2 , k n C 2 ( n ) e k n h m + δ D 2 ( n ) e k n h m + δ τ si ( m , n )
σ s i ( m , n ) = 2 τ b s η s i ( m , n )
τ si ( m , n ) = 2 τ b s ξ si ( m , n )
According to the boundary condition Equation (24),
δ m C 3 i ( m ) = C s i ( m ) = n = 1 , 2 , · k n A 2 ( n ) e k n h m + δ B 2 ( n ) e k n h m + δ σ s i ( m , n ) + n = 1 , 2 , · k n C 2 ( n ) e k n h m + δ D 2 ( n ) e k n h m + δ τ s i ( m , n )
Summarizing the above equations, they can be concentrated as Equation (A27).
Afterward, these equations can be converted to matrix form. There are three parameters in the matrix m , n , i , where m is the harmonic coefficient of the stator slot. The larger m is, the more clearly the interaction between the slots will be revealed. n is the harmonic coefficient of the PM, and i is the number of slots, so the dimension of the matrix is m × n . In order to solve the convenient, the following formula is n = 1, m = 25, as in Equation (A28). Therefore, when n is taken in the range that it can be, it is possible to solve for the possible coefficients.
A 1 ( n ) + B 1 ( n ) = 0 C 1 ( n ) + D 1 ( n ) = 0 A 1 ( n ) e K n h m + B 1 ( n ) e K n h m A 2 ( n ) e K n h m B 2 ( n ) e K n h m = 0 C 1 ( n ) e K n h m + D 1 ( n ) e K n h m C 2 ( n ) e K n h m D 2 ( n ) e K n h m 0 μ r A 1 ( n ) e K n h m μ r B 1 ( n ) e K n h m A 2 ( n ) e K n h m + B 2 ( n ) e K n h m = 1 K n M x n c μ r C 1 ( n ) e K n h m μ r D 1 ( n ) e K n h m C 2 ( n ) e K n h m + D 2 ( n ) e K n h m = 1 K n M x n s A 2 ( n ) e K n ( h m + δ ) + B 2 ( n ) e K n ( h m + δ ) i m D 3 i ( m ) e F m ( h m + δ ) η s i ( m , n ) = 0 C 2 ( n ) e K n ( h m + δ ) + D 2 ( n ) e K n ( h m + δ ) i m D 3 i ( m ) e F m ( h m + δ ) δ s i ( m , n ) = 0 2 τ b s i n K n A 2 ( n ) e K n y 3 B 2 ( n ) e K n ( h m + δ ) η s i ( m , n ) + C 2 ( n ) e K n ( h m + δ ) D 2 ( n ) e K n ( h m + δ ) δ s i ( m , n ) + F m D 3 ( m ) e F m ( h m + δ ) = 0
1 1 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 E 3 n n E 4 n n 0 0 E 5 n n E 6 n n 0 0 0 0 0 E 3 n n E 4 n n 0 0 E 5 n n E 6 n n 0 E 7 n n E 8 n n 0 0 E 9 n n E 10 n n 0 0 0 0 0 E 7 n n E 8 n n 0 0 E 9 n n E 10 n n 0 0 0 0 0 E 11 n n E 12 n n 0 0 E 17 S n n 0 0 0 0 0 0 E 11 n n E 12 n n E 18 S n n 0 0 0 0 E 13 S n n E 14 S n n E 15 S n n E 16 S n n W m m A 1 B 1 C 1 D 1 A 2 B 2 C 2 D 2 D 3 = 0 0 0 0 B b 5 B b 6 0 0 0
E 3 n n ( n , n ) = e K n h m , E 4 n n ( n , n ) = e K n h m E 5 n n ( n , n ) = e K n h m , E 6 n n ( n , n ) = e K n h m E 7 n n ( n , n ) = μ r e K n h m , E 8 n n ( n , n ) = μ r e K n h m E 9 n n ( n , n ) = e K n h m , E 10 n n ( n , n ) = e K n h m E 11 n n ( n , n ) = e K n ( h m + δ ) , E 12 n n ( n , n ) = e K n ( h m + δ )
E 13 S m n = i n 2 T p P b K n e K n ( h m + δ ) η s i ( m , n ) E 14 S m n = i n 2 T p P b K n e K n ( h m + δ ) η s i ( m , n )
E 15 S m n = i n 2 T p P b K n e K n ( h m + δ ) δ s i ( m , n ) E 16 S m n = i n 2 T p P b K n e K n ( h m + δ ) δ s i ( m , n ) E 17 S m n = i m e F m ( h m + δ ) η s i ( m , n ) E 18 S m n = i m e F m ( h m + δ ) δ s i ( m , n )
W m m m i , m i = F m e F m ( h m + δ ) B b 5 = 1 K n M x n c , B b 6 = 1 K n M x n s . D 3 = D 3 i i = 1 , 2 , z n
where z n is the number of stator slots. The equations can be solved by matrix operations, which is a systematic and complex project due to its implementation by computer programming and will not be expanded on here.

References

  1. Mihalachi, M.; Leidhold, R.; Mutschler, P. Motion Control for Long Primary Linear Drives Used in Material Handling. In Proceedings of the 14th International Power Electronics and Motion Control Conference EPE-PEMC 2010, Ohrid, Macedonia, 6–8 September 2010; pp. T5-94–T5-101. [Google Scholar]
  2. Ranky, P.G. MagneMotion’s linear synchronous motor (LSM) driven assembly automation and material handling system designs. Assem. Autom. 2007, 27, 97–102. [Google Scholar] [CrossRef]
  3. Zhang, T.; Mei, X.; Du, X. Absolute Position Acquisition for Linear Synchronous Motor With Passive Mover. IEEE Access 2021, 9, 100757–100768. [Google Scholar] [CrossRef]
  4. Wang, M.; Li, L.; Pan, D. Detent Force Compensation for PMLSM Systems Based on Structural Design and Control Method Combination. IEEE Trans. Ind. Electron. 2015, 62, 6845–6854. [Google Scholar] [CrossRef]
  5. Huang, Q.; Huang, X.; Yu, W.; Ren, C. Modeling and Optimization of Detent Force for Flexible Permanent Magnet Synchronous Linear Motor Considering Spacing Between Primaries. In Proceedings of the 2021 13th International Symposium on Linear Drives for Industry Applications (LDIA), Wuhan, China, 1–3 July 2021. [Google Scholar]
  6. Shi, H.; Niguchi, N.; Hirata, K. Design and analysis of surface permanent magnet vernier linear motor based on air gap magnetic flux density distribution. In Proceedings of the 2019 12th International Symposium on Linear Drives for Industry Applications (LDIA), Neuchatel, Switzerland, 1–3 July 2019; IEEE: Piscataway, NJ, USA, 2019; pp. 1–5. [Google Scholar]
  7. Huang, X.; Tianpeng, J.; 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]
  8. Zahid, A.; Khan, F.; Ahmad, N.; Sami, I.; Alkhammash, H.I. Design and Analysis of Dual Mover Multi-Tooth Permanent Magnet Flux Switching Machine for Ropeless Elevator Applications. Actuators 2021, 10, 81. [Google Scholar] [CrossRef]
  9. Huang, X.; Liang, J.; Zhenyu, Q.; Li, J. An Iterative Estimation Algorithm of Prepositioning Focusing on the Detent Force in the Permanent Magnet Linear Synchronous Motor System. IEEE Trans. Ind. Electron. 2020, 67, 8252–8261. [Google Scholar] [CrossRef]
  10. Benavides Oswald, R. Investigation of Control Methods for Segmented Long Stator Linear Drives. Ph.D. Thesis, Technische Universität, Berlin, Germany, 2008. [Google Scholar]
  11. Cui, L.; Zhang, H.; Jiang, D. Research on High Efficiency V/f Control of Segment Winding Permanent Magnet Linear Synchronous Motor. IEEE Access 2019, 7, 138904–138914. [Google Scholar] [CrossRef]
  12. Lines, C.R. Modelling and Control of Segmented Long-Stator Permanent-Magnet Linear Synchronous Motors. Ph.D. Thesis, University of the Witwatersrand, Faculty of Engineering and the Built, Johannesburg, South Africa, 2013. [Google Scholar]
  13. Creppe, R.C.; Ulson, J.C.; Rodrigues, J.F. Influence of Design Parameters on Linear Induction Motor End Effect. IEEE Trans. Energy Convers. 2008, 23, 358–362. [Google Scholar] [CrossRef]
  14. Lu, J.; Ma, W. Research on End Effect of Linear Induction Machine for High-Speed Industrial Transportation. IEEE Trans. Plasma Sci. 2011, 39, 116–120. [Google Scholar] [CrossRef]
  15. Boff, B.H.B.; Eckert, P.R.; Amara, Y. A Comprehensive Review on the End Effects of Linear Permanent Magnet Machines. In Proceedings of the 2021 13th International Symposium on Linear Drives for Industry Applications (LDIA), Wuhan, China, 1–3 July 2021; IEEE: Piscataway, NJ, USA, 2021; pp. 1–6. [Google Scholar]
  16. Benavides, R.; Mutschler, P. Detent force compensation in segmented long stator permanent magnet linear drives using Finite Element Models. In Proceedings of the European Conference on Power Electronics & Applications, Aalborg, Denmark, 2–5 September 2007. [Google Scholar]
  17. Cui, F.; Sun, Z.; Xu, W.; Qian, H.; Cao, C. Optimization Analysis of Long Primary Permanent Magnet Linear Synchronous Motor. IEEE Trans. Appl. Supercond. 2021, 31, 1–4. [Google Scholar] [CrossRef]
  18. Heidary, M.; Naderi, P.; Shiri, A. Modeling and analysis of a multi-segmented linear permanent-magnet synchronous machine using a parametric magnetic equivalent circuit. Electr. Eng. 2022, 104, 705–715. [Google Scholar] [CrossRef]
  19. Bianchi, N. Analytical field computation of a tubular permanent-magnet linear motor. IEEE Trans. Magn. 2000, 36, 3798–3801. [Google Scholar] [CrossRef]
  20. Yudin, M.; Sevostyanov, N.; Tsareva, D. Simulation of Electromagnetic Fields in Layered Media Based on Semi-Analytical Integral Transformations. Conf. Proc. 2021, 2021, 1–6. [Google Scholar]
  21. Huangfu, Y.; Wang, S.; Qiu, J.; Zhang, H.; Wang, G.; Zhu, J. Transient Performance Analysis of Induction Motor Using Field-Circuit Coupled Finite-Element Method. IEEE Trans. Magn. 2014, 50, 873–876. [Google Scholar] [CrossRef]
  22. Maroufian, S.S. Winding Function Modeling of Synchronous Reluctance Machines. Ph.D. Thesis, Concordia University, Montreal, QC, Canada, 2018. [Google Scholar]
  23. Lakshmi, N.S.; Allirani, S. Modelling and Simulation of Permanent Magnet Synchronous Motor for Performance Enhancement Using ANSYS Maxwell; Springer: Berlin/Heidelberg, Germany, 2021. [Google Scholar]
  24. Liu, Z.; Li, J. Accurate Prediction of Magnetic Field and Magnetic Forces in Permanent Magnet Motors Using an Analytical Solution. IEEE Trans. Energy Convers. 2008, 23, 717–726. [Google Scholar] [CrossRef]
  25. Hu, H.; Zhao, J.; Liu, X.; Guo, Y.; Zhu, J. No-Load Magnetic Field and Cogging Force Calculation in Linear Permanent-Magnet Synchronous Machines With Semiclosed Slots. IEEE Trans. Ind. Electron. 2017, 64, 5564–5575. [Google Scholar] [CrossRef] [Green Version]
  26. Lu, Q.; Wu, B.; Yao, Y.; Shen, Y.; Jiang, Q. Analytical model of permanent magnet linear synchronous machines considering end effect and slotting effect. IEEE Trans. Energy Convers. 2019, 35, 139–148. [Google Scholar] [CrossRef]
  27. Chen, F.; Zhang, C.; Chen, J.; Yang, G. Accurate Subdomain Model for Computing Magnetic Field of Short Moving-Magnet Linear Motor With Halbach Array. IEEE Trans. Magn. 2020, 56, 8200509. [Google Scholar] [CrossRef]
  28. Markovic, M.; Jufer, M.; Perriard, Y. Reducing the cogging torque in brushless DC motors by using conformal mappings. IEEE Trans. Magn. 2004, 40, 451–455. [Google Scholar] [CrossRef]
  29. O’Connell, T.C.; Krein, P.T. A Schwarz–Christoffel-Based Analytical Method for Electric Machine Field Analysis. IEEE Trans. Energy Convers. 2009, 24, 565–577. [Google Scholar] [CrossRef]
  30. Liyi, L.; Hong, J.; Lu, Z.; Ying, L.; Song, Y.; Rizhong, L.; Xiaopeng, L. Fields and Inductances of the Sectioned Permanent-Magnet Synchronous Linear Machine Used in the EMALS. IEEE Trans. Plasma Sci. 2011, 39, 87–93. [Google Scholar] [CrossRef]
  31. Tursini, M.; Villani, M.; Tullio, A.D.; Fabri, G.; Collazzo, F.P. Nonlinear Model Suitable for the Offline Cosimulation of Fault-Tolerant PM Motors Drives. IEEE Trans. Ind. Appl. 2017, 53, 3719–3729. [Google Scholar] [CrossRef]
  32. Benavides, R.; Mutschler, P. Compensation of Disturbances in Segmented Long Stator Linear Drives using Finite Element Models. In Proceedings of the International Symposium on Industrial Electronics, Montreal, QC, Canada, 9–13 July 2006. [Google Scholar]
  33. Mihalachi, M.; Leidhold, R.; Mutschler, P. Linear drive system for combined transportation and processing of materials. In Proceedings of the Conference of the Industrial Electronics Society, Porto, Portugal, 3–5 November 2009. [Google Scholar]
  34. Gysen, B.L.J.; Meessen, K.J.; Paulides, J.J.H.; Lomonova, E.A. General Formulation of the Electromagnetic Field Distribution in Machines and Devices Using Fourier Analysis. IEEE Trans. Magn. 2010, 46, 39–52. [Google Scholar] [CrossRef] [Green Version]
  35. Liu, Z.; Li, J. Analytical Solution of Air-Gap Field in Permanent-Magnet Motors Taking Into Account the Effect of Pole Transition Over Slots. IEEE Trans. Magn. 2007, 43, 3872–3883. [Google Scholar] [CrossRef]
  36. Mao, Y.; Sun, Z.; Zhou, W.; Cui, F. Equivalent Circuit Model of Long Primary Permanent Magnet Linear Synchronous Motor Considering End Effect. In Proceedings of the 2021 13th International Symposium on Linear Drives for Industry Applications (LDIA), Wuhan, China, 1–3 July 2021. [Google Scholar]
  37. Virtic, P.; Stumberger, B. Analytical Analysis of Magnetic Field and Force Calculation in a Slotless-Type Permanent Magnet Linear Synchronous Machine; Verification with Numerical Analysis. In Proceedings of the International Electric Machines and Drives Conference, Antalya, Turkey, 3–5 May 2007. [Google Scholar]
Figure 1. WS-PMLSM for transmission lines.
Figure 1. WS-PMLSM for transmission lines.
Symmetry 14 01049 g001
Figure 2. Structure of the WS-PMLSM.
Figure 2. Structure of the WS-PMLSM.
Symmetry 14 01049 g002
Figure 3. WS-PMLSM exact sub-domain model.
Figure 3. WS-PMLSM exact sub-domain model.
Symmetry 14 01049 g003
Figure 4. Magnetization intensity distribution function.
Figure 4. Magnetization intensity distribution function.
Symmetry 14 01049 g004
Figure 5. Primary side end effect model diagram.
Figure 5. Primary side end effect model diagram.
Symmetry 14 01049 g005
Figure 6. Coordinate plane of z w t .
Figure 6. Coordinate plane of z w t .
Symmetry 14 01049 g006
Figure 7. End effect of the air gap ratio permeability.
Figure 7. End effect of the air gap ratio permeability.
Symmetry 14 01049 g007
Figure 8. Analysis model of PM secondary end effects.
Figure 8. Analysis model of PM secondary end effects.
Symmetry 14 01049 g008
Figure 9. Experimental prototype of the proposed system.
Figure 9. Experimental prototype of the proposed system.
Symmetry 14 01049 g009
Figure 10. Flux simulation of WS-PMLSM using JMAG.
Figure 10. Flux simulation of WS-PMLSM using JMAG.
Symmetry 14 01049 g010
Figure 11. Permanent magnet surface flux density distribution.
Figure 11. Permanent magnet surface flux density distribution.
Symmetry 14 01049 g011
Figure 12. Air gap center flux density distribution.
Figure 12. Air gap center flux density distribution.
Symmetry 14 01049 g012
Figure 13. Slot surface flux density distribution.
Figure 13. Slot surface flux density distribution.
Symmetry 14 01049 g013
Figure 14. PM end density distribution.
Figure 14. PM end density distribution.
Symmetry 14 01049 g014
Figure 15. Flux with symmetrical voltage excited.
Figure 15. Flux with symmetrical voltage excited.
Symmetry 14 01049 g015
Figure 16. Detent force characteristics.
Figure 16. Detent force characteristics.
Symmetry 14 01049 g016
Table 1. Structure parameters of PMSLM.
Table 1. Structure parameters of PMSLM.
ParametersSymbolValue (Unit)
Length of the stator L s t a t o r 200 mm
Length of the mover L m o v e r 200 mm
Air-gap length δ 1.5 mm
Height of slot opening H s 15 mm
Height of PM H m 4.0 mm
Height of the back core of PM H z c 4.0 mm
Height of the stator core H s t 25 mm
Stator tooth width b s 4 mm
Slot pitch t s 8 mm
Pole pitch τ 50 mm
Relative permeability of PM μ r 1.04
PM remanence B r 1.12 T
Number of armature teeth Z p 25
Number of pole pairsp1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, T.; Mei, X. Research on Detent Force Characteristics of Winding Segmented Permanent Magnet Linear Synchronous Motor Based on Analytical Model. Symmetry 2022, 14, 1049. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14051049

AMA Style

Zhang T, Mei X. Research on Detent Force Characteristics of Winding Segmented Permanent Magnet Linear Synchronous Motor Based on Analytical Model. Symmetry. 2022; 14(5):1049. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14051049

Chicago/Turabian Style

Zhang, Tuanshan, and Xuesong Mei. 2022. "Research on Detent Force Characteristics of Winding Segmented Permanent Magnet Linear Synchronous Motor Based on Analytical Model" Symmetry 14, no. 5: 1049. https://0-doi-org.brum.beds.ac.uk/10.3390/sym14051049

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