Dynamic stochastic finite element method using time‐dependent generalized polynomial chaos. (20th August 2020)

• Record Type: Journal Article
• Authors: Lacour, Maxime; Bal, Guillaume; Abrahamson, Norman
• Is Part Of: International journal for numerical and analytical methods in geomechanics. Volume 45:Number 3(2021)
• Page Start: 293
• Page End: 306
• Publication Date: 2020-08-20
• Subjects: dynamic polynomial chaos -- stochastic finite elements; Soil mechanics -- Mathematics -- Periodicals; Rock mechanics -- Mathematics -- Periodicals; 624.1510151
• DOI: 10.1002/nag.3137
• Languages: English
• ISSNs: 0363-9061

Abstract

Summary: We present an intrusive formulation for the dynamic stochastic finite‐element method to propagate the epistemic uncertainty in material properties into a finite‐element system over time. The stochastic finite‐element method, originally developed for the static case, uses generalized polynomial chaos (gPC) expansions to represent the uncertainty in both material/load fields and displacement fields and solves for the unknown PC coefficients of displacement at each degree of freedom of the finite‐element system. In this case, the gPC basis used to represent the solution is optimal and can be kept the same throughout the static simulation; however, when integrating a stochastic system over time, it is proven that using gPC tends to break down for integrations over long times. The reason is that the solution's complexity increases in time, and the set of polynomials used in the gPC expansion to represent the solution, therefore, does not stay optimal. In this formulation, new stochastic variables and orthogonal polynomials are constructed as time progresses. These variables are obtained as the Kahrunen‐Loeve expansion of the finite‐element solution at the times of update, thus, optimally representing the distribution of the solution using a minimal set of orthogonal polynomials. The result from the method is a time‐domain polynomial chaos representation of the entire finite‐element solution. A fast post‐processing phase can be used to (a) obtain the probability distribution of the solution at each degree of freedom and (b) generate any number of time series realizations of the solution, which correspond to the same time series that would be obtained from a set of simulations based on direct Monte‐Carlo simulations of the uncertain material properties.