Effect of prestrain on fatigue crack propagation in mild steel containing non-metallic inclusions

The effect of tensile prestrain on fatigue crack propagation behaviour of commercial mild steel with significant amount of stringer inclusions has been studied. In prestrained materials the usual stable stage II crack growth region is preceded by a phase wherein a retardation in crack growth rate occurs. No such behaviour is observed in annealed material. The amount of retardation is found to increase with increase in prestrain. A mechanism for the observed retardation in crack growth rate is also presented.

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.

Fatigue crack propagation model and size effect in concrete using dimensional analysis

It is well known that fatigue in concrete causes excessive deformations and cracking leading to structural failures. Due to quasi-brittle nature of concrete and formation of a fracture process zone, the rate of fatigue crack growth depends on a number of parameters, such as, the tensile strength, fracture toughness, loading ratio and most importantly the structural size. In this work, an analytical model is proposed for estimating the fatigue crack growth in concrete by using the concepts of dimensional analysis and including the above parameters. Knowing the governed and the governing parameters of the physical problem and by using the concepts of self-similarity, a relationship is obtained between different parameters involved. It is shown that the proposed fatigue law is able to capture the size effect in plain concrete and agrees well with different experimental results. Through a sensitivity analysis, it is shown that the structural size plays a dominant role followed by loading ratio and the initial crack length in fatigue crack propagation. (C) 2010 Elsevier Ltd. All rights reserved.

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.

Fatigue crack propagation at concrete-concrete bi-material interfaces

This paper presents experimental and analytical studies on fatigue crack propagation in concrete-concrete cold jointed interface specimens. Beams of different sizes having jointed interface between two concretes with different elastic properties are tested under fatigue loading. The acoustic emission technique is used for monitoring the fatigue crack growth. It is observed that the interface having a higher moduli mismatch tends to behave in a brittle manner. The CMOD compliances at different loading cycles are measured and the equivalent crack lengths are determined from a finite element analysis. An analytical model for crack growth rate is proposed using the concepts of the dimensional analysis. (C) 2014 Elsevier Ltd. All rights reserved.

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.

An analysis on overall crack-number-density of short-fatigue-cracks

The evolution of dispersed short-fatigue-cracks is analysed based on the equilibrium of crack-number-density (CND). By separating the mean value and the stochastic fluctuation of local CND, the equilibrium equation of overall CND is derived. Comparing with the mean-field equilibrium equation, the equilibrium equation of overall CND has different forms in the expression of crack-nucleation-rate or crack-growth-rate. The simulation results are compared with experimental measurements showing the stochastic analyses provide consistent tendency with experiments. The discrepancy in simulation results between overall CND and mean-field CND is discussed.

An analysis of collective damage for short fatigue cracks based on equilibrium of crack numerical density

Collective damage of short fatigue cracks was analyzed in the light of equilibrium of crack numerical density. With the estimation of crack growth rate and crack nucleation rate, the solution of the equilibrium equation was studied to reveal the distinct feature of saturation distribution for crack numerical density. The critical time that characterized the transition of short and long-crack regimes was estimated, in which the influences of grain size and grain-boundary obstacle effect were investigated. Furthermore, the total number of cracks and the first order of damage moment were discussed.

Ice induced low temperature fatigue crack propagation in offshore structural steel A131

To investigate the low temperature fatigue crack propagation behavior of offshore structural steel A131 under random ice loading, three ice failure modes that are commonly present in the Bohai Gulf are simulated according to the vibration stress responses induced by real ice loading. The test data are processed by a universal software FCPUSL developed on the basis of the theory of fatigue crack propagation and statistics. The fundamental parameter controlling the fatigue crack propagation induced by random ice loading is determined to be the amplitude root mean square stress intensity factor K-arm. The test results are presented on the crack propagation diagram where the crack growth rate da/dN is described as the function of K-arm. It is evident that the ice failure modes have great influence on the fatigue crack propagation behavior of the steel in ice-induced vibration. However, some of the experimental phenomena and test results are hard to be physically explained at present. The work in this paper is an initial attempt to investigate the cause of collapse of offshore structures due to ice loading.

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.

A probabilistic model for fatigue crack propagation analysis

A simple probabilistic model for predicting crack growth behavior under random loading is presented. In the model, the parameters c and m in the Paris-Erdogan Equation are taken as random variables, and their stochastic characteristic values are obtained through fatigue crack propagation tests on an offshore structural steel under constant amplitude loading. Furthermore, by using the Monte Carlo simulation technique, the fatigue crack propagation life to reach a given crack length is predicted. The tests are conducted to verify the applicability of the theoretical prediction of the fatigue crack propagation.

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.