BEM for crack-hole problems in thermopiezoelectric materials

The boundary element formulation for analysing interaction between a hole and multiple cracks in piezoelectric materials is presented. Using Green's function for hole problems and variational principle, a boundary element model (BEM) for a 2-D thermopiezoelectric solid with cracks and holes has been developed and used to calculate stress intensity factors of the crack-hole problem. In BEM, the boundary condition on the hole circumference is satisfied a priori by Green's function, and is not involved in the boundary element equations. The method is applicable to multiple-crack problems in both finite and infinite solids. Numerical results for stress and electric displacement intensity factors at a particular crack tip in a crack-hole system of piezoelectric materials are presented to illustrate the application of the proposed formulation.     

Cracks emanating from circular hole or square hole in rectangular plate in tension

A numerical analysis of cracks emanating from a circular hole (Fig. 1) or a square hole (Fig. 2) in rectangular plate in tension is performed by means of the displacement discontinuity method with crack-tip elements (a boundary element method) presented recently by the author. Detail solutions of the stress intensity factors (SIFs) of the two plane elastic crack problems are given, which can reveal the effect of geometric parameters of the cracked bodies on the SIFs. By comparing the SIFs of the two crack problems with those of the center crack in rectangular plate in tension (Fig. 3), in addition, an effect of the circular hole or the square hole on the SIFs of the center crack is discussed in detail. The numerical results reported here also illustrate that the boundary element method is simple, yet accurate for calculating the SIFs of complex crack problems in finite plate.     

BOUNDARY ELEMENT ANALYSIS OF CRACK INTERSECTION, INITIATION AND GROWTH AT PIN LOADED HOLES

Abstract The dual boundary element method for the analysis of cracks in linear elastic materials has been previously generalised by the authors to allow for automatic remeshing when crack tips intersect other cracks or boundaries, and initiation and growth of small cracks at positions of high stress concentration. The new cracks are assumed to result from sudden events such as an overload or the subsequent stress redistribution when cracks intersect other cracks or holes. In this paper a crack at the edge of one hole in a row of pin-loaded holes is investigated; various values are considered for the stress at which new cracks may initiate. Two rows of aligned or staggered holes are examined also. The spacing between the holes was typical for lines of holes in overlap joints in plates. For the same load transfer between the plates, new cracks are initiated less readily and grow more slowly for a double row of aligned holes than for a single row or for staggered rows.     

Cracks emanating from a square hole in rectangular plate in tension

A numerical analysis of cracks emanating from a square hole in a rectangular plate in tension is performed using a hybrid displacement discontinuity method (a boundary element method). Detailed solutions of the stress intensity factors (SIFs) of the plane elastic crack problem are given, which can reveal the effect of geometric parameters of the cracked body on the SIFs. By comparing the calculated SIFs of the plane elastic crack problem with those of the centre crack in a rectangular plate in tension, in addition, an amplifying effect of the square hole on the SIFs is found. The numerical results reported here also prove that the boundary element method is simple, yet accurate, for calculating the SIFs of complex crack problems in finite plate.     

Stress intensity factors and crack extension in a cracked pressurised cylinder

This paper considers the plane elastic problem corresponding to single or multiple radial cracks emanating from the internal boundary of a circular ring, under uniform external tension and internal pressure. The stress intensity factors are calculated by using the dual boundary element method with the J-integral technique. Accurate data are found for varying crack depths over a representative range of wall ratios for fracture mechanics applications to pressurised circular cylinders. The interaction of multiple cracks and crack extension are investigated in the case of an internal pressure loading condition. The analysis shows that, for a multi-cracked pressurised cylinder, it is sufficient to calculate the stress intensity of the main crack in isolation for the purposes of safety assessment.     

The problem of internal cracks in an infinite strip having a circular hole

Summary. The elastostatic plane problem of an infinite strip having a circular hole and containing two symmetrically located internal cracks perpendicular to the boundary is formulated in terms of triply coupled integral equations. The solution of the problem is obtained for various crack geometries and for uniaxial tension applied to the strip away from the crack region. Quantities of physical interest are displayed in graphical forms.     

T-stress Solutions for Cracks Emanating from a Circular Hole in a Finite Plate

The boundary element method is employed to obtain T-stress solutions for cracks emanating from a circular hole in finite rectangular plates. Numerical values of the T-stress are obtained using the M-contour integral approach. A range of crack lengths are analyzed for two hole sizes, and the cases of a single crack and double-cracks emanating from the hole in the plate under both uniform remote tension and simple bending are considered. For completeness, stress intensity factor solutions are also presented. These results will be useful for failure assessments using two-parameter linear elastic fracture mechanics.     

Spannungsfaktoren eines Risses in der Umgebung eines Kreisloches und ihr Einfluß auf das Bruchverhalten

Stress Intensity Factors in the Neighbourhood of a Circular Hole and their Influence on Crack Behavior . A photoelastic method was developed to determine the stress intensity factors K_I and K_II for cracks subjected to mixed-mode loading. The constants in the near field expansion about the crack tip were computed using a non-linear optimization program to give a best fit to the observed isochromatics. Copying the latter onto an equal density film increased their sharpness and, thus, the accuracy of the determination. The method was applied to cracks lying perpendicular to the external stress near a circular hole in a plate under uniaxial tension and the results used to describe the paths of cracks in the neighbourhood of a hole.     

Analysis of crack interactions at adjacent holes and onset of multi-site fatigue damage in aging airframes

Focusing on the geometry of one hot spot in airframes, this paper discusses the onset of the interaction of two collinear cracks at adjacent holes and defines the onset as a criterion for multi-site fatigue damage failure. The finite element method is used to calculate the stress intensity factors at the tips of two collinear cracks at adjacent holes growing towards each other. The stress intensity factor is found to increase rapidly at the onset of interaction. Since a rapid increase in stress intensity factor results in a rapid and unstable growth of the crack, the onset of the interaction is proposed as the point where the multi-site fatigue damage starts. A criterion to avoid multi-site fatigue damage locally is then established based on the separation distance of two crack tips at the onset of the interaction. To speed up the simulation of crack growth under multi-site fatigue damage with the finite element method, a semi-empirical criterion is derived to determine the time at which the stress intensity factors at the tips of the cracks correlate. The numerical examples show that the proposed criterion saves simulation time while incurring negligible relative error in the computation of the final crack length.     

Fatigue growth prediction of center slant crack in rectangular plate

This paper presents a numerical prediction model of mixed‐mode crack fatigue growth in a plane elastic plate. It involves a formulations of fatigue growth of multiple crack tips under mixed‐mode loading and a displacement discontinuity method with crack‐tip elements (a boundary element method) proposed recently by Yan is extended to analyse the fatigue growth process of multiple crack tips. Due to an intrinsic feature of the boundary element method, a general growth problem of multiple cracks can be solved in a single‐region formulation. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not necessary. Crack extension is conveniently modelled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. At the same time, the element characters of some related elements are adjusted according to the manner in which the boundary element method is implemented. As an example, the present numerical approach is used to analyse the fatigue growth of a centre slant crack in a rectangular plate. The numerical results illustrate the validation of the numerical prediction model and can reveal the effect of the geometry of the cracked plate on the fatigue growth.     

Dual boundary element method modelling of aircraft structural joints with multiple site damage

The computation of crack growth from a bolt or rivet hole in a structural joint practically requires that the geometry be approximated to some degree. In this paper a simplified quasi-2D stress analysis method, using the boundary element method is presented, where the load transfer rate and the contact stresses at the hole edge for the full 3D geometry are fairly well approximated. Coupled with a dual boundary element formulation for the crack propagation problem, this model is used to evaluate stress intensity factors for through cracks emanating from holes in several double shear lap joint configurations. As the calculated stress intensity factors compare well with experimental data, this procedure is considered to approximate satisfactorily the load transfer rate and the contact stresses at the hole edge of the full 3D geometry, when secondary bending is not a factor.     

A new boundary element method for mixed boundary value problems involving cracks and holes: Interactions between rigid inclusions and cracks

This paper derives a new boundary integral equation (BIE) formulation for plane elastic bodies containing cracks and holes and subjected to mixed displacement/ traction boundary conditions, and proposes a new boundary element method (BEM) based upon this formulation. The basic unknown in the formulation is a complex boundary function H(t), which is a linear combination of the boundary traction and boundary displacement density. The present BIE formulation can be related directly to Muskhelishvili's formalism. Singular interpolation functions of order r ^–1/2 (where r is the distance measured from the crack tip) are introduced such that singular integrand involved at the element level can be integrated analytically. By applying the BEM, the interaction between a rigid circular inclusion and a crack is investigated in details. Our results for the stress intensity factor are comparable with those given by Erdogan and Gupta (1975) and Gharpuray et al. (1990) for a crack emanating from a stiff inclusion, and with those by Erdogan et al. (1974) for a crack in the neighborhood of a stiff inclusion.     

Effect of defective holes on crack using the boundary element method

Boundary value problem of a plate with crack and defect such as the circular and/or elliptical holes is a multi-connected domain problem, this kind of problem is suitable for solving by boundary element method with its higher precision. The sub-region method is used in the paper, a center cracked plate subjected to remote tensile and shear loading is studied numerically. The effect of the circular hole on Mode-I and Mode-II stress intensity factors is studied assuming the plane strain condition, the results are more precise than that of using FEM. The effect of the elliptical hole on Mode-I stress intensity factors is studied as well and some significant results are obtained.     

Analysis of three-dimensional interface cracks using enriched finite elements

Many important interface crack problems are inherently three-dimensional in nature, e.g., debonding of laminated structures at corners and holes. In an effort to accurately analyze three-dimensional interface fracture problems, an efficient computational technique was developed that utilizes enriched crack tip elements containing the correct interface crack tip asymptotic behavior. In the enriched element formulation, the stress intensity factors K _I, K _II, and K _III are treated as additional degrees of freedom and are obtained directly during the finite element solution phase. In this study, the results that should be of greatest interest are obtained for semi-circular surface and quarter-circular corner cracks. Solutions are generated for uniform remote tension and uniform thermal loading, over a wide range of bimaterial combinations. Of particular interest are the free surface effects, and the influence of Dundurs’ material parameters on the strain energy release rate magnitudes and corresponding phase angles.     

An empirical formula for stress intensity factors of cracks emanating from a circular hole in a rectangular plate in tension

An empirical formula of the SIFs of cracks emanating from a circular hole in a rectangular plate in tension is presented and examined. It is found that the empirical formula is simple, yet accurate for evaluating the SIFs of the crack problem.     

Numerical analysis of the inclusion-crack interactions using an integral equation

 A general-purpose integral formulation is proposed for the analysis of the interaction between inclusions and cracks embedded in an elastic isotropic homogeneous infinite medium subjected to a remote loading. This formulation is tailored for the inclusions of arbitrary shapes with the presence of cracks. The discretization is limited to the inclusions (with continuous quadratic triangular and quadrilateral elements) and the cracks (using discontinuous quadratic elements). For the calculation of the stress intensity factors at the crack tips, special crack tip elements are used to model the variation of the displacements near the crack tips. Maximum circumferential stress criterion is adopted to determine the crack propagating direction. Numerical results of benchmark examples are compared with other available methods. Received: 8 January 2002 / Accepted: 24 September 2002     

Cracks propagation and interaction in an orthotropic elastic material: Analytical and numerical methods

An elastic orthotropic material containing a crack in Mode I is considered to formulate a new analytical model. The boundary conditions for the crack existence in the material lead to the solution of the homogeneous Riemann–Hilbert problems. The mathematical model was elaborated for a single and two collinear cracks of different lengths and distance for Mode I in order to investigate cracks interaction problem. Using the theory of Cauchy’s integral and the numerical analysis, the fields in the vicinity of the crack tips were determined.Finite Element Method was applied to compare the mathematical analytical solution and to determine the fields in the vicinity of the crack tips. The critical values of applied stress which caused cracks propagation were evaluated. The interaction of cracks in an orthotropic aramid-epoxy material was studied in details. Comparison of both approaches to crack propagation leads to the conclusion that the new analytical model is correct and can be applied to more complex cracks geometries, including inclined cracks.     

The Interaction between Microcapsules with Different Sizes and Propagating Cracks

The microcapsule-contained self-healing materials are appealing since they can heal the cracks automatically and be effective for a long time. Although many experiments have been carried out, the influence of the size of microcapsules on the self-healing effect is still not well investigated. This study uses the two-dimensional discrete element method (DEM) to investigate the interaction between one microcapsule and one microcrack. The influence of the size of microcapsules is considered. The potential healing time and the influence of the initial damage are studied. The results indicate that the coalescence crack is affected by the size of holes. The elastic modulus, the compressive strength and the coalescence stress decrease with the rising radius of holes. The initial damage in experiments should be greater than 95% of the compressive strength to enhance the self-healing effect. The large microcapsules require slight initial damage. Both a new type of displacement field near the crack and a new category of coalescence crack are observed. The influence of sizes of holes on the cracking behavior of concrete with a circular hole and a pre-existing crack is clarified.     

Stress intensity factors for cracks in thin plates

This paper presents stress intensity factor solutions for several crack configurations in plates. The loadings considered include internal pressure, and also combined bending and tension. The dual boundary element method is used to model the plate and mixed mode stress intensity factors are evaluated by a crack surface displacement extrapolation technique and the J-integral technique. Several cases including centre crack, edge crack and cracks emanating from a hole in finite width plates are presented.     

Interaction between a cracked hole and a line crack under uniform heat flux

This article deals with the interaction between a cracked hole and a line crack under uniform heat flux. Using the principle of superposition, the original problem is converted into three particular cracked hole problems: the first one is the problem of the hole with an edge crack under uniform heat flux, the second and third ones are the problems of the hole under distributed temperature and edge dislocations, respectively, along the line crack surface. Singular integral equations satisfying adiabatic and traction free conditions on the crack surface are obtained for the solution of the second and third problems. The solution of the first problem, as well as the fundamental solutions of the second and third, is obtained by the complex variable method along with the rational mapping function approach. Stress intensity factors (SIFs) at all three crack tips are calculated. Interestingly, the results show that the interaction between the cracked hole and the line crack under uniform heat flux can lead to the vanishing of the SIFs at the hole edge crack tip. The fact has never been seen for the case of a cracked hole and a line crack under remote uniform tension.