Perturbation methods in evolutionary spectral analysis for linear dynamics and equivalent statistical linearization

Abstract The analysis of large-scale structures subject to transient random loads, coherent in space and time, is a classic problem encountered in earthquake and wind engineering. The simulation-based framework is usually seen as the most convenient approach for both linear and nonlinear dynamics. However, the generation of statistically consistent samples of an excitation field remains a heavy computational task. In light of this, perturbation techniques are applied to develop and improve evolutionary spectral analysis. Advantageously performed in a standard modal basis, this evolutionary spectral analysis for linear structures requires the computation of the modal impulse response matrix. However, this matrix has no general closed-form expression in the presence of modal coupling. We propose therefore to model it by an asymptotic approximation, obtained by the inverse Fourier transform of an asymptotic expansion of the modal transfer matrix of the structure. This latter expansion considers the modal coupling as a perturbation of a main decoupled system. This strategy leads to an expansion known in a closed-form. Finally, the semi-group property allows the use of an efficient recurrence relation to approximate the modal evolutionary transfer matrix, i.e. the evolutionary extension of the transfer matrix. The asymptotic expansion-based method and the recurrence relation are then applied to nonlinear transient dynamics by using Gaussian equivalent linearization. This extension is formalized by a multiple timescales approach, allowing to consider a linearized structure, namely a time variant system, as piecewise linear time invariant depending on a statistical timescale. The proposed developments are finally illustrated on realistic civil engineering applications

A Method for the Combination of Stochastic Time Varying Load Effects

The problem of evaluating the probability that a structure becomes unsafe under a combination of loads, over a given time period, is addressed. The loads and load effects are modeled as either pulse (static problem) processes with random occurrence time, intensity and a specified shape or intermittent continuous (dynamic problem) processes which are zero mean Gaussian processes superimposed 'on a pulse process. The load coincidence method is extended to problems with both nonlinear limit states and dynamic responses, including the case of correlated dynamic responses. The technique of linearization of a nonlinear limit state commonly used in a time-invariant problem is investigated for timevarying combination problems, with emphasis on selecting the linearization point. Results are compared with other methods, namely the method based on upcrossing rate, simpler combination rules such as Square Root of Sum of Squares and Turkstra's rule. Correlated effects among dynamic loads are examined to see how results differ from correlated static loads and to demonstrate which types of load dependencies are most important, i.e., affect' the exceedance probabilities the most. Application of the load coincidence method to code development is briefly discussed.National Science Foundation Grants CME 79-18053 and CEE 82-0759

Mean-Variance Cointegration and the Expectations Hypothesis

The present paper sheds further light on a well-known (alleged) violation of the expec- tations hypothesis of the term structure (EHT) - the frequent finding of unit roots in interest rate spreads. We show that the EHT implies (i) that the nonstationarity stems from the holding premium, which is hence (ii) cointegrated with the spread. In a stochas- tic discount factor framework we model the premium as being driven by the integrated variance of excess returns. Introducing the concept of mean-variance cointegration we actually find cointegration relations between spreads and premia in US data.Expectations Hypothesis, Holding Premium, Persistence, Cointegration, GARCH

Frequency domain analysis method of nonstationary random vibration based on evolutionary spectral representation

PurposeA novel frequency domain approach, which combines the pseudo excitation method modified by the authors and multi-domain Fourier transform (PEM-FT), is proposed for analyzing nonstationary random vibration in this paper.Design/methodology/approachFor a structure subjected to a nonstationary random excitation, the closed-form solution of evolutionary power spectral density of the response is derived in frequency domain.FindingsThe deterministic process and random process in an evolutionary spectrum are separated effectively using this method during the analysis of nonstationary random vibration of a linear damped system, only modulation function of the system needs to be estimated, which brings about a large saving in computational time.Originality/valueThe method is general and highly flexible as it can deal with various damping types and nonstationary random excitations with different modulation functions.</jats:sec

A new SSI algorithm for LPTV systems: Application to a hinged-bladed helicopter

Many systems such as turbo-generators, wind turbines and helicopters show intrinsic time-periodic behaviors. Usually, these structures are considered to be faithfully modeled as linear time-invariant (LTI). In some cases where the rotor is anisotropic, this modeling does not hold and the equations of motion lead necessarily to a linear periodically time- varying (referred to as LPTV in the control and digital signal field or LTP in the mechanical and nonlinear dynamics world) model. Classical modal analysis methodologies based on the classical time-invariant eigenstructure (frequencies and damping ratios) of the system no more apply. This is the case in particular for subspace methods. For such time-periodic systems, the modal analysis can be described by characteristic exponents called Floquet multipliers. The aim of this paper is to suggest a new subspace-based algorithm that is able to extract these multipliers and the corresponding frequencies and damping ratios. The algorithm is then tested on a numerical model of a hinged-bladed helicopter on the ground

Response statistics and failure probability determination of nonlinear stochastic structural dynamical systems

Novel approximation techniques are proposed for the analysis and evaluation of nonlinear dynamical systems in the field of stochastic dynamics. Efficient determination of response statistics and reliability estimates for nonlinear systems remains challenging, especially those with singular matrices or endowed with fractional derivative elements. This thesis addresses the challenges of three main topics. The first topic relates to the determination of response statistics of multi-degree-of-freedom nonlinear systems with singular matrices subject to combined deterministic and stochastic excitations. Notably, singular matrices can appear in the governing equations of motion of engineering systems for various reasons, such as due to a redundant coordinates modeling or due to modeling with additional constraint equations. Moreover, it is common for nonlinear systems to experience both stochastic and deterministic excitations simultaneously. In this context, first, a novel solution framework is developed for determining the response of such systems subject to combined deterministic and stochastic excitation of the stationary kind. This is achieved by using the harmonic balance method and the generalized statistical linearization method. An over-determined system of equations is generated and solved by resorting to generalized matrix inverse theory. Subsequently, the developed framework is appropriately extended to systems subject to a mixture of deterministic and stochastic excitations of the non-stationary kind. The generalized statistical linearization method is used to handle the nonlinear subsystem subject to non-stationary stochastic excitation, which, in conjunction with a state space formulation, forms a matrix differential equation governing the stochastic response. Then, the developed equations are solved by numerical methods. The accuracy for the proposed techniques has been demonstrated by considering nonlinear structural systems with redundant coordinates modeling, as well as a piezoelectric vibration energy harvesting device have been employed in the relevant application part. The second topic relates to code-compliant stochastic dynamic analysis of nonlinear structural systems with fractional derivative elements. First, a novel approximation method is proposed to efficiently determine the peak response of nonlinear structural systems with fractional derivative elements subject to excitation compatible with a given seismic design spectrum. The proposed methods involve deriving an excitation evolutionary power spectrum that matches the design spectrum in a stochastic sense. The peak response is approximated by utilizing equivalent linear elements, in conjunction with code-compliant design spectra, hopefully rendering it favorable to engineers of practice. Nonlinear structural systems endowed with fractional derivative terms in the governing equations of motion have been considered. A particular attribute pertains to utilizing localized time-dependent equivalent linear elements, which is superior to classical approaches utilizing standard time-invariant statistical linearization method. Then, the approximation method is extended to perform stochastic incremental dynamical analysis for nonlinear structural systems with fractional derivative elements exposed to stochastic excitations aligned with contemporary aseismic codes. The proposed method is achieved by resorting to the combination of stochastic averaging and statistical linearization methods, resulting in an efficient and comprehensive way to obtain the response displacement probability density function. A stochastic incremental dynamical analysis surface is generated instead of the traditional curves, leading to a reliable higher order statistics of the system response. Lastly, the problem of the first excursion probability of nonlinear dynamic systems subject to imprecisely defined stochastic Gaussian loads is considered. This involves solving a nested double-loop problem, generally intractable without resorting to surrogate modeling schemes. To overcome these challenges, this thesis first proposes a generalized operator norm framework based on statistical linearization method. Its efficiency is achieved by breaking the double loop and determining the values of the epistemic uncertain parameters that produce bounds on the probability of failure a priori. The proposed framework can significantly reduce the computational burden and provide a reliable estimate of the probability of failure

Solving, Estimating and Selecting Nonlinear Dynamic Economic Models without the Curse of Dimensionality

A welfare analysis of a risky policy is impossible within a linear or linearized model and its certainty equivalence property. The presented algorithms are designed as a toolbox for a general model class. The computational challenges are considerable and I concentrate on the numerics and statistics for a simple model of dynamic consumption and labor choice. I calculate the optimal policy and estimate the posterior density of structural parameters and the marginal likelihood within a nonlinear state space model. My approach is even in an interpreted language twenty time faster than the only alternative compiled approach. The model is estimated on simulated data in order to test the routines against known true parameters. The policy function is approximated by Smolyak Chebyshev polynomials and the rational expectation integral by Smolyak Gaussian quadrature. The Smolyak operator is used to extend univariate approximation and integration operators to many dimensions. It reduces the curse of dimensionality from exponential to polynomial growth. The likelihood integrals are evaluated by a Gaussian quadrature and Gaussian quadrature particle filter. The bootstrap or sequential importance resampling particle filter is used as an accuracy benchmark. The posterior is estimated by the Gaussian filter and a Metropolis- Hastings algorithm. I propose a genetic extension of the standard Metropolis-Hastings algorithm by parallel random walk sequences. This improves the robustness of start values and the global maximization properties. Moreover it simplifies a cluster implementation and the random walk variances decision is reduced to only two parameters so that almost no trial sequences are needed. Finally the marginal likelihood is calculated as a criterion for nonnested and quasi-true models in order to select between the nonlinear estimates and a first order perturbation solution combined with the Kalman filter.stochastic dynamic general equilibrium model, Chebyshev polynomials, Smolyak operator, nonlinear state space filter, Curse of Dimensionality, posterior of structural parameters, marginal likelihood

Reliability problems in eartquake engineering

This monograph deals with the problem of reliability analysis in the field of Earthquake Engineering. Chapter 1 is devoted to a summary of the most widely used reliability methods, with emphasis on Monte Carlo and solver surrogate techniques used in the subsequent chapters. Chapter 2 presents a discussion of the Monte Carlo from the viewpoint of Information Theory. Then, a discussion is made in Chapter 3 on the selection of random variables in Earthquake Engineering. Next, some practical methods for computing failure probabilities under seismic loads are reported in Chapter 4. Finally, a method for reliability-based design optimization under seismic loads is presented in Chapter 5

Deterministic continutation of stochastic metastable equilibria via Lyapunov equations and ellipsoids

Numerical continuation methods for deterministic dynamical systems have been one of the most successful tools in applied dynamical systems theory. Continuation techniques have been employed in all branches of the natural sciences as well as in engineering to analyze ordinary, partial and delay differential equations. Here we show that the deterministic continuation algorithm for equilibrium points can be extended to track information about metastable equilibrium points of stochastic differential equations (SDEs). We stress that we do not develop a new technical tool but that we combine results and methods from probability theory, dynamical systems, numerical analysis, optimization and control theory into an algorithm that augments classical equilibrium continuation methods. In particular, we use ellipsoids defining regions of high concentration of sample paths. It is shown that these ellipsoids and the distances between them can be efficiently calculated using iterative methods that take advantage of the numerical continuation framework. We apply our method to a bistable neural competition model and a classical predator-prey system. Furthermore, we show how global assumptions on the flow can be incorporated - if they are available - by relating numerical continuation, Kramers' formula and Rayleigh iteration.Comment: 29 pages, 7 figures [Fig.7 reduced in quality due to arXiv size restrictions]; v2 - added Section 9 on Kramers' formula, additional computations, corrected typos, improved explanation