Next Article in Journal
Optimal Placement and Sizing of DGs in Distribution Networks Using MLPSO Algorithm
Next Article in Special Issue
Fault Diagnosis and Tolerant Control of Three-Level Neutral-Point Clamped Inverters in Motor Drives
Previous Article in Journal
Performance Review of Three Car Integrated ABS Types: Development of a Tire Independent Wheel Speed Control
Previous Article in Special Issue
Design of High-Speed Permanent Magnet Motor Considering Rotor Radial Force and Motor Losses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Optimal Torque Feedforward and Modal Current Feedback Control for Low Inductance PM Motors

Chair of Mechatronics, Otto von Guericke University, 39106 Magdeburg, Germany
*
Author to whom correspondence should be addressed.
Energies 2020, 13(23), 6184; https://doi.org/10.3390/en13236184
Submission received: 7 October 2020 / Revised: 10 November 2020 / Accepted: 22 November 2020 / Published: 25 November 2020
(This article belongs to the Special Issue Design and Control of Electrical Motor Drives)

Abstract

:
Small sized electric motors providing high specific torque and power are required for many mobile applications. Air gap windings technology allows to create innovative lightweight and high-power electric motors that show low phase inductances. Low inductance leads to a small motor time constant, which enables fast current and torque control, but requires a high switching frequency and short sampling time to keep current ripples and losses in an acceptable range. This paper proposes an optimal torque feedforward control method, minimizing either torque ripples or motor losses, combined with a very robust and computation-efficient modal current feedback control. Compared to well-known control methods based on the Clarke-Park Transformations, the proposed strategy reduces torque ripples and motor losses significantly and offers a very fast implementation on standard microcontrollers with high robustness, e.g., against measurement errors of rotor angle. To verify the accuracy of the proposed control method, an experimental setup was used including a wheel hub motor built with a slotless air gap winding of low inductance, a standard microcontroller and GaN (Gallium Nitride) Power Devices allowing for high PWM switching frequencies. The proposed control method was validated first by correlation of simulation and experimental results and second by comparison to conventional field-oriented control.

1. Introduction

Manufacturers of electric vehicles or electric planes require electric motors with small size, lightweight construction and maximum efficiency, providing at the same time high torque and power. For low-power applications (<200 W), electrical motor weight can be reduced by air gap winding designs offered, e.g., by Faulhaber and Maxxon [1]. A new slotless air gap winding design for a 15-inch wheel-hub motor providing high power of 40 kW at a speed of 1500 rpm and a very low weight of 20 kg is presented in References [2,3,4]. Wheel hub motors take a new step in technology of electric vehicles as they allow to create completely new designs of mobile applications without a lot of the mechanics used today and a lot more room for passengers and baggage. High specific power and torque densities are the main indicators of the quality of wheel hub motors. This way, the very lightweight wheel hub motors based on air gap windings [2,3,4] are most promising.
Usually, motors are based on slotted windings, which significantly increase the motor weight and reduce specific torque and power density. Slotless air gap winding design reduces the weight of electric motors, as it avoids a big part of stator back-iron needed for conventional slotted machines. Additionally, a very high torque can be generated without any torque ripple induced by slotted geometry. However, low number of turns and low iron volume in air gap winding designs lead to very low phase inductances of only several µH, which makes control a challenge [4,5]. For this reason, a high pulse width modulation (PWM) switching frequency of about 100 kHz is needed to keep current ripples, and accordingly torque ripples and associated losses, in a demanded range. In addition, the motor control must be adapted to low phase inductance and high PWM frequency for microcontrollers with a very small sampling step size.
Depending on the phase current waveform, conventional control methods for permanent magnet (PM) motors can be divided into two groups: based on the block commutation [5,6,7,8,9] and based on the Clarke-Park Transformations [10,11,12,13,14,15,16]. The first one is usually used for low-cost and low-power applications because of low manufacturing costs and simple design. On the other hand, this leads to a high torque ripple waveform caused by non-ideal commutation currents [5,6,7,8,9]. This disadvantage is even more significant for PM motors with low inductance due to the inverse relationship between phase inductances and current ripples. Trigonometric Clarke-Park Transformations transfer phase currents into a torque generating Q-part and a field generating D-part (DQ control), dependent on rotor position [10,11,12,13,14,15,16]. This gives the opportunity for smoother phase commutations, which decreases current and torque ripple but assumes pure sinusoidal waveforms for B-field (magnetic flux density) and phase currents. For high sampling frequencies, calculation of transformations is a high burden for the microcontroller. If the B-field consists of significantly higher harmonics, they should be compensated by adding additionally higher harmonic parts to the phase currents. This method may be unacceptable for high-frequency control, and only iterative approaches are known [10,11,12,13,14].
This article proposes a combined optimal torque feedforward and modal current feedback control method for PM motors with very low phase inductance. Previously, in References [17,18], a high-frequency Combined Optimal Torque and Modal Current Control (OTMIC) was proposed for low inductance PM motors focused on the minimization of motor losses. In this article, this method for loss minimization is described as mathematically complete. Further, it is extended to torque ripple minimization to give a powerful optimization method flexibly adaptable to different goals. Limitations of References [17,18] given by growing size of symbolic expressions will be avoided by a new matrix formulation of the algorithm that allows for adaptive real-time applications. It is possible to switch online between both optimization strategies to fulfill changing demands depending on certain operation points (e.g., low speed → minimal torque ripple, high speed → minimal loss). In case of motors based on air gap winding designs, motor B-field is nearly independent of speed, torque and phase currents. This conclusion was used to explicitly take into consideration all harmonics of B-field for torque control. In an off-line optimization step, optimal current values for a set of electric angles are stored in a table for on-line usage under hard real-time conditions. By means of a very robust modal current control approach, which works independently on rotor angle, the 3-phase current control problem can be reduced to a 2-phase modal one, based on a linear model of the PM machine with constant coefficients. It becomes evident that the non-controllable state (one of the 3-phase currents) can be eliminated without any pre-conditions, e.g., on B-field waveform, and all statements on stability are global using this approach. Additionally, modal transformation matrices are constant, especially independent of electrical angle and any physical parameter. This reduces computation effort significantly compared, e.g., to Park-Clarke Transformation, where transformation matrices change with the angle. Errors of angle measurement or estimation have only a very limited effect on current control. Thus, independence of motor parameters increases the robustness of the transformation due to any physical changes, further allowing for a very fast software implementation on a low-cost microcontroller with an acceptable sampling step size. Optimization as well as control problem can be solved symbolically, allowing for very fast and adaptive implementations.
In the second section of this article, a new approach of PM motor control based on the mathematical model of slotless air gap winding PM motor with very low phase inductance is presented. The proposed method combines modal feedback control to guarantee the desired stability, dynamics and robustness margins with optimal feedforward control to minimize important motor features like ohmic losses or torque ripple. The last section shows the experimental and simulation results of the proposed method and comparison of the proposed control with conventional Field-Oriented Control (FOC) based on DQ Clarke-Park Transformations (DQ control). A PM wheel hub motor with slotless air gap windings and very small inductances (near 1.5 μH per phase) was used as a test object.

2. Design of the Proposed Combined Optimal Feedforward Torque and Modal Current Feedback Control for Low Inductance PM Motors

2.1. Operation Principle of PM Motor with Very Low Phase Inductance

The dynamic model of a PM motor with slotless air gap winding is obtained considering the following assumptions [3]:
1.
Rotor consists of surface-mounted permanent magnets.
2.
No cogging torque due to slotless design.
3.
Harmonic magnetic field density generated by permanent magnets.
4.
B-field harmonics and inductance are constant, e.g., independent from speed and phase currents due to very small current-induced magnetic fields of the windings.
5.
Stator has symmetric 3-phase star-connected windings commutated by a 3-phase inverter (B6-bridge).
A model is shown in Figure 1, where R —phase resistance, L ,   M —phase self and mutual inductance, e _ = [ e a e b e c ] T —phase back-EMFs (back electromotive force), u _ = [ u a u b u c ] T , i _ = [ i a i b i c ] T —phase voltages and currents, u S —voltage at star point and u D C ,   i D C —input direct voltage, current to PM motor inverter [2,3,4].
Considering the Y-connection, Kirchhoff’s first Law
i a + i b + i c = 0
and Kirchhoff’s second Law
[ u a u b u c ] = [ R 0 0 0 R 0 0 0 R ] [ i a i b i c ] + [ L M M M L M M M L ] [ d i a / d t d i b / d t d i c / d t ] + [ e a e b e c ] + [ 1 1 1 ] u S
or more compact
u _ = R _   i _ + L _   d i _ / d t + e _ + u _ S
state the mathematical model of a 3-phase PM motor. The 3-phase back-EMF values
e _ = ω   k M   B _ ( φ ) ,
depend on ω —the mechanical angular velocity of the rotor, k M —a motor geometric constant and B _ = [ B a ( φ ) B b ( φ ) B c ( φ ) ] T   —the vector of averaged B-fields acting on air gap winding phases [3], which itself depends on φ —the electrical angle. Obviously, B a ( φ ) = B ( φ ) , B b ( φ ) = B ( φ 2 π 3 ) and B c ( φ ) = B ( φ 4 π 3 ) are given by the harmonic magnetic field density generated by permanent magnets, calculated by:
B ( φ ) = k b k sin ( k   φ ) .
as a function of electric angle and odd harmonic coefficients, b k . Either FE analysis of the underlying magnetic circuit or line-fitting of measured back-EMF values give the B-field harmonic coefficients with very good accuracy [2,3,4].
Electrical torque is generated by Lorentz force [2,3,4]:
T E = k M [ B a B b B c ] [ i a i b i c ] = k M B _ T i _
acting on the airgap winding. Loss torque:
T L = d   ω + c
collects linear mechanical friction together with back-iron and airgap winding eddy losses into the linear friction coefficient d and mechanical coulomb friction together with back-iron hysteresis losses into the constant friction coefficient c . It will be noted that loss torque T L and resulting loss power depends only on speed and is independent of motor current. Electrical torque T E reduced by loss torque T L gives motor output torque:
T M = T E T L .
Equations (1)–(8) describe the motor model in terms of back-EMF, output torque and phase currents based on information of angular position φ , angular velocity ω , B-field harmonics and motor parameters k M , d and c . Torque Equation (6) shows that the adaptation of phase current waveforms taking into account the B-field harmonics is the key to motor torque control [2,3,4].
Phase currents i ( φ ) , like B-field, can also be defined as harmonic functions of electrical angle φ :
i ( φ ) = k a k sin ( k   φ ) .
Cosine terms are omitted. If included, they would vanish during optimization. From (1), we get additional constraints:
a 3 k = 0
for current coefficients. Average electrical torque:
T A = 1 2 π 0 2 π T E   d φ = 3 2   k M k a k b k
can then be expressed by Fourier coefficients of phase current and B-Field. Using Equations (5) and (8)–(11), in this paper will be proposed a new method to generate optimal phase currents to generate a given torque T R e f for an arbitrary harmonic B-field, while minimizing either the motor losses or alternatively while minimizing the output torque ripples.

2.2. Optimal Control for Minimal Motor Losses

Generation of a given torque T R e f with maximum efficiency or, equivalently expressed, with minimum losses P L is a typical design requirement of an electrical motor. Averaged total motor losses consist of ohmic and non-ohmic losses [3]:
P L = 1 2 π 0 2 π ( R ( i a 2 + i b 2 + i c 2 ) + d   ω 2 + c   ω ) d φ   P L = 3 R 2 k a k 2 + d · ω ¯   2 + c · ω ¯ ,
where non-ohmic losses depend only on the average angular speed in a pole pair segment:
ω ¯ = 0 2 π ω d φ and ω ¯   2 = 0 2 π ω 2 d φ
and ohmic losses can be expressed by Fourier coefficients of phase currents according to (9). Given demands lead to the mathematical optimization problem for minimization of motor losses P L while generating a demanded torque T R e f :
M i n = P L with   constraint T A = T R e f
for unknown parameters a k of the phase current harmonics based on the given harmonics b k of the B-field. The solution of this quadratic optimization problem can be determined analytically, e.g., by combining cost function and constraints in a Lagrange function, calculating derivatives and solving first-order conditions for unknown coefficients:
a m = 2 3 T R e f k M b m b ¯ 2
with the sum of squares of selected harmonics of B-field:
b ¯ 2 = m b m 2
and the selection indices:
m = 1 , 5 , 7 , 11 , 13 , 17 , 19 , 23
Remaining parameters a k are zero. Selection indices show that some odd harmonics of phase current {3, 9, 15, …} will not contribute to torque. The optimal current waveforms can be calculated from (15) for each electrical angle taking into account the harmonics of B-field, motor constant k M and given reference torque T R e f . This approach achieves that each demanded torque can be generated with minimum motor losses for arbitrary B-field and regardless of motor speed. Air gap winding PM motor design allows to have nearly ideally constant B-field harmonics with minor sensitivity to temperature changes. If necessary, e.g., for application to other motor types, the effect of temperature on B-field and phase resistances can be easily incorporated into the model and adapted online.

2.3. Optimal Control for Minimal Torque Ripples

Generation of a given torque T R e f with minimum torque ripples is another typical design requirement of an electrical motor, especially if vibrations and noise are critical. Torque ripple is defined as mean square deviation between electrical torque T E and average torque T A :
T R 2 = 1 2 π 0 2 π ( T E T A ) 2   d φ
In this case, optimal current can be determined solving the optimization problem:
M i n = T R 2 with   constraint   T A = T R e f
Inserting current definition (9) into the integral (18), the optimization problem stated in (19) again can be solved analytically by formulating Lagrange function, calculating derivatives and solving the first-order conditions. Besides this straightforward path, there exists another approach, which leads to more simple expressions and allows further insight into the solution. Starting from (5), (6) harmonics
t k = 1 π 0 2 π T E cos ( 6   k   φ ) d φ
of electrical torque T E can be determined and collected in vector
t _ = A _   a _
with matrix
A _ = 2 3 k M [ 2 b 1 2 b 5 2 b 7 2 b 11 2 b 13 2 b 17 2 b 19 2 b 23 b 7 b 5 b 11 b 1 b 13 + b 1 b 17 + b 5 b 19 + b 7 b 23 + b 11 b 13 b 17 b 13 b 11 b 17 b 7 b 19 b 5 b 23 b 1 b 1 b 5 b 7 b 11 b 19 b 17 b 23 b 13 b 11 b 7 b 5 b 1 b 1 b 5 b 23 b 19 b 17 b 13 b 11 b 7 b 5 b 1 0 0 b 23 b 19 b 17 b 13 b 11 b 7 0 0 0 0 b 23 b 19 b 17 b 13 0 0 0 0 0 0 b 23 b 19 ]
and vector of current harmonics
a _ = [ a 1 , a 5 , a 7 , a 11 , a 13 , a 17 , a 19 , a 23 ]
for an example with 23 harmonics. Only average value at zero and multiples of the sixth harmonic appear. All other coefficients are zero. Matrix A _ is square and has full rank. Defining the torque harmonics
t _ = [ 2   T R e f , 0 , 0 , 0 , 0 , 0 , 0 , 0 ]
in a way that the zero coefficient is given by the demanded torque T R e f and all other torque harmonics vanish, delivers the solution for the above-mentioned optimization problem by solving the linear system of equations defined by (21). For practical applications, a much smaller number of harmonics than given in (23) will be sufficient. The larger number presented here helps to demonstrate the structure of the matrix. It is important to state that also the torque ripples minimization leads to an analytical solution for optimal currents that is well suited for online and adaptive applications.

2.4. Modal Current Control

So far, we have two feedforward strategies to implement either loss minimal or torque ripples minimal operation of the motor. Both strategies rely on current control, which is able to deal with higher current harmonics, which means also with higher frequencies. A modal current control method, exploiting all structural benefits of given control problem and allowing for very high-speed control, will be presented in more detail in this section. The complete model of the PM motor with air gap winding is defined by (1)–(3). Node Equation (1) shows that only 2 of the 3 phase currents are independent. Taking this dependency into account, (3) can be written as
G _ ( u _ e _ ) = R _   i _ + L _ d i _ / d t
with
G _ = 1 3 [ 2 1 1 1 2 1 1 1 2 ] .
For simplicity, back-EMF e _ can be neglected and later added to input voltage u _ as a feedforward compensation or simply treated as a disturbance for the controller. This coupled system for the 3 phase currents i _ can be completely decoupled with a modal transformation
[ J Σ J 1 J 2 ] = 1 3 [ 1 1 1 1 1 2 1 2 1 ] W _ [ i a i b i c ]
applied to system (25) in a way that modal inductance matrix
Λ _ =   W _   L _   W _ 1 = [ L 2 M 0 0 0 L + M 0 0 0 L + M ]
gets diagonal. Phase currents i _ are transformed to modal currents J _ =   [ J Σ J 1 J 2 ] T . Considering (1) and the first row of the transformation matrix defined in (27), modal current J Σ will be zero. Thus, only 2 out of 3 modal currents associated with the modal inductances L + M are relevant. The decoupled modal system
1 3 [ 0 0 0 1 1 2 1 2 1 ] [ u a u b u c ] = [ R 0 0 0 R 0 0 0 R ] [ J Σ J 1 J 2 ] + [ L 2 M 0 0 0 L + M 0 0 0 L + M ] [ d J Σ / d t d J 1 / d t d J 2 / d t ]
shows that modal current J Σ is not controllable as the associated row of the modal input matrix is completely zero. The remaining 2 modal currents, J 1 and J 2 , come from 2 identical first-order systems with time constant L + M R . Eliminating non-controllable current J Σ gives the reduced controllable system
1 3 [ 1 1 2 1 2 1 ] Γ _ [ u a u b u c ] = [ R 0 0 R ] [ J 1 J 2 ] + [ L + M 0 0 L + M ] [ d J 1 / d t d J 2 / d t ]
Thus, the 3-phase current control problem is reduced to the control of two completely identical but decoupled first-order systems, associated with modal currents J 1 and J 2 . The redundant total current J Σ displays any errors in current control loop and can be used as an indicator of the quality of the system online. Modal currents J 1 and J 2 are not connected and can be controlled completely independent by modal input voltages
V _ = [ V 1 V 2 ] =   Γ _   ·   u _ ,
as shown in Figure 2. This simple separation based on the modal transformation given in (27) does not use any trigonometric calculations based on the rotor angle as in conventional Clarke-Park Transformations, is independent of motor parameters and is well defined for any shape of B-field.
According to (30). motor transfer function is given by
G M ( s ) = 1 / R 1 + L + M R s .
For high-speed control of low inductance PM motors with high number of poles, current sensor dynamics must be included in the control loop with a first-order sensor model
G S ( s ) = 1 1 + T S s
and sensor time constant T S . As a microcontroller-driven B6-bridge generates PWM pulses at frequency f P W M = 1 / Δ t to control motor input voltage and ADC (Analog-to-Digital Converter) samples current sensor signals with the same rate, time continuous transfer function
G M S ( s ) = G M ( s ) G S ( s )
has to be discretized with that frequency to get the discrete transfer function
D M S ( z ) = ( β 1 + ( 1 α ) δ ) z + ( ( α 1 ) δ α ) β + α R ( δ 1 ) ( z α ) ( z β )
with sampling time Δ t , discrete eigenvalues of motor α = e R Δ t L + M , sensor β = e Δ t T S and motor/sensor speed ratio δ = L + M R   T S . Very small motor time constant due to low phase inductance allows to achieve very fast current and electrical torque response with requested time constant T R e q , which expands the possibilities of design for more complex control systems like dynamic stability control or autonomous driving. To get a controlled system of first order
D ( z ) = 1 z R z z R
with requested eigenvalues z R = e Δ t T R e q , it is necessary to use a discrete PID (Proportional-Integral-Derivative) controller with filtered derivative action:
D C ( z ) = K P + K I z 1 + K D N D + 1 z 1 .
Control parameters can be determined by comparison of coefficients of discrete closed loop
D ( z ) = D C ( z ) D M S ( z ) 1 + D C ( z ) D M S ( z )
and reference transfer function given in (36), as
K P = ( 1 α + δ ( β 1 ) ) R ( z R 1 ) ( δ 1 ) ( β 1 ) ( α 1 ) ;   K I = R ( z R 1 ) ;                                                                                 K D = R δ ( z R 1 ) ( β α ) 2 ( δ 1 ) 2 ( β 1 ) 2 ( α 1 ) 2 ; N D = β 1 + δ ( 1 α ) ( δ 1 ) ( β 1 ) ( α 1 ) .              
A saturation of output and I-part of the PID controller guarantees safe operation in the case of limited supply voltages. Outputs from PID controllers D C ( z ) are two modal control voltages, V 1 and V 2 . According to (31), V 1 and V 2 must be transformed back to three control voltages u _ by means of a Moore Penrose pseudo-inverse transformation matrix Γ _ +
[ u a u b u c ] = [ 1 1 0 1 1 0 ] Γ _ + [ V 1 V 2 ] .
Modal current control ensures system stability regardless of fluctuations in system parameters caused by changes of motor temperature and other factors. Moreover, it can be implemented for PM synchronous motors with any waveform of B-field. However, it should be noted, that the motor model should satisfy the mathematical model described with Equations (1)–(3), which can be very good when performed for PM motors with air gap windings and also other PM motors.

2.5. Control Diagram of OTMIC

The proposed control diagram of Combined Optimal Torque and Modal Current Control (OTMIC) is shown in Figure 3. The control algorithm uses a torque reference T R e f and measured electrical angle φ to estimate optimal phase current waveforms for motor loss or torque ripples minimization, i _ R e f = [ i a R e f ( φ ) i b R e f ( φ ) i c R e f ( φ ) ] T , based on one of the optimization methods described in Section 3.1 and Section 3.2. The optimal currents i _ R e f are used as a reference for measured phase currents i _ . Current errors Δ i _ as differences between reference and measured currents are fed to the Modal Current Control described in Section 2.4 and used to generate the control voltages u _ , which set PM motor in motion with the optimal torque waveform via a SVPWM (Space Vector Pulse Width Modulation) and a B6-bridge. Architecture of OTMIC control allows to simplify the microcontroller torque/current control loop by using constant matrices with simple coefficients for modal transformation.
In addition, optimal reference currents according to (9) and (15)–(17) for minimal Loss Control or from (9) and (21)–(24) for minimal Torque Ripple Control can be calculated off-line for a set of angles and stored in a table. These simplifications allow to implement the control loop with very small step size of Δ t < 10   μ s on a common low-cost microcontroller.

3. Units Simulation and Experimental Results

3.1. Motivation and Experimental Setup

A wheel hub motor for a hub-less scooter, shown in Figure 4 on the experimental setup (left above) and in the dissembled state (right) for better demonstration of the air gab winding, will be used for verification of OTMIC control method. The motor has a very low total weight of 2.7 kg, a high nominal power of 4.8 kW and a nominal torque of 80 Nm. More motor data is given in Table 1, where B odd Harmonics are 1st, 3rd, 5th and 7th harmonic numbers. The experimental platform shown in Figure 4 was used to verify the proposed OTMIC method with either minimized motor losses or minimized torque ripples for low inductance PM motors.
Due to elimination of stator iron, motor inductance is very low (near 1.5 μH per phase). Thus, motor needs an adapted control with high switching frequency of about 100 kHz to keep current ripples and total losses in an acceptable range. Control parameters calculated according to (39) are also given in Table 1.
To drive the motor, six GS61008P MOSFETs (metal–oxide–silicon transistors) from company GaN Systems, able to switch up to 90 A at a rate up to 300 kHz, were connected as a B6-bridge. Due to limited MOSFETs cooling of experimental boards, phase current was restricted to 20 A for the experiments. Hall-effect-based sensors from Allegro ACS730 were used for phase current measurement with a very small sensor time constant of 1 μs. Rotor angle is delivered by an encoder using 5 k increments. 100 kHz PWM switching frequency to power the B6-bridge requires a very small step size in the microcontroller (lower than 10 µs). To achieve this requirement, an ARM Cortex-M4F has been used running at 200 MHz. It allows to implement the control loop shown in Figure 3 in 8.3 µs. In comparison to the proposed method, a conventional DQ control needs a minimum of 9.4 µs, which underlines the effectiveness of the proposed method. This comparison can be seen in Figure 5.
This means that OTMIC control can be calculated faster and in addition, gives the optimal torque in contrast to FOC (DQ) based on the Clarke-Parke Transformations.

3.2. Verification of the Proposed Method

The simulation results presented in this study have been performed through Matlab/Simulink environment based on the mathematical model of the PM motor described in Section 2 and the control model presented in Section 3. Conventional DQ field-oriented vector control used as a reference to the proposed OTMIC was also implemented in the experiment to control the wheel hub motor, with the parameters in Table 1.
The motor dynamic behavior can be verified by means of step response of electrical torque T E . As shown in Figure 6, measured and simulated torque of proposed OTMIC and conventional DQ methods were carried out according to zero velocity, mechanically fixed rotor and 5 Nm reference torque T R e f (green). Because torque sensor bandwidth was too low for this experiment, electrical torque generated in air gap winding (grey) was calculated using (6) based on the values of measured currents and the electric angle with a given motor geometric constant k M , and compared with the simulation for OTMIC (blue) and DQ (red) control methods. As can be seen form Figure 6, the difference between both methods is minimal due to the very low motor time constant. The exciting fast torque step response within 35 µs underlines the very high dynamics of this type of motor. The slight difference between simulation and measurement is due to parasitic inductances and capacitances in power lines and supply. This can be easily considered in the motor model but was not done in this experiment to show the robustness of the control.
For a more detailed demonstration of the basic principles of the proposed control algorithm, Figure 7 shows the results of modeling the ideal optimal current waveforms to minimize the torque ripples (OTMIC TRO) and losses (OTMIC LO) for two waveforms of B-Field (Case A with odd harmonic coefficients b k from Table 1 and Case B with b k from Reference [19], which result from a change in the air gap between stator and rotor). For a more detailed analysis, the torque RMS (root mean square) and efficiency value are also presented for a sinusoidal current of conventional DQ control based on the Clarke-Park Transformations. Figure 7 shows the simulation results of an ideal system without any errors in the phase current or rotor position measurements. As can be seen, the optimal phase current waveforms are able to the minimize torque ripples or maximize motor efficiency based on the Equations (6) and (12). The advantages of the OTMIC TRO are seen very well in Figure 7. In an ideal system, OTMIC TRO eliminants any torque fluctuations for any waveform of the B-field. OTMIC with loss optimization (LO) changes the phase current waveform and improves loss rates and efficiency of the system with increased torque ripples. The advantages of the OTMIC LO in comparison with the standard DQ method are visible only in the presence of significant parasitic harmonics in the form of a B-field. Simply said, the ideal sinusoidal B-field leads to the same torque ripples and motor loss for every method.
Simulated and measured phase currents and torque results during normal operation of proposed OTMIC control with minimization of torque ripple (TRO Figure 8a—simulation, Figure 9a—measurement) and losses (LO Figure 8b and Figure 9b) as well as a conventional vector DQ control (DQ Figure 8c and Figure 9c) correspond to the references: ω = 8   rad / s and T Ref = 10   Nm . The behavior of phase currents and torque for both OTMIC optimizations given in Figure 9a,b has a slight asymmetry caused by speed fluctuations of load machine and errors in angle measurements, which could not be completely avoided in experimental setup. The current and torque harmonic distribution for the aforementioned control approaches are shown in Figure 8d for simulation and in Figure 9d for measured results. The amplitude of each harmonic component is normalized relative to that of a fundamental one. The FFT (Fast Fourier transform) analyses of simulation currents and torques underline the main idea of the proposed method. As shown in Figure 8d (left), the phase current waveforms of OTMIC have 5th and 7th harmonics expected to conventional control. On the other hand, the torque waveform of OTMIC with torque ripple optimization TRO (blue) is deprived of such harmonics, which leads to the very small torque fluctuations as opposed to DQ control, as shown in Figure 8d (right).
The harmonic spectrum of measured phase current with OTMIC control shown in Figure 9d (left) has parasitical 2nd, 4th and 6th harmonics due to inaccuracies in experimental setup, mainly caused by interactions of control system and load machine, but they are less than 1% of that of the fundamental one. The necessary 5th and 7th harmonics have a slight deviation of about 0.5% of the fundamental one compared to simulation results. The current waveform of conventional DQ control should be quite similar to a sinusoidal, as shown in the simulation results in Figure 8c and harmonic analyses in Figure 8d (left). However, the harmonic spectrum of measured phase currents by the conventional method consist of high 3rd and 4th harmonics due to the high sensitivity of the standard method based on the trigonometric transformations to errors in the control system, especially errors in angle measurement.
The measured torque harmonic distribution shown in Figure 9d (right) consists of a large number of harmonics due to speed errors in rotating shaft and additional torque produced by load machine. However, the 6th and 12th harmonics retain their tendency and can be compared with the simulation. In addition, the FFT analyses of measured torques in Figure 9d (right) shows that the proposed OTMIC control with torque ripple minimization is able to compensate additional harmonics in B-field. Summarized, the measured results show very good correspondence to theoretical calculations. As opposed to OTMIC, the Clarke-Park Transformations need a very accurate rotor angle measurement, which can be seen in the larger asymmetries given in Figure 9c. Thus, the proposed OTMIC method is less sensitive to speed fluctuations and angle errors due to the absence of trigonometric transformations of phase currents based on the electric angle.
Figure 10 compares the measured values of motor efficiency for OTMIC control with LO and TRO and conventional DQ control. Motor efficiency was measured with given constant velocity of 30 rad/s, variable torque from 1 to 10 Nm and output power from 30 to 300 W. Results show that OTMIC control with loss minimization increases the system efficiency by 0.7% in compassion to the conventional method, while delivering the same torque. However, for a correct analysis of the proposed control in case of a motor loss and efficiency, any errors in the measurement system should be avoided and motor torque and current should be increased to the rated values, which could not be fully realized in our experimental setup. More detailed analyses of the system efficiency should be investigated in the future work.
Simulation and experimental results show that the OTMIC control system is less sensitive to errors in the electrical angle and phase currents measurements compared to FOC. In addition, the proposed method worked very robustly independent of deviations of motor parameters like phase resistance due to the temperature and the speed of fluctuations due to the load motor.

4. Conclusions

This paper proposed a Combined Optimal Torque feedforward and Modal Current feedback Control architecture (OTMIC) for PM motors with very low inductance and higher harmonic B-field, which minimizes motor losses or torque ripples/motor noise by means of offline calculation of optimal current waveform. Optimal three-phase currents are used as reference values to the measured currents. The optimization algorithm considers all relevant coefficients of B-field harmonic distribution to calculate the appropriate current waveform off-line and employs an extremely fast table-lookup for real-time implementation in the microcontroller. Modal current feedback control provides stability, high system dynamics without deviation of a given reference and high robustness against variation of motor parameters or errors in rotor angle measurement.
Compared to classical field-oriented control methods based on the cumbersome and complex Clarke-Park Transformations, the proposed control is simpler for implementation on the microcontroller and allows to minimize either torque ripples or motor losses. As a result, the proposed OTMIC control can reach its full potential in high-frequency systems, where simplicity of the control algorithm and microcontroller calculation time matters. Due to the transformations used in the conventional DQ method, the phase currents’ waveforms have to be nearly ideally sinusoidal. This means that conventional methods to achieve the specified requirements of motor torque or total losses demand the special concepts and additional real-time calculations to adapt the phase current waveform to the non-sinusoidal waveform of the magnetic field. It complicates the control system and increases the calculation time in the microcontroller. Moreover, in real high-frequency control systems with errors in current and angle measurements, this problem can be more significant. Experimental results show that the Clarke-Park Transformations for electric drives with a very small motor time constant and high number of poles introduces significant current deviations from the values given. In this case, the proposed OTMIC control method used for PM synchronous motors with a very low phase inductance and tested experimentally shows a very good match between simulated and measured waveforms of phase current and torque. The proposed control can also be applied to other electric motors such as brushless DC (BLDC), asynchronous or PM motors, with any waveform of B-field waveform.

Author Contributions

Conceptualization, R.K. and D.G.; methodology, R.K.; software, R.K. and D.G.; validation, D.G.; formal analysis, R.K. and D.G.; writing—original draft preparation, R.K. and D.G.; writing—review and editing, R.K. and D.G.; funding acquisition, R.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the European Commission under Grant ERDF (European Regional Development Fund) and the County of Saxony-Anhalt.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ragot, P.; Markovic, M.; Perriard, Y. Analytical Determination of the Phase Inductances of a Brushless DC Motor with Faulhaber Winding. IEEE Trans. Ind. Appl. 2010, 46, 1360–1366. [Google Scholar] [CrossRef]
  2. Borchardt, N.; Heinemann, W.; Kasper, R. Design of a wheel hub motor with air gap winding and simultaneous utilization of all magnetic poles. In Proceedings of the 2012 IEEE International Electric Vehicle Conference, Greenville, SC, USA, 4–8 March 2012; pp. 1–7. [Google Scholar]
  3. Borchardt, N.; Kasper, R. Nonlinear design optimization of electric machines by using parametric Fourier coefficients of air gap flux density. In Proceedings of the IEEE International Conference on Advanced Intelligent Mechatronics (AIM), Banff, AB, Canada, 12–15 July 2016; pp. 645–650. [Google Scholar]
  4. Kasper, R.; Borchardt, N. Boosting Power Density of Electric Machines by Combining Two Different Winding Types. In Proceedings of the 7th IFAC Symposium on Mechatronic Systems, Leicestershire, UK, 5–8 September 2016; pp. 322–329. [Google Scholar]
  5. Tan, B.; Hua, Z.; Zhang, L.; Fang, C. A New Approach of Minimizing Commutation Torque Ripple for BLDCM. Energies 2017, 10, 173. [Google Scholar] [CrossRef] [Green Version]
  6. Li, X.; Jiang, G.; Chen, W.; Shi, T.; Zhang, G.; Geng, Q. Commutation Torque Ripple Suppression Strategy of Brushless DC Motor Considering Back Electromotive Force Variation. Energies 2019, 12, 1932. [Google Scholar] [CrossRef] [Green Version]
  7. Li, H.; Zheng, S.; Ren, H. Self-Correction of Commutation Point for High-Speed Sensorless BLDC Motor with Low Inductance and Nonideal Back EMF. IEEE Trans. Power Electron. 2017, 32, 642–651. [Google Scholar] [CrossRef]
  8. Li, X.; Xia, C.; Cao, Y.; Chen, W.; Shi, T. Commutation Torque Ripple Reduction Strategy of Z-Source Inverter Fed Brushless DC Motor. IEEE Trans. Power Electron. 2016, 31, 7677–7690. [Google Scholar] [CrossRef]
  9. Yoon, K.-Y.; Baek, S.-W. Robust Design Optimization with Penalty Function for Electric Oil Pumps with BLDC Motors. Energies 2019, 12, 153. [Google Scholar] [CrossRef] [Green Version]
  10. Song, Q.; Li, Y.; Jia, C. A Novel Direct Torque Control Method Based on Asymmetric Boundary Layer Sliding Mode Control for PMSM. Energies 2018, 11, 657. [Google Scholar] [CrossRef] [Green Version]
  11. Qian, J.; Ji, C.; Pan, N.; Wu, J. Improved Sliding Mode Control for Permanent Magnet Synchronous Motor Speed Regulation System. Appl. Sci. 2018, 8, 2491. [Google Scholar] [CrossRef] [Green Version]
  12. Mohd Zaihidee, F.; Mekhilef, S.; Mubin, M. Robust Speed Control of PMSM Using Sliding Mode Control (SMC)—A Review. Energies 2019, 12, 1669. [Google Scholar] [CrossRef] [Green Version]
  13. Wang, Q.; Yu, H.; Wang, M.; Qi, X. A Novel Adaptive Neuro-Control Approach for Permanent Magnet Synchronous Motor Speed Control. Energies 2018, 11, 2355. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, X.; Yu, H.; Yu, J.; Zhao, L. Combined Speed and Current Terminal Sliding Mode Control With Nonlinear Disturbance Observer for PMSM Drive. IEEE Access 2018, 6, 29594–29601. [Google Scholar] [CrossRef]
  15. Zeb, K.; Din, W.U.; Khan, M.A.; Khan, A.; Younas, U.; Busarello, T.D.C.; Kim, H.J. Dynamic Simulations of Adaptive Design Approaches to Control the Speed of an Induction Machine Considering Parameter Uncertainties and External Perturbations. Energies 2018, 11, 2339. [Google Scholar] [CrossRef] [Green Version]
  16. Kang, L.; Cheng, J.; Hu, B.; Luo, X.; Zhang, J. A Simplified Optimal-Switching-Sequence MPC with Finite-Control-Set Moving Horizon Optimization for Grid-Connected Inverter. Electronics 2019, 8, 457. [Google Scholar] [CrossRef] [Green Version]
  17. Kasper, R.; Golovakha, D.; Süberkrüb, F. Combined Optimal Torque and Modal Current Control for Low Inductance PM Motor. In Proceedings of the IEEE International Conference on Mechatronics (ICM), Ilmenau, Germany, 18–20 March 2019; pp. 491–497. [Google Scholar]
  18. Golovakha, D. Combined Optimal Torque and Modal Current Control for Low Inductance PM Motor. Ph.D. Thesis, Otto von Guericke University, Magdeburg, Germany, 2020. [Google Scholar]
  19. Borchardt, N.; Kasper, R. Analytical magnetic circuit design optimization of electrical machines with air gap winding using a Halbach array. In Proceedings of the International Electric Machines and Drives Conference (IEMDC), Miami, FL, USA, 21–24 May 2017; pp. 1–7. [Google Scholar]
Figure 1. B6-bridge-fed 3-phase permanent magnet (PM) drive connections.
Figure 1. B6-bridge-fed 3-phase permanent magnet (PM) drive connections.
Energies 13 06184 g001
Figure 2. Modal control loop for modal currents J 1 , 2 .
Figure 2. Modal control loop for modal currents J 1 , 2 .
Energies 13 06184 g002
Figure 3. Proposed Combined Optimal Torque and Modal Current Control (OTMIC) diagram of PM motor drive system.
Figure 3. Proposed Combined Optimal Torque and Modal Current Control (OTMIC) diagram of PM motor drive system.
Energies 13 06184 g003
Figure 4. Experimental setup for OTMIC verification [18].
Figure 4. Experimental setup for OTMIC verification [18].
Energies 13 06184 g004
Figure 5. Step response of simulated and measured torque for OTMIC and FOC (DQ).
Figure 5. Step response of simulated and measured torque for OTMIC and FOC (DQ).
Energies 13 06184 g005
Figure 6. Step response of simulated and measured torque for OTMIC and DQ.
Figure 6. Step response of simulated and measured torque for OTMIC and DQ.
Energies 13 06184 g006
Figure 7. Working principles of the proposed OTMIC control algorithm for two B-Field waveforms.
Figure 7. Working principles of the proposed OTMIC control algorithm for two B-Field waveforms.
Energies 13 06184 g007
Figure 8. Simulation results of stator current and torque of: (a) OTMIC with torque ripples optimization, (b) OTMIC with losses optimization, (c) conventional DQ control and (d) harmonic spectrum of phase current (left) and torque (right) for every method [18].
Figure 8. Simulation results of stator current and torque of: (a) OTMIC with torque ripples optimization, (b) OTMIC with losses optimization, (c) conventional DQ control and (d) harmonic spectrum of phase current (left) and torque (right) for every method [18].
Energies 13 06184 g008
Figure 9. Experimental results of stator current and torque of: (a) OTMIC with torque ripples optimization, (b) OTMIC with losses optimization, (c) conventional DQ control and (d) harmonic spectrum of phase current (left) and torque (right) for every method [18].
Figure 9. Experimental results of stator current and torque of: (a) OTMIC with torque ripples optimization, (b) OTMIC with losses optimization, (c) conventional DQ control and (d) harmonic spectrum of phase current (left) and torque (right) for every method [18].
Energies 13 06184 g009
Figure 10. Motor efficiency for OTMIC and DQ control.
Figure 10. Motor efficiency for OTMIC and DQ control.
Energies 13 06184 g010
Table 1. Air gab winding wheel-hub motor parameters [17,18].
Table 1. Air gab winding wheel-hub motor parameters [17,18].
SymbolDescriptionValueUOM
Motor Parameters
u D C DC voltage48 V
L s + M s phase inductance1.5 μ H
R s phase resistance0.026 Ω
t M motor time constant58 μ s
p number of poles94
k M motor constant0.304Nm/(TA)
b _ k B odd harmonics[1.15 0.2 0.06 0.01] T
c constant of coulomb friction and hysteresis loss0.0832 N m
d constant of mechanical friction and eddy loss0.0008 N m · s / r a d
Control Parameters
K P Proportional term0.166 V / A
K I Integral term0.026 V / A
K D Derivative term0.013 V / A
N D PID filter term0.9854
T S sensor time constant1 μ s
Δ t sampling time10 μ s
f P W M switching frequency100 k H z
T R e q requested time constant20 μ s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kasper, R.; Golovakha, D. Combined Optimal Torque Feedforward and Modal Current Feedback Control for Low Inductance PM Motors. Energies 2020, 13, 6184. https://doi.org/10.3390/en13236184

AMA Style

Kasper R, Golovakha D. Combined Optimal Torque Feedforward and Modal Current Feedback Control for Low Inductance PM Motors. Energies. 2020; 13(23):6184. https://doi.org/10.3390/en13236184

Chicago/Turabian Style

Kasper, Roland, and Dmytro Golovakha. 2020. "Combined Optimal Torque Feedforward and Modal Current Feedback Control for Low Inductance PM Motors" Energies 13, no. 23: 6184. https://doi.org/10.3390/en13236184

APA Style

Kasper, R., & Golovakha, D. (2020). Combined Optimal Torque Feedforward and Modal Current Feedback Control for Low Inductance PM Motors. Energies, 13(23), 6184. https://doi.org/10.3390/en13236184

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