Spatial-Mode-Shape Identification using a Continuous Wavelet Transform

Size: px
Start display at page:

Download "Spatial-Mode-Shape Identification using a Continuous Wavelet Transform"

Transcription

1 Strojniški vetnik - Journal of Mechanical Engineering 55(009)5, StartPage-EndPage Paper received: UDC 53/533 Paper accepted: x Spatial-Mode-Shape Identification uing a Continuou Wavelet Tranform Martin Čenik - Janko Slavič - Miha Boltežar * Univerity of Ljubljana, Faculty of Mechanical Engineering, Ljubljana, Slovenia Thi paper preent an experimental modal analyi of a damped multi-degree-of-freedom mechanical ytem uing a continuou wavelet tranform. An approximation of the wavelet tranform of the impule repone function i deduced, which erve a a bai for the extraction of the natural frequencie, the damping ratio and the mode hape. Due to an approximation with a finite Taylor erie, a computational error in the identified ocillatory amplitude occur and i oberved for the imulated ytem repone. The preented approach of modal identification i applied to real mechanical ytem, uch a a teel beam and the horizontal tail of an ultralight aircraft. Uing the propoed meaurement methodology, it i poible to recontruct the patial mode hape of any dynamic linear ytem with an arbitrary geometry. 009 Journal of Mechanical Engineering. All right reerved. Keyword: modal parameter, continuou wavelet tranform, patial vibration 0 INTRODUCTION One of the key tep in the development of a numerical model i it validation, where the modal analyi of a numerical model and the experimental modal analyi (EMA) of the real ytem are compared. With the reult of the comparion it i poible to determine the adequatene of the numerical model. The mot common way to undertake a modal analyi i by meauring the frequency repone-function (FRF) [] and []. Thi i done by ignal tranformation from the time domain to the frequency domain uing a Fourier tranformation. With the ue of a fat Fourier tranform (FFT) algorithm, the calculation of the FRF i a imple and ueful method ued to perform an EMA [3] to [5]. The decribed method i uitable for analyzing tationary ignal, ince it average the ignal in time domain. For an analyi of non-tationary ignal it i convenient to ue tranformation in the time or time-frequency domain, for intance, leatquare complex exponential (LSCE) [], continuou wavelet tranform (CWT) [6], etc. The impule repone of a mechanical ytem i a typical non-tationary proce where the ocillation amplitude decreae exponentially with time. Generally, the modal identification, baed on non-tationary ignal, i uually done with time-domain method [], e.g. LSCE for impule repone function or Ibrahim time domain (ITD) for free decay repone. The idea of applying a CWT to a dynamic ytem repone for damping identification wa firt introduced by Stazewki [7]. While Stazewki ued the Morlet wavelet function, Slavič et al. [8, 9, 0] developed their idea uing the Gabor wavelet function and focuing it on the edge effect and the relatively hort ignal. The propertie of the Gabor wavelet have been tudied in detail by Simonovki and Boltežar [], who alo ued it for fault detection in DC electro-motor []. Le and Argoul [3] extended the work of precendent reearcher and identified all the modal parameter by applying a CWT to the theoretical impule repone of a patial 4-degreeof-freedom (DOF) model with the Morlet, Cauchy and harmonic wavelet function. They alo explored the problem of the edge-effect and time-frequency localization. Furthermore, Lardie and Gouttebroze [4] applied a CWT to analyze real ignal meaured on a tower excited by the wind and determined it natural frequencie, damping ratio and one-dimenional mode hape. They compared the reult, obtained with a CWT to thoe obtained by an autocorrelation method. Their work wa extended by Huang and Su [5], who applied a CWT to dicrete equation of motion in order to identify the modal parameter of a teel frame ubjected to an earthquake excitation. A a reult, they obtained the frame * Corr. Author' Addre: Univerity of Ljubljana, Faculty of Mechanical Engineering, Aškerčeva 6, Ljubljana, Slovenia, miha.boltezar@f.uni-lj.i

2 Strojniški vetnik - Journal of Mechanical Engineering 55(009)5, StartPage-EndPage mode hape by a determination of the planar diplacement in a horizontal plane. To continue their work, Čenik et al. [6] tried to extract the patial mode hape of a mechanical ytem with an arbitrary geometry. Thi work repreent a continuation and enhancement of that reearch. In thi paper, a procedure for the identification of the modal parameter i introduced, baed on a CWT of the impule repone function. Due to the low-order approximation, a computation error of the ocillatory amplitude occur when dealing with ytem with high damping, which i preented in a numerical example. The method i alo teted on a real ytem of a teel beam The method i alo teted on a real ytem of a teel beam where the etimated mode hape are compared with the theoretical one. With the identification of the patial mode hape of the horizontal tail of an ultralight aircraft, the feaibility of the developed method i proven. CONTINUOUS WAVELET TRANSFORM The introductory part of CWT theory will be ummarized according to Mallat [6] and Tchamitchian et al. [7]. The continuou wavelet tranformation of the time function f(t), that atifie f(t) L (R) i defined a + * ψ u, () Wf ( u, ) = f ( t) ( t)d t, where ψ(t) repreent the wavelet function, the upercript * denote the complex conjugate, u i the tranlation parameter related to time and i the cale parameter that erve a an invere of the frequency. Fig. how the influence of the parameter u and on the wavelet function. The function f(t) L (R)may be treated a a mother wavelet function when atifying ocillatory, energy preervation and admiibility condition [6]. A modified wavelet function ψ u, (t), tranlated in the time domain and dilated in the frequency domain alo ha to atify the energy-preervation condition, and o we obtain t u ψu, () t = ψ, () ˆ i ω u ψ, ( ) ˆ u ω = e ψ ( ω ), (3) where ψˆu, ( ω ) denote the integral Fourier tranformation of the wavelet function ψ u, (t) in the frequency domain. The impule repone of a multi-degreeof-freedom (MDOF) ytem i expected to be a ummation of N mutually independent function pair [7] N i i () t () = A () t e ϕ, (4) f t i= i where A i (t) denote the amplitude and ϕ i (t) the phae (9). Therefore, the linearity property of the CWT [6] can be ued N N W αi fi ( u, ) = αi( Wfi)( u, ). (5) i= i=. Gabor Wavelet Function The Gabor mother wavelet function i defined a [6] t i (6) ψgabor ( t, σ, η) = e σ e η t. 4 πσ ( ) By introducing the tranlation and dilation parameter into (6), we obtain a family of Gabor wavelet function, defined in the time and frequency domain a [] ( t u) t u iη σ ψ, (,, u t σ η) = e e, 4, ( πσ ) ( ) ( ) η ω σ 4 i (7) ω u ψˆ u ω, σ, η = 4 π σ e e. (8) Fig.. ( ) mother ψ Gabor ( tu, = 0, = ) and ( ) modified ψ ( = = ) Gabor, u, tu, 4, 0.5 Čenik, M. - Slavič, J - Boltežar, M.

3 Strojniški vetnik - Journal of Mechanical Engineering 55(009)5, Startpage-EndPage. Approximation of the CWT of the Aymptotic Function The impule repone of a vicouly damped ingle-degree-of-freedom (SDOF) ytem can be written a [8] ζω0 t f () t = Ae co ( ω0d t+ ϕ0 ), (9) where ζ denote the damping ratio, ϕ 0 i the phae hift, ω 0 i the undamped angular frequency and ω 0d = ω 0 ζ i the damped angular frequency. In order to eparate the amplitude and the phae information of the repone function it can be written in a more general form a an analytic function [7] i f () t fa() t = Af () t e ϕ. (0) Similarly, a wavelet function may alo be given in the form of an analytic function a i () t a () t A () t e ϕ ψ = ψ ψ. With regard to the relation Wf ( u, ) = Wf ( u, ) [6] we can deduce that a + t u Wf ( u, ) = Af () t Aψ e i ϕf () t t u ϕψ d. t () In order to olve the equation (), a few approximation have to be made: the function S f (t) i approximated with a zero-order Taylor erie (i.e., a contant value A), and the function ϕ f (t) i approximated with a firt-order Taylor erie. By uing thee approximation, a general olution of the CWT can be deduced [8] ( ) ( f ( ) f ( )) ( f '', f ' ) i ϕ f ( u) Wf ( u, ) = Af ( u) e ψˆ ϕf '( u), η, σ + Er A '' u, ϕ ' u, where Er A ( u) ( u) () ϕ denote the error function due to the ue of the finite Taylor erie. On the bai of the linearity property of the CWT (5) we can apply an approximation () to multicomponent function, uch a the impule repone of the MDOF ytem (4) in order to identify the modal parameter of a particular natural frequency..3 Modal Parameter Identification Due to the nature of a CWT it i poible to tranform the meaured ignal in form of a time erie to the frequency domain without the lo of time-domain information. Furthermore, a CWT offer very good reitance to the noie in the meaured ignal. Regarding the nature of the CWT, it i poible to oberve the development of natural frequencie with time. In order to apply the CWT to modal analyi it i eential to analytically decribe the ditinctive propertie of the impule repone wavelet tranformation. The deduction of the impule-repone CWT will be preented according to Stazewki [7] and Slavič [8], who analyzed the Morlet and Gabor wavelet, repectively. From equation () we can ee that the eential contribution to the magnitude i provided by the tationary point of the argument of the integrand [6], i.e., the point t, uch that t u ϕf '( t) = ϕ ψ '. (3) The canonical pair ( u, ), which atify the condition t (u,)=u define the ridge of the wavelet tranform [7]. The ridge can be detected with different approache, and Stazewki preented three: the cro-ection method, the amplitude method and the phae method [7]. In thi work, the phae method, which i baed on a frequency match (3), i ued. In order to determine the ridge an initial etimation of the intantaneou frequency ω ha to be made. It i then entered into the recurive equation φ ( ui+, ηωi) φ( ui, ηωi) ωi+ =, (4) ui+ ui which gradually converge to an intantaneou frequency in the ignal (Fig.). The value of the wavelet tranform retricted to it ridge are called the keleton of the wavelet tranform; they erve for the modalparameter identification (6, 7). When conidering an impule repone function (9) an aumption of contant frequency can be made. Spatial-Mode-Shape Identification uing a Continuou Wavelet Tranform 3

4 Strojniški vetnik - Journal of Mechanical Engineering 55(009)5, StartPage-EndPage Fig.. Ridge detection uing the phae method with the recurive equation (4) Due to the linear relation between frequency and cale for the Gabor wavelet φ'(u,)=η/ we can conclude that the cale variable on the ridge (u)i contant (u)= 0. By introducing the Gabor wavelet function into equation (), regarding the frequency match and neglecting the error function, it follow that [6] 4 ζω0d u Wf ( u, 0) = Ae ( 4 πσ 0). (5) Finally, we can deduce equation for the direct identification of the damping ratio ζ and the initial amplitude A ( ln Wf ( u, 0 )) ζ = u cont., (6) ω0d ζω0d u e Wf ( u, 0 ) A = cont. (7) 4 4πσ η ω ( ( 0d ) ) are coniderably fater than the amplitude variation A A f f '( t) () t () t << ϕ '. (8) f When conidering ignal by which the aymptotic condition i not entirely atified, e.g., a highly damped impule repone, an error occur during the amplitude A identification. The preented error repreent a deviation of the recontructed amplitude A recontr from the real amplitude in the ytem repone A real. It can be oberved from a numerical example of a highly damped impule repone and depend on the damping ratio ζ and the parameter σ (Fig.3 and Fig.4); however, it i independent of the real amplitude value A real. Conequently, the ratio between the real amplitude A real and the recontructed amplitude A recontr which erve for the mode-hape recontruction, i contant A real A real = = cont. (9) A recontr A recontr Thi finding i important when a recontruction of the mode hape i made only by comparing the amplitude of the diplacement at different meaurement poition (Section.). Fig. 3. Recontructed initial amplitude at variable parameter σ and A real (t=0)=..4 A( t) Variation Influence on the Identified Amplitude An approximation for the theoretical CWT of the impule repone i deduced by approximating an impule repone (9) a a finite Taylor erie. An approximation of A(t) = cont. i exact only for aymptotic ignal, like the undamped repone, where the phae variation Fig. 4. Recontructed initial amplitude at Areal ( t = 0) = and A real ( t = 0) =. 4 Čenik, M. - Slavič, J - Boltežar, M.

5 Strojniški vetnik - Journal of Mechanical Engineering 55(009)5, Startpage-EndPage.5 Spatial Mode Shape Geometry In order to recontruct the patial mode hape the patial vibration have to be meaured with accelerometer at pre-defined meaurement poition. When performing a modal analyi on imple geometry uch a a teel beam (Section.), the repone accelerometer can be oriented in the ame direction and no iue arie from the geometry point of view. When analyzing a ytem with three-dimenional geometry, the orientation of the repone accelerometer in the ame direction a the reference one i not alway poible (Fig. 5). To overcome thi problem, the local coordinate ytem of each accelerometer i referenced to the global coordinate ytem of the meaured ample. One poibility of decribing thi linkage i with the Euler angle (Eq. 0 to ). Fig. 5. Problem of patial vibration meaurement To recontruct a ingle patial mode hape, the identified diplacement need to be tranformed from the local coordinate ytem (a meaured) into the global coordinate ytem with the ue of the rotational matrice 0 0 Rx, α = 0 co( α ) in( α ), (0) 0 in( α) co( α) ( ) ( ) co β 0 in β Ry, β = 0 0, () in( β ) 0 co( β ) ( ) ( ) ( γ) ( γ) co γ in γ 0 Rz, γ = in co 0, () 0 0 where α, β and γ denote the Euler rotation angle about the x, y or z axi, repectively..6 Identification of the Spatial Mode Shape The procedure to perform the experiment and recontruct the ingle patial mode hape i preented in the following tep: a) Experiment. Definition of the excitation poition.. Definition of the global coordinate ytem, n repone poition with n local coordinate ytem on the meaured ample geometry. 3. Selection of the reference poition and the reference direction (thi hould be the poition/direction where high ocillation of the identified mode are expected). 4. Excitation of the ytem with the impule hammer, acceleration meaurement on the reference poition/direction and on the repone poition() in the choen direction(). b) Recontruction 5. Identification of the ridge (4) and of the initial amplitude A ref j and A rep i,j (7) for the reference and repone ignal, repectively; j denote the conecutive number of the impule excitation and i i the number of the repone poition (i=,,n). 6. Computation of the normalized repone amplitude A norm i with repect to the reference amplitude A ref j with the ratio A norm i = A rep i,j / A ref j. 7. Tranformation of the normalized diplacement A norm i from the local to the global coordinate ytem (Eq. 0 to ). 8. Graphical preentation of the patialmode hape. Spatial-Mode-Shape Identification uing a Continuou Wavelet Tranform 5

6 Strojniški vetnik - Journal of Mechanical Engineering 55(009)5, StartPage-EndPage EXPERIMENT For a relatively imple teel beam, the preented approach for modal identification i experimentally compared to the beam theoretical mode hape. Furthermore, to demontrate an experimental identification of the patial modal analyi, the horizontal tail of an ultralight aircraft i preented. identification of global mode hape only negligible error are expected at the validation experiment. Eleven meauring poition were et along the beam: poition wa defined a the "reference" poition, and the remaining poition were defined a the "repone" poition, a hown in Fig. 6.. Experimental Validation In order to perform a pecific numerical manipulation, cutom program and algorithm were developed by the author. Cutom-made program offer a higher flexibility and provide the author with an inight into the computation proce. To check the correctne of the developed program a imple experiment with well-defined theoretical mode hape wa carried out. With algorithm that are verified on imple model, it i poible to perform a reliable modal analyi on more complex patial tructure. A imple validation experiment for method verification wa performed on a free-free upported teel beam with dimenion of 5x40x000 mm, excited by an impule hammer at a pre-defined "excitation" poition. The experiment wa made with everal impule excitation; during each excitation the repone accelerometer wa at a different poition, a decribed in the following. For the ake of experimental implicity, only two accelerometer were ued for the meaurement. One accelerometer wa denoted a the "reference" enor and wa alway placed in the "reference" poition. The econd accelerometer, denoted a the "repone" accelerometer wa placed at a different "repone" poition during each impule excitation. When carrying out the decribed experimental procedure a precaution hould be taken; due to changing of poition of the repone accelerometer a ma matrix of the whole beam ytem i changing. A a reult ome error of the recontructed mode hape and belonging natural frequencie could occur, epecially when uing accelerometer with coniderable ma and when identifying high-order mode hape with local nature. With application of a lightweight accelerometer with a ma of g and at Fig. 6. Meauring-poition placement on the analyzed beam After the whole et of meaurement wa carried out, we could obtain the impule repone from any pre-defined poition along the beam. Fig. 7. The validation-experiment configuration et-up With the ue of a CWT on the meaured repone it i poible to identify the modal parameter according to the theoretical backgrounddecribed in Section.3. The equation (4) i ued to define the time development of an oberved natural frequency. Once the value of the natural frequency i known, equation (6) and (7) are ued to define the damping ratio ζ of the ytem and the initial amplitude A at the location of the meaurement. Baed on the ratio between the diplacement at the reference and elected poition (Section.6), the firt nine normalized mode hape were recontructed and compared to the theoretical mode hape by computing the Modal Aurance Criterion MAC (Φ X, Φ A ) [], hown in Figure 8 and tated in Table. Table how the aociated natural frequencie f i and the damping ratio ζ i. 6 Čenik, M. - Slavič, J - Boltežar, M.

7 Strojniški vetnik - Journal of Mechanical Engineering 55(009)5, Startpage-EndPage Table. MAC(Φ X, Φ A ) matrix Mode MAC (Φ X, Φ A ) Table. Beam natural frequencie and damping ratio i f i [Hz] ζ i Fig. 8. Modal aurance criterion matrix. Horizontal Tail Experiment Fig. 9. Horizontal tail in iometric view An experiment imilar to the one decribed in Section.. wa carried out on the horizontal tail of an ultralight aircraft. Since the tail geometry i decribed uing patial curve (Fig. 9), a patial vibration meaurement had to be performed. In order to reduce the experimental complexity and the number of nonlinearitie in the ytem, the tail wa upported free-free. Twenty-eight meauring poition were defined on the tail urface, a hown in Fig.0. Altogether, 8 local coordinate ytem were defined in uch a way that the z axi wa normal to the urface and the x axi followed obviou line in the tail geometry. A a reference, meaurement poition wa choen (Fig. 0 and ), oriented in the local direction z. Fig. 0. Meauring poition placement and the global coordinate ytem of the horizontal tail In each meaurement two ignal were ampled: - the acceleration at the reference poition in the reference direction, - the acceleration at the repone poition in one of the defined local axe. Altogether, 84 eparate meaurement were made (8 meauring poition in 3 direction). Spatial-Mode-Shape Identification uing a Continuou Wavelet Tranform 7

8 Strojniški vetnik - Journal of Mechanical Engineering 55(009)5, StartPage-EndPage b) Second mode hape Fig.. A horizontal tail with accelerometer placed at the reference poition (left) and at the i-th repone poition (right) In order to tranform the identified diplacement from the local to the global coordinate ytem, pecial oftware wa developed by the author, to which an arbitrary geometry can be imported. The geometrical characteritic are imported into the oftware in the form of a finite-element-method (FEM) meh that conit of volume element. The uer define the poition and the orientation of the local coordinate ytem and the correponding identified diplacement of the ingle mode hape. The oftware generate a cript which can be run in ANSYS and give a recontructed mode hape in the form of a tatic problem olution. The identified natural frequencie and the damping ratio are hown in Table, and the experimentally recontructed normalized mode hape are hown in Fig.. Table. Horizontal tail natural frequencie and damping ratio i f i [Hz] ζ i c) Third mode hape Fig.. Recontructed normalized horizontal tail mode hape 3 DISCUSSION From the validation experiment, which erved a an ideal, real ytem, it i obviou that the identified mode hape agree with theoretical mode hape to a high degree, although ome change in ma matrix occured due to the diplacement of a repone accelerometer. With the imple beam experiment the applicability of the continuou wavelet tranform to determine the modal parameter i confirmed. A meaurement approach with the introduction of a reference accelerometer wa hown to be uitable. According to the experiment realization and the reult, performed and obtained in Section., the preented approach provide a imple and feaible method for an experimental determination of the modal parameter. Although a mechanical ytem, uch a a horizontal tail, i not an ideal ytem without nonlinearitie, the extraction of it modal parameter with a CWT wa hown to be a robut and effective method. 4 CONCLUSIONS a) Firt mode hape In thi paper a modal parameter identification of a dynamic ytem uing a continuou wavelet tranformation of ytem impule repone wa preented. The computation error due to the finite Taylor erie wa indicated 8 Čenik, M. - Slavič, J - Boltežar, M.

9 Strojniški vetnik - Journal of Mechanical Engineering 55(009)5, Startpage-EndPage and it influence on the mode-hape recontruction wa hown. Steel beam mode hape were recontructed and compared to the theoretical mode hape. The patial mode hape of a free-free upported horizontal tail of an ultralight aircraft were recontructed. A comparion among the theoretical and experimental mode hape of the beam confirmed the uitability of the propoed method for modehape recontruction. The main contribution of thi paper i an enhancement of known method for EMA uing a CWT for one-dimenional cae with a new method for the recontruction of patial mode hape for an arbitrary geometry. In thi cae it i a horizontal tail. The experimental reult can be ued for the validation of the numerical model. 5 ACKNOWLEDGEMENT The author would like to thank Pipitrel d.o.o. for providing the horizontal tail of the ultralight aircraft. 6 REFERENCES [] Maia, N.M.M., Silva, J.M.M., Theoretical and experimental modal analyi, Reearch Studie Pre, Someret, 997. [] Ewin, D.J., Modal Teting: Theory and Practice, Reearch Studie Pre, Lechtworth, 986. [3] Čermelj, P., Boltežar, M., An indirect approach to invetigating the dynamic of a tructure containing ball bearing, Journal of Sound and Vibration, 004, vol. 76, no. /, p [4] Čermelj, P., Boltežar, M., Modelling localied nonlinearitie uing the harmonic nonlinear uper model, Journal of Sound and Vibration, 006, vol. 98, no. 4/5, p [5] Boltežar, M., Čermelj, P., Dynamic of complex tructure - valid modeling of indutrial product, In: Korelc, J., Zupan, D. (Ed.): Kuhljevi dnevi 007, Snovik, p. 7-4, 007. [6] Mallat, S., A wavelet tour of ignal proceing, Academic Pre, 999, econd edition. [7] Stazewki, W.J., Identification of damping in mdof ytem uing time-cale decompoition, Journal of Sound and Vibration, 997, vol. 03, p [8] Slavič, J., Simonovki, I., Boltežar, M., Damping identification uing continuou wavelet tranform: application to real data, Journal of Sound and Vibration, 003, vol. 6, p [9] Slavič, J., Boltežar, M., Enhanced identification of damping uing continuou wavelet tranform, Journal of Mechanical Engineering, 00, vol. 48, no., p [0] Boltežar, M., Slavič, J., Enhancement to the continuou wavelet tranform for damping identification on hort ignal, Mechanical ytem and ignal proceing, 004, vol. 8, no. 5, p [] Simonovki, I., Boltežar, M., The norm and variance of the gabor, morlet and general harmonic wavelet function, Journal of Sound and Vibration, 003, vol. 64, p [] Boltežar, M., Simonovki, I., Furlan, M., Fault detection in DC electro motor uing the continuou wavelet tranform, Meccanica, 003, vol. 38, no., p [3] Le, T.P., Argoul, P., Continuou wavelet tranform for modal identification uing free decay repone, Journal of Sound and Vibration, 004, vol. 77, p [4] Lardie, J., Gouttebroze, S., Identification of modal parameter uing free decay repone, International Journal of Mechanical Science, 00, vol. 44, p [5] Huang, C.S., Su, W.C., Identification of modal parameter of a time invariant linear ytem by continuou wavelet tranform, Mechanical Sytem and Signal Proceing, 006, vol., p [6] Čenik, M., Slavič, J., Boltežar, M., Modal parameter identification uing continuou wavelet tranform, In: Boltežar, M., Slavič, J. (Ed.): Kuhljevi dnevi 008, Cerklje na Gorenjkem, p. 4-48, 008. [7] Tchamitcian, P., Torreani, B., Aymptotic wavelet and gabor analyi: Extraction of intantaneou frequencie, IEEE tranaction on information technology, 99, vol. 38, p Spatial-Mode-Shape Identification uing a Continuou Wavelet Tranform 9

Determination of the local contrast of interference fringe patterns using continuous wavelet transform

Determination of the local contrast of interference fringe patterns using continuous wavelet transform Determination of the local contrat of interference fringe pattern uing continuou wavelet tranform Jong Kwang Hyok, Kim Chol Su Intitute of Optic, Department of Phyic, Kim Il Sung Univerity, Pyongyang,

More information

CHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL

CHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL 98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i

More information

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3

More information

Finding the location of switched capacitor banks in distribution systems based on wavelet transform

Finding the location of switched capacitor banks in distribution systems based on wavelet transform UPEC00 3t Aug - 3rd Sept 00 Finding the location of witched capacitor bank in ditribution ytem baed on wavelet tranform Bahram nohad Shahid Chamran Univerity in Ahvaz bahramnohad@yahoo.com Mehrdad keramatzadeh

More information

THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER

THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER Proceeding of IMAC XXXI Conference & Expoition on Structural Dynamic February -4 Garden Grove CA USA THE EXPERIMENTAL PERFORMANCE OF A NONLINEAR DYNAMIC VIBRATION ABSORBER Yung-Sheng Hu Neil S Ferguon

More information

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get Lecture 25 Introduction to Some Matlab c2d Code in Relation to Sampled Sytem here are many way to convert a continuou time function, { h( t) ; t [0, )} into a dicrete time function { h ( k) ; k {0,,, }}

More information

SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD

SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem

More information

Quantifying And Specifying The Dynamic Response Of Flowmeters

Quantifying And Specifying The Dynamic Response Of Flowmeters White Paper Quantifying And Specifying The Dynamic Repone Of Flowmeter DP Flow ABSTRACT The dynamic repone characteritic of flowmeter are often incompletely or incorrectly pecified. Thi i often the reult

More information

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281 72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition

More information

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,

More information

Modelling and Wavelet-Based Identification of 3-DOF Vehicle Suspension System

Modelling and Wavelet-Based Identification of 3-DOF Vehicle Suspension System Journal of Software Engineering and Application, 011, 4, 67-681 doi:10.436/jea.011.41079 Publihed Online December 011 (http://www.scirp.org/journal/jea) Modelling and Wavelet-Baed Identification of 3-DOF

More information

Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis

Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis Advanced Digital ignal Proceing Prof. Nizamettin AYDIN naydin@yildiz.edu.tr Time-Frequency Analyi http://www.yildiz.edu.tr/~naydin 2 tationary/nontationary ignal Time-Frequency Analyi Fourier Tranform

More information

The Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor

The Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor The Influence of the Load Condition upon the Radial Ditribution of Electromagnetic Vibration and Noie in a Three-Phae Squirrel-Cage Induction Motor Yuta Sato 1, Iao Hirotuka 1, Kazuo Tuboi 1, Maanori Nakamura

More information

USING NONLINEAR CONTROL ALGORITHMS TO IMPROVE THE QUALITY OF SHAKING TABLE TESTS

USING NONLINEAR CONTROL ALGORITHMS TO IMPROVE THE QUALITY OF SHAKING TABLE TESTS October 12-17, 28, Beijing, China USING NONLINEAR CONTR ALGORITHMS TO IMPROVE THE QUALITY OF SHAKING TABLE TESTS T.Y. Yang 1 and A. Schellenberg 2 1 Pot Doctoral Scholar, Dept. of Civil and Env. Eng.,

More information

A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: CORRESPONDENCE: ABSTRACT

A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: CORRESPONDENCE: ABSTRACT A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: Zenon Medina-Cetina International Centre for Geohazard / Norwegian Geotechnical Intitute Roger

More information

Evolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis

Evolutionary Algorithms Based Fixed Order Robust Controller Design and Robustness Performance Analysis Proceeding of 01 4th International Conference on Machine Learning and Computing IPCSIT vol. 5 (01) (01) IACSIT Pre, Singapore Evolutionary Algorithm Baed Fixed Order Robut Controller Deign and Robutne

More information

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,

More information

Supplementary Figures

Supplementary Figures Supplementary Figure Supplementary Figure S1: Extraction of the SOF. The tandard deviation of meaured V xy at aturated tate (between 2.4 ka/m and 12 ka/m), V 2 d Vxy( H, j, hm ) Vxy( H, j, hm ) 2. The

More information

Bogoliubov Transformation in Classical Mechanics

Bogoliubov Transformation in Classical Mechanics Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How

More information

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement

More information

Fluid-structure coupling analysis and simulation of viscosity effect. on Coriolis mass flowmeter

Fluid-structure coupling analysis and simulation of viscosity effect. on Coriolis mass flowmeter APCOM & ISCM 11-14 th December, 2013, Singapore luid-tructure coupling analyi and imulation of vicoity effect on Corioli ma flowmeter *Luo Rongmo, and Wu Jian National Metrology Centre, A*STAR, 1 Science

More information

Finite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization

Finite Element Analysis of a Fiber Bragg Grating Accelerometer for Performance Optimization Finite Element Analyi of a Fiber Bragg Grating Accelerometer for Performance Optimization N. Baumallick*, P. Biwa, K. Dagupta and S. Bandyopadhyay Fiber Optic Laboratory, Central Gla and Ceramic Reearch

More information

LTV System Modelling

LTV System Modelling Helinki Univerit of Technolog S-72.333 Potgraduate Coure in Radiocommunication Fall 2000 LTV Stem Modelling Heikki Lorentz Sonera Entrum O heikki.lorentz@onera.fi Januar 23 rd 200 Content. Introduction

More information

S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS

S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS S_LOOP: SINGLE-LOOP FEEDBACK CONTROL SYSTEM ANALYSIS by Michelle Gretzinger, Daniel Zyngier and Thoma Marlin INTRODUCTION One of the challenge to the engineer learning proce control i relating theoretical

More information

Simple Observer Based Synchronization of Lorenz System with Parametric Uncertainty

Simple Observer Based Synchronization of Lorenz System with Parametric Uncertainty IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih

More information

Mechanics. Free rotational oscillations. LD Physics Leaflets P Measuring with a hand-held stop-clock. Oscillations Torsion pendulum

Mechanics. Free rotational oscillations. LD Physics Leaflets P Measuring with a hand-held stop-clock. Oscillations Torsion pendulum Mechanic Ocillation Torion pendulum LD Phyic Leaflet P.5.. Free rotational ocillation Meauring with a hand-held top-clock Object of the experiment g Meauring the amplitude of rotational ocillation a function

More information

NONLINEAR CONTROLLER DESIGN FOR A SHELL AND TUBE HEAT EXCHANGER AN EXPERIMENTATION APPROACH

NONLINEAR CONTROLLER DESIGN FOR A SHELL AND TUBE HEAT EXCHANGER AN EXPERIMENTATION APPROACH International Journal of Electrical, Electronic and Data Communication, ISSN: 232-284 Volume-3, Iue-8, Aug.-25 NONLINEAR CONTROLLER DESIGN FOR A SHELL AND TUBE HEAT EXCHANGER AN EXPERIMENTATION APPROACH

More information

Lecture 10 Filtering: Applied Concepts

Lecture 10 Filtering: Applied Concepts Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering

More information

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou

More information

Hybrid Projective Dislocated Synchronization of Liu Chaotic System Based on Parameters Identification

Hybrid Projective Dislocated Synchronization of Liu Chaotic System Based on Parameters Identification www.ccenet.org/ma Modern Applied Science Vol. 6, No. ; February Hybrid Projective Dilocated Synchronization of Liu Chaotic Sytem Baed on Parameter Identification Yanfei Chen College of Science, Guilin

More information

Social Studies 201 Notes for March 18, 2005

Social Studies 201 Notes for March 18, 2005 1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

Estimating floor acceleration in nonlinear multi-story moment-resisting frames

Estimating floor acceleration in nonlinear multi-story moment-resisting frames Etimating floor acceleration in nonlinear multi-tory moment-reiting frame R. Karami Mohammadi Aitant Profeor, Civil Engineering Department, K.N.Tooi Univerity M. Mohammadi M.Sc. Student, Civil Engineering

More information

The Hassenpflug Matrix Tensor Notation

The Hassenpflug Matrix Tensor Notation The Haenpflug Matrix Tenor Notation D.N.J. El Dept of Mech Mechatron Eng Univ of Stellenboch, South Africa e-mail: dnjel@un.ac.za 2009/09/01 Abtract Thi i a ample document to illutrate the typeetting of

More information

An estimation approach for autotuning of event-based PI control systems

An estimation approach for autotuning of event-based PI control systems Acta de la XXXIX Jornada de Automática, Badajoz, 5-7 de Septiembre de 08 An etimation approach for autotuning of event-baed PI control ytem Joé Sánchez Moreno, María Guinaldo Loada, Sebatián Dormido Departamento

More information

POWER SYSTEM SMALL SIGNAL STABILITY ANALYSIS BASED ON TEST SIGNAL

POWER SYSTEM SMALL SIGNAL STABILITY ANALYSIS BASED ON TEST SIGNAL POWE YEM MALL INAL ABILIY ANALYI BAE ON E INAL Zheng Xu, Wei hao, Changchun Zhou Zheang Univerity, Hangzhou, 37 PChina Email: hvdc@ceezueducn Abtract - In thi paper, a method baed on ome tet ignal (et

More information

Digital Control System

Digital Control System Digital Control Sytem Summary # he -tranform play an important role in digital control and dicrete ignal proceing. he -tranform i defined a F () f(k) k () A. Example Conider the following equence: f(k)

More information

Sampling and the Discrete Fourier Transform

Sampling and the Discrete Fourier Transform Sampling and the Dicrete Fourier Tranform Sampling Method Sampling i mot commonly done with two device, the ample-and-hold (S/H) and the analog-to-digital-converter (ADC) The S/H acquire a CT ignal at

More information

A BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH

A BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Mathematical and Computational Application Vol. 11 No. pp. 181-191 006. Aociation for Scientific Reearch A BATCH-ARRIVA QEE WITH MTIPE SERVERS AND FZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Jau-Chuan

More information

Practical issues in signal processing for structural flexibility identification

Practical issues in signal processing for structural flexibility identification See dicuion, tat, and author profile for thi publication at: http://www.reearchgate.net/publication/276384859 Practical iue in ignal proceing for tructural flexibility identification Article in SMART STRUCTURES

More information

Analysis of Step Response, Impulse and Ramp Response in the Continuous Stirred Tank Reactor System

Analysis of Step Response, Impulse and Ramp Response in the Continuous Stirred Tank Reactor System ISSN: 454-50 Volume 0 - Iue 05 May 07 PP. 7-78 Analyi of Step Repone, Impule and Ramp Repone in the ontinuou Stirred Tank Reactor Sytem * Zohreh Khohraftar, Pirouz Derakhhi, (Department of hemitry, Science

More information

The continuous time random walk (CTRW) was introduced by Montroll and Weiss 1.

The continuous time random walk (CTRW) was introduced by Montroll and Weiss 1. 1 I. CONTINUOUS TIME RANDOM WALK The continuou time random walk (CTRW) wa introduced by Montroll and Wei 1. Unlike dicrete time random walk treated o far, in the CTRW the number of jump n made by the walker

More information

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified

More information

Digital Control System

Digital Control System Digital Control Sytem - A D D A Micro ADC DAC Proceor Correction Element Proce Clock Meaurement A: Analog D: Digital Continuou Controller and Digital Control Rt - c Plant yt Continuou Controller Digital

More information

Reliability Analysis of Embedded System with Different Modes of Failure Emphasizing Reboot Delay

Reliability Analysis of Embedded System with Different Modes of Failure Emphasizing Reboot Delay International Journal of Applied Science and Engineering 3., 4: 449-47 Reliability Analyi of Embedded Sytem with Different Mode of Failure Emphaizing Reboot Delay Deepak Kumar* and S. B. Singh Department

More information

Lecture Notes II. As the reactor is well-mixed, the outlet stream concentration and temperature are identical with those in the tank.

Lecture Notes II. As the reactor is well-mixed, the outlet stream concentration and temperature are identical with those in the tank. Lecture Note II Example 6 Continuou Stirred-Tank Reactor (CSTR) Chemical reactor together with ma tranfer procee contitute an important part of chemical technologie. From a control point of view, reactor

More information

III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES

III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SUBSTANCES III.9. THE HYSTERESIS CYCLE OF FERROELECTRIC SBSTANCES. Work purpoe The analyi of the behaviour of a ferroelectric ubtance placed in an eternal electric field; the dependence of the electrical polariation

More information

CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS

CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Copyright 22 IFAC 5th Triennial World Congre, Barcelona, Spain CHEAP CONTROL PERFORMANCE LIMITATIONS OF INPUT CONSTRAINED LINEAR SYSTEMS Tritan Pérez Graham C. Goodwin Maria M. Serón Department of Electrical

More information

Robust Decentralized Design of H -based Frequency Stabilizer of SMES

Robust Decentralized Design of H -based Frequency Stabilizer of SMES International Energy Journal: Vol. 6, No., Part, June 005-59 Robut Decentralized Deign of H -baed Frequency Stabilizer of SMES www.erd.ait.ac.th/reric C. Vorakulpipat *, M. Leelajindakrirerk *, and I.

More information

EE 4443/5329. LAB 3: Control of Industrial Systems. Simulation and Hardware Control (PID Design) The Inverted Pendulum. (ECP Systems-Model: 505)

EE 4443/5329. LAB 3: Control of Industrial Systems. Simulation and Hardware Control (PID Design) The Inverted Pendulum. (ECP Systems-Model: 505) EE 4443/5329 LAB 3: Control of Indutrial Sytem Simulation and Hardware Control (PID Deign) The Inverted Pendulum (ECP Sytem-Model: 505) Compiled by: Nitin Swamy Email: nwamy@lakehore.uta.edu Email: okuljaca@lakehore.uta.edu

More information

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall Beyond Significance Teting ( nd Edition), Rex B. Kline Suggeted Anwer To Exercie Chapter. The tatitic meaure variability among core at the cae level. In a normal ditribution, about 68% of the core fall

More information

Lateral vibration of footbridges under crowd-loading: Continuous crowd modeling approach

Lateral vibration of footbridges under crowd-loading: Continuous crowd modeling approach ateral vibration of footbridge under crowd-loading: Continuou crowd modeling approach Joanna Bodgi, a, Silvano Erlicher,b and Pierre Argoul,c Intitut NAVIER, ENPC, 6 et 8 av. B. Pacal, Cité Decarte, Champ

More information

Jan Purczyński, Kamila Bednarz-Okrzyńska Estimation of the shape parameter of GED distribution for a small sample size

Jan Purczyński, Kamila Bednarz-Okrzyńska Estimation of the shape parameter of GED distribution for a small sample size Jan Purczyńki, Kamila Bednarz-Okrzyńka Etimation of the hape parameter of GED ditribution for a mall ample ize Folia Oeconomica Stetinenia 4()/, 35-46 04 Folia Oeconomica Stetinenia DOI: 0.478/foli-04-003

More information

Seismic responses analysis of long continuous rigid-framed bridge subjected to multi-support excitations

Seismic responses analysis of long continuous rigid-framed bridge subjected to multi-support excitations Seimic repone analyi of long continuou rigid-framed ridge ujected to multi-upport excitation Jitao Li 1 and Qinghan Yang 2 1 Potgraduate, School of Civil Engineering, Beijing Jiaotong Univerity, Beijing.

More information

Codes Correcting Two Deletions

Codes Correcting Two Deletions 1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of

More information

Design By Emulation (Indirect Method)

Design By Emulation (Indirect Method) Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal

More information

Chapter 1 Basic Description of Laser Diode Dynamics by Spatially Averaged Rate Equations: Conditions of Validity

Chapter 1 Basic Description of Laser Diode Dynamics by Spatially Averaged Rate Equations: Conditions of Validity Chapter 1 Baic Decription of Laer Diode Dynamic by Spatially Averaged Rate Equation: Condition of Validity A laer diode i a device in which an electric current input i converted to an output of photon.

More information

Recent progress in fire-structure analysis

Recent progress in fire-structure analysis EJSE Special Iue: Selected Key Note paper from MDCMS 1 1t International Conference on Modern Deign, Contruction and Maintenance of Structure - Hanoi, Vietnam, December 2007 Recent progre in fire-tructure

More information

Molecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions

Molecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,

More information

RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS

RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS www.arpapre.com/volume/vol29iue1/ijrras_29_1_01.pdf RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS Sevcan Demir Atalay 1,* & Özge Elmataş Gültekin

More information

Social Studies 201 Notes for November 14, 2003

Social Studies 201 Notes for November 14, 2003 1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

APPLICATION OF THE SINGLE IMPACT MICROINDENTATION FOR NON- DESTRUCTIVE TESTING OF THE FRACTURE TOUGHNESS OF NONMETALLIC AND POLYMERIC MATERIALS

APPLICATION OF THE SINGLE IMPACT MICROINDENTATION FOR NON- DESTRUCTIVE TESTING OF THE FRACTURE TOUGHNESS OF NONMETALLIC AND POLYMERIC MATERIALS APPLICATION OF THE SINGLE IMPACT MICROINDENTATION FOR NON- DESTRUCTIVE TESTING OF THE FRACTURE TOUGHNESS OF NONMETALLIC AND POLYMERIC MATERIALS REN A. P. INSTITUTE OF APPLIED PHYSICS OF THE NATIONAL ACADEMY

More information

MODERN CONTROL SYSTEMS

MODERN CONTROL SYSTEMS MODERN CONTROL SYSTEMS Lecture 1 Root Locu Emam Fathy Department of Electrical and Control Engineering email: emfmz@aat.edu http://www.aat.edu/cv.php?dip_unit=346&er=68525 1 Introduction What i root locu?

More information

EE/ME/AE324: Dynamical Systems. Chapter 8: Transfer Function Analysis

EE/ME/AE324: Dynamical Systems. Chapter 8: Transfer Function Analysis EE/ME/AE34: Dynamical Sytem Chapter 8: Tranfer Function Analyi The Sytem Tranfer Function Conider the ytem decribed by the nth-order I/O eqn.: ( n) ( n 1) ( m) y + a y + + a y = b u + + bu n 1 0 m 0 Taking

More information

LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM. P.Dickinson, A.T.Shenton

LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM. P.Dickinson, A.T.Shenton LOW ORDER MIMO CONTROLLER DESIGN FOR AN ENGINE DISTURBANCE REJECTION PROBLEM P.Dickinon, A.T.Shenton Department of Engineering, The Univerity of Liverpool, Liverpool L69 3GH, UK Abtract: Thi paper compare

More information

Dynamic Simulation of a Three-Phase Induction Motor Using Matlab Simulink

Dynamic Simulation of a Three-Phase Induction Motor Using Matlab Simulink Dynamic Simulation of a ThreePhae Induction Motor Uing Matlab Simulink Adel Aktaibi & Daw Ghanim, graduate tudent member, IEEE, M. A. Rahman, life fellow, IEEE, Faculty of Engineering and Applied Science,

More information

SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48)

SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lectures 41-48) Chapter 5 SERIES COMPENSATION: VOLTAGE COMPENSATION USING DVR (Lecture 41-48) 5.1 Introduction Power ytem hould enure good quality of electric power upply, which mean voltage and current waveform hould

More information

CONTROL OF INTEGRATING PROCESS WITH DEAD TIME USING AUTO-TUNING APPROACH

CONTROL OF INTEGRATING PROCESS WITH DEAD TIME USING AUTO-TUNING APPROACH Brazilian Journal of Chemical Engineering ISSN 004-6632 Printed in Brazil www.abeq.org.br/bjche Vol. 26, No. 0, pp. 89-98, January - March, 2009 CONROL OF INEGRAING PROCESS WIH DEAD IME USING AUO-UNING

More information

SIMPLIFIED MODEL FOR EPICYCLIC GEAR INERTIAL CHARACTERISTICS

SIMPLIFIED MODEL FOR EPICYCLIC GEAR INERTIAL CHARACTERISTICS UNIVERSITY OF PITESTI SCIENTIFIC BULLETIN FACULTY OF ECHANICS AND TECHNOLOGY AUTOOTIVE erie, year XVII, no. ( 3 ) SIPLIFIED ODEL FOR EPICYCLIC GEAR INERTIAL CHARACTERISTICS Ciobotaru, Ticuşor *, Feraru,

More information

Application of Gradient Projection for Sparse Reconstruction to Compressed Sensing for Image Reconstruction of Electrical Capacitance Tomography

Application of Gradient Projection for Sparse Reconstruction to Compressed Sensing for Image Reconstruction of Electrical Capacitance Tomography Journal of Electrical and Electronic Engineering 08; 6(): 46-5 http://www.ciencepublihinggroup.com/j/jeee doi: 0.648/j.jeee.08060. ISS: 39-63 (Print); ISS: 39-605 (Online) Application of Gradient Projection

More information

Lecture 8: Period Finding: Simon s Problem over Z N

Lecture 8: Period Finding: Simon s Problem over Z N Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing

More information

One Class of Splitting Iterative Schemes

One Class of Splitting Iterative Schemes One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi

More information

Optimal Coordination of Samples in Business Surveys

Optimal Coordination of Samples in Business Surveys Paper preented at the ICES-III, June 8-, 007, Montreal, Quebec, Canada Optimal Coordination of Sample in Buine Survey enka Mach, Ioana Şchiopu-Kratina, Philip T Rei, Jean-Marc Fillion Statitic Canada New

More information

FLEXOELECTRIC SIGNALS ON RINGS IN TRANSVERSE MOTIONS

FLEXOELECTRIC SIGNALS ON RINGS IN TRANSVERSE MOTIONS Proceeding of the ASME 011 International Deign Engineering Technical Conference & Computer and Information in Engineering Conference IDETC/CIE 011 Augut 8-31, 011, Wahington, DC, USA DETC011-4819 FLEXOELECTRIC

More information

Lecture 7: Testing Distributions

Lecture 7: Testing Distributions CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting

More information

ABSTRACT- In this paper, a Shunt active power filter (SAPF) is developed without considering any harmonic detection

ABSTRACT- In this paper, a Shunt active power filter (SAPF) is developed without considering any harmonic detection Special Iue of International Journal of Advance in Applied Science and Engineering (IJAEAS) ISSN (P): 2348-1811; ISSN (E): 2348-182X Vol. 4, Iue 1,2, March 2017, 34-39 IIST SHUNT ACTIVE POWER FILTER PERFORMANCE

More information

Numerical Simulations of Coriolis Flow Meters for Low Reynolds Number Flows

Numerical Simulations of Coriolis Flow Meters for Low Reynolds Number Flows MAPAN - Journal Numerical of Metrology Simulation Society of of Corioli India, Vol. Flow 26, Meter No. 3, 2011; for Low pp. Reynold 225-235 Number Flow ORIGINAL ARTICLE Numerical Simulation of Corioli

More information

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang

More information

Module 4: Time Response of discrete time systems Lecture Note 1

Module 4: Time Response of discrete time systems Lecture Note 1 Digital Control Module 4 Lecture Module 4: ime Repone of dicrete time ytem Lecture Note ime Repone of dicrete time ytem Abolute tability i a baic requirement of all control ytem. Apart from that, good

More information

FOURIER-BASED METHODS FOR THE SPECTRAL ANALYSIS OF MUSICAL SOUNDS. Sylvain Marchand

FOURIER-BASED METHODS FOR THE SPECTRAL ANALYSIS OF MUSICAL SOUNDS. Sylvain Marchand FOUIE-BAS METHOS FO THE SPECTAL ANALYSIS OF MUSICAL SOUNS Sylvain Marchand Univerity of Bret, Lab-STICC CNS UM 628, 292 Bret, Brittany, France ABSTACT When dealing with muical ound, the hort-time Fourier

More information

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Nonlinear Single-Particle Dynamics in High Energy Accelerators Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction

More information

Copyright 1967, by the author(s). All rights reserved.

Copyright 1967, by the author(s). All rights reserved. Copyright 1967, by the author(). All right reerved. Permiion to make digital or hard copie of all or part of thi work for peronal or claroom ue i granted without fee provided that copie are not made or

More information

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -

More information

Critical Height of Slopes in Homogeneous Soil: the Variational Solution

Critical Height of Slopes in Homogeneous Soil: the Variational Solution Critical Height of Slope in Homogeneou Soil: the Variational Solution Chen, Rong State Key Laboratory of Coatal and Offhore Engineering & Intitute of Geotechnical Engineering, Dalian Univerity of Technology,

More information

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax:

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax: Control Sytem Engineering ( Chapter 7. Steady-State Error Prof. Kwang-Chun Ho kwangho@hanung.ac.kr Tel: 0-760-453 Fax:0-760-4435 Introduction In thi leon, you will learn the following : How to find the

More information

5.5 Application of Frequency Response: Signal Filters

5.5 Application of Frequency Response: Signal Filters 44 Dynamic Sytem Second order lowpa filter having tranfer function H()=H ()H () u H () H () y Firt order lowpa filter Figure 5.5: Contruction of a econd order low-pa filter by combining two firt order

More information

Non-stationary Random Vibration Analysis of Vehicle with Fractional

Non-stationary Random Vibration Analysis of Vehicle with Fractional 13 th National Conference on Mechanim and Machine (NaCoMM7), IISc, Bangalore, India, December 1-13, 7 NaCoMM-7-77 Non-tationary Random Vibration Analyi of Vehicle with Fractional Damping Yunqing Zhang

More information

NAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE

NAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE POLITONG SHANGHAI BASIC AUTOMATIC CONTROL June Academic Year / Exam grade NAME (pinyin/italian)... MATRICULATION NUMBER... SIGNATURE Ue only thee page (including the bac) for anwer. Do not ue additional

More information

Module: 8 Lecture: 1

Module: 8 Lecture: 1 Moule: 8 Lecture: 1 Energy iipate by amping Uually amping i preent in all ocillatory ytem. It effect i to remove energy from the ytem. Energy in a vibrating ytem i either iipate into heat oun or raiate

More information

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam

Solutions. Digital Control Systems ( ) 120 minutes examination time + 15 minutes reading time at the beginning of the exam BSc - Sample Examination Digital Control Sytem (5-588-) Prof. L. Guzzella Solution Exam Duration: Number of Quetion: Rating: Permitted aid: minute examination time + 5 minute reading time at the beginning

More information

SENSORLESS DIRECT TORQUE CONTROL OF RUSHLESS AC MACHINE USING LUENBERGER OBSERVER

SENSORLESS DIRECT TORQUE CONTROL OF RUSHLESS AC MACHINE USING LUENBERGER OBSERVER 5-8 JATIT. All right reerved. SENSORLESS DIRECT TORQUE CONTROL OF RUSHLESS AC MACHINE USING LUENBERGER OBSERVER 1 Khalil. Nabti, K. Abed, H. Benalla 1, Student, Prof. Department of Electrical Engineering,

More information

Alternate Dispersion Measures in Replicated Factorial Experiments

Alternate Dispersion Measures in Replicated Factorial Experiments Alternate Diperion Meaure in Replicated Factorial Experiment Neal A. Mackertich The Raytheon Company, Sudbury MA 02421 Jame C. Benneyan Northeatern Univerity, Boton MA 02115 Peter D. Krau The Raytheon

More information

Fast Convolutional Sparse Coding (FCSC)

Fast Convolutional Sparse Coding (FCSC) Fat Convolutional Spare Coding (FCSC) Bailey ong Department of Computer Science Univerity of California, Irvine bhkong@ic.uci.edu Charle C. Fowlke Department of Computer Science Univerity of California,

More information

Adaptive Control of Level in Water Tank: Simulation Study

Adaptive Control of Level in Water Tank: Simulation Study Adaptive Control of Level in Water Tank: Simulation Study Jiri Vojteek and Petr Dotal Abtract An adaptive control i a popular, o called modern control method which could be ued for variou type of ytem

More information

ISSN: [Basnet* et al., 6(3): March, 2017] Impact Factor: 4.116

ISSN: [Basnet* et al., 6(3): March, 2017] Impact Factor: 4.116 IJESR INERNAIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH ECHNOLOGY DIREC ORQUE CONROLLED INDUCION MOOR DRIVE FOR ORQUE RIPPLE REDUCION Bigyan Banet Department of Electrical Engineering, ribhuvan Univerity,

More information

Damage Diagnosis of Beam-like Structures Based on Sensitivities of Principal Component Analysis Results

Damage Diagnosis of Beam-like Structures Based on Sensitivities of Principal Component Analysis Results Damage Diagnoi of Beam-like Structure Baed on Senitivitie of Principal Component Analyi Reult NGUYEN Viet Ha, GOLNVAL Jean-Claude Univerity of Liege, Aeropace & Mechanical Engineering Department, Structural

More information

Study on the effect of vent on the electroacoustic absorber

Study on the effect of vent on the electroacoustic absorber INTER-NOISE 016 Study on the effect of vent on the electroacoutic aborber Youngeun Cho 1 ; Semyung Wang 1, Gwangju Intitute of Science and Technology, Republic of Korea ABSTRACT Electroacoutic aborber

More information

Clustering Methods without Given Number of Clusters

Clustering Methods without Given Number of Clusters Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,

More information

Equivalent Strain in Simple Shear Deformations

Equivalent Strain in Simple Shear Deformations Equivalent Strain in Simple Shear Deformation Yan Beygelzimer Donetk Intitute of Phyic and Engineering The National Academy of Science of Ukraine Abtract We how that imple hear and pure hear form two group

More information

OBSERVER-BASED REDUCED ORDER CONTROLLER DESIGN FOR THE STABILIZATION OF LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

OBSERVER-BASED REDUCED ORDER CONTROLLER DESIGN FOR THE STABILIZATION OF LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS International Journal o Computer Science, Engineering and Inormation Technology (IJCSEIT, Vol.1, No.5, December 2011 OBSERVER-BASED REDUCED ORDER CONTROLLER DESIGN FOR THE STABILIZATION OF LARGE SCALE

More information

ANALYTICAL BEARING MODEL FOR ANALYSIS OF INNER LOAD DISTRIBUTION AND ESTIMATION OF OPERATIONAL LUBRICATION REGIME

ANALYTICAL BEARING MODEL FOR ANALYSIS OF INNER LOAD DISTRIBUTION AND ESTIMATION OF OPERATIONAL LUBRICATION REGIME 58 th ICMD 017 6-8 September 017, Prague, Czech Republic ANALYTICAL BEARING MODEL FOR ANALYSIS OF INNER LOAD DISTRIBUTION AND ESTIMATION OF OPERATIONAL LUBRICATION REGIME Jakub CHMELAŘ 1, Votěch DYNYBYL

More information