Correlation of EM Torque Settling Time with Shaft Velocity Response Time

## 1. Introduction

As already stated for the reasons given in a previous chapter a good continuous time model, of low complexity, of a BLMD system is essential to adequately describe mathematically the PWM inverter switching process with dead time and subsequent binary waveform generation in terms of the switching instant occurrences for accurate computer aided design (CAD) of embedded BLMD model simulation in proposed electric vehicle (EV) propulsion systems. In this chapter a complete software model of the BLMD system as a set of difference equations representing subsystem functionality, the organization of these subsystem activities into flowchart form and the processing details of these modular activities as software function calls in C-language for simulation purposes (Guinee, 2003) is presented.

Furthermore in this the second chapter, concerning BLMD model fidelity for EV applications, BLMD model simulation accuracy for embedded EV CAD is next checked for a range of restraining shaft load torques via numerical simulation and then extensively compared and benchmarked for accuracy against theoretical estimates using known manufacturer’s catalogued specifications and motor drive constants (Guinee, 2003).

Model simulation accuracy is further substantiated and validated through evaluation of the shaft velocity step response rise time when cross checked against (i) experimental test data and (ii) that evaluated from the catalogued performance index relating to the brushless motor dynamic factor (Guinee, 2003). Numerical simulation with outer velocity loop closure is used to demonstrate the accuracy of the completed BLMD reference model, based on established model confidence in torque control mode, in ASD configuration when compared with experimental test data.

In addition to the BLMD model structure presented in the previous chapter for actual drive emulation two innovative measures which relate to increased drive performance are also provided. These novel techniques (Guinee, 2003), which include

inverter dead time cancellation and

motor stator winding impedance angle compensation,

are encapsulated within the BLMD model framework and simulated for validation purposes and prediction of enhanced drive performance in EV systems. An approximate analysis is given to support the approach taken and verify the performance outcome in each case.

In the first of these BLMD performance enhancements a novel compensation method has already been presented in the first chapter to offset the torque reduction effects of inverter delay during BLMD operation. This simple expedient relies on the zener diode clamping of the triangular carrier voltage during the carrier waveform comparison with the modulating current control signal in the comparator modulator to nullify power transistor turnon delay. This approach obviates the need for separate compensation timing circuitry in each phase as required in other schemes. The accuracy of this methodology is supported by current feedback, EM torque generation and shaft velocity trace simulation when compared with similar traces from the BLMD benchmark reference model with the effects of the inverter basedrive trigger delay neglected.

The second proposed innovative improvement, presented in this chapter, relates to the progressive introduction of commutation phase lead with increased shaft speed as BLMD impedance angle compensation which forces the impedance angle to the same value as the internal power factor angle. This effect maintains zero load angle between the stator winding terminal voltage and the back emf. It also results in rated load torque delivery at lower shaft speeds with minimal rise time, overshoot and settling time in the generated torque for a range of torque demand input values. This novel technique greatly enhances the dynamic performance of the embedded BLMD prime mover in EV applications without overstressing mechanical assembly components during periods of rapid acceleration and deceleration. The incorporation of this novel impedance angle compensation technique thus minimizes component wear-out such as gear boxes, transmission shafts and wheel velocity joints and consequently enhances overall EV reliability improvement. BLMD simulation is provided in torque control mode at rated torque load conditions, for the actual drive system represented, with and without impedance angle compensation to gauge model performance accuracy over a range of torque demand step input values.

## 2. BLMD model structure and program sequence of activities

The BLMD model structure is composed of interconnected subsystems with feedback as shown in Figure 1, of varying complexity according to physical principles. Consequently it can be described by a discrete time configuration of first order digital filter realizations for linear elements cascaded with difference equations representing nonlinear PWM inverter behaviour into a complete software model for simulation purposes as illustrated in Figure 2. The BLMD model program is organized into a sequence of software activities, coded in C-language as function calls, representing the functionality of various subsystem modules shown as the flowchart in Figure 3. All subsystem output (o/p) variable quantities in the cascaded activity chain are assumed to remain constant, once computed irrespective of feedback linkage, throughout the remainder of the time step interval t_{k} based on the simulation sampling rates (1/t_{k}) chosen from considerations given in section 3.1 of the previous chapter. The essential features of the BLMD model program in Figure 3 can be explained by means of the linked modular software configuration encoded as the functional block sequence in Figure 2 along with the appropriate C-language code segments illustrated in Figure 4.

### 2.1. BLMD model simulation

Numerical simulation commences with the declaration of known BLMD system parameters followed by a declaration with initialization of variables and three phase (3) arrays for global usage, over the program linked function call sequence, as outlined in the code blocks shown in Figure 4-1.

All first order linear system discretization is accomplished by complex variable substitution of the Euler backward rectangular rule using the Z transform. The alternative filter discretization process using Tustin’s bilinear method (Franklin et al, 1980) or the trapezoidal integration rule (Balabanian, 1969) can also be used but with negligible observable differences at the step size t chosen. Resultant digital filter implementation for simulation purposes is facilitated by transfer of the appropriate filter time constant and gain coefficients using the C-language ‘structure’ mechanism in the function call * setup_fo_filt* ( ) illustrated in Figure 4-2.

Before proceeding with model program execution in the ^{th} time step (t_{k-1}
_{k}) all global variables and arrays are captured in the function call * capture_filt_out* ( ) for later reinstatement, during accurate resolution of the width modulated pulse edge transition time via the regula-falsi method, with restore_filt_out ( ) in Figure 4-3.

The instruction code group * run_to_pwmsw* ( ) processes the sequence of BLMD software activities up to the comparator modulator o/p using the following list of function calls in Figures 4-1and 4-4.

* tdemf* ( ) filters the i/p torque demand signal

*with o/p*torq_dem

*.*ftorq_dem

* mot_commutator* ( ) establishes the 3-phase reference

*[j], from the computed shaft rotor displacement*psink

^{k}

_{m}, for 3 phase stator winding voltage commutation with modulated amplitude.

* tor_sink*[j] based on the filtered torque demand

*.*ftorq_dem

* idemf* ( ) filters the i/p torque related current command signal

*[j] with o/p*tor_sink

*[j].*idemk

* cur_cont* ( ) The lag compensator ‘optimizes’ the current error as the difference between the current command

*[j] and the filtered stator winding current feedback*idemk

*[j] in each phase of the 3 current control loop. The o/p*fifbk

*[j] from each of the high gain controllers is amplitude limited to the saturation voltage levels V*vmpwmk

_{z}(~10v) by zener diodes.

* pwm_mod* ( ) produces a width modulated o/p pulse sequence

*[j] for each phase in accordance with the amplitude comparison of the modulating control signal*pwmtrk

*[j] o/p and the triangular dither signal*vmpwmk

*. The modulator has a gain K*osc

_{mod}with complementary outputs, for basedrive operation, that are hard limited to V

_{z}(~10v) by zener diodes.

A test is used to interrogate the o/p status of the simulated comparator modulator by monitoring any observational sign change in the o/p polarity (V_{Z}), which is indicative of a modulated pulse edge transition, after execution of the software code module * pwm_mod* ( ) as shown in Figure 4-5.

The o/p status of the comparator modulator is examined by comparing the trapped value * cpwmtrk*[j] at the beginning of the time step t

_{k-1}with the new o/p

*[j] at t*pwmtrk

_{k}for each phase j and signalling any change via the

*in the function call*pwm_test_flag

*( ) detailed in Figure 4-6. If a crossover event occurs during simulation then the transition interval t*test_pwm_xover

_{X}= (t

_{X}-t

_{k-1}), denoted by

*, is determined by the regula falsi method in (LXVII) in the previous chapter.*min-time

A check is made to see if this value occurs within the imposed tolerance limit (* tol**

*) at the end of the time step interval t*delt

_{k}, denoted by

*, in which it is discarded in the affirmative without test flag registration. If multiple crossover events occur within the simulation interval, corresponding to different phases, then all transition times with the appropriate phase tag number are logged in the respective*delt

*[j] and*sw_time

*[j] arrays along with the signaled transition count via the PWM test flag. A check is also performed for identical multiple switch transition times without an increase in test flag count. The test routine is completed by arranging the multiple switch times, with corresponding phase listing, in increasing order of magnitude for subsequent detailed PWM simulation in the function call*test_flag

*( ). Accurate internal simulation of a modulated pulse transition, indigenous to the time step, commences with restoration of the captured global variables preceding the time step and temporary storage of the original step size (*inter_pwm_simulation

*) and time delay (*zeit

*) settings, relevant to the*t_del

*( ) signal source, for later retrieval. The new time step (*dither

*) is initially set to the smallest switching interval*delt

t

^{1}

_{X}(

*[1]) for discretization of all first order linear subsystems using the call function*sw_time

*( ) as per the C-code module in Figure 4-7.*setup_fo_filt

The necessary delay offset _{d} (* tdel*) is determined by back tracking (t

_{k}-

t

^{1}

_{X}) from

t

_{k}for proper time registration in the execution of the carrier function

*( ) and rerun of the call sequence*dither

*( ) followed by the block function call*run_to_pwmsw

*( ). The post PWM simulation call list contains the additional BLMD model basedrive switching features as an embedded layer in the nested*run_post_pwmsw

*( ) and associated switch event signalling*base_drive

*( ) program routines. The complete BLMD model program is subsequently exercised for othermultiple switch time intervals t*test_drk_xover

^{i>1}

_{x}with updated linear system discretization and adjusted delay offset. Termination of the remainder of the original time step simulation is accomplished by setting the integration interval

*equal to the time step residue (t*delt

_{k}-t

^{max}

_{X}) followed by the call sequences

*( ) and*run_to_pwmsw

*( ) and exiting to the main program with a reinstatement of original time settings*run_post_pwmsw

Numerical BLMD model simulation proceeds to the next program step in the flowchart cycle shown in Figure 3, by processing the call sequence * run_post_pwmsw* ( ), with the execution of the switch event routine

*( ) associated with the basedrive turn-on/off as a consequence of the PWM process. The relevant global variables and arrays associated with this call sequence run are trapped by the command*base_drive

*( ) as a precursor to basedrive simulation. The ‘lockout’ circuit routine, illustrated in Figure 4-8, consists of the integrating capacitor action when the PWM comparator o/p*capture-drk_out

v

^{k}

_{sm}[

*] exceeds*j

v

^{k-1}

_{lj}[j] and charge dumping when

v

^{k}

_{sm}[

*] <*j

v

^{k-1}

_{lj}[

*] as shown in Figure 3 for the basedrive BDJ with a similar microprogram description for complementary basedrive*j

V

_{th}(0.0) setting determines the inverter turn-on time. The charge dump action by the shunt diode across the delay timing resistor is virtually instantaneous when the switched comparator PWM output v

^{k}

_{sm}[j]=K

_{mod}(v

^{k}

_{c}[j]-v

^{k}

_{tri}) is hard limited to -v

_{z}.

This effect results in swift basedrive turn-off with zero delay when referenced to the trailing edge of the PWM o/p. However the capacitor discharge can be gradual, when the PWM o/p is soft switched (v_{z} > v^{k}_{sm}[j] > -v_{z}), due to the limited magnitude of the product combination of modulator gain _{osc} (~68) and error response v^{k}_{c}[j] of the current loop controller which is implicitly dependent on the filtered current feedback response ^{k}_{f}[j] for fixed current demand ^{k}_{d}[j]. The gradual reduction in capacitor voltage protracts the basedrive switch-off time, when referenced to the initial point of the logic “1-to-0” transition, associated with the PWM trailing edge. This delay has to be accounted for in an accurate inverter software model description with a search of the basedrive turn-off in addition to the turn-on times associated with exponential voltage growth.

The BLMD program test function * test_drk_xover* (&

*), which is shown in Figure 4-10and is very similar to*drk_test_flag

*( ) in code content, checks for basedrive on/off firing signal occurrence within a simulation time step interval. This search is complemented with the evaluation of associated multiple phase activation times t*test_pwm_xover

^{i}

_{x}for both normal BDJ and complementary

^{i}

_{x}are determined by piecewise linear approximation using (LXXVI) in the previous chapter, ranked in increasing order of magnitude and phase tagged via the global storage arrays

*[j] and*drk_sw_time

*[j] for subsequent use in detailed basedrive simulation. Accurate simulation of the basedrive trigger timing signals for subsequent inverter operation is achieved using the software routine call*drk_phase_flag

*(&*inter_base_drk_loop_sim

*) which is shown in Figure 4-11and has similar execution features to*drk_test_flag

*( ).*inter_pwm_simulation

The function call begins with the reinstatement of the global arrays at the beginning (t_{k-1}) of the time step using * restore_drk_out* ( ), illustrated in Figure 4-9, and temporary storage (

*) of the original step size t. The routine proceeds with linear system discretization appropriate to and with execution of the*zeit

*( ) function and the subsequent call sequence run_post_drksw ( ), listed in Figure 4-12, for progressive substitution of multiple differential switch times as the temporary variable*base_drive

*.*delt

This simulation call is completed with restoration of the original time step size followed by first order system discretization with a return to the main BLMD program to begin the new time step t_{k} t_{k+1}. The function call group run_post_drksw ( ), summoned during main program execution in the flowchart of Figure 3, processes the following sequence of modular software activities illustrated in Figures 4-12and 4-13 pertaining to BLMD system electrodynamic operation with inverter interaction.

* pwm_inv* ( ) generates the 3 inverter output HT binary voltage

*[j] v*ph_htk

^{k}

_{g}[j] in response to the PWM basedrive gating signals

_{th}establish the conduction states

V

_{O}, of the complementary power transistor pair T

_{J+}and T

_{J-}in each leg J of the 3 inverter shown in Figure 15 in the previous chapter. The tristate flag condition in conjunction with the sustained stator winding current flow

*[j]*ph_curk

i

^{k}

_{s}[j] through the free wheeling shunt diodes establish the inverter o/p binary voltage as 0 or

V

_{dc}. The neutral star point voltage

*(*Vng

v

_{sg}) of the stator winding is determined from (LVIII) in the previous chapter for subsequent evaluation of the phase voltages

*[j] v*ph_voltk

^{k}

_{s}[j] via (LIX) in the previous chapter.

If one or more phase currents are zero during the condition V_{O} = 0 a * tf* is incremented by unity for each null occurrence and the relevant phases are tagged using the

*[j] array for later identification in the computation of the relevant motor winding phase voltages. A test flag count of 2 or 3 indicates that the phase voltages are identical to the back*test_cur

*voltages*EMF

*[j] v*back_emfk

^{k}

_{e}[j] from the previous simulation step using expressions (XXIII) and (LVII) in the previous chapter. If the test flag (

*) value is unity the associated phase current*u

*[*ph_curk

*] is zero and the following relations result*test_flag

which enable the phase voltages to be calculated from the 3 inverter o/p.

* winding* ( )determines the stator winding current

*[j] from a knowledge of the phase voltages and motor parameters {r*ph_curk

_{s}, L

_{s}} encoded in the discrete filter representation of the winding electrical behaviour.

* convert_torqunit* ( )computes the developed electromagnetic torque, manifested in the winding current flow, via the 3 commutation vector

*[j] and motor torque constant K*psink

_{t}for dynamic operation.

* mot_shaft_vel* ( )evaluates the rotor shaft velocity

mot_shaft_vel

^{k}

_{m}from the net torque

tot_torq

^{k}

_{e}available, with the effects of load torque retardation

torq_load

^{k}

_{l}considered, using a discrete filter representation of the rotor dynamics as per (IL) in the previous chapter with parameters {

J

_{m},

J

_{l},

B

_{m}}. The back emf can be determined via the motor voltage constant Ke, along with the mechanical (

*) and electrical (*mech_power

*) power delivery, once the shaft velocity is known.*elec_power

* mot_shaft_pos* ( )evaluates the motor shaft position

*from the rotor velocity.*mot_shaft_pos

## 3. BLMD model simulation with restraining shaft load torque

The effect of a fixed applied shaft load _{l} on BLMD model behaviour can be monitored via its simulation characteristics, in torque control mode with and without impedance angle compensation considered, for a range of torque demand step input stimuli capable of matching the posed restraining torque. A suitable choice for the target load magnitude is based on the manufacturers continuous rated stall torque of 5Nm for the particular motor type specified in Table 1 of the previous chapter. The set of motor shaft velocity step response characteristics _{d} range is displayed in Figures 6 and 7 without impedance angle compensation.

The response time _{res} required for the motor drive to reach the maximum shaft velocity of

All the torque response characteristics, with the exception of that at _{d} = 4v, exhibit overshoot before settling to the required value of ~5.3Nm to overhaul the fixed restraining load torque (5Nm) and frictional effects. The degree of overshoot increases in proportion to the torque demand i/p, as exhibited in Figure 10 for the average peak EM torque _{ep} responsible for overshoot _{ep}, accompanied by a corresponding reduction in settling time as shown in Figure 11.

### 3.1. Theoretical consideration of motor accelerative dynamical performance

The reduction in settling time is paralleled by the shaft velocity response time improvement in reaching rated motor speed. It is evident from inspection of the velocity and torque simulation traces that a direct correlation exists between the EM torque settling time and motor shaft velocity response time as indicated in Table 1.

Torque Load =5Nm_{l} | “Inertial” Time Constant =_{m}/_{m} =0.318_{m} | Tran. Gain (Fig.10) K _{(} =1.28 | ||||

Torq-Dem_{d}(Volts) | Av. Peak EM Torq ep (Nm) | Torque Overshoot -l_{ep}(Nm) | Torq. Settling Time (sec)_{setl}Figs. 8/9 | Shaft Velocity Rise Time (sec)_{res }Figs. 6/7 | Theoretical Rise Time (sec)_{r}Eqn. (IV) | Rise Time T _{(}(sec) via Dyn-Fac. Eqn (VI) |

5 | 6.2 | 1.2 | ~0.13 | ~0.13 | 0.131 | 0.107 |

6 | 7.45 | 2.45 | ~0.06 | ~0.06 | 0.057 | 0.0524 |

7 | 8.98 | 3.98 | ~0.04 | ~0.037 | 0.034 | 0.0323 |

8 | 10.29 | 5.29 | ~0.03 | ~0.027 | 0.025 | 0.0243 |

9 | 11.634 | 6.634 | ~0.024 | ~0.022 | 0.02 | 0.02 |

The shaft velocity step response rise time, as defined in Figure 6, can be obtained directly from the solution of the transfer function (XCIX) from the previous chapter in the time domain with a step input approximation for the average peak torque overshoot

with time constant

The step response time, for the shaft velocity under load conditions to reach maximum speed

The estimated rise times are in excellent agreement with the approximate settling and response times obtained from the BLMD model simulation traces. An alternative crude estimate of the response time can be obtained from the motor “dynamic factor”

for average peak torque endurance as the acceleration time

from standstill to maximum speed assuming a shaft velocity linear transient response which is valid for torque demand values in excess of 5 volts.

These response estimates, given in above Table 1, are in good agreement with those already obtained except for that at _{d} = 5v where the rise time is longer with exponential speed ramp-up.

### 3.2. Torque demand BLMD model response - internal node simulation

The simulated back-EMF along with the stator impedance voltage drop are illustrated in Figures 12and 13 for two relatively close values of torque demand i/p. In the former case the torque demand i/p of 4volts results in sufficient motor torque to meet the imposed shaft load constraint (5Nm) without reaching rated speed and saturation (10v) of the current compensator o/p trace shown in Figure 14. The corresponding reaction EMF exceeds the winding impedance voltage V_{Z} and is almost in phase with the stator current, which is proportional to V_{Z}, at the particular low motor speed reached. The torque demand i/p of 5v in the latter case results in the onset of a clipped current controller o/p in Figure 15 due to saturation (10) at rated motor speed _{rmax}.

The back-EMF generated at this speed greatly exceeds the winding impedance voltage, as in the former case, and leads the stator current necessary to surmount the torque load by the internal power factor (PF) angle (~27°) with a correspondingly low power factor (~0.7).

The stator winding currents corresponding to the inputs _{d} >5v.

The simulated motive power characteristic with the steady state threshold value of ~2.3kW necessary to sustain shaft motion, for _{d} =5v with restraining load torque and friction losses is shown in Figure 18 at base speed _{rmax} 420 rad.sec^{-1}.

This can be rationalized from the power budget required to sustain the load torque at rated speed via (LXXXVIII) in the previous chapter as

The excess coupling field power required to surmount mechanical shaft friction losses is shown simulated in Figure 19 with a steady state estimate of ~200 watts.

Torq_Dem d Step i/p | Shaft_Vel rmax rad.sec ^{-1} | Elec_Power P_{e} volts (XLVII) – Prev. Chap | Back_EMF V_{ej}volts | Imped_Vol V_{Z}_{ }volts(XC) – Previous Chap | Ph_Cur I _{js }amps | ||

5v | 419.2 | 2301 | 94.82 | 44.87 | 9.09 | ||

6v | 420.3 | 2305 | 95.24 | 48 | 9.7 | ||

7v | 418.9 | 2298 | 94.78 | 50.89 | 10.32 | ||

8v | 410.3 | 2251 | 92.05 | 52.42 | 10.84 | ||

9v | 405.5 | 2224 | 91.5 | 53.66 | 11.23 | ||

Derived Phase Quantities as per Figure 42 in previous chapter | |||||||

Torq_Dem d volts | Int. PF Ang I (XCII) – Prev. Chap. | _{I} Estimatevia Figure 13 | Ph_Vol V _{js}(XCIII) – Prev. Chap | Imp_Ang Z (LXXXIV) – Prev. Chap. | Load Ang_{T}(XCV) – Prev. Chap | PF Ang_{I} + _{T} | |

5 | 27.13 | 27.13 | 126.44v | 81.26 | 15.52 | 42.65 | |

6 | 33.7 | 32.16 | 131.75v | 81.28 | 15.6 | 49.3 | |

7 | 38.43 | 38.74 | 136.56v | 81.26 | 14.68 | 53.11 | |

8 | 41.24 | 42.58 | 136.5v | 81.08 | 13.83 | 55.07 | |

9 | 43.8 | 45.12 | 138.1v | 80.97 | 13.58 | 57.38 |

The effect of shaft load on the BLMD model simulation characteristics for _{d} >5v is summarized in above Table 2 for steady state conditions with the aid of the general phasor diagram in Figure 42 of the previous chapter.

It is evident from the table that the back EMF has reached its peak rms value with the onset of maximum shaft velocity, for all values of _{d} >5V, with

Furthermore the impedance voltage drop V_{z} in (XC) of the previous chapter is limited to a very small increase with torque demand current I_{dj} listed in Table 2 and is shown almost stabilized to a constant value in Figure 20. This voltage clamping effect, due to current compensator o/p saturation in response to tracking current feedback, is controlled to achieve the desired rms level of clipped current flow in the stator winding as shown in Figure 17 to satisfy torque load requirements. The rms winding current flow necessary at unity internal power factor to meet steady state toque load and friction demands at ~5.4Nm in Figures 8 and 9 can be determined from (XLV) in the previous chapter as

This is almost identical to the rms values obtained from BLMD model simulations in Table 2, which are consistent with increased torque current demand, when internal power factor self adjustments are accounted for as in

The internal power factor angles, listed in Table 2 and displayed in Figure 21, are deduced for _{d} >5v from the mechanical power transfer by substituting the rms quantities obtained from back EMF and winding current simulations in expression (XCII) of the previous chapter. These angles, which increase with torque demand i/p, can be alternatively calculated from the simulated winding current response using (X) with knowledge of_{z}, which is fixed at rated machine speed _{rmax}, is determined from (LXXXIV) as ~81.2 in Table 2.

The rms winding voltage V_{js} is obtained in its pure spectral form, instead of the PWM version furnished by the current controlled inverter, upon application of (XCIII) to the known rms phasor quantities given in Table 2 for different values of _{d} >5V.

Knowledge of the relevant phasor magnitudes with corresponding phase angles enable the load angle _{T} to be determined from (XCV) of the previous chapter for given shaft load conditions. This is approximately fixed, at ~15 as indicated in Table 2 with about 2 variation, over the torque demand i/p range as shown in Figure 21. The resulting power factor angle listed in Table 2 increases with _{I}, for fixed load angle over the torque demand i/p range as shown, in a way that is commensurate in (X) with motor current requirements towards sustaining shaft load torque with a decreasing power factor as illustrated in Figure 22.

### 3.3. BLMD model simulation with novel impedance angle compensation

The effect of motor impedance angle compensation (MIAC), manifested as commutation phase lead angle incorporated into the BLMD model in (XCVIII) of the last chapter as

The variation of peak torque overshoot with i/p demand, displayed as the mutual characteristic in Figure 25, is linear with a transfer gain that is lower than that without MIAC in Figure 10. Consequently the maximum peak torque delivery, for a given i/p demand to sustain shaft load requirements, is lower in amplitude and of shorter overshoot pulse duration as seen in Figure 24 when compared with that without MIAC in Figures 8 and 9. Furthermore the persistence of torque overshoot is lower with a much reduced settling time (<0.015 sec), in reaching steady state sustained load conditions in all cases albeit at lower acceleration and much smaller drive speeds, thereby exerting less mechanical stress on the drive shaft components and minimizing shaft flexure in EV propulsion applications.

The shaft velocity characteristics also indicate a much lower steady state motor run speed, with MIAC deployed, which never reaches velocity saturation _{dmax} =10V. This speed reduction is singly due to the maintenance of an almost zero load angle _{T} shown in Figure 27, between the motor terminal V_{js} and back EMF V_{ej} rms voltage phasors in Figure 45 of the previous chapter, by commutation phase angle advance for optimal torque production as indicated from the BLMD simulation results in Table 3.

This phase compensation technique results in back EMF and winding impedance voltage V_{z} phasors that appear approximately equal in magnitude over the allowable torque demand input range as shown in Figure 28. Furthermore the internal power factor angle _{I} is forced to adopt approximately the same value as the machine impedance angle _{z} as indicated in Table 3, by the phase advance measure _{z} in the current commutation circuit, with a consequent collinear alignment of phasors V_{ej} and V_{z} in Figure 45. This collinear arrangement can only be sustained at a particular machine speed that is dependent on the torque demand i/p which determines the subsequent winding current flow and thus the necessary impedance angle for alignment. This reasoning can be deduced as follows by noting that for a given torque load _{l} the rms winding current flow is linear with torque demand i/p as per Table 3 and Figure 29.

Torq_Dem d Step i/p | Shaft_Vel rmax rad.sec ^{-1} | Elec_Power P_{e} (XLVII) in Prev. Chap. | Back_EMF Vej volts | Imped_Vol V_{Z }volts – (XC) in Prev. Chap. | Ph_Cur Ijs amps |

4v | 18.6 | 94.44 | 4.17 | 6.06 | 7.76 |

5v | 48.95 | 257.2 | 11.5 | 9.23 | 9.7 |

6v | 70.87 | 363.67 | 16.01 | 13.05 | 11.71 |

7v | 87.9 | 452.6 | 19.95 | 17.33 | 13.66 |

8v | 102.9 | 531.2 | 23.28 | 22.18 | 15.7 |

9v | 116.3 | 602.2 | 26.3 | 27.45 | 17.74 |

Torq_Dem d volts | Int. PF Ang I (XCII) in Prev. Chap. | Ph_Vol V_{js}(XCIII) in Prev. Chap. | Imp_Ang Z (LXXXIV) in Prev. Chap. | Load Ang _{T}(XCV) in Prev. Chap. | PFAng_{I} + _{T} |

4 | 13.75 | 10.23 | 16.1 | 1.39 | 15.14 |

5 | 36.13 | 20.73 | 37.22 | 0.51 | 36.64 |

6 | 49.71 | 29.06 | 47.72 | -1.0 | 48.71 |

7 | 56.38 | 37.27 | 53.76 | -1.22 | 55.16 |

8 | 61.02 | 45.44 | 57.95 | -1.5 | 59.52 |

9 | 64.52 | 53.73 | 61.01 | -1.79 | 62.73 |

#### 3.3.1. MIAC substantiation via theoretical analysis and validation

The internal power factor angle _{I} can be determined theoretically for fixed winding current flow corresponding to a given torque demand i/p using (Eq. 9) and (Eq. 10), assuming negligible dynamic friction at the shaft speeds concerned with

The motor terminal voltage i/p V_{js} from previous chapter can be optimized with respect to the motor impedance angle _{z}, which is unknown, in terms of the rms phasor quantities V_{ej}, V_{z} and the fixed internal power angle _{I} from (Eq. 11) by letting

This procedure results in the impedance angle _{z} in terms of the known angle _{I} as

which is unknown as both V_{ej} and V_{z} depend on the motor shaft velocity _{r}. The shaft velocity can now be determined from (LXXXIV) from previous chapter using expression (Eq. 13) as

Torq_Dem d Step i/p | Ph_Cur Ijs Table 3 | Int_PF_{I}Eqn. (XI) | Shaft_Vel r Eqn (XV) | Back_EMF_{ej}Eqn (VIII) | Imp_Vol _{Z}Eqn. (XC) in Prev. Chap. | Ph_Vol_{js}Eqn. (XIV) |

4v | 7.76 A | 15.37 | 17.71 rad/s | 3.94v | 6.04v | 9.98v |

5v | 9.70 A | 39.52 | 53.17 rad/s | 11.84v | 9.43v | 21.27v |

6v | 11.71 A | 50.28 | 77.55 rad/s | 17.27v | 13.74v | 31.01v |

7v | 13.66 A | 56.79 | 98.43 rad/s | 21.92v | 18.71v | 40.63v |

8v | 15.70 A | 61.54 | 118.87 rad/s | 26.48v | 24.71v | 51.19v |

9v | 17.74 A | 65.05 | 138.49 rad/s | 30.85v | 31.54v | 62.39v |

This value of _{r} can be used to theoretically generate the rms voltage phasors V_{ej}, V_{z} and V_{js} using expressions (VIII), (XC) and (XCIII) in the previous chapter respectively from a knowledge of the motor winding phasor current I_{js} as per Table 1V over the i/p torque demand range range_{T} from (XCV) in the previous chapter due to the phasor collinearity and thus improved torque control via the PWM voltage supplied by the current controlled inverter. The power factor angle , internal power factor angle _{I} and machine impedance angle _{z} variations with torque demand i/p, which are displayed in Figure 27 using estimates extracted from BLMD model simulation in Table 3 for_{T} 0).

The internal power factor _{I} adjustment, from the mirrored motor current increase in Figure 29, constrained by a fixed shaft load in (Eq. 10). Impedance angle compensation results in a improved motor power factor as shown in Figure 31 than that without MIAC over the torque demand i/p range _{l}.

Motor speed reduction is also mirrored with a decrease of the shaft velocity step response rise time as shown Figure 32 with maximum values falling below the velocity time response floor of the uncompensated BLMD model. This results in constant motor speed operation, though small by comparison to that without phase angle advance, well below the rated value in torque control mode with smooth torque delivery to satisfy load requirements.

The simulated motor winding impedance and back EMF voltages for mid (5V) and full range (9V) torque demand input values, which result in developed torque capable of surmounting the fixed restraining shaft load (5Nm), are displayed in Figures 33 and 34. Both sets of characteristics exhibit comparable amplitudes appropriate to the level of torque demand i/p, with speed related motor current phase lags _{I} as per Table 3, that are much lower than those without MIAC in Figure 13. The impedance and back EMF voltages are interrelated which can be shown as follows by starting with expression (XC) for V_{z} and using (Eq. 9) and (Eq.10) giving

This can be rewritten by using (LXXXIV) in the previous chapter with optimized value of _{I} in (Eq. 13) as

The shaft velocity _{r}, linking the back EMF, can be replaced in (Eq. 18) by using (Eq. 13) yielding

where

and

from substitution of parameters in Table 1 of the previous chapter and _{l} = 5Nm. The impedance voltage in (Eq. 17) is expressed as a quadratic equation in terms of the back EMF with points of equality corresponding to

which are visible in Figure 28 as points of intersection of the two voltage traces.

These crossover points divide the rms V_{z} amplitude variation along with V_{ej} in Figure 28 into three distinct regions, over the usable torque demand i/p range as per Table 1V, with

These regions can also be inferred from the voltage amplitude traces in Figure 33 and 34 where V_{ej} exceeds V_{z} in the former case with

The simulated stator winding current along with current feedback response and current controller o/p are displayed in Figures 35, 36 and 37 respectively for

## 4. BLMD reference model simulation in velocity control mode

In this section the BLMD reference model performance as an ASD emulator is examined and compared with experimental step response data for shaft inertial load conditions with J_{l} ~ 3J_{m}. Adjustable speed drive operation, with embedded inner PWM current control, is effected by closing the outer velocity loop via a two term PI term controller G_{v} as shown in Figures 1 and 5. The analog velocity controller shown in Figure 41, which has an inbuilt velocity offset adjustment and speed gain control adjustment K_{s} for the chosen BLMD system modelled here (Moog GmbH, 1989) has a transfer function

with proportional and integral compensation gain settings K_{p} and K_{I} respectively. The inclusion of this outer loop velocity compensator, in addition to the inner torque control current loop, results in a complete holistic BLMD reference model that can now be used for ASD simulation and performance evaluation in embedded applications. Proportional and integral control is easily incorporated in C-language routine during BLMD simulation of velocity closed loop operation as a digital filter code module via (LXV) of the previous chapter using the backward Euler method in (LXIV) of the previous chapter. The resulting ASD model was exercised at low and high shaft velocities corresponding to 36.4% and 73.6% of rated motor speed n_{o} in Table 1 of the previous chapter and compared with experimental test data at critical internal nodes in Figure 1 for model validation and simulation accuracy.

The current controller G_{I} step response simulation traces for a velocity step command input V_{} of 2volts, corresponding to 36% of rated motor speed (~4000rpm), are exhibited in Figures 42 to 44 for linear pulsewidth modulator operation.

The accuracy of these simulation traces, which capture the essence of the velocity transient step response V_{r} overshoot in Figure 45, is characterized by a large waveform correlation coefficient of fit in Table 5 which provides a good indication of the model fidelity when matched with experimental data.

Target Data Length =3000_{D} | Data Sampling Rate: 12.5kHz | BLMD Simulation Time Step: 1s |

Waveform Correlation Analysis for Total Inertial Shaft Load= _{Tot } + _{l}. =12.3kg.cm_{m}^{2} | ||

ASD Waveform | Velocity Command i/p V _{( }= 2V | Velocity Command i/p V _{( }= 4V |

Current Demand _{dj} | 96.8% | 92.85% |

Current Feedback (FC)_{fj} | 97.26% | 93.27% |

Current Compensator output _{cj} | 59.81% | 45.46% |

Motor Shaft Velocity output _{(r} | 99.8% | 99.68% |

The simulated current demand and feedback waveforms, which have high matching coefficients with test data, exhibit an amplitude modulated step response with velocity transient overshoot and ringing, before eventually setting to negligible constant amplitude traces with fixed frequency commensurate with reached shaft speed _{r} demanded (V_{r} ~2V) in Figure 45.

The compensated velocity error output for 2Volts operation V_{r} in the BLMD network structure in Figure 5 is equivalent to the filtered torque demand _{df}, as the velocity control effort V_{e} shown in Figure 46, applied to the inner closed loop for motor current control and BLMD output torque regulation. This optimized velocity error V_{e} in Figure 46 is a short duration pulse for reasons of fast BLMD shaft velocity risetime T_{res} and short setting time T_{setl} as required in high performance ASD industrial applications.

The presence of overshoot in the BLMD velocity step response in Figure 45 is due to the non-optimal tuning of the velocity controller PI parameters required to ensure stiff dynamical operation for the total drive shaft inertial load J_{Tot} = J_{m}+J_{l} in question.

Examination of the ASD velocity step response trace simulations over a range of shaft inertial load multiples of the rotor value J_{m} in Figure 47 reveal that the PI parameters have been optimized only at zero load with J_{Tot} = J_{m} for good drive dynamic transient performance with little overshoot. The effect of increased shaft inertia on the velocity control effort V_{e} in Figure 48, for a velocity command input of 2 volts, is a greater sustained oscillation accompanied by longer settling times T_{setl} manifested in the simulated ASD velocity step response as shown in Figure 49. This behaviour is mirrored by an increased overshoot, as defined in Figure 47, in the BLMD shaft velocity step response with shaft load inertia as shown in Figure 50.

The effect of increased load inertia on ASD dynamic performance also translates into slower rise times T_{res} as shown in Figure 51 for a non optimally tuned velocity controller.

Further ASD step response simulation and comparison with experimental measurements in Figures 52 to 55, for a 4 volts velocity command input which corresponds to 74% of rated motor speed n_{0} with resulting saturated pulsewidth modulator operation in Figure 54 for the load inertia considered (~12.3kg.cm^{2}), reveal good BLMD model accuracy.

The quality of ASD simulation trace match with test data is indicated by the high value of the waveform correlation coefficients given in Table 5.

The velocity loop derived torque command input stimulus V_{e} in Figure 56 has a pulse duration that is much shorter than the time constant (_{ml} ~ J_{tot}/B_{m}) of the load dynamics, eventhough the pulse amplitude is of sufficient strength to force the shaft velocity to the value demanded (V_{}~4volts). The endurance ^{*} of the velocity control effort in Figure 56, associated with pronounced PWM saturation in Figure 54, is a measure of the maximum sustained EM torque necessary to accelerate the total BLMD inertial load masses to the appropriate shaft velocity demanded by the ASD command setting V_{}. This velocity error pulse has amplitude that is clipped to a maximum saturation limit of 10 volts at the three phase current generator input, which limits the size of the torque loop input stimulus, in the derivation of the BLMD current command signals.

Examination of the family of characteristics pertaining to velocity control effort over a range of motor shaft inertial loads in Figure 57 indicate peak saturation over long pulse intervals ^{*} proportional to the inertial masses as in Figure 58 to be accelerated to the required speed V_{r}. This velocity error saturation is absent in the characteristics displayed in Figure 48 for 2volt ASD operation and results in linear PWM operation with a BLMD acceleration torque delivery commensurate with the velocity effort.

The variation of the simulated ASD velocity step response overshoot and settling time with shaft inertial load for a 4volt velocity command corresponding to 75% of rated shaft speed in Figures 50 and 49 respectively appear to be lower than those for 2 volt ASD operation at 36% rated speed. This is due to the saturation effect of the velocity error in ASD torque generation, which limits the peak amplitude swing of the oscillatory step response on velocity overshoot in Figure 59. However the rise time is longer in this instance as indicated through ASD simulation in Figure 51 for higher velocity command input (~4V) and increases almost linearly with shaft inertial load. The bigger the load inertia being handled during normal ASD operation the greater the corrective action required in the control effort to limit shaft velocity overshoot with inertial deceleration and settling time with reached velocity experienced with a non optimal speed controller during large step changes in velocity command input.

Unless the inertial load is known apriori in this scenario the variable PI parameters cannot be optimally selected for fast risetime and minimum overshoot except through offline manual tuning approximation procedures during the installation and commissioning phase (Moog GmbH, 1988) of the embedded drive in industrial applications. Inertial load parameter extraction in new ASD industrial applications, using the BLMD reference model with step response testing in an offline identification strategy for autotuning purposes, is difficult without knowledge of the initial PI settings of the velocity controller for zero load conditions before drive hookup to the embedded application. In this instance the procedure of accurate inertia parameter extraction using the ASD model in any identification strategy is complicated by the control action of outer velocity loop closure on the wideband inner torque loop when the variable PI settings are unknown. The ASD parameter identification problem is this case has to incorporate evaluation of the existing PI term settings in addition to the inertial load parameter in the new application in order to optimally design the embedded drive velocity controller. However if the drive is configured in torque control mode, thereby eliminating the velocity controller during step response testing in the commissioning phase, inertial load parameters are more amenable to extraction using low torque command input stimuli (<1V) for linear BLMD behaviour with small perturbation of the drive dynamics. The identification problem in this case has to focus only on the extraction of inertial load parameter J_{Tot} and also on the friction coefficient B_{m} if required.

## 5. Conclusions

Further BLMD simulation for various load torque settings, based on model confidence, yield results that compare favourably with those obtained from theoretical considerations using catalogued data for the actual drive concerned. The introduction of stator impedance angle compensation at high shaft speed results in improved motor power factor and better BLMD steady state performance. This is verified theoretically and illustrated through model simulation. Detailed BLMD simulation, configuration as an ASD with velocity feedback, is provided at internal observation nodes and checked against measured data at low and high command speed settings for confirmation of model accuracy and validation purposes.

## Acknowledgments

The author wishes to acknowledge

Eolas – The Irish Science and Technology Agency – for research funding.

Moog Ireland Ltd for brushless motor drive equipment for research purposes.

## References

- 1.
Balabanian N. Bickart T. A. 1969 Electrical Network Theory J.Wiley & Sons - 2.
J.D.,Franklin G. F. Powell J. 1980 - 3.
Guinee R. A. 2003 - 4.
Guinee R. A. Lyden C. 1999 Accurate Modelling And Simulation Of A DC Brushless Motor Drive System For High Performance Industrial Applications IEEE ISCAS’99 - 5.
Moog GmbH; 1988 - 6.
Moog GmbH; 1989 Moog Brushless Technology User Manual:D31X-XXX Motors,T158 01 X Controllers,T157-001 Power Supply,, D-7030 Bblingen.