Plasticity induced fatigue crack growth retardation model for steel elements reinforced by composite patch

Fatigue tests on notched steel plates reinforced by composite patch showed that the application of carbon fiber reinforced polymers (CFRP) strips with pretension of the overlays prior to bonding. This resulted in a significant amount of additional fatigue life. In particular, the pre-tension produces a compressive field in the steel plate which reduces the stress ratio that enhances crack growth retardation. The fatigue crack propagation rate is postulated to be a function of the effective strain energy density factor range. Fatigue crack growth data showed that standard crack growth retardation model cannot be used to evaluate the minimum effective stress. Hence, an ad hoc plasticity model is introduced and validated using experimental results. The proposed technique is an extension of the well know Newman’s model. The bridging effect due to the reinforcing strips is analytically modeled in order to estimate the reduction of crack opening displacement and finally the magnification of the crack growth retardation. Numerical and experimental results match well and show a significant influence of the pre-tension level on the expected fatigue crack growth rate of a reinforced steel plate.     

A cyclic plasticity model for single crystals

The Armstrong–Frederick type kinematic hardening rule was invoked to capture the Bauschinger effect of the cyclic plastic deformation of a single crystal. The yield criterion and flow rule were built on individual slip systems. Material memory was introduced to describe strain range dependent cyclic hardening. The experimental results of copper single crystals were used to evaluate the cyclic plasticity model. It was found that the model was able to accurately describe the cyclic plastic deformation and properly reflect the dislocation substructure evolution. The well-known three distinctive regimes in the cyclic stress–strain curve of the copper single crystals oriented for single slip can be reproduced by using the model. The model can predict the enhanced hardening for crystals oriented for multislip, showing the model's ability to describe anisotropic cyclic plasticity. For a given loading history, the model was able to capture not only the saturated stress–strain response but also the detailed transient stress–strain evolution. The model was used to predict the cyclic plasticity under a high–low loading sequence. Both the stress–strain responses and the microstructural evolution can be appropriately described through the slip system activation.     

A continuum model for intermittent deformation of single crystal micropillars

Strain bursts are often observed during compression tests of single crystal micropillars. In this work, we formulate a new continuum model that accounts for the strain bursts within the framework of crystal plasticity. The strain bursts are separated from the loading stage (nearly elastic loading) by introducing a dimensionless constant in the continuum model, and are detected by load serrations. The boundary conditions in the context of micropillar compression are studied and they are shown to be changing and unpredictable as plastic deformation proceeds. To evaluate the validity of our model, finite element simulations of the uniaxial compression tests on nickel micropillars are performed. Our simulations produce clearly visible strain bursts during the plastic flow and the produced intermittent flows are comparable with the experimental observations. For the bulk crystal, a series of strain bursts is identified in the course of plastic flow, despite an apparently smooth stress–strain response. We also show that the intermittent flow is intensified in the micrometer-scale due to both increasing numbers of the successive strain bursts and increasing amplitude of the strain burst, when the specimen size decreases. Finally, we show that the occurrences of the strain bursts are always associated with negative values of the second-order work.     

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.     

Thermodynamic-based cohesive zone healing model for self-healing materials

This work presents a thermodynamic-based cohesive zone framework to model healing in materials that tend to self-heal. The nominal, healing and effective configurations of continuum damage-healing mechanics are extended to represent cohesive zone configurations. To incorporate healing in a cohesive zone model, the principle of virtual power is used to derive the local static/dynamic macroforce balance and the boundary traction as well as the damage and healing microforce balances. A thermodynamic framework for constitutive modeling of damage and healing mechanisms of cracks is used to derive the evolution equations for the damage and healing internal state variables. The effects of temperature, resting time, crack closure, history of healing and damage, and level of damage on the healing behavior of the cohesive zone are incorporated. The proposed model promises solid basis for understanding the self-healing phenomena in self-healing materials.     

Experimental studies on the effect of overload on fatigue crack growth

The fatigue crack growth behavior resulting from a single overload is investigated. In order to clarify the mechanism of overload on fatigue crack growth, the processes of crack closure and opening and their stress levels are monitored by strain gages placed on the back surface of specimens, and the fracture surface morphologies are examined by the microfractography. Experimental results may be used to explain quantitatively the mechanisms of retardation and delayed retardation after a single overload.     

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.     

Perturbation approaches of a planar crack in linear elastic fracture mechanics: A review

One current challenge of linear elastic fracture mechanics (LEFM) is to take into account the non-linearities induced by the crack front deformations. For this, a suitable approach is the crack front perturbation method initiated by Rice (1985). It allows to update the stress intensity factors (SIFs) when the crack front of a planar crack is perturbed in its plane. This approach and its later extensions to more complex cases are recalled in this review. Applications concerning the deformation of the crack front when it propagates quasistatically in a homogeneous or heterogeneous media have been considered in brittle fracture, fatigue or subcritical propagation. The crack shapes corresponding to uniform SIF have been derived: cracks with straight or circular fronts, but also when bifurcations exist, with wavy front. For an initial straight crack, it has been shown that, in homogeneous media, in the quasistatic case, perturbations of all lengthscales progressively disappear unless disordered fracture properties yields Family and Vicsek (1985) roughness of the crack front. Extension of those perturbation approaches to more realistic geometries and to coalescence of cracks is also envisaged.     

A model and numerical study for coiling of Kelvin-type viscoelastic filament

Coiling problems of elastic rope and viscous jet have been fully studied both experimentally and theoretically, but very few studies exist for viscoelastic material. In this paper, a system of one-dimensional two-point boundary nonlinear equations for Kelvin material coiling is presented. The equations are solved numerically by continuation method. It is found that the coiling frequency depends on the dimensionless retardation time and other continuation parameters involving inertia and gravity. The multivaluedness of the solution is observed in the numerical study.     

A generalized Paris’ law for fatigue crack growth

An extension of the celebrated Paris law for crack propagation is given to take into account some of the deviations from the power-law regime in a simple manner using the Wöhler SN curve of the material, suggesting a more general “unified law”. In particular, using recent proposals by the first author, the stress intensity factor K(a) is replaced with a suitable mean over a material/structural parameter length scale Δa, the “fracture quantum”. In practice, for a Griffith crack, this is seen to correspond to increasing the effective crack length of Δa, similarly to the Dugdale strip-yield models. However, instead of including explicitly information on cyclic plastic yield, short-crack behavior, crack closure, and all other detailed information needed to eventually explain the SN curve of the material, we include directly the SN curve constants as material property. The idea comes as a natural extension of the recent successful proposals by the first author to the static failure and to the infinite life envelopes. Here, we suggest a dependence of this fracture “quantum” on the applied stress range level such that the correct convergence towards the Wöhler-like regime is obtained. Hence, the final law includes both Wöhler's and Paris’ material constants, and can be seen as either a generalized Wöhler's SN curve law in the presence of a crack or a generalized Paris’ law for cracks of any size.     

An anisotropic micromechanics model for predicting the rafting direction in Ni-based single crystal superalloys

An anisotropic micromechanics model based on the equivalent inclusion method is developed to investigate the rafting direction of Ni-based single crystal superalloys. The micromechanical model considers actual cubic structure and orthogonal anisotropy properties. The von Mises stress, elastic strain energy density, and hydrostatic pressure in different inclusions of micromechanical model are calculated when applying a tensile or compressive loading along the [001] direction. The calculated results can successfully predict the rafting direction for alloys exhibiting a positive or a negative mismatch, which are in agreement with pervious experimental and theoretical studies. Moreover, the elastic constant differences and mismatch degree of the matrix and precipitate phases and their influences on the rafting direction are carefully discussed.     

Limitation on crack speed for a strip-craze model applied to PMMA

A strip-craze model is proposed to study crack propagation in polymers. A nonlinear differential equation is derived to govern the dynamic process of crack propagation. The viscous feature of the material in the craze zone is taken into account by means of an experimentally determined relationship between the craze stress and crack speed. By fitting experimental data of PMMA into the model, some parameters including the strip-craze length are deduced. A non-singular stress is introduced to control the crack propagation with a strip craze at its tip. Variations of the crack length and the crack speed with time are computed and their dependence on the non-singular stress is investigated. For PMMA, three stages of crack propagation are identified in terms of initial non-singular stress σ_ns0. When σ_ns0<60 MPa, the crack speed mm/s and the crack is basically stationary; when 60 <σ_ns0<95 MPa, then mm/s the crack is in slow propagation; when σ_ns0>95 MPa, then mm/s and the crack is in rapid propagation. The proposed model is applicable only in slow crack propagation.     

Geometrical disorder of the fronts of a tunnel-crack propagating in shear in some heterogeneous medium

A statistical analysis of the deformation of the fronts of a tensile tunnel-crack propagating in fatigue in some medium with spatially varying Paris constant was recently performed, with special emphasis on the evolution of the power spectra and correlation functions of the fluctuations of the fronts around reference straight lines. This study is extended here to coplanar propagation (along a weak plane) of a tunnel-crack loaded in mode 2+3. The results are rather similar to those previously obtained for mode 1. In particular, just like for tensile loadings, there is an effect of gradual selection in time of Fourier components of the fluctuations of the fronts of large wavelength. One novelty, however, is that for shear loadings, the fronts no longer tend to become symmetrical in time, so that correlations between crack front fluctuations at two points are higher for points located on the same front than for points located on distinct ones.     

A refined cohesive zone model that accounts for inertia of cohesive zone of a moving crack

Existing cohesive zone models assume that actual fracture zone of non-zero mass can be modeled by a line segment (cohesive zone) with no mass and inertia. In the present work, a simplified mass-spring model is presented to study inertia effect of cohesive zone on a mode-I steady-state moving crack. It is showed that fracture energy predicted by the present model increases dramatically when a finite limiting crack speed is approached. Reasonable agreement with known experiments indicates that the present model has the potential to catch the inertia effect of cohesive zone which has been ignored in existing cohesive zone models and better simulate dynamic fracture at high crack speed.     

A yield function for single crystals containing voids

A yield function for single crystals containing voids has been developed based on a variational approach. This first yield function is phenomenologically extended by modifying the dependence on the mean stress and introducing three adjustable parameters. Unit cell finite element calculations are performed for various stress triaxiality ratios, main loading directions and porosity levels in the case of a perfectly plastic FCC single crystal. The three model parameters are adjusted on the unit cell calculations so that a very good agreement between simulation results and the proposed model is obtained.     

Macroscopic and microscopic examination of the relationship between crack velocity and path and Rayleigh surface wave speed in single crystal silicon

The speed of Rayleigh surface waves, denoted C_R, is the accepted upper limit for Mode I crack velocity in monolithic solids. In the current contribution, we discuss several critical issues associated with the velocity of Rayleigh surface waves and crack velocity in single crystal (SC) brittle solids, and the global and local influence of C_R on crack path selection in particular.Recent cleavage experiments in SC silicon showed that crack velocity at certain cleavage planes and crystallographic orientations cannot exceed a small fraction of C_R, and thereafter the crack deflects to other cleavage planes. Indeed, C_R defined by the continuum mechanics ignores atomistic phenomena occurring during rapid crack propagation, and therefore is limited in predicting the crack velocity. Examination of these anomalies shows that this limitation lies in microstructural lattice arrangement and in anisotropic phonon radiation during rapid crack propagation. Globally, C_R has no influence on the crack deflection phenomenon. However, the misfit in C_R between the original plane of propagation and the deflected plane generates local instabilities along the deflection zone.     

A single segment differential tube model with interchain tube pressure effect

We develop a single segment differential tube model including interchain tube pressure effect (ITPE) [G. Marrucci, G. Ianniruberto, Interchain pressure effect in extensional flows of entangled polymers, Macromolecules 36 (2004) 3934–3942], able to describe the non-linear behaviour of entangled linear polymers. The model accounts for the effect of flow on the tube length and diameter. It is presented in two versions, depending on which tube dimension is assumed to deform affinely. The classical relaxation mechanisms, i.e., reptation, stretch dynamics, convective constraint release (CCR), as well as finite extensibility, are incorporated in a simple manner; hence the model allows an explicit comparison of the relative importance of various effects. A striking result is the insignificance of finite extensibility and the detrimental influence of CCR for moderately entangled systems when ITPE is taken into account. For highly entangled systems, CCR regains importance to avoid the well-known shear stress instability. The proposed model is able to make quantitative predictions of steady elongational and shear data for monodisperse melts, while transient values are less accurate but within experimental errors.     

A novel yield function for single crystals based on combined constraints optimization

A novel yield function representing the overall plastic deformation in a single crystal is developed using the concept of optimization. Based on the principle of maximum dissipation during a plastic deformation, the problem of single crystal plasticity is first considered as a constrained optimization problem in which constraints are yield functions for slip systems. To overcome the singularity that usually arises in solving the above problem, a mathematical technique is used to replace the above constrained optimization problem with an equivalent problem which has only one constraint. This single constraint optimization problem, the so-called combined constraints crystal plasticity (CCCP) model, is implemented into a finite element code and the results of modeling the uniaxial tensions of the single crystal copper along different crystallographic directions and also hydroforming of aluminum tubes proved the capability of the proposed CCCP model in accurately predicting the deformation in polycrystalline materials.     

A micro-mechanically based continuum damage model for fatigue life prediction of filled rubbers

This paper presents a continuum damage model based on two mechanisms: decohesion between fillers and matrix at a micro-scale followed by a crack nucleation at a macro-scale. That scenario was developed considering SEM observations and an original experimental procedure based on simple shear and tension specimens. Damage accumulation is related to fatigue life using the continuum damage mechanics (CDM). The material behavior is investigated using the statistical framework introduced by Martinez et al. (2011). A Finite Element implementation is proposed and some numerical examples are provided.