Effect of plastic zone on the fatigue crack propagation behavior between two fatigue cracks

In this study, the fatigue crack propagation behavior in the stress interaction field between two different fatigue cracks is studied by experiment and finite element analysis. In the experiment, the offset distance between two cracks and the applied stress are varied to create different stress interaction fields. The size of the plastic zone area is used to examine the crack propagation path and rate. Three types of crack propagation in the interaction field were found by experiment, and the crack propagation behavior of two cracks was significantly changed as different stresses were applied. The size of the plastic zone obtained by finite element analysis can be used to explain crack propagation behavior qualitatively.     

Morphology of horizontal cracks in swelling soils

Horizontally formed cracks in the network of cracked swelling soils tend to influence water and solute transport. An approach is suggested for estimating the mean width and volume of horizontal cracks. It is assumed that the nearly horizontal cracks appear as a result of inhomogeneous soil subsidence caused by rapid drying and shrinkage of thin layers at the walls of vertical cracks. Compared with the moist soil matrix, at the same soil depth, horizontal cracks originate as ruptures in stretched layers of the drying walls of vertical cracks. A characteristic of the average inhomogeneity of soil subsidence, i.e., the mean potential relative subsidence (MPRS) depending on the soil depth is defined. It is calculated on the basis of linear shrinkage in the clay soil matrix and at the walls of vertical cracks of different depths, and on two geometrical parameters of crack networks. They are namely the maximum crack depth and the thickness of the upper intensive-cracking layer. The absolute value and sign of the derivative of the MPRS function with respect to soil depth determine the specific volume of horizontal cracks (horizontal-cracks porosity), and their mean width as functions of depth. Model predictions are obtained using published data on variation of linear shrinkage with depth in 19 soil profiles. For lack of data specific to horizontal-crack characteristics model, predictions were compared with data on vertical cracks and subsidence at the soil surface. Satisfactory agreement was obtained for all soil depths up to the maximum crack depth.     

Measurement of energy release rate and energy flux of rapidly bifurcating crack in Homalite 100 and Araldite B by high-speed holographic microscopy

High-speed holographic microscopy is applied to take three successive photographs of fast propagating cracks in Homalite 100 or in Araldite B at the moment of bifurcation. Crack speed at bifurcation is about 540 m/s on Homalite 100, and about 450 m/s on Araldite B. From the photographs, crack speeds immediately before and after bifurcation are obtained, and it is found that discontinuous change of crack speed does not exist at the moment of bifurcation in the case of Homalite 100, but exists in the case of Araldite B. From the photographs, crack opening displacement (COD) is also measured along the cracks as a function of distance r from the crack tips. The measurement results show that the CODs are proportional to √r before bifurcation. After bifurcation, the CODs of mother cracks are proportional to √r, though the CODs of branch cracks are not always proportional to √r. The energy release rate is obtained from the measured CODs, and it is found that energy release rate is continuous at bifurcation point in both cases of Homalite 100 and Araldite B. Energy flux that shows the energy flow toward a crack tip is also obtained.     

The nonlocal solution of two parallel cracks in functionally graded materials subjected to harmonic anti-plane shear waves

In this paper, the dynamic interaction of two parallel cracks in functionally graded materials (FGMs) is investigated by means of the non-local theory. To make the analysis tractable, the shear modulus and the material density are assumed to vary exponentially with the coordinate vertical to the crack. To reduce mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the dynamic problem to obtain stress fields near the crack tips. By use of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations, in which the unknown variables are the jumps of displacements across the crack surfaces. To solve the dual integral equations, the jumps of displacements across the crack surfaces are expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present at the crack tips. The non-local elastic solutions yield a finite hoop stress at the crack tips. The present result provides theoretical references helpful for evaluating relevant strength and preventing material failure of FGMs with initial cracks. The magnitude of the finite stress field depends on relevant parameters, such as the crack length, the distance between two parallel cracks, the parameter describing the FGMs, the frequency of the incident waves and the lattice parameter of materials. The project supported by the National Natural Science Foundation of China (90405016, 10572044) and the Specialized Research Fund for the Doctoral Program of Higher Education (20040213034). The English text was polished by Yunming Chen.     

On the plastic zone size of solids containing doubly periodic rectangular-shaped arrays of cracks under longitudinal shear

Multiple cracks interaction plays an important role in fracture behavior of materials. A number of studies have been devoted to analytical and numerical analyses of the doubly periodic arrays of cracks. A very natural and highly accurate solution procedure is proposed to describe the interaction effect among the doubly periodic rectangular-shaped arrays of cracks. The proposed solution is implemented in the framework of continuously distributed dislocation model and singular integral equation approach. The accuracy of this solution is proved through a comparison of results from the present simulation and known closed form solutions. Further, the interaction effects among the periodic cracks on the plastic zone size and crack tip opening displacement are studied. It is found that the interaction distance among the vertical and horizontal periodic cracks is quite different.     

Effective permeability of cracked unsaturated porous materials

This study focuses on a theoretical estimation of the effective permeability of unsaturated cracked porous media. The closed-form flow solution around and in a superconductive crack, embedded in an infinite porous matrix under a far-field condition, is recalled first. Then the solution of flow around a completely unsaturated (empty) crack that is considered as an obstruction against the flow is determined. The flow solution for partially saturated crack in special configurations is obtained by superposition of the two basic solutions for superconductive and empty cracks. The contribution of an unsaturated crack, with a given saturation degree, to the effective permeability is estimated by using dilute upscaling scheme. Numerical results obtained by Finite Elements Method, are in good agreement with the theoretical results for weak crack densities but show the additional effect of cracks interaction for higher densities.     

A thermomechanical cohesive zone model for bridged delamination cracks

The coupled thermomechanical numerical analysis of composite laminates with bridged delamination cracks loaded by a temperature gradient is described. The numerical approach presented is based on the framework of a cohesive zone model. A traction-separation law is presented which accounts for breakdown of the micromechanisms responsible for load transfer across bridged delamination cracks. The load transfer behavior is coupled to heat conduction across the bridged delamination crack. The coupled crack-bridging model is implemented into a finite element framework as a thermomechanical cohesive zone model (CZM). The fundamental response of the thermomechanical CZM is described. Subsequently, bridged delamination cracks of fixed lengths are studied. Values of the crack tip energy release rate and of the crack heat flux are computed to characterize the loading of the structure. Specimen geometries are considered that lead to crack opening through bending deformation and buckling delamination. The influence of critical mechanical and thermal parameters of the bridging zone on the thermomechanical delamination behavior is discussed. Bridging fibers not only contribute to crack conductance, but by keeping the crack opening small they allow heat flux across the delamination crack to be sustained longer, and thereby contribute to reduced levels of thermal stresses. The micro-mechanism based cohesive zone model allows the assessment of the effectiveness of the individual mechanisms contributing to the thermomechanical crack bridging embedded into the structural analysis.     

Evolution of fatigue crack corrosion from surface irregularities

A moving boundary model is presented for crack nucleation and growth from surface flaws. It concerns with chemical attack that results in material dissolution. A controlling mechanism for evolution is the rupture of a brittle corrosion-protective film that is built up along the corroding surface. The evolution rate is a function of the degree of protective film damage caused by the surface straining. The problem is formulated for an elastic body containing a single and double pits. Low-frequency cyclic loading is considered. Numerical solution is proposed. The behaviours of a growing crack and of two competing cracks are described. Stages of incubation, blunting and steady-state growth characterise a single crack evolution. The steady-state growth rate is found independent of the initial geometry. Stages of independent growth, interactive growth and arrest of one crack characterise the evolution of two competing cracks. The lengths of the arrested cracks are presented as functions of the ratio between the pit depth for a series of different distances between the pits. It is emphasized that the solutions correspond to a homogeneous material. Further work is required to account for the material microstructure.     

Dynamic crack deflection and penetration at interfaces in homogeneous materials: experimental studies and model predictions

We examine the deflection/penetration behavior of dynamic mode-I cracks propagating at various speeds towards inclined weak planes/interfaces of various strengths in otherwise homogeneous isotropic plates. A dynamic wedge-loading mechanism is used to control the incoming crack speeds, and high-speed photography and dynamic photoelasticity are used to observe, in real-time, the failure mode transition mechanism at the interfaces. Simple dynamic fracture mechanics concepts used in conjunction with a postulated energy criterion are applied to examine the crack deflection/penetration behavior and, for the case of interfacial deflection, to predict the crack tip speed of the deflected crack. It is found that if the interfacial angle and strength are such as to trap an incident dynamic mode-I crack within the interface, a failure mode transition occurs. This transition is characterized by a distinct, observable and predicted speed jump as well as a dramatic crack speed increase as the crack transitions from a purely mode-I crack to an unstable mixed-mode interfacial crack.     

DYNAMIC BEHAVIOR OF TWO COLLINEAR PERMEABLE CRACKS IN A PIEZOELECTRIC LAYER BONDED TO TWO HALF SPACES

IntroductionIn the fracture mechanics studies for piezoelectric materials,differently electricboundary conditions at the crack surfaces have been proposed by many researchers.Forexample,for the sake of analytical simplification,the assumption that the cra…     

Coupled field state around three parallel non-symmetric cracks in a piezoelectric/piezomagnetic material plane

In this paper, the behavior of three parallel non-symmetric permeable cracks in a piezoelectric/piezomagnetic material plane subjected to anti-plane shear stress loading was studied by the Schmidt method. The problem was formulated through Fourier transform into three pairs of dual integral equations, in which unknown variables are jumps of displacements across the crack surfaces. To solve the dual integral equations, the jumps of displacements across the crack surfaces were directly expanded as a series of Jacobi polynomials. Finally, the relations among the electric displacement, the magnetic flux and the stress fields near the crack tips can be obtained. The results show that the stress, the electric displacement and the magnetic flux intensity factors at the crack tips depend on the lengths and spacing of cracks. It was also revealed that the crack shielding effect is present in piezoelectric/piezomagnetic materials.     

Multi-interface cracks in a piezoelectric layer bonded to two half-piezoelectric spaces under anti-plane shear loading

The behavior of four parallel symmetry permeable interface cracks in a piezoelectric layer bonded to two half-piezoelectric spaces under anti-plane shear loading is investigated. By using the Fourier transform, the problem can be solved with the help of two pairs of triple integral equations. These equations are solved by the Schmidt method. This process is quite different from that papers adopted previously. The normalized stress and electrical displacement intensity factors are determined for different geometric and property parameters for permeable crack surface conditions. Numerical examples are provided to show the effect of the geometry of the interacting cracks, the thickness and the materials constants of the piezoelectric layer upon the stress and electric displacement intensity factors of the cracks. It is found that the electric displacement intensity factors for the permeable crack surface conditions are much smaller than the results for the impermeable crack surface conditions.     

Steady-state mode I cracks in a viscoelastic triangular lattice

We construct exact solutions for Mode I steady-state cracks in an ideally brittle viscoelastic triangular lattice model. Our analytic solutions for the infinite lattice are compared to numerical results for finite width systems. The issues we address include the crack velocity versus driving curve as well as the onset of additional bond breaking, signaling the emergence of complex spatio-temporal behavior. Somewhat surprisingly, the critical velocity for this transition becomes a decreasing function of the dissipation for sufficiently large values thereof. Lastly, we briefly discuss the possible relevance of our findings for experiments on mode I crack instabilities.     

CTOA and crack-tunneling measurements in thin sheet 2024-T3 aluminum alloy

The stable tearing behavior of 2.3-mm thick sheets of 2024-T3 aluminum alloy was experimentally investigated for middle crack tension specimens having initial flaws that were: (a) flat fatigue cracks (low fatigue stress) and (b) 45-deg through-thickness slant cracks (high fatigue stress). The critical CTOA values during stable tearing were measured by two independent methods, optical microscopy and digital-image correlation. Results from the two methods agree well.The CTOA measurements and observations of the fracture surfaces have shown that the initial stable tearing behavior of low and high fatigue stress tests is significantly different. The cracks in the low fatigue stress tests underwent a transition from flat-to-slant crack growth, during which the CTOA values were high and significant crack tunneling occurred. After crack growth equal to about the thickness (a>B), CTOA reached a constant value of 6 deg and after crack growth equal to about twice the thickness (a>2B), crack tunneling stabilized. The cracks in the high fatigue stress tests reach the same constant CTOA value after crack growth equal to about the thickness, but produced only slightly higher CTOA values during initial crack growth. The amount of tunneling in the high fatigue stress tests was about the same as that in the low fatigue stress tests after the flat-to-slant transition.This study indicates that stress history has an influence on the initial portion of the stable tearing behavior. The initial high CTOA values, in the low fatigue crack tests, coincided with large three-dimensional crack front shape changes due to a variation in the through-thickness crack-tip constraint. The measured CTOA reached a constant value of 6 deg for crack growth of about the specimen thickness. This coincided with the onset of 45-deg slant crack growth and a stabilized, slightly tunneled (about 20 percent of the thickness) crack-front shape. For crack growth on the 45-deg slant, the crack front and local field variables are still highly three dimensional.     

New stress intensity factor solutions of a periodic array of cracks in residual stress field

Many important applications of crack mechanics involve self-equilibrating residual or thermal stress fields. For these types of problems, the traditional fracture mechanics approach based on the superposition principle has ignored the effect of crack surface contact when the crack-tip propagates into the residual compressive region. Contact between the crack faces and the wedging action are responsible for subsequent crack-tip reopening, which often leads to a much larger mode I stress intensity factor. In this study, an analytical approach is used to study the effect of crack face contact for a period array of collinear cracks embedded in several typical residual stress fields. It is found that the nonlinear contact between crack surfaces dominates the cracking behavior in residual/thermal stress fields, which is responsible for crack coalescence.     

Interactive effect of multiple crack growth on fatigue

Interaction of multiple cracks are found in mud when drying and crazing, thin film varnishes or coatings of aeroengine turbine blades. A two-dimensional multiple crack interaction model is developed to simulate the growth of interacting parallel surface cracks. Density and the initial distribution of the microcracks are accounted for in analyzing the growth of a crack to a pre-determined length. Analytical predictions are discussed with reference to experimental observations of fatigue cracks on coated turbine blades. Introduction of a large density of similar cracks can enhance the fatigue life of structural components.     

Fracture of Edge Cracked Layered Plates Subjected to In-plane Bending

Layered structures are used in protection systems such as personal and heavy armor, windshields and also in functionally graded thermal barriers. The focus of this study is to understand the behavior of cracks in such systems, especially when the crack orientation is such that there are property changes along the crack front. Layered plates were prepared by bonding together sheets of epoxy and Polymethylmethacrylate (PMMA) using an epoxy adhesive. Among the two, epoxy has higher elastic modulus and lower fracture toughness compared to PMMA. Two different layer configurations; a single sheet of epoxy bonded to a single PMMA sheet (two layer) and a single sheet of epoxy sandwiched between two PMMA sheets (three layer) were considered. Single edge notched specimens were loaded in three point bending and the thickness averaged stress intensity factor (SIF) was estimated through photoelasticity. Subsequently, the behavior of crack-extension in these plates was also investigated. In both configurations, crack growth initiated in the epoxy layer first and extended stably before the start of crack extension in the PMMA layer. Once the crack extension started in the PMMA layer, the plate looses its structural integrity. It was observed that the onset of crack extension in the epoxy layer can be predicted using the thickness averaged SIF. A method using analysis of patched cracks is presented for estimating the load at which the plate completely looses its load carrying capacity. The estimates from this method match well with experimental results.     

Mixed mode axisymmetric annular cracks in a finite layer: Off angle crack initiation

The solutions of axisymmetric Volterra type climb and glide edge dislocations are obtained in a layer by means of the Hankel transforms. Utilizing the same procedure, Green’s function solution is obtained for a layer under self-equilibration normal ring traction. The distributed dislocation technique is used to construct integral equations for a system of co-axial annular cracks where the layer is under axisymmetric normal loads. These equations are solved numerically to obtain dislocation density on the cracks surfaces. The results are employed to determine stress intensity factors for annular and penny-shaped cracks and the interaction between two co-axial penny-shaped cracks is studied. Moreover, the stress intensity factors of the interacting cracks are determined such that they can be further used in conjunction with strain energy density (SED) failure criterion to obtain the possible direction of crack initiation that may not be apparent under mixed mode conditions.     

Analysis of cracked functionally graded piezoelectric strip

The fracture behavior of a cracked strip under antiplane mechanical and inplane electrical loading is studied. A functionally graded piezoelectric strip with exponential material gradation is under consideration. The mechanical and electrical loading is combined via loading coupling factor. The problem of a graded piezoelectric strip containing a screw dislocation is solved. This solution results in stress and electric displacement components with Cauchy singularity. Based on the solution achieved for the dislocation, the distributed dislocation technique (DDT) is utilized to form any geometry of multiple cracks and analyze the behavior of a cracked strip under antiplane mechanical and inplane electrical loading. This technique is capable of the analysis of a strip with a system of interacting cracks. Several examples including strips with single crack, two straight cracks and two curved cracks are presented.