Theoretical analysis of crack front instability in mode I+III

This paper focusses on the theoretical prediction of the widely observed crack front instability in mode I+III, that causes both the crack surface and crack front to deviate from planar and straight shapes, respectively. This problem is addressed within the classical framework of fracture mechanics, where the crack front evolution is governed by conditions of constant energy-release-rate (Griffith criterion) and vanishing stress intensity factor of mode II (principle of local symmetry) along the front. The formulation of the linear stability problem for the evolution of small perturbations of the crack front exploits previous results of Movchan et al. (1998) (suitably extended) and Gao and Rice (1986), which are used to derive expressions for the variations of the stress intensity factors along the front resulting from both in-plane and out-of-plane perturbations. We find exact eigenmode solutions to this problem, which correspond to perturbations of the crack front that are shaped as elliptic helices with their axis coinciding with the unperturbed straight front and an amplitude exponentially growing or decaying along the propagation direction. Exponential growth corresponding to unstable propagation occurs when the ratio of the unperturbed mode III to mode I stress intensity factors exceeds some “threshold” depending on Poisson's ratio. Moreover, the growth rate of helical perturbations is inversely proportional to their wavelength along the front. This growth rate therefore diverges when this wavelength goes to zero, which emphasizes the need for some “regularization” of crack propagation laws at very short scales. This divergence also reveals an interesting similarity between crack front instability in mode I+III and well-known growth front instabilities of interfaces governed by a Laplacian or diffusion field.     

压电压磁材料中正n边形孔边裂纹分析

基于电磁复合材料力学,运用Stroh型公式和复变函数方法,针对压电压磁材料中含正n边形孔边裂纹反平面问题进行了研究。利用Schwarz-Christoffel变换技术,结合Cauchy积分公式和留数定理,导出了磁电全非渗透型边界条件下任意正n边形裂纹尖端场强度因子和能量释放率的解析解。当缺失磁场时,所得解退化为已有结果,以此验证方法的有效性。通过数值算例,对比分析了n=3,n=4和n=5三种特殊情形对应的孔口边长、裂纹长度和受到的力、电和磁载荷对等效场强度因子和无量纲能量释放率的影响规律。研究结果发现,正n边形孔洞的尺寸和裂纹长度均会促进裂纹扩展,且前者的影响更显著一些;正n边形边的数量增加会阻止裂纹的扩展;在磁电全非渗透型边界条件下,机械载荷始终促进裂纹的扩展,电位移载荷可以促进或抑制裂纹的扩展,磁载荷对裂纹的扩展贡献较少。本研究结果适用于任意正n边形孔边裂纹求解问题,为压电压磁材料元器件的优化设计和断裂特性分析提供了新思路。     

Elliptical crack normal to functionally graded interface of bonded solids

This paper presents an analysis of an elliptical crack that is perpendicular to a functionally graded interfacial zone between two fully bonded solids. The functionally graded interfacial zone is treated as a non-homogeneous solid layer with its elastic modulus varying in the thickness direction. A generalized Kelvin solution based boundary element method is employed for the calculation of the stress intensity factors associated with the three-dimensional crack problem. The elliptical crack surface is subject to either uniform normal traction or uniform shear traction. The stress intensity factors are examined by taking into account the effects of the non-homogeneity parameter and thickness of the functionally graded interfacial zone, as well as the crack distance to the zone. The SIF values are further incorporated into the S-criterion for prediction of crack growth. The paper presents the most possible direction and location of the elliptical crack growth under an inclined tensile (or compressive) load. The paper further presents results of the critical external loads that would cause the elliptical crack to grow at the most possible location and along the most possible direction. The paper also examines the effects of external load direction and material and geometrical parameters on the critical loads.     

基于裂纹闭合效应的高载迟滞预测模型修正

裂纹闭合效应通常是导致I 型裂纹扩展在高载作用下发生迟滞效应的主要因素之一．本文采用汽车薄板QSTE340TM 材料，针对不同应力比，高载比条件下疲劳裂纹扩展行为进行了实验研究．论文通过断面分析，针对各参数对裂纹闭合效应的具体影响进行了分析讨论，认为裂纹作用区域随裂纹扩展而动态变化，从而提出了一种对有效应力强度因子幅的修正方法．通过在原有模型中引入幂函数形式的动态变量α，表征裂纹闭合效应的作用比例随裂纹长度的动态变化，取得了较好的效果．     

一种经验型疲劳小裂纹扩展模型

金属材料疲劳寿命由裂纹萌生和裂纹扩展寿命两部分组成,其中对于萌生寿命中的小裂纹分析是精确描述裂纹萌生寿命的关键．而小裂纹在扩展过程中由于尺寸相对较小,导致传统线弹性断裂力学预测方法失效,需要对其进行改进,考虑裂纹尖端塑性区引起的残余压应力对小裂纹扩展速度的影响．本文针对此问题进行了初步分析,通过对塑性区引起的残余应力的量化,结合小裂纹门槛值特性,提出了一种经验型修正的小裂纹扩展模型,用于定量预测裂纹的萌生寿命．使用铝合金6082-T6缺口试样进行了疲劳实验,并与理论结果进行了对比,验证了所提模型的有效性．     

Fatigue growth modeling of cracks emanating from a circular hole in infinite plate

This paper presents a numerical approach of fatigue growth analysis of cracks emanating from a hole in infinite elastic plate subjected to remote loads. It involves a generation of Bueckner’s principle and a hybrid displacement discontinuity method (a boundary element method) proposed recently by the senior author of the paper. Because of an intrinsic feature of the boundary element method, a general crack growth problem can be solved in a single region formulation. In the numerical simulation, for each increment of crack extension, remeshing of existing boundaries is not necessary. Crack extension is modeled conveniently by adding new boundary elements on the incremental crack extension to the previous crack boundaries. As an example, fatigue growth process of an inclined crack in an infinite plate under uniaxial cycle load is modeled to illustrate the effectiveness of the numerical approach. In addition, fatigue growth of cracks emanating from a circular hole in infinite elastic plate subjected to remote loads is investigated by using the numerical approach. Many numerical results are given     

Stable crack growth along a brittle\ductile interface—II. Small scale yielding solutions and interfacial toughness predictions

The problem of a crack growing steadily and quasi-statically along a brittle\ductile interface under plane strain, mixed mode, and small scale yielding conditions is considered. The ductile material is assumed to be characterized by the J₂-flow theory of plasticity with linear strain hardening, while the brittle material is assumed to be linear elastic. A displacement-based finite element method, exploiting the convective nature of the problem, is utilized to solve the relevant boundary value problem. In Part I of this work, the corresponding asymptotic problem was solved. This paper addresses the full-field problem in order to validate the asymptotic solutions, and to explore the physical implications of the results. The numerical full-field results are found to be in good agreement with the analytical asymptotic solutions. In particular, the full-field results strongly suggest that the stress fields in the vicinity of the crack tip are variable-separable of the power singular type; and also that the mode mix of the near-tip stress fields is, to a large extent, independent of the applied elastic mode mix. The amplitude (the plastic stress intensity factor) and the regions of validity of the asymptotic fields are estimated from the full-field results, and are observed to be strongly dependent on the applied mode mix. The remote elastic loading fields appear to influence the near-tip fields, primarily, through the plastic stress intensity factor. The present work also explores the suggestion made by Bose and Ponte Castaneda, 1992 that the solutions to the small scale yielding problem may be used in the context of a standard crack growth criterion, requiring that continued growth take place with a fixed near-tip crack opening profile, to obtain theoretical predictions for the dependence of interfacial toughness on the applied mode mix. Based on the numerical results, predictions for mixed mode toughness of the brittle\ductile interface are reported. The results, which are in qualitative agreement with available experimental data and also with some recent theoretical results, predict a strong dependence of interfacial toughness on mode mix. This suggests that ductility provides the main operating mechanism for explaining the dependence of interfacial toughness on the mode mix of the applied loading fields, during steady crack growth.     

巴西圆盘试样中心斜裂纹疲劳扩展轨迹的边界元模拟

提出了一个平面弹性体多裂纹疲劳扩展模型. 它主要涉及到复合型加载情况下多裂纹尖端疲劳扩展的数学模型及杂交位移不连续法(一种边界元法). 在数值模拟中, 对每一裂纹扩展增量分析时,在其先前的边界上增添裂纹扩展增量, 且只对新增添的裂纹扩展增量划分单元, 同时, 按照这种边界元法的实施方法对一些单元特征进行调整, 就可以方便地模拟裂纹扩展. 用这种数值方法模拟了巴西圆盘试样中心斜裂纹疲劳扩展轨迹,数值结果说明了预报模型的有效性, 揭示了裂纹体几何对疲劳扩展的影响.      

Retardation of a crack with connections between the faces using an induced thermoelastic stress field

This paper considers local temperature variations near the tip of a crack in the presence of regions in which the crack faces interact. It is assumed that these regions are adjacent to the crack tip and are comparable in size to the crack size. The problem of local temperature variations consists of delay or retardation of crack growth. For a crack with connections between the crack faces subjected to external tensile loads, an induced thermoelastic stress field, and the stresses at the connections preventing crack opening, the boundary-value problem of the equilibrium of the crack reduces to a system of nonlinear singular integrodifferential equations with a Cauchy kernel. The normal and tangential stresses at the connections are found by solving this system of equations. The stress intensity factors are calculated. The energy characteristics of cracks with tip regions are considered. The limiting equilibrium condition for cracks with tip regions is formulated using the criterion of limiting stretching of the connections.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 133–143, January–February, 2005     

弹性T项对裂尖参数和裂纹扩展路径稳定性的作用

研究了弹性T项在主裂纹与近裂尖微空洞干涉问题中对裂尖参数和裂纹扩展路径稳定性的影响作用.利用“伪力”法,并通过解决主裂纹与近裂尖微空洞干涉问题,对远场纯1型载荷和远场弹性T项载荷下,该干涉问题中弹性T项的影响作用进行了分析从数值结果可以看出：由于空洞的存在;释放了弹性T项所引起的应力,弹性T项对裂尖参数;应力强度因子和J积分都有直接显著的影响,因而,它对该载荷下的裂纹扩展路径的稳定性有控制作用。     

Optimal control of crack growth in elastic body with inclusions

In the paper, an optimal control problem of crack growth is considered, allowing to choose the most safe inclusions in elastic bodies from the standpoint of their influence on a crack propagation.     

Cohesive zone laws for fatigue crack growth: Numerical field projection of the micromechanical damage process in an elasto-plastic medium

Cohesive zone failure models are widely used to simulate fatigue crack propagation under cyclic loading, but the model parameters are phenomenological and are not closely tied to the underlying micromechanics of the problem. In this paper, we will inversely extract the cohesive zone laws for fatigue crack growth in an elasto-plastic ductile solid using a field projection method (FPM), which projects the equivalent tractions and separations at the cohesive crack-tip from field information outside the process zone. In our small-scale yielding model, a single row of discrete voids is deployed directly ahead of a crack in an elasto-plastic medium subjected to cyclic mode I K-field loading. Damage accumulation under cyclic loading is captured by the growth of voids within the micro-voiding zone ahead of the crack, while the evolution of the cohesive zone law representing the micro-voiding zone is inversely extracted via the FPM. We show that the field-projected cohesive zone law captures the essential micromechanisms of fatigue crack growth in the ductile medium: from loading and unloading hysteresis caused by void growth and plastic hardening, to the softening damage locus associated with crack propagation via a void by void growth mechanism. The results demonstrate the effectiveness of the FPM in obtaining a micromechanics-based cohesive zone law in-place of phenomenological models, which opens the way for a unified treatment of fatigue crack problems.     

Numerical simulation of debonding of adhesively bonded joint

This paper describes a numerical method to simulate the debonding of adhesively bonded joints. Assuming that the adhesive thickness and the adhesive Young’s modulus are small with respect to the characteristic length of the joint and to the Young’s modulus of the adherents, a simplified model is derived in the case of large displacements using the asymptotic expansion technique. Then, the problem of the crack growth is stated, in the case of a stable growth, as the search of the local minima of the total energy of the joint, sum of the mechanical energy and the Griffith’s fracture energy. This is made using the Newton’s method. To this end, the expressions of the first and second derivatives of the mechanical energy with respect to a crack front displacement are derived analytically. Finally, numerical examples are presented, highlighting the unstable character of the crack growth at initiation.     

Stable growth of a crack interacting with a circular inclusion

The problem of failure of a plate containing a circular inclusion and a crack is studied. The crack is oriented along a diameter of the inclusion and the plate is subjected to a remote uniaxial stress perpendicular to the crack axis. The process of slow stable crack growth from initiation to termination is studied by the strain energy density theory. The crack growth is simulated by predicting finite increments of crack extension when material elements near the crack tip absorb a critical amount of strain energy density level, . Unstable crack growth occurs when the strain energy density factor S reaches a critical value where r_c is the critical size of the final crack increment prior to instability. The stress at crack initiation and the critical stress and crack length at failure are determined. The influence of the mechanical properties of the plate and the inclusion, the relative position of the inclusion and the crack and the crack length on the characteristic quantities of stable crack growth is analyzed. The dependence of the stable crack growth process on the loading rate is also investigated. Results are displayed in graphical form.     

A non-intrusive methodology for the representation of crack growth stochastic processes

In this paper, we present a methodology to pursue the uncertainty quantification of the stochastic process that represents the crack growth problem. The main idea of this methodology is to discretize the crack growth process in a sequence of random variables and then, approximate each of them using a stochastic polynomial approach. This methodology is non-intrusive, i.e. it is based on the representation of random variables using stochastic polynomials, whose coefficients are evaluated using a least squares method and only a few realizations of the stochastic process. The Paris–Erdogan law was used as crack growth model in order to focus the reader's attention on the uncertainty quantification methodology. We modeled the parameters of the Paris–Erdogan law as random variables, i.e. the initial crack length and the coefficients of the Paris–Erdogan model are treated as random variables. Two numerical examples are presented in order to shown the effectiveness and accuracy of the proposed methodology. From the results of these examples, it is shown that the proposed methodology is able to successfully approximate the stochastic process that represents the crack growth for the Paris–Erdogan model, with a much lower computational cost than the MCS. The main limitation of the proposed approach is that, in the form it was presented, it is not able to handle random processes as input parameters.     

Initiation and growth characterization of corner cracks near circular hole

The finite element analysis of crack problems often incorporates the asymptotic character of the local solution into the formulation. Embedment of stress or strain singularities can impose serious restrictions on the outcome and inconsistencies in predicting crack and/or growth. These restrictions are discussed in connection with the problem of two diametrically opposite corner cracks near a circular hole subjected to remote uniform tension. Enforced in the numerical treatment is the 1/r character of the strain energy density function local to the corner crack border where r is the radial distance measured from the crack front. The tendency for the corner crack to become a through crack is predicted by assuming that each point of the crack border extends by an amount proportional to the strain energy density factor. The path would correspond to the loci of minimum strain energy density function. Numerical results are displayed graphically and discussed in connection with crack initiation and non-self-similar crack growth.     

Optimization formulation of the evolutionary problem of crack propagation under quasibrittle fracture

For the evolutionary problem describing crack propagation in a solid with allowance for the irreversible work of plastic deformation due to the crack propagation, a general optimization formulation is proposed and investigated. For the optimum crack, data on the H²-smoothnesses of the displacement field in the solid and, hence, on the finiteness of the stress at the crack tip, are obtained. The solvability of the optimization problem (i.e., the existence of an optimum crack) is proved for a curvilinear crack propagation path specified a priori. For the particular case of a straight path, a generalized criterion of crack growth is proposed. The question of the choice of a crack propagation path is discussed and a comparison with existing fracture criteria is made. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 5, pp. 107–118, September–October, 2006.     

An investigation of hydraulic fracturing. Part I: one-phase flow model

The prediction of the growth of a hydraulic fracture in an oil bearing formation based on the injection rate of fluid is valuable in applications of the waterflood technique in secondary oil recovery. In this paper, the problem of hydraulic fracture growth is studied under the assumption of uniform distribution of pressure in the fracture and unidirectional permeating flow in an infinitely large isothermal linearly elastic porous medium saturated with a one-phase incompressible fluid. The condition of plane strain is imposed in the study. A comparison of the constant fracture toughness criterion based on the asymptotic value for large crack growth with the crack tip ductility criterion for an ideally plastic solid under plane strain and small-scale yielding conditions indicates that the effect of ductility of rock on the crack growth is so small that the steady state value of the energy release rate can be reached within a short period of crack growth. Thus we can employ the constant fracture toughness criterion in our study. The analysis includes the effects of both fracture volume increase and leak-off of fluid from the surface of the fracture. A nonlinear singular integro-differential equation can be formulated for the quasi-static hydraulic fracture growth under a prescribed injection rate. It is solved numerically by a modified fourth order Runge-Kutta method.     

A MOVING CRACK IN A NONHOMOGENEOUS MATERIAL STRIP

This paper considers an anti-plane moving crack in a nonhomogeneous material strip of finite thickness. The shear modulus and the mass density of the strip are considered for a class of functional forms for which the equilibrium equation has analytical solutions. The problem is solved by means of the singular integral equation technique. The stress field near the crack tip is obtained. The results are plotted to show the effect of the material non-homogeneity and crack moving velocity on the crack tip field. Crack bifurcation behaviour is also discussed. The paper points out that use of an appropriate fracture criterion is essential for studying the stability of a moving crack in nonhomogeneous materials. The prediction whether the unstable crack growth will be enhanced or retarded is strongly dependent on the type of the fracture criterion used. Based on the analysis, it seems that the maximum 'anti-plane shear' stress around the crack tip is a suitable failure criterion for moving cracks in nonhomogeneous materials.