Application of a Stochastic Finite Element Procedure to Reliability Analysis

Transcription

1 Application of a Stochastic Finite Element Procedure to Reliability Analysis Bruno Sudret & Marc Berveiller Electricité de France, R&D Division, Site des Renardières, F Moret-sur-Loing Maurice Lemaire Institut Français de Mécanique Avancée, LaRAMA, Campus des Cézeaux, BP65, F Aubière Keywords: Stochastic Finite Element, Polynomial Chaos, Structural Reliability, FORM/SORM, Parametric Study ABSTRACT: The stochastic finite element procedure presented in this paper consists in representing in a probabilistic form the response of a system whose material properties and loading are random Each random variable is expanded into Hermite polynomial series The response is expanded onto the so-called polynomial chaos The coefficients of the expansion of the response (eg the nodal displacement vector) are obtained by a Galerkin-type method in the space of probability This stochastic finite element procedure is applied to the parametric study of the reliability of a foundation on a two-layer elastic soil mass The method reveals accurate and more efficient than a direct coupling between a finite element code and a probabilistic code 1 INTRODUCTION The Spectral Stochastic Finite Element Method (SSFEM) was introduced in the early 90 s for solving stochastic boundary value problems in which the spatial variability of a material property (eg Young s modulus) was modeled as a random field (Ghanem and Spanos, 1991) The approach, originally developed for a linear elastic problem using Gaussian random fields, was later extended to heat transfer problems (Ghanem, 1999b; Xiu and Karniadakis, 00) and lognormal fields (Ghanem, 1999b) Recently, applications in fluid mechanics (Le Maître et al, 001), modal analysis (Van Den Nieuwenhof, 003), soil structure interaction (Ghiocel and Ghanem, 00; Clouteau and Lafargue, 003) have also been proposed In most papers, the end result is the calculation of the two first moments of the response (ie mean and standard deviation) Application to structural reliability as a postprocessing of SSFEM has been demonstrated in Sudret and Der Kiureghian (000); Sudret and Der Kiureghian (00) In all the above references, the basic assumption is the description of the system randomness in terms of one or two random fields, which reduces the spectrum of applications The aim of this paper is to develop a stochastic finite element procedure for problems involving any number of random variables of any type (ie non necessarily Gaussian and lognormal) describing uncertainties in the material properties as well as the loading The proposed stochastic finite element method (SFEM) relies upon three ingredients : 1

2 expansion of the input random variables in terms of Hermite polynomials of standard normal variables, expansion of the response onto the so-called polynomial chaos, post-processing of the results for reliability analysis The method is applied to study the reliability of a foundation on a two-layer elastic soil mass Comparisons with a classical finite element reliability analysis (direct coupling between the FORM method and a finite element code) are given in terms of accuracy and efficiency EXPANSION OF RANDOM VARIABLES ONTO THE POLYNOMIAL CHAOS Let us denote by L (Θ,F,P) the Hilbert space of random variables with finite variance Let us consider a random variable X with prescribed probability density function f X (x) Classical results Malliavin (1997) allow to expand X as a Hermite polynomial series expansion : X = a i H i (ξ) (1) where ξ denotes a standard normal variable, {H i,i = 0,, } are Hermite polynomials defined by : ( ) H i (x) = ( 1) i e x d dx i e x () and {a i, i = 0,, } are coefficients to be evaluated Two methods are now presented for this purpose 1 Projection method This method was used by Puig et al (00); Xiu and Karniadakis (00) Due to the orthogonality of the Hermite polynomials with respect to the Gaussian measure, it comes from Eq(1) : E[XH i (ξ)] = a i E[Hi (ξ)] (3) where E[H i (ξ)] = i! Let F X denote the cumulative distribution function (CDF) of the random variable to approximate and Φ the standard normal CDF By using the transformation to the standard normal space X ξ : F X (X) = Φ(ξ), one can write : Thus : a i = 1 i! E[X(ξ)H i(ξ)] = 1 i! X(ξ) = F 1 X (Φ(ξ)) (4) R F 1 X (Φ(t))H i(t)ϕ(t)dt (5) where ϕ is the standard normal probability density function (PDF) When X is a normal or lognormal random variable, coefficients {a i,i = 0,, } can be evaluated analytically : X N(µ,σ) a 0 = µ, a 1 = σ, a i = 0 for i X LN(λ,ζ) a i = ζi i! exp[λ + 1 ζ ] for i 0 (6) For other types of distribution, the quadrature method may be used for evaluating the integral in Eq(5) (Berveiller and Sudret, 003)

3 Collocation method This method was introduced by Isukapalli (1999) It is based on a least square minimization of the discrepancy between the input variable X and its truncated approximation X : X = p a i H i (ξ) (7) Let {ξ (1),,ξ (n) } be n outcomes of ξ From Eq(4), we obtain n outcomes {X (1),,X (n) } The least square method consists in minimizing the following quantity with respect to {a i,i = 0,,p} : n n p X = (X (i) X (i) ) = X (i) a j H j (ξ (i) ) (8) This leads to the following linear system yielding the expansion coefficients {a i,i = 0,,p} : n n n H 0 (ξ (i) )H 0 (ξ (i) ) H 0 (ξ (i) )H p (ξ (i) ) X (i) H 0 (ξ (i) ) a 0 = n n H p (ξ (i) )H 0 (ξ (i) ) H p (ξ (i) )H p (ξ (i) ) a p n X (i) H p (ξ (i) ) (9) Both methods are illustrated in Figure 1 in case of a lognormal distribution LN(06501, 0936) with mean and standard deviation theoretical projection collocation Method a 0 a 1 a a 3 Projection Collocation FIG 1: Theoretical and approximated PDF for a lognormal distribution LN(06501, 0936) - Coefficients of the expansion (p = 3) 3 STOCHASTIC FINITE ELEMENT METHOD IN LINEAR ELASTICITY Using classical notations (Zienkiewicz and Taylor, 000), the finite element method for static problems in linear elasticity yields a linear system of size N ddl N ddl where N ddl denotes the number of degrees of freedom of the structure : K U = F (10) where K is the global stiffness matrix, U is the basic response quantity (vector of nodal displacements) and F is the vector of nodal forces 3

4 In SFEM, due to the introduction of input random variables, the basic response quantity is a random vector of nodal displacements U(θ) Each component is a random variable expanded onto the so-called polynomial chaos : U(θ) = U j Ψ j ({ξ k (θ)} M k=1 ) (11) where {ξ k (θ)} M k=1 denotes the set of standard normal variables appearing in the discretization of all input random variables and {Ψ j ({ξ k } M k=1 )} are multidimensional Hermite polynomials that form an orthogonal basis of L (Θ,F,P) As the polynomial chaos contains the one-dimensional Hermite polynomials in each random variable, the expressions of each input random variable may be injected in the polynomial chaos expression : X i = x i jh j (ξ i ) x i jψ j ({ξ k } M k=1 ) (1) In the sequel, the dependency of Ψ j in {ξ k } M k=1 will be omitted for the sake of clarity 31 Taking into account randomness in material properties In the deterministic case, the global stiffness matrix reads : K = k e = B T D BdΩ e (13) e e Ω e where B is the matrix that relates the components of strain to the element nodal displacement, D is the elasticity matrix and is the assembly procedure over all elements When material properties e are described by means of random variables, the elasticity matrix hence the global stiffness matrix become random The latter may be expanded onto the polynomial chaos as follows : where K = K j = E[KΨ j ] = e K j Ψ j (14) Ω e B T E[DΨ j ] BdΩ e (15) Note that B is a deterministic matrix while D is random In case of an isotropic elastic material with random independent Young s modulus E and Poisson s ratio ν, the latter may be written as (plane strain problems) : D = E( λ(ν)d 1 + µ(ν)d ) (16) where and D 1 = λ(ν) = ν (1 + ν)(1 ν) , µ(ν) =, D = 1 (1 + ν) (17) (18) 4

5 Random Young s modulus E is expanded as in Eq(1) Functions of the random Poisson s ratio { λ(ν), µ(ν)} may be expanded in the same way, using either the projection or the collocation method : E = e i H i (ξ E ) λ(ν) = λ j H j (ξ ν ) µ(ν) = µ j H j (ξ ν ) (19) Note that the same standard normal variable ξ ν is used to expand both functions λ(ν) and µ(ν) By substituting for Eq(19) in (16), one finally gets : D = e i λ j H i (ξ E )H j (ξ ν )D 1 + e i µ j H i (ξ E )H j (ξ ν )D (0) Products H i (ξ E )H j (ξ ν ) may be injected into the polynomial chaos Ψ k (ξ E,ξ ν ), finally yielding : D = (α k D 1 + β k D )Ψ k (ξ E,ξ ν ) (1) k=0 If the structure under consideration is made of several materials, the above procedure is applied using different random variables in each element group having the same material properties 3 Taking into account randomness in loading The vector of nodal forces may be written : N q F = q i F i () where N q is the number of load cases, {q i } Nq denote random loading parameters and F i load pattern vectors corresponding to a unit value of q i Note that this formulation equally applies to pinpoint forces, pressure or initial stresses Coefficients q i can be expanded onto the polynomial chaos (see Eq(1)) : q i = qj i H j(ξ i ) = q j i Ψ j (3) Finally the random vector of nodal forces reads : F = N q q j i Ψ jf i = F j Ψ j (4) 33 Global linear system By using (11), (14) and (4), the discretized stochastic equilibrium equation reads : ( ) K i Ψ i U j Ψ j = F j Ψ j (5) After a truncature of the series appearing in Eq(5) at order P, the residual in the equilibrium equation is : ( P 1 ) P 1 P 1 ɛ P = K i Ψ i U j Ψ j F j Ψ j (6) 5

6 Coefficients {U 0,,U P 1 } are computed using the Galerkin method minimizing the residual defined above, which is equivalent to requiring that this residual be orthogonal to the space spanned by {Ψ j } P 1 (Ghanem and Spanos, 1991) : E[ɛ P Ψ k ] = 0, k = {0,,P 1} (7) This leads to the linear system : K 0,0 K 0,P 1 K 1,0 K 1,P 1 U 0 U 1 = F 0 F 1 (8) K P 1,0 K P 1,P 1 U P 1 F P 1 P 1 where K j,k = K i d ijk and : i!j!k! ( ) ( ) ( ), if (i + j + k) even and k [ i j,i + j] d ijk = E[Ψ i Ψ j Ψ k ] = i+j k! j+k i! k+i j! 0 otherwise (9) This system can be solve directly (Ghanem and Kruger, 1996) or using a hierarchical approach (Ghanem, 1999a, 000) after a direct resolution at a lower order 4 APPLICATION TO RELIABILITY ANALYSIS Solving the linear system (8) yields the expansion coefficients of the vector of nodal displacements U = U j Ψ j It is easy to show that any response quantity (eg strain or stress P 1 component) may be also expanded onto the polynomial chaos (Berveiller and Sudret, 003) Thus the mechanical response of the system S (ie the set of all nodal displacements, strain or stress components) may be P 1 written as S = S j Ψ j In reliability analysis, the failure criterion of a structure is defined in terms of a limit state function g(x, S(X)) which may depend both on basic random variables X and response quantities S(X) After performing a SFEM analysis as described above, it is clear that any limit state function is analytical and defined in terms of standard normal variables : P 1 g(x,s(x)) g {ξ k } M k=1, S j Ψ j ({ξ k } M k=1 ) (30) Thus the reliability problem, which is already formulated in the standard normal space, may be solved by any available method including Monte-Carlo Simulation, FORM/SORM, Importance Sampling, etc (Ditlevsen and Madsen, 1996) Application of this approach (SFEM+FORM) in the context of spatial variability in an elastic mechanical problem can be found in Sudret and Der Kiureghian (00) Application to one (resp two)-layer elastic foundation with random elastic properties and loading can be found in Sudret et al (003) Both references show the accuracy of the approach for reasonable 6

7 reliability indices up to β = 4 5 (less than 5% discrepancy in β with respect to the reference solution) In this paper we further develop the approach to carry out efficient parametric studies Suppose the reliability of an elastic structure is investigated as a function of a load parameter λ, for instance : g(x,λ) = S λs(x) (31) where S is an admissible threshold and S is a response quantity (eg nodal displacement) depending on the basic random variable X and the load parameter λ In order to get the probability of failure with respect to λ (fragility curve), two strategies may be adopted : for each value of λ, a finite element reliability analysis is carried out by coupling a finite element code and a reliability code (Der Kiureghian and Ke, 1988; Lemaire and Mohamed, 000) A parametric study with respect to λ may be carried out, which enhances the efficiency in terms of numbers of calls to the finite element code a single SFEM analysis corresponding to λ = 1 is performed The limit state function (31) rewrites : P g SFEM (X,λ) = S 1 λ S j Ψ j ({ξ k } M k=1 ) (3) The successive reliability analysis corresponding to Eq(3) are then almost costless since this limit state function is analytical Both approach are now compared on an application example in geotechnical engineering 5 APPLICATION EXAMPLE : SERVICEABILITY OF A FOUNDATION Let us consider an elastic soil mass made of two layers of different isotropic linear elastic materials lying on a rigid substratum A foundation on this soil mass is modeled as a uniform pressure λ applied over a length B of the free surface (figure 5) Due to the symmetry, half of the structure is modeled by finite element The mesh of half of the foundation (due to symmetry) comprises 99 nodes and 80 4-node elements which allows a 1%-accurate evaluation of the maximal settlement compared to a reference solution The model parameters are listed in Table 1 There are four random variables, namely the Young s modulus and Poisson s ratio of each layer The reliability of the foundation with respect to the maximum admissible settlement ū is investigated Hence the parametric limit state function : g(x,λ) = ū λu A (E 1,ν 1,E,ν ) (33) where u A (E 1,ν 1,E,ν ) is the maximal displacement obtained for a unit load λ = 1 The reliability index is computed by FORM for different values of λ using the two strategies mentioned in Section 4, namely : a direct coupling between the finite element code Code AsterR [http ://wwwcode-asterorg] and the probabilistic code PROBAN (parametric study) a single SFEM analysis followed by a parametric FORM reliability analysis from Eq(33) These tools are implemented in a Matlab package Results are reported in Table 3 and plotted in Figure 3 Column # corresponds to the direct coupling reference solution, columns #3-5 corresponds to different strategies of resolution in the SFEM analysis, namely a complete resolution at order or 3 and a hierarchical resolution at order 3 with pre-resolution at order (denoted by p = 3()) It appears that for reliability indices lower than 5, we have excellent results for the three different methods used For reliability indices greater than 5, the direct method with an order gives poor results, the two other results are equivalent More 7

8 FIG : Mesh of the foundation TAB 1: Parameters of the foundation Parameter Notation Type Mean Coef of Var Upper layer Young s Modulus E 1 Lognormal 50 MPa 0% Lower layer Young s Modulus E Lognormal 70 MPa 0% Upper layer Poisson s Ratio ν 1 Beta [0,05] 03 0% Lower layer Poisson s Ratio ν Beta [0,05] 03 0% Width B Deterministic 10 m - Mesh width L Deterministic 60 m - Upper layer soil thickness t 1 Deterministic 775 m - Lower layer soil thickness t Deterministic 5 m - Threshold ū Deterministic 01 m - generally, other results not reported in this paper show that the hierarchical resolution at order 3 with a pre-resolution at order is the best compromise between accuracy and efficiency Table reports the computer processing time for the complete parametric study in each case, where the unit time corresponds to one single deterministic analysis (ie p = 0 in the SFEM context) It appears that the direct coupling is about as computationally expansive as the SFEM with p = 3 However the same accuracy is obtained for the SFEM (p = 3()) hierarchical solution at about one tenth of the cost TAB : Computer processing time required by the direct coupling and SFEM at various orders Direct Coupling SFEM p = 0 Coupling p = 0 p = p = 3 p = 3() Time direct p= direct p=3 hierarchical p=3() Coupling direct p= direct p=3 hierarchical p=3() Coupling FIG 3: Left : Reliability index β vs load P - Right : Failure probability P f vs P (fragility curve) 8

9 TAB 3: Reliability index for the settlement of a foundation with different loads λ λ (MPa) Direct Coupling SFEM β DirectCoupling β SFEMp=3 p = p = 3 p = 3() CONCLUSION This paper presented a stochastic finite element method able to solve linear mechanical problems with random materials properties and loading Basic random variables can be in any number and of any type The paper focuses on the application of SFEM to solving parametric reliability problems The approach appears accurate for reliability indices up to β = 5, at a cost which is one order of magnitude smaller than the direct coupling in the application example Note that the SFEM also provides straightforwardly the statistical moments of the response (Sudret et al, 003) and that any other reliability analysis (eg corresponding to other limit state functions) is costless since the complete probabilistic content of the response has been characterized once for all in Eq(11) Work is currently in progress to obtain the expansion coefficients in Eq(11) in a different way, which will open SFEM to non linear problems Références Berveiller, M and B Sudret (003) Thèse sur les éléments finis stochastiques Rapport d avancement 1 Technical Report HT-6/03/039/A, EDF R&D internal report (in French) Clouteau, D and R Lafargue (003) An iterative solver for stochastic soil-structure interaction In P Spanos and G Deodatis (Eds), Computational Stochastic Mechanics (CSM-4), pp Der Kiureghian, A and J-B Ke (1988) The stochastic finite element method in structural reliability Prob Eng Mech 3(), Ditlevsen, O and H Madsen (1996) Structural reliability methods J Wiley and Sons, Chichester Ghanem, R-G (1999a) Ingredients for a general purpose stochastic finite elements implementation Comp Meth Appl Mech Eng 168, Ghanem, R-G (1999b) Stochastic finite elements with multiple random non-gaussian properties J Eng Mech, ASCE 15(1), 6 40 Ghanem, R-G (000) Iterative solution of systems of linear equations arising in the context of stochastic finite elements Advances in Engineering Software 31, Ghanem, R-G and R Kruger (1996) Numerical solution of spectral stochastic finite element systems Comp Meth Appl Mech Eng 19(3),

10 Ghanem, R-G and P-D Spanos (1991) Stochastic finite elements - A spectral approach Springer Verlag Ghiocel, D and R Ghanem (00) Stochastic finite element analysis of seismic soil-structure interaction J Eng Mech 18, Isukapalli, S S (1999) Uncertainty Analysis of Transport-Transformation Models Ph D thesis, The State University of New Jersey Le Maître, O, O M Knio, H N Najm, and R Ghanem (001) A stochastic projection method for fluid flow J Comp Phys 173, Lemaire, M and A Mohamed (000) Finite element and reliability : A happy marriage? In Reliability and Optimization of Structural Systems, Keynote Lecture at the 9th IFIP WG 75 Working Conference, pp 3 14 Nowak, A and Szerszen, M (Eds) Malliavin, P (1997) Stochastic Analysis Springer Puig, B, F Poirion, and C Soize (00) Non-Gaussian simulation using Hermite polynomial expansion : convergences Prob Eng Mech 17, Sudret, B, M Berveiller, and M Lemaire (003) Eléments finis stochastiques spectraux, nouvelles perspectives In Proc 16ème Congrès Français de Mécanique - Nice (in French) Sudret, B and A Der Kiureghian (000) Stochastic finite elements and reliability : A state-of-the-art report Technical Report n o UCB/SEMM-000/08, University of California, Berkeley 173 pages Sudret, B and A Der Kiureghian (00) Comparison of finite element reliability methods Prob Eng Mech 17, Van Den Nieuwenhof, B (003) Stochastic Finite Elements for elastodynamics : random fields and shape uncertainty modeling using direct and modal perturbation-based aproaches Ph D thesis, Université Catholique de Louvain Xiu, D and G E Karniadakis (00) The Wiener-Askey polynomial chaos for stochastic differentiel equations J Sci Comput 4(), Zienkiewicz, O-C and R-L Taylor (000) The finite element method Butterworth Heinemann, 5 th edition 10