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 dislocation and its applications to interface cracks in anisotropic bimaterials

Interfacial dislocations and cracks in anisotropic bimaterials are considered. The displacement and the stress fields due to an interfacial dislocation are obtained in a real and simple form. Explicit solutions to the traction along the interface and the crack opening displacement for a Griffith interface crack are derived. Possible definitions of stress intensity factors are given which reduce to the classical definition for a crack in a homogeneous medium. It is found that a planar interface between dissimilar anisotropic solids is completely characterized by no more than 9 independent parameters. Some invariant properties of the dislocation and crack solutions under coordinate transformation are also discussed.     

Green's functions for anisotropic piezoelectric bimaterials with a slipping interface

An explicit full-field expression of the Green's functions for anisotropic piezoelectric bimaterials with a slipping interface is derived. When the electro-elastic singularity reduces to a pure dislocation in displacement and electric potential, interaction energy between the dislocation and the bimaterials is obtained explicitly while the generalized force on the dislocation is given in a real form which is also valid for degenerate materials. The investigation demonstrates that the boundary conditions at the slipping interface between two piezoelectric materials will exert a prominent influence on the mobility of the dislocation. Project supported by the National Natural Science Foundation of China (No. 59635140).     

Fracture Properties of Concrete–Concrete Interfaces Using Digital Image Correlation

The mode I and mode II fracture toughness and the critical strain energy release rate for different concrete–concrete jointed interfaces are experimentally determined using the Digital Image Correlation technique. Concrete beams having different compressive strength materials on either side of a centrally placed vertical interface are prepared and tested under three-point bending in a closed loop servo-controlled testing machine under crack mouth opening displacement control. Digital images are captured before loading (undeformed state) and at different instances of loading. These images are analyzed using correlation techniques to compute the surface displacements, strain components, crack opening and sliding displacements, load-point displacement, crack length and crack tip location. It is seen that the CMOD and vertical load-point displacement computed using DIC analysis matches well with those measured experimentally.     

各向异性双材料界面裂纹扩展裂尖场

应用界面断裂力学理论和Stroh方法,研究了广义平面变形下动态裂纹沿着各向异性双材料界面扩展时的裂尖奇异应力及动态应力强度因子.双材料界面的动态裂尖区域特性主要由两个实矩阵W和D确定,且裂尖奇异应力和动态应力强度因子可以由包含这两个矩阵的柯西奇异积分方程确定,同时给出了动态应力强度因子和能量释放率的显示表达式.算例得出当裂纹以小速度扩展时,裂尖振荡因子ε与静态时几乎相同,当界面裂纹扩展速度接近瑞利波速时,ε趋于无穷大;同时得出应力强度因子及能量释放率随裂纹扩展速度的变化关系.     

Energy release rate and contact zone in a cohesive and an interface crack by hypersingular integral equations

Solution of Cauchy-type singular integral equations permits the evaluation of the fracture parameters at the crack tips very accurately. However, it does not permit the determination of the crack opening and sliding displacements while ensuring no crack surface interpenetration unless the location of the contact zone is known a priori. In order to circumvent this shortcoming, this study presents a solution method based on the Hadamard-type singular integral equations to obtain the crack opening and sliding displacements directly while enforcing the appropriate conditions to prevent interpenetration. Furthermore, the crack opening displacements are physically more meaningful and readily validated against the finite element analysis predictions. The numerical solutions of the hypersingular integral equations provide not only crack opening and sliding displacements directly but also the stress intensity factors and energy release rates. Also, the behavior of the energy release rate is examined as the cohesive crack located parallel to the interface approaches the interface from either the soft or the stiff side of the interface. The limiting value of the energy release rate is established by considering an interface crack. As the cohesive crack approaches the interface from either side of the interface, the energy release rate approaches to that of the interface crack. However, the length of contact zone between the cohesive crack surfaces under uniform shear loading does not approach to that of the interface crack.     

垂直磁压电材料界面三维裂纹的超奇异积分法

应用有限部积分概念和广义位移基本解，垂直于磁压电双材料界面三维复合型裂纹问题被转 化为求解一组以裂纹表面广义位移间断为未知函数的超奇异积分方程问题. 进而，通过主部 分析法精确地求得裂纹尖端光滑点附近的奇性应力场解析表达式. 然后，通过将裂纹表面 位移间断未知函数表达为位移间断基本密度函数与多项式之积，使用有限部积分法对超奇异 积分方程组建立了数值方法. 最后，通过典型算例计算，讨论了广义应力强度因子的变化规 律.     

Analysis method of planar interface cracks of arbitrary shape in three-dimensional transversely isotropic magnetoelectroelastic bimaterials

Using the fundamental solutions for three-dimensional transversely isotropic magnetoelectroelastic bimaterials, the extended displacements at any point for an internal crack parallel to the interface in a magnetoelectroelastic bimaterial are expressed in terms of the extended displacement discontinuities across the crack surfaces. The hyper-singular boundary integral–differential equations of the extended displacement discontinuities are obtained for planar interface cracks of arbitrary shape under impermeable and permeable boundary conditions in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. An analysis method is proposed based on the analogy between the obtained boundary integral–differential equations and those for interface cracks in purely elastic media. The singular indexes and the singular behaviors of near crack-tip fields are studied. Three new extended stress intensity factors at crack tip related to the extended stresses are defined for interface cracks in three-dimensional transversely isotropic magnetoelectroelastic bimaterials. A penny-shaped interface crack in magnetoelectroelastic bimaterials is studied by using the proposed method.The results show that the extended stresses near the border of an impermeable interface crack possess the well-known oscillating singularity r^?1/2±iε or the non-oscillating singularity r^?1/2±κ. Three-dimensional transversely isotropic magnetoelectroelastic bimaterials are categorized into two groups, i.e., ε-group with non-zero value of ε and κ-group with non-zero value of κ. The two indexes ε and κ do not coexist for one bimaterial. However, the extended stresses near the border of a permeable interface crack have only oscillating singularity and depend only on the mechanical loadings.     

Debonding of graded coatings under in-plane compression

Graded materials are multiphase composites with continuously varying thermophysical properties. The concept provides material scientists and engineers with an important tool to develop new materials tailored for some specific applications. One such application of this new class of materials is as top coats or interfacial regions in thermal barrier systems. A widely observed failure mode in these layered materials is known to be interfacial cracking that leads to spallation. In many cases it is the buckling instability of coating under mechanically or thermally induced compressive stresses that triggers spallation. Under in-plane loading since the linear elastic small deformation theory gives only a trivial solution, in this study the plane strain interface crack problem for a graded coating bonded to a homogeneous substrate is formulated by using a kinematically nonlinear continuum theory. Both the instability and the postbuckling problems are considered. The main objective of the study is the investigation of the influence of material nonhomogeneity, kinematic nonlinearity and plate approximation on the critical instability load and on such fracture mechanics parameters as strain energy release rate, stress intensity factors and crack opening displacements.     

Green’s function of bimaterials comprising all cases of material degeneracy

For bimaterials with planar interfaces subjected to a line force and dislocation, Green’s functions are determined for all types of anisotropic materials including the nondegenerate, degenerate and extra-degenerate cases. The changes in Green’s function caused by material degeneracy are twofold: (i) implicit changes, attributable to material effects only and characterized by high-order eigenvectors and their intrinsic coupling in the higher-order eigensolutions; (ii) explicit changes, influenced by boundary and interface conditions, that cause additional terms in Green’s function. Material degeneracy affects the angular variation of the singular stress field, which may have significant implication on the failure prediction of strongly anisotropic materials. For all material types, Green’s functions are obtained for bimaterials with a planar interface, and for multi-material wedges subjected to a line force and dislocation at the vertex. The results are expressed in a concise notation in terms of the complete set of eigenvectors and kernel matrices of analytic functions.     

十次对称二维准晶中的圆弧形裂纹

采用复变函数法研究了在远场承受均布声子场机械载荷作用下十次对称二维准晶中穿透周期轴的一个圆弧形裂纹；报导了表征声子场和相位子场的4个复应力函数的解析表达式，并由复应力函数获得了两个裂尖处声子场和相位子场应力强度因子、裂面张开位移以及能量释放率的表达式．     

A DECAGONAL QUASICRYSTAL WITH AN ARC-SHAPED CRACK

A decagonal quasicrystal,which is weakened by an arc-shaped crack penetratingthrough the solid in the period direction,and which is subjected to remote uniform phonon stresses,is investigated by applying the complex variable method which is just developed by the authors.It is found that the phonon and phason stresses near the crack tips exhibit inverse square rootsingularities.The four complex stress functions characterizing the phonon and phason fields arederived.Explicit expressions for the phonon and phason stress intensity factors,crack openingdisplacements and energy release rate are also presented.     

On the crack-tip stress singularity of interfacial cracks in transversely isotropic piezoelectric bimaterials

In this paper, characteristics of the interface crack-tip stress and electric displacement fields in transversely isotropic piezoelectric bimaterials are studied. The authors have proven, within the framework of the generalized Stroh formalism for piezoelectric bimaterials, that there is no coexistence of the parameters (oscillating) and κ (non-oscillating) in the interface crack-tip generalized stress field for all transversely isotropic piezoelectric bimaterials. This leads to the classification of piezoelectric bimaterials into one group that exhibits the oscillating property in the interface crack-tip generalized stress field and the other that does not. Fifteen (15) pair-combinations of six (6) piezoelectric materials PZT-4, PZT-5H, PZT-6B, PZT-7A, P-7, and BaTiO₃, which are commonly used in practice, are numerically analyzed in this study, and the results backup the above theoretical conclusions. Moreover, the associated eigenvectors for such material systems (with either =0 or κ=0) are also obtained numerically, and the result show that there still exist four linear independent associate eigenvectors for each bimaterial.     

Line-integral representations for the elastic displacements,stresses and interaction energy of arbitrary dislocation loops in transversely isotropic bimaterials

The elastic displacements, stresses and interaction energy of arbitrarily shaped dislocation loops with general Burgers vectors in transversely isotropic bimaterials (i.e. joined half-spaces) are expressed in terms of simple line integrals for the first time. These expressions are very similar to their isotropic full-space counterparts in the literature and can be easily incorporated into three-dimensional (3D) dislocation dynamics (DD) simulations for hexagonal crystals with interfaces/surfaces. All possible degenerate cases, e.g. isotropic bimaterials and isotropic half-space, are considered in detail. The singularities intrinsic to the classical continuum theory of dislocations are removed by spreading the Burgers vector anisotropically around every point on the dislocation line according to three particular spreading functions. This non-singular treatment guarantees the equivalence among different versions of the energy formulae and their consistency with the stress formula presented in this paper. Several numerical examples are provided as verification of the derived dislocation solutions, which further show significant influence of material anisotropy and bimaterial interface on the elastic fields and interaction energy of dislocation loops.     

Computing lower and upper bounds on stress intensity factors in bimaterials

This paper presents a finite element approach for finding complementary bounds of stress intensity factors (SIFs) in bimaterials. The SIF is formulated as an explicit computable linear function of displacements by means of the two-point extrapolation method. An appropriate and computable form of the SIF plays a crucial role in the dual problem involved in the computing procedure of both lower and upper bounds. In our discussions, computable forms of stress intensity factors, K₀ and K_r, are derived, which are related to the energy release rate, and the ratio of the open mode and shear mode SIFs, respectively. Based on a posteriori finite element error estimation, a bounding procedure is used to compute the bounds on the two stress intensity factors. Finally, bounds on the SIFs in a bimaterial interface crack problem are provided to verify the procedure.     

Arbitrarily oriented crack near interface in piezoelectric bimaterials

Arbitrarily oriented crack near interface in piezoelectric bimaterials is considered. After deriving the fundamental solution for an edge dislocation near the interface, the present problem can be expressed as a system of singular integral equations by modeling the crack as continuously distributed edge dislocations. In the paper, the dislocations are described by a density function defined on the crack line. By solving the singular integral equations numerically, the dislocation density function is determined. Then, the stress intensity factors (SIFs) and the electric displacement intensity factor (EDIF) at the crack tips are evaluated. Subsequently, the influences of the interface on crack tip SIFs, EDIF, and the mechanical strain energy release rate (MSERR) are investigated. The J-integral analysis in piezoelectric bimaterals is also performed. It is found that the path-independent of J₁-integral and the path-dependent of J₂-integral found in no-piezoelectric bimaterials are still valid in piezoelectric bimaterials.     

Theoretical and Experimental Investigations of the Stress Concentration in Nonthin Cylindrical Shells with Rectangular Openings

The problem on the stress state of cylindrical shells with rectangular openings under axial compression is considered based on the theory of shells of average thickness. For a shell with two diametrical openings, numerical calculations are carried out by the finite-difference method for anisotropic and isotropic materials. The anisotropy of the material and the sizes of the openings are shown to affect the displacements and stresses on the opening periphery. An experiment is implemented for a cylindrical shell with rectangular openings under axial compression. The experimental and numerical data are compared     

2维各向异性弹性力学中的Stroh公式和某些不变量

本文简要综述了各向异性弹性体2维变形的Stroh公式。该公式数学上优美而应用上有力。它的最近进展使我们将已有的复数形式的解转变成实数形式,进而可去解决一些迄今尚未解出的问题。结果,发现了许多有趣的而在物理上难以预料的现象。例如,对于界面包含x_3轴的各向异性双介质体,如果线力和线位错都沿x_3轴,那末包含x_3轴的任意径向平面上的面拉力将不随径向平面的选择而变化。而且,在任何固定点上的应力与界面的方位无关。要揭示的其他一些不变性质是:双介质材料中的虚拟力,界面裂纹上的应力奇异,界面裂纹表面上的振荡,等等;另外,在具有一个有限裂纹,或一个椭圆孔,或一个刚性椭圆核的无穷大均匀介质中的某些应力和位移也是不变量。这些不变量的大部分尚未有它们的物理解释。在实数形式的解中,Barnett-Lothe张量常常出现。在绕x_3轴转动时,这些张量的分量的某些组合是不变的。这些组合将在本文中推导并列出。最后我们给出了某些各向异性弹性材料,它们的Barnett-Lothe张量是不变量,即这些张量的所有分量在绕x_3轴转动时是不变的。     

Crack deflection at an interface between dissimilar elastic materials: Role of residual stresses

A crack intersecting an interface between two dissimilar materials may advance by either penetrating through the interface or deflecting into the interface. The competition between deflection and penetration can be assessed by comparison of two ratios: (i) the ratio of the energy release rates for interface cracking and crack penetration; and (ii) the ratio of interface to material fracture energies. Residual stresses caused by thermal expansion misfit can influence the energy release rates of both the deflected and penetrating crack. This paper analyses the role of residual stresses. The results reveal that expansion misfit can be profoundly important in systems with planar interfaces (such as layered materials, thin film structures, etc.), but generally can be expected to be of little significance in fiber composites. This paper corrects an earlier result for the ratio of the energy release rate for the doubly deflected crack to that for the penetrating crack in the absence of residual stress.     

SINGULARITIES AND EXPLICIT SOLUTIONS OF COLLINEAR CRACK PROBLEM IN ANISOTROPIC MEDIA

The singularities of collinear cracks both in anisotropic single mediumand at the interface of anisotropie bimaterials are studied by combining Stroh formalismand the analytic function method.The formulae for calculating the field potential andstress intensity factor(SIF)are obtained.It is found that the field potentials areexplicitly related to material matrix L and the in-plane and anti-plane fields can beseparately calculated when orthotropic bimaterials are considered.