Evaluation of the three-dimensional elastic Green's function in anisotropic cubic media

By using the Fourier transforms method, the three-dimensional Green's function solution for a unit force applied in an infinite cubic material is evaluated in this paper. Although the elastic behavior of a cubic material can be characterized by only three elastic constants, the explicit solutions of Green's function for a cubic material are not available in the literatures. The central problem for explicitly solving the elastic Green's function of anisotropic materials depends upon the roots of a sextic algebraic equation, which results from the inverse Fourier transforms and is composed of the material constants and position vector parameters. The close form expression of Green's function is presented here in terms of roots of the sextic equation. The sextic equation for an anisotropic cubic material is discussed thoroughly and specific results are given for possible explicit solutions.     

Efficient computation of the Green's function and its derivatives for three-dimensional anisotropic elasticity in BEM analysis

An alternative scheme to compute the Green's function and its derivatives for three dimensional generally anisotropic elastic solids is presented in this paper. These items are essential in the formulation of the boundary element method (BEM); their evaluation has remained a subject of interest because of the mathematical complexity. The Green's function considered here is the one introduced by Ting and Lee [Q. J. Mech. Appl. Math. 1997; 50: 407–26] which is of real-variable, explicit form expressed in terms of Stroh's eigenvalues. It has received attention in BEM only quite recently. By taking advantage of the periodic nature of the spherical angles when it is expressed in the spherical coordinate system, it is proposed that this Green's function be represented by a double Fourier series. The Fourier coefficients are determined numerically only once for a given anisotropic material; this is independent of the number of field points in the BEM analysis. Derivatives of the Green's function can be performed by direct spatial differentiation of the Fourier series. The resulting formulations are more concise and simpler than those derived analytically in closed form in previous studies. Numerical examples are presented to demonstrate the veracity and superior efficiency of the scheme, particularly when the number of field points is very large, as is typically the case when analyzing practical three dimensional engineering problems.     

Green's function for anisotropic dielectric halfspace

A spatial representation of the Green's function for anisotropic halfspace is considered. This representation can be used for analyzing an electric field connected with planar metal structures applied in surface-acoustic-wave devices     

A numerical Green's function for multiple cracks in anisotropic bodies

The numerical construction of a Green's function for multiple interacting planar cracks in an anisotropic elastic space is considered. The numerical Green's function can be used to obtain a special boundary-integral method for an important class of two-dimensional elastostatic problems involving planar cracks in an anisotropic body.     

Green's function for two-and-a-half dimensional elastodynamic problems in a half-space

 This paper presents an analytical solution, together with explicit expressions, for the steady state response of a homogeneous three-dimensional half-space subjected to a spatially sinusoidal, harmonic line load. These equations are of great importance in the formulation of three-dimensional elastodynamic problems in a half-space by means of integral transform methods and/or boundary elements. The final expressions are validated here by comparing the results with those obtained with the boundary element method (BEM) solution, for which the free surface of the ground is discretized with boundary elements.     

Efficient implementation of anisotropic three dimensional boundary-integral equation stress analysis

The boundary-integral equation medthod is particularly well suited for solution of stress concentration and elastic fracture mechanics problems. The method was not previously applicable to anisotropic three dimensional problems because no efficient technique existed for calculation of the required point load solution for an infinite body. A technique has been developed to evaluate numerically the anisotropic point load soultions, and used to generate data bases for various materials. An intrpolation technique is used to evaluate the point load solutions efficiently within a higher order boundary-intgral equation code. The effectiveness of the technique is verified by solution of problems involving both uniaxial stress states and stress concentrations.     

Decoupling the nonlocal elasticity equations for three dimensional vibration analysis of nano-plates

In this paper, a three dimensional vibration analysis of nano-plates is studied by decoupling the field equations of Eringen theory. Considering the small scale effect, the three dimensional equations of nonlocal elasticity are obtained. At first, three decoupled equations in terms of displacement components and three decoupled equations in terms of rotation components are obtained. In order to find the solution for a nano-plate based on the presented formulation, one of the three equations in terms of displacement components and corresponding rotation equation should be solved independently. Using some relations, the other two displacement components can be obtained in terms of the mentioned displacement and rotation component. A Navier-type method for finding the exact three dimensional solution of a nano-plate is presented using the Fourier series technique. Exact natural frequencies of nano-plates are presented and compared with the results of nonlocal first order and third order shear deformation theories.     

Real-space Green's function of an isolated point charge in an unbounded anisotropic dielectric

In modern electronic devices the finite size of the substrates requires three-dimensional analysis. To this end, the boundary element method (BEM) may be utilized. The BEM involves problem-specific Green's functions (GFs), which are to be constructed first. This work has derived a GF which is the electric potential response to an isolated point charge in an unbounded anisotropic dielectric.     

Green's function method in the radiative transfer problem. II. Spatially heterogeneous anisotropic surface

The most recent theoretical studies have shown that three-dimensional (3-D) radiation effects play an important role in the optical remote sensing of atmospheric aerosol and land surface reflectance. These effects may contribute notably to the error budget of retrievals in a broad range of sensor resolutions, introducing systematic biases in the land surface albedo data sets that emerge from the existing global observation systems. At the same time, 3-D effects are either inadequately addressed or completely ignored in data processing algorithms. Thus there is a need for further development of the radiative transfer theory that can rigorously treat both 3-D and surface anisotropy effects and yet be flexible enough to permit the development of fast forward and inversion algorithms. We describe a new theoretical solution to the 3-D radiative transfer problem with an arbitrary nonhomogeneous non-Lambertian surface. This solution is based on an exact semianalytical solution derived in operator form by the Green's function method. The numerical implementation is based on several parameterizations that accelerate the solution dramatically while keeping its accuracy within several percent under most general conditions.     

A new boundary element method formulation for three dimensional problems in linear elasticity

Summary A new boundary element method (BEM) formulation for planar problems of linear elasticity has been proposed recently [6]. This formulation uses a kernel which has a weaker singularity relative to the corresponding kernel in the standard formulation. The most important advantage of the new formulation, relative to the standard one, is that it delivers stresses accurately at internal points that are extremely close to the boundary of a body. A corresponding BEM formulation for three dimensional problems of linear elasticity is presented in this paper. This formulation is derived through the use of Stokes' theorem and has kernels which are only 1/r singular (wherer is the distance between a source and a field point) for the displacement equation. The standard BEM formulation for three-dimensional elasticity problems has a kernel which is 1/r ² singular.With 2 Figures     

Comments on real-space Green's function of an isolated point-chargein an unbounded anisotropic medium

The electric potential of an isolated point-charge in an unbounded anisotropic dielectric can be derived by an approach which merely involves a coordinate transformation and the knowledge of the solution to the corresponding isotropic problem     

Spectral decomposition of anisotropic elasticity

Summary A newly developed approach, based on the spectral decomposition principle, which is especially useful in crystallography, is applied in this paper. The compliance fourth-rank tensor of crystalline media belonging to the monoclinic system is spectrally decomposed, its eigenvalues are evaluated, together which its elementary idempotent tensors, which expand uniquely the fourth-rank tensor space into orthogonal subspaces. Next, the compliance tensor is spectrally analysed for anisotropic media of the orthorhombic, tetragonal, hexagonal and cubic crystal systems, by regarding these decompositions as particular cases of the spectral decomposition for monolinic media. Consequently, the characteristic values and the idempotent fourth-rank tensors are derived, as well as the stress and strain second-rank eigentensors for all the above mentioned symmetries.     

A numerical integration Green's function model for ultrasonic field profiles in mildly anisotropic media

A numerical integration technique utilizing a point source Green's function is introduced to analyze the wave behavior in transversely isotropic-type anisotropic media allowing us to make fast and accurate computations of the acoustic field. The centrifugally cast stainless steel (CCSS) used in nuclear power plants is chosen as a sample medium because of its columnar grain character leading to material anisotropy. A representative number of field profiles are computed and plotted to illustrate the quasi-longitudinal, quasi-transverse, and horizontally-polarized shear wave propagation in a transversely-isotropic medium. Phenomena such as beam skewing, beam splitting, beam focusing, unsymmetrical beams, and other anisotropic effects, some of which are already known from earlier experimental observations, emerge as a computational result of the introduced technique.     

Screw dislocation in nonlocal anisotropic elasticity

Based on Eringen’s model of nonlocal anisotropic elasticity, new solutions for the stress fields of screw dislocations in anisotropic materials are derived. In the theory of nonlocal anisotropic elasticity the anisotropy is twofold. The anisotropic material behavior is not only included in the anisotropy of the elastic stiffness properties, but also in the anisotropy of the nonlocality which is expressed by the anisotropy of the length scale parameters, which is incorporated in the anisotropy of the nonlocal kernel function. Particularly, a new two-dimensional anisotropic kernel which is the Green function of a linear differential operator with three length scale parameters is derived analytically. New solutions for the stresses of straight screw dislocations in anisotropic (monoclinic and hexagonal) materials are found. The stresses do not have singularities and possess interesting features of anisotropy, which are presented and discussed.     

Quasiclassical Green's function in the BCS pairing theory

We study a generalization of the quasiclassical Green's function which allows us to include the transfer of momentum to the particles. In this approach, we may handle Galilei transformations, rotations, and gauge transformations in a systematic way. As an example, we calculate the quasiparticle flow pattern which arises during the motion of the orbital vector in the ABM phase, and discuss the meaning of the intrinsic angular momentum of the Cooper pairs. Finally, we consider charged particles in a magnetic field, and derive a Boltzmann equation for a superconductor which applies to the Hall effect in the case of moving vortices.On leave from Institut für Théorie der Kondensierten Materie, Universität Karlsruhe.     

Efficient and accurate multi-layered elastostatic Green's functions via the bi-material Green's function

The three-dimensional elastostatic Green's functions for multi-layered half-spaces are expressed in terms of a cylindrical system of vector functions using the propagator matrix method. To obtain efficient and accurate results near the singularity at the interface between layers, the singularity extraction method is employed in the transform domain, based on the singular solutions for the bi-material full space expressed here in vector cylindrical function form. Numerical trials demonstrate that this approach improves on previous work, which relied on singularity extraction based on the well-known singular (Kelvin) solution for the homogenous full-space. A boundary element implementation illustrates the practical potential of this work to boundary value problems.     

An alternating iterative algorithm for the Cauchy problem in anisotropic elasticity

The alternating iterative algorithm proposed by Kozlov et al. [An iterative method for solving the Cauchy problem for elliptic equations. USSR Comput Math Math Phys 1991;31:45–52] for obtaining approximate solutions to the Cauchy problem in two-dimensional anisotropic elasticity is analysed and numerically implemented using the boundary element method (BEM). The ill-posedness of this inverse boundary value problem is overcome by employing an efficient regularising stopping. The numerical results confirm that the iterative BEM produces a convergent and stable numerical solution with respect to increasing the number of boundary elements and decreasing the amount of noise added into the input data.     

Fundamental solution of the incompatibility problem in three-dimensional infinite anisotropic elasticity theory

Summary The decomposition of the strain field into its deformation and its incompatible parts which is unique in an infinite elastic medium is used to derive a fourth-rank Green's function tensorP _abij (fundamental solution) of the incompatibility problem of linear anisotropic elasticity theory. This fundamental solution is homogeneous of degree — 1 and non-unique. For the calculation of internal stresses only one of the possible fundamental solutions is needed. In order to solve the integration problem, it is convenient to decompose the fundamental solution into two parts. The one part consists of a particular solution of the incompatibility condition (Eq. (3)) and is obtained in a straightforward manner. The other part results from a deformation part of the strain. The integration of this part is worked out by using the theory of residues. In the case of hexagonal crystal symmetry an explicit solution is presented.     

An anisotropic theory of elasticity for continuum damage mechanics

This paper presents the development of an anisotropic elastic damage theory. This is achieved by deriving a modified damage effect tensor M(D) for the effective stress equations capable of including the effect of anisotropic material damage. The modified tensor removes the restriction of a priori knowledge of the directions of principal stresses imposed by a damage effect tensor developed earlier and can now be made for general practical engineering applications of failure analysis. Reduction of the proposed tensor to a scalar for isotropic damage is shown to be possible when it is expressed not only in the principal directions but also in any arbitrary coordinate system, a necessary condition to verify the validity of the proposed tensor. Uniaxial tension and pure torsion are chosen to illustrate the application of the theory as well as associated damage variables that may be experimentally determined using laboratory size specimens. The measured damage variables confirm the presence of anisotropic damage from an initially isotropic material specimen and the magnitude is more pronounced at higher stresses and strains.

---

Résumé On présente un développement d'une théorie sur l'endommagement élastique anisotrope en déduisant un tenseur modifié décrivant l'effet de l'endommagement pour un système d'équations de contraintes effectives susceptible d'inclure l'effet d'un endommagement dans un matériau anisotrope. Le tenseur modifié supprime la restriction de la connaissance a priori des directions des contraintes principales imposées par un tenseur d'effet d'endommagement développé précédemment; il peut à présent entrer dans les applications pratiques en construction de l'analyse des ruptures.On montre qu'il est possible de réduire le tenseur proposé à une valeur scalaire dans le cas d'un dommage isotrope, dès lors qu'il est exprimé non seulement suivant les directions principales, mais dans un système de coordonnées arbitraires, ce qui est une condition nécessaire pour en vérifier la validité.On choisit une traction multiaxiale et une torsion pure pour illustrer l'application de la théorie ainsi que des variables d'endommagement associées, susceptibles d'être déterminées expérimentalement à l'aide d'éprouvettes de laboratoire.Les variables d'endommagement mesurées confirment la présence d'un dommage anisotrope dans le cas d'une éprouvette d'un matériau initialement isotrope; son amplitude est plus prononcée à des contraintes ou des déformations plus importantes.

     

Surface Green's function of a piezoelectric half-space

The computation of the two-dimensional harmonic spatial-domain Green's function at the surface of a piezoelectric half-space is discussed. Starting from the known form of the Green's function expressed in the spectral domain, the singular contributions are isolated and treated separately. It is found that the surface acoustic wave contributions (i.e., poles in the spectral Green's function) give rise to an anisotropic generalization of the Hankel function H0(2), the spatial Green's function for the scalar two-dimensional wave equation. The asymptotic behavior at infinity and at the origin (for the electrostatic contribution) also are explicitly treated. The remaining nonsingular part of the spectral Green's function is obtained numerically by a combination of fast Fourier transform and quadrature. Illustrations are given in the case of a substrate of Y-cut lithium niobate.