Showing papers in "Nonlinear Dynamics in 2005"

Spectral Properties of Dynamical Systems, Model Reduction and Decompositions

Abstract: In this paper we discuss two issues related to model reduction of deterministic or stochastic processes. The first is the relationship of the spectral properties of the dynamics on the attractor of the original, high-dimensional dynamical system with the properties and possibilities for model reduction. We review some elements of the spectral theory of dynamical systems. We apply this theory to obtain a decomposition of the process that utilizes spectral properties of the linear Koopman operator associated with the asymptotic dynamics on the attractor. This allows us to extract the almost periodic part of the evolving process. The remainder of the process has continuous spectrum. The second topic we discuss is that of model validation, where the original, possibly high-dimensional dynamics and the dynamics of the reduced model – that can be deterministic or stochastic – are compared in some norm. Using the “statistical Takens theorem” proven in (Mezic, I. and Banaszuk, A. Physica D, 2004) we argue that comparison of average energy contained in the finite-dimensional projection is one in the hierarchy of functionals of the field that need to be checked in order to assess the accuracy of the projection.

The Method of Proper Orthogonal Decomposition for Dynamical Characterization and Order Reduction of Mechanical Systems: An Overview

Abstract: Modal analysis is used extensively for understanding the dynamic behavior of structures. However, a major concern for structural dynamicists is that its validity is limited to linear structures. New developments have been proposed in order to examine nonlinear systems, among which the theory based on nonlinear normal modes is indubitably the most appealing. In this paper, a different approach is adopted, and proper orthogonal decomposition is considered. The modes extracted from the decomposition may serve two purposes, namely order reduction by projecting high-dimensional data into a lower-dimensional space and feature extraction by revealing relevant but unexpected structure hidden in the data. The utility of the method for dynamic characterization and order reduction of linear and nonlinear mechanical systems is demonstrated in this study.

Reduced-Order Models for MEMS Applications

Abstract: We review the development of reduced-order models for MEMS devices. Based on their implementation procedures, we classify these reduced-order models into two broad categories: node and domain methods. Node methods use lower-order approximations of the system matrices found by evaluating the system equations at each node in the discretization mesh. Domain-based methods rely on modal analysis and the Galerkin method to rewrite the system equations in terms of domain-wide modes (eigenfunctions). We summarize the major contributions in the field and discuss the advantages and disadvantages of each implementation. We then present reduced-order models for microbeams and rectangular and circular microplates. Finally, we present reduced-order approaches to model squeeze-film and thermoelastic damping in MEMS and present analytical expressions for the damping coefficients. We validate these models by comparing their results with available theoretical and experimental results.

Low-Dimensional Modelling of Turbulence Using the Proper Orthogonal Decomposition: A Tutorial

Abstract: The proper orthogonal decomposition identifies basis functions or modes which optimally capture the average energy content from numerical or experimental data. By projecting the Navier–Stokes equations onto these modes and truncating, one can obtain low-dimensional ordinary differential equation models for fluid flows. In this paper we present a tutorial on the construction of such models. In addition to providing a general overview of the procedure, we describe two different ways to numerically calculate the modes, show how symmetry considerations can be exploited to simplify and understand them, comment on how parameter variations are captured naturally in such models, and describe a generalization of the procedure involving projection onto uncoupled modes that allow streamwise and cross-stream components to evolve independently. We illustrate for the example of plane Couette flow in a minimal flow unit – a domain whose spanwise and streamwise extent is just sufficient to maintain turbulence.

Identification of Nonlinear Aeroelastic Systems Based on the Volterra Theory: Progress and Opportunities

Abstract: The identification of nonlinear aeroelastic systems based on the Volterra theory of nonlinear systems is presented. Recent applications of the theory to problems in computational and experimental aeroelasticity are reviewed. Computational results include the development of computationally efficient reduced-order models (ROMs) using an Euler/Navier–Stokes flow solver and the analytical derivation of Volterra kernels for a nonlinear aeroelastic system. Experimental results include the identification of aerodynamic impulse responses, the application of higher-order spectra (HOS) to wind-tunnel flutter data, and the identification of nonlinear aeroelastic phenomena from flight flutter test data of the active aeroelastic wing (AAW) aircraft.

Homotopy solutions for a generalized second-grade fluid past a porous plate

Abstract: The flow of a second-grade fluid past a porous plate subject to either suction or blowing at the plate has been studied. A modified model of second-grade fluid that has shear-dependent viscosity and can predict the normal stress difference is used. The differential equations governing the flow are solved using homotopy analysis method (HAM). Expressions for the velocity have been constructed and discussed with the help of graphs. Analysis of the obtained results showed that the flow is appreciably influenced by the material and normal stress coefficient. Several results of interest are deduced as the particular cases of the presented analysis.

On the Hysteretic Bouc–Wen Model

Abstract: This paper deals with the problem of identifying the parameters of the hysteretic Bouc–Wen model. In the existing literature, the methods devoted to this problem rely mainly on numerical simulations and do not have, to a very large extent, a rigorous mathematical justification. Our method consists in exciting the hysteretic system with a periodic input and obtain the desired parameters from the resulting limit cycle. The identification method that we propose has a rigorous mathematical basis as it based on the analytic description of the limit cycle, and, unlike existing identification methods for the Bouc–Wen model, gives guaranteed relative errors between the (unknown) exact model parameters and their corresponding estimates. We also prove that this method is robust with respect to constant and T-periodic disturbances commonly present in any laboratory experiment. A numerical simulation example illustrates the use of our identification method.

Dimension Reduction of Dynamical Systems: Methods, Models, Applications

Abstract: After presenting some basic introductory ideas concerning dimension reduction and reduced order modelling, an overview of the contents of the papers collected in this Special Issue of Nonlinear Dynamics is given.

Peristaltic Flow of a Magnetohydrodynamic Johnson–Segalman Fluid

Abstract: The problem of peristaltic transport of non-Newtonian fluid represented by the constitutive equation for a Johnson–Segalman fluid is analyzed for the case of a planar channel. The fluid is electrically conducting. The walls of the channel are electrically insulated and are transversely displaced by an infinite, harmonic travelling wave of long wavelength. The general solution of the non-linear equation resulting from the momentum equation is constructed for all values of Weissenberg number. The perturbation solution is also obtained. Some graphs are plotted for interesting physical parameters and discussed.

Resolving Controversies in the Application of the Method of Multiple Scales and the Generalized Method of Averaging

Abstract: I compare application of the method of multiple scales with reconstitution and the generalized method of averaging for determining higher-order approximations of three single-degree-of-freedom systems and a two-degree-of-freedom system. Three implementations of the method of multiple scales are considered, namely, application of the method to the system equations expressed as second-order equations, as first-order equations, and in complex-variable form. I show that all of these methods produce the same modulation equations. I address the problem of determining higher-order approximate solutions of the Duffing equation in the case of primary resonance. I show that the conclusions of Rahman and Burton that the method of multiple scales, the generalized method of averaging, and Lie series and transforms might lead to incorrect results, in that spurious solutions occur and the obtained frequency–response curves bear little resemblance to the actual response, is the result of their using parameter values for which the neglected terms are the same order as the retained terms. I show also that spurious solutions cannot be avoided, in general, in any consistent expansion and their presence does not constitute a limitation of the methods. In particular, I show that, for the Duffing equation, the second-order frequency–response equation does not possess spurious solutions for the case of hardening nonlinearity, but possesses spurious solutions for the case of softening nonlinearity. For sufficiently small nonlinearity, the spurious solutions are far removed from the actual response. But as the strength of the nonlinearity increases, these solutions move closer to the backbone and eventually distort it. This is not a drawback of the perturbation methods but an indication of an application of the analysis for parameter values outside the range of validity of the expansion. Also, I address the problem of obtaining non-Hamiltonian modulation equations in the application of the method of multiple scales to multi-degree-of-freedom Hamiltonian systems written as second-order equations in time and how this problem can be overcome by attacking the state-space form of the governing equations. Moreover, I show that application of a variation of the method of Rahman and Burton to multi-degree-of-freedom systems leads to results that do not agree with those obtained with the generalized method of averaging.

On the Computer Formulations of the Wheel/Rail Contact Problem

Abstract: In this investigation, four nonlinear dynamic formulations that can be used in the analysis of the wheel/rail contact are presented, compared and their performance is evaluated. Two of these formulations employ nonlinear algebraic kinematic constraint equations to describe the contact between the wheel and the rail (constraint approach), while in the other two formulations the contact force is modeled using a compliant force element (elastic approach). The goal of the four formulations is to provide accurate nonlinear modeling of the contact between the wheel and the rail, which is crucial to the success of any computational algorithm used in the dynamic analysis of railroad vehicle systems. In the formulations based on the elastic approach, the wheel has six degrees of freedom with respect to the rail, and the normal contact forces are defined as function of the penetration using Hertz’s contact theory or using assumed stiffness and damping coefficients. The first elastic method is based on a search for the contact locations using discrete nodal points. As previously presented in the literature, this method can lead to impulsive forces due to the abrupt change in the location of the contact point from one time step to the next. This difficulty is avoided in the second elastic approach in which the contact points are determined by solving a set of algebraic equations. In the formulations based on the constraint approach, on the other hand, the case of a non-conformal contact is assumed, and nonlinear kinematic contact constraint equations are used to impose the contact conditions at the position, velocity and acceleration levels. This approach leads to a model, in which the wheel has five degrees of freedom with respect to the rail. In the constraint approach, the wheel penetration and lift are not permitted, and the normal contact forces are calculated using the technique of Lagrange multipliers and the augmented form of the system dynamic equations. Two equivalent constraint formulations that employ two different solution procedures are discussed in this investigation. The first method leads to a larger system of equations by augmenting all the contact constraint equations to the dynamic equations of motion, while in the second method an embedding procedure is used to obtain a reduced system of equations from which the surface parameter accelerations are systematically eliminated. Numerical results are presented in order to examine the performance of various methods discussed in this study.

Global Bifurcations and Chaotic Dynamics in Nonlinear Nonplanar Oscillations of a Parametrically Excited Cantilever Beam

Abstract: This paper presents the analysis of the global bifurcations and chaotic dynamics for the nonlinear nonplanar oscillations of a cantilever beam subjected to a harmonic axial excitation and transverse excitations at the free end. The governing nonlinear equations of nonplanar motion with parametric and external excitations are obtained. The Galerkin procedure is applied to the partial differential governing equation to obtain a two-degree-of-freedom nonlinear system with parametric and forcing excitations. The resonant case considered here is 2:1 internal resonance, principal parametric resonance-1/2 subharmonic resonance for the in-plane mode and fundamental parametric resonance–primary resonance for the out-of-plane mode. The parametrically and externally excited system is transformed to the averaged equations by using the method of multiple scales. From the averaged equation obtained here, the theory of normal form is applied to find the explicit formulas of normal forms associated with a double zero and a pair of pure imaginary eigenvalues. Based on the normal form obtained above, a global perturbation method is utilized to analyze the global bifurcations and chaotic dynamics in the nonlinear nonplanar oscillations of the cantilever beam. The global bifurcation analysis indicates that there exist the heteroclinic bifurcations and the Silnikov type single-pulse homoclinic orbit in the averaged equation for the nonlinear nonplanar oscillations of the cantilever beam. These results show that the chaotic motions can occur in the nonlinear nonplanar oscillations of the cantilever beam. Numerical simulations verify the analytical predictions.

Lyapunov-Based Controller for the Inverted Pendulum Cart System

Abstract: A nonlinear control force is presented to stabilize the under-actuated inverted pendulum mounted on a cart. The control strategy is based on partial feedback linearization, in a first stage, to linearize only the actuated coordinate of the inverted pendulum, and then, a suitable Lyapunov function is formed to obtain a stabilizing feedback controller. The obtained closed-loop system is locally asymptotically stable around its unstable equilibrium point. Additionally, it has a very large attraction domain.

Periodic and Chaotic Motions of a Rotor-Active Magnetic Bearing with Quadratic and Cubic Terms and Time-Varying Stiffness

Abstract: In this paper, we use the asymptotic perturbation method to investigate nonlinear oscillations and chaotic dynamics in a rotor-active magnetic bearings (AMB) system with 8-pole legs and the time-varying stiffness. The stiffness in the AMB is considered as the time varying in a periodic form. Because of considering the weight of the rotor, the formulation on the electromagnetic force resultants includes the quadratic and cubic nonlinearities. The resulting dimensionless equations of motion for the rotor-AMB system with the time-varying stiffness in the horizontal and vertical directions are a two-degree-of-freedom nonlinear system with quadratic and cubic nonlinearities and parametric excitation. The asymptotic perturbation method is used to obtain the averaged equations in the case of primary parametric resonance and 1/2 subharmonic resonance. It is found that there exist period-3, period-4, period-6, period-7, period-8, quasiperiodic and chaotic modulated amplitude oscillations in the rotor-AMB system with the time-varying stiffness. It is seen from the numerical results that there are the phenomena of the multiple solutions and the soft-spring type and the hardening-spring type in nonlinear frequency-response curves for the rotor-AMB system. The parametric excitation, or the time-varying stiffness produced by the PD controller is considered to be a controlling force which can control the chaotic response in the rotor-AMB system to a period n motion.

Ši’lnikov Chaos in the Generalized Lorenz Canonical Form of Dynamical Systems

Abstract: This paper studies the generalized Lorenz canonical form of dynamical systems introduced by Celikovský and Chen [International Journal of Bifurcation and Chaos12(8), 2002, 1789]. It proves the existence of a heteroclinic orbit of the canonical form and the convergence of the corresponding series expansion. The Si’lnikov criterion along with some technical conditions guarantee that the canonical form has Smale horseshoes and horseshoe chaos. As a consequence, it also proves that both the classical Lorenz system and the Chen system have Si’lnikov chaos. When the system is changed into another ordinary differential equation through a nonsingular one-parameter linear transformation, the exact range of existence of Si’lnikov chaos with respect to the parameter can be specified. Numerical simulation verifies the theoretical results and analysis.

An Internal Damping Model for the Absolute Nodal Coordinate Formulation

Abstract: Introducing internal damping in multibody system simulations is important as real-life systems usually exhibit this type of energy dissipation mechanism. When using an inertial coordinate method such as the absolute nodal coordinate formulation, damping forces must be carefully formulated in order not to damp rigid body motion, as both this and deformation are described by the same set of absolute nodal coordinates. This paper presents an internal damping model based on linear viscoelasticity for the absolute nodal coordinate formulation. A practical procedure for estimating the parameters that govern the dissipation of energy is proposed. The absence of energy dissipation under rigid body motion is demonstrated both analytically and numerically. Geometric nonlinearity is accounted for as deformations and deformation rates are evaluated by using the Green–Lagrange strain–displacement relationship. In addition, the resulting damping forces are functions of some constant matrices that can be calculated in advance, thereby avoiding the integration over the element volume each time the damping force vector is evaluated.

Fuzzy Logic Control of Seat Vibrations of a Non-Linear Full Vehicle Model

Abstract: In this paper, the dynamic behavior of a non-linear eight degrees of freedom vehicle model having active suspensions and a fuzzy logic (FL) controlled passenger seat is examined. The non-linearity occurs due to dry friction on the dampers. Three cases of control strategies are taken into account. In the first case, only the passenger seat is controlled. In the second case, only the vehicle body is controlled. In the third case, both the vehicle body and the passenger seat are fully controlled at the same time. The time responses of the non-linear vehicle model due to road disturbance and the frequency responses are obtained for each control strategy. At the end, the performances of these strategies are compared.

Inverse Dynamics of Servo-Constraints Based on the Generalized Inverse

Abstract: The acceleration form of constraint equations is utilized in this paper to solve for the inverse dynamics of servo-constraints. A condition for the existence of control forces that enforce servo-constraints is derived. For overactuated dynamical systems, the generalized Moore-Penrose inverse of the constraint matrix is used to parameterize the solutions for these control forces in terms of free parameters that can be chosen to satisfy certain requirements or optimize certain criterions. In particular, these free parameters can be chosen to minimize the Gibbsian (i.e., the acceleration energy of the dynamical system), resulting in “ideal” control forces (those satisfying the principle of virtual work when the virtual displacements satisfy the servo-constraint equations). To achieve this, the nonminimal nonholonomic form recently derived by the authors in the context of Kane’s method is used to determine the accelerations of the system, and hence to determine the forces to be generated by the redundant manipulators. Finally, an extension to inverse dynamics of servo-constraints involving control variables is made. The procedures are illustrated by two examples.

Three-dimensional large deformation analysis of the multibody pantograph/catenary systems

Abstract: To accurately model the nonlinear behavior of the pantograph/catenary systems, it is necessary to take into consideration the effect of the large deformation of the catenary and its interaction with the nonlinear pantograph system dynamics. The large deformation of the catenary is modeled in this investigation using the three-dimensional finite element absolute nodal coordinate formulation. To model the interaction between the pantograph and the catenary, a sliding joint that allows for the motion of the pan-head on the catenary cable is formulated. To this end, a non-generalized arc-length parameter is introduced in order to be able to accurately predict the location of the point of contact between the pan-head and the catenary. The resulting system of differential and algebraic equations formulated in terms of reference coordinates, finite element absolute nodal coordinates, and non-generalized arc-length and contact surface parameters are solved using computational multibody system algorithms. A detailed three-dimensional multibody railroad vehicle model is developed to demonstrate the use of the formulation presented in this paper. In this model, the interaction between the wheel and the rail is considered. For future research, a method is proposed to deal with the problem of the loss of contact between the pan-head and the catenary cable.

Energy Transfers in a System of Two Coupled Oscillators with Essential Nonlinearity: 1:1 Resonance Manifold and Transient Bridging Orbits

Abstract: The purpose of this study is to highlight and explain the vigorous energy transfers that may take place in a linear oscillator weakly coupled to an essentially nonlinear attachment, termed a nonlinear energy sink. Although these energy exchanges are encountered during the transient dynamics of the damped system, it is shown that the dynamics can be interpreted mainly in terms of the periodic orbits of the underlying Hamiltonian system. To this end, a frequency-energy plot gathering the periodic orbits of the system is constructed which demonstrates that, thanks to a 1:1 resonance capture, energy can be irreversibly and almost completely transferred from the linear oscillator to the nonlinear attachment. Furthermore, it is observed that this nonlinear energy pumping is triggered by the excitation of transient bridging orbits compatible with the nonlinear attachment being initially at rest, a common feature in most practical applications. A parametric study of the energy exchanges is also performed to understand the influence of the parameters of the nonlinear energy sink. Finally, the results of experimental measurements supporting the theoretical developments are discussed.

Dynamics of a Gear System with Faults in Meshing Stiffness

Abstract: Gear box dynamics is characterised by a periodically changing stiffness. In real gear systems, a backlash also exists that can lead to a loss in contact between the teeth. Due to this loss of contact the gear has piecewise linear stiffness characteristics, and the gears can vibrate regularly and chaotically. In this paper we examine the effect of tooth shape imperfections and defects. Using standard methods for nonlinear systems we examine the dynamics of gear systems with various faults in meshing stiffness.

Milling Bifurcations from Structural Asymmetry and Nonlinear Regeneration

Abstract: This paper investigates multiple modeling choices for analyzing the rich and complex dynamics of high-speed milling processes. Various models are introduced to capture the effects of asymmetric structural modes and the influence of nonlinear regeneration in a discontinuous cutting force model. Stability is determined from the development of a dynamic map for the resulting variational system. The general case of asymmetric structural elements is investigated with a fixed frame and rotating frame model to show differences in the predicted unstable regions due to parametric excitation. Analytical and numerical investigations are confirmed through a series of experimental cutting tests. The principal results are additional unstable regions, hysteresis in the bifurcation diagrams, and the presence of coexisting periodic and quasiperiodic attractors which is confirmed through experimentation.

Energy Pumping with Various Nonlinear Structures: Numerical Evidences

Abstract: Numerical investigations are carried out on a linear structure, weakly coupled to a small nonlinear attachment. The essential nonlinearity of the attachment enables it to resonate with any of the linearized modes of the structure leading to energy pumping, i.e. passive, one-way, irreversible transfer of energy from the structure to the attachment. Different nonlinear structures (piecewise linear system, chaotic system) and efficiency of energy pumping are studied in each case in order to be able to apply it to civil engineering. As a specific application, attenuation of vibrations of a building is studied with two building models. In particular, the case of stochastic excitations is analyzed to examine if it is possible to process energy pumping when a seism occurs and an indicator of efficiency has been introduced.

Oblique Frictional Impact of a Bar: Analysis and Comparison of Different Impact Laws

Abstract: In this paper a basic, easily to multi-contact problems extendable, non-smooth approach is applied to analyze a bar striking an inelastic half-space. Coulomb contact is assumed and modeled by using set-valued Newtonian impact laws in normal as well as in tangential direction. The resulting linear complementarity problem contains all possible impact states and provides an instantaneous collision operator that respects all inequality constraints. This operator depends on the orientation of the bar and determines uniquely the post-impact velocities as functions of the pre-impact state. Different types of solutions may occur, including “stick’’ and “slip’’. In this context, stick and slip have to be understood as the two cases characterized by the tangential impulsive force as an element of either the set-valued or of the single-valued domain of the friction law. Depending on the choice of parameters, sign reversal of the tangential contact velocity is possible. For certain inertia properties and initial conditions, the collision operator yields an impact, even for initially vanishing normal contact velocity. This phenomenon is well known as the Painleve paradox. The results obtained by this fully non-smooth rigid body approach are compared with those of other impact models, such as a lumped mass model with compliance elements, and a collision operator used for particle interactions in flows.

2:1:1 Resonance in the Quasi-Periodic Mathieu Equation

Abstract: We present a small e perturbation analysis of the quasi-periodic Mathieu equation $$ \ddot x + (\delta + \epsilon \cos t + \epsilon \cos \omega t) x=0 $$ in the neighborhood of the point δ = 0.25 and ω = 0.5. We use multiple scales including terms of O(e2) with three time scales. We obtain an asymptotic expansion for an associated instability region. Comparison with numerical integration shows good agreement for e = 0.1. Then we use the algebraic form of the perturbation solution to approximate scaling factors which are conjectured to determine the size of instability regions as we go from one resonance to another in the δ−ω parameter plane.

Component mode synthesis using nonlinear normal modes

Abstract: This paper describes a methodology for developing reduced-order dynamic models of structural systems that are composed of an assembly of nonlinear component structures. The approach is a nonlinear extension of the fixed-interface component mode synthesis (CMS) technique developed for linear structures by Hurty and modified by Craig and Bampton. Specifically, the case of nonlinear substructures is handled by using fixed-interface nonlinear normal modes (NNMs). These normal modes are constructed for the various substructures using an invariant manifold approach, and are then coupled through the traditional linear constraint modes (i.e., the static deformation shapes produced by unit interface displacements). A class of systems is used to demonstrate the concept and show the effectiveness of the proposed procedure. Simulation results show that the reduced-order model (ROM) obtained from the proposed procedure outperforms the ROM obtained from the classical fixed-interface linear CMS approach as applied to a nonlinear structure. The proposed method is readily applicable to large-scale nonlinear structural systems that are based on finite-element models.

Regenerative Tool Chatter Near a Codimension 2 Hopf Point Using Multiple Scales

Abstract: We study a well-known regenerative machine tool vibration model (a delay differential equation) near a codimension 2 Hopf bifurcation point. The method of multiple scales is used directly, bypassing a center manifold reduction. We use a nonstandard choice of expansion parameter that helps understand practically relevant aspects of the dynamics for not-too-small amplitudes. Analytical expressions are then obtained for the double Hopf points. Both sub- and supercritical bifurcations are predicted to occur near the reference point; and analytical conditions on the parameter variations for each type of bifurcation to occur are obtained as well. Analytical approximations are supported by numerics.

Continuous-Time Bilinear System Identification

Abstract: The objective of this paper is to describe a new method for identification of a continuous-time multi-input and multi-output bilinear system. The approach is to make judicious use of the linear-model properties of the bilinear system when subjected to a constant input. Two steps are required in the identification process. The first step is to use a set of pulse responses resulting from a constant input of one sample period to identify the state matrix, the output matrix, and the direct transmission matrix. The second step is to use another set of pulse responses with the same constant input over multiple sample periods to identify the input matrix and the coefficient matrices associated with the coupling terms between the state and the inputs. Numerical examples are given to illustrate the concept and the computational algorithm for the identification method.