Abstract: An Adomian decomposition based mathematical framework to derive the mean square responses of nonlinear structural systems subjected to stochastic excitation is presented. The exact mean square response estimation of certain class of nonlinear (opens in a new window)stochastic systems is achieved using Fokker–Planck–Kolmogorov (FPK) equations resulting in analytical expressions or using (opens in a new window)Monte Carlo simulations. However, for most of the (opens in a new window)nonlinear systems, the response estimation using (opens in a new window)Monte Carlo simulations is computationally expensive, and, also, obtaining solution of FPK equation is mathematically exhaustive owing to the requirement to solve a stochastic (opens in a new window)partial differential equation. In this context, the present work proposes an Adomian decomposition based formalism to derive semi-analytical expressions for the second order response statistics. Further, a derivative matching based moment approximation technique is employed to reduce the higher order moments in (opens in a new window)nonlinear systems into functions of lower order moments without resorting to any sort of linearization. Three case studies consisting of Duffing oscillator with negative stiffness, Rayleigh Van-der Pol oscillator and a Pendulum (opens in a new window)tuned mass damper inerter system with linear (opens in a new window)auxiliary spring–damper arrangement subjected to white noise excitation are undertaken. The accuracy of the closed form expressions derived using the proposed framework is established by comparing the mean square responses of the systems with the exact solutions. The results demonstrate the robustness of the proposed framework for accurate statistical analysis of (opens in a new window)nonlinear systems under stochastic excitation.