3D dynamic stress intensity factors around two parallel square cracks in an infinite elastic medium under impact load

Summary  Transient stresses around two parallel cracks in an infinite elastic medium are investigated in the present paper. The shape of the cracks is assumed to be square. Incoming shock stress waves impinge upon the two cracks normal to tzheir surfaces. The mixed boundary value equations with respect to stresses and displacements are reduced to two sets of dual integral equations in the Laplace transform domain using the Fourier transform technique. These equations are solved by expanding the differences in the crack surface displacements in a double series of a function that is equal to zero outside the cracks. Unknown coefficients in the series are calculated using the Schmidt method. Stress intensity factors defined in the Laplace transform domain are inverted numerically to the physical space. Numerical calculations are carried out for transient dynamic stress intensity factors under the assumption that the shape of the upper crack is identical to that of the lower crack. Received 2 February 2000; accepted for publication 10 May 2000     

DYNAMIC STRESS INTENSITY FACTORS AROUND TWO CRACKS NEAR AN INTERFACE OF TWO DISSIMILAR ELASTIC HALF－PLANES UNDER IN－PLANE SHEAR IMPACT LOAD

ＤＹＮＡＭＩＣＳＴＲＥＳＳＩＮＴＥＮＳＩＴＹＦＡＣＴＯＲＳＡＲＯＵＮＤＴＷＯＣＲＡＣＫＳＮＥＡＲＡＮＩＮＴＥＲＦＡＣＥＯＦＴＷＯＤＩＳＳＩＭＩＬＡＲＥＬＡＳＴＩＣＨＡＬＦ－ＰＬＡＮＥＳＵＮＤＥＲＩＮ－ＰＬＡＮＥＳＨＥＡＲＩＭＰＡＣＴＬＯＡＤＱｉａｎＲ...     

Transient response of a magnetoelectroelastic solid with two collinear dielectric cracks under impacts

Dynamic analysis of two collinear electro-magnetically dielectric cracks in a piezoelectromagnetic material is made under in-plane magneto-electro-mechanical impacts. Generalized semi-permeable crack-face boundary conditions are proposed to simulate realistic opening cracks with dielectric. Ideal boundary conditions of a combination of electrically permeable or impermeable and magnetically permeable or impermeable assumptions are several limiting cases of the semi-permeable dielectric crack. Utilizing the Laplace and Fourier transforms, the mixed initial-boundary-value problem is reduced to solving singular integral equations with Cauchy kernel. Dynamic intensity factors of stress, electric displacement, magnetic induction and crack opening displacement (COD) near the inner and outer crack tips are determined in the Laplace transform domain. Numerical results for a special magnetoelectroelastic solid are calculated to show the influences of the dielectric permittivity and magnetic permeability inside the cracks on the crack-face electric displacement and magnetic induction. By means of a numerical inversion of the Laplace transform, the variations of the normalized intensity factors of stress and COD are discussed against applied magnetoelectric impact loadings and the geometry of the cracks for fully impermeable, vacuum, fully permeable cracks and shown in graphics.     

Dynamic stress intensity factors for two parallel interface cracks between a nonhomogeneous bonding layer and two dissimilar elastic half-planes subject to an impact load

Some composite materials are constructed of two dissimilar half-planes bonded by a nonhomogeneous elastic layer. In the present study, a crack is situated at the interface between the upper half-plane and the bonding layer of such a material, and another crack is located at the interface between the lower half-plane and the bonding layer. The material properties of the bonding layer vary continuously from those of the lower half-plane to those of the upper half-plane. Incoming shock stress waves impinge upon the two interface cracks normal to their surfaces. Fourier transformations were used to reduce the boundary conditions for the cracks to two pairs of dual integral equations in the Laplace domain. To solve these equations, the differences in the crack surface displacements were expanded in a series of functions that are zero-valued outside the cracks. The unknown coefficients in the series were solved using the Schmidt method so as to satisfy the conditions inside the cracks. The stress intensity factors were defined in the Laplace domain and were inverted numerically to physical space. Dynamic stress intensity factors were calculated numerically for selected crack configurations.     

Dynamic interaction of multiple arbitrarily oriented planar cracks in a piezoelectric space: A semi-analytic solution

The problem of determining the electro-elastic fields around arbitrarily oriented planar cracks in an infinite piezoelectric space is considered. The cracks which are acted upon by a transient load are either electrically impermeable or permeable. A semi-analytic method based on the theory of exponential Fourier transformation is proposed for solving the problem in the Laplace transform domain. The Laplace transforms of the jumps in the displacements and electric potential across opposite crack faces are determined by solving a system of hypersingular integral equations. Once these displacement and electric potential jumps are obtained, the displacements and electric potential and other physical quantities of interest, such as the crack tip stress and electric displacement intensity factors, can be computed with the help of a suitable algorithm for inverting Laplace transforms. The stress and electric displacement intensity factors are computed for some specific cases of the problem.     

SH波对一维六方准晶中直裂纹的散射

分析了SH波对一维六方准晶中直裂纹的散射问题。利用积分变换技术,结合Copson方法,通过求解对偶积分方程,得到声子场和相位子场应力、位移及裂纹尖端动应力强度因子的解析表达式。通过数值算例讨论了裂纹长度、入射角和入射波频率对标准动应力强度因子的影响,此研究在工程材料应用中有一定的参考价值。     

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.     

A new method for solving Biot＇s consolidation of a finite soil layer in the cylindrical coordinate system

A new method is developed to solve Biot＇s consolidation of a finite soil layer in the cylindrical coordinate system. Based on the governing equations of Biot＇s consolidation and the technique of Laplace transform, Fourier expansions and Hankel transform with respect to time t, coordinate θ and coordinate r, respectively, a relationship of displacements, stresses, excess pore water pressure and flux is established between the ground surface （z = 0） and an arbitrary depth z in the Laplace and Hankel transform domain. By referring to proper boundary conditions of the finite soil layer, the solutions for displacements, stresses, excess pore water pressure and flux of any point in the transform domain can be obtained. The actual solutions in the physical domain can be acquired by inverting the Laplace and the Hankel transforms.     

The stress intensity factors near the tips of two collinear cracks in composite under in-plane shear

In this paper, the stress-intensity factors for two collinear cracks in a composite bonded by an isotropic and an anisotropic half-plane were calculated. The cracks are paralell to the interface, and the crack surfaces are loaded by uniform shear stresses. By using Fourier transform, the mixed boundary value problem is reduced to a set of singular integral equations. For solving the integral equations, the crack surface displacements are expanded in triangular series and the unknown coefficients in the series are determined by the Schmidt method. The stress intensity factors for the cracks in the boron-fibre plastics and aluminium joined composite and in carbon-fibre reinforced plastics were calculated numerically.     

Opening-Function Simulation of the Three-Dimensional Nonstationary Interaction of Cracks in an Elastic Body

Integral relations between three-dimensional dynamic displacements (stresses) in an infinite elastic body with arbitrarily located plane cracks and discontinuities in the displacements of the opposite crack faces are presented. The influence of opening cracks on each other is considered in the problem on crack faces loaded by pulse forces. This problem is reduced to a system of boundary integral equations of the wave-potential type in a time domain. The dynamic mode I stress intensity factors are determined for two coplanar elliptic cracks under forces in the form of the Heaviside function     

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.     

饱和地基任意形状衬砌内部均布突加荷载作用下的动力响应

地下衬砌结构经常会受到内部动荷载的作用,内荷载引起的衬砌结构的动应力集中备受关注.论文利用波函数展开法和Laplace变换,推导了饱和土中突加荷载作用下衬砌结构和土体的位移、应力、孔压表达式.应用复变函数和保角变换,将任意形状边界映射为圆形边界,利用饱和土和衬砌结构的连续条件和边界条件,求得了任意形状衬砌结构的动力响应...     

Scattering of harmonic anti-plane shear waves by an interface crack in magneto-electro-elastic composites

IntroductionCompositematerialconsistingofapiezoelectricphaseandapiezomagneticphasehasdrawnsignificantinterestinrecentyears,duetotherapiddevelopmentinadaptivematerialsystems .Itshowsaremarkablylargemagnetoelectriccoefficient,thecouplingcoefficientbetweenst…     

Full field analysis of an anti-plane interface crack subjected to dynamic body forces

In this study, the transient full field response of an interface crack between two different media subjected to dynamic body force at one material is investigated. For time t < 0, the bimaterial medium is stress free and at rest. At t = 0, a concentrated anti-plane dynamic point loading is applied at the medium as shown in Fig. 1. The total wave field is due to the effect of this point loading and the scattering of the incident waves by the interface crack. An alternative methodology that is different from the conventional superposition method is used to construct the reflected, refracted and diffracted wave fields. A useful fundamental solution is proposed in this study and the full field solution is determined by superposition of the fundamental solution in the Laplace transform domain. The proposed fundamental problem is the problem of applying an exponentially distributed traction (in the Laplace transform domain) on the interfacial crack faces. The Cagniard–de Hoop method of Laplace inversion is used to obtain the transient solution in time domain. Exact transient closed form solutions for stresses and stress intensity factors are obtained. Numerical results for the time history of stresses and stress intensity factors during the transient process are discussed in detail.     

Dynamic analysis of a cracked magnetoelectroelastic medium under antiplane mechanical and inplane electric and magnetic impacts

The transient analysis of a magnetoelectroelastic medium containing a crack is made under antiplane mechanical and inplane electric and magnetic impacts. The crack is assumed to penetrate through the solid along the poling direction. By using the Fourier and Laplace transforms, the associated mixed boundary value problem is reduced to a Fredholm integral equation of the second kind, which is solved numerically. By means of a numerical inversion of the Laplace transform, dynamic field intensity factors are obtained in the time domain. Numerical results are presented graphically to show the effects of the material properties and applied electric and magnetic impacts on the dynamic intensity factors of COD and stress, and dynamic energy density factors. The results indicate that except for the intensity factors of electric displacement and magnetic induction, other field intensity factors exhibit apparent transient feature. Moreover, they depend strongly on mechanical input as well as electric and magnetic impacts.     

Basic solution of two parallel non-symmetric permeable cracks in piezoelectric materials

The behavior of two parallel non-symmetric cracks in piezoelectric materials subjected to the anti-plane shear loading was studied by the Schmidt method for the permeable crack electric boundary conditions. Through the Fourier transform, the present problem can be solved with two pairs of dual integral equations ip which the unknown variables are the jumps of displacements across crack surfaces. To solve the dual integral equations, the jumps of displacements across crack surfaces were directly expanded in a series of Jacobi polynomials. Finally, the relations between electric displacement intensity factors and stress intensity factors at crack tips can be obtained. Numerical examples are provided to show the effect of the distance between two cracks upon stress and electric displacement intensity factors at crack tips. Contrary to the impermeable crack surface condition solution, it is found that electric displacement intensity factors for the permeable crack surface conditions are much smaller than those for the impermeable crack surface conditions. At the same time, it can be found that the crack shielding effect is also present in the piezoelectric materials.     

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

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

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.     

Stress intensity factors for a moving cylindrical crack in a nonhomogeneous cylindrical layer in composite materials

Stresses are determined for a finite cylindrical crack that is propagating with a constant velocity in a nonhomogeneous cylindrical elastic layer, sandwiched between an infinite elastic medium and a circular elastic cylinder made from another material. The Galilean transformation is employed to express the wave equations in terms of coordinates that are attached to the moving crack. An internal gas pressure is then applied to the crack surfaces. The solution is derived by dividing the nonhomogeneous interfacial layer into several homogeneous cylindrical layers with different material properties. The boundary conditions are reduced to two pairs of dual integral equations. These equations are solved by expanding the differences in the crack surface displacements into a series of functions that are equal to zero outside the crack. The Schmidt method is then used to solve for the unknown coefficients in the series. Numerical calculations for the stress intensity factors were performed for speeds and composite material combinations.     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程,结合广义Betti-Rayleigh互易等式与时域下的边界积分方程,推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解,将经过主部分析的积分核函数分解为静态和动态部分,其中动态积分核不具有奇异性。在裂纹前沿附近单元,采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换,采用配置点法计算超奇异积分,获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序,得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律,数值结果稳定且收敛速度快。