Aerodynamic Performance and Wind-Induced Responses of Large Wind Turbine Systems with Meso-Scale Typhoon Effects

The theoretical system of existing civil engineering typhoon models is too simplified and the simulation accuracy is very low. Therefore, in this work a meso-scale weather forecast model (WRF) based on the non-static Euler equation model was introduced to simulate typhoon “Nuri” with high spatial and temporal resolution, focusing on the comparison of wind direction and wind intensity characteristics before, during and after the landing of the typhoon. Moreover, the effectiveness of the meso-scale typhoon “Nuri” simulation was verified by a comparison between the track of the typhoon center based on minimum sea level pressure and the measured track. In this paper, the aerodynamic performance of large wind turbines under typhoon loads is studied using WRF and CFD nesting technology. A 5 MW wind turbine located in a wind power plant on the southeast coast of China has been chosen as the research object. The average and fluctuating wind pressure distributions as well as airflow around the tower body and eddy distribution on blade and tower surface were compared. A dynamic and time-historical analysis of wind-induced responses under different stop positions was implemented by considering the finite element complete transient method. The influence of the stop position on the wind-induced responses and wind fluttering factor of the system were analyzed. Finally, under a typhoon process, the most unfavorable stop position of the large wind turbine was concluded. The results demonstrated that the internal force and wind fluttering factor of the tower body increased significantly under the typhoon effect. The wind-induced response of the blade closest to the tower body was affected mostly. The wind fluttering factor of this blade was increased by 35%. It was concluded from the analysis that the large wind turbine was stopped during the typhoon. The most unfavorable stop position was at the complete overlapping of the lower blade and the tower body (Condition 1). The safety redundancy reached the maximum when the upper blade overlapped with the tower body completely (Condition 5). Therefore, it is suggested that during typhoons the blade of the wind turbine be rotated to Condition 5.


Introduction
The microturbulence structure of meso-scale typhoons shows significant differences with that of normal winds. Therefore, it is difficult to simulate the wind field characteristics of typhoons near the ground directly through wind tunnel tests or micro-scale CFD. Wind load is the key load of large dynamics and thermodynamics was introduced to simulate typhoon "Nuri" with high spatial and temporal resolution. The characteristics of the 3D near-ground wind field are analyzed emphatically. The validity of the meso-scale typhoon simulation was verified by comparing the simulated and measured paths of the typhoon center. In this paper, a 5 MW wind turbine located in a wind power plant on the southeast coast of China has been chosen as the study object. Near-ground wind field information throughout the landing of the typhoon was acquired and fitted by the nonlinear least squares method. The fitting results provide the initial conditions for the computational fluid mechanics (CFD) model. On this basis, the average and fluctuating wind pressure distributions as well as airflow around the tower body and eddy distribution on the wind turbine system under different stop positions were compared. The effects of stop position on the wind-induced response and wind fluttering factor of the system were discussed. The most unfavorable stop position for the large wind turbine during the typhoon was revealed.

WRF Mode Meshing
The 12# Typhoon "Nuri" in 2008 was chosen as the simulation object. Typhoon "Nuri" brought heavy winds and rains to Hong Kong, Guangdong and Guangxi, and caused serious geological disasters and waterlogging. In order to carry out effective numerical simulation of the typhoon path, pressure field and wind speed field, the simulation includes the whole progress of the typhoon from Saigon, Hong Kong to Zhongshan, Guangdong.
Firstly, meshing shall be performed: the WRF mode was the completely compressible non-static mode. The Arakawa-C grid was used to divide the horizontal computing domain in (Figure 1a), while the vertical grids applied the Euler mass coordinate η along the topology (Figure 1b). Arakawa-C grid can represent both scalar and vector magnitudes [25,26], which is beneficial to improve accuracy in high-resolution simulations. However, the positions on the C-type grid were different. The u and v components of the horizontal wind speed were set on the orthogonal boundary of the quadrilateral cell point region. They were respectively perpendicular to the boundary of the vertical and longitudinal elements. In addition, scalars such as temperature, humidity and pressure were set in the center of the quadrilateral cell. The velocity of Arakawa-C grid and other variables were calculated on the single grid distance, which had good dispersion properties and conservation.
Energies 2018, 11, x FOR PEER REVIEW 3 of 24 In order to systematically study the aerodynamic performance and wind effect characteristics of large-scale wind turbine systems under the typhoon process, a meso-scale WRF model based on fluid dynamics and thermodynamics was introduced to simulate typhoon "Nuri" with high spatial and temporal resolution. The characteristics of the 3D near-ground wind field are analyzed emphatically. The validity of the meso-scale typhoon simulation was verified by comparing the simulated and measured paths of the typhoon center. In this paper, a 5 MW wind turbine located in a wind power plant on the southeast coast of China has been chosen as the study object. Near-ground wind field information throughout the landing of the typhoon was acquired and fitted by the nonlinear least squares method. The fitting results provide the initial conditions for the computational fluid mechanics (CFD) model. On this basis, the average and fluctuating wind pressure distributions as well as airflow around the tower body and eddy distribution on the wind turbine system under different stop positions were compared. The effects of stop position on the wind-induced response and wind fluttering factor of the system were discussed. The most unfavorable stop position for the large wind turbine during the typhoon was revealed.

WRF Mode Meshing
The 12# Typhoon "Nuri" in 2008 was chosen as the simulation object. Typhoon "Nuri" brought heavy winds and rains to Hong Kong, Guangdong and Guangxi, and caused serious geological disasters and waterlogging. In order to carry out effective numerical simulation of the typhoon path, pressure field and wind speed field, the simulation includes the whole progress of the typhoon from Saigon, Hong Kong to Zhongshan, Guangdong.
Firstly, meshing shall be performed: the WRF mode was the completely compressible non-static mode. The Arakawa-C grid was used to divide the horizontal computing domain in (Figure 1a), while the vertical grids applied the Euler mass coordinate η along the topology (Figure 1b). Arakawa-C grid can represent both scalar and vector magnitudes [25,26], which is beneficial to improve accuracy in high-resolution simulations. However, the positions on the C-type grid were different. The u and v components of the horizontal wind speed were set on the orthogonal boundary of the quadrilateral cell point region. They were respectively perpendicular to the boundary of the vertical and longitudinal elements. In addition, scalars such as temperature, humidity and pressure were set in the center of the quadrilateral cell. The velocity of Arakawa-C grid and other variables were calculated on the single grid distance, which had good dispersion properties and conservation.  The triple nesting program was applied for computation of the typhoon wind field of WRF mode. The horizontal resolution was set 13.5 km, 4.5 km and 1.5 km, respectively. The latitude and longitude coordinates of the center of the simulation area were 114.1°E, 22.5°N. The spacing of the outermost grid D01 was 13.5 km and the number of grids was 211 × 211. The grid spacing of the The triple nesting program was applied for computation of the typhoon wind field of WRF mode. The horizontal resolution was set 13.5 km, 4.5 km and 1.5 km, respectively. The latitude and longitude coordinates of the center of the simulation area were 114.1 • E, 22.5 • N. The spacing of the outermost grid D01 was 13.5 km and the number of grids was 211 × 211. The grid spacing of the second layer D02 was 4.5 km and the number of grids was 217 × 217. The third layer mesh spacing D03 in the inner Energies 2019, 12, 3696 4 of 24 layer was 1.5 km and the number of grids was 241 × 241. The map projection applied the Lambert program [5]. In the vertical direction, it is divided into 37 layers, with upper sparse layer and lower dense layer, among which the lower 1000 m is divided into 19 layers. The simulation domain is shown in Figure 2.
Energies 2018, 11, x FOR PEER REVIEW 4 of 24 second layer D02 was 4.5 km and the number of grids was 217 × 217. The third layer mesh spacing D03 in the inner layer was 1.5 km and the number of grids was 241 × 241. The map projection applied the Lambert program [5]. In the vertical direction, it is divided into 37 layers, with upper sparse layer and lower dense layer, among which the lower 1000 m is divided into 19 layers. The simulation domain is shown in Figure 2.

Selection of the Parameterization Scheme
Generally, the hub height of wind turbines is lower than 150 m and at the lower part of the boundary layer thickness, which is significantly affected by the flow quality, humidity, heat transfer and other parameters in the boundary layer [25][26][27]. Therefore, the setting of parametric physics scheme has a direct impact on the boundary layer simulation accuracy of typhoon fields. The parameter settings of each physical scheme of the WRF mode in this paper are shown in Table 1.

Validity Verification and Results Analysis
Air pressure, precipitation and temperature nephograms after the second landing of typhoon "Nuri" are shown in Figure 3. It can be seen from the figure that: (1) the typhoon center was a low-pressure center. After the landing of the typhoon, the underlying surface of the land was coarser than that of the ocean, which led to its heat source being cut off. The updraft was weakened and the surrounding airflow continued to converge towards the center of the typhoon. The simulation

Selection of the Parameterization Scheme
Generally, the hub height of wind turbines is lower than 150 m and at the lower part of the boundary layer thickness, which is significantly affected by the flow quality, humidity, heat transfer and other parameters in the boundary layer [25][26][27]. Therefore, the setting of parametric physics scheme has a direct impact on the boundary layer simulation accuracy of typhoon fields. The parameter settings of each physical scheme of the WRF mode in this paper are shown in Table 1.

Validity Verification and Results Analysis
Air pressure, precipitation and temperature nephograms after the second landing of typhoon "Nuri" are shown in Figure 3. It can be seen from the figure that: (1) the typhoon center was a low-pressure center. After the landing of the typhoon, the underlying surface of the land was coarser than that of the ocean, which led to its heat source being cut off. The updraft was weakened and the surrounding airflow continued to converge towards the center of the typhoon. The simulation results were compared with measured results of the typhoon website of National Meteorological Center of CMA, finding a high agreement in the air pressure distribution. (2) Cloud system precipitation near the typhoon was long-term convective precipitation, and the rainfall in this region was large and related to typhoon intensity and center distance. The Noah road progress program could better simulate the surface flux and airflow convergence field. (3) Affected by the downward transfer of heat from the upper atmosphere, the MYJ boundary layer scheme was used to calculate the turbulent pulsation flux of heat and momentum in the boundary layer in this region with high accuracy. precipitation near the typhoon was long-term convective precipitation, and the rainfall in this region was large and related to typhoon intensity and center distance. The Noah road progress program could better simulate the surface flux and airflow convergence field. (3) Affected by the downward transfer of heat from the upper atmosphere, the MYJ boundary layer scheme was used to calculate the turbulent pulsation flux of heat and momentum in the boundary layer in this region with high accuracy.  By comparison with Figure 4, we can see that there's no strong convection and a partial northeaster before the landing of the typhoon. However, the wind direction of the typhoon turned into a southeaster gradually at the landing point of the typhoon. After landing, the airflow continues to converge towards the center of the typhoon, and the wind velocity streamlines do not show a regular trend. Nevertheless, it gradually turned into a partial northeaster. The abovementioned wind speed streamlines agreed highly with the actual development trend of the typhoon path.
The simulated path was shown by the coordinates of the minimum sea level pressure which were output by the mode every 6 hours in Figures 5 and 6. The blue line was the measured path and the red line was the simulated path. According to the observations: (1) the invasion of clouds around the typhoon increased the wind speed, which increased continuously as the typhoon closed. The the wind speed was higher closer to the center of the typhoon. (2) Typhoon "Nuri" generally moved toward the northwest. There's a certain deviation between the simulated path and measured path, but the simulated typhoon path was consistent with the measured one. (3) The evolution of the typhoon intensity was mainly characterized by the lowest sea level pressure. In the simulation period, typhoon "Nuri" was still in a very mature stage and the minimum sea level pressure was mainly maintained at about 970 hPa. It was slightly lower than the measured result in the late simulation stage. Moreover, the simulated typhoon intensity and variation trend were similar to the measured results.   By comparison with Figure 4, we can see that there's no strong convection and a partial northeaster before the landing of the typhoon. However, the wind direction of the typhoon turned into a southeaster gradually at the landing point of the typhoon. After landing, the airflow continues to converge towards the center of the typhoon, and the wind velocity streamlines do not show a regular trend. Nevertheless, it gradually turned into a partial northeaster. The abovementioned wind speed streamlines agreed highly with the actual development trend of the typhoon path.
The simulated path was shown by the coordinates of the minimum sea level pressure which were output by the mode every 6 hours in Figures 5 and 6. The blue line was the measured path and the red line was the simulated path. According to the observations: (1) the invasion of clouds around the typhoon increased the wind speed, which increased continuously as the typhoon closed. The the wind speed was higher closer to the center of the typhoon. (2) Typhoon "Nuri" generally moved toward the northwest. There's a certain deviation between the simulated path and measured path, but the simulated typhoon path was consistent with the measured one. (3) The evolution of the typhoon intensity was mainly characterized by the lowest sea level pressure. In the simulation period, typhoon "Nuri" was still in a very mature stage and the minimum sea level pressure was mainly maintained at about 970 hPa. It was slightly lower than the measured result in the late simulation stage. Moreover, the simulated typhoon intensity and variation trend were similar to the measured results.
was large and related to typhoon intensity and center distance. The Noah road progress program could better simulate the surface flux and airflow convergence field. (3) Affected by the downward transfer of heat from the upper atmosphere, the MYJ boundary layer scheme was used to calculate the turbulent pulsation flux of heat and momentum in the boundary layer in this region with high accuracy.
(a) Air pressure (b) Precipitation (c) Potential temperature Figure 3. Air pressure, precipitation and potential temperature nephograms at landing.
By comparison with Figure 4, we can see that there's no strong convection and a partial northeaster before the landing of the typhoon. However, the wind direction of the typhoon turned into a southeaster gradually at the landing point of the typhoon. After landing, the airflow continues to converge towards the center of the typhoon, and the wind velocity streamlines do not show a regular trend. Nevertheless, it gradually turned into a partial northeaster. The abovementioned wind speed streamlines agreed highly with the actual development trend of the typhoon path.
The simulated path was shown by the coordinates of the minimum sea level pressure which were output by the mode every 6 hours in Figures 5 and 6. The blue line was the measured path and the red line was the simulated path. According to the observations: (1) the invasion of clouds around the typhoon increased the wind speed, which increased continuously as the typhoon closed. The the wind speed was higher closer to the center of the typhoon. (2) Typhoon "Nuri" generally moved toward the northwest. There's a certain deviation between the simulated path and measured path, but the simulated typhoon path was consistent with the measured one. (3) The evolution of the typhoon intensity was mainly characterized by the lowest sea level pressure. In the simulation period, typhoon "Nuri" was still in a very mature stage and the minimum sea level pressure was mainly maintained at about 970 hPa. It was slightly lower than the measured result in the late simulation stage. Moreover, the simulated typhoon intensity and variation trend were similar to the measured results.

Brief Introduction to the Wind Turbine Model
This paper mainly studies the aerodynamic performance and wind-induced response of a wind turbine system under eight different shutdown positions after the wind turbine stops during atyphoon. In addition, in order to exclude the influence of other interferences and ensure that there is only one variable, the influence of yaw angle and pitch angle under typhoon conditions will not be considered. The chosen 5 MW wind turbine tower was 124 m high, with an upper diameter of 3.0 m and a bottom diameter of 3.5 m. The tower body was a structure of varying length and thickness. The upper wall thickness was 0.06 m and the bottom wall thickness was 0.15 m. The nacelle dimension was 18 m (length) × 6 m (width) × 6 m (height). The cut-out wind speed of the wind wheel was 25 m/s. The blades were distributed uniformly along the circumferential direction at an angle of 120°. The main design parameters of blade and the whole machine are listed in Table 2. Tower, cabin and blade models were built successively, and then the 3D solid model of the large-scale wind turbine was formed through Boolean operations.
According to the relative position between blade and tower body, considering the rotation periodicity of the 3-blade system, the situation of complete overlapping between blade and the tower body was defined as Condition 1 and the remaining conditions were defined every 15° of clockwise rotation, forming eight conditions. In addition, the blade which was overlapped with the tower body in Condition 1 was denoted as Blade 3, and the remaining blades making a clockwise rotation were Blade 1 and Blade 2. The single-period rotation directions of the wind turbine blades under different conditions are shown in Figure 7.

Brief Introduction to the Wind Turbine Model
This paper mainly studies the aerodynamic performance and wind-induced response of a wind turbine system under eight different shutdown positions after the wind turbine stops during atyphoon. In addition, in order to exclude the influence of other interferences and ensure that there is only one variable, the influence of yaw angle and pitch angle under typhoon conditions will not be considered. The chosen 5 MW wind turbine tower was 124 m high, with an upper diameter of 3.0 m and a bottom diameter of 3.5 m. The tower body was a structure of varying length and thickness. The upper wall thickness was 0.06 m and the bottom wall thickness was 0.15 m. The nacelle dimension was 18 m (length) × 6 m (width) × 6 m (height). The cut-out wind speed of the wind wheel was 25 m/s. The blades were distributed uniformly along the circumferential direction at an angle of 120°. The main design parameters of blade and the whole machine are listed in Table 2. Tower, cabin and blade models were built successively, and then the 3D solid model of the large-scale wind turbine was formed through Boolean operations.
According to the relative position between blade and tower body, considering the rotation periodicity of the 3-blade system, the situation of complete overlapping between blade and the tower body was defined as Condition 1 and the remaining conditions were defined every 15° of clockwise rotation, forming eight conditions. In addition, the blade which was overlapped with the tower body in Condition 1 was denoted as Blade 3, and the remaining blades making a clockwise rotation were Blade 1 and Blade 2. The single-period rotation directions of the wind turbine blades under different conditions are shown in Figure 7.

Brief Introduction to the Wind Turbine Model
This paper mainly studies the aerodynamic performance and wind-induced response of a wind turbine system under eight different shutdown positions after the wind turbine stops during atyphoon. In addition, in order to exclude the influence of other interferences and ensure that there is only one variable, the influence of yaw angle and pitch angle under typhoon conditions will not be considered. The chosen 5 MW wind turbine tower was 124 m high, with an upper diameter of 3.0 m and a bottom diameter of 3.5 m. The tower body was a structure of varying length and thickness. The upper wall thickness was 0.06 m and the bottom wall thickness was 0.15 m. The nacelle dimension was 18 m (length) × 6 m (width) × 6 m (height). The cut-out wind speed of the wind wheel was 25 m/s. The blades were distributed uniformly along the circumferential direction at an angle of 120 • . The main design parameters of blade and the whole machine are listed in Table 2. Tower, cabin and blade models were built successively, and then the 3D solid model of the large-scale wind turbine was formed through Boolean operations.
According to the relative position between blade and tower body, considering the rotation periodicity of the 3-blade system, the situation of complete overlapping between blade and the tower body was defined as Condition 1 and the remaining conditions were defined every 15 • of clockwise rotation, forming eight conditions. In addition, the blade which was overlapped with the tower body in Condition 1 was denoted as Blade 3, and the remaining blades making a clockwise rotation were Blade 1 and Blade 2. The single-period rotation directions of the wind turbine blades under different conditions are shown in Figure 7.

Computational Domain and Meshing
In order to ensure the full development of flows, the calculation domain was set as 12D×5D×5D (flow direction X × spanwise Y × vertical Z, and D is the diameter of wind wheel). The wind turbine was located 3D away from the entrance of the computing domain. Because the blade surface is twisted and complex, the mixed mesh was adopted. In addition, the whole computing domain was divided into internal and external parts: the core area adopted tetrahedral grid, and the local grid around the wind turbine was encrypted; the peripheral area adopted high-quality hexahedral grid, with a total number of 13.4 million. The specific division of the computing domain and grid are shown in Figure 8.

Micro-Scale Typhoon Simulation
According to the analysis, the typhoon wind speed distribution close to the near-ground surface was relatively smooth at different moments, while the wind speed distribution at high

Computational Domain and Meshing
In order to ensure the full development of flows, the calculation domain was set as 12D × 5D × 5D (flow direction X × spanwise Y × vertical Z, and D is the diameter of wind wheel). The wind turbine was located 3D away from the entrance of the computing domain. Because the blade surface is twisted and complex, the mixed mesh was adopted. In addition, the whole computing domain was divided into internal and external parts: the core area adopted tetrahedral grid, and the local grid around the wind turbine was encrypted; the peripheral area adopted high-quality hexahedral grid, with a total number of 13.4 million. The specific division of the computing domain and grid are shown in Figure 8.

Micro-Scale Typhoon Simulation
According to the analysis, the typhoon wind speed distribution close to the near-ground surface was relatively smooth at different moments, while the wind speed distribution at high altitude was relatively disordered. The typhoon wind speed at 20:00 on 22nd August was significantly lower than other moments. The reason was that the typhoon had landed at this moment and the underlying surface changed from the ocean surface to land surface, thus dissipating the energy gradually. altitude was relatively disordered. The typhoon wind speed at 20:00 on 22nd August was significantly lower than other moments. The reason was that the typhoon had landed at this moment and the underlying surface changed from the ocean surface to land surface, thus dissipating the energy gradually.  According Figure 9, where the near surface wind speed in the simulated central area and the typhoon near surface wind speed fitting curve obtained by fitting according to the nonlinear least square principle are shown, it can be observed that the near ground typhoon field is well fitted (simulation accuracy was 94.78%). The typhoon wind speed at the height of 10 m was relatively high and it increased slowly with the altitude. At the same time, in order to qualitatively and quantitatively compare the difference between the benign wind and the typhoon wind field, the typhoon field and the normal wind field have the same wind speed when the Class A geomorphologic gradient height is 300 m in the definition specification (GB50009-2012) [28], and the near-ground typhoon profile is used as the wind speed boundary conditions in the subsequent micro-scale wind field numerical simulation.     According Figure 9, where the near surface wind speed in the simulated central area and the typhoon near surface wind speed fitting curve obtained by fitting according to the nonlinear least square principle are shown, it can be observed that the near ground typhoon field is well fitted (simulation accuracy was 94.78%). The typhoon wind speed at the height of 10 m was relatively high and it increased slowly with the altitude. At the same time, in order to qualitatively and quantitatively compare the difference between the benign wind and the typhoon wind field, the typhoon field and the normal wind field have the same wind speed when the Class A geomorphologic gradient height is 300 m in the definition specification (GB50009-2012) [28], and the near-ground typhoon profile is used as the wind speed boundary conditions in the subsequent micro-scale wind field numerical simulation.
altitude was relatively disordered. The typhoon wind speed at 20:00 on 22nd August was significantly lower than other moments. The reason was that the typhoon had landed at this moment and the underlying surface changed from the ocean surface to land surface, thus dissipating the energy gradually.  According Figure 9, where the near surface wind speed in the simulated central area and the typhoon near surface wind speed fitting curve obtained by fitting according to the nonlinear least square principle are shown, it can be observed that the near ground typhoon field is well fitted (simulation accuracy was 94.78%). The typhoon wind speed at the height of 10 m was relatively high and it increased slowly with the altitude. At the same time, in order to qualitatively and quantitatively compare the difference between the benign wind and the typhoon wind field, the typhoon field and the normal wind field have the same wind speed when the Class A geomorphologic gradient height is 300 m in the definition specification (GB50009-2012) [28], and the near-ground typhoon profile is used as the wind speed boundary conditions in the subsequent micro-scale wind field numerical simulation.

Large Eddy Simulation (LES) Governing Equation
In the study of structural wind resistance, fluid was considered as viscous incompressible and the transient N-S equation was processed by space averaging, thus getting the LES governing equation [29,30]: where ρ is the air density. t is time. v is the kinematic coefficient of air viscosity. µ i and µ j are speeds along three directions after smoothing (i = 1, 2, 3 and j = 1, 2, 3). x i and x j are coordinate vectors along three directions in the space. p is the pressure after filtering. τ ij is the open term in the N-S equation after space averaging, that is, the subgrid-scale stress: According to the subgrid model based on eddy viscosity hypothesis proposed by Smagorinsky and Boussinesq hypothesis, the subgrid stress can be expressed as: where S ij was the solvable scale strain rate tensor. τ kk was the isotropic part of subgrid-scale stress, including the stress term after smoothing. δ ij was the Kronecker delta function. µ t was the subgrid-scale turbulence eddy viscosity, which generally applied the Smagorinsky hypothesis: where C S is the Smagorinsky constant. It generally ranges from 0.1~0.23 and it is determined as 0.1 in this paper. The strain rate tensor is S = 2S ij S ij . ∆ is the spatial grid scale, where ∆ x , ∆ y and ∆ z are grid sizes along the x, y and z directions. This is the standard Smagorinsky subgrid-scale model which can simulate micro-turbulence flow accurately. Considering the high-rise special structure (wind turbine system) with high accuracy of simulation, it is necessary to use the above model to carry out ultra-high precision numerical simulation in order to achieve better results. Based on the abovementioned subgrid-scale model and fitted typhoon profile, the User Defined Function file was used to integrate the model into the Fluent software ( Figure 10) to realize LES of the wind turbine system under meso-and micro-scale nesting conditions and different stop positions.

Large Eddy Simulation (LES) Governing Equation
In the study of structural wind resistance, fluid was considered as viscous incompressible and the transient N-S equation was processed by space averaging, thus getting the LES governing equation [29,30]: where ρ is the air density. t is time. v is the kinematic coefficient of air viscosity. i μ and j μ are speeds along three directions after smoothing (i = 1, 2, 3 and j = 1, 2, 3). i x and j x are coordinate vectors along three directions in the space. p is the pressure after filtering. ij τ is the open term in the N-S equation after space averaging, that is, the subgrid-scale stress: According to the subgrid model based on eddy viscosity hypothesis proposed by Smagorinsky and Boussinesq hypothesis, the subgrid stress can be expressed as: where ij S was the solvable scale strain rate tensor. kk τ was the isotropic part of subgrid-scale stress, including the stress term after smoothing. ij δ was the Kronecker delta function. t μ was the subgrid-scale turbulence eddy viscosity, which generally applied the Smagorinsky hypothesis: where CS is the Smagorinsky constant. It generally ranges from 0.1~0.23 and it is determined as 0.1 in this paper. The strain rate tensor is Δ is the spatial grid scale, Δ are grid sizes along the x, y and z directions. This is the standard Smagorinsky subgrid-scale model which can simulate micro-turbulence flow accurately. Considering the high-rise special structure (wind turbine system) with high accuracy of simulation, it is necessary to use the above model to carry out ultra-high precision numerical simulation in order to achieve better results. Based on the abovementioned subgrid-scale model and fitted typhoon profile, the User Defined Function file was used to integrate the model into the Fluent software ( Figure 10) to realize LES of the wind turbine system under meso-and micro-scale nesting conditions and different stop positions.

Boundary Conditions and Parameter Setting
The inlet and outlet of the calculation domain were the velocity inlet and pressure boundary outlet, respectively. The top and side surfaces applied the symmetry boundaries of free slip wall, while the ground and structural surface used the non-slip wall. The fluctuating wind field was defined by UDF. Meanwhile, the exponential wind profile and turbulence profile of class A geomorphic atmospheric boundary layer were set. The ground roughness in the wind profile was set 0.076. The theoretical formula of the fluctuating wind load generation based on Fluent [31] was: x is the coordinate vector, t is time and k m is the wave number. ξ ∈ N(0, 1). I s is the turbulence integral scale, used to adjust spatial correlation. a is evenly distributed between 0 and 1. In this calculation case, k max = 500 and N = 100.
3D single-precision and a separate solver was adopted for numerical calculation, and incompressible flow field was used for air and the wind field. The boundary conditions of calculation domain were set as shown in Figure 11. Based on the SIMPLEC algorithm of Large Eddy Simulation, the pressure velocity nesting method had good convergence and was suitable for large eddy simulation calculation with small time step. Grid tilt correction was set in the calculation process to improve the mixed grid calculation effect, and the computational residual of the N-S equation was set as 1 × 10 −6 .

Boundary Conditions and Parameter Setting
The inlet and outlet of the calculation domain were the velocity inlet and pressure boundary outlet, respectively. The top and side surfaces applied the symmetry boundaries of free slip wall, while the ground and structural surface used the non-slip wall. The fluctuating wind field was defined by UDF. Meanwhile, the exponential wind profile and turbulence profile of class A geomorphic atmospheric boundary layer were set. The ground roughness in the wind profile was set 0.076. The theoretical formula of the fluctuating wind load generation based on Fluent [31] was: x is the coordinate vector, t is time and km is the wave number.
Is is the turbulence integral scale, used to adjust spatial correlation. a is evenly distributed between 0 and 1. In this calculation case, kmax = 500 and N = 100.
3D single-precision and a separate solver was adopted for numerical calculation, and incompressible flow field was used for air and the wind field. The boundary conditions of calculation domain were set as shown in Figure 11. Based on the SIMPLEC algorithm of Large Eddy Simulation, the pressure velocity nesting method had good convergence and was suitable for large eddy simulation calculation with small time step. Grid tilt correction was set in the calculation process to improve the mixed grid calculation effect, and the computational residual of the N-S equation was set as 1 × 10 −6 .

Average Wind Pressure Coefficient
In Figure 12, it is shown that: (1) the blade disturbance to the tower body was dominated by the proximity effect. When the incoming blade overlapped with the tower body completely (Condition 1), a narrow long negative wind pressure coefficient zone was developed on the windward side,

Average Wind Pressure Coefficient
In Figure 12, it is shown that: (1) the blade disturbance to the tower body was dominated by the proximity effect. When the incoming blade overlapped with the tower body completely (Condition 1), a narrow long negative wind pressure coefficient zone was developed on the windward side, which covered an area of 0.5S (S = windward/leeside area of the tower body). The maximum negative pressure was -0.8. With the reduction of shielding area, the negative pressure zone shrank gradually.
(2) Besides, a large negative pressure zone occurred on the leeside of tower body under Condition 1. This region shrank gradually with the weakening of the shielding effect.
Energies 2018, 11, x FOR PEER REVIEW 11 of 24 which covered an area of 0.5S (S = windward/leeside area of the tower body). The maximum negative pressure was -0.8. With the reduction of shielding area, the negative pressure zone shrank gradually.
(2) Besides, a large negative pressure zone occurred on the leeside of tower body under Condition 1. This region shrank gradually with the weakening of the shielding effect.
(a) Windward side (b) Leeside  The contour map of wind pressure coefficients of three blades under eight conditions under typhoon loads is shown in Figure 13.  It shows that: (1) the height of Blade 3 increased gradually in the clockwise rotation, accompanied with gradual increase of positive pressure on the windward side of blade. This was because the wind speed underwent an exponential increase with the height and the wind pressure coefficient was relatively small at a relatively low position. A similar trend was observed in Blade 1 and Blade 2.
(2) Regions with different reduction amplitudes of negative pressure were developed on the leeside of blades. Blade 3 which was close to the tower structure presented tendency changes. The negative pressure reduction region was expanded gradually with the increase of relative distance. (3) Under Condition 1, the absolute value of the wind pressure coefficient at root the Blade 3 on the windward side and at about 1/3 on the leeside along the span direction were 0.45 and 0.4, respectively. The wind pressure coefficients at these two positions were significantly higher than those under the remaining conditions. In addition, the absolute value of wind pressure coefficient decreased gradually with the narrowing of the shielding area. However, the wind pressure coefficients of Blade 1 and Blade 2 presented exponential changes with height. Blades shielded the tower body, and the tower blade also would disturb the blades.

Coefficient of Fluctuating Wind Pressure
According to the relative position of tower body and blades, the tower body was divided into a non-disturbed zone (0-64 m) and a disturbed zone (64-124 m). The time-history curves of the wind pressure coefficient on the leeside at 60 m and 110 m height of the tower body are shown. Figure 14 shows that when the tower body was not shielded by blades, there's a small disparity between the wind pressure coefficients at 60 m and 110 m. However, such difference increased when the blades began to shield the two zones of the tower body, reaching the peak under Condition 1. Later, the difference was narrowed with the reduction of shielding area.

Coefficient of Fluctuating Wind Pressure
According to the relative position of tower body and blades, the tower body was divided into a non-disturbed zone (0-64 m) and a disturbed zone (64-124 m). The time-history curves of the wind pressure coefficient on the leeside at 60 m and 110 m height of the tower body are shown. Figure 14 shows that when the tower body was not shielded by blades, there's a small disparity between the wind pressure coefficients at 60 m and 110 m. However, such difference increased when the blades began to shield the two zones of the tower body, reaching the peak under Condition 1. Later, the difference was narrowed with the reduction of shielding area.  From the Figures 15 and 16, It can be seen that: (1) The absolute value of the wind pressure coefficient in the disturbed zone of the windward side and lateral surface is significantly greater than the non-disturbed section. The mean wind pressure coefficient at 60 m height of the windward side in Condition 1 was higher than the numerical value at 110 m, indicating that blade shielding could influence the wind pressure coefficient of the tower body markedly. Additionally, the RMS of wind pressure on the windward side and lateral surface in the disturbed zone under no blade shielding effect was higher than that on the undisturbed zone.    Extreme wind pressure is crucial to study the mechanical performance of the tower body and blade in the wind turbine system. The calculation formula (GB50009-2012) was: where Ĉ pi , C pi and pi σ are extreme, mean and RMS of wind pressure coefficient at the measuring point i. g is the peak factor at the measuring point i, which is determined 2.5 (GB50009-2012) in this paper. Distributions of extreme wind pressures on the tower body and blade under different conditions are shown in Figure 17. The results demonstrate that: (1) the extreme wind pressure on the windward side was similar at different heights, but it was significantly different on the lateral surface. The split inflow at the blade disturbed stresses on the tower body during the falling process, thus influencing extreme wind pressure distribution in the disturbed zone significantly. This phenomenon was weakened gradually with the narrowing of the shielded area.     Extreme wind pressure is crucial to study the mechanical performance of the tower body and blade in the wind turbine system. The calculation formula (GB50009-2012) was: where Ĉ pi , C pi and pi σ are extreme, mean and RMS of wind pressure coefficient at the measuring point i. g is the peak factor at the measuring point i, which is determined 2.5 (GB50009-2012) in this paper.
Distributions of extreme wind pressures on the tower body and blade under different conditions are shown in Figure 17. The results demonstrate that: (1) the extreme wind pressure on the windward side was similar at different heights, but it was significantly different on the lateral surface. The split inflow at the blade disturbed stresses on the tower body during the falling process, thus influencing extreme wind pressure distribution in the disturbed zone significantly. This phenomenon was weakened gradually with the narrowing of the shielded area.  Extreme wind pressure is crucial to study the mechanical performance of the tower body and blade in the wind turbine system. The calculation formula (GB50009-2012) was: whereĈ pi , C pi and σ pi are extreme, mean and RMS of wind pressure coefficient at the measuring point i. g is the peak factor at the measuring point i, which is determined 2.5 (GB50009-2012) in this paper. Distributions of extreme wind pressures on the tower body and blade under different conditions are shown in Figure 17. The results demonstrate that: (1) the extreme wind pressure on the windward side was similar at different heights, but it was significantly different on the lateral surface. The split inflow at the blade disturbed stresses on the tower body during the falling process, thus influencing extreme wind pressure distribution in the disturbed zone significantly. This phenomenon was weakened gradually with the narrowing of the shielded area. By comparing previous research results under normal wind [22,23], the maximum absolute value (2.82) of average wind pressure under typhoon conditions was achieved at the lateral surface at the tower top, which was 0.27 lower than that under normal wind. The maximum RMS of wind pressure (0.82) was on the windward side at the tower top, which was 0.238 higher than that under normal wind.

Characteristics of Airflow around the Tower Body
Wind speed streamlines in the undisturbed zone and disturbed zone under different conditions under typhoon loads are shown in Figures 18 and 19. Clearly, the inflow at 60 m was not hindered by the upstream blade. It was split at 0° under different conditions and the inflow split point was basically consistent. Backflow and different scales of eddies were generated on the leeside of the tower body. Since the disturbed zone was shielded by the upstream blade to different extents, the airflow around the tower body was significantly different. Airflow around the tower body was affected the most under Condition 1. Later, such influences were weakened gradually with the narrowing of shielding area and disappeared completely when the Blade 2 and Blade 3 were symmetrically placed around the tower body (Condition 5).  By comparing previous research results under normal wind [22,23], the maximum absolute value (2.82) of average wind pressure under typhoon conditions was achieved at the lateral surface at the tower top, which was 0.27 lower than that under normal wind. The maximum RMS of wind pressure (0.82) was on the windward side at the tower top, which was 0.238 higher than that under normal wind.

Characteristics of Airflow around the Tower Body
Wind speed streamlines in the undisturbed zone and disturbed zone under different conditions under typhoon loads are shown in Figures 18 and 19. Clearly, the inflow at 60 m was not hindered by the upstream blade. It was split at 0 • under different conditions and the inflow split point was basically consistent. Backflow and different scales of eddies were generated on the leeside of the tower body. Since the disturbed zone was shielded by the upstream blade to different extents, the airflow around the tower body was significantly different. Airflow around the tower body was affected the most under Condition 1. Later, such influences were weakened gradually with the narrowing of shielding area and disappeared completely when the Blade 2 and Blade 3 were symmetrically placed around the tower body (Condition 5). By comparing previous research results under normal wind [22,23], the maximum absolute value (2.82) of average wind pressure under typhoon conditions was achieved at the lateral surface at the tower top, which was 0.27 lower than that under normal wind. The maximum RMS of wind pressure (0.82) was on the windward side at the tower top, which was 0.238 higher than that under normal wind.

Characteristics of Airflow around the Tower Body
Wind speed streamlines in the undisturbed zone and disturbed zone under different conditions under typhoon loads are shown in Figures 18 and 19. Clearly, the inflow at 60 m was not hindered by the upstream blade. It was split at 0° under different conditions and the inflow split point was basically consistent. Backflow and different scales of eddies were generated on the leeside of the tower body. Since the disturbed zone was shielded by the upstream blade to different extents, the airflow around the tower body was significantly different. Airflow around the tower body was affected the most under Condition 1. Later, such influences were weakened gradually with the narrowing of shielding area and disappeared completely when the Blade 2 and Blade 3 were symmetrically placed around the tower body (Condition 5).

Vorticity Distribution
Vorticity reflects the strength and direction of eddy. Vorticity distributions at the disturbed zone under different conditions are shown in Figure 20. The vorticity distribution differed significantly under different stop positions. The existence of the tower body changed the trailing vortex development at the blades and the trailing vortex of the blades destroyed the separated flow of eddy. The tower body and blades influenced each other mutually. The vorticity distribution between the blades and tower body was dominated by the proximity effect. The existence of the blade wake caused a large-scaled vorticity increment zone. The vorticity increment zone expands with the decrease of relative distance between the blade and tower body. The disturbance redundancy between the tower body and blades under Conditions 4-6 was the most significant, so the interaction could be overlooked.

Vorticity Distribution
Vorticity reflects the strength and direction of eddy. Vorticity distributions at the disturbed zone under different conditions are shown in Figure 20. The vorticity distribution differed significantly under different stop positions. The existence of the tower body changed the trailing vortex development at the blades and the trailing vortex of the blades destroyed the separated flow of eddy. The tower body and blades influenced each other mutually. The vorticity distribution between the blades and tower body was dominated by the proximity effect. The existence of the blade wake caused a large-scaled vorticity increment zone. The vorticity increment zone expands with the decrease of relative distance between the blade and tower body. The disturbance redundancy between the tower body and blades under Conditions 4-6 was the most significant, so the interaction could be overlooked.  By comparing previous research results under normal wind conditions [22,23], the shielding effect of blades to the tower body played the amplification effect under typhoon loads. The shielding effect of blades to the tower body became more significant and the vorticity increment zone decreased with the weakening of shielding effect.

Finite Element Modeling and Dynamic Characteristics
The wind turbine tower-blade integrated finite element model under different conditions was established based on ANSYS. The tower body and blade was simulated by Shell63 units and the cabin was simulated by Beam188 units. The internal structure was neglected. The round raft basis unit was Solid65 with a diameter of 30 m and a height of 3 m. The foundation bottom was solidified. The foundation and basis were simulated by Combin14 spring unit.
The first 100 orders of inherent frequency distribution are shown in Figure 21. It shows the numerical value of the inherent frequency under Condition 1, but only expressed increase/decrease amplitudes under the remaining conditions relative to the value under Condition 1. By comparison we know that: (1) the frequency distribution of the wind turbine was similar under different conditions. The low-order frequency of the structure under Condition 1 was significantly higher than those under the remaining conditions. However, the high-order frequency was slightly lower than those under conditions 2~8. The blade rotation influenced the low-order frequency of the wind turbine system significantly, but influenced the high-order frequency slightly. (2) Under different conditions, there's a small frequency difference among some orders and the frequency increase/decrease of the structure concentrated on the horizontal level of 0 Hz. The fundamental frequency of the structure under the Condition 1 was 0.14 Hz. It decreased gradually with the increase of relative distance between the blade and tower body. By comparing previous research results under normal wind conditions [22,23], the shielding effect of blades to the tower body played the amplification effect under typhoon loads. The shielding effect of blades to the tower body became more significant and the vorticity increment zone decreased with the weakening of shielding effect.

Finite Element Modeling and Dynamic Characteristics
The wind turbine tower-blade integrated finite element model under different conditions was established based on ANSYS. The tower body and blade was simulated by Shell63 units and the cabin was simulated by Beam188 units. The internal structure was neglected. The round raft basis unit was Solid65 with a diameter of 30 m and a height of 3 m. The foundation bottom was solidified. The foundation and basis were simulated by Combin14 spring unit.
The first 100 orders of inherent frequency distribution are shown in Figure 21. It shows the numerical value of the inherent frequency under Condition 1, but only expressed increase/decrease amplitudes under the remaining conditions relative to the value under Condition 1. By comparison we know that: (1) the frequency distribution of the wind turbine was similar under different conditions. The low-order frequency of the structure under Condition 1 was significantly higher than those under the remaining conditions. However, the high-order frequency was slightly lower than those under conditions 2~8. The blade rotation influenced the low-order frequency of the wind turbine system significantly, but influenced the high-order frequency slightly. (2) Under different conditions, there's a small frequency difference among some orders and the frequency increase/decrease of the structure concentrated on the horizontal level of 0 Hz. The fundamental frequency of the structure under the Condition 1 was 0.14 Hz. It decreased gradually with the increase of relative distance between the blade and tower body. To sum up, the stop position of blades wouldn't influence the frequency of the wind turbine significantly. Therefore, the list of typical vibration modes of wind turbines was given by taking working condition 1 as an example. It can be seen from Table 3 that the low-order mode of vibration was mainly manifested by blade driving front-to-rear swing and left-to-right swing of the cabin and tower body. The swing became more violent with the increase of fundamental frequency. The structural deformation and instability of the tower and the blade appeared in the higher mode. To sum up, the stop position of blades wouldn't influence the frequency of the wind turbine significantly. Therefore, the list of typical vibration modes of wind turbines was given by taking working condition 1 as an example. It can be seen from Table 3 that the low-order mode of vibration was mainly manifested by blade driving front-to-rear swing and left-to-right swing of the cabin and tower body. The swing became more violent with the increase of fundamental frequency. The structural deformation and instability of the tower and the blade appeared in the higher mode.

Responses of the Tower Body
The negative value means downwind pointing with respect to the tower body and the positive value means pointing out from the tower against the wind direction in Figure 22. It can be known from the comparison that: (1) Figure 23. The RMS on windward and leeward of the tower are large. RMS of radial displacement on leeside of the tower body was high when the inflow split and fell at blades. Besides, the relative position between blade root and tower body would influence RMS of leeside significantly. With the reduction of relative distance, the RMS increased gradually and reached the significance level under Condition 6. The maximum value was close to −0.1.

Responses of the Tower Body
The negative value means downwind pointing with respect to the tower body and the positive value means pointing out from the tower against the wind direction in Figure 22. It can be known from the comparison that: (1)  Distributions of RMS of radial displacement under different conditions are shown in Figure 23. The RMS on windward and leeward of the tower are large. RMS of radial displacement on leeside of the tower body was high when the inflow split and fell at blades. Besides, the relative position between blade root and tower body would influence RMS of leeside significantly. With the reduction of relative distance, the RMS increased gradually and reached the significance level under Condition 6. The maximum value was close to −0.1.

Responses of the Tower Body
The negative value means downwind pointing with respect to the tower body and the positive value means pointing out from the tower against the wind direction in Figure 22. It can be known from the comparison that: (1)  Distributions of RMS of radial displacement under different conditions are shown in Figure 23. The RMS on windward and leeward of the tower are large. RMS of radial displacement on leeside of the tower body was high when the inflow split and fell at blades. Besides, the relative position between blade root and tower body would influence RMS of leeside significantly. With the reduction of relative distance, the RMS increased gradually and reached the significance level under Condition 6. The maximum value was close to −0.1.

Responses of the Tower Body
The negative value means downwind pointing with respect to the tower body and the positive value means pointing out from the tower against the wind direction in Figure 22. It can be known from the comparison that: (1)  Distributions of RMS of radial displacement under different conditions are shown in Figure 23. The RMS on windward and leeward of the tower are large. RMS of radial displacement on leeside of the tower body was high when the inflow split and fell at blades. Besides, the relative position between blade root and tower body would influence RMS of leeside significantly. With the reduction of relative distance, the RMS increased gradually and reached the significance level under Condition 6. The maximum value was close to −0.1.

Responses of the Tower Body
The negative value means downwind pointing with respect to the tower body and the positive value means pointing out from the tower against the wind direction in Figure 22. It can be known from the comparison that: (1)  Distributions of RMS of radial displacement under different conditions are shown in Figure 23. The RMS on windward and leeward of the tower are large. RMS of radial displacement on leeside of the tower body was high when the inflow split and fell at blades. Besides, the relative position between blade root and tower body would influence RMS of leeside significantly. With the reduction of relative distance, the RMS increased gradually and reached the significance level under Condition 6. The maximum value was close to −0.1. Distributions of RMS of radial displacement under different conditions are shown in Figure 23. The RMS on windward and leeward of the tower are large. RMS of radial displacement on leeside of the tower body was high when the inflow split and fell at blades. Besides, the relative position between blade root and tower body would influence RMS of leeside significantly. With the reduction of relative distance, the RMS increased gradually and reached the significance level under Condition 6. The maximum value was close to −0.1.                     Distributions of mean and RMS of shearing force of blades under different conditions are shown in Figures 28 and 29. The results demonstrated that: (1) at the roots of the three blades, the mean shearing force increased sharply in the beginning and then dropped suddenly. It tended to be stable with the increase of span length. The RMS of shearing force increased sharply at first, then became stable and finally increased dramatically with the increase span length. RMS at three blade tips was significantly higher than those at blade roots.      Distributions of mean and RMS of shearing force of blades under different conditions are shown in Figures 28 and 29. The results demonstrated that: (1) at the roots of the three blades, the mean shearing force increased sharply in the beginning and then dropped suddenly. It tended to be stable with the increase of span length. The RMS of shearing force increased sharply at first, then became stable and finally increased dramatically with the increase span length. RMS at three blade tips was significantly higher than those at blade roots. (2) Under Condition 5 and Condition 6, the mean of shearing force at root of Blade 1 increased dramatically. Under Condition 1 and Condition 2, the mean of shearing force of Blade 2 increased dramatically, indicating that shearing force at blade roots was positively related with blade height. Strengthening the blade roots in the structural design is suggested.  Distributions of mean and RMS of shearing force of blades under different conditions are shown in Figures 28 and 29. The results demonstrated that: (1) at the roots of the three blades, the mean shearing force increased sharply in the beginning and then dropped suddenly. It tended to be stable with the increase of span length. The RMS of shearing force increased sharply at first, then became stable and finally increased dramatically with the increase span length. RMS at three blade tips was significantly higher than those at blade roots. (2) Under Condition 5 and Condition 6, the mean of shearing force at root of Blade 1 increased dramatically. Under Condition 1 and Condition 2, the mean of shearing force of Blade 2 increased dramatically, indicating that shearing force at blade roots was positively related with blade height. Strengthening the blade roots in the structural design is suggested. (3) The variation trend of Blade 3 was not observed in Blade 1 and Blade 2. Under Condition 1 and Condition 2, shearing force increased significantly. This was closely related with the tower-blade interaction. According to previous research results under normal wind conditions [22,23] the typhoon intensified the wind-induced responses of the tower body. The maximum mean radial displacement was increased by about 35%. Mean and RMS of blade responses were influenced slightly. It was suggested to strengthen the tower body in regions with frequent occurrences of typhoons.

Wind Fluttering Factor
The wind fluttering factor distributions of the tower body under different conditions when using internal stresses and displacement as the equivalent goals are shown in Figure 30. It was found that: (1) the wind fluttering factor of tower body under typhoon loads was positively related with height. Under different conditions, the wind fluttering factor when using radial displacement and the meridian axial force at 0° as the equivalent goal was mainly in the range of 2.0-4.0. The wind fluttering factor when using the meridian axial force at 0° as the equivalent goal was slightly larger than that when using the radial displacement as the equivalent goal. Distortion was developed at the tower top. (2) The wind fluttering factor was similar between Condition 2 and Conditions 8, as well as between the Condition 2 and Condition 8. Given two equivalent goals, the wind fluttering factor increased as the shielding effect weakened. The wind fluttering factor under Condition 4 and Condition 6 was significantly higher than those under rest conditions. (3) The wind fluttering factor distribution was gentle under the Condition 1, mainly concentrated in the 2.0-2.5 range. According to previous research results under normal wind conditions [22,23] the typhoon intensified the wind-induced responses of the tower body. The maximum mean radial displacement was increased by about 35%. Mean and RMS of blade responses were influenced slightly. It was suggested to strengthen the tower body in regions with frequent occurrences of typhoons.

Wind Fluttering Factor
The wind fluttering factor distributions of the tower body under different conditions when using internal stresses and displacement as the equivalent goals are shown in Figure 30. It was found that: (1) the wind fluttering factor of tower body under typhoon loads was positively related with height. Under different conditions, the wind fluttering factor when using radial displacement and the meridian axial force at 0° as the equivalent goal was mainly in the range of 2.0-4.0. The wind fluttering factor when using the meridian axial force at 0° as the equivalent goal was slightly larger than that when using the radial displacement as the equivalent goal. Distortion was developed at the tower top. (2) The wind fluttering factor was similar between Condition 2 and Conditions 8, as well as between the Condition 2 and Condition 8. Given two equivalent goals, the wind fluttering factor increased as the shielding effect weakened. The wind fluttering factor under Condition 4 and Condition 6 was significantly higher than those under rest conditions. (3) The wind fluttering factor distribution was gentle under the Condition 1, mainly concentrated in the 2.0-2.5 range. According to previous research results under normal wind conditions [22,23] the typhoon intensified the wind-induced responses of the tower body. The maximum mean radial displacement was increased by about 35%. Mean and RMS of blade responses were influenced slightly. It was suggested to strengthen the tower body in regions with frequent occurrences of typhoons.

Wind Fluttering Factor
The wind fluttering factor distributions of the tower body under different conditions when using internal stresses and displacement as the equivalent goals are shown in Figure 30. It was found that: (1) the wind fluttering factor of tower body under typhoon loads was positively related with height. Under different conditions, the wind fluttering factor when using radial displacement and the meridian axial force at 0 • as the equivalent goal was mainly in the range of 2.0-4.0. The wind fluttering factor when using the meridian axial force at 0 • as the equivalent goal was slightly larger than that when using the radial displacement as the equivalent goal. Distortion was developed at the tower top.
(2) The wind fluttering factor was similar between Condition 2 and Conditions 8, as well as between the Condition 2 and Condition 8. Given two equivalent goals, the wind fluttering factor increased as the shielding effect weakened. The wind fluttering factor under Condition 4 and Condition 6 was significantly higher than those under rest conditions. (3) The wind fluttering factor distribution was gentle under the Condition 1, mainly concentrated in the 2.0-2.5 range.
Distributions of wind fluttering factor of blades under different conditions are shown in Figures 31 and 32. Downwind displacement and shearing force as equivalent goals were used as the equivalent goal. It can be known by comparison that: (1) the wind fluttering factor of blade was significantly higher than that of the tower body. Given different equivalent goals, the wind fluttering factor mainly distributed in the range of 2.0-5.0. The wind fluttering factor when using downwind displacement as the equivalent factor was slightly larger than that when using shearing force as the equivalent factor. (2) By using the shearing force as the equivalent goal, the wind fluttering factor when there's a sharp increase of shearing force peak at blade root was relatively high. When the shearing force peak was relatively weak, the wind fluttering factor at blade tip was larger than that at blade root. Distributions of wind fluttering factor of blades under different conditions are shown in Figures 31 and 32. Downwind displacement and shearing force as equivalent goals were used as the equivalent goal. It can be known by comparison that: (1) the wind fluttering factor of blade was significantly higher than that of the tower body. Given different equivalent goals, the wind fluttering factor mainly distributed in the range of 2.0-5.0. The wind fluttering factor when using downwind displacement as the equivalent factor was slightly larger than that when using shearing force as the equivalent factor. (2) By using the shearing force as the equivalent goal, the wind fluttering factor when there's a sharp increase of shearing force peak at blade root was relatively high. When the shearing force peak was relatively weak, the wind fluttering factor at blade tip was larger than that at blade root. (3) The wind fluttering factor when using downwind displacement as the equivalent goal developed serious distortion, especially under blade shielding effects. In conditions without blade shielding effect, the distribution of wind fluttering factor was relatively gentle and the numerical value was in the range of 2.2-3.8.

Conclusions
In this paper, the aerodynamic performance of a large wind turbine under typhoon loads wasstudied by WRF and CFD nesting technology. The dynamic characteristics and wind-induced responses of the structural system under different stop positions of blades are discussed by

Conclusions
In this paper, the aerodynamic performance of a large wind turbine under typhoon loads wasstudied by WRF and CFD nesting technology. The dynamic characteristics and wind-induced responses of the structural system under different stop positions of blades are discussed by

Conclusions
In this paper, the aerodynamic performance of a large wind turbine under typhoon loads wasstudied by WRF and CFD nesting technology. The dynamic characteristics and wind-induced responses of the structural system under different stop positions of blades are discussed by combining the finite element method. A series of conclusions were drawn: (1) By comparing typhoon wind speed, wind intensity, wind direction and wind profile, it is concluded that the meso-scale WRF model can effectively simulate the near surface typhoon wind field. Based on the minimum least square method, the profile index of typhoon "Nuri" was fitted at 0.076. In this paper, the downscaling method can effectively simulate the 3D typhoon field of this kind of large-scale wind turbine system, which verifies the validity of the meso-and micro-scale nested simulation method, and provides the load input for the subsequent random wind pressure characteristics and wind vibration dynamic analysis. (2) The simulation results show that the wake zone is expanded gradually and the vorticity increment zone increases significantly under typhoon loads. Such a phenomenon is weakened with the reduction of the shielding area. The fluctuating and extreme wind pressures on blades and tower body are increase dramatically under typhoon loads. The maximum increase reaches 29%, which is appeared on the tower top of windward side. The fluctuating wind pressure is increased gradually as the blade shielding effect weakens. (3) Internal stresses and wind fluttering factor of blades and tower body are increased significantly under typhoon loads. When the tower body is not shielded by blades, the wind fluttering factor is decreased significantly. When the tower body is shielded by blades, the wind fluttering factor is increased significantly. The increase of radial displacement of the tower body is intensified as the shielding effect weakens. The maximum increase reaches 35%. (4) Based on our comprehensive analysis, when the large wind turbine is stopped under typhoon loads, the most unfavorable condition is when blade overlaps with the tower body completely (Condition 1). The safety redundancy reaches the maximum when the upper blade overlaps with the tower body completely (Condition 5).