Strategies for finding the design point in non-linear finite element reliability analysis

• Published in: Probabilistic engineering mechanics Vol. 21; no. 2; pp. 133 - 147
• Main Authors: Haukaas, Terje, Kiureghian, Armen Der
• Format: Journal Article
• Language: English
• Published: Oxford Elsevier Ltd 01-04-2006; Elsevier Science
• Subjects: Exact sciences and technology; Finite elements; Fracture mechanics (crack, fatigue, damage...); Fundamental areas of phenomenology (including applications); Inelasticity (thermoplasticity, viscoplasticity...); Non-linear analysis; Optimization algorithms; Physics; Response sensitivity; Smooth material models; Solid mechanics; Structural and continuum mechanics; Structural reliability; Vibration, mechanical wave, dynamic stability (aeroelasticity, vibration control...); Structural reliability; Optimization algorithms; Finite elements; Smooth material models; Response sensitivity; Non-linear analysis; Discontinuity; Probabilistic approach; Rupture; Iterative method; Modeling; Search algorithm; Constrained optimization; Finite element method; Inelasticity; Plasticity; Non linear effect; Numerical convergence; Importance sampling
• ISSN: 0266-8920; 1878-4275
• DOI: 10.1016/j.probengmech.2005.07.005

An essential step in FORM, SORM and importance sampling reliability methods is the determination of the so-called design point. This point is the solution of a constrained optimization problem in the outcome space of the random variables, which is commonly solved by an iterative, gradient-based search algorithm. In solving this problem in the context of non-linear finite element reliability analysis, two serious impediments are encountered: (a) for certain material models, the constraint function may have a discontinuous gradient, leading to failure of the search algorithm to converge. (b) The search algorithm may generate trial points too far in the failure domain, where the finite element code fails to produce a result due to lack of numerical convergence. In this paper, remedying strategies are developed for both impediments. The first impediment is addressed by using smooth or smoothed material models, including a smoothed bi-linear model, a Bouc–Wen model and a generalized plasticity model. This is complemented by a proof that sudden elastic unloading does not give rise to gradient discontinuities. The second impediment is addressed by modifying or introducing search algorithms that prevent the trial points from overshooting into the failure domain. Numerical examples are used to demonstrate the two impediments and effectiveness of the proposed remedies.