Microsoft Word - A_33_Hatos_R.doc HUNGARIAN JOURNAL OF INDUSTRIAL CHEMISTRY VESZPRÉM Vol. 39(1) pp. 157-161 (2011) PARAMETER SENSITIVITY ANALYSIS OF AN INDUCTION MOTOR P. HATOS, A. FODOR , A. MAGYAR University of Pannonia, Department of Electrical Engineering and Information Systems, Veszprém, HUNGARY E-mail: foa@almos.uni-pannon.hu A simple dynamic model of an induction motor is presented in this paper based on engineering principles that describe the mechanical phenomena together with the electrical model. The investigated state space model consists of nonlinear state equations and linear output equations. The model has been verified under the usual controlled operating conditions when the speed is controlled. The effect of load on the controlled induction motor has been analyzed by simulation. The sensitivity analysis of the induction motor and the bridge of the inverter have been applied to determine the model parameters to be estimated. Keywords: induction machine, dynamic state space model, parameter sensitivity analysis Introduction The induction motors are the most commonly used electrical rotating machines in several industrial applications including the automotive industry, too. In the modern adjustable speed induction motor drives inverters are used to drive the three-phase motor as variable frequency voltage or current sources. Whatever the size and the application area, these motors share the most important dynamic properties, and their dynamic models have a similar structure. Therefore the final aim of our study is to design a controller that can control the speed and the torque of the induction motor. Because of the specialties and great practical importance of the induction motor in industrial applications, their modelling for control purposes is well investigated in the literature. Besides of the basic textbooks (see e.g. [1-3]), there are several papers that describe the modelling and use the developed models for the design of different types of controllers: vector control [1] and [4], sensor less vector control [5] and direct torque control (DTC) [6]. The aim of this paper is to build a simple dynamical model of the induction motor together with the three- phase inverter and analyze the models sensitivity of its parameters. The result of this analysis will be the basis of a subsequent parameter estimation step. The state space model has been implemented in Matlab/Simulink environment which enables us to analyze the parametric sensitivity based on simulation experiments. The model of the induction motor In this section the state- space model for an induction motor is developed. Modelling assumptions For constructing the induction motor model the following assumptions are made: ● symmetrical three phase windings, ● the slotting effect and the copper losses are neglected, ● the permeability of the iron parts is assumed to be infinite with linear magnetic properties, ● flux density is radial in the air gap, ● the spatial distribution of fluxes and apertures wave are considered to be sinusoidal, ● the spatial distribution of the stator fluxes and apertures wave are considered to be sinusoidal. According to the above modeling conditions the mathematical description of the induction motor is developed through the space vector theory. If the voltage of the stator is presumed to be the input excitation of the machine, then the spatial distribution along the stator of the x phase voltage can be described by the complex vector Usx(t). We can determine the orientation of the voltage vector Us the direction of the respective phase axis and the voltage polarity. (2.1) (2.2) 158 (2.3) , (2.4) where a = ej120° In equation (2.1) 2/3 is a normalizing factor. The flux density distribution can be obtained by integrating the current density wave along the cylinder of the stator. The flux linkage wave as a system variable, because it contains detailed information about the winding geometry. The rotating flux density wave induces voltages in the individual stator windings. Thus stator voltage Us(t) can be represented in the overall distributed voltages in all phase windings: (2.5) (2.6) (2.7) (2.8) Considering the stator of the induction machine as the primer side of the transformer, then using the Kirchoff’s voltage law the following equation can be written: (2.9) Figure 1: The equivalent circuit of the induction motor As for the secondary side of the transformer, it can be deduced that the same relationship is true for the rotor side space vectors: (2.10) Equations (2.9) and (2.10) describe the electromagnetic interaction as the connection of first order dynamical subsystems. Since four complex variables (is(t), ir(t), Ψs(t), Ψr(t)) are presented in these two equations, (2.1) and (2.5) flux equations are needed to complete the relationship between them. (2.11) (2.12) where angle ρ(t) defines the position of the rotor compared to the axis of the stator, while and are the three-phase inductances and ls, lr are the inductances of a stator and a rotor phase winding, Lm = 3/2*lm is the mutual inductance between the stator and the rotor. By applying the following substitutions: (2.13) (2.14) then the following equations are obtained with the flux connections in the model: (2.15) (2.16) The mechanical energy Pmech(t) of the system can be defined as: (2.17) where the mechanical energy Wmech(t) in case of rotating motor can be given by: (2.18) On the other hand, there is another expression for the mechanical energy: (2.19) where is the input electric power, is the resistive power loss, and is the air gap power. Using the above equations it can be concluded that: (2.20) The transformer can be decomposed into d-axis and q-axis. Park’s transformation converts the equations to a simplified and more tractable form. Figure 2: The equivalent circuit of the d axis of the induction motor Figure 3: The equivalent circuit of the q axis of the induction motor The actual terminal voltage v of the windings can be written in the following form 159 (2.21) where ij are the currents, Rj are the winding resistances, and Ψj are the flux linkages. Assume, that the positive directions of the stator currents point out of the induction motor terminals. By considering the d-axis and the q-axis of the induction motor, the following equations can be written: (2.22) (2.23) (2.24) (2.25) , (2.26) where ω is the reference frame angular velocity and ωr is the electrical angular velocity. (2.27) (2.28) (2.28) (2.29) The above model can be summarized in a state- space model by expressing the fluxes from the voltage equations. Parameter sensitivity analysis Thirteen parameters of the state space model of the induction motor and the bridge have been selected for sensitivity analysis (collected in Table 1), and the sensitivity of the state variables: voltage, phase a current, speed, electric torque, and outputs has been investigated for all of them by means of Matlab/Simulink dynamical simulation. Some simulation results are shown in Figs 4-7. The blue signal represents the simulation result with the nominal parameter values and the red signal represents the simulation result with the modified parameter values. Fig. 4 shows the model responses for changing a critically sensitive parameter (stator self inductance LS). It is apparent, that the speed diverges even for a 10% change of the parameter value. It can be seen that the speed of the motor becomes minus infinity and the electronic torque is zero. The case of a sensitive parameter (stator resistance RS) can be seen on Fig. 6. Table 1: The parameters of the induction motor and the bridge Parameter Initial value Dimension Name of the parameter Rs 0.435 Ohm Stator resistance Lls 0.002 H Stator leakage inductance Rr 0.816 Ohm Rotor resistance Llr 0.002 H Rotor leakage inductance M 0.0693 H Mutual inductance P 2 - Number of pole pairs In 0.089 kg·m2 Inertia of the motor Ed 1000 V Voltage of the inverter Rsn 10 5 Ohm Resistance of the snubber circuit Csn 10 10 F Capacitor of the snubber circuit Rbr 10 -3 Ohm Resistance of the bridge Ls M+Lls H Stator self inductance Lr M+Llr H Rotor self inductance 160 Figure 4: The -10% changing of parameter Ls Figure 5: The +50% changing of parameter M Figure 6: The -90% changing of parameter Rs Figure 7: The -90% changing of the resistance of the snubber circuit 161 As a result, the model parameters have been partitioned to four groups: ● Critically sensitive: The self inductance of the rotor widings (Lr) and the self inductance of the stator widings (Ls) ● Sensitive: The inertia (In) and the resistance of the stator (Rs) ● Less sensitive: The stator leakage inductance (Lls) and the resistance of the rotor (Rr), the rotor leakage inductance (Llr) and the Mutual inductance (M) ● Not sensitive: The Resistance of the bridge (RBr ), the capacitor of the snubber circuit (CSn) and the resistance of the snubber circuit (Rsn). Conclusions and future works Based on the results presented here, it is possible to select the candidate parameters for model parameter estimation based on real data that is a further aim of the authors, the four parameters are rotor self inductance (Lr), stator self inductance (Ls), Inertia and rotor resistance (Rs). The final aim of is to develop a simple yet detailed state space model of the induction motor for control purposes which gives us the possibility to develop and analyze different control strategies for the induction motor. ACKNOWLEDGEMENT We acknowledge the financial support of this work for the Hungarian State and the European Union under the TAMOP-4.2.1/B-09/1/KONV-2010-0003 project. REFERENCES 1. P. VAS: Artifical-intelligence-Based Electrical Machines and Drives, Oxford University Press, (1999) 2. P. VAS: Sensorless Vector and Direct Torque Control, Oxford University Press, (1998) 3. L. ZHENG, J. E. FLETCHER, B. W. WILLIAMS, X. HE: Dual-Plane Vector Control of a Five-Phase Induction Machine for an Improved Flux Pattern, IEEE Transaction on Industrial Electronics, 55(5), (2008), 1996–2005 4. E. LEVI: Impact of Iron Loss on Behavior of Vector Controlled Induction Machines, IEEE Transaction on Industry Applications, 31(6), (1995), 1287–1296 5. M. HASEGAWA, K. MATSUI: Robust Adaptive Full- Order Observer Design with Novel Adaptive Scheme for Speed Sensorless Vector Controlled Induction Motors, IEEE-IECON, (2002) 6. T. GEYER, G. PAPAFOTIU, M. MORARI: Model Predictive Direct Torque Control—Part I: Concept, Algorithm and Analysis, IEEE Transaction on Industrial Electronics, 56(6), (2009), 1894–1905 7. H. M. EMARA, W. ELSHAMY, A. BAHGAT: Parameter Identification of Induction Motor Using Modified Particle Swarm Optimization Algorithm, IEEE International Symposium on Industrial Electronics 2008, 841–847 << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice