DESIGN OF MECHANICAL SYSTEMS HAVING MAXIMALLY FLAT RESPONSE AT LOW FREQUENCIES

Similar documents
8.1 Force Laws Hooke s Law

About the definition of parameters and regimes of active two-port networks with variable loads on the basis of projective geometry

Course support: Control Engineering (TBKRT05E)

EE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng

Mutual capacitor and its applications

Fourier Series Summary (From Salivahanan et al, 2002)

Chapter 1: Basics of Vibrations for Simple Mechanical Systems

Electric Power System Transient Stability Analysis Methods

Q ESTIMATION WITHIN A FORMATION PROGRAM q_estimation

Chapter 10 Objectives

A New Algorithm for Reactive Electric Power Measurement

Reducing Vibration and Providing Robustness with Multi-Input Shapers

Chapter 2: Introduction to Damping in Free and Forced Vibrations

Chapter 10: Sinusoidal Steady-State Analysis

List Scheduling and LPT Oliver Braun (09/05/2017)

TELLEGEN S THEOREM APPLIED TO MECHANICAL, FLUID AND THERMAL SYSTEMS

Chapter 6 1-D Continuous Groups

Scattering and bound states

ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER

2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all

General Properties of Radiation Detectors Supplements

CHAPTER 1: INTRODUCTION

New upper bound for the B-spline basis condition number II. K. Scherer. Institut fur Angewandte Mathematik, Universitat Bonn, Bonn, Germany.

13 Harmonic oscillator revisited: Dirac s approach and introduction to Second Quantization

ma x = -bv x + F rod.

. The univariate situation. It is well-known for a long tie that denoinators of Pade approxiants can be considered as orthogonal polynoials with respe

Force and dynamics with a spring, analytic approach

Feature Extraction Techniques

NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT

Uniform Approximation and Bernstein Polynomials with Coefficients in the Unit Interval

Virtual physics laboratory for distance learning developed in the frame of the VccSSe European project

Using a De-Convolution Window for Operating Modal Analysis

Chapter 10: Sinusoidal Steady-State Analysis

Automated Frequency Domain Decomposition for Operational Modal Analysis

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Mechanical Engineering 2.010: Systems Modeling and Dynamics III. Final Examination Review Problems

Estimation of ADC Nonlinearities from the Measurement in Input Voltage Intervals

THE KALMAN FILTER: A LOOK BEHIND THE SCENE

Now multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get:

Kernel Methods and Support Vector Machines

ACTIVE VIBRATION CONTROL FOR STRUCTURE HAVING NON- LINEAR BEHAVIOR UNDER EARTHQUAKE EXCITATION

Introduction to Discrete Optimization

Work, Energy and Momentum

Study Committee B5 Colloquium 2005 September Calgary, CANADA

The accelerated expansion of the universe is explained by quantum field theory.

Envelope frequency Response Function Analysis of Mechanical Structures with Uncertain Modal Damping Characteristics

ECEG 351 Electronics II Spring 2017

Optimization of ripple filter for pencil beam scanning

A Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair

a a a a a a a m a b a b

Chapter 2. Small-Signal Model Parameter Extraction Method

Extension of CSRSM for the Parametric Study of the Face Stability of Pressurized Tunnels

DETECTION OF NONLINEARITY IN VIBRATIONAL SYSTEMS USING THE SECOND TIME DERIVATIVE OF ABSOLUTE ACCELERATION

Lecture #8-3 Oscillations, Simple Harmonic Motion

&& && F( u)! "{ f (x)} = f (x)e # j 2$ u x. f (x)! " #1. F(u,v) = f (x, y) e. f (x, y) = 2D Fourier Transform. Fourier Transform - review.

Chapter 28: Alternating Current

ELG3311: Assignment 3

EFFECT OF SURFACE ASPERITY TRUNCATION ON THERMAL CONTACT CONDUCTANCE

2.003 Engineering Dynamics Problem Set 2 Solutions

Fast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials

RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS

th Annual IEEE Power Electronics Specialists Conference Aachen, Germany, Parallel Connection of Piezoelectric Transformers

In the session you will be divided into groups and perform four separate experiments:

Ufuk Demirci* and Feza Kerestecioglu**

Complexity reduction in low-delay Farrowstructure-based. filters utilizing linear-phase subfilters

Determination of Active and Reactive Power in Multi-Phase Systems through Analytical Signals Associated Current and Voltage Signals

Multi-Scale/Multi-Resolution: Wavelet Transform

Four-vector, Dirac spinor representation and Lorentz Transformations

Supplementary Information for Design of Bending Multi-Layer Electroactive Polymer Actuators

Understanding Machine Learning Solution Manual

Lecture 16: Scattering States and the Step Potential. 1 The Step Potential 1. 4 Wavepackets in the step potential 6

Generalized Sampling Theorem for Bandpass Signals

m potential kinetic forms of energy.

Proc. of the IEEE/OES Seventh Working Conference on Current Measurement Technology UNCERTAINTIES IN SEASONDE CURRENT VELOCITIES

Physics 2107 Oscillations using Springs Experiment 2

ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics

Nonlinear Stabilization of a Spherical Particle Trapped in an Optical Tweezer

Energy and Momentum: The Ballistic Pendulum

Kinetic Theory of Gases: Elementary Ideas

Q5 We know that a mass at the end of a spring when displaced will perform simple m harmonic oscillations with a period given by T = 2!

In this chapter we will start the discussion on wave phenomena. We will study the following topics:

Solving initial value problems by residual power series method

13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices

COS 424: Interacting with Data. Written Exercises

Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi.

III.H Zeroth Order Hydrodynamics

EE 330 Lecture 31. Basic amplifier architectures. Common Emitter/Source Common Collector/Drain Common Base/Gate

The Wilson Model of Cortical Neurons Richard B. Wells

i ij j ( ) sin cos x y z x x x interchangeably.)

Chapter 8 Deflection. Structural Mechanics 2 Dept of Architecture

THERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS = (1) ln = COMPOSITES & POLYCON 2009

Kinetic Theory of Gases: Elementary Ideas

JOURNAL OF PHYSICAL AND CHEMICAL SCIENCES

On Conditions for Linearity of Optimal Estimation

1 Analysis of heat transfer in a single-phase transformer

SOLUTIONS. PROBLEM 1. The Hamiltonian of the particle in the gravitational field can be written as, x 0, + U(x), U(x) =

IDAN Shock Mount Isolation Vibration Study November 1, The operation of shock and vibration isolation base plate

On the characterization of non-linear diffusion equations. An application in soil mechanics

A Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine. (1900 words)

Optical Properties of Plasmas of High-Z Elements

BEF BEF Chapter 2. Outline BASIC PRINCIPLES 09/10/2013. Introduction. Phasor Representation. Complex Power Triangle.

Transcription:

DESIGN OF MECHANICAL SYSTEMS HAVING MAXIMALLY FLAT RESPONSE AT LOW FREQUENCIES V.Raachran, Ravi P.Raachran C.S.Gargour Departent of Electrical Coputer Engineering, Concordia University, Montreal, QC, CANADA, HG M8. Departent of Electrical Coputer Engineering, Rowan University, Glassboro, New ersey, U.S.A., 8. École de Technologie Supérieure, University of Quebec, Montreal, CANADA, HC. Abstract This paper discusses a ethod of designing translational rotational echanical systes having axially flat response at the lower end of the frequency spectru. Introduction One of the ethods of analyzing a syste is to obtain an analogous electrical networ then use nown analysis techniques to obtain the required inforation []. The syste ay be echanical (translational rotational, fluid, theral or of any other type. The analysis can be carried out either in the tie-doain or in the frequency-doain. Also, Laplacetransfor techniques can be conveniently used in all cases, when the syste considered is linear []. It is the purpose of this paper to show that the sae concept can be applied in the design of such systes. Specifically, since a considerable aount of literature exists in the design of electrical low-pass filters, it can be ade use of in the design of siilar echanical systes also. Magnitude of the Transfer function One starts with the agnitude of the transfer function H ( j V out Ω ( V n in Ω where V in is the input velocity, V out is the output velocity, Ω is the noralized frequency ω, ω c ω is the operating frequency, ω c is the cut-off frequency, n is the order of the syste. Since the velocity is an acrossvariable, the analogous variable in the electrical syste is the voltage therefore Eq.( can be considered to be a voltage transfer function. This perits one to use the theory of electrical filters to study the properties of Eq.( its design. Soe iportant properties of the agnitude transfer function considered. Property : At zero frequency, that is, at Ω, the response is always unity. Property : At unity noralized frequency, that is, at Ω, (which is the cut-off frequency, the response is

always, irrespective of the order of the syste. Property : It is readily verified that the first (n- derivatives at Ω will always be zero. This gives rise to the axially flat response in the neighbourhood of Ω. These three properties do not give any inforation regarding the selection of order n. Property : The agnitude response is onotonically decreasing. Property 5 : It can be shown that d H ( jω dω Ω -.55n ( This gives the slope of the agnitude response at the cut-off frequency. This is a possibility for the specification of the syste. Note that n has to be an integer. If n coes out to be a fraction, it is converted into the next higher integer. Alternatively, the order n could be deterined by specifying the agnitude response at a frequency Ω >. In this case also, if n coes out to be a fraction, it is converted into the next higher integer. Table I gives the agnitude of the transfer functions for different orders for Ω.5 Ω.. Deterination of the transfer function After the deterination of n, the order of the syste, the next step is to generate the appropriate transfer function. This is carried out with the help of the coplex frequency given by S σ jω ( S can also be considered as a Laplace-transfor variable. This perits one to write Table I Magnitudes of the transfer function for different orders at Ω.5 Ω. radians. Order Ω.5 Ω. Magnitude in db Magnitude in db.557-5..7-7..6-7.8.5 -..8 -.9. -8..98 -.5.69 -. 5.6-7.68.5 -. 6.875 -.6. -.5 H j Ω H(S. H(-S ( Ω -S ( ( S n It is evident that when the denoinator is equated to zero, it contains n roots S, S,...., S n-, S n, each one having agnitude of unity. In other words, all of the lie on the unit circle in the coplex-plane S. Since the function has to be stable, H(S shall contain all its poles in the left-half of the S -plane, or H ( S (5 D( S where D(S is a strictly Hurwitz polynoial. It is readily shown that the roots of the equation (- n S n (6a are given by π (- j n S e,,,,n (6b Fro these roots, D(S a a S a S a S.. (6c can be readily obtained these are given upto n 6 in Table. For details, please see []

Table Coefficients of the denoinator polynoial of the transfer function upto order 6. Order a a a a, a. a a a a a a a a.66 a. 5 a 5 a a a.668 a a 5.668 6 a 6 a a 5 a.867 a a.9785 a 9.6 Ipleentation of the generated Transfer Function We are ainly interested in a echanical syste, which results in the generated transfer function. One possibility is to design a passive electrical networ then obtain the corresponding echanical syste by appropriate odeling. One such networ is shown in Fig.. Vin LC - networ R L Vout ' ' Fig.: The passive electrical networ considered. The LC-networ consists of only inductors capacitors the utual inductances are absent. (There are other possibilities they are not considered here. It can be readily shown that ( H S V - y out (7 Vin y R L where y input adittance at - with - short-circuited R L reoved, -y feedbac transfer adittance between - - with R L reoved - shortcircuited. It can be shown that y y are the ratio of even to odd polynoials [], because they are obtained for LCnetwors only. Writing D(S M(S N(S (8 where M(S Even part of D(S, N(S Odd part of D(S. It can be shown that -y N( S (9 M(S -y ( N( S with R L Ω (taen without any loss of generality. M(S Exping around S N( S into a continued fraction expansion, the required LC-networ is obtained. Once the electrical networ is realized, the corresponding echanical syste is obtained by replacing the resistor by a daper, the inductor by a spring the capacitor by a ass (for translational systes a flywheel (for rotational systes. This is illustrated by the following two exaples where the orders considered are n n respectively. Exaple Let the order n be. The transfer function will be

H ( S S S S This yields y S S S y S S y has to be exped into continued fractions this yields y ( S S S S This results in the electrical networ shown in Fig.(a..5H.5H F ' ' Fig.(a: A third-order low-pass filter having axially flat response around ω. Ω The corresponding noralized echanical translation syste is shown in Fig.(b the noralized echanical rotational syste is shown in Fig.(c. n n n b where n, n, n b. Fig.(b: A third-order noralized echanical translational syste corresponding to the electrical networ of Fig.(a. n n n B where n, n, n B. Fig.(c: A third-order noralized echanical rotational syste corresponding to the electrical networ of Fig.(a. Exaple Let the order n be. The transfer function will be H ( S S a S a S a S where a a.66, a.. This yields - y.6s.66s S.S y.6s.66s The continued fraction of y around S yields y C S L S C S LS where L.57, C.5776, L.89, C.868. This results in the electrical networ shown in Fig.(a.

H L H L C F ' ' Fig.(a: A fourth-order low-pass filter having axially flat response around ω. F C Ω The corresponding noralized echanical translational syste is shown in Fig.(b the noralized echanical rotational syste is shown in Fig.(c. n n n n b where n.98797, n.658, n.86, n.5776, b. Fig.(b: A fourth-order noralized echanical translational syste corresponding to the electrical networ of Fig.(a. n n n n B where n.988, n.658, n.86, n.5776, B. Fig.(c: A fourth-order noralized echanical rotational syste corresponding to the electrical networ of Fig.(a. Scaling Having obtained the noralized values in the design, one has to denoralize the in order to obtain the actual values. This process is nown as scaling. There are two types of scaling, naely (a Ipedance scaling (b Frequency scaling. Both these can be carried out independently or siultaneously, depending on the requireents. They are discussed below: (a Ipedance scaling : This is necessitated by the fact that b or B is taen to be unity initially. By ipedance scaling, the ipedance of every coponent has to be scaled by the sae quantity. However, the transfer function reains unaltered. Let be the ipedance scaling factor. It is readily seen that b B becoe b B S respectively. The ipedances S becoes S S respectively. These yield the results that the spring constants becoes respectively. Considering the asses, the ipedances S becoe respectively, S S S giving us the new values as after ipedance scaling. These new values ay be further subjected to frequency scaling discussed below. (b Frequency scaling : This is required because the operating frequency has been noralized to unity. We can s readily put S, where β is the β frequency scaling factor. The quantities b B reain unchanged. By a siilar treatent as given in the case of 5

Ipedance scaling, it can be shown that the quantities,, becoe β, β, β β respectively. (c Cobined Ipedance Frequency scaling : Both the scaling operations can be cobined the final results are given in Table III. The quantities in colun are given after they are subjected to frequency scaling only. The quantities in colun can be obtained fro those in colun are obtained by putting those in colun are obtained by putting β in those given in colun. Table Effect of scaling on the coponent values Quantity b B Ipedance scaling Frequency scaling b b B B β β β β Cobined scalings b B β β β β translational syste considered in Fig.(a with b N.s/ a cutoff frequency of 5 radians/second. This gives β 5. The various denoralized coponent values coe out to be: n N/, n 9 8 N/ n g. If the fourth 9 order syste considered in Fig.(a is considered with the sae scaling factors, the various denoralized will be n 5.98 N/, n 8.88 N/, n.559 g n.5 g. Exaple In this exaple, we shall consider the denoralization of the rotational syste considered in Fig.(c. Let B N..s the cut-off frequency be 5 radians/second. The various coponent values coe out to be n N./radians, n N./radians, n 6 g.. When the fourth order rotational syste given in Fig.(c is considered, the various denoralized coponent values coe out to be n.796 N./radians, n.776 N./radians, n.5 g. n 7.86 g.. The responses for the orders are shown in Fig.. By using these scalings, we shall now deterine the actual values to be used for the two exaples considered above. Exaple In this exaple, we shall consider the scaling of a third order 6

Fig.: Magnitude responses for n n. Conclusions In this paper, it is shown that the design ethod used in electrical wave filters can be effectively used in the design of echanical low-pass systes. Specifically, axially-flat response near the zero frequency has been considered the corresponding echanical syste has been designed. Though lower order systes are considered, the sae concept can be used to design higher order systes as well. It is envisaged that the sae concept can be used when other types of responses are considered. References [].F.Lindsay V.Raachran, Modeling Analysis of Linear Physical Systes, Weber systes Inc., 99. [] V.Raachran Ravi P.Raachran, Tellegan s Theore applied to Mechanical, Fluid Theral Systes, Annual Conference Exposition of Aerican Society for Engineering Education, New Mexico, Session 79, une -7,. [] V.Raachran, Course notes for Modeling Analysis of Linear Physical Systes, Concordia University, Montreal, QC, CANADA, 999. 7