Stochastic meshless analysis of elastic–plastic cracked structures

This paper presents a stochastic mesh-free method for probabilistic fracture-mechanics analysis of nonlinear cracked structures. The method involves enriched element-free Galerkin formulation for calculating the J-integral; statistical models of uncertainties in load, material properties, and crack geometry; and the first-order reliability method (FORM) for predicting probabilistic fracture response and reliability of cracked structures. The sensitivity of fracture parameters with respect to crack size, required for probabilistic analysis, is calculated using a virtual crack extension technique. Numerical examples based on mode-I fracture problems have been presented to illustrate the proposed method. The results from sensitivity analysis indicate that the maximum difference between sensitivity of the J-integral calculated using the proposed method and reference solutions obtained by the finite-difference method is about three percent. The results from reliability analysis show that the probability of fracture initiation using the proposed sensitivity and meshless-based FORM are very accurate when compared with either the finite-element-based Monte Carlo simulation or finite-element-based FORM. Since all gradients are calculated analytically, the reliability analysis of cracks can be performed efficiently using meshless methods. The authors would like to acknowledge the financial support of the U.S. National Science Foundation (NSF) under Award No. CMS-9900196. The NSF program director was Dr. Ken Chong.     

Probabilistic fracture mechanics by Galerkin meshless methods – part II: reliability analysis

 This is the second in a series of two papers generated from a study on probabilistic meshless analysis of cracks. In this paper, a stochastic meshless method is presented for probabilistic fracture-mechanics analysis of linear-elastic cracked structures. The method involves an element-free Galerkin method for calculating fracture response characteristics; statistical models of uncertainties in load, material properties, and crack geometry; and the first-order reliability method for predicting probabilistic fracture response and reliability of cracked structures. The sensitivity of fracture parameters with respect to crack size, required for probabilistic analysis, is calculated using a virtual crack extension technique described in the companion paper [1]. Numerical examples based on mode-I and mixed-mode problems are presented to illustrate the proposed method. The results show that the predicted probability of fracture initiation based on the proposed formulation of the sensitivity of fracture parameter is accurate in comparison with the Monte Carlo simulation results. Since all gradients are calculated analytically, reliability analysis of cracks can be performed efficiently using meshless methods. Received 20 February 2001 / Accepted 19 December 2001     

A methodology for computing nonlinear fracture parameters for a bulging crack in a pressurised aircraft fuselage

A computational methodology for obtaining nonlinear fracture parameters which account for the effects of plasticity at the tips of a bulging crack in a pressurised aircraft fuselage is developed. The methodology involves a hierarchical three stage analysis (global, intermediate, and local) of the cracked fuselage, with the crack incorporated into the model at each stage. The global analysis is performed using a linear elastic shell finite element model in which the stiffeners are treated as beam elements. The geometrically nonlinear nature of the bulging phenomenon is emulated in the intermediate analysis using a geometrically nonlinear shell finite element model. The local analysis is a three-dimensional solid finite element model of the cracked skin using a hypoelastic-plastic rate formulation. Kinematic boundary conditions for each stage are obtained from the preceding stage in the hierarchy using a general mesh independent mechanism. The T ^*integral, which accounts for both large deformations and plasticity, is taken to be the fracture parameter characterising the severity of the conditions at the crack tip, and is evaluated from the local analysis using the Equivalent Domain Integral (EDI) method. The implementation of the EDI technique for finite deformations in shell space is also outlined. The methodology is applied to a number of example problems for which correction factors relating the nonlinear T ^* values to those obtained from a linear elastic stiffened shell analysis are computed. The issue of flapping is addressed by investigating the behaviour of the longitudinal stress parallel to the crack for various cases.This research is supported under a grant to the Center of Excellence for Computational Modeling of Aircraft Structures at Georgia Institute of Technology, from the Federal Aviation Administration. The authors also thank Dr. R. Singh, Mr. V Nagaswamy, Mr. L. Wang for many helpful discussions. The first author also wishes to express his heartfelt thanks to S. for helping him through some turbulent times of this life     

Probabilistic fracture mechanics analysis of linear-elastic cracked structures by spline fictitious boundary element method

The inherent uncertainties in crack geometry, material properties and loadings have large influence on fracture response characteristics of cracked structures. This paper presents the probabilistic fracture mechanics analysis of linear-elastic cracked structures subjected to mixed-mode loading conditions using the spline fictitious boundary element method (SFBEM). The response surface method (RSM) is used to predict the fracture probability of the cracked structure. To determine the unknown coefficients of the response surface function, the SFBEM based on the Erdogan fundamental solutions for infinite cracked plates is adopted to perform deterministic analyses of stress intensity factors (SIFs) corresponding to different test points with given parameters. Numerical examples based on mode-I and mixed-mode crack problems are presented to illustrate the present method. The results show that the predicted failure probability obtained by the present approach is accurate in comparison with the Monte Carlo simulation (MCS) results. Since a much lesser number of numerical tests are required in RSM as compared with that needed in MCS, and since the SFBEM based on the Erdogan fundamental solutions has been used to conduct the numerical tests, reliability analysis of cracked structures can be performed efficiently using the present method.     

A generalized stochastic perturbation technique for plasticity problems

The main aim of this paper is to present an algorithm and the solution to the nonlinear plasticity problems with random parameters. This methodology is based on the finite element method covering physical and geometrical nonlinearities and, on the other hand, on the generalized nth order stochastic perturbation method. The perturbation approach resulting from the Taylor series expansion with uncertain parameters is provided in two different ways: (i) via the straightforward differentiation of the initial incremental equation and (ii) using the modified response surface method. This methodology is illustrated with the analysis of the elasto-plastic plane truss with random Young’s modulus leading to the determination of the probabilistic moments by the hybrid stochastic symbolic-finite element method computations.     

Reliability analysis of a mechanical contact between deformable solids

The contact problems in mechanics are very complex. The behavior of a structure, where a contact occurs, depends on some assumptions on the design parameters and is in general highly nonlinear. Several deterministic numerical methods have been developed on finite element codes and give accurate solutions of contact problems. However some design variables could be random and the deterministic results could be unacceptable. Thus a new reliability analysis of a mechanical contact is presented. A reliability–mechanical combination based on an augmented Lagrangian method and a response surface method is proposed to compute the failure probability of this nonlinear problem. Then the advantages of this probabilistic approach are exposed and commentated on an Hertz numerical application.     

Two-parameter fracture assessment of surface cracked cylindrical shells during collapse

The present study addresses the use of CTOD and T-stress in fracture assessments of surface cracked shell structures. A new software is developed for this purpose, denoted LINK_pipe. It is based on a combination of a quadrilateral assumed natural deviatoric strain thin shell finite element and an improved linespring finite element. Plasticity is accounted for using stress resultants. A power law hardening model is used for shell and linespring materials. A co-rotational formulation is employed to represent nonlinear geometry effects. With this, one can carry out nonlinear fracture mechanics assessments in structures that show instabilities due buckling (local/global), ovalisation and large rigid body motion. Many constraint-measuring parameters have been proposed, with the Q-parameter or the T-stress being the most popular ones. Solid finite element meshing for complex structures such as pipes containing semi-elliptical surface cracks in order to compute Q is at present not a feasible approach. However, shell structures are most conveniently meshed with shell finite elements, and the linespring finite element is a natural way of accounting for surface cracks. The T-stress is readily obtained from the linespring membrane force and bending moment along the surface crack. In this study we present a new approach to analyse cracked shell structures subjected to large geometric changes. By numerical examples it is shown how geometric instabilities and fracture compete as governing failure mode.     

On the dual iterative stochastic perturbation‐based finite element method in solid mechanics with Gaussian uncertainties

The main idea is a dual mathematical formulation and computational implementation of the iterative stochastic perturbation‐based finite element method for both linear and nonlinear problems in solid mechanics. A general‐order Taylor expansion with random coefficients serves here for the iterative determination of the basic probabilistic characteristics, where linearization procedure widely applicable in stochastic perturbation technique is replaced with the iterative one. The expected values and, in turn, the variances are derived first, and then, they are substituted into the equations for higher central probabilistic moments and additional probabilistic characteristics. The additional formulas for up to the fourth‐order probabilistic characteristics are derived thanks to the 10th‐order Taylor expansion. Computational implementation of this idea in the stochastic finite element method is provided by using the direct differentiation method and, independently, the response function method with polynomial basis. Numerical experiments include the simple tension of the elastic bar, nonlinear elasto‐plastic analysis of the aluminum 2D truss, and solution to the homogenization problem of periodic fiber‐reinforced composite with random elastic properties. The expected values, coefficients of variation, skewness, and kurtosis of the structural response determined via this iterative scheme are contrasted with these estimated by the Monte Carlo simulation as well as with the results of the semi‐analytical probabilistic technique following the response function method itself. Although the entire methodology is illustrated here by using the Gaussian variables where all odd‐order terms simply vanish, it can be extended towards non‐Gaussian processes as well and completed with all the perturbation orders. Copyright © 2015 John Wiley & Sons, Ltd.     

Probabilistic fracture mechanics by the boundary element method

In this paper, a new boundary element formulation is presented for probabilistic analysis of crack problems in context of linear elastic fracture mechanics. The method involves an implicit differentiation method for calculating fracture response derivatives with respect to random parameters and applies first-order reliability method to predicting the reliability of cracked structures. Direct differentiation is used to obtain the derivatives of fracture parameters with respect to crack size which is required for probabilistic analysis. Numerical examples are presented to illustrate the proposed method. The accuracy of the proposed method is demonstrated by comparing with either analytical solution or other established method.     

Element‐local level set method for three‐dimensional dynamic crack growth

An approximate level set method for three‐dimensional crack propagation is presented. In this method, the discontinuity surface in each cracked element is defined by element‐local level sets (ELLSs). The local level sets are generated by a fitting procedure that meets the fracture directionality and its continuity with the adjacent element crack surfaces in a least‐square sense. A simple iterative procedure is introduced to improve the consistency of the generated element crack surface with those of the adjacent cracked elements. The discrete discontinuity is treated by the phantom node method which is a simplified version of the extended finite element method (XFEM). The ELLS method and the phantom node technology are combined for the solution of dynamic fracture problems. Numerical examples for three‐dimensional dynamic crack propagation are provided to demonstrate the effectiveness and robustness of the proposed method. Copyright © 2009 John Wiley & Sons, Ltd.     

T-stress solutions for two-dimensional crack problems in anisotropic elasticity using the boundary element method

The importance of a two‐parameter approach in the fracture mechanics analysis of many cracked components is increasingly being recognized in engineering industry. In addition to the stress intensity factor, the T stress is the second parameter considered in fracture assessments. In this paper, the path‐independent mutual M‐integral method to evaluate the T stress is extended to treat plane, generally anisotropic cracked bodies. It is implemented into the boundary element method for two‐dimensional elasticity. Examples are presented to demonstrate the veracity of the formulations developed and its applicability. The numerical solutions obtained show that material anisotropy can have a significant effect on the T stress for a given cracked geometry.     

A dimensional decomposition method for stochastic fracture mechanics

This paper presents a new dimensional decomposition method for obtaining probabilistic characteristics of crack-driving forces and reliability analysis of general cracked structures subject to random loads, material properties, and crack geometry. The method involves a novel function decomposition permitting lower-variate approximations of a crack-driving force or a performance function, Lagrange interpolations for representing lower-variate component functions, and Monte Carlo simulation. The effort required by the proposed method can be viewed as performing deterministic fracture analyses at selected input defined by sample points. Compared with commonly-used first- and second-order reliability methods, no derivatives of fracture response are required by the new method developed. Results of three numerical examples involving both linear-elastic and nonlinear fracture mechanics of cracked structures indicate that the decomposition method provides accurate and computationally efficient estimates of probability density of the J-integral and probability of fracture initiation for various cases including material gradation characteristics and magnitudes of applied stresses and loads.     

工程弹塑性断裂随机分析的多维插值方法

为解决用确定性有限元或已有无量纲断裂参数表值时进行弹塑性断裂随机分析的问题，本文研究了利用数值插值方法进行随机参数计算的可行性。以单边裂纹板为例，根据确定性有限元计算的无量纲断裂参数表值，通过多维插值计算获得断裂参数及其对基本随机变量变化率，应用“弹塑性断裂随机分析的工程方法”进行弹塑性断裂随机分析。算例表明在敏感区用插值方法计算的结果小于用随机有限元方法计算的结果。     

Damage tolerance reliability analysis of automotive spot-welded joints

This paper develops a damage tolerance reliability analysis methodology for automotive spot-welded joints under multi-axial and variable amplitude loading history. The total fatigue life of a spot weld is divided into two parts, crack initiation and crack propagation. The multi-axial loading history is obtained from transient response finite element analysis of a vehicle model. A three-dimensional finite element model of a simplified joint with four spot welds is developed for static stress/strain analysis. A probabilistic Miner's rule is combined with a randomized strain-life curve family and the stress/strain analysis result to develop a strain-based probabilistic fatigue crack initiation life prediction for spot welds. Afterwards, the fatigue crack inside the base material sheet is modeled as a surface crack. Then a probabilistic crack growth model is combined with the stress analysis result to develop a probabilistic fatigue crack growth life prediction for spot welds. Both methods are implemented with MSC/NASTRAN and MSC/FATIGUE software, and are useful for reliability assessment of automotive spot-welded joints against fatigue and fracture.     

油管螺纹应力应变场的有限元分析与检测

利用ANSYS大型分析软件,建立了属于表面非线性和材料非线性相偶合问题的油管接头的弹塑性轴对称接触有限元模型,对油管接头在不同载荷工况和不同配合下的应力应变场进行了数值计算。又用电测法对油管接头的油管内壁进行了实际测量。对比两种不同方法得到的结果,验证了轴对称模型有限元法的适用范围。可靠地给出了工程中最危险截面应力应变场的各种变化规律,为研究油管的疲劳断裂问题提供了重要的理论分析依据和数据。     

On semi‐analytical probabilistic finite element method for homogenization of the periodic fiber‐reinforced composites

The main aim of this paper is a development of the semi‐analytical probabilistic version of the finite element method (FEM) related to the homogenization problem. This approach is based on the global version of the response function method and symbolic integral calculation of basic probabilistic moments of the homogenized tensor and is applied in conjunction with the effective modules method. It originates from the generalized stochastic perturbation‐based FEM, where Taylor expansion with random parameters is not necessary now and is simply replaced with the integration of the response functions. The hybrid computational implementation of the system MAPLE with homogenization‐oriented FEM code MCCEFF is invented to provide probabilistic analysis of the homogenized elasticity tensor for the periodic fiber‐reinforced composites. Although numerical illustration deals with a homogenization of a composite with material properties defined as Gaussian random variables, other composite parameters as well as other probabilistic distributions may be taken into account. The methodology is independent of the boundary value problem considered and may be useful for general numerical solutions using finite or boundary elements, finite differences or volumes as well as for meshless numerical strategies. Copyright © 2011 John Wiley & Sons, Ltd.     

Anomalous behavior of bilinear quadrilateral finite elements for modeling cohesive cracks with XFEM/GFEM

The performance of partition‐of‐unity based methods such as the generalized finite element method or the extended finite element method is studied for the simulation of cohesive cracking. The focus of investigation is on the performance of bilinear quadrilateral finite elements using these methods. In particular, the approximation of the displacement jump field, representing cohesive cracks, by extended finite element method/generalized finite element method and its effect on the overall behavior at element and structural level is investigated. A single element test is performed with two different integration schemes, namely the Newton‐Cotes/Lobatto and the Gauss integration schemes, for the cracked interface contribution. It was found that cohesive crack segments subjected to a nonuniform opening in unstructured meshes (or an inclined crack in a structured finite element mesh) result in an unrealistic crack opening. The reasons for such behavior and its effect on the response at element level are discussed. Furthermore, a mesh refinement study is performed to analyze the overall response of a cohesively cracked body in a finite element analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

Stochastic analysis on flexural behavior of reinforced concrete beams based on piecewise response surface scheme

To save the computational efforts of Monte Carlo simulation together with nonlinear finite element method, the analysis framework combining the response surface method and Monte Carlo simulation is usually adopted to investigate the stochastic nonlinear behavior of structures. It is found that the traditional response surface method could not describe stochastic behavior of cracked concrete beams. In order to overcome the discontinuity caused by cracking nonlinearity, a scheme by introducing a piecewise response surface is proposed in this paper. This scheme is evaluated by the stochastic analysis of several concrete beams. The comparison between the proposed method and some traditional methods, including Monte Carlo simulation and Monte Carlo simulation with traditional response surface, shows that the proposed method could well depict the stochastic behavior of concrete beams. Finally, the probability density evolution of concrete beams from short-term deflection to a long-term deflection is analyzed.     

A probabilistic fracture mechanics model including 3D ductile tearing of bi-axially loaded pipes with surface cracks

This paper presents a probabilistic fracture mechanics model established from three-dimensional FEM analyses of surface cracked pipes subjected to tension load in combination with internal pressure. The models are particularly interesting for offshore pipelines under operational conditions or during laying, where inelastic deformations may occur. In the numerical models, the plastic deformations, including ductile tearing effects, are accounted for by use of the Gurson-Tvergaard-Needleman model. This model is calibrated to represent a typical X65 pipeline steel behaviour under ductile crack growth and collapse. Several parameters are taken into account, such as crack depth, crack length and material hardening. Another important topic is the examination of the influence of bi-axial loading due to internal pressure on capacity. From the results of the deterministic analyses a probabilistic fracture mechanics model is established using the response surface methodology. Two failure criteria are examined to represent the structural capacity. Based on the established model, we illustrate the methodology by examples employing the two different failure criteria solved with first and second order reliability methods.