A moving Griffith crack at the interface of two dissimilar elastic layers

The anti-plane shear problem of a Griffith crack traveling with a constant velocity at the interface of two dissimilar isotropic elastic layers is considered. Integral transform method is used to reduce the problem to the solution of a singular integral equation which is further reduced, by using Chebyshev polynomials, to a system of algebraic equations. The results for the particular cases of a moving Griffith crack at the interface of a layer and a half-space and two half-spaces are derived. Numerical results for the stress intensity factor are displayed graphically.     

Axisymmetrically loaded thin circular plate in adhesive contact with an elastic half-space

The interaction problem of adhesive contact between an axisymmetrically loaded thin circular plate and an isotropic or transversely isotropic elastic half-space is reduced to a problem of solving a pair of coupled integral equations for the unknown normal and shearing interfacial tractions. The system of integral equations is solved numerically and some results are presented.     

Axisymmetric time-harmonic response of a transversely isotropic substrate–coating system

By virtue of a method of displacement potentials, an analytical treatment of the response of a transversely isotropic substrate–coating system subjected to axisymmetric time-harmonic excitations is presented. In determination of the corresponding elastic fields, infinite line integrals with singular complex kernels are encountered. Branch points, cuts, and poles along the path of integration are accounted for exactly, and the physical phenomena pertinent to wave propagation in the medium are also highlighted. For evaluation of the integrals at the singular points, an accurate analytical residual theory is presented. Comparisons with the existing numerical solutions for a two-layered transversely isotropic medium under static surface load, and a transversely isotropic half-space subjected to buried time-harmonic load are made to confirm the accuracy of the present solutions. Selected numerical results for displacements and stresses are presented to portray the dependence of the response of the substrate–coating system on the frequency of excitation and the role of coating layer.     

Boundary integral equation method for conductive cracks in two and three-dimensional transversely isotropic piezoelectric media

Using the fundamental solutions and the Somigliana identity of piezoelectric medium, the boundary integral equations are obtained for a conductive planar crack of arbitrary shape in three-dimensional transversely isotropic piezoelectric medium. The singular behaviors near the crack edge are studied by boundary integral equation approach, and the intensity factors are derived in terms of the displacement discontinuity and the electric displacement boundary value sum near the crack edge on crack faces. The boundary integral equations for two dimensional crack problems are deduced as a special case of infinite strip planar crack. Based on the analogy of the obtained boundary integral equations and those for cracks in conventional isotropic elastic material and for contact problem of half-space under the action of a rigid punch, an analysis method is proposed. As an example, the solution to conductive Griffith crack is derived.     

An integral formulation for steady-state elastoplastic contact over a coated half-plane

 A boundary-domain integral equation for a coated half-space (elastically isotropic homogeneous substratum, possibly anisotropic coating layer) is developed. The half-space fundamental solution is used, so that the discretization is limited to the potential contact zone (boundary elements), the potentially plastic part of the substratum and the coating layer (domain integration cells). Steady-state elastoplastic analysis is implemented within this framework, for plane-strain conditions, for solving rolling and/or sliding contact problems, where at the moment the contact load comes from either a purely elastic contact analysis or is of Hertz type. The constitutive integration is of implicit type. In order to improve accuracy and computational efficiency, infinite elements are used. Comparison of numerical results with other sources, when available, is satisfactory. The present formulation is also used to compute the contact pressure for an isotropic (or anisotropic) coating on an isotropic homogeneous half-space indented by an elastic punch. Received 29 May 2001     

Axial soil-pile interaction in a transversely isotropic half-space

A rigorous integral equation formulation is presented for the axisymmetric load-transfer analysis of a thin-walled pile embedded in a transversely isotropic half-space under axial load. By virtue of a set of ring-load Green’s functions for the pile and one for the half-space, the problem is shown to be reducible to a pair of Fredholm integral equations. Through a mathematical analysis of an auxiliary pair of Cauchy integral equations, the inherent singularities of the contact stress distributions are rendered explicit. With a direct incorporation of the singular nature of the resultant load transfers, the numerical solution of the integral equations is shown to be possible by an interpolation of regular functions only. Typical results for various material and geometrical conditions are presented, including a comparison with past classical solutions, to illustrate the effects of transverse anisotropy on the load-transfer process.     

Spherical indentation of a transversely isotropic elastic half-space reinforced with a thin layer

The frictionless spherical indentation test is considered for a transversely isotropic elastic half-space reinforced with a thin layer whose flexural stiffness is negligible compared to its tensile stiffness. It is assumed that the deformation of the reinforcing layer can be treated as the generalized plane stress state. Closed-form analytical approximate equations for the maximum contact pressure, contact radius, and contact force are presented. The isotropic case is considered in detail.     

Effect of elastic property inhomogeneity for a ceramic substrate on separation conditions for a glass coating

Results are presented for a study of the contact interaction of a glass coating with a ceramic substrate. The study is carried out within the framework of classical linear elasticity theory: the glass coating is described by the Kirchhoff-Love thin plate theory, and the substrate is modeled by an elastic isotropic half-space with a variable Poisson's ratio over the depth. The possibility of separation of the edge of the coating from the substrate is considered. Solution of the corresponding edge problem of linear elasticity theory is reduced to a Fredholm integral equation of the second kind by methods of the theory of integral transformation. In order to determine the functional dependence of substrate elastic property inhomogeneity a procedure is described for numerical determination of values of parameters for this dependence with which the edge of the coating does not separate.Translated from Problemy Prochnosti, No. 2, pp. 36–39, February, 1991.     

重力对具有表面层的半空间中Rayleigh波的影响

重力对具有表面层的半空间中Rayleigh波具有重要的影响,假设表面层和半空间均为均匀各向同性介质,首先利用考虑重力作用的运动方程得出了重力影响下具有表面层的半空间中Rayleigh波的弥散方程。该方程经退化后得出了不受重力作用时的弥散方程,且与忽略重力时得到的方程完全一致。然后利用数值方法得到了重力影响下的Rayleigh波弥散曲线,分析了泊松比和表面层的厚度的影响,结果表明泊松比和表面层的厚度对弥散曲线具有明显影响。     

Anti-plane problem of periodic interface cracks in a functionally graded coating-substrate structure

In this paper, the interface cracking between a functionally graded material (FGM) and an elastic substrate is analyzed under antiplane shear loads. Two crack configurations are considered, namely a FGM bonded to an elastic substrate containing a single crack and a periodic array of interface cracks, respectively. Standard integral-transform techniques are employed to reduce the single crack problem to the solution of an integral equation with a Cauchy-type singular kernel. However, for the periodic cracks problem, application of finite Fourier transform techniques reduces the solution of the mixed-boundary value problem for a typical strip to triple series equations, then to a singular integral equation with a Hilbert-type singular kernel. The resulting singular integral equation is solved numerically. The results for the cases of single crack and periodic cracks are presented and compared. Effects of crack spacing, material properties and FGM nonhomogeneity on stress intensity factors are investigated in detail.     

饱和弹性半空间地基与中厚圆板的动力相互作用

利用Fourier-Bessel级数,对横观各向同性饱和弹性半空间地基与中厚圆板的动力相互作用问题进行了系统地分析。板的位移函数、荷载、地基反力及板下地基表面的沉降均被展开为二重Fourier-Bessel级数,这些级数中的待定系数由板的边界条件、板的控制方程及板-地基的相容条件加以确定,从而将饱和弹性半空间地基与中厚圆板的动力相互作用问题转化为数值积分和代数方程组的求解问题。数值计算表明,该级数解答具有较快的收敛速度。可以将该方法推广到弹性半空间与梁、矩形板相互作用问题的分析。     

位于弹性半空间上的理想流体层动力反应—平面P波入射

根据弹性介质与理想流体交界面的边界条件以及理想流体层自由表面的边界条件,基于位移势函数,推导了平面P波从弹性半空间入射到与理想流体层的交界面时,确定交界面波的透射与反射系数的理论方程。根据该方程,通过数值算例,分析了弹性半空间的弹性模量和密度、平面P波的入射角以及入射频率对理想流体层动压力的影响。     

Acoustic microscopy measurement of elastic constants and massdensity

A method is presented to determine the elastic constants and the mass density of isotropic and anisotropic solids and anisotropic thin films. The velocity and attenuation of leaky surface acoustic waves (SAWs) have been obtained for specified propagation directions from V(z) curves measured by line-focus acoustic microscopy (LFAM). The experimentally obtained velocities have been compared to velocities obtained from a measurement model for the V(z) curve which simulates the experiment. Since the measured and simulated V(z) curves have the same systemic errors, the material constants are free of such errors. For an isotropic solid, Young's modulus E, the shear modulus G and the mass density ρ have been determined from the leaky Rayleigh wave velocity and attenuation, measured by LFAM, and a longitudinal wave velocity measured by a pulse-echo transit-time technique. For a cubic-crystalline solid, the ratios of the elastic constants to the mass density (c₁₁ /ρ, c₁₂/ρ, c₄₄/ρ) have been determined from the directional variation of measured SAW velocities, using a preliminary estimate of ρ. The mass density ρ has subsequently been determined by additionally using the attenuation of leaky SAWs in crystal symmetry directions. For a cubic-crystalline thin film deposited on a substrate, the elastic constants and the mass density (c₁₁, c₁₂, c₄₄, ρ) of the film have been determined from the directional variation of the measured SAW velocities, and a comparison of the corresponding attenuation coefficient with the measured attenuation coefficient has been used to verify the results     

Elastodynamic analysis of nonplanar cracks in a half-space

Summary The interaction of time harmonic antiplane shear waves with nonplanar cracks embedded in an elastic half-space is studied. Based on the qualitatively similar features of crack and dislocation, with the aid of image method, the problem can be formulated in terms of a system of singular integral equations for the density functions and phase lags of vibrating screw dislocations. The integral equations, with the dominant singular part of Hadamard's type, can be solved by Galerkin's numerical scheme. Resonance vibrations of the layer between the cracks and the free surface are observed, which substantially give rise to high elevation of local stresses. The calculations show that near-field stresses due to scattering by a single crack and two cracks are quite different. The interaction between two cracks is discussed in detail. Furthermore, by assuming one of the crack tips to be nearly in contact with the free surface, the problem can be regarded as the diffraction of elastic waves by edge cracks. Numerical results are presented for the elastodynamic stress intensity factors as a function of the wave number, the incident angle, and the relative position of the cracks and the free surface.     

Effect of periodic surface roughness on V(z) curves for the line-focus acoustic microscope

A specimen with a periodic surface profile is considered to estimate the effect of surface roughness on the V(z) curve for the line-focus acoustic microscope. The Fourier optics approach is used to obtain the response of the lens and the Rayleigh-Fourier method is used to obtain the reflection coefficients for plane wave incidence from the fluid side on the periodic surface. An integral expression is obtained to calculate V(z) curves for periodic surface profiles. The V(z) curves are used to calculate the leaky Rayleigh wave velocities by applying the fast Fourier technique. Numerical results are presented to display the effect of sinusoidal surface roughness on the V(z) curves and the corresponding leaky Rayleigh wave velocities.     

Axisymmetric solution for a semi-infinite elastic circular cylinder of one material indenting an elastic half-space of another material

The axisymmetric bonded contact problem of a semi-infinite right circular cylinder of one elastic material indenting a half-space of a different elastic material is reduced to a system of singular integral equations of the second kind. The kernels of the integral equations are found to contain Cauchy and generalized Cauchy-type singularities. The index of the singularity for various material parameters combinations is determined by solving a characteristic determinant, which is obtained by considering the dominant part of the kernels. Using a modification of the method employed in [1], the system of singular integral equations is reduced to a system of simultaneous algebraic equations. The latter may then be solved numerically as in [1].     

位于弹性半空间上的理想流体层动力反应--平面SV波入射

根据弹性固体与理想流体动力学方程,导出了固-液介质交界面上波的透射与反射系数的理论公式,分析了平面SV波从弹性半空间入射到与理想流体层的交界面时,弹性半空间的弹性模量和密度,平面SV波的入射角和频率,以及流体的体积模量和密度对理想流体层动压力的影响。     

Study on SH-SAW in imperfectly bonded piezoelectric structures loaded with viscous liquid

We analytically investigate shear horizontal surface acoustic wave (SH-SAW) propagation in layered piezoelectric structures loaded with viscous liquid, which involves a thin piezoelectric layer imperfectly bonded to an unbounded elastic substrate. The coupling wave equations are obtained based on the linear piezoelectric theory. The governing equations are solved by means of the analytical method with consideration of electrically open and shorted cases, respectively. The dispersive relations are obtained, and the effects of the imperfect constant on the properties of waves are presented and discussed. From the numerical results, we can find that the phase velocity decreases with the increase of the interface parameter n, and for a specified viscosity, the attenuation increases with the interface parameter. The results show that the effects of the imperfect constant on the properties of SH-SAW are remarkable.     

Dynamic SIF of interface crack between two dissimilar viscoelastic bodies under impact loading*

In this paper, the transient dynamic stress intensity factor (SIF) is determined for an interface crack between two dissimilar half-infinite isotropic viscoelastic bodies under impact loading. An anti-plane step loading is assumed to act suddenly on the surface of interface crack of finite length. The stress field incurred near the crack tip is analyzed. The integral transformation method and singular integral equation approach are used to get the solution. By virtue of the integral transformation method, the viscoelastic mixed boundary problem is reduced to a set of dual integral equations of crack open displacement function in the transformation domain. The dual integral equations can be further transformed into the first kind of Cauchy-type singular integral equation (SIE) by introduction of crack dislocation density function. A piecewise continuous function approach is adopted to get the numerical solution of SIE. Finally, numerical inverse integral transformation is performed and the dynamic SIF in transformation domain is recovered to that in time domain. The dynamic SIF during a small time-interval is evaluated, and the effects of the viscoelastic material parameters on dynamic SIF are analyzed.     

Thermal loading of a thin layer with circular debonding over a substrate

By using techniques appropriate to mixed boundary value problems, this study addresses the determination of stress intensity factors for a circular interface debonding between a thin layer and a substrate subjected to nearly uniform temperature change. The solution method involves three-dimensional equilibrium equations of thermo-elasticity under axisymmetry conditions. The stress intensity factors are obtained by solving the resulting pair of coupled singular integral equations numerically. This revised version was published online in August 2006 with corrections to the Cover Date.