Effect of T-stress on crack growth under mixed mode I–III loading

For a crack subjected to combined mode I and III loading the influence of a T-stress is analyzed, with focus on crack growth. The solid is a ductile metal modelled as elastic–plastic, and the fracture process is represented in terms of a cohesive zone model. The analyzes are carried out for conditions of small scale yielding, with the elastic solution applied as boundary conditions on the outer edge of the region analyzed. For several combinations of the stress intensity factors K_I and K_III and the T-stress crack growth resistance curves are calculated numerically in order to determine the fracture toughness. In all situations it is found that a negative T-stress adds to the fracture toughness, whereas a positive T-stress has rather little effect. For given values of K_I and T the minimum fracture toughness corresponds to K_III = 0.     

Discrete dislocation simulations of the development of a continuum plastic zone ahead of a stationary Mode III crack

We present results from discrete dislocation simulations showing the development of the plastic zone in front of a Mode III crack under constant load. We find that the equilibrated zone is circular, in agreement with continuum mechanics predictions of the elastic-perfectly plastic Mode III crack. The size of the equilibrated zone scales as the square of the applied load (K_III), also in agreement with the continuum results. The zone approaches saturation exponentially, with a time that scales as K_III²/σ_p³, where σ_p is the Peierls stress. These results delineate conditions under which the classical, continuum predictions of elastic-plastic fracture mechanics are applicable.     

On the effects of characteristic lengths in bending and torsion on Mode III crack in couple stress elasticity

The problem of a stationary semi-infinite crack in an elastic solid with microstructures subject to remote classical K_III field is investigated in the present work. The material behavior is described by the indeterminate theory of couple stress elasticity developed by Koiter. This constitutive model includes the characteristic lengths in bending and torsion and thus it is able to account for the underlying microstructure of the material as well as for the strong size effects arising at small scales. The stress and displacement fields turn out to be strongly influenced by the ratio between the characteristic lengths. Moreover, the symmetric stress field turns out to be finite at the crack tip, whereas the skew-symmetric stress field displays a strong singularity. Ahead of the crack tip within a zone smaller than the characteristic length in torsion, the total shear stress and reduced tractions occur with the opposite sign with respect to the classical LEFM solution, due to the relative rotation of the microstructural particles currently at the crack tip. The asymptotic fields dominate within this zone, which however has limited physical relevance and becomes vanishing small for a characteristic length in torsion of zero. In this limiting case the full-field solution recovers the classical K_III field with square-root stress singularity. Outside the zone where the total shear stress is negative, the full-field solution exhibits a bounded maximum for the total shear stress ahead of the crack tip, whose magnitude can be adopted as a measure of the critical stress level for crack advancing. The corresponding fracture criterion defines a critical stress intensity factor, which increases with the characteristic length in torsion. Moreover, the occurrence of a sharp crack profile denotes that the crack becomes stiffer with respect to the classical elastic response, thus revealing that the presence of microstructures may shield the crack tip from fracture.     

Nonlinear interface shear fracture of end notched flexure specimens

The nonlinear analytical solutions of an end notched flexure adhesive joint or fracture test specimen with identical or dissimilar adherends are investigated. In the current study, a cohesive zone model (with arbitrary nonlinear cohesive laws) based analytical solution is obtained for the interface shear fracture of an end notched flexure (ENF) specimen with sufficiently long bond length. It is found that the scatter and inconsistency in calculating Mode II toughness may be significantly reduced by this model. The present work indicates that the Mode II toughness G_IIc under pure shear cracking condition is indeed very weakly dependent on the initial crack length. And this conclusion is well supported by the experimental results found in the literature. The parametric studies show that the interface shear strength is the most dominant parameter on the critical load. It is also interesting to note that with very short initial crack length and identical interface shear strength, higher Mode II toughness indeed cannot increase the critical load. Unlike the high insensitivity of critical load to the detailed shape of the cohesive law for Mode I peel fracture, the shape of the cohesive law becomes relatively important for the critical load of joints under pure Mode II fracture conditions, especially for joints with short initial crack length. The current study may help researchers deepen the understanding of interface shear fracture and clarify some previous concepts on this fracture mode.     

Effect of specimen size and crack depth on 3D crack-front constraint for SENB specimens

Three-dimensional (3D) elastic–plastic finite element analyses (FEA) are performed to study constraint effect on the crack-front stress fields for single-edge notched bend (SENB) specimens. Both rectangular and square cross-section of the specimens with a deep crack of a/W=0.5 are considered to investigate the effect of specimen size. A square-cross-section specimen with a shallow crack of a/W=0.15 is also considered to examine the effect of crack depth. Stresses from FEA at the crack front on different planes of the specimen are compared with those determined by the J–A₂ three-term solution. Results show that in-plane stress fields can be characterized by the three-term solution throughout the thickness even in the region near the free surface. Cleavage fracture toughness data is compared to predict the effects of specimen size and crack depth on fracture behavior. It is found that the distributions of crack opening stress are nearly the same for the SENB specimens at the critical J which is consistent with the RKR model. Furthermore our results indicate that there is a distinct relationship between the crack-front constraint and the cleavage fracture toughness. By introducing the failure curves, the minimum fracture toughness and scatter band can be well captured using the J–A₂ approach.     

On determination of mode II fracture toughness using semi-circular bend specimen

The cracked semi-circular specimen subjected to three-point bending has been recognized as an appropriate test specimen for conducting mode I, mode II and mixed mode I/II fracture tests in brittle materials. The manufacturing and pre-cracking of the specimen are simple. No complicated loading fixture is also required for a fracture test. However, almost all of the theoretical criteria available for mixed mode brittle fracture fail to predict the experimentally determined mode II fracture toughness obtained from the semi-circular bend (SCB) specimen. In this paper, a modified maximum tangential stress criterion is used for calculating mode II fracture toughness K_IIc in terms of mode I fracture toughness K_Ic. The modified criterion is used for predicting the reported values of mode II fracture toughness for two brittle materials: a rock material (Johnstone) and a brittle polymer (PMMA). It is shown that the modified criterion provides very good predictions for experimental results.     

A novel method in determination of dynamic fracture toughness under mixed mode I/II impact loading

The objective of this paper is to propose a novel methodology for determining dynamic fracture toughness (DFT) of materials under mixed mode I/II impact loading. Previous experimental investigations on mixed mode fracture have been largely limited to qusi-static conditions, due to difficulties in the generation of mixed mode dynamic loading and the precise control of mode mixity at crack tip, in absence of sophisticated experimental techniques. In this study, a hybrid experimental–numerical approach is employed to measure mixed mode DFT of 40Cr high strength steel, with the aid of the split Hopkinson tension bar (SHTB) apparatus and finite element analysis (FEA). A fixture device and a series of tensile specimens with an inclined center crack are designed for the tests to generate the components of mode I and mode II dynamic stress intensity factors (DSIF). Through the change of the crack inclination angle β (=90°, 60°, 45°, and 30°), the K_II/K_I ratio is successfully controlled in the range from 0 to 1.14. A mixed mode I/II dynamic fracture plane, which can also exhibit the information of crack inclination angle and loading rate at the same time, is obtained based on the experimental results. A safety zone is determined in this plane according to the characteristic line. Through observation of the fracture surfaces, different fracture mechanisms are found for pure mode I and mixed mode fractures.     

Role of inclusion stiffness and interfacial strength on dynamic matrix crack growth: An experimental study

Experimental simulations of dynamic crack growth past inclusions of two different elastic moduli, stiff (glass) and compliant (polyurethane) relative to the matrix (epoxy), are carried out in a 2D setting. Full-field surface deformations are mapped in the crack–inclusion vicinity optically. The crack growth behavior as a function of inclusion–matrix interfacial strength and the inclusion location relative to the crack is studied under stress-wave loading conditions. An ultra high-speed rotating mirror-type digital camera is used to record random speckle patterns in the crack–inclusion vicinity to quantify in-plane displacement fields. The crack-tip deformation histories from the time of impact until complete fracture are mapped and fracture parameters are extracted. The crack front is arrested by the symmetrically located compliant inclusion for about half the duration needed for complete fracture event. The dynamically propagating crack is attracted and trapped by the weakly bonded inclusion interface for both stiff and compliant symmetrically located inclusion cases, whereas it is deflected away by the strongly bonded stiff inclusion and attracted by strongly bonded compliant inclusion when located eccentrically. The crack is arrested by a strongly bonded compliant inclusion for a significant fraction of the total dynamic event and is longer than the one for the weakly bonded counterpart. The compliant inclusion cases show higher fracture toughness than the stiff inclusion cases. Measured crack-tip mode-mixities correlate well with the observed crack attraction and repulsion mechanisms. Macroscopic examination of fracture surfaces reveals much higher surface roughness and ruggedness after crack–inclusion interaction for compliant inclusion than the stiff one. Implications of these observations on the dynamic fracture behavior of micron size A-glass and polyamide (PA6) particle filled epoxy is demonstrated. Filled-epoxy with 3% V_f of PA6 filler is shown to produce the same dynamic fracture toughness enhancement as the one due to 10% V_f glass.     

Effect of pure mode I,II or III loading or mode mixity on crack growth in a homogeneous solid

Crack growth is analyzed numerically under combined mode I, II and III loading, or under loading in one of these modes alone. The solid is a ductile metal modelled as elastic–plastic, and the fracture process is represented in terms of a cohesive zone model. The analyses are carried out for conditions of small-scale yielding, with the elastic mixed mode solution applied as boundary conditions on the outer edge of the region analyzed. For pure mode I loading crack growth continued far beyond the maximum fracture toughness shows that the predicted subsequent steady-state toughness is well below the maximum. The reason for this is discussed in terms of the local stress and strain fields around the tip. For pure mode II or mode III loading it is shown that there is no maximum before the steady-state. Also results for different mixed mode conditions are presented and discussed in relation to the results for loading in only one mode. Most of the results are based on assuming that the peak tractions for tangential separation are equal to that for normal separation, but it is shown that a relatively smaller peak traction for tangential separation may significantly affect the predictions.     

Initiation and arrest of cracks on two pipeline steels at low temperature

Dynamic fracture toughness at initiationK _1d and fracture toughness at arrestK _1a were measured on two pipeline steel grades. Dynamic fracture toughness was measured at a very high loading rate with the help of split Hopkinson pressure bars. The values ofK _1d andK _1a are compared. The purpose of this work is to examine the possibilities of using dynamic fracture toughness at crack initiation as a lower bound of crack arrest toughness. This work has practical applications because crack arrest tests are difficult to perform, give scattered results and are costly and time consuming. This procedure shows that it is possible to economize and rationalize using intelligent technology.     

K Ic and J-resistance-curve measurements on Nevada tuff

Linear-elastic fracture mechanics and J-integral test methods were used to determine the static fracture behavior of a porous volcanic tuff. Notched and prefatigued specimens of two different sizes were tested in three-point bending. The fracture toughness (K _Ic ) and theJ-resistance (J vs. crack growth) curve for each specimen was determined. The results indicate that there is good agreement between the fracture parameters determined by the two methods; however, there is some dependence of the results on specimen size. Possible reasons for this specimen-size dependence are discussed.     

Critique on fracture criteria in mixed-mode-loading

The paper presents an improved formulation of mixed-mode fracture criteria. Both the fracture toughness K_IC and K_IIC are taken into account. Consequently, the crack propagation direction and the fracture envelope in the and plane depends on the toughness ratio . In the plane only one envelope for each criterion is obtained.     

Dynamic mechanical and fracture properties of an infiltrated TiC-1080 steel cermet

The dynamic mechanical and fracture properties of a TiC porous network infiltrated with1080 steel are reported. Following infiltration, the cermet is subjected to various heat treatments that affect essentially the steel matrix. Dynamic compression tests show that the heat treatments increase the fracture strength of the cermet. The quasi-static fracture toughness (K_Ic) is also increased by the heat treatments. The dynamic (initiation) fracture toughness (K_Id) is substantially higher (by about a factor of 3) than its static counterpart. Failure mechanisms consist mainly of cleavage of the TiC and matrix grains, along with minor interfacial decohesion. However, dynamic loading induces substantial damage around the crack tip, consisting essentially of cleavage of TiC grains. Microcrak toughnening is believed to be responsible for the high dynamic toughness of the material. The critical microstructural fracture event is thus identified as the spreading of TiC cleavage microcracks into the neighboring steel grains.     

Rate sensitivity of mixed mode interface toughness of dissimilar metallic materials: Studied at steady state

Crack propagation in metallic materials produces plastic dissipation when material in front for the crack tip enters the active plastic zone traveling with the tip, and later ends up being part of the residual plastic strain wake. Thus, the macroscopic work required to advance the crack is typically much larger than the work needed in the near tip fracture process. For rate sensitive materials, the amount of plastic dissipation typically depends on the rate at which the material is deformed. A dependency on the crack velocity should therefore be expected. The objective of this paper is to study the macroscopic toughness of crack advance along an interface joining two dissimilar rate dependent materials, characterized by an elastic-viscoplastic material model that approaches the response of a J₂-flow material in the rate independent limit. The emphasis here is on the rate sensitivity of the macroscopic fracture toughness under mixed Mode I/II loading. Moreover, special cases of joined similar rate dependent materials, as well as dissimilar materials where one substrate remains either elastic or approaches the rate independent limit is also included. The numerical analysis is carried out using the SSV model [Suo, Z., Shih, C., Varias, A., 1993. A theory for cleavage cracking in the presence of plastic flow. Acta Metall. Mater. 41, 1551–1557] embedded in a steady state finite element formulation, here assuming plane strain conditions and small-scale yielding. Results are presented for a wide range of material parameters, including noteworthy observations of a characteristic crack velocity at which the macroscopic toughness becomes independent of the material rate sensitivity. The potential of this phenomenon is elaborated on from a modeling point of view.     

A study of hysteresis in theK ID -å relation

A photoelastic study of high speed crack propagation in Homalite 100 was conducted to measure hysteresis in the constitutive relation forK _ID -å. The fracture specimen was designed to obtain both crack acceleration and deceleration during a single crack extension. Additional loads perpendicular to the cráck-propagation path were applied at secondary locations to accentuate the magnitude of acceleration-deceleration observed in this specimen. The photoelastic data were analyzed using dynamic stress field equations in conjunction with the least-squares over-deterministic method to obtain the crack-propagation fracture toughness,K _ID . Crack velocity, å, was determined numerically by differentiating a polynomial fitted to the crack length-time data in a leastsquares sense. Experimental results indicate that Homalite 100 does not exhibit significant hysteresis in theK _ID -å relation.     

In-Situ SEM J-Testing and Crack-Tip Micromechanics of Zircaloy-4 by Three-Point Bending

A three-point bend fixture has been designed, fabricated, and utilized to demonstrate the feasibility of performing in-situ J-testing at ambient and elevated temperatures inside a scanning electron microscope (SEM). Using the three-point bend test technique, in-situ SEM J-testing has been performed to measure the crack mouth opening displacement and crack extension as a function of the applied load in order to generate J-R curves for Zircaloy-4 at 25°C and 316°C. Once the J-R curve is determined, an equivalent K_J-resistance (K_J-R) curve is computed on the basis of a relationship between the J-integral (J) and the stress intensity factor (K). The J-R and K_J-R curves of Zircaloy-4 exhibit a rising R-curve behavior, while the elastic K-R curve underestimates the fracture resistance of Zircaloy-4 once substantial crack extension has occurred. For the specimen dimensions considered, the J-R curves generated by in-situ SEM J-tests are not sensitive to the specimen geometry and measure the actual fracture resistance of the material. Furthermore, the onset of crack extension is dictated by the emission of one or more slipbands from the crack tip, and a change in the crack-tip displacement field, followed by void formation along the slipband, and linkage of the voids with the main crack.     

Anomalies associated with energy release parameters for cracks in piezoelectric materials

For a central crack in a piezoelectric plate, the mode-I stress intensity factor (K_I), electric displacement intensity factor (K_D), energy release rates (G, G_M) and energy density factor (S) are obtained from the finite element results. For the impermeable crack, the numerical results of K_I and K_D are coupled; this error is contrary to the uncoupled analytical solutions. The error has little effect on the total energy release rate G and energy density factor S, but in some cases, large errors in the mechanical energy release rate G_M are observed. G is global while SED is local. Also G is negative which defies physics where energy cannot be created while crack attempts to extend as implied by G. Computations should be made for the J-integral and also show that J becomes negative. What this shows is that the global fracture energy criterion is not suitable to address the local release of energy because it includes the overall energy which are irrelevant to fracture initiation being a local behavior. In addition, the case study shows that the energy density theory is the better fracture criterion for the piezoelectric material. According to the results of S, it retards the crack growth when the external electric field and piezoelectric poling are on opposite directions. This conclusion agrees with analytical and experimental evidence in the past references.     

A finite element investigation of the effect of crack tip constraint on hole growth under mode i and mixed mode loading

In this work, the effect of constraint on hole growth near a notch tip in a ductile material under mode I and mixed mode loading (involving modes I and II) is investigated. To this end, a 2-D plane strain, modified boundary layer formulation is employed in which the mixed mode elastic K–T field is prescribed as remote boundary conditions. A finite element procedure that accounts for finite deformations and rotations is used along with an appropriate version of J₂ flow theory of plasticity with small elastic strains. Several analyses are carried out corresponding to different values of T-stress and remote elastic mode-mixity. The interaction between the notch and hole is studied by examining the distribution of hydrostatic stress and equivalent plastic strain in the ligament between the notch tip and the hole, as well as the growth of the hole. The implications of the above results on ductile fracture initiation due to micro-void coalescence are discussed.     

An experimental investigation of crack-path directional stability

Using the criterion that a crack will extend perpendicular to the maximum circumferential stress,σ _θ, we show that the directional stability of crack growth is governed by the location of microcrack initiation ahead of the crack tip. At distances greater than a geometrical radiusr _o, the maximum value ofσ _θ deviates from the position of symmetry. Thus, if we assume that the physical processes involved in fracture lead to crack initiation at a distancer _c ahead of the crack tip, the criterion for directional stability isr _o>r _c. Experimental and theoretical values ofr _o verify that, asr _o becomes small, the crack's directional stability deteriorates. Observing that a lengthwise compressive stress increasesr _o, a center-cracked specimen was developed which allows the application of controlled lengthwise compression independently of the opening-mode load. A detailed photoelastic analysis of the specimen has provided the value ofr _o as a function of the crack length. The value ofr _o is then compared with the expected microcrack initiation distances in ductile fracture. By applying sufficient lengthwise compression, we are able to make the crack grow straight and obtain numerous data points from this specimen which would otherwise be directionally unstable. The results indicate that, as the total lengthwise tensile stress at the crack tip increases, the fracture toughness also increases. Using this information we can then adjustK _Ic for zero lengthwise loading and obtain a geometry independent fracture toughness.     

Asymptotic fields in steady crack growth with linear strain-hardening

The asymptotic stress and velocity fields of a crack propagating steadily and quasi-statically into an elastic-plastic material are presented. The material is characterized by J₂-flow theory with linear strain- hardening. The possibility of reloading on the crack flanks is taken into account. The cases of anti-plane strain (mode III), plane strain (modes I and II), and plane stress (modes I and II) are considered. Numerical results are given for the strength of the singularity and for the distribution of the stress and velocity fields in the plastic loading, elastic unloading and plastic reloading regions, as functions of the strain-hardening parameter. An attempt is made to make a connection with the perfectly-plastic solutions in the limit of vanishing strain-hardening.