Probabilistic analysis of the steady state of weakly perturbed linear oscillators subject to a class of Gaussian inputs

• Published in: Chaos, solitons and fractals Vol. 187; p. 115451
• Main Authors: Cortés, J.-C., Romero, J.-V., Roselló, M.-D., Valencia Sullca, J.F.
• Format: Journal Article
• Language: English
• Published: Elsevier Ltd 01.10.2024
• Subjects: Equivalent linearization; Nonlinear oscillator; Perturbation technique; Principle of maximum entropy; Nonlinear oscillator; Principle of maximum entropy; Perturbation technique; Equivalent linearization
• ISSN: 0960-0779; 1873-2887
• DOI: 10.1016/j.chaos.2024.115451

This paper aims to probabilistically study a class of nonlinear oscillator subject to weak perturbations and driven by stationary zero-mean Gaussian stochastic processes. For the sake of generality in the analysis, we assume that the perturbed term is a polynomial of arbitrary degree in the spatial position, that contains, as a particular case, the important case of the Duffing equation. We then take advantage of the so-called stochastic equivalent linearization technique to construct an equivalent linear model so that its behavior consistently approximates, in the mean-square sense, that of the nonlinear oscillator. This approximation allows us to take extensive advantage of the probabilistic properties of the solution of the linear model and its first mean-square derivative to construct reliable approximations of the main statistical moments of the steady state. From this key information, we then apply the principle of maximum entropy to construct approximations of the probability density function of the steady state. We illustrate the superiority of the equivalent linearization technique over the perturbation method through some examples. •A family of nonlinear stochastic oscillators is studied.•Stochastic equivalent linearization and maximum entropy methods are combined.•Mean, variance and all higher moments of the solution are calculated.•An approximation of the stationary probability density function is computed.•Gaussian excitations are considered in examples.