Research Article | Open Access

Fangyuan Zhou, Xuezhang Wen, Hongping Zhu, "Effect of Soil Structure Interaction on Torsional Response of Structure Supported by Asymmetric Soil Foundation System", *Shock and Vibration*, vol. 2016, Article ID 2647675, 14 pages, 2016. https://doi.org/10.1155/2016/2647675

# Effect of Soil Structure Interaction on Torsional Response of Structure Supported by Asymmetric Soil Foundation System

**Academic Editor:**Nerio Tullini

#### Abstract

The torsional response of a structure supported by asymmetric foundation was investigated in this study. Several types of the asymmetric soil foundation system were employed to analyze the effect of soil structure interaction on torsional response of the superstructure. It can be concluded from the study that torsional response would be generated for a structure supported by asymmetric soil foundation system under horizontal seismic excitation, and the generated torsional response of the superstructure changed with the degree of the asymmetry of the foundation.

#### 1. Introduction

In recent decades, several reports have pointed out that structure failures due to torsional response were common during previous earthquake [1, 2]. Meanwhile, through numerical analysis of the structure under seismic excitation, researchers found that asymmetric structures were subjected to more damage than symmetric structures. Chandler and Hutchinson [3] presented a detailed parametric study of the coupled lateral and torsional responses of a partially symmetric single storey building subjected to both harmonic and earthquake base excitation. It was concluded that torsional coupling induced a significant amplification of earthquake forces which should be accounted for in design. de La Llera and Chopra [4–6] calculated the value of the accidental eccentricity with the equivalent lateral force procedure by studying the dynamic response of single and multistory buildings subjected to torsional ground motion. Nagarajaiah et al. [7, 8], using multistory models, concluded that eccentricities in the isolation system and superstructure both contributed to the torsional response. Similarly Jangid and Datta [9, 10] found that significant eccentricity of the superstructure could reduce the effectiveness of the base isolation system. In Tena-Colunga et al. [11–13], parametric study was carried out with eccentricities in both the isolation system and superstructure, which concluded that eccentricity in the isolation system might lead to a torsional response. Most of these studies assumed base fixity in the structural models, neglecting the soil structure interaction (SSI).

Some studies showed that soil structure interaction considerably influenced the dynamic responses of the structures subjected to earthquake loading [14–16]. Although currently there exist adequate computational capabilities for building design, SSI analyses were rarely carried out for the seismic design of building structures [17].

Neglecting SSI in the seismic response analysis of building structures may result in a false torsional response of the structure. In many cases, fixed-base analyses cannot predict accurately the torsional response of building structures. As we know the soil properties and the foundation type can change the dynamic response of the superstructure. Thus, the assumption that the base of the structure is fixed is not realistic. Therefore, several researchers attempted to incorporate the flexibility of foundation in asymmetric system models. Among them, Balendra et al. [18] used simple springs to approximate the foundation impedance in an asymmetric multistory structure. Chandler and Hutchinson [19, 20] analyzed the seismic responses subjected to different types of earthquake loading considering soil structure interaction. Sivakumaran et al. [21] and Sivakumaran and Balendra [22] inspected the dynamic response of asymmetric multistorey structure-foundation systems using frequency independent foundation impedance.

Previous studies mainly focused on the torsional responses of structures induced by the asymmetry in superstructures. Little attention has been paid to the torsional responses of structures induced by asymmetric soil foundation system. In this study, the torsional response of a structure supported on asymmetric soil foundation system was investigated. Several types of asymmetric soil foundation were employed to analyze the effect of soil structure interaction on the torsional response of the superstructure.

#### 2. Analytical Theory

The analysis procedure in this paper was based on the following basic assumptions:(1)the soil was assumed to be elastic medium;(2)the foundation remains in close contact with the surrounding soil; no slippage was allowed at the soil foundation interface.For a structure (Figure 1) under seismic excitation, the equation of motion of the structure considering soil structure interaction can be depicted as formula (1) (the calculation was performed in the frequency domain, and the term is omitted in the formula):where . , , and are the stiffness matrices, mass matrices, and displacement vector, respectively. Subscripts and represent the superstructure and foundation, respectively. represents the damping ratio. and stand for the force vector acting at the superstructure and the foundation, respectively. For the seismic wave input problem, the force vector acting at the superstructure .

For the superstructure, the well-known beam model is one of the most commonly used models for conducting comprehensive parametric studies. In the MDOF structure models used in this study, each floor was assumed to be a lumped mass connected by a beam element. The stiffness matrices of the foundation and excavated soil were obtained by finite element method.

As shown in Figure 2, the structure-foundation-soil system was separated into three components, that is, structure-foundation, the free field, and excavated soil, using the concept of the substructure method [23].

The equation of motion of the free field can be depicted as the formulawhere character represents the free field, stands for the force vector acting at the nodes of the free field.

The stiffness matrix of the free field can be calculated with the following equations:In (4), Green’s function solution for a concentrated load and for a uniform distributed load acting on a disk (when the excitation node and reception node coincide with each other) was applied to form the flexibility matrix of the free field [24] (see Appendix A). At the same time, a paraxial boundary condition (see Appendix A) [25, 26] was used to simulate the bottom boundary of the layered soil.

The equation of motion of excavated soil can be depicted as the formula where . Superscript stands for the excavated soil which has same geometric dimension as the foundation and same physical parameters as the free field.

For the soil structure interaction system,Substituting formulae (1)~(5) into formula (6) producesFormula (7) is the equation of motion for a structure-foundation-soil system under seismic excitation. Through formula (7) the dynamic response of the structure under seismic excitation can be obtained.

#### 3. Numerical Study

In order to investigate the torsional response of the structure influenced by asymmetric soil foundation system, influences of the following asymmetric soil foundation systems were analyzed:(1)Embedded foundation supported by piles with nonunique pile lengths.(2)L shape and T shape embedded foundation supported by piles.(3)Foundation with an adjacent foundation.For all the cases discussed in this section, normal incident wave was used as an excitation source; the frequency of the excitation wave ranged from 0.1 Hz to 10 Hz.

A 10-storey structure was used in this section, the first, second, and third horizontal natural frequencies of the superstructure were 1.4, 4.2, and 6.8 Hz, respectively. The first, second, and third torsional natural frequencies of the superstructures were 1.9, 5.7 and 9.4 Hz, respectively.

##### 3.1. Embedded Foundation Supported by Piles with Nonunique Pile Lengths

Plane sketch of the employed embedded foundation supported by piles of the superstructure was depicted in Figure 3. The dimension of the embedded foundation was 72 m in length () and 54 m in width () and buried depth of the foundation () was 4 m. Piles were defined by diameter (1.5 m), center-to-center spacing between adjacent piles (9.0 m), and length . The material parameters of the pile were as follows: density = 2.4 t/m^{3}; moduli of elasticity of the pile = 2.1 × 10^{4} MPa; Poisson ratio of the pile = 1/6. The soil condition surrounding the foundation was shown in Figure 4. To ensure the accuracy of the thin layer method at low frequency, the relationship between the wavelength and the depth of the soil should be . Meanwhile, to ensure the accuracy of the thin layer method at high frequency, the relationship between the wavelength and the depth of the soil layer should be [27]. In this paper, the frequency domain of the excitation varied from 0.1 Hz to 10 Hz; thus, the soil profile shown in Figure 4 was employed.

In order to investigate the impact of irregular soil foundation system on the torsional response of the superstructure, the following cases as shown in Figure 5 were investigated:(1)Case : embedded foundation supported by piles with unique pile length.(2)Cases 2~3: embedded foundation supported by piles with nonunique piles lengths (nonunique piles were distributed along -axis).(3)Cases 4~5: embedded foundation supported by piles with nonunique piles lengths (nonunique piles were distributed along -axis).The following ratio was employed to investigate the torsional response of the superstructure: . Here is the torsional response of the structure; is the lateral free field displacement.

**(a) Embedded foundation supported by piles with unique pile length**

**(b) Embedded foundation supported by piles with nonunique pile length (nonunique piles distributed along -axis)**

**(c) Embedded foundation supported by piles with nonunique pile length (nonunique piles distributed along -axis)**

Torsional response of the roof floor (RF) of the superstructure for Cases 1~5 is depicted in Figure 6. It can be observed that when the excited direction was perpendicular to -axis, the torsional responses of the roof floor in Case and Case were generated, and the torsional response of the roof floor in Case and Case occurred when the excited direction was parallel to the -axis. However, the torsional response of the roof floor in Case cannot be observed under both directions of horizontal excitation. Through Figures 6(a) and 6(b), it can be seen that the torsional response of the roof floor increased with the length difference of the piles.

**(a) Excited perpendicular to direction**

**(b) Excited parallel to direction**

The above phenomenon can be explained as follows. When the embedded foundation was supported by piles with different pile lengths the mass center and the rigidity center of the soil foundation system did not coincide with each other; thus eccentricity existed in the soil foundation system. The foundation input motion of the embedded foundation supported by piles in Cases 1~5 is depicted in Figure 7. It can be concluded that in Case (symmetric soil foundation system) torsional foundation input motion was not generated, but in Cases 2~5 the torsional foundation input motion was generated due to the asymmetric soil foundation system. Therefore, the asymmetric soil foundation system can induce the torsional foundation input motion; thus for the structure supported by asymmetric soil foundation system, the torsional response of the superstructure is of concern.

**(a) Excited perpendicular to direction**

**(b) Excited parallel to direction**

##### 3.2. L Shape and T Shape Embedded Foundation Supported by Piles

Figure 8 shows the L shape (Cases 6~7) and T shape (Cases 8~10) embedded foundations supported by piles. The diameter and length of the piles used in this section were 1.5 m and 16 m, respectively. The distance between the adjacent piles was 9 m. The parameters of the superstructure and the soil were the same as those in Section 3.1. The plan dimension of the L shape and T shape embedded foundation supported by piles can be obtained through the distance between the two adjacent piles as shown in Figure 8. The torsional response of Cases 6~7 is compared with that of Case depicted in Figure 9(a). It can be observed from Figure 9(a) that for a 10-storey structure supported on L shape embedded foundation supported by piles, torsional response of the roof floor was generated regardless of the horizontal excitation direction, and it increased with the irregular degree of the L shape foundation. The torsional responses in Cases 8~10 are compared with that of Case depicted in Figure 9(b). From Figure 9(b), it can be concluded that the torsional response phenomenon can be observed in Cases 8~10 when the excitation direction was parallel to -axis. However, if the excitation direction was perpendicular to -axis, there was torsional response in Case but none in Case and Case because both the mass center and the rigidity center of the soil foundation system were on the -axis.

**(a) Case : L shape**

**(b) Case : L shape**

**(c) Case : T shape**

**(d) Case : T shape**

**(e) Case : T shape**

**(a) L shape**

**(b) T shape**

##### 3.3. Foundation with an Adjacent Foundation

The influence of an adjacent structure (named AS) on the torsional response of the targeted structure (named TS) is addressed in this section. The parameters of the soil were the same as those depicted in Figure 4. The dimension of the embedded foundation was 60 m in length () and 60 m in width () and the buried depth of the foundation () was 4 m. The material parameters of the pile were the same as those in Section 3.1. The diameter and the length of the piles were 1.5 m and 16 m, respectively. The distance between the two adjacent piles was 12 m. Figure 10(a) shows the sketch paragraph of a 10-storey structure (TS) adjacent to a 10-storey structure (AS) with varying distance () between the two adjacent foundations. Figure 10(b) shows the sketch paragraph of a 10-storey building (TS) adjacent to a 10-storey structure (AS) with varying placed angle () between the two adjacent foundations. The following cases are discussed: Case : m, ; Case : m, ; Case : m, ; Case : m, ; Case : m, .

**(a) Variation of distance**

**(b) Variation of angle (plan sketch)**

Torsional responses of the roof floor in Cases 11~13 compared with that in Case (excited direction was perpendicular to the -axis) are depicted in Figure 11(a). It can be observed that torsional response of the superstructure could be induced by the existence of the adjacent structure, and it increased with the decrease of the distance . The above phenomenon can be explained as follows. As the movement of the soil around the foundation of the targeted structure can be restricted by the presence of an adjacent structure, motion characteristic of the soil around the targeted structure at different spatial location may be impacted in different extent. Thus, the torsional response of the roof floor of the targeted structure can be observed under horizontal seismic excitation.

**(a) Variation of distance (excited perpendicular to -axis)**

**(b) Variation of angle**

In order to investigate the effect of the placed angle of the adjacent structure to the torsional response of the existing structure, comparison results between Case , Case , and Case are depicted in Figure 11(b). It can be observed that when the excitation direction was parallel to the -axis, the torsional responses of the roof floor of the TS in Cases and were generated, and if the excitation direction was perpendicular to the -axis, the torsional responses of the roof floor of TS in Cases and occurred. Thus, it can be concluded that the torsional response of TS can be influenced by the placed angle of the adjacent structure (AS).

Figure 12 shows the response of a structure adjacent to another structure with varying storey numbers. The first horizontal natural frequencies of the four-, six-, and eight-floor structures were 3.2, 2.3, and 1.7 Hz, respectively. And the first torsional natural frequencies of the four-, six-, and eight-floor structures were 4.5, 3.2, and 2.3 Hz, respectively. Symbol 4_6 was used to represent a 4-storey structure (TS) adjacent to a 6-storey structure (AS). Figure 13(a) shows the torsional response of the TS affected by the adjacent structure with varying storey numbers in Case ; it can be concluded that the torsional response of the roof of the structure can be affected by the adjacent structure, and when the targeted structure was adjacent to a structure with same structural and foundation configuration (Case (4_4)), the torsional dynamic response of the targeted structure may be more severe than the other cases (Case (4_6) and Case (4_8)). Meanwhile, through Figure 13(a), it can also be observed that around the horizontal natural frequency of the adjacent structure the torsional response of the targeted structure was generated. In order to investigate the torsional dynamic response of the targeted structure affected by the horizontal natural frequency of the adjacent structure, two cases shown in Figures 12(b)~12(c) are discussed here. As it has been described that the first horizontal natural frequencies of the four- and six-floor structures coincided with the first torsional natural frequencies of the six- and eight-floor structures. Figures 13(b)~13(c) show the torsional response of the structure affected by the adjacent structure; when the horizontal natural frequency of the adjacent structure coincided with the torsional natural frequency of the targeted structure, the torsional dynamic response of the targeted structure was more significant as can be seen in Figures 13(b) and 13(c).

**(a) Case 16: 4-floor structure (TS) adjacent to structures**

**(b) Case 17: 6-floor structure (TS) adjacent to structures**

**(c) Case 18: 8-floor structure (TS) adjacent to structures**

**(a) Case : excited perpendicular to direction**

**(b) Case : excited perpendicular to direction**

**(c) Case : excited perpendicular to direction**

#### 4. Conclusion

The effects of soil structure interaction on torsional response of structure supported on asymmetric foundation under horizontal seismic excitation were assessed in this study. The results of the study may lead to the following conclusions:(1)The torsional response can be generated for the structure supported on asymmetric foundation under horizontal seismic excitation due to the eccentricity in the soil foundation system.(2)For the structure supported by embedded-pile foundation with different pile lengths, the torsional response of the roof floor will be induced under horizontal seismic excitation as the mass center and the rigidity center of the soil foundation system with nonunique pile lengths do not coincide with each other, and the induced torsional response increases with the increase in the difference of pile lengths. Torsional response of the roof floor can also be observed for the structure supported by L shape or T shape embedded-pile foundation.(3)As the movement of the soil around the foundation of the targeted structure can be restricted by the presence of an adjacent structure, motion characteristics of the soil around the targeted structure at different spatial location may be impacted in different extent by the adjacent structure. Thus, the torsional response of the roof floor of the targeted structure can be observed under horizontal seismic excitation. Especially when the horizontal natural frequency of the adjacent structure coincides with the torsional natural frequency of the targeted structure, the torsional dynamic response of the targeted structure can be affected significantly.

#### Appendix

#### A. Formulation of Green’s Function Used in the Thin Layer Method and Expression of the Paraxial Boundary Condition

The displacement solution of a point located at surface for the load acting at surface (see Figure 14) includes(1)point load solution ();(2)disk load solution (, pile tip);(3)ring load solution (, pile shaft).

##### A.1. Point Load Solution

Let , , and displace components at elevation for a point load acting at elevation . Green’s function for point load is where is the second Hankel function of order . is the th eigenvalue of characteristic equation (A.7). is the th eigenvalue of characteristic equation (A.8). , are the corresponding eigenvector of . is the corresponding eigenvector of . is the th component of the eigenvector . is the th component of the eigenvector . is th component of the eigenvector .

##### A.2. Disk Load Solution

Let , , and displace components at elevation for a disk load acting at elevation . Green’s function for disk load is thenwhere is the radius of the disk.

##### A.3. Ring Load Solution

Consider the following:where is the radius of the ring.

##### A.4. Characteristic Equation of the Thin Layer Model as Shown in Figure 14

The characteristic equations of the thin layer model (shown in Figure 1) for Rayleigh (in-plane) and Love (anti-plane) waves are expressed by (A.7) and (A.8), respectively:where , , , , and are the matrices assembled by , , , , and , respectively, in the manner presented in where is the matrix corresponding to the th soil layer and is matrix corresponding to the boundary condition. Consider the following:where is the excited frequency of the load and , , , and are the Lame constant, shear elastic modulus, thickness, and density of the th soil layer, respectively.

and () represent the eigenvalue and eigenvector, respectively, of characteristic equation (A.7).

and () represent the eigenvalue and eigenvector, respectively, of the characteristic equation (A.8).

##### A.5. Paraxial Boundary

(1)The boundary condition of (A.8) is as follows:(2)The boundary condition of (A.7) is shown as follows: where , , , , , and are the Lame constant, shear elastic modulus, density, Poisson ratio, -wave velocity, and -wave velocity of the boundary, respectively.

#### Conflict of Interests

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

#### Acknowledgments

The authors would like to gratefully acknowledge the financial support for this research provided by the National Natural Science Foundation of China (Grant no. 51308243), the National Basic Research Program (973) of China (Grant no. 2011CB013800), and the Fundamental Research Funds for the Central Universities: HUST (Grant no. 2015TS122).

#### References

- E. Rosenblueth and R. Meli,
*The 1985 Earthquake: Causes and Effects in Mexico City*, vol. 5 of*Concrete International*, 1986. - S. L. Wood, R. Stark, and S. A. Greer, “Collapse of eight-story RC building during 1985 Chile earthquake,”
*Journal of Structural Engineering*, vol. 117, no. 2, pp. 600–619, 1991. View at: Publisher Site | Google Scholar - A. M. Chandler and G. L. Hutchinson, “Torsional coupling effects in the earthquake response of asymmetric buildings,”
*Engineering Structures*, vol. 8, no. 4, pp. 222–236, 1986. View at: Publisher Site | Google Scholar - J. C. de la Llera and A. K. Chopra, “Accidental torsion in buildings due to stiffness uncertainty,”
*Earthquake Engineering & Structural Dynamics*, vol. 23, no. 2, pp. 117–136, 1994. View at: Publisher Site | Google Scholar - J. C. de la Llera and A. K. Chopra, “Accidental torsion in buildings due to base rotational excitation,”
*Earthquake Engineering & Structural Dynamics*, vol. 23, no. 9, pp. 1003–1021, 1994. View at: Publisher Site | Google Scholar - J. C. de la Llera and A. K. Chopra, “Evaluation of code accidental-torsion provisions from building records,”
*Journal of Structural Engineering*, vol. 120, no. 2, pp. 597–616, 1994. View at: Publisher Site | Google Scholar - S. Nagarajaiah, A. M. Reinhorn, and M. C. Constantinou, “Torsion in base-isolated structures with elastomeric isolation systems,”
*Journal of Structural Engineering*, vol. 119, no. 10, pp. 2932–2951, 1993. View at: Publisher Site | Google Scholar - S. Nagarajaiah, A. M. Reinhorn, and M. C. Constantinou, “Torsional coupling in sliding base-isolated structures,”
*Journal of Structural Engineering*, vol. 119, no. 1, pp. 130–149, 1993. View at: Publisher Site | Google Scholar - R. S. Jangid and T. K. Datta, “Nonlinear response of torsionally coupled base isolated structure,”
*ASCE Journal of Structural Engineering*, vol. 120, no. 1, pp. 1–22, 1994. View at: Google Scholar - R. S. Jangid and T. K. Datta, “Seismic response of torsionally coupled structure with elasto-plastic base isolation,”
*Engineering Structures*, vol. 16, no. 4, pp. 256–262, 1994. View at: Publisher Site | Google Scholar - A. Tena-Colunga and C. Zambrana-Rojas, “Dynamic torsional amplifications of base-isolated structures with an eccentric isolation system,”
*Engineering Structures*, vol. 28, no. 1, pp. 72–83, 2006. View at: Publisher Site | Google Scholar - A. Tena-Colunga and L. Gómez-Soberón, “Torsional response of base-isolated structures due to asymmetries in the superstructure,”
*Engineering Structures*, vol. 24, no. 12, pp. 1587–1599, 2002. View at: Publisher Site | Google Scholar - A. Tena-Colunga and J. L. Escamilla-Cruz, “Torsional amplifications in asymmetric base-isolated structures,”
*Engineering Structures*, vol. 29, no. 2, pp. 237–247, 2007. View at: Publisher Site | Google Scholar - E. Çelebi and A. N. Gündüz, “An efficient seismic analysis procedure for torsionally coupled multistory buildings including soil–structure interaction,”
*Turkish Journal of Engineering and Environmental Sciences*, vol. 29, no. 3, pp. 143–157, 2005. View at: Google Scholar - M. Nasser, “Seismic response of R/C frames considering dynamic soil-structure interaction,” in
*Proceedings of the 18th International Conference on the Application of Computer Science and Mathematics in Architecture and Civil Engineering*, Weimar, Germany, July 2009. View at: Google Scholar - C.-C. Lin, J.-F. Wang, and C.-H. Tsai, “Dynamic parameter identification for irregular buildings considering soil-structure interaction effects,”
*Earthquake Spectra*, vol. 24, no. 3, pp. 641–666, 2008. View at: Publisher Site | Google Scholar - S. Mohasseb and B. Abdollahi, “Soil-structure interaction analyses using cone models,”
*Journal of Seismology and Earthquake Engineering*, vol. 10, no. 4, pp. 167–174, 2009. View at: Google Scholar - T. Balendra, C. W. Tat, and S. L. Lee, “Modal damping for torsionally coupled buildings on elastic foundation,”
*Earthquake Engineering & Structural Dynamics*, vol. 10, pp. 735–756, 1982. View at: Google Scholar - A. M. Chandler and G. L. Hutchinson, “Code design provisions for torsionally coupled buildings on elastic foundation,”
*Earthquake Engineering & Structural Dynamics*, vol. 15, no. 4, pp. 517–536, 1987. View at: Publisher Site | Google Scholar - A. M. Chandler and G. L. Hutchinson, “Parametric earthquake response of torsionally coupled buildings with foundation interaction,”
*Soil Dynamics and Earthquake Engineering*, vol. 6, no. 3, pp. 138–148, 1987. View at: Publisher Site | Google Scholar - K. S. Sivakumaran, M.-S. Lin, and P. Karasudhi, “Seismic analysis of asymmetric building-foundation systems,”
*Computers and Structures*, vol. 43, no. 6, pp. 1091–1103, 1992. View at: Publisher Site | Google Scholar - K. S. Sivakumaran and T. Balendra, “Seismic analysis of asymmetric multistorey buildings including foundation interaction and P-Δ effects,”
*Engineering Structures*, vol. 16, no. 8, pp. 609–625, 1994. View at: Publisher Site | Google Scholar - J. Lysmer, R. M. Tabatabaie, F. Tajirian, S. Vahdani, and F. Ostadan, “SASSI—a system for analysis of soil-structure interaction,”
*Report*UCB/GT/81-02, University of California, Berkeley, Calif, USA, 1981. View at: Google Scholar - E. Kausel, “An explicit solution for the green function for dynamic loads in layered media,” Tech. Rep. R81-13, Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, Mass, USA, 1981. View at: Google Scholar
- S. Hull,
*Dynamic loads in layered halfspaces [Ph.D. thesis]*, Massachusetts Institute of Technology, 1985. - S. Hull and E. Kausel, “Dynamic loads in layered halfspace,” in
*Proceedings of the 5th ASCE/EMD Specialty Conference*, pp. 201–204, Laramie, Wyo, USA, August 1984. View at: Google Scholar - Architectural Institute of Japan,
*An Introduction to Dynamic Soil-Structure Interaction*, Gihodoshppan, Tokyo, Japan, 1996.

#### Copyright

Copyright © 2016 Fangyuan Zhou 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.