## Article

# Reliability analysis of a large curved-roof structure considering wind and snow coupled effects

Correspondence to: Dr. Yongbo Peng, Shanghai Institute of Disaster Prevention and Relief, Tongji University, 1239 Siping Road, Shanghai 200092, China. E-mail:

**Received:**9 Apr 2022 |

**Revised:**22 Nov 2022 |

**Accepted:**5 Dec 2022 |

**Published:**20 Dec 2022

**Academic Editor:**Han-Ping Hong |

**Copy Editor:**Fangling Lan |

**Production Editor:**Fangling Lan

**Open Access**This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

## Abstract

Ensuring the safety of large curved-roof structures subjected to blizzards in cold regions remains challenging. This paper assesses the reliability of a large curved-roof structure considering the coupled effects of wind and snow. To this end, a two-way coupled simulation considering the interaction between snow particles and turbulent wind is first carried out using a computational fluid dynamics scheme and a spectral representation method based on the wavenumber-frequency joint power spectrum and stochastic harmonic function; thereby, the uneven snow distribution and stochastic wind field on a large curved roof are revealed. The probability density evolution method is then adopted to conduct a stochastic response and reliability analysis of the structure. The simulation results reveal that a credible uneven distribution of snow on the large curved roof can be obtained using the suggested snowdrift simulation scheme, and the refined spectral representation method can efficiently and accurately simulate the two-dimensional fluctuating wind field of the large curved-roof structure. It is shown that under stochastic wind, snow and self-weight loads for a 100-year return period, the large curved-roof structure has a higher risk of member yielding than of node deflection. Meanwhile, the axial stress acting on the structural member is far less than the maximum allowable stress and the structure is thus deemed to be sufficiently safe.

## Keywords

*,*two-way coupled simulation

*,*uneven snow distribution

*,*stochastic wind field

*,*probability density evolution method

*,*large curved-roof structure

## INTRODUCTION

Urbanization has received much attention from researchers worldwide in recent decades. With an increasing urban population, there has been considerable development of non-residential buildings. Many public buildings have been constructed with a large spatial configuration; e.g., industrial plants, exhibition halls, airport terminals, indoor gymnasiums and concert halls.

Compared with traditional buildings having normal spans, modern large-span spatial structures have more spacious and more open interior spaces and provide more possibilities for indoor activities. However, the increasing structural span means that the design of large-span spatial structures has to consider two issues. (ⅰ) The sensitivity of large structures to the gravity load is high, which is of particular concern for large-span spatial structures located in cold regions, where the snow load is a control load of the structural design ^{[1]}; and (ⅱ) Large structures have high structural flexibility and low natural vibration frequency and are thus prone to wind-induced damage under strong winds or winds with dominating low-frequency fluctuating components ^{[2]}. In fact, many instances of snow-induced collapse and wind-induced damage have been recorded for large-span spatial structures. For instance, in August 2005, the hurricane Katrina struck the city of New Orleans in the southern United States (with an average wind speed of 40.2 m/s), tearing and blowing over the ethylene propylene diene monomer rubber roof of the Louisiana Superdome. In April 2016, winds gusting at magnitude 13 (with an average wind speed of 39.9 m/s) hit the Changzhou area in Guangzhou, China, overturning the iron roof of the Guangwai Gymnasium in University Town. In January 2008, rarely heavy snowfall fell in Shanghai, China, and damaged and collapsed many houses, parking garages and warehouses. In February 2006, a blizzard hit Moscow, Russia, and snow accumulated continually on the roof of the Basmany Market on the outskirts of the city over a period of days and finally collapsed the roof, killing more than 50 people. In December 2010, strong winds and blizzards struck Minneapolis in the United States, resulting in the accumulation of more than 61 cm of snow on the roof of the home stadium of the Vikings football team, eventually destroying the roof.

The above-mentioned structural failure events suggest that there has been insufficient research on wind-induced damage, snow-induced damage and their coupled effects on large-span spatial structures. It is thus important to structural engineering practice to determine the global reliability of large-span spatial structures considering the coupled effects of wind and snow. Relevant studies on the distribution of snow and wind loads on a roof ^{[3]}, statistical analysis of structural responses ^{[4]} and reliability assessment ^{[5]} have been conducted. However, although the associated work has made great progress, there remain critical issues to address. (ⅰ) Most numerical simulations of the snow distribution on building roofs have been carried out on small models and only occasionally calibrated with the results of wind tunnel tests or field measurements ^{[6, 7]}, and the numerical simulation of prototype models of building roofs has not received enough attention; (ⅱ) Previous wind-induced snowdrift models are mostly based on the assumption of one-way coupled effects ^{[3, 8]}, ignoring the reaction of moving snow particles to turbulent wind and its further effect on snow particle motion, and the simulation of snow- and wind-load distributions considering the two-way coupled effects of snow particles and turbulent wind needs to be addressed; and (ⅲ) The vibration response analysis of large curved-roof structures subjected to snow and wind loads in previous studies has been largely based on the equivalent static wind load or frequency-domain analysis ^{[9, 10]}, the argument for assessing the structural performance has been a deterministic index or second-order statistical moment, and the accurate reliability of structures has not been secured owing to a lack of sufficient probabilistic information of the responses.

The present study thus conducts the reliability assessment of a large curved-roof structure through time-history analysis considering the coupled effects of wind and snow. As a practical example of structural design, a large curved-roof structure located on the outskirts of a city in a cold region of China considering the design return period of 100 years is used. The contributions of the study are as follows.

• Wind and snow pressure coefficients of a large curved-roof structure are determined through computational fluid dynamics simulations considering two-way coupled effects of the snow particles and turbulent wind.

• A two-dimensional fluctuating wind field for time history analysis of a large curved-roof structure is simulated by adopting a refined spectral representation method combining the wavenumber-frequency joint power spectrum and stochastic harmonic function.

• A stochastic response analysis and performance evaluation of a large curved-roof structure subjected to an uneven snow distribution and fluctuating wind field are carried out using the probability density evolution method.

• The reliability of the large curved-roof structure under stochastic wind, snow and self-weight loads with a 100-year return period is assessed.

The rest of this paper is organized as follows. Section 2 simulates the wind pressure and snow pressure coefficients by introducing the two-way coupling of snow particles and turbulent wind and also simulates the fluctuating wind field using a spectral representation method based on the wavenumber-frequency joint power spectrum and stochastic harmonic function. Section 3 presents the results of stochastic response analysis of the structure based on the probability density evolution method (PDEM). Section 4 conducts a dynamic reliability assessment of the large curved-roof structure. Concluding remarks are given in Section 5.

## SIMULATION OF THE SNOWDRIFT DISTRIBUTION AND FLUCTUATING WIND FIELD

### Overview of a large curved-roof structure

The considered large curved-roof structure serves as a terminal on the outskirts of a city in a cold region of China. A Google map of the terminal and its adjacent buildings is shown in Figure 1A. The main section of the terminal has a longitudinal span of 204 m (with a longitudinal dimension of 216 m), a transverse span of 75 m and a maximum overhang width of 7.5 m (with a transverse dimension of 90 m), which represents a typical large-span spatial grid structure; see the plan of the terminal in Figure 1B. The building has three stories, with the first and second stories having a reinforced-concrete frame structure and a height of 7.8 and 8.2 m, respectively, whereas the third story is a steel frame structure with a bidirectional curved and double-layer reticulated shell roof, with the floor height being 21 m; see the elevation of the terminal in Figure 1C. The roof comprises plane trusses, restraint trusses, steel pipe columns and steel pipe diagonal braces and varies in elevation by approximately 5 m. All connections of the steel pipe columns and steel pipe diagonal braces with the bottom concrete frame column and the top steel roof are cast steel nodes, which are approximately hinged. The region of the terminal has a temperate continental climate, including a long cold winter with severe wind and snow conditions. According to *Chinese Load Code for the Design of Building Structures* (GB50009-2012), the average wind pressure and average snow pressure for a 100-year return period are 0.7 and 1.0 kN/m

### Simulation of wind pressure and snow pressure coefficients

Wind-induced snow drift is a phenomenon of interactions between the near-ground wind and snow particles at the surface of ground snow and floating in the air. The classical scheme of wind-induced snow drift simulation thus includes the simulation of the near-ground wind field and the simulation of snow drift motion.

The Reynolds averaged Navier-Stokes (RANS) equations are often used to model the near-ground wind field by considering flow properties such as the wind velocity, wind pressure, flow density, wall shear stress and turbulent characteristics. To consider the effect of snow particle motion on the air turbulence characteristics and the air phase-averaged momentum, an extended two-way coupled scheme for wind-induced snow drift simulation is used. In this scheme, different damping source terms are introduced into the RANS equations and the transport equations of the turbulent kinetic energy and energy dissipation rate of the standard ^{[11]}:

where ^{[12]}. ^{[11]}:

Here,

The simulated snowdrift is often quantified as the flux of snow erosion/deposition, which is mainly contributed by snow particles in both the saltation and suspension layers ^{[13]}:

The flux of snow erosion and the flux of snow deposition are expressed as ^{[14, 15]}

where ^{[16]}; ^{[13]}:

where

In representing the uneven distribution of wind and snow loads, the wind pressure and snow pressure coefficients are defined as ^{[17, 18]}

where

The finite volume method is applied to solve the partial differential equations in Equations 1-3 and 10 using the finite element software ANSYS 14.0 FLUENT. The airshed modeling and grid division of the flow domain of the large curved roof subject to wind blowing in the most unfavorable direction (i.e., along the transverse span of the terminal building because this span is just one-third of the longitudinal span) are shown in Figure 2. The distances from the centroid of the horizontal projection plane of the terminal building to the wind inlet, outlet and top surface are set at 10

Figure 2. Modeling and mesh of the flow domain of the large curved roof subject to wind blowing in the most unfavorable direction. (A) Modeling of the flow domain; (B) mesh of the flow domain.

Figure 3 presents the coefficients of wind pressure and snow pressure on the large curved roof. For ease of analysis, the roof is divided into 18 sub-regions of the same size (36 m

Figure 3. Coefficients of wind pressure and snow pressure on the large curved roof. (A) Wind pressure coefficient; (B) snow pressure coefficient.

Figure 4. Distribution of coefficients of wind pressure and snow pressure in sub-regions of the large curved roof. For ease of analysis, the roof is divided into 18 sub-regions of the same size (36 m

The two-way coupled simulation and one-way coupled simulation were compared in the authors' previous study ^{[19]}. That study revealed that the two-way coupled simulation was more accurate relative to field observations (having 7.9% lower simulation error) and more efficient (having 37.1% lower computational cost) than the one-way coupled simulation.

### Simulation of the fluctuating wind field

The wind pressure coefficient for the roof shown in Figure 4A is attained in steady computational fluid dynamics (CFD) simulation, and it is essentially an average value of the real wind pressure field. In reality, carrying out a transient CFD simulation is time-consuming. A practical approach is to decompose the wind field into a deterministic average wind component and a stochastic fluctuating wind component through Reynolds decomposition. The deterministic average wind component can be simulated using the RANS equations (see Equations 1-3), whereas the simulation of the stochastic fluctuating wind component often resorts to a spectral representation.

The elevation variation of the roof is far less (by two orders of magnitude) than the plan size, and the three-dimensional spatial fluctuating wind field can thus be approximately presented as a two-dimensional fluctuating wind field. Moreover, to reduce the computational cost, the same sub-regions of the roof used for the wind pressure coefficient distribution are adopted, where the center point is taken as the representative point and other points in the same sub-region are assumed to have the same fluctuating wind as the representative point. Considering the low efficiency of the conventional scheme of spectral representation methods in wind field simulation, the spectral representation method based on the wavenumber-frequency joint power spectrum and stochastic harmonic function is used in this study ^{[20, 21]}:

where ^{[21]}

where

This study uses the Davenport spectrum and coherence function model to derive the wavenumber-frequency joint power spectrum as ^{[21]}

where ^{[22]}.

To validate the accuracy of the wavenumber-frequency joint power spectrum and stochastic harmonic function based spectral representation method, the simulated fluctuating wind velocities at three representative points on the roof, namely the center points of the sub-regions Z9 (labeled P1), Z10 (labeled P2) and Z4 (labeled P3), are considered. Representative time histories of the fluctuating wind velocities at P1, P2 and P3, and their auto-power spectra and coherence functions are shown in Figures 5 and 6. It is seen that the auto-power spectra of the simulated fluctuating wind velocities at the different spatial points and their coherence functions are in good agreement with the targets; see the Davenport spectrum and the related coherence function, which indicates that the wavenumber-frequency joint power spectrum and stochastic harmonic function based spectral representation method can be used to accurately simulate the two-dimensional fluctuating wind field for a large curved-roof structure.

Figure 5. Representative time histories of fluctuating wind velocities at P1, P2 and P3. For ease of analysis, the roof is divided into 18 sub-regions (6 sub-regions

Figure 6. Auto-power spectra and coherence functions of fluctuating wind velocities at P1, P2 and P3. For ease of analysis, the roof is divided into 18 sub-regions (6 sub-regions

## STOCHASTIC RESPONSE ANALYSIS OF A LARGE CURVED-ROOF STRUCTURE

### PDEM

There is randomness associated with the fluctuating wind field, and the dynamical system of the large curved-roof structure subjected to the coupled effects of the wind and snow is thus a typical random vibration problem, which needs to be solved using probabilistic methods. In this study, the PDEM, which has been widely adopted in recent years, is used to efficiently and accurately estimate the statistical information of the structural response ^{[23]}.

Without loss of generality, we write the equation of motion of an

where

As regards the physical quantities ^{[23]}

Specifically, as

where

It is noted that the GDEE reveals the intrinsic connections between a stochastic dynamical system and its deterministic counterpart. To solve the GDEE, an initial condition is given by

where

Solving Equation 25 under the initial condition given as Equation 26, the instantaneous probability density function (PDF) of

where

The numerical procedure for solving the GDEE is as follows.

*Step 1*: Select the representative point set P @^{[24]}.

*Step 2*: For the prescribed

*Step 3*: Introduce

*Step 4*: Obtain the numerical result of

### Finite element modeling and analysis

Using the finite element analysis platform SAP2000 v19, the deterministic structural dynamic analysis involved in *Step 2* in Section 3.1 is carried out. In the concrete frame-steel roof hybrid structure, the lower concrete structure has higher rigidity than the top large-span steel roof structure, and the structural internal force and displacement are generally small under wind and snow loads. To simplify the calculation, therefore, only the steel roof at the top of the terminal building is modeled using finite elements, and the bottom of the steel pipe column and steel diagonal bracing adopts hinge supports, of which the number is 58. For the steel pipe columns, steel diagonal bracings, upper chords, lower chords and diagonal rods of the steel roof, frame elements are used, and the section size of each element is selected according to the actual design size. The steel roof support shell adopts massless shell elements for roof load transfer. By doing so, the finite element model of the large curved-roof structure includes a total of 4241 nodes, 16, 073 frame elements and 3083 massless shell elements, with a total structural mass of 2924.6 t. The three-dimensional finite element model is shown in Figure 7A.

Figure 7. Three-dimensional finite element model and the first mode of the structure. (A) View of the finite element model; (B) first mode of the structure.

To clarify the basic dynamic characteristics of the large curved-roof structure, structural modal analysis is carried out using the Rayleigh-Ritz vector method. It is found that the first four orders of the modal frequencies are 0.854, 0.767, 0.553 and 0.463 s, and they sequentially correspond to the modal characteristics of horizontal vibration along the *Chinese Technical Regulation for Spatial Grid Structures* (JGJ7-2010); therefore, the wind-induced vibration problem should be considered in the structural analysis. Additionally, the natural vibration periods of the structure are distributed in a narrow range; i.e., the structure has densely distributed frequencies. Thereby, the contribution of each mode to the dynamic response of the structure is equally important, and ignoring the contribution of high-order modes will inevitably lead to large computational errors. Moreover, the overall vertical stiffness of the structure is high, and it is predicted that the vertical deflection of nodes of the roof under the loads will be smaller than that for a flexible roof of the same size.

As a practical situation, the large curved-roof structure is simultaneously subjected to the wind load comprising average and fluctuating wind components, the snow load unevenly distributed on the roof and the self-weight load. Structural analysis should consider a certain combination of the effects of these loads. In this study, they are in weighted linear combination according to *Chinese Load Code for the Design of Building Structures* (GB50009-2012). Figure 8A and B show the steady member axial stress and node vertical deflection of the large curved-roof structure under representative wind, snow and self-weight loads (at a time instant of 300 s).

Figure 8. Steady member axial stress and node deflection of the large curved-roof structure under representative wind, snow and self-weight loads (at a time instant of 300 s). (A) Member axial stress of the structure (MPa); (B) node vertical deflection of the structure (mm).

It is seen that under the representative wind, snow and self-weight loads, the axial stress acting on members is in the range of -198 to 154 MPa, and the lower roof chord and roof inclined rod experience higher stresses as same as the range in the structural level, whereas the upper roof chord and roof support column experience lower stresses in the range of -141 to 106 MPa. Under the representative wind, snow and self-weight loads, the node vertical deflection is in the range of -71.5 to 16.5 mm and the deflection is a maximum at the nodes in the mid-span area of the roof, whereas the overhanging part of the roof has a small upward deflection. Therefore, the member axial stress of the large curved-roof structure reaches a maximum of 198 MPa, and the node vertical deflection reaches a maximum of 71.5 mm. Considering that the strength of the structural member is 345 MPa, the maximum member axial stress is approximately 0.574 times the maximum allowable stress. The threshold of the node vertical deflection is 1/250 of the horizontal span of the structure according to the design provision (i.e., 300 mm), and the maximum node vertical deflection is approximately 0.238 times the threshold. The member axial stress is closer to its limit value than the node vertical deflection, and the risk of the member yielding is thus higher than that of node deflection, which is in agreement with the results of structural modal analysis.

### Statistical analysis of the stochastic responses of structures

In this study, there is a random vibration problem due to the randomness associated with the fluctuating wind field. Using the PDEM, a total of 300 samples of the fluctuating wind field are simulated and used for the stochastic response analysis and reliability assessment of the large curved-roof structure. For the purpose of illustration, the stochastic responses of sub-region Z10 are addressed. Time histories of the means and standard deviations of the axial stress of Member #14294 and the vertical deflection of Node #3019 in this sub-region are presented in Figure 9A-D, respectively. The simulated fluctuating wind field is a stationary state, and the time histories of statistical moments within just the first 300 s are thus presented to clearly show the initial non-stationary stage of the structural responses.

Figure 9. Time histories of the means and standard deviations of the axial stress acting on Member #14294 and vertical deflection of Node #3019 in sub-region Z10. For ease of analysis, the roof is divided into 18 sub-regions (6 sub-regions

It is seen that the means and standard deviations of the member axial stress and node vertical deflection quickly reach a steady state; i.e., within 50 periods of the basic natural vibration of the structure. Additionally, the average wind component of the stochastic wind load has a suction effect on the roof, which partially offsets the structural response relating to the self-weight load and the unevenly distributed snow load, and the average wind component thus reduces the overall response of the structure. Meanwhile, the fluctuating wind component has an appreciable dynamic effect on the responses of the structure. Furthermore, the probability densities of the axial stress of Member #14294 at 100, 200, and 300 s and in the time interval [150, 160] s are presented in Figure 10. It is seen that the probability density of the structural response in the stage of the steady state only slightly changes and approximately follows a normal distribution. This can be explained in that the simulation of the stochastic wind field is a Gaussian process and the dynamic analysis using finite element software reveals the linear state of the structure.

Figure 10. Probability densities of the axial stress of Member #14294. (A) Probability densities at typical instants; (B) probability densities in the time interval [150, 160] s.

Figure 10 also shows that although the mean of the member axial stress is far less than the maximum allowable stress, the member still faces the possibility of yielding because the axial stress varies appreciably. It is thus necessary to carry out reliability assessment in terms of the member axial stress in determining the safety of the structure.

## DYNAMIC RELIABILITY ANALYSIS OF A LARGE CURVED-ROOF STRUCTURE

### PDEM-based reliability method

Without loss of generality, the associated reliability of the first-passage problem of structural dynamical systems is defined as

where Pr

When the structural system has multiple modalities of failure, the first-passage problem is expressed as

where

We construct an equivalent extreme-value event as a random variable:

The reliability function in Equation 29 can then be defined alternatively as

Equation (32) can be rewritten as

where

Therefore, the introduction of the equivalent extreme-value event criterion (see Equation 31) reduces the high-dimension integral associated with the reliability of systems into the one-dimension integral ^{[25]}.

To derive the PDF of the equivalent extreme-value event, a pseudo-random process is first defined as

which satisfies the conditions

where

Theoretically, arbitrary forms of the pseudo-random process

It is observed that the probabilistic information of the equivalent extreme-value event

### Reliability assessment of a structure

Using the PDEM and the associated reliability method, the PDFs and cumulative density functions (CDFs) of the equivalent extreme-value axial stress of members in the sub-regions Z1, Z8, Z10 and Z17 are derived, as shown in Figure 11. It is seen that all the curves related to the sub-regions under investigation are far below the maximum allowable stress (i.e., 345 MPa), indicating safety in these regions. It is also seen that among the four representative sub-regions under investigation, the sub-region Z10 has the largest member axial stress, followed by Z1, Z8 and Z17 in order. In view of Figure 4, the sub-region Z1 features lower wind suction and higher snow pressure and is thus expected to have appreciable member axial stress, whereas the other three sub-regions feature almost similar wind suction and snow pressure and thus, compared with Z8 and Z17, the sub-region Z10 has a greater contribution from the fluctuating wind component to the member axial stress because this area is located in the rolling region of the curved roof. Furthermore, all the sub-regions are analyzed and their reliabilities are shown in Figure 12. Additionally, the PDF and CDF of the equivalent extreme-value axial stress of members of the roof structure are presented in Figure 13, Showing a safe structure in the global sense.

Figure 11. PDFs and CDFs of the axial stress acting on members in sub-regions Z1, Z8, Z10 and Z17. For ease of analysis, the roof is divided into 18 sub-regions (6 sub-regions

Figure 12. Reliabilities of axial stress of members in sub-regions of the roof structure. For ease of analysis, the roof is divided into 18 sub-regions of the same size (36 m

Figure 13. PDF and CDF of equivalent extreme-value axial stress of members of the roof structure. (A) PDF of equivalent extreme-value axial stress; (B) CDF of equivalent extreme-value axial stress. PDF: Probability density function; CDF: cumulative density function.

In summary, under stochastic wind, snow and self-weight loads with a 100-year return period, the equivalent extreme-value axial stresses of members in each sub-region and in the global structural system are both far less than the maximum allowable stress, and the large curved roof structure is deemed to be sufficiently safe.

## CONCLUSIONS

Aiming at the reliability assessment of large curved-roof structures considering the coupled effects of wind and snow with a 100-year return period, this study conducted pertinent works, including a snowdrift simulation using computational fluid dynamics, a fluctuating wind field simulation using the wavenumber-frequency joint power spectrum and stochastic harmonic function based spectral representation method, and a stochastic response and dynamic reliability analysis of a structure using the probability density evolution method. The main remarks and conclusions of the study are listed as follows.

• Snowdrift simulation considering the two-way coupling of snow particles and turbulent wind provides a credible uneven distribution of snow on a large curved-roof structure.

• The wavenumber-frequency joint power spectrum and stochastic harmonic function based spectral representation method can be used to efficiently and accurately simulate the two-dimensional fluctuating wind field of a large curved-roof structure.

• The overall vertical stiffness of the structure is large, and the vertical deflection of nodes of the presented roof under loads is thus smaller than that of a flexible roof of the same size, and there is a higher risk of structural member yielding than structural node deflection.

• The average wind component of the stochastic wind load has a suction effect on the roof, which partially offsets the structural response for the self-weight load and the unevenly-distributed snow load and reduces the overall response of the structure. Meanwhile, the fluctuating wind component has an appreciable dynamic effect on the responses of the structure.

• The probability density of the structural response in the stage of a steady state changes only slightly and approximately follows a normal distribution because the simulation of the stochastic wind field is a Gaussian process and the dynamic analysis using finite element software reveals the linear state of the structure.

• Under stochastic wind, snow and self-weight loads having a 100-year return period, the member axial stress of the structural system is far less than the maximum allowable stress, and the large curved-roof structure is thus deemed sufficiently safe.

This study simulated snowdrift considering the two-way coupling of snow particles and turbulent wind. However, only steady simulations based on non-thermal and constant meteorological conditions were involved, and the simulation results may differ from reality. Therefore, unsteady simulations considering thermal effects such as roof heating and varying meteorological conditions including the wind velocity, wind direction, temperature and humidity need to be conducted in the near future. Additionally, a field measurement of the snow distribution on a large curved roof will be performed to validate the simulation results of the present study.

## DECLARATIONS

### Acknowledgments

The support of the National Key R & D Program of China (Grant No. 2017YFC0803300) is highly appreciated.

### Authors' contributions

Conceptualization, methodology, validation, formal analysis, writing - review & editing, resources, supervision, project administration: Peng Y

Software, investigation, data curation: Zhao W, Zhou J

Methodology, writing - original draft, visualization: Zhao W

### Availability of data and materials

Some or all data and materials that support the findings of this study are available from the corresponding author upon reasonable request.

### Financial support and sponsorship

National Key R & D Program of China (Grant No. 2017YFC0803300).

### Conflicts of interest

All authors declared that there are no conflicts of interest.

### Ethical approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Copyright

© The Author(s) 2022.

## REFERENCES

1. Vazna RV, Zarrin M. Sensitivity analysis of double layer Diamatic dome space structure collapse behavior. *Eng Struct* 2020;212:110511.

2. Wang XJ, Li QS, Yan BW, Li JC. Field measurements of wind effects on a low-rise building with roof overhang during typhoons. *J Wind Eng Ind Aerodyn* 2018;176:143-57.

3. Zhou XY, Li XF. Simulation of snow drifting on roof surface of terminal building of an airport. Available from: https://www.researchgate.net/publication/285439767_Simulation_of_Snow_Drifting_on_Roof_Surface_of_Terminal_Building_of_an_Airport . [Last accessed on 5 Dec 2022].

4. Zhang GL, Zhang QW, Fan F, Shen SZ. Numerical simulations of snowdrift characteristics on multi-span arch roofs. *J Wind Eng Ind Aerodyn* 2021;212:104593.

5. Aghamohammadi R, Nasrollahzadeh K, Mofidi A, Gosling P. Reliability-based assessment of bond strength models for near-surface mounted FRP bars and strips to concrete. *Compos Struct* 2021;272:114132.

6. Zhang GL, Zhang QW, Fan F, Shen SZ. Numerical simulations of development of snowdrifts on long-span spherical roofs. *Cold Reg Sci Technol* 2021;182:103211.

7. Zhang GL, Zhang QW, Fan F, Shen SZ. Field measurements of snowdrift characteristics on reduced scale building roofs based on the size effect study. *Structures* 2021;32:2020-31.

8. Liu DT, Wang B, Li YL, Liu S. A source term model for drifting snow based on the assumption of local equilibrium saltation. *Cold Reg Sci Technol* 2021;181:103175.

9. Yang QS, Chen B, Wu Y, Tamura Y. Wind-induced response and equivalent static wind load of large-size curved-roof structures by combined Ritz-proper orthogonal decomposition method. *J Struct Eng* 2013;139:997-1008.

10. Uematsu Y, Watanabe K, Sasaki A, Yamada M, Hongo T. Wind-induced dynamic response and resultant load estimation of a circular flat roof. *J Wind Eng Ind Aerodyn* 1999;83:251-61.

11. Okaze T, Takano Y, Mochida A, et al. Development of a new k-epsilon model to reproduce the aerodynamic effects of snow particles on a flow field. *J Wind Eng Ind Aerodyn* 2015;144:118-24.

12. Alfonsi G. Reynolds-Averaged Navier-Stokes equations for turbulence modeling. *Appl Mechan Rev* 2009;62:040802.

13. Uematsu T, Nakata T, Takeuchi K, et al. Three-dimensional numerical simulation of snowdrift. *Cold Reg Sci Technol* 1991;20:65-73.

14. Anderson RS, Haff PK. Wind modification and bed response during saltation of sand in air. *Acta Mechan* 1991:21-51.

15. Naaim M, Naaim-Bouvet F, Martinez H. Numerical simulation of drifting snow: erosion and deposition models. *Ann Glaciol* 1998;26:191-6.

16. Kind RJ. Snow drifting. In Gray DM, Male DH, editors. Handbook of Snow, principles, processes, management, and use. Oxford: Pergamon Press; 1981, pp. 338-59.

17. Zhou XY, Zhang Y, Kang L, et al. CFD simulation of snow redistribution on gable roofs: impact of roof slope. *J Wind Eng Ind Aerodyn* 2019;185:16-32.

18. Sun X, He R, Wu Y. Numerical simulation of snowdrift on a membrane roof and the mechanical performance under snow loads. *Cold Reg Sci Technol* 2018;150:15-24.

19. Peng YB, Zhao WJ, Li S, Zhou J. Prediction of wind-induced snow redistribution on large-span building roof based on two-way coupled simulation. *Nat Hazards Rev* 2022. under review.

20. Song YP, Chen JB, Peng YB Spanos PD, Li J. Simulation of nonhomogeneous fluctuating wind speed field in two-spatial dimensions via an evolutionary wavenumber-frequency joint power spectrum. *J Wind Eng Ind Aerodyn* 2018;179:250-9.

21. Chen JB, Song YP, Peng YB, et al. Simulation of homogeneous fluctuating wind field in two spatial dimensions via a joint wave number-frequency power spectrum. *J Eng Mechan* 2018;144:04018100.

22. Peng YB, Wang SF, Li J. Field measurement and investigation of spatial coherence for near-surface strong winds in Southeast China. *J Wind Eng Ind Aerodyn* 2018;172:423-40.

23. Li J, Chen JB. Stochastic dynamics of structures. Singapore: John Wiley & Sons; 2009.

24. Chen JB, Ghanem R, Li J. Partition of the probability space in probability density evolution analysis of non-linear stochastic structures. * Probabilistic Eng Mech* 2009;24:27-42.

## Cite This Article

**OAE Style**

Peng Y, Zhao W, Zhou J. Reliability analysis of a large curved-roof structure considering wind and snow coupled effects. *Dis Prev Res* 2022;1:8. http://dx.doi.org/10.20517/dpr.2022.02

**AMA Style**

Peng Y, Zhao W, Zhou J. Reliability analysis of a large curved-roof structure considering wind and snow coupled effects. *Disaster Prevention and Resilience*. 2022; 1(2): 8. http://dx.doi.org/10.20517/dpr.2022.02

**Chicago/Turabian Style**

Peng, Yongbo, Weijie Zhao, Jian Zhou. 2022. "Reliability analysis of a large curved-roof structure considering wind and snow coupled effects" *Disaster Prevention and Resilience*. 1, no.2: 8. http://dx.doi.org/10.20517/dpr.2022.02

**ACS Style**

Peng, Y.; Zhao W.; Zhou J. Reliability analysis of a large curved-roof structure considering wind and snow coupled effects. *Dis. Prev. Res.* **2022**, *1*, 8. http://dx.doi.org/10.20517/dpr.2022.02

## Related

## Comments

Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at support@oaepublish.com.

**12**clicks

**0**comments

**9**likes

*|*risk assessment

*|*resilience

*|*

^{0}