**Elements for numerical simulation of wind time series**

**I. Fernández ^{1}*, V. Elena Parnás ***

* Instituto Superior Politécnico José Antonio Echeverría (CUJAE). CUBA

**ABSTRACT**

Dynamic characteristics of structures determine the methods to be used for their analysis and design against environmental loads such as wind. In some structures, the effect of the energy contained in gusts sequences must be attended, emphasizing the fluctuating components of the gusts that resonate with them. In these cases, dynamic methods in time domain are frequently used. Numerical simulation of wind loads is a tool that has been widely developed and applied. The aim of this study is to perform a literature review to determine the necessary parameters for simulating wind time series as well as the existing numerical methods for this purpose, emphasizing in tropical cyclone-prone region.

**Keywords:** Wind velocity, mean velocity, turbulence, numerical simulation, time history series

**1. Introduction**

The wind is represented by its movement in two different time scales: 1) movements related to the variations of the planetary weather systems and 2) gusts. Van der Hoven (Van der Hoven, 1957) explained this behavior by identifying two phases within the wind power spectrum. The first phase is associated to the first and highest peak within the spectrum, which he calls macro or mesometeorological peak, which shows an energy concentration period of approximately four days corresponding to the typical transit time of a fully developed weather system. The second important peak, called micrometeorological, contains energies with periodic fluctuations around 1 minute up to tenths of seconds, and it is caused by turbulences caused by topographic effects: land roughness or obstacles around the site. The frequency range between the two peaks lacks energy, and it is known as a spectral gap. The greatest relevance of this gap is that enables to represent the wind in terms of a mean speed depending on the height above the ground level, which reflects the synoptic variations only, with the superposition of time-dependent fluctuations in three directions, the so-called gusts.

The type of analyzed structure also determines how that load is going to be assigned on the structure, either as a static force or as a time-dependent dynamic force, given the fluctuating nature of the action and, consequently, of the response. Some structures respond much less to the intensity of an individual gust than to the energy contained in gust sequences, particularly regarding the fluctuating components of the gusts that resonate with them. Therefore, dynamic methods in the time or frequency domains are recommended for the analysis and design of these structural types; and for those having also a strong non-linear behavior, the time domain is more reliable.

The numerical simulation of synthetic time series is a tool that has been widely developed and applied in the generation of wind time functions that can subsequently be applied as load in the structures. Two wind components, mean and fluctuation, have to be studied for the simulation, and their elements have to be very carefully defined to accurately reflect the large and small-scale movements of the wind. The objective of this paper is to define a theoretical framework to determine the simulation parameters and methods of wind speed time series, differentiated according to the two components, with emphasis on hurricane-prone regions.

**1.1 Mean wind speed. Vertical profiles**

The wind’s mean speed component is mainly conditioned by temperature and frictional interactions with the ground, due to the roughness and orography, which try to slow down the airflow by forming a vertical profile.

Several analytic models can describe the vertical profile of speeds and the three most known and used are: the logarithm law, the power law, and the Deaves and Harris model. According to Tamura and Kareem (Y. Tamura and Kareem, 2013), the log law is the most adequate for strong wind systems; however, it offers some mathematical complications that engineers, for practical reasons, prefer to avoid by using the power law. Studies by Yaojun et al (Yaojun Ge, Xinyang Jin and Cao, 2010) compare the Chinese standard with a total of 15 codes and design standards addressing wind loads in several countries of the region, referred by them as “Asian Pacific” economies. They reaffirm that 87% of the studied codes prefer and use the power law (Canada, China, Hong Kong, India, Indonesia, Japan, Korea, Malaysia, Philippines, Taiwan, Thailand, United States and Vietnam).

The adjustment of the behavior of the mean wind to previous profiles in tropical cyclone-prone regions has been studied. Among the main questionings of the scientific community is the comparison of profiles close to the cyclone eyewall with those of the outer regions of the center vortex of the eye. In 2009, Vickery et al. (Vickery, Wadhera, Powell and Chen, 2009) studied this radial dependence and determined that, in lower areas close to the ground level (above 100 m), the variation of the wind’s mean speed is approximately logarithmic, and as the height increases, the applicability of the log law is broken and wind velocities start to reduce; a similar behavior is observed in the inner and outer areas of the cyclonic vortex. The authors propose a new empirical model, thereby recognizing that it properly adjusts to the sea surface, but they recommend extending the studies in order to apply them inland. In 2012 and 2013, Giammanco et al. (Ian M. Giammanco, Schröder and Powell, 2012; Ian M. Giammanco and Schröder, 2013) obtained results that are quite adjustable to those of Vickery et al. (Vickery et al., 2009). Cao et al. (Yukio Tamura, Cao and Giang, 2012) studied the vertical profiles during the path of the Maemi typhoon across Japan, concluding that the height variation of the mean wind speed in the regions close to the surface can be obtained through the classic log law or power law. They also confirmed that no significate differences were observed in the profile’s characteristics on the inner or outer area of the eye of the typhoon. Based on researches of Vickery et al. (Vickery et al., 2009) and Giammanco et al. (Ian M. Giammanco et al., 2012; Ian M. Giammanco and Schröder, 2013), Tse et al. (Tse, Li, Chan, Mok and Weerasuriya, 2013) use the power law, the log law and the empirical profile obtained by Vickery et al. (Vickery et al., 2009) to adjust the mean profiles of typhoons observed from 2007 to 2009 in Hong Kong, at averages of 1 h and 10 min. According to their results, they do not recommend using the empirical law of Vickery et al. (Vickery et al., 2009) for 10-min averages and heights over 300 m; on the contrary, both the power law and the log law gave adequate results for these same conditions.

**1.2 Wind speed fluctuations**

The atmospheric turbulence , is a stochastic process involving three variables (*3V*) (*u, v, w*) and four dimensions (*4D*) (*x,y ,z; t*). It is usually assumed as a stationary random field, with a media of zero, with Gaussian behavior and time-dependent; this review accepts this behavior. It can be assumed that stochastic components (*u,v,w*) are statistically independent from each other, changing from a 3V-4D process towards a 3(1V-4D). Through the discretization of the space domain in N points, sites where the wind has an impact on the structure, each component of the turbulence (*u,v,w*) can be considered a stochastic vector. Each component vector is a stochastic field (1V-1D) (representing one of the components of the wind speed fluctuation in a specific location, dependent on the deterministic parameter t, which is correlated to other components of the same vector). Thereby, it is possible to change from a 1V-4D process to three NV-1D process, each one corresponding to one of the space components of the turbulence. This method conduces to a simplified model, which disregards any link with another point close to the one considered; rather than taking into account the space distribution of the wind field, from a full probabilistic characterization of the stochastic field reached by the complex form of the cross power spectral density (CPSD). In general, to determine the dynamic longitudinal response of the action of the wind on the structures, the components that are orthogonal to the direction of the fluid (x) may be disregarded; thus, the CPSD is expressed for two points located on the plane perpendicular to the action of the wind (*y-z*) through Equation 1, assuming that the imaginary part of the cross spectrum is insignificant (Ubertini and Giuliano, 2010).

**(1)**

Where (-f_{12}(ω)), is the coherence function, which quantifies the cross correlation of the analog components of the turbulence at different points of the stochastic field.

One of the alternative formules for f_{12}(ω) can be obtained by equation 2 (Cottone and Di Paola, 2011; Guoqing Huang, Liao and Li, 2013; Koulatsou, Petrini, Vernardos and Gantes, 2013).

**(2)**

Where is the mean speed and k_{y} y k_{z} is the mean speed and k_{y} and k_{z} are experimentally determined decay coefficients (Di Paola and Gullo, 2001). The latter differ for each turbulence component and depend on various factors such as height, surface roughness, atmospheric stability, among others. According to Chen et al. (Chen, Li and Xiang, 2014), factors normally range from 7 to 10, and according to Xu et al, they are between 6 and 26 (Xu et al., 2014). Solari and Piccardo (G. Solari and Piccardo, 2001) analyzed 39 reference sources to study the variability of decay coefficients, which allowed creating the summary Table 1 for different directions of the wind action and the turbulence speed components.

**Table 1.** Coefficients proposed by Solari and Piccardo (G. Solari and Piccardo, 2001). x, y, z: longitudinal, lateral and vertical directions of the wind action

With regard to these coefficients in hurricane-prone regions, there are some works based on field measurements. The study of Tamura et al. (Yukio Tamura et al., 2012) is one of them, which shows values that do not differ greatly from those proposed by Solari and Piccardo (G. Solari and Piccardo, 2001)

In Equation 1, S_{u1}(ω) and S_{u2}(ω) are spectral density functions determined at each point. Meteorology and wind energy use different mathematical formulas for them.

In 1948, as a result of a large experimentation developed in wind tunnels with isotropic and homogenous turbulent airflow, von Kármán (von Kármán, 1948) proposed the first recognized formula for S_{u}(ω), (Equation 3), specifically formulated for the longitudinal component of fluctuations.

**(3)**

Where z the height measured from the ground to the point, L_{u} is the integral scale of the turbulence component, σ_{u2} is the variance of the component, and u^{2} is the square of the friction velocity.

Based on the expression of von Kármán, other ones have derived. The Davenport spectrum (Davenport, 1961), the Kaimal-type spectrum (Kaimal, Wyngaard, Izumi and Coté, 1972) and the spectrum proposed by Solari and Piccardo (G. Solari and Piccardo, 2001) are examples thereof. Among them, the spectra of Solari and Piccardo (Ciampoli, Petrini and Augusti, 2011; Cottone and Di Paola, 2011; Ubertini and Giuliano, 2010), Kaimal (Hernandez, Bernal and Caracoglia, 2013; J. Li, Li and Chen, 2011; Pinheiro, 2010) and von Kármán (An, Quan and Gu, 2012; Fu, Wu, Xu, Li and Xiao, 2012; Q. S. Li, Lunhai Zhi and Fei Hu, 2010) are the most commonly used by the international scientific community. The standard of Australia and New Zealand (AS/NZS1170.2-2011, 2011) and ISO (International Organization for Standardization, 2009) use the von Kármán spectrum. And according to Tamura and Kareem (Y. Tamura and Kareem, 2013), and endorsed by a group of works which monitored structures subjected to the action of hurricanes (Fu et al., 2012; Q. S. Li et al., 2010; Yukio Tamura et al., 2012), this is the most adequate spectrum for describing the wind turbulence during this type of storms.

**1.3 Turbulence parameters associated to the spectrum**

Certain parameters derived from the spectrum in equation 3 deserve special attention, due to a bibliographic divergence regarding their determination. These parameters are mainly: the integral scale and the turbulence intensity, defining the terms variance, friction velocity and turbulence intensity factor for the latter.

The turbulence integral length scale, L_{u}, defines the position of the spectral content of the turbulence. It is a measure of the average size of vortices in the fluid inside the atmospheric boundary layer. It depends essentially on the roughness length, Z_{o}, the height from the ground z, the mean speed, and the parameters of Coriolis, which depend on the latitude of the observation site (WMO, 2010). Z_{o} is the height in relation to the ground, where the wind speed is theoretically zero when the profile is assumed as logarithmic. There are two types of approaches to determine this parameter inside the atmospheric boundary layer (L. Li, Kareem, Xiao, Song and Zhou, 2015; Q. S. Li et al., 2010): 1) the morphometric methods, which use the characteristics of the surface to determine the roughness and 2) the micrometeorological methods. The work of Li et al. (L. Li et al., 2015) reviews some of the main works dealing with each of these methods.

There are several ways of obtaining L_{u} that can be consulted in the references (L. Li, Xiao, Kareem, Song and Qin, 2012; Wang, Hu and Cheng, 2011). Solari and Piccardo (G. Solari and Piccardo, 2001), based on the study of approximately 20 bibliographical sources, propose equation 4 to determine the turbulence integral scale. Equation 4 has been used by several authors who recognize its validity (Carassale and Solari, 2006; Cottone and Di Paola, 2011; Torrielli, 2011).

**(4)**

Some works obtain L_{u} a partir de mediciones de campo en presencia de tifones (Shiau, 2000; Tamura et al., 2012; Wang et al., 2011). They report that L_{u} increase as the mean speed increases, being higher than the integral scales for the rest of the turbulence components.

According to a study by Li et al. (L. Li et al., 2005), L_{u} are higher in typhoons than in hurricanes, and they associate this to the fact that basins and latitudes originating the typhoons and hurricanes are a potential source of variations. The mentioned paper makes a comparison between the experimentally measured values and the results derived from applying the formulation proposed by the US standard ASCE7-10, 2010, given by Equation 5.

**(5)**

Where *l*, is equal to 97.54, 152.4 and 198.12 and equal to 1⁄3, 1⁄5 and 1⁄8 for exposure categories B, C and D respectively. They conclude that the formulation of ASCE presents good approximations for speeds lower than 25 m/s and underestimates the value according to the real measurements when the speed is higher than 30 m/s.

The Japanese code AIJ-RLB of 2004 recommends the use of Equation 6 to determine the integral scale (Fu et al., 2012; Q. S. Li et al., 2010).

**(6)**

According to the work of Li et al.(Q.S. Li, Lunhai Zhi and Fei Hu, 2010), which compares the results of finding L_{u} in field measurements of an instrumented structure in the city of Beijing and the empirical formulations proposed by AIJ-RLB (6) and ASCE (5), it is valid to apply the formulation by AIJ-RLB, due to its good adjustments. An et al. (An et al., 2012) and Fu et al. (Fu et al., 2012) differ from this criterion by concluding that the expression of the Japanese code does not offer good approximations of the turbulence scale, compared with the measurements made during the path of the Muifa and Megi typhoons across the Chinese cities of Shanghai and Guangzhou, respectively.

Table 2 summarizes some of the formulations proposed by four of the main wind codes, in addition to the Japanese and the US codes.

The variance of the longitudinal turbulence component σ_{u}2 = β_{u*}2 where β is a non-dimensional coefficient called turbulence intensity factor, depends on the same parameters as the integral scale. In order to obtain it, the terms turbulence intensity factor and friction velocity must be specified.

According to a large number of studies, which Solari and Piccardo (G. Solari and Piccardo, 2001) quote from Bietry, Sacré and Simiu of 1978, the former propose equation 7 to determine the value of the turbulence intensity factor in the longitudinal direction of the wind, which has been adopted as valid in several published researches (Ciampoli et al., 2011; Cottone and Di Paola, 2011; Ubertini and Giuliano, 2010).

**(7)**

During the path of tropical cyclones, the turbulence mix is strong, which entails higher turbulence intensity factors. Li et al. (L. Li et al., 2012) propose Equation 8 to calculate the β, factor, obtained from measurements of the Hagupit typhoon, which struck the southern coast of China in 2008. This paper compares the results with those of Masters et al. (Masters, Tieleman and Balderrama, 2010), concluding that both studies are quite consistent.

**(8)**

In a more updated paper using a larger number of hurricanes, Li et al. (L. Li et al., 2015) propose a better adjustment of equation 8 through Equation 9.

**(9)**

As for the determination of u_{*}, friction velocity, two approaches are frequently used: 1) extrapolations from the logarithmic profile and 2) the expression based on the turbulent exchange motion (L. Li et al., 2015; Wang et al., 2011). The second method is used when anemometers are recording measurements in the three directions.

Concerning the turbulence intensity, *I*=σ_{u}⁄ it can be perceived that it tends to decrease as the mean speed increases, and in the hurricane-prone regions there is a tendency of greater intensity than in the extra-tropical regions.

While Fu et al. (Fu et al., 2012) studied the results from the instrumentation of a building in Guangzhou, China during the Megi typhoon, they concluded that the turbulence intensity values were consistent with those resulting from applying the Solari and Piccardo model to determine the variance (G. Solari and Piccardo, 2001), as well as the results of Li et al. (Q. S. Li et al., 2010) regarding observations made in Beijing.

The work of Li et al. (L. Li et al., 2015) compares the experimentally obtained turbulence intensity during the path of typhoons and hurricanes, with the expression proposed by ASCE7-10 (2010), given by Equation 10.

**(10)**

**Table 2.** Formulations for the turbulence integral scale according to certain standards

Where c is 0.30, 0.20 and 0.15 for exposure categories B, C and D of the same standard. Based on the comparison, they concluded that ASCE7-10 (2010) offered good approximations up to 60 m above the ground for categories D and C, while for B it overestimated the value. These results differ from those obtained by the previous work of Li et al. (Q. S. Li et al., 2010), which compares real-scale measurements with empirical values resulting from the formulations proposed by ASCE (10) and AIJ-RLB (11). The comparison was made for the B land category, according to ASCE, and from a height of 47 m above the ground level.

I_{u}=0.1(Zg/Z)^{α+0.05}**(11)**

α is 0.27 for a urban land.

The authors of that study concluded that the experimental values did not properly adjust to any of the empirical profiles of the ASCE or AIJ-RLB codes; thus, their application to high structures is not valid.

**1.4 Simulation of wind speed time series for gaussian processes**

There are several simulation methods for generating stationary stochastic fields with Gaussian behavior. Neural networks (Wu and Kareem, 2011) and wave transformations (Benowitz and Deodatis, 2015) are currently applied alternatives. However, the most used methods are classified according to the following diagram in Figure 1.

(AR) and (MA) models are particular (ARMA) cases (Rossi, Lazzari and Vitaliani, 2004). The work of Samaras et al. (Samaras, Shinzuka and Tsurui, 1985) was one of the first to introduce (ARMA) models in 1985. A special use of these methods is observed in modelling of complex structures and large simulation domains (Di Paola and Gullo, 2001; J. Li and Li, 2012; Rossi et al., 2004). The (ARIMA) models are an extension of the (ARMA) models, and they are used mainly in the non-stationary behavior of wind series. The main strength of the (ARIMA) models lies in their capacity to reveal temporary interdependencies in the series (Zhihua Zhang and Moore, 2015).

The second traditional approach is to use the superposition of trigonometric functions with random phase angles, “spectral representation methods”. Basic ideas were originated in the 50’s (Chen et al., 2014), but its implementation in structural dynamics begun to be studied in the 70’s, with Shinozuka and Jan as precursors (Shinozuka and Jan, 1972). These authors developed their simulation formulas aimed at stochastic processes, including multivariable, multi-dimensional and non-stationary cases; but when the number of variables increased, it turned out to be computationally inefficient. In 1972, with the purpose of improving the computational efficiency, Yang (Yang, 1972) introduced that the sum of the trigonometric functions can be made through the Fast Fourier Transform (FFT) and proposes a formula to simulate envelope processes. In 1996, Deodatis (Deodatis, 1996) developed an algorithm that improved the simulation formulas established by Shinozuka, because it guaranteed the generation of ergodic series by optimizing the implementation of the FFT. However, the use of the FFT increases the computer storage demand, which is a particular problem for simulations of multivariable and multi-dimensional processes that require long data-processing times (Chen et al., 2014; J. Li, Li, He and Shen, 2015). Therefore, in order to improve the computational efficiency, techniques have arisen that require the decomposition of XPSD matrix.

In the 90’s, Li and Kareem (Yosun Li and Kareem, 1991) introduced the concept of XPSD decomposition, with the aim of simulating random stationary processes, with further application possibilities in multivariable and non-stationary random processes (Y. Li and Kareem, 1997). Two decomposition techniques are mainly used: 1) the Cholesky Decomposition and 2) the Proper Orthogonal Decomposition (POD) (Chen et al., 2014; Ubertini and Giuliano, 2010).

The improved procedure, resulting from combining the studies of Shinozuka et al. (Shinozuka and Jan, 1972) and Deodatis (Deodatis, 1996), applies the Cholesky decomposition of the power spectral density matrix.

Solari et al. (Giovanni Solari, Carassale and Tubino, 2007) review the use of POD. According to this work, its use in structural dynamics dates back to the 80’s. The works of Di Paola (Di Paola and Gullo, 2001) and Carassale and Solari (Carassale and Solari, 2006) are some of the most quoted by the scientific community regarding the application of this technique.

Some studies have used this technique to establish a comparison between the results’ precision and the computational efficiency in the same problem. In 2010, Ubertini and Giuliano (Ubertini and Giuliano, 2010) made two numerical examples. In the first, they simulated the wind field to study the behavior of a suspension bridge tower (in 9 nodes); therefore, they analyzed the AR model efficiency, of the Shinozuka and Deodatis method, and the POD, with the improvements proposed by Carassale and Solari (Carassale and Solari, 2006). The second example simulated the wind field for the entire bridge (83 nodes), using only the AR model and the POD. They concluded that the Shinozuka and Deodatis method is the most accurate, but in relation to the first example, it demands 10 times more computing memory and takes approximately 8 times longer. In the first case, the variations of the other two methods were not significant among them. As for the second, they concluded that the POD is more accurate and computationally efficient in relation to time, but the AR consumes 5 times less computing memory.

Concerning the analysis of flexible structures such as towers and tall buildings (Serrano, Mora and Salazar, 2014), transmission towers (Mercanti, Pizzutti, Aguirre, Fank and Möller, 2011; T. Zhang, Lin and Bai, 2013; Zhuoqun Zhang, Li, Li, Wang and Tian, 2013) or guyed towers (Bu, Law and Zhu, 2012; Clobes and Peil, 2011; Gani and Légeron, 2010), and in agreement with these results, a group of authors has been identified who prefer the use of the spectral representation method, through the Cholesky decomposition. However, if they need to simulate large wind fields in connection with structures like suspension or cable-stayed bridges, they select techniques such as POD (Carassale and Solari, 2006; Y. F. Huang and Yang, 2011) or linear filter models (Proppe and Zhang, 2015).

Based on the methods presented, several works have refined them (Chen et al., 2014; Guoqing Huang et al., 2013; J. Li et al., 2011) or combined them. Regarding the latter, the papers of Huang (G. Huang, 2015; Guoqing Huang, 2014) addressing the simulation of non-stationary processes are an example thereof.

**Figure 1.** Diagram of most used methods for simulating Gaussian processes

**2. Conclusions**

From the bibliographic study dealing with the simulation of wind speed time series, it is concluded that, regarding vertical profiles of mean wind speed, both the power law and the log law can be applied to simulate the behavior of the mean wind speed component under storm conditions. As for the formulations to obtain the spectrum, the expression of van Kármán seems to be the most adequate to describe the wind turbulence during the path of strong storms, such as hurricanes. Due to the wide range of criteria associated to the determination of the turbulence scale, the model of Solari and Piccardo (G. Solari and Piccardo, 2001) is proposed. Moreover, concerning the parameter of turbulence intensity factor, a comparison between the application of formulas 7 and 9 is proposed.

In relation to the synthetic generation methods of wind speed time series, it is concluded that digital filtering methods for generating wind time series does not require large storage memories, unlike those using the FFT, but they are less accurate than the spectral representation method. From the two main XPSD decomposition procedures, in the spectral representation method, the Cholesky decomposition leads to higher computational costs related to memory and time, but it is more accurate in predicting the response induced by the wind. This is an important aspect for the analysis of flexible and strongly non-linear structures such as transmission line towers, guyed towers and cable-stayed bridges. The Proper Orthogonal Decomposition (POD) is adequate when it is necessary to densify the discretization of the speed fields.

**3. References**

** An Y., Quan Y. y Gu, M. (2012),** Field Measurement of Wind Characteristics of Typhoon Muifa on the Shanghai World Financial Center. International Journal of Distributed Sensor Networks, 2012, 11. doi: 10.1155/2012/893739

* AS/NZS1170.2-2011. (2011),* Australian/New Zeland Standart: Structural Design actions, Part 2: Wind Actions.

** ASCE7-10. (2010),** American Society of Civil Engineers, Minimum Design Loads for Buildings and Other Structures.

** Benowitz B. A. y Deodatis, G. (2015)**. Simulation of wind velocities on long span structures: A novel stochastic wave based model. Journal of Wind Engineering and Industrial Aerodynamics, 147, 154-163, http://dx.doi.org/10.1016/j.jweia.2015.10.004

** Bu J. Q., Law S. y Zhu X. (2012),** Time-varying Wind Speed Identification from Structural Responses by Modified Iterative Method. Journal of Civil Engineering and Construction, 1(1).

** Carassale L. y Solari G. (2006),** Monte Carlo simulation of wind velocity fields on complex structures. Journal of Wind Engineering and Industrial Aerodynamics, 94(5), 323-339, http://dx.doi.org/10.1016/j.jweia.2006.01.004

** Chen N., Li Y. y Xiang H. (2014),** A new simulation algorithm of multivariate short-term stochastic wind velocity field based on inverse fast Fourier transform. Engineering Structures, 80, 251-259, http://dx.doi.org/10.1016/j.engstruct.2014.09.012

** Ciampoli M., Petrini F. y Augusti G. (2011),** Performance-based wind engineering: towards a general procedure. Structural Safety, 33(6), 367-378, http://dx.doi.org/10.1016/j.strusafe.2011.07.001

** Clobes M. y Peil U. (2011),** Unsteady buffeting wind loads in the time domain and their effect on the life-cycle prediction of guyed masts. Structure and Infrastructure Engineering, 7(1-2), 187-196, http://dx.doi.org/10.1080/15732471003588783

* Cottone G. y Di Paola M. (2011),* Fractional spectral moments for digital simulation of multivariate wind velocity fields. Journal of Wind Engineering and Industrial Aerodynamics, 99(6–7), 741-747, http://dx.doi.org/10.1016/j.jweia.2011.03.006

** Davenport A. G. (1961),** The spectrum of horizontal gustiness near the ground in high winds. Q.J.R. Meteorol. Soc., 87, 194–211, http://dx.doi.org/10.1002/qj.49708737208

* Deodatis G. (1996),* Simulation of Ergodic Multivariate Stochastic Processes. Journal of Engineering Mechanics, 122(8), 778-787, http://dx.doi.org/10.1061/(ASCE)0733-9399(1996)122:8(778)

** Di Paola M. y Gullo, I. (2001),** Digital generation of multivariate wind field processes. Probabilistic Engineering Mechanics, 16(1), http://dx.doi.org/1-10, 10.1016/S0266-8920(99)00032-6

** EN1991-1-4. (2004),** Eurocode 1: Actions on structures — General actions — Part 1-4: Wind actions.

** Fu J.Y., Wu J.R., Xu A., Li Q. S. y Xiao Y. Q. (2012),** Full-scale measurements of wind effects on Guangzhou West Tower. Engineering Structures 35, 120-139, http://dx.doi.org/10.1016/j.engstruct.2011.10.022

** Gani F. y Légeron F. (2010),** Dynamic response of transmission lines guyed towers under wind loading. Canadian Journal of Civil Engineering, 37(3), 450-465, http://dx.doi.org/10.1139/l09-160

** Giammanco I. M., Schroeder J. L. y Powell M. D. (2012),** Observed characteristics of tropical cyclone vertical wind profiles. Wind and Structures, 15(1), 65-86, http://dx.doi.org/10.12989/was.2012.15.1.065

** Giammanco I. M. y Schroeder J. L. (2013),** GPS Dropwindsonde and WSR-88D Observations of Tropical Cyclone Vertical Wind Profiles and Their Characteristics. Weather and Forecasting, 28(1), 77-99, http://dx.doi.org/10.1175/WAF-D-11-00155.1

** Hernandez E. M., Bernal D. y Caracoglia L. (2013),** On-line monitoring of wind-induced stresses and fatigue damage in instrumented structures. Structural Control and Health Monitoring, 20(10), 1291-1302, http://dx.doi.org/10.1002/stc.1536

** Huang G. (2015),** Application of Proper Orthogonal Decomposition in Fast Fourier Transform—Assisted Multivariate Nonstationary Process Simulation. Journal of Engineering Mechanics, 141(7), 04015015, http://dx.doi.org/10.1061/(ASCE)EM.1943-7889.0000923

** Huang G. (2014),** An efficient simulation approach for multivariate nonstationary process: Hybrid of wavelet and spectral representation method. Probabilistic Engineering Mechanics, 37, 74-83, http://dx.doi.org/10.1016/j.probengmech.2014.06.001

* Huang G., Liao H. y Li M. (2013),* New formulation of Cholesky decomposition and applications in stochastic simulation. Probabilistic Engineering Mechanics, 34, 40-47, http://dx.doi.org/10.1016/j.probengmech.2013.04.003

* Huang Y. F. y Yang J. X. (2011),* Simulation of Wind Field on Long-Span Cable-Stayed Bridge Based on SPT. Applied Mechanics and Materials, 90-93, 2451-2455, http://dx.doi.org/ 10.4028/www.scientific.net/AMM.90-93.2451

* ISO 4354. (2009),* Wind action on structures. Geneva. Switzerland.

* IWC, Indian Wind Code (2012),* Wind Loads on Buildings and Structures - Proposed draft & Commentary. Document No: IITK GSDMA-Wind 02-V 50 IS: 875 (Part 3). India.

* Kaimal J. C., Wyngaard J. C., Izumi Y. y Coté O. R. (1972),* Spectral characteristics of surface layer turbulence. Quart J Roy Meteorol Soc, 98, 563-589.

* Koulatsou K., Petrini F., Vernardos S. y Gantes C. J. (2013),* Artificial Time Histories of Wind Actions For Structural Analysis of Wind Turbines. Paper presented at the 2nd International Balkans Conference on Challenges of Civil Engineering, Epoka University, Tirana, Albania. Available in: http://dspace.epoka.edu.al/handle/1/1189

* Li J. y Li C. (2012),* Simulation of non-Gaussian stochastic process with target power spectral density and lower-order moments. Journal of Engineering Mechanics, 138(5), 391-404, http://dx.doi.org/10.1061/(ASCE)EM.1943-7889.0000349

** Li J., Li C. y Chen S. (2011),** Spline-Interpolation-Based FFT Approach to Fast Simulation of Multivariate Stochastic Processes. Mathematical Problems in Engineering, 2011, 24, http://dx.doi.org/10.1155/2011/842183

** Li J., Li C., He L. y Shen J. (2015),** Extended modulating functions for simulation of wind velocities with weak and strong nonstationarity. Renewable Energy, 83, 384-397, http://dx.doi.org/10.1016/j.renene.2015.04.044

** Li L., Kareem A., Xiao Y., Song L. y Zhou C. (2015),** A comparative study of field measurements of the turbulence characteristics of typhoon and hurricane winds. Journal of Wind Engineering and Industrial Aerodynamics, 140, 49-66, http://dx.doi.org/10.1016/j.jweia.2014.12.008

** Li L., Xiao Y., Kareem A., Song L. y Qin P. (2012),** Modeling typhoon wind power spectra near sea surface based on measurements in the South China sea. Journal of Wind Engineering and Industrial Aerodynamics, 104–106, 565-576, http://dx.doi.org/10.1016/j.jweia.2012.04.005

** Li Q. S., Zhi L. y Hu F. (2010),** Boundary layer wind structure from observations on a 325 m tower. Journal of Wind Engineering and Industrial Aerodynamics, 98(12), 818-832, http://dx.doi.org/10.1016/j.jweia.2010.08.001

** Li Y. y Kareem A. (1997),** Simulation of Multivariate Nonstationary Random Processes: Hybrid DFT and Digital Filtering Approach. Journal of Engineering Mechanics, 123(12), 1302-1310, http://dx.doi.org/10.1061/(ASCE)0733-9399(1997)123:12(1302)

** Li Y. y Kareem A. (1991),** Simulation of Multivariate Nonstationary Random Processes by FFT. Journal of Engineering Mechanics, 117(5), 1037-1058, http://dx.doi.org/10.1061/(ASCE)0733-9399(1991)117:5(1037)

** Masters F. J., Tieleman H. W. y Balderrama J. A. (2010),** Surface wind measurements in three Gulf Coast hurricanes of 2005. Journal of Wind Engineering and Industrial Aerodynamics, 98(10–11), 533-547, http://dx.doi.org/10.1016/j.jweia.2010.04.003

** Mercanti J. A, Pizzutti H. D., Aguirre M. A, Fank P. Y. y Möller O. (2011),** Acción del viento como proceso estocástico aplicado a líneas de transmisión de energía eléctrica. Mecánica Computacional, 30(35). Available in: http://www.cimec.org.ar/ojs/index.php/mc/article/viewFile/3945/3862

**Pinheiro M. A. S. (2010),** Simulação numérica da ação de vento turbulento considerando correlação espacial. Paper presented at the Proceedings of the 9th Brazilian Conference on Dynamics Control and their Applications. Serra Negra, SP, Brazil. Available in: http://www.sbmac.org.br/dincon/trabalhos/PDF/methods/69117.pdf

** Proppe C. y Zhang X. (2015),** Influence of Uncertainties on Crosswind Stability of Vehicles. Procedia IUTAM, 13, 98-107, http://dx.doi.org/10.1016/j.piutam.2015.01.006

** Rossi R., Lazzari M. y Vitaliani R. (2004),** Wind field simulation for structural engineering purposes. International Journal for Numerical Methods in Engineering, 61(5), 738-763, http://dx.doi.org/10.1002/nme.1083

** Samaras E., Shinzuka M. y Tsurui A. (1985),** ARMA Representation of Random Processes. Journal of Engineering Mechanics, 111(3), 449-461, http://dx.doi.org/10.1061/(ASCE)0733-9399(1985)111:3(449)

** Serrano O. J. P., Mora E. B. y Salazar A. R. (2014),** Simulación de la componente longitudinal del viento por representación espectral y el análisis dinámico en edificios de cortante. Paper presented at the XIX Congreso Nacional de Ingeniería Estructural, Puerto Vallarta, Jalisco, México. Available in: http://www.smie.org.mx/SMIE_Articulos/co/co_18/te_05/ar_05.pdf

** Shiau B.S. (2000),** Velocity spectra and turbulence statistics at the northeastern coast of Taiwan under high-wind conditions. Journal of Wind Engineering and Industrial Aerodynamics, 88(2–3), 139-151, http://dx.doi.org/10.1016/S0167-6105(00)00045-3

** Shinozuka M. y Jan C. M. (1972),** Digital simulation of random processes and its applications. Journal of Sound and Vibration, 25(1), 111-128, https://doi.org/10.1016/0022-460X(72)90600-1

**Solari G. y Piccardo G. (2001),** Probabilistic 3-D turbulence modeling for gust buffeting of structures. Probabilistic Engineering Mechanics, 16(1), 73-86, http://dx.doi.org/10.1016/S0266-8920(00)00010-2

** Solari G., Carassale L. y Tubino F. (2007),** Proper orthogonal decomposition in wind engineering. Part 1: A state-of-the-art and some prospects. Wind and Structures, 10(2), 153-176, http://dx.doi.org/10.12989/was.2007.10.2.153

** Tamura Y. y Kareem A. (2013),** Advanced Structural Wind Engineering. Japan: Springer. 414 p.p. 6-11, 16-18, 263-268. ISBN 978-4-431-54337-4 (eBook). http://dx.doi.org/10.1007/978-4-431-54337-4

** Tamura Y., Cao S. y Giang L. T. (2012),** Wind characteristics of strong tropical cyclones. Paper presented at the VI National Conference on Wind Engineering, Roorkee at New Delhi, India.

** Torrielli Alessio (2011),** Long term simulation and reliability analysis of structures subjected to the wind action. (PhD), Università degli Studi di Genova, Italia. Available in: http://www.scuoladottorato-sti.unige.it

** Tse K. T., Li S. W., Chan P. W., Mok H. Y. y Weerasuriya A. U. (2013),** Wind profile observations in tropical cyclone events using wind-profilers and doppler SODARs. Journal of Wind Engineering and Industrial Aerodynamics, 115(0), 93-103, http://dx.doi.org/10.1016/j.jweia.2013.01.003

** Ubertini F. y Giuliano F. (2010),** Computer Simulation of Stochastic Wind Velocity Fields for Structural Response Analysis: Comparisons and Applications. Advances in Civil Engineering, 2010. http://dx.doi.org/10.1155/2010/749578

** Van der Hoven Isaac (1957),** Power spectrum of horizontal wind speed in the frequency range from 0.0007 to 900 cycles per hour. Journal of Meteorology, 14(2), 160-164. http://dx.doi.org/10.1175/1520-0469(1957)014<0160:PSOHWS>2.0.CO;2

** Vickery P. J., Wadhera D., Powell M. D. y Chen Y. (2009),** A Hurricane Boundary Layer and Wind Field Model for Use in Engineering Applications. Journal of Applied Meteorology and Climatology, 48(2), 381-405, http://dx.doi.org/10.1175/2008jamc1841.1

** Von Kármán T. (1948),** Progress in the statistical theory of turbulence. Paper presented at the Proceedings of the National Academy of Sciences of the United States of America, 34(11), 530-539. Available in: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC1079162/pdf/pnas01712-0042.pdf

** Wang B., Hu F. y Cheng X. (2011),** Wind gust and turbulence statistics of typhoons in South China. Acta Meteorologica Sinica, 25(1), 113-127, http://dx.doi.org/10.1007/s13351-011-0009-8

** WMO (2010),** Guidelines for converting between various wind averaging periods in tropical cyclone conditions WMO/TD-No.1555. Geneva, Switzerland World Meteorological Organization.

** Wu T. y Kareem A. (2011),** Modeling hysteretic nonlinear behavior of bridge aerodynamics via cellular automata nested neural network. Journal of Wind Engineering and Industrial Aerodynamics, 99(4), 378-388, http://dx.doi.org/10.1016/j.jweia.2010.12.011

** Xu K., Zhao L., Cao S. Y., Yoshida A., Tamura Y. y Ge Y. J. (2014),** Investigation of Spatial Coherences of Aerodynamic Loads on a Streamlined Bridge Deck in an Actively-Controlled Wind Tunnel. Advances in Structural Engineering, 17(1), 53-66, http://dx.doi.org/10.1260/1369-4332.17.1.53

** Yang J. N. (1972),** Simulation of random envelope processes. Journal of Sound and Vibration, 21(1), 73-85, http://dx.doi.org/10.1016/0022-460X(72)90207-6

** Yaojun G., Xinyang J. y Cao S. (2010),** Comparison of APEC wind loading codification and revision of Chinese National Code. Paper presented at the Proceedings of 6th workshop on regional harmonization of wind loading and wind environmental specifications in Asia-Pacific economies, Korea. Available in: http://wind.arch.t-kougei.ac.jp/APECWW/Report/2010/CHN.pdf

** Zhang T., Lin Y. Y. y Bai H. F. (2013a),** Numerical Analysis of Transmission Towers-lines Construction on Wind Forces. Journal of Applied Sciences, 13(9), 1587-1591, http://dx.doi.org/10.3923/jas.2013.1587.1591

** Zhang Z. y Moore J. C. (2015),** Chapter 8 - Autoregressive Moving Average Models. In Zhihua Zhang John C. Moore (Ed.), Mathematical and Physical Fundamentals of Climate Change (239-290). Boston: Elsevier.

** Zhang Z., Li H., Li G., Wang W. y Tian L. (2013b),** The Numerical Analysis of Transmission Tower-Line System Wind-Induced Collapsed Performance. Mathematical Problems in Engineering, 2013, 11, http://dx.doi.org/10.1155/2013/413275

Ingeniera Civil del Instituto Superior Politécnico José Antonio Echeverrría (CUJAE). Máster en Reducción de Desastres, Mención estructural CUJAE. Actual estudiante de doctorado. Profesora del departamento de Estructuras de la Facultad de Ingeniería Civil del CUJAE

E-mail: ingridfl@civil.cujae.edu.cu

Fecha de Recepción: 13/12/2016 Fecha de Aceptación: 22/02/2017

### Refbacks

- There are currently no refbacks.

Copyright (c) 2017 Revista Ingeniería de Construcción

This work is licensed under a Creative Commons Attribution 4.0 International License.