The effect of material combinations and relative crack size to the stress intensity factors at the crack tip of a bi-material bonded strip

Several types of singular stress fields may appear at the corner where an interface between two bonded materials intersects a traction-free edge depending on the material combinations. Since the failure of the multi-layer systems often originates at the free-edge corner, the analysis of the edge interface crack is the most fundamental to simulate crack extension. In this study, the stress intensity factors for an edge interfacial crack in a bi-material bonded strip subjected to longitudinal tensile stress are evaluated for various combinations of materials using the finite element method. Then, the stress intensity factors are calculated systematically with varying the relative crack sizes from shallow to very deep cracks. Finally, the variations of stress intensity factors of a bi-material bonded strip are discussed with varying the relative crack size and material combinations. This investigation may contribute to a better understanding of the initiation and propagation of the interfacial cracks.     

Computation of the stress intensity factors for repaired cracks with bonded composite patch in mode I and mixed mode

In this study, the finite element method is used to analyse the behaviour of repaired cracks with bonded composite patches in mode I and mixed mode by computing the stress intensity factors at the crack tip. The effects of the patch size and the adhesive properties on the stress intensity factors variation were highlighted. The plot of the stress intensity factors according to the crack length in mode I, shows that the stress intensity factor exhibits an asymptotic behaviour as the crack length increases. In mixed mode, the obtained results show that the Mode I stress intensity factor is more affected by the presence of the patch than that of mode II.     

Static finite element stress intensity factors for annular cracks

Results of finite element static stress intensity factor calculations for an annular crack around a spherical inclusion (void) are presented and compared with those from approximate analytical methods.     

Stress intensity factors and crack opening displacements for slanted axial through-wall cracks in pressurized pipes

This paper proposes elastic stress intensity factors and crack opening displacements (CODs) for a slanted axial through-wall cracked cylinder under an internal pressure based on detailed three-dimensional (3D) elastic finite element (FE) analyses. The FE model and analysis procedure were validated against existing solutions for both elastic stress intensity factor and COD of an idealized axial through-wall cracked cylinder. To cover a practical range, four different values of the ratio of the mean radius of cylinder to the thickness ( R _m/ t ) were selected. Furthermore, four different values of the normalized crack length and five different values of the ratio of the crack length at the inner surface to the crack length at the outer surface representing the slant angle were selected. Based on the elastic FE results, the stress intensity factors along the crack front and CODs through the thickness at the centre of the crack were provided. These values were also tabulated for three selected points, that is, the inner and outer surfaces and at the mid-thickness. The present results can be used to evaluate the crack growth rate and leak rate of a slanted axial through-wall crack due to stress corrosion cracking and fatigue. Moreover, the present results can be used to perform a detailed leak-before-break analysis considering more realistic crack shape development.     

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

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

Stress intensity factor analysis for an interface crack between dissimilar isotropic materials under thermal stress

Thermal stresses, one of the main causes of interfacial failure between dissimilar materials, arise from different coefficients of linear thermal expansion. Two efficient numerical procedures in conjunction with the finite element method (FEM) for the stress intensity factor (SIF) analysis of interface cracks under thermal stresses are presented. The virtual crack extension method and the crack closure integral method are modified using the superposition method. The SIF analyses of some interface crack problems under mechanical and thermal loads are demonstrated. Very accurate mode separated SIFs are obtained using these methods.     

Multiple penny-shaped cracks interaction in a finite body and their effect on stress intensity factor

In this paper we investigate the stress intensity factors (SIFs) of multiple penny-shaped cracks in an elastic solid cylinder under mode I (axial tension) loading. The cracks are located symmetrically and in parallel to one another in the isotropic cylinder. The fractal-like finite element method (FFEM) is employed to study the interaction of multiple cracks and to demonstrate the efficiency of the FFEM for multiple crack problems. The results show that the SIF values of the inner cracks, which are denoted as crack number 1,2,3,…,(n+1)/2 of a stack of n parallel cracks, are lower than the SIF values of a single crack by between 16% and 48%. Also, the outermost crack, that is the crack closest to the boundaries of a multiple cracked body, has the highest SIF values and is, therefore, likely to fail first.     

Stress intensity factors of an inclined elliptical crack near a bimaterial interface

In this paper the stress intensity factors are discussed for an inclined elliptical crack near a bimaterial interface. The solution utilizes the body force method and requires Green’s functions for perfectly bonded semi-infinite bodies. The formulation leads to a system of hypersingular integral equation whose unknowns are three modes of crack opening displacements. In the numerical calculation, unknown body force densities are approximated by using fundamental density functions and polynomials. The results show that the present method yields smooth variations of stress intensity factors along the crack front accurately. Distributions of stress intensity factors are presented in tables and figures with varying the shape of crack, distance from the interface, and elastic modulus ratio. It is found that the inclined crack can be evaluated by the models of vertical and parallel cracks within the error of 24% even for the cracks very close to the interface.     

Stress intensity factors for through-thickness cracks in a wide plate: Derivation and application to arbitrary weld residual stress fields

Exact closed-form stress intensity factor (SIF) solutions have been developed for mode-I, II and III through-thickness cracks in an infinite plate. Centre-crack problems have been analysed comprehensively in the literature, but the focus has been on the effect of simple loading about the crack centre. In the current work, the formula of Sih-Paris-Erdogan was extended to consider the difference in SIF on the left and right crack tips under an asymmetric stress field. Mathematical manipulations were performed to derive exact stress magnification factors for SIF computations and simultaneously circumvent the problem of crack-tip stress singularity. The solutions so obtained are applied to derive the residual SIFs that would act on a crack growing under the influence of the residual stress fields associated with VPPA (variable polarity plasma arc) and friction stir welds, using measured residual stress profiles.     

Stress intensity factor analysis of a three-dimensional interface crack between dissimilar anisotropic materials

A new numerical method to calculate the stress intensity factors (SIFs) of a three-dimensional interface crack between dissimilar anisotropic materials was developed. In this study, the M-integral method was employed for mode separation of the SIFs. The moving least-square method was utilized to calculate the M-integral. Using the M-integral with the moving least-square method, SIFs can be automatically calculated with only the nodal displacements from the finite element method (FEM). Here, SIFs analyses of some typical three-dimensional problems are demonstrated. Excellent agreement was achieved between the numerical results obtained by the present method and the corresponding results proposed by other researchers. In addition, the SIFs of a single-edge crack, a through crack, and a semi-circular crack between two anisotropic solids in three-dimensional structures were analyzed.     

Stress intensity factor analyses of interface cracks between dissimilar anisotropic materials using the finite element method

Delamination along an interface between dissimilar materials is the primary cause of failure in microstructures like electronic packages, micro-electro-mechanical systems (MEMS), and so on. Fracture mechanics is a powerful tool for the evaluation of delamination. However, many materials used in microstructures such as composite materials and single crystals are anisotropic materials. Stress intensity factors of an interface crack between dissimilar anisotropic materials, which were proposed by Hwu, are useful for evaluating the reliability of microstructures. However, numerical methods that can analyze the stress intensity factors of an interface crack between anisotropic materials have not been developed. We propose herein a new numerical method for the analysis of an interface crack between dissimilar anisotropic materials. The stress intensity factors of an interface crack are based on the generalized plane strain condition. The energy release rate is obtained by the virtual crack extension method in conjunction with the finite element method for the generalized plane strain condition. The energy release rate is separated into individual modes of the stress intensity factors K_I, K_II, and K_III, using the principal of superposition. The target problem to be solved is superposed on the asymptotic solution of displacement in the vicinity of an interface crack tip, which is described using the Stroh formalism. Analyses of the stress intensity factors of center interface cracks between semi-infinite dissimilar anisotropic media subjected to concentrated self-balanced loads on the center of crack surfaces and to uniform loads are demonstrated. The present method accurately provides mode-separated stress intensity factors using relatively coarse meshes for the finite element method.     

Stress intensity factors for slanted through-wall cracks based on elastic finite element analyses

Based on detailed 3‐dimensional (3‐D) elastic finite element (FE) analyses, the present paper provides stress intensity factors (SIFs) for plates with slanted through‐wall crack (TWC) and cylinders with slanted circumferential TWC. Regarding loading conditions, axial tension was considered for the plates, whereas axial tension, global bending and internal pressure were considered for the cylinders. To cover a practical range, the geometric variables affecting the SIF were systematically varied. Based on FE results, SIFs along the crack front, including the inner and outer surface points, were provided. The present results can be used to evaluate the fatigue crack growth or stress corrosion cracking behaviour of a slanted TWC and furthermore to perform detailed Leak‐Before‐Break analysis considering a more realistic crack shape.     

Stress intensity factors for internal circular cracks in fibers under tensile loading

Stress intensity factors for inner circular cracks placed eccentrically in a fiber with round cross section were computed and are presented in this paper in both analytical and graphical form. The crack plane was perpendicular to the fiber axis and remote tensile loading was assumed. The stress intensity factors were numerically computed using the finite element method. Mesh objectivity and some other aspects of computational precision are considered. The asymptotic behaviour when the crack size and the ligament depth vanish were considered in order to formulate accurate interpolation expressions.     

Stress intensity factor for a crack normal to an interface between two orthotropic materials

In this paper, the problem of a crack normal to an interface in two joined orthotropic plates is studied as a plane problem. Body force method is used to investigate dependence of the stress intensity factor on the elastic constants: E _x1, E _y1, G _xy1, V _xy1 for material 1 and E _x2, E _y2, G _xy2, V _xy2 for material 2. A particular attention is paid to simplifying kernel functions, which is used in the body force method, so that all the elastic constants involved can be represented by three new parameters: H ₁, H _2I, H ₃ for the mode I deformation and H ₁, H _2II, H ₃ for the mode II deformation. From the kernel function so obtained it is found that the effects of the eight elastic constants on the stress intensity factors can be expressed by the three material parameters, H ₁, H _2I, H ₃ and H ₁, H _2II, H ₃, respectively for K _I and K _II. Furthermore, it is also found that the dependence of K _I on H ₁, H _2I, H ₃ is exactly the same as the dependence of K _II on H ₁, H _2II, H ₃. This revised version was published online in August 2006 with corrections to the Cover Date.     

Numerical analysis of the stress intensity factors for repaired cracks from a notch with bonded composite semicircular patch

In this paper, we investigated the crack growth behaviour of cracked thin aluminium plate repaired with bonded composite patch. The finite element method is used to study the performance of the bonded composite reinforcement or repair for reducing the stress concentration at a semicircular lateral notch and repairing cracks emanating from this kind of notch. The effects of the adhesive properties and the patch size on the stress intensity factor variation at the crack tip in mode I were highlighted. The obtained results show that the stress concentration factor at the semicircular notch root and the stress intensity factor of a crack emanating from notch are reduced with the increase of the diameter and the number of the semicircular patch. The maximal reduction of stress intensity factor is about 42% and 54%, respectively, for single and double patch. However, the gain in the patch thickness increases with the increase of the crack length and it decreases when the patch thickness increases. The adhesive properties must be optimised in order to increase the performance of the patch repair or reinforcement.     

The effects of disbonds on the pure mode II stress intensity factor of aluminium plate reinforced with bonded composite materials

Adhesively bonded composite patch repair has been widely used to restore or extend the service life of cracked structural components due to its effectiveness to mechanical repair technique. In this work, the finite element method is applied to analyse the performance of the bonded composite patch for repairing cracks emanating from semicircular notch root in pure mode II. The stress intensity factor was computed at the crack tip repaired using a boron/epoxy patch for different orientation of fibers, taking into account the disbond. In this case, the increase of a patch thickness reduces the negative effects of disbond. When this effect is significant between the patch and the plate, it reduces the repair effectiveness. The maximum reduction obtained by using a boron/epoxy of fibers in the x-direction is of the order of 20% more important compared to a patch having its fibers in the y-direction. The stress intensity factor exhibits an asymptotic behaviour as the disbond size increases.     

Stress intensity factors of semi-elliptical surface cracks in pressure vessels by global-local finite element methodology

In the present study the problem of calculation of the stress intensity factors (SIF) of semi-elliptical cracks located in the stress concentration areas of a pressure vessel is numerically solved by advanced global-local finite element (FE) analysis. The common characteristic of the cases solved is that the stress field at the crack area varies along the axial, the circumferential, as well as, the through-the-thickness directions. SIF solutions for such problems are not available, neither analytically, nor numerically, as the currently existing solutions in the literature (numerical results, Newman-Raju empirical equations, weight function solutions, etc.) are only valid for uniform stress distribution along the axial and circumferential directions of the pressure vessel and allow variation only through-the-thickness. The crack locations considered are the intersection of the cylinder to a nozzle and the connection of the cylinder with its hemi-spherical ends. The stress intensity factors are presented in a suitable table format for various geometrical configurations of both the pressure vessel and the semi-elliptical crack, thus providing a useful tool for the fracture mechanics design of cracked pressure vessels. The modeling details of the sub-structuring methodology, employed in the analysis, are extensively discussed and the numerical approach is proven to be very efficient for the SIF calculation of pressure vessel semi-elliptical cracks.     

Shear loaded interface crack under the influence of friction: a finite difference solution

Formulation of the elastic two‐dimensional problem of contact with friction is presented. Two‐dimensional equilibrium equations and boundary conditions in an orthogonal curvilinear co‐ordinate system are written explicitly. The above formulation is solved with the aid of the finite difference technique. An iterative algorithm which does not require load increments is employed for solving interface fracture problems with contact and friction subjected to a monotonically increasing load. The J‐integral is extended for problems in which there is friction along the crack faces. Stress intensity factors are calculated by means of the J‐integral, as well as an asymptotic expansion of the tangential shift. Two problems are analysed: (1) a crack in homogeneous material in the presence of friction involving stationary contact; and (2) an interface crack in the presence of friction involving receding contact. Results are compared to those found by analytical and semi‐analytical methods which are presented in the literature, as well as to those obtained by means of the finite element method. The accuracy of the results establishes the reliability of the finite difference analysis, as well as the post‐processors. In addition, a problem involving stick conditions is considered. It is observed that with increasing friction, the normal gaps and tangential shifts decrease. The size of the contact zone increases and values of the stress intensity factor decrease. Copyright © 2004 John Wiley & Sons, Ltd.     

Stress intensity factors in spot welds

This paper looks at stress intensity factors of cracks in resistance spot welded joints. Stress intensity factors have been used in the past to predict fatigue crack propagation life of resistance spot welds. However, the stress intensity factors from all previous work was based on assumed initial notch cracks at the nugget, parallel to the sheets. Physical evidence shows, however, that fatigue cracks in spot welds propagate through the thickness of the sheets rather than through the nugget. In this work, stress intensity factors of assumed notch cracks and through thickness cracks in tensile shear (TS) and modified coach peel (MCP) specimens were determined by the finite element method. The finite element results from the assumed notch cracks were compared with the results in the literature and were found to be in agreement with the results from Zhang’s equations [Int. J. Fract. 88 (1997) 167]. The stress intensity factors of assumed notch cracks were found to be different from those of through thickness cracks. To date, no analytic equations for stress intensity factors of through thickness cracks in spot welds have been published. In the current work, simple equations are proposed to estimate the K_I and K_II values of through thickness cracks in TS and MCP specimens.     

Stress intensity factors for penny-shaped cracks perpendicular to graded interfacial zone of bonded bi-materials

This paper examines the stress intensity factors that are associated with a penny-shaped crack perpendicular to the interface of a bi-material bonded with a graded interfacial zone. Elastic modulus of the graded interfacial zone is assumed to be an exponential function of the depth. The stress intensity factors are calculated numerically using a so-called generalized Kelvin solution based boundary element method. Three cases of normal or shear tractions acting on the crack surfaces are examined. Values of the stress intensity factors are examined by taking into account the effects of the following four parameters: (a) the crack front position; (b) the non-homogeneity parameter of the graded interfacial zone; (c) the crack distance to the graded interfacial zone; and (d) the graded interfacial zone thickness. The numerical results are compared well with existing solutions under some degenerated conditions. These results are useful to furthering our knowledge on fracture behavior of bi-material systems with or without a graded interfacial zone.