Explicit Form of Exact Analytical Solution for Calculating Ground Displacement and Stress Induced by Shallow Tunneling and Its Application

Guo, Caixia; Han, Kaihang; Kong, Heng; Shi, Leilei

doi:https://doi.org/10.1155/2019/5739123

Advances in Civil Engineering

On this page

Abstract Introduction Conclusion Appendix Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 5739123 | https://doi.org/10.1155/2019/5739123

Explicit Form of Exact Analytical Solution for Calculating Ground Displacement and Stress Induced by Shallow Tunneling and Its Application

Caixia Guo,¹Kaihang Han,^2,3Heng Kong,⁴and Leilei Shi⁵

Academic Editor: Claudio Tamagnini

Received28 Apr 2018

Accepted23 Oct 2018

Published28 Mar 2019

Abstract

In urban environment, it is often unavoidable for shallow tunnels to be constructed adjacent to existing pile foundations. To obtain the ground displacements and stresses induced by shallow tunneling and existing pile foundation loads, the key procedure involves superimposing the analytical solution for shallow tunneling in green-field with the analytical solution for existing structure loads. In green-field, the complex variable method provides exact analytical solutions of ground displacements and stresses caused by shallow tunneling. However, the exact analytical solutions are not directly expressed as explicit functions of the coordinates (x, y) in the physical plane (called implicit form of exact analytical solutions), whereas the displacements and stresses induced by existing structure loads are explicit functions of the coordinates (x, y) in the physical plane, which makes it difficult to superpose the displacements and stresses induced by existing structure loads. In this paper, explicit form of exact analytical solutions of displacements and stresses induced by shallow tunneling in green-field is obtained by using the inverse conformal transformation and the Cauchy–Riemann equations. Comparison with implicit form of exact analytical solutions shows that the explicit form of exact analytical solutions is intuitional and easily used by engineers, and moreover, the calculation amount is much smaller than that for the implicit form of exact analytical solutions. Then, an application involving superimposing the explicit form of exact analytical solutions with Mindlin’s solution is implemented to analyze the secondary stress field and the related potential plastic zone caused by shallow tunneling adjacent to pile foundations. Moreover, the influences of pile foundation parameters on the ranges and shapes of the potential plastic zones induced by nearby tunneling are analyzed.

1. Introduction

Due to continuous expansion of cities and decreases in available land, the demand for urban tunnels (roads or metro) has increased sharply in recent years. At the same time, the foundations of nearby existing structures are unavoidably disturbed by this process and affect it, especially in densely constructed areas located adjacent to the tunnel construction site. The prediction and mitigation of damage caused by construction-induced ground movement represents a major factor in the design of tunnels. This is an especially important problem for shallow tunnels excavated in soft soils, where expensive remedial measures such as compensation grouting or structural underpinning may need to be considered prior to construction. To obtain the ground displacements and stresses induced by shallow tunneling and existing structure loads, the key procedure involves superimposing the analytical solution for tunneling in green-field with the analytical solution for existing structure loads.

To estimate the ground displacements and stresses caused by shallow tunneling in green-field, a number of empirical and analytical solutions have been developed. In engineering practice, these ground surface settlements are often described by Peck’s empirical formula [1, 2], which provides a Gaussian distribution curve for ground settlement and is based on field observations and intuitive deductions, thus it is lacking in theory and ambiguous in its range of applicability. In addition, there are four main analytical methods, namely, the virtual image technique [3–5], the complex variable method [6, 7], the general series form stress function in polar coordinates [8–11], and the stochastic medium theory [12]. Pinto and Whittle [13] and Xiang [14] presented detailed reviews and comparisons of the empirical and analytical solutions of ground displacements and stresses for shallow tunnels in green-field. However, those solutions mentioned above are in general all based upon the premise that the ground is free from existing structure loads, and they cannot directly predict ground displacements and stresses induced by the combination of shallow tunneling and existing structure loads.

For the condition of shallow tunneling adjacent to an existing pile foundation, Xiang and Feng [15] proposed a superposition method for predicting the potential plastic zone of shallow tunneling adjacent to a pile foundation in soils. Although their results achieve good effects, there is still one approximation. According to Verruijt [6] and Pinto and Whittle [13], the tunneling-induced stresses adopted by Xiang and Feng [15] are an approximate solution that implicitly ignores the finite dimensions of the tunnel itself, and the results from approximate solution would lead to certain differences, especially for tunnels too close to the ground surface. This approximation would make the superposition method proposed by Xiang and Feng [15] unsuitable for tunnels too close to the ground surface. Referring to Xiang [14] and Pinto and Whittle [13], the complex variable method not only provides exact analytical solutions of ground displacements and stresses caused by shallow tunneling in green-field but also makes the exact analytical solutions suitable for various types of boundary conditions. Therefore, to obtain more accurate results for secondary stress field and the related potential plastic zone, the exact analytical solutions based on the complex variable method are used for the superposition method in this paper. However, there is still a problem that exact analytical solutions are not directly expressed as explicit functions of the coordinates (x, y) in the physical plane (called implicit form of exact analytical solutions), whereas the displacements and stresses induced by existing structure loads are explicit functions of the coordinates (x, y) in the physical plane, which makes it difficult to superpose the displacements and stresses induced by existing structure loads.

The present paper extends and partially revises past research work [16]. The implicit form of exact analytical solutions, which are based on the complex variable method, is derived into explicit form of exact analytical solutions by using the inverse conformal transformation and the Cauchy–Riemann equations. Then, an application is conducted to analyze the secondary stress field and the related potential plastic zone caused by shallow tunneling adjacent to pile foundations. Moreover, the influence of pile foundation parameters (pile length, load magnitude, and pile offsets) on the ranges and shapes of the potential plastic zones induced by nearby tunneling is analyzed.

2. Explicit Form of Exact Analytical Solutions of Ground Displacements and Stresses Induced by Shallow Tunneling

2.1. Implicit Form of Exact Analytical Solutions

The problem considers an elastic half-plane with a circular tunnel (Figure 1) [6]. The radius of the tunnel is expressed by r, the depth of its center below the free surface by h, and the cover by d. The ground surface boundary is free of stress and the boundary of the tunnel undergoes a given distribution of displacements (for instance, two typical boundary conditions: uniform radial displacement (the ground loss problem) u₀ and ovalization u_d). In the complex variable method used by Verruijt [6], it is assumed that the original domain in the z-plane (physical plane) is mapped conformally onto an annular region bounded by the circles and , where , on the auxiliary domain in the -plane (mapped plane) by the following conformal transformation:where is given by

In z-plane, the solutions are expressed in terms of two analytic functions and . The stresses and displacements are related to these functions as follows:where is the shear modulus of the elastic material, is related to Poisson’s ratio ν by , is the imaginary constant, and overbar is a complex conjugate.

By virtue of the conformal transformation function , the functions and can be considered as functions of as follows:

As the functions and are analytical, they can be expanded in Laurent series in the -space as follows:where the coefficients , and are determined with recursive relations derived from the boundary conditions. First, the coefficients , and are found with the specific displacement boundary condition at the boundary of the tunnel (for instance, convergence and ovalization) by using related coefficients in Fourier series terms.where the related coefficients are in Fourier series terms (Appendix A).

Second, the coefficients , and can be obtained with stress-free boundary condition at the ground surface by using the following recursive relations:

The analytical solutions of ground displacements and stresses induced by shallow tunneling are obtained with the analytic functions (equations (3a)–(3c)). Moreover, Pinto [17] claims that 10–15 terms of analytic functions (equations (3a)–(3c)) are sufficient to achieve accurate solutions for both the convergence and ovalization modes of deformation.

However, the solutions proposed by Verruijt [6] (equations (3a)–(3c)) are expressed as two analytical functions ( and ) and are not directly expressed as explicit functions of the coordinates (x, y) in z-space (referred to as “implicit form of exact analytical solutions”), which makes the extensive application of implicit exact analytical solutions inconvenient. On the one hand, implicit form of exact analytical solutions is somewhat complicated and not intuitive for engineers. Implicit form of exact analytical solutions cannot be directly used by engineers like Peck empirical formulas if they have no certain professional theoretical knowledge of complex variable method. On the other hand, the superposition applicability of implicit form of exact analytical solutions is poor when solving the ground displacements and stresses induced by shallow tunneling adjacent to existing structures (surface building or pile foundation).

2.2. Explicit Form of Exact Analytical Solution

In this section, the main work is to derive the explicit form of exact analytical solution of the stress and the displacement expressed by the coordinates (x, y) in z-space (physical plane).

2.2.1. Series Form of and with the Functions of Coordinate (x, y) in z-Space

The corresponding inverse conformal transformation of equation (1) is

Considering , equation (8) can be transformed intowhere

To obtain the power of (i.e., and ), it is better to transform equation (9) into the triangle form of as follows:

According to the related operation rules of the power of complex numbers,

On the basis of a consideration of symmetry, Verruijt [6] assumed that all the coefficients are purely imaginary:where , and are the real parts of purely imaginary , and , respectively.

With equations (12a), (12b), and (13), equations (5a) and (5b) can be transformed into

Equations (14a) and (14b) can be simplified to the form of real and imaginary parts of and as follows:where

With equations (4a) and (4b), equations (15a) and (15b) can be written into equations (17a) and (17b):

2.2.2. Derivations of Analytic Functions Based on the Cauchy–Riemann Equations

if the function is analytical, the derivation of analytical function can be calculated with

Because and are analytical, the derivations of and can be obtained:

2.2.3. Explicit Analytical Solutions of Ground Displacements and Stresses

By taking equations (19a)–(19c) into equations (3a)–(3c), the explicit exact analytical solutions of ground displacements and stresses can be obtained:

The explicit form of exact analytical solutions (equations (20a)–(20e)) relates to only the original domain in the z-plane and is directly expressed as functions of coordinates (x, y) in z-space. The superscript ST indicates that the reason of induced ground displacements and stresses is shallow tunneling.

3. Application of Explicit Form of Exact Analytical Solutions: Secondary Stress Field and Potential Plastic Zone of Shallow Tunneling Adjacent to a Pile Foundation

3.1. Conceptual Model and Calculation Procedure

To predict the degree and extent of tunneling effects on a pile foundation, Xiang and Feng [15] proposed a theoretical superposition method for predicting the potential plastic zone of shallow tunneling adjacent to a pile foundation in soils. The practical problem of a shallow tunneling project adjacent to a pile foundation is simplified into the mechanics model shown in Figure 2.

In this paper, to obtain more accurate results for the secondary stress field and the related potential plastic zone, a similar calculation procedure is adopted. The differences between the calculation procedure for the superposition method proposed by Xiang and Feng [15] and that for the new method used in this paper are detailed below. It is worth noting that Mindlin’s solution, which is used to calculate the ground stresses due to pile foundation loads, is an exact solution and is preserved in the new superposition method. The cause of the inaccuracy of the superposition method proposed by Xiang and Feng [15] is mainly due to the approximate solution [5] that is used to calculate the ground displacements induced by shallow tunneling in green-field. Therefore, the improvement of this paper is to adopt an exact solution (explicit form of exact analytical solutions) for calculating tunneling-induced stress instead of the approximate solution used by Xiang and Feng [15].

Through the above explicit derivation, the explicit form of exact analytical solutions are expressed as functions of the coordinates (x, y), which can be directly superimposed with the stress induced by pile foundation loads. By superimposing these several equations, the envelope of the potential plastic zone induced by tunneling adjacent to the pile is easily obtained, which highlights the convenience of explicit form of exact analytical solutions. Then, by using the software MATLAB, the theoretical procedure described above is conducted successfully.

Referring to Xiang and Feng [15], the assumed parameters in all the presented calculations are as follows: silt clayey soil (E = 10 MPa, = 0.25, c = 30 kPa, and φ = 30); a circular tunnel R = 3 m, h = 8 m, 10 m, uniform convergence u₀ = 30 mm, oval deformation u_d = 10 mm, s = 78 kN/m, and .

3.2. Comparison with the Results of Xiang and Feng

To detect any differences in the plastic zone in green-field or plastic zone with pile load calculated by the different methods, a comparison between the results of this paper and those of Xiang and Feng [15] is performed for two relative tunnel depths (the ratio of cover depth to tunnel diameter) (Figures 3 and 4).

(a)

(b)

(a)

(b)

In Figure 3, comparisons of the potential plastic zones in green-field from the two superposition methods are presented for two relative tunnel depths. It is shown that the horizontal range between the two plastic zones in green-field with a relative depth of 0.83 differs by 38.7 cm in Figure 3(a), whereas the horizontal range with a relative depth of 1.17 differs by 27.8 cm in Figure 3(b).

In Figure 4, comparisons of the potential plastic zones with pile load from the two superposition methods are presented for two relative tunnel depths. There exists a more marked difference between the ranges and shapes of the two plastic zones with pile load with a relative depth of 0.83 in Figure 4(a). The difference between the range and shape of two plastic zones with pile load is greater with a relative depth of 0.83 than with a relative depth of 1.17 in Figure 4(b).

As a conclusion, the plastic zones around relatively deep tunnels adjacent to a pile foundation obtained from the two superposition methods are similar to each other, whereas the plastic zones around relatively shallow tunnels adjacent to a pile foundation from the two superposition methods are obviously different from each other. In other words, the superposition method proposed by Xiang and Feng [15] can obtain satisfactory results for deep tunnels, whereas the superposition method used in this paper is suitable for both deep tunnels and shallow tunnels because it is an exact analytical solution.

3.3. Influences of Different Pile Foundation Parameters

The influences of different pile foundation parameters (pile length, load magnitude, and pile offsets) on the ranges and shapes of the potential plastic zones induced by nearby tunneling are analyzed for a relative depth of 1.17 in Figure 5. The results indicate that the whole potential plastic zones induced by tunneling and pile foundation loads would coalesce when the pile is located close enough to the tunnel, whereas the two plastic zones are separated from each other.

(a)

(b)

(c)

3.4. Influences of Different Tunnel Boundary Conditions

It should be noted that the tunnel boundary conditions determine the calculation results of the ranges and shapes of the potential plastic zones caused by shallow tunneling adjacent to a pile foundation in soils. To compare the calculation results with those by Xiang and Feng [15] under the same conditions, the same boundary conditions of a shallow tunnel are adopted, which are a combination of uniform convergence and ovalization proposed by Verruijt and Booker [5] and shown in Figure 6. Pinto and Whittle [13] summarized three types of shallow tunnel boundary conditions. They claimed that the vertical translation of a shallow tunnel should be incorporated with uniform convergence and ovalization when considering the buoyancy effect, usually induced because the weight of the tunnel is usually less than the weight of the excavated soil. And by incorporating the vertical translation with the uniform convergence and ovalization of the tunnel, Park [10] introduced four types of complex boundary conditions, and those complex boundary conditions are more in accordance with engineering practice, as shown in Figures 7 and 8. In fact, the boundary condition B.C.-1 is the same as the boundary condition discussed by Verruijt and Booker [5], whereas the boundary conditions B.C.-2, B.C.-3, and B.C.-4 are different and contain vertical translation of a shallow tunnel.

(a)

(b)

As mentioned above, the part of the superposition method proposed by Xiang and Feng [15] that is used to calculate the ground displacement induced by tunneling in green-field adopts the approximate solution proposed by Verruijt and Booker [5]; thus, the superposition method proposed by Xiang and Feng [15] can adopt only the combination of uniform convergence and ovalization as the boundary conditions of a tunnel. However, the superposition method used in this paper is based on the explicit form of exact analytical solutions, which allows the solution to be suitable for various types of boundary conditions. In other words, the superposition method used in this paper can use not only the boundary conditions proposed by Verruijt and Booker [5] as the boundary conditions of a shallow tunnel but also the complex boundary conditions proposed by Park [10] as the boundary conditions of a shallow tunnel. The related coefficients A_k in Fourier series terms [18] in equations (6a)–(6d) are shown in Appendix A.

The influences of different tunnel boundary conditions (for instance, B.C.-3 and B.C.-4) on the ranges and shapes of the potential plastic zones are indicated in Figure 6. It can be concluded from comparison with Figures 3–5 that the plastic zones with boundary conditions (B.C.-3 and B.C.-4) are different from the plastic zones with the boundary conditions proposed by Verruijt and Booker [5], especially in the lower part of the plastic zone. In other words, it would overestimate the plastic zones if the buoyancy effect of a shallow tunnel is not considered. The main reason for this difference is that vertical translation is not considered in the boundary conditions proposed by Verruijt and Booker [5].

4. Conclusion

In green-field, the complex variable method provides exact analytical solutions of ground displacements and stresses caused by shallow tunneling. However, the exact analytical solutions [6] are not directly expressed as explicit functions of the coordinates (x, y) in the physical plane (called implicit form of exact analytical solutions), whereas the displacements and stresses induced by existing structure loads are explicit functions of the coordinates (x, y) in the physical plane, which makes it difficult to superpose the displacements and stresses induced by existing structure loads. This paper transforms implicit form of exact analytical solutions into explicit form of exact analytical solutions, which improves the superposition applicability of exact analytical solution with the analytical solution for the existing structure load. With the explicit form of exact analytical solution, the secondary stress field and the related potential plastic zone caused by tunneling adjacent to pile foundations are obtained. The main conclusions are presented as follows:(1)With the inverse conformal transformation, the series forms of complex potential functions in z-plane are obtained. By taking the derivative of the analytic functions proposed by Verruijt [5] with the Cauchy–Riemann equations, the explicit form of exact analytical solutions of displacement and stresses induced by shallow tunneling is obtained. The explicit form of exact analytical solutions is intuitional and easily used by engineers, and the calculation amount is smaller than that for the implicit analytical solutions through comparison with the implicit form of exact analytical solutions.(2)An application involving superimposing the explicit form of exact analytical solutions with Mindlin’s solution [1] is implemented to analyze the secondary stress field and the related potential plastic zone caused by tunneling adjacent to pile foundations. A comparison between the results of this paper and one of the existing approaches proposed by Xiang and Feng [15], which is an approximate solution, is performed. The plastic zones around a relatively deep tunnel adjacent to a pile foundation obtained from the two superposition methods are similar to each other, whereas the plastic zones around a relatively shallow tunnel adjacent to a pile foundation from the two superposition methods are obviously different from each other. In other words, the superposition method proposed by Xiang and Feng [15] can obtain satisfactory results for deep tunnels, whereas the superposition method used in this paper is suitable for both deep tunnels and shallow tunnels because it is an exact analytical solution.(3)The influences of different pile foundation parameters (pile length, load magnitude, and pile offset) on the ranges and shapes of the potential plastic zones induced by nearby tunneling are also analyzed. The results indicate that the whole tunneling-induced potential plastic zones induced by tunneling and pile foundation loads around the tunnel and around the pile would coalesce when the pile is located close enough to the tunnel, whereas if the pile is far enough away from the tunnel, the two plastic zones due to tunneling-induced stress changes are separated from each other.(4)The superposition method used in this paper can use not only the simple boundary conditions proposed by Verruijt and Booker [5] but also the complex boundary conditions proposed by Park [10]. Therefore, the explicit form of exact analytical solution for calculating ground displacement and stress induced by shallow tunneling proposed in this paper has more extensive adaptability, so it can solve the more complex problems of shallow tunnels. For example, these solutions can be used to analyze differences in the plastic zone under the influence of the buoyancy effect. The results show that it would overestimate the plastic zones if the buoyancy effect of a shallow tunnel is not considered, especially in the lower part of the plastic zone.

Appendix

A. Fourier Coefficients for Boundary Deformations of a Tunnel

(1)Uniform convergence [6](2)Ovalization [13](3)B.C.-2 [10, 18](4)B.C.-3 [10, 18](5)B.C.-4 [10, 18]

Data Availability

An executable file to create the potential plastic zone caused by shallow tunneling adjacent to pile foundations and the executable file used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors acknowledge the financial support provided by the Fundamental Research Funds for the Central Universities of China (Grant no. 2015YJS128) and the National Key Research Program of China (2017YFC0805000).

References

R. J. Mair and R. N Taylor, “Theme lecture: bored tunneling in the urban environment,” in Proceedings of the Fourteenth International Conference on Soil Mechanics and Foundation Engineering, pp. 2353–2385, Balkema, Hamburg, Germany, September 1997.
View at: Google Scholar
R. B. Peck, “Deep excavations and tunneling in soft ground,” in Proceedings of the 7th International Conference on Soil Mechanics and Foundation Engineering, pp. 225–290, Mexico City, Mexico, 1969.
View at: Google Scholar
N. Loganathan and H. G. Poulos, “Analytical prediction for tunneling-induced ground movements in clays,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 124, no. 9, pp. 846–856, 1998.
View at: Publisher Site | Google Scholar
C. Sagaseta, “Analysis of undraind soil deformation due to ground loss,” Géotechnique, vol. 37, no. 3, pp. 301–320, 1987.
View at: Publisher Site | Google Scholar
A. Verruijt and J. R. Booker, “Surface settlements due to deformation of a tunnel in an elastic half plane,” Geotechnique, vol. 46, no. 4, pp. 751–756, 1996.
View at: Publisher Site | Google Scholar
A. Verruijt, “A complex variable solution for a deforming circular tunnel in an elastic half-plane,” International Journal for Numerical and Analytical Methods in Geomechanics, vol. 21, no. 2, pp. 77–89, 1997.
View at: Publisher Site | Google Scholar
A. Verruijt, “Deformations of an elastic half plane with a circular cavity,” International Journal of Solids and Structures, vol. 35, no. 21, pp. 2795–2804, 1998.
View at: Publisher Site | Google Scholar
A. Bobet, “Analytical solutions for shallow tunnels in saturated ground,” Journal of Engineering Mechanics, vol. 127, no. 12, pp. 1258–1266, 2001.
View at: Publisher Site | Google Scholar
W. I. Chou and A. Bobet, “Prediction of ground deformations in shallow tunnels in clay,” Tunnelling and Underground Space Technology, vol. 17, no. l, pp. 3–19, 2002.
View at: Publisher Site | Google Scholar
K. H. Park, “Elastic solution for tunneling-induced ground movements in clays,” International Journal of Geomechanics, vol. 4, no. 4, pp. 310–318, 2004.
View at: Publisher Site | Google Scholar
K.-H. Park, “Analytical solution for tunnelling-induced ground movement in clays,” Tunnelling and Underground Space Technology, vol. 20, no. 3, pp. 249–261, 2005.
View at: Publisher Site | Google Scholar
J. S. Yang, B. C. Liu, and M. C. Wang, “Modeling of tunneling-induced ground surface movements using stochastic medium theory,” Tunnelling and Underground Space Technology, vol. 19, no. 2, pp. 113–123, 2004.
View at: Publisher Site | Google Scholar
F. Pinto and A. J. Whittle, “Ground movements due to shallow tunnels in soft ground. I: analytical solutions,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 140, no. 4, Article ID 04013040, 2013.
View at: Publisher Site | Google Scholar
Y. Y. Xiang, “An evaluation of theoretical solutions for prediction of ground movements due to shallow tunneling in soils,” in Proceedings of Underground Construction and Ground Movement, vol. 155, pp. 296–303, ASCE-Geotechnical Special Publication, Shanghai, China, June 2006.
View at: Publisher Site | Google Scholar
Y. Xiang and S. Feng, “Theoretical prediction of the potential plastic zone of shallow tunneling in vicinity of pile foundation in soils,” Tunnelling and Underground Space Technology, vol. 38, pp. 115–121, 2013.
View at: Publisher Site | Google Scholar
K. H. Han, C. P. Zhang, and M. S. Wang, “Explicit analytical solutions for stress and displacement of surrounding rock in shallow tunnels,” Chinese Journal of Geotechnical Engineering, vol. 36, no. 12, pp. 2253–2259, 2014, in Chinese.
View at: Google Scholar
F. Pinto, “Analytical methods to interpret ground deformations due to soft ground tunneling,” Deptartment of Civil and Environmental Engineering, Massachusetts Institute of Technology (MIT), Cambridge, MA, USA, 1999, M.S. thesis.
View at: Google Scholar
L.-Z. Wang, L.-L. Li, and X.-J. Lv, “Complex variable solutions for tunneling-induced ground movement,” International Journal of Geomechanics, vol. 9, no. 2, pp. 63–72, 2009.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2019 Caixia Guo 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

917

Downloads

1009

Citations