Energy-Rate Method and Stability Chart of Parametric Vibrating Systems

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Energy-Rate Method and Stability Chart of Parametric Vibrating Systems"

Transcription

1 Reza N. Jazar Department of Mechanical Engineering Manhattan College New York, NY, A M. Mahinfalah Department of Mechanical Engineering Milwaukee chool of Engineering Milwaukee, WI, A N. Mahmoudian and Department of Aerospace Engineering, Virginia Institute of Technology Blacksurg, VI, A Energy-Rate Method and taility Chart of Parametric Virating ystems The Energy-Rate method is an applied method to determine the transient curves and staility chart for the parametric equations. This method is ased on the first integral of the energy of the systems. Energy-Rate method finds the values of parameters of the system equations in such a way that a periodic response can e produced. In this study, the Energy-Rate method is applied to the following forced Mathieu equation: && y + hy& + ( 1 + cos( rt) ) y = sin ( rt) This equation governs the lateral viration of a microcanilever resonator in linear domain. Its staility chart in the -r plane shows a complicated map, which cannot e detected y perturation methods. Keywords: energy-rate method, staility chart, parametric virations M. A. Rastgaar Department of Mechanical Engineering Virginia Institute of Technology Blacksurg, VI, A Introduction 1 In a recent paper, an energy-ased method, called Energy-Rate method has een developed to investigate the staility chart of the parametric differential equations. The method has een introduced and applied to the Mathieu equation. && x+ a x x cos( t) = Application of Energy-Rate method has shown that there are new periodic curves in stale domain of a- plane that have not een detectale y conventional methods (Jazar, 4) (Butcher and inha, 1998). In this study, the Energy-Rate method will e utilized to detect the staility diagram of a forced parametric equation, and a function to e used as a gauge function for determining relative staility in parametric space. In general, a solution is unstale if it grows unoundedly as time approaches to infinity, while a solution is stale if it remains ounded as time approaches to infinity. A periodic solution is neutral and is considered as a special case of a stale solution (Hayashi, 1964). The term parametric equations or parametrically-excited oscillation is generally used for phenomena in which a parameter of the system is time-varying, usually periodically. In early experiments, such systems showed resonance when the variale parameter has a frequency equal to the doule of the natural frequency of the linearized system. Existence of periodic viration with frequencies equal to multiple of the principal resonance frequency is an essential characteristic of parametric equations. Parametric equations can e seen in many ranches of applied science. They have een the suject of a vast numer of investigations since the eginning of the last century (Whittaker and Watson, 197) (Van der Pol and trutt, 198). Investigation of parametric phenomena has een started since 1831 when Faraday produced wave motion in water y virating an immersed plate (Iwanowski, 1965). hortly afterwards, the importance of parametric phenomena was revealed y many scientists in the late 19th and the eginning of th century (Richards, 1983). The Mathieu equation expressed in the form: dx && x+ a x x cos( t) = x& = (1) dt is a special case of Hill equation, && x+ f(t) x =, is the simplest and the first well defined parametric differential equation, in which a and are constant parameters. Depending on the value of a and, Mathieu equation has ounded as well as unounded solutions. Because the Mathieu equation is linear, it s staility is not dependent on initial conditions. The values of the pairs (a, ) that produces periodic responses at the oundaries of stale and unstale regions, make a set of continuous curves in the parametric plane, as shown in Figure 1 (McLachlan1947). taility charts of timevarying systems, including Mathieu equation, have een investigated using a variety of methods and concepts, such as Lyapunov exponents, Poincaré maps, Liapunov-Floquet transformations, and perturative Hamiltonian normal forms (inha and Butcher, 1997) (Guttalu and Flashner, 1996) (Butcher and inha, 1995). Paper accepted April, 8. Technical Editor: Domingos Alves Rade. 18 / Vol. XXX, No. 3, July-eptemer 8 ABCM

2 Energy-Rate Method and taility Chart of T Γ= 1 1 Edt & x g ( x,x,t) dt T = T & &. (6) T ce se 1 ce1 se ce se 3 ce 3 se 4 ce 4 se 5 nstale tale Figure 1. Mathieu staility chart ased on the numerical values, generated y (McLachlan, 1947). Energy-Rate Method The Energy-Rate method can e applied to the following general parametric differential equation: ( ) ( ) && x+ f x + g x, x&, t = () where, f(x) is a single variale function, and g ( x,x,t & ) is a smooth periodic time varying function with the following conditions. g (,,t) =. (3) ( + ) = ( ) g x,x,t & T g x,x,t & (4) Furthermore, the functions f and g can depend on finite sets of parameters. Equation () may e assumed as a model of a unit mass attached to a spring, acted upon y a non-conservative force g ( x,x,t & ). Introducing the kinetic, potential, and mechanical energies for the system, T( x& ) = x& /, V( x) = f( x) dx, and E = T( x& ) + V( x), we may transform Equation () to an integral of energy as follows: & d d 1 = = & + ( ) = ( ) dt dt & & &. (5) E ( E) x f( x) x dt x g x, x, t & & &, which is equal to the time derivative of the mechanical energy of the system, represents the instantaneous rate of generated or asored energy y the applied g x,x,t &. Instantaneous generation or asorption of energy The function E = x g( x, x, t) force ( ) depends on the value of parameters, time, x, and x&. If the overall value of the rate of energy in an excitation period is negative for a set of parameters and a nonzero response x(t), the response of the system shrinks along the path of x(τ), t <τ< t + T, and the amplitude of oscillation decreases. On the other hand, if the overall value of the rate of energy in an excitation period is positive for a set of parameters and a nonzero response x(t), then the response of the system expands along the path of x(τ), t <τ< t + T, and the amplitude of oscillation increases. We define an Energy-Rate function as an average of the integral of the energy rate: ce 5 se 6 The Energy-Rate function is zero for free viration of a conservative system, as well as the steady state (T/n)-periodic response of the system (). Equation () includes two parameters, a and, so the parameter space is two-dimensional. The numer of parameters determines the dimension of the parameter space. Determination of staility diagram is carried-out y finding the pairs of (a, ) located on the oundary of stale regions. For a steady state T-periodic response, we pick a pair of (a, ) and integrate Equation () numerically to evaluate Γ. If Γ> then, (a, ) elongs to an unstale region where energy is eing inserted to the system. However, if Γ<, then (a, ) elongs to a stale region in which energy is eing extracted from the system. On the common oundary of these two regions, Γ=, and (a, ) elongs to a transition curve. Let us fix and search for a on a transition curve such that its left hand side is stale and its right hand side is unstale. o, if a point (a, ) shows that Γ is less than zero, the point is in the stale region. In such a case, increasing a will increase Γ. On the other hand, if Γ is greater than zero, the point is in the unstale region so, decreasing a will decrease Γ. The Energy-Rate function Γ can e used as a scale for converging a to an appropriate value, such that Γ ecomes zero. Applying this procedure, we find a point on a transient curve with a stale region on the left and an unstale region on the right. Varying y an increment and repeating the calculation, we will find another periodic point. The transient curve can e found when is varied enough to cover the entire domain of interest. It is often useful to have a prior understanding of ehavior of the system. Reversing the strategy and searching for transient curves whose right hand side is stale and left hand side is unstale, completes the staility diagram of the system. This procedure can e arranged in an algorithm to e set up in a computer program, as follows: tep 1 - set a tep - set, equal to some aritrary small value tep 3 - solve the differential equation numerically tep 4 - evaluate Γ tep 5 - decrease / increase a, if Γ > / Γ <, y a small increment tep 6 - the increment of a must e decreased if Γ (a i ) Γ (a i-1)< tep 7 - save a and when Γ <<1 tep 8 - while < final, increase and go to tep 3 tep 9 - reverse the decision in tep 5 and go to tep 1 An alternative method is to set a domain of interest in the parameter plane and then, evaluate Γ for the domain using a fine resolution. A two-dimensional matrix will e found for a twoparameter equation. Each element of the matrix would e the value of Γ for the associated pair of parameters. sing this matrix, a threedimensional surface can e plotted to illustrate the Energy-Rate function. This plot shows the relative staility of different points in the parameter plane. Intersection of the surface with plane Γ= designates the transition curves, while the domain elow it indicate stale regions, and the domain aove this plane indicate instale regions. The accuracy of this analysis depends on the resolution in the parametric plane. There are some advantages in Energy-Rate method over classical perturation methods. First, Energy-Rate method can e applied to nonlinear parametric equations as well as linear ones. J. of the Braz. oc. of Mech. ci. & Eng. Copyright 8 y ABCM July-eptemer 8, Vol. XXX, No. 3 / 183

3 econd, it can e applied to oth, small and large values of parameters. Third, the values of the parameters for a periodic response can e found faster and more accurately than in other methods. Fourth, application of the Energy-Rate method for detection of periodic solutions can e adjusted to detect any periodic response, not necessarily on the oundary of stale and instale regions. Fifth, the Energy-Rate function can e used as a gauge function to compare the relative staility of different points of parameter space. The greatest disadvantage of perturation methods is that they can usually e applied to very small values of the parameters. Also, the accuracy of the results of most asymptotic methods cannot e increased y increasing the degree of polynomials or the numer of terms in the pertured solution. In the next section we apply the Energy-Rate method to the Mathieu equation and to forced Mathieu type equation to show that the Energy-Rate method provides much more accurate staility diagram as compared to perturation methods. Γ Figure. Plot of π for Mathieu equation as a function of a for =1. taility Diagram of Mathieu Equation To show the applicaility of the Energy-Rate method, Mathieu equation is examined as the first example. Mathieu equation is a well investigated equation, and ecause of its linearity, there exist some effective methods for finding its staility diagram. Transition curves of Mathieu equation are π and π-periodic. The analytic solution of the equation on transient curves are availale and are called cosine and sine elliptic functions (Richards, 1983). The first three even and odd cosine and sine elliptic functions, ce n, se n, ce n+1, and se n+1 are shown in Figure 1. The unstale tongues are ounded with two transition curves starting from integer roots of a n =, n=, 1,, 3,. It can e shown that although the odd order cosine and sine elliptic functions are not symmetric with respect to, the overall staility diagram looks symmetric. ince the parametric resonance of Mathieu equation happens in periods T and T, the ound of integral may e defined from zero to T to include T and T periodic solutions. The Energy-Rate function for Mathieu Equation (1) would e a two-parameter equation. Γ Figure 3(a). taility surface of Mathieu Equation (1). a π 1 Γ =Γ ( a,) = ( x x cos( t )) dt π & (7) Consider a case in which =1, then Γ is only a function of a. This function is plotted in Figure using 1 points in the interval <a<45. It can e seen that the curve Γ intersects the a axis exactly at the same points where the line =1 intersects the transition curves of the Mathieu staility diagram (see Figure 1). Evaluating Γ over a domain of a and generates a three dimensional staility surface. Figure 3(a) displays the staility surface of Mathieu equation over the area < a< 1 and < < 1. Based on Energy-Rate method, wherever Γ is negative, Mathieu Equation (1) is stale, and wherever Γ is positive, it is unstale. The staility surface has een found y integrating over T=π. Therefore, the intersection of the staility surface with the plane Γ=, indicates π or π-periodic curves, ecause a π-periodic function is π-periodic as well. Figure 3() depicts a top view of the staility surface. It also illustrates the intersection curves of the staility surface with plane Γ= a Figure 3(). π-periodicity curves of Mathieu Equation (1). If a periodic curve separates stale and unstale regions, it is called a transition curve, otherwise it is a splitting curve. A splitting curve is a periodic curve emedded in a stale region. Figure 4 depicts the π and π-periodic transition curves of Mathieu equation, along with a splitting curve in each stale region (Jazar, 4). 184 / Vol. XXX, No. 3, July-eptemer 8 ABCM

4 Energy-Rate Method and taility Chart of h=1. h=.5 h=. h= a Figure 4. π-periodicity curves of Mathieu Equation (1), intersection of the staility surface and plane Γ =. r Figure 5. Transient curves and instaility region for Equation (1). taility Diagram of A Forced Mathieu Type Equation As a second example, forced viration of a Mathieu type system will e studied in this section. Consider the following forced and parametrically excited equation: ( ( )) ( ) && y+ hy& cos rt y = sin rt. (1) which is the nondimensionalized governing equation of a linearized model of a micro-electro-mechanical system (Luo, ) (Jazar, 6). The system is a cantilever resonator activated y an alternating electric field. Equation (1) is dependent on three independent parameters. Therefore the parameter space of the system is a three-dimensional r--h space. Assuming h=, the condition for parametric resonance, which states that n times of the forcing frequency is equal to twice of the natural frequency of the unforced system, would e 1 = nr,n N. Therefore, transition curves, if there is any, must start from r = 1/n,n N on the line =. The periodic curves corresponding to n=1 makes the principal instaility domain, if they are transition curves. Assuming a periodic solution, y = Y n n sin( nrt +ϕ), and = applying multiple scale perturation method provides the following relationship for the periodic condition with Y =. ( ( ) ) 1 r = h 4± h 4 4h 1 + 4, n N. () n Traditionally, the periodic condition is used to determine the periodic curves in parameter plane r-. The periodic curves need to e checked to determine if they are transition or splitting curves. The periodic curves () corresponding to n=1 are plotted in Figure 5 for different values of h. It is expected that two periodic curves come out of r=1 and make the principal instaility tongue. Time response of the undamped system for selected points (r, ) = (.9,.1), and (1.,.1), are found and shown in Figure 5. It is concluded that Equation () are transition curves, and the region surrounded y curves coming out of r=1 is unstale. Figure (5) also illustrates that the instaility domain shrinks y increasing the dimensionless damping, h. Applying Poincare-Lindstad method provides the following equations for no damping transient curves will e found 4 1 = 1 r + O(r ) 1 r = + O( r ) (3) which are in agreement with Equation () and Figure 5. It must e mentioned that not only the curves derived from perturation analysis are not necessarily transition curves, the staility ehavior of the system is also questionale when parameters are close to a curve. Furthermore, we know that although Equations () and (3) might e used to determine the periodic curves and staility regions in r- space around =, they cannot e used to determine the gloal staility regions of the system. Defining the Energy-Rate function as follow: π /r r Γ= y( sin ( rt) hy y cos( rt) ) dt π &. (4) and applying Energy-Rate method provides an exact staility diagram. To present a etter view of staility diagram, we study the system in - plane instead of r- plane. However, it is also possile to apply a transformation rt =τ, and convert Equation (1) into the form: where, ( ( )) ( ) y + Hy + A+ Bcos τ y = Bsin τ y = dy/rdτ, H A = ( 1 ) /r, = h/ r, B = /r. Therefore, the Energy-Rate function may e defined as ( ( ) ( )) * πγ = π y Bsin rt Hy Bycos τ dτ. A three dimensional illustration of staility surface πγ /r is depicted in Figure 6(a) to provide a relative staility view. The staility surface is made y 5 1 points. The intersection of the staility surface with Γ= specifies the transition curves of the system that are shown in Figure 6(). Figure 6(c) shows a top view of the staility surface along with the staility diagram. The staility characteristic of stale and unstale regions of Figure 6() is clear J. of the Braz. oc. of Mech. ci. & Eng. Copyright 8 y ABCM July-eptemer 8, Vol. XXX, No. 3 / 185

5 y investigating Figures 6(a) or 6(c). It can also e determined y numerical solution of the equation for a sample point within each region. A deep and an elevated regions are associated to strong stale and an unstale region: a stale region for large values of around r=1, and an unstale region corresponding to large values of and. There is also an unstale region connected to =1 which is the principal instaility region that could e predicted y perturation methods approximately. Figure 6(c). Top view of the staility surface for the parametric Equation (1) along with transition curves. πγ/r Figure 6(a). taility surface for the parametric Equation (1). πγ/r Figure 7(a). Magnification of staility surface at the principal instaility region of the parametric Equation (1). Figure 6(). taility diagram for the parametric Equation (1). Because of the importance of the principal instaility region, a magnified view of the staility surface for <<.4 and 1<<.5 is shown in Figure 7(a). A etter view of the transition curves is depicted in Figure 7(). A top view of the staility surface along with the transition curves are also shown in Figure 7(c) to highlight the staility characteristic of Figure 7(). In order to etter understand the ehaviour of the system when parameters are around the principal instaility region, time responses of the system for a few selected points in different regions are shown in Figure 8. The examined points are marked in Figure 7(). Figure 7(). Transition curves around the principal instaility region of the parametric Equation (1). 186 / Vol. XXX, No. 3, July-eptemer 8 ABCM

6 Energy-Rate Method and taility Chart of Perturation method Energy-Rate method Figure 7(c). Top view of the staility surface of Equation (8) along with transition curves around the principal instaility region. Figure 9. Comparison of the result of perturation and Energy-Rate method close to the principal instaility region of Equation (8). Comparison of the results of perturation analysis in Equation (9) or (1) with the results of Energy-Rate method is shown in Figure 9. From a design point of view, the strong stale region related to r 1 and large values of are important. This is proaly the est area to set the parameters of the system governed y Equation (8). To clarify the oundary of the stale and unstale regions in this area, a three dimensional illustration of staility surface, and a two dimensional staility diagram are shown in Figures 1-1. πγ/r Y =.4 =1.5 h=. =.13 =1.8 Figure 1. Magnification of staility surface around the strong stale region of the parametric Equation (8). t =.37 =.4 =.15 =.45 =.5 =1. Figure 8. Time response of the parametric Equation (8) for some points of the parameter plane. Figure 11. Magnification of transition curves around the strong stale region of the parametric Equation (8). J. of the Braz. oc. of Mech. ci. & Eng. Copyright 8 y ABCM July-eptemer 8, Vol. XXX, No. 3 / 187

7 are also found y applying the Energy-Rate method. Threedimensional illustration of staility surface shows a sense of relative staility of parameter plane. The intersection of a zero plane with staility surface determines the periodic and transition curves of the system. The procedure is sometimes time consuming, ut it generates a more efficient and more accurate staility diagram. Figure 1. Top view of the staility surface of the parametric Equation (8) along with transition curves around the strong stale region. Conclusion The Energy-Rate function, Γ= x g( x,x,t) dt /T & & is used to T generate a three-dimensional staility surface for two parametric equations example. The value of the function can e determined numerically and compared to zero for a parametric equation such as && x+ f ( x) + g( x, x&, t, a, ) =. The Energy-Rate method may e used to investigate the staility of the parametric equations. Procedure of searching for periodic curves in the parameter space is presented as an algorithm. To validate the algorithm, it was applied to the Mathieu equation, and its staility chart was developed. The results are in agreement with a previously reported results. In the second example, the method was applied to a forced and parametrically excited equation, which is ased on the model of a physical system. Multiple time scale, and Poincaré-Lindstad perturation methods were applied and periodic curves for small values of parameters were found. The staility surface and staility diagram of the system References Butcher E.A., and inha.c., 1998, ymolic Computation of Local taility and Bifurcation urfaces for Nonlinear Time-Periodic ystems, Nonlinear Dynamics, 17, 1-1. Butcher, E.A., and inha,.c., 1995, On the Analysis of Time - Periodic Nonlinear Hamiltonian Dynamical ystems, AME Design Engineering Technical Conferences DE 84-1, Vol. 3, part A, Guttalu, R.., and Flashner, H., 1996, taility Analysis of Periodic ystems y Truncated Point Mapping, Journal of ound and Viration, 189, Hayashi, C., 1964, Nonlinear Oscillations in Physical ystems, Princeton. Iwanowski, E.R.M., 1965, On the Parametric Response of tructures, Applied Mechanics Review, 18(9), Jazar, Reza N., 4, taility Chart of Parametric Virating ystems, sing Energy Method, International Journal of Nonlinear Mechanics, 39(8), pp Jazar, Reza N., 6, Mathematical Modeling and imulation of Thermoelastic Effects in Flexural Microcantilever Resonators Dynamics, accepted for pulication in the Journal of Viration and Control. Luo, A. C. J.,, Chaotic Motion in the Generic eparatrix Band of a Mathieu-Duffing Oscillator with a Twin-Well Potential, Journal of ound and Viration, 48(3), McLachlan, N. W., 1947, Theory and Application of Mathieu Functions, Clarendon, Oxford. P., Reprinted y Dover, New York, Richards, J. A., 1983, Analysis Periodically Time Varying ystems, pringer-verlag, New York. inha.c., and Butcher E.A, 1997, ymolic Computation of Fundamental olution Matrices for Linear Time-Periodic Dynamical ystems, Journal of ound and Viration, 6(1), Van der Pol, B.,, and trutt, M.J.O., 198, On the taility of the solution of Mathieu s Equation, The Philosophical Magazine, vol. 5, p. 18. Whittaker, E.T., and Watson, G.N., 197, A Course of Modern Analysis, Camridge niv. Press, 4th ed. 188 / Vol. XXX, No. 3, July-eptemer 8 ABCM

Relative Stability of Mathieu Equation in Third Zone

Relative Stability of Mathieu Equation in Third Zone Relative tability of Mathieu Equation in Third Zone N. Mahmoudian, Ph.D. tudent, Nina.Mahmoudian@ndsu.nodak.edu Tel: 701.231.8303, Fax: 701.231-8913 G. Nakhaie Jazar, Assistant Professor, Reza.N.Jazar@ndsu.nodak.edu

More information

Melnikov s Method Applied to a Multi-DOF Ship Model

Melnikov s Method Applied to a Multi-DOF Ship Model Proceedings of the th International Ship Staility Workshop Melnikov s Method Applied to a Multi-DOF Ship Model Wan Wu, Leigh S. McCue, Department of Aerospace and Ocean Engineering, Virginia Polytechnic

More information

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)).

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)). Difference Equations to Differential Equations Section 8.5 Applications: Pendulums Mass-Spring Systems In this section we will investigate two applications of our work in Section 8.4. First, we will consider

More information

Effect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium

Effect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium IOSR Journal of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 39-765X. Volume, Issue Ver. III (Jan. - Fe. 06), 5-65 www.iosrjournals.org Effect of Uniform Horizontal Magnetic Field on Thermal Instaility

More information

A Stable and Convergent Finite Difference Scheme for 2D Incompressible Nonlinear Viscoelastic Fluid Dynamics Problem

A Stable and Convergent Finite Difference Scheme for 2D Incompressible Nonlinear Viscoelastic Fluid Dynamics Problem Applied and Computational Mathematics 2018; (1): 11-18 http://www.sciencepulishinggroup.com/j/acm doi: 10.11648/j.acm.2018001.12 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Stale and Convergent

More information

2. Determine whether the following pair of functions are linearly dependent, or linearly independent:

2. Determine whether the following pair of functions are linearly dependent, or linearly independent: Topics to be covered on the exam include: Recognizing, and verifying solutions to homogeneous second-order linear differential equations, and their corresponding Initial Value Problems Recognizing and

More information

ASEISMIC DESIGN OF TALL STRUCTURES USING VARIABLE FREQUENCY PENDULUM OSCILLATOR

ASEISMIC DESIGN OF TALL STRUCTURES USING VARIABLE FREQUENCY PENDULUM OSCILLATOR ASEISMIC DESIGN OF TALL STRUCTURES USING VARIABLE FREQUENCY PENDULUM OSCILLATOR M PRANESH And Ravi SINHA SUMMARY Tuned Mass Dampers (TMD) provide an effective technique for viration control of flexile

More information

2:2:1 Resonance in the Quasiperiodic Mathieu Equation

2:2:1 Resonance in the Quasiperiodic Mathieu Equation Nonlinear Dynamics 31: 367 374, 003. 003 Kluwer Academic Publishers. Printed in the Netherlands. ::1 Resonance in the Quasiperiodic Mathieu Equation RICHARD RAND Department of Theoretical and Applied Mechanics,

More information

Responses of two nonlinear microbial models to warming and increased carbon input

Responses of two nonlinear microbial models to warming and increased carbon input Biogeosciences, 13, 887 92, 216 www.iogeosciences.net/13/887/216/ doi:1.5194/g-13-887-216 Author(s) 216. CC Attriution 3. License. Responses of two nonlinear microial models to warming and increased caron

More information

SIMILARITY METHODS IN ELASTO-PLASTIC BEAM BENDING

SIMILARITY METHODS IN ELASTO-PLASTIC BEAM BENDING Similarity methods in elasto-plastic eam ending XIII International Conference on Computational Plasticity Fundamentals and Applications COMPLAS XIII E Oñate, DRJ Owen, D Peric and M Chiumenti (Eds) SIMILARIT

More information

Hopf bifurcation and Turing instability in the reaction diffusion Holling Tanner predator prey model

Hopf bifurcation and Turing instability in the reaction diffusion Holling Tanner predator prey model IMA Journal of Applied Mathematics 203 78, 287 306 doi:0.093/imamat/hxr050 Advance Access pulication on Novemer 7, 20 Hopf ifurcation and Turing instaility in the reaction diffusion Holling Tanner predator

More information

Robot Position from Wheel Odometry

Robot Position from Wheel Odometry Root Position from Wheel Odometry Christopher Marshall 26 Fe 2008 Astract This document develops equations of motion for root position as a function of the distance traveled y each wheel as a function

More information

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point

More information

CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD

CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen

More information

ROUNDOFF ERRORS; BACKWARD STABILITY

ROUNDOFF ERRORS; BACKWARD STABILITY SECTION.5 ROUNDOFF ERRORS; BACKWARD STABILITY ROUNDOFF ERROR -- error due to the finite representation (usually in floatingpoint form) of real (and complex) numers in digital computers. FLOATING-POINT

More information

VARIABLE COEFFICIENT OF FRICTION: AN EFFECTIVE VFPI PARAMETER TO CONTROL NEAR-FAULT GROUND MOTIONS

VARIABLE COEFFICIENT OF FRICTION: AN EFFECTIVE VFPI PARAMETER TO CONTROL NEAR-FAULT GROUND MOTIONS ISET Journal of Earthquake Technology, Paper No. 5, Vol. 49, No. 3-4, Sept.-Dec. 0, pp. 73 87 VARIABLE COEFFICIENT OF FRICTION: AN EFFECTIVE VFPI PARAMETER TO CONTROL NEAR-FAULT GROUND MOTIONS Girish Malu*

More information

MATH 251 Final Examination May 4, 2015 FORM A. Name: Student Number: Section:

MATH 251 Final Examination May 4, 2015 FORM A. Name: Student Number: Section: MATH 251 Final Examination May 4, 2015 FORM A Name: Student Number: Section: This exam has 16 questions for a total of 150 points. In order to obtain full credit for partial credit problems, all work must

More information

Design and Analysis of Silicon Resonant Accelerometer

Design and Analysis of Silicon Resonant Accelerometer Research Journal of Applied Sciences, Engineering and Technology 5(3): 970-974, 2013 ISSN: 2040-7459; E-ISSN: 2040-7467 Maxwell Scientific Organization, 2013 Sumitted: June 22, 2012 Accepted: July 23,

More information

ME 323 Examination #2

ME 323 Examination #2 ME 33 Eamination # SOUTION Novemer 14, 17 ROEM NO. 1 3 points ma. The cantilever eam D of the ending stiffness is sujected to a concentrated moment M at C. The eam is also supported y a roller at. Using

More information

4.2 Homogeneous Linear Equations

4.2 Homogeneous Linear Equations 4.2 Homogeneous Linear Equations Homogeneous Linear Equations with Constant Coefficients Consider the first-order linear differential equation with constant coefficients a 0 and b. If f(t) = 0 then this

More information

Order Reduction of Parametrically Excited Linear and Nonlinear Structural Systems

Order Reduction of Parametrically Excited Linear and Nonlinear Structural Systems Venkatesh Deshmukh 1 Eric A. Butcher e-mail: ffeab@uaf.ediu Department of Mechanical Engineering, University of Alaska Fairbanks, Fairbanks, AK 99775 S. C. Sinha Nonlinear Systems Research Laboratory,

More information

MULTIPLE POSITIVE SOLUTIONS FOR A KIRCHHOFF TYPE PROBLEM

MULTIPLE POSITIVE SOLUTIONS FOR A KIRCHHOFF TYPE PROBLEM Electronic Journal of Differential Equations, Vol. 205 (205, No. 29, pp. 0. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu MULTIPLE POSITIVE SOLUTIONS

More information

Construction of Manifolds of Positive Ricci Curvature with Big Volume and Large Betti Numbers

Construction of Manifolds of Positive Ricci Curvature with Big Volume and Large Betti Numbers Comparison Geometry MSRI Pulications Volume 30, 1997 Construction of Manifolds of Positive Ricci Curvature with Big Volume and Large Betti Numers G. PERELMAN Astract. It is shown that a connected sum of

More information

A NEW METHOD OF INVERSION OF THE LAPLACE TRANSFORM* ATHANASIOS PAPOULIS. Polytechnic Institute of Brooklyn

A NEW METHOD OF INVERSION OF THE LAPLACE TRANSFORM* ATHANASIOS PAPOULIS. Polytechnic Institute of Brooklyn 405 A NEW METHOD OF INVERSION OF THE LAPLACE TRANSFORM* BY ATHANASIOS PAPOULIS Polytechnic Institute of Brooklyn Introduction. In determining a function r(t) from its Laplace transform R(p) R(p) = [ e"'r(t)

More information

Vector Spaces. EXAMPLE: Let R n be the set of all n 1 matrices. x 1 x 2. x n

Vector Spaces. EXAMPLE: Let R n be the set of all n 1 matrices. x 1 x 2. x n Vector Spaces DEFINITION: A vector space is a nonempty set V of ojects, called vectors, on which are defined two operations, called addition and multiplication y scalars (real numers), suject to the following

More information

AN IMPROVED EFFECTIVE WIDTH METHOD BASED ON THE THEORY OF PLASTICITY

AN IMPROVED EFFECTIVE WIDTH METHOD BASED ON THE THEORY OF PLASTICITY Advanced Steel Construction Vol. 6, No., pp. 55-547 () 55 AN IMPROVED EFFECTIVE WIDTH METHOD BASED ON THE THEORY OF PLASTICITY Thomas Hansen *, Jesper Gath and M.P. Nielsen ALECTIA A/S, Teknikeryen 34,

More information

AP Calculus Testbank (Chapter 9) (Mr. Surowski)

AP Calculus Testbank (Chapter 9) (Mr. Surowski) AP Calculus Testbank (Chapter 9) (Mr. Surowski) Part I. Multiple-Choice Questions n 1 1. The series will converge, provided that n 1+p + n + 1 (A) p > 1 (B) p > 2 (C) p >.5 (D) p 0 2. The series

More information

Vibrations Qualifying Exam Study Material

Vibrations Qualifying Exam Study Material Vibrations Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering vibrations topics. These topics are listed below for clarification. Not all instructors

More information

1 Systems of Differential Equations

1 Systems of Differential Equations March, 20 7- Systems of Differential Equations Let U e an open suset of R n, I e an open interval in R and : I R n R n e a function from I R n to R n The equation ẋ = ft, x is called a first order ordinary

More information

Limit cycle oscillations at resonances

Limit cycle oscillations at resonances IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Limit cycle oscillations at resonances To cite this article: K Hellevik and O T Gudmestad 2017 IOP Conf. Ser.: Mater. Sci. Eng.

More information

DYNAMICAL MODELING AN INTERDISCIPLINARY APPROACH OF THE ECONOMIC REALITY

DYNAMICAL MODELING AN INTERDISCIPLINARY APPROACH OF THE ECONOMIC REALITY DYNAMICAL MODELING AN INTERDISCIPLINARY APPROACH OF THE ECONOMIC REALITY Assoc. Prof. Laura Ungureanu, Spiru Haret University Craiova, Romania Assoc. Prof. Ion Viorel Matei Spiru Haret University Craiova,

More information

Chapter 1. Harmonic Oscillator. 1.1 Energy Analysis

Chapter 1. Harmonic Oscillator. 1.1 Energy Analysis Chapter 1 Harmonic Oscillator Figure 1.1 illustrates the prototypical harmonic oscillator, the mass-spring system. A mass is attached to one end of a spring. The other end of the spring is attached to

More information

Composite Plates Under Concentrated Load on One Edge and Uniform Load on the Opposite Edge

Composite Plates Under Concentrated Load on One Edge and Uniform Load on the Opposite Edge Mechanics of Advanced Materials and Structures, 7:96, Copyright Taylor & Francis Group, LLC ISSN: 57-69 print / 57-65 online DOI:.8/5769955658 Composite Plates Under Concentrated Load on One Edge and Uniform

More information

Transduction Based on Changes in the Energy Stored in an Electrical Field

Transduction Based on Changes in the Energy Stored in an Electrical Field Lecture 6- Transduction Based on Changes in the Energy Stored in an Electrical Field Actuator Examples Microgrippers Normal force driving In-plane force driving» Comb-drive device F = εav d 1 ε oε F rwv

More information

arxiv: v1 [cs.fl] 24 Nov 2017

arxiv: v1 [cs.fl] 24 Nov 2017 (Biased) Majority Rule Cellular Automata Bernd Gärtner and Ahad N. Zehmakan Department of Computer Science, ETH Zurich arxiv:1711.1090v1 [cs.fl] 4 Nov 017 Astract Consider a graph G = (V, E) and a random

More information

Applications of Second-Order Differential Equations

Applications of Second-Order Differential Equations Applications of Second-Order Differential Equations ymy/013 Building Intuition Even though there are an infinite number of differential equations, they all share common characteristics that allow intuition

More information

Math Assignment 5

Math Assignment 5 Math 2280 - Assignment 5 Dylan Zwick Fall 2013 Section 3.4-1, 5, 18, 21 Section 3.5-1, 11, 23, 28, 35, 47, 56 Section 3.6-1, 2, 9, 17, 24 1 Section 3.4 - Mechanical Vibrations 3.4.1 - Determine the period

More information

The logistic difference equation and the route to chaotic behaviour

The logistic difference equation and the route to chaotic behaviour The logistic difference equation and the route to chaotic behaviour Level 1 module in Modelling course in population and evolutionary biology (701-1418-00) Module author: Sebastian Bonhoeffer Course director:

More information

The Sommerfeld Polynomial Method: Harmonic Oscillator Example

The Sommerfeld Polynomial Method: Harmonic Oscillator Example Chemistry 460 Fall 2017 Dr. Jean M. Standard October 2, 2017 The Sommerfeld Polynomial Method: Harmonic Oscillator Example Scaling the Harmonic Oscillator Equation Recall the basic definitions of the harmonic

More information

transformation of the reactants through highly reactive intermediates (active centers) to final products.

transformation of the reactants through highly reactive intermediates (active centers) to final products. LETUE. ETIN MEHNISMS ND EVLUTIN F TE FMS We have already seen that single reactions, in general, do not occur in one step; ut that it is a particular sequence o elementary reactions (reaction mechanism)

More information

Chaotic motion. Phys 750 Lecture 9

Chaotic motion. Phys 750 Lecture 9 Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to

More information

Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps

Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit

More information

Kink, singular soliton and periodic solutions to class of nonlinear equations

Kink, singular soliton and periodic solutions to class of nonlinear equations Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 10 Issue 1 (June 015 pp. 1 - Applications and Applied Mathematics: An International Journal (AAM Kink singular soliton and periodic

More information

Partial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators

Partial Differential Equations Separation of Variables. 1 Partial Differential Equations and Operators PDE-SEP-HEAT-1 Partial Differential Equations Separation of Variables 1 Partial Differential Equations and Operators et C = C(R 2 ) be the collection of infinitely differentiable functions from the plane

More information

arxiv: v1 [nlin.ao] 16 Jul 2017

arxiv: v1 [nlin.ao] 16 Jul 2017 arxiv:1707.05676v1 [nlin.ao] 16 Jul 017 Andronov Hopf ifurcation with and without parameter in a cuic memristor oscillator with a line of equiliria Ivan A. Korneev 1 1, a) and Vladimir V. Semenov Department

More information

Response of A Hard Duffing Oscillator to Harmonic Excitation An Overview

Response of A Hard Duffing Oscillator to Harmonic Excitation An Overview INDIN INSTITUTE OF TECHNOLOGY, KHRGPUR 710, DECEMBER 8-0, 00 1 Response of Hard Duffing Oscillator to Harmonic Excitation n Overview.K. Mallik Department of Mechanical Engineering Indian Institute of Technology

More information

STATIONARY AND TRANSIENT VIBRATIONS OF DAVI VIBROINSULATION SYSTEM.

STATIONARY AND TRANSIENT VIBRATIONS OF DAVI VIBROINSULATION SYSTEM. STTIONRY ND TRNSIENT VIBRTIONS OF DVI VIBROINSULTION SYSTEM. Jerzy Michalczyk, ul.zagorze c, 3 398 Cracow, Poland Piotr Czuak, ul. Szujskiego /4, 3-3 Cracow, Poland address for correspondence Grzegorz

More information

Copyright 2012 Tech Science Press CMES, vol.84, no.5, pp , 2012

Copyright 2012 Tech Science Press CMES, vol.84, no.5, pp , 2012 Copyright 01 Tech Science Press CMES, vol84, no5, pp459-497, 01 A Simple Collocation Scheme for Obtaining the Periodic Solutions of the Duffing Equation, and its Equivalence to the High Dimensional Harmonic

More information

Solution of Excited Non-Linear Oscillators under Damping Effects Using the Modified Differential Transform Method

Solution of Excited Non-Linear Oscillators under Damping Effects Using the Modified Differential Transform Method Article Solution of Excited Non-Linear Oscillators under Damping Effects Using the Modified Differential Transform Method H. M. Abdelhafez Department of Physics and Engineering Mathematics, Faculty of

More information

Tarek S. Ellaithy 1, Alexander Prexl 1, Nicolas Daynac 2. (1) DEA Egypt Branches, (2) Eliis SAS, France. Introduction

Tarek S. Ellaithy 1, Alexander Prexl 1, Nicolas Daynac 2. (1) DEA Egypt Branches, (2) Eliis SAS, France. Introduction Methodology and results of Nile Delta Messinian-Pliocene Holistic 3D seismic interpretation A case study from the Disouq area, Onshore Nile Delta, Egypt Tarek S. Ellaithy 1, Alexander Prexl 1, Nicolas

More information

status solidi Stability of elastic material with negative stiffness and negative Poisson s ratio

status solidi Stability of elastic material with negative stiffness and negative Poisson s ratio physica pss www.pss-.com status solidi asic solid state physics Editor s Choice Staility of elastic material with negative stiffness and negative Poisson s ratio Shang Xinchun 1, and Roderic S. Lakes,

More information

OPTIMIZATION OF A BRAYTON CYCLE AT MAXIMUM POWER OPERATING CONDITION AND AT MAXIMUM THERMAL EFFICIENCY OPERATING CONDITION

OPTIMIZATION OF A BRAYTON CYCLE AT MAXIMUM POWER OPERATING CONDITION AND AT MAXIMUM THERMAL EFFICIENCY OPERATING CONDITION SSN 76-580 nd nternational ongress of Mechanical Engineering Novemer -7, 0, Rieirão Preto, SP, Brazil opyright 0 y ABM OPMZAON OF A BRAYON YE A MAXMUM POER OPERANG ONDON AND A MAXMUM ERMA EFFENY OPERANG

More information

Dynamics of a mass-spring-pendulum system with vastly different frequencies

Dynamics of a mass-spring-pendulum system with vastly different frequencies Dynamics of a mass-spring-pendulum system with vastly different frequencies Hiba Sheheitli, hs497@cornell.edu Richard H. Rand, rhr2@cornell.edu Cornell University, Ithaca, NY, USA Abstract. We investigate

More information

Report E-Project Henriette Laabsch Toni Luhdo Steffen Mitzscherling Jens Paasche Thomas Pache

Report E-Project Henriette Laabsch Toni Luhdo Steffen Mitzscherling Jens Paasche Thomas Pache Potsdam, August 006 Report E-Project Henriette Laabsch 7685 Toni Luhdo 7589 Steffen Mitzscherling 7540 Jens Paasche 7575 Thomas Pache 754 Introduction From 7 th February till 3 rd March, we had our laboratory

More information

Upper Bounds for Stern s Diatomic Sequence and Related Sequences

Upper Bounds for Stern s Diatomic Sequence and Related Sequences Upper Bounds for Stern s Diatomic Sequence and Related Sequences Colin Defant Department of Mathematics University of Florida, U.S.A. cdefant@ufl.edu Sumitted: Jun 18, 01; Accepted: Oct, 016; Pulished:

More information

TOPOGRAPHY INFLUENCE ON THE LAKE EQUATIONS IN BOUNDED DOMAINS

TOPOGRAPHY INFLUENCE ON THE LAKE EQUATIONS IN BOUNDED DOMAINS TOPOGRAPHY INFLUENCE ON THE LAKE EQUATIONS IN BOUNDED DOMAINS CHRISTOPHE LACAVE, TOAN T. NGUYEN, AND BENOIT PAUSADER Astract. We investigate the influence of the topography on the lake equations which

More information

Design Principles of Seismic Isolation

Design Principles of Seismic Isolation Design Principles of Seismic Isolation 3 George C. Lee and Zach Liang Multidisciplinary Center for Earthquake Engineering Research, University at Buffalo, State University of New York USA 1. Introduction

More information

Weak bidders prefer first-price (sealed-bid) auctions. (This holds both ex-ante, and once the bidders have learned their types)

Weak bidders prefer first-price (sealed-bid) auctions. (This holds both ex-ante, and once the bidders have learned their types) Econ 805 Advanced Micro Theory I Dan Quint Fall 2007 Lecture 9 Oct 4 2007 Last week, we egan relaxing the assumptions of the symmetric independent private values model. We examined private-value auctions

More information

NONLINEAR SYSTEMS. Step response of a servo using an ideal relay. Stability of a nonlinear device using Liapunov s method

NONLINEAR SYSTEMS. Step response of a servo using an ideal relay. Stability of a nonlinear device using Liapunov s method NONLINEAR SYSTEMS Describing function of a nonlinear device ( output is the cube of its input) Stability of a system described by its state equation and assumed Liapunov function Step response of a servo

More information

On Universality of Blow-up Profile for L 2 critical nonlinear Schrödinger Equation

On Universality of Blow-up Profile for L 2 critical nonlinear Schrödinger Equation On Universality of Blow-up Profile for L critical nonlinear Schrödinger Equation Frank Merle,, Pierre Raphael Université de Cergy Pontoise Institut Universitaire de France Astract We consider finite time

More information

Notes for Expansions/Series and Differential Equations

Notes for Expansions/Series and Differential Equations Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated

More information

One dimensional Maps

One dimensional Maps Chapter 4 One dimensional Maps The ordinary differential equation studied in chapters 1-3 provide a close link to actual physical systems it is easy to believe these equations provide at least an approximate

More information

DRIVEN and COUPLED OSCILLATORS. I Parametric forcing The pendulum in 1:2 resonance Santiago de Compostela

DRIVEN and COUPLED OSCILLATORS. I Parametric forcing The pendulum in 1:2 resonance Santiago de Compostela DRIVEN and COUPLED OSCILLATORS I Parametric forcing The pendulum in 1:2 resonance Santiago de Compostela II Coupled oscillators Resonance tongues Huygens s synchronisation III Coupled cell system with

More information

IN this paper we study a discrete optimization problem. Constrained Shortest Link-Disjoint Paths Selection: A Network Programming Based Approach

IN this paper we study a discrete optimization problem. Constrained Shortest Link-Disjoint Paths Selection: A Network Programming Based Approach Constrained Shortest Link-Disjoint Paths Selection: A Network Programming Based Approach Ying Xiao, Student Memer, IEEE, Krishnaiyan Thulasiraman, Fellow, IEEE, and Guoliang Xue, Senior Memer, IEEE Astract

More information

2 discretized variales approach those of the original continuous variales. Such an assumption is valid when continuous variales are represented as oat

2 discretized variales approach those of the original continuous variales. Such an assumption is valid when continuous variales are represented as oat Chapter 1 CONSTRAINED GENETIC ALGORITHMS AND THEIR APPLICATIONS IN NONLINEAR CONSTRAINED OPTIMIZATION Benjamin W. Wah and Yi-Xin Chen Department of Electrical and Computer Engineering and the Coordinated

More information

The Fundamentals of Duct System Design

The Fundamentals of Duct System Design Design Advisory #1: CAS-DA1-2003 The Fundamentals of Duct System Design Wester's Dictionary defines "fundamentalism" as "a movement or point of view marked y a rigid adherence to fundamental or asic principles".

More information

Physics 235 Chapter 4. Chapter 4 Non-Linear Oscillations and Chaos

Physics 235 Chapter 4. Chapter 4 Non-Linear Oscillations and Chaos Chapter 4 Non-Linear Oscillations and Chaos Non-Linear Differential Equations Up to now we have considered differential equations with terms that are proportional to the acceleration, the velocity, and

More information

Practice Problems For Test 3

Practice Problems For Test 3 Practice Problems For Test 3 Power Series Preliminary Material. Find the interval of convergence of the following. Be sure to determine the convergence at the endpoints. (a) ( ) k (x ) k (x 3) k= k (b)

More information

Abelian Theorem for Generalized Mellin-Whittaker Transform

Abelian Theorem for Generalized Mellin-Whittaker Transform IOSR Journal of Mathematics (IOSRJM) ISSN: 2278-5728 Volume 1, Issue 4 (July-Aug 212), PP 19-23 Aelian Theorem for Generalized Mellin-Whittaker Transform A. S. Gudadhe 1, R. V. Kene 2 1, 2 (Head Department

More information

Chapter #6 EEE State Space Analysis and Controller Design

Chapter #6 EEE State Space Analysis and Controller Design hapter #6 EEE3-83 State Space nalsis and ontroller Design State feedack Linear Quadratic Regulator (LQR) Open and closed loop estimators Reduced Order Estimators (not assessed) racking Module Leader: Dr

More information

1.17 Triangle Numbers

1.17 Triangle Numbers .7 riangle Numers he n triangle numer is nn ( ). he first few are, 3,, 0, 5,, he difference etween each numer and e next goes up y each time. he formula nn ( ) gives e ( n ) triangle numer for n. he n

More information

Chaos and R-L diode Circuit

Chaos and R-L diode Circuit Chaos and R-L diode Circuit Rabia Aslam Chaudary Roll no: 2012-10-0011 LUMS School of Science and Engineering Thursday, December 20, 2010 1 Abstract In this experiment, we will use an R-L diode circuit

More information

physica status solidi REPRINT basic solid state physics Negative compressibility, negative Poisson s ratio, and stability

physica status solidi REPRINT basic solid state physics Negative compressibility, negative Poisson s ratio, and stability physica pss www.pss-.com asic solid state physics Negative compressiility, negative Poisson s ratio, and staility Rod Lakes 1 and K. W. Wojciechowski 1 University of Wisconsin, Madison, WI 53706, USA Institute

More information

DYNAMICS OF AN UMBILICAL CABLE FOR UNDERGROUND BORE-WELL APPLICATIONS

DYNAMICS OF AN UMBILICAL CABLE FOR UNDERGROUND BORE-WELL APPLICATIONS DYNAMICS OF AN UMBILICAL CABLE FOR UNDERGROUND BORE-WELL APPLICATIONS 1 KANNAN SANKAR & 2 JAYA KUMAR.V.K Engineering and Industrial Services, TATA Consultancy Services, Chennai, India E-mail: sankar.kannan@tcs.com,

More information

NONLINEAR VIBRATIONS

NONLINEAR VIBRATIONS NONLINEAR VIBRATIONS Prof. S. K. Dwivedy Mechanical Engineering Department Indian Institute of Technology Guwahati dwivedy@iitg.ernet.in Joint initiative of IITs and IISc Funded by MHRD Page 1 of 31 Module

More information

Parametric optimization of a neutral system with two delays and PD-controller

Parametric optimization of a neutral system with two delays and PD-controller 10.2478/acsc-2013-0008 Archives of Control Sciences Volume 23LIX, 2013 No. 2, pages 131 143 Parametric optimization of a neutral system with two elays an PD-controller JÓZEF DUDA In this paper a parametric

More information

1. Define the following terms (1 point each): alternative hypothesis

1. Define the following terms (1 point each): alternative hypothesis 1 1. Define the following terms (1 point each): alternative hypothesis One of three hypotheses indicating that the parameter is not zero; one states the parameter is not equal to zero, one states the parameter

More information

Nonlinear System Analysis

Nonlinear System Analysis Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control

More information

Air and Heat Flow through Large Vertical Openings

Air and Heat Flow through Large Vertical Openings Air and Heat Flow through Large Vertical Openings J L M Hensen *, J van der Maas #, A Roos * * Eindhoven University of Technology # Ecole Polytechnique Federale de Lausanne After a short description of

More information

Perturbation Theory of Dynamical Systems

Perturbation Theory of Dynamical Systems Perturbation Theory of Dynamical Systems arxiv:math.ho/0111178 v1 15 Nov 2001 Nils Berglund Department of Mathematics ETH Zürich 8092 Zürich Switzerland Lecture Notes Summer Semester 2001 Version: November

More information

Universit at Fachbereich Stuttgart Mathematik Preprint 2005/006

Universit at Fachbereich Stuttgart Mathematik Preprint 2005/006 Universität Stuttgart Fachereich Mathematik Spatially nondecaying solutions of 2D Navier-Stokes equation in a strip Sergey Zelik Preprint 25/6 Universität Stuttgart Fachereich Mathematik Spatially nondecaying

More information

SKATEBOARD: A HUMAN CONTROLLED NON-HOLONOMIC SYSTEM

SKATEBOARD: A HUMAN CONTROLLED NON-HOLONOMIC SYSTEM Proceedings of the ASME 205 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 205 August 2-5, 205, Boston, Massachusetts, USA DETC205-4752

More information

On an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University

On an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University On an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University njrose@math.ncsu.edu 1. INTRODUCTION. The classical eigenvalue problem for the Legendre Polynomials

More information

d mx kx dt ( t) Note first that imposed variation in the mass term is easily dealt with, by simply redefining the time dt mt ( ).

d mx kx dt ( t) Note first that imposed variation in the mass term is easily dealt with, by simply redefining the time dt mt ( ). . Parametric Resonance Michael Fowler Introduction (Following Landau para 7 A one-dimensional simple harmonic oscillator, a mass on a spring, d mx kx ( = has two parameters, m and k. For some systems,

More information

David A. Pape Department of Engineering and Technology Central Michigan University Mt Pleasant, Michigan

David A. Pape Department of Engineering and Technology Central Michigan University Mt Pleasant, Michigan Session: ENG 03-091 Deflection Solutions for Edge Stiffened Plates David A. Pape Department of Engineering and Technology Central Michigan University Mt Pleasant, Michigan david.pape@cmich.edu Angela J.

More information

p" N p" 2 p" 1 -ij" 1 =0 (W + ) ... q" 1 q" 2 q'' N =0 (W - )

p N p 2 p 1 -ij 1 =0 (W + ) ... q 1 q 2 q'' N =0 (W - ) Few-Body Systems 0, 6 (5) Few- Body Systems cfl y Springer-Verlag 5 PrintedinAustria Dynamical Hierarchy in Transition States of Reactions C.B. Li ; Λ, A. Shojiguchi ; ΛΛ,M.Toda ; ΛΛΛ,T. Komatsuzaki ;3;

More information

Improving convergence of incremental harmonic balance method using homotopy analysis method

Improving convergence of incremental harmonic balance method using homotopy analysis method Acta Mech Sin (2009) 25:707 712 DOI 10.1007/s10409-009-0256-4 RESEARCH PAPER Improving convergence of incremental harmonic balance method using homotopy analysis method Yanmao Chen Jike Liu Received: 10

More information

The Important State Coordinates of a Nonlinear System

The Important State Coordinates of a Nonlinear System The Important State Coordinates of a Nonlinear System Arthur J. Krener 1 University of California, Davis, CA and Naval Postgraduate School, Monterey, CA ajkrener@ucdavis.edu Summary. We offer an alternative

More information

Time Response of Systems

Time Response of Systems Chapter 0 Time Response of Systems 0. Some Standard Time Responses Let us try to get some impulse time responses just by inspection: Poles F (s) f(t) s-plane Time response p =0 s p =0,p 2 =0 s 2 t p =

More information

MA 266 Review Topics - Exam # 2 (updated)

MA 266 Review Topics - Exam # 2 (updated) MA 66 Reiew Topics - Exam # updated Spring First Order Differential Equations Separable, st Order Linear, Homogeneous, Exact Second Order Linear Homogeneous with Equations Constant Coefficients The differential

More information

Integer Sequences and Circle Chains Inside a Circular Segment

Integer Sequences and Circle Chains Inside a Circular Segment Forum Geometricorum Volume 18 018) 47 55. FORUM GEOM ISSN 1534-1178 Integer Sequences and Circle Chains Inside a Circular Segment Giovanni Lucca Astract. We derive the conditions for inscriing, inside

More information

Fourier Series. Green - underdamped

Fourier Series. Green - underdamped Harmonic Oscillator Fourier Series Undamped: Green - underdamped Overdamped: Critical: Underdamped: Driven: Calculus of Variations b f {y, y'; x}dx is stationary when f y d f = 0 dx y' a Note that y is

More information

2.7 - Materials selection and shape. Outline. The shape factor, and shape limits. Material indices that include shape

2.7 - Materials selection and shape. Outline. The shape factor, and shape limits. Material indices that include shape - shy Method.7 - Materials selection and shape Outline Shape efficiency The shape factor, and shape limits Material indices that include shape Graphical ways of dealing with shape Resources: M.. shy, Materials

More information

Topic 1: Plotting Data Basic introduction to entering and plotting data, displaying error bars and appropriate labels.

Topic 1: Plotting Data Basic introduction to entering and plotting data, displaying error bars and appropriate labels. PHY 530 Module C1 Spread Sheet Physics Topic 1: Plotting Data Basic introduction to entering and plotting data, displaying error bars and appropriate labels. Topic 2: Plotting functions How to generate

More information

Alabama Course of Study: Mathematics Precalculus

Alabama Course of Study: Mathematics Precalculus A Correlation of : Graphical, Numerical, Algebraic 8 th Edition to the INTRODUCTION This document demonstrates how : Graphical, Numerical, Algebraic, 8th Edition 2011, (Demana, et al.), meets the indicators

More information

Francisco Paulo Lépore Neto. Marcelo Braga dos Santos. Introduction 1. Nomenclature. Experimental Apparatus and Formulation

Francisco Paulo Lépore Neto. Marcelo Braga dos Santos. Introduction 1. Nomenclature. Experimental Apparatus and Formulation Francisco Paulo Lépore Neto and Marcelo Braga dos Santos Francisco Paulo Lépore Neto fplepore@mecanica.ufu.br Federal University of Uberlandia School of Mechanical Engineering 38408-902 Uberlandia, MG,

More information

Only this exam and a pen or pencil should be on your desk.

Only this exam and a pen or pencil should be on your desk. Lin. Alg. & Diff. Eq., Spring 16 Student ID Circle your section: 31 MWF 8-9A 11 LATIMER LIANG 33 MWF 9-1A 11 LATIMER SHAPIRO 36 MWF 1-11A 37 CORY SHAPIRO 37 MWF 11-1P 736 EVANS WORMLEIGHTON 39 MWF -5P

More information

On the principal eigenvalue of elliptic operators in R N and applications

On the principal eigenvalue of elliptic operators in R N and applications On the principal eigenvalue of elliptic operators in R N and applications Henri Berestycki a and Luca Rossi a, a EHESS, CAMS, 54 Boulevard Raspail, F-75006 Paris, France Università La Sapienza Roma I,

More information

19 Tree Construction using Singular Value Decomposition

19 Tree Construction using Singular Value Decomposition 19 Tree Construction using Singular Value Decomposition Nicholas Eriksson We present a new, statistically consistent algorithm for phylogenetic tree construction that uses the algeraic theory of statistical

More information

Notes on Special Functions

Notes on Special Functions Spring 25 1 Notes on Special Functions Francis J. Narcowich Department of Mathematics Texas A&M University College Station, TX 77843-3368 Introduction These notes are for our classes on special functions.

More information