Antiplane deformation of a partially bonded elliptical inclusion

A new method that introduces two holomorphic potential functions (the two-phase potentials) is applied to analyze the antiplane deformation of an elliptical inhomogeneity partially-bonded to an infinite matrix. Elastic fields are obtained when either the matrix is subject to a uniform longitudinal shear or the inhomogeneity undergoes a uniform shear transformation. The stress field possesses the square-root singularity of a Mode III interface crack, which, in the special case of a rigid line inhomogeneity, changes in order, as the crack tip approaches the inhomogeneity end. In the latter situation the crack-tip elastic fields are linear in two real stress intensity factors related to a strong and a weak singularity of the stress field.     

增强相/夹杂内螺位错偶极子与含共焦钝裂纹椭圆夹杂的干涉效应

运用弹性力学的复势方法,研究了纵向剪切下增强相/夹杂内螺型位错偶极子与含共焦钝裂纹椭圆夹杂的干涉效应,得到了该问题复势函数的封闭形式解答,由此推导出了夹杂区域的应力场、作用在螺型位错偶极子中心的像力和像力偶矩以及裂纹尖端应力强度因子级数形式解。并分析了位错偶极子倾角 、钝裂纹尺寸和材料常数对位错像力、像力偶矩以及应力强度因子的影响。数值计算结果表明：位错像力、像力偶矩以及应力强度因子均随位错偶极子倾角做周期变化;夹杂内部的椭圆钝裂纹明显增强了硬基体对位错的排斥,减弱了软基体对位错的吸引,且对于硬夹杂,位错出现了一个不稳定平衡位置,该平衡位置随钝裂纹曲率的增大不断向界面靠近;变化 值将出现改变位错偶极子对应力强度因子作用方向的临界值。     

Stress intensity factor of an elliptical crack as influenced by a spherical inhomogeneity

The interaction between an elliptical crack and a spherical inhomogeneity embedded in a three-dimensional solid subject to uniaxial tension is investigated. Both the inhomogeneity and the solid are isotropic but have different elastic moduli. The Eshelby's equivalent inclusion method is applied together with the principle of superposition. An approximate solution for the stress intensity factor is obtained by an approach that expands the distance between the center of the crack and inhomogeneity in series. The local stress field can be increased or decreased depending on the relative modulus of the spherical inhomogeneity and matrix. If the inhomogeneity modulus is larger than that of the matrix, a reduction in the stress intensity factor prevails. Displayed numerically are results to exhibit the influence of inhomogeneity and its distance to the crack.     

Edge dislocation interacting with an interfacial crack along a circular inhomogeneity

The elastic interaction of an edge dislocation, which is located either outside or inside a circular inhomogeneity, with an interfacial crack is dealt with. Using Riemann–Schwarz’s symmetry principle integrated with the analysis of singularity of the complex potentials, the closed form solutions for the elastic fields in the matrix and inhomogeneity regions are derived explicitly. The image force on the dislocation is then determined by using the Peach–Keohler formula. The influence of the crack geometry and material mismatch on the dislocation force is evaluated and discussed when the dislocation is located in the matrix. It is shown that the interfacial crack has significant effect on the equilibrium position of the edge dislocation near a circular interface. The results also reveal a strong dependency of the dislocation force on the mismatch of the shear moduli and Poisson’s ratios between the matrix and inhomogeneity.     

位错偶极子和Ⅲ型界面裂纹的干涉效应

研究了多晶体材料中螺型位错偶极子和界面裂纹的弹性干涉作用.利用复变函数方法,得到了该问题复势函数的封闭形式解答.求出了由位错偶极子诱导的应力场和裂纹尖端应力强度应子,分析了偶极子的方向,偶臂和位置以及材料失配对应力强度因子的影响.推导了作用在螺型位错偶极子中心的像力和力偶矩,并讨论了界面裂纹几何条件和不同材料特征组合对位错偶极子平衡位置的影响规律.结果表明,裂纹尖端的螺型位错偶极子对应力强度因子会产生强烈的屏蔽或反屏蔽效应.同时,界面裂纹对螺型位错偶极子在材料中运动有很强的扰动作用.     

On the anti-plane shear of an elliptic nano inhomogeneity

The elastic field of an elliptic nano inhomogeneity embedded in an infinite matrix under anti-plane shear is studied with the complex variable method. The interface stress effects of the nano inhomogeneity are accounted for with the Gurtin–Murdoch model. The conformal mapping method is then applied to solve the formulated boundary value problem. The obtained numerical results are compared with the existing closed form solutions for a circular nano inhomogeneity and a traditional elliptic inhomogeneity under anti-plane. It shows that the proposed semi-analytic method is effective and accurate. The stress fields inside the inhomogeneity and matrix are then systematically studied for different interfacial and geometrical parameters. It is found that the stress field inside the elliptic nano inhomogeneity is no longer uniform due to the interface effects. The shear stress distributions inside the inhomogeneity and matrix are size dependent when the size of the inhomogeneity is on the order of nanometers. The numerical results also show that the interface effects are highly influenced by the local curvature of the interface. The elastic field around an elliptic nano hole is also investigated in this paper. It is found that the traction free boundary condition breaks down at the elliptic nano hole surface. As the aspect ratio of the elliptic hole increases, it can be seen as a Mode-III blunt crack. Even for long blunt cracks, the surface effects can still be significant around the blunt crack tip. Finally, the equivalence between the uniform eigenstrain inside the inhomogeneity and the remote loading is discussed.     

含圆形嵌体弹性平面中径向裂纹问题的超奇异积分方程方法

根据含圆形嵌体平面问题在极坐标下的弹性力学基本解,使用Betti互换定理,在有限部积分意义下将问题归结为两个以裂纹岸位移间断为基本未知量、对于Ⅰ型和Ⅱ型问题相互独立的超奇异积分方程,对含圆形嵌体弹性平面中的径向裂纹问题进行了研究.根据有限部积分原理,建立了问题的数值算法.计算结果表明,嵌体半径、裂纹位置及材料剪切弹性模量等都对裂纹应力强度因子具有较为明显的影响.     

On the nucleation of a Zener crack from a wedge disclination dipole in the presence of a circular inhomogeneity

The nucleation of a Mode-I Zener crack from a wedge disclination dipole in the presence of a circular inhomogeneity is investigated. It is assumed that the disclination dipole and the nucleated Zener crack are along the radial direction of the inhomogeneity. Two cases are studied herein, i.e., the positive or negative wedge disclination of the dipole locates nearer to the inhomogeneity respectively. In order to investigate how various factors such as the elastic mismatch between the inhomogeneity and the matrix influence the nucleation of the Zener crack, the Stress Intensity Factor (SIF) at the sharp tip of the Zener crack is determined for different sets of geometric and material parameters with the distributed dislocation technique. The formulated singular integral equations are then solved numerically. Our results indicate that a nearby ‘hard’ inhomogeneity (having a higher shear modulus than the matrix) is beneficial to the crack nucleation for the first case (the positive disclination locates nearer to the inhomogeneity) while it retards the crack nucleation for the second case (the negative disclination locates nearer to the inhomogeneity). A nearby ‘soft’ inhomogeneity is helpful to the crack nucleation for the second case while it has inverse effects on the crack nucleation for the first case. This phenomenon can be explained with the concept of material force. The characteristics of the crack nucleation and the effects of the disclination strength, the distance between the inhomogeneity and the dipole, the disclination dipole arm length and inhomogeneity size on the crack nucleation are also systematically studied. The obtained results are helpful to characterize and enhance the strength of precipitate alloys and particle reinforced composites.     

界面裂纹问题中的弹性T项和应力强度因子

研究两相材料有限板含单边界面裂纹的断裂力学特性，对不同的材料组合用广义变分法分析了不同尺寸试件和裂纹长度下的应力强度因子和弹性T项，讨论了材料特性对应力强度因子和弹性T项的作用．分析了试件尺寸和裂纹长度对应力强度因子和弹性T项的影响．     

Dynamic interaction between a matrix crack and a circular inhomogeneity with a distinct interphase

This article provides a theoretical treatment of the dynamic interaction between a matrix crack and an arbitrarily located circular inhomogeneity with a distinct interphase under antiplane loading. The matrix⧹inhomogeneity interphase is characterized by a linear spring model. The theoretical formulations governing the steady state problem are based upon the use of integral transform techniques, Bessel function expansions and a Pseudo-incident wave technique. The closed form expression for the resulting stress intensity factor at the matrix crack is obtained by solving the appropriate singular integral equations using Chebyshev polynomials. Typical examples are provided to show the effect of the location of the inhomogeneity, the material combination and the interface property upon the dynamic stress intensity factor of the matrix crack.     

A general solution method for an anti-plane shear crack dynamically accelerating along a bimaterial interface

We develop a general solution method for a dynamically accelerating crack under anti-plane shear conditions along the interface between two different homogeneous isotropic elastic materials. The crack is initially at rest, and after loading is applied the crack-tip speed which may accelerate up to the shear wave speed of the more compliant material. The analysis includes an exact, closed-form expression for the stress intensity factor for an arbitrary time-dependent crack-face traction, as well as expressions for computing the crack-face displacements and the stress in front of the crack. We also present some numerical examples for fixed loads and for loads moving with the crack tip, using a stress intensity factor fracture criterion, in order to examine the predicted effect of material mismatch on interfacial fracture.     

On a semi-infinite crack penetrating a piezoelectric circular inhomogeneity with a viscous interface

We investigate a semi-infinite crack penetrating a piezoelectric circular inhomogeneity bonded to an infinite piezoelectric matrix through a linear viscous interface. The tip of the crack is at the center of the circular inhomogeneity. By means of the complex variable and conformal mapping methods, exact closed-form solutions in terms of elementary functions are derived for the following three loading cases: (i) nominal Mode-III stress and electric displacement intensity factors at infinity; (ii) a piezoelectric screw dislocation located in the unbounded matrix; and (iii) a piezoelectric screw dislocation located in the inhomogeneity. The time-dependent electroelastic field in the cracked composite system is obtained. Particularly the time-dependent stress and electric displacement intensity factors at the crack tip, jumps in the displacement and electric potential across the crack surfaces, displacement jump across the viscous interface, and image force acting on the piezoelectric screw dislocation are all derived. It is found that the value of the relaxation (or characteristic) time for this cracked composite system is just twice as that for the same fibrous composite system without crack. Finally, we extend the methods to the more general scenario where a semi-infinite wedge crack is within the inhomogeneity/matrix composite system with a viscous interface.     

双材料界面裂纹应力强度因子的边界元分析

采用双材料基本解建立边界元法基本方程,计算双材料界面裂纹尖端附近的应用力和位移场。不离散界面,并设置面力奇异四分之一点裂尖单元以提高计算精度。数值结果表明,本文的方法具有较高的精度和效率。     

Three coplanar Griffith cracks in an infinite elastic strip

We consider the problem of determining the stress intensity factors and crack energy in an infinitely long isotropic, homogeneous elastic strip containing three coplanar Griffith cracks. The three coplanar Griffith cracks are situated symmetrically on a line perpendicular to the edges of the strip. We assume that the cracks are opened by an internal pressure and the edges of the strip are fixed. By using the theory of Fourier series we reduce the problem to solving a set of quadruple trigonometrical series equations with a cosine kernel. Closed form solution is obtained for the quadruple series equations. Closed form analytical expressions are derived for the stress intensity factors, the shape of the deformed cracks and the crack energy. Solutions to some particular problems are derived as limiting cases.     

Computing bounds to mixed-mode stress intensity factors in elasticity

A bounding procedure combined with an effective error bound method for linear functionals of the displacements and a simple two points displacement extrapolation method is presented to compute the lower and upper bounds to the stress intensity factors in elastic fracture problems. First, the displacements of two nodes (or node pairs) on the crack edges are used to construct the linear extrapolation to obtain the stress intensity factors at the crack tip, so that stress intensity factors are explicitly expressed as linear functionals of the displacements. Then, a posteriori bounding method is utilized to compute the bounds to the stress intensity factors. Finally, the bounding procedure is verified by a mixed-mode homogenous elastic fracture problem and a bimaterial interface crack problem.     

Evaluation of fracture mechanics parameters for free edges in multi-layered structures with weak singularities

This study focuses on the stress intensity factors for free edges in multi-layered structural components. The effects of elastic constants of various material combinations on the weak singularity at free edges are analyzed. Using the H-integral approach, the effects of elastic mismatch parameters, the bond area and the thickness of the thin metal layer on the stress intensity factor are quantified and the results are compared with detailed finite element solutions. A good agreement between numerical predictions obtained from the H-integral and the detailed FE results is achieved, showing the applicability of this approach. Similar to a crack problem, the singular elastic field dominates in an annular region adjacent to the notch tip. The relationship between the valid range of the K-dominated field and the thin-film thickness is then demonstrated. Furthermore, the competition of crack initiation between the free edge interface (180° opening angle) and a 90° notch interface in a generic specimen is investigated, in order to find out which is the prevailing failure mode. Comparison between isotropic Si and anisotropic Si substrate is also illustrated. Anisotropy of the Si substrate has a significant influence on the stress intensity factor when combined with an Au or Al metal layer but not with a Cu layer. Additionally, standardized numerical formulae of the dimensionless stress intensity factor have been derived to guide the engineering application.     

Electroelastic interaction between a piezoelectric screw dislocation and a confocal elliptic blunt crack with elliptical inhomogeneity

The electroelastic interaction between a piezoelectric screw dislocation and an elliptical inhomogeneity containing a confocal blunt crack under infinite longitudinal shear and in-plane electric field is investigated. Using the sectionally holomorphic function theory, Cauchy singular integral, singularity analysis of complex functions and theory of Rieman boundary problem, the explicit series solution of stress field is obtained when the screw dislocation is located in inhomogeneity. The intervention law of the interaction between blunt crack and screw dislocation in inhomogeneity is discussed. The analytical expressions of generalized stress and strain field of inhomogeneity are calculated, while the image force, field intensity factors of blunt crack are also presented. Moreover, a new matrix expression of the energy release rate and generalized strain energy density (SED) are deduced. With the size variation of blunt crack, the results can be reduced to the case of the interaction between a piezoelectric screw dislocation and a line crack in inhomogeneity. Numerical analysis are then conducted to reveal the effects of the dislocation location, the size of inhomogeneity and blunt crack and the applied load on the image force, energy release rate and strain energy density. The influence of dislocation on energy release rate and strain energy density is also revealed.     

The deformation of the front of a 3D interface crack propagating quasistatically in a medium with random fracture properties

One studies the evolution in time of the deformation of the front of a semi-infinite 3D interface crack propagating quasistatically in an infinite heterogeneous elastic body. The fracture properties are assumed to be lower on the interface than in the materials so that crack propagation is channelled along the interface, and to vary randomly within the crack plane. The work is based on earlier formulae which provide the first-order change of the stress intensity factors along the front of a semi-infinite interface crack arising from some small but otherwise arbitrary in-plane perturbation of this front. The main object of study is the long-time behavior of various statistical measures of the deformation of the crack front. Special attention is paid to the influences of the mismatch of elastic properties, the type of propagation law (fatigue or brittle fracture) and the stable or unstable character of 2D crack propagation (depending on the loading) upon the development of this deformation.     

Stress intensity factor analysis of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress

A numerical method using a path-independent H-integral based on the conservation integral was developed to analyze the singular stress field of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress. In the present method, the shape of the corner front is smooth. According to the theory of linear elasticity, asymptotic stress near the tip of a sharp interfacial corner is generally singular as a result of a mismatch of the materials’ elastic constants. The eigenvalues and the eigenfunctions are obtained using the Williams eigenfunction method, which depends on the anisotropic materials’ properties and the geometry of an interfacial corner. The order of the singularity related to the eigenvalue is real, complex or power-logarithmic. The amplitudes of the singular stress terms can be calculated using the H-integral. The stress and displacement around an interfacial corner for the H-integral are obtained using finite element analysis. In this study, a proposed definition of the stress intensity factors of an interfacial corner, which includes those of an interfacial crack and a homogeneous crack, is used to evaluate the singular stress fields. Asymptotic solutions of stress and displacement around an interfacial corner front are uniquely obtained using these stress intensity factors. To prove the accuracy of the present method, several different kinds of examples are shown such as interfacial corners or cracks in three-dimensional structures.     

Interaction of an Ellipsoidal Inclusion with an Elliptic Crack in an Elastic Material under Triaxial Tension

The interaction of an elastic ellipsoidal inclusion with an elliptic crack in an infinite elastic medium under triaxial loading is analyzed. The stress state in the elastic space is represented as a superposition of the principal state and perturbed states, which are due to the presence and interaction of the inclusion and the crack. The analytical solution of the problem is found using the method of equivalent inclusion, the potential of an inhomogeneous ellipsoid, and a system of harmonic functions for an elliptic crack. The effect of triaxial loading on the stress intensity factors is analyzed