The computation of stress intensity factors in dissimilar materials

A reciprocal work contour integral method for calculating stress intensity factors is extended to treat the problem of two bonded dissimilar materials containing a crack along the bond. The method is based on Betti's Reciprocal work theorem from which the singular stress intensities at the crack tip may be evaluated in terms of an integral involving tractions and displacements on a contour remote from the crack tip.     

概率断裂分析的随机边界元法

应用随机边界元法分析材料弹性常数的随机性和裂纹面随机性对应力强度因子的影响。文中首先简介了随机边界元法，给出了具有随机材料或几何参数的弹性体的边界位移或面力的协方差，进而给出了材料参数和裂纹面随机时应力强度因子均值和方差的计算公式。算例中详细讨论了杨氏模量、泊松比及裂纹面的随机性对应力强度因子的影响。     

Exact solution for mixed boundary value problems at anisotropic piezoelectric bimaterial interface and unification of various interface defects

The special mixed boundary value problem in which a debonded conducting rigid line inclusion is embedded at the interface of two piezoelectric half planes is solved analytically by employing the 8-D Stroh formalism. Different from existing interface insulating crack model and interface conducting rigid line inclusion model, the presently analyzed model is based on the assumption that all of the physical quantities, i.e., tractions, displacements, normal component of electric displacements and electric potential, are discontinuous across the interface defect. Explicit solutions for stress singularities at the tips of debonded conducting rigid line inclusion are obtained. Closed form solutions for the distribution of tractions on the interface, surface opening displacements and jump in electric potential on the debonded inclusion are also obtained, in addition real form solutions for these physical quantities are derived. Various forms of interface defect problems encountered in practice are solved within a unified framework and the stress singularities induced by those interface defects are discussed in detail. Particularly, we find that the analysis of interface cracks between the embedded electrode layer and piezoelectric ceramics can also be carried out within the unified framework.     

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.     

The cantilever beam under tension,bending or flexure at infinity

A novel technique, the method of projection, is applied to the plane strain problems of determining the tractions, and stress intensity factors, at the fixed end of a cantilever beam under tension, bending or flexure at infinity. The method represents a useful alternative to the integral equation method of Erdogan, Gupta and Cook, and possesses certain advantages; in particular it is much easier to extend the present method to the more difficult dynamics case. An unusual feature of the method is that the required tractions are expanded as a series whose terms have the natural role of displacements rather than stresses.     

V形切口应力强度因子的一种边界元分析方法

将V形切口结构分成围绕切口尖端的小扇形和剩余结构两部分. 尖端处扇形域应力场表示成关于尖端距离$\rho$的渐近级数展开式,从线弹性理论方程推导出了一组分析平面V形切口奇异性的常微分方程特征值问题,通过求解特征方程,得到前若干个奇性指数和相应的特征向量. 再将切口尖端的位移和应力表示为有限个奇性阶和特征向量的组合. 然后用边界元法分析挖去小扇形后的剩余结构. 将位移和应力的线性组合与边界积分方程联立,求解获得切口根部区域的应力场、应力幅值系数和整体结构的位移和应力. 从而准确计算出平面V形切口的奇异应力场和应力强度因子.      

Transient dynamic analysis of interface cracks in anisotropic bimaterials by the scaled boundary finite-element method

The transient response of finite bimaterial plates with interface cracks is analyzed directly in the time domain by using the scaled boundary finite-element method. A bimaterial plate is divided into a few subdomains. Only the boundaries of the subdomains are discretized with line elements leading to great flexibility in mesh generation. The displacement and stress fields are expressed as a series solution which separates the singular stress term from other high-order terms. The oscillatory stress singularity in the radial direction emanating from the scaling center is represented analytically. The complex dynamic stress intensity factors are evaluated directly from either the stresses or the crack opening displacements of the singular stress term. Numerical examples of cracked anisotropic bimaterial plates are presented to verify the accuracy of the present technique and to provide additions to the very limited number of reference solutions in the literature.     

Dynamic interaction between a sub-interface crack and the interface in a bi-material: time-domain BEM analysis

Summary Transient response of a sub-interface crack in a bi-material is studied with emphasis on the dynamic interaction between the crack and the interface, by combining the traditional time-domain displacement boundary element method (BEM) and the non-hypersingular traction BEM. Computations are performed for an unbounded bi-material with a crack subjected to impact tensile loading on its faces or incident impact waves and a bounded rectangular bi-material plate under remote impact tensile loading. Numerical results of the dynamic stress intensity factors (DSIFs) and dynamic interface tractions are presented for various material combinations and crack locations. It is shown that pronounced increases in DSIFs and the interface tractions may be caused in some cases because of the dynamic interaction between the crack and the interface.This work was initialized during the second author's stay at Institute of Mechanics, TU Darmstadt, Germany under the support of the Alexander von Humboldt Foundation. Discussion on the BEM formulation with Dr. Seelig is gratefully acknowledged. The first two authors are also grateful for the partial support by the China National Natural Science Foundation under Grant No. 10025211 and the NJTU Scientific Paper Fund (PD195).     

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

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

Fracture mechanics analysis of curved stiffened panels using BEM

A dual boundary element method is developed for a analysis of reinforced cracked shallow shells. Boundary integral equations are derived from the Betti’s reciprocal theorem for a cracked shallow shell with transverse frames and longitudinal stiffeners. The effect of frames and stiffeners are treated as a distribution of line body forces. The radial basis function is used to transform domain integrals to boundary integrals. Stress intensity factors are evaluated from crack opening displacements. The effect of curvature on the stress intensity factors is illustrated by numerical examples. Three examples are presented to demonstrate the accuracy of this method compared with solutions obtained using the finite element method.     

Periodic array of traction-free interface cracks subjected to far-field uniform heat flow

A periodic array of interface cracks is subjected to a uniform heat flow in the far field. The crack opening displacements and complex stress intensity factors are determined, using analytic function theory, for the case that the upper half-space is less distortive than the lower half-space.     

Coulomb traction on a penny-shaped crack in a three dimensional piezoelectric body

The axisymmetric problem of a penny-shaped crack embedded in an infinite three-dimensional (3D) piezoelectric body is considered. A general formulation of Coulomb traction on the crack surfaces can be obtained based on thermodynamical considerations of electromechanical systems. Three-dimensional electroelastic solutions are derived by the classical complex potential theory when Coulomb traction is taken into account and the poling direction of piezoelectric body is perpendicular to the crack surfaces. Numerical results show that the magnitude of Coulomb tractions can be large, especially when a large electric field in connection with a small mechanical load is applied. Unlike the traditional traction-free crack model, Coulomb tractions induced by an applied electric field influence the Mode I stress intensity factor for a penny-shaped crack in 3D piezoelectric body. Moreover, compared to the current model, the traditional traction-free crack model always overestimates the effect of the applied electric load on the field intensity factors and energy release rates, which has consequences for 3D piezoelectric fracture mechanics.     

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.     

含界面裂纹的不可压缩双金属胶接件的解析变分－守恒积分解法

本文根据平面问题的复变函数理论推导了含界面裂纹双金属胶接件满足微分方程、开裂界面边界条件与未开裂界面连续条件的应力与位移本征函数展开式，并建立了不可压缩双金属界面裂纹的复合型守恒积分及其与应力强度因子之关系，进而利用分区广义变分原理满足其余边界条件确定包含应力强度因子在内的展开式系数，得到守恒积分并求出应力强度因子．数值计算表明，沿不同回路的在恒积分具有很好的守恒性而且由这两种方法所得应力强度因子具有很好的一致性．     

Penny-shaped crack at the interface between elastic half-spaces under the action of a shear wave

The three-dimensional problem of elasticity for a bimaterial body with a penny-shaped crack at the interface under the action of a normal harmonic shear wave is solved by the boundary-element method. The distribution of displacements of crack faces and tractions and displacements at the interface is analyzed     

黏弹性界面断裂分析的增量“加料”有限元法

提出了一种适用于黏弹性界面裂纹问题的增量“加料” 有限元方法. 利用弹性界面裂纹尖端位移场的解答,通过对应原理和拉普拉斯逆变换近似方法,得到了黏弹性界面裂纹的尖端位移场. 用该位移场构造了黏弹性界面裂纹“加料” 单元和过渡单元位移模式,推导了增量“加料” 有限元方程,求解有限元方程可获得应力强度因子和应变能释放率等断裂参量. 建立了典型黏弹性界面裂纹平面问题“加料” 有限元模型,计算结果表明,对于弹性/黏弹性界面裂纹和黏弹性/黏弹性界面裂纹,该方法都能得到相当精确地断裂参量,并能很好地反映蠕变和松弛特性,可推广应用于黏弹性界面断裂问题的计算分析.      

A Self-Similar Crack Expansion method for three-dimensional cracks in an infinite or semi-infinite medium

The Self-Similar Crack Expansion (SSCE) method is used to calculate stress intensity factors for three-dimensional cracks in an infinite medium or semi-infinite medium by the boundary integral element technique, whereby, the stress intensity factors at crack tips are determined by calculating the crack-opening displacements over the crack surface. For elements on the crack surface, regular integrals and singular integrals are precisely evaluated based on closed form expressions, which improves the accuracy. Examples show that this method yields very accurate results for stress intensity factors of penny-shaped cracks and elliptical cracks in the full space, with errors of less than 1% as compared with analytical solutions. The stress intensity factors of subsurface cracks are in good agreement with other analytical solutions.     

Boundary Value Problems of Plane Elasticity Involving Orientations of Displacements and Tractions

This paper presents unconventional formulations of boundary problems of plane elasticity formulated in terms of orientations of tractions and displacements on a closed contour separating internal and external domains as the boundary conditions. These are combined with the conditions of continuity of tractions or displacements across the boundary. Therefore the magnitudes of neither tractions nor displacements are assumed on the contour. Four boundary value problems for both external and internal domains are investigated by analyzing the solvability of the corresponding singular integral equations. It is shown that all considered problems can have non-unique solutions expressed as linear combinations of particular solutions and, hence, contain free arbitrary parameters, the number of which is finite and determined by the contour orientations of tractions and/or displacements     

界面裂纹问题中的权函数方法

本文将Ｐａｒｉｓ等确定均匀材料中裂纹尖端应力强度因子的权函数方法推广应用到界面裂纹问题，给出了界面裂纹尖端附近或无限大体半无限界面裂纹问题的权函数的显式表达式。利用此权函数表达式可以很简便地求解界面裂纹尖端附近一些外来作用引起的应力强度因子，比如任意分布力、相变应变、位错和热等。作为一个算例，本文计算了界面一侧一个刃型位错引起的应力强度因子。     

INTERACTION OF A DISLOCATION DIPOLE WITH A MODE III INTERFACE CRACK

Interaction between a screw dislocation dipole and a mode III interface crack is investigated. By using the complex variable method, the closed form solutions for complex potentials are obtained when a screw dislocation dipole lies inside a medium. The stress fields and the stress intensity factors at the tip of the interface crack produced by the screw dislocation dipole are given. The influence of the orientation, the dipole arm and the location of the screw dislocation dipole as well as the material mismatch on the stress intensity factors is discussed. The image force and the image torque acting on the screw dislocation dipole center are also calculated. The mechanical equilibrium position of the screw dislocation dipole is examined for various material property combinations and crack geometries. The results indicate that the shielding or anti-shielding effect on the stress intensity factor increases abruptly when the dislocation dipole approaches the tip of the crack. Additionally, the disturbation of the interface crack on the motion of the dislocation dipole is also significant.