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.     

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.     

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.     

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 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.     

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.     

Two-dimensional Eshelby’s problem for piezoelectric materials with a parabolic boundary

We use complex variable techniques to obtain analytic solutions of Eshelby’s problem consisting of an inclusion of arbitrary shape in an anisotropic piezoelectric plane with a parabolic boundary. The region of the physical plane below the parabola is mapped onto the lower half of the image plane. The problem is then more conveniently studied in the image plane rather than in the physical plane. The critical step in our approach lies in the construction of certain auxiliary functions in the image plane which allow for the technique of analytic continuation to be applied to an inclusion of arbitrary shape.     

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.     

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.     

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.     

Prandtl’s formulation for the Saint–Venant’s torsion of homogeneous piezoelectric beams

The Saint–Venant torsional problem for homogeneous, monoclinic piezoelectric beams is formulated in terms of Prandtl’s stress function and electric displacement potential function. The analytical approach presented in this paper generalizes the known formulation of Prandtl’s solution which refers to homogeneous elastic beams. The Prandtl’s stress function and electric displacement potential function satisfy the so called coupled Dirichlet problem (CDP) in the cross-sectional domain. A direct and a variational formulation are developed. Exact analytical solutions for solid elliptical cross-section and hollow circular cross-section and an approximate solution based on a variational formulation for thin-walled closed cross-section are presented.     

Anti-plane shear Green’s functions for an isotropic elastic half-space with a material surface

This paper presents analytical Green’s function solutions for an isotropic elastic half-space subject to anti-plane shear deformation. The boundary of the half-space is modeled as a material surface, for which the Gurtin–Murdoch theory for surface elasticity is employed. By using Fourier cosine transform, analytical solutions for a point force applied both in the interior or on the boundary of the half-space are derived in terms of two particular integrals. Through simple numerical examples, it is shown that the surface elasticity has an important influence on the elastic field in the half-space. The present Green’s functions can be used in boundary element method analysis of more complicated problems.     

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.     

Minimization of semicoercive functions: a generalization of Fichera’s existence theorem for the Signorini problem

The existence theorem of Fichera for the minimum problem of semicoercive quadratic functions in a Hilbert space is extended to a more general class of convex and lower semicontinuous functions. For unbounded domains, the behavior at infinity is controlled by a lemma which states that every unbounded sequence with bounded energy has a subsequence whose directions converge to a direction of recession of the function. Thanks to this result, semicoerciveness plus the assumption that the effective domain is boundedly generated, that is, admits a Motzkin decomposition, become sufficient conditions for existence. In particular, for functions with a smooth quadratic part, a generalization of the existence condition given by Fichera’s theorem is proved.     

Analytical solution for squeeze film damping of MEMS perforated circular plates using Green’s function

Squeezed air film between two closely spaced vibrating microstructures is the important source of energy dissipation and has profound effects on the dynamics of microelectromechanical systems (MEMS). Perforations in the design are one of the methods to model these damping effects. The literature reveals that the analytical modeling of squeeze film damping of perforated circular microplates is less explored; however, these microplates are also an imperative part of the numerous MEMS devices. Here, we derive an analytical model of transverse and rocking motions of a perforated circular microplate. A modified Reynolds equation that incorporates compressibility and rarefaction effects is utilized in the analysis. Pressure distribution under the vibrating microplate is derived by using Green’s function and also derived by finite element method (FEM) to visualize the pressure distribution under perforated and non-perforated areas of the microplate. The analytical damping results are validated with previous renowned analytical models and also with the FEM results. The outcomes confirm the potential of the present analytical model to accurately predict the squeeze film damping parameters.     

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.     

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.     

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.     

On Green’s function for a bimaterial elastic half-plane

The problem of a point force acting in a composite, two-dimensional, isotropic elastic half-plane is considered. An exact solution is obtained, using Mellin transforms and the Melan solution for a point force in a homogeneous half-plane.