Some problems in the antiplane shear deformation of bi-material wedges

The antiplane shear deformation of a bi-material wedge with finite radius is studied in this paper. Depending upon the boundary condition prescribed on the circular segment of the wedge, traction or displacement, two problems are analyzed. In each problem two different cases of boundary conditions on the radial edges of the composite wedge are considered. The radial boundary data are: traction–displacement and traction–traction. The solution of governing differential equations is accomplished by means of finite Mellin transforms. The closed form solutions are obtained for displacement and stress fields in the entire domain. The geometric singularities of stress fields are observed to be dependent on material property, in general. However, in the special case of equal apex angles in the traction–traction problem, this dependency ceases to exist and the geometric singularity shows dependency only upon the apex angle. A result which is in agreement with that cited in the literature for bi-material wedges with infinite radii. In part II of the paper, Antiplane shear deformation of bi-material circular media containing an interfacial edge crack is considered. As a special case of bi-material wedges studied in part I of the paper, explicit expressions are derived for the stress intensity factor at the tip of an edge crack lying at the interface of the bi-material media. It is seen that in general, the stress intensity factor is a function of material property. However, in special cases of traction–traction problem, i.e., similar materials and also equal apex angles, the stress intensity factor becomes independent of material property and the result coincides with the results in the literature.     

Analysis of an Anisotropic Finite Wedge under Antiplane Deformation

The antiplane deformation of an anisotropic wedge with finite radius is considered in this paper within the classical linear theory of elasticity. The traction-free condition is imposed on the circular segment of the wedge. Three different cases of boundary conditions on the radial edges are considered, which are: traction-displacement, displacement-displacement and traction-traction. The solution to the governing differential equation of the problem is accomplished in the complex plane by relating the displacement field to a complex function. Several complex transformations are defined on this complex function and its first and second derivatives to formulate the problem in each of the three cases of the problem corresponding to the radial boundary conditions, separately. These transformations are then related to integral transforms which are complex analogies to the standard finite Mellin transforms of the first and second kinds. Closed form expressions are obtained for the displacement and stress fields in the entire domain. In all cases, explicit expressions for the strength of singularity are derived. These expressions show the dependence of the order of stress singularity on the wedge angle and material constants. In the displacement-displacement case, depending upon the applied displacement, a new type of stress singularity has been observed at the wedge apex. This revised version was published online in August 2006 with corrections to the Cover Date.     

Electroelastic analysis of a piezoelectric finite wedge with mixed type boundary conditions under a pair of concentrated shear forces and free charges

The antiplane problem of a piezoelectric finite wedge with mixed type boundary conditions (i.e. traction free-electrically grounded or clamped-electrically insulated) assumed on the circular edge is studied in this paper. Using the first and second kinds of the finite Mellin transforms, the stress and electrical displacement in the piezoelectric finite wedge subjected to a pair of concentrated shear forces and free charges are derived analytically. The singularity orders and all generalized intensity factors can also be obtained. After being reduced to the crack problem, the energy density theory is applied to study the effects of the boundary conditions on the circular edge and the wedge radius on the fracture behaviors of the cracked piezoelectric medium. The problem can be further degenerated to several simple problems and the results agree well with those of previous studies.     

Analysis of anisotropic sector with a radial crack under anti-plane shear loading

In this paper, the anti-plane shear deformation of an anisotropic sector with a radial crack is investigated. The traction–traction boundary conditions are imposed on the radial edges and the traction-free condition is considered on the circular segment of the sector. A novel mathematical technique is employed for the solution of the problem. This technique consists of the use of some recently proposed finite complex transforms (Shahani, 1999), which have complex analogies to the standard finite Mellin transforms of the first and second kinds. However, it is essential to state the traction-free condition of the crack faces in the form of a singular integral equation which is done in this paper by describing an exact analytical method. The resultant dual integral equations are solved numerically to determine the stress intensity factors at the crack tips. In the special cases, the obtained results coincide with those cited in the literature.     

Numerical methods for solving singular integral equations obtained by fracture mechanical analysis of cracked wedge简

The numerical solutions to the singular integral equations obtained by the fracture mechanical analyses of a cracked wedge under three different conditions are considered. The three considered conditions are: (i) a radial crack on a wedge with a non-finite radius under the traction-traction boundary condition, (ii) a radial crack on a wedge with a finite radius under the traction-traction boundary condition, and (iii) a radial crack on a finite radius wedge under the traction-displacement boundary condition. According to the boundary conditions, the extracted singular integral equations have different forms. Numerical methods are used to solve the obtained coupled singular integral equations, where the Gauss-Legendre and the Gauss-Chebyshev polynomials are used to approximate the responses of the singular integral equations. The results are presented in figures and compared with those obtained by the analytical response. The results show that the obtained Gauss-Chebyshev polynomial response is closer to the analytical response.     

Antiplane electro-mechanical field of a piezoelectric finite wedge under shear loading and at fixed-grounded boundary conditions

This paper presents the general solutions of antiplane electro-mechanical field solutions for a piezoelectric finite wedge subjected to a pair of concentrated forces and free charges. The boundary conditions on the circular segment are considered as fixed and grounded. Employing the finite Mellin transform method, the stress and electrical displacement at all fields of the piezoelectric finite wedge are derived analytically. In addition, the singularity orders and intensity factors of stress and electrical displacement can also be obtained. These parameters can be applied to examine the fracture behavior of the wedge structure. After being reduced to the problem of an antiplane edge crack or an infinite wedge in a piezoelectric medium, the results compare well with those of previous studies.     

Boundary element analysis of interaction between an elastic rectangular inclusion and a crack

The interaction between an elastic rectangular inclusion and a kinked crack inan infinite elastic body was considered by using boundary element method. The new complexboundary integral equations were derived. By introducing a complex unknown function H(t)related to the interface displacement density and traction and applying integration by parts,the traction continuous condition was satisfied automatically. Only one complex boundaryintegral equation was obtained on interface and involves only singularity of order l/ r. Toverify the validity and effectiveness of the present boundary element method, some typicalexamples were calculated. The obtained results show that the crack stress intensity factorsdecrease as the shear modulus of inclusion increases. Thus, the crack propagation is easiernear a softer inclusion and the harder inclusion is helpful for crack arrest.     

界面下主微裂纹间干涉的研究

涉及两相正交各向异性体界面干涉问题的研究,多裂纹问题被分解为只含单裂纹的子问题,利用位错理论和裂面应力自由条件,列出一组可数值求解位错密度函数的奇异积分方程,从耐 注得应力强度因子。     

Cracking of an orthotropic substrate reinforced by an orthotropic plate

A solution method is derived to determine the stress intensity factors for both an internal crack and an edge crack in an orthotropic substrate that is reinforced on its boundary by a finite-length orthotropic plate. The method utilizes the Green’s functions for a pair of dislocations and a concentrated force on the boundary while invoking the concept of superposition. Enforcing the traction-free boundary condition along the crack surfaces and the continuity of displacement gradients along the plate/substrate interface results in a coupled system of singular integral equations. An asymptotic analysis of the kernels in these equations for the region of the junction point between the plate corner and the substrate boundary reveals the strength of the singularity in the case of an edge crack. The numerical solution of the integral equations provides results for the stress intensity factors for both an internal crack and an edge crack perpendicular to the substrate boundary and aligned with one of the corners of the plate. The present results have been validated against previously published stress intensity factors for an internal crack and an edge crack in an isotropic substrate.     

A contour integral method for dynamic stress intensity factors

The contour integral method previously used to determine static stress intensity factors is applied to dynamic crack problems. The required derivatives of the traction in the reference problem are obtained numerically by the displacement discontinuity method. Stress intensity factors are determined by an integral around a contour which contains a crack tip. If the contour is chosen as the outer boundary of the body, the stress intensity factor is obtained from the boundary values of traction and displacement. The advantage of this path-independent integral is that it yields directly both the opening-mode and sliding-mode stress intensity factors for a straight crack. For dynamic problems, Laplace transforms are used and the dynamic stress intensity factors in the time domain are determined by Durbin's inversion method. An indirect boundary element method, incorporating both displacement discontinuity and fictitious load techniques, is used to determine the boundary or contour values of traction and displacement numerically.     

Singular stress fields of angle-ply and monoclinic bimaterial wedges

The characteristic equations for the order of stress singularity of anisotropic bimaterial wedges subjected to traction boundary conditions are investigated. For an angle-ply bimaterial wedge, both fully bonded and frictional interfaces are considered, whereas for a monoclinic bimaterial wedge, a frictional interface is considered. Here, the Stroh formalism and the separation of variables technique are used. In general, the order of stress singularity can be real or complex, but for the special geometry of a crack along the frictional interface of a monoclinic composite, it is always real. Explicit characteristic equations for the order of singularity are presented for an aligned orthotropic composite with a frictional interface. Numerical results are given for an angle-ply bimaterial wedge and a monoclinic bimaterial wedge consisting of a graphite/epoxy fiber-reinforced composite.     

同震断层三维裂纹扩展波动时域超奇异积分法

应用波动时域超奇异积分法将P波、S波和磁电热弹多场耦合作用下同震断层任意形状三维裂纹扩展问题转化为求解以广义位移间断率为未知函数的超奇异积分方程组问题;定义了广义应力强度因子,得到裂纹前沿广义奇异应力增量解析表达式;应用波动时域有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以三维矩形裂纹扩展问题为例,通过典型算例,研究了广义应力强度因子随裂纹位置变化规律;分析了同震断层裂纹扩展中力、磁、电场辐射规律.      

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

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

A Long Crack Penetrating a Circular Inhomogeneity

This paper presents the effects of elastic mismatch and crack-tip position on the stress intensity factors of a long crack penetrating a circular inhomogeneity. The analysis relies on closed-form solutions, derived using complex variable techniques, for the stresses and displacements produced by dislocations positioned inside and outside the inhomogeneity. Dislocation distributions are introduced to express the traction boundary condition along the crack surfaces as a system of singular integral equations, whose solution is obtained through a numerical procedure. It is shown that if the elastic mismatch is interpreted correctly, then the stress intensity factors of this micromechanical model are very good approximations to those computed using a Monte Carlo finite element model of a long crack in a polycrystalline plate with compliant grain boundaries.     

Higher order fields for damaged nonlinear antiplane shear notch, crack and inclusion problems

Damaged nonlinear antiplane shear problems with a variety of singularities are studied analytically. A deformation plasticity theory coupled with damage is employed in analysis. The effect of microscopic damage is considered in terms of continuum damage mechanics approach. An exact solution for the general damaged nonlinear singular antiplane shear problem is derived in the stress plane by means of a hodograph transformation, then corresponding higher order asymptotic solutions are obtained by reversing the stress plane solution to the physical plane. As example, traction free sharp notch and crack, rigid sharp wedge and flat inclusion, and mixed boundary sharp notch problems are investigated, respectively. Consequently, higher order fields are obtained, in which analytical expressions of the dominant and second order singularity exponents and angular distribution functions of the near tip fields are derived. Effects of the damage and hardening exponents of materials and the geometric angle of notch/wedge on the near tip quantities are discussed in detail. It is found that damage leads to a weaker dominant singularity of stress, but to little stronger singularities of the dominant and second order terms of strain compared to that for undamaged material. It is also seen that damage has important effect on the angular distribution functions of the near tip stress and strain fields. As special cases, higher order analytical solutions of the crack and rigid flat inclusion tip fields are obtained, respectively, by reducing the notch/wedge tip solutions. Effects of damage and hardening exponents on the dominant and second order terms in the solutions of the crack and inclusion tip fields are discussed.     

Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

INTERFACIAL CRACK ANALYSIS IN THREE-DIMENSIONAL TRANSVERSELY ISOTROPIC BI-MATERIALS BY BOUNDARY INTEGRAL EQUATION METHOD

The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropic materials.     

三维体内含垂直于双材料界面混合型裂纹的超奇异积分解法

利用双材料位移基本解和Somigliana公式,将三维体内含垂直于双材料界面混合型裂纹问题归结为求解一组超奇异积分方程。使用主部分析法,通过对裂纹前沿应力奇性的分析,得到用裂纹面位移间断表示的应力强度因子的计算公式,进而利用超奇异积分方程未知解的理论分析结果和有限部积分理论,给出了超奇异积分方程的数值求解方法。最后,对典型算例的应力强度因子做了计算,并讨论了应力强度因子数值结果的收敛性及其随各参数变化的规律。     

压电材料平面裂纹尖端场的杂交应力有限元分析

基于复势理论和杂交变分原理建立了一种适用于力电耦合分析的杂交应力有限元模 型. 给出了建立刚度矩阵的主要公式和推导过程，单元内的位移场和应力场采用满足平 衡方程的复变函数级数解，假设的复变函数级数解事先精确满足裂纹的无应力和电位移法向 分量为零的条件，单元外边界的位移场假设按抛物线变化, 单元的刚度矩阵采用Gauss积分的方法得出. 通过对力电耦合裂尖场的数值计算验证了程序 的正确性和单元的有效性，同时也用所得结果校验了理论解.     

三维间断位移法及强奇异和超奇异积分的处理方法

从积分方程Somigliana等式出发,导出三维状态下单位位错集度的基本解.在此基础上,建立了边界积分方程,并给出了其离散形式.对强奇异和超奇异积分,采用了Hadamard定义的有限部分积分来处理.最后,给出了计算裂纹应力强度因子的算例,并与解析解进行了比较,证实了该方法的有效性.