Elastodynamic Analysis of a Hollow Cylinder with Decagonal Quasicrystal Properties: Meshless Implementation of Local Integral Equations

A meshless approximation and local integral equation (LIE) formulation are proposed for elastodynamic analysis of a hollow cylinder made of quasicrystal materials with decagonal quasicrystal properties. The cylinder is assumed to be under shock loading. Therefore, the general transient elastodynamic problem is considered for coupled phonon and phason displacements and stresses. The equations of motion in the theory of compatible elastodynamics of wave type for phonons and wave-telegraph type for phasons are employed and can be easily modified to the elasto-hydro dynamic equations for quasicrystals (QCs). The angular dependence of the tensor of phonon–phason coupling coefficients handicaps utilization of polar coordinates, when the governing equations would be given by partial differential equations with variable coefficients. Despite the symmetry of the geometrical shape, the local weak formulation and meshless approximation are developed in the Cartesian coordinate system. The response of the cylinder in terms of both phonon and phason stress fields is obtained and studied in detail.


Introduction
The observation of five-fold rotational symmetry by the electron diffraction patterns with bright diffraction spots in a rapid cooled Al-Mn alloy in April 1982 by Shechtman [1] yields surprising discoveries on the structure of matter and its properties.Since the five-fold rotational symmetry was in contradiction with the basic laws of the symmetry of crystals and the pattern also did not exhibit translational periodicity, this discovery finally resulted in redefinition of the concept of crystals, and the order is now associated not only with periodicity but also with quasi-periodicity.Soon after that discovery, five-fold rotational symmetry and three-dimensional icosahedral quasicrystals were found.Now, we also know about pentagonal, octagonal, decagonal and twelve-fold decagonal quasicrystals, which are two-dimensional quasicrystals with quasi-periodic atomic arrangement along two directions and periodic along the third direction, which is just the direction of 5-, 8-, 10-, and 12-fold rotational symmetry axes, respectively.The planar quasicrystals are different from the two-dimensional ones, since the former belong to a two-dimensional structure (without the third dimension) with a quasiperiodic atomic arrangement.In one-dimensional quasicrystals, the atomic arrangement is quasiperiodic along one axis, and periodic along the plane perpendicular to the axis.The icosahedral and decagonal quasicrystals present the majority of observed quasicrystals.It should be stressed that the unusual structure of quasicrystals leads to a series of properties different from those of crystals.The thermal as well as electrical conductivity of quasicrystals is lower than that of the conventional metals.Due to the lack of periodicity, the Bloch theorem and Brillouin zone concept are inapplicable and the electronic structure as well as electronic energy spectra of quasicrystals exhibit different properties than crystals.Moreover, no conservative slip along crystallographic planes can be allowed in quasicrystals, which causes quasicrystals to be hard and brittle.Quasicrystals are deformable under applied external sources, and the physical basis of the elasticity of quasicrystals is substantially different from the classical elasticity [2,3].The Landau theory is mostly used on elementary excitations and symmetry-breaking of condensed matter.The other approach is based on the unit cell description based on the Penrose tiling.
The numerical simulations play an important role in the investigation of phenomena in new materials as well as in design of structural elements.Some exact solutions have been presented for one-dimensional cases in beams and plates made of quasicrystals [4,5].Yang et al. [5] proposed an exact solution for analysis of plates subjected to surface loading and made of one-dimensional orthorhombic quasicrystal.They obtained the displacement and stress fields in a sandwich plate with simply supported boundary conditions, which was assumed to be made of quasicrystals and crystals with two different stacking sequences.
One of the most important analyses in engineering applications is elastodynamic analysis of structures.The structures made of quasicrystals (QCs) may be under dynamic and transient loadings.Thus, the general dynamic theory for QCs should be used to obtain the realistic dynamic behaviors of QCs.There are some theories for the description of dynamic behaviors of QCs, which were presented in previously published papers.Recently, the importance of them are addressed in a research paper by Agiasofitou et al. [6], of which the most important of them are the elastodynamic properties of wave types for QCs presented by Fan et al. [2], Shi [7] and Fan [3].A minimal model of the phonon-phason elastodynamic, which is a combination of models based on wave type and hydrodynamics of QCs, was proposed by Rochal et al. [8].This model was adopted by many researchers as the elasto-hydrodynamics of QCs [9].In another work, the icosahedral quasicrystals with inclusion of the pinning effect is dynamically analyzed using a minimal model of the phonon-phason dynamics considering resonant attenuation of low-frequency acoustic waves in the temperature range corresponding to thermal activation of phasons by Kozinkina et al. [10].Recently, an elastodynamic model of the wave-telegraph type was proposed for the description of dynamic behaviors of QCs by Agiasofitou et al. [6].
There are some computational methods for solving the governing equations based on various elastodynamic models of QCs.For example, recently, the meshless local Petrov-Galerkin (MLPG) was successfully used by Sladek et al. [11] for bending analysis in 1D orthorhombic QCs under static and dynamic loadings.They used the Bak and elasto-hydrodynamic models for phason governing equation in the elastodynamic case of their problem.In addition, the MLPG method was employed for analysis of initial-boundary value crack problems in decagonal QCs by Sladek et al. [12].There are some analytical techniques that can be applied to elastodynamic analysis of QCs, and were successfully used in other elastodynamic or thermoelasticity analyses.Hosseini et al. [13] proposed an analytical method based on series expansions for solutions of elastodynamic problems in functionally graded thick hollow cylinders under shock loading.In another work, they used the proposed analytical method for coupled thermoelasticity analysis based on Green-Naghdi theory (without energy dissipation) [14].Wang et al. [15] presented analytical solutions for some defect problems in 1D hexagonal and 2D octagonal quasicrystals.
A fundamental solution was proposed for the time-dependent differential equations of anisotropic elasticity in 3D quasicrystals using the Fourier transform with respect to the space variables and matrix transformations by Cerdik Yaslan [16].
In this paper, we present an effective numerical method for dynamic analysis of cylinder made of 2D decagonal QCs while assuming translational symmetry along the axis of the cylinder.The conventional wave propagation equation is used for description of phonon waves while the dynamics of phasons is described by the wave-telegraph equation.Recall that the wave-telegraph equation can be easily modified to a diffusion type equation used in the elasto-hydrodynamic model of QCs.The governing equations of motions and constitutive equations for QCs cylinder are formulated and the utilization of the polar coordinate system is discussed.Then, the equations of motions in terms of momenta and stresses are considered in a local weak sense and the Local Integral Equations (LIE) are derived.The meshless approach is used for approximation of both the phonon and phason displacements, and the stress tensors are approximated by using the constitutive relationships and derivatives of approximated displacements.The dynamic response of the QC cylinder is studied in detail while presenting the comparison of results by using various formulations for QCs dynamics.Furthermore, the influence of the quasi periodicity (quasi periodic arrangement of atoms in plane orthogonal to the cylinder axis) on the angular dependence of the response of the QC hollow cylinder on the angularly symmetric loading is also studied.

Mathematical Formulations
Similar to other aperiodic crystals, quasi-periodicity induces new degrees of freedom.The Landau phenomenological theory on the symmetry-breaking of condensed matter [17] played the central role in the development of elasticity theory of quasicrystals and has been widely acknowledged.In many physical systems (classical or quantum), the motion presents the discrete spectrum (corresponding to an eigenvalue problem of a certain operator from the mathematical point of view), the lowest energy level state is called ground state, and that beyond the ground state is called the excited state.According to the Debye and Born theory of crystal lattices, the propagation of the lattice vibration under the long-wavelength approximation can be seen as continuum elastic waves.The quanta of lattice vibration (the smallest portion of energy of the elastic wave) is phonon (elastic waves belong to acoustic waves).Unlike the photon, the phonon is not an elementary particle, but, in the sense of quantization, the phonon exhibits natural similarity to elementary particles (phonons are quasi-particles emerging from the quantization of elastic waves in linearized elastic setting).It is known that standard liquid exhibits arbitrary translational as well as arbitrary rotational symmetry, and the periodicity of crystals breaks these symmetries of liquid.According to the Landau theory, the phonon is elementary excitation resulting from symmetry-breaking in the transition from liquid phase to crystal.In transition from crystals to quasicrystals, the break of periodic symmetry of crystals results in new elementary excitation called phasons.Thus, quasi-periodicity (similar to other aperiodic crystals) induces new degrees of freedom (phason displacements), but the phason modes in quasicrystals present different nature, i.e., the motion of atoms exhibits discontinuous jump rather than the long-wavelength propagation.The phason displacements describe relative atomic shifts, which assure quasi-periodicity of the atomic distribution.Hence, the phason field is insensitive to translations of coordinate frames in the physical space.Since we shall deal with quasicrystals within the phenomenological theory of continuum mechanics, the adjectives phonon and phason will be utilized just to distinguish between the displacements known in classical elasticity (phonon displacements) and new degrees of freedom induced by the quasi-periodicity in quasicrystalline materials (phason displacements).
Within linear elasticity, elastic free energy of quasicrystals consists of three types of quadratic terms: phonon-phonon, phason-phason and phonon-phason coupling terms [3].Then, the constitutive laws are given as: where σ ij and H ij are the phonon and phason stresses, respectively, while the corresponding phonon strains ε kl and phason strains ω kl are defined in terms of gradients of the phonon displacements u i (x, t) and phason fields w i (x, t) as: Recall that ω kl is not symmetric, while ε kl is a symmetric tensor (the rigid body rotations have no contribution to deformation energy).The tensors of material coefficients and coupling coefficients are given as [3]: where c 11 , c 12 , K 1 , K 2 , and R are five constants characterizing the material properties of quasi-crystals, and ε 3lk is the permutation symbol with ε 123 = 1.Remember the following symmetry relationships, which are satisfied by the tensors of material coefficients: In problems for a hollow cylinder with translational symmetry along the cylinder axis, the plane elasticity conditions are applicable and polar coordinates are appropriate for modelling of geometry.Performing the transformation from the Cartesian to polar components of tensors of material coefficients, we get the constitutive laws in polar coordinate system It is seen that the tensor of coupling coefficients is position dependent in a plane perpendicular to the periodic symmetrical axis of quasicrystals.Thus, the angular symmetry is inapplicable in such a quasicrystal cylinder even if the loading and boundary conditions were angularly symmetric.Moreover, the expressions for the phonon as well as phason strains and stresses are more harmful in the polar coordinate system than in the Cartesian coordinate system, and the governing equations in polar coordinates would be given by partial differential equations with variable coefficients.Therefore, the utilization of the polar coordinate system seems to be inappropriate in the considered problem despite the angular symmetry of the geometrical shape.Nevertheless, the constitutive laws are symmetric with respect to rotations on the angle equal to π as can be seen from the representations of the tensors of material coefficients in Equations ( 5) and ( 6) in the polar coordinate system.
The equations of motion in the theory of compatible elastodynamics of wave-telegraph type of quasicrystals in terms of the momenta and stresses take the form [4] .
where f i , g i are the conventional (phonon) body force density and a generalized (phason) body force density, respectively, while g w j is the frictional force (dissipative force) vector.Recall the principal difference between f i and g i .There are no bulk external fields acting on phason degrees of freedom, except we presume phason inertia [18][19][20].However, it is proven that bulk phason self-action which has both the dissipative (frictional force) and conservative components with the latter being dependent on the phason field [20].The phason self-action is inside the material elements in contrast to external actions.In the isotropic case, D ij = Dδ ij and D > 0 is the phason friction coefficient.Introduction of the constitutive laws into Equations ( 7) and ( 8) yields the equations of motion for phonon and phason displacements.Both of these equations of motion are hyperbolic partial differential equation (PDE) with the phonon modes being described by a wave type equation while the phason modes by a telegraph type PDE.The latter equation leads to waves damped in time and finite propagation velocity.Since these two equations of motion are coupled, the phonon dynamics are affected by phason dynamics and vice versa.Recall that this model also includes two particular models developed before the Agiasofitou and Lazar models occur [6]: (i) the elastodynamic model of wave type ( ρ = ρ, D = 0) could describe phason waves with finite velocity but without attenuation; (ii) the hydrodynamic and/or the elasto-hydrodynamic model ( ρ = 0, D = 0) support a diffusion type equation for the dynamics of the phason field.However, as is well known, the diffusion equation supports exclusively an infinite velocity of propagation.The physically inadmissible velocity of information propagation can be overcome by replacing the parabolic PDE by the hyperbolic one.
In order to model waves damped exponentially in time and propagating with finite velocity as observed experimentally [21,22], Agiasofitou and Lazar [6] introduced a telegraph type equation for phason displacements.As shown above, utilization of the wave-telegraph equation enables us to step towards both the wave type and diffusion type description of phason dynamics by proper choice of the coefficients ρ and D.
Since initial-boundary value problems for the considered system of PDE are rather complicated for finding analytical solutions, we propose combining the Laplace transformation technique (with respect to the time variable) with a spatial discretization and approximation.Moreover, the application of the discretization to governing equations for displacements seems to be awkward, and we apply it directly to Equations ( 7) and ( 8) in a weak sense.

Weak Formulation of Local Integral Equations (LIEs)
Eliminating the time variable in governing Equations ( 7) and ( 8) by using the Laplace transformation, the governing equations for Laplace transforms become: If we assumed homogeneous initial conditions for both the phonon and phason displacement fields and vanishing body forces, the right-hand sides in Equations ( 9) and (10) would be equal to zero.
Let us consider the set of coupled governing Equations ( 9) and (10) in a weak sense on local subdomains considered around interior points x c ∈ Ω c ⊂ Ω. Taking the weight function as the Heaviside function supported on Ω c , we obtain the system of local integral equations (LIE): in which we have utilized the Gauss divergence theorem and ∂Ω c is the boundary of Ω c .The derived local integral equations (LIEs) (11) and ( 12) are written in general vector form, and one can represent the vectors in an arbitrary coordinate system.The Cartesian coordinate system appears to be the most appropriate since the basis vectors are independent on the integration variables.Then, the component representation is obtained directly from Equations ( 11) and ( 12) as: with i, j ∈ {1, 2}.
In order to facilitate the integrations, it is appropriate to choose the local subdomain as a circle around the point Then, the integration over ∂Ω c can be performed as the integration with respect to the angular variable.Now, in view of Equations ( 1)-( 4), the LIE can be rewritten as: with τ i being the unit tangent vector related to the unit outward normal vector as τ i = ε 3ji n j .
Thus, what we need is to utilize certain approximation for the phonon as well as phason displacements and their gradients.

Meshless Approximation for Phonon and Phason Displacements
To the spatial variations of phonon and phason displacements u k (x, s), w k (x, s), the meshless approximation over a number of randomly located nodes x c , c = 1, 2, ..., n, is applied using the radial basis functions (RBF).Thus, the considered approximations take the form: where R c (x) is the radial basis functions associated with the node x c , and U c k (s), W c k (s) are expansion coefficients for phonon and phason displacements, respectively.The radial basis functions employed in this article are multiquadrics The shape parameter c is proportional to the minimum distance between nodes, while m is chosen to be m = −2.Taking the approximations as Equations ( 17) and (18) at nodal points, one obtains the following systems of equations: which can be rewritten in the matrix notation as Hence, the expansion coefficients can be expressed in terms of nodal unknowns as: Substituting Equation (23) into Equations ( 17) and ( 18), we receive where is the shape function associated with node x c and u c k (s), w c k (s) are the nodal values of the Laplace transforms of phonon and phason displacements, respectively.
Considering the problem without body forces (phonon body forces and conservative bulk phason self-action) and with homogeneous initial conditions, the discretized local integral equations (LIEs) (15) and ( 16) collocated at an interior node x c ∈ Ω c ⊂ Ω can be written as: The global boundary of the analyzed domain is decomposed into two parts Γ = Γ e ∪ Γ n , where Γ e and Γ n are parts on which essential and natural boundary conditions are prescribed, respectively.
For x b ∈ Γ e , we collocate the meshless approximation of phonon and phason displacements, i.e., For x b ∈ Γ n , we consider the discretized LIEs ( 15) and ( 16) collocate at x b : where t i (x, s) and q i (x, s) are the Laplace transforms of prescribed values of phonon and phason tractions defined as t i (x, t) = σ ij (x, t)n j (x), q i (x, t) = H ij (x, t)n j (x), respectively.Equations ( 27)-(32) give us the complete set of algebraic equations for computation of unknown nodal values of the Laplace transforms of phonon and phason displacements.Having solved this system for unknowns, one can obtain the time-dependent solutions by applying the inverse Laplace-transformation to the transformed quantities such as the Talbot technique.
Recall that the proposed formulation and calculations are performed in the Cartesian representation of the phonon and phason stresses and displacements.The representation in the polar coordinate system can be obtained by using the Mohr's circle technique.

Numerical Example and Discussions:
To show some numerical results of the problem, an infinite length sector of a hollow cylinder is assumed with the following mechanical properties and inner radius r a = 1 m, outer radius r b = 1 m and θ min ≤ θ ≤ θ max , θ min = 0, θ max = π for numerical analysis (see Figure 1).Firstly, to verify the presented method and data, the material of cylinder is assumed to be isotropic homogenous with the following mechanical properties: The considered boundary conditions for phonon modes are shown partially in Figure 1 and vanishing initial conditions are assumed.Since the phonon-phason coupling is absent, there is no phason response in the considered example and the dynamic behavior of cylinder is described as classical elastodynamic.The obtained results by the MLPG method for phonon radial stress along π Firstly, to verify the presented method and data, the material of cylinder is assumed to be isotropic homogenous with the following mechanical properties: The considered boundary conditions for phonon modes are shown partially in Figure 1 and vanishing initial conditions are assumed.Since the phonon-phason coupling is absent, there is no phason response in the considered example and the dynamic behavior of cylinder is described as classical elastodynamic.The obtained results by the MLPG method for phonon radial stress along the radial direction, and θ = π 2 at several time instants are compared with those obtained by finite element method (FEM), which are illustrated in Figure 2. The comparison shows a good agreement between the results.isotropic homogenous with the following mechanical properties: The considered boundary conditions for phonon modes are shown partially in Figure 1 and vanishing initial conditions are assumed.Since the phonon-phason coupling is absent, there is no phason response in the considered example and the dynamic behavior of cylinder is described as classical elastodynamic.The obtained results by the MLPG method for phonon radial stress along the radial direction, and 2 π θ = at several time instants are compared with those obtained by finite element method (FEM), which are illustrated in Figure 2. The comparison shows a good agreement between the results.The following mechanical properties correspond to Al-Ni-Co quasicrystals [11,12] and are considered for study on the transient response of QCs cylinder to Heaviside impact loading: The following mechanical properties correspond to Al-Ni-Co quasicrystals [11,12] and are considered for study on the transient response of QCs cylinder to Heaviside impact loading: Since the phason friction coefficient is assumed to be zero, the employed compatible elastodynamics of wave-telegraph type of quasicrystals degenerates to the elastodynamic model of wave type.The boundary of the analyzed domain is composed of four parts ∂Ω = ∂Ω a ∪ ∂Ω b ∪ ∂Ω m ∪ ∂Ω M as shown in Figure 1.In the shape of Cartesian components, the following boundary conditions are considered to study the dynamic behavior of both phonon and phason stress fields (see Figure 1): where H (t) is a Heaviside unit step function, n i (x) is the unit outward normal vector on the boundary, τ i (x) is the unit tangent vector on the boundary, and σ 0 = −1 Pa.The response to the shock loading by normal stresses applied on the inner cylindrical surface of the considered sector of the hollow cylinder is illustrated via the spatial distributions at several time instants as well as via the time evolution of phonon and phason fields at certain points of the analyzed domain.
Figure 3 shows the distribution of phonon radial stresses σ rr along the ray (r, θ = π/2) at several time instants, with the time step being assumed to be ∆t = 0.86 × 10 −5 s.The propagation of phonon radial stress wave fronts can be tracked in Figure 3.Likewise, the wave front of phonon hoop stress σ θ θ is propagated along the ray (r, θ = π/2) with a finite speed, as can be seen in Figure 4.
instants as well as via the time evolution of phonon and phason fields at certain points of the analyzed domain.
Figure 3 shows the distribution of phonon radial stresses r r σ along the ray ( , / 2) r θ = π at several time instants, with the time step being assumed to be 5 0.86 10 s t − Δ = × . The propagation of phonon radial stress wave fronts can be tracked in Figure 3.Likewise, the wave front of phonon hoop stress θθ σ is propagated along the ray ( , / 2) r θ = π with a speed, as can be seen in Figure 4.   Owing to the coupling, the phason stress field is affected by the considered excitation in phonon stress field.Figure 5 depicts the wave propagation of the phason radial stress field r r H along ray ( , / 6) r θ = π with a finite speed.A similar behavior can be observed also for the phason hoop stress field H θθ , as shown in Figure 6.Figures 5 and 6 demonstrate that the presented MLPG method has the capability to simulate the dynamic behavior of both phonon and phason fields and the interaction among them.Owing to the coupling, the phason stress field is affected by the considered excitation in phonon stress field.Figure 5 depicts the wave propagation of the phason radial stress field H r r along ray (r, θ = π/6) with a finite speed.A similar behavior can be observed also for the phason hoop stress field H θθ , as shown in Figure 6.Figures 5 and 6 demonstrate that the presented MLPG method has the capability to simulate the dynamic behavior of both phonon and phason fields and the interaction among them.
phonon stress field.Figure 5 depicts the wave propagation of the phason radial stress field r r H along ray ( , / 6) r θ = π with a finite speed.A similar behavior can be observed also for the phason hoop stress field H θθ , as shown in Figure 6.Figures 5 and 6 demonstrate that the presented MLPG method has the capability to simulate the dynamic behavior of both phonon and phason fields and the interaction among them.The problem is not angularly symmetric because of the phonon-phason coupling.Since the classical (phonon) stress loading is applied on the inner surface of a half of the infinite hollow cylinder, and the coupling coefficients are 100 times smaller than the classical elasticity coefficients, the response of phason fields is 100 times smaller than the response of phonon fields.Therefore the response of phonon fields is almost independent of the angular variable (see Figure 7), though such a dependence is evident in the response of the phason fields (see Figure 8).The problem is not angularly symmetric because of the phonon-phason coupling.Since the classical (phonon) stress loading is applied on the inner surface of a half of the infinite hollow cylinder, and the coupling coefficients are 100 times smaller than the classical elasticity coefficients, the response of phason fields is 100 times smaller than the response of phonon fields.Therefore the response of phonon fields is almost independent of the angular variable (see Figure 7), though such a dependence is evident in the response of the phason fields (see Figure 8).classical (phonon) stress loading is applied on the inner surface of a half of the infinite hollow cylinder, and the coupling coefficients are 100 times smaller than the classical elasticity coefficients, the response of phason fields is 100 times smaller than the response of phonon fields.Therefore the response of phonon fields is almost independent of the angular variable (see Figure 7), though such a dependence is evident in the response of the phason fields (see Figure 8).The distributions of both phonon and phason radial displacements, r u and r w , along the ray ( , / 2) r θ = π can be observed in Figures 9 and 10.The wave fronts of both the phonon and phason radial displacements can be tracked at various time instants because of a finite speed of wave propagation.Again, the coupling between the phonon and phason fields is confirmed.The distributions of both phonon and phason radial displacements, u r and w r , along the ray (r, θ = π/2) can be observed in Figures 9 and 10.The wave fronts of both the phonon and phason radial displacements can be tracked at various time instants because of a finite speed of wave propagation.Again, the coupling between the phonon and phason fields is confirmed.The distributions of both phonon and phason radial displacements, r u and r w , along the ray ( , / 2) r θ = π can be observed in Figures 9 and 10.The wave fronts of both the phonon and phason radial displacements can be tracked at various time instants because of a finite speed of wave propagation.Again, the coupling between the phonon and phason fields is confirmed.The effects of the value of R on the results for the phason fields can be observed in Figures 11   and 12 in which the distributions of phason radial and hoop stresses along the radial direction are drawn at 2 π θ = in a half of the hollow cylinder for various values of R , respectively.As it can be seen in Figures 11 and 12, the values of phason stress fields are increased by increasing the value of R .The effects of the value of R on the results for the phason fields can be observed in Figures 11  and 12 in which the distributions of phason radial and hoop stresses along the radial direction are drawn at θ = π 2 in a half of the hollow cylinder for various values of R, respectively.As it can be seen in Figures 11 and 12, the values of phason stress fields are increased by increasing the value of R.
Crystals 2016, 6, 94 14 of 16 and 12 in which the distributions of phason radial and hoop stresses along the radial direction are drawn at 2 π θ = in a half of the hollow cylinder for various values of R , respectively.As it can be seen in Figures 11 and 12, the values of phason stress fields are increased by increasing the value of R .

Conclusions
The elastodynamic model of a wave-telegraph type is used to study the dynamic response of decagonal quasicrystals.It is shown that angular symmetry is inapplicable in the case of a hollow cylinder even if symmetric boundary conditions are assumed.The tensor of material coefficients obeying coupling between the phonon and phason fields is dependent on the angular coordinate in polar coordinate system.Since the governing equations in such a coordinate system would be given by the PDE with variable coefficients, the utilization of the polar coordinates appears to be inconvenient.Therefore, the formulation is presented in the Cartesian coordinate system.The local Petrov-Galerkin (MLPG) method together with the Laplace transform technique with respect to time variables are employed for numerical solutions of boundary value problems.For the spatial variation of field variables, we have used the meshless approximation based on the radial basis functions, and the Talbot technique is employed for numerical inversion of the Laplace transformation.In a numerical example, a finite sector of an infinitely long thick hollow cylinder is considered under a shock loading with holding angularly symmetric boundary conditions.
Attention is paid to a detailed explanation of coupling between the phonon and phason fieldswave propagation of both the phonon and phason stress and displacement waves with a finite speed.Since the phason response to force loading is approximately 1000 times smaller than the phonon response, the latter is affected negligibly by phason fields, despite the fact that phononphason coupling is considered.It has been shown that the proposed method is suitable for numerical

Conclusions
The elastodynamic model of a wave-telegraph type is used to study the dynamic response of decagonal quasicrystals.It is shown that angular symmetry is inapplicable in the case of a hollow cylinder even if symmetric boundary conditions are assumed.The tensor of material coefficients obeying coupling between the phonon and phason fields is dependent on the angular coordinate in polar coordinate system.Since the governing equations in such a coordinate system would be given by the PDE with variable coefficients, the utilization of the polar coordinates appears to be inconvenient.Therefore, the formulation is presented in the Cartesian coordinate system.The local Petrov-Galerkin (MLPG) method together with the Laplace transform technique with respect to time variables are employed for numerical solutions of boundary value problems.For the spatial variation of field variables, we have used the meshless approximation based on the radial basis functions, and the Talbot technique is employed for numerical inversion of the Laplace transformation.In a numerical example, a finite sector of an infinitely long thick hollow cylinder is considered under a shock loading with holding angularly symmetric boundary conditions.
Attention is paid to a detailed explanation of coupling between the phonon and phason fieldswave propagation of both the phonon and phason stress and displacement waves with a finite speed.Since the phason response to force loading is approximately 1000 times smaller than the phonon response, the latter is affected negligibly by phason fields, despite the fact that phonon-phason coupling is considered.It has been shown that the proposed method is suitable for numerical simulations in rather complex problems.There are no restrictions on the geometry of the analyzed domain as well as on applied loads, and arbitrary material coefficients are allowed within the considered general elastodynamic model of the wave-telegraph type in decagonal quasicrystals.The presented methodology can be easily extended to other kinds of quasicrystals.Mariano and Planas [20] highlighted the role of the conservative bulk phason self-action in mechanical description of quasicrystals.In the presented analyses, the conservative bulk self-action is omitted.The next study will concern verification of such a role by numerical simulations.
where the tilde denotes the prescribed values by boundary conditions.If the collocation point lies on the boundary, x b ∈ Γ n , we define local subdomain as Ω b = Ω b ∩ Ω, where Ω b is the circular domain centered at x b .Then, the boundary of local sub-domain consists of two parts ∂Ω b = Γ b ∪ L b , where Crystals 2016, 6, 94 8 of 16

Figure 1 .
Figure 1.Illustration of the analyzed domain and boundary conditions.

Figure 1 .
Figure 1.Illustration of the analyzed domain and boundary conditions.

Figure 2 .
Figure 2. Comparison of phonon radial stresses r r σ distributions along the radial direction at

Figure 2 .
Figure 2. Comparison of phonon radial stresses σ r r distributions along the radial direction at θ = π 2

Figure 3 .
Figure 3. Distribution of phonon radial stresses r r σ in QC cylinder along the radial direction for

Figure 4 .
Figure 4. Distribution of phonon hoop stresses θθ σ in QC cylinder along the radial direction for

Figure 5 .
Figure 5. Distribution of the phason radial stresses

Figure 6 .
Figure 6.Distribution of the phason hoop stresses H θ θ in QC cylinder along the radial direction for θ = π 2 .

Figure 7 .
Figure 7. Distribution of phonon radial stresses r r σ in QC cylinder along the radial direction for

Figure 7 . 16 Figure 8 .
Figure 7. Distribution of phonon radial stresses σ r r in QC cylinder along the radial direction for various values of θ.Crystals 2016, 6, 94 13 of 16

2 Figure 8 .
Figure 8. Distribution of the phason radial stresses H r r in QC cylinder along the radial direction for various values of θ.

Figure 9 .
Figure 9. Distribution of phonon radial displacements r u in QC cylinder along the radial direction

ΔFigure 10 .
Figure 10.Distribution of phason radial displacements w r in QC cylinder along the radial direction for θ = π 2 .

Figure 11 . 16 Figure 12 .
Figure 11.Distribution of phason radial stresses H r r in QC cylinder along the radial direction for various values of R. Crystals 2016, 6, 94 15 of 16

Figure 12 .
Figure 12.Distribution of phason hoop stress H θθ in QC cylinder along the radial direction for various values of R.