Boundary Value Problems of Potential Functions in Decagonal Quasicrystals

Li, Wu; Xin, Hao; Fan, Tianyou

doi:https://doi.org/10.1155/2014/863047

Mathematical Problems in Engineering

On this page

Abstract Introduction Conclusion and Discussion Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2014 | Article ID 863047 | https://doi.org/10.1155/2014/863047

Boundary Value Problems of Potential Functions in Decagonal Quasicrystals

Wu Li,¹Hao Xin,²and Tianyou Fan³

Academic Editor: Evangelos J. Sapountzakis

Received16 Apr 2014

Revised29 Aug 2014

Accepted30 Aug 2014

Published14 Oct 2014

Abstract

A unified form of potential functions in decagonal quasicrystals (QCs) and conformal mappings are applied in a novel way to solve the boundary value problems emanating from the generalized theory of elasticity for decagonal QCs. By executing the reduction of boundary value problem to function equations, two crack problems are investigated. In the first one, an approximate analysis for bending specimen with a crack is given. In the other, a finite width strip with single edge crack of decagonal QCs is analytically estimated. Using the basic idea underlying Dugdale’s crack model, the extent of cohesive force zone in each of the two cases is analytically derived.

1. Introduction

Soon after the discovery of QCs [1], the boundary value problems of elasticity and defects of QCs have been put forward. In 2011, some scholars found 12-fold and 18-fold diffraction symmetric QCs in colloid, which may arouse great interest [2]. The unusual characters of the QCs are originated from their special atomic constitution. Due to the quasiperiodicity, QCs have a special type of elastic degrees of freedom, termed as phason degrees of freedom apart from phonon degrees which exist too for crystals, which have been studied extensively and widely [3–13]. On the basis of the physical framework and by extending the methodology of mathematical physics and classical elasticity, the mathematical theory of elasticity of QCs has been developed [11, 12, 14, 15]. Recent studies on the elasto-/hydrodynamics of QCs have made preliminary but significant progress [8, 10, 13, 16–18]. Wang et al. provide some useful insight into the interface problems of QCs [19]. Feuerbacher et al. have studied plastic deformation of QCs in microscopic mechanism, [20–29]. Li et al. have studied thermoelastic general solutions of one-dimensional hexagonal quasicrystals [30, 31]. More recently, Radi and Mariano investigated the straight cracks and dislocations in two-dimensional QCs by using the Stroh formalism and obtained the closed-form solution to the balance equations in terms of phonon and phason fields [32, 33]. The Stroh formalism for two-dimensional deformations of the icosahedral QCs was then studied by Li and Liu [34]. Guo et al. have studied a finite crack in a one-dimensional hexagonal quasicrystal strip [35]. As a material, QC is deformable under applied forces, thermal loads, and certain internal effects. The theoretical physicists have proposed various descriptions of its elasticity. The majority agrees that the Landau density wave theory is the physical basis of elasticity of QCs. Essentially, the description suggested that there are two displacement fields and in a QC: the former is similar to that in crystals, named as the phonon field; the latter is a displacement field, named as the phase field. The total displacement field in a QC is expressed by where denotes the direct sum; is in the physical space or the parallel space ; is in the complement space or the perpendicular space , which is an internal space. Under given loading conditions, they are giving rise to two elastic strain tensors and , respectively [36]: In the case of small deformation of the decagonal QCs, the phonon and phason stresses obey generalized Hooke’s law, namely, where , are phonon and phason displacements, and phonon stresses and strains, and phason stresses and strains, , , and the phonon, phason, and phonon-phason coupling elastic constants, respectively. Now we give the boundary and initial conditions; according to practical case the boundary consists of two parts and ; that is, ; at the tractions are given and at the displacements are prescribed. For the former case, where represents the outward normal unit vector at any point of the boundary , and represent the traction and generalized traction vectors, which are given functions at the boundary . And for the latter case, where and are known functions at the boundary .

In over 200 individual QCs observed to date there is over 1/3 belonging to two-dimensional decagonal QCs. So this kind of QCs is very important from fundamental point and applications. The present work is on a bending specimen with a crack and a finite width strip with a single edge crack, each of which is embedded in a two-dimensional decagonal quasicrystalline media. The complex potential theory method initially proposed by Muskhelishvili [37] is generalized to elliptic notch problem of the bending specimen. This promoted method postulated that it is valid only for rational conformal mappings. For finite width strip with a single edge crack, which belongs to an ideal plastic deformation problem, a transcendental function conformal mapping is accomplished to transform the region at the physical plane onto the upper-half plane of the mapping plane. In the first one, the complex potential functions are obtained, which can express the stresses and displacements. The stress intensity factors and the size of cohesive force zone can be derived analytically in these two problems. The practice shows that the procedure is powerful for complicated boundary value problem of the equations. Under these treatments, the original problems can be transformed into boundary value problems.

2. Mathematical Physics Theory of Elasticity of Decagonal QCs

First some preliminary results given using the complex potential theory will be recalled for a crack in a decagonal quasicrystalline medium.

Assume that it is perpendicular to the periodic symmetry axis (e.g., axis ); they are giving rise to two displacement fields for phonon and for phason in the deformation of the decagonal QCs. Further we consider here only for plane problem; that is, it is assumed that this leads to . In addition, in the absence of body forces, the equilibrium equations are as follows: When the phason field vanishes, these balance equations can be reduced to the balance equations of classical elasticity [11–15]. The explanation for readers may be seen in [38]. Equations (2), (3), and (7) are the basic equations describing elasticity of decagonal QCs under plane deformation; this is a set of 18 field equations. By using the displacement potential or stress potentials (see Fan [36]), the equation set can be reduced to a unique equation of higher order; that is, the governing equation is quadruple harmonic equation for two-dimensional decagonal QCs: The general solution of the quadruple harmonic equation can be written in terms of four analytic functions of a single complex variable : where the bar denotes the complex conjugate hereinafter; that is, . These analytic functions will be determined by boundary value conditions of practical problems. It is easy to prove that has no contribution to stress and displacement fields, so we can take that The phonon stresses , phason stresses , and the phonon displacements and phason displacements of QCs can be expressed as follows [36]:in which the constants are defined by

Next we define new symbols for convenience [36]:in which the prime, two primes, three primes, and superscript () denote the first to fourth order differentiation to variable .

In the following we consider only the stress boundary value problem; that is, the specific form of the boundary conditions at the boundary curve for decagonal QCs can be expressed by [36]Equation (14b) can be considered as the boundary conditions at the boundary curve in phason field, where and are the generalized stress in phason field. The right-hand side term of (14a) and (14b) is the resultant force of the boundary in phonon and phason fields. From (14a) and after some derivation the phonon stress boundary condition can be reduced to the equivalent form: Meanwhile, we have the equivalent form of phason stress boundary condition: Boundary conditions include both opening and shear modes (with respect to a crack); the stress intensity factors have a universal form; that is, Here and are the “stress intensity factors,” is the stress potential function, and is the conformal mapping corresponding to loading conditions, and is located in the tip of the crack. We choose to focus on the case of uniform remote traction-free crack boundaries. This choice, although the most general, is of great interest and will serve to obtain results. In linear elasticity, many other solutions may be obtained by superposition. Thus the boundary conditions at infinity, for the two in-plane symmetry modes of fracture, are presented as In addition, the free boundary conditions on the crack are expressed as where is the arc-length parametrization of the crack boundary and denotes the outward normal direction. In order to compute the full stress field one should first formulate the boundary conditions in terms of the analytic functions and remove the gauge freedom in (13a)–(13e). Indeed, we do not have enough boundary conditions to determine the stress potential uniquely. So we can allow in (13a)–(13e) arbitrary transformations of the form where is a real constant and , , and are complex constants. This choice provides seven degrees of freedom in the definition of the stress potential. Four of them are removed according to single-valued invariance of the displacement and stress and the free boundary conditions on the crack. It is important to stress that whatever the choice of the function the stresses of phonon field are unaffected, so it is not discussed. Meanwhile, the loadings on QCs considered are counterweight system. To proceed we represent and in Laurent expansion form: in which , , , , and are principal stresses at infinity, and is the angle of and -axis. This form is in agreement with the boundary conditions at infinity that disallow higher order terms in . Hence, the form of generalized solution using the complex potential has been established.

The existence and uniqueness of a solution to the boundary value problem originated from quasicrystal linear elasticity have been pointed out by Fan on the basis of (1)–(7); details are in [36].

3. Boundary Value Problems of Decagonal QCs in Virtue of Dugdale’s Manner

In the linear elasticity, the stresses at the tip of a crack in a body are singular and the stresses in the neighbourhood of the crack tip may exceed the yield stress of the material constituting the body whereas inelastic deformations are nowhere permissible. To avoid this paradox, Dugdale suggested a crack model which allows the occurrence of yielding and gives an estimate for the spread of plasticity in terms of external applied loading [39]. The cohesive zone is defined as a transition region in which the surface traction smoothly changes from zero at the crack tip to a certain magnitude at the cohesive zone tip, where the displacement discontinuity across the zone disappears. The cohesive zone is a mathematical extension of the crack and can be viewed physically as the fracture process zone. Near the crack tip there is a so-called cohesive force zone whose value is unknown temporarily; this zone is subjected to stress (phonon field), in which represents the stress yield strength of the material. In essence, the Dugdale model may be regarded as a superposition of two linear elastic crack solutions: firstly, the stress field resulting from the uniform traction which tends to open the crack and secondly the surface stress over the yield zones, which is equal in magnitude to the material’s yield stress and act to prevent this opening action such that the stresses are bounded at the crack tip. Based on the concept of the Dugdale model [39], the unknown strip yield zone sizes are determined by the condition of vanishing stress singularity at the crack tip; namely, where is the stress intensity factor of fictitious crack tip under uniform pressure stress . And, represents the stress intensity factor due to the shear yield strength , which can be determined by the experimental results.

3.1. A Bending Specimen of Point Group 10 mm Decagonal QCs with a Crack

Consider a bending specimen made of point group 10 mm decagonal QCs, weakened by a crack, as shown in Figure 1.

We do not have yet ideas of loading devices for the phason loads, so tractions in the phason field are assumed equal to zero on the boundary of a QC, from a theoretical point of view. Then the following boundary conditions should be satisfied:in which denotes the bending moment subjected to the upper and lower surface of the specimen, and is the thickness and the width of the specimen. We assume that the plastic zone around crack tip is of a narrow band type, with length denoted by : its value to be determined. The difficulty of the problem at hands rests essentially on the complicated boundary conditions. The exact analytic solution is not available. We here can obtain an approximate solution. For this purpose we decompose the problem into two secondary problems, one corresponding to the boundary conditions:and the other is related to the boundary conditions:The boundary value problem (8) ((23a)–(23c)) is a pure bending problem of cracked specimen. Muskhelishvili studied a bending specimen with an elliptic notch, but the major axis of the ellipse is parallel to the -axis; the solution is no usage to the present problem; the methodology of Muskhelishvili is still effective when the condition is satisfied [37]. By making these simplifications, the boundary value problems also cannot be solved at the physical plane (i.e., the -plane). We need to employ a conformal mapping to transform the problem for solving at mapping plane. Here we must take an approximation and assume that the size of crack is quite smaller than that of the specimen, and we can use the conformal mapping to transform the region of the specimen onto the interior of the unit circle at the -plane; we find the approximate solution at the conformal mapping plane by solving (15a) and the boundary conditions (23a)–(23c) such as In the equation some terms being independent from crack stress state are not included.

The result related to the boundary value problem (24a)–(24d) is obtained; then, where represents the value of corresponding point at the unit circle , represents the conjugate of , and corresponds to the point at the physical plane and represents the conjugate of . It is evident that and . The solution (27) is constructed for the infinite plane with a Dugdale crack, which approximately holds for the specimen with finite width . Since the complex analysis is well known, the detail of above calculation is omitted. In light of (26) and (25), meanwhile noting that (16), we have From (27) and (25) one finds that Because at the tip of the cohesive force zone there is no stress singularity; that is, the stress intensity factor at the point should be zero; that is, the requirement of (21) leads to The size of the plastic zone is determined already, though the relation is a nonlinear equation obviously. Due to and , the formula (30) is highly approximated. To obtain the stresses and displacements, the solution of and is enough. All can be determined on basis of (11a)–(11i) and (13a)–(13e).

3.2. A Finite Width Strip of Decagonal QCs with a Single Edge Crack

Consider a strip with a single crack embedded in decagonal QCs, in which stands for the width. Figure 2 shows a scheme of the crack containing a cohesive zone ahead of the crack tip.

A transcendental function conformal mapping constructed corresponding to the transformation between the region in the physical plane and the upper half-plane in -plane reads It can transform the cracked material shown in Figure 1 onto the upper-half plane of the mapping plane. Under the mapping (31), the crack tip is mapped to , while are mapped to . As shown in Figure 2, the boundary conditions of this problem of QCs can be written as follows:

Meanwhile, based on the above figure of a finite width strip with single edge crack of two-dimensional QCs and the boundary conditions, the functional equations can be established by means of (15a) and its conjugated equation: in which represents the value of at the real axis in the mapping plane, and

By solving the foregoing equations, one gets the result Finally, substituting (36) and (31) into (16) yields the exact stress intensity factors: Similarly, if the uniform pressure on the surface of the crack is replaced by a shear force , we can find the stress intensity factors: The stress intensity factor due to the yield strength can be expressed by in which can be determined like . In the same way the requirement of (21) leads to the size of cohesive force zone: in which

4. Conclusion and Discussion

This work employs the complex potentials and conformal mapping techniques to analyze the equilibrium boundary value problems for elastic quasicrystals. In virtue of the generalized solutions, stress and displacement boundary value problems of decagonal QCs are transferred to boundary value problems of analytic functions in a given region. Two separated problems emerge and can be analyzed. The stress intensity factors can be expressed in exact analytic forms. Although there exists phason field in QCs, (28), (29) and (37), (38) show that the intensity factors of phonon and phason fields are independent of material constants. Since both phonon and phason stress and displacement fields are dependent on all of the material constants, this result indicates the important effect of phason and phonon-phason coupling to the fracture parameters. The Dugdale model is a simple and efficient method to study the ideal elastic-plastic cracks. The results of these problems studied in this paper show that phason field and phonon-phason coupling have no effect on the plastic zone size, which is useful for plastic deformation of decagonal QCs. All the results can be reduced to the corresponding results of linear elastic theory when the phason elastic constants are approaching to zero. This work reveals the physical sense of the results relative to phason and the difference between the mechanical behaviors of the crack problem in crystals and decagonal QCs. It is presented that the mechanic problems can be transferred to boundary value problems of analytic functions.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant no. 11402158) and the Qualified Personnel Foundation of Taiyuan University of Technology (Grant no. tyut-rc201358a).

References

D. Shechtman, I. Blech, D. Gratias, and J. W. Cahn, “Metallic phase with long-range orientational order and no translational symmetry,” Physical Review Letters, vol. 53, no. 20, pp. 1951–1953, 1984.
View at: Publisher Site | Google Scholar
S. Fischer, A. Exner, K. Zielske et al., “Colloidal quasicrystals with 12-fold and 18-fold diffraction symmetry,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, no. 5, pp. 1810–1814, 2011.
View at: Publisher Site | Google Scholar
R. Mikulla, J. Stadler, and H. R. Trebin, “Crack propagation in quasicrystals,” Physical Review Letters, vol. 81, pp. 3163–3166, 1998.
View at: Google Scholar
J. Bohsung and H. R. Trebin, “Disclinations in quasicrystals,” Physical Review Letters, vol. 58, no. 12, pp. 1204–1207, 1987.
View at: Publisher Site | Google Scholar | MathSciNet
K. Edagawa, “Phonon-phason coupling in decagonal quasicrystals,” Philosophical Magazine, vol. 87, no. 18–21, pp. 2789–2798, 2007.
View at: Publisher Site | Google Scholar
K. Edagawa, Dislocation in Solids, Elsevier, Amsterdam, The Netherlands, 2007.
C. Walz, Zur Hydrodynamik in Quasikristallen, Diplomarbeit, Universitaet Stuttgart, 2003.
H. Lipp, M. Engel, S. Sonntag, and H.-R. Trebin, “Phason dynamics in one-dimensional lattices,” Physical Review B, vol. 81, no. 6, Article ID 064302, 2010.
View at: Publisher Site | Google Scholar
M. Engel, M. Umezaki, H.-R. Trebin, and T. Odagaki, “Dynamics of particle flips in two-dimensional quasicrystals,” Physical Review B—Condensed Matter and Materials Physics, vol. 82, no. 13, Article ID 134206, 2010.
View at: Publisher Site | Google Scholar
F. Rösch, C. Rudhart, J. Roth, H.-R. Trebin, and P. Gumbsch, “Dynamic fracture of icosahedral model quasicrystals: a molecular dynamics study,” Physical Review B, vol. 72, no. 1, Article ID 014128, 2005.
View at: Publisher Site | Google Scholar
P. Bak, “Phenomenological theory of icosahedral incommensurate (“Quasiperiodic”) order in Mn-Al Alloys,” Physical Review Letters, vol. 54, no. 14, pp. 1517–1519, 1985.
View at: Publisher Site | Google Scholar
P. Bak, “Symmetry, stability, and elastic properties of icosahedral incommensurate crystals,” Physical Review B: Condensed Matter. Third Series, vol. 32, no. 9, pp. 5764–5772, 1985.
View at: Publisher Site | Google Scholar | MathSciNet
T. C. Lubensky, S. Ramaswamy, and J. Toner, “Hydrodynamics of icosahedral quasicrystals,” Physical Review B, vol. 32, no. 11, pp. 7444–7452, 1985.
View at: Publisher Site | Google Scholar
D.-H. Ding, W. Yang, C. Hu, and R. Wang, “Generalized elasticity theory of quasicrystals,” Physical Review B, vol. 48, no. 10, pp. 7003–7010, 1993.
View at: Publisher Site | Google Scholar
C. Hu, R. Wang, and D. Ding, “Symmetry groups, physical property tensors, elasticity and dislocations in quasicrystals,” Reports on Progress in Physics, vol. 63, no. 1, pp. 1–39, 2000.
View at: Publisher Site | Google Scholar | MathSciNet
W. Li and T. Y. Fan, “Dynamic analysis of cracked octagonal quasicrystals,” Mathematical Problems in Engineering, vol. 2011, Article ID 567489, 18 pages, 2011.
View at: Publisher Site | Google Scholar | MathSciNet
S. Francoual, F. Livet, M. De Boissieu et al., “Dynamics of phason fluctuations in the i-AlPdMn quasicrystal,” Physical Review Letters, vol. 91, Article ID 225501, 2003.
View at: Google Scholar
G. Coddens, “On the problem of the relation between phason elasticity and phason dynamics in quasicrystals,” European Physical Journal B, vol. 54, no. 1, pp. 37–65, 2006.
View at: Publisher Site | Google Scholar
J. B. Wang, L. Mancini, R. H. Wang, and J. Gastaldi, “Phonon- and phason-type spherical inclusions in icosahedral quasicrystals,” Journal of Physics Condensed Matter, vol. 15, no. 24, pp. L363–L370, 2003.
View at: Publisher Site | Google Scholar
M. Feuerbacher, M. Bartsch, B. Grushko, U. Messerschmidt, and K. Urban, “Plastic deformation of decagonal Al-Ni-Co quasicrystals,” Philosophical Magazine Letters, vol. 76, no. 6, pp. 369–376, 1997.
View at: Publisher Site | Google Scholar
M. Feuerbacher, B. Baufeld, R. Rosefeld et al., “Intrinsic deformation properties of icosahedral Al-Pd-Mn single quasicrystals,” Philosophical Magazine Letters, vol. 71, pp. 91–98, 1995.
View at: Google Scholar
M. Wollgarten, M. Beyss, K. Urban, H. Liebertz, and U. Köster, “Direct evidence for plastic deformation of quasicrystals by means of a dislocation mechanism,” Physical Review Letters, vol. 71, no. 4, pp. 549–552, 1993.
View at: Publisher Site | Google Scholar
R. Rosenfeld and M. Feuerbacher, “Study of plastically deformed icosahedral Al-Pd-Mn single quasicrystals by transmission electron microscopy,” Philosophical Magazine Letters, vol. 72, pp. 375–784, 1995.
View at: Google Scholar
P. Schall, M. Feubrbacher, M. Bartsch, U. Messerschmidt, and K. Urban, “Dislocation density evolution upon plastic deformation of Al-Pd-Mn single quasicrystals,” Philosophical Magazine Letters, vol. 79, no. 10, pp. 785–796, 1999.
View at: Publisher Site | Google Scholar
B. Geyer, M. Bartsch, M. Feuerbacher, K. Urban, and U. Messerschmidt, “Plastic deformation of icosahedral Al-Pd-Mn single quasicrystals I. Experimental results,” Philosophical Magazine A, vol. 80, no. 5, pp. 1151–1163, 2000.
View at: Google Scholar
U. Messerschmidt, M. Bartsch, M. Feuerbacher, B. Geyer, and K. Urban, “Friction mechanism of dislocation motion in icosahedral Al-Pd-Mn quasicrystals,” Philosophical Magazine A, vol. 79, no. 9, pp. 2123–2135, 1999.
View at: Google Scholar
P. Ebert, M. Feuerbacher, N. Tamura, M. Wollgarten, and K. Urban, “Evidence for a cluster-based structure of AIPdMn single quasicrystals,” Physical Review Letters, vol. 77, no. 18, pp. 3827–3830, 1996.
View at: Publisher Site | Google Scholar
H. Takakura, M. Shiono, T. J. Sato, A. Yamamoto, and A. P. Tsai, “Ab initio structure determination of icosahedral Zn-Mg-Ho quasicrystals by density modification method,” Physical Review Letters, vol. 86, no. 2, pp. 236–239, 2001.
View at: Publisher Site | Google Scholar
D. Caillard, G. Vanderschaeve, L. Bresson, and D. Gratias, “Transmission electron microscopy study of dislocations and extended defects in as-grown icosahedral Al-Pd-Mn single grains,” Philosophical Magazine A, vol. 80, no. 1, pp. 237–253, 2000.
View at: Google Scholar
X. Y. Li and P. D. Li, “Three-dimensional thermo-elastic general solutions of one-dimensional hexagonal quasi-crystal and fundamental solutions,” Physics Letters A: General, Atomic and Solid State Physics, vol. 376, no. 26-27, pp. 2004–2009, 2012.
View at: Publisher Site | Google Scholar
X.-Y. Li, “Fundamental solutions of penny-shaped and half-infinite plane cracks embedded in an infinite space of one-dimensional hexagonal quasi-crystal under thermal loading,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 469, no. 2154, Article ID 20130023, 2013.
View at: Publisher Site | Google Scholar
E. Radi and P. M. Mariano, “Stationary straight cracks in quasicrystals,” International Journal of Fracture, vol. 166, pp. 102–120, 2010.
View at: Google Scholar
E. Radi and P. M. Mariano, “Steady-state propagation of dislocations in quasi-crystals,” Proceedings of The Royal Society of London A, vol. 467, no. 2136, pp. 3490–3508, 2011.
View at: Publisher Site | Google Scholar | MathSciNet
L. H. Li and G. T. Liu, “Stroh formalism for icosahedral quasicrystal and its application,” Physics Letters A, vol. 376, no. 8-9, pp. 987–990, 2012.
View at: Publisher Site | Google Scholar
J. H. Guo, J. Yu, and Y. M. Xing, “Anti-plane analysis on a finite crack in a one-dimensional hexagonal quasicrystal strip,” Mechanics Research Communications, vol. 52, pp. 40–45, 2013.
View at: Publisher Site | Google Scholar
T. Y. Fan, Mathematical Theory of Elasticity of Quasicrystals and Its Applications, Springer, Heideberg, Germany, 2010.
N. I. Muskhelishvili, Some Basic Problems of Mathematical Theory of Elasticity, Noordhoff, Groningen, The Netherlands, 1963.
View at: MathSciNet
P. M. Mariano and J. Planas, “Phason self-actions in quasicrystals,” Physica D, vol. 249, pp. 46–57, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
D. S. Dugdale, “Yielding of steel sheets containing slits,” Journal of the Mechanics and Physics of Solids, vol. 8, no. 2, pp. 100–104, 1960.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2014 Wu Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1143

Downloads

1033

Citations