Abstract
The response of a structural system to white noise excitation (deltacorrelated) constitutes a Markov vector process whose transitional probability density function (TPDF) is governed by both the forward Fokker-Planck and backward Kolmogorov equations. Numerical solution of these equations by finite element and finite difference methods for dynamical systems of engineering interest has been hindered by the problem of dimensionality. In this paper numerical solution of the stationary and transient form of the Fokker-Planck (FP) equation corresponding to two state nonlinear systems is obtained by standard sequential finite element method (FEM) using C0 shape function and Crank-Nicholson time integration scheme. The method is applied to Van-der-Pol and Duffing oscillators providing good agreement between results obtained by it and exact results. An extension of the finite difference discretization scheme developed by Spencer, Bergman and Wojtkiewicz is also presented. This paper presents an extension of the finite difference method for the solution of FP equation up to four dimensions. The difficulties associated in extending these methods to higher dimensional systems are discussed.
Similar content being viewed by others
References
Caughey T K 1971 Nonlinear theory of random vibrations. InAdvances in applied mechanics (ed.) C-S Yih (New York: Academic Press) 11: 209–253
Caughey T K, Dienes J K 1962 The behaviour of linear systems with white noise input.J. Math. Phys. 32: 2476–2479
Crandall S H 1963 Perturbation techniques for random vibration of nonlinear systems.J. Acoust. Soc. Am. 35: 1700–1705
Johnson E A, Wojtkiewicz S F, Bergman L A, Spencer B F Jr 1997a Observations with regard to massively parallel computation for Monte Carlo simulation of stochastic dynamical systems.Int. J. Nonlinear Mech. 32(4): 721–734
Johnson E A, Wojtkiewicz S F, Spencer B F Jr, Bergman LA1 997b Finite element and finite difference solutions to the transient Fokker—Planck equation. DESY-97-161: 290–306
Khasminskii R Z 1996 A limit theorem for the solution of differential equations with random righthand sides.Theor. Probab. Appl. 11: 390–405
Kunert A 1991 Efficient numerical solution of multi-dimensional Fokker-Planck equations with chaotic and nonlinear random vibration. InVibration analysis -Analytical and computational (eds.) T C Huanget al DE-37: 57–60
Langley R S 1985 A finite element method for the statistics of nonlinear random vibration.J. Sound Vibr. 101:41–54
Langtangeh HP 1991 A general numerical solution method for Fokker-Planck equations with applications to structural reliability.Probab. Eng. Mech. 6: 33–48
Lin Y K 1967Probabilistic theory of structural dynamics (New York: McGraw Hill)
Naess A, Johnson J M 1992 Response statistics of nonlinear dynamics systems by path integration. InProc. IUTAM Symp. on Nonlinear Stochastic Mechanics, (eds) N Bellomo, F Casciatti (Berlin: Springer-Verlag) pp 401–414
Risken H 1989The Fokker-Planck equation: Methods of solution and applications (New York: Springer-Verlag)
Roberts J B, Spanos P D 1990Random vibration and statistical linearization (Wiley: New York)
Schueller G I, Pradlwarter H J, Pandey M D 1993 Methods for reliability assessment of nonlinear system under stochastic dynamic loading—A review. InStructural dynamics: Proc. Second European Conf. on Structural Dynamics: EURODYN’93 (eds) T Moanet al 1993 (Rotterdam: Balkema) pp 751–759
Soong T 1994The Fokker-Planck equation for stochastic dynamical system and its explicit steady state solution (Singapore: World Scientific)
Spencer Jr B F, Bergman L A 1993 On the numerical solution of the Fokker equations for nonlinear stochastic systems.Nonlinear Dynamics 4: 357–372
Sun J Q, Hsu C S 1988 First-passage time probability of nonlinear stochastic systems by generalized cell mapping method.J. Sound Vibr. 124: 233–248
Sun J-Q, Hsu C S 1990 The generalized cell mapping method in nonlinear random vibration based upon short-time Gaussian approximation.ASME J. Appl. Mech. 57: 1018–1025
Wang M C, Uhlenbeck G 1945 On the theory of Brownian motion II.Rev. Modern Phys. 17: 323–342 (reprinted 1954 inSelected papers on noise and stochastic processes (ed.) N Wax (New York: Dover)
Wen Y K 1976 Method for random vibration of hysteretic systems.ASCE J. Engineering Mech. 102(EM2): 249–263
Wojtkiewicz S F, Bergman L A, Spencer B F Jr 1994a Robust numerical solution of the Fokker-Planck-Kolmogorov equation for two-dimensional stochastic dynamical systems. Technical Report AAE 94-08, Department of Aeronautical and Astronautical Engineering, University of Illinois at Urbana-Champaign, Urbana-Champaign, IL
Wojtkiewicz S F, Bergman L A, Spencer B F Jr 1994b High fidelity numerical solutions of the Fokker-Planck equation.Proc. ICOSSAR ’97, the 7th Int. Conf on Structural Safety and Reliability, Kyoto, Japan
Wojtkiewicz S F, Spencer B F Jr, Bergman L A 1996 On the cumulant-neglect closure method in stochastic dynamics.Int. J. Nonlinear Mech. 31: 657–684
Wong E, Zakai M 1965 On the relation between ordinary and stochastic differential equation.Int. J. Eng. Sci. 31965: 213–229
Wu W F, Lin Y K 1984 Cumulant-neglect closure for nonlinear oscillators under parametric and external excitations.Int. J. Nonlinear Mech. 19: 349–362
Author information
Authors and Affiliations
Additional information
This paper is dedicated to Prof R N Iyengar of the Indian Institute of Science on the occasion of his formal retirement.
Rights and permissions
About this article
Cite this article
Kumar, P., Narayanan, S. Solution of Fokker-Planck equation by finite element and finite difference methods for nonlinear systems. Sadhana 31, 445–461 (2006). https://doi.org/10.1007/BF02716786
Issue Date:
DOI: https://doi.org/10.1007/BF02716786