TY - JOUR

T1 - Extracting stochastic stress intensity factors using generalized polynomial chaos

AU - Omer, Netta

AU - Yosibash, Zohar

N1 - Publisher Copyright:
© 2018 Elsevier Ltd

PY - 2019/2/1

Y1 - 2019/2/1

N2 - Realistic material properties, as the Young modulus E and Poisson ratio ν (isotropic materials), are measured by experimental observations and are inherently stochastic. Having their stochastic representation E(ξ) or ν(ξ) where ξ is a random variable, we formulate the elastic solution of the stochastic elasticity system in the vicinity of a crack tip. We show that the stochastic asymptotic displacements are of the form. u(r,θ;ξ)=A01(ξ)φ(01)(θ;ξ)+A02(ξ)φ(02)(θ;ξ)+A1(ξ)r1/2φ(1)(θ;ξ)+A2(ξ)r1/2φ(2)(θ;ξ)+O(r)with deterministic eigenvalues and either deterministic or stochastic eigenfunctions φ(i)(θ;ξ) and coefficients Ai. However, the stresses are represented by an asymptotic series with a stochastic behavior manifested only in the SIF: σ(r,θ;ξ)=[Formula presented]ϕ(1)(θ)+[Formula presented]ϕ(2)(θ)+O(1)We present explicitly whether the expressions in series expansions are stochastic or not, depending both on the material properties and on the boundary conditions far from the crack faces. The generalized polynomial chaos (GPC) method is used thereafter to compute the stochastic expressions from deterministic finite element solutions. As an example we consider either a stochastic Young modulus or Poisson ratio to be given as random variable with a normal distribution: E(ξ)=E0+E1ξ,orν(ξ)=ν0+ν1ξ,ξ∼N(0,σ2)Numerical examples are presented in which we compute φi(θ;ξ) and thereafter Ai(ξ) and KI(ξ) from deterministic finite element analyses using the GPC. Monte Carlo simulations are used to demonstrate the efficiency of the proposed methods. Numerical examples are provided that show the efficiency and accuracy of the proposed methods.

AB - Realistic material properties, as the Young modulus E and Poisson ratio ν (isotropic materials), are measured by experimental observations and are inherently stochastic. Having their stochastic representation E(ξ) or ν(ξ) where ξ is a random variable, we formulate the elastic solution of the stochastic elasticity system in the vicinity of a crack tip. We show that the stochastic asymptotic displacements are of the form. u(r,θ;ξ)=A01(ξ)φ(01)(θ;ξ)+A02(ξ)φ(02)(θ;ξ)+A1(ξ)r1/2φ(1)(θ;ξ)+A2(ξ)r1/2φ(2)(θ;ξ)+O(r)with deterministic eigenvalues and either deterministic or stochastic eigenfunctions φ(i)(θ;ξ) and coefficients Ai. However, the stresses are represented by an asymptotic series with a stochastic behavior manifested only in the SIF: σ(r,θ;ξ)=[Formula presented]ϕ(1)(θ)+[Formula presented]ϕ(2)(θ)+O(1)We present explicitly whether the expressions in series expansions are stochastic or not, depending both on the material properties and on the boundary conditions far from the crack faces. The generalized polynomial chaos (GPC) method is used thereafter to compute the stochastic expressions from deterministic finite element solutions. As an example we consider either a stochastic Young modulus or Poisson ratio to be given as random variable with a normal distribution: E(ξ)=E0+E1ξ,orν(ξ)=ν0+ν1ξ,ξ∼N(0,σ2)Numerical examples are presented in which we compute φi(θ;ξ) and thereafter Ai(ξ) and KI(ξ) from deterministic finite element analyses using the GPC. Monte Carlo simulations are used to demonstrate the efficiency of the proposed methods. Numerical examples are provided that show the efficiency and accuracy of the proposed methods.

KW - Generalized polynomial chaos

KW - Stochastic stress intensity factors

UR - http://www.scopus.com/inward/record.url?scp=85058422076&partnerID=8YFLogxK

U2 - 10.1016/j.engfracmech.2018.12.001

DO - 10.1016/j.engfracmech.2018.12.001

M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???

AN - SCOPUS:85058422076

SN - 0013-7944

VL - 206

SP - 375

EP - 391

JO - Engineering Fracture Mechanics

JF - Engineering Fracture Mechanics

ER -