Poroelastodynamic potential method for transversely isotropic fluid-saturated poroelastic media

By introduction of two scalar potentials, an analytical method is developed for the solution of poroelastodynamic boundary value problems in transversely isotropic fluid-saturated poroelastic media. The governing equations of motion are considered in the framework of Biot's complete model without any assumption or simplification. As a case of application, solutions in three dimensions for a transversely isotropic fluid saturated porous half space loaded by an arbitrary distribution of time harmonic tractions at the free surface is derived. The free surface of the half space may be considered either permeable or impermeable. As a particular solution, Green's functions for uniform vertical and horizontal circular patch loads are presented as semi-infinite integrals which may be evaluated by means of an appropriate numerical method proposed. The accuracy of the solutions is verified both analytically and numerically against the preceding solutions. Some numerical results are also presented to clarify the influence of different degrees of anisotropy and frequency of excitation on the response of the medium.     

Displacement potentials for functionally graded piezoelectric solids

Two new displacement potential functions are introduced for the general solution of a three-dimensional piezoelasticity problem for functionally graded transversely isotropic piezoelectric solids. The material properties vary continuously along the axis of symmetry of the medium. The four coupled equilibrium equations in terms of displacements and electric potential are reduced to two decoupled sixth- and second-order linear partial differential equations for the potential functions. The obtained results are verified with two limiting cases: (i) a functionally graded transversely isotropic medium, and (ii) a homogeneous transversely isotropic piezoelectric solid. The simplified relations corresponding to the special case of similar variation of material properties are also given. Furthermore, the special cases of axisymmetric problems, exponentially graded piezoelectric media and transversely isotropic piezoelectric media with power law variation are discussed in detail.     

Green's function approach to study the propagation of SH‐wave in piezoelectric layer influenced by a point source

Green's function technique serves as a powerful tool to find the particle displacements due to SH‐wave propagation in layer of a shape different from the space between two parallel planes. Therefore, the present paper undertook to study the propagation of SH‐wave in a transversely isotropic piezoelectric layer under the influence of a point source and overlying a heterogeneous substrate using Green's function technique. The coupled electromechanical field equations are solved with the aid of Green's function technique. Expression for displacements in both layer and substrate, scalar potential and finally the dispersion relation is obtained analytically for the case when wave propagates along the direction of layering. Numerical computations are carried out and demonstrated with the aid of graphs for six different piezoelectric materials namely PZT‐5H ceramics, Barium titanate (BaTiO₃) ceramics, Silicon dioxide (SiO₂) glass, Borosilicate glass, Cobalt Iron Oxide (CoFe₂O₄), and Aluminum Nitride (AlN). The effects of heterogeneity, piezoelectric and dielectric constants on the dispersion curve are highlighted. Moreover, comparative study is carried out taking the phase velocity for different piezoelectric materials on one hand and isotropic case on the other. Dispersion relation is reduced to well‐known classical Love wave equation with a view to illuminate the authenticity of problem. Copyright © 2017 John Wiley & Sons, Ltd.     

具有横观各向同性球形压电夹杂的压电复合材料的有效性质

压电复合材料的有效性质是一个十分重要的问题.在分析了一无限大均匀线性非压电介质中含有一个横观各向同性球形压电夹杂问题的基础上,导出了夹杂内的约束应变场和约束电场的解析表达式,进而得到了含横观各向同性球形压电夹杂的压电复合材料的有效电弹模量的稀疏解.本文结果可以用于压电复合材料或智能材料和智能结构的设计与分析.     

横观各向同性弹性层点力解

本文根据弹性层状结构的传递矩阵法思想，由横观各向同性弹性力学基本方程，导出了含应力和位移两类变量的混合方程，利用Ｆｏｕｒｉｅｒ变换和文献［７］的位移函数通解，以及计算机代数软件，得到了横观各向同性层的点力解，这个点力解可直接退化到各同性情形的解．     

3D dynamic responses of a multi-layered transversely isotropic saturated half-space under concentrated forces and pore pressure

Dynamic Green's function plays an important role in the study of various wave radiation, scattering and soil-structure interaction problems. However, little research has been done on the response of transversely isotropic saturated layered media. In this paper, the 3D dynamic responses of a multi-layered transversely isotropic saturated half-space subjected to concentrated forces and pore pressure are investigated. First, utilizing Fourier expansion in circumferential direction accompanied by Hankel integral transform in radial direction, the wave equations for transversely isotropic saturated medium in cylindrical coordinate system are solved. Next, with the aid of the exact dynamic stiffness matrix for in-plane and out-of-plane motions, the solutions for multi-layered transversely isotropic saturated half-space under concentrated forces and pore pressure are obtained by direct stiffness method. A FORTRAN computer code is developed to achieve numerical evaluation of the proposed method, and its accuracy is validated through comparison with existing solutions that are special cases of the more general problems addressed. In addition, selected numerical results for a homogeneous and a layered material model are performed to illustrate the effects of material anisotropy, load frequency, drainage condition and layering on the dynamic responses. The presented solutions form a complete set of Green's functions for concentrated forces (including horizontal load in x(y)-direction, vertical load in z-direction) as well as pore pressure, which lays the foundation for further exploring wave propagation of complex local site in a layered transversely isotropic saturated half-space by using the BEMs.     

A general solution for three dimensional steady-state transversely isotropic hygro-thermo-magneto-piezoelectric media

In this article, we extract the general solution of three dimensional (3D) equations using potential theory method (PTM) for steady-state, transversely isotropic, hygro-thermo-magneto-piezoelectric media (HTMPM). The governing equations are simplified by introducing the displacement functions. A general solution is completely determined by advantage of the superposition principle and operator theory, which is connected in terms of two functions, fulfilling a second-order and twelfth-order homogeneous partial differential equation (PDE), separately. With the help of Almansi’s theorem, the general solution can be further shortened, which is stated by seven harmonic functions only. The acquired general solutions are straightforward structure and helpful in boundary value problems of HTMPM. Further, we apply the 3D fundamental solutions inside an infinite and on the surface of semi-infinite of a steady point heat source united with a steady point moisture source transversely isotropic HTMPM. Comprehensive and exact solutions are given in the form of elementary functions, which appear as a standard for various types of approximate solutions and numerical codes. Some numerical simulation is conducted based on the obtained general solutions.     

横观各向同性含液饱和多孔介质中应力波传播的特征分析

根据广义特征理论,对横观各向同性含液饱和多孔介质中应力波传播特性进行了特征分析．给出了特征曲面的微分方程以及沿次特征线的相容条件,得到了波阵面的解析表达式．详细地讨论了应力波在横观各向同性含液饱和多孔介质中传播时,其速度曲面和波阵面的形状及性质．分析结果亦表明,纯固体中应力波传播的特征方程,是含液饱和多孔介质中应力波特征方程的特例．     

横观各向同性饱和弹性多孔介质非轴对称动力响应

应用Fourier展开和Hankel变换求解了简谐激励下横观各向同性饱和弹性多孔介质的非轴对称Biot波动方程,得到了一般解。用一般解给出了多孔介质总应力分量的表达式。最后对求解横观各向同性饱和弹性多孔介质非轴对称动力响应边值问题的方法作了系统说明,并且给出了数值分析特例。     

Strip electro-mechanical yielding model for piezoelectric plate cut along two equal collinear cracks

The multiple-crack problems for piezoelectric ceramics till now have not yet address the crack opening arrest problem. The present work addresses this paucity. A 2-D strip-electro-mechanical yielding model is proposed for a transversely isotropic piezoelectric media weakened by two internal equal collinear straight cracks. The infinite boundary is prescribed with combined uniform constant in-plane mechanical and electrical loads. Developed mechanical and electric strip zones are arrested by prescribing over their rims uniform, normal, cohesive yield point stress and saturation limit electric displacement. Two cases are considered when saturation zone is bigger than developed yield zone and vice versa. Stroh formulation together with complex variable technique is employed to obtain the solution. Closed form expressions are derived for saturation zone length, yield zone length, crack opening displacement (COD), crack opening potential jump (COP) and energy release rate (ERR). An illustrative numerical study is prescribed to determine the effect of various parameters on the crack growth arrest and presented graphically. The results reveal that the model is capable of crack arrest under small-scale mechanical and electric yielding.     

横观各向同性饱和地基的三维动力响应

首先引入位移函数,将直角坐标系下横观各向同性饱和土Biot波动方程转化为2个解耦的六阶和二阶控制方程;然后基于双重Fourier变换,求解了Biot波动方程,得到以土骨架位移和孔隙水压力为基本未知量的积分形式的一般解,并用一般解给出了饱和土总应力分量的表达式．在此基础上系统研究了横观各向同性饱和半空间体的稳态动力响应问题,考虑表面排水和不排水两种情况,得到了半空间体在任意分布的表面谐振荷载作用下,表面位移的稳态动力响应,文末给出了算例．     

层状横观各向同性饱和地基上圆板的非轴对称振动

研究了层状横观各向同性饱和地基上弹性圆板的非轴对称振动问题．首先,通过方位角的Fourier变换,将圆柱坐标系下横观各向同性饱和土的三维动力方程转化为一阶常微分方程组,基于径向Hankel变换,建立问题的状态方程,求解状态方程后得到传递矩阵;其次,利用传递矩阵,结合层状饱和地基的边界条件、排水条件及层间接触和连续条件,给出了任意简谐激振力作用下层状横观各向同性饱和地基动力响应的通解;然后,按混合边值问题建立层状饱和地基上弹性圆板非轴对称振动的对偶积分方程,并将对偶积分方程化为易于数值计算的第二类Fredholm积分方程,并给出了算例．     

二维各向异性压电介质机电耦合场的基本解

本文研究各向异性压电介质的机电耦台问题．应用平面波分解法和留数定理，首次得到了线力和线电荷作用下一般二维各向异性压电介质机电耦合场的基本解．本文的解适用于平面问题、反平面问题以及平面和反平面相互耦合问题．作为特例，文中给出了横观各向同性压电介质的基本解．     

Analysis solution method for 3D planar crack problems of two-dimensional hexagonal quasicrystals with thermal effects

An analysis solution method (ASM) is proposed for analyzing arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystal (QC) media. The extended displacement discontinuity (EDD) boundary integral equations governing three-dimensional (3D) crack problems are transferred to simplified integral-differential forms by introducing some complex quantities. The proposed ASM is based on the analogy between these EDD boundary equations for 3D planar cracks problems of 2D hexagonal QCs and those in isotropic thermoelastic materials. Mixed model crack problems under combined normal, tangential and thermal loadings are considered in 2D hexagonal QC media. By virtue of ASM, the solutions to 3D planar crack problems under various types of loadings for 2D hexagonal QCs are formulated through comparison to the corresponding solutions of isotropic thermoelastic materials which have been studied intensively and extensively. As an application, analytical solutions of a penny-shaped crack subjected uniform distributed combined loadings are obtained. Especially, the analytical solutions to a penny-shaped crack subjected to the anti-symmetric uniform thermal loading are first derived for 2D hexagonal QCs. Numerical solutions obtained by EDD boundary element method provide a way to verify the validity of the presented formulation. The influences of phonon-phason coupling effect on fracture parameters of 2D hexagonal QCs are assessed.     

Transient response of a transversely isotropic multilayered half-space due to a vertical loading

This paper investigates the transient response of a transversely isotropic multilayered half-space under vertical loadings. With the aid of a Laplace–Hankel transform, the global stiffness matrix for a multilayered half-space is acquired by assembling the analytical layer-element of each layer medium. The solutions for the displacements in the time domain are obtained by using the global stiffness matrix equations and a numerical inversion procedure. The accuracy of the proposed method is verified through comparisons with existing solutions for displacements induced by a step and rectangular pulse loading. In addition, selected numerical results for displacements induced by the buried loading are presented to illustrate the effect of transient loading type and material anisotropy on the transient response.     

3D Green’s functions for a transversely isotropic thermoelastic cone

We use the compact harmonic general solutions of transversely isotropic thermoelastic materials to construct the three-dimensional Green’s functions of a steady point heat source on the apex of a transversely isotropic thermoelastic cone by three newly introduced harmonic functions. All components of thermoelastic field are expressed in terms of elementary functions and are convenient to use. When the apex angle 2α equals to π, the solution reduce to the important solution of semi-infinite body with a surface point heat source. Numerical results are given graphically by contours.     

Time-harmonic response of transversely isotropic and layered poroelastic half-spaces under general buried loads

A semi-analytical method is developed for solving the dynamic response of transversely isotropic, multilayered, and poroelastic half-spaces with different surface hydraulic conditions and subjected to time-harmonic vertical and horizontal loads buried in the layered half-space. The coupled governing equations of motion are presented in details in terms of the Biot's poroelastodynamic theory via the (u,p) formulation. The cylindrical system of vector functions is introduced to express the unknown primary quantities so that the coupled governing partial differential equations can be reduced and separated into two sets of first-order ordinary differential equations (i.e., the LM- and N-types). A recursive relation for the expansion coefficients among different layers is established by virtue of the stable and efficient dual variable and position method. Making use of the boundary and interface conditions, the fundamental solutions are obtained in terms of the vector-function system. The corresponding physical-domain solutions are then derived via an accurate semi-infinite integral algorithm. The developed fundamental solutions are carefully checked with existing solutions, and numerical examples are further presented to demonstrate the effect of material anisotropy, loading depth, material layering, and surface hydraulic condition on the dynamic response, which should be useful to design engineers. These solutions could be further served as benchmarks for future numerical methods.     

集中载荷作用下层合厚圆板的轴对称弯曲

从三维弹性力学基本方程出发,建立了横观各向同性层合圆板轴对称弯曲问题的状态方程,并将板面的集中载荷展成付里叶贝塞尔级数,从而给出问题的解析解,此解满足弹性力学全部方程,计及了所有独立的弹性常数,并满足层间连续性条件。     

Least-squares mixed finite elements with applications to anisotropic elasticity and viscoplasticity

The objective of this work is to discuss a least-squares finite element method with applications to physically nonlinear and anisotropic constitutive equations at small strains. The L₂-norm minimization of the residuals of the given first order system of differential equations leads to a functional, which is a two field formulation in the displacements and the stresses, see e.g. Cai & Starke [1]. These functionals provide the foundation for the formulations of the related least-squares mixed finite elements. A main focus of the presentation lies on the extension of plane elasticity to anisotropic or nonlinear material behavior. In this context transversely isotropic elasticity and viscoplasticity is considered. Finally a numerical example is presented. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)