Fast matrix exponent for deterministic or random excitations
Author(s)
Date Issued
2001
Publisher
John Wiley & Sons
Journal
International Journal for Numerical Methods in Engineering
Volume
50
Issue
2
Start page
377
End page
394
Abstract
The solution of ż=Az is z(t)=exp(At)z0=Etz0, z0=z(0). Since z(2t)=E2tz0=Et2z0, z(4t)=E4tz0=E2t2z0, etc., one function evaluation can double the time step. For an n-degree-of-freedoms system, A is a 2n matrix of the nth-order mass, damping and stiffness matrices M, C and K. If the forcing term is given as piecewise combinations of the elementary functions, the force response can be obtained analytically. The mean-square response P to a white noise random force with intensity W(t) is governed by the Lyapunov differential equation: Ṗ=AP+PAT+W. The solution of the homogeneous Lyapunov equation is P(t)=exp(At) P0 exp(ATt), P0=P(0). One function evaluation can also double the time step. If W(t) is given as piecewise polynomials, the mean-square response can also be obtained analytically. In fact, exp(At) consists of the impulsive- and step-response functions and requires no special treatment. The method is extended further to coloured noise. In particular, for a linear system initially at rest under white noise excitation, the classical non-stationary response is resulted immediately without integration. The method is further extended to modulated noise excitations. The method gives analytical mean-square response matrices for lightly damped or heavily damped systems without using modal expansion. No integration over the frequency is required for the mean-square response. Four examples are given. The first one shows that the method include the result of Caughy and Stumpf as a particular case. The second one deals with non-white excitation. The third finds the transient stress intensity factor of a gun barrel and the fourth finds the means-square response matrix of a simply supported beam by finite element method.
SFU Affiliated Publication
No
Availability at SFU Library
No database links found.