This article discusses the estimation of Lyapunov characteristic exponents related to the solution of non-linear multibody dynamics problems expressed in the form of systems of differential–algebraic equations. Current literature on the topic projects the differential–algebraic equations in the state space of the underlying strangeness free ordinary differential equations. This article shows how Lyapunov exponents can be successfully estimated with the discrete QR method directly applied to the matrix pencil associated to the differential–algebraic equations, using algorithmically dissipative multistep methods. Numerical results show the soundness of the proposed approach.
DieciLVan VleckESEstepDTavernerS. Lyapunov and other spectra: a surveyCollected lectures on the preservation of stability under discretization, chapter 112002Philadelphia, PASIAM197–218.
2.
DieciLRussellRDVan VleckES. On the computation of Lyapunov exponents for continuous dynamical systemsSIAM J Numer Anal1997; 341402–423.
Van VleckESDieciLJollyM. Codes for approximating Lyapunov exponentsASME IDETC/CIE 200930 August–2 September 2009San Diego, CADETC2009–86530.
5.
ShimadaINagashimaT. A numerical approach to ergodic problem of dissipative dynamical systemsProg Theor Phys1979; 6161605–1616.
6.
BenettinGGalganiLGiorgilliA. Lyapunov characteristic exponents for smooth dynamical systems and for Hamiltonian systems; a method for computing all of them. Part 1: TheoryMeccanica1980; 1519–20.
7.
ChristiansenFRughHH. Computing Lyapunov spectra with continuous Gram-Schmidt orthonormalizationNonlinearity1997; 1051063–1072.
8.
RosensteinMTCollinsJJDe LucaCJ. A practical method for calculating largest Lyapunov exponents from small data setsPhysica D1993; 651–2117–134.
WolfASwiftJBSwinneyHL. Determining Lyapunov exponents from a time seriesPhysica D1985; 163285–317.
11.
SanoMSawadaY. Measurement of the Lyapunov spectrum from a chaotic time seriesPhys Rev Lett1985; 55: 1082–1085.
12.
BroomheadDSKingGP. Extracting qualitative dynamics from experimental dataPhysica D1986; 202–3217–236.
13.
BrownRBryantPAbarbanelHDI. Computing the Lyapunov spectrum of a dynamical system from an observed time seriesPhys Rev A1991; 43: 2787–2806.
14.
AbarbanelHDIBrownRKennelMB. Local Lyapunov exponents computed from observed dataJ Nonlinear Sci1992; 23343–365.
15.
KantzH. A robust method to estimate the maximal Lyapunov exponent of a time seriesPhys Lett A1994; 185177–87.
16.
ParlitzU. Identification of true and spurious Lyapunov exponents from time seriesInt J Bifurcation Chaos1992; 21155–165.
17.
MüllerPC. Calculation of Lyapunov exponents for dynamic systems with discontinuitiesChaos, Solitons Fractals1995; 591671–1681.
18.
AsokanthanSFWangX-H. Characterization of torsional instabilities in a Hooke's joint driven system via maximal Lyapunov exponentsJ Sound Vib1996; 194183–91.
19.
GilatRAboudiJ. Parametric stability of non-linearly elastic composite plates by Lyapunov exponentsJ Sound Vib2000; 2354627–637.
20.
YangCWuQ. On stabilization of bipedal robots during disturbed standing using the concept of Lyapunov exponentsRobotica2006; 24: 621–624.
21.
YangCWuQ. On stability analysis via Lyapunov exponents calculated from a time series using nonlinear mapping – a case studyNonlinear Dyn2010; 59: 239–257.
22.
EberhardPSchiehlenW. Computational dynamics of multibody systems: history, formalisms, and applicationsJ Comput Nonlinear Dyn2006; 1: 3–12.
23.
CongNDNamH. Lyapunov's inequality for linear differential algebraic equationActa Mathematica Vietnamica2003; 28173–88.
24.
CongNDNamH. Lyapunov regularity of linear differential algebraic equations of index 1Acta Math Vietnamica2004; 2911–21.
25.
AwrejcewiczJKudraG. Stability analysis and Lyapunov exponents of a multi-body mechanical system with rigid unilateral constraintsNonlinear Anal2005; 63: e909–e918.
26.
FuSWangQ. Estimating the largest Lyapunov exponent in a multibody system with dry friction by using chaos synchronizationActa Mech Sin2006; 22: 277–283.
27.
BaumgarteJ. Stabilization of constraints and integrals of motion in dynamical systemsComput Methods Appl Mech Eng1972; 1: 1–36.
28.
LinhVHMehrmannV. Lyapunov, Bohl and Sacker-Sell spectral intervals for differential-algebraic equationsJ Dyn Differ Equ2009; 21: 153–194.
29.
OseledecVI. A multiplicative ergodic theorem: Lyapunov characteristic numbers for dynamical systemsTrans Moscow Math Soc1968; 19: 197–231.
30.
MärzR. Notes on linearization of differential-algebraic equations and on optimization with differential-algebraic constraints. Technical Report 11-16, Department of Mathematics, Humboldt University2011.
31.
MasaratiP. Direct eigenanalysis of constrained system dynamicsProc IMechE Part K: J Multi-body Dynamics2009; 2234335–342.
32.
LambertJD. Numerical methods for ordinary differential systems1991Chichester, UKJohn Wiley & Sons Ltd.
33.
StewartGW. On the sensitivity of the eigenvalue problem Ax = λBxSIAM J Numer Anal1972; 94669–686.
34.
MolerCBStewartGW. An algorithm for generalized matrix eigenvalue problemsSIAM J Numer Anal1973; 102241–256.
35.
RipepiMMasaratiP. Reduced order models using generalized eigenanalysisProc IMechE Part K: J Multi-body Dynamics2011; 225: 1–14.
36.
LeimkuhlerBPetzoldLRGearCW. Approximation methods for the consistent initialization of differential-algebraic equationsSIAM J Numer Anal1991; 281205–226.
DormandJRPrincePJ. A family of embedded Runge-Kutta formulaeJ Comput Appl Math1980; 6119–26.
39.
BoschiVDe SalvoRMasaratiP. Seismic attenuation system synthesis by reduced order models from multibody analysis25–28 June 2007Multibody dynamics 2007, ECCOMAS thematic conferenceMilano.
40.
AsokanthanSFHwangMC. Torsional instabilities in a system incorporating a Hooke's jointJ Vibr Acoust1996; 1183368–374.
41.
MasaratiPLanzMMantegazzaP. Multistep integration of ordinary, stiff and differential-algebraic problems for multibody dynamics applicationsXVI Congresso Nazionale AIDAA24–28 September 2001Palermo71.1–71.10.
