Next Article in Journal
A Study on the Influence of Suspended Matter by the Foundation Construction of Different Offshore Wind Turbines in the Taiwan Sea Area
Next Article in Special Issue
Control of Runoff Peak Flow for Urban Flooding Mitigation
Previous Article in Journal
Influence of Floodplain Flooding on Channel Siltation Adjustment under the Effect of Vegetation on a Meandering Riverine Beach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Seismic Fragility Evaluation of Simply Supported Aqueduct Accounting for Water Stop’s Leakage Risk

College of Water Resources and Architectural Engineering, Northwest A&F University, Xianyang 712100, China
*
Authors to whom correspondence should be addressed.
Water 2021, 13(10), 1404; https://doi.org/10.3390/w13101404
Submission received: 16 April 2021 / Revised: 9 May 2021 / Accepted: 11 May 2021 / Published: 18 May 2021
(This article belongs to the Special Issue Frontiers in Urban Water Infrastructure)

Abstract

:
Due to the demands of booming Chinese cities and the increase in urban residents, the safety of aqueduct water transportation structures is noteworthy. A lot of old aqueducts were built in the 1990s and even earlier in the last century and may become vulnerable to potential earthquakes. This paper deals with an evaluation of an aqueduct’s seismic vulnerability accounting for leakage risk. Based on the Hua Shigou aqueduct in Ningxia, a probabilistic investigation was carried out to obtain the seismic fragility using Latin hypercube sampling. In the numerical study, the superstructure and substructure of the aqueduct were modeled as beam elements, and the lumped mass method was adopted to simulate the fluid–structure interaction. The rubber water stop’s mechanical performance was studied, and its damage states were proposed. Parametric numerical models were then subjected to a set of ground motions according to incremental dynamic analysis (IDA), which contained probabilistic parameters such as water, concrete strength, and bearing performance degradation. Both the system and component levels of the old aqueduct’s seismic fragility curves were obtained. It was found that the probability of the water stop’s leakage risk is significantly elevated with the increase in ground motion.

1. Introduction

The South-to-North Water Diversion Project is a national strategic project in China [1]. Its final water transfer scale is estimated to be 44.8 billion m3, which will greatly alleviate the imbalance of water resources between the southern and western areas of China. The Middle Route of the South-to-North Water Diversion Project (MRP), part of the South-to-North Water Diversion Project, is the largest water conveyance project in China and 1432 km in length to transfer [2]. A total of 27 large aqueducts have been constructed in the MRP to cross rivers and valleys.
However, performance-based seismic theory was not considered in the design stage of the aqueduct back in 2000, and it may be more vulnerable to earthquakes. The structure of an aqueduct is very similar to that of a bridge [3]. Aqueducts and bridges have similar super- and substructures. However, compared with a bridge, the coupling effect of fluid in an aqueduct should be considered [4]. Therefore, fluid–structure interaction (FSI) should be taken into account. In terms of the seismic fragility assessment methods of structures in a water tower and liquid storage tank such as the FSI effect, generally speaking, there are two main analysis methods. The first method is the lumped mass method, which is also called the added mass method. In 1957, Housner [5,6] studied the interaction between shaking water and a water tank and proposed the “equivalent concentrated method.” He simplified the mass of fluid into two points. The ACI standard [7] was developed based on this study. The lumped mass method was adopted to study dynamic buckling, seismic protection, and the response of elastic soil under horizontal ground excitation of a liquid storage tank in fluid–structure coupling studies [8,9,10,11]. However, this method is mainly used in rigid tanks and cannot evaluate aqueduct sidewalls. The second method is the nonlinear analysis including coupled Lagrangian–Eulerian (CLE), arbitrary Lagrangian–Eulerian (ALE), Smoothed Particle Hydrodynamics (SPH), and so on. D’Alessandro et al. [12] discussed the linear sloshing theory and nonlinear sloshing approach when studying the fluid sloshing coupling simulation of tank vehicle dynamics modeling. Rawat et al. [13] used coupled structural–acoustic (CSA) analysis and CLE to compare the response of water sloshing and hydrodynamic pressure in the two different analysis methods. Zhao et al. [14] applied some simplifications. Masses and shell elements were used to simulate the internal components and equipment of NPP due to its complicated shapes and arrangements. When studying the seismic performance and optimal design of a water tank with an inner ring baffle in a nuclear power plant under seismic load, water and air were simulated by the Eulerian element, while in a shield building, the water tank and ring baffles were simulated by the shell element. Ozdemir et al. [15] used ALE to simulate nonlinear fluid–structure coupling to investigate the application of anchored and unanchored tanks in seismic analysis. In general, the lumped mass method pays more attention to the overall response, while the nonlinear analysis places more emphasis on the local response of the liquid storage tank.
Seismic fragility methods have been widely used in bridge engineering [16]. Moreover, there have been a large number of studies on the seismic fragility of liquid storage tanks, industrial storage pipes, and silos [17,18,19,20,21,22]. In terms of aqueducts, Valeti et al. [4] studied the seismic response of an elevated aqueduct with a hydrodynamic force and the nonlinear interaction of soil–structure, and the effect of convective masses was found to be significant. Prashar [23] conducted an aqueduct earthquake disaster assessment. Ali Rafiee et al. [24] designed a numerical simulation of a masonry aqueduct under earthquake conditions. Liu et al. [25] studied the shock absorption and seismic effect of rubber support for large aqueducts.
As mentioned above, there has been little investigation on aqueducts’ seismic fragility. Because aqueducts are used to transport water, leakage risk should be considered. However, to the author’s knowledge, limited research about water stops has been reported. Considering that aqueducts and bridges are similar in structure, this paper carried out an analysis of aqueduct seismic fragility by referring to bridges. The fluid–structure coupling effect in the aqueduct was involved in the analysis model, and the limit state of the water stop deformation was proposed and considered in the seismic vulnerability analysis.

2. The Steps of Seismic Fragility Analysis

(1)
The appropriate numerical model was built.
(2)
According to the fortification intensity and site type of the selected project, the design response spectrum was obtained, which in this study is under Chinese Code GB_51247-2018 [26]. The first nine groups of ground motions were then selected according to the designed response spectrum from the PEER Ground Motion Database [27].
(3)
The Latin hypercube sampling (LHS) method was implemented to reduce the correlation between variables. This paper used LHS to extract six groups of three variables.
(4)
By running two groups of data analysis and making a comparison of the results, the most vulnerable ground motion direction was identified, and the transverse direction was found to be in the governing ground motion direction.
(5)
Each group of models was analyzed with nine groups of ground motion, and the obtained data were extracted and processed.
(6)
The damage index of the aqueduct pier and the water stop was proposed and determined.
(7)
The fragility curve was drawn, and then the possibility of the exceedance of a limit state under a specific intensity measure was discussed.

3. Numerical Modeling

3.1. Introduction of the Aqueduct

The water supply project of the north main canal is an important water resource optimal allocation project in Ningxia Province, of which the investigated Hua Shigou aqueduct was part. It is of great significance to improve the ecological environment at the eastern foot of the Helan Mountains, promote the development of characteristic agriculture, and stabilize the industrial and residential water supply in the west of Yinchuan. Hua Shigou aqueduct’s total length is 252 m and is divided into 21 spans. The length of each span was 12 m. This study selected a two-span model (shown in Figure 1 and Figure 2). The water transmission project consisted of three sections, among which the north main canal was 47.28 km in total length, including irrigated an area of 25,616.6 acres. Although freezing may occur in extremely low temperatures of Ningxia in winter, some studies have shown that the covered aqueduct will not freeze when the fluid depth is more than 2 m, which we did not cover in this paper.
The water stop is the waterproof structural component between two girders of aqueducts and mainly consists of two types: a rubber water stop and a plastic water stop. A rubber water stop was widely used in the aqueduct, and it can fit the expansion and deformation of concrete in the form of rubber material elasticity. According to structure deformation, water pressure, and waterproof grade, the mechanical property of the water stop was determined.
The geological conditions are slightly dense to medium dense breccia in the surface layer and argillaceous sandstone in the bottom layer. The strata is a nonliquefied soil layer.

3.2. Methods of Modeling and Analysis

3.2.1. Lumped Mass Method

FSI is simplified as the lumped mass model according to Housner’s theory [6]. The water filled in the aqueducts could be divided into two parts: the impulsive mass and convective mass (shown in Figure 3). For the rectangular tank with a width and water depth of 2 L and H, respectively, the tank with an oscillating water surface is shown in Figure 3a. The equivalent dynamic system is shown in Figure 3b. In Housner’s theory [6], the convective mass and impulsive mass are calculated as Equations (2) and (3). M is the total mass in the tank. The height of M0 and M1 are H0 and H1, respectively.
M = M 0 + M 1
M 0 = M tan h 1.7 L / H 1.7 L / H
M 1 = M 0.83 tan h 1.6 H / L 1.6 H / L
K 1 = 3 × M 1 2 M × gH L 2
In this paper, the connection between the impulsive mass and convective mass was simplified as a connection link with a transverse spring K1.

3.2.2. Modeling and Analysis

General FE software SAP2000 was used for incremental dynamic analysis (IDA). The superstructure of the aqueduct and the columns were discrete, as the frame element and pile–soil structural response was simplified as a spring element. The plate rubber bearing was simulated as a multilinear elastic link element for a nonlinear setting [28]. The water stop was simulated by a nonlinear gap element to capture its unique function. The water was considered as lumped mass, as stated previously. The nearby span weight was also modeled as a point load on the column to accurately reflect the aqueduct’s dynamic characteristic. The FE model, connecting links, and a model 3D view are illustrated in Figure 4.

4. Ground Motions and LHS

4.1. Selection of Ground Motions

The site condition of the aqueduct is II. The basic seismic peak ground acceleration (PGA) is 0.20 g. The characteristic period of the response spectrum of the basic ground motion acceleration is 0.45 s. According to the fortification intensity and site type of the project, the design response spectrum was obtained, as shown in Figure 5 [26].
The first nine groups of ground motions selected by the utility of the NGA-West2 On-Line Ground-Motion Database Tool are shown in Figure 6.
In each group, the ground motion with the maximum of PGA was selected and illustrated in Table 1.

4.2. Latin Hypercube Sampling

During the service of the aqueduct, the strength of concrete and the friction coefficient of the rubber bearings are affected by environmental conditions such as freeze-frozen temperature, ultraviolet light, etc. The amount of water is also a variable, mainly due to seasonal variation. In probabilistic theory, the variables are found to be conformed to certain probability distributions. The strength of C30 concrete (ƒ′c) in this study followed a normal distribution of mean = 20.1 MPa and standard deviation = 2.5 MPa [29]. The flow rate conformed to a normal distribution of flow, mean = 10.4 m3/s, and standard deviation = 2.6 m3/s [30]. The weight of water per segment was assumed to follow a normal distribution with a mean = 56.6 kN/m and a standard deviation = 14.1 kN/m. The friction coefficient of rubber bearing is a uniform distribution from 0.02 to 0.05 [31]. These variables are summarized in Table 2.
Insignificant differences between groups of the same variable would make the theory lack universality. To avoid it, the Latin hypercube sampling (LHS) method was implemented to reduce the correlation of variables between groups. In this paper, six groups of variables were obtained by using LHS, as shown in Table 3.
In the scaling approach, all the selected ground motions were incrementally scaled up and down to produce different levels of earthquake intensities. The IDA method was adopted to adjust the amplitude of ground motion Numbers 2–7. The PGA of nine selected groups of ground motion and the amplitude modulation (AM) acceleration are shown in Table 4.

5. Seismic Fragility Evaluation of Aqueducts

5.1. Definition of the Water Stop’s Limit State

Water stops are key components of aqueducts. The leakage of the water stop during work will not only affect normal water delivery but also accelerate the corrosion of reinforcement and soften the foundation soil. This will reduce the durability and stability of the structure and seriously endanger the safety of the structure. Especially in the case of a lagged inspection, an earthquake may lead to the aqueduct overturning because the soil around the pier foundation becomes soft. Therefore, it is very necessary to study the vulnerability of a water stop under the action of an earthquake.
Under the action of a transverse earthquake, the shear deformation of the water stop is illustrated in Figure 7.
The damage index is defined as Equation (5):
tan γ = Δ L
where γ is the rotation of shear deformation. Δ is the aqueduct transverse relative displacement. L is the width of the connection gap of the water stop.
The rubber water stop is a sealed rubber element, and its shear failure strain is up to more than 300%. Categorizing by function, the water stop’s failure modes include relaxation creep failure, stiffness failure, and stability failure. Categorizing by material properties, the water stop’s failure modes include fatigue failure, limit failure, and bond failure. Comprehensively considering the bonding, vulcanization process, functionality, and ultimate damage of the water stop, the damage index is generalized in Table 5.

5.2. Definition of the Pier Limit State

At present, the damage indexes of bridges are mainly characterized by the curvature or drift rate of the pier. From the angle of bridge engineering, the drift rate of the pier is applicable and accurate during the fragility analysis. Thus, the drift rate of the pier is used as the pier damage index and defined as Equation (6):
d 0 = d c / H
where H is the pier height, and dc is the maximum displacement of the pier top.
The crack of the protective layer of the pier, the yield of the tensile reinforcement, the limit of the tensile reinforcement, and the collapse of the pier column were selected as the boundaries of the damage indices. The damage indices are listed in Table 6.

5.3. Seismic Fragility Evaluation

Seismic fragility represents the conditional possibility of the exceedance of a limit state under a specific intensity measure. It was expressed as the probability of exceeding a specified performance level for the structure in question by convolving the probability distributions for demand, capacity, and ground motion intensity hazard [32]. The lack of experience in aqueduct damage data, the subjective definition of the limit state, and the advancement of computational tools prove the superiority of this analytical method. It permits predicting the maximum deformation behavior and the collapse risk of an aqueduct under dynamic loadings properly.
To be specific, the exceeding probability Pf = P(ai, Ek) is defined as the probability that the randomly selected aqueduct samples are in the limit state of Ek under PGA = ai. The conditional probability of a structure exceeding a specific damage level for a given earthquake intensity can be generally expressed as Equation (7).
P ƒ = P m d m c 1 IM = Φ ln m d / m c β d 2 + β c 2
where md and mc are the demand and capacity of the aqueduct. IM is the PGA (g); βd is the logarithmic standard deviation of structural demand. βc is the logarithmic standard deviation of structural capacity. βc us 0.6 in this paper. Φ(x) is the standard normal cumulative distribution function. After solving and obtaining x, the corresponding transcendence probability Pf can be obtained by checking the table according to the size of x.
The median engineering demand parameters md and the ground IMs are assumed to follow a logarithmic correlation, as expressed in Equation (8) [33].
ln ( m d ) = ln α + b ln IM
where a and b are constant coefficients that can be predicted from the regression analysis.
Compared with the bridge, the water load of the aqueduct lifts the center of gravity of the structure, which theoretically makes the aqueduct structure more susceptible to the transverse force, especially the hydrodynamic influence. Therefore, it was important to ensure whether the horizontal or vertical action was more vulnerable in preanalysis. Itwas found that the model is more vulnerable in the transverse direction after comparing the bending moment of the pier column, which also agreed with the dynamic characteristic of the bridge. A large number of seismic time-history analyses were carried out to investigate the transverse seismic vulnerability of the aqueduct. Using regression parameters, the fitting results are expressed by Equations (9) and (10):
ln tan γ = 0.82053 × ln PGA + 0.205969
ln d 0 = 0.974 × ln PGA 0.5781
Each of the depicted points in Figure 8a shows the logarithmic regression relationship between tanγ and PGA. Each of the depicted points in Figure 8b shows the logarithmic regression relationship between d0 and PGA.
After that, Equations (11) and (12) can be derived from combining Equations (8)–(10).
P ƒ = P m d m c 1 P G A = Φ ln tan Ƴ / tan Ƴ β C = Φ ln 7.844 × P G A t a n γ 0.82053 0.72744
P ƒ = P m d m c 1 P G A = Φ ln tan Ƴ / tan Ƴ β C = Φ ln 7.844 × P G A t a n γ 0.82053 0.72744
In the end, according to the standard normal distribution table, the fragility curves of the water stop and the pier were drawn, as shown in Figure 9a,b.

6. Discussion

It can be seen from Figure 9a that the performance index of the water stop is put forward based on the definition of the shear failure of the rubber bearing. The analysis results show that with the increase in ground motion, the probability of damage to the water stop significantly increases. In the range of PGA < 0.5 g, the water stop is in an intact state, and slight damage is a high probability, while in the range of 0.15–0.3 g, the water stop is prone to be in a slight or moderate damage state. The water stop is gradually transformed to complete damage status when the PGA exceeds 0.3 g and is more than 50% to be completely destroyed. It can be found from Figure 9a that the water stop has good ductility.
Through the analysis of 90 samples, the relationship between the upper and lower limits of the dimensionless drift rate of the pier under the damage limit state and PGA was obtained. The results are shown in Figure 10.
Comparisons of the fragility curves in terms of the pier are presented in Figure 10. In the range of PGA < 0.15 g, the aqueduct is most likely to be in an intact or slight damage state. In the PGA range of 0.15–0.2g, the aqueduct is most likely to be in a state of slight damage or complete damage, while when the PGA range exceeds 0.2 g, the probability of complete damage is more than 50% and increasing rapidly.
As observed from Figure 11, the aqueduct at present is vulnerable. Seismic rehabilitation is necessary under the current circumstances, and the seismic vulnerability of the water stop is not as significant as the pier. However, the seismic vulnerability of the water stop may be influenced by the performance of the pier and the height difference between adjacent piers.

7. Conclusions

This paper presents a seismic fragility evaluation of an old aqueduct considering leakage risk. Based on the Hua Shigou aqueduct in Ningxia Province, the numerical analysis for risk evaluation was evolved from a probabilistic model, and specific water stop damage states were proposed. By a cluster of IDA, the seismic fragility of the old aqueduct was obtained. The conclusions below are drawn:
  • An analytical model incorporating water variant, concrete strength, and bearing performance degradation is proposed, in which a water stop is considered a nonlinear gap element. This analytical model can be largely used for the vulnerability analysis of general aqueducts.
  • The damage state of the rubber water stop is proposed based on the rubber pad theory. The damage limit state of the water stop is defined by the tangent of the shear angle. In the range of tanγ < 0.5, the damage status is defined as intact; in the range of 0.5–1, the damage status is defined as slight damage; in the range of 1–2, the damage status is defined as moderate damage; in the range of 2–3, the damage status is defined as extensive damage; and when tanγ exceeds 3, the damage status is defined as complete damage.
  • With the increase in ground motion, the probability of the water stop’s leakage risk is significantly increased. However, the water stop has better ductility compared with the pier according to the vulnerability curve shown in Figure 11. It is necessary to evaluate the seismic fragility of old aqueducts with long service history, and seismic rehabilitation could be needed for these aqueducts based on the vulnerability curve.
The effect of the height difference between adjacent piers and the local damage of prestressed concrete superstructures under fluid–structure coupling will be studied in future work.

Author Contributions

Conceptualization, Z.X. and A.Z.; methodology, C.L.; investigation, Z.X. and C.L.; resources, A.Z.; writing—original draft preparation, Z.X. and C.L.; writing—review and editing, Z.X.; visualization, H.Z. and J.L.; supervision, Z.X.; funding acquisition, Z.X. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the National Natural Science Foundation of China (Grant No. 51908422) and the Elite Scholar Program of Northwest A&F University (Grant No. Z111022001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, C.; Zheng, H. South-to-north Water Transfer Schemes for China. Int. J. Water Resour. Dev. 2002, 18, 453–471. [Google Scholar] [CrossRef]
  2. Wu, Y.; Li, L.; Liu, Z.; Chen, X.; Huang, H. Real-Time Control of the Middle Route of South-to-North Water Diversion Project. Water 2021, 13, 97. [Google Scholar] [CrossRef]
  3. Xiong, Z.H.; Zhang, A.J.; Liu, Y.Q. Structural type and load effects research of large span aqueduct–bridge co-built structure. J. China Inst. Water Resour. Hydropower Res 2021. in press (In Chinese) [Google Scholar] [CrossRef]
  4. Valeti, B.; Ray-Chaudhuri, S.; Raychowdhury, P. Seismic response of an elevated aqueduct considering hydrodynamic and soil-structure interactions. Int. J. Adv. Struct. Eng. 2016, 8, 53–71. [Google Scholar] [CrossRef] [Green Version]
  5. Housner, G.W. Dynamic pressures on accelerated fluid containers. Bull. Seismol. Soc. Am. 1957, 47, 15–35. [Google Scholar]
  6. Housner, G.W. The dynamic behavior of water tanks. Bull. Seismol. Soc. Am. 1963, 53, 381–389. [Google Scholar]
  7. ACI Committee. Seismic Design of Liquid-Containing Concrete Structures (ACI 350.3-01) and Commentary (ACI 350.3 R-01); An ACI Standard; American Concrete Institute: Farmington Hills, MI, USA, 2001. [Google Scholar]
  8. Virella, J.C.; Godoy, L.A.; Suárez, L.E. Dynamic buckling of anchored steel tanks subjected to horizontal earthquake excitation. J. Constr. Steel Res. 2006, 62, 521–531. [Google Scholar] [CrossRef]
  9. La Salandra, V.; Wenzel, M.; Bursi, O.S.; Carta, G.; Movchan, A.B. Conception of a 3D metamaterial-based foundation for static and seismic protection of fuel storage tanks. Front. Mater. 2017, 4, 30. [Google Scholar] [CrossRef] [Green Version]
  10. Terenzi, G.; Rossi, E. Seismic analysis and retrofit of the oldest R/C elevated water tank in Florence. Bull. Earthq. Eng. 2018, 16, 3081–3102. [Google Scholar] [CrossRef]
  11. Meng, X.; Li, X.; Xu, X.; Zhang, J.; Zhou, W.; Zhou, D. Earthquake Response of Cylindrical Storage Tanks on an Elastic Soil. J. Vib. Eng. Technol. 2019, 7, 433–444. [Google Scholar] [CrossRef]
  12. D’alessandro, V. Modeling of tank vehicle dynamics by fluid sloshing coupled simulation. Ph.D. Dissertation, Politecnico di Milano, Milan, Lombardia, Italy, March 2012. [Google Scholar]
  13. Rawat, A.; Mittal, V.; Chakraborty, T.; Matsagar, V. Earthquake induced sloshing and hydrodynamic pressures in rigid liquid storage tanks analyzed by coupled acoustic-structural and Euler-Lagrange methods. Thin Walled Struct. 2019, 134, 333–346. [Google Scholar] [CrossRef]
  14. Zhao, C.; Chen, J.; Wang, J.; Yu, N.; Xu, Q. Seismic mitigation performance and optimization design of NPP water tank with internal ring baffles under earthquake loads. Nucl. Eng. Des. 2017, 318, 182–201. [Google Scholar] [CrossRef]
  15. Ozdemir, Z.; Souli, M.; Fahjan, Y.M. Application of nonlinear fluid–structure interaction methods to seismic analysis of anchored and unanchored tanks. Eng. Struct. 2010, 32, 409–423. [Google Scholar] [CrossRef]
  16. Gibson, R.; Mistry, A.; Go, J.; Lubkowski, Z. Fluid Structure Interaction vs. Lumped Mass Analogue for Storage Tank Seismic Assessment. In Proceedings of the Earthquake Risk and Engineering towards a Resilient World, Cambridg, UK, 9–10 July 2015. [Google Scholar]
  17. Xiong, Z.H.; Feng, Y.C.; Song, S.L.; Wang, J.B. Optimization Design of Large Span Cable-Stayed Bridge in High Seismic Risk Zone. Appl. Mech. Mater. 2013, 353–356, 2015–2019. [Google Scholar] [CrossRef]
  18. Phan, H.N.; Paolacci, F.; Corritore, D.; Alessandri, S. Seismic vulnerability analysis of storage tanks for oil and gas industry. Sci. Technol. 2018, 8, 55–65. [Google Scholar] [CrossRef] [Green Version]
  19. Bakalis, K.; Vamvatsikos, D.; Fragiadakis, M. Seismic risk assessment of liquid storage tanks via a nonlinear surrogate model. Earthq. Eng. Struct. Dyn. 2017, 46, 2851–2868. [Google Scholar] [CrossRef]
  20. Kanyilmaz, A.; Castiglioni, C.A. Reducing the seismic vulnerability of existing elevated silos by means of base isolation devices. Eng. Struct. 2017, 143, 477–497. [Google Scholar] [CrossRef]
  21. Tsipianitis, A.; Tsompanakis, Y. Impact of damping modeling on the seismic response of base-isolated liquid storage tanks. Soil Dyn. Earthq. Eng. 2019, 121, 281–292. [Google Scholar] [CrossRef]
  22. Faggella, M.; LaGuardia, R.; Gigliotti, R.; Morelli, F.; Braga, F.; Salvatore, W.; Papadrakakis, M. Performance-Based Nonlinear response History Analysis Framework for the “Proindustry” Project Case Studies. In Proceedings of the VII European Congress on Computational Methods in Applied Sciences and Engineering (ECCOMAS Congress 2016), Cret Island, Greece, 5–10 June 2016; Volume 3, pp. 5912–5925. [Google Scholar] [CrossRef]
  23. Prashar, Y.; Irias, X.J.; Shewbridge, S.E. East Bay Municipal Utility Districts Aqueduct System in the Sacramento-San Joaquin Delta: Hazard Evaluation. In Proceedings of the Pipelines Specialty Conference 2009, San Diego, CA, USA, 15–19 August 2009; pp. 1263–1275. [Google Scholar] [CrossRef]
  24. Rafiee, A.; Vinches, M.; Bohatier, C. Modelling and analysis of the Nîmes arena and the Arles aqueduct subjected to a seismic loading, using the Non-Smooth Contact Dynamics method. Eng. Struct. 2008, 30, 3457–3467. [Google Scholar] [CrossRef]
  25. Liu, Y.; Dang, K.; Dong, J. Finite element analysis of the aseismicity of a large aqueduct. Soil Dyn. Earthq. Eng. 2017, 94, 102–108. [Google Scholar] [CrossRef]
  26. Ministry of Water Resources of the People’s Republic of China. Standard for Seismic Design of Hydraulic Structures (GB 51247-2018); Ministry of Housing and Urban-Rural Development of the People’s Republic of China (MOHURD): Beijing, China, 2018.
  27. The PEER Ground Motion Database. Available online: https://ngawest2.berkeley.edu (accessed on 13 May 2021).
  28. Pan, Y.; Agrawal, A.K.; Ghosn, M.; Alampalli, S. Seismic Fragility of Multispan Simply Supported Steel Highway Bridges in New York State. I: Bridge Modeling, Parametric Analysis, and Retrofit Design. J. Bridg. Eng. 2010, 15, 448–461. [Google Scholar] [CrossRef]
  29. Ge, S.J.; Xiong, Z.H.; Zhai, M.G.; Pan, C.P. Research on Seismic Fragility of Medium and Small Sized Concrete Continuous Girder Bridge. J. Highw. Transp. Res. Dev. 2013, 30, 60–65. [Google Scholar] [CrossRef]
  30. Song, T.H.; Meng, C.Q.; Wang, Y.Q.; Li, Y.M.; Zhong, D.Y. Further Discussion on the Necessity of the Construction of Daliushu Water Control Project. Yellow River 2020, 42, 42–47. [Google Scholar] [CrossRef]
  31. Zhang, J. Bridge Bearing, 3rd ed.; China Railway Publishing House: Beijing, China, 2010. [Google Scholar]
  32. Cornell, C.A.; Jalayer, F.; Hamburger, R.O.; Foutch, D.A. Probabilistic Basis for 2000 SAC Federal Emergency Management Agency Steel Moment Frame Guidelines. J. Struct. Eng. 2002, 128, 526–533. [Google Scholar] [CrossRef] [Green Version]
  33. Gardoni, P.; Mosalam, K.M.; Der Kiureghian, A. Probabilistic seismic demand models and fragility estimates for RC bridges. J. Earthq. Eng. 2003, 7, 79–106. [Google Scholar] [CrossRef]
Figure 1. Elevation of the aqueduct and the layout of the water stop (unit: mm).
Figure 1. Elevation of the aqueduct and the layout of the water stop (unit: mm).
Water 13 01404 g001
Figure 2. Schematic diagram of the aqueduct in the transverse direction (unit: mm).
Figure 2. Schematic diagram of the aqueduct in the transverse direction (unit: mm).
Water 13 01404 g002
Figure 3. (a) Sloshing diagram of a tank; (b) equivalent dynamic system for a water tank.
Figure 3. (a) Sloshing diagram of a tank; (b) equivalent dynamic system for a water tank.
Water 13 01404 g003
Figure 4. (a) FE models and connecting links; (b) 3D model view.
Figure 4. (a) FE models and connecting links; (b) 3D model view.
Water 13 01404 g004
Figure 5. The design response spectrum.
Figure 5. The design response spectrum.
Water 13 01404 g005
Figure 6. The response spectrum from the PEER Ground Motion Database.
Figure 6. The response spectrum from the PEER Ground Motion Database.
Water 13 01404 g006
Figure 7. Shear deformation of the water stop.
Figure 7. Shear deformation of the water stop.
Water 13 01404 g007
Figure 8. Logarithmic regression relationship: (a) tanγ and PGA(g); (b) d0 and PGA (g).
Figure 8. Logarithmic regression relationship: (a) tanγ and PGA(g); (b) d0 and PGA (g).
Water 13 01404 g008
Figure 9. (a) Fragility curves of the water stop; (b) fragility curves of the pier.
Figure 9. (a) Fragility curves of the water stop; (b) fragility curves of the pier.
Water 13 01404 g009
Figure 10. Fragility curves of aqueducts: (a) slight damage; (b) moderate damage; (c) extensive damage; (d) complete damage.
Figure 10. Fragility curves of aqueducts: (a) slight damage; (b) moderate damage; (c) extensive damage; (d) complete damage.
Water 13 01404 g010aWater 13 01404 g010b
Figure 11. Comparison of fragility curves of the pier and the water stop.
Figure 11. Comparison of fragility curves of the pier and the water stop.
Water 13 01404 g011
Table 1. Ground motion records characteristics.
Table 1. Ground motion records characteristics.
No.EarthquakeStationYearMagnitudeHorizontal Acc.
1Imperial Valley-02El Centro Array #919406.95RSN6_IMPVALL.I_I-ELC180
2Kern CountySanta Barbara Courthouse19527.36RSN14_KERN_SBA132
3Kern CountyTaft Lincoln School19527.36RSN15_KERN_TAF021
4Northern Calif-03Ferndale City Hall19546.5RSN20_NCALIF.FH_H-FRN044
5San FernandoLA–Hollywood Stor FF19716.61RSN68_SFERN_PEL090
6San FernandoLake Hughes #119716.61RSN70_SFERN_L01021
7San FernandoPalmdale Fire Station19716.61RSN78_SFERN_PDL120
8San FernandoPasadena–CIT Athenaeum19716.61RSN79_SFERN_PAS090
9Managua_Nicaragua-02Managua_ ESSO19725.2RSN96_MANAGUA_B-ESO180
Table 2. Variable distribution of the sampling results.
Table 2. Variable distribution of the sampling results.
VariableProbability DistributionμσRange of Variations
Water weightNormal distribution56.6 kN/m14.1 kN/m0.2–113 kN
ƒ′cNormal distribution20.1 MPa2.5 MPa15–25 MPa
Bearing friction coefficientUniform distribution 0.02–0.05
Table 3. Latin hypercube sampling results.
Table 3. Latin hypercube sampling results.
Group No.Water Weight (kN)ƒ′c (MPa)Friction Coefficient of Bearing
158.9720.70.0468
222.122.30.0492
383.916.50.0395
448.923.90.0300
539.918.40.0356
686.220.30.0221
Table 4. Ground motion parameter.
Table 4. Ground motion parameter.
Ground Motion No.123456789
Response acceleration, gmax0.800.350.600.480.410.410.380.410.80
Maximum acceleration, g0.280.130.160.160.220.150.110.110.22
AM acceleration IDA, g 0.390.320.480.660.750.77
Table 5. Definition of different damage states for the water stop.
Table 5. Definition of different damage states for the water stop.
Damage Limit StateDamage Index
Intacttanγ < 0.5
Slight damage0.5 ≤ tanγ < 1
Moderate damage1 ≤ tanγ < 2
Extensive damage2 ≤ tanγ < 3
Complete damagetanγ ≥ 3
Table 6. Definition of different damage states of the pier.
Table 6. Definition of different damage states of the pier.
Damage Limit StateDamage Index
Intactd0 < 0.0246%
Slight damage0.0246% ≤ d0 < 0.084%
Moderate damage0.084% ≤ d0 < 0.09%
Extensive damage0.09% ≤ d0 < 0.113%
Complete damaged0 ≥ 0.113%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xiong, Z.; Liu, C.; Zhang, A.; Zhu, H.; Li, J. Seismic Fragility Evaluation of Simply Supported Aqueduct Accounting for Water Stop’s Leakage Risk. Water 2021, 13, 1404. https://doi.org/10.3390/w13101404

AMA Style

Xiong Z, Liu C, Zhang A, Zhu H, Li J. Seismic Fragility Evaluation of Simply Supported Aqueduct Accounting for Water Stop’s Leakage Risk. Water. 2021; 13(10):1404. https://doi.org/10.3390/w13101404

Chicago/Turabian Style

Xiong, Zhihua, Chen Liu, Aijun Zhang, Houda Zhu, and Jiawen Li. 2021. "Seismic Fragility Evaluation of Simply Supported Aqueduct Accounting for Water Stop’s Leakage Risk" Water 13, no. 10: 1404. https://doi.org/10.3390/w13101404

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop