Transient analysis of a propagating crack with finite length subjected to a horizontally polarized shear wave

In this study, the transient response of a finite crack subjected to an incident horizontally polarized shear wave and then propagated with a constant speed in an unbounded elastic solid is investigated. Initially, the finite crack with crack length l is stress-free and at rest. At time t = 0, an incident horizontally polarized shear wave strikes at one of the crack tips and will arrive at the other tip at a later time. Then, two crack tips propagate along the crack tip line with different velocities as the corresponding stress intensity factors reach their fracture toughness. The correspondent configuration is shown in Fig. 1

1. Download : Download high-res image (31KB)
2. Download : Download full-size image

Fig. 1. Configuration and coordinate systems of a finite crack in an unbounded medium.

. In analyzing this problem, diffracted waves generated by two propagating crack tips must be taken into account and it makes the analysis extremely difficult. In order to solve this problem, the transform formula in the Laplace transform domain between moving and stationary coordinates is first established. Complete solutions are determined by superposition of proposed fundamental solutions in the Laplace transform domain. The fundamental solutions to be used are from the problems of applying exponentially distributed traction and screw dislocation on crack faces and along the crack tip line, respectively. The exact transient solutions of dynamic stress intensity factor for the first few diffracted waves that arrive at two crack tips are obtained and expressed in compact formulations. Numerical calculations of dynamic stress intensity factors for both tips are evaluated and the results are discussed in detail.     

Dynamic stress intensity factors around a cylindrical crack in an infinite elastic medium subject to impact load

This paper investigates transient stresses around a cylindrical crack in an infinite elastic medium subject to impact loads. Incoming stress waves resulting from the impact load impinge on the crack in a direction perpendicular to the crack axis. In the Laplace transform domain, by means of the Fourier transform technique, the mixed boundary value equations with respect to stresses and displacements were reduced to two sets of dual integral equations. To solve the equations, the differences in the crack surface displacements were expanded in a series of functions that are zero outside the crack. The boundary conditions for the crack were satisfied by means of the Schmidt method. Stress intensity factors were defined in the Laplace transform domain and were numerically inverted to physical space. Numerical calculations were carried out for the dynamic stress intensity factors corresponding to some typical shapes assumed for the cylindrical crack.     

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     

The extended Durbin method and its application for piezoelectric wave propagation problems

Based on the well-known Durbin method, an efficient numerical method was developed for the inversion of the two-sided Laplace transform. The accuracy of the method was verified using examples. As an application of the method, transient elastic waves propagating in a two-layered piezoelectric medium subjected to anti-plane concentrated loading and in-plane electric displacement loading were investigated. One-sided and two-sided Laplace transforms were applied to determine the shear stresses and electric displacements in the double Laplace transform domain. Subsequently, the Durbin method for one-sided Laplace transform inversion and the extended Durbin method for two-sided Laplace transform inversion were used to implement the numerical inversions. Additionally, the numerical results of the transient stresses and electric displacements were evaluated and discussed. It showed that the arrival time of transient waves satisfies physical phenomena, and the transient solution oscillates near the static solution and rapidly approximates the static solution.     

A crack in a piezoelectric layer bonded to a dissimilar elastic layer under transient load

Summary  A piezoelectric layer bonded to the surface of an elastic structure is considered. The piezoelectric and the elastic layers are infinite along the x-axis and have finite thickness in the y-direction. The polarization direction of the piezoelectric material is along the y-axis. By means of the method of singular integral equations, the solution in a Laplace transform plane is demonstrated. Laplace inversion yields the results in the time domain. Numerical values of the crack tip fields under in-plane transient electromechanical loading are obtained. The influence of layers thickness on stress and electric displacement intensity factors is investigated. Received 16 March 2000; accepted for publication 16 August 2000     

热对流作用下筒壁涂层的边裂行为

边裂(边缘开裂)是涂层热致损伤的主要模式之一. 边缘裂纹穿透涂层后,常导致界面脱粘从而驱使涂层与基体剥离,最终丧失对基体的保护作用. 本文以热应力强度因子表征边缘裂纹的扩展驱动力,研究筒壁涂层在热对流作用下的边裂行为. 首先,利用拉普拉斯变换法,得到了瞬态温度场及热应力场的封闭解. 其次,运用Fett等的三参数法确定了筒壁涂层边缘裂纹的权函数. 最后,基于叠加原理和权函数方法计算了边缘裂纹的热应力强度因子. 探讨了无量纲时间、边缘裂纹深度、基体/涂层厚度比、热对流强度等参数对热应力强度因子的影响规律. 结果表明:热应力强度因子的峰值既非发生在热载荷初始时刻,也非发生在热稳态时刻,而出现在时间历程的中间时刻;增大热对流强度不仅可提高热应力强度因子的峰值,而且使峰值提前出现;其他条件相同时,热应力强度因子随着边缘裂纹长度的增大而降低;增大涂层厚度或减小基体厚度可增强涂层抵抗瞬态热载荷的能力.      

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.     

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.     

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.     

Solutions for anti-plane problems of a composite material with a crack terminating at the interface

Previous researchers have investigated anti-plane problems of a crack terminating at an interface within an infinite domain. This study is devoted to the theoretical analysis of these crack terminating problems for a circular composite material with a finite radius. Detailed solutions for a composite with a crack terminating at the interface under anti-plane loading were derived by employing a double transform that consists of the finite Mellin and Laplace transforms. Using the proposed method, related problems of anti-plane shearing can be solved in a straightforward manner. According to these solutions, stress intensity factors are obtained and discussed. Based on the problems studied with the maximum shear stress criterion, the crack may be interfacial debonded or may penetrate through the interface, but it cannot be reflected back.     

Thermoelastic transient response of an infinitely long annular multilayered cylinder

Thermoelastic transient response of multilayered annular cylinders of infinite lengths subjected to known temperature at traction-free inner and outer surfaces are considered. A method based on the Laplace transformation and finite difference method has been developed to analyze the thermoelasticity problem. Using the Laplace transform with respect to time, the general solutions of the governing equation are obtained in transform domain. The solution is obtained by using the matrix similarity transformation and inverse Laplace transform. Solutions for the temperature and thermal stress distributions in a transient state were obtained. It was found that the temperature distribution, the displacement and the thermal stresses change slightly as time increases. There is no limit of number of annular layers of the cylinder in the presented computational procedures.     

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

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

Hybrid graded element model for transient heat conduction in functionally graded materials

This paper presents a hybrid graded element model for the transient heat conduction problem in functionally graded materials (FGMs). First, a Laplace transform approach is used to handle the time variable. Then, a fundamental solution in Laplace space for FGMs is constructed. Next, a hybrid graded element is formulated based on the obtained fundamental solution and a frame field. As a result, the graded properties of FGMs are naturally reflected by using the fundamental solution to interpolate the intra-element field. Further, Stefest’s algorithm is employed to convert the results in Laplace space back into the time-space domain. Finally, the performance of the proposed method is assessed by several benchmark examples. The results demonstrate well the efficiency and accuracy of the proposed method.     

The transient fracture behavior for a functionally graded layered structure subjected to an in-plane impact load

The transient fracture behavior of a functionally graded layered structure subjected to an in-plane impact load is investigated. The studied structure is composed of two homogeneous layers and a functionally graded interlayer with a crack perpendicular to the boundaries. The impact load is applied on the face of the crack. Fourier transform and Laplace transform methods are used to formulate the present problem in terms of a singular integral equation in Laplace transform domain. Considering variations of parameters such as the nonhomogeneity constant, the thickness ratio and the crack length, the dynamic stress intensity factors (DSIFs) in time domain are studied and some meaningful conclusions are obtained.The project supported by the National Science Foundation for Excellent Young Investigators (10325208), the National Natural Science Foundation of China (10432030) and the China Postdoctoral Science Foundation (2004036018)The English text was polished by Ron Marshall.     

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.     

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

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

Diffraction of a plane stress wave by a semi-infinite crack in a general anisotropic elastic material

The problem of a semi-infinite crack subjected to an incident stress wave in a general anisotropic elastic solid is considered. The plane wave impinges the crack at a general oblique angle and is of any of the three types propagating in that direction. A related problem of a semi-infinite crack loaded by a pair of concentrated forces moving along the crack surfaces is also considered. In contrast to the conventional approach by Laplace transforms, a Stroh-like formalism is employed to construct the solution directly in the time domain. The solution is shown to depend on a Wiener–Hopf factorization of a symmetric matrix. Closed-form solution of the stress intensity factors is derived. A remarkably simple expression for the energy release rate is obtained for normal incidence.     

Dynamic investigation of a functionally graded layered structure with a crack crossing the interface

The dynamic response of a functionally graded layered structure with a crack crossing the interface is analyzed. The in-plane impact loading condition is considered. By using the Laplace and Fourier integral transforms, singular integral equation method and residue theory, the present problem is reduced to a singular integral equation in the Laplace transform domain. The influences of Young’s modulus ratio, thickness ratio, and crack length and location on the dynamic stress intensity factors (DSIFs) are investigated. Particularly, the DSIFs corresponding to different crack locations are shown in the case when the crack center moves from one layer to another layer through the interface. The peak and static values and overshoot characteristics of the DSIFs are analyzed. It is found that these values typically exhibit kinking behavior when the crack tips arrive at the interface. This study is different from previous other investigations in the following respects: (1) the dynamic response of a crack crossing the interface of a functionally graded structure is studied analytically, which has hardly been done in the past and (2) the present model can be reduced to some important problems, such as a functionally graded coating-substrate structure with a crack in the graded coating or homogeneous substrate or one intersecting the interface.     

Theoretical analysis,numerical calculation and experimental measurement of transient elastic waves in strips subjected to dynamic loadings

In this study, the transient response of an elastic strip subjected to dynamic in-plane loadings on the surface is investigated in detail. One of the objectives of this study is to develop an effective analytical method for determining transient solutions in a strip. By applying Laplace transform, the analytical solution in the transformed domain is derived and expressed in matrix form. The solution is then decomposed into infinite wave groups in which the multiple reflected waves with the same reflection are involved. Each multi-reflected wave can be identified by a coding method and be verified by the theory of generalized ray. The inverse transform is performed by using the well-known Cagniard method. The transient solutions in time domain for stresses and displacements are expressed in a closed form and are discussed in detail by an example. The experimental results show that the early time transient responses of displacements on the surface agree very well with the numerical calculations based on the theoretical solutions.     

New,real fundamental solutions to the transient thermal contact problem in a piezoelectric strip under the coupling actions of a rigid punch and a convective heat supply

A thermo-electro-mechanical contact analysis has been performed for a finite piezoelectric strip, which is subjected to the joint actions of a rigid, flat punch and a transient convective heat supply. The Laplace transform and Fourier sine and cosine transforms were applied in solving the governing equations. A detailed analysis of the characteristic roots of the corresponding characteristic equation was made. Real fundamental solutions were derived, which can readily lead to real solutions to the thermo-electro-mechanical quantities. A Cauchy-type singular integral equation was obtained for the stated problem and then solved numerically. Closed form solutions of a special case were obtained. To obtain the accurate solution in the time domain, an effective numerical inversion algorithm of the Laplace transform was applied. Detailed analyses were performed to reveal the variation law of temperature, contact stress beneath the punch, stress intensity factor at the punch edge and strain with time. Parametric studies were performed to discover the effects of the layer thickness on the distribution of temperature, contact stress beneath the punch and stress intensity factor at the punch edge.