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.     

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.     

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.     

Constant moving crack in a magnetoelectroelastic material under anti-plane shear loading

Analytical solutions for an anti-plane Griffith moving crack inside an infinite magnetoelectroelastic medium under the conditions of permeable crack faces are formulated using integral transform method. The far-field anti-plane mechanical shear and in-plane electrical and magnetic loadings are applied to the magnetoelectroelastic material. Expressions for stresses, electric displacements and magnetic inductions in the vicinity of the crack tip are derived. Field intensity factors for magnetoelectroelastic material are obtained. The stresses, electric displacements and magnetic inductions at the crack tip show inverse square root singularities. The moving speed of the crack have influence on the dynamic electric displacement intensity factor (DEDIF) and the dynamic magnetic induction intensity factor (DMIIF), while the dynamic stress intensity factor (DSIF) does not depend on the velocity of the moving crack. When the crack is moving at very lower or very higher speeds, the crack will propagate along its original plane; while in the range of M_c1 < M < M_c2, the propagation of the crack possibly brings about the branch phenomena in magnetoelectroelastic media.     

Anisotropic plastic stress fields at a slowly propagating crack tip

Under the condition that any perfectly plastic stress components at a crack tip are nothing but the functions of 0 only making use of equilibrium equations. Hill anisotropic yield condition and unloading stress-strain relations, in this paper, we derive the general analytical expressions of anisotropic plastic stress fields at the slowly steady propagating tips of plane and anti-plane strain. Applying these general analytical expressions to the concrete cracks, the analytical expressions of anisotropic plastic stress fields at the-slowly steady propagating tips of Mode I and Mode III cracks are obtained. For the isotropic plastic material, the anisotropic plastic stress fields at a slowly propagating crack tip become the perfectly plastic stress fields.     

Stress intensity factor and crack velocity relationship for polyester/TiO2 nanocomposites

The dynamic fracture behavior of polyester/TiO₂ nanocomposites has been characterized and compared with that of the matrix material. A relationship between the dynamic stress intensity factor,K _I and the crack tip velocity,å, has been established. Dynamic photoelasticity coupled with high-speed photography has been used to obtain crack tip velocities and dynamic stress fields around the propagating cracks. Birefringent coatings were used to conduct the photoelastic study due to the opaqueness of the nanocomposites. Single-edge notch tension and modified compact tension specimens were used to obtain a broad range of crack velocities. Fractographic analysis was conducted to understand the fracture process. The results showed that crack arrest toughness in nanocomposites was 60% greater than in the matrix material. Crack propagation velocities prior to branching in nanocomposites were found to be 50% greater than those in polyester.     

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°.     

Brittle fracture dynamics with arbitrary paths I. Kinking of a dynamic crack in general antiplane loading

We present a new method for determining the elasto-dynamic stress fields associated with the propagation of anti-plane kinked or branched cracks. Our approach allows the exact calculation of the corresponding dynamic stress intensity factors. The latter are very important quantities in dynamic brittle fracture mechanics, since they determine the crack path and eventual branching instabilities. As a first illustration, we consider a semi-infinite anti-plane straight crack, initially propagating at a given time-dependent velocity, that changes instantaneously both its direction and its speed of propagation. We will give the explicit dependence of the stress intensity factor just after kinking as a function of the stress intensity factor just before kinking, the kinking angle and the instantaneous velocity of the crack tip.     

加载角度对层理页岩裂纹扩展影响的实验研究

采用分离式霍普金森压杆（SHPB）系统对页岩进行冲击实验,研究层理角度对页岩动态断裂过程的影响,在裂尖设置裂纹扩展计,借助高速摄影和数字图像相关（DIC）技术对页岩中心切槽半圆盘弯曲（NSCB）试件断裂的全过程进行研究,得到了不同加载角度下页岩的动态起裂韧度、裂纹扩展速度、断裂过程中应变场和水平位移场的变化规律。实验发现:不同加载角度下,页岩的动态起裂韧度具有显著的各向异性,加载角度与动态起裂韧度呈正相关;加载角度对试样的裂纹扩展速度具有显著影响,与裂纹扩展速度呈负相关;当冲击速度较低时,切槽方向是裂纹扩展的优势方向,而当冲击速度较高时,试样会产生沿层理弱面的次生裂纹,次生裂纹对试样的断裂具有显著影响。     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

Effect of orthotropy on crack propagation

The tendency of moving cracks to spread along the preferred directions of material anisotropy is treated. Depending on the velocity of crack propagation, the change of material properties in orthogonal planes is shown to affect the bifurcation characteristics. The problem is reduced to a system of dual integral equations that can be solved in a standard fashion. Of particular interest is the dynamic stress field near the tip of a moving crack in an orthotropic material. Although the 1√r stress singularity is preserved with r being the radial distance measured from the crack tip, the angular variations of the stresses are dependent on crack speed and material anisotropy. The possibility of crack bifurcation is examined by application of the strain energy density criterion for several composite systems. Crack branching is found to be enhanced by material anisotropy, a phenomenon that is not uncommon in composite materials.     

Rayleigh waves emitted by a propagating crack in a strain-rate dependent elastic medium

A simple model was proposed for the interpretation of the non-circular form of the Rayleigh wavefronts emitted by a fast running crack in a plate. The surface deformation around the crack tip, due to the high stress concentration there, propagated as a surface wave after fracture of this zone. On the other hand, the moving singularity of the crack tip created a dynamic stress field of varying intensity with time all over the specimen. This dynamic stress field resulted in a significant change of the mechanical properties of a strain-rate dependent material and therefore it influenced the velocity of propagation of fracture-Rayleigh wavefronts. An analysis of this varying dynamic strain field explained the non-circular form of Rayleigh waves, accompanying the propagating crack. For the experimental evaluation of the K₁-factor the method of dynamic caustics was used in conjunction with the high-speed photography technique.     

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.     

Asymptotic crack-tip fields for dynamic crack growth in compressive power-law hardening materials

The effect of material compressibility on the stress and strain fields for a mode-I crack propagating steadily in a power-law hardening material is investigated under plane strain conditions. The plastic deformation of materials is characterized by the J₂ flow theory within the framework of isotropic hardening and infinitesimal displacement gradient. The asymptotic solutions developed by the present authors [Zhu, X.K., Hwang K.C., 2002. Dynamic crack-tip field for tensile cracks propagating in power-law hardening materials. International Journal of Fracture 115, 323–342] for incompressible hardening materials are extended in this work to the compressible hardening materials. The results show that all stresses, strains, and particle velocities in the asymptotic fields are fully continuous and bounded without elastic unloading near the dynamic crack tip. The stress field contains two free parameters σ_eq0 and s₃₃₀ that cannot be determined in the asymptotic analysis, and can be determined from the full-field solutions. For the given values of σ_eq0 and s₃₃₀, all field quantities around the crack tip are determined through numerical integration, and then the effects of the hardening exponent n, the Poisson ratio ν, and the crack growth speed M on the asymptotic fields are studied. The approximate behaviors of the proposed solutions are discussed in the limit of ν → 0.5 or n → ∞.     

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

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

Mixed Mode Dynamic Fracture in Particulate Reinforced Functionally Graded Materials

A detailed analytical and experimental investigation is presented to understand the dynamic fracture behavior of functionally graded materials (FGMs) under mode I and mixed mode loading conditions. Crack-tip stress, strain and displacement fields for a mixed mode crack propagating at an angle from the direction of property gradation were obtained through an asymptotic analysis coupled with a displacement potential approach. This was followed by a comprehensive series of experiments to gain further insight into the behavior of propagating cracks in FGMs. Dynamic photoelasticity coupled with high-speed photography was used to obtain crack tip velocities and dynamic stress fields around the propagating cracks. Birefringent coatings were used to conduct the photoelastic study due to the opaqueness of the FGMs. Dynamic fracture experiments were performed using different specimen geometries to develop a dynamic constitutive fracture relationship between the mode I dynamic stress intensity factor (K _ID ) and crack-tip velocity ( ) for FGMs with the crack moving in the direction of increasing fracture toughness. A similar -K _ID relation was also obtained for matrix material (polyester) for comparison purposes. The results obtained show that crack propagation velocities in FGMs were about 80% higher than the polyester matrix. Crack arrest toughness was found to be about 10% lower than the value of local fracture toughness in FGMs.     

Asymptotic fields around an interfacial crack with a cohesive zone ahead of the crack tip

This paper considers an interfacial crack with a cohesive zone ahead of the crack tip in a linearly elastic isotropic bi-material and derives the mixed-mode asymptotic stress and displacement fields around the crack and cohesive zone under plane deformation conditions (plane stress or plane strain). The field solution is obtained using elliptic coordinates and complex functions and can be represented in terms of a complete set of complex eigenfunction terms. The imaginary portion of the eigenvalues is characterized by a bi-material mismatch parameter ε = arctanh(β)/π, where β is a Dundurs parameter, and the resulting fields do not contain stress singularity. The behaviors of “Mode I” type and “Mode II” type fields based on dominant eigenfunction terms are discussed in detail. For completeness, the counterpart for the Mode III solution is included in an appendix.     

Evaluation of crack initiation angle under mixed mode loading at diverse strain rates

Crack initiation angle, under mixed mode loading at several strain rates, is analysed using an experimental–numerical approach. The physical phenomenon for the problem at hand is influenced by the local and global conditions. One of such factors is the strain rate at the crack tip. For this purpose, PMMA plates with centred angled cracks under mixed mode loading were tested. The strain rate at the neighbourhood of the crack tip before crack propagation was evaluated. Considering that this material is strain rate sensitive, the numerical models were calibrated with the modulus of elasticity measured in tension tests at the observed strain rates. Numerical evaluations were performed with the finite element method in conjunction with the volume energy density criterion. An improvement in the evaluation of the crack propagation angle was observed. In order to complete the analysis, the crack initiation angle was also evaluated with the strain energy density factor S, considering the mechanical properties of PMMA, as evaluated at the observed strain rates, and the stress intensity factors k₁ and k₂. Results are in agreement with those observed experimentally.     

Inclined semi-elliptical crack for predicting crack growth direction based on apparent stress intensity factors

Linear elastic criterion of the inclined semi-elliptical crack growth direction is elaborated on the basis of the strain energy density theory. Stress and displacement fields are presented for higher order terms asymptotic expansion. Solutions for elastic stress intensity factors are accounting for the function describing of the crack tip fields near the free surface of plate. The mixed mode behavior of crack growth direction angle along the semi-elliptical crack front for different combination of biaxial loading, inclination crack angle and surface flaw geometry is determined.     

Effects of characteristic material lengths on mode III crack propagation in couple stress elastic–plastic materials

The asymptotic fields near the tip of a crack steadily propagating in a ductile material under Mode III loading conditions are investigated by adopting an incremental version of the indeterminate theory of couple stress plasticity displaying linear and isotropic strain hardening. The adopted constitutive model is able to account for the microstructure of the material by incorporating two distinct material characteristic lengths. It can also capture the strong size effects arising at small scales, which results from the underlying microstructures. According to the asymptotic crack tip fields for a stationary crack provided by the indeterminate theory of couple stress elasticity, the effects of microstructure mainly consist in a switch in the sign of tractions and displacement and in a substantial increase in the singularity of tractions ahead of the crack-tip, with respect to the classical solution of LEFM and EPFM. The increase in the stress singularity also occurs for small values of the strain hardening coefficient and is essentially due to the skew-symmetric stress field, since the symmetric stress field turns out to be non-singular. Moreover, the obtained results show that the ratio η introduced by Koiter has a limited effect on the strength of the stress singularity. However, it displays a strong influence on the angular distribution of the asymptotic crack tip fields.