Next Article in Journal
Microstructure Evolution and Mechanical Properties of TiB2/Al–Cu–Mn–Cd Composite with the Adoption of Two-Stage Solution and Aging Treatment
Next Article in Special Issue
Three-Dimensional Axisymmetric Analysis of Annular One-Dimensional Hexagonal Piezoelectric Quasicrystal Actuator/Sensor with Different Configurations
Previous Article in Journal
Friedel Oscillations and He-He Interactions in Mo
Previous Article in Special Issue
A Closed-Form Solution to the Mechanism of Interface Crack Formation with One Contact Area in Decagonal Quasicrystal Bi-Materials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three-Dimensional and Two-Dimensional Green Tensors of Piezoelectric Quasicrystals

Institute of Mechanics, Darmstadt University of Technology, D-64287 Darmstadt, Germany
*
Author to whom correspondence should be addressed.
Crystals 2024, 14(10), 835; https://doi.org/10.3390/cryst14100835
Submission received: 4 September 2024 / Revised: 22 September 2024 / Accepted: 23 September 2024 / Published: 26 September 2024
(This article belongs to the Special Issue Structures, Properties and Applications of Quasicrystals)

Abstract

:
In this work, within the framework of the linear piezoelectricity theory of quasicrystals, the three-dimensional and two-dimensional Green tensors for arbitrary piezoelectric quasicrystals are derived. In the piezoelectricity of quasicrystals, where phonon, phason and electric fields exist, we introduce the corresponding multifields by developing a hyperspace notation for piezoelectric quasicrystals. Using Fourier transform and the multifield formalism, the three-dimensional Green tensor for piezoelectric quasicrystals as well as its spatial gradient necessary for applications, are derived. The solutions for the “displacement”, “distortion” and “stress” multifields in the presence of a “force” multifield in a piezoelectric quasicrystal as well as the solution of the generalised Kelvin problem, are given. In addition, the two-dimensional Green tensors of piezoelectric quasicrystals as well as of quasicrystals, are determined.

1. Introduction

Nowadays, piezoelectric materials are an intriguing class of materials with broad use in technology, including energy conversion and sensing as well as in energy harvesting and transducer devices (see, e.g., [1]). For a piezoelectric material, the elastic and electric response is coupled and anisotropic (see, e.g., [2,3,4]). The Green functions or fundamental solutions of partial differential equations play an essential role in the solution of many problems in various scientific areas, like mathematics, physics, engineering science and material science. In general, the importance of Green functions is at least twofold. A Green function is the fundamental solution of a partial differential equation and can be used to clarify various approximative methods, such as finite element and boundary element methods for the study of defects like cracks, dislocations, disclinations, point defects and inclusions. The three-dimensional elastic Green tensor was derived by Lord Kelvin [5] for isotropic materials and by Lifshitz and Rosenzweig [6], Kröner [7] and Willis [8] for hexagonal materials. Lifshitz and Rosenzweig [6] and Synge [9] (see also [10,11]) derived the three-dimensional elastic Green tensor for arbitrary anisotropic materials. The two-dimensional Green tensor of anisotropic elasticity was derived by Wang and Achenbach [12], Wang [13] and Lazar [14] using Radon transform and Fourier transform, respectively. In analogy to the three-dimensional Green tensor of anisotropic elasticity given by Synge [9], the three-dimensional Green tensor for arbitrary piezoelectric materials was derived by Deeg [15], Chen [16], Dunn [17], Akamatsu and Tanuma [18], Pan and Tonon [19], Li and Wang [20] and Nowacki [21]. The three-dimensional Green tensor for hexagonal piezoelectric crystals with point group 6   m m has been calculated by Michelitsch [22]. The two-dimensional Green tensor of piezoelectricity has been given by Wang [23].
To date, quasicrystals represent an interesting class of novel materials due to their particular properties, such as physical, structural, chemical, electronic and magnetic properties among others. Quasicrystals were discovered by Shechtman in 1982 and reported by Shechtman et al. [24] in 1984. In 2011, Shechtman received the Nobel Prize in Chemistry for the discovery of quasicrystals. Quasicrystals belong to aperiodic crystals and possess long-range orientational order but no translational symmetry. The basis of the continuum theory of solid quasicrystals is set up by two elementary excitations: the phonons and the phasons [25,26]. Quasicrystals possess piezoelectric effects, attracting interest for their applications in technology, for example, in electronic systems as sensor and actuator devices by which an electric voltage can induce an elastic deformation and vice versa. There is considerable interest in the subject of piezoelectric quasicrystals, which is an ongoing important and active research field (see, e.g., [27,28,29]).
Concerning Green functions in quasicrystals, De and Pelcovits [30] found the two-dimensional Green functions for pentagonal quasicrystals and Ding et al. [31] found the explicit expressions of two-dimensional Green tensors for various forms of planar quasicrystals. Bachteler and Trebin [32] gave an approximative solution for the elastic Green tensor of three-dimensional icosahedral quasicrystals, assuming that the coupling between phonons and phasons is small. Lazar and Agiasofitou [33] derived the three-dimensional elastic Green tensor for arbitrary quasicrystals.
In the theory of piezoelectric quasicrystals, the phonon fields as well as the phason fields are coupled with the electric field (see, e.g., [34,35]). It is clear that, due to this fact, the piezoelectric behaviour of quasicrystals is much more complicated than that of conventional piezoelectric crystals. An improvement to the existing constitutive modelling of piezoelectric quasicrystals has been recently given by Agiasofitou and Lazar [36] with a focus on the tensor of phason piezoelectric moduli, which is proved to be fully asymmetric without any major or minor symmetry. Moreover, it has been shown that piezoelectric effects in quasicrystals can exist only due to phonons or only due to phasons or due to both phonon and phason fields, depending, in general, on the Laue class as well as on the point group (see [35,36]). Using group representation theory, Hu et al. [35] studied piezoelectric effects in all two-dimensional quasicrystals with crystallographic and non-crystallographic symmetries, and in three-dimensional icosahedral and cubic quasicrystals. Agiasofitou and Lazar [36] have investigated piezoelectric effects in all one-dimensional quasicrystals with crystallographic and non-crystallographic symmetries. The interested reader can find concrete results concerning which crystallographic point groups possess piezoelectric effects and which do not, and if these are induced only by phonons or only by phasons or due to both phonon and phason fields in [36] and the references therein.
Regarding Green functions in piezoelectric quasicrystals, Green functions have been given only for particular cases of one-dimensional piezoelectric quasicrystals in the literature (see, e.g., [37,38,39]). In this work, we derive exact analytical expressions for the three-dimensional and two-dimensional Green tensor functions for arbitrary piezoelectric quasicrystals, that means for one-dimensional, two-dimensional and three-dimensional piezoelectric quasicrystals.
The paper is organized as follows. In Section 2, the basic framework of the generalised linear piezoelectricity theory of quasicrystals is presented. In Section 3, the hyperspace notation of piezoelectric quasicrystals is introduced, being an essential tool for the derivation of the three-dimensional Green tensor for arbitrary piezoelectric quasicrystals, which is done in Section 4 using the Fourier transform and the multifield formalism. Additionally, Section 4 includes other fundamental properties, such as the solution of the generalised Kelvin problem for piezoelectric quasicrystals and the derivation of the two-dimensional Green tensors for arbitrary piezoelectric quasicrystals as well as for arbitrary quasicrystals. Conclusions are given in Section 5. In Appendix A, the major symmetry of the multifield tensor of the material moduli in the hyperspace notation is proved. In Appendix B, the multifield tensor of the material moduli in the hyperspace notation is explicitly given for one-dimensional, two-dimensional and three-dimensional piezoelectric quasicrystals.

2. Basic Framework of Piezoelectric Quasicrystals

In this section, the basic framework of the generalised linear piezoelectricity theory of quasicrystals [36], which is the generalisation of the linear (compatible) elasticity of quasicrystals [40] towards the piezoelectricity of quasicrystals, is given.
An ( n 3 ) -dimensional quasicrystal can be generated by the projection of an n-dimensional periodic structure to the three-dimensional physical space ( n = 4 , 5 , 6 ). The n-dimensional hyperspace E n can be decomposed into the direct sum of two orthogonal subspaces,
E n = E 3 E ( n 3 ) ,
where E 3 is the three-dimensional physical or parallel space of the phonon fields and E ( n 3 ) is the ( n 3 ) -dimensional perpendicular space of the phason fields. Here, ⊕ denotes the direct sum. For n = 4 , 5 , 6 , we speak of one-dimensional, two-dimensional and three-dimensional quasicrystals with the dimension of the corresponding hyperspace being 4D, 5D, and 6D, respectively. Here, indices in the parallel space are denoted by small letters i , j , k , l , with i , j , k , l = 1 , 2 , 3 , and indices in the perpendicular space are denoted by Greek letters α , γ , with α , γ = 3 for one-dimensional quasicrystals (with quasiperiodicity in the z-direction), α , γ = 1 , 2 for two-dimensional quasicrystals, and α , γ = 1 , 2 , 3 for three-dimensional quasicrystals. Throughout the text, phonon fields will be denoted by ( · ) and phason fields by ( · ) . All field quantities (phonon and phason fields) depend on the so-called material space coordinates (or spatial coordinates) x R 3 . Note that in the linear theory of quasicrystals, the material space coincides with the parallel space.
In the theory of the generalised compatible elasticity of quasicrystals, the (elastic) phonon and phason distortion tensors, β k l and β γ l , are defined as the gradient of the phonon displacement vector  u k and of the phason displacement vector  u γ , respectively,
β k l = u k , l , β γ l = u γ , l
and they fulfill the compatibility conditions
ϵ i j l β k l , j = 0 , ϵ i j l β γ l , j = 0 ,
where ϵ i j l is the three-dimensional Levi–Civita tensor. The subscript comma denotes the partial differentiation j = / x j with respect to the spatial coordinates x j . The phonon strain tensor is given by
e k l = 1 2 ( u k , l + u l , k ) .
Using the two different coordinate systems and indices, it can be seen that a phason strain tensor cannot be defined (see the relative discussion in [36]). On the other hand, the electric field strength vector  E l is defined as the negative gradient of the electrostatic potential  φ
E l = φ , l
and it satisfies the Bianchi identity
ϵ i j l E l , j = 0 ,
which is the compatibility condition for the electric field. Equation (6) is nothing but the electrostatic version of the Faraday law [41].
For piezoelectric quasicrystals, the (force) equilibrium conditions read [36,42]
σ i j , j + f i = 0 ,
σ α j , j + f α = 0 ,
and the Gauss law of electrostatics (see, e.g., [41]) reads
D j , j = q .
Here, σ i j and σ α j are the phonon and phason stress tensors, respectively, D j is the electric displacement vector or electric excitation, f i is the conventional phonon body force density, f α is the phason body force density (see, e.g., [33]) and q denotes the body charge density. Note that the phonon stress tensor is symmetric, σ i j = σ j i , whereas the phason stress tensor σ α j is an asymmetric two-point tensor, σ α j σ j α (see also [33,42]).
Following the improved constitutive modelling for linear piezoelectric quasicrystals given recently by Agiasofitou and Lazar [36], the constitutive relations are
σ i j = C i j k l e k l + D i j γ l β γ l e l i j E l ,
σ α j = D k l α j e k l + E α j γ l β γ l f l α j E l ,
D j = e j k l e k l + f j γ l β γ l + ε j l E l .
In Equations (10)–(12), it can be seen that the electric field strength vector contributes to the phonon and phason stress tensors and that the elastic phonon strain tensor and phason distortion tensor contribute to the electric displacement vector for a linear piezoelectric quasicrystal. The constitutive tensors, which are the physical property tensors, possess the following symmetries [36]. The tensor of the elastic moduli of phonons  C i j k l possesses the major symmetry and the minor symmetries
C i j k l = C k l i j , C i j k l = C i j l k = C j i k l .
The tensor of the elastic moduli of phasons  E α j γ l possesses only a major symmetry
E α j γ l = E γ l α j ,
whereas the tensor of the elastic moduli of the phonon–phason coupling  D i j γ l possesses only a minor symmetry with respect to the first two indices
D i j γ l = D j i γ l .
In addition, the tensor of the phonon piezoelectric moduli  e k i j possesses a minor symmetry with respect to the last two indices
e k i j = e k j i ,
the tensor of the phason piezoelectric moduli  f j γ l is a fully asymmetric tensor possessing neither a major symmetry nor a minor symmetry with respect to the last two indices
f j γ l f j l γ
and the tensor of the dielectric moduli  ε i j exhibits the major symmetry
ε i j = ε j i .
Substituting Equations (10)–(12) into Equations (7)–(9), and using Equations (2), (4) and (5), the coupled equations of equilibrium for the phonon and phason displacement fields and the electrostatic potential are given by
C i j k l u k , l j + D i j γ l u γ , l j + e l i j φ , l j = f i ,
D k l α j u k , l j + E α j γ l u γ , l j + f l α j φ , l j = f α ,
e j k l u k , l j + f j γ l u γ , l j ε j l φ , l j = q ,
which are inhomogeneous partial differential equations with the body forces and body charges as source terms.
The electric enthalpy density  Ψ for piezoelectric quasicrystals is given by [36]
Ψ ( e i j , β α j , E j ) = 1 2 σ i j e i j + 1 2 σ α j β α j 1 2 D j E j ,
which, using Equations (10)–(12), takes the form
Ψ = 1 2 C i j k l e i j e k l + 1 2 D i j γ l e i j β γ l 1 2 e l i j E l e i j + 1 2 D k l α j e k l β α j + 1 2 E α j γ l β α j β γ l 1 2 f l α j E l β α j 1 2 e j k l E j e k l 1 2 f j γ l E j β γ l 1 2 ε j l E j E l .
The form of the electric enthalpy density (23) is useful in order to write the constitutive tensors of the material moduli as a multifield tensor in the hyperspace notation, which follows in the next section.

3. Hyperspace Notation and Multifields for Piezoelectric Quasicrystals

Throughout this paper, we introduce the hyperspace notation appropriate for the description of piezoelectric quasicrystals. The hyperspace notation for piezoelectric quasicrystals is the unification of the notation introduced by Barnett and Lothe [43] for piezoelectric materials and of the hyperspace notation introduced by Lazar and Agiasofitou [33] for quasicrystals towards piezoelectric quasicrystals.
For a piezoelectric quasicrystal, the corresponding hyperspace is an ( n + 1 ) -dimensional hyperspace E ( n + 1 ) which is decomposed into three orthogonal subspaces, E 3 , E ( n 3 ) , and E φ 1 ,
E ( n + 1 ) = E 3 E ( n 3 ) E φ 1 ,
where the additional space E φ 1 is the one-dimensional space of the electrostatic potential. For n = 4 , 5 , 6 , we speak of one-dimensional, two-dimensional and three-dimensional piezoelectric quasicrystals, and the dimension of the corresponding hyperspace being 5D, 6D, and 7D, respectively. The obtained fields which “live” in the hyperspace E ( n + 1 ) are called multifields. An index associated with a multifield is denoted by a capital letter, e.g., I , K = 1 , , n + 1 . In the hyperspace notation, we define the following multifields:
  • The “displacement” multifield vector  U K = ( u k , u γ , φ ) E 3 E ( n 3 ) E φ 1 or U K E ( n + 1 )
    U K = u k , K = 1 , 2 , 3 , u γ , K = 4 , , n , φ , K = n + 1 .
  • The “distortion” multifield tensor  B K l E 3 E ( n 3 ) E φ 1 E 3 or B K l E ( n + 1 ) E 3
    B K l = β k l , K = 1 , 2 , 3 , β γ l , K = 4 , , n , E l , K = n + 1 ,
    which has one multifield index and one index in the parallel space. Above, ⊗ denotes the tensor product. It is easy to see that
    B K l = U K , l .
  • The “stress” multifield tensor  Σ I j E 3 E ( n 3 ) E φ 1 E 3 or Σ I j E ( n + 1 ) E 3
    Σ I j = σ i j , I = 1 , 2 , 3 , σ α j , I = 4 , , n , D j , I = n + 1 ,
    which has one multifield index and one index in the parallel space.
  • The “body force density” multifield vector  F I = ( f i , f α , q ) E 3 E ( n 3 ) E φ 1 or F I E ( n + 1 )
    F I = f i , I = 1 , 2 , 3 , f α , I = 4 , , n , q , I = n + 1 .
  • The multifield tensor of the material moduli  C I j K l E 3 E ( n 3 ) E φ 1 E 3 E 3 E ( n 3 ) E φ 1 E 3 or C I j K l E ( n + 1 ) E 3 E ( n + 1 ) E 3
    C I j K l = C i j k l , I = 1 , 2 , 3 ; K = 1 , 2 , 3 , D i j γ l , I = 1 , 2 , 3 ; K = 4 , , n , e l i j , I = 1 , 2 , 3 ; K = n + 1 , D k l α j , I = 4 , , n ; K = 1 , 2 , 3 , E α j γ l , I = 4 , , n ; K = 4 , , n , f l α j , I = 4 , , n ; K = n + 1 , e j k l , I = n + 1 ; K = 1 , 2 , 3 , f j γ l , I = n + 1 ; K = 4 , , n , ε j l , I = n + 1 ; K = n + 1 .
    The tensor C I j K l retains a major symmetry in the hyperspace notation (see Appendix A)
    C I j K l = C K l I j .
    In matrix form, Equation (30) reads
    C I j K l = C i j k l D i j γ l e l i j D k l α j E α j γ l f l α j e j k l f j γ l ε j l .
    The multifield tensor of the material moduli C I j K l in the hyperspace notation is explicitly given for one-dimensional, two-dimensional and three-dimensional piezoelectric quasicrystals in the Appendix B.
Strictly speaking, the tensors B K l , Σ I j and C I j K l appearing in Equations (26), (28) and (30), respectively, are double tensor fields [44] or two-point tensors [45] since they have indices in two spaces, namely the hyperspace and the parallel space. Specifically, B K l and Σ I j are two-point tensors of rank two, and C I j K l is a two-point tensor of rank four.
Using the hyperspace notation for piezoelectric quasicrystals, the basic equations can be now written in a simple “compact” form, which is easy to handle. The compatibility conditions (3) and (6) are simply written
ϵ i j l B K l , j = 0 .
The constitutive relations (10)–(12) can be written in the form
Σ I j = C I j K l B K l
and the equilibrium conditions (7)–(9) read
Σ I j , j + F I = 0 .
By substituting Equation (34) into Equation (35) and using Equation (27), the coupled equations of equilibrium for the “displacement” multifield vector U K read as follows
C I j K l U K , l j + F I = 0 .
Equation (36) is the “compact” form of the coupled equations of the equilibrium for the phonon and phason displacement fields and the electrostatic potential, Equations (19)–(21), using the hyperspace notation, and it is a Navier-type partial differential equation for the “displacement” multifield vector U K . Finally, the electric enthalpy density (23) in the hyperspace notation reads as
Ψ = 1 2 B I j C I j K l B K l ,
where the major symmetry (31) can be easily seen.

4. Fundamental Properties of Piezoelectric Quasicrystals

In this section, we derive the three-dimensional Green tensor for homogeneous piezoelectric quasicrystals, which is used to calculate the “displacement”, “distortion” and “stress” multifields in the presence of a “force” multifield in a piezoelectric quasicrystal. By considering a Kelvin-type “force” multifield, the generalised Kelvin problem for piezoelectric quasicrystals is solved. The section ends with the derivation of the two-dimensional Green tensor for homogeneous piezoelectric quasicrystals. It should be mentioned that the Green functions and the corresponding Kelvin solutions serve, in general, as a pertinent basis of important numerical methods, such as the boundary element method [46].

4.1. Three-Dimensional Green Tensor of Piezoelectric Quasicrystals

The method of Green functions (see, e.g., [47]) is usually used to solve linear inhomogeneous partial differential equations like Equation (36). The Green tensor G K M ( R ) , K , M = 1 , , n + 1 of the three-dimensional Navier-type Equation (36), describing the electro-elastic fields in piezoelectric quasicrystals, is defined by
C I j K l G K M , l j ( R ) + δ I M δ ( R ) = 0 ,
where I , K , M = 1 , , n + 1 , j , l = 1 , 2 , 3 , δ I M is the ( n + 1 ) × ( n + 1 ) unit matrix in the hyperspace, R = x x , and δ ( R ) is the three-dimensional Dirac delta function. G K M ( R ) represents the “displacement” multifield in the hyperspace in the K direction at the point x , arising from a unit point “force” multifield in the M direction applied at the point x . Therefore, G K M ( R ) is also a multifield Green tensor. The Green tensor G K M ( R ) satisfies the symmetry relations
G K M ( R ) = G M K ( R ) = G K M ( R ) ,
that is, G K M ( R ) is symmetric and is an even function of R . The fact that the Green tensor G K M ( R ) is a symmetric tensor of rank two in the hyperspace is a direct consequence of the major symmetry (31).
Using the three-dimensional Fourier integral representation of the Green tensor
G K M ( R ) = 1 ( 2 π ) 3 G K M ( k ) e i k · R d k
and of the Dirac delta function
δ ( R ) = 1 ( 2 π ) 3 e i k · R d k ,
Equation (38) can be transformed to an algebraic equation in the Fourier space
C I j K l k j k l G K M = δ I M .
If we introduce the unit vector in the three-dimensional Fourier space with k R 3
κ = k / | k |
and the notation
( κ C κ ) I K = κ j C I j K l κ l = κ j C i j k l κ l κ j D i j γ l κ l κ j e l i j κ l κ j D k l α j κ l κ j E α j γ l κ l κ j f l α j κ l κ j e j k l κ l κ j f j γ l κ l κ j ε j l κ l ,
the multifield Green tensor can be written in the Fourier space as follows:
G K M ( k ) = 1 k 2 ( κ C κ ) K M 1 ,
which is a homogeneous function of k of degree 2 . The matrix ( κ C κ ) K M 1 is the inverse of ( κ C κ ) K M and is given by
( κ C κ ) K M 1 = A K M ( κ ) D ( κ ) ,
where D ( κ ) and A K M ( κ ) are the determinant and the adjoint of the matrix ( κ C κ ) K M , respectively. Substituting Equation (45) into Equation (40), the three-dimensional Fourier integral can be reduced to a line integral along the unit circle in the plane orthogonal to R (see, e.g., [9,11,33])
G K M ( k ) e i k · R d k = π R 0 2 π G K M ( n ) d ϕ .
Here, n is a unit vector “scanning” the circle of integration and remaining orthogonal to R thus n · R = 0 and R = | R | . Note that n is a function of ϕ . Combining Equations (40), (45) and (47), the ( n + 1 ) × ( n + 1 ) Green tensor G K M ( R ) of piezoelectric quasicrystals can be written in the form
G K M ( R ) = 1 8 π 2 R 0 2 π ( n C n ) K M 1 d ϕ ,
which represents the three-dimensional Green tensor of piezoelectric quasicrystals. One can see that the Green tensor (48) is the generalisation of the three-dimensional Green tensors of the anisotropic elasticity (see, e.g., [9,10,11,48,49]), of piezoelectricity (see, e.g., [19,20,21]) and of the elasticity of quasicrystals [33] towards the piezoelectricity of quasicrystals. The integral in Equation (48) can be computed by standard numerical methods if C I j K m is given.
The Green tensor G K M ( R ) of piezoelectric quasicrystals can be decomposed into its phonon, phason and electric parts, and their coupling
G K M ( R ) = G k m ( R ) G k μ ( R ) G k n + 1 ( R ) G γ m ( R ) G γ μ ( R ) G k n + 1 ( R ) G n + 1 m ( R ) G n + 1 μ ( R ) G n + 1 n + 1 ( R )
with the symbol denoting the electric contribution, and it has the following physical interpretation:
  • G k m ( R ) is the phonon displacement at x in the e k caused by a unit phonon point force at x in the e m direction, k , m = 1 , 2 , 3 ;
  • G k μ ( R ) is the phonon displacement at x in the e k direction caused by a unit phason point force at x in the e μ direction, k = 1 , 2 , 3 , μ = 4 , , n ;
  • G k n + 1 ( R ) is the phonon displacement at x in the e k direction caused by a unit electric point charge at x , k = 1 , 2 , 3 ;
  • G γ m ( R ) is the phason displacement at x in the e γ direction caused by a unit phonon point force at x in the e m direction, m = 1 , 2 , 3 , γ = 4 , , n ;
  • G γ μ ( R ) is the phason displacement at x in the e γ direction caused by a unit phason point force at x in the e μ direction, γ , μ = 4 , , n ;
  • G γ n + 1 ( R ) is the phason displacement at x in the e γ direction caused by a unit electric point charge at x , γ = 4 , , n ;
  • G n + 1 m ( R ) is the electrostatic potential at x caused by a unit phonon point force at x in the e m direction, m = 1 , 2 , 3 ;
  • G n + 1 μ ( R ) is the electrostatic potential at x caused by a unit phason point force at x in the e μ direction, μ = 4 , , n ;
  • G n + 1 n + 1 ( R ) is the electrostatic potential at x caused by a unit electric point charge at x .

4.2. “Displacement”, “Distortion” and “Stress” Multifields in the Presence of a “Force” Multifield

The three-dimensional Green tensor (48) of piezoelectric quasicrystals can be applied to calculate the phonon, phason and electric fields caused by phonon and phason body forces and a body charge in homogeneous piezoelectric quasicrystalline materials.
For an arbitrary “body force density” multifield vector F M (see Equation (29)), the solution of the partial differential Equation (36) can be written as
U K = G K M F M ,
where ∗ denotes the three-dimensional spatial convolution. Using the Green tensor (48), Equation (50) gives the “displacement” multifield vector caused by a “body force density” multifield in a piezoelectric quasicrystal
U K ( x ) = 1 8 π 2 V F M ( x ) R 0 2 π ( n C n ) K M 1 d ϕ d V .
It is worth noticing that Equation (51) represents the generalisation towards the piezoelectricity of the quasicrystals in the following problems:
  • The solution of the displacement for an arbitrary body force density known from anisotropic elasticity [9];
  • The solution of the electrostatic potential for an arbitrary body charge density known from electrostatics [41];
  • The solution of the “extended” displacement for an arbitrary “extended” force density known from piezoelectricity [21].
Combining Equations (27) and (50), the “distortion” multifield tensor reads in terms of the gradient of the Green tensor (48)
B K l = G K M , l F M ,
where the spatial gradient of the Green tensor is given by (see also [33])
G K M , l ( R ) = 1 8 π 2 R 2 0 2 π τ l ( n C n ) K M 1 n l F K M d ϕ
with
F K M = ( n C n ) K N 1 ( n C τ ) N P + ( τ C n ) N P ( n C n ) P M 1
and τ = R / R is the unit vector along R . Substituting Equation (53) into Equation (52), we obtain the “distortion” multifield tensor caused by a“body force density” multifield in a piezoelectric quasicrystal
B K l ( x ) = 1 8 π 2 V F M ( x ) R 2 0 2 π τ l ( n C n ) K M 1 n l F K M d ϕ d V .
Using the constitutive relation (34), the corresponding“stress”multifield tensor for piezoelectric quasicrystals reads
Σ I j ( x ) = 1 8 π 2 V F M ( x ) R 2 0 2 π C I j K l τ l ( n C n ) K M 1 n l F K M d ϕ d V .

Generalised Kelvin Problem in Piezoelectric Quasicrystals

Here, we examine the generalised Kelvin problem towards piezoelectric quasicrystals. We specify a Kelvin-type “force” density multifield vector with strength f ^ M
F M ( x ) = f ^ M δ ( x ) ,
which is a point “force” density multifield in three dimensions. Then, Equation (50) becomes
U K ( x ) = f ^ M G K M
and Equation (51) reduces to
U K ( x ) = f ^ M 8 π 2 r 0 2 π ( n C n ) K M 1 d ϕ ,
where r = | x | = x 2 + y 2 + z 2 . Equation (59) is the Kelvin-type “force” multifield solution for the “displacement” multifield vector for piezoelectric quasicrystals. In this case, Equation (55) reduces to
B K l ( x ) = f ^ M 8 π 2 r 2 0 2 π τ l ( n C n ) K M 1 n l F K M d ϕ ,
which is the Kelvin-type “force” multifield solution for the “distortion” multifield tensor for piezoelectric quasicrystals. Analogously, Equation (56) yields to
Σ I j ( x ) = f ^ M 8 π 2 r 2 0 2 π C I j K l τ l ( n C n ) K M 1 n l F K M d ϕ ,
which is the Kelvin-type “force” multifield solution for the “stress” multifield tensor for piezoelectric quasicrystals.
It is interesting to examine that for some particular cases, the solution (59) leads to well-known results from classical anisotropic elasticity, isotropic elasticity and electrostatics. Indeed:
  • If the phason modes and the electric effects are absent, that is, considering only classical anisotropic elasticity, then Equation (59) reads
    u k ( x ) = f ^ m 8 π 2 r 0 2 π ( n C n ) k m 1 d ϕ ,
    which is the displacement vector caused by a point force with strength f ^ m in anisotropic elasticity (see, e.g., [9]).
    Next, if we consider additionally that the medium is isotropic with the tensor of elastic moduli to be given by
    C i j k l = 2 μ ν 1 2 ν δ i j δ k l + μ δ i k δ j l + δ i l δ j k ,
    where μ is the shear modulus and ν is the Poisson ratio, then Equation (62) reduces to
    u k ( x ) = f ^ m 16 π μ ( 1 ν ) r ( 3 4 ν ) δ k m + x k x m r 2 ,
    which is the well-known Kelvin solution of a point force with strength f ^ m in an infinite isotropic body (see, e.g., [5,50]).
  • If only electric effects are considered, then Equation (59) gives
    φ ( x ) = q ^ 8 π 2 r 0 2 π ( n ε n ) 1 d ϕ ,
    which is an integral representation of the electrostatic potential φ caused by a point charge with strength q ^ in anisotropic electrostatics.
    Next, if we consider the tensor of the dielectric permittivity ε i j = δ i j ε 0 with ε 0 the permittivity of the free space, then Equation (65) reduces to the well-known formula of the electrostatic potential of a point charge with strength q ^ (see, e.g., [41])
    φ ( x ) = q ^ 4 π ε 0 1 r .

4.3. Two-Dimensional Green Tensor of Piezoelectric Quasicrystals

In this section, the two-dimensional Green tensor for arbitrary piezoelectric quasicrystals is derived. Moreover, the two-dimensional Green tensor for arbitrary quasicrystals in the framework of elasticity is given here for the first time.
The Green tensor G K M ( R ) with K , M = 1 , , n + 1 of the two-dimensional Navier-type equation is defined by
C I j K l G K M , l j ( R ) + δ I M δ ( R ) = 0 ,
where I , K , M = 1 , , n + 1 , j , l = 1 , 2 , R = x x R 2 , and δ ( R ) is the two-dimensional Dirac delta function. The two-dimensional Fourier transform of the Green tensor in Equation (67) reads ( k R 2 )
G K M ( k ) = 1 k 2 ( κ C κ ) K M 1 ,
where κ denotes the unit vector in the k 1 k 2 -plane of the two-dimensional Fourier space ( κ j with j = 1 , 2 ) defined by κ = k / k with k = | k | . Here, the matrix ( κ C κ ) K M 1 is the inverse of ( κ C κ ) K M given in Equation (44) with j , l = 1 , 2 .
The two-dimensional Green tensor in the two-dimensional real space is obtained by the inverse Fourier transform of Equation (68) as follows:
G K M ( R ) = 1 ( 2 π ) 2 R 2 ( κ C κ ) K M 1 cos ( k · R ) k 2 d k = 1 ( 2 π ) 2 0 2 π ( κ C κ ) K M 1 0 cos ( k κ · R ) k d k d ϕ .
In Equation (69), d k = k d k d ϕ indicates the two-dimensional volume element in Fourier space in polar coordinates, and ϕ ( 0 < ϕ 2 π ) is an appropriate polar angle scanning a unit circle κ 2 = 1 . The integration in k is performed using the relation
0 cos ( k κ · R ) k d k = 1 2 1 k e i k κ · R d k
and the principal value integral [47]
P 1 k e i k t d k = 2 γ E 2 ln | t | ,
where γ E denotes the Euler constant and P means the principal value. The constant term can be neglected for the Green tensor (see, e.g., [47]). Using Equations (70) and (71), the two-dimensional Green tensor (69) can be expressed as integral over the unit circle in Fourier space
G K M ( R ) = 1 ( 2 π ) 2 0 2 π ( κ C κ ) K M 1 ln | κ · R | d ϕ .
Because κ ( ϕ ) varies with ϕ , the integrand ( κ C κ ) K M 1 is a function of the integration variable ϕ and n κ ( ϕ ) , and the ( n + 1 ) × ( n + 1 ) Green tensor G K M ( R ) becomes
G K M ( R ) = 1 ( 2 π ) 2 0 2 π ( n C n ) K M 1 ln | n · R | d ϕ ,
which represents the two-dimensional Green tensor of piezoelectric quasicrystals.
If the piezoelectric effects are not considered, that is, if only the elasticity of quasicrystals is considered, then the piezoelectric ( n + 1 ) × ( n + 1 ) Green tensor G K M ( R ) , K , M = 1 , , n + 1 , Equation (73), reduces to the elastic n × n Green tensor of quasicrystals
G K M ( R ) = 1 ( 2 π ) 2 0 2 π ( n C n ) K M 1 ln | n · R | d ϕ
with K , M = 1 , , n . Equation (74) represents the two-dimensional Green tensor of quasicrystals.
In the limit to piezoelectricity, that is, if the phason modes are absent, then the Green tensor (73) reduces to the two-dimensional Green tensor of piezoelectricity
G K M ( R ) = 1 ( 2 π ) 2 0 2 π ( n C n ) K M 1 ln | n · R | d ϕ
with K , M = 1 , , 4 , which is in agreement with the two-dimensional Green tensor of piezoelectricity given by Wang [23].
In the limit to anisotropic elasticity, that is, if the phason modes and electric fields are absent, then the Green tensor (73) reduces to the two-dimensional Green tensor of anisotropic elasticity
G k m ( R ) = 1 ( 2 π ) 2 0 2 π ( n C n ) k m 1 ln | n · R | d ϕ
with k , m = 1 , 2 , 3 . Note that the Green tensor (76) is in agreement with the two-dimensional Green tensor of anisotropic elasticity given by Wang and Achenbach [12], Wang [13] and Lazar [14].

5. Conclusions

In this work, the three-dimensional Green tensor and its spatial gradient have been calculated for arbitrary piezoelectric quasicrystals. With this aim, the hyperspace notation for piezoelectric quasicrystals has been developed, and the corresponding multifields have been defined. Using the three-dimensional Green tensor, we have given a solution for the “displacement”, “distortion” and “stress” multifields in the presence of a “force” multifield in a piezoelectric quasicrystal as well as the solution of the generalised Kelvin problem. These solutions allow for the consideration of all elastic, dielectric and piezoelectric material properties in piezoelectric quasicrystals. Moreover, the two-dimensional Green tensors of piezoelectric quasicrystals and of quasicrystals have been calculated. The derived Green tensors may be used for the study of dislocation, fracture, interface and inclusion problems in piezoelectric quasicrystals.

Author Contributions

Conceptualization, M.L. and E.A.; methodology, M.L. and E.A.; writing—original draft preparation, M.L. and E.A.; writing—review and editing, M.L. and E.A.; funding acquisition, M.L. and E.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deutsche Forschungsgemeinschaft, Grant Numbers LA1974/4-2, AG322/1-2.

Data Availability Statement

The data are contained within the article.

Acknowledgments

Markus Lazar and Eleni Agiasofitou gratefully acknowledge the grants from the Deutsche Forschungsgemeinschaft (Grant numbers LA1974/4-2 and AG322/1-2).

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. The Symmetry of the Multifield Tensor of the Material Moduli CIjKl

In this Appendix, the major symmetry (see Equation (31)) of the multifield tensor of the material moduli C I j K l in the hyperspace notation is proved. Indeed,
C K l I j = C k l i j , K = 1 , 2 , 3 ; I = 1 , 2 , 3 , D k l α j , K = 1 , 2 , 3 ; I = 4 , , n , e j k l , K = 1 , 2 , 3 ; I = n + 1 , D i j γ l , K = 4 , , n ; I = 1 , 2 , 3 , E γ l α j , K = 4 , , n ; I = 4 , , n , f j γ l , K = 4 , , n ; I = n + 1 , e l i j , K = n + 1 ; I = 1 , 2 , 3 , f l α j , K = n + 1 ; I = 4 , , n , ε l j , K = n + 1 ; I = n + 1 ,
which corresponds to the electric enthalpy density
Ψ = 1 2 C k l i j e k l e i j + 1 2 D k l α j e k l β α j 1 2 e j k l E j e k l + 1 2 D i j γ l e i j β γ l + 1 2 E γ l α j β γ l β α j 1 2 f j γ l E j β γ l 1 2 e l i j E l e i j 1 2 f l α j E l β α j 1 2 ε l j E l E j .
Making use only of the major symmetries C i j k l = C k l i j , E γ l α j = E α j γ l and ε l j = ε j l , it is easy to see that
C K l I j = C i j k l , I = 1 , 2 , 3 ; K = 1 , 2 , 3 , D k l α j , I = 4 , , n ; K = 1 , 2 , 3 , e j k l , I = n + 1 ; K = 1 , 2 , 3 , D i j γ l , I = 1 , 2 , 3 ; K = 4 , , n , E α j γ l , I = 4 , , n ; K = 4 , , n , f j γ l , I = n + 1 ; K = 4 , , n , e l i j , I = 1 , 2 , 3 ; K = n + 1 , f l α j , I = 4 , , n ; K = n + 1 , ε j l , I = n + 1 ; K = n + 1 , = C I j K l .

Appendix B. The Multifield Tensor of the Material Moduli CIjKl

In this Appendix, we give the multifield tensor of the material moduli C I j K l in the hyperspace notation explicitly for one-dimensional, two-dimensional and three-dimensional piezoelectric quasicrystals.
(i)
For one-dimensional piezoelectric quasicrystals
C I j K l = C i j k l , I = 1 , 2 , 3 ; K = 1 , 2 , 3 , D i j γ l , I = 1 , 2 , 3 ; K = 4 , e l i j , I = 1 , 2 , 3 ; K = 5 , D k l α j , I = 4 ; K = 1 , 2 , 3 , E α j γ l , I = 4 ; K = 4 , f l α j , I = 4 ; K = 5 , e j k l , I = 5 ; K = 1 , 2 , 3 , f j γ l , I = 5 ; K = 4 , ε j l , I = 5 ; K = 5 ,
where i , j , k , l = 1 , 2 , 3 and α , γ = 3 .
(ii)
For two-dimensional piezoelectric quasicrystals,
C I j K l = C i j k l , I = 1 , 2 , 3 ; K = 1 , 2 , 3 , D i j γ l , I = 1 , 2 , 3 ; K = 4 , 5 , e l i j , I = 1 , 2 , 3 ; K = 6 , D k l α j , I = 4 , 5 ; K = 1 , 2 , 3 , E α j γ l , I = 4 , 5 ; K = 4 , 5 , f l α j , I = 4 , 5 ; K = 6 , e j k l , I = 6 ; K = 1 , 2 , 3 , f j γ l , I = 6 ; K = 4 , 5 , ε j l , I = 6 ; K = 6 ,
where i , j , k , l = 1 , 2 , 3 and α , γ = 1 , 2 .
(iii)
For three-dimensional piezoelectric quasicrystals,
C I j K l = C i j k l , I = 1 , 2 , 3 ; K = 1 , 2 , 3 , D i j γ l , I = 1 , 2 , 3 ; K = 4 , 5 , 6 , e l i j , I = 1 , 2 , 3 ; K = 7 , D k l α j , I = 4 , 5 , 6 ; K = 1 , 2 , 3 , E α j γ l , I = 4 , 5 , 6 ; K = 4 , 5 , 6 , f l α j , I = 4 , 5 , 6 ; K = 7 , e j k l , I = 7 ; K = 1 , 2 , 3 , f j γ l , I = 7 ; K = 4 , 5 , 6 , ε j l , I = 7 ; K = 7 ,
where i , j , k , l = 1 , 2 , 3 and α , γ = 1 , 2 , 3 .

References

  1. Cook-Chennault, K.A.; Thambi, N.; Sastry, A.M. Powering MEMS portable devices—A review of non-regenerative and regenerative power supply systems with special emphasis on piezoelectric energy harvesting systems. Smart Mater. Struct. 2008, 17, 043001. [Google Scholar] [CrossRef]
  2. Schouten, J.A. Tensor Analysis for Physicists; Oxford University Press: Oxford, UK, 1951. [Google Scholar]
  3. Nye, J.F. Physical Properties of Crystals; Oxford University Press: Oxford, UK, 1957. [Google Scholar]
  4. Landau, L.D.; Lifshitz, E.M. Electrodynamics of Continuous Media; Pergamon Press: Oxford, UK, 1960. [Google Scholar]
  5. Kelvin, L. Mathematical and Physical Papers; Cambridge University Press: Cambridge, UK, 1882; Volume 1, p. 97. [Google Scholar]
  6. Lifshitz, I.M.; Rosenzweig, L.N. On the construction of the Green’s tensor for the basic equation of the theory of elasticity of an anisotropic medium. Zh. Eksper. Teor. Fiz. 1947, 17, 783–791. [Google Scholar]
  7. Kröner, E. Das Fundamentalintegral der anisotropen elastischen Differentialgleichungen. Z. Phys. 1953, 136, 402–410. [Google Scholar] [CrossRef]
  8. Willis, J.R. The elastic interaction energy of dislocation loops in anisotropic media. Quart. J. Mech. Appl. Math. 1965, 18, 419–433. [Google Scholar] [CrossRef]
  9. Synge, J.L. The Hypercircle in Mathematical Physics; Cambridge University Press: Cambridge, UK, 1957. [Google Scholar]
  10. Barnett, D.M. The precise evaluation of derivatives of the anisotropic elastic Green’s functions. Phys. Status Solidi B 1972, 49, 741–748. [Google Scholar] [CrossRef]
  11. Teodosiu, C. Elastic Models of Crystal Defects; Springer: Berlin/Heidelberg, Germany, 1982. [Google Scholar]
  12. Wang, C.-Y.; Achenbach, J.D. Elastodynamic fundamental solutions for anisotropic solids. Geophys. J. Int. 1994, 118, 384–392. [Google Scholar] [CrossRef]
  13. Wang, C.-Y. Two-dimensional elastostatic Green’s functions for general anisotropic solids and generalization of Stroh’s formalism. Int. J. Solids Struct. 1994, 31, 2591–2597. [Google Scholar] [CrossRef]
  14. Lazar, M. Displacements and stress functions of straight dislocations and line forces in anisotropic elasticity: A new derivation and its relation to the integral formalism. Symmetry 2021, 13, 1721. [Google Scholar] [CrossRef]
  15. Deeg, W.F.J. The analysis of dislocation, crack, and inclusion problems in piezoelectric solids. Ph.D. Thesis, Stanford University, Stanford, CA, USA, 1980. [Google Scholar]
  16. Chen, T. Green’s functions and the non-uniform transformation problem in a piezo electric medium. Mech. Res. Commun. 1993, 20, 271–278. [Google Scholar] [CrossRef]
  17. Dunn, M.L. Electroelastic Green’s functions for transversely isotropic piezoelectric media and their application to the solution of inclusion and inhomogeneity problems. Int. J. Eng. Sci. 1994, 32, 119–131. [Google Scholar] [CrossRef]
  18. Akamatsu, M.; Tanuma, K. Green’s function of anisotropic piezoelectricity. Proc. R. Soc. Lond. A 1997, 453, 473–487. [Google Scholar] [CrossRef]
  19. Pan, E.; Tonon, F. Three-dimensional Green’s functions in anisotropic piezoelectric solids. Int. J. Solids Struct. 2000, 37, 943–958. [Google Scholar] [CrossRef]
  20. Li, X.; Wang, M. Three-dimensional Green’s functions for infinite anisotropic piezoelectric media. Int. J. Solids Struct. 2007, 44, 1680–1684. [Google Scholar] [CrossRef]
  21. Nowacki, J. Static and Dynamic Coupled Fields in Bodies with Piezoeffects or Polarization Gradient; Lecture Notes in Applied and Computational Mechanics; Springer: Berlin/Heidelberg, Germany, 2006. [Google Scholar]
  22. Michelitsch, T. Calculation of the electroelastic Green’s function of the hexagonal infinite medium. Z. Phys. B 1997, 104, 497–503. [Google Scholar] [CrossRef]
  23. Wang, C.-Y. Green’s functions and general formalism for 2D piezoelectricity. Appl. Math. Lett. 1996, 9, 1–7. [Google Scholar] [CrossRef]
  24. Shechtman, D.; Blech, I.; Gratias, D.; Cahn, J.W. Metallic phase with long-range orientational order and no translational symmetry. Phys. Rev. Lett. 1984, 53, 1951–1953. [Google Scholar] [CrossRef]
  25. Bak, P. Symmetry, stability, and elastic properties of icosahedral incommensurate crystals. Phys. Rev. B 1985, 32, 5764–5772. [Google Scholar] [CrossRef]
  26. Levine, D.; Lubensky, T.C.; Ostlund, S.; Ramaswamy, S.; Steinhardt, P.J.; Toner, J. Elasticity and dislocations in pentagonal and icosahedral quasicrystals. Phys. Rev. Lett. 1985, 54, 1520–1523. [Google Scholar] [CrossRef]
  27. Hu, K.; Meguid, S.A.; Wang, L.; Jin, H. Electro-elastic field of a piezoelectric quasicrystal medium containing two cylindrical inclusions. Acta Mech. 2021, 232, 2513–2533. [Google Scholar] [CrossRef]
  28. Wang, Y.; Liu, C.; Zhang, L.; Pan, E.; Gao, Y. Mechanical analysis of functionally graded multilayered two-dimensional decagonal piezoelectric quasicrystal laminates with imperfect interfaces. Crystals 2024, 14, 170. [Google Scholar] [CrossRef]
  29. Wang, Y.; Liu, C.; Zhu, Z.; Zhang, L.; Gao, Y. Analysis of multilayered two-dimensional decagonal piezoelectric quasicrystal beams with mixed boundary conditions. Eur. J. Mech. A-Solids 2024, 106, 105333. [Google Scholar] [CrossRef]
  30. De, P.; Pelcovits, R.A. Linear elasticity theory of pentagonal quasicrystals. Phys. Rev. B 1987, 35, 8609–8620. [Google Scholar] [CrossRef] [PubMed]
  31. Ding, D.-H.; Wang, R.; Yang, W.; Hu, C. General expressions for the elastic displacement fields induced by dislocations in quasicrystals. J. Phys. Condens. Matter 1995, 7, 5423–5436. [Google Scholar] [CrossRef]
  32. Bachteler, J.; Trebin, H.-R. Elastic Green’s function of icosahedral quasicrystals. Eur. Phys. J. B 1998, 4, 299–306. [Google Scholar] [CrossRef]
  33. Lazar, M.; Agiasofitou, E. Fundamentals in generalized elasticity and dislocation theory of quasicrystals: Green tensor, dislocation key-formulas and dislocation loops. Philos. Mag. 2014, 94, 4080–4101. [Google Scholar] [CrossRef]
  34. Yang, W.; Wang, R.; Ding, D.-H.; Hu, C. Elastic strains induced by electric fields in quasicrystals. J. Phys. Condens. Matter 1995, 7, L499–L502. [Google Scholar] [CrossRef]
  35. Hu, C.; Wang, R.; Ding, D.-H.; Yang, W. Piezoelectric effects in quasicrystals. Phys. Rev. B 1997, 56, 2463–2468. [Google Scholar] [CrossRef]
  36. Agiasofitou, E.; Lazar, M. On the constitutive modelling of piezoelectric quasicrystals. Crystals 2023, 13, 1652. [Google Scholar] [CrossRef]
  37. Li, X.Y.; Li, P.D.; Wu, T.H.; Shi, M.X.; Zhu, Z.W. Three-dimensional fundamental solutions for one-dimensional hexagonal quasicrystal with piezoelectric effect. Phys. Lett. A 2014, 378, 826–834. [Google Scholar] [CrossRef]
  38. Fan, C.Y.; Li, Y.; Xu, G.T.; Zhao, M.H. Fundamental solutions and analysis of three-dimensional cracks in one-dimensional hexagonal piezoelectric quasicrystals. Mech. Res. Commun. 2016, 74, 39–44. [Google Scholar] [CrossRef]
  39. Zhang, L.; Wu, D.; Xu, W.; Yang, L.; Ricoeur, A.; Wanf, Z.; Gao, Y. Green’s functions of one-dimensional quasicrystal bi-material with piezoelectric effect. Phys. Lett. A 2016, 380, 3222–3228. [Google Scholar] [CrossRef]
  40. Ding, D.-H.; Wang, W.; Hu, C.; Yang, R. Generalized elasticity theory of quasicrystals. Phys. Rev. B 1993, 48, 7003–7010. [Google Scholar] [CrossRef] [PubMed]
  41. Jackson, J.D. Classical Electrodynamics, 2nd ed.; Wiley: New York, NY, USA, 1975. [Google Scholar]
  42. Lazar, M.; Agiasofitou, E. Eshelbian mechanics of novel materials: Quasicrystals. J. Micromech. Mol. Phys. 2016, 1, 1640008. [Google Scholar] [CrossRef]
  43. Barnett, D.M.; Lothe, J. Dislocations and line charges in anisotropic piezoelectric insulators. Phys. Status Solidi B 1975, 67, 105–111. [Google Scholar] [CrossRef]
  44. Ericksen, J.L. Appendix, Tensor Fields. In Encyclopedia of Physics Vol. III/1; Flügge, S., Ed.; Springer: Berlin/Heidelberg, Germany, 1960; pp. 794–858. [Google Scholar]
  45. Marsden, J.E.; Hughes, T.J.S. Mathematical Foundations of Elasticity; Prentice Hall: Upper Saddle River, NJ, USA, 1982. [Google Scholar]
  46. Becker, A.A. The Boundary Element Method in Engineering: A Complete Course; Mcgraw-Hill: New York, NY, USA, 1992. [Google Scholar]
  47. Vladimirov, V.S. Equations of Mathematical Physics; Marcel Dekker, Inc.: New York, NY, USA, 1971. [Google Scholar]
  48. Lazar, M.; Kirchner, H.O.K. Dislocation loops in anisotropic elasticity: Displacement field, stress function tensor and interaction energy. Philos. Mag. 2013, 93, 174–185. [Google Scholar] [CrossRef]
  49. Li, S.; Wang, G. Introduction to Micromechanics and Nanomechanics; World Scientific: Singapore, 2008. [Google Scholar]
  50. Gurtin, M.E. The linear theory of elasticity. In Handbuch der Physik VIa/2; Flügge, S., Ed.; Springer: Berlin/Heidelberg, Germany, 1972; pp. 1–345. [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lazar, M.; Agiasofitou, E. Three-Dimensional and Two-Dimensional Green Tensors of Piezoelectric Quasicrystals. Crystals 2024, 14, 835. https://doi.org/10.3390/cryst14100835

AMA Style

Lazar M, Agiasofitou E. Three-Dimensional and Two-Dimensional Green Tensors of Piezoelectric Quasicrystals. Crystals. 2024; 14(10):835. https://doi.org/10.3390/cryst14100835

Chicago/Turabian Style

Lazar, Markus, and Eleni Agiasofitou. 2024. "Three-Dimensional and Two-Dimensional Green Tensors of Piezoelectric Quasicrystals" Crystals 14, no. 10: 835. https://doi.org/10.3390/cryst14100835

APA Style

Lazar, M., & Agiasofitou, E. (2024). Three-Dimensional and Two-Dimensional Green Tensors of Piezoelectric Quasicrystals. Crystals, 14(10), 835. https://doi.org/10.3390/cryst14100835

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop