Next Article in Journal
Influence of Fluid Viscosity on Cavitation Characteristics of a Helico-Axial Multiphase Pump (HAMP)
Previous Article in Journal
Three Terminal Perovskite/Silicon Solar Cell with Bipolar Transistor Architecture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Speed Control of Magnetic Drive-Trains with Pole-Slipping Amelioration

1
Guangdong University of Petrochemical Technology (GDUPT), Maoming 525000, China
2
The Joint Research Lab between GDUPT and University of Lincoln, Lincoln LN6 7TS, UK
3
School of Engineering, University of Lincoln, Lincoln LN6 7TS, UK
*
Authors to whom correspondence should be addressed.
Energies 2022, 15(21), 8148; https://doi.org/10.3390/en15218148
Submission received: 22 September 2022 / Revised: 26 October 2022 / Accepted: 28 October 2022 / Published: 1 November 2022
(This article belongs to the Topic Advanced Electrical Machines and Drives Technologies)

Abstract

:
The paper introduces new techniques to reduce the potential for pole-slipping induced by control systems and presents a low-cost pole-slipping detection and recovery scheme for magnetic drive-trains (MDTs). For the first time, the paper shows that a combination of electromagnetic and load-torque excitations which individually are not greater than the maximum coupling torque can initiate pole-slipping. For applications where acceleration feedback is unavailable, the motor-side inertia is virtually increased with a tracking differentiator to provide feedback of acceleration. Subsequently, controller design and parameter optimization are discussed. Experimental measurements on a custom test facility verify the presented principles that low-bandwidth controller designs with low inertia ratios can accommodate a wider range of on-load startup torque and load-torque disturbances without pole-slipping. To address overload issues, a pole-slipping detection method based on the kurtosis of electromagnetic torque and a recovery strategy based on converting the state of pole-slipping into that of on-load startup are presented. Experimental results demonstrate that detecting slip anomalies without load-side information, and recovery from pole-slipping without auxiliary mechanical devices are both feasible.

1. Introduction

Due to the physical isolation of motor- and load-side shafts when employing magnetic couplings or magnetic gears, magnetic drive-trains (MDTs) can overcome the inherent problems of traditional couplings/gearboxes, such as overload failures, jamming, mechanic fatigue, and regular maintenance requirements [1]. Since the first modern magnetic gear solution was presented in 2001 [2], much published work has discussed various topologies, computational models, and design aspects to improve torque densities and gearing ratios for magnetic-torque transmission devices [3,4,5]. However, studies around the high-performance control of MDTs, with pole-slipping amelioration, remain limited.
MDTs, which present with variable nonlinear stiffness, can be regarded as a variant of two-mass systems. Montague [1] linearized the model of MDTs at the no-load condition to derive a constant torsional stiffness; MDTs are then modelled as classic two-mass systems. Recently reported work [6] shows that linearsing the MDT model at 75% of the maximum coupling torque is more appropriate. Using the linearized model, Refs. [1,6] adopt basic pseudo derivative feedback (PDF) and 2-degree-of-freedom PI as the control strategy, respectively. By contrast, Refs. [7,8,9] employ observers to achieve a full-state feedback control of MDTs. Among these studies, only [8] introduces a proportional gain to reduce the probability of controller-induced slip; however, as in [1,7,9], specifications such as the maximum starting load-torque or load-torque disturbance that designed controllers can accommodate, have not been specified in [8].
Pole-slipping detection and recovery are crucial for the control design of MDTs. Montague [1,10] observed that pole-slipping imposes a modulation onto motor-side speed and that pole-slipping can then be observed through the speed oscillations of motor-side components. However, [1,10] do not give a detection criterion due to the amplitude of oscillations being affected by inertia ratios and design coupling torque. Ref. [11] calculates the relative mechanical displacement angle to identify pole-slipping, but it requires sensors to be installed on both motor- and load-side shafts. If load-side information is obtained through estimations, high-resolution sensors (to obtain motor-side information) and high-performance computational units are required to be used.
Ref. [1] proposes a slip recovery strategy based on reconfiguring PDF controllers. The pole-slipping recovery approach presented in [11] requires measuring or observing the load-side position of MDTs. The pole-slipping recovery experimental results given by [1,11] show that the speed of the load-side is almost zero during the overload pole-slipping. This means the MDTs often require extra mechanic devices to isolate the shafts during overload conditions and then re-couple them when recovery begins.
If speed oscillations can be suppressed, the possibility for pole-slipping can be reduced. Methodologies [12,13,14] adopted to inherently reduce speed oscillations of two-inertia systems indicate that effective oscillation suppression can be achieved by adjusting the inertia ratio with additional acceleration or torsional torque feedback loops. In contrast, speed oscillation attenuation with inertia ratio regulation has yet to be studied in the control methods proposed for MDTs.
Moreover, previously reported work [1,6,11,15] has only considered pole-slipping in response to overload or aggressive control effort induced by references or load-torque disturbances. However, load-startup failures caused by the mismatch of load- and motor-side inertias have yet to be discussed. The primary aim of this paper is therefore to provide a mathematical analysis for pole-slipping and then use the developed principles to improve the operational robustness of MDTs. The benefits presented in this paper include.
(1)
Through a formal mathematical derivation, the paper firstly highlights that a MDT can be led into a pole-slipping regime with electromagnetic and load-torque excitations that are not individually greater than the maximum coupling torque. Moreover, the analysis also shows that the robustness of MDTs to actuation induced pole-slipping can be improved by decreasing the inertia ratio and reducing the bandwidth of designed controllers.
(2)
The kurtosis of electromagnetic torque is introduced to provide a pole-slipping detection criterion that doesn’t rely on the information of load-side position/speed. Moreover, a novel pole-slipping recovery is presented by transforming the pole-slipping regime into that of an on-load startup problem. The pole-slipping remedial strategy ensures MDTs recover from failures without the need for additional mechanical intervention/devices, as have been required in previous studies, to facilitate recovery.

2. Dynamics of MDTs

2.1. Experimental Rig

A magnetically-coupled drive-train test facility has been setup as shown in Figure 1 to provide a practical focus to the investigation. The motor-side position sensor is a 12 bit incremental encoder. The magnetic coupling has 5 pole pairs (p), and the maximum (pull-out) coupling torque ( T G ) is 1.6 N·m. The test rig has the load- and motor-side friction coefficients ( B L and B M ) of 0.003 N·m/(rad/s) and the load- and motor-side inertias ( J L and J M ) are 0.001 kg · m 2 .
For the driving unit, the motor-side speed ω M and position θ M can be measured, but the load-side speed ω L and position θ L are assumed unavailable for control purposes. This test facility emulates the most common application scenarios for MDTs, where ω L and θ L are ideally not measured. Control effort is electromagnetic torque ( T e m ), and hence the prime mover is operating in a torque control mode. The driven unit is used to emulate load ( T L ) variations.

2.2. Dynamic Model

Defining T C as the torsional torque developed between the two magnetically-coupled shafts, T C is then given by
T C = T G sin p θ D
where θ D = θ M θ L is the relative mechanical displacement angle.
If damping torque is not taken into consideration, MDTs are described by the following set of equations [1].
θ ˙ M = ω M ω ˙ M = 1 J M T e m T C B M ω M θ ˙ L = ω L ω ˙ L = 1 J L T C T L B L ω L .
θ M , ω M , θ L , ω L are considered as state variables, T e m and T L as inputs, and θ D as the output. Linearizing (2) at an operating point, and letting K lin represent the resulting torsional stiffness K lin = T G p cos p θ D , the transfer functions from T e m and T L to θ D are then, respectively, given by
θ D ( s ) T e m ( s ) = J L s + B L D ( s ) ,
θ D ( s ) T L ( s ) = J M s + B M D ( s )
where D ( s ) = J M J L s 3 + J M B L + B M J L s 2 + B L B M + ( J L + J M ) K lin s + ( B M + B L ) K lin .
The anti-resonance frequency ( ω a ) and natural resonance frequency ( ω n ) are then given by
ω a = K lin J L , ω n = K lin J L ( 1 + R )
where R = J L / J M . Table 1 enumerates ω a and ω n , which ranges from no-load to the 99% of the pull-out torque.
If B L B M and ( B M + B L ) K lin are negligible, then (3) and (4) simplify to
θ D ( s ) T e m ( s ) = 1 / J M s 2 + 2 ζ ω n s + ω n 2 ,
θ D ( s ) T L ( s ) = 1 / J L s 2 + 2 ζ ω n s + ω n 2
where ζ = J M B L + B M J L 2 J M J L ω n .

2.3. On-Load Startup and Criteria for MDTs to Maintain Normal Running

On-load startup is the most demanding operating condition that can result in actuation induced pole-slipping, where the relative angle of both shafts is simultaneously developed by the transients of electromagnetic torque demands and load-torque.
If T e m and T L are considered as step inputs ( T e m ( s ) = T e m / s , and T L ( s ) = T L / s ), then θ D resulting from T e m and T L can be, respectively, rewritten based on (6) and (7), and are given by
θ D ( s ) | T e m = T e m / J M s s 2 + 2 ζ ω n s + ω n 2 ,
θ D ( s ) | T L = T L / J L s s 2 + 2 ζ ω n s + ω n 2 .
Inverse Laplace transforming of (8) and (9) yields
θ D ( t ) | T e m = T e m ω n 2 J M 1 ( 1 + ζ ω n ω d ) e ζ ω n t sin ω d t ,
θ D ( t ) | T L = T L ω n 2 J L 1 ( 1 + ζ ω n ω d ) e ζ ω n t sin ω d t
where ω d = 1 ζ 2 ω n . Since ζ is small, ω n ω d . Referring to (10) and (11), the maximum increment of θ D results from when T e m and T L occur simultaneously and approximates to
θ D m = T e m R + T L ω n 2 J L 1 + e 3 π 2 ζ .
It can be seen from (12) that increasing J M improves the on-load startup capability provided that J L is fixed ( R = J L / J M ). Moreover, if T e m output exhibits first-order dynamics compared to T L which is usually assumed to be a step input, θ D m would take a much lower value. This also means that low-bandwidth controllers can potentially accommodate a higher on-load startup torque and load-torque disturbances than high-bandwidth designs.
A MDT will enter a pole-slip operating regime if (13) or (14) holds.
| θ D | > π 2 p ,
T e m > T G + T f m , or T L > T G + T f l
where T f m and T f l represent motor- and load-side friction torque. Now, assuming θ D 0 is the relative angle of both shafts while MDTs are running under normal conditions, then a criterion for a MDT to avoid pole-slipping is given by
θ D 0 + θ D m < π 2 p .
The experimental study of [6] shows that selecting 75% of the pull-out torque as the operating point can accommodate a wider range of load-torque disturbance. Here, we will show how to use Equations (12), (13), and (15) to estimate the maximum on-load startup torque and the permitted load-torque disturbance that a MDT can accommodate through two example scenarios.
Example 1: supposing the MDT rig is operating at 800 r/min under no-load conditions ( θ D 0 = arcsin ( 0.18 ) / p ), 75% of the pull-out torque is selected as the operating point ( ω n = 103.0 rad/s), and the increment of T e m can immediately follows T L to find the maximum load-torque disturbance that the MTD rig accommodate without pole-slipping.
According to (12) and (15), we have
T L < π 2 p θ D 0 4 ω n 2 J L .
Substituting the values of θ D 0 , ω n , and J L into (16) yields T L = 0.43 T G . However, due to the delay induced by current/speed control loops, the permitted load-torque disturbance can conceivably be greater than 0.43 T G .
Example 2: Assuming 75% of the pull-out torque is selected as the operating point, and T e m is equal to T G , find the maximum on-load starting torque ( T L m ) that the facility can accommodate without pole-slipping.
According to (12) and (13), T L m for the rig discussed in the paper is 0.66 T G . If friction torque is taking into consideration, then, T L m = 0.48 T G .

3. Controller Design and Pole-Slipping Amelioration

3.1. Virtually Increasing J M Using Acceleration Feedback

From the analysis above, we know that increasing J M (effectively decreasing R, see (12)) decreases the possibility of pole-slipping related to on-load startup. Figure 2 shows a model of the MDT with motor-side acceleration feedback ( K a is the feedback coefficient). Clearly, with acceleration feedback, the equivalent motor-side inertia J ^ M is given by
J ^ M = J M + K a .
Accordingly, the inertia ratio is adjusted to
R = J L J ^ M .
Letting T e m represent the actual torque developed by the drive, then
T e m = T e m K a ω ˙ M .
Equation (19) indicates the aggressive control effort T e m can be modified by K a , and this mechanism can be used to prevent θ D exceeding its constraint under on-load startup and inertia mismatch.
The drive for a practical facility usually integrates period or frequency measurement methods to calculate speed signals with feedback from encoders. However, most drives do not provide acceleration information. Calculating acceleration from speed measurements at a particular sampling frequency will inevitably amplify noise components.
A tracking differentiator (TD) [16] can effectively suppress the amplification of derivative noise while quickly tracking the derivative of input signals without relying on model information. Here, the TD is adopted to obtain an estimation of acceleration.
Supposing v 1 tracks the motor-side speed ω M , v 2 is the derivative of v 1 , then the TD is described by
v 1 ˙ = v 2 v 2 ˙ = fst ( v 1 , v 2 , ω M , r , h )
where r dictates tracking speed, and h is a parameter that determines the filtering effect ( h T s , T s is the sampling period), and fst ( v 1 , v 2 , ω M , r , h ) is given by
fst = r a , a d r sign ( a ) , a > d
where d = r h , and a is given by
a = v 2 + y / h , y d 0 v 2 + 1 2 ( a 0 d ) sign ( y ) , y > d 0
with d 0 = d h , y = v 1 ω M + h v 2 , a 0 = d 2 + 8 r y .

3.2. Controller Design and Parameter Optimization

The paper now describes the design of a generalized controller (fractional-order controller) represented by
F C ( s ) = K p + K i s λ + K d s μ
where K p , K i , K d represent the proportional, integral, and derivative gains and λ , μ are the order of integrals and derivative items. Of note is that traditional integer PID controllers can be regarded as fractional-order PID controllers (FOPID) with λ = μ = 1 . Hence, the following design process is also applicable for the design of integer-order controllers. Figure 3 shows the principle of optimizing parameters for speed controller using integral of time and absolute error as a performance metric. Specifically, the objective function f ( x , t ) is considered:
f ( x , t ) = 0 t | e ( t ) | d t
where x = K p , K i , K d , λ , μ and e ( t ) is the error of the control loop at time t.
The purpose of parameter optimization is to find an x subject to constraints to ensure f ( x , t ) has a minimum value, described by
min x f ( x , t ) such that c ( x ) < 0 l b x u b
where c ( x ) represents general nonlinear inequality constraints; and l b and u b are the lower- and upper-boundaries for x, respectively. Normally, c ( x ) can be constructed with a combined specifications of gain crossover frequency, phase margin, gain margin, sensitivity function, and complementary sensitivity function.
From the design perspective of two-mass systems, the gain cross-over frequency ( ω g c ) should be set much smaller than ω a to obtain responses similar to that of a low-pass filter, i.e. ideally with no overshoot. It can be seen from Table 1 that ω a = 59.1 rad/s while the MDT rig is running at 90% of the pull-out load-torque. Moreover, the analysis of Section 2.3 shows that a low-bandwidth design can accommodate a wider range of load-torque disturbances. Here, ω g c is set to 5 rad/s.
The sensitivity function S j ω represents the transfer function from the disturbance (d) to ω M and the speed reference (r) to e ( S j ω = 1 / 1 + F C ( s ) G T e m ω M ( s ) , see Figure 3 ) . Equations (10) and (11) show that θ D under excitations takes the form of a step plus a damped oscillation which is regarded as disturbances in this design. To obtain a small control error, | S j ω | is also required to be set small at low frequencies since references usually take the form of relatively low-frequency inputs. Hence, the sensitivity specification for controller design can therefore be given by
ω s = 5 rad / s , | S j ω | 20 dB , ω ω s .
Note: here, ω s is assumed to be a value approximately equal to ω g c .
The complementary sensitivity function T j ω relates speed measurement noise to motor-side speed ω M . Typically, the cut-off frequency of shaft speed filters is usually designed to be 3∼5 ω n ( ω n = 126 rad/s while the MDT rig is running at no-load conditions, see Table 1) to meet dynamic processing requirements.The specifications for the complementary sensitivity function are then given by
ω t = 5 ω n , | T j ω | 20 dB , ω ω t .
Previously reported work [17] shows that satisfactory performance can be obtained by setting the gain margin G M > 6 dB and the phase margin 30 0 P M 60 0 . Here, P M and G M are chosen to be
PM = 45 0 , GM = 10 dB .
The parameters for the TD are r = 500 and h = 0.004 (the sampling frequency for running speed signals is set to 250 Hz). In this case, K a is set to be 0.1 (see Figure 2). The optimization is carried out using the open-source software FOMCON [18], and the optimized parameters for the design controller are K p = 3.97 , K i = 3.17 , K d = 3.8 × 10 5 , λ = 1 , and μ = 0.12 . The optimized controller can be further approximated to a classical PI structure as 3.97 s + 3.17 / s .

3.3. A Criterion for Pole-Slipping Detection

Here, an overload pole-slipping detection methodology is considered. MDTs will enter a pole-slipping regime when the relative mechanical displacement angle between the magnetically-coupled shafts exceeds the designed range. During pole-slipping, ω M will rise for a short period and then return to the reference speed since the speed controller remains in operation. The control effort T e m will therefore converge to be roughly equal to the friction torque in steady state ( T e m T f m ).
The abrupt change of T e m can be measured by kurtosis, which is a dimensionless indicator widely used for evaluating the characteristics of vibrations. The kurtosis of T e m is derived from [19]
K u r = 1 n i = 1 n x i x ¯ 4 1 n i = 1 n x i x ¯ 2 2
where x i represents the control effort T e m at time i, and x ¯ is the mean value of x i .
A pole-slipping detection criteria can then be described as when
T e m T f m , K u r > 5 .

3.4. A Strategy for Pole-Slipping Recovery

It can be seen from the paragraph below Equation (2) that the equivalent torsional stiffness ( K l i n ) varies with the displacement angle ( θ D ) between the magnetically-coupled shafts. Once the load-side shaft synchronizes with the motor-side, the displacement angle is constant, resulting in a stable coupling torque ( T c = K l i n θ D ). Hence, the variance of ω M during synchronization is much smaller than the commensurate value when the motor- and load-side speeds are different. The proposed pole-slipping recovery strategy considered here is based on the assumption that T e m can be greater than T G , and T L can restore from the state T L > T G to T L < T G . The proposed recovery strategy is summarised as follows:
  • At the onset of pole-slipping, reverse the running direction of the motor-side shaft and control the speed of it up to the maximum specification speed, then;
  • Control the electromagnetic torque to be zero ( T e m = 0 ); | ω M | will consequently then decrease;
  • As | ω M | decreases, the motor-side shaft will re-synchronize with the load-side shaft (the variance of ω M is now used as a synchronization detection indicator). The PI plus TD configuration can then be reconfigured to control the MDT to the reference speed.

4. Simulation and Experimental Results

To verify the theoretical analysis in Section 2.3, this Section will initially consider three cases of on-load startup (first case: T L = 0.4 T G , second case: T L = 0.5 T G , third case: T L = 0.6 T G ) through simulation. For the discussed cases, the electromagnetic torque is assumed to be saturated immediately ( T e m = T G ). Note that in this paper, T e m and T L are normalized at T G for both simulation and experimental trials.
Figure 4a shows the simulation results for R = 1 . It can be seen that slip between two shafts of the MDT occurs while T L is 0.5 T G or above. Although increasing J M results in a slower transient response, Figure 4b shows that the on-load starting torque can be up to 0.5 T G by modifying J M to 2 J L . However, as shown in Figure 4c, by setting J M to 0.5 J L , the two sides of the MDT are out of synchronization in all studied cases. It is clear that decreasing the inertia ratio through increasing J M improves on-load startup capability. In practice, low-cost drives often cannot provide sufficient bandwidth for their current control loops. Hence, the output electromagnetic torque and load-toque are assumed to exhibit first-order dynamics. Under these circumstances, from Equation (12) and subsequent discussion, the permitted on-load startup torque can be a higher value compared with the simulation studies.
Figure 5 shows tshows the measured load rejection dynamics for the optimized controller. As shown in Figure 5a, the speed of load-side coupling follows the motor-side without speed oscillations even when the MDT rig is subject to 2 significant load disturbances. Figure 5b shows that by adjusting the inertia ratio with K a = 0.1 offers first-order filter dynamics for the reactive control effort T e m compared to the step load-torque disturbance occurring at 13 s. This characteristic greatly improves the range of load-torque disturbance and on-load startup torque that the MDT rig can accommodate while referring to (10), (11), and (12). It can also be seen from Figure 5b that the peak control effort resulting from the speed reference is about 0.7 T G . Assuming K a = 0 (remove the acceleration feedback) and T e m can immediately follow T L without overshoot, the maximum on-load startup torque would be around 0.48 T G if friction torque is included (as shown in Figure 5b, T f m = T f l 0.18 T G at around 800 r/min).
The on-load startup torque is now set to 0.63 T G . Figure 6 shows that virtually increasing J M through acceleration feedback overcomes the pole-slipping problem associated with on-load startup and that a starting load-torque that is greater than 0.63 T G can lead the MDT to pole-slipping ( T e m slightly exceeds T G at around 1.8 s; however, the effective coupling torque is still a slightly lower than T G due to the existence of motor-side friction). Of note is that if load-side friction torque ( T f l ) is considered, the maximum on-load starting torque is 0.81 T G in this case.
It can also be seen from the experimental trials that the inertia ratio is set to 0.09 ( K a = 0.1 ) to prevent pole-slipping resulting from on-load startup and inertia mismatch. By contrast, previously reported work [12,13,20] shows that effective oscillation suppression is achieved by setting the resonance ratio ( ω n / ω a ) of two-mass systems to be 2 5 ; thus, R ranges between 1 and 4 (see (5)). Although MDTsare traditionally regarded as two-mass systems, it can be seen from the simulations and experimental trials that the recommended inertia ratio for MDTs should generally be set much lower than that of classical two-inertia systems to prevent MDTs from entering pole-slipping initiated by the combined effects of T e m and T L .
The efficacy of pole-slipping detection and recovery strategies are now investigated. Figure 7b shows that the kurtosis of T e m is lower than three, when the MDT is operating in normal conditions. However, the kurtosis rises dramatically to 8.71 during the period of 10.42 s to 12.64 s (pole-slipping occurs at 12.1 s). It can also be seen from Figure 7a that T e m T f m at around 12.5 s. Hence, the proposed method can effectively detect the transition to pole-slipping.
The experimental results shown in Figure 8 provide an example of the recovery strategy. As shown in the figure, pole-slipping is induced at 12 s. At the beginning of pole-slipping, the load-side shaft reverses its direction due to T L being temporarily higher than the maximum coupling torque and T e m . Following that, the load-side shaft loses synchronization with the primary-side, running at around 1600 r/min during pole-slipping. After detecting pole-slipping with the kurtosis methodology discussed above, the motor-side drive reverses its running direction from 14 s and controls ω M to 2000 r/min, a setting higher than ω L . It can be seen from Figure 8 that the motor-side speed varies around 2000 r/min from at 16 s and that the variance of ω M increases considerably when both side speeds are different. From 25 s to 30 s, T e m is set to 0, | ω M | then decreases. Following that, the motor-side shaft synchronizes with the load-side shaft from roughly 27 s. Now, the coupling torque is relative constant; therefore, speed oscillations are reduced. This is a situation that can be regarded as on-load startup. The reconfigured scheme can then be used to recover the MDT rig from pole-slipping. It can be seen from the figure that the MDT recovers after 30 s. According to Table 2, the variance of ω M during the period of 19–21 s ( ω L = 1600 r/min, ω M = 2000 r/min) is far higher than the variance calculated from 9 s–11 s or 27 s–29 s ( ω L = ω M ). Hence, the variance of ω M is a suitable indicator for synchronisation.

5. Conclusions

In this paper, pole-slipping anomalies induced by electromagnetic inputs and load-torque disturbances are investigated. Experimental trials have been used to verify that a low-bandwidth controller design incorporating a tracking differentiator to virtually increase motor-side inertia significantly improves the on-load startup performance and load-torque disturbance handling range for the MDT where acceleration feedback is unavailable. MDTs are traditionally classified as two-inertia systems; however, the experimental trials show that an appropriate inertia ratio for MDTs to maintain normal operation is much lower than that of traditional dual-inertia control systems. Experimental measurements have confirmed that the proposed amelioration methods can effectively detect pole-slipping without using load-side position/speed measurements and can overcome the requirement for including additional mechanic devices to facilitate recovery.

Author Contributions

Conceptualization, X.L., C.B. and T.S.; formal analysis, X.L.; project administration, C.B.; software, X.L.; supervision, C.B.; experiment investigation, X.L., C.B. and T.S.; writing—original draft, X.L.; writing—review and editing, C.B.. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Argyrios Zolotas for inspiring to employ fractional-order control as the control scheme.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ω a    Anti-resonance frequency (rad/s)
ω d    The frequency of damped sine signal (rad/s)
ω L    Load-side angular velocity (rad/s)
ω M    Motor-side angular velocity (rad/s)
ω n    Natural resonance frequency (rad/s)
θ L    Load-side angle (rad)
θ M    Motor-side angle (rad)
θ D     θ M θ L (rad)
ζ    Damping coefficient (kg/s), and  ζ = J M B L + B M J L 2 J M J L ω n
λ    Integer order
μ    Derivative order
B L    Load-side friction coefficient N·m/(rad/s)
B M    Motor-side friction coefficient N·m/(rad/s)
J L    Load-side inertia ( kg · m 2 )
J M    Motor-side inertia ( kg · m 2 )
J ^ M    The equivalent value of motor-side inertia ( kg · m 2 )
K a    Feedback coefficient
K d    Derivative gain
K i    Integral gain
K p    Proportional gain
K lin    Linearized torsional stiffness ( N · m/rad)
p   Pole-pairs
R   Inertia ratio ( J L / J M )
T C    Coupling torque ( N · m)
T e m    Electromagnetic torque ( N · m)
T f l    Load-side friction torque ( N · m)
T f m    Motor-side friction torque ( N · m)
T L    Load torque ( N · m)
T G    The maximum coupling torque ( N · m)

References

  1. Montague, R.; Bingham, C.; Atallah, K. Servo control of magnetic gears. IEEE/Asme Trans. Mechatron. 2011, 17, 269–278. [Google Scholar] [CrossRef]
  2. Atallah, K.; Howe, D. A novel high-performance magnetic gear. IEEE Trans. Magn. 2001, 37, 2844–2846. [Google Scholar] [CrossRef] [Green Version]
  3. Wang, Y.; Filippini, M.; Bianchi, N.; Alotto, P. A Review on Magnetic Gears: Topologies, Computational Models and Design Aspects. IEEE Trans. Ind. Appl. 2019, 55, 4557–4566. [Google Scholar] [CrossRef]
  4. Rasmussen, P.O.; Andersen, T.O.; Jorgensen, F.T.; Nielsen, O. Development of a high-performance magnetic gear. IEEE Trans. Ind. Appl. 2005, 41, 764–770. [Google Scholar] [CrossRef]
  5. Jorgensen, F.T.; Andersen, T.O.; Rasmussen, P.O. The Cycloid Permanent Magnetic Gear. IEEE Trans. Ind. Appl. 2008, 44, 1659–1665. [Google Scholar] [CrossRef]
  6. Liao, X.; Bingham, C.; Zolotas, A.; Zhang, Q.; Smith, T. Servo Control of Drive-Trains Incorporating Magnetic Couplings. IEEE Trans. Ind. Appl. 2022, 58, 3674–3684. [Google Scholar] [CrossRef]
  7. Montague, R.; Bingham, C.; Atallah, K. Dual-observer-based position-servo control of a magnetic gear. IET Electr. Power Appl. 2011, 5, 708–714. [Google Scholar] [CrossRef]
  8. Bouheraoua, M.; Wang, J.; Atallah, K. Rotor position estimation of a pseudo direct-drive pm machine using extended kalman filter. IEEE Trans. Ind. Appl. 2016, 53, 1088–1095. [Google Scholar] [CrossRef] [Green Version]
  9. Bouheraoua, M.; Wang, J.; Atallah, K. Design and implementation of an observer-based state feedback controller for a pseudo direct drive. IET Electr. Power Appl. 2013, 7, 643–653. [Google Scholar] [CrossRef]
  10. Montague, R.; Bingham, C.; Atallah, K. Magnetic gear dynamics for servo control. In Proceedings of the MELECON 2010–2010 15th IEEE Mediterranean Electrotechnical Conference, Valletta, Malta, 26–28 April 2010; pp. 1192–1197. [Google Scholar]
  11. Bouheraoua, M.; Wang, J.; Atallah, K. Slip recovery and prevention in pseudo direct drive permanent-magnet machines. IEEE Trans. Ind. Appl. 2014, 51, 2291–2299. [Google Scholar] [CrossRef]
  12. O’Sullivan, T.M.; Bingham, C.M.; Schofield, N. High-performance control of dual-inertia servo-drive systems using low-cost integrated SAW torque transducers. IEEE Trans. Ind. Electron. 2006, 53, 1226–1237. [Google Scholar] [CrossRef] [Green Version]
  13. Saarakkala, S.E.; Hinkkanen, M. State-space speed control of two-mass mechanical systems: Analytical tuning and experimental evaluation. IEEE Trans. Ind. Appl. 2014, 50, 3428–3437. [Google Scholar] [CrossRef] [Green Version]
  14. Szabat, K.; Orlowska-Kowalska, T. Vibration Suppression in a Two-Mass Drive System Using PI Speed Controller and Additional Feedbacks—Comparative Study. IEEE Trans. Ind. Electron. 2007, 54, 1193–1206. [Google Scholar] [CrossRef]
  15. Montague, R.G.; Bingham, C.; Atallah, K. Magnetic gear pole-slip prevention using explicit model predictive control. IEEE/ASME Trans. Mechatronics 2012, 18, 1535–1543. [Google Scholar] [CrossRef]
  16. Han, J. From PID to active disturbance rejection control. IEEE Trans. Ind. Electron. 2009, 56, 900–906. [Google Scholar] [CrossRef]
  17. Ogata, K. Modern Control Engineering; Prentice Hall: Upper Saddle River, NJ, USA, 2010; Volume 5. [Google Scholar]
  18. Tepljakov, A. FOMCON: Fractional-order modeling and control toolbox. In Fractional-Order Modeling and Control of Dynamic Systems; Springer: Berlin/Heidelberg, Germany, 2017; pp. 107–129. [Google Scholar]
  19. Immovilli, F.; Cocconcelli, M.; Bellini, A.; Rubini, R. Detection of generalized-roughness bearing fault by spectral-kurtosis energy of vibration or current signals. IEEE Trans. Ind. Electron. 2009, 56, 4710–4717. [Google Scholar] [CrossRef]
  20. Hori, Y.; Sawada, H.; Chun, Y. Slow resonance ratio control for vibration suppression and disturbance rejection in torsional system. IEEE Trans. Ind. Electron. 1999, 46, 162–168. [Google Scholar] [CrossRef]
Figure 1. Experimental test rig (a) motor- and load-side actuation (b) magnetic coupling.
Figure 1. Experimental test rig (a) motor- and load-side actuation (b) magnetic coupling.
Energies 15 08148 g001
Figure 2. Block diagram of MDTs with acceleration feedback.
Figure 2. Block diagram of MDTs with acceleration feedback.
Energies 15 08148 g002
Figure 3. Block diagram to optimize controller parameters.
Figure 3. Block diagram to optimize controller parameters.
Energies 15 08148 g003
Figure 4. Simulation results of on-load startup: (a) R = 1 ( J M = 0.001 kg·m 2 ); (b) R = 0.5 ( J M = 0.002 kg·m 2 ); (c) R = 2 ( J M = 0.0005 kg·m 2 ).
Figure 4. Simulation results of on-load startup: (a) R = 1 ( J M = 0.001 kg·m 2 ); (b) R = 0.5 ( J M = 0.002 kg·m 2 ); (c) R = 2 ( J M = 0.0005 kg·m 2 ).
Energies 15 08148 g004
Figure 5. Experimental responses using the optimized controller with K a = 0.1 (speed reference ω r e f is imposed at 1 s; T L is set to 0%, 63%, and 0% of the pull-out torque at 1–13 s, 13–28 s, and 28–40 s, respectively): (a) shaft speed of on both sides of the magnetic coupling; (b) T e m and T L .
Figure 5. Experimental responses using the optimized controller with K a = 0.1 (speed reference ω r e f is imposed at 1 s; T L is set to 0%, 63%, and 0% of the pull-out torque at 1–13 s, 13–28 s, and 28–40 s, respectively): (a) shaft speed of on both sides of the magnetic coupling; (b) T e m and T L .
Energies 15 08148 g005
Figure 6. Experimental measurement while starting with 63% of the pull-out torque ( ω r e f and T L are imposed at 1 s, simultaneously): (a) speed of both side shafts; (b) T e m and T L .
Figure 6. Experimental measurement while starting with 63% of the pull-out torque ( ω r e f and T L are imposed at 1 s, simultaneously): (a) speed of both side shafts; (b) T e m and T L .
Energies 15 08148 g006
Figure 7. Pole-slipping detection using kurtosis: (a) the output of T e m with respect to T L ; (b) kurtosis of T e m (pole-slipping occurs at 12.1 s, and kurtosis is calculated individually with 200 samples of T e m output).
Figure 7. Pole-slipping detection using kurtosis: (a) the output of T e m with respect to T L ; (b) kurtosis of T e m (pole-slipping occurs at 12.1 s, and kurtosis is calculated individually with 200 samples of T e m output).
Energies 15 08148 g007
Figure 8. Experimentally measured speed of both coupled shafts components before, during, and after the overload condition.
Figure 8. Experimentally measured speed of both coupled shafts components before, during, and after the overload condition.
Energies 15 08148 g008
Table 1. ω a and ω n of the MDT rig subject to various percentages of pull-out torque conditions.
Table 1. ω a and ω n of the MDT rig subject to various percentages of pull-out torque conditions.
Percentage
(%)
050759099
ω n (rad/s)126.0119.0103.071.148.4
ω a (rad/s)89.483.272.759.133.6
Table 2. Variance of ω M .
Table 2. Variance of ω M .
Time (s)9–1111–1319–2125–2727–29
σ 2.01 × 10 6 0.880.020.71 4.70 × 10 5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liao, X.; Bingham, C.; Smith, T. Speed Control of Magnetic Drive-Trains with Pole-Slipping Amelioration. Energies 2022, 15, 8148. https://doi.org/10.3390/en15218148

AMA Style

Liao X, Bingham C, Smith T. Speed Control of Magnetic Drive-Trains with Pole-Slipping Amelioration. Energies. 2022; 15(21):8148. https://doi.org/10.3390/en15218148

Chicago/Turabian Style

Liao, Xiaowen, Chris Bingham, and Tim Smith. 2022. "Speed Control of Magnetic Drive-Trains with Pole-Slipping Amelioration" Energies 15, no. 21: 8148. https://doi.org/10.3390/en15218148

APA Style

Liao, X., Bingham, C., & Smith, T. (2022). Speed Control of Magnetic Drive-Trains with Pole-Slipping Amelioration. Energies, 15(21), 8148. https://doi.org/10.3390/en15218148

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