Mixed mode (I and II) crack tip fields in bulk metallic glasses

Stationary crack tip fields in bulk metallic glasses under mixed mode (I and II) loading are studied through detailed finite element simulations assuming plane strain, small scale yielding conditions. The influence of internal friction or pressure sensitivity on the plastic zones, notch deformation, stress and plastic strain fields is examined for different mode mixities. Under mixed mode loading, the notch deforms into a shape such that one part of its surface sharpens while the other part blunts. Increase in mode II component of loading dramatically enhances the normalized plastic zone size, lowers the stresses but significantly elevates the plastic strain levels near the notch tip. Higher internal friction reduces the peak tangential stress but increases the plastic strain and stretching near the blunted part of the notch. The simulated shear bands are straight and extend over a long distance ahead of the notch tip under mode II dominant loading. The possible variations of fracture toughness with mode mixity corresponding to failure by brittle micro-cracking and ductile shear banding are predicted employing two simple fracture criteria. The salient results from finite element simulations are validated by comparison with those from mixed mode (I and II) fracture experiments on a Zr-based bulk metallic glass.     

Plasticity aspects of interfacial fracture for polymer/glass specimens

Experimental results suggest that the interfacial fracture resistance is minimal for approximate near tip Mode I accompanied by positive and negative near tip Mode II. Finite-strain FE analysis is made for an elastic–plastic medium bonded to an ideally elastic medium with an interface crack. Small-scale plasticity conditions are invoked and examined in relation to the elastic–plastic stress distribution along the bond line. Plasticity engenders a tendency to turn near tip biaxiality towards pure Mode I regardless of the mixed-mode loading. High levels of hydrostatic stress are attained. For different mode mixities of the applied load, the dependence of the elastic–plastic normal bond stress on load level is examined. It is found that under positive Mode II loading, the normal bond stress σ_yy tends to saturate as the load level rises. This does not occur for Mode I and negative Mode II loading. In addition, deformation patterns inside the plastic zone are examined for mixed-mode situations. A displacement criterion based on the normal bond crack opening suggests a dependence of the critical load level on the extent of mixed mode. Under positive mode II fracture, traces of the ductile material are found at the top of the elastic substrate. Some of these conclusions appear to be consistent with the fracture patterns observed for LD-polyethylene/glass interfacial mixed-mode fracture.     

一种基于SHTB的Ⅱ型动态断裂实验技术

冲击剪切载荷作用下动态断裂韧性的测定是材料力学性能和断裂行为研究中重要组成部分.为了测定材料的Ⅱ型动态断裂韧性,许多学者采用不同的试样与实验方法进行了实验,但限于实验条件,裂纹断裂模式往往是I+Ⅱ复合型,而不是纯Ⅱ型,因而不能准确测得材料的Ⅱ型动态断裂韧性.鉴于此,本文基于分离式霍普金森拉杆(split Hopkinson tension bar,SHTB)实验技术,提出一种改进的紧凑拉伸剪切(modified compact tension shear,MCTS)试样,通过夹具对MCTS试样施加约束,从而保证试样按照纯Ⅱ型模式断裂.采用实验-数值方法对MCTS试样动态加载过程进行分析,将实验测得的波形输入有限元软件ANSYS-LSDYNA,得到了裂纹尖端应力强度因子-时间曲线,并与紧凑拉伸剪切(compact tension shear,CTS)试样进行了对比.同时采用数字图像相关法进行了实验,验证了有限元分析结果.结果表明,MCTS试样在整个加载过程中K_I K_Ⅱ,裂纹没有张开;而CTS试样在同样的加载过程中K_IK_Ⅱ,出现裂纹张开现象.这说明MCTS试样能够准确地测定材料的Ⅱ型动态断裂韧性,为材料动态力学测试提供了一种有效的实验技术.     

Effect of temperature on the fracture modes of BHW-35 steel under mode-II loading

Finite element method (FEM) has been used to analyze the stress and strain fields and the stress tri-axial levels around the tip of the crack under mode- II loading. The results show that: under mode- II loading, the direction of the maximum tensile stress and that of the maximum tri-axial levels (R _o ) exist at an angle of –75. 3° from the original crack plane; the maximum shear stress andR _o = 0 exist along the original crack plane.Mode- II loading experiment using BHW-35 steel at different temperatures show that there are two kinds of fracture mode, opening mode (or tensile mode) and sliding mode (or shear mode). A decrease in temperature causes the fracture mode to change from shear mode to tensile mode. For BHW-35 steel, this critical temperature is about –90 C. Actually, under any kind of loading mode (mode I . mode II , mode III or mixed mode), there always exist several kinds of potenital fracture modes (for example, opening mode, sliding mode, tearing mode or mixed mode). The effect of temperature under mode- II loading is actually related to the change of the elastic-plastic properties of the material.     

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.     

K III -Deformation modes in internal oblique cracks under plane-stress conditions

The exact shape of the deformed internal oblique crack in an infinite elastic plate under conditions of plane stress was studied. The Muskhelishvili potential function yielded the exact stress and displacement field around the crack. While in previous papers by the author the study was mainly concerned with the definition of the in-plane shape of the deformed crack and important properties were disclosed concerning especially the contribution of the shear loading of the crack, in this paper the out-of-plane component of displacements is determined and its influence on the exact shape of the deformed crack in space is presented. It is shown that, as soon as shear displacements appear in the cracked plate under plane stress, out-of-plane shear displacements are a compulsory consequence for plane-stress conditions of the plate. The elliptic form of the deformed internal crack was twisted out of the plane of the plate with its zero twisting displacement near the new crack tip of its deformed shape corresponding to the vertex of the ellipse. The points of maximum and minimum out-of-plane displacements were placed close to the vertices of the ellipse at polar angles, θ, depending only on the eccentricity of the ellipse and displaced always on both sides of the vertices. The compulsory coexistence of mode II and mode III deformations makes the internal crack in plane stress to present a complicated pattern of deformation at its deformed crack tips. All of these results are amply supported by experimental evidence with caustics, which always show either a simple mode I pattern or a complex mode II and III pattern as soon as shear interferes in the mode of deformation of the plate.     

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.     

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅰ)—THE ANALYSES OF CONSTRAINT PARAMETERS AND FRACTURE PARAMETERS

In the present paper, three dimensional analyses of some general constraint parameters and fracture parameters near the crack tip. of Mode I CT specimens in two different thicknesses are carried out by employing ADINA program. The results reveal that the constraints along the thickness direction are obviously separated into two parts: the keeping similar high constraint field (Z₁) and rapid reducing constraints one (Z₂). The two fields are experimentally confiremed to correspond to the smooth region and the shear lip on the fracture face respectively. So the three dimensional stress structure of Mode I specimens can be derived through discussing the two fields respectively. The distribution of the Crack Tip Opening Displacement (CTOD) along the thickness direction and the three dimensional distribution of the void growth ratio (V_g) near the crack tip are also obtained. The two fracture parameters are in similar trends along the thickness direction, and both of them can reflect the effect of thickness and that of the loading level to a certain degree.     

Fracture Properties of Concrete–Concrete Interfaces Using Digital Image Correlation

The mode I and mode II fracture toughness and the critical strain energy release rate for different concrete–concrete jointed interfaces are experimentally determined using the Digital Image Correlation technique. Concrete beams having different compressive strength materials on either side of a centrally placed vertical interface are prepared and tested under three-point bending in a closed loop servo-controlled testing machine under crack mouth opening displacement control. Digital images are captured before loading (undeformed state) and at different instances of loading. These images are analyzed using correlation techniques to compute the surface displacements, strain components, crack opening and sliding displacements, load-point displacement, crack length and crack tip location. It is seen that the CMOD and vertical load-point displacement computed using DIC analysis matches well with those measured experimentally.     

A finite element study of the near crack tip deformation of a ductile material under mixed mode loading

In ductile fracture, voids near a crack tip play an important role. From this point of view, a large deformation finite element analysis has been made to study the deformation, stress and strain, and void ratio near the crack tip under mixed mode plane strain loading conditions, employing Gurson's constitutive equation which has taken into account the effects of void nucleation and growth. The results show that: (i) one corner of the crack tip sharpens while the other corner blunts, (ii) the stress and strain distributions except for the near crack tip region, can be superimposed by normalizing distance from the crack tip by a crack tip deformation length, i.e., a steady-state solution under a mixed mode condition has been obtained, (iii) the field near a crack tip can be divided into four characteristic fields (K field, HRR field, blunted crack tip field, and damaged region), and (iv) the strain and void volume fraction become concentrated in the sharpened part of a crack tip with increasing Mode II component.     

Experimental and Theoretical Assessment of Brittle Fracture in Engineering Components Containing a Sharp V-Notch

A criterion was proposed to predict brittle fracture in engineering components containing sharp V-shaped notches and subjected to mixed mode I/II loading. The criterion, called SV-MTS, was developed based on the maximum tangential stress (MTS) criterion proposed originally for analyzing crack problems. The curves which are obtained from the SV-MTS criterion could be used conveniently to predict the fracture resistance and also the notch bifurcation angle in sharp V-notched components under pure mode II and also mixed mode loading. To evaluate the validity of the proposed criterion, a set of fracture tests were conducted on a new test specimen, called sharp V-notched Brazilian disc (SV-BD), under mixed mode loading conditions. It is shown that the experimental results obtained from PMMA specimens are in very good agreement with the curves of SV-MTS criterion.     

LC4—M材料复合型韧断主要影响因素的分析

通过分析LC4－M铝合金材料在不同复合比载荷下的断裂试验结果,参照常规断裂现象,修正了一般断裂试验中认定裂纹启裂方向的方法,结合不同复合比下裂尖附近应力三维、主应力方向的计算分析,得到：在裂端的钝化变形区域、应力三维度的极大值处,对应于裂纹的启裂位置,即使在高韧性材料中发生剪切断裂的情况下也是如此;裂纹的启裂方向在拉伸断裂时与启裂点最大拉应力方向有关,在剪切断裂时启裂点最大剪应力方向有关,引起两种形式断裂的主要因数和破坏机理有很大不同。     

Numerical simulation of fracture mode transition in ductile plates

Fracture mode of ductile solids can vary depending on the history of stress state the material experienced. For example, ductile plates under remote in-plane loading are often found to rupture in mode I or mixed mode I/III. The distinct crack patterns are observed in many different metals and alloys, but until now the underlying physical principles, though highly debated, remain unresolved. Here we show that the existing theories are not capable of capturing the mixed mode I/III due to a missing ingredient in the constitutive equations. We introduce an azimuthal dependent fracture envelope and illustrate that two competing fracture mechanisms, governed by the pressure and the Lode angle of the stress tensor, respectively, exist ahead of the crack tip. Using the continuum damage plasticity model, we demonstrate that the distinctive features of the two crack propagation modes in ductile plates can be reproduced using three dimensional finite element simulations. The magnitude of the tunneling effect and the apparent crack growth resistance are calculated and agree with experimental observations. The finite element mesh size dependences of the fracture mode and the apparent crack growth resistance are also investigated.     

Mode II fracture behavior of a Zr-based bulk metallic glass

Shear band formation and fracture are characterized during mode II loading of a Zr-based bulk metallic glass. The measured mode II fracture toughness, K_IIc=75±4 MPa√m, exceeds the reported mode I fracture toughness by ∼4 times, suggesting that normal or mean stresses play a significant role in the deformation process at the crack tip. This effect is explained in light of a mean stress modified free volume model for shear localization in metallic glasses. Thermal imaging of deformation at the mode II crack tip further reveals that shear bands initiate, arrest, and reactivate along the same path, indicating that flow in the shear band leads to permanent changes in the glass structure that retain a memory of the shear band path. The measured temperature increase within the shear band is a fraction of a degree. However, heat dissipation models indicate that the temperature could have exceeded the glass transition temperature for less than 1 ms immediately after the shear band formed. It is shown that this time scale is sufficient for mechanical relaxation slightly above the glass transition temperature.     

An exact transient analysis of plane wave diffraction by a crack in an orthotropic or transversely isotropic solid

A transient plane strain analysis of diffraction of plane waves by a semi-infinite crack in an unbounded orthotropic or transversely isotropic solid is performed. The waves approach the crack at a general oblique angle, and are of two types, a normal stress pulse and a shear stress pulse, i.e. a P- and an SV-wave, respectively, in the isotropic limit. A class of materials that includes this limit and beryl, cobalt, ice, magnesium and titanium is chosen for illustration, and exact solutions are obtained for the initial/mixed boundary value problems.In contrast to related work, a factorization in the Laplace transform space is used to simplify the solution forms and the Wiener-Hopf component of the solution process, and to yield a more compact expression for the Rayleigh wave speed. Calculations for this speed, the two allowable, direction-dependent, plane wave speeds, and quantities related to the Mode I and Mode II dynamic stress intensity factors are given for the five anisotropic materials mentioned.     

Ti-6Al-4V在高应变率下的动态剪切特性及失效机理

采用基于霍普金森压杆的新型加载技术对Ti-6Al-4V材料的动态剪切特性及失效机理进行了测试研究。获得了Ti-6Al-4V材料在超过10⁴ s^-1应变率下的剪应力-剪应变曲线及失效参数。研究发现,材料的流动应力存在明显的应变率强化效应;随着应变率的增加,材料的失效应力逐渐增大,而失效应变逐渐减小。采用ABAQUS/Explicit对加载过程进行了数值模拟。结果显示,剪切区材料基本处于平面剪切状态,应力应变场分布较为均匀,计算得到的剪应力-剪应变曲线与实验结果吻合较好。经断口分析可知,随着应变率的升高,Ti-6Al-4V的失效机理存在由韧窝、拉伸韧窝至台阶及河流花样的演化过程,材料的失效模式主要表现为韧性断裂。

     

THE ANALYSES OF THE THREE-DIMENSIONAL STRESS STRUCTURE NEAR THE CRACK TIP OF MODE I CT SPECIMENS IN ELASTICPLASTIC STATE(Ⅰ)—THE ANALYSES OF CONSTRAINT PARAMETERS AND FRACTURE PARAMETERS

In the present paper,three dimensional analyses of some general constraintparameters and fracture parameters near the crack tip of Mode I CT specimens in twodifferent thicknesses are carried out by employing ADINA program.The results revealthat the constraints along the thickness direction are obviously separated into twoparts:the keeping similar high constraint field(Z_1)and rapid reducing constraintsone(Z_2).The two fields are experimentally confiremed to correspond to the smoothregion and the shear lip on the fracture face respectively.So the three dimensionalstress structure of Mode I specimens can be derived through discussing the two fieldsrespectively.The distribution of the Crack Tip Opening Displacement(CTOD)alongthe thickness direction and the three dimensional distribution of the void growth ratio(V_9 )near the crack tip are also obtained.The two fracture parameters are in similartrends along the thickness direction,and both of them can reflect the effect ofthickness and that of the loading l     

Failure modes of specimens containing surface flaws under cyclic torsion

Turbine-generator shafts are often subjected to a complex transient torsional loading. Such a torsional loading may initiate yielding at the outer radius of the shaft or in the fillets. The methods for predicting turbine-shaft fatigue life due to transient loading depent upon the mode of crack growth from an undetected crack. The most common location for the existence of a crack is the fillets or shoulders of the shaft. Specimens were designed from AISI 4340 steel with two diametrically opposed flat surfaces. Initial defect orientations of 0 deg, 45 deg and 90 deg with respect to the sepcimen axis on the fillet were studied. The specimens were subjected to cyclic torsion with zero mean torque and with a torque amplitude necessary to cause yielding at the outer radius of the specimen. When initial defects were aligned with a plane of maximum shear stress (0 deg and 90 deg), the cracks propagated along that plane. For 45-deg defects (aligned to a plane of maximum tensile stress) the crack still propagated along the plane of maximum shear. However, the number of cycles to initiate and to propagate the crack to failure for 45-deg defects were (two to three times) larger than those for 0-deg and 90-deg defects. Mode II and Mode III crack-growth rates were measured from specimens containing 0-deg and 90-deg defects. It was found that the crack-growth rate in Mode II was higher than in Mode III. However, all the specimens failed due to reduction of the net cross section, mostly attributed to Mode III crack growth. Similar results were obtained from specimens of turbine-shaft material (A469 steel), and 2024 aluminum with different rolling directions. Fatigue-crack-growth rates in Mode III were measured from circumferentially notched bar. They were found to be a unique function of ΔK _III alternating stress intensity in Mode III. It was found that the mechanism of crack growth is produced by the formation and linkage of elongated cavities at the crack tip.     

T应力对光弹性条纹影响的理论分析

以裂尖弹性应力场的多参数模型为基础,研究I型、II型以及I/II混合型裂纹参数对光弹性条纹的影响. T应力的存在和符号影响着等色线条纹环的半径大小和旋转方向,对于纯I型或II型裂纹而言,条纹环的旋转角度只与T应力有关;而对于I/II混合型裂纹,条纹环旋转角度与K_{\rm I}, K_{\rm II}和T应力有关. T应力的存在使得I型裂纹在裂尖±π/3方向上出现2个各向同性点(T应力为正时),使得II型裂纹在裂尖后的裂纹面上出现1个各向同性点. 对于I/II混合型裂纹而言,当T应力为正时等倾线出现距裂尖半径不等的3个各向同性点;反之, T应力为负时在裂尖后只存在1个各向同性点,这些各向同性点分别与I型和II型裂纹情况具有相同的规律.