Derivatives of the three-dimensional Green’s functions for anisotropic materials

A concise formulation is presented for the derivatives of Green’s functions of three-dimensional generally anisotropic elastic materials. Direct calculation for derivatives of the Green’s function on the Cartesian coordinate system is a common practice, which, however, usually leads to a complicated course. In this paper the Green’s function derived by Ting and Lee [Ting, T.C.T., Lee, V.G., 1997. The three-dimensional elastostatic Green’s function for general anisotropic linear elastic solids. The Quarterly Journal of Mechanics and Applied Mathematics 50 (3) 407–426] is extended to obtain the derivatives. Using a spherical coordinate system, the Green’s function can be shown as the composition of two independent functions, one depends only on the radial distance of the field point to the origin and the other is in spherical angles. The method of derivation is based on the total differential scheme and then takes its partial differentiation accordingly. With the application of Cauchy residue theorem, the contour integral can be evaluated in terms of the Stroh eigenvalues of a sextic equation. For the degenerate case, evaluation of residues at multiple poles is also given. Applications of the present result are made to examine the Green’s functions and stress components for isotropic and transversely isotropic materials. The results are in exact agreement with existing solutions.     

Three-dimensional Green’s functions for composite laminates

The three-dimensional Green’s functions due to a point force in composite laminates are solved by using generalized Stroh formalism and two-dimensional Fourier transforms. Each layer of the composite is generally anisotropic and linearly elastic. The interfaces between different layers are parallel to the top and bottom surfaces of the composite and are perfectly bonded. The Green’s functions of point forces applied at the free surface, interface, and in the interior of a layer are derived in the Fourier transformed domain respectively. The surfaces are imposed by a proportional spring-type boundary condition. The spring-type condition may be reduced to traction-free, displacement-fixed, and mirror-symmetric conditions. Numerical examples are given to demonstrate the validity and elegance of the present formulation of three-dimensional point-force Green’s functions for composite laminates.     

Time-harmonic Green’s functions for anisotropic magnetoelectroelasticity

Two-dimensional (2-D) and three-dimensional (3-D) time-harmonic Green’s functions for linear magnetoelectroelastic solids are derived in this paper by means of Radon-transform. Displacement field and electric and magnetic potentials in a fully anisotropic magnetoelectroelastic infinite solid due to a time-harmonic point force, point charge and magnetic monopole are obtained in form of line integrals over a unit circle in 2-D case and surface integrals over a unit sphere in 3-D case. This dynamic fundamental solution is then split into the sum of regular dynamic plus singular terms. The singular terms coincide with the Green’s functions for the static problem and may be further reduced to closed form expressions. The proposed Green’s functions can be used in the corresponding boundary element method (BEM) formulation.     

Explicit expression of derivatives of elastic Green’s functions for general anisotropic materials

The analytical expressions of Green’s function and their derivatives for three-dimensional anisotropic materials are presented here. By following the Fourier integral solutions developed by Barnett [Phys. Stat. Sol. (b) 49 (1972) 741], we characterize the contour integral formulations for the derivatives into three types of integrals H, M, and N. With Cauchy’s residues theorem and the roots of a sextic equation from Stroh eigenrelation, these integrals can be solved explicitly in terms of the Stroh eigenvalues P_i (i=1,2,3) on the oblique plane whose normal is the position vector. The results of Green’s functions and stress distributions for a transversely isotropic material are discussed in this paper.     

Dynamic Green’s Functions for an Anisotropic Multilayered Poroelastic Half-Space

The dynamic responses of an anisotropic multilayered poroelastic half-space to a point load or a fluid source are studied based on Stroh formalism and Fourier transforms. Taking the boundary conditions and the continuity of the materials into consideration, the three-dimensional Green’s functions of generalized concentrated forces (force and fluid source) applied at the free surface, interface and in the interior of a layer are derived in the Fourier transformed domain, respectively. The actual solutions in the frequency domain can further be acquired by inverting the Fourier transform. Finally, numerical examples are carried out to verify the presented theory and discuss the Green’s fields due to three cases of a concentrated force or a fluid source applied at three different locations for an anisotropic multilayered poroelastic half-space.

     

Superposing scheme for the three-dimensional Green’s functions of an anisotropic half-space

This paper presents a method of superposition for the half-space Green’s functions of a generally anisotropic material subjected to an interior point loading. The mathematical concept is based on the addition of a complementary term to the Green’s function in an anisotropic infinite domain. With the two-dimensional Fourier transformation, the complementary term is derived by solving the generalized Stroh eigenrelation and satisfying the boundary conditions on the free surface with the use of Green’s functions in the full-space case. The inverse Fourier transform leads to the contour integrals, which can be evaluated with the application of Cauchy residue theorem. Application of the present results is made to obtain analytical expression for the orthotropic materials which were not reported previously. The closed-form solutions for the transversely isotropic and isotropic materials derived directly from the solutions as being a special case are also given in this paper.     

High-order derivatives of Green’s functions in magneto-electro-elastic materials

Three-dimensional Green’s functions and their arbitrary order derivatives in general anisotropic magneto-electro-elastic materials are derived by using Fourier transform. They are analytical solutions expressed in line integral forms, and can be evaluated by a standard numerical integration method. With this method, we can obtain results with high accuracy. Besides, a numerical finite difference method is also given to evaluate the second-order derivatives quickly. When setting the appropriate material coefficients to zero, the piezoelectric, piezomagnetic, and purely anisotropic elastic Green’s functions and their derivatives can all be obtained from the current solutions.     

Wave propagation in two-dimensional anisotropic acoustic metamaterials of K4 topology

An acoustic metamaterial is envisaged as a synthesised phononic material the mechanical behaviour of which is determined by its unit cell. The present study investigates one aspect of mechanical behaviour, namely the band structure, in two-dimensional (2D) anisotropic acoustic metamaterials encompassing locally resonant mass-in-mass units connected by massless springs in a K4 topology. The 2D lattice problem is formulated in the direct space (r-space) and the equations of motion are derived using the principle of least action (Hamilton’s principle). Only proportional anisotropy and attenuation-free shock wave propagation have been considered. Floquet–Bloch’s principle is applied, therefore a generic unit cell is studied. The unit cell can represent the entire lattice regardless of its position. It is transformed from the direct lattice in r-space onto its reciprocal lattice conjugate in Fourier space (k-space) and point symmetry operations are applied to Wigner–Seitz primitive cell to derive the first irreducible Brillouin Zone (BZ). The edges of the first irreducible Brillouin Zone in the k-space have then been traversed to generate the full band structure. It was found that the phenomenon of frequency filtering exists and the pass and stop bands are extracted. A follow-up parametric study appreciated the degree and direction of influence of each parameter on the band structure.     

Stroh formalism in evaluation of 3D Green’s function in thermomagnetoelectroelastic anisotropic medium

The paper presents studies on the Green’s function for thermomagnetoelectroelastic medium and its reduction to the contour integral. Based on the previous studies the thermomagnetoelectroelastic Green’s function is presented as a surface integral over a half-sphere. The latter is then reduced to the double integral, which inner integral is evaluated explicitly using the complex variable calculus and the Stroh formalism. Thus, the Green’s function is reduced to the contour integral. Since the latter is evaluated over the period of the integrand, the paper proposes to use trapezoid rule for its numerical evaluation with exponential convergence. Several numerical examples are presented, which shows efficiency of the proposed approach for evaluation of Green’s function in thermomagnetoelectroelastic anisotropic solids.     

Steady and transient Green’s functions for anisotropic conduction in an exponentially graded solid

The problem of the determination of Green’s function in conduction for a rectilinearly anisotropic solid with an exponential grading along a certain direction is studied. Domains of an unbounded space and a half-space, either three-dimensional or two-dimensional, are considered. Along the boundary of the domain, homogeneous boundary conditions of the first and second kinds are imposed. We find interestingly that, under this specific type of grading, the Green’s functions permit an algebraic decomposition, which will in turn greatly simplify the formulation. The method of Fourier transform is employed for the Green’s function for a half-space or a half-plane. Although the derivation process is quite tedious, we show analytically that the inverse transform can be found exactly and their resulting expressions are surprisingly neat and compact. In addition, both steady-state and transient-state field solutions are considered. By taking Laplace transform with respect to the time variable, we show that the mathematical frameworks for the steady-state and transient-state Green’s functions are entirely analogous. Thereby, the transient-state Green’s function is readily obtained by taking Laplace inverse transform back to the time domain. These derived fundamental solutions will serve as benchmark results for modeling some inhomogeneous materials. In the absence of grading term, we have verified analytically that our solutions agree exactly with previously known Green’s functions for homogeneous media.     

3D frictional contact of anisotropic solids using BEM

A numerical method to study three-dimensional (3D) contact problems in solids with anisotropic elastic behavior is developed in this work. This formulation is based on the Boundary Element Method (BEM) for computing the elastic influence coefficients and on projection functions over the augmented Lagrangian for contact restrictions fulfillment. The constitutive equations of the potential contact zone are Signorini’s contact conditions and Coulomb’s law of friction. The formulation uses a recently introduced explicit approach for fundamental solutions evaluation, which are valid for general anisotropic behavior meanwhile mathematical degeneracies are allowed. The accuracy and robustness of the proposed method is illustrated by solving some examples previously presented in the literature. This approach is further applied to study the influence of solids anisotropy on the contact problem.     

Boundary element method for thermoelastic analysis of three-dimensional transversely isotropic solids

Thermal effects are well known to manifest themselves as additional volume integral terms in the direct formulation of the boundary integral equation (BIE) for linear elastic solids when using the boundary element method (BEM). This domain integral has been successfully transformed in an exact manner to surface ones only in isotropy and in 2D anisotropy, thereby restoring the BEM as a truly boundary solution technique. The difficulties with extending it to 3D general anisotropic solids lie in the mathematical complexity of the Green’s function and its derivatives for such materials. These quantities are required items in the BEM formulation. In this paper, the exact, analytical transformation of the volume integral associated with thermal effects to surface ones is achieved for a transversely isotropic material using a similar approach which the authors have previously employed for the same task in BEM for 2D general anisotropy. A numerical scheme, however, needs to be employed to evaluate some of the new terms introduced in the surface integrals that arise from this process here. The mathematical soundness of the formulation is demonstrated by a few examples; the numerical results obtained are checked by alternative means, including those obtained from the commercial FEM code, ANSYS.     

Decoupling coefficients of dilatational wave for Biot’s dynamic equation and its Green’s functions in frequency domain

Green’s functions for Biot’s dynamic equation in the frequency domain can be a highly useful tool for the investigation of dynamic responses of a saturated porous medium. Its applications are found in soil dynamics, seismology, earthquake engineering, rock mechanics, geophysics, and acoustics. However, the mathematical work for deriving it can be daunting. Green’s functions have been presented utilizing an analogy between the dynamic thermoelasticity and the dynamic poroelasticity in the frequency domain using the u-p formulation. In this work, a special term “decoupling coefficient” for the decomposition of the fast and slow dilatational waves is proposed and expressed to present a new methodology for deriving the poroelastodynamic Green’s functions. The correctness of the solution is demonstrated by numerically comparing the current solution with Cheng’s previous solution. The separation of the two waves in the present methodology allows the more accurate evaluation of Green’s functions, particularly the solution of the slow dilatational wave. This can be advantageous for the numerical implementation of the boundary element method (BEM) and other applications.     

Extended meshfree analysis of transverse and inplane loading of a laminated anisotropic plate of general planform geometry

An extended meshfree method is presented for the analysis of a laminated anisotropic plate under elastostatic loading. The plate may be of any planform shape with its thickness profile composed of perfectly bonded uniform thickness layers of distinct anisotropic materials. Both transverse and inplane loads are considered using a first order shear deformation theory for flexural behavior and generalized plane stress for the membrane behavior. In this extended meshfree method, a rectangular domain is initially considered with the plate of arbitrary geometry inscribed within it. A particular solution in the form of an analytic generalized Navier solution (a compound double Fourier series) is used to capture the response due to the loading within the rectangular domain. Then, a homogeneous solution by meshfree analysis is added to treat the augmented boundary conditions on the actual contour of the plate. These augmented conditions are composed of the prescribed values and that of the particular solution evaluated around the plate’s contour.Concentrated transverse and inplane loads in the form of uniform loads over a very small patch are considered with this generalized Navier solution representation. When a meshfree portion is added to account for the boundary conditions, such solutions constitute the Green’s functions for the plate. The viability of these double Fourier series representations is shown by the convergence rates for the kinematic and force/moment fields. An additional example of a two layer ±30° angleply circular plate is given to illustrate the capability of this extended meshfree method.     

On Green’s correlation of Stokes’ equation

The derivation of Green’s correlation naturally arises when identifying a linear propagation medium with uncontrolled random sources or ambient noise. As expected, this involves convolution of the well known Green’s function with its time-reversed version. The purpose of this paper is to derive a general expression of Green’s correlation function of a linear visco-acoustic propagation medium, in which the pressure field satisfies Stokes’ equation. From the expression obtained for a visco-acoustic medium, the Ward identity that was recently obtained for unbounded media is extended to the case of bounded propagation media. This extension appears necessary as the unbounded model is not valid in many practical cases, as for acoustic rooms for example. It is illustrated with both simulations and real-world aerial acoustics experimental data recorded in a closed room and in the framework of passive identification. In these experiments, Green’s correlation is estimated by the classical coda-based approach, and the performances are studied in this new context.     

Structure and properties of the solution space of general anisotropic laminates

The algebraic structure of the solution space of all types of anisotropic laminates is determined. The full space is shown to be the direct sum of a number of orthogonal eigenspaces, one for each simple or multiple eigenvalue, whose dimension equals the multiplicity. There are eight different types of eigenvalues, which combine to yield eleven distinct types of laminates with peculiar representations of the general solution. All such representations are explicitly obtained, along with the pseudo-metrics based on the binary product of the eigenvectors. This leads to the projection operators in the solution space, spectral sums and intrinsic tensors analogous to the Stroh–Barnett–Lothe tensors in 2-D elasticity. The present theoretical results are obtained by adopting a mixed formulation involving the deflection function and Airy’s stress function, and by using new laminate elasticity matrices different from the conventional stiffness matrices A, B and D. The new formulation also discloses an isomorphism relating each anisotropic laminate to an image laminate, such that every equilibrium solution of the former directly yields an image solution of the latter by interchanging the kinematical and kinetic variables and the in-plane and out-of-plane variables. This implies, in particular, that the classical bending theory of homogeneous plates and symmetric laminates is not a distinct subject, despite its historical development and pedagogical recognition, but is mathematically identical to the plane stress problem of anisotropic elasticity.     

On linking n-dimensional anisotropic and isotropic Green’s functions for infinite space,half-space,bimaterial, and multilayer for conduction

We establish exact mathematical links between the n-dimensional anisotropic and isotropic Green’s functions for diffusion phenomena for an infinite space, a half-space, a bimaterial and a multilayered space. The purpose of this work is not to attempt to present a solution procedure, but to focus on the general conditions and situations in which the anisotropic physical problems can be directly linked with the Green’s functions of a similar configuration with isotropic constituents. We show that, for Green’s functions of an infinite and a half-space and for all two-dimensional configurations, the exact correspondences between the anisotropic and isotropic ones can always be established without any regard to the constituent conductivities or any other information. And thus knowing the isotropic Green’s functions will readily provide explicit expressions for anisotropic Green’s functions upon back transformation. For three- and higher-dimensional bimaterials and layered spaces, the correspondence can also be found but the constituent conductivities need to satisfy further algebraic constraints. When these constraints are fully satisfied, then the anisotropic Green’s functions can also be obtained from those of the isotropic ones, or at least in principle.     

Three-dimensional Green’s functions for infinite anisotropic piezoelectric media

A direct, effective and concise method is adopted in this paper to find out the Green’s functions for infinite anisotropic piezoelectric media. The partial differential equations satisfied by the Green’s functions turn into a set of inhomogeneous algebraic equations after by using Fourier transform. Then inverse transform the solutions of the algebraic equations, the Green’s functions can be expressed by contour integral. Finally, the explicit expression can be obtained for the Green’s functions by using residual theory. The method demonstrated in this paper is easier to follow by people without knowledge of Radon transform, which has been used to obtain the Green’s functions by others.     

Thermo-elastic crack branching in general anisotropic media

Thermal fields may exist in addition to mechanical loading, for example, due to short term exposure to fire. In this paper, the branching of cracks in the presence of combined thermal and mechanical loads is investigated for general anisotropic media by employing the theory of Stroh’s dislocation formalism, extended to thermo-elasticity in matrix notation. A general solution to the thermo-elastic crack problem for an anisotropic material under arbitrary loading is obtained in a compact form. Green’s functions are also presented for a thermal dislocation (heat vortex) and a conventional dislocation (or, referred as mechanical dislocation), which are formulated considering the cuts located at an arbitrary angle with respect to the x₁ axis of the coordinate system (x₁, x₂, x₃). Using the derived compact expressions, the interaction between the crack and the dislocation is studied and a closed form solution for this interaction is obtained. The branching portion of the thermo-elastic crack is modelled as a continuous distribution of dislocations. This problem is then converted into a set of singular integral equations. Numerical results are presented to illustrate the possible effects of thermal loading on the propagation of the branched crack.