On crack path stability in a layered material

Crack paths in an elastic layer on top of a substrate are considered. Crack growth is initiated from an edge crack in the layer. The plane of the initially straight crack forms an angle to the free surface. The load consists of a pair of forces applied at the crack mouth and parallel to the interface. Crack paths are calculated using a boundary element method. Crack growth is assumed to proceed along a path for which the mode II stress intensity factor vanishes. The inclination and the length of the initial crack are varied. The effect of two different substrates on the crack path evolution is demonstrated. A crack path initially leading perpendicularly to the interface is shown to be directionally unstable for a rigid substrate. Irrespective of its initial angle, the crack does not reach the interface, but reaches the free surface if the layer is infinitely long. At finite layer length the crack reaches the upper free surface if the initial crack inclination to the surface is small enough. For an inextendable flexible substrate, on the other hand, the crack reaches the interface if its initial inclination is large enough. For the flexible substrate an unstable path parallel with the sides of an infinitely long layer is identified. The results are compared with experimental results and discussed in view of characterisation of directionally unstable crack paths. The energy release rate for an inclined edge crack is determined analytically.

A local-global scheme for tracking crack path in three-dimensional solids

This paper aims to contribute to the three-dimensional generalization of numerical prediction of crack propagation through the formulation of finite elements with embedded discontinuities. The analysis of crack propagation in two-dimensional problems yields lines of discontinuity that can be tracked in a relatively simple way through the sequential construction of straight line segments oriented according to the direction of failure within each finite element in the solid. In three-dimensional analysis, the construction of the discontinuity path is more complex because it requires the creation of plane surfaces within each element, which must be continuous between the elements. In the method proposed by Chaves (2003) the crack is determined by solving a problem analogous to the heat conduction problem, established from local failure orientations, based on the stress state of the mechanical problem. To minimize the computational effort, in this paper a new strategy is proposed whereby the analysis for tracking the discontinuity path is restricted to the domain formed by some elements near the crack surface that develops along the loading process. The proposed methodology is validated by performing three-dimensional analyses of basic problems of experimental fractures and comparing their results with those reported in the literature.

Short fatigue crack path determinants in polycrystalline Ni-base superalloys

Fifty seven short fatigue cracks in the Ni-base superalloy AP1 have been examined, to ascertain how the paths taken by growing fatigue cracks are determined. The observations were made on the surface of a smooth specimen, and on the exposed fracture surfaces. Three dimensional reconstructions of the vulnerable microstructures in the vicinity of the cracks were produced. Initiation occurred in mode II, with the lines of intersection of the initiation sites with the specimen top surface orientated at approximately 45° to the tensile axis. These initiation sites developed in slip bands which crossed a large grain and at least one other grain via a grain boundary with a low angle of misorientation. 'River markings' on one of the initiation facets, indicated that the crack first opened from the top centre of the initiation grain. Subsequent to initiation, the growth paths of these cracks are related to the misorientations of the grains and the progress of the crack front.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

2D evaluation of crack openings using smeared and embedded crack models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Microstructural effects on fatigue crack growth behavior of a microalloyed steel

Thermal transformations on microalloyed steels can produce multiphase microstructures with different amounts of ferrite, martensite, bainite and retained austenite. These different phases, with distinct morphologies, are determinant of the mechanical behavior of the steel and can, for instance, affect the crack path or promote crack shielding, thus resulting in changes on its propagation rate under cyclic loading. The aim of the present work is to evaluate the effects of microstructure on the tensile strength and fatigue crack growth (FCG) behaviour of a 0.08%C-1,5%Mn (wt. pct.) microalloyed steel, recently developed by a Brazilian steel maker under the designation of RD480. This steel is being considered as a promising alternative to replace low carbon steel in wheel components for the automotive industry. Various microstructural conditions were obtained by means of heat treatments followed by water quench, in which the material samples were kept at the temperatures of 800, 950 and 1200 °C. In order to describe the FCG behavior, two models were tested: the conventional Paris equation and a new exponential equation developed for materials showing non-linear FCG behavior. The results allowed correlating the tensile properties and crack growth resistance to the microstructural features. It is also shown that the Region II FCG curves of the dual and multiphase microstructural conditions present crack growth transitions that are better modeled by dividing them in two parts. The fracture surfaces of the fatigued samples were observed via scanning electron microscopy in order to reveal the fracture mechanisms presented by the various material conditions. © 2010 Published by Elsevier Ltd.

Laser shock peening treatment to control and moderate fatigue crack growth in aircraft structure based on residual stress engineering approach

Laser Shock Peening (LSP) is a surface enhancement treatment which induces a significant layer of beneficial compressive residual stresses of up to several mm underneath the surface of metal components in order to improve the detrimental effects of the crack growth behavior rate in it. The aim of this thesis is to predict the crack growth behavior in metallic specimens with one or more stripes which define the compressive residual stress area induced by the Laser Shock Peening treatment. The process was applied as crack retardation stripes perpendicular to the crack propagation direction with the object of slowing down the crack when approaching the peened stripes. The finite element method has been applied to simulate the redistribution of stresses in a cracked model when it is subjected to a tension load and to a compressive residual stress field, and to evaluate the Stress Intensity Factor (SIF) in this condition. Finally, the Afgrow software is used to predict the crack growth behavior of the component following the Laser Shock Peening treatment and to detect the improvement in the fatigue life comparing it to the baseline specimen. An educational internship at the “Research & Technologies Germany – Hamburg” department of AIRBUS helped to achieve knowledge and experience to write this thesis. The main tasks of the thesis are the following: •To up to date Literature Survey related to “Laser Shock Peening in Metallic Structures” •To validate the FE model developed against experimental measurements at coupon level •To develop design of crack growth slowdown in Centered Cracked Tension specimens based on residual stress engineering approach using laser peened strip transversal to the crack path •To evaluate the Stress Intensity Factor values for Centered Cracked Tension specimens after the Laser Shock Peening treatment via Finite Element Analysis •To predict the crack growth behavior in Centered Cracked Tension specimens using as input the SIF values evaluated with the FE simulations •To validate the results by means of experimental tests

Effect of Laser Shock Peening on Fatigue Crack Propagation of Aeronautical Structures

Laser Shock Peening (LSP) is a surface enhancement treatment which induces a significant layer of beneficial compressive residual stresses up to several mm underneath the surface of metal components in order to improve the detrimental effects of crack growth behavior rate in it. The aim of this thesis is to predict the crack growth behavior of thin Aluminum specimens with one or more LSP stripes defining a compressive residual stress area. The LSP treatment has been applied as crack retardation stripes perpendicular to the crack growing direction, with the objective of slowing down the crack when approaching the LSP patterns. Different finite element approaches have been implemented to predict the residual stress field left by the laser treatment, mostly by means of the commercial software Abaqus/Explicit. The Afgrow software has been used to predict the crack growth behavior of the component following the laser peening treatment and to detect the improvement in fatigue life comparing to the specimen baseline. Furthermore, an analytical model has been implemented on the Matlab software to make more accurate predictions on fatigue life of the treated components. An educational internship at the Research and Technologies Germany- Hamburg department of Airbus helped to achieve knowledge and experience to write this thesis. The main tasks of the thesis are the following: -To up to date Literature Survey related to laser shock peening in metallic structures -To validate the FE models developed against experimental measurements at coupon level -To develop design of crack growth slow down in centered and edge cracked tension specimens based on residual stress engineering approach using laser peened patterns transversal to the crack path -To predict crack growth behavior of thin aluminum panels -To validate numerical and analytical results by means of experimental tests.

The effect of frequency and microstructure on corrosion fatigue crack propagation in high strength aluminium alloys

Fatigue crack growth in high strength aluminium alloy 7150 commercial plate material has been studied in both laboratory air and acidified aqueous salt solution. The aggressive aqueous environment enhanced fatigue crack growth rates by up to an order in magnitude compared to laboratory air. The enhancement in fatigue crack growth rate was accompanied by evidence of embrittlement in the crack path, involving both brittle intergranular and transgranular failure modes. Both the enhancement of fatigue crack growth rates and the extent of intergranular growth modes are dependent on cyclic frequency which, along with the absence of a similar frequency effect in a spray-formed version of the material with a significantly different grain structure, supports a mechanism of grain boundary hydrogen diffusion for intergranular corrosion fatigue crack growth. The convergence of corrosion fatigue crack growth rates at high ΔK in both spray-formed and conventional plate materials coincides with the operation of identical transgranular corrosion fatigue modes dependent on strain-controlled hydrogen diffusion ahead of the crack tip. © 1997 Acta Metallurgica Inc.

Modelling short crack growth behaviour in nickel-base superalloys

This paper provides a description of the features and mechanisms of facetted short crack growth in Ni-base superalloys, and briefly reviews existing short crack growth models in terms of their application to Ni-base alloys. The concept of “soft barriers” is introduced to produce a new two-phase model for local microstructural effects on short crack growth in Waspaloy. This is derived from detailed observations of crack growth through individual grains. The model differs from all previous approaches in highlighting the importance of crack path perturbations within grains. Potential applications of the model in alloy development are discussed.

A novel two-stage discrete crack method based on the screened Poisson equation and local mesh refinement

We propose an alternative crack propagation algo- rithm which effectively circumvents the variable transfer procedure adopted with classical mesh adaptation algo- rithms. The present alternative consists of two stages: a mesh-creation stage where a local damage model is employed with the objective of defining a crack-conforming mesh and a subsequent analysis stage with a localization limiter in the form of a modified screened Poisson equation which is exempt of crack path calculations. In the second stage, the crack naturally occurs within the refined region. A staggered scheme for standard equilibrium and screened Poisson equa- tions is used in this second stage. Element subdivision is based on edge split operations using a constitutive quantity (damage). To assess the robustness and accuracy of this algo- rithm, we use five quasi-brittle benchmarks, all successfully solved.

Finite element simulations of notch tip fields in magnesium single crystals

Recent experiments using three point bend specimens of Mg single crystals have revealed that tensile twins of {10 (1) over bar2}-type form profusely near a notch tip and enhance the fracture toughness through large plastic dissipation. In this work, 3D finite element simulations of these experiments are carried out using a crystal plasticity framework which includes slip and twinning to gain insights on the mechanics of fracture. The predicted load-displacement curves, slip and tensile twinning activities from finite element analysis corroborate well with the experimental observations. The numerical results are used to explore the 3D nature of the crack tip stress, plastic slip and twin volume fraction distributions near the notch root. The occurrence of tensile twinning is rationalized from the variation of normal stress ahead of the notch tip. Further, deflection of the crack path at twin-twin intersections observed in the experiments is examined from an energy standpoint by modeling discrete twins close to the notch root.

Collective evolution characteristics and computer simulation of short fatigue cracks

Fatigue testing was conducted using a kind of triangular isostress specimen to obtain the short-fatigue-crack behaviour of a weld low-carbon steel. The experimental results show that short cracks continuously initiate at slip bands within ferrite grain domains and the crack number per unit area gradually increases with increasing number of fatigue cycles. The dispersed short cracks possess an orientation preference, which is associated with the crystalline orientation of the relevant slip system. Based on the observed collective characteristics, computer modelling was carried out to simulate the evolution process of initiation, propagation and coalescence of short cracks. The simulation provides progressive displays which imitate the appearance of experimental observations. The results of simulation indicate that the crack path possesses a stable value of fractal dimension whereas the critical value of percolation covers a wide datum band, suggesting that the collective evolution process of short cracks is sensitive to the pattern of crack site distribution.