Compressional waves in fluid-saturated porous solid containing a penny-shaped crack

Diffraction of normal compression waves by a penny-shaped crack in a fluid-saturated porous medium is investigated. Two wave types are considered, namely, compressional wave of the first kind, and the second kind. The former, also known as fast wave, propagates primarily through the solid, whereas the latter or slow wave, propagates mainly in the fluid. Each wave propagates in the medium along with induced wave of the same type in the companion constituent of the material. Application of Biot’s theory in conjunction with integral transform technique reduces the problem to a mixed boundary-value problem whose solution is in turn governed by a Fredholm integral equation of the second kind. Near-field and far-field solutions are obtained in terms of the dynamic stress-intensity factor and the scattering cross section, respectively. They are of particular importance to the linear elastic fracture mechanics (LEFM) and in the scattering theory of elastic waves. The mode I stress-intensity factors are computed numerically for a set of selected material property values, and shown graphically for various mass density and viscosity-to-permeability ratios. The obtained results reveal significant impact of the presence of pore fluid upon the stress-intensity factors, both magnitudes and frequencies at their peak values. The influence of the fluid is also observed from the calculated scattering cross sections of the scattered far-field. Accuracy of the present solution procedure is verified by comparing the numerical results with existing results in the limiting case of dry elastic materials.     

Elastic wave scattering from a penny-shaped interface crack in coated materials

This paper is concerned with the elastic wave scattering induced by a penny-shaped interface crack in coated materials. Using the integral transform, the problem of wave scattering is reduced to a set of singular integral equations in matrix form. The singular integral equations are solved by the asymptotic analysis and contour integral technique, and the expressions for the stress and displacement as well as the dynamic stress intensity factors (SIFs) are obtained. Using numerical analysis, this approach is verified by the finite element method (FEM), and the numerical results agree well with the theoretical results. For various crack sizes and material combinations, the relations between the SIFs and the incident frequency are analyzed, and the amplitudes of the crack opening displacements (CODs) are plotted versus incident wavenumber. The investigation provides a theoretical basis for the dynamic failure analysis and nondestructive evaluation of coated materials.     

Analysis of cracked magnetoelectroelastic composites under time-harmonic loading

This paper presents a numerical model for the analysis of cracked magnetoelectroelastic materials subjected to in-plane mechanical, electric and magnetic dynamic time-harmonic loading. A traction boundary integral equation formulation is applied to solve the problem in combination with recently obtained time-harmonic Green’s functions (Rojas-Diaz et al., 2008). The hypersingular boundary integral equations appearing in the formulation are first regularized via a simple change of variables that permits to isolate the singularities. Relevant fracture parameters, namely stress intensity factors, electric displacement intensity factor and magnetic induction intensity factor are directly evaluated as functions of the computed nodal opening displacements and the electric and magnetic potentials jumps across the crack faces. The method is checked by comparing numerical results against existing solutions for piezoelectric solids. Finally, numerical results for scattering of plane waves in a magnetoelectroelastic material by different crack configurations are presented for the first time. The obtained results are analyzed to evaluate the dependence of the fracture parameters on the coupled magnetoelectromechanical load, the crack geometry and the characteristics of the incident wave motion.     

Coupled hydro-mechanical effects in a poro-hyperelastic material

Fluid-saturated materials are encountered in several areas of engineering and biological applications. Geologic media saturated with water, oil and gas and biological materials such as bone saturated with synovial fluid, soft tissues containing blood and plasma and synthetic materials impregnated with energy absorbing fluids are some examples. In many instances such materials can be examined quite successfully by appeal to classical theories of poroelasticity where the skeletal deformations can be modelled as linear elastic. In the case of soft biological tissues and even highly compressible organic geological materials, the porous skeleton can experience large strains and, unlike rubberlike materials, the fluid plays an important role in maintaining the large strain capability of the material. In some instances, the removal of the fluid can render the geological or biological material void of any hyperelastic effects. While the fluid component can be present at various scales and forms, a useful first approximation would be to treat the material as hyperelastic where the fabric can experience large strains consistent with a hyperelastic material and an independent scalar pressure describes the pore fluid response. The flow of fluid within the porous skeleton is defined by Darcy's law for an isotropic material, which is formulated in terms of the relative velocity between the pore fluid and the porous skeleton. It is assumed that the form of Darcy's law remains unchanged during the large strain behaviour. This approach basically extends Biot's theory of classical poroelasticity to include finite deformations. The developments are used to examine the poro-hyperelastic behaviour of certain one-dimensional problems.     

SH波在压电材料条中垂直界面裂纹处的散射

研究了SH波在压电材料条中裂纹处的散射.压电材料条两侧涂有相同梯度参数的两个半无限大功能梯度材料,裂纹垂直于界面.通过Fourier变换,利用边界条件把问题转化为柯西核奇异积分方程,然后利用Chebyshev多项式对奇异积分方程进行数值求解.通过数值计算,分析讨论了压电条的几何参数和SH波频率对标准动应力强度因子的影响.     

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

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

Scattering of plane monochromatic waves from a heterogeneous inclusion of arbitrary shape in a poroelastic medium: An efficient numerical solution

Scattering of plane longitudinal monochromatic waves from a heterogeneous inclusion of arbitrary shape in an infinite poroelastic medium is considered. Wave propagation in the medium is described by Biot’s equations of poroelasticity. The scattering problem is formulated in terms of the volume integral equations for displacements of the solid skeleton and fluid pressure in the pore space in the region occupied by the inclusion. An efficient numerical method is applied to solve these equations. In the method, Gaussian approximating functions are used for discretization of the problem. For regular node grids, the matrix of the discretized problem has Toeplitz’s properties, and the Fast Fourier Transform technique can be used for the calculation of matrix–vector products. The latter accelerates substantially the process of iterative solution of the discretized problem. For material parameters of typical sedimentary rocks, the system of differential equations of poroelasticity contains a differential operator with a small parameter. As the result, the wave field in the inclusion region is split up into a slowly changing part, and boundary layer functions concentrated near the inclusion interface. The method of matched asymptotic expansions is used for the numerical solution in this case. For a spherical inclusion, the results of the numerical and matched asymptotic expansion methods are compared with a semi-analytical series solution. For a non-spherical heterogeneous inclusion, an example of the numerical solution is presented.     

The constitutive equations of finite strain poroelasticity in the light of a micro-macro approach

After recalling the constitutive equations of finite strain poroelasticity formulated at the macroscopic level, we adopt a microscopic point of view which consists of describing the fluid-saturated porous medium at a space scale on which the fluid and solid phases are geometrically distinct. The constitutive equations of poroelasticity are recovered from the analysis conducted on a representative elementary volume of porous material open to fluid mass exchange. The procedure relies upon the solution of a boundary value problem defined on the solid domain of the representative volume undergoing large elastic strains. The macroscopic potential, computed as the integral of the free energy density over the solid domain, is shown to depend on the macroscopic deformation gradient and the porous space volume as relevant variables. The corresponding stress-type variables obtained through the differentiation of this potential turn out to be the macroscopic Boussinesq stress tensor and the pore pressure. Furthermore, such a procedure makes it possible to establish the necessary and sufficient conditions to ensure the validity of an ‘effective stress’ formulation of the constitutive equations of finite strain poroelasticity. Such conditions are notably satisfied in the important case of an incompressible solid matrix.     

Wave scattering from an interface crack in multilayered piezoelectric plate

This paper focuses on the theoretical basis for the study of wave scattering from an interface crack in multilayered piezoelectric media. The materials are taken to be anisotropic with arbitrary symmetry. Based on the Fourier transform technique together with the aid of the stiffness matrix approach, the boundary value problem of wave scattering is reduced to solving a system of Cauchy-type singular equations. The intensity factors and crack opening displacements are defined in terms of the solutions of the corresponding integral equations for any incident frequencies and incident angles. Numerical results are presented. The effects of incident frequencies and crack location on both the major and coupling intensity factors are illustrated. The influence of the piezoelectricity is also shown.     

Dynamic bulk and shear moduli due to grain-scale local fluid flow in fluid-saturated cracked poroelastic rocks: Theoretical model

Grain-scale local fluid flow is an important loss mechanism for attenuating waves in cracked fluid-saturated poroelastic rocks. In this study, a dynamic elastic modulus model is developed to quantify local flow effect on wave attenuation and velocity dispersion in porous isotropic rocks. The Eshelby transform technique, inclusion-based effective medium model (the Mori–Tanaka scheme), fluid dynamics and mass conservation principle are combined to analyze pore-fluid pressure relaxation and its influences on overall elastic properties. The derivation gives fully analytic, frequency-dependent effective bulk and shear moduli of a fluid-saturated porous rock. It is shown that the derived bulk and shear moduli rigorously satisfy the Biot-Gassmann relationship of poroelasticity in the low-frequency limit, while they are consistent with isolated-pore effective medium theory in the high-frequency limit. In particular, a simplified model is proposed to quantify the squirt-flow dispersion for frequencies lower than stiff-pore relaxation frequency. The main advantage of the proposed model over previous models is its ability to predict the dispersion due to squirt flow between pores and cracks with distributed aspect ratio instead of flow in a simply conceptual double-porosity structure. Independent input parameters include pore aspect ratio distribution, fluid bulk modulus and viscosity, and bulk and shear moduli of the solid grain. Physical assumptions made in this model include (1) pores are inter-connected and (2) crack thickness is smaller than the viscous skin depth. This study is restricted to linear elastic, well-consolidated granular rocks.     

A 2.5-D dynamic model for a saturated porous medium. Part II: Boundary element method

The 3-D boundary integral equation is derived in terms of the reciprocal work theorem and used along with the 2.5-D Green’s function developed in Part I [Lu, J.F., Jeng, D.S., Williams, S., submitted for publication. A 2.5-D dynamic model for a saturated porous medium: Part I. Green’s function. Int. J. Solids Struct.] to develop the 2.5-D boundary integral equation for a saturated porous medium. The 2.5-D boundary integral equations for the wave scattering problem and the moving load problem are established. The Cauchy type singularity of the 2.5-D boundary integral equation is eliminated through introduction of an auxiliary problem and the treatment of the weakly singular kernel is also addressed. Discretisation of the 2.5-D boundary integral equation is achieved using boundary iso-parametric elements. The discrete wavenumber domain solution is obtained via the 2.5-D boundary element method, and the space domain solution is recovered using the inverse Fourier transform. To validate the new methodology, numerical results of this paper are compared with those obtained using an analytical approach; also, some numerical results and corresponding analysis are presented.     

薄膜涂层材料中币形界面裂纹的弹性波散射

研究了薄膜涂层材料中币形界面裂纹的弹性波散射问题,建立了含有币形界面裂纹的覆层半空间模型,采用Hankel积分变换,将裂纹对弹性波散射的问题转化为求解矩阵形式的奇异积分方程。结合渐近分析和围道积分技术得到积分方程的解,进一步推导了散射波的应力场和位移场,以及动应力强度因子的理论计算公式。在数值算例中,分析了不同材料组合和裂纹尺寸情况下动应力强度因子与入射波频率的关系,并给出了裂纹张开位移的结果。为薄膜涂层材料的动态破坏分析提供了一定的理论基础。     

A boundary element solution to scattering of elastic SH wave by arbitrarily shaped inclusions embedded in an infinite anisotropic medium

Scattering problems for inhomogeneous bodies are investigated by the integral equation method. The boundary integral equation (BIE) for the scattered displacement field associated with finite inhomogeneities in an anisotropic medium are derived with the help of the generalized Green's identity. The discretization of BIE is based upon the constant element, linear element and quadratic element. Several numerical examples for calculating the scattering displacement, stress and scattering cross section from a cylinder, an interface crack, and two elliptic cylinders are given. Results show that the present method can be advantageously applied to a wide range of scattering problems of elastic waves.     

Propagation of axial shear magneto–electro-elastic waves in piezoelectric–piezomagnetic composites with randomly distributed cylindrical inhomogeneities

The integral equations of the scattering problem for piezoelectric–piezomagnetic composites with an inhomogeneity are derived. In the long-wave limit, the solutions of these integral equations for the composites containing a single inhomogeneous fiber are solved in close forms. The total scattering cross-section for the one-fiber composites is also obtained. By the so-called effective field method, the multi-fiber scattering problem is simplified to the one-fiber scattering problem, and the analytical expressions of magneto–electro-elastic fields for the multi-fiber composites are obtained in the long-wave limit. These solved magneto–electro-elastic fields are then used to solve the expressions of the static effective moduli, effective wave velocity and attenuation factor of piezoelectric–piezomagnetic composites with randomly distributed cylindrical inhomogeneities. Through numerical examples, it concludes that, if the random set of fiber cross-sections is homogeneous and isotropic, the effective field method is coincident with the Mori–Tanaka mean field method when the static effective moduli of piezoelectric–piezomagnetic composites are looked for. Moreover, the rules of the effective wave velocity versus the volume fraction of fibers are investigated for specific materials.     

Flexural waves scattering by a through crack in a magnetically saturated Mindlin plate under a uniform magnetic field

This paper deals with the scattering of time harmonic flexural waves by a through crack in a magnetically saturated plate under a uniform magnetic field normal to the plate surfaces. The analysis is based on Mindlin's plate theory of magneto-elastic interactions under a magnetic field. An incident wave giving rise to moments symmetric about the crack plane is applied. Fourier transforms are used to reduce the mixed boundary value problem to one involving the numerical solution of a Fredholm integral equation. The dynamic moment intensity factor versus frequency is computed and the influence of the magnetic field on the normalized values is displayed graphically.     

Two non-elliptical inhomogeneities with internal uniform stresses interacting with a mode-III crack

We use the method of Green's functions to analyze an inverse problem in which we aim to identify the shapes of two non-elliptical elastic inhomogeneities, embedded in an infinite matrix subjected to uniform remote stress, which enclose uniform stress distributions despite their interaction with a finite mode-III crack. The problem is reduced to an equivalent Cauchy singular integral equation, which is solved numerically using the Gauss–Chebyshev integration formula. The shapes of the two inhomogeneities and the corresponding location of the crack can then be determined by identifying a conformal mapping composed in part of a real density function obtained from the solution of the aforementioned singular integral equation. Several examples are given to demonstrate the solution.     

SCATTERING OF ANTI-PLANE SHEAR WAVES IN A FUNCTIONALLY GRADED MATERIAL STRIP WITH AN OFF-CENTER VERTICAL CRACK

The scattering problem of anti-plane shear waves in a functionally graded material strip with an off-center crack is investigated by use of Schmidt method. The crack is vertically to the edge of the strip. By using the Fourier transform, the problem can be solved with the help of a pair of dual integral equations that the unknown variable is the jump of the displacement across the crack surfaces. To solve the dual integral equations, the jump of the displacement across the crack surfaces was expanded in a series of Jacobi polynomials. Numerical examples were provided to show the effects of the parameter describing the functionally graded materials, the position of the crack and the frequency of the incident waves upon the stress intensity factors of the crack.     

功能梯度压电带中裂纹对SH波的散射

研究功能梯度压电带中裂纹对SH波的散射问题,为了便于分析,材料性质假定为指数模型,并假设裂纹面上的边界条件为电渗透型的.根据压电理论得到压电体的状态方程,利用Fourier积分变换,问题转化为对偶积分方程的求解.用Copson方法求解积分方程.求得了裂纹尖端动应力强度因子、电位移强度因子的解析表达式,最后数值结果显示了标准动应力强度因子与入射波数、材料参数、带宽、波数以及入射角之间的关系.     

含直对称裂纹的一维六方压电准晶对SH波的散射Scattering of SH wave by a linear symmetric crack at one-dimensional piezoelectric hexagonal quasicrystals

利用积分变换技术,结合Copson方法,研究了含直线型对称裂纹的一维六方压电准晶对SH波的散射问题。通过求解对偶积分方程,得到声子场、相位子场应力、位移及电场电位移分量的解析解。定义了裂纹尖端应力强度因子及电位移强度因子,给出了电非渗透性条件下应力强度因子及电位移强度因子的解析解。此研究结果对压电准晶材料的工程应用有一定的理论价值。     

Application of boundary integral equations for analyzing the dynamics of elastic,viscoelastic, and poroelastic bodies

Two approaches (classical and nonclassical) of the boundary integral equation method for solving three-dimensional dynamical boundary value problems of elasticity, viscoelasticity, and poroelasticity are considered. The boundary integral equation model is used for porous materials. The Kelvin–Voigt model and the weakly singular hereditary Abel kernel are used to describe the viscoelastic properties. Both approaches permit solving the dynamic problems exactly not only in the isotropic but also in the anisotropic case. The boundary integral equation solution scheme is constructed on the basis of the boundary element technique. The numerical results obtained by the classical and nonclassical approaches are compared.