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.     

Analytical layer-element solutions of Biot’s consolidation with anisotropic permeability and incompressible fluid and solid constituents

Based on the governing equations of 2D plane-strain Biot’s consolidation, the relationship between generalized displacements and stresses of a single soil layer with anisotropic permeability and incompressible fluid and solid constituents is described by an analytical layer-element, which is deduced in the Laplace–Fourier transform domain by using the eigenvalue approach. Taking the boundary conditions and the continuity of the soil layers into consideration, a global stiffness matrix is subsequently assembled and solved. As to the 3D case, the same derivation is employed after the application of a decoupling transformation. The actual solutions in the physical domain can further be acquired by inverting the Laplace–Fourier transform. Finally, numerical examples are carried out to verify the presented theory and discuss the influence of the anisotropic permeability on the consolidation behavior.     

3D dynamic Green’s functions in a multilayered poroelastic half-space

The complete 3D dynamic Green’s functions in the multilayered poroelastic media are presented in this study. A method of potentials in cylindrical coordinate system is applied first to decouple the Biot’s wave equations into four scalar Helmholtz equations, and then, general solutions to 3D wave propagation problems are obtained. After that, a three vector base and the propagator matrix method are introduced to treat 3D wave propagation problems in the stratified poroelastic half-space disturbed by buried sources. It is known that the original propagator algorithm has the loss-of-precision problem when the waves become evanescent. At present, an orthogonalization procedure is inserted into the matrix propagation loop to avoid the numerical difficulty of the original propagator algorithm. At last, the validity of the present approach for accurate and efficient calculating 3D dynamic Green’s functions of a multilayered poroelastic half-space is confirmed by comparing the numerical results with the known exact analytical solutions of a uniform poroelastic half-space.     

Fast integral equation solver for Maxwell’s equations in layered media with FMM for Bessel functions

The paper presents a new fast integral equation solver for Maxwell’s equations in 3-D layered media. First, the spectral domain dyadic Green’s function is derived, and the 0-th and the 1-st order Hankel transforms or Sommerfeld-type integrals are used to recover all components of the dyadic Green’s function in real space. The Hankel transforms are performed with the adaptive generalized Gaussian quadrature points and window functions to minimize the computational cost. Subsequently, a fast integral equation solver with O(N _z ² N _x N _y log(N _x N _y )) in layered media is developed by rewriting the layered media integral operator in terms of Hankel transforms and using the new fast multipole method for the n-th order Bessel function in 2-D. Computational cost and parallel efficiency of the new algorithm are presented.     

Vibration of a pre-stressed plate on a transversely isotropic multilayered half-plane due to a moving load

This paper presents an analytical research on the dynamic interaction problem between a pre-stressed plate and a transversely isotropic multilayered half-plane subjected to a moving load. The pre-stressed plate is governed by the Kirchhoff plate theory, and the transversely isotropic multilayered half-plane is solved by the analytical layer-element method. Combining the frictionless contact and displacement compatibility conditions between the plate and the soil, the contact stress and the deflection of plate in the Fourier transform domain are derived. With the aid of the inverse Fourier transform, the actual solutions can be further achieved. Numerical examples are given to illustrate the influence of load speed, the rigidity of plate, the axial force applied on the plate and the stratified character of the soil.     

Fractional non-axisymmetric consolidation of stratified cross-anisotropic visco-poroelastic media

This paper investigates the fractional non-axisymmetric consolidation of stratified cross-anisotropic visco-poroelastic media. Firstly, based on the Laplace–Hankel transform, the ordinary differential equations of cross-anisotropic poroelastic media are derived, which are extended to those of fractional visco-poroelastic media by introducing the fractional Merchant model and the elastic-viscoelastic correspondence principle. Then, the extended equations are solved by the extended precise integration method (PIM), and the proposed theory is compared with the results reported in the references. Finally, selected numerical examples are performed to investigate the effects of the order of fractional derivative, viscoelasticity, cross-anisotropy and stratification on the non-axisymmetric consolidation.     

State vector solutions for nonaxisymmetric problem of multilayered half space piezoelectric medium

A state space formulation is established for the nonaxisymmetric space problem of transversely isotropic piezoelectric media in a system of cylindrical coordinate by introducing the state vector. Using the Hankel transform and the Fourier series, the state vector equations are transformed into ordinary differential equations. By the use of the matrix methods, the analytical solutions of a single piezoelectric layer are presented in the form of the product of initial state variables and transfer matrix. The applications of state vector solutions are discussed. An analytical solution for a semiinfinite piezoelectric medium subjected to the vertical point forceP _z, horizontal point forceP _x along x-direction and point electric charge Q at the origin of the surface is presented. According to the continuity conditions at the interfaces, the general solution formulation forN-layered transversely isotropic piezoelectric media is given. Project supported by the National Natural Science Foundation of China (Grant No. 59648001).     

Three-dimensional dynamic Green’s functions in transversely isotropic tri-materials

An analytical derivation of the elastodynamic fundamental solutions for a transversely isotropic tri-material full-space is presented by means of a complete representation using two displacement potentials. The complete set of three-dimensional point-load, patch-load, and ring-load Green’s functions for stresses and displacements are given, for the first time, in the complex-plane line-integral representations. The formulation includes a complete set of transformed stress-potential and displacement-potential relations in the framework of Fourier expansions and Hankel integral transforms, that is useful in a variety of elastodynamic as well as elastostatic problems. For the numerical computation of the integrals, a robust and effective methodology is laid out. Selected numerical results for point-load and patch-load Green’s functions are presented to portray the dependence of the response on layering, the frequency of excitation, and type of loading.     

A new off-step high order approximation for the solution of three-space dimensional nonlinear wave equations

In this paper, we propose a new high accuracy numerical method of O(k² + k²h² + h⁴) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form u_tt = A(x,y,z,t)u_xx + B(x,y,z,t)u_yy + C(x,y,z,t)u_zz + g(x,y,z,t,u,u_x,u_y,u_z,u_t), 0 < x,y,z < 1,t > 0 subject to given appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in  and . We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method.     

Analytical and numerical aspects of a generalization of the complementary error function

In this paper we discuss analytical and numerical properties of the function , with α,β,Rz>0, which can be viewed as a generalization of the complementary error function, and in fact also as a generalization of the Kummer U-function. The function V_ν,μ(α, β, z) is used for certain values of the parameters as an approximate in a singular perturbation problem. We consider the relation with other special functions and give asymptotic expansions as well as recurrence relations. Several methods for its numerical evaluation and examples are given.     

Asymptotics of cylindrical functions in the complex domain: II

We obtain asymptotic formulas uniform with respect to the index p > 0 for the Hankel functions H _p ^(j) (z) (j = 1, 2) for large |z| in the complex domain. These formulas generalize those known for the real argument.     

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.     

Asymptotics of cylindrical functions in the complex domain: I

We obtain asymptotic formulas uniform with respect to the index p > 0 for the Hankel functions H _p ^(j)(z) (j = 1, 2) for large |z| in the complex domain. These formulas generalize those well known for the real argument.     

Landau’s theorem for functions with logharmonic Laplacian

In this paper, we show the existence of Landau constant for functions with logharmonic Laplacian of the form F(z) = ∣z∣²L(z) + K(z), ∣z∣ < 1, where L is logharmonic and K is harmonic. Moreover, the problem of minimizing the area is solved     

Linear systems and determinantal random point fields

In random matrix theory, determinantal random point fields describe the distribution of eigenvalues of self-adjoint matrices from the generalized unitary ensemble. This paper considers symmetric Hamiltonian systems and determines the properties of kernels and associated determinantal random point fields that arise from them; this extends work of Tracy and Widom. The inverse spectral problem for self-adjoint Hankel operators gives sufficient conditions for a self-adjoint operator to be the Hankel operator on L²(0,∞) from a linear system in continuous time; thus this paper expresses certain kernels as squares of Hankel operators. For suitable linear systems (−A,B,C) with one-dimensional input and output spaces, there exists a Hankel operator Γ with kernel ?(x)(s+t)=Ce−(2x+s+t)AB such that gx(z)=det(I+(z−1)ΓΓ†) is the generating function of a determinantal random point field on (0,∞). The inverse scattering transform for the Zakharov-Shabat system involves a Gelfand-Levitan integral equation such that the trace of the diagonal of the solution gives . When A?0 is a finite matrix and B=C†, there exists a determinantal random point field such that the largest point has a generalised logistic distribution.     

A stable algorithm for numerical evaluation of Hankel transforms using Haar wavelets

The purpose of the paper is to propose a stable algorithm for the numerical evaluation of the Hankel transform F _n (y) of order n of a function f(x) using Haar wavelets. The integrand \(\sqrt x f(x)\) is replaced by its wavelet decomposition. Thus representing F _n (y) as a series with coefficients depending strongly on the local behavior of the function \(\sqrt x f(x)\), thereby getting an efficient and stable algorithm for their numerical evaluation. Numerical evaluations of test functions with known analytical Hankel transforms illustrate the proposed algorithm.     

Symbolic and numerical computation on Bessel functions of complex argument and large magnitude

The Lanczos τ-method, with perturbations proportional to Faber polynomials, is employed to approximate the Bessel functions of the first kind J_v(z) and the second kind Y_v(z), the Hankel functions of the first kind H_v⁽¹⁾(z) and the second kind H_v⁽²⁾(z) of integer order v for specific outer regions of the complex plane, i.e. ¦z¦ ⩾ R for some R. The scaled symbolic representation of the Faber polynomials and the appropriate automated τ-method approximation are introduced. Both symbolic and numerical computation are discussed. In addition, numerical experiments are employed to test the proposed τ-method. Computed accuracy for J₀(z) and Y₀(z) for ¦z¦ ⩾ 8 are presented. The results are compared with those obtained from the truncated Chebyshev series approximations and with those of the τ-method approximations on the inner disk ¦z¦ ⩽ 8. Some concluding remarks and suggestions on future research are given.     

Meromorphic Functions That Share Fixed-Points

In this paper, we prove the following result: Let f(z) and g(z) be two nonconstant meromorphic(entire) functions, n ≥ 11(n ≥ 6) a positive integer. If fⁿ(z)f′(z) and gⁿ(z)g′(z) have the same fixed-points, then either f(z) = c₁e^cz2, g(z) = c₂e^− cz2, where c₁, c₂, and c are three constants satisfying 4(c₁c₂)^n + 1c² = −1, or f(z) ≡ tg(z) for a constant t such that t^n + 1 = 1.     

Algorithms for the computation of Hankel functions of complex order

Hankel functions of complex order and real argument arise in the study of wave propagation and many other applications. Hankel functions are computed using, for example, Chebyshev expansions, recursion relations and numerical integration of the integral representation. In practice, approximation of these functions is required when the order and the argumentz are large.When andz are large, the Chebyshev series expansion of the Hankel function is of limited use. The situation is remedied by the use of appropriate asymptotic expansions. These expansions are normally expressed in terms of coefficients which are defined recursively involving derivatives and integrals of polynomials. The applicability of these expansions in both numerical and symbolic software is discussed with illustrative examples.