Analysis of a three-dimensional dissimilar material joint with one real singularity using a conservative integral

In the present study, a conservative integral based on the Betti reciprocal principle is formulated to determine the intensity of singularity at a vertex of the interface in three-dimensional dissimilar material joints with one real singularity. Eigenanalysis formulated using a three-dimensional finite element method (FEM) is used to calculate the order of stress singularity, angular functions of displacements and stresses. Models with various element sizes and various integral areas are used to investigate the effect of the integration area on the accuracy of the results. The results are compared with those obtained from the boundary element method (BEM) using a curve-fitting technique to calculate the intensity of singularity. In addition, models of various lengths and various material combinations are used to investigate the stress singularity characteristics in three-dimensional dissimilar material joints. The results of the present study indicate that the conservative integral can be used to determine the intensity of singularity in three-dimensional bi-material joints. The accuracy of the results can be improved by mesh refinement. Finally, the relationships among the intensity of singularity, the order of stress singularity and the model geometry are discussed.     

三维双材料结构的应力奇异性分析

应用有限单元法子模型技术，对具有不同界面角的三维双材料结构的应力奇异性进行了分析。结果表明，应用子模型技术估算三维双材料结构的应力奇异性指数是有效的。然后分析了界面端线和界面端点处附近奇异性指数，得到了一些重要而有趣的结果。最后对消除三维双材料结构应力奇异性的几何条件进行了讨论。     

多重应力奇异性及其强度系数的数值分析方法

以具有两个应力奇异性次数的平面问题为例,提出了一种利用普通的数值分析结果确定奇异点附近多重应力奇异性的各阶次数以及相应的应力强度系数的数值分析方法,计算实例表明,本方法可以精确地求得各阶应力奇异性的次数,并且可以很方便地应用外插法确定出对应的应力强度系数。     

结合材料界面端的三维应力奇异性

本文利用特殊有限元方法，开发了一个用来求解结合材料界面端三维应力奇异性问题的数值分析程序。该方法只需对界面端的角度方向进行离散即可求得应力奇异性。结合材料的应力奇异性取决于两种材料的材料常数和界面端形状。选用三个材料参数作为变量，用来研究结合材料三维应力奇异性随材料常数的变化规律。文中计算了几种重要而且常见的情况，并以此为基础建立了数据库。同时，还分析了应力奇异性随界面端形状的变化规律，并得到了应力函数的分布图。     

单纤维压出实验模型的力学分析

单纤维压出实验是一种细观实验方法,用于测量复合材料的界面强度。本文采用子域法处理多域问题,利用轴对称边界元程序,对单纤维压出实验模型的应力传递进行了分析,并与常用的剪切滞后理论的分析结果进行了比较,发现在大部分纤维埋置区域,两种分析吻合很好,在纤维与基体的界面端部,边界元的分析结果显示界面应力在该区域有奇异现象,而剪滞理论则无法反映应力奇异现象。应力奇异现象的存在使得我们对该实验的界面强度判据需要     

The Analysis of Dynamic Stress Intensity Factor for Semi-Circular Surface Crack Using Time-Domain BEM Formulation

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｒａｎｓｉｅｎｔａｎａｌｙｓｉｓｆｏｒ 3＿Ｄｅｌａｓｔｏｄｙｎａｍｉｃｃｒａｃｋｐｒｏｂｌｅｍｓｉｓａｎｉｍｐｏｒｔａｎｔｃｏｎｔｅｎｔｉｎｆｒａｃｔｕｒｅｍｅｃｈａｎｉｃｓ.Ｂｅｃａｕｓｅｏｆｔｈｅｄｉｆｆｉｃｕｌｔｉｅｓｏｆｔｈｅｐｒｏｂｌｅｍｓ,ｔｈｅｃａｌｃｕｌａｔｉｏｎｏｆＤＳＩＦｆｏｒ 3＿Ｄｃｒａｃｋｉｓａｐｒｏｂｌｅｍｕｎｓｏｌｖｅｄｃｏｍｐｌｅｔｅｌｙｕｐｔｏｔｈｅｐｒｅｓｅｎｔ.Ｃｈｅｎ[1]ｆｉｒｓｔｌｙｕｓｅｄｔｈｅｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅ…     

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).     

Stress singularity analysis around the singular point on the stress singularity line in three-dimensional joints

The characteristics of the stress fields around a singular point on the stress singularity line of dissimilar materials in three-dimensional joints are investigated using BEM. Contour for the order of stress singularity around the point is mapped on Dundurs’ parameters plane using eigen value analysis by FEM. The results in 3D joints are compared with those in 2D joints having the same cross section and material combination. The order of stress singularity around the singular point on the stress singularity line in 3D joints is almost identical with that in 2D joints in the singularity region. However, the zero boundary of singularity in 3D joints is slightly different from that in 2D joints. Furthermore, the multiple root of p = 1 exists in the eigen value analysis by FEM. Therefore, logarithmic singularity possibly occurs around the singular point on the stress singularity line. Then, the stress distributions around this point are expressed by the combination of the r^λ term and logarithmic singularity terms. Finally, the characteristics of the stress intensity factors of the r^λ term and logarithmic singularity terms around the singular points are investigated.     

黏弹性材料界面裂纹应力场奇性分析

研究两半无限大黏弹性体间Griffith界面裂纹在简谐载荷作用下裂纹尖端动应力场的奇异特性.通过引入裂纹张开位移和裂纹位错密度函数,相应的混合边值问题归结为一组耦合的奇异积分方程.渐近分析表明裂尖动应力场的奇异特征完全包含在奇异积分方程的基本解中.通过对基本解的深入分析发现黏弹性材料界面裂纹裂尖动应力场具有与材料参数和外载荷频率相关的振荡奇异特性.以标准线性固体黏弹材料为例讨论了材料参数和载荷频率对奇性指数和振荡指数的影响.     

Elastic-viscoplastic field at mixed-mode interface crack-tip under compression and shear

For a compression-shear mixed mode interface crack, it is difficult to solve the stress and strain fields considering the material viscosity, the crack-tip singularity, the frictional effect, and the mixed loading level. In this paper, a mechanical model of the dynamic propagation interface crack for the compression-shear mixed mode is proposed using an elastic-viscoplastic constitutive model. The governing equations of propagation crack interface at the crack-tip are given. The numerical analysis is performed for the interface crack of the compression-shear mixed mode by introducing a displacement function and some boundary conditions. The distributed regularities of stress field of the interface crack-tip are discussed with several special parameters. The final results show that the viscosity effect and the frictional contact effect on the crack surface and the mixed-load parameter are important factors in studying the mixed mode interface crack- tip fields. These fields are controlled by the viscosity coefficient, the Mach number, and the singularity exponent.     

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the interface crack tip for mode Ⅱ of double dissimilar orthotropic composite materials are studied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized bi-harmonic equations,the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions,a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about himaterial engineering parameters. According to the uniqueness theorem of limit,both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same,the stress singularity exponents,stress intensity factors and stresses for mode Ⅱ crack of the orthotropic single material are obtained.     

移动接触下求解二维闭合裂纹的双重迭代边界元法

使用子域边界元法对受移动接触弹性体作用下的二维闭合裂纹问题进行了数值计算。由于两弹性体的接触界面和裂纹表面的接触范围的大小和接触状态事先是未知的 ,对此 ,在两个接触表面同时采用迭代的方法进行了求解。在裂纹的每个裂尖上都采用了四分之一的奇异单元以保证裂尖位移场和应力场奇异性的满足。用我们编制的二维裂纹问题程序对一些中心裂纹问题进行了计算 ,计算结果与经典断裂力学的理论值比较吻合。在无摩擦的条件下 ,对一些具有不同角度且受移动接触弹性体作用下的闭合裂纹问题进行了数值计算 ,得到了一些耦合作用下的应力强度因子的计算结果     

Corner stress singularities in an FGM thin plate

Describing the behaviors of stress singularities correctly is essential for obtaining accurate numerical solutions of complicated problems with stress singularities. This analysis derives asymptotic solutions for functionally graded material (FGM) thin plates with geometrically induced stress singularities. The classical thin plate theory is used to establish the equilibrium equations for FGM thin plates. It is assumed that the Young’s modulus varies along the thickness and Poisson’s ratio is constant. The eigenfunction expansion method is employed to the equilibrium equations in terms of displacement components for an asymptotic analysis in the vicinity of a sharp corner. The characteristic equations for determining the stress singularity order at the corner vertex and the corresponding corner functions are explicitly given for different combinations of boundary conditions along the radial edges forming the sharp corner. The non-homogeneous elasticity properties are present only in the characteristic equations corresponding to boundary conditions involving simple support. Finally, the effects of material non-homogeneity following a power law on the stress singularity orders are thoroughly examined by showing the minimum real values of the roots of the characteristic equations varying with the material properties and vertex angle.     

Application of the reciprocity theorem to scattering of surface waves by a cavity

Scattering of surface waves by a cylindrical cavity at the surface of a homogenous, isotropic, linearly elastic half-space is analyzed in this paper. In the usual manner, the scattered field is shown to be equivalent to the radiation from a distribution of tractions, obtained from the incident wave on the surface of the cavity. For the approximation used in this paper, these tractions are shifted to tractions applied to the projection of the cavity on the surface of the half-space. The radiation of surface waves from a normal and a tangential line load, recently determined by the use of the reciprocity theorem, is employed to obtain the field scattered by the cavity from the superposition of displacements due to the distributed surface tractions. The vertical displacement at some distance from the cavity is compared with the solution of the scattering problem obtained by the boundary element method (BEM) for various depths and widths of the cavity. Comparisons between the analytical and BEM results are graphically displayed. The limitations of the approximate approach are discussed based on the comparisons with the BEM results.     

The collinear crack problem for an orthotropic functionally graded coating–substrate structure

A theoretical treatment of antiplane crack problem of two collinear cracks on the two sides of and perpendicular to the interface between a functionally graded orthotropic strip bonded to an orthotropic homogeneous substrate is put forward. Various internal cracks and crack terminating at the interface and crack crossing the interface configurations are investigated, respectively. The problem is formulated in terms of a singular integral equation with the crack face displacement as the unknown variable. The asymptotic stress field near the tip of a crack crossing the interface is examined, and it is shown that, unlike the corresponding stress field in piecewise homogeneous materials, in this case, the “kink” in material property at the interface does not introduce any singularity. Numerical calculations are carried out, and the influences of the orthotropy and nonhomogeneous parameters and crack interactions on the mode III stress intensity factors are investigated.     

Interface crack problems for mode II of double dissimilar orthotropic composite materials

The fracture problems near the similar orthotropic composite materials are interface crack tip for mode Ⅱ of double disstudied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized hi-harmonic equations, the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions, a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about bimaterial engineering parameters. According to the uniqueness theorem of limit, both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same, the stress singularity exponents, stress intensity factors and stresses for mode II crack of the orthotropic single material are obtained.     

纤维增强复合材料圆柱型界面裂纹分析

以裂纹面上的位错函数为未知量将圆柱型界面裂纹问题化成一组奇异积分方程的求解问题．应用Ｍｕｓｋｈｅｌｉｓｈｖｉｌｉ的奇异积分方程理论，分析了圆柱型界面裂纹尖端应力场．针对裂纹尖端分别存在和不存在接触区两种情况，确定了裂纹尖端应力场的奇异性．利用数值方法计算了圆柱型界面裂纹尖端接触区尺寸对剪应力强度因子的影响．     

An experimental study of stress singularities at a sharp corner in a contact problem

The photoelastic method was used to investigate the nature of the local stress field at a sharp corner of a wedge that was compressed against a larger body. Planar wedge specimens made of photoelastic material were compressed against a half plane (larger body) of identical material at various load levels. Several wedge angles were studied. The nature of the singular stress field postulated by linear elastic analysis was verified and the strength of the singularity was obtained by plotting the variation of fringe order as a function of radial distance from the sharp corner on a logarithmic scale. The experimental results were found to be in good agreement with the theoretical predictions. The effect of interface friction and the effect of rounding off the sharp edge are brietly discussed.     

有限尺寸下双材料界面附近裂纹位置对裂纹尖端的影响

随着复合材料的应用和发展,不同材料组成的界面结构越来越受到人们的重视。界面层两侧材料的性能相异会引起材料界面端奇异性,同时界面和界面附近存在裂纹会引起裂尖处的应力奇异性。因此双材料界面附近的力学分析是比较复杂的。本文建立双材料直角界面模型,在材料界面附近预设初始裂纹,计算了有限材料尺寸对界面应力场及其附近裂纹应力强度因子的影响。运用弹性力学中的 Goursat 公式求得直角界面端在有限尺寸下的应力场以及其应力强度系数。通过叠加原理和格林函数法进一步得到在直角界面端附近的裂纹尖端应力强度因子。计算结果表明,在适当范围内改变材料内裂纹与界面之间的距离,界面附近裂纹尖端的应力强度因子随着裂纹与界面距离的增加而减少,并且逐渐趋于稳定。分析结果可以为预测双材料结构复合材料界面失效位置提供参考。     

Singular stress field near the edge of interface of bonded dissimilar materials with an interlayer

For bonded dissimilar materials, the free-edge stress singularity usually prevails near the intersection of the free-surface and the interface. When two materials are bonded by using an adhesive, an interlayer develops between the two bonded materials. When a ceramic and a metal are bonded, the residual stress develops because of difference in the coefficient of thermal expansion. An interlayer may be inserted between the two materials to defuse the residual stress. Stress field near the intersection of the interface and free-surface in the presence of the interlayer is then very important for evaluating the strength of bonded dissimilar materials.In this study, stress distributions on the interface of bonded dissimilar materials with an interlayer were calculated by using the boundary element method to investigate the effect of the interlayer on the stress distribution. The relation between the free-edge singular stress fields of bonded dissimilar materials with and without an interlayer was investigated numerically. It was found that the influence of the interlayer on the stress distributions was confined within a small area of the order of interlayer thickness around the intersection of the interface and the free-surface when the interlayer was very thin. The stress distribution near the intersection of the interface and the free-surface was controlled by the free-edge stress singularity of the bonded dissimilar materials without the interlayer. In this case, the interlayer can be called free-edge singularity-controlled interlayer. If a stress distribution on the interface is known for one thickness of an interlayer h, stress distributions on the interface for other values of h can be estimated.