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.     

Green’s functions of two-dimensional piezoelectric quasicrystal in half-space and bimaterials

In this paper, we obtain Green’s functions of two-dimensional (2D) piezoelectric quasicrystal (PQC) in half-space and bimaterials. Based on the elastic theory of QCs, the Stroh formalism is used to derive the general solutions of displacements and stresses. Then, we obtain the analytical solutions of half-space and bimaterial Green’s functions. Besides, the interfacial Green’s function for bimaterials is also obtained in the analytical form. Before numerical studies, a comparative study is carried out to validate the present solutions. Typical numerical examples are performed to investigate the effects of multi-physics loadings such as the line force, the line dislocation, the line charge, and the phason line force. As a result, the coupling effect among the phonon field, the phason field, and the electric field is prominent, and the butterfly-shaped contours are characteristic in 2D PQCs. In addition, the changes of material parameters cause variations in physical quantities to a certain degree.     

Green’s functions for anti-plane deformations of a circular arc-crack at the interface of piezoelectric materials

Summary The anti-plane deformation problem of an interfacial debounding crack between a circular piezoelectric inclusion and a piezoelectric matrix is investigated by means of the complex variables method. For a line load applied within the matrix or inside the inclusion, Greens functions are presented for the complex potentials, intensity factors and electric fields on the crack faces, respectively, in closed and explicit form. The solutions are valid for both permeable and impermeable crack models. It is shown that, in the general case of permeable cracks, the electric field singularity is always proportional to the stress singularity.The first author (C.F.Gao) would like to express his gratitude for the support of the Alexander von Humboldt Foundation (Germany).     

Green’s functions of a surface-stiffened transversely isotropic half-space

Green’s functions of a transversely isotropic half-space overlaid by a thin coating layer are analytically obtained. The surface coating is modeled by a Kirchhoff thin plate perfectly bonded to the half-space. With the aid of superposition technique and employing appropriate displacement potential functions, the Green’s functions are expressed in two parts; (i) a closed-form part corresponding to the transversely isotropic half-space with surface kinematic constraints, and (ii) a numerically evaluated part reflecting the interaction between the half-space and the plate in the form of semi-infinite integrals. Some limiting cases of the problem such as surface-stiffened isotropic half-space, Boussinesq and Cerruti loadings, and extremely flexible and rigid plates are also studied. For the classical Cerruti problem in transversely isotropic materials, the effects of incompressibility are highlighted. Numerical results are provided to show the effects of material anisotropy, relative stiffness factor, and load buried depth. The obtained Green’s functions play a key role in treating further mixed-boundary-value problems in surface stiffened transversely isotropic half-spaces.     

Green’s functions of an exponentially graded transversely isotropic half-space

By virtue of a complete set of displacement potential functions and Hankel transform, the analytical expressions of Green’s function of an exponentially graded elastic transversely isotropic half-space is presented. The given solution is analytically in exact agreement with the existing solution for a homogeneous transversely isotropic half-space. Employing a robust asymptotic decomposition technique, the Green’s function is decomposed to the closed-form Green’s function corresponding to the homogeneous transversely isotropic half-space and grading term with strong decaying integrands. This representation is very useful for numerical methods which are based on boundary-integral formulations such as boundary-element method since the numerically evaluated part is not responsible for the singularity. The high accuracy of the proposed numerical scheme is confirmed by some numerical examples.     

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.     

Static equilibrium of a thermoelastic half-plane: Green’s functions and solutions in integrals

This article presents in a closed form new influence functions of a unit point heat source on the displacements for three boundary value problems of thermoelasticity for a half-plane. We also obtain the corresponding new integral formulas of Green’s and Poisson’s types that directly determine the thermoelastic displacements and stresses in the form of integrals of the products of specified internal heat sources or prescribed boundary temperature and constructed already thermoelastic influence functions (kernels). All these results are presented in terms of elementary functions in the form of three theorems. Based on these theorems and on derived early by author the general Green-type integral formula, we obtain in elementary functions new solutions to two particular boundary value problems of thermoelasticity for half-plane. The graphical presentation of the temperature and thermal stresses of one concrete boundary value problems of thermoelasticity for half-plane also is included. The proposed method of constructing thermoelastic Green’s functions and integral formulas is applicable not only for a half-plane, but also for many other two- and three-dimensional canonical domains of different orthogonal coordinate systems.     

Green’s functions for a loaded rolling tyre

A new formulation to determine the unit impulse response (Green’s) functions of a loaded rotating tyre in the vehicle-fixed (Eulerian) reference frame for tyre/road noise predictions is presented. The proposed formulation makes use of the set of eigenfrequencies and eigenmodes for the statically loaded tyre obtained from a finite element (FE) model of the tyre. A closed-form expression for the Green’s functions of a rotating tyre in the Eulerian reference system as a function of the eigenfrequencies and eigenmodes of the statically loaded tyre is found. Non-linear effects during loading are accounted for in the FE model, while the frequency shift due to the rotational velocity is included in the calculation of the Green’s functions. In the literature on tyre/road noise these functions are generally used to determine the tyre response during tyre/road contact calculations. The presented formulation opens the possibility to solve the contact problem directly in the Eulerian reference frame and to include local tyre softening due to non-linear effects while keeping the computational advantage of describing the tyre dynamics as a set of impulse response functions. The advantage of obtaining the Green’s functions in the Eulerian reference system is that only the Green’s functions corresponding to the potential contact zone need to be determined, which significantly reduces the computational cost of solving the tyre/road contact and since the mesh is fixed in space, a finer mesh can be used for the potential contact zone, improving the accuracy of the contact force calculations. Although these effects might be less pronounced if a more accurate tyre model is used, it is found that using the Green’s functions of the loaded tyre in a contact force calculation leads to smaller forces than in the unloaded case, lower frequencies are present in the response and they decrease faster as the rotational velocity increases.     

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.     

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.     

Three-dimensional Green’s functions for transversely isotropic thermoelastic bimaterials

Green’s functions for transversely isotropic thermoelastic biomaterials are established in the paper. We first express the compact general solutions of transversely isotropic thermoelastic material in terms of harmonic functions and introduce six new harmonic functions. The three-dimensional Green’s function having a concentrated heat source in steady state is completely solved using these new harmonic functions. The analytical results show some new phenomena of temperature and stress distributions at the interface. The temperature contours are normal to the interface for the isotropic material but not for the orthotropic one. The normal stress contours are parallel to the interface at the boundary in the isotropic region only and shear failure is most likely at the heat source due to the highly degenerated direction of shear stress contours.     

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.     

Dynamic Green’s functions of a two-phase saturated medium subjected to a concentrated force

The Green’s functions of a two-phase saturated medium subjected to a concentrated force are known to play an important role in seismology, earthquake engineering, soil dynamics, geophysics, and dynamic foundation theory. This paper presents a physical method for obtaining the dynamic Green’s functions of a two-phase saturated medium for materials considered to be isotropic and for low frequencies. First, the pore-fluid pressure in a two-phase saturated medium is divided into two parts: flow pressure and deformation pressure. Next, based on the compatibility condition of Biot’s equation and the property of the δ-function, the problem of coupled_fast and slow dilational waves is solved using the decomposition condition of the potential dilation field. The Green’s function for a concentrated force is then obtained by solving Biot’s complex modular equations, and their physical characteristics are discussed. The behavior of Green’s functions for the solid and fluid phases of a δ-impulsive force is investigated, from which the Green’s functions for a unit Heaviside force are also obtained by time integration. Finally, the present Green’s functions for a unit Heaviside force are compared with those obtained by a purely mathematical method; the two differ in form, but the numerical results are identical. The physical meaning of the expressions of Green’s functions obtained in this paper is evident. Therefore, the results may benefit future research on the dynamic responses of a two-phase saturated medium.     

Three-dimensional dynamic Green’s functions for a multilayered transversely isotropic half-space

With the aid of a method of displacement potentials, an efficient and accurate analytical derivation of the three-dimensional dynamic Green’s functions for a transversely isotropic multilayered half-space is presented. Constituted by proper algebraic factorizations, a set of generalized transmission–reflection matrices and internal source fields that are free of any numerically unstable exponential terms are proposed for effective computations of the potential solution. Three-dimensional point-load Green’s functions for stresses and displacements are given, for the first time, in the complex-plane line-integral representations. The present formulations and solutions are analytically in exact agreement with the existing solutions given by Pak and Guzina (2002) for the isotropic case. For the numerical computation of the integrals, a robust and effective methodology which gives the necessary account of the presence of singularities including branch points and poles on the path of integration is laid out. A comparison with the existing numerical solutions for multilayered isotropic half-space is made to confirm the accuracy of the numerical solutions.     

Green’s functions for plane problem under various boundary conditions and applications

Green’s functions of a point dislocation as well as a concentrated force for the plane problem of an infinite plane containing an arbitrarily shaped hole under stress, displacement, and mixed boundary conditions are stated. The Green’s functions are obtained in closed forms by using the complex stress function method along with the rational mapping function technique, which makes it possible to deal with relatively arbitrary configurations. The stress functions for these problems consist of two parts: a principal part containing singular and multi-valued terms, and a complementary part containing only holomorphic terms. These Green’s functions can be derived without carrying out any integration. The applications of the Green’s functions are demonstrated in studying the interaction of debonding and cracking from an inclusion with a line crack in an infinite plane subjected to remote uniform tension. The Green’s functions should have many other potential applications such as in boundary element method analysis. The boundary integral equations can be simplified by using the Green’s functions as the kernels.     

Static and Dynamic Green’s Functions in Peridynamics

We derive the static and dynamic Green’s functions for one-, two- and three-dimensional infinite domains within the formalism of peridynamics, making use of Fourier transforms and Laplace transforms. Noting that the one-dimensional and three-dimensional cases have been previously studied by other researchers, in this paper, we develop a method to obtain convergent solutions from the divergent integrals, so that the Green’s functions can be uniformly expressed as conventional solutions plus Dirac functions, and convergent nonlocal integrals. Thus, the Green’s functions for the two-dimensional domain are newly obtained, and those for the one and three dimensions are expressed in forms different from the previous expressions in the literature. We also prove that the peridynamic Green’s functions always degenerate into the corresponding classical counterparts of linear elasticity as the nonlocal length tends to zero. The static solutions for a single point load and the dynamic solutions for a time-dependent point load are analyzed. It is analytically shown that for static loading, the nonlocal effect is limited to the neighborhood of the loading point, and the displacement field far away from the loading point approaches the classical solution. For dynamic loading, due to peridynamic nonlinear dispersion relations, the propagation of waves given by the peridynamic solutions is dispersive. The Green’s functions may be used to solve other more complicated problems, and applied to systems that have long-range interactions between material points.     

Image singularities of Green’s functions for isotropic elastic bimaterials subjected to concentrated forces and dislocations

Green’s functions for isotropic materials in the two-dimensional problem for elastic bimaterials with perfectly bonded interface are reexamined in the present study. Although the Green’s function for an isotropic elastic bimaterial subjected to a line force or a line dislocation has been discussed by many authors, the physical meaning and the structure of the solution are not clear. In this investigation, the Green’s function for an elastic bimaterial is shown to consist of eight Green’s functions for a homogeneous infinite plane. One of the novel features is that Green’s functions for bimaterials can be expressed directly by knowing Green’s functions for the infinite plane. If the applied load is located in material 1, the solution for the half-plane of material 1 is constructed with the help of five Green’s functions corresponding to the infinite plane. However, the solution for the half-plane of material 2 only consists of three Green’s functions for the infinite plane. One of the five Green’s functions of material 1 and all the three Green’s functions of material 2 have their singularities located in the half-plane where the load is applied, and the other four image singularities of material 1 are located outside the half-plane at the same distance from the interface as that of the applied load. The nature and magnitude of the image singularities for both materials are presented explicitly from the principle of superposition, and classified according to different loads. It is known that for the problem of anisotropic bimaterials subjected to concentrated forces and dislocations, the image singularities are simply concentrated forces and dislocations with the stress singularity of order O(1/r). However, higher orders (O(1/r²) and O(1/r³)) of stress singularities are found to exist in this study for isotropic bimaterials. The highest order of the stress singularity is O(1/r³) for the image singularities of material 1, and is O(1/r²) for material 2. Using the present solution, Green’s functions associated with the problems of elastic half-plane with free and rigidly fixed boundaries, for homogeneous isotropic elastic solid, are obtained as special cases.     

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.     

Local existence for Frémond’s model of damage in elastic materials

We consider a dissipative model recently proposed by M. Frémond to describe the evolution of damage in elastic materials. The corresponding PDEs system consists of an elliptic equation for the displacements with a degenerating elastic coefficient coupled with a variational dissipative inclusion governing the evolution of damage. We prove a local-in-time existence and uniqueness result for an associated initial and boundary value problem, namely considering the evolution in some subinterval where the damage is not complete. The existence result is obtained by a truncation technique combined with suitable a priori estimates. Finally, we give an analogous local-in-time existence and uniqueness result for the case in which we introduce viscosity into the relation for macroscopic displacements such that the macroscopic equilibrium equation is of parabolic type.Received: 31 July 2002, Accepted: 9 August 2003, Published online: 21 November 2003Correspondence to: E. Bonetti