Research Article | Open Access
Keunju Kim, Boo Hyun Nam, Heejung Youn, "Effect of Cyclic Loading on the Lateral Behavior of Offshore Monopiles Using the Strain Wedge Model", Mathematical Problems in Engineering, vol. 2015, Article ID 485319, 12 pages, 2015. https://doi.org/10.1155/2015/485319
Effect of Cyclic Loading on the Lateral Behavior of Offshore Monopiles Using the Strain Wedge Model
This paper presents the effect of cyclic loading on the lateral behavior of monopiles in terms of load-displacement curves, deflection curves, and - curves along the pile. A commercial software, Strain Wedge Model (SWM), was employed, simulating a 7.5 m in diameter and 60 m long steel monopile embedded into quartz sands. In order to account for the effect of cyclic loading, accumulated strains were calculated based on the results of drained cyclic triaxial compression tests, and the accumulated strains were combined with static strains representing input strains into the SWM. The input strains were estimated for different numbers of cycles ranging from 1 to 105 and 3 different cyclic lateral loads (25%, 50%, and 75% of static capacity). The lateral displacement at pile head was found to increase with increasing number of cycles and increasing cyclic lateral loads. In order to model these deformations resulting from cyclic loading, the initial stiffness of the - curves has to be significantly reduced.
Monopiles have been frequently used to resist offshore environmental loads such as wind, wave, tidal, and ice loads. Unlike onshore loading conditions, the offshore environmental loads have two main differentiated characteristics; the magnitude of lateral load is significant and the load is cyclic. Cyclic loading may cause serious problems in offshore structures and surrounding soils depending on the soil type and properties. During the design lifetime, the offshore structure is known to undergo about 108 lateral cyclic loading events of varying amplitude . During this great number of cycles, the offshore structures need to endure fatigue, and surrounding soils should not alter much to satisfy the lateral capacity and displacement criteria. Engineering properties of surrounding soils can deteriorate because of a possible accumulation of excess pore water pressure induced by cyclic loading. For cohesive soils, the deterioration of the clay structure or the forming of a gap between soil and pile could be the main concerns causing degradation of engineering properties. Furthermore, an excessive number of cycles may induce friction fatigue along the pile regardless of soil type, which is another critical factor causing the deterioration of side resistance. For cohesionless soils, the excess pore water pressure is unlikely to develop as the soils are normally in drained conditions; rather both the strength and stiffness will increase due to densification of surrounding soils. Still, the lateral displacement and rotation at pile head accumulate possibly leading to a serviceability problem. This is believed to be considerable especially in the case of a one way cyclic loading.
The - curves have been widely used in practice to analyze the lateral behavior of piles. Many - curves have been proposed chiefly considering soil types (sand or clay), loading types (static or cyclic), and ground water conditions (soil under water or above water). In conventional design process, the soils supporting piles are replaced by a set of nonlinear springs, and the pile is assumed to act as a beam supported by the springs. A soil spring acts as a reaction, , per lateral displacement, , at a pile elevation. The concept was first introduced by Winkler  and has been updated by many researchers. Reese et al.  suggested a parabola - curve and an empirical adjustment factor with different definitions for static and cyclic loading. Later, simplified approaches to establish the - curves have been reported .
As the aforementioned - curves were developed not for offshore wind turbines but for offshore oil platforms, there have been concerns that the offshore monopiles with their large lateral load compared to dead weight may behave differently. Moreover, the available - curves are validated for small piles with diameter of 1-2 m under small number of cycles only, while offshore monopiles have 5–7 m diameter and are subject to much more cycles. The existing - curves may not be applicable to the offshore monopiles without appropriate validation. Achmus et al.  instead developed design charts for offshore wind turbine foundation installed in sandy soils based on a stiffness degradation model. Cyclic triaxial tests and series of numerical analyses were carried out in order to investigate the long-term behavior under cyclic loading. Wichtmann et al. [5, 6] suggested a High Cycle Accumulation (HCA) model, which captures the effect of cyclic loading, and performed series of finite element analyses. The proposed HCA model was based on series of drained cyclic triaxial compression and extension tests [7, 8]. Instead of using complicated FEM, practical software analysing the lateral behavior of a pile by means of the Strain Wedge Model (SWM) was applied in  to quantify the deformation of an offshore wind turbine, illustrating an increase in the lateral displacement with increasing number of cycles. Recently, nonlinear - curves in cohesive soils were proposed accounting for soil-pile interaction, degradation in soil stiffness and strength due to cyclic loading, soil-pile gap forming, and radiation damping .
The numerical analyses were often supported by small or full scale test results. Roesen et al. [11, 12] investigated the soil-pile interaction of a monopile embedded in sands using small scale tests. The monopile was subject to a high number of one way cyclic loads (up to 60,000 cycles) with varying amplitude, mean load, and loading period. It was reported that the increase of the accumulated rotation with growing number of cycles could be approximated by a power function. Moller and Christiansen  conducted small scale tests to evaluate the existing - methods for static and cyclic loading and compared them with numerical analysis results.
This paper investigates the lateral behavior of a monopile subjected to a lateral cyclic loading, focusing particularly on the influences of the number of cycles and magnitude of cyclic lateral load. The effect of cyclic loading was assessed in terms of accumulated strains, which were later modified into the input strains used in the SWM . The load-displacement curves, deflection curves, and cyclic - curves were obtained from SWM; and the lateral behavior of the monopile was analyzed.
2. Material Properties and Loading Conditions
Figure 1 illustrates the offshore monopile embedded in quartz sands used in this study along with the pile and soil properties. The monopile, which has to resist relatively large lateral load, was a 60 m long and 90 mm thick steel tube, the diameter of which was 7.5 m. The pile was assumed to be linear elastic with Young’s modulus of 210 GPa; the soils were assumed to be homogeneous through the layer. The material properties of the soil were adopted from a previous research, which investigated the effect of cyclic loading on the accumulated displacement of pile embedded into quartz sand . The investigated sands have a peak internal friction angle of 37.5°, effective unit weight of 11 kN/m3, and a cohesion of 0 kPa. The used soil was dense; thus, the relative density of the soil was assumed to be 75% throughout the layer, and the maximum and minimum void ratio were assumed to be 0.874 and 0.577 as used in the literature .
The information of cyclic lateral load is difficult to quantify because the load is site specific and variable in direction and magnitude. Therefore, the static lateral capacity of the monopile was calculated using the load-displacement curve at seabed from the SWM analysis, and the cyclic lateral loads were chosen as 25%, 50%, and 75% of the static capacity. The static lateral capacity was determined to be 22.04 MN at 38.1 mm of lateral displacement at seabed . Consequently, the used cyclic lateral loads, , were 5.51, 11.02, and 16.53 MN, respectively. In the simulations, after applying a certain number of cyclic loads, the pile head was loaded with the lateral load of 22.04 MN to obtain the load-displacement relationship at pile head, and the curves were obtained at 5, 15, 25, and 35 m depth along the pile below seabed.
3. Input Parameter Calculation
3.1. SWM Analysis
The Strain Wedge Model was initially suggested by Norris  and has been used to predict the behavior of flexible piles under lateral loading. Figure 2 presents the configuration of the Strain Wedge Model in uniform soil. The soil resistance against lateral loading is developed by a 3-dimensional passive wedge of soil at the front of the pile, and the configuration of the passive wedge is determined by the mobilized friction angle and pile diameter: , is the pile diameter, is the height of the passive wedge, is the variation in horizontal stress at wedge face, and is the shear stress at pile side. The advantage of the SWM is the transfer of the stress-strain-strength behavior of soil observed in element tests to the 1-dimensional Beam on Elastic Foundation (BEF) parameters. While the conventional - curves were experimentally derived from strains developed along the laterally loaded pile, the SWM provides a theoretical link between soil properties and lateral behavior of pile. The horizontal strain () in the passive wedge is used to determine the deflection of the pile , and the variation of horizontal stress at the passive wedge is used to determine the soil resistance associated with the BEF. Details of the theory and the assumptions used in the SWM are provided, for example, in .
Figure 3 shows a flow chart illustrating the numerical process proposed in this paper to analyze the monopile under cyclic lateral loading. The key input parameters into the SWM include the static strain at 50% peak strength (), dimension and Young’s modulus of the steel pile, and the material properties of the surrounding soils. The SWM requires the strain at 50% peak strength () in order to assess the lateral behavior of the pile under the static loading condition. On the other hand, the strain under the cyclic loading condition accumulates as the number of cycles increases (), which means the strain will not recover to the original state of the pile upon unloading . Therefore, the strain under cyclic loading needs to be modified into an input strain () considering key factors such as the magnitude of cyclic lateral loads, average mean stress, void ratio, and number of cyclic loads. The proposed methodology to assess the key factors will be presented in later sections.
3.2. Input Strain under Cyclic Loading
Figure 4 shows the concept of calculating the accumulated strain developed by the applied number of cycles in the plane. At the end of cycles, the plastic strain has accumulated to some extent, which is denoted by . During a subsequent monotonic loading, the deviatoric stress of the soil would start from the accumulated strain at zero deviatoric stress to peak strength (as shown by line ①). The probable strain at failure after the considered number of cycles was calculated by summing twice the accumulated strain () and twice the static strain at 50% peak strength (). Then the strain at failure under cyclic loading was modified into input strain at 50% peak strength (), which means the corresponding line for a monotonic loading incorporating the cumulative effects follows line ②. Even though this approach does not perfectly capture the real behavior of soils and the real accumulated strain, the calculation procedure seems to be simple and reasonable. The input strains used in this study were calculated by
The static strain () varies with the void ratio (i.e., relative density) and coefficient of uniformity of the sand as shown in Figure 5. With the chosen relative density (75%) and coefficient of uniformity (1.8) from previous research , the static strain () of the studied soil was calculated to be 0.42% at the reference confining pressure of 42.5 kPa. The static strain needs modification to account for the effect of confining pressure, which varies along the pile length. Thus, the static strain under the reference confining pressure (42.5 kPa) is modified into an input value for the effective overburden pressure at desired depth using (2) proposed in . Consider the following:
3.3. Accumulated Strain
The accumulated strains as a result of a large number of cycles were calculated based on test data from cyclic drained triaxial tests in compression and could be assessed using the following proposed equation : where is a coefficient for strain amplitude, is a coefficient for average mean stress (), is a coefficient for average stress ratio (), is a coefficient for void ratio, is a coefficient for number of cycles, and is a factor associated with polarization change. Equation (3) was originally developed for the norm of the strain tensor, but within this paper the assumption is made for simplification that the total strain predicted by (3) equals the horizontal strain in the soil at the front of the pile. Table 1 summarizes the partial functions for the coefficients and the suggested material constants used in this study.
Figure 6 illustrates the stress path in the - plane applied in a cyclic triaxial test in compression. When cyclic lateral load is applied to the pile head, the soils at the pile front are compressed and released repeatedly, which can be simulated by cyclic triaxial test in compression. In fact, the SWM recommends using triaxial compression test results ; thus the results from cyclic triaxial test in compression were adopted in this study. Increase in axial (horizontal) stress () causes an increasing deviatoric stress from the initial in situ stress on the – line to the desired deviatoric stress amplitude above the – line, and then the axial stress decreased causing a decrease in deviatoric stress to the desired amplitude under the – line. In one way loading condition, the cyclic lateral load () induces the deviatoric stress amplitude () above and under the – line, indicating that the increase in the axial stress () is twice the deviatoric stress amplitude. It should be noted that Ashour et al.  recommend using conventional triaxial compression test where the confining pressure is isotropic. However, the cyclic triaxial test has anisotropic confining pressure before applying cyclic axial stress. The anisotropic confining pressure as shown in the - plane is appropriate for in situ stress condition, and the cyclic lateral load induces axial stress always greater than zero which occurs in one way cyclic loading condition.
The coefficient for the amplitude () could be calculated from strain amplitude () which decreased with the number of cycles in the cyclic triaxial tests with constant stress amplitude. Figure 7 shows the development of strain amplitude with the number of cycles for various deviatoric stress amplitudes at an average stress with kPa and . In order to calculate deviatoric stress amplitude along the pile resulting from cyclic lateral load, SWM analyses were conducted using static strains. From the static analyses, the - curves and pile deflection curves were obtained. The - curves provide the soil resistance per unit length at a specific lateral displacement , and the lateral displacement can be obtained from pile deflection curves. The estimated soil resistance per unit length was divided by the pile diameter (7.5 m), resulting in the increase in cyclic axial stress . The deviatoric stress amplitude is half the cyclic axial stress. In the same way, the deviatoric stress amplitudes along the pile at 5, 15, 25, and 35 m depth were obtained.
Using the obtained deviatoric stress amplitudes along the pile, the strain amplitudes were estimated using Figure 7. Linear interpolation technique was used when exact deviatoric stress amplitude is unavailable. As the strain amplitude was measured with an average mean pressure of 200 kPa, the strain amplitude needs correction to account for the effect of pressure dependency of the secant stiffness. With simple assumption that Young’s modulus is proportional to the square root of average mean stress (4), the strain amplitudes at different mean stresses (different depth) were calculated. Using the correction factor, , provided in (5), the strain amplitudes for every mean stress were calculated as shown in (6). Applying the correction factors to the strain amplitudes, the pressure dependency of strain amplitudes was considered along the pile.
Table 2 provides the strain amplitudes for different depths and number of cycles (1 and 105) under 3 cyclic lateral loads. The cyclic lateral loads were selected at 25%, 50%, and 75% of the static capacity. As shown in the table, the strain amplitudes were found to be larger for higher cyclic lateral loads, implying that the largest cyclic lateral load (75% of the static capacity) would lead to the largest accumulated strains, and consequently to the largest lateral displacement:
The coefficient for the number of cycles () was simply calculated by inserting the number of cycles under consideration together with the suggested material constants into the partial function provided in Table 1. It was reported that the accumulated strain increased with decreasing average mean stress, and the coefficient for average mean stress () was slightly affected by the number of cycles . However, the recommended material constants neglected the effect of the number of cycles and were adopted as such in this study. The average stress ratio (the ratio of the average deviatoric stress to the average mean stress) as well as the void ratio was assumed to be constant; accordingly, the coefficients, and , were calculated to be 1.53 and 0.13. The coefficient was 1.0 since polarization changes were not considered in this study.
Tables 3 and 4 present the calculated coefficients for the influences of the strain amplitude, number of cycles, average stress ratio, void ratio, polarization, and average mean stress. The homogeneous subground was divided into 10 m thick layers, and the coefficients were calculated for each layer. It is shown that the coefficients for strain amplitude significantly vary with depths but only slightly with the number of cycles and that the coefficients for the number of cycles vary from 1.49 () to 57.51 () for 1 and 105 cycles, respectively. Table 5 tabulates the calculated accumulated strains and input strains at and 105 accounting for all the coefficients provided in Tables 3 and 4. It is clearly shown that most accumulated strains develop at a depth shallower than 25 m regardless of the number of cycles. This agrees well with the fact that the strain amplitudes are very small at the greater depths for the given cyclic lateral loads. The magnitude of input strains is found to be larger for larger cyclic lateral loads and larger numbers of cycles.
4. Numerical Results and Discussion
4.1. SWM Analysis versus FEM
Achmus et al.  performed the Finite Element Method (FEM) with the stiffness degradation model to investigate the accumulated displacement due to cyclic loading. In order to check the applicability of the SWM analysis for cyclic lateral loading, the lateral displacement of pile at seabed level was compared with the FEM results (Figure 8). The used soil and pile properties are identical to those in Figure 1, and the used cyclic lateral load was 15.0 MN. The number of cycles ranges from 1 to 104. It is shown that the displacement at seabed is similar for both analyses when , but the SWM analysis tends to underestimate the displacement when the number of cycles gets higher. For = 104, the SWM analysis predicts 46% increase in lateral displacement whereas the FEM predicts 80% increase. Both FEM and SWM analyses predict the increase in accumulated lateral displacement with increasing number of cycles but underestimate the displacements compared to the predictions proposed by other researchers [18, 19].
4.2. Load-Displacement Curve
Figure 9 shows the load-displacement curves at the pile head for different cyclic lateral loads and numbers of cycles. For the given lateral load of 22.04 MN, the pile was found to be more displaced as the number of cycles increased. When cyclic lateral load was assumed to be 25% of the static capacity, the displacement for 1 cycle of lateral load was 183.11 mm but increased up to 201.93 mm for 105 cycles. The 105 cycles are likely to induce approximately 111% additional displacement for the same magnitude of lateral load (22.04 MN). The discrepancy between the displacements at and 105 appeared to significantly grow when the largest magnitude in cyclic lateral load was considered (75% of the static capacity). With the largest cyclic load, the lateral displacement for 105 cycles was calculated to be approximately 140% of the displacement under static loading. It was found that an increase both in the magnitude of the cyclic lateral load and in the number of cycles significantly increases the accumulated lateral displacement.
(a) 25% of
(b) 50% of
(c) 75% of
Figure 10 summarizes the effect of the number of cycles and the cyclic lateral load on the lateral displacement at pile head. The lateral displacement appeared to increase exponentially with increasing logarithm of the number of cycles, leading to about 40% increase of the displacement when comparing with cycles. This implies that the lateral displacement properly predicted for static lateral load would be inappropriate for cyclic loading conditions. It should be noted that the partial safety factor for displacement is usually 1.0 in offshore pile design.
4.3. Pile Deflection Curve
Figure 11 presents deflection curves along the pile. Similar to the load-displacement curves, the pile deflection increased with increasing number of cycles and growing cyclic lateral load. Even though the pile does not deflect significantly below 30 m depth, the deflection near the ground surface considerably differs depending on loading conditions. For the smallest cyclic lateral load, the pile deflections calculated for the various numbers of cycles did not significantly differ; however, the deflection was clearly dependent on the number of cycles for the biggest lateral load.
(a) 25% of
(b) 50% of
(c) 75% of
4.4. Cyclic - Curves
Figures 12-13 present - curves at 2 and 4 m depth considering the effect of cyclic loading, indicating that the consideration of the cumulative deformations due to cyclic loading leads to a considerable decrease of the initial stiffness of the - curves. Even though SWM provides - curves at any desired locations, two shallow depths were particularly chosen because soil properties near the ground surface play a key role in the lateral behavior of the monopile. As shown in the figures, the initial stiffness of the - curves sharply decreased with increasing number of cycles and increasing cyclic lateral load. The stiffness decrease becomes more prominent with increasing cyclic lateral load. Such trends in the - curves do not differ with depth.
(a) 25% of
(b) 50% of
(c) 75% of
(a) 25% of
(b) 50% of
(c) 75% of
In this study, the lateral behavior of a monopile for an offshore wind turbine subject to cyclic loading was investigated using the SWM analysis in combination with the data from drained cyclic triaxial compression test results on quartz sand. From the SWM analyses, the load-displacement curves, deflection curves, and - curves were obtained and compared with each other. The following conclusions were drawn:(1)A new approach was proposed to investigate the effect of cyclic lateral loads on the lateral behavior of monopiles. The approach employed the input strains (), which were calculated based on the accumulated strains caused by cyclic lateral loads () and the static strains ().(2)The displacement at pile head, representing the sum of accumulated displacements due to a cyclic loading and a subsequent monotonic loading towards a maximum load larger than the amplitude of the cycles, was found to increase exponentially as the logarithm of the number of cycles increased. For a given lateral load of 22.04 MN (lateral static capacity of the pile under consideration), the increase in lateral displacement at pile head was 111% if 105 cycles with an amplitude corresponding to 25% of the static capacity had been previously applied. The increase of the lateral displacement due to preceding cycles becomes as much as 140% at an amplitude corresponding to 75% of the static capacity.(3)The - curves obtained at 2 and 4 m depths indicated that the initial stiffness of the - curves considerably decreased with increasing number of cycles. The decrease became prominent when the cyclic lateral load was 75% of the static capacity. However, it should be kept in mind that this decrease of stiffness is some kind of “calculation” in order to consider the accumulated deformations. The real stiffness of the sand surrounding the pile usually increases due to compaction.(4)The application of existing - curves without considering cyclic loading would lead to an underestimation of the lateral displacement at pile head. Modified - curves accounting for the magnitude of cyclic loading and number of cycles are required to be used in the BEF analysis.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was partially supported by National Research Foundation of Korea (NRF) funded by Ministry of Science, ICT & Future Planning (NRF-2013R1A1A1011983), and Korea Institute of Ocean Science and Technology (20113010020010-11-3-510).
- M. Achmus, Y.-S. Kuo, and K. Abdel-Rahman, “Behavior of monopile foundations under cyclic lateral load,” Computers and Geotechnics, vol. 36, no. 5, pp. 725–735, 2009.
- E. Winkler, Theory of Elasticity and Strength, Domimicus, Prague, Czech Republic, 1867.
- L. C. Reese, W. R. Cox, and F. D. Koop, “Analysis of laterally loaded piles in sand,” in Proceedings of the Offshore Technology in Civil Engineering: sHall of Fame Papers from the Early Years, pp. 95–105, Houston, Tex, USA, May 1974.
- D. Bogard and H. Matlock, “Simplified calculation of p-y curves for laterally loaded piles in sand,” Tech. Rep., The Earth Technology, Houston, Tex, USA, 1980.
- T. Wichtmann, A. Niemunis, and T. Triantafyllidis, “Prediction of long-term deformations for monopile foundations of offshore wind power plants,” in Proceedings of the 11th Baltic Sea Geotechnical Conference: Geotechnics in Maritime Engineering, Gdańsk, Poland, September 2008.
- T. Wichtmann, A. Niemunis, and T. Triantafyllidis, “Towards the FE prediction of permanent deformations of offshore wind power plant foundations using a high-cycle accumulation model,” in Proceedings of the 2nd International Symposium on Frontiers in Offshore Geotechnics (ISFOG '10), pp. 635–640, November 2010.
- T. Wichtmann, A. Niemunis, and T. Triantafyllidis, “Strain accumulation in sand due to cyclic loading: drained triaxial tests,” Soil Dynamics and Earthquake Engineering, vol. 25, no. 12, pp. 967–979, 2005.
- T. Wichtmann, A. Niemunis, and T. Triantafyllidis, “Strain accumulation in sand due to cyclic loading: drained cyclic tests with triaxial extension,” Soil Dynamics and Earthquake Engineering, vol. 27, no. 1, pp. 42–48, 2007.
- P. Hinz, K. Lesny, and W. Richwien, “Prediction of monopile deformation under high cyclic lateral loading,” in Proceedings of the 8th German Wind Energy Conference, DEWEK, Bremen, Germany, 2006.
- M. Heidari, M. Jahanandish, H. E. Naggar, and A. Ghahramani, “Nonlinear cyclic behavior of laterally loaded pile in cohesive soil,” Canadian Geotechnical Journal, vol. 51, no. 2, pp. 129–143, 2014.
- H. R. Roesen, L. B. Ibsen, and L. V. Andersen, “Experimental testing of monopiles in sand subjected to one-way long-term cyclic lateral loading,” in Proceedings of the 18th International Conference on Soil Mechanics and Geotechnical Engineering, 2013.
- H. R. Roesen, L. B. Ibsen, M. Hansen, T. K. Wolf, and K. L. Rasmussen, “Laboratory testing of cyclic laterally loaded pile in cohesionless soil,” in Proceedings of the 23rd International Offshore and Polar Engineering Conference (ISOPE '13), pp. 594–601, July 2013.
- I. F. Moller and T. H. Christiansen, Laterally loaded monopile in dry and saturated sand—static and cyclic loading [M.S. thesis], Aalborg University, Aalborg, Denmark, 2011.
- L. C. C. Geopile, “SWM 6.2 version 6.2 : Laterally and axially loaded single pile and pile group,” 2000.
- AASHTO, AASHTO LRFD Bridge Design Specifications, American Association of Highway and Transportation Officials, Washington, DC, USA, 2007.
- G. Norris, “Theoretically based BEF laterally loaded pile analysis,” in Proceedings of the 3rd International Conference on Numerical Methods in Offshore Piling, pp. 361–386, Nantes, France, 1986.
- M. Ashour, G. Norris, and P. Pilling, “Lateral loading of a pile in layered soil using the strain wedge model,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 124, no. 4, pp. 303–314, 1998.
- R. L. Little and J. L. Briaud, “Full scale cyclic lateral load tests on six single piles in sand,” Miscellaneous Paper GL-88-27, Geotechnical Division, Texas A&M Unuiversity, College Station, Tex, USA, 1988.
- A. Hettler, Verschiebungen Starrer und Elastischer Gründungskörper in Sand bei Monotoner und Zyklischer Belastung, Institut für Bodenmechanik und Felsmechanik der Universität Fridericiana, Karlsruhe, Germany, 1981.
Copyright © 2015 Keunju Kim 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.