Electromagnetic Modeling of the Propagation Characteristics of Satellite Communications Through Composite Precipitation Layers Science and Technology, 7 (2002) 109-121 © 2002 Sultan Qaboos University MATLAB and PSPICE Dynamic Model of Axial- Field Permanent-Magnet Motor A.H. Al-Badi and A. Gastli Department of Electrical Engineering, College of Engineering, Sultan Qaboos University, P.O.Box 33, Al Khod 123, Muscat, Sultanate of Oman, Email:albadi@squ.edu.om, Gastli@squ.edu.om. MATLABو PSPICE النموذج الديناميكي للمحرك الكهربائي ذي المجال المغناطيسي المحوري باستخدام عبداهللا البادي و عادل القصطلي يمكن استخدام هذا . تشرح هذه الورقة النموذج الديناميكي للمحرك الكهربائي ذي المجال المغناطيسي المحوري : خالصة تم مقارنة النموذج الديناميكي بالنتائج العملية . العابرة عند عمل هذه اآللة كمحرك أو كمولد كهربائي النموذج لدراسة االستجابة .ووجد تطابق بين النتائج مما يثبت صحة النموذج الديناميكي المستخدم ABSTRACT: In this paper a dynamic model of the axial-field, permanent-magnet machine will be discussed. The model may be used to study the transient response with the machine operating as a motor or as a generator. Experimental results are compared with Matlab/Simulink and Pspice predictions for a number of conditions, including the complete starting transient. Excellent correlation is obtained which validated both modeling approaches. KEYWORDS: Torus Motor, Dynamic Model, Pspice and Matlab / Simulink. 1. Introduction O ne of the important types of axial-field permanent-magnet machines is the Torus machine. A Torus machine is a double-sided, axial-field, disc-type, permanent-magnet, brushless d.c. motor. It is a compact lightweight machine, which offers significant advantages in application where space is restricted, due to its extremely short axial length. It has a simple toroidal, strip- wound stator core that carries a slotless toroidal winding. The rotor comprises two discs carrying six axially polarised magnets. The disc rotors act naturally as fans, so even with high electric loading good cooling of the stator winding is obtained. High magnetic loading may be achieved using sintered Nd-Fe-B for the magnets. The Torus has very small stator leakage inductances and relatively small mutual inductances, because of the large effective air-gap compared with a conventional machine with slotted stator. The axial length of this machine can be made small, for a given power rating, because the active conductor length is the radial length of the stator winding, which falls between the rotor magnets. The axially directed end-winding lengths are relatively short, yielding low resistance compared with the end-winding resistance of a radial-field machine. Because of the geometry of the machine, the self-inductances and the mutual inductances between the stator phases are constant (independent of the rotor angular position). Some potential applications of Torus machines include auxiliary power unit, wind-power generator, electric vehicle, high speed generator driven by a gas turbine, adjustable-speed pump drive, lawnmower motor, etc. Due to the complexity of the machine-drive systems and high switching frequency in the power electronics, the mathematical analysis and the experimental development are becoming difficult to perform as well costly and time consuming. The use of simulation is therefore strongly 109 AL-BADI and GASTLI recommended so that the performance can be predicted and optimised before physically building the system. Pspice and Matlab are widely used by researchers and designers for the design and analysis of electric circuits and systems. They have been intensively used in the simulation of the performance of electrical machines and power electronic circuits (Al-Badi and Green, 1998, Chalmers et al 1997, Gastli et al 1998, Gastli and Akherraz, 1998, Green and Al-Badi, 1997). 2. Mathematical Model 2.1 Circuit Equation Model The Torus motor has permanent magnets on the rotor and three stator windings. The equations of the motor can be expressed as follows:           +                     +                     =           c b a c b a cbcac cbbab cabaa c b a c b a e e e I I I LMM MLM MML p I I I R R R V V V ,, ,, ,, 00 00 00 (1) where, R is the stator winding resistance (0.133 ). Ω Va, Vb, Vc are the phase voltages. Ia, Ib, Ic are the phase winding currents. ea, eb, ec are the back emfs. La, Lb, Lc are the winding self inductances (200 µ H). Ma,b, Ma,c, Mb,c are the winding mutual inductances between the three phases (73 µ H). p is the mathematical differential operator d /dt. 2.2 Electromagnetic Torque Model The electromagnetic torque Te and the equations describing the mechanical response are: mccbbaae IeIeIeT ω/)( ++= (2) )( mlem BTTJp ωω −−= (3) mmp ωθ = mr Pp ωθ = where, Te is the electromagnetic torque. Tl is the load torque. J is the moment of inertia of the rotating mass (0.065 Kg m2). P is the pole pair number. B is the damping constant [0.00467 Nm/(rad/s]. p is the differential operator. mω is the rotor speed in mechanical rad/s. mθ is the angle defining the rotor position in mechanical radians. rθ is the angle defining the rotor position in electrical radians. 3. The Back EMF Model The back emf waveform is trapezoidal in shape, hence the transformation to d, q axis reference frame can not be easily accomplished. Therefore, the a, b, c phase model will be used. 110 MATLAB AND PSPICE DYNAMIC MODEL For the test machine with high-field magnets and air-gap winding, the effective air-gap length is constant. The measured back emf for a particular speed (104.72 rad/s) is represented in Figure 1. A Fourier analysis was performed on the waveform and values of the Fourier coefficients up to 37th harmonic are considered (Al-Badi and Green, 1998). 0 10 20 30 40 50 Time, [ms] V ol ta ge , [ V ] 20 10 0 - 10 - 20 Figure 1. The back emf. The initial model included coefficients up to 71st harmonic but this was found to have a negligible effect on the shape of the waveform. The Fourier series for the back emf of phase a is given as: ( ) ( ) ([ ]37373311 1 37sin3sinsin φ−ω+φ−ω+φ−ω ω ω = teteteea ) (4) Similar expressions can be obtained for the other two phases after including the appropriate phase shifts (-120o for phase b and -240o for phase c). The Fourier coefficients e1, e3,…,e37 are known for the particular speed 1ω . Fourier coefficients can be calculated at any speed using the following back emf-speed relationship: ω ω'' = i i e e (5) 4. Pspice Model for the Motor Drive System Figures 2, 3 and 4 show the 3-phase Torus motor system model, torque model, speed and position models respectively, that were developed using Pspice. The description of the full model was reported in (Al-Badi and Green, 1998). The resulting model is used to perform a comprehensive analysis of the system including the determination of the motor phase currents, line voltages, speed and electromagnetic torque. 5. Matlab/Simulink Model A computer model of the 3-phase Torus motor system was also developed using the Matlab/Simulink software. Figure 5 shows the whole drive system with the inverter and the motor. The description of the full model is reported hereafter. The resulting model was used to perform a detailed analysis of the system including the determination of the motor phase currents and line voltages. 111 AL-BADI and GASTLI ea eb ec R R R L a L b Lc M a,b dib/dt M a,c dic/dt M b,cdic/dt M b,adia/dt M c,bdib/dt M c,adia/dt V s Control circuit Rotor position sensor S1 S3 S5 S 4 S6 S 2 a b c Ia Figure 2. Circuit diagram for the drive system. mcIce ω/ mbIbe ω/ Tl R=106 maIae ω/ mBω Figure 3. Torque circuit, 1Nm/V. Jd m dtω / R=109 C=J Farad L=1 H P mω R=10-9 Figure 4. Speed circuit ( mω ), 1 rad/s/V and position circuit (θr ), electrical. 112 MATLAB AND PSPICE DYNAMIC MODEL +- v Phase A Phase B Phase C Wm Te Torus Machine Torque va wm t Speed Pulses +Vdc A B C -Vdc MOS_Inverter wm Pulses MOSFET Triggering Pulse Controller Te wm Load Plot Results Initialize Simulation Clock -VdcVdc 3_phase Diode Bridge Rectifier Figure 5. Torus motor drive system implemented with Simulink 5.1 Machine Model The Torus machine model described by (1) can be easily implemented with Simulink blocks after separating the machine parameters matrix and the back emf model. 5.1.1 Machine Parameters Figure 6 shows the model of the Torus machine. The machine windings consist of three-phase R-L circuits with self-inductances and mutual inductances. This three-phase circuit can be implemented with the Mutual Inductance block build in the Power System Blocksets library. The Mutual Inductance block implements a magnetic coupling between three separate windings. The self-resistance and inductance of each winding are specified on the first three entries of the dialog box and the mutual resistance and inductance in the last entry. 5.1.2 Back emf The back emf was modeled for the three phases using (4) and including the appropriate phase shift between phases. The implementation is done with Simulink Function block. Figure 7 shows the block diagram of the implemented model. Note that, for each phase, (4) was divided into four parts: part one includes all the harmonics from 1st to 9th, part two includes the sum of all harmonics from 10th to 19th, part three includes all harmonics from 20th to 29th and part 4 includes all harmonics from 30th to 37th. ωt is implemented using the integral block and integrating the electric angular velocity ω. 5.2 Inverter Model The inverter used was made of six power switches controlled with six synchronized pulse generators. The inverter configuration is shown in Figure 8. The dc voltage supply to the inverter is generated using the three-phase ac-dc converter (rectifier) shown in Figure 9. 5.2.1 Switch Model The power electronics switch was implemented using the MOSFET model, which is a block build in the Matlab Power System Bolckset. The MOSFET is a semiconductor device controllable by the gate signal (g > 0) if its drain current Id is positive (Id > 0). The MOSFET is connected in parallel with an internal diode that turns on when the MOSFET is reverse biased (Vds < 0). The model is simulated as a series combination of a variable resistor (Rt) and inductor (Lon) in series with a switch controlled by a logical signal (g > 0 or g = 0). The MOSFET turns on when the drain-source voltage is positive and a positive signal is applied at the gate input (i.e. g = 1). The 113 AL-BADI and GASTLI MOSFET block also contains a series Rs-Cs snubber circuit, which is connected in parallel with the MOSFET. 1 Te + - v + i - + i - + i - Vab f(u) Torque Te ic ib vab ia Mux Ia V + - V + - V + - wm e_c e_b e_a Back emf 4 W m 3 Phase C 2 Phase B 1 Phase A Figure 6. Torus machine model implemented with Simulink. 3 e_a 2 e_b 1 e_c f(u) ec3 f(u) ec2 f(u) ec1 f(u) ec f(u) eb3 f(u) eb2 f(u) eb1 f(u) eb f(u) ea3 f(u) ea2 f(u) ea1 f(u) ea ec eb ea Or Mux s 1 P 1 wm Figure 7. Back emf model implemented with Simulink. 114 MATLAB AND PSPICE DYNAMIC MODEL 4 -Vdc 3C 2B 1A m1 d g s m Mosfet6 d g s m Mosfet5 d g s m Mosfet4 d g s m Mosfet3 d g s m Mosfet2 d g s m Mosfet1 Demux 2 +Vdc 1 Pulses Figure 8. Inverter model implemented with Simulink 1 Vdc v 3 v 2 v 1 Vd(V) + - v Vd Vdc L3 L2 L1 a k m Diode6 a k m Diode5 a k m Diode4 a k m Diode3 a k m Diode2 a k m Diode1 C 1 -Vdc Figure 9. AC-DC converter circuit. 5.2.2 Gate Signal Model The gate pulses for the six transistors were generated by six-pulse generator blocks taken from the built in Simulink library. The pulses period, duty cycle, and phases can be easily adjusted and set to desired values. Figure 10 shows the block diagram of the model as it was implemented with Simulink blocks. 1 Pulses teta su3 sd3 teta su2 sd2 teta su1 sd1 Mux s 1 P 1 wm Figure 10. Gate signal model implemented with Simulink 115 AL-BADI and GASTLI 6. Models Validation The validity of both models was verified by applying them to the prediction of the performance of the Torus motor drive system and comparing the simulation results with the corresponding experimental data. 6.1 At Steady-State The dynamic model was first validated under steady state conditions by comparing simulated and measured phase current waveforms. Figures 11 and 12 show the simulated and measured phase current and line voltage waveforms at 40.5 V, 3.78 Nm. It indicated good correlation between the measured and simulated results. 6.2 At Starting A 30 V rectified d.c. supply is applied across the inverter and then the switches are controlled on and off. The motor phase current and speed during starting and subsequent acceleration up to a final steady speed are shown in Figures 13 and 14. Figure 13 shows the experimental and simulated results of the current waveform during starting. Figure 14 shows the experimental and simulated speed build up during starting. Closer examination of the current waveform shows that the peak current is determined almost by the circuit resistance because the inductance is very small, also evident in the finite rise in the phase current. The Figure shows also a spike superimposed on the top of each current pulse, owing to the switching which takes place every 60 degrees in other phases and to normal ripple in the rectified d.c. supply. The frequency of the first component is speed-dependent while the second is six times the a.c. supply frequency. Comparison between measured and simulated results of the phase currents and speed waveforms shows good correlation between these results. 7. Conclusion Dynamic performance prediction of an existing Torus motor drive system using Pspice and Matlab/Simulink models has been presented and compared with the actual drive performance. The two modeling approaches have been verified and validated with the experiments. It is possible to conclude that both approaches provide satisfactory tools for the investigation of the axial-field machine drive system performance. They can be helpful in the design and analysis of any axial-field machine drive system, which includes the machine, the load, the converter and the controller. 116 MATLAB AND PSPICE DYNAMIC MODEL 0 4 8 12 16 20 Time, [ms] 10 0 -10 C ur re nt , [ A ] (a) Measured 40 45 50 55 60 Time, [ms] - 10 0 10 C ur re nt , [ A ] (b) Simulated, Pspice 0 . 0 3 8 0 . 0 4 0 . 0 4 2 0 . 0 4 4 0 . 0 4 6 0 . 0 4 8 0 . 0 5 0 . 0 5 2 0 . 0 5 4 0 . 0 5 6 0 . 0 5 - 1 5 - 1 0 - 5 0 5 1 0 1 5 T i m e , [ s e c ] Li ne C ur re nt , [ A ] (C) Simulated, Matlab Figure 11. Phase current waveforms with Tl =3.78 Nm, 40.5 V, speed=1000 r.p.m. 117 AL-BADI and GASTLI 0 4 8 12 16 20 Time, [ms] - 50 - 25 0 25 50 V ol ta ge , [ V ] (a) Measured 40 45 50 55 60 Time, [ms] -50 -25 0 25 50 V ol ta ge , [ V ] (b) Simulated, PSpice 0 . 0 3 8 0 . 0 4 0 . 0 4 2 0 . 0 4 4 0 . 0 4 6 0 . 0 4 8 0 . 0 5 0 . 0 5 2 0 . 0 5 4 0 . 0 5 6 0 . 0 5 8 - 5 0 - 4 0 - 3 0 - 2 0 - 1 0 0 1 0 2 0 3 0 4 0 5 0 T i m e , [ s e c ] L in e V ol ta ge , [ V ] (c) Simulated, Matlab Figure 12. Line voltage waveforms with Tl =3.78 Nm, 40.5 V, speed=1000 r.p.m. 118 MATLAB AND PSPICE DYNAMIC MODEL 0 1 2 3 T ime, s -10 0 10 20 -20 Current, A (a) Measured 0 1 2 3 -2 5 -1 5 -5 5 1 5 2 5 T im e , s C u rre n t, A (b) Simulated, PSpice 0 0.5 1 1.5 2 2.5 3 -30 -20 -10 0 10 20 30 Tim e, s C urrent, A (c) Simulated, Matlab Figure 13. Phase current transient. 119 AL-BADI and GASTLI 0 1 2 3 T i m e , s 4 0 0 2 0 0 0 6 0 0 S p e e d , r e v / m i n (a) Measured 0 1 2 3 0 2 0 0 4 0 0 6 0 0 T i m e , s S p e e d , r e v / m i n (b) Simulated, Pspice 0 1 2 3 0 2 0 0 4 0 0 6 0 0 T i m e , s S p e e d , r e v / m i n (c) Simulated, Matlab Figure 14. Speed transient. 120 MATLAB AND PSPICE DYNAMIC MODEL References AL-BADI, A.H and GREEN, A.M. 1998. A transient model of the Torus electrical machine using PSPICE. 33rd Universities Power Engineering Conference, September 1998, Edinburgh, pp.254-257. CHALMERS, B.J., GREEN, A.M., REECE, A.B. and AL-BADI, A.H. 1997. Modelling and simulation of the Torus generator. IEE Proc. B, Vol.144, 6: 446-452. GASTLI, A., AKHERRAZ, M. and GAMMAL, M. 1998. Matlab/ Simulink/ANN Based Modeling and Simulation of A Stand-Alone Self-Excited Induction Generator. International Conference on Communication, Computer and Power, ICCCP’98, December 7-10 1998, Muscat, Sultanate of Oman, Sponsored by the IEEE and IEE, pp. 93-98. GASTLI, A., AKHERRAZ, M. 1998. Matlab/Simulink Simulation of Variable-Speed Busbar- Connected Self-Excited Induction Generator. International Conference on Communication, Computer and Power, ICCCP’98, December 7-10 1998, Muscat, Sultanate of Oman, Sponsored by the IEEE and IEE, pp. 258-263. GREEN, A.M. and AL-BADI, A.H. 1997. Modelling and simulation of a Torus generator using PSPICE. 32nd Universities Power Engineering Conference, September 1997, Manchester, pp.806-809. Received 24 June 2001 Accepted 13 November 2001 121 MATLAB and PSPICE Dynamic Model of Axial-Field Permanent-Magnet Motor A.H. Al-Badi and A. Gastli Department of Electrical Engineering, College of Engineering, Sultan Qaboos University, P.O.Box 33, Al Khod 123, Muscat, Sultanate of Oman, Email:albadi@squ.edu.om, Gastli@squ.edu.om. Mathematical Model 2.1 Circuit Equation Model 2.2 Electromagnetic Torque Model The Back EMF Model Pspice Model for the Motor Drive System Matlab/Simulink Model 5.1 Machine Model 5.1.1 Machine Parameters 5.1.2 Back emf 5.2 Inverter Model 5.2.1 Switch Model 5.2.2 Gate Signal Model Models Validation 6.1 At Steady-State 6.2 At Starting Conclusion References