Neumann enriched polynomial chaos approach for stochastic finite element problems

An enrichment scheme based upon the Neumann expansion method is proposed to augment the deterministic coefficient vectors associated with the polynomial chaos expansion method. The proposed approach relies upon a split of the random variables into two statistically independent sets. The principal variability of the system is captured by propagating a limited number of random variables through a low-ordered polynomial chaos expansion method. The remaining random variables are propagated by a Neumann expansion method. In turn, the random variables associated with the Neumann expansion method are utilised to enrich the polynomial chaos approach. The effect of this enrichment is explicitly captured in a new augmented definition of the coefficients of the polynomial chaos expansion. This approach allows one to consider a larger number of random variables within the scope of spectral stochastic finite element analysis in a computationally efficient manner. Closed-form expressions for the first two response moments are provided. The proposed enrichment method is used to analyse two numerical examples: the bending of a cantilever beam and the flow through porous media. Both systems contain distributed stochastic properties. The results are compared with those obtained using direct Monte Carlo simulations and using the classical polynomial chaos expansion approach.     

A non‐intrusive model reduction approach for polynomial chaos expansion using proper orthogonal decomposition

In this paper, a non‐intrusive stochastic model reduction scheme is developed for polynomial chaos representation using proper orthogonal decomposition. The main idea is to extract the optimal orthogonal basis via inexpensive calculations on a coarse mesh and then use them for the fine‐scale analysis. To validate the developed reduced‐order model, the method is implemented to: (1) the stochastic steady‐state heat diffusion in a square slab; (2) the incompressible, two‐dimensional laminar boundary‐layer over a flat plate with uncertainties in free‐stream velocity and physical properties; and (3) the highly nonlinear Ackley function with uncertain coefficients. For the heat diffusion problem, the thermal conductivity of the slab is assumed to be a stochastic field with known exponential covariance function and approximated via the Karhunen–Loève expansion. In all three test cases, the input random parameters are assumed to be uniformly distributed, and a polynomial chaos expansion is found using the regression method. The Sobol's quasi‐random sequence is used to generate the sample points. The numerical results of the three test cases show that the non‐intrusive model reduction scheme is able to produce satisfactory results for the statistical quantities of interest. It is found that the developed non‐intrusive model reduction scheme is computationally more efficient than the classical polynomial chaos expansion for uncertainty quantification of stochastic problems. The performance of the developed scheme becomes more apparent for the problems with larger stochastic dimensions and those requiring higher polynomial order for the stochastic discretization. Copyright © 2015 John Wiley & Sons, Ltd.     

Error estimation in a stochastic finite element method in electrokinetics

Input data to a numerical model are not necessarily well known. Uncertainties may exist both in material properties and in the geometry of the device. They can be due, for instance, to ageing or imperfections in the manufacturing process. Input data can be modelled as random variables leading to a stochastic model. In electromagnetism, this leads to solution of a stochastic partial differential equation system. The solution can be approximated by a linear combination of basis functions rising from the tensorial product of the basis functions used to discretize the space (nodal shape function for example) and basis functions used to discretize the random dimension (a polynomial chaos expansion for example). Some methods (SSFEM, collocation) have been proposed in the literature to calculate such approximation. The issue is then how to compare the different approaches in an objective way. One solution is to use an appropriate a posteriori numerical error estimator. In this paper, we present an error estimator based on the constitutive relation error in electrokinetics, which allows the calculation of the distance between an average solution and the unknown exact solution. The method of calculation of the error is detailed in this paper from two solutions that satisfy the two equilibrium equations. In an example, we compare two different approximations (Legendre and Hermite polynomial chaos expansions) for the random dimension using the proposed error estimator. In addition, we show how to choose the appropriate order for the polynomial chaos expansion for the proposed error estimator. Copyright © 2009 John Wiley & Sons, Ltd.     

Reduced chaos expansions with random coefficientsin reduced‐dimensional stochastic modeling of coupled problems

We address the curse of dimensionality in methods for solving stochastic coupled problems with an emphasis on stochastic expansion methods such as those involving polynomial chaos expansions. The proposed method entails a partitioned iterative solution algorithm that relies on a reduced‐dimensional representation of information exchanged between subproblems to allow each subproblem to be solved within its own stochastic dimension while interacting with a reduced projection of the other subproblems. The proposed method extends previous work by the authors by introducing a reduced chaos expansion with random coefficients. The representation of the exchanged information by using this reduced chaos expansion with random coefficients enables an expeditious construction of doubly stochastic polynomial chaos expansions that separate the effect of uncertainty local to a subproblem from the effect of statistically independent uncertainty coming from other subproblems through the coupling. After laying out the theoretical framework, we apply the proposed method to a multiphysics problem from nuclear engineering. Copyright © 2013 John Wiley & Sons, Ltd.     

Stochastic fracture mechanics using polynomial chaos

Crack propagation in metals has long been recognized as a stochastic process. As a consequence, crack propagation rates have been modeled as random variables or as random processes of the continuous. On the other hand, polynomial chaos is a known powerful tool to represent general second order random variables or processes. Hence, it is natural to use polynomial chaos to represent random crack propagation data: nevertheless, no such application has been found in the published literature. In the present article, the large replicate experimental results of Virkler et al. and Ghonem and Dore are used to illustrate how polynomial chaos can be used to obtain accurate representations of random crack propagation data. Hermite polynomials indexed in stationary Gaussian stochastic processes are used to represent the logarithm of crack propagation rates as a function of the logarithm of stress intensity factor ranges. As a result, crack propagation rates become log-normally distributed, as observed from experimental data. The Karhunen–Loève expansion is used to represent the Gaussian process in the polynomial chaos basis. The analytical polynomial chaos representations derived herein are shown to be very accurate, and can be employed in predicting the reliability of structural components subject to fatigue.     

Identification of composite uncertain material parameters from experimental modal data

Stochastic analysis of structures using probability methods requires the statistical knowledge of uncertain material parameters. This is often quite easier to identify these statistics indirectly from structure response by solving an inverse stochastic problem. In this paper, a robust and efficient inverse stochastic method based on the non-sampling generalized polynomial chaos method is presented for identifying uncertain elastic parameters from experimental modal data. A data set on natural frequencies is collected from experimental modal analysis for sample orthotropic plates. The Pearson model is used to identify the distribution functions of the measured natural frequencies. This realization is then employed to construct the random orthogonal basis for each vibration mode. The uncertain parameters are represented by polynomial chaos expansions with unknown coefficients and the same random orthogonal basis as the vibration modes. The coefficients are identified via a stochastic inverse problem. The results show good agreement with experimental data.     

Variance-based simplex stochastic collocation with model order reduction for high-dimensional systems

In this work, an adaptive simplex stochastic collocation method is introduced in which sample refinement is informed by variability in the solution of the system. The proposed method is based on the concept of multi-element stochastic collocation methods and is capable of dealing with very high-dimensional models whose solutions are expressed as a vector, a matrix, or a tensor. The method leverages random samples to create a multi-element polynomial chaos surrogate model that incorporates local anisotropy in the refinement, informed by the variance of the estimated solution. This feature makes it beneficial for strongly nonlinear and/or discontinuous problems with correlated non-Gaussian uncertainties. To solve large systems, a reduced-order model (ROM) of the high-dimensional response is identified using singular value decomposition (higher-order SVD for matrix/tensor solutions) and polynomial chaos is used to interpolate the ROM. The method is applied to several stochastic systems of varying type of response (scalar/vector/matrix) and it shows considerable improvement in performance compared to existing simplex stochastic collocation methods and adaptive sparse grid collocation methods.     

Stochastic convergence acceleration through basis enrichment of polynomial chaos expansions

Given their mathematical structure, methods for computational stochastic analysis based on orthogonal approximations and projection schemes are well positioned to draw on developments from deterministic approximation theory. This is demonstrated in the present paper by extending basis enrichment from deterministic analysis to stochastic procedures involving the polynomial chaos decomposition. This enrichment is observed to have a significant effect on the efficiency and performance of these stochastic approximations in the presence of non‐continuous dependence of the solution on the stochastic parameters. In particular, given the polynomial structure of these approximations, the severe degradation in performance observed in the neighbourhood of such discontinuities is effectively mitigated. An enrichment of the polynomial chaos decomposition is proposed in this paper that can capture the behaviour of such non‐smooth functions by integrating a priori knowledge about their behaviour. The proposed enrichment scheme is applied to a random eigenvalue problem where the smoothness of the functional dependence between the random eigenvalues and the random system parameters is controlled by the spacing between the eigenvalues. It is observed that through judicious selection of enrichment functions, the spectrum of such a random system can be more efficiently characterized, even for systems with closely spaced eigenvalues. Copyright © 2007 John Wiley & Sons, Ltd.     

Strain and stress computations in stochastic finite element methods

This paper focuses on the computation of statistical moments of strains and stresses in a random system model where uncertainty is modeled by a stochastic finite element method based on the polynomial chaos expansion. It identifies the cases where this objective can be achieved by analytical means using the orthogonality property of the chaos polynomials and those where it requires a numerical integration technique. To this effect, the applicability and efficiency of several numerical integration schemes are considered. These include the Gauss–Hermite quadrature with the direct tensor product—also known as the Kronecker product—Smolyak's approximation of such a tensor product, Monte Carlo sampling, and the Latin Hypercube sampling method. An algorithm for reducing the dimensionality of integration under a direct tensor product is also explored for optimizing the computational cost and complexity. The convergence rate and algorithmic complexity of all of these methods are discussed and illustrated with the non‐deterministic linear stress analysis of a plate. Copyright © 2007 John Wiley & Sons, Ltd.     

Adaptive polynomial chaos expansions applied to statistics of extremes in nonlinear random vibration

This paper presents a new module towards the development of efficient computational stochastic mechanics. Specifically, the possibility of an adaptive polynomial chaos expansion is investigated. Adaptivity in this context refers to retaining, through an iterative procedure, only those terms in a representation of the solution process that are significant to the numerical evaluation of the solution. The technique can be applied to the calculation of statistics of extremes for nongaussian processes. The only assumption involved is that these processes be the response of a nonlinear oscillator excited by a general stochastic process. The proposed technique is an extension of a technique developed by the second author for the solution of general nonlinear random vibration problems. Accordingly, the response process is represented using its Karhunen-Loeve expansion. This expansion allows for the optimal encapsulation of the information contained in the stochastic process into a set of discrete random variables. The response process is then expanded using the polynomial chaos basis, which is a complete orthogonal set in the space of second-order random variables. The time dependent coefficients in this expansion are then computed by using a Galerkin projection scheme which minimizes the approximation error involved in using a finite-dimensional subspace. These coefficients completely characterize the solution process, and the accuracy of the approximation can be assessed by comparing the contribution of successive coefficients. A significant contribution of this paper is the development and implimentation of adaptive schemes for the polynomial chaos expansion. These schemes permit the inclusion of only those terms in the expansion that have a significant contribution.     

Collocation-based stochastic finite element analysis for random field problems

A stochastic response surface method (SRSM) which has been previously proposed for problems dealing only with random variables is extended in this paper for problems in which physical properties exhibit spatial random variation and may be modeled as random fields. The formalism of the extended SRSM is similar to the spectral stochastic finite element method (SSFEM) in the sense that both of them utilize Karhunen–Loeve (K–L) expansion to represent the input, and polynomial chaos expansion to represent the output. However, the coefficients in the polynomial chaos expansion are calculated using a probabilistic collocation approach in SRSM. This strategy helps us to decouple the finite element and stochastic computations, and the finite element code can be treated as a black box, as in the case of a commercial code. The collocation-based SRSM approach is compared in this paper with an existing analytical SSFEM approach, which uses a Galerkin-based weighted residual formulation, and with a black-box SSFEM approach, which uses Latin Hypercube sampling for the design of experiments. Numerical examples are used to illustrate the features of the extended SRSM and to compare its efficiency and accuracy with the existing analytical and black-box versions of SSFEM.     

Robust topology optimization of phononic crystals with random field uncertainty

The uncertain spatial variation of material properties can remarkably affect the band gap characteristics of phononic crystals (PnCs). It is necessary to consider this issue when designing and manufacturing PnC materials/structures. This paper investigates a robust topology optimization method for designing the microstructures of PnCs by considering random‐field material properties. Herein, the spatial distribution of the material properties is first represented by a random field and then discretized into uncorrelated stochastic variables with the expansion optimal linear estimation method; stochastic band gap analysis is then conducted with polynomial chaos expansion. Furthermore, a robust topology optimization formulation of PnCs is proposed on the basis of the relative elemental density, where a weighted objective function handles the compromise of the mean value and standard deviation of the PnC band gap. The band gap response is analyzed, employing the finite element method for each sample of polynomial chaos expansion. In this context, the sensitivities of the stochastic band gap behaviors to the design variables are also derived. Numerical examples demonstrate that the proposed method can generate meaningful optimal topologies of PnCs with a relatively large width and less sensitive band gap. Additionally, the effects of the weight factors in the objective function and the variation coefficient of material properties are discussed.     

Hybridization of stochastic reduced basis methods with polynomial chaos expansions

We propose a hybrid formulation combining stochastic reduced basis methods with polynomial chaos expansions for solving linear random algebraic equations arising from discretization of stochastic partial differential equations. Our objective is to generalize stochastic reduced basis projection schemes to non-Gaussian uncertainty models and simplify the implementation of higher-order approximations. We employ basis vectors spanning the preconditioned stochastic Krylov subspace to represent the solution process. In the present formulation, the polynomial chaos decomposition technique is used to represent the stochastic basis vectors in terms of multidimensional Hermite polynomials. The Galerkin projection scheme is then employed to compute the undetermined coefficients in the reduced basis approximation. We present numerical studies on a linear structural problem where the Youngs modulus is represented using Gaussian as well as lognormal models to illustrate the performance of the hybrid stochastic reduced basis projection scheme. Comparison studies with the spectral stochastic finite element method suggest that the proposed hybrid formulation gives results of comparable accuracy at a lower computational cost.     

An iterative polynomial chaos approach toward stochastic elastostatic structural analysis with non-Gaussian randomness

Stochastic analysis of structure with non-Gaussian material property and loading in the framework of polynomial chaos (PC) is considered. A new approach for the solution of stochastic mechanics problem with random coefficient is presented. The major focus of the method is to consider reduced size of expansion in an iterative manner to overcome the problem of large system matrix in conventional PC expansion. The iterative method is based on orthogonal expansion of stochastic responses and generation of an iterative PC based on the responses of the previous iteration. The polynomials are evaluated using Gram-Schmidt orthogonalization process. The numbers of random variables in PC expansion are reduced by considering only the dominant components of the response characteristics, which is evaluated using Karhunen-Loève (KL) expansion. In case of random material field problem, the KL expansion is used to discretize and simulate the non-Gaussian random field. Independent component analysis (ICA) is carried out on the non-Gaussian KL random variables to minimize statistical dependence. The usefulness of the proposed method in terms of accuracy and computational efficiency is examined. From the numerical analysis of three different types of structural mechanics problems, the proposed iterative method is observed to be computationally more efficient and accurate than conventional PC method for solution of linear elastostatic structural mechanics problems.     

Stochastic projection schemes for deterministic linear elliptic partial differential equations on random domains

We present stochastic projection schemes for approximating the solution of a class of deterministic linear elliptic partial differential equations defined on random domains. The key idea is to carry out spatial discretization using a combination of finite element methods and stochastic mesh representations. We prove a result to establish the conditions that the input uncertainty model must satisfy to ensure the validity of the stochastic mesh representation and hence the well posedness of the problem. Finite element spatial discretization of the governing equations using a stochastic mesh representation results in a linear random algebraic system of equations in a polynomial chaos basis whose coefficients of expansion can be non‐intrusively computed either at the element or the global level. The resulting randomly parametrized algebraic equations are solved using stochastic projection schemes to approximate the response statistics. The proposed approach is demonstrated for modeling diffusion in a square domain with a rough wall and heat transfer analysis of a three‐dimensional gas turbine blade model with uncertainty in the cooling core geometry. The numerical results are compared against Monte–Carlo simulations, and it is shown that the proposed approach provides high‐quality approximations for the first two statistical moments at modest computational effort. Copyright © 2010 John Wiley & Sons, Ltd.     

基于同伦分析方法的随机结构静力响应求解

该文提出了一种基于同伦分析方法的求解含随机参数结构的静力响应的新方法。该方法将随机静力平衡方程重新进行同伦构造,利用含随机变量和趋近函数的同伦级数展式来表示结构的随机静力位移响应,该同伦级数的各阶确定性系数和趋近函数可通过对一系列的变形方程求解得到。由于趋近函数的引入,该同伦级数解相较于传统的摄动法有更大的收敛范围,对于含较大变异性随机参数的结构也能获得不错的求解精度。同时,该文提出了一种降维策略来提高该方法的计算效率。数值算例表明,与目前广泛应用的广义正交多项式展开法（GPC）相比,从计算精度上看,该文方法的3阶展开与GPC 2阶展开相当,该文方法的6阶展开与GPC 4阶展开相当,而计算时间上前者均明显少于后者。此外,该文方法也可以方便地应用到随机结构的几何非线性分析当中,并具有较好的计算精度和计算效率。     

Time response of structures with uncertain parameters under stochastic inputs based on coupled polynomial chaos expansion and component mode synthesis methods

This article presents a numerical procedure to compute the stochastic dynamic response of large finite element models with uncertain parameters based on polynomial chaos and component mode synthesis methods. Polynomial chaos expansions at higher orders are used to derive the statistical solution of the dynamic response as well as the Monte Carlo simulation procedure. Based on various component mode synthesis methods, the size of the model is reduced. These methods are coupled with polynomial chaos expansion and the explicit mathematical formulations are given. Numerical results illustrating the accuracy and efficiency of the proposed coupled methodological procedures are presented.     

A reduced polynomial chaos expansion method for the stochastic finite element analysis

The stochastic finite element analysis of elliptic type partial differential equations is considered. A reduced method of the spectral stochastic finite element method using polynomial chaos is proposed. The method is based on the spectral decomposition of the deterministic system matrix. The reduction is achieved by retaining only the dominant eigenvalues and eigenvectors. The response of the reduced system is expanded as a series of Hermite polynomials, and a Galerkin error minimization approach is applied to obtain the deterministic coefficients of the expansion. The moments and probability density function of the solution are obtained by a process similar to the classical spectral stochastic finite element method. The method is illustrated using three carefully selected numerical examples, namely, bending of a stochastic beam, flow through porous media with stochastic permeability and transverse bending of a plate with stochastic properties. The results obtained from the proposed method are compared with classical polynomial chaos and direct Monte Carlo simulation results.     

Hierarchical surrogate model with dimensionality reduction technique for high-dimensional uncertainty propagation

In this article, hierarchical surrogate model combined with dimensionality reduction technique is investigated for uncertainty propagation of high-dimensional problems. In the proposed method, a low-fidelity sparse polynomial chaos expansion model is first constructed to capture the global trend of model response and exploit a low-dimensional active subspace (AS). Then a high-fidelity (HF) stochastic Kriging model is built on the reduced space by mapping the original high-dimensional input onto the identified AS. The effective dimensionality of the AS is estimated by maximum likelihood estimation technique. Finally, an accurate HF surrogate model is obtained for uncertainty propagation of high-dimensional stochastic problems. The proposed method is validated by two challenging high-dimensional stochastic examples, and the results demonstrate that our method is effective for high-dimensional uncertainty propagation.     

An adaptive algorithm to build up sparse polynomial chaos expansions for stochastic finite element analysis

Polynomial chaos (PC) expansions are used in stochastic finite element analysis to represent the random model response by a set of coefficients in a suitable (so-called polynomial chaos) basis. The number of terms to be computed grows dramatically with the size of the input random vector, which makes the computational cost of classical solution schemes (may it be intrusive (i.e.of Galerkin type) or non-intrusive) unaffordable when the deterministic finite element model is expensive to evaluate.