A thermodynamic framework for fatigue crack growth in concrete

The phenomenon of fatigue is commonly observed in majority of concrete structures and it is important to mathematically model it in order to predict their remaining life. An energy approach is adopted in this research by using the framework of thermodynamics wherein the dissipative phenomenon is described by a dissipation potential. An analytical expression is derived for the dissipation potential using the concepts of dimensional analysis and self-similarity to describe a fatigue crack propagation model for concrete. This is validated using available experimental results. Through a sensitivity analysis, the hierarchy of importance of different parameters is highlighted.

Analysis of fatigue crack growth in reinforced concrete beams

Mechanical behavior of reinforced concrete members is influenced by the action of unknown crack bridging reactions of rebars. Under cyclic loading, due to progressive growth of cracks, this bridging action contributes to the overall strength, stiffness and hysteretic behavior of the member. In this work, fatigue behavior of reinforced concrete beams are studied using a crack propagation law, developed using dimensional analysis for plain concrete with the effect of reinforcement being simulated through constraint exerted on the crack opening. The parameters considered in the model are fracture toughness, crack length, loading ratio and structural size. A numerical procedure is followed to compute fatigue life of RC beams and the dissipated energy in the steel reinforcement due to the shake down phenomenon under cyclic loading. Through a sensitivity study, it is concluded that the structural size is the most sensitive parameter in the fatigue crack propagation phenomenon. Furthermore, the residual moment carrying capacity of an RC member is determined as a function of crack extension by including the bond-slip mechanism.

Monte-Carlo simulation of surface crack growth rate for offshore structural steel E36-Z35

The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.

Dislocations emission and crack extension at the atomistic crack tip in body-centered-cubic metal Mo

The behaviors of a crack in body-centered-cubic metal Mo under different loading modes were studied using the molecular dynamics method. Dislocation emission was observed near the crack tip in response to mode II loading with theta = 0 degrees in which theta is the inclination angle of the slip plane with respect to the crack plane, and two full dislocations were observed at the stress level of K-II = 1.17 MPa m(1/2) without any evidence of crack extension. Within the range of 0 degrees less than or equal to theta less than or equal to 45 degrees, crack extension was observed in response to mode I loading, and the effect of crystal orientation on the crack propagation was studied, The crack propagated along the [111] slip direction without any evidence of dislocations emission.

I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

Dynamic behaviors of mode III interfacial crack under a constant loading rate

In an earlier study on intersonic crack propagation, Gao et al. (J. Mech. Phys. Solids 49: 2113-2132, 2001) described molecular dynamics simulations and continuum analysis of the dynamic behaviors of a mode II dominated crack moving along a weak plane under a constant loading rate. The crack was observed to initiate its motion at a critical time after the onset of loading, at which it is rapidly accelerated to the Rayleigh wave speed and propagates at this speed for a finite time interval until an intersonic daughter crack is nucleated at a peak stress at a finite distance ahead of the original crack tip. The present article aims to analyze this behavior for a mode III crack moving along a bi-material interface subject to a constant loading rate. We begin with a crack in an initially stress-free bi-material subject to a steadily increasing stress. The crack initiates its motion at a critical time governed by the Griffith criterion. After crack initiation, two scenarios of crack propagation are investigated: the first one is that the crack moves at a constant subsonic velocity; the second one is that the crack moves at the lower shear wave speed of the two materials. In the first scenario, the shear stress ahead of the crack tip is singular with exponent -1/2, as expected; in the second scenario, the stress singularity vanishes but a peak stress is found to emerge at a distance ahead of the moving crack tip. In the latter case, a daughter crack supersonic with respect to the softer medium can be expected to emerge ahead of the initial crack once the peak stress reaches the cohesive strength of the interface.

Element-Free Galerkin Modelling of Crack Migration in Composite Laminates

The development of a virtual testing environment, as a cost-effective industrial design tool in the design and analysis of composite structures, requires the need to create models efficiently, as well as accelerate the analysis by reducing the number of degrees of freedom, while still satisfying the need for accurately tracking the evolution of a debond, delamination or crack front. The eventual aim is to simulate both damage initiation and propagation in components with realistic geometrical features, where crack propagation paths are not trivial. Meshless approaches, and the Element-Free Galerkin (EFG) method, are particularly suitable for problems involving changes in topology and have been successfully applied to simulate damage in homogeneous materials and concrete. In this work, the method is utilized to model initiation and mixed-mode propagation of cracks in composite laminates, and to simulate experimentally-observed crack migration which is difficult to model using standard finite element analysis. N

Describing fatigue crack growth and load ratio effects in Al 2524 T3 alloy with an enhanced exponential model

The fatigue crack behavior in metals and alloys under constant amplitude test conditions is usually described by relationships between the crack growth rate da/dN and the stress intensity factor range Delta K. In the present work, an enhanced two-parameter exponential equation of fatigue crack growth was introduced in order to describe sub-critical crack propagation behavior of Al 2524-T3 alloy, commonly used in aircraft engineering applications. It was demonstrated that besides adequately correlating the load ratio effects, the exponential model also accounts for the slight deviations from linearity shown by the experimental curves. A comparison with Elber, Kujawski and "Unified Approach" models allowed for verifying the better performance, when confronted to the other tested models, presented by the exponential model. (C) 2012 Elsevier Ltd. All rights reserved.

Fatigue crack growth of Al 7475 T761 under constant and variable amplitude loading

Experimental programs in constant and variable amplitude loading were performed to obtain a x N curves and to study retardation in fatigue crack growth due to overloads. The main aim of this research program was to analyse the effect of overload ratio and number of overload peaks. The effect of underloads, before and after the overload blocks was also studied. The generalised equation of Paris-Erdogan type was used for modelling of obtained data on crack propagation under constant amplitude load.

Embedded crack elements with non-uniform discontinuity modes

The consequences of the use of embedded crack finite elements with uniform discontinuity modes (opening and sliding) to simulate crack propagation in concrete are investigated. It is shown the circumstances in which the consideration of uniform discontinuity modes is not suitable to accurately model the kinematics induced by the crack and must be avoided. It is also proposed a technique to embed cracks with non-uniform discontinuity modes into standard displacement-based finite elements to overcome the shortcomings of the uniform discontinuity modes approach.

Mechanical strength and subcritical crack growth under wet cyclic loading of glass-infiltrated dental ceramics

Objectives. Evaluate the flexural strength (sigma) and subcritical crack growth (SCG) under cyclic loading of glass-infiltrated alumina-based (IA, In-Ceram Alumina) and zirconia-reinforced (IZ, In-Ceram Zirconia) ceramics, testing the hypothesis that wet environment influences the SCG of both ceramics when submitted to cyclic loading.Methods. Bar-shaped specimens of IA (n = 45) and IZ ( n = 45) were fabricated and loaded in three-point bending (3P) in 37 degrees C artificial saliva (IA(3P) and IZ(3P)) and cyclic fatigued (F) in dry (D) and wet (W) conditions (IA(FD), IA(FW), IZ(FD), IZ(FW)). The initial sigma and the number of cycles to fracture were obtained from 3P and F tests, respectively. Data was examined using Weibull statistics. The SCG behavior was described in terms of crack velocity as a function of maximum stress intensity factor (K(Imax)).Results. The Weibull moduli (m = 8) were similar for both ceramics. The characteristic strength (sigma(0)) of IA and IZ was and 466 MPa 550 MPa, respectively. The wet environment significantly increased the SCG of IZ, whereas a less evident effect was observed for IA. In general, both ceramics were prone to SCG, with crack propagation occurring at K(I) as low as 43-48% of their critical K(I). The highest sigma of IZ should lead to longer lifetimes for similar loading conditions.Significance. Water combined with cyclic loading causes pronounced SCG in IZ and IA materials. The lifetime of dental restorations based on these ceramics is expected to increase by reducing their direct exposure to wet conditions and/or by using high content zirconia ceramics with higher strength. (C) 2010 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Modelling of fatigue crack growth following overloads

In variable-amplitude loading there are interaction effects between the loading history and the crack propagation rate. The most important of these effects is the retardation in the crack propagation, which may raise the life of the cracked structureconsiderably. The main objective of this research is to analyse and quantify the retardation of crack propagation in a thin plate of the high-resistance aluminium alloy 2024-T3, comparing the results obtained from the mathematical models proposed to account for the retardation effect. The specimens were tested under high-low loading sequences, for different crack sizes and overload ratios. A simplified model was developed, based on crack closure theory, that could represent the crack behaviour during retardation very well. © 1991.

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.

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.

Fatigue crack growth rate in mode I of a carbon fiber 5HS weave composite laminate processed via RTM

Delamination or crack propagation between plies is a critical issue for structural composites. In viewing this issue and the large application of woven fabrics in structural applications, especially the ones that requires high drapeability to be preformed in a RTM mold cavity such as the asymmetric ones, e.g HS series, this research aimed in dynamically testing the carbon fiber 5HS/RTM6 epoxy composites under opening mode using DCB set up in order to investigate the crack growth rate behavior in an irregular surface produced by the fabric waviness. The evaluation of the energy involved in each crack increment was based on the Irwin-Kies equation using compliance beam theory. The tests were conducted at constant stress ratio of R=0.1 with displacement control, frequency of 10 Hz, in accordance to ASTM E647-00 for measurement of crack growth rate. The results showed large scatter when compared to unidirectional carbon fiber composites due to damage accumulation at the fill tows.