## Abstract

Owing to their specific structure, which can neither be classified as crystalline nor amorphous, quasi-crystals (QCs) exhibit properties that are interesting to both material science and mathematical physics or continuum mechanics. Within the framework of a mathematical theory of elasticity, one major focus is on features evolving from the coupling of phonon and phason fields, which is not observed in classical crystalline or amorphous materials. This paper deals with the problems of combinations of point phonon forces and point phason forces, which are applied to the interior of infinite solids and bimaterial solids of two-dimensional hexagonal QCs. By using the general solution of QCs, a series of displacement functions is adopted to obtain the analytical results when the two half-spaces are supposed to be ideally bonded or to be in smooth contact. In the final expressions, we provide three-dimensional Green’s functions for infinite bimaterial QC solids in the closed form, which are very convenient to be used in the study of dislocations, cracks and inhomogeneities of the new solid phase. Furthermore, the paper is concluded by a discussion of some special cases, in which Green’s functions for infinite transversely isotropic solids and Green’s functions for a half-space with free or fixed boundary are given.

## 1. Introduction

A two-dimensional quasi-crystal (QC) is defined as a three-dimensional body in which the atom arrangement is quasi-periodic in a plane and periodic in the orthogonal direction. Therefore, there are two kinds of displacement fields in the elastic theory of QCs. One is a phonon displacement field, corresponding to the displacement field of classical crystals. The other is the phason displacement field, which is diffusive due to the elementary excitation associated with the phason mode and describes the local rearrangements of the unit cells. Both types of fields are coupled with each other. Accordingly, there are phonon and phason stresses and strains, the phason appearing just in QC elasticity.

From the theoretical point of view, QCs are an interesting subject of investigation since they exhibit a clearly defined anisotropy and the phenomenon of field coupling. It is a fundamental concern to extend the theory of elasticity towards this extraordinary class of solids. From the technological point of view, there is a variety of one-and two-dimensional quasi-crystalline materials (Gao & Ricoeur 2011, submitted). For engineering applications, a theoretical basis for stress and strain analyses has to be elaborated. Here, Green’s functions supply a perfect background since arbitrary loading situations can be reproduced by superposition. Hence, even cracks and other defects can be included in a stress analysis based on Green’s functions. This is particularly important to assess fracture loads and the lifetime of engineering structures composed of quasi-crystalline materials.

QCs have become the focus of theoretical and experimental studies in the physics of condensed matter since the first discovery of the icosahedral QC in Al–Mn alloys (Levine & Steinhardt 1984; Shechtman *et al.* 1984). The physical properties, including elasticity and defects of QCs have been intensively investigated in experimental and theoretical analyses (Wollgarten *et al.* 1993; Athanasiou *et al.* 2002; Park *et al.* 2005*a*,*b*). In particular, the field of linear elastic theory of QCs has been formulated for many years (Levine *et al.* 1985; Socolar 1989; Ding *et al.* 1993; Hu *et al.* 1996). Great progress has been made in the fields of the mechanics involving the elasticity and defects, see review articles for details (Hu *et al.* 2000; Fan & Mai 2004).

Many engineering structures are made by binding together two or more materials with different physical properties. The elastic body joined with two dissimilar materials with point forces applied at an arbitrary point is fundamental to the development of three-dimensional elastic theory and is of vital significance in the structural design. For isotropic materials, Vijayakumar & Cormack (1987) and Huang & Wang (1991) studied the fundamental solutions for the case when a point force is applied at one of the two bonded semi-infinite solids. With regard to transversely isotropic elastic materials, Pan & Chou (1979) conducted a systematic study on three-dimensional Green’s functions for point forces applied to the interior of a two-phase infinite space. In the Fourier transformed domain, Ting (1996) obtained Green’s functions of point forces in a two-phase anisotropic elastic solid. In the case of transversely isotropic piezoelectric materials, Ding *et al.* (1997*a*) developed effectively a novel method determining the closed-form point force and point charge solutions for two-phase piezoelectric media using the analysis technique of the image source. In a similar way, they considered plane problems (Ding *et al.* 1997*c*). For plane problems of cubic QCs with imperfect interface, Gao & Ricoeur (2010) obtained two-dimensional Green’s functions for line forces applied at the interior of a two-phase infinite plane. The significance of fundamental solutions and Green’s functions in constructing solutions to various kinds of boundary-value problems has been well recognized in the mechanics literature (Eshelby 1957; Mura 1987; Ting 1996, 2000).

Obviously, three-dimensional Green’s functions play an important role in the analysis, because they not only have theoretical merits themselves, but also can be benchmarks to clarify various approximate methods, such as the finite element method and the boundary element method as well as in the study of cracks, defects and inclusions. However, relevant three-dimensional Green’s functions for two-dimensional QCs have not been attempted. The purpose of this paper is to develop the previous work (Gao & Ricoeur 2010), and to study the problem of point forces applied to the interior of an infinite bimaterial two-dimensional hexagonal QC. To achieve this, the general solution of two-dimensional hexagonal QC (Gao & Zhao 2009) is used to uncouple the system of equations of equilibrium. Similar to the derivation of Green’s functions for transversely isotropic elastic and piezoelectric materials (Ding *et al.* 1997*a*,*b*,*c*), Green’s functions for infinite QC solids and two half-spaces with point forces are obtained in the closed-form by generalizing the analysis technique of the image source to two-dimensional QC media.

## 2. Basic equations and the general solution

In a fixed rectangular coordinate system (*x*_{1}, *x*_{2}, *x*_{3}), a two-dimensional QC refers to a three-dimensional solid structure with quasi-periodic arrangement in a plane (*x*_{1}–*x*_{2} plane) and periodic arrangement in the orthogonal direction (*x*_{3}-direction). From two-dimensional QC elastic theory (Hu *et al.* 1996), in the absence of body forces, the general equations governing the three-dimensional theory of two-dimensional hexagonal QC can be written as:
2.1
2.2
2.3where the subscripts *i*, *p*, *q*=1, 2, 3, *j*=4, 5 and *α*, *β*=1, 2 are used throughout this paper. *u*_{p} and *w*_{q} denote phonon and phason displacements in the physical and perpendicular spaces, respectively; *σ*_{pq} and *ε*_{pq} are phonon stresses and strains, respectively; *H*_{αq} and *w*_{αq} represent phason stresses and strains, respectively; *C*_{kl}, *K*_{k} and *R*_{k} stand for the elastic constants in the phonon, phason and phonon–phason coupling fields, respectively, with the relationships 2*C*_{66}=*C*_{11}−*C*_{12}, *K*_{6}=*K*_{1}−*K*_{2}−*K*_{3}, 2*R*_{6}=*R*_{1}−*R*_{2}.

For the sake of compactness, the notations *U*_{1β}=*u*_{β}, *U*_{2β}=*w*_{β}, *T*_{11}=*σ*_{11}, *T*_{21}=*H*_{11}, *T*_{12}=*σ*_{22}, *T*_{22}=*H*_{22}, *T*_{13}=*σ*_{23}, *T*_{23}=*H*_{23}, *T*_{14}=*σ*_{13} and *T*_{24}=*H*_{13} are applied in the present paper. According to the general solution of two-dimensional hexagonal QC with distinct eigenvalues (Gao & Zhao 2009), the components of displacements take the form:
2.4The constants *k*_{αi}, *k*_{αj} and *k*_{3i} are defined as
where *δ*_{kl} is the Kronecker delta symbol, and the following summation convention has been used throughout this paper: the Einstein summation over repeated lower case indices is applied, while upper case indices take on the same numbers as the corresponding lower case ones but are not summed. Besides, the potential functions *ψ*_{i} and *ψ*_{j} satisfy the equations
2.5in which is the planar Laplacian; are three eigenvalues of the following cubic algebraic equation of *s*^{2}, *as*^{6}−*bs*^{4}+*cs*^{2}−*d*=0; and are two eigenvalues of the following quadratic algebraic equation of *s*^{2}, , where
The components of stresses obtained from equations (2.3) and (2.4) can be shown to be
2.6where
For the sake of brevity and conciseness, in the following two sections the fundamental solutions for an infinite QC solid and an infinite bimaterial QC solid will be given only for the case of distinct eigenvalues. When equal eigenvalues appear, the fundamental solutions can be obtained by using a similar analysis technique, although for these cases the general solution will take a more complicated form (Gao & Zhao 2009).

## 3. Three-dimensional Green’s functions for infinite QC solids

A two-dimensional hexagonal QC full-space is considered as the domain of the problem. As indicated in figure 1, a point phonon force with components *Q*_{1}, *P*_{1} and *F* in *x*_{1}, *x*_{2} and *x*_{3}-directions, respectively, and a point phason force with components *Q*_{2}, *P*_{2} in *x*_{1} and *x*_{2}-directions, respectively are applied simultaneously at an arbitrary point in the space. Without loss in generality, we take this point as the origin of Cartesian coordinates. Based on the theorem of superposition the problem can be divided into two sub-problems: the problem of a point phonon force *F* in the *x*_{3}-direction and the problem of a combination of point phonon and phason forces *Q*_{α} in the *x*_{1}-direction or *P*_{α} in the *x*_{2}-direction. As a classical elastic problem, the fundamental solutions of point phonon forces applied in the *x*_{3}-direction were derived based on the general solution (Gao & Zhao 2009).

### (a) Solution for a point force *F* in the *x*_{3}-direction

This is an axisymmetric problem discussed in Gao & Zhao (2009). The potential functions have the form as
3.1where
*A*_{i} are undetermined constants and sign() is the signum function.

Since the general solution expressed by the potential functions *ψ*_{i} and *ψ*_{j} satisfies the basic equations of QC elasticity, such as the deformation geometry equations, the equilibrium equations and the constitutive equations, the boundary-value problems of QC elasticity are transformed into finding a solution that satisfies the continuity conditions and the prescribed boundary conditions. Substitution of equations (3.1) into equations (2.4) and (2.6) yields the expressions for the phonon and phason fields below
3.2and
3.3

Two kinds of continuity conditions should be considered in the boundary-value problems to determine the free constants of the potential functions. One kind of continuity condition is the mathematical continuity which means that *ψ*_{i} and *ψ*_{j} are continuous functions. The other kind of continuity condition is the mechanical continuity which assumes that the displacements and stresses are continuous. This solution from equations (3.2) and (3.3) shows that the displacement component *u*_{3} and the stress components *T*_{α3}, *T*_{α4} and *σ*_{33} are continuous functions except at the origin. However, the continuity of the other displacement and stress components and the potential functions *ψ*_{i} in the plane *x*_{3}=0 is unclear and needs further examination. Since the components *U*_{αβ}, *T*_{αβ}, *σ*_{12}, *H*_{12}, *H*_{21} and *ψ*_{i} are odd functions of *x*_{3}, they must vanish across the plane *x*_{3}=0. This implies that
3.4In view of equations (3.4), the second and third equations of equations (3.4) are satisfied automatically. Besides, the equilibrium condition of a layer between any two planes, such as *x*_{3}=±*h*, must meet the requirement stipulated as follows:
3.5Substituting equations (3.3) into equation (3.5), one obtains
3.6Thus, the three unknown constants *A*_{i} can be calculated from the first equation of equations (3.4) and (3.6), written down matrix notation yielding
3.7aor in another form
3.7b

### (b) Solution for point forces *Q*_{α} in the *x*_{1}-direction or *P*_{α} in the *x*_{2}-direction

When the point forces *Q*_{α} are applied at the origin along the *x*_{1}-direction, it can be assumed that
3.8where *B*_{i} and *B*_{j} are constants to be determined. Following the same manipulation as in the previous case, the expressions for the phonon and phason fields are summarized as follows:
3.9
3.10

As before, the continuity of the components *u*_{3}, *T*_{α3} and *T*_{α4} on *x*_{3}=0 demands
3.11

Also, the equilibrium condition of the layer bounded by the planes *x*_{3}=±*h* becomes
3.12

Substituting the expression for *T*_{α4} in equations (3.10) into equation (3.12) yields
3.13From the integral result in equation (3.13), we know that the integral in equation (3.12) is independent of the value of *h*. The five unknown constants *B*_{i} and*B*_{j} can be determined from equations (3.11) and (3.13) as follows:
3.14where
3.15

When the problem of infinite QC solids subjected to point forces *P*_{α} in the *x*_{2}-direction is considered, we should replace *x*_{1} by *x*_{2} and *x*_{2} by −*x*_{1} in equations (3.8), respectively, and replace *Q*_{α} by *P*_{α} in equations (3.14). Therefore, if the problem of combination of point forces *Q*_{α} in the *x*_{1}-direction and point forces *P*_{α} in the *x*_{2}-direction is considered, the potential functions *ψ*_{i} and *ψ*_{j} take the form
3.16

By using the general solution in equations (2.4) and heuristic functions with simple forms, the expressions of the displacements and stresses can be obtained. It is important to note that these expressions are continuous except at the origin. The only approximation is introduced by the approximate specification of the boundary conditions of the layer, i.e. the boundary conditions in equations (3.5) and (3.12) are specified in terms of the stress resultants, instead of the stress distribution in the vicinity of the origin. Therefore, in the cases where Saint-Venant’s principle holds, Green’s function solutions should be very accurate ones.

### (c) The degenerated form of infinite quasi-crystal solids

Determination of the independent elastic constants *C*_{ij}, *K*_{i} and *R*_{i} for different kinds of QCs depends on their symmetries with the group representation theory (Ding *et al.* 1993; Hu *et al.* 1996, 2000). It is noted that, although *C*_{ij} in QCs can be measured by some experimental methods, *K*_{i} are difficult to measure (Tanaka *et al.* 1996). Significant progress in this area has been made by Jeong & Steinhardt (1993), who evaluated *K*_{i} of decagonal QCs by Monte Carlo simulation. The values of *K*_{i} are of the same order of magnitude as *C*_{ij} obtained by resonant ultrasound spectroscopy (Chernikov *et al.* 1998). There are no data available for *R*_{i} which, based on the estimation of some experts (Edagawa 2007; Takeuchi & Edagawa 2007) working in the field of QCs, hold lower values than *K*_{i}.

However, the relevant data such as *K*_{i} and*R*_{i} associated with the present paper are still lacking. Alternatively, we will discuss a degenerated form of infinite QC solids to investigate its validity, i.e. a two-dimensional QC body reduces to a transversely isotropic elastic body. In this case, it can be shown that *R*_{i}=0. Hence the governing equations (2.1)–(2.3) reduce to two groups of equations for uncoupled phonon and phason field problems, respectively. Then, the cubic algebraic equation *as*^{6}−*bs*^{4}+*cs*^{2}−*d*=0 can be reformulated as
3.17Let be the roots of the first multiplier of equation (3.17), and be the root of the second multiplier of equation (3.17) with no loss of generality. The quadratic algebraic equation is rewritten as
3.18We assume and . Then it can be seen that and relate only to elastic constants in the phonon field, while and associate only with elastic constants in the phason field.

For the transversely isotropic elastic body, the constants *k*_{3i}, *k*_{2i} and *k*_{2j} take the following form inserting from equation (3.17):
3.19and *k*_{23}, *k*_{25}≠0 which associate with and . Since the analysis in the following calculation does not involve *k*_{33}, *k*_{23}, *k*_{24} and *k*_{25} except the requirement *k*_{23}, *k*_{25}≠0, it suffices to discuss only *k*_{3α} and .

For the case of the point force *F* in the *x*_{3}-direction, the solution of equations (3.7*a*,*b*) reduces to
3.20Substitution of equation (3.20) into the degenerated equation of equations (3.2) leads to the expressions of displacements in transversely isotropic elastic solids
3.21From equations (3.2) and (3.7*a*,*b*), the expressions of displacements in two-dimensional QC solids are
3.22

In comparison with Green’s function of transversely isotropic solids, the configuration and deformation in two-dimensional QC solids are strongly related with not only the phonon elastic constants, but also phason and phonon–phason coupling elastic constants, even if the point phonon force is applied alone. Owing to the introduction of the phason field, a theoretical description of the deformed state of QCs requires a combined consideration of interrelated phonon and phason fields, so the elasticity of QCs is much more complex than that of conventional crystals. Green’s function solutions provide important information for studying the mechanical behaviours of the new solid phase and understanding clearly the interplay of the interaction between the phonon and phason activity.

For the case of the point forces *Q*_{1} in the *x*_{1}-direction, the solution of equations (3.14) becomes
3.23

The above results in equations (3.20) and (3.23) are the same as the corresponding results deduced by Ding *et al*. (1997*b*). To illustrate the applications of Green’s functions, two numerical examples, isotropic elastic solids with linear distributed loads on one surface and transversely isotropic elastic solids containing a spheroidal cavity under uniform tension along the *x*_{3}-direction, were examined by the boundary element method (Ding *et al.* 1997*b*). Comparing the results with the known solutions, these examples show that the exact or accurate solutions may be obtained by applying Green’s function. When the relevant material data are available, the same numerical examples of two-dimensional QCs can be also obtained by performing similar derivations.

As Ding *et al.* (1997*b*) mentioned, Green’s function solutions for transversely isotropic elastic solids are valid for two cases in which two eigenvalues are distinct or equal to each other. Thus, these are the unified solutions for both isotropic elasticity and transversely isotropic elasticity. This property means that there is no need to worry about the eigenvalue conditions. This is particularly important for materials with very close eigenvalues. By setting *R*_{i}=0, the coefficients in equations (3.7*a*,*b*) and (3.14) are available for both the two-dimensional QC and degenerated cases (transversely isotropic elastic solids), so Green’s function solutions are expressed in the united form for both two-dimensional QC elasticity and pure elasticity. Therefore, the same codes can be used to perform the numerical computation for different cases implementing an appropriate algorithm.

## 4. Three-dimensional Green’s functions for infinite QC solids composed oftwo half-spaces

An infinite solid composed of two two-dimensional hexagonal QC half-spaces with different material constants is shown in figure 2. Let the lower half-space *x*_{3}≥0 be occupied by material 1 and the upper half-space *x*_{3}≤0 be occupied by material 2. The interface of the two half-spaces is parallel to their quasi-periodic plane and the point forces are applied at the point (0, 0, *h*).

If the two half-spaces are rigidly bonded together along the interface *x*_{3}=0 such that the components of stresses and displacements are continuous across the interface, we have the following boundary conditions:
4.1where superscripts ()^{−} and ()^{+} denote the variables in the lower half-space *x*_{3}≥0 and the upper half-space *x*_{3}≤0, respectively. On the other hand, if the half-spaces are in smooth contact, i.e. in a complete contact with relative frictionless movement of both half-spaces along the *x*_{1}–*x*_{2} plane, then
4.2

In the following discussion, the components of displacements in the lower half-space can be decomposed into two parts:
4.3where *U*_{αβ} and *u*_{3} are the displacements owing to Green’s functions for infinite QC solids, and *U*′_{αβ} and *u*′_{3} are the displacements which make equations (4.3) satisfy the boundary conditions (4.1) or (4.2) on the interface. In the lower half-space, *U*_{αβ} and *u*_{3} can be obtained simply by replacing *x*_{3} by *x*_{3}−*h* in equations (3.2), (3.3), (3.9) and (3.10).

### (a) Point force *F* in the *x*_{3}-direction

When only a point phonon force *F* is applied at the origin along the *x*_{3}-axis, for the case of the lower half-space, namely, *x*_{3}≥0, we assume that
4.4where *m*=1, 2, 3, and
and
*A*_{i} have been obtained in equations (3.7*a*,*b*), and are nine constants to be determined. On substituting equations (4.4) into equations (2.4) and (2.6), we can obtain the corresponding components of displacements and stresses as follows:
4.5and
4.6

For the case of the upper half-space, i.e. *x*_{3}≤0, and take the form
4.7where
and are also nine constants to be determined. By using equations (2.4), (2.6) and (4.7), we can obtain the displacements and stresses in the upper half-space
4.8
4.9

The undetermined constants can be obtained using the contact conditions on the interface. When the two half-spaces are ideally bonded, the conditions in equations (4.1) lead to 4.10When the two half-spaces are in smooth contact as defined by equations (4.2), we get 4.11

It can be seen that there are totally 18 independent linear algebraic equations involved in equations (4.10) or (4.11), from which the 18 unknown constants and can be uniquely determined. However, the solutions will be lengthy and we do not list them here. When calculating the displacements and stresses numerically, it is easier to solve the equations listed in equations (4.10) or (4.11) directly. According to the interface conditions (4.1), only some components (*U*_{αβ}, *u*_{3}, *T*_{α3}, *T*_{α4} and *σ*_{33}) are continuous across the interface, while the others (*T*_{αβ}, *σ*_{12}, *H*_{12} and *H*_{21}) are discontinuous if the two bonded half-spaces are really dissimilar. This continuous property is in line with the physical understanding of two half-spaces with perfect interface. Therefore, more work has to be done if these Green’s function solutions are to be adopted to analyse the coupling fields at interfaces. One example is using Green’s function as fundamental solutions in integral kernels of the boundary element method.

### (b) Point forces *Q*_{α} in the *x*_{1}-direction or *P*_{α} in the *x*_{2}-direction

When point forces *Q*_{α} are applied at the origin along the *x*_{1}-direction, for material 1 in *x*_{3}≥0, we take and as
4.12where *n*=4, 5. and have the same structures as and . *B*_{i} and *B*_{j} have been obtained in equations (3.14), yet and are 13 undetermined constants. Substitution of equations (4.12) into equations (2.4) and (2.6) leads to
4.13
4.14For material 2 in *x*_{3}≤0, we assume that
4.15where and are 13 constants to be determined again. has the same structure as . With the use of equations (4.15), we rewrite equations (2.4) and (2.6) as
4.16
4.17

If the two half-spaces are ideally bonded, similarly, the conditions listed in equations (4.1) lead to 4.18Again, when the two half-spaces are in smooth contact on their interface, the conditions specified by equations (4.2) lead to 4.19

It can be seen again that there are totally 26 independent linear algebraic equations involved in equations (4.18) or (4.19), from which the 26 unknown constants , , and can be solved.

If the problem of the QC half-spaces subjected to point forces *P*_{α} in the *x*_{2}-direction is considered, in a similar manner to infinite solids we replace *x*_{1} by *x*_{2} and *x*_{2} by −*x*_{1} in equations (4.12) and (4.15), respectively, and replace *Q*_{α} by *P*_{α} in equations (4.18) and (4.19). Thus, the potential functions , , and for the problem of combination of point forces *Q*_{α} in the *x*_{1}-direction and point forces *P*_{α} in the *x*_{2}-direction are as follows:
4.20
4.21In equations (4.20) and (4.21), *l*_{αi} and *l*_{αj} have been determined in equations (3.15). Besides, the algebraic equations to determine , , and are the same as equations (4.18) or (4.19), only , , and in these equations are replaced by the corresponding constants , , and , respectively. When the QC materials of the two half-spaces are the same, that is to say, , , , , and in these equations are equal to zero, the results obtained above reduce directly to those in the previous section for infinite QC solids, provided that *h*=0.

## 5. Three-dimensional Green’s functions for half-space QC solids

For classical elastic problems, Lorentz (1907) developed the fundamental solutions for the displacement boundary-value problem of the elastic half-space. Phanthien (1983) showed that Lorentz’s results could be further extended to point forces applied in the interior of a half-space with a fixed boundary. By superposing a complementary part of the solution to Kelvin’s full-space function, Mindlin (1936) gave the fundamental solutions of the same problem under traction-free boundary conditions. These two solutions construct the basis of solving elastic problems in the half-space. According to the above results, we can directly extend Mindlin’s and Lorentz’s results for half-space QC solids with either free or fixed boundaries.

For the generalized Lorentz problem of half-space QC solids, the boundary condition at the surface is that the displacement components in the lower half-space are zero, i.e. . As a result, the right-hand sides of the first two equations of equations (4.10) vanish, and we can get nine algebraic equations to determine the constants . 5.1Similarly, the constants and can be determined from the first three equations of equations (4.18), so they can be conveniently solved as follows: 5.2

For the generalized Mindlin problem of half-space QC solids, the boundary condition at the surface of the lower half-space is traction-free, i.e. . As a result, the right-hand sides of the last two equations of equations (4.10) and the last three equations of equations (4.18) also vanish, so we can get the corresponding algebraic equations to determine the constants , and . It is useful to list them here:
5.3
5.4When the point forces are applied to the surface of the half-space, i.e. *h*=0, the generalized Mindlin problem degenerates to the generalized Boussinesq problem (Gurtin 1972).

Up to here, these three generalized problems of half-space QC solids can be solved directly from the fundamental solutions of two infinite half-spaces. Results show that the fundamental solutions developed in this paper are reliable and help to describe problems in an incisive way, so they can serve as a basis for further applications.

## 6. Conclusions

The present paper studies the problems of combination of point phonon forces and point phason forces applied in infinite spaces and bimaterials, which consist of two half-spaces of dissimilar two-dimensional QC materials bonded together. By using the general solution of two-dimensional hexagonal QC, a series of displacement functions is established, and the fundamental solutions of infinite spaces and two half-spaces with point forces applied in the interior of QC solids are obtained. Furthermore, the present solutions are reduced to the solutions of infinite transversely isotropic solids, the fundamental solutions of the generalized Lorentz problem, the generalized Mindlin problem and the generalized Boussinesq problem of half-space QC solids with either free or fixed boundary. The method for deriving the displacements and stresses presented here has some merits. One of the key features of the method is that the physical quantities can be readily calculated without the need of performing any transformation operations.

Three-dimensional Green’s functions can be extended to a wide range of engineering problems, for instance interlaminar stresses and surface responses in a laminated structure composed of different material layers (Gao & Noda 2005), interfacial cracks and contact in earthquake/rock engineering (Gao & Wang 2001), and inverse evaluation of materials properties using experimental approaches. Moreover, the solutions can be directly applied to the study of strained quantum dot semiconductor devices (Pan & Yang 2001), since a semiconductor device commonly appears as an anisotropic heterostructure, i.e. a layered structure.

The performance of QC materials is influenced by the presence of defects such as dislocations, inclusions, cracks, etc. Such solutions are very convenient to be used in the study of point defects and inhomogeneities in the materials. These provide important information for further studying the deformation and fracture of the new solid phase and understanding clearly the interplay of the interaction between the phonon and phason activity. They also play an important role in numerical simulations such as the finite element method and the boundary element method. The straightforward but tedious formulation process can be simplified if a computer symbolic manipulator software package, such as Mathematica or Maple, is used.

## Acknowledgements

The work is supported by the National Natural Science Foundation of China (no. 10702077) and the Alexander von Humboldt Foundation in Germany.

- Received January 10, 2011.
- Accepted March 8, 2011.

- This journal is © 2011 The Royal Society