Dynamic growth of tensile cracks by ductile and brittle fracture mechanisms in a viscoplastic material

In this paper, a finite element analysis of steady-state dynamic crack growth under Mode I, plane strain, small-scale yielding conditions is performed in a rate dependent plastic material characterized by the over-stress model. The main objective of the paper is to obtain theoretically the dependence of dynamic fracture toughness on crack speed. Crack propagation due to a ductile (micro-void) mechanism or a brittle (cleavage) mechanism, as well as transition from one mode to another are considered. The conversion from ductile to brittle has been observed experimentally but has received very little attention using analytical methods. Local fracture criteria based on strains and stresses are used to describe ductile and brittle fracture mechanisms. The results obtained in this paper are in general agreement with micro-structural observations of mode conversion during fracture initiation. Finally, the particular roles played by material rate sensitivity and inertia are examined in some detail.

Effect of strain rate on the high-temperature low-cycle fatigue properties of a nimonic PE-16 superalloy

Strain-rate effects on the low-cycle fatigue (LCF) behavior of a NIMONIC PE-16 superalloy have been evaluated in the temperature range of 523 to 923 K. Total-strain-controlled fatigue tests were per-formed at a strain amplitude of +/-0.6 pct on samples possessing two different prior microstructures: microstructure A, in the solution-annealed condition (free of gamma' and carbides); and microstructure B, in a double-aged condition with gamma' of 18-nm diameter and M23C6 carbides. The cyclic stress response behavior of the alloy was found to depend on the prior microstructure, testing temperature, and strain rate. A softening regime was found to be associated with shearing of ordered gamma' that were either formed during testing or present in the prior microstructure. Various manifestations of dynamic strain aging (DSA) included negative strain rate-stress response, serrations on the stress-strain hysteresis loops, and increased work-hardening rate. The calculated activation energy matched well with that for self-diffusion of Al and Ti in the matrix. Fatigue life increased with an increase in strain rate from 3 x 10(-5) to 3 x 10(-3) s-1, but decreased with further increases in strain rate. At 723 and 823 K and low strain rates, DSA influenced the deformation and fracture behavior of the alloy. Dynamic strain aging increased the strain localization in planar slip bands, and impingement of these bands caused internal grain-boundary cracks and reduced fatigue life. However, at 923 K and low strain rates, fatigue crack initiation and propagation were accelerated by high-temperature oxidation, and the reduced fatigue life was attributed to oxidation-fatigue interaction. Fatigue life was maximum at the intermediate strain rates, where strain localization was lower. Strain localization as a function of strain rate and temperature was quantified by optical and scanning electron microscopy and correlated with fatigue life.

Evaluation of Fracture Parameters by Double-G, Double-K Models and Crack Extension Resistance for High Strength and Ultra High Strength Concrete Beams

This paper presents the advanced analytical methodologies such as Double- G and Double - K models for fracture analysis of concrete specimens made up of high strength concrete (HSC, HSC1) and ultra high strength concrete. Brief details about characterization and experimentation of HSC, HSC1 and UHSC have been provided. Double-G model is based on energy concept and couples the Griffith's brittle fracture theory with the bridging softening property of concrete. The double-K fracture model is based on stress intensity factor approach. Various fracture parameters such as cohesive fracture toughness (4), unstable fracture toughness (K-Ic(c)), unstable fracture toughness (K-Ic(un)) and initiation fracture toughness (K-Ic(ini)) have been evaluated based on linear elastic fracture mechanics and nonlinear fracture mechanics principles. Double-G and double-K method uses the secant compliance at the peak point of measured P-CMOD curves for determining the effective crack length. Bi-linear tension softening model has been employed to account for cohesive stresses ahead of the crack tip. From the studies, it is observed that the fracture parameters obtained by using double - G and double - K models are in good agreement with each other. Crack extension resistance has been estimated by using the fracture parameters obtained through double - K model. It is observed that the values of the crack extension resistance at the critical unstable point are almost equal to the values of the unstable fracture toughness K-Ic(un) of the materials. The computed fracture parameters will be useful for crack growth study, remaining life and residual strength evaluation of concrete structural components.

Failure modes and fatigue behavior of conventional and refilled friction stir spot welds in AA 6061-T6 sheets

The fatigue behavior of conventional friction stir spot welding (FSSW) and friction stir spot welding refilled by the friction forming process (FSSW-FFP) in aluminum 6061-T6 lap shear specimens, are investigated based on the experimental observations. Optical micrographs of the welds after fatigue failure in both the cases are examined to study the fatigue crack propagation and failure modes. Experimental results indicate that the fatigue strength of the FSSW-FFP weld samples is higher than that of the conventional FSSW samples at all loads. Fracture surfaces are analyzed in detail using the scanning electron microscope. (C) 2013 Elsevier Ltd. All rights reserved.

Static shear strength and fatigue life of refill friction stir spot welded AA 6061-T6 sheets

This study was aimed at evaluating the static shear strength and fatigue properties of the newly developed refilled friction stir spot welded AA 6061-T6 joints. The keyhole, the process disadvantage of conventional friction stir spot welding, was refilled successfully, using an additional filler plate, with specially designed tools. Two different tool profiles, namely, convex and concave, were used for the refilling process. Sound and defect free joints were obtained by the refilling process. Joints refilled with convex tools showed better static shear strength than those with the concave ones. The variation of microhardness in different regions of the weld was analysed. Fatigue tests were conducted on the lap shear specimens at a stress ratio of R=0.1. The optical micrographs of the welds after fatigue failure in both the conventional and refilled processes were examined to study the fatigue crack propagation and failure modes.

Crack stability in edge-notched clamped beam specimens: modeling and experiments

Stability of a fracture toughness testing geometry is important to determine the crack trajectory and R-curve behavior of the specimen. Few configurations provide for inherent geometric stability, especially when the specimen being tested is brittle. We propose a new geometrical construction called the single edge notched clamped bend specimen (SENCB), a modified form of three point bending, yielding stable cracking under load control. It is shown to be particularly suitable for small-scale structures which cannot be made free-standing, (e.g., thin films, coatings). The SENCB is elastically clamped at the two ends to its parent material. A notch is inserted at the bottom center and loaded in bending, to fracture. Numerical simulations are carried out through extended finite element method to derive the geometrical factor f(a/W) and for different beam dimensions. Experimental corroborations of the FEM results are carried out on both micro-scale and macro-scale brittle specimens. A plot of vs a/W, is shown to rise initially and fall off, beyond a critical a/W ratio. The difference between conventional SENB and SENCB is highlighted in terms of and FEM simulated stress contours across the beam cross-section. The `s of bulk NiAl and Si determined experimentally are shown to match closely with literature values. Crack stability and R-curve effect is demonstrated in a PtNiAl bond coat sample and compared with predicted crack trajectories from the simulations. The stability of SENCB is shown for a critical range of a/W ratios, proving that it can be used to get controlled crack growth even in brittle samples under load control.

Crack propagation in graphene

The crack initiation and growth mechanisms in an 2D graphene lattice structure are studied based on molecular dynamics simulations. Crack growth in an initial edge crack model in the arm-chair and the zig-zag lattice configurations of graphene are considered. Influence of the time steps on the post yielding behaviour of graphene is studied. Based on the results, a time step of 0.1 fs is recommended for consistent and accurate simulation of crack propagation. Effect of temperature on the crack propagation in graphene is also studied, considering adiabatic and isothermal conditions. Total energy and stress fields are analyzed. A systematic study of the bond stretching and bond reorientation phenomena is performed, which shows that the crack propagates after significant bond elongation and rotation in graphene. Variation of the crack speed with the change in crack length is estimated. (C) 2015 AIP Publishing LLC.

Simulations and Experiments of Stochastic Characteristics for Collective Short Fatigue Cracks in Steels

Stochastic characteristics prevail in the process of short fatigue crack progression. This paper presents a method taking into account the balance of crack number density to describe the stochastic behaviour of short crack collective evolution. The results from the simulation illustrate the stochastic development of short cracks. The experiments on two types of steels show the random distribution for collective short cracks with the number of cracks and the maximum crack length as a function of different locations on specimen surface. The experiments also give the variation of total number of short cracks with fatigue cycles. The test results are consistent with numerical simulations.

Finite Element Solutions for Plane Strain Mode I Crack with Strain Grading Effect

In this paper, a new phenomenological theory with strain gradient effects is proposed to account for the size dependence of plastic deformation at micro- and submicro-length scales. The theory fits within the framework of general couple stress theory and three rotational degrees of freedom omega(i) are introduced in addition to the conventional three translational degrees of freedom mu(i). omega(i) is called micro-rotation and is the sum of material rotation plus the particles' relative rotation. While the new theory is used to analyze the crack tip field or the indentation problems, the stretch gradient is considered through a new hardening law. The key features of the theory are that the rotation gradient influences the material character through the interaction between the Cauchy stresses and the couple stresses; the term of stretch gradient is represented as an internal variable to increase the tangent modulus. In fact the present new strain gradient theory is the combination of the strain gradient theory proposed by Chen and Wang (Int. J. Plast., in press) and the hardening law given by Chen and Wang (Acta Mater. 48 (2000a) 3997). In this paper we focus on the finite element method to investigate material fracture for an elastic-power law hardening solid. With remotely imposed classical K fields, the full field solutions are obtained numerically. It is found that the size of the strain gradient dominance zone is characterized by the intrinsic material length l(1). Outside the strain gradient dominance zone, the computed stress field tends to be a classical plasticity field and then K field. The singularity of stresses ahead of the crack tip is higher than that of the classical field and tends to the square root singularity, which has important consequences for crack growth in materials by decohesion at the atomic scale. (C) 2002 Elsevier Science Ltd. All rights reserved.

Interface crack problems with strain gradient effects

In this paper, the strain gradient theory proposed by Chen and Wang (2001 a, 2002b) is used to analyze an interface crack tip field at micron scales. Numerical results show that at a distance much larger than the dislocation spacing the classical continuum plasticity is applicable; but the stress level with the strain gradient effect is significantly higher than that in classical plasticity immediately ahead of the crack tip. The singularity of stresses in the strain gradient theory is higher than that in HRR field and it slightly exceeds or equals to the square root singularity and has no relation with the material hardening exponents. Several kinds of interface crack fields are calculated and compared. The interface crack tip field between an elastic-plastic material and a rigid substrate is different from that between two elastic-plastic solids. This study provides explanations for the crack growth in materials by decohesion at the atomic scale.

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.

Crack propagation in structures subjected to periodic excitation

In the present paper, a simple mechanical model is developed to predict the dynamic response of a cracked structure subjected to periodic excitation, which has been used to identify the physical mechanisms in leading the growth or arrest of cracking. The structure under consideration consists of a beam with a crack along the axis, and thus, the crack may open in Mode I and in the axial direction propagate when the beam vibrates. In this paper, the system is modeled as a cantilever beam lying on a partial elastic foundation, where the portion of the beam on the foundation represents the intact portion of the beam. Modal analysis is employed to obtain a closed form solution for the structural response. Crack propagation is studied by allowing the elastic foundation to shorten (mimicking crack growth) if a displacement criterion, based on the material toughness, is met. As the crack propagates, the structural model is updated using the new foundation length and the response continues. From this work, two mechanisms for crack arrest are identified. It is also shown that the crack propagation is strongly influenced by the transient response of the structure.

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.

Experimental relaxation behavior of chopped-strand mat glass-fiber reinforced polyester beams in 3 point bending

A dimensionless relation of the form for collating fatigue crack starting growth data is proposed in which Δkth represents the stress intensity factor range at the threshold. Based on experimental results, this relation attains the value of 0.6 for a fatigue crack to start growth in the Austenitic stainless steel investigated in this work. Metallurgical examinations were also carried out to show a transgranular shear mode of cyclic cleavage and plastic shear.

Ductile crack propagation with the strain-energy density criterion

The strain energy density criterion is used to characterize subcritical crack growth in a thin aluminum alloy sheet undergoing general yielding. A finite element analysis which incorporates both material and geometrical nonlinear behaviors of the cracked sheets is developed to predict fracture loads at varying crack growth increments. The predicted results are in excellent agreement with those measured experimentally, thus confirming the validity of the strain energy density criterion for characterizing ductile crack propagation.