功能梯度夹层多个环形界面裂纹扭转冲击

研究位于功能梯度层和外部均匀材料之间多个环形界面裂纹的扭转冲击问题,功能梯度材料 (FGM)粘结在两种不同的弹性材料之间,功能梯度层和外部材料之间环形界面裂纹的数目是任意的．引进积分变换和位错密度函数将问题化为求解Laplace域里标准的Cauchy奇异积分方程,进而化为求解代数方程;应用Laplace数值反演技术,计算时域里的动应力强度因子(DSIF)．考查了结构几何尺度和材料特性对裂尖动态断裂特性的影响．数值结果表明,DSIF存在一个主峰,到达主峰后,在其相应的静态值附近波动并最终趋于稳定;增加FGM的梯度能减小DSIF的峰值．     

Dynamic stress intensity factors of two collinear mode-III cracks perpendicular to and on the two sides of a bi-FGM weak-discontinuous interface

The mechanical model was established for the anti-plane dynamic fracture problem for two collinear cracks on the two sides of and perpendicular to a weak-discontinuous interface between two materials with smoothly graded elastic properties, as opposed to a sharp interface with discontinuously changing elastic properties. The problem was reduced as a system of Cauchy singular integral equations of the first kind by Laplace and Fourier integral transforms. The integral equations were solved by Erdogan's collocation method and the dynamic stress intensity factors in the time domain were obtained through Laplace numerical inversion proposed by Miller and Guy. The influences of geometrical and physical parameters on the dynamic stress intensity factors were illustrated and discussed, based on which some conclusions were drawn: (a) to increase the thickness of the FGM strip on either side of the interface will be beneficial to reducing the DSIF of a crack perpendicular to a bi-FGM interface and embedded at the center of one of the FGM strips; (b) To increase the rigidity of the FGM strip where the crack is located will increase the DSIF. However, when the material in one side of the interface is more rigid, the DSIF of the interface-perpendicular embedded crack in the other side will be reduced; (c) To decrease the weak-discontinuity of a bi-FGM interface will not necessarily reduce the stress intensity factor of a crack perpendicular to it, which is different from the case of interfacial crack; (d) For two collinear cracks with equal half-length, when the distance between the two inner tips is less than about three times of the half-length, the interaction of them is intensified, however, when the distance is greater than this the interaction becomes weak.     

The Anti-plane Dynamic Fracture Analysis of Functionally Graded Materials with Arbitrary Spatial Variations of Material Properties

Anti-plane dynamic fracture analysis is presented for functionally graded materials (FGM) with arbitrary spatial variations of material properties. The FGM with the material properties varying continuously in an arbitrary manner is modeled as a multi-layered medium with the elastic modulus and mass density varying linearly in each sub-layer and continuous at the interfaces between two adjacent sub-layers. With this linearly inhomogeneous multi-layered model, the problem of a crack in a graded interfacial zone bonded to two homogeneous half-spaces or in a coating bonded to a homogeneous half-space subjected to the anti-plane shear impact load is investigated. Laplace and Fourier transforms and transfer matrix are applied to reduce the associated mixed boundary value problem to a Cauchy singular integral equation which is solved numerically in the Laplace transformed domain. The dynamic stress intensity factors (DSIF) are obtained by using the numerical technique of Laplace inversion.     

Fracture analysis on a piezoelectric sensor with a viscoelastic interface

Fracture analysis is performed on a layered piezoelectric sensor possessing a Kelvin-type viscoelastic interface. An electrically permeable anti-plane crack is situated in the piezoelectric layer and perpendicular to the interface. The crack problem is solved by the methods of integral transform and Cauchy singular integral equation. The variations of the dynamic stress intensity factor (DSIF) vs. physical and geometrical parameters are investigated. At the beginning of creep and relaxation, larger viscosity coefficient always induces smaller DSIF. With time elapsing, the effect of viscosity coefficient becomes weaker and weaker. When time approaches infinity, the viscous effect disappears, and the DSIF converges to a value corresponding to the case of an elastic interface. The effect of the viscoelastic interface on the fracture behavior of the piezoelectric layer also depends on the substrate thickness. To some extent, thicker substrate may intensify the effect of the interface.     

Dynamic stress intensity factor of a crack perpendicular to the weak-discontinuous interface in a nonhomogeneous coating–substrate structure

A mechanical model was established for the antiplane dynamic fracture problem of a functionally graded coating–substrate structure with a coating crack perpendicular to the weak-discontinuous interface. The problem was reduced to a Cauchy singular integral equation by the methods of Laplace and Fourier integral transforms. Erdogan’s collocation method and the Laplace numerical inversion proposed by Miller and Guy were used to calculate the dynamic stress intensity factors. Three conclusions were drawn through parametric studies: (a) unlike the conclusion drawn for an interfacial crack, reducing the weak discontinuity of the interface will not necessarily decrease the dynamic stress intensity factor (DSIF) of the coating crack perpendicular to the interface; (b) increasing the stiffness of the substrate when that of the coating is fixed, or decreasing the stiffness of coating when that of the substrate is fixed, will be beneficial for the reduction of the DSIF of a coating crack perpendicular to the interface; and (c) the free surface has a greater influence on the DSIF than the interface does, and the effect of the interface on the DSIF is greater than that of the material stiffness in the crack-tip region.     

Elastodynamic analysis of a crack at an arbitrary angle to the graded interfacial zone in bonded half-planes under antiplane shear impact

A solution is provided for the elastodynamic problem of a crack at an arbitrary angle to the graded interfacial zone in bonded media under the action of antiplane shear impact. The interfacial zone is modeled by a nonhomogeneous interlayer with the spatially varying shear modulus and mass density in terms of power functions between the two dissimilar, homogeneous half-planes. Based on the use of Laplace and Fourier integral transforms and the coordinate transformations of basic field variables, formulation of the transient crack problem is reduced to solving a Cauchy-type singular integral equation in the Laplace transform domain. The crack-tip response in the physical domain is recovered via the inverse Laplace transform and the values of dynamic mode III stress intensity factors are obtained as a function of time. A comprehensive parametric study is then presented of the effects of crack obliquity on the overshoot behavior of the transient crack-tip response, by plotting the peak values of the dynamic stress intensity factors versus the crack orientation angle for various material and geometric combinations of the bonded system.     

Antiplane mechanical and inplane electric time-dependent load applied to two coplanar cracks in piezoelectric ceramic material

This work is concerned with the dynamic response of two coplanar cracks in a piezoelectric ceramic under antiplane mechanical and inplane electric time-dependent load. The cracks are assumed to act either as an insulator or as a conductor. Laplace and Fourier transforms are used to reduce the mixed boundary value problems to Cauchy-type singular integral equations in Laplace transform domain. A numerical Laplace inversion algorithm is used to determine the dynamic stress and electric displacement factors that depend on time and geometry. A normalized equivalent parameter describing the ratio of the equivalent magnitude of electric load to that of mechanical load is introduced in the numerical computation of the dynamic stress intensity factor (DSIF) which has a similar trend as that for the pure elastic material. The results show that the dynamic electric field will impede or enhance crack propagation in a piezoelectric ceramic material at different stages of the dynamic electromechanical load. Moreover, the electromechanical response is greatly affected by the ratio of the crack length to the ligament between the cracks. The stress and electric displacement intensity factor can be combined by the energy density factor or function to address the fracture of piezoelectric materials under the combined influence of electromechanical loading.     

Torsional impact of transversely isotropic solid with functionally graded shear moduli and a penny-shaped crack

The torsional impact response of a penny-shaped crack in an unbounded transversely isotropic solid is considered. The shear moduli are assumed to be functionally graded such that the mathematics is tractable. Laplace transform and Hankel transform are used to reduce the problem to solving a Fredholm integral equation. The crack tip stress fields are obtained. Investigated are the influence of material nonhomogeneity and orthotropy on the dynamic stress intensity factor. The peak value of the dynamic stress intensity factor can be suppressed by increasing the shear moduli's gradient and/or increasing the shear modulus in a direction perpendicular to the crack surface.     

Dynamic fracture analysis of adhesive bonded material under normal loading

Dynamic fracture behavior of a Griffith crack along the interface of an adhesive bonded material under normal loading is studied. The singular integral equations are obtained by employing integral transformation and introducing dislocation density functions. By adopting Gauss-Jacobi integration formula, the problem is reduced to the solution of algebraic equations, and by collocation dots method. their solutions can be obtained Based on the parametric discussions presented in the paper, the following conclusions can be drawn： （1） Mode I dynamic stress intensity factor （DSIF） increases with increasing initial crack length and decreasing visco-elastic layer thickness, revealing distinct size effect; （2） The influence of the visco-elastic adhesive relaxation time on the DSIF should not be ignored.     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。     

Impact failure prediction of Mode III crack in orthotropic functionally graded strip

Mode III impact of a crack in an orthotropic functionally graded strip is investigated. The shear moduli in two directions of the material are assumed to vary proportionately with gradient. Laplace transform and Fourier cosine transform are 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. Energy density factor criterion is applied to obtain the maximum of minimum energy density and direction of crack initiation. Numerical results are given graphically. The effects of orthotropy, nonhomogeneity and height of the strip on the energy density factor are discussed.     

Fracture analysis of a weak-discontinuous interface in a symmetrical functionally gradient composite strip loaded by anti-plane impact

The anti-plane impact fracture analysis was performed for a weak-discontinuous interface in a symmetrical functionally gradient composite strip. A new bi-parameter exponential function was introduced to simulate the continuous variation of material properties. Using Laplace and Fourier integral transforms, we reduced the problem to a dual integral equation and obtained asymptotic analytical solution of crack-tip stress field. Based on the numerical solution of the second kind of Fredholm integral equation transformed from the dual integral equation, the effects of the two non-homogeneity parameters on DSIF were discussed. It was indicated that the relative stiffness of the interface and the general stiffness of the whole structure are two important factors affecting the impact fracture behavior of the weak-discontinuous interface. The greater the relative stiffness of the interface is, the higher the value of the dynamic stress intensity factor will be. The greater the general stiffness of the whole structure is, the shorter the time for DSIF to arrive at the peak value and then to stabilize to the steady one. If the general stiffness of the whole structure is great enough, there will be an oscillation between the peak and steady values of DSIF.     

Transient dynamic response of a coated piezoelectric strip with a vertical crack

In this study, the dynamic response of a coated piezoelectric strip containing a crack vertical to the interfaces under normal impact load is considered. Based on the superposition principle and the integral transform techniques, the solution in the Laplace transformed plane is obtained in terms of a singular integral equation. The order of stress singularity around the tip of the terminated crack is also obtained. The singular integral equation is solved by using the Gauss–Jacobi integration formula, and the numerical Laplace inversion is then carried out to obtain the resulting dynamic stress and electric displacement intensities. The effects of the material properties and the geometric parameters on the dynamic stress intensity factor and the dynamic energy density factors are shown graphically.     

Dynamic fracture analysis of an annular interfacial crack between dissimilar magnetoelectroelastic layers

Transient response of an annular interfacial crack between dissimilar magnetoelectroelastic layers under impacts is investigated. On the crack surface, magnetoelectrically impermeable boundary condition is adopted. Using Laplace and Hankel transform techniques, the mixed boundary value problem is reduced to a system of singular integral equations. The integral equations are further reduced to a system of algebraic equations with the aid of Jacobi polynomials. The dynamic field intensity factor and dynamic energy release rate are determined. Numerical results reveal the effects of electric or magnetic loadings and material parameters of composite on crack propagation and growth.     

A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.     

Dynamic analysis of a penny-shaped dielectric crack in a magnetoelectroelastic solid under impacts

The transient response of a magneto-electro-elastic material with a penny-shaped dielectric crack subjected to in-plane magneto-electro-mechanical impacts is made. To simulate an opening crack with a dielectric interior, the crack-face electromagnetic boundary conditions are supposed to depend on the crack opening displacement and the jumps of electric and magnetic potentials across the crack. Four ideal crack-face electromagnetic boundary conditions involving a combination of electrically permeable or impermeable and magnetically permeable or impermeable assumptions can be reduced. The Laplace and Hankel transform techniques are further utilized to solve the mixed initial-boundary-value problem. Three coupling Fredholm integral equations are obtained and solved by the composite Simpson's rule. Dynamic field intensity factors of stress, electric displacement, magnetic induction, crack opening displacement (COD), electric potential and magnetic potential are given in the Laplace transform domain. By means of a numerical inversion of the Laplace transform, numerical results are calculated to show the variations of the physical parameters of concern versus the normalized time in graphics. The effects of applied electric and magnetic loads on the dynamic intensity factors of stress and COD, and the dynamic energy release rate for a BaTiO₃-CoFe₂O₄ composite with a penny-shaped vacuum crack are discussed in detail.     

Impact response of a piezoelectric layered composite plate with a crack

The dynamic response of a central crack in a piezoelectric layered composite plate under normal impact is analyzed. The crack is oriented normally to the interfaces. The Laplace and Fourier transform techniques are used to formulate the problem in terms of a singular integral equation. The order of stress singularity around the tip of the terminated crack is also obtained. Numerical calculations are carried out, and the main results presented are the variations of the dynamic stress intensity factor and the dynamic energy density factor versus time as functions of the geometric parameters and the piezoelectric material properties of the layered composite plate.     

动态载荷下功能梯度复合材料的圆币形裂纹问题

研究了动态载荷下功能梯度材料中的圆币形裂纹问题．假设材料为横观各向同性,并且含有多个垂直于厚度方向的裂纹,材料参数沿轴向（与裂纹面垂直的方向）为变化的,沿该方向将材料划分为许多单层,各单层材料参数为常数,利用Hankel变换祛,在Laplace域内推导出了控制问题的对偶积分方程组．利用Laplace数值反演,得出了裂纹尖端的动态应力强度因子和能量释放率．研究了含两个裂纹的功能梯度接头结构,分析了材料非均匀性参数对应力强度因子和能量释放率的影响．     

SCATTERING OF ANTI-PLANE SHEAR WAVES BY A SINGLE CRACK IN AN UNBOUNDED TRANSVERSELY ISOTROPIC ELECTRO-MAGNETO-ELASTIC MEDIUM

A theoretical treatment of the scattering of anti-plane shear (SH) waves is provided by a single crack in an unbounded transversely isotropic electro-magneto-elastic medium. Based on the differential equations of equilibrium, electric displacement and magnetic induction intensity differential equations, the governing equations for SH waves were obtained. By means of a linear transform, the governing equations were reduced to one Helmholtz and two Laplace equations. The Cauchy singular integral equations were gained by making use of Fourier transform and adopting electro-magneto imperme ableboundary conditions. The closed form expression for the resulting stress intensity factor at the crack was achieved by solving the appropriate singular integral equations using Chebyshev polynomial. Typical examples are provided to show the loading frequency upon the local stress fields around the crack tips. The study reveals the importance of the electro-magneto-mechanical coupling terms upon the resulting dynamic stress intensity factor.     

A PENNY-SHAPED CRACK IN A MAGNETOELECTROELASTIC LAYER UNDER RADIAL SHEAR IMPACT LOADING

This paper analyzes the dynamic magnetoelectroelastic behavior induced by a penny- shaped crack in a magnetoelectroelastic layer.The crack surfaces are subjected to only radial shear impact loading.The Laplace and Hankel transform techniques are employed to reduce the prob- lem to solving a Fredholm integral equation.The dynamic stress intensity factor is obtained and numerically calculated for different layer heights.And the corresponding static solution is given by simple analysis.It is seen that the dynamic stress intensity factor for cracks in a magnetoelec- troelastic layer has the same expression as that in a purely elastic material.And the influences of layer height on both the dynamic and static stress intensity factors are insignificant as h/a>2.