The construction of non-linear normal modes for systems with internal resonance

A numerical method, based on the invariant manifold approach, is presented for constructing non-linear normal modes for systems with internal resonances. In order to parameterize the non-linear normal modes of interest, multiple pairs of system state variables involved in the internal resonance are kept as ‘seeds’ for the construction of the multi-mode invariant manifold. All the remaining degrees of freedom are then constrained to these ‘seed’, or master, variables, resulting in a system of non-linear partial differential equations that govern the constraint relationships, and these are solved numerically. The computationally-intensive solution procedure uses a combination of finite difference schemes and Galerkin-based expansion approaches. It is illustrated using two examples, both of which focus on the construction of two-mode models. The first example is based on the analysis of a simple three-degree-of-freedom example system, and is used to demonstrate the approach. An invariant manifold that captures two non-linear normal modes is constructed, resulting in a reduced order model that accurately captures the system dynamics. The methodology is then applied to a larger order system, specifically, an 18-degree-of-freedom rotating beam model that features a three-to-one internal resonance between the first two flapping modes. The accuracy of the non-linear two-mode reduced order model is verified by comparing time-domain simulations of the two DOF model and the full system equations of motion.     

Normal modes for piecewise linear vibratory systems

A method to construct the normal modes for a class of piecewise linear vibratory systems is developed in this study. The approach utilizes the concepts of Poincaré maps and invariant manifolds from the theory of dynamical systems. In contrast to conventional methods for smooth systems, which expand normal modes in a series form around an equilibrium point of interest, the present method expands the normal modes in a series form of polar coordinates in a neighborhood of an invariant disk of the system. It is found that the normal modes, modal dynamics and frequency-amplitude dependence relationship are all of piecewise type. A two degree of freedom example is used to demonstrate the method.     

Nonlinear Dynamic Analysis of MEMS Switches by Nonlinear Modal Analysis

A nonlinear modal analysis approach based on the invariant manifoldmethod proposed earlier by Boivin et al. [10] is applied in this paperto perform the dynamic analysis of a micro switch. The micro switch ismodeled as a clamped-clamped microbeam subjected to a transverseelectrostatic force. Two kinds of nonlinearities are encountered in thenonlinear system: geometric nonlinearity of the microbeam associatedwith large deflection, and nonlinear coupling between two energydomains. Using Galerkin method, the nonlinear partial differentialgoverning equation is decoupled into a set of nonlinear ordinarydifferential equations. Based on the invariant manifold method, theassociated nonlinear modal shapes, and modal motion governing equationsare obtained. The equation of motion restricted to these manifolds,which provide the dynamics of the associated normal modes, are solved bythe approach of nonlinear normal forms. Nonlinearities and the pull-inphenomena are examined. The numerical results are compared with thoseobtained from the finite difference method. The estimate for the pull-involtage of the micro device is also presented.     

Advanced Proper Orthogonal Decomposition Tools: Using Reduced Order Models to Identify Normal Modes of Vibration and Slow Invariant Manifolds in the Dynamics of Planar Nonlinear Rods

Reduced order models for the dynamics of geometrically exact planar rods are derived by projecting the nonlinear equations of motion onto a subspace spanned by a set of proper orthogonal modes. These optimal modes are identified by a proper orthogonal decomposition processing of high-resolution finite element dynamics. A three-degree-of-freedom reduced system is derived to study distinct categories of motions dominated by a single POD mode. The modal analysis of the reduced system characterizes in a unique fashion for these motions, since its linear natural frequencies are near to the natural frequencies of the full-order system. For free motions characterized by a single POD mode, the eigen-vector matrix of the derived reduced system coincides with the principal POD-directions. This property reflects the existence of a normal mode of vibration, which appears to be close to a slow invariant manifold. Its shape is captured by that of the dominant POD mode. The modal analysis of the POD-based reduced order system provides a potentially valuable tool to characterize the spatio-temporal complexity of the dynamics in order to elucidate connections between proper orthogonal modes and nonlinear normal modes of vibration.     

Damping for large-amplitude vibrations of plates and curved panels,Part 1: Modeling and experiments

Theoretical and experimental non-linear vibrations of thin rectangular plates and curved panels subjected to out-of-plane harmonic excitation are investigated. Experiments have been performed on isotropic and laminated sandwich plates and panels with supported and free boundary conditions. A sophisticated measuring technique has been developed to characterize the non-linear behavior experimentally by using a Laser Doppler Vibrometer and a stepped-sine testing procedure. The theoretical approach is based on Donnell's non-linear shell theory (since the tested plates are very thin) but retaining in-plane inertia, taking into account the effect of geometric imperfections. A unified energy approach has been utilized to obtain the discretized non-linear equations of motion by using the linear natural modes of vibration. Moreover, a pseudo arc-length continuation and collocation scheme has been used to obtain the periodic solutions and perform bifurcation analysis. Comparisons between numerical simulations and the experiments show good qualitative and quantitative agreement. It is found that, in order to simulate large-amplitude vibrations, a damping value much larger than the linear modal damping should be considered. This indicates a very large and non-linear increase of damping with the increase of the excitation and vibration amplitude for plates and curved panels with different shape, boundary conditions and materials.     

An invariant manifold approach for studying waves in a one-dimensional array of non-linear oscillators

We consider a one-dimensional linear spring-mass array coupled to a one-dimensional array of uncoupled pendula. The principal aim of this study is to investigate the non-linear dynamics of this large-scale system in the limit of weak non-linearities, i.e. when the (fast) non-linear pendulum effects are small compared to the underlying (slow) linear dynamics of the linear spring-mass chain. We approach the dynamics in the context of invariant manifolds of motion. In particular, we prove the existence of an invariant manifold containing the (predominantly) slow dynamics of the system, with the fast pendulum dynamics providing small perturbations to the motions on the invariant manifold. By restricting the motion on the slow invariant manifold and performing asymptotic analysis we prove that the non-linear large-scale system possesses propagation and attenuation zones (PZs and AZs) in the frequency domain, similarly to the corresponding zones of the linearized system. Inside PZs non-linear travelling wave solutions exist, whereas in AZs only attenuating waves are permissible.     

Resonant non-linear normal modes. Part I: analytical treatment for structural one-dimensional systems

Approximations of the resonant non-linear normal modes of a general class of weakly non-linear one-dimensional continuous systems with quadratic and cubic geometric non-linearities are constructed for the cases of two-to-one, one-to-one, and three-to-one internal resonances. Two analytical approaches are employed: the full-basis Galerkin discretization approach and the direct treatment, both based on use of the method of multiple scales as reduction technique. The procedures yield the uniform expansions of the displacement field and the normal forms governing the slow modulations of the amplitudes and phases of the modes. The non-linear interaction coefficients appearing in the normal forms are obtained in the form of infinite series with the discretization approach or as modal projections of second-order spatial functions with the direct approach. A systematic discussion on the existence and stability of coupled/uncoupled non-linear normal modes is presented. Closed-form conditions for non-linear orthogonality of the modes, in a global and local sense, are discussed. A mechanical interpretation of these conditions in terms of virtual works is also provided.     

A novel statistical linearization solution for randomly excited coupled bending-torsional beams resting on non-linear supports

A novel statistical linearization technique is developed for computing stationary response statistics of randomly excited coupled bending-torsional beams resting on non-linear elastic supports. The key point of the proposed technique consists in representing the non-linear coupled response in terms of constrained linear modes. The resulting set of non-linear equations governing the modal amplitudes is then replaced by an equivalent linear one via a classical statistical error minimization procedure, which provides algebraic non-linear equations for the second-order statistics of the beam response, readily solved by a simple iterative scheme. Data from Monte Carlo simulations, generated by a pertinent boundary integral method in conjunction with a Newmark numerical integration scheme, are used as benchmark solutions to check accuracy and reliability of the proposed statistical linearization technique.

     

Evaluation of multi-modes for finite-element models: systems tuned into 1:2 internal resonance

A non-linear multi-mode of vibration arises from the coupling of two or more normal modes of a non-linear system under free-vibration. The ensuing motion takes place on a 2M-dimensional invariant manifold in the phase space of the system, M being the number of coupled linear modes; the manifold contains a stable equilibrium point of interest, and at that point is tangent to the 2M-dimensional eigenspace of the system linearised about that equilibrium point, which characterises the corresponding M linear modes. On this manifold, M pairs of state variables govern the dynamics of the system; that is, the system behaves like an M-degree-of-freedom oscillator. Non-linear multi-modes may therefore come about when the system exhibits non-linear coupling among generalised co-ordinates. That is the case, for instance, of internal resonance of the 1:2 or 1:3 types, for systems with quadratic or cubic non-linearities, respectively, in which a four-dimensional manifold should be determined. Evaluation of non-linear multi-modes poses huge computational challenges, which is the explanation for very limited reports on the subject in the literature so far. The authors developed a procedure to determine the non-linear multi-modes for finite-element models of plane frames, using the method of multiple scales. This paper refers to the case of quadratic non-linearities. The results obtained by the proposed technique are in good agreement with those coming out from direct integration of the equations of motion in the time domain and also with those few available in the literature.     

Nonlinear normal modes and their superposition in a two degrees of freedom asymmetric system with cubic nonlinearities

This paper investigates nonlinear normal modes and their superposition in a two degrees of freedom asymmetric system with cubic nonlinearities for all nonsingular conditions, based on the invariant subspace in nonlinear normal modes for the nonlinear equations of motion. The focus of attention is to consider relation between the validity of superposition and the static bifurcation of modal dynamics. The numerical results show that the validity has something to do not only with its local restriction, but also with the static bifurcation of modal dynamics. Project Supported by the National Natural Science Foundation and PSF of China     

Non-linear modal analysis for beams subjected to axial loads: Analytical and finite-element solutions

A rigorous derivation of non-linear equations governing the dynamics of an axially loaded beam is given with a clear focus to develop robust low-dimensional models. Two important loading scenarios were considered, where a structure is subjected to a uniformly distributed axial and a thrust force. These loads are to mimic the main forces acting on an offshore riser, for which an analytical methodology has been developed and applied. In particular, non-linear normal modes (NNMs) and non-linear multi-modes (NMMs) have been constructed by using the method of multiple scales. This is to effectively analyse the transversal vibration responses by monitoring the modal responses and mode interactions. The developed analytical models have been crosschecked against the results from FEM simulation. The FEM model having 26 elements and 77 degrees-of-freedom gave similar results as the low-dimensional (one degree-of-freedom) non-linear oscillator, which was developed by constructing a so-called invariant manifold. The comparisons of the dynamical responses were made in terms of time histories, phase portraits and mode shapes.     

Numerical and experimental dynamic characteristics of thin-film membranes

Presented is a total-Lagrangian displacement-based non-linear finite-element model of thin-film membranes for static and dynamic large-displacement analyses. The membrane theory fully accounts for geometric non-linearities. Fully non-linear static analysis followed by linear modal analysis is performed for an inflated circular cylindrical Kapton membrane tube under different pressures, and for a rectangular membrane under different tension loads at four corners. Finite-element results show that shell modes dominate the dynamics of the inflated tube when the inflation pressure is low, and that vibration modes localized along four edges dominate the dynamics of the rectangular membrane. Numerical dynamic characteristics of the two membrane structures were experimentally verified using a Polytec PI PSV-200 scanning laser vibrometer and an EAGLE-500 8-camera motion analysis system.     

Non-linear modal properties of non-shallow cables

A non-linear mechanical model of non-shallow linearly elastic suspended cables is employed to investigate the non-linear modal characteristics of the free planar motions. An asymptotic analysis of the equations of motion is carried out directly on the partial-differential equations overcoming the drawbacks of a discretization process. The direct asymptotic treatment delivers the approximation of the individual non-linear normal modes. General properties about the non-linearity of the in-plane modes of different type—geometric, elasto-static and elasto-dynamic—are unfolded. The spatial corrections to the considered linear mode shape caused by the quadratic geometric forces are investigated for modes belonging to the three mentioned classes. Moreover, the convergence of Galerkin reduced-order models is discussed and the influence of passive modes is highlighted.     

Non-linear interactions in the flexible multi-body dynamics of cable-supported bridge cross-sections

An elastic section model is proposed to analyze some characteristic issues of the cable-supported bridge dynamics through an equivalent planar multi-body system. The quadratic non-linearities of the four-degree-of-freedom model essentially describe the geometric coupling which may strongly characterize the dynamic interactions of the bridge deck and a pair of identical suspension cables (hangers or stays). The linear modal solution shows that the flexural and torsional modes of the deck (global modes) typically co-exist with symmetric or anti-symmetric modes of the cables (local modes). The combinations of parameters which realize remarkable 2:1:1 internal resonance conditions among one of the global modes (with higher natural frequency) and two local modes (with lower and close natural frequencies) are obtained by virtue of a multiparameter perturbation method. The non-linear response of the resonant systems shows that the global deck motion – directly forced at primary resonance by an external harmonic load – can parametrically excite the local cable motion, when the deck vibration amplitude overcomes the critical value at which a period-doubling bifurcation occurs. The relevant effects of both viscous damping and internal detuning on the instability boundaries are parametrically investigated. All the internal resonance conditions as well as the critical vibration amplitudes are expressed as an explicit, though asymptotically approximate, function of the structural parameters.     

Component Mode Synthesis Using Nonlinear Normal Modes

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.     

Stability and modal analysis of shock/boundary layer interactions

The dynamics of oblique shock wave/turbulent boundary layer interactions is analyzed by mining a large-eddy simulation (LES) database for various strengths of the incoming shock. The flow dynamics is first analyzed by means of dynamic mode decomposition (DMD), which highlights the simultaneous occurrence of two types of flow modes, namely a low-frequency type associated with breathing motion of the separation bubble, accompanied by flapping motion of the reflected shock, and a high-frequency type associated with the propagation of instability waves past the interaction zone. Global linear stability analysis performed on the mean LES flow fields yields a single unstable zero-frequency mode, plus a variety of marginally stable low-frequency modes whose stability margin decreases with the strength of the interaction. The least stable linear modes are grouped into two classes, one of which bears striking resemblance to the breathing mode recovered from DMD and another class associated with revolving motion within the separation bubble. The results of the modal and linear stability analysis support the notion that low-frequency dynamics is intrinsic to the interaction zone, but some continuous forcing from the upstream boundary layer may be required to keep the system near a limit cycle. This can be modeled as a weakly damped oscillator with forcing, as in the early empirical model by Plotkin (AIAA J 13:1036–1040, 1975).     

Two-to-one resonant multi-modal dynamics of horizontal/inclined cables. Part I: Theoretical formulation and model validation

This paper is first of the two papers dealing with analytical investigation of resonant multi-modal dynamics due to 2:1 internal resonances in the finite-amplitude free vibrations of horizontal/inclined cables. Part I deals with theoretical formulation and validation of the general cable model. Approximate nonlinear partial differential equations of 3-D coupled motion of small sagged cables – which account for both spatio-temporal variation of nonlinear dynamic tension and system asymmetry due to inclined sagged configurations – are presented. A multi-dimensional Galerkin expansion of the solution of nonplanar/planar motion is performed, yielding a complete set of system quadratic/cubic coefficients. With the aim of parametrically studying the behavior of horizontal/inclined cables in Part II [25], a second-order asymptotic analysis under planar 2:1 resonance is accomplished by the method of multiple scales. On accounting for higher-order effects of quadratic/cubic nonlinearities, approximate closed-form solutions of nonlinear amplitudes, frequencies and dynamic configurations of resonant nonlinear normal modes reveal the dependence of cable response on resonant/nonresonant modal contributions. Depending on simplifying kinematic modeling and assigned system parameters, approximate horizontal/inclined cable models are thoroughly validated by numerically evaluating statics and non-planar/planar linear/non-linear dynamics against those of the exact model. Moreover, the modal coupling role and contribution of system longitudinal dynamics are discussed for horizontal cables, showing some meaningful effects due to kinematic condensation.     

A Variational Approach to Dynamics of Flexible Multibody Systems

Abstract

This paper presents a variational formulation of constrained dynamics of flexible multibody systems, using a vector-variational calculus approach. Body reference frames are used to define global position and orientation of individual bodies in the system, located and oriented by position of its origin and Euler parameters, respectively. Small strain linear elastic deformation of individual components, relative to their body reference frames, is defined by linear combinations of deformation modes that are induced by constraint reaction forces and normal modes of vibration. A library of kinematic couplings between flexible and/or rigid bodies is defined and analyzed. Variational equations of motion for multibody systems are obtained and reduced to mixed differential-algebraic equations of motion. A space structure that must deform during deployment is analyzed, to illustrate use of the methods developed     

Two-to-one resonant multi-modal dynamics of horizontal/inclined cables. Part II: Internal resonance activation,reduced-order models and nonlinear normal modes

Resonant multi-modal dynamics due to planar 2:1 internal resonances in the non-linear, finite-amplitude, free vibrations of horizontal/inclined cables are parametrically investigated based on the second-order multiple scales solution in Part I [1] (in press). The already validated kinematically non-condensed cable model accounts for the effects of both non-linear dynamic extensibility and system asymmetry due to inclined sagged configurations. Actual activation of 2:1 resonances is discussed, enlightening on a remarkable qualitative difference of horizontal/inclined cables as regards non-linear orthogonality properties of normal modes. Based on the analysis of modal contribution and solution convergence of various resonant cables, hints are obtained on proper reduced-order model selections from the asymptotic solution accounting for higher-order effects of quadratic nonlinearities. The dependence of resonant dynamics on coupled vibration amplitudes, and the significant effects of cable sag, inclination and extensibility on system non-linear behavior are highlighted, along with meaningful contributions of longitudinal dynamics. The spatio-temporal variation of non-linear dynamic configurations and dynamic tensions associated with 2:1 resonant non-linear normal modes is illustrated. Overall, the analytical predictions are validated by finite difference-based numerical investigations of the original partial-differential equations of motion.     

Partially invariant solutions in gas dynamics and implicit equations

This paper studies a nonbarochronic, regular, partially invariant solution (submodel) of rank one and defect two to the equations of gas dynamics which describes spatial unsteady gas motion. The equations of gas dynamics are reduced to an implicit ordinary differential equation of the first order for an auxiliary function and to an integrable system. A complete classification of the irregular singular points of the key equation according to a parameter characterizing the gas flow is given, and transformations of the irregular singular points with variation in the parameter are obtained. Qualitative properties of the solution are investigated and physically interpreted in terms of gas motion. It is shown that there are two modes of motion, one of which is supersonic, and in the second modes, a continuous transition through the speed of sound is possible.