Fracture analysis of interfacial crack by global-local finite element

An eigenfunction expansion is used to formulate the global element on the crack tip. The global-local finite element method employs both conventional finite element and classical Rayleigh-Ritz kinematic approach. The hybrid Ritz method not only preserves the finite element modelling capability but adds the advantage of using prior information regarding the anticipate behaviour of the particular problem. Thus, it is able to achieve better accuracy with fewer elements in comparison with conventional finite element. Several examples relative to crack problems are presented to demonstrate the global-local finite element method.     

A finite element analysis of cracked square plates and bars under antiplane loading

ABSTRACT Finite element analyses were carried out on cracked 20 mm square plates and bars ranging in thickness from 2.5 mm to a length of 60 mm. The crack extended from the middle of one side of the square to its centre, and was modelled as a narrow, parallel‐sided notch with a semicircular tip. An antiplane loading was applied to the side containing the crack. An infinitely long bar under the antiplane loading used is in pure Mode III. It was found that the central portions of 40, 56 mm and 60 mm long bars were in pure Mode III, and also that K_III was approximately constant. These central portions were therefore representative of an infinitely long bar. Towards the ends of a bar K_III decreased. At the ends of a bar corner point effects meant that Mode II stress intensity factors and displacements were induced in the corner region. The size of the corner region was independent of bar length. In the 2.5, 5 and 10 mm thick plates out of plane bending means that the antiplane loading became a mixed Mode II and Mode III loading. At a centre line K_II is zero by symmetry. Behaviour in the corner region was a function of plate thickness. For both plates and bars, as has been predicted theoretically, the ratio K_II/K_III tends to a constant value as a surface is approached. For a thickness of 20 mm, that is a 20‐mm cube, behaviour represents a transition between plate and bar behaviour.     

Fractal finite element method based shape sensitivity analysis of multiple crack system

This paper presents fractal finite element based continuum shape sensitivity analysis for a multiple crack system in a homogeneous, isotropic, and two dimensional linear-elastic body subjected to mixed-mode (modes I and II) loading conditions. The salient feature of this method is that the stress intensity factors and their derivatives for the multiple crack system can be obtained efficiently since it only requires an evaluation of the same set of fractal finite element matrix equations with a different fictitious load. Three numerical examples are presented to calculate the first-order derivative of the stress intensity factors or energy release rates.     

Fully automatic modelling of mixed-mode crack propagation using scaled boundary finite element method

The newly-developed scaled boundary finite element method (SBFEM) is able to calculate stress intensity factors directly because the singularity in stress solutions at crack tips is analytically represented. By taking this advantage, a mixed-mode crack propagation model based on linear elastic fracture mechanics (LEFM) was developed in this study. A domain is first divided into a few subdomains. Because the dimensions and shapes of subdomains can be flexibly varied and only the domain boundaries or common edges between subdomains are discretised in the SBFEM, a remeshing procedure as simple as in boundary element methods was developed with minimum mesh changes whereas the generality and flexibility of the FEM is well maintained. Fully-automatic modelling of mixed-mode crack propagation is then achieved by combining the remeshing procedure with a propagation criterion. Three mixed-mode examples were modelled. Comparisons of the numerical results with those from available publications show that the developed model is capable of predicting crack trajectories and load-displacement relations accurately and efficiently.     

非均匀材料界面裂纹的Cell-Based光滑有限元法

为提高非均匀材料界面裂纹尖端断裂参数的求解精度,基于非均匀材料界面断裂力学、Cell-Based光滑有限元(Cell-SFEM)和非均匀材料的互交作用积分法,提出了求解非均匀材料界面裂纹尖端断裂参数的CellBased光滑有限元法,推导了基于Cell-Based光滑有限元法的非均匀材料的互交作用积分法,对非均匀材料间的界面裂纹尖端处正则应力强度因子进行了求解,并与参考解进行了比较,讨论了互交积分区域大小和光滑子元个数与正则应力强度因子的关系。数值算例结果表明:本方法具有很高的计算精度,对积分区域大小不敏感,可为设计、制造抗破坏非均匀材料提供依据。     

Fracture mechanics analysis using the wavelet Galerkin method and extended finite element method

This paper presents fracture mechanics analysis using the wavelet Galerkin method and extended finite element method. The wavelet Galerkin method is a new methodology to solve partial differential equations where scaling/wavelet functions are used as basis functions. In solid/structural analyses, the analysis domain is divided into equally spaced structured cells and scaling functions are periodically placed throughout the domain. To improve accuracy, wavelet functions are superposed on the scaling functions within a region having a high stress concentration, such as near a hole or notch. Thus, the method can be considered a refinement technique in fixed‐grid approaches. However, because the basis functions are assumed to be continuous in applications of the wavelet Galerkin method, there are difficulties in treating displacement discontinuities across the crack surface. In the present research, we introduce enrichment functions in the wavelet Galerkin formulation to take into account the discontinuous displacements and high stress concentration around the crack tip by applying the concept of the extended finite element method. This paper presents the mathematical formulation and numerical implementation of the proposed technique. As numerical examples, stress intensity factor evaluations and crack propagation analyses for two‐dimensional cracks are presented. Copyright © 2012 John Wiley & Sons, Ltd.     

An over‐deterministic method for calculation of coefficients of crack tip asymptotic field from finite element analysis

An over‐deterministic method has been employed for calculating the stress intensity factors (SIFs) as well as the coefficients of the higher‐order terms in the Williams series expansions in cracked bodies, using the conventional finite element analysis. For a large number of nodes around the crack tip, an over‐determined set of simultaneous linear equations is obtained, and using the fundamental concepts of the least‐squares method, the coefficients of the Williams expansion can be calculated for pure mode I, pure mode II and mixed mode I/II conditions. A convergence study has been conducted to examine the effects of the number of nodes used, the number of terms in Williams expansion and the distance of the selected nodes from the crack tip, on the accuracy of the results. It is shown that the simple method presented in this paper, yields accurate results even for coarse finite element meshes or in the absence of singular elements. The accuracy of SIFs and the coefficients of higher‐order terms are validated by using the available results in the literature.     

An embedded cohesive crack model for finite element analysis of brickwork masonry fracture

This paper presents a numerical procedure for fracture of brickwork masonry based on the strong discontinuity approach. The model is an extension of the cohesive model prepared by the authors for concrete, and takes into account the anisotropy of the material. A simple central-force model is used for the stress versus crack opening curve. The additional degrees of freedom defining the crack opening are determined at the crack level, thus avoiding the need of performing a static condensation at the element level. The need for a tracking algorithm is avoided by using a consistent procedure for the selection of the separated nodes. Such a model is then implemented into a commercial code by means of a user subroutine, consequently being contrasted with experimental results. Fracture properties of masonry are independently measured for two directions on the composed masonry, and then input in the numerical model. This numerical procedure accurately predicts the experimental mixed-mode fracture records for different orientations of the brick layers on masonry panels.     

Stress intensity factors of a central slant crack with frictional surfaces in plates with biaxial loading

The stress intensity factor is a traditional topic in mechanics and there have been many solutions for many different cases. The closed frictional crack problem has been modeled in the rock mechanics field where fractures are mostly under compression. Further, the effect of finite plate dimensions under biaxial loading has not been considered in the literature. The key contribution of the present paper is to evaluate the effect of the crack length to plate width ratio on the mode I and II stress intensity factors (SIF) of a central slant crack with frictional surfaces in plates with biaxial loading of different patterns, i.e. tension-tension, tension-compression, compression-tension or compression-compression. A plane strain elastic two-dimensional finite element analysis was adopted. Crack length to plate width ratios equal to 0.1, 0.3 and 0.5 with biaxial ratios from –1 to 1, crack angles from 0° to 90° and friction coefficients from 0 to 1 were considered. Contact regimes and the effect of the crack length to plate width ratio were found dependent on biaxial ratio and pattern, friction coefficient and crack angle.     

Fatigue crack closure determination by means of finite element analysis

Plasticity-induced crack closure is an observed phenomenon during fatigue crack growth. However, accurate determination of fatigue crack closure has been a complex task for years. It has been approached by means of experimental and numerical methods. The finite element method (FEM) has been the principal numerical tool employed. In this paper the results of a broad study of fatigue crack closure in plane stress and plane strain by means of FEM are presented. The effect of three principal factors has been analysed in depth, the maximum load, the crack length and the stress ratio. It has been found that the results are independent of maximum load and the crack length, and there exists a direct influence of the stress ratio. This relation has been numerically correlated and compared with experimental results. Differences have also been established between opening and closure points and between the different criteria employed to compute crack closure.     

An embedded cohesive crack model for finite element analysis of mixed mode fracture of concrete*

An embedded cohesive crack model is proposed for the analysis of the mixed mode fracture of concrete in the framework of the Finite Element Method. Different models, based on the strong discontinuity approach, have been proposed in the last decade to simulate the fracture of concrete and other quasi‐brittle materials. This paper presents a simple embedded crack model based on the cohesive crack approach. The predominant local mode I crack growth of the cohesive materials is utilized and the cohesive softening curve (stress vs. crack opening) is implemented by means of a central force traction vector. The model only requires the elastic constants and the mode I softening curve. The need for a tracking algorithm is avoided using a consistent procedure for the selection of the separated nodes. Numerical simulations of well‐known experiments are presented to show the ability of the proposed model to simulate the mixed mode fracture of concrete.     

Elastodynamic analysis of the Finite punch and finite crack problems in orthotropic materials

The transient elastodynamic response of the finite punch and finite crack problems in orthotropic materials is examined. Solution for the stress intensity factor history around the punch corner and crack tip is found. Laplace and Fourier transforms together with the Wiener–Hopf technique are employed to solve the equations of motion in terms of displacements. A detailed analysis is made in the simplified case when a flat rigid punch indents an elastic orthotropic half-plane, the punch approaches with a constant velocity normally to the boundary of the half-plane. An asymptotic expression for the singular stress near the punch corner is analyzed leading to an explicit expression for the dynamic stress intensity factor which is valid for the time the dilatational wave takes to travel twice the punch width. In the crack problem, a finite crack is considered in an infinite orthotropic plane. The crack faces are loaded by impact uniform pressure in mode I. An expression for the dynamic stress intensity factor is found which is valid while the dilatational wave travels the crack length twice. Results for orthotropic materials are shown to converge to known solutions for isotropic materials derived independently.     

Stress intensity factors and weight functions for deep semi-elliptical surface cracks in finite-thickness plates

ABSTRACT Three-dimensional finite element analyses have been conducted to calculate the stress intensity factors for deep semi-elliptical cracks in flat plates. The stress intensity factors are presented for the deepest and surface points on semi-elliptic cracks with a/t -values of 0.9 and 0.95 and aspect ratios ( a/c ) from 0.05 to 2. Uniform, linear, parabolic or cubic stress distributions were applied to the crack face. The results for uniform and linear stress distributions were combined with corresponding results for surface cracks with a/t = 0.6 and 0.8 to derive weight functions over the range 0.05 ≤  a/c  ≤ 2.0 and 0.6 ≤  a/t  ≤ 0.95. The weight functions were then verified against finite element data for parabolic or cubic stress distributions. Excellent agreements are achieved for both the deepest and surface points. The present results complement stress intensity factors and weight functions for surface cracks in finite thickness plate developed previously.     

Analytical solutions for two penny-shaped crack problems in thermo-piezoelectric materials and their finite element comparisons

Thermally loaded penny-shaped cracks in thermopiezoelectric materials are investigated in this paper. The analytical solutions for the penny-shaped cracks subjected to uniform temperature and steady heat flow are discussed. Comparisons are made between the stress-intensity factors derived by the analytical solutions and the numerical results using different finite element techniques.     

Determination of optimum effective parameters on thermal buckling of hybrid composite plates with quasi-square cut-out using a genetic algorithm

ABSTRACT

The thermal buckling load on perforated composite plates is affected by several parameters, including design variables such as cut-out orientation, fibre angle, bluntness of cut-out corners, cut-out size to plate size ratio and stacking sequence. This study investigates the effect of these parameters on the thermal buckling load of a composite plate with a quasi-square cut-out. Optimal values of the parameters are determined using a genetic algorithm to achieve the maximum buckling load. The composite used herein is a four-layer laminated composite plate. The stacking sequences of the plate are also studied. Stability equations are obtained using first order shear deformation theory. The results showed that a plate with a quasi-square cut-out is more resistant to thermal buckling than one with a circular cut-out; thermal buckling of a composite plate is dependent on various parameters, and the maximum thermal buckling load can be achieved by appropriate selection of these parameters.     

Stress intensity factors for surface fatigue crack in a round bar under cyclic axial loading

In this paper, the surface fatigue crack growth shape for an initial straight-fronted edge crack in an elastic bar of circular cross-section is determined through experiments under pure fatigue axial loading. Three different initial notch depths are discussed. The relations of the aspect ratio (b/c) and relative crack depth (b/D) are obtained, and it is shown that there is a great difference in the growth of cracks with different initial front shapes and crack depths. Further, using the three-dimensional finite element method, the stress intensity factors (SIFs) are determined under remote uniform tension loading. Since the relationship of b/c and b/D changes during the fatigue crack growth, the SIFs are determined for different surface crack configurations.     

Influence of minimum element size to determine crack closure stress by the finite element method

Measuring opening or closure stress is a complex process that influences the low accuracy of obtained data. Finite element models have been one of the available ways to deal with this problem. The difficulty of modelling the whole process of crack growth (due to the great number of cycles implied) as the great complexity of the phenomenon itself (with a high plastic strain concentrated in a small area, with elevated stress gradients) has made the results to be quite varied, being influenced by a great number of modelling parameters. Of those parameters, the minimum size of the element used to mesh the area around the crack tip vicinity presents a great influence on the results.In this work, a detailed analysis of the influence of this parameter in the results in terms of closure or opening stress is presented. The effect that different meshing criteria can have on the result is complex and it has been necessary to reduce the element size around the crack tip to a size that had not been reached before. Procedures and modelling criteria stricter than the ones shown in the current bibliography are proposed. A methodology for the correct interpretation of the results is also established.     

Dynamic analysis of fixed cracks in composites by the extended finite element method

This paper is dedicated to simulation of dynamic analysis of fixed cracks in orthotropic media using an extended finite element method. This work is in fact an extension to dynamic problems of the recently developed orthotropic extended finite element method for fracture analysis of composites. In this method, the Heaviside and near-tip enrichment functions are used in the framework of the partition of unity for modeling crack discontinuity and crack-tip singularities within the classical finite element method. In this procedure, elements that include a crack are not required to conform to crack edges. Therefore, mesh generation can be performed without any need to comply to crack edges and the method is capable of modeling the crack propagation without any remeshing. To determine the fracture properties, mixed-mode dynamic stress intensity factors (DSIFs) are evaluated by means of domain separation integral (J^′-integral) method. Results of the proposed method are compared with other available analytical and computational results.     

Multiple cohesive crack growth in brittle materials by the extended Voronoi cell finite element model

This paper is aimed at modeling the propagation of multiple cohesive cracks by the extended Voronoi cell finite element model or X-VCFEM. In addition to polynomial terms, the stress functions in X-VCFEM include branch functions in conjunction with level set methods and multi-resolution wavelet functions in the vicinity of crack tips. The wavelet basis functions are adaptively enriched to accurately capture crack-tip stress concentrations. Cracks are modeled by an extrinsic cohesive zone model in this paper. The incremental crack propagation direction and length are adaptively determined by a cohesive fracture energy based criterion. Numerical examples are solved and compared with existing solutions in the literature to validate the effectiveness of X-VCFEM. The effect of cohesive zone parameters on crack propagation is studied. Additionally, the effects of morphological distributions such as length, orientation and dispersion on crack propagation are studied.