Dynamic crack tip fields and dynamic crack propagation characteristics of anisotropic material

Based on mechanics of anisotropic material, the dynamic crack propagation problem of I/II mixed mode crack in an infinite anisotropic body is investigated. Expressions of dynamic stress intensity factors for modes I and II crack are obtained. Components of dynamic stress and dynamic displacements around the crack tip are derived. The strain energy density theory is used to predict the dynamic crack extension angle. The critical strain energy density is determined by the strength parameters of anisotropic materials. The obtained dynamic crack tip fields are unified and applicable to the analysis of the crack tip fields of anisotropic material, orthotropic material and isotropic material under dynamic or static load. The obtained results show Crack propagation characteristics are represented by the mechanical properties of anisotropic material, i.e., crack propagation velocity M and fiber direction α. In particular, the fiber direction α and the crack propagation velocity M give greater influence on the variations of the stress fields and displacement fields. Fracture angle is found to depend not only on the crack propagation but also on the anisotropic character of the material.     

爆炸载荷下板条边界斜裂纹的动态扩展行为

为了研究爆炸应力波作用下板条边界斜裂纹的动态扩展行为,首先分析了爆炸应力波在含边界斜裂纹板条中的传播,其次采用动态焦散线实验方法,进行了爆炸载荷下板条边界斜裂纹扩展规律的实验研究.研究结果表明,爆炸应力波作用下,板条试件边界斜裂纹的扩展过程中,裂纹扩展速度、扩展加速度和裂尖动态应力强度因子随时间波动变化,扩展速度最大值...     

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

Multiscale modeling of crack initiation and propagation at the nanoscale

Fracture occurs on multiple interacting length scales; atoms separate on the atomic scale while plasticity develops on the microscale. A dynamic multiscale approach (CADD: coupled atomistics and discrete dislocations) is employed to investigate an edge-cracked specimen of single-crystal nickel, Ni, (brittle failure) and aluminum, Al, (ductile failure) subjected to mode-I loading. The dynamic model couples continuum finite elements to a fully atomistic region, with key advantages such as the ability to accommodate discrete dislocations in the continuum region and an algorithm for automatically detecting dislocations as they move from the atomistic region to the continuum region and then correctly “converting” the atomistic dislocations into discrete dislocations, or vice-versa. An ad hoc computational technique is also applied to dissipate localized waves formed during crack advance in the atomistic zone, whereby an embedded damping zone at the atomistic/continuum interface effectively eliminates the spurious reflection of high-frequency phonons, while allowing low-frequency phonons to pass into the continuum region.The simulations accurately capture the essential physics of the crack propagation in a Ni specimen at different temperatures, including the formation of nano-voids and the sudden acceleration of the crack tip to a velocity close to the material Rayleigh wave speed. The nanoscale brittle fracture happens through the crack growth in the form of nano-void nucleation, growth and coalescence ahead of the crack tip, and as such resembles fracture at the microscale. When the crack tip behaves in a ductile manner, the crack does not advance rapidly after the pre-opening process but is blunted by dislocation generation from its tip. The effect of temperature on crack speed is found to be perceptible in both ductile and brittle specimens.     

冲击载荷作用下圆孔缺陷对裂纹动态扩展行为的影响规律

空腔和裂纹缺陷通常共存于深部地下岩体中,它们共同影响着岩体的结构安全性与稳定性。为了探究动力扰动载荷下圆形空腔对裂隙岩体内裂纹扩展行为的影响规律,提出了不同圆孔倾角的直裂纹空腔圆弧开口试件（circular opening specimen with straight crack cavity, COSSCC）,利用自制大型落锤冲击实验装置进行动态加载实验,同时采用裂纹扩展计系统测试了裂纹的动态起裂时刻与裂纹扩展速度等各种断裂力学参数,随后采用有限差分软件Autodyn进行裂纹扩展路径与圆孔周围应力场的数值分析,并采用有限元软件Abaqus计算裂纹的动态起裂韧度与裂纹扩展过程中的动态扩展韧度。结果表明:（1）当圆孔倾角θ小于10°时,裂纹扩展路径会偏折并穿过圆孔表面;当圆孔倾角θ为20°与30°时,裂纹扩展路径向圆孔方向发生偏折但不会穿过圆孔,圆孔具有明显的裂纹扩展引导作用; 当圆孔倾角θ为40°与50°时,裂纹扩展路径不会发生偏折,圆孔引导作用明显减弱。（2）当裂纹扩展路径达到圆孔空腔附近时,裂纹尖端的拉伸应力区与圆孔边缘的拉伸应力区发生重合,此时裂纹扩展速度显著增大,裂纹动态断裂韧度显著减小。（3）裂纹的偏折方向与裂纹尖端最大周向应力的方向基本一致。（4）裂纹动态断裂韧度始终小于裂纹起裂韧度,且裂纹动态断裂韧度与裂纹动态扩展速度呈负相关关系。裂纹动态扩展速度越大,裂纹动态断裂韧度越小。     

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.     

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.     

Analysis of dynamic growth of a tensile crack in an elastic-plastic material

The influence of inertia on the stress and deformation fields near the tip of a crack growing in an elastic-plastic material is studied. The material is characterized by the von Mises yield criterion and J₂ flow theory of plasticity. The crack grows steadily under plane strain conditions in the tensile opening mode. Features of the stress and deformation state at points near the moving crack tip are described for elastic-perfectly plastic response and for several crack propagation speeds. It is found that inertia has a significant effect on the elastic-plastic response of material particles near the crack tip, and that elastic unloading may occur behind the crack tip for higher speeds. The relationship between the applied crack driving force, represented by a remote stress intensity factor, and the crack tip speed is examined on the basis of a critical crack tip opening angle growth criterion. The calculated result is compared with dynamic fracture toughness versus crack speed data for a 4340 steel.     

Effects of fractal crack

Experimental results indicate that propagation paths of cracks in geomaterials are often irregular, producing rough fracture surfaces which are fractal. In this paper, crack tip motion along a fractal crack trace is discussed. A fractal kinking model of the crack extension path is established to describe irregular crack growth. The length, velocity and kinking effects of the fractal crack are analysed. A formula is derived to describe the effects of fractal crack propagation on the dynamic stress intensity factor and on crack velocity. Finally, expressions of stress and displacement fields near the fractal crack tip are given.     

Quasicontinuum simulation of single crystal nano-plate with a mixed-mode crack

The quasicontinuum (QC) method is employed to simulate a nickel single crystal nano-plate with a mixed-mode crack. Atomic stresses near the crack tip are fitted according to the elastoplastic fracture mechanics equations. It is found that the atomic stress fields neighboring the crack tip are also singular and controlled by the atomic stress intensity factors. And then the critical energy release rates for brittle and ductile fracture are computed and compared in order to predict crack propagation or dislocation emission. Four possible slip directions at the crack tip are pointed out. Finally, the slip direction around the crack tip is determined by the shear stress and it is well consistent with the atomic pictures from the QC simulation.     

Stress-strain field near crack tip and calculation of critical stress of crack propagation in a pure bending beam of rectangular section with one-sided mode i crack

This paper studies the stress-strain field near crack tip in a pure bending beam of rectangular section with one-sided mode I crack by the analytic method of Ref. [1], then it gives the stress and strain components at the crack tip when the crack propagates and further it obtains the formulas of calculating the elastic deformed area width, the deformed intensity area width and the equation groups of calculating the critical stress of crack propagation, last the equation group of calculating critical stress of crack propagation is verified by calculating instance. The maximum error is 0.18%. First Received May 7, 1994.     

Brittle fracture dynamics with arbitrary paths. II. Dynamic crack branching under general antiplane loading

The dynamic propagation of a bifurcated crack under antiplane loading is considered. The dependence of the stress intensity factor just after branching is given as a function of the stress intensity factor just before branching, the branching angle and the instantaneous velocity of the crack tip. The jump in the dynamic energy release rate due to the branching process is also computed. Similar to the single crack case, a growth criterion for a branched crack is applied. It is based on the equality between the energy flux into each propagating tip and the surface energy which is added as a result of this propagation. It is shown that the minimum speed of the initial single crack which allows branching is equal to 0.39c, where c is the shear wave speed. At the branching threshold, the corresponding bifurcated cracks start their propagation at a vanishing speed with a branching angle of approximately 40°.     

Remarks on the energy release rate for an antiplane moving crack in couple stress elasticity

This paper is concerned with the steady-state propagation of an antiplane semi-infinite crack in couple stress elastic materials. A distributed loading applied at the crack faces and moving with the same velocity of the crack tip is considered, and the influence of the loading profile variations and microstructural effects on the dynamic energy release rate is investigated. The behavior of both energy release rate and maximum total shear stress when the crack tip speed approaches the critical speed (either that of the shear waves or that of the localized surface waves) is studied. The limit case corresponding to vanishing characteristic scale lengths is addressed both numerically and analytically by means of a comparison with classical elasticity results.     

Continuing crack-tip deformation and fracture for plane-strain crack growth in elastic-plastic solids

Analysis of the deformation field consistent with a Prandtl stress distribution travelling with an advancing plane-strain crack reveals the functional form of the near tip crack profile in an elastic-plastic solid. The crack opening δ is shown to have the form δ ～ r In (const./r) at a distance r from the tip. This observation coupled with data generated from finite element investigations of growing cracks in small-scale yielding permits the construction of a relation characterizing the deformation at an extending crack tip. A ductile crack-growth criterion consisting of the attainment of a critical opening at a small characteristic material distance from the tip is adopted. Predictions of the stability of a growing crack for both small-scale yielding specimens and those subject to general yielding are discussed.     

基于DIC的爆炸加载下脆性材料裂纹扩展规律的试验研究

通过数字图像相关法（DIC），应用PMMA对爆炸加载条件下脆性材料的裂纹扩展规律进行了试验研究。基于对称性试验模型，实现了裂纹尖端位置和应变场信息的同步记录。以此为基础，通过对比分析获知，主应变场应变值最大点不能作为裂纹尖端的判断依据。并以动态裂纹扩展速度为参量，应用断裂动力学和最小二乘牛顿迭代法，计算出了考虑惯性效应的Ⅰ-Ⅱ混合型裂纹的应力强度因子:K_Ⅰ和K_Ⅱ值会随着裂纹扩展方向改变而发生突变；K_Ⅰ最大值为2.63 MPa·m1/2，最小值为0.89 MPa·m1/2；其整体变化趋势表明，爆炸加载条件下脆性材料裂纹扩展随能量积聚和释放呈循环阶梯式递减发展。     

Molecular-dynamics simulation-based cohesive zone representation of intergranular fracture processes in aluminum

A traction-displacement relationship that may be embedded into a cohesive zone model for microscale problems of intergranular fracture is extracted from atomistic molecular-dynamics (MD) simulations. An MD model for crack propagation under steady-state conditions is developed to analyze intergranular fracture along a flat Σ99 [1 1 0] symmetric tilt grain boundary in aluminum. Under hydrostatic tensile load, the simulation reveals asymmetric crack propagation in the two opposite directions along the grain boundary. In one direction, the crack propagates in a brittle manner by cleavage with very little or no dislocation emission, and in the other direction, the propagation is ductile through the mechanism of deformation twinning. This behavior is consistent with the Rice criterion for cleavage vs. dislocation blunting transition at the crack tip. The preference for twinning to dislocation slip is in agreement with the predictions of the Tadmor and Hai criterion. A comparison with finite element calculations shows that while the stress field around the brittle crack tip follows the expected elastic solution for the given boundary conditions of the model, the stress field around the twinning crack tip has a strong plastic contribution. Through the definition of a Cohesive-Zone-Volume-Element—an atomistic analog to a continuum cohesive zone model element—the results from the MD simulation are recast to obtain an average continuum traction-displacement relationship to represent cohesive zone interaction along a characteristic length of the grain boundary interface for the cases of ductile and brittle decohesion.     

Dynamic stress intensity factor KⅢ and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials,the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated.Stress,strain and displacement around the crack tip are expressed as an analytical complex function,which can be represented in power series.Constant coefficients of series are determined by boundary conditions.Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained.Components of dynamic stress,dynamic strain and dynamic displacement around the crack tip are derived.Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials,i.e.,crack propagation velocity M and the parameter α.The faster the crack velocity is,the greater the maximums of stress components and dynamic displacement components around the crack tip are.In particular,the parameter α affects stress and dynamic displacement around the crack tip.     

On the dynamics of crack growth in glass

An experimental investigation of crack growth has been performed using plates of glass subjected to impulsive loading by means of small explosive charges. For comparison, some experiments were also performed with two “glass-like” polymers. The formation of the resultant crack patterns in glass was observed by the use of high-speed photography during the period following the detonation. Information has been obtained about the physical features of crack growth with regard to crack velocity, bifurcation, and the interaction between growing cracks and propagating stress waves. The existence of a terminal velocity of crack propagation is considered and, in particular, the degree to which crack bifurcation provides an explanation for this phenomenon is discussed.     

Equilibrium of a pressurized plastic fluid in a wellbore crack

This paper investigates equilibrium of a pressurized plastic fluid invading a tensile wellbore crack in a linear elastic, permeable rock. The crack is initially filled by pore fluid at ambient pressure, that is immiscibly displaced by the plastic fluid invading from the wellbore. The plastic fluid comes to rest to form a “plug” within the elastically deformed crack when the limit equilibrium between the shear stresses generated at the fracture walls and the pressure drop between the wellbore wall and the crack tip is reached. The model assumes that the leak-off of the plastic fluid into the rock is negligible, while the displaced pore fluid in the crack tip region is freely exchanged with the surrounding permeable rock to maintain the ambient pressure level. When the crack length ? is small or large compared to the wellbore radius R, the problem reduces to that of a pressurized edge or Griffith’s crack, respectively, subjected to a uniform far-field confining stress. In these two end-member cases, the normalized solution for the net pressure distribution, the plug length, and the stress intensity factor at the crack tip is obtained as a function of two numbers – the normalized net fluid pressure at the crack inlet and at the crack tip (partial plugs only) – that embody the solution’s dependence on the wellbore and the far field loading, the fluid yield strength, and the rock modulus. In the general case of an intermediate crack length (? ～ R), the normalized solution is a function of two additional parameters, the length-to-radius ratio and a normalized measure of the far field stress anisotropy, respectively, which accurate approximation is devised from an end-member solution using a rescaling argument. The equilibrium plug solutions are used to evaluate the breakdown pressure, the critical wellbore pressure at which the crack propagation condition is first met, and to analyze the stability of the ensuing crack propagation.