Dynamic anti-plane analysis for two symmetrically interfacial cracks near circular cavity in piezoelectric bi-materials

The present paper is exposed theoretically to the influence on the dynamic stress intensity factor （DSIF） in the piezoelectric bi-materials model with two symmet- rically permeable interracial cracks near the edges of a circular cavity, subjected to the dynamic incident anti-plane shearing wave （SH-wave）. An available theoretical method to dynamic analysis in the related research field is provided. The formulations are based on Green＇s function method. The DSIFs at the inner and outer tips of the left crack are obtained by solving the boundary value problems with the conjunction and crack- simulation technique. The numerical results are obtained by the FORTRAN language program and plotted to show the influence of the variations of the physical parameters, the structural geometry, and the wave frequencies of incident wave on the dimensionless DSIFs. Comparisons with previous work and between the inner and outer tips are con- cluded.     

Special approach for the determination of fracture mechanical parameters at an interface crack tip

An interface crack of finite length is considered between two semi-infinite planes with an artificial contact zone at one of the two crack tips. A transcendental equation and certain simple asymptotic formulas are established for the real contact zone (in the Comninou-Dundurs sense) in terms of the stress intensity factors (SIFs) of the considered model. In these terms analytical expressions are also provided for the energy release rate and for the SIF of the classical interface crack model with an oscillating singularity at the crack tip. The appropriate length of the artifical contact zone is shown to be attainable on the basis of the analysis of the stresses at the crack tip. The use of the proposed model is suggested for integrity assessment of inhomogeneous structural elements of composites containing interface cracks. Received 26 March 1997; accepted for publication 12 September 1997     

节理岩体爆生裂纹扩展动态焦散线模型实验研究

应用动态焦散线测试系统,模拟含节理岩体断裂爆破过程,进行了PMMA模型透射式动态焦散线实验,着重研究了爆炸初始裂纹与节理面夹角不同的情况下,裂纹尖端动态强度因子的变化规律、裂纹穿过节理面的扩展规律、以及炮孔与节理距离不同时裂纹穿过节理扩展的规律。实验结果分析表明,爆生裂纹穿过节理面时,裂纹尖端的动态强度因子和裂纹扩展速度显著下降,穿过节理面后,强度因子又增强;裂纹穿过节理面时,裂纹会沿节理面偏离一段距离后沿原方向继续扩展。实验结果还表明,炮孔与节理间距适宜的情况下,裂纹才会穿过节理继续扩展,间距太小和过大都不利于裂纹的扩展。     

Dynamic self-similar debonding of interface at very high velocity: Anti-plane case

This paper analyzes the anti-plane problem of dynamic self-similar debonding of interface at very high velocity. The debonding is modeled as an interface crack propagating self-similarly from zero-length. The extending speed is assumed to be transonic or supersonic. We first consider the dynamic debonding under moving concentrated loads. The moving dislocation model of self-similar propagation of an interface crack is used to formulate the problem to a singular integral equation which is solved analytically. The singularity of stresses near the crack tip is discussed and the dynamic stress intensity factors are presented. Finally the solution of dynamic debonding underx ²-type loads is obtained by using the superposition method.     

功能梯度双材料弱/微间断界面的冲击断裂分析

提出强间断、弱间断、微间断和全连续界面的概念与分类,建立功能梯度弹性双材料弱间断 界面冲击断裂问题的力学模型,采用积分变换法推导问题的Cauchy奇异积分方程,并用配 点法求得数值解. 分析表明,弱/微间断性对于FGMs界面裂纹应力强度因子有着重要影响, 而且微间断性是优于弱间断性的一种界面力学性能连接关系. 以FGMs界面某一侧 的力学性能函数在界面处的Taylor展开式的低阶项作为界面另一侧的力学性能函数,便可 以使FGMs界面成为``微间断'界面. 界面的一阶微间断对应力强度因子的减小作用较为明 显,而高阶(二阶及以上)微间断对应力强度因子的影响较小. 减小界面的弱间断程度或使 FGMs界面具备``微间断性',都将利于提高功能梯度双材料界面抗冲击断裂能力,在一定 程度上达到界面增韧的目的.     

Mechanical modeling and transient anti-plane fracture analysis for the graded inter-diffusion regions in a bonded structure

An inter-diffusion interface model (IDIM) is put forward for a bonded structure. Laplace and Fourier integral transforms are applied to reduce the transient anti-plane fracture problem of the structure as a Cauchy singular integral equation. Lobatto-Chebyshev collocation method and Laplace numerical inversion transform are employed to evaluate transient stress intensity factors (TSIFs). The effects of geometrical and physical parameters on TSIFs are studied. Specially discussed are the influences of the weak/micro-discontinuity of the interfaces. Comparison between IDIM and the graded interlayer model indicates that if the inter-diffusion between the two original materials is prominent, the former should be applied instead of the latter in fracture analyses of bonded structures.     

Analysis of stress intensity factor in orthotropic bi-material mixed interface crack

Adopting the complex function approach, the paper studies the stress intensity factor in orthotropic bi-material interface cracks under mixed loads. With consideration of the boundary conditions, a new stress function is introduced to transform the problem of bi-material interface crack into a boundary value problem of partial differential equations. Two sets of non-homogeneous linear equations with 16 unknowns are constructed. By solving the equations, the expressions for the real bi-material elastic constant εt and the real stress singularity exponents λt are obtained with the bi-material engineering parameters satisfying certain conditions. By the uniqueness theorem of limit,undetermined coefficients are determined, and thus the bi-material stress intensity factor in mixed cracks is obtained. The bi-material stress intensity factor characterizes features of mixed cracks. When orthotropic bi-materials are of the same material, the degenerate solution to the stress intensity factor in mixed bi-material interface cracks is in complete agreement with the present classic conclusion. The relationship between the bi-material stress intensity factor and the ratio of bi-material shear modulus and the relationship between the bi-material stress intensity factor and the ratio of bi-material Young's modulus are given in the numerical analysis.     

The dynamic stress intensity factor analysis of adhesively bonded material interface crack with damage under shear loading

This paper studies the dynamic stress intensity factor (DSIF) at the interface in an adhesive joint under shear loading. Material damage is considered. By introducing the dislocation density function and using the integral transform, the problem is reduced to algebraic equations and can be solved with the collocation dots method in the Laplace domain. Time response of DSIF is calculated with the inverse Laplace integral transform. The results show that the mode Ⅱ DSIF increases with the shear relaxation parameter, shear module and Poisson ratio, while decreases with the swell relaxation parameter. Damage shielding only occurs at the initial stage of crack propagation. The singular index of crack tip is -0.5 and independent on the material parameters, damage conditions of materials, and time. The oscillatory index is controlled by viscoelastic material parameters.     

Scattering of love waves by an interface crack between a piezoelectric layer and an elastic substrate

The scattering of Love waves by an interface crack between a piezoelectric layer and an elastic substrate is investigated by using the integral transform and singular integral equation techniques. The dynamic stress intensity factors of the left and the right crack tips are determined. It is found from numerical calculation that the dynamic response of the system depends significantly on the crack configuration, the material combination and the propagating direction of the incident wave. It is expected that specifying an appropriate material combination may retard the growth of the crack for a certain crack configuration. Project supported by the National Natural Science Foundation of China (No. 19891180), the Fundamental Research Foundation of Tsinghua University (JZ 2000.007) and the Fund of the Education Ministry of China.     

A study on stress intensity factors and singular stress fields of three-dimensional interface crack

By using the finite-part integral concepts and limit technique, the hypersingular integrodifferential equations of three-dimensional (3D) planar interface crack were obtained; then the dominant-part analysis of 2D hypersingular integral was further used to investigate the stress fields near the crack front theoretically, and the accurate formulae were obtained for the singular stress fields and the complex stress intensity factors. After that, a numerical method is proposed to solve the hypersingular integrodifferential equations of 3D planar interface crack, and the problem of elliptical planar crack is then considered to show the application of the method. The numerical results obtained are satisfactory. Project supported by the Foundation of Solid Mechanics Open Research Laboratory of State Education Commission at Tongji University and the National Natural Science Foundation.     

Effects of the weak/micro-discontinuity of interface on the fracture behavior of a functionally graded coating with an inclined crack

The mechanical model was established for the anti-plane fracture problem of a functionally graded coating–substrate system with a coating crack inclined to the weak/micro-discontinuous interface. The Cauchy singular integral equation for the crack was derived using Fourier integral transform, and the Lobatto–Chebyshev collocation method put up by Erdogan and Gupta was used to get its numerical solution. Finally, the effects of the weak/micro-discontinuity of the interface on SIFs were analyzed, the “affected regions” corresponding to the two crack tips have been obtained and their engineering significance was discussed. It was indicated that, for the crack tip in the corresponding “affected region”, to reduce the weak-discontinuity of the interface and to make the interface micro-discontinuous are the two effective ways to reduce the SIF, and the latter way always has more remarkable SIF-reduction effect. For the crack tip outside the “affected region”, its SIF is mainly influenced by material stiffness, and to prevent such a tip from growing toward the interface “softer coating and stiffer substrate” is a more advantageous combination than “stiffer coating and softer substrate”.     

冲击载荷作用下裂纹动态响应的数值模拟

对垂直、剪切以及斜向等各种冲击载荷作用下裂纹的动态响应进行了数值模拟,得到了一系列随时间变化的动态应力场以及应变场图;根据其定义,计算出了相应的动态应力强度因子,进而分析了斜向载荷作用下裂纹起裂情况,并对最优断裂问题进行了阐述。     

Torsional impact response of a penny-shaped crack in a functional graded strip

The torsional impact response of a penny-shaped crack in a nonhomogeneous strip is considered. The shear modulus is assumed to be functionally graded such that the mathematics is tractable. Laplace and Hankel transforms were used to reduce the problem to solving a Fredholm integral equation. The crack tip stress field is obtained by considering the asymptotic behavior of Bessel function. Explicit expressions of both the dynamic stress intensity factor and the energy density factor were derived. And it is shown that, as crack driving force, they are equivalent for the present crack problem. Investigated are the effects of material nonhomogeneity and strip‘s highness on the dynamic fracture behavior.Numerical results reveal that the peak of the dynamic stress intensity factor can be suppressed by increasing the nonhomogeneity parameter of the shear modulus, and that the dynamic behavior varies little with the adjusting of the strip‘ s highness.     

Fracture analysis on a layered multiferroic smart structure containing non-periodic interfacial cracks under magnetic/electric and mechanical loadings

The problems of periodic cracks have been widely solved in existing literature of fracture mechanics. Although the assumption of periodicity facilitates the theoretical derivation, practical multiple cracks in most cases are not periodically distributed in materials and structures. The present article performs fracture analysis on a piezoelectric/piezomagnetic/piezoelectric (PE/PM/PE) tri-layered smart structure containing nonperiodic interfacial cracks under magnetic/electric and mechanical loadings. The methods of dislocation simulation, Green's function and Cauchy singular integral equations are employed to solve the problem. Stress intensity factors (SIFs) are numerically calculated, and parametric studies are conducted to demonstrate the variation laws of SIFs versus geometrical parameters. The obtained conclusions may provide references for the optimal design of PE/PM/PE smart structures.     

Scattering of SH-wave by cracks originating at an elliptic hole and dynamic stress intensity factor

The method of complex function and the method of Green‘s function are used to investigate the problem of SH-wave scattering by radial cracks of any limited length along the radius originating at the boundary of an elliptical hole, and the solution of dynamicstress intensity factor at the crack tip was given. A Green‘s function was constructed for the problem, which is a basic solution of displacement field for an elastic half space containing a half elliptical gap impacted by anti-plane harmonic linear source force at any point of its horizontal boundary. With division of a crack technique, a series of integral equations can be established on the conditions of continuity and the solution of dynamic stress intensity factor can be obtained. The influence of an elliptical hole on the dynamic stress intensity factor at the crack tip was discussed.     

Effects of electric field on crack growth for a penny-shaped dielectric crack in a piezoelectric layer

The assumptions of impermeable and permeable cracks give rise to significant errors in determining electro-elastic behavior of a cracked piezoelectric material. The former simply imposes that the permittivity or electric displacement of the crack interior vanishes, and the latter neglects also the effects of the dielectric of an opening crack interior. Considering the presence of the dielectric of an opening crack interior and the permeability of the crack surfaces for electric field, this paper analyzes electro-elastic behavior induced by a penny-shaped dielectric crack in a piezoelectric ceramic layer. In the cases of prescribed displacement or prescribed stress at the layer surfaces, the Hankel transform technique is employed to reduce the problem to Fredholm integral equations with a parameter dependent nonlinearly on the unknown functions. For an infinite piezoelectric space, a closed-form solution can be derived explicitly, while for a piezoelectric layer, an iterative technique is suggested to solve the resulting nonlinear equations. Field intensity factors are obtained in terms of the solution of the equations. Numerical results of the crack opening displacement intensity factors are presented for a cracked PZT-5H layer and the effect of applied electric field on crack growth are examined for both cases. The results indicate that the fracture toughness of a piezoelectric ceramic is affected by the direction of applied electric fields, dependent on the elastic boundary conditions.     

Three-dimensional interfacial fracture analysis of a one-dimensional hexagonal quasicrystal coating

In this paper, the three-dimensional (3D) interfacial fracture is analyzed in a one-dimensional (1D) hexagonal quasicrystal (QC) coating structure under mechanical loading. A planar interface crack with arbitrary shape is studied by a displacement discontinuity method. Fundamental solutions of interfacial concentrated displacement discontinuities are obtained by the Hankel transform technique, and the corresponding boundary integral-differential equations are constructed with the superposition principle. Green’s functions of constant interfacial displacement discontinuities within a rectangular element are derived, and a boundary element method is proposed for numerical simulation. The singularity of stresses near the crack front is investigated, and the stress intensity factors (SIFs) as well as energy release rates (ERRs) are determined. Finally, relevant influencing factors on the fracture behavior are discussed.     

残余应力对压力容器破坏模式的影响分析

破前漏(简称LBB)是压力容器、核电站设备结构设计与评价中的一个重要准则.表面裂纹准静态扩展的几何形貌变化规律的预测是破前漏(LBB)评判十分重要的课题之一.本文对特定焊接残余应力场加载作用下,含三维表面裂纹的压力容器模型,用有限元软件(ABAQUS)进行了表面裂纹准静态扩展模拟计算,得到在此残余应力场作用下应力强度因子沿裂纹前缘的分布规律.结合外载引起的应力强度因子,就可以判别裂纹的扩展形貌,从而判断结构是否满足LBB要求.     

Mode III fracture of an arbitrary oriented crack in two dimensional functionally graded material

In this paper, a two dimensional functionally graded material (2D-FGM) under an anti-plane load with an internal crack is considered. The crack is oriented in an arbitrary direction. The material properties are assumed to vary exponentially in two planar directions. The problem is analyzed and solved by two different methods namely Fourier integral transforms with singular integral equation technique, and then by the finite element method. The effects of crack orientation, material non-homogeneity, and other parameters on the value of stress intensity factor (SIF) are studied. Finally, the obtained results for Mode III stress intensity factor of different methods are compared.     

Anti-plane shear of an arbitrary oriented crack in a functionally graded strip bonded with two dissimilar half-planes

An internal crack located within a functionally graded material (FGM) strip bonded with two dissimilar half-planes and under an anti-plane load is considered. The crack is oriented in an arbitrary direction. The material properties of strip are assumed to vary exponentially in the thickness direction and two half-planes are assumed to be isotropic. Governing differential equations are derived and to reduce the difficulty of the problem dealing with solution of a system of singular integral equations Fourier integral transform is employed. Semi closed form solution for the stress distribution in the medium is obtained and mode III stress intensity factor (SIF), at the crack tip is calculated and its validity was verified. Finally, the effects of nonhomogeneous material parameter and crack orientation on the stress intensity factor are studied.