Towards a Technique for Nonlinear Modal Analysis

Author(s):  
Simon A. Neild ◽  
Andrea Cammarano ◽  
David J. Wagg

In this paper we discuss a theoretical technique for decomposing multi-degree-of-freedom weakly nonlinear systems into a simpler form — an approach which has parallels with the well know method for linear modal analysis. The key outcome is that the system resonances, both linear and nonlinear are revealed by the transformation process. For each resonance, parameters can be obtained which characterise the backbone curves, and higher harmonic components of the response. The underlying mathematical technique is based on a near identity normal form transformation. This is an established technique for analysing weakly nonlinear vibrating systems, but in this approach we use a variation of the method for systems of equations written in second-order form. This is a much more natural approach for structural dynamics where the governing equations of motion are written in this form as standard practice. In fact the first step in the method is to carry out a linear modal transformation using linear modes as would typically done for a linear system. The near identity transform is then applied as a second step in the process and one which identifies the nonlinear resonances in the system being considered. For an example system with cubic nonlinearities, we show how the resulting transformed equations can be used to obtain a time independent representation of the system response. We will discuss how the analysis can be carried out with applied forcing, and how the approximations about response frequencies, made during the near-identity transformation, affect the accuracy of the technique. In fact we show that the second-order normal form approach can actually improve the predictions of sub- and super-harmonic responses. Finally we comment on how this theoretical technique could be used as part of a modal testing approach in future work.

Author(s):  
Zhenfang Xin ◽  
S. A. Neild ◽  
D. J. Wagg

The normal form technique is an established method for analysing weakly nonlinear vibrating systems. It involves applying a simplifying nonlinear transform to the first-order representation of the equations of motion. In this paper we consider the normal form technique applied to a forced nonlinear system with the equations of motion expressed in second-order form. Specifically we consider the selection of the linearised natural frequencies on the accuracy of the normal form prediction of sub- and superharmonic responses. Using the second-order formulation offers specific advantages in terms of modeling lightly damped nonlinear dynamic response. In the second-order version of the normal form, one of the approximations made during the process is to assume that the linear natural frequency for each mode may be replaced with the response frequencies. Here we will show that this step, far from reducing the accuracy of the technique, does not affect the accuracy of the predicted response at the forcing frequency and actually improves the predictions of sub and superharmonic responses. To gain insight into why this is the case, we consider the Duffing oscillator. The results show that the second-order approach gives an intuitive model of the nonlinear dynamic response which can be applied to engineering applications with weakly nonlinear characteristics.


Energies ◽  
2020 ◽  
Vol 13 (5) ◽  
pp. 1249 ◽  
Author(s):  
Nnaemeka Sunday Ugwuanyi ◽  
Xavier Kestelyn ◽  
Bogdan Marinescu ◽  
Olivier Thomas

Increasing nonlinearity in today’s grid challenges the conventional small-signal (modal) analysis (SSA) tools. For instance, the interactions among modes, which are not captured by SSA, may play significant roles in a stressed power system. Consequently, alternative nonlinear modal analysis tools, notably Normal Form (NF) and Modal Series (MS) methods are being explored. However, they are computation-intensive due to numerous polynomial coefficients required. This paper proposes a fast NF technique for power system modal interaction investigation, which uses characteristics of system modes to carefully select relevant terms to be considered in the analysis. The Coefficients related to these terms are selectively computed and the resulting approximate model is computationally reduced compared to the one in which all the coefficients are computed. This leads to a very rapid nonlinear modal analysis of the power systems. The reduced model is used to study interactions of modes in a two-area power system where the tested scenarios give same results as the full model, with about 70% reduction in computation time.


1994 ◽  
Vol 116 (1) ◽  
pp. 291-297 ◽  
Author(s):  
A. I. Mahyuddin ◽  
A. Midha ◽  
A. K. Bajaj

A method to study parametric stability of flexible cam-follower systems is developed. This method is applied to an automotive valve train which is modeled as a single-degree-of-freedom vibration system. The inclusion of the transverse and rotational flexibilities of the camshaft results in a system governed by a second-order, linear, ordinary differential equation with time-dependent coefficients. This class of equations, known as Hill’s equations, merits special notice in determination of the system response and stability. The analysis includes development of the equivalent model of the system, derivation of its equation of motion, and a method to evaluate its parametric stability based on Floquet theory. A closed-form numerical algorithm, developed to compute the periodic response of systems governed by second-order, linear, ordinary differential equations of motion with time-dependent coefficients, is utilized. The results of this study are presented in a companion paper in the forms of parametric stability charts and three-dimensional stability and response charts.


Author(s):  
Simon A. Neild ◽  
David J. Wagg

Vibration problems are naturally formulated with second-order equations of motion. When the vibration problem is nonlinear in nature, using normal form analysis currently requires that the second-order equations of motion be put into first-order form. In this paper, we demonstrate that normal form analysis can be carried out on the second-order equations of motion. In addition, for forced, damped, nonlinear vibration problems, we show that the invariance properties of the first- and second-order transforms differ. As a result, using the second-order approach leads to a simplified formulation for forced, damped, nonlinear vibration problems.


1993 ◽  
Vol 115 (2) ◽  
pp. 314-323 ◽  
Author(s):  
A. G. Jablokow ◽  
S. Nagarajan ◽  
D. A. Turcic

This paper presents a modal analysis solution technique for the matrix equations of motion of elastic mechanism systems including the rigid-body elastic motion coupling terms and general damping. In many cases, researchers have neglected these terms because they complicate and diminish the efficiency of the solution process. This has been justified by assuming that the terms are small and do not affect the system response. The results obtained using the techniques developed within show this is true for some, but not all mechanisms. The solution technique adds the rigid-body elastic motion coupling terms to the system mass, damping, and stiffness matrices thus allowing them to affect the system response appropriately. The resulting nonsymmetric matrices are then rewritten in first order form allowing general damping to be included in the analysis. Modal analysis techniques are utilized to solve for the steady-state response of the elastic mechanism system. Thus the methods developed in this work provide a technique for including the rigid-body elastic motion coupling terms and general damping in the equations of motion while maintaining the advantage of using efficient modal analysis techniques of finding the response of the system. A number of examples are presented that establish the validity of this approach to the solution of the matrix equations of motion for elastic mechanism systems. The results show that the rigid-body elastic motion coupling terms can become significant at higher operating speeds.


2020 ◽  
Vol 35 (4) ◽  
pp. 3247-3257 ◽  
Author(s):  
Nnaemeka Sunday Ugwuanyi ◽  
Xavier Kestelyn ◽  
Olivier Thomas ◽  
Bogdan Marinescu ◽  
Arturo Roman Messina

1995 ◽  
Vol 117 (1) ◽  
pp. 38-45 ◽  
Author(s):  
J. D. Hicks ◽  
A. W. Troesch ◽  
C. Jiang

The high speeds, small trim angles, and shallow drafts of planing hulls produce large changes in vessel wetted surface which, in turn, lead to significant hydrodynamic and dynamic nonlinearities. Due to the complex nonlinearities of this type of craft, naval architects and planing boat designers tend to rely upon experimental tests or simulation for guidance. In order for simulation to be an effective design tool, a fundamental understanding of the system’s dynamic characteristics is required. This paper describes a developing methodology by which the necessary insight may be obtained. A demonstration of the combined use of modern methods of dynamical system analysis with simulation is given in the evaluation of the vertical motions of a typical planing hull. Extending the work of Troesch and Hicks (1992) and Troesch and Falzarano (1993), the complete nonlinear hydrodynamic force and moment equations of Zarnick (1978) are expanded in a multi-variable Taylor series. As a result, the nonlinear integro-differential equations of motion are replaced by a set of highly coupled, ordinary differential equations with constant coefficients, valid through third order. Closed-form, analytic expressions are available for the coefficients (Hicks, 1993). Numerical examples for all first-order and some second-order terms are presented. Once completely determined, the coefficient matrices will serve as input to path following or continuation methods (e.g., Seydel, 1988) where heave and pitch magnification curves can be generated, allowing the entire system response to be viewed. The branching behavior of the solutions resulting from a variation of the center of gravity is examined in detail. These studies of the second-order accurate model show the potential of the method to identify areas of critical dynamic response, which in turn can be verified and explored further through the use of the simulator.


1992 ◽  
Vol 114 (3) ◽  
pp. 369-374 ◽  
Author(s):  
S. M. Shahruz ◽  
G. Langari

The simplest technique of decoupling the normalized equations of motion of a linear nonclassically damped second-order system is to neglect the off-diagonal elements of the normalized damping matrix. In this paper, the error introduced in the system response due to this decoupling technique is studied. Conditions are derived under which the solution of an approximately decoupled system is “close” to the exact solution of the system.


2021 ◽  
Vol 502 (3) ◽  
pp. 3976-3992
Author(s):  
Mónica Hernández-Sánchez ◽  
Francisco-Shu Kitaura ◽  
Metin Ata ◽  
Claudio Dalla Vecchia

ABSTRACT We investigate higher order symplectic integration strategies within Bayesian cosmic density field reconstruction methods. In particular, we study the fourth-order discretization of Hamiltonian equations of motion (EoM). This is achieved by recursively applying the basic second-order leap-frog scheme (considering the single evaluation of the EoM) in a combination of even numbers of forward time integration steps with a single intermediate backward step. This largely reduces the number of evaluations and random gradient computations, as required in the usual second-order case for high-dimensional cases. We restrict this study to the lognormal-Poisson model, applied to a full volume halo catalogue in real space on a cubical mesh of 1250 h−1 Mpc side and 2563 cells. Hence, we neglect selection effects, redshift space distortions, and displacements. We note that those observational and cosmic evolution effects can be accounted for in subsequent Gibbs-sampling steps within the COSMIC BIRTH algorithm. We find that going from the usual second to fourth order in the leap-frog scheme shortens the burn-in phase by a factor of at least ∼30. This implies that 75–90 independent samples are obtained while the fastest second-order method converges. After convergence, the correlation lengths indicate an improvement factor of about 3.0 fewer gradient computations for meshes of 2563 cells. In the considered cosmological scenario, the traditional leap-frog scheme turns out to outperform higher order integration schemes only when considering lower dimensional problems, e.g. meshes with 643 cells. This gain in computational efficiency can help to go towards a full Bayesian analysis of the cosmological large-scale structure for upcoming galaxy surveys.


Sign in / Sign up

Export Citation Format

Share Document