New partial differential equations governing the joint, response-excitation, probability distributions of nonlinear systems, under general stochastic excitation

In the present work the problem of determining the probabilistic structure of the dynamical response of nonlinear systems Subjected to general, external, stochastic excitation is considered. The starting point Of Our approach is a Hopf-type equation, governing the evolution of the joint, response-excitation, characteristic functional. Exploiting this equation, we derive new linear partial differential equations governing the joint, response-excitation, characteristic (or probability density) function, which can be considered as an extension of the well-known Fokker-Planck-Kolmogorov equation to the case of a general, correlated excitation and, thus, non-Markovian response character. These new equations are Supplemented by initial conditions and a marginal compatibility condition (with respect to the known probability distribution of the excitation), which is of non-local character. The validity of this new equation is also checked by showing its equivalence with the infinite system of moment equations. The method is applicable to any differential system, in state-space form, exhibiting polynomial nonlinearities. In this paper the method is illustrated through a detailed analysis of a simple, first-order, scalar equation, with a cubic nonlinearity. It is also shown that various versions of Fokker-Planck-Kolmogorov equation, corresponding to the case of independent-increment excitations, can be derived by using the same approach. A numerical method for the solution of these new equations is introduced and illustrated through its application to the simple model problem. It is based on the representation of the joint probability density (or characteristic) function by means of a convex superposition of kernel functions, which permits us to satisfy a priori the non-local marginal compatibility condition. On the basis of this representation, the partial differential equation is eventually transformed to a system of ordinary differential equations for the kernel parameters. Extension to general, multidimensional, dynamical systems exhibiting any polynomial nonlinearity will be presented in a forthcoming paper. (C) 2008 Elsevier Ltd. All rights reserved.