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
VL - 206
SP - 375
EP - 391
JO - Engineering Fracture Mechanics
JF - Engineering Fracture Mechanics
SN - 0013-7944
ER -