Dynamic Damage Mechanism and Seismic Fragility Analysis of an Aqueduct Structure

The Concrete Damaged Plasticity (CDP) constitutive is introduced to study the dynamic failure mechanism and the law of damage development to the aqueduct structure during the seismic duration using a large-scale aqueduct structure from the South-to-North Water Division Project (SNWDP) as a research object. Incremental dynamic analysis (IDA) and multiple stripe analysis (MSA) seismic fragility methods are introduced. The spectral acceleration is used as the scale of ground motion record intensity measure (IM), and the aqueduct pier top offset ratio quantifies the limit of structural damage measure (DM). The aqueduct structure’s seismic fragility evaluation curves are constructed with indicators of different seismic intensity measures to depict the damage characteristics of aqueduct structures under different seismic intensities through probability. The results show that penetrating damage is most likely to occur on both sides of the pier cap and around the pier shaft in the event of a rare earthquake, followed by the top of the aqueduct body, which requires the greatest care during an earthquake. The results of two fragility analysis methodologies reveal that the fragility curves are very similar. The aqueduct structure’s first limit state level (LS1) is quite steep and near the vertical line, indicating that maintaining the excellent condition without damage in the seismic analysis will be challenging. Except for individual results, the overall fragility results are in good agreement, and the curve change rule is the same. The exceedance probability in the case of any ground motion record IM may be estimated using only two factors when using the MSA approach, and the computation efficiency is higher. The study of seismic fragility analysis methods in this paper can provide a reference for the seismic safety evaluation of aqueducts and similar structures.


Introduction
The aqueduct is a water conveyance structure commonly used for irrigation, water transfer, and supply in water conservancy projects. The construction of the aqueduct structure can not only alleviate the severe water shortage disaster for human beings, but also contributes to the local economic, science and technology, and social development, achieving the main sustainable development goals (SDG) that are employed by the United Nations (UN) Agenda, such as the Valens Aqueduct (Turkey) [1]. The structural type of aqueduct is mainly divided into rectangular, semicircular, U-shaped, and so on [2,3]. Its dynamic response law under seismic loads has always been at the center of. Zhang et al. [4] considered the concrete elastoplastic random damage constitutive relationship and the random seismic excitation model to carry out a random nonlinear seismic response analysis of the double trough aqueduct structure. The results showed that the response value of the mid-span bottom plate and the pier bottom of the aqueduct is relatively large, which belongs to the weak position of earthquake resistance. Ying et al. [5] studied the foundation soil property on the dynamic characteristics and responses of rectangular aqueducts. The results show that the increase of the soil stiffness causes amplification in impulsive displacement. Ding et al. [6] studied the influence of aqueduct height on the structure and proved that the higher aqueduct has a greater structural dynamic response. Some researchers have looked into the dynamic failure process of aqueduct structures under earthquake ground motion, taking into account the effect of fluid-structure interactions [7][8][9][10]. Li et al. [11] proposed a simplified beam-water coupled system to investigate the seismic ground motion response to long spans of large-scale aqueduct structures. The traditional seismic dynamic analysis focuses on the dynamic response of the aqueduct structure itself under the action of ground motions, which can directly reflect the actual conditions of the aqueduct structure under seismic excitation. Most of the large-scale aqueduct structures studied by researchers today are in areas prone to earthquakes. Research on the damage and collapse of aqueduct structures under powerful earthquakes is critical to the structural safety assessment of engineering procedures to assure the safety of aqueducts and normal regulation of water resources during strong earthquakes.
Performance-based seismic engineering theory collectively refers to the damage analysis and seismic response analysis in the structural earthquake-resistant evaluation as a seismic fragility analysis [12][13][14]. Bertero et al. [15] first proposed the idea of an incremental time-history analysis, which Vamvatsikos and Cornell [16,17] perfected as the earliest method of incremental dynamic analysis. The aqueduct structure is similar to bridge structures in terms of engineering design. The pier shaft is responsible for sustaining the superstructure and affecting the overall stability of both structures. Pietro et al. [18,19] proposed an efficient procedure for the collapse mechanism evaluation of existing reinforced concrete motorway bridges under horizontal loads, providing useful information for the maintenance of existing bridges in terms of repair and strengthening interventions. In recent years, various researchers have applied IDA to the construction of frames, tunnels, bridges, and other fields [20][21][22][23]. However, the probability of structures exceeding the damage index under the action of earthquakes may not strictly increase with the ground motion record IM in a large number of numerical simulations. This statistical probability method of incremental dynamic analysis may not accurately give the corresponding value when the structure reaches a limit state. Therefore, the seismic fragility assessment may be limited. In 2003, Jalayer [24] proposed the MSA, a fragility analysis method using discrete IM. The fragility analysis with discrete intensity parameters is performed through the MSA method based on the mean value and logarithmic standard deviation of the fragility function, estimated through the maximum likelihood estimation (MLE). There is no need to adjust the ground motion intensity parameters to the state where the structural response goes beyond the limit. Jalayer and Cornell [25] compared the IDA and MSA methods, and the results showed that the MSA method could provide the precise mean and standard deviation required for a fragility assessment. On performing a structural dynamic analysis under a set number of seismic waves, Baker and Jack [26] discovered that the MSA technique could better predict the fragility equation parameters than the IDA method, with a shorter analysis time and higher computation efficiency. Yuan et al. [27] and Mahmoodi et al. [28] used the MSA method to analyze the fragility of wind turbines and dams, providing a new path for the safety assessment of engineering structures. Giammaria et al. [29] adapted the MSA method to analyze the seismic performance assessment of a liquid-storage tank installed in an industrial steel moment-resisting frame building. The dynamic response of the tank was systematically studied, and has a certain reference value for this kind of structure. So far, the MSA method with good applicability has not been applied to the structural engineering of aqueducts.
This paper aimed at studying the dynamic failure mechanism and the rules of damage for concrete aqueduct structures subject to seismic ground motion. For this purpose, the CDP constitutive is introduced to investigate the dynamic damage mechanism and damage development law of a concrete aqueduct structure under the influence of a rare earthquake Appl. Sci. 2021, 11, 11709 3 of 22 and to identify weak spots where concrete fracturing could occur. Depending on the spectral similarity of the response spectrum of the aqueduct, the natural seismic waves that satisfy the calculation conditions and artificial waves that are in good agreement with the response spectrum are selected for the hunt and fill amplitude modulation and nonlinear dynamic analysis. To construct a fragility curve that reflects the deformation capacity and describes the specific damage degree of the aqueduct structure under the exceedance probability condition, the spectral acceleration is used as the ground motion record intensity measure, and structural response indicators suitable for evaluating weak parts of the aqueduct structure against earthquakes are used. The calculation results based on different fragility methods are discussed, and various limit states are evaluated.

Concrete Damage Plasticity (CDP) Constitutive Model
The CDP constitutive model assumes that material is damaged by isotropic compression and tension. The damage factor d (Damage) is presented based on the plastic model to represent the irrecoverability of materials produced by damage to simulate the actual mechanical behavior of concrete [30,31]. It ranges from 0 (the material is not harmed) to 1 (the material is entirely destroyed). Based on the energy equivalent hypothesis, this study adapts Sidoroff's [32] universal compression and tensile damage factor calculation approach and integrates it directly into the ABAQUS program. In Figure 1, taking C40 concrete as an example, the tensile-compression damage plastic parameters are calculated. ural seismic waves that satisfy the calculation conditions and artificial wa good agreement with the response spectrum are selected for the hunt and modulation and nonlinear dynamic analysis. To construct a fragility cur the deformation capacity and describes the specific damage degree of structure under the exceedance probability condition, the spectral accelera the ground motion record intensity measure, and structural response ind for evaluating weak parts of the aqueduct structure against earthquakes calculation results based on different fragility methods are discussed, an states are evaluated.

Concrete Damage Plasticity (CDP) Constitutive Model
The CDP constitutive model assumes that material is damaged by pression and tension. The damage factor d (Damage) is presented based model to represent the irrecoverability of materials produced by damage actual mechanical behavior of concrete [30,31]. It ranges from 0 (the m harmed) to 1 (the material is entirely destroyed). Based on the energy pothesis, this study adapts Sidoroff's [32] universal compression and t factor calculation approach and integrates it directly into the ABAQUS pr ure 1, taking C40 concrete as an example, the tensile-compression dam rameters are calculated. According to the plastic damage model theory in ABAQUS, the tota ε is composed of elastic strain rate el ε , and equivalent plastic strain ra determined as follows el pl ε ε ε = + When there is no damage to the concrete, the stress-strain relations crete is According to the plastic damage model theory in ABAQUS, the total strain tensor ε is composed of elastic strain rate ε el , and equivalent plastic strain rate ε pl can be determined as follows When there is no damage to the concrete, the stress-strain relationship of the concrete is in which σ is the total stress and D el is the elastic stiffness matrix. When the concrete is damaged, the material damage factor is introduced to describe the degradation of stiffness, In this formula, σ is the effective stress, and d is the damage factor, which varies from 0 (no damage) to 1 (complete damage).
The complex damage mechanism of concrete under cyclic alternating stress is related to the cracking, consolidation, and changing of early cracks. When concrete is changed from tension to compression, the stiffness of the concrete recovers locally. To reflect the phenomenon, it is assumed that the relationship between the damage factor and tensile and compressive damage variables (d t and d c ) under an alternating load is where s t and s c are the relationship of stress action under stiffness recovery, which can be determined as follows where r * (σ) is the weight factor in the relationship of principal stress under a multiaxial state, which is in which w t and w c are stiffness recovery weighting factors related to material properties. σ i (i = 1, 2, 3) is the main stress component; • can be expressed as x = (|x| + x)/2 Figure 1 shows the stiffness recovery curve of the concrete damage model under a uniaxial alternating load. The weight factors under tension and compression are ω t = 0 (compression → tension) and ω c = 1 (tension → compression), respectively.
The different evolution laws of material strength under tension and compression load are considered in the yield function of the CDP model, and the expression of effective stress is where α and γ are dimensionless constants;σ max is the maximum effective principal stress; the effective hydrostatic pressure and Mises equivalent stress are p and q; σ b0 is the equal axial initial yield stress, and σ c0 is the non-equiaxial initial yield stress; K c is the parameter of the plane yield curve of deviatoric stress; I is the unit matrix quantity; and S is the effective stress deviation. The CDP model adopts the flow potential G expression of the Drucker-Prager hyperbola as λ is the plasticity factor, and G is the Drucker-Prager hyperbolic formula, which can be expressed as G = (∈ σ t0 tan ψ) 2 + q 2 − p tan ψ (10) in which ψ is the expansion angle in the meridian plane under a high confining pressure; σ t0 is the ultimate uniaxial tensile strength; and ∈ is the variable of hyperbolic eccentricity, which is used to describe the rate at which the potential function approaches its asymptote.

Case Model
This paper takes an aqueduct of the SNWDP in China as an example. The body of the aqueduct is rectangular with a length of 60 m, with a bottom width of 15 m and a height of 8.95 m. The pier cap is 17.10 m long, 5.20 m wide, and 2.50 m high. The pier shaft is 5.50 m in height; and the pile cap is 16.10 m in length, 10.20 m in width, and 3.50 m in height. There were 71,172 three-dimensional solid elements used to build the body of the aqueduct and the foundation model. The model mainly uses the eight-node hexagonal Solid65 element, and the pull beam is simulated by the Beam 3 element (Figure 2). The main materials are shown in Table 1. In the dynamic calculation, an artificial viscous-spring boundary is applied to the foundation to simulate the seismic wave propagation process and reflect the energy dispersion at the boundary. The far boundary of the foundation is at a distance from the aqueduct of approximately 1.5~2 times the aqueduct in all directions. The dynamic response is analyzed by the commercial finite element software ABAQUS. In this study, a dynamic implicit analysis is performed.
For the seismic analysis of the aqueduct structure, natural frequencies and modes are usually used to study its vibration characteristics. In the solutions to most seismic problems, the first sloshing mode of a structure makes a dominant contribution to all the other modes [2]. The natural frequency of the aqueduct selected in this paper is 1.162 Hz when the design water level is 7.50 m. Table 2 shows the first eight vibration modes of the aqueduct structure. Compared to the natural frequencies of other aqueducts of similar height, width, and length under wet mode conditions, this aqueduct's natural frequency and mode of vibration are similar to other aqueducts [33,34]. For example, the first natural frequency of an aqueduct (50 m × 14 m × 7 m) is 1.199 Hz [35]. We refer to the results of a doctoral dissertation. The aqueduct structure is similar to it in this paper (single 30 m × 20 m × 8 m) [35]. In the dynamic model experiment, the first frequency of modal identification is about 0.9~1.2 Hz, and the finite element calculation result is 0.973 Hz. Therefore, the natural frequency results of the aqueduct structure in this paper can be used to validate the accuracy of the finite element model. In this study, Rayleigh damping is used for the analysis, and the damping ratio corresponding to the structural frequency is taken as 5%, α = 0.366, and β = 0.007.

Ground Motion Input
The engineering foundation site classification is Class II, and the seismic intensity scale is VIII. To better understand the structural vibration response, damage development process, and damage degree caused by strong earthquakes, the data from the strong earthquake database of the Pacific Earthquake Engineering Research (PEER) Centre was used as a reference to select waves with a magnitude of 6~8, with shear wave velocity of 500~750 m/s, and epicentral distance of 0~105 km. The natural wave Northridge-01 (Whittier-S. Alta Dr), which was close to the characteristic site period of 0.4 s, has been chosen, following the frequency band fixation method proposed by the Applied Technology Council (ATC-3) in 1978 [36]. The basic parameters of selected natural waves are shown in Table 3. The horizontal peak acceleration of severe earthquakes is 0.19 g, and the ratio to the vertical peak acceleration is 1:0.65. The seismic wave input interval is 0.01 s, and the total duration is 25 s. Two acceleration time−history curves are shown in Figure 3. The additional mass approach proposed by Westergaard in 1933 during a seismic study on the dam body is also used when considering the force of the seismic load on the water body [37].

Material
Elastic Modulus (GPa)    quency and mode of vibration are similar to other aqueducts [33,34]. For example, the first natural frequency of an aqueduct (50 m × 14 m × 7 m) is 1.199 Hz [35]. We refer to the results of a doctoral dissertation. The aqueduct structure is similar to it in this paper (single 30 m × 20 m × 8 m) [35]. In the dynamic model experiment, the first frequency of modal identification is about 0.9~1.2 Hz, and the finite element calculation result is 0.973 Hz. Therefore, the natural frequency results of the aqueduct structure in this paper can be used to validate the accuracy of the finite element model. In this study, Rayleigh damping is used for the analysis, and the damping ratio corresponding to the structural frequency is taken as 5%, =0.366 α , and =0.007 β .

Ground Motion Input
The engineering foundation site classification is Class II, and the seismic intensity scale is Ⅷ. To better understand the structural vibration response, damage development process, and damage degree caused by strong earthquakes, the data from the strong earthquake database of the Pacific Earthquake Engineering Research (PEER) Centre was used as a reference to select waves with a magnitude of 6~8, with shear wave velocity of 500~750 m/s, and epicentral distance of 0~105 km. The natural wave Northridge-01 (Whittier-S. Alta Dr), which was close to the characteristic site period of 0.4 s, has been chosen, following the frequency band fixation method proposed by the quency and mode of vibration are similar to other aqueducts [33,34]. For example, the first natural frequency of an aqueduct (50 m × 14 m × 7 m) is 1.199 Hz [35]. We refer to the results of a doctoral dissertation. The aqueduct structure is similar to it in this paper (single 30 m × 20 m × 8 m) [35]. In the dynamic model experiment, the first frequency of modal identification is about 0.9~1.2 Hz, and the finite element calculation result is 0.973 Hz. Therefore, the natural frequency results of the aqueduct structure in this paper can be used to validate the accuracy of the finite element model. In this study, Rayleigh damping is used for the analysis, and the damping ratio corresponding to the structural frequency is taken as 5%, =0.366 α , and =0.007 β .

Ground Motion Input
The engineering foundation site classification is Class II, and the seismic intensity scale is Ⅷ. To better understand the structural vibration response, damage development process, and damage degree caused by strong earthquakes, the data from the strong earthquake database of the Pacific Earthquake Engineering Research (PEER) Centre was used as a reference to select waves with a magnitude of 6~8, with shear wave velocity of 500~750 m/s, and epicentral distance of 0~105 km. The natural wave Northridge-01 (Whittier-S. Alta Dr), which was close to the characteristic site period of 0.4 s, has been chosen, following the frequency band fixation method proposed by the 1.271 Hz Vertical 2 quency and mode of vibration are similar to other aqueducts [33,34]. For example, the first natural frequency of an aqueduct (50 m × 14 m × 7 m) is 1.199 Hz [35]. We refer to the results of a doctoral dissertation. The aqueduct structure is similar to it in this paper (single 30 m × 20 m × 8 m) [35]. In the dynamic model experiment, the first frequency of modal identification is about 0.9~1.2 Hz, and the finite element calculation result is 0.973 Hz. Therefore, the natural frequency results of the aqueduct structure in this paper can be used to validate the accuracy of the finite element model. In this study, Rayleigh damping is used for the analysis, and the damping ratio corresponding to the structural frequency is taken as 5%, =0.366 α , and =0.007 β .

Ground Motion Input
The engineering foundation site classification is Class II, and the seismic intensity scale is Ⅷ. To better understand the structural vibration response, damage development process, and damage degree caused by strong earthquakes, the data from the strong earthquake database of the Pacific Earthquake Engineering Research (PEER) Centre was used as a reference to select waves with a magnitude of 6~8, with shear wave velocity of 500~750 m/s, and epicentral distance of 0~105 km. The natural wave Northridge-01 (Whittier-S. Alta Dr), which was close to the characteristic site period of 0.4 s, has been chosen, following the frequency band fixation method proposed by the quency and mode of vibration are similar to other aqueducts [33,34]. For example, the first natural frequency of an aqueduct (50 m × 14 m × 7 m) is 1.199 Hz [35]. We refer to the results of a doctoral dissertation. The aqueduct structure is similar to it in this paper (single 30 m × 20 m × 8 m) [35]. In the dynamic model experiment, the first frequency of modal identification is about 0.9~1.2 Hz, and the finite element calculation result is 0.973 Hz. Therefore, the natural frequency results of the aqueduct structure in this paper can be used to validate the accuracy of the finite element model. In this study, Rayleigh damping is used for the analysis, and the damping ratio corresponding to the structural frequency is taken as 5%, =0.366 α , and =0.007 β .

Ground Motion Input
The engineering foundation site classification is Class II, and the seismic intensity scale is Ⅷ. To better understand the structural vibration response, damage development process, and damage degree caused by strong earthquakes, the data from the strong earthquake database of the Pacific Earthquake Engineering Research (PEER) Centre was used as a reference to select waves with a magnitude of 6~8, with shear wave velocity of 500~750 m/s, and epicentral distance of 0~105 km. The natural wave Northridge-01 (Whittier-S. Alta Dr), which was close to the characteristic site period of 0.4 s, has been chosen, following the frequency band fixation method proposed by the 1.301 Hz Vertical 3 lar height, width, and length under wet mode conditions, this aqueduct's natural fre-quency and mode of vibration are similar to other aqueducts [33,34]. For example, the first natural frequency of an aqueduct (50 m × 14 m × 7 m) is 1.199 Hz [35]. We refer to the results of a doctoral dissertation. The aqueduct structure is similar to it in this paper (single 30 m × 20 m × 8 m) [35]. In the dynamic model experiment, the first frequency of modal identification is about 0.9~1.2 Hz, and the finite element calculation result is 0.973 Hz. Therefore, the natural frequency results of the aqueduct structure in this paper can be used to validate the accuracy of the finite element model. In this study, Rayleigh damping is used for the analysis, and the damping ratio corresponding to the structural frequency is taken as 5%, =0.366 α , and =0.007 β .

Ground Motion Input
The engineering foundation site classification is Class II, and the seismic intensity scale is Ⅷ. To better understand the structural vibration response, damage development process, and damage degree caused by strong earthquakes, the data from the strong earthquake database of the Pacific Earthquake Engineering Research (PEER) Centre was used as a reference to select waves with a magnitude of 6~8, with shear wave velocity of 500~750 m/s, and epicentral distance of 0~105 km. The natural wave Northridge-01 (Whittier-S. Alta Dr), which was close to the characteristic site period of 0.4 s, has been chosen, following the frequency band fixation method proposed by the quency and mode of vibration are similar to other aqueducts [33,34]. For example, the first natural frequency of an aqueduct (50 m × 14 m × 7 m) is 1.199 Hz [35]. We refer to the results of a doctoral dissertation. The aqueduct structure is similar to it in this paper (single 30 m × 20 m × 8 m) [35]. In the dynamic model experiment, the first frequency of modal identification is about 0.9~1.2 Hz, and the finite element calculation result is 0.973 Hz. Therefore, the natural frequency results of the aqueduct structure in this paper can be used to validate the accuracy of the finite element model. In this study, Rayleigh damping is used for the analysis, and the damping ratio corresponding to the structural frequency is taken as 5%, =0.366 α , and =0.007 β .

Ground Motion Input
The engineering foundation site classification is Class II, and the seismic intensity scale is Ⅷ. To better understand the structural vibration response, damage development process, and damage degree caused by strong earthquakes, the data from the strong earthquake database of the Pacific Earthquake Engineering Research (PEER) Centre was used as a reference to select waves with a magnitude of 6~8, with shear wave velocity of 500~750 m/s, and epicentral distance of 0~105 km. The natural wave Northridge-01 (Whittier-S. Alta Dr), which was close to the characteristic site period of 0.4 s, has been chosen, following the frequency band fixation method proposed by the 1.304 Hz Transverse 4 quency and mode of vibration are similar to other aqueducts [33,34]. For example, the first natural frequency of an aqueduct (50 m × 14 m × 7 m) is 1.199 Hz [35]. We refer to the results of a doctoral dissertation. The aqueduct structure is similar to it in this paper (single 30 m × 20 m × 8 m) [35]. In the dynamic model experiment, the first frequency of modal identification is about 0.9~1.2 Hz, and the finite element calculation result is 0.973 Hz. Therefore, the natural frequency results of the aqueduct structure in this paper can be used to validate the accuracy of the finite element model. In this study, Rayleigh damping is used for the analysis, and the damping ratio corresponding to the structural frequency is taken as 5%, =0.366 α , and =0.007 β .

Ground Motion Input
The engineering foundation site classification is Class II, and the seismic intensity scale is Ⅷ. To better understand the structural vibration response, damage development process, and damage degree caused by strong earthquakes, the data from the strong earthquake database of the Pacific Earthquake Engineering Research (PEER) Centre was used as a reference to select waves with a magnitude of 6~8, with shear wave velocity of 500~750 m/s, and epicentral distance of 0~105 km. The natural wave Northridge-01 (Whittier-S. Alta Dr), which was close to the characteristic site period of 0.4 s, has been chosen, following the frequency band fixation method proposed by the 1.214 Hz Transverse 8 lar height, width, and length under wet mode conditions, this aqueduct's natural fre-quency and mode of vibration are similar to other aqueducts [33,34]. For example, the first natural frequency of an aqueduct (50 m × 14 m × 7 m) is 1.199 Hz [35]. We refer to the results of a doctoral dissertation. The aqueduct structure is similar to it in this paper (single 30 m × 20 m × 8 m) [35]. In the dynamic model experiment, the first frequency of modal identification is about 0.9~1.2 Hz, and the finite element calculation result is 0.973 Hz. Therefore, the natural frequency results of the aqueduct structure in this paper can be used to validate the accuracy of the finite element model. In this study, Rayleigh damping is used for the analysis, and the damping ratio corresponding to the structural frequency is taken as 5%, =0.366 α , and =0.007 β .

Ground Motion Input
The engineering foundation site classification is Class II, and the seismic intensity scale is Ⅷ. To better understand the structural vibration response, damage development process, and damage degree caused by strong earthquakes, the data from the strong earthquake database of the Pacific Earthquake Engineering Research (PEER) Centre was used as a reference to select waves with a magnitude of 6~8, with shear wave velocity of 500~750 m/s, and epicentral distance of 0~105 km. The natural wave Northridge-01 (Whittier-S. Alta Dr), which was close to the characteristic site period of 0.4 s, has been chosen, following the frequency band fixation method proposed by the  Applied Technology Council (ATC-3) in 1978 [36]. The basic parameters of selected natural waves are shown in Table 3. The horizontal peak acceleration of severe earthquakes is 0.19 g, and the ratio to the vertical peak acceleration is 1:0.65. The seismic wave input interval is 0.01 s, and the total duration is 25 s. Two acceleration time−history curves are shown in Figure 3. The additional mass approach proposed by Westergaard in 1933 during a seismic study on the dam body is also used when considering the force of the seismic load on the water body [37].   Figure 4 depicts the progression of the aqueduct structure's dynamic damage as a result of significant earthquakes. The pier shaft and pier cap were severely damaged because they functioned as the principal supports. Penetrating damage occurred on both sides of the pier cap and around the pier shaft, and minor damage near the cushion cap, as indicated by the damage development range. Because of the aqueduct's high rigidity, damage only emerged near the top. The damage development pattern shows that the seismic excitation was weak at the early stages of the earthquake, and the structure was in the elastic stage, implying that damage has little impact on the result. The damage first showed on both sides of the pier shaft and at the bottom of the pier cap when the earthquake struck 3.70 s. It is due to the introduction of CDP constitutive, which describes the phenomenon of the concrete structure after damage. The structure was un-   However, due to the effect of foundations, the aqueduct body and cushion cap showed damage later than the pier shafts and pier caps, and the damage expanded slower than the aqueduct body. The damage ranges gradually increased throughout the entire aqueduct during the period 11.00 s~25.00 s, under the continuous action of the seismic force, as shown in Figure 4c,d. At the later stage of the earthquake, the damage ranges eventually stabilized when the earthquake duration curve shifted.

Damage Mechanism of Aqueduct Structures under Ground Motion
Changes in acceleration cause the structural stress state to vary during the seismic wave motion process, an important component of dynamic load time history calculation. The feature points on one side of the structure are chosen for analysis to examine the development process and laws of aqueduct structural damage during seismic wave motion. Figure 5 depicts the evolution and changes in damage value at each point. The slope of the damage development process curve of the pier cap is the largest, and it is also the position where the damage occurs first. The damage range quickly increases over the entire aqueduct construction as the seismic force continues, and the pier shaft and pier cap are damaged almost simultaneously, and entirely. Compared to the pier shaft and pier cap, the cushion cap is damaged later and has a smaller damage range.
The damage to the aqueduct structure proliferates, as evidenced by the damage development curves of characteristic points. The structural damage aggravates rapidly when the ground motion acceleration reaches its peak, indicating that the peak acceleration of ground motion can considerably impact the damage to the structure, and the seismic damage is severe. When comparing the development trend of damage to various elements, pier caps and pier shafts exceed the aqueduct body and cushion caps significantly. Since pier shafts and pier caps are responsible for sustaining the upper structure and affecting its overall stability, it is hypothesized that the rapid progression of damage to pier shafts and pier caps is due to their increased displacement under seismic forces. The above assessments show that the most susceptible parts of the aqueduct are the pier tops and pier shafts, followed by the top of the aqueduct body, which should be the focus of concern during a strong earthquake. The entire aqueduct structure may collapse if the damage and cracks are severe.
As can be seen from the above analysis results, each element of the aqueduct structure is damaged to varying degrees when subjected to the action of a rare earthquake. The arrival time and intensity of earthquakes will occur randomly in real-world engineering. Earthquakes of various magnitudes have also struck places where earlier seismic disasters have happened. The occurrence of these earthquakes of varying magnitudes is also likely to cause structural damage. So far, the fragility mechanism of aqueduct structures at various earthquake intensities has remained mostly unexplored. In the following section, we proceed to analyse the seismic fragility of an aqueduct using the IDA and MSA fragility methodologies, respectively, to verify the accuracy of fragility research results and provide a foundation for aqueduct seismic research and post-earthquake decision-making. Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 24 (a) Changes in acceleration cause the structural stress state to vary during the seismic wave motion process, an important component of dynamic load time history calculation. The feature points on one side of the structure are chosen for analysis to examine the development process and laws of aqueduct structural damage during seismic wave motion. Figure 5 depicts the evolution and changes in damage value at each point. The slope of the damage development process curve of the pier cap is the largest, and it also the position where the damage occurs first. The damage range quickly increase over the entire aqueduct construction as the seismic force continues, and the pier sha and pier cap are damaged almost simultaneously, and entirely. Compared to the pi shaft and pier cap, the cushion cap is damaged later and has a smaller damage range. The damage to the aqueduct structure proliferates, as evidenced by the damag development curves of characteristic points. The structural damage aggravates rapid when the ground motion acceleration reaches its peak, indicating that the peak acceler tion of ground motion can considerably impact the damage to the structure, and th seismic damage is severe. When comparing the development trend of damage to variou elements, pier caps and pier shafts exceed the aqueduct body and cushion caps signif cantly. Since pier shafts and pier caps are responsible for sustaining the upper structu and affecting its overall stability, it is hypothesized that the rapid progression of damag to pier shafts and pier caps is due to their increased displacement under seismic force The above assessments show that the most susceptible parts of the aqueduct are the pi tops and pier shafts, followed by the top of the aqueduct body, which should be the f cus of concern during a strong earthquake. The entire aqueduct structure may collapse the damage and cracks are severe.
As can be seen from the above analysis results, each element of the aqueduct stru ture is damaged to varying degrees when subjected to the action of a rare earthquak The arrival time and intensity of earthquakes will occur randomly in real-world eng neering. Earthquakes of various magnitudes have also struck places where earlier sei mic disasters have happened. The occurrence of these earthquakes of varying magn tudes is also likely to cause structural damage. So far, the fragility mechanism of aqu duct structures at various earthquake intensities has remained mostly unexplored. In th following section, we proceed to analyse the seismic fragility of an aqueduct using th IDA and MSA fragility methodologies, respectively, to verify the accuracy of fragili research results and provide a foundation for aqueduct seismic research an post-earthquake decision-making.

Methods to Develop Fragility Curves
To understand the importance of using the fragility functions and how scholars can be used to assess various structures against earthquakes, this section presents the fragility function types and the difference between them.
The fragility curves are an important tool to assess seismic risk. Every building or structure has its fragility curve. Based on the literature review [38], four methods to develop fragility curves can be identified, namely: (1) expert-based method; (2) empirical; (3) analytical; (4) hybrid.
(1) Expert-Based Method: The expert-based method is the oldest and simplest one among those mentioned above. This method depends on the use of questionnaires, the experiences of experts, and the number of experts consulted. In general, their judgments or opinions may contain uncertainties and may be less accurate, thus affecting the quality of the result.
(2) Empirical: The fragility curve developed by this empirical method is based on previous earthquake events. The advantage of this approach is that it shows the actual fragility to represent a realistic picture.
(3) Analytical: This approach is developed using simulated data from a time history analysis of the structural model for real. The analytical method is the most popular method in developing fragility seismic curves because this approach has less bias. The weakness of the analytical method is its requirement to produce a realistic model, that is, if the model is improperly designed or unrealistic, then it may result in inaccurate estimation that can affect the fragility curve.
(4) Hybrid: the hybrid approach is a calibrated empirical and analytical method, and is conducted by integrating numerical method to solve a numerical structural model equation. The disadvantage of the hybrid method is that it requires the combination of experiment and analysis.

Basic Principles of Incremental Dynamic Analysis (IDA)
The numerical method is commonly used in the seismic fragility analysis of a structure. That is, if a house fails to withstand an earthquake and collapses totally, the fragility probability is 1; if the house remains intact, the likelihood is 0. The relationship between the structural engineering demand measure (DM) and the seismic intensity measure (IM) in fragility analysis is Assuming the median D of the seismic demand parameter and that the seismic intensity measure is subject to an exponential relationship, which is Take logarithms of both sides of the above formula to get, where a and b are regression coefficients, Letλ d = LnD, then In the formula, λ d is the logarithmic mean of DM; and β d is the logarithmic standard deviation of DM.
Define the capacity parameter C of the structure with the logarithmic mean λ d and logarithmic standard deviation β d ; then: That is, P f = P(C − D < 0). If Z = C − D, C and D are both independent random variables and subject to normal distribution so that Z is also subject to normal distribution, and its mean is λ z = λ c − λ d and standard deviation β z = β 2 c + β 2 d . The exceedance probability P f for a certain performance level of a structure under a specific action of earthquakes can be expressed as:

Basic Principles of Multiple Strip Analysis (MSA)
According to reference [26], θ and σ are obtained through the maximum likelihood estimation (MLE) method to represent the logarithmic mean and logarithmic standard deviation of each performance state. The two-parameter (σ and θ) log-normal distribution function for the seismic fragility is: where I M = x is the seismic intensity measure, P(C|I M = x) , represents the probability of the structure reaching a certain performance criterion when the input seismic intensity I M = x, and Φ(x) represents the standard normal distribution function. With the two parameters θ and σ, the normal distribution function with the seismic intensity measure I M as an independent variable can be obtained for the fragility analysis According to the MLE principle and formula (18), the parameters θ and σ are estimated thus: Figure 6 illustrates the process for generating fragility curves. The engineering demand measure and the seismic intensity measure for aqueduct structure are estimated from the IDA, and the double parameter values of performance level at each limit state are calculated from the MSA.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 13 of 24 Figure 6. Illustration of the steps for generating fragility curves.

Selection of Ground Motions
Luco et al. [39] discovered that 10~20 seismic wave records might be used to estimate the structural earthquake-resistant capabilities of mid-to high-rise buildings correctly. Refer to Section 2.3 of this study for seismic wave selection criteria, and use the duration to attain 1/10 of the peak acceleration as the threshold [40]. Fifteen natural seismic waves from the database were chosen, along with two artificial waves that matched the engineering site spectrum. Table 4 lists the 17 seismic records, whereas Figure 7 depicts the response spectrum. The red line represents the target response spectrum.

Selection of Ground Motions
Luco et al. [39] discovered that 10~20 seismic wave records might be used to estimate the structural earthquake-resistant capabilities of mid-to high-rise buildings correctly. Refer to Section 2.3 of this study for seismic wave selection criteria, and use the duration to attain 1/10 of the peak acceleration as the threshold [40]. Fifteen natural seismic waves from the database were chosen, along with two artificial waves that matched the engineering site spectrum. Table 4 lists the 17 seismic records, whereas Figure 7 depicts the response spectrum. The red line represents the target response spectrum.

Selection of Ground Motions
Luco et al. [39] discovered that 10~20 seismic wave records might be used to estimate the structural earthquake-resistant capabilities of mid-to high-rise buildings correctly. Refer to Section 2.3 of this study for seismic wave selection criteria, and use the duration to attain 1/10 of the peak acceleration as the threshold [40]. Fifteen natural seismic waves from the database were chosen, along with two artificial waves that matched the engineering site spectrum. Table 4 lists the 17 seismic records, whereas Figure 7 depicts the response spectrum. The red line represents the target response spectrum.

Determination of Structural Performance Indicators
Studies have shown that PGA, PGV, and spectral acceleration S a (T 1 , 5%) can serve as common ground motion intensity measures [41][42][43]. Yu [44] chose 60 ground motion IM related to structural performance to construct ground motion parameter evaluation models. He believes that the analysis results involving the spectral acceleration could accurately describe the probability statistic relationship between IM and DM and avoid dispersion between analysis results. The aqueduct structure is similar to bridge structures in terms of engineering design. The pier shaft is responsible for sustaining the superstructure and affecting the overall stability of both structures. Bridge damage can be measured using a variety of indicators, including pier top displacement, pier top offset ratio, cracking degree, curvature ductility ratio, and displacement ductility ratio. Dutta and Mander [45] proposed the pier top offset ratio during their research, using it to indicate the pier shaft's deformation capacity. It is the ratio of the maximum displacement of the pier top to the pier shaft height during the earthquake. The limit values are shown in Table 5. The joint rotation features, structural damage degree, and deformation performance in the dynamic analysis are all directly related to displacement in the structure. Therefore, referring to the four limit states divided by FEMA356 and combined with the existing achievement of previous scholars [45][46][47], this paper divides structural earthquake-resistant performance into four grades: normal occupancy (NO), immediate occupancy (IO), life safety (LF), and collapse prevention (CP). We use the S a (T 1 , 5%) as ground motion IM and the aqueduct pier shaft top offset ratio to quantify the structural damage index limit to reflect the specific damage degree to a structure in case of exceedance probability. The limit state values of the quantitative indicators of the structural performance level can be found in Table 5. Minor repair is needed for structural and non-structural members.
The structure remains stable and has enough capacity.
The structure does not collapse and the damage is acceptable.
Note: dr = U max /H, dr is the aqueduct pier top offset ratio, U max is the maximum displacement of aqueduct pier shaft top, and H is the height of aqueduct pier shaft, H = 5.5 m.

Criteria for Amplitude Modulation of Seismic Waves
In this paper, the hunt and fill amplitude modulation criterion proposed by Vamvatasikos and Jalayer is used to balance the calculation accuracy and efficiency [48]. Using the natural wave Northridge-01 (Whittier-S. Alta Dr) as an example, we used the damping ratio of 5% to obtain the S a (T 1 , 5%) = 0.005 g at the first amplitude modulation, and a step increment of 0.05 g. The structure's performance evolution law under ground motion excitation was determined progressively. Table 6 depicts the exact amplitude modulation process.  Figure 8 depicts the results of a nonlinear analysis of the seismic dynamic response of an aqueduct structure under earthquake activity. When the IM remains constant, the data points in the graph reflect the maximum displacement of the pier shaft top of an aqueduct construction against different seismic waves. The quantitative indications of the four limit structural damage states are represented by the red dotted line. LS1 represents the limit state of the aqueduct structure changing from no damage to normal occupancy, LS2 represents the limit state of the aqueduct structure changing from normal occupancy to immediate occupancy, LS3 represents the limit state of the aqueduct structure changing from immediate occupancy to life safety, and LS4 represents the limit state of the aqueduct structure from life safety to collapse prevention.

IDA Method-Based Probability Demand Analysis
Appl. Sci. 2021, 11, x FOR PEER REVIEW 1 Figure 8 depicts the results of a nonlinear analysis of the seismic dynamic res of an aqueduct structure under earthquake activity. When the IM remains constan data points in the graph reflect the maximum displacement of the pier shaft top aqueduct construction against different seismic waves. The quantitative indicatio the four limit structural damage states are represented by the red dotted line. LS1 r sents the limit state of the aqueduct structure changing from no damage to norm cupancy, LS2 represents the limit state of the aqueduct structure changing from n occupancy to immediate occupancy, LS3 represents the limit state of the aqu structure changing from immediate occupancy to life safety, and LS4 represents the state of the aqueduct structure from life safety to collapse prevention. Since a single IDA data cannot accurately predict the limit behavior of the stru this paper selects a series of seismic records to perform the incremental dynamic a sis and obtain multiple IDA data. The 17 seismic ground records IDA data gene during ground motion excitation were summarized in the same IM-DM [(Umax-5%))] coordinate system, as shown in Figure 8. It has been conducted using the method as obtaining the IDA curve during ground motion excitation for the n wave Northridge-01 (Whittier-S. Alta Dr). The response of the same structure to d ent seismic records input is discrete to some extent since the shape of the IDA cu directly related to the selected seismic records. Therefore, it is necessary to summ these different curves with the median and the 16% and 84% percentile curves to r the difference. The three IDA percentile curves obtained indicate that among t seismic records, the 16%, 50%, and 84% seismic ground motion records exceede corresponding structural limit, as shown in Figure 9. The corresponding capacity v of each limit state of the structure are shown in Table 7. Since a single IDA data cannot accurately predict the limit behavior of the structure, this paper selects a series of seismic records to perform the incremental dynamic analysis and obtain multiple IDA data. The 17 seismic ground records IDA data generated during ground motion excitation were summarized in the same IM-DM [(Umax-S a (T 1 , 5%))] coordinate system, as shown in Figure 8. It has been conducted using the same method as obtaining the IDA curve during ground motion excitation for the natural wave Northridge-01 (Whittier-S. Alta Dr). The response of the same structure to different seismic records input is discrete to some extent since the shape of the IDA curve is directly related to the selected seismic records. Therefore, it is necessary to summarize these different curves with the median and the 16% and 84% percentile curves to reduce the difference. The three IDA percentile curves obtained indicate that among the 17 seismic records, the 16%, 50%, and 84% seismic ground motion records exceeded the corresponding structural limit, as shown in Figure 9. The corresponding capacity values of each limit state of the structure are shown in Table 7.

IDA Method-Based Probability Demand Analysis
Taking the spectral acceleration S a (T 1 , 5%) as the independent variable IM and the maximum displacement of aqueduct pier shaft top as the dependent variable, the functional relation between them is as follows: Ln(U max ) = a + bLn(S a (T 1 , 5%)) (22) In the formula, a and b are regression coefficients. We performed a logarithmic regression analysis with the data in Figure 10, with a correlation coefficient R = 0.998. According to the seismic design standard HAZUS99 for building fragility curve parameters [49], when the fragility curve takes the spectral acceleration S a (T 1 , 5%) as the independent variable,  Taking the spectral acceleration Sa(T1, 5%) as the independent variable IM and the maximum displacement of aqueduct pier shaft top as the dependent variable, the functional relation between them is as follows: Ln U a bLn S T = + (22) In the formula, a and b are regression coefficients. We performed a logarithmic regression analysis with the data in Figure 10, with a correlation coefficient R = 0.998. According to the seismic design standard HAZUS99 for building fragility curve parameters [49], when the fragility curve takes the spectral acceleration Sa(T1, 5%) as the independent variable,    Figure 11 shows the seismic fragility curves of the aqueduct structure calculated based on the IDA method. For the limit state LS1, the IDA method obtains that the ex  Figure 11 shows the seismic fragility curves of the aqueduct structure calculated based on the IDA method. For the limit state LS1, the IDA method obtains that the exceedance probability of the structure is 64.18% at 0.11 g and 96.77% at 0.23 g. In the limit state LS2, the exceedance probability of the structure at 0.23 g is 62.50%, and the exceedance probability at 0.51 g is 97.10%. According to the IDA method, when the structure is 0.51 g, the exceedance probability in LS3 is close to 80%. It implies that if the IM is greater than 0.50 g, the influence of weak structure components, such as the pier body and pier cap, should be adequately evaluated, as they are highly likely to affect structural safety. When the value exceeds 1.00 g, the probability is gradually close to 100%. In the LS4, the results show that the exceedance probability value of the structure is close to 90% at 1.10 g. The fragility analysis curve of the LS1 of the aqueduct structure under IDA conditions is highly steep and close to the vertical line, making it very simple to surpass the structural damage index of the LS1, according to studies on the characteristics of fragility curves. The LS2 has a somewhat steep fragility curve, while LS3 and LS4 gradually slow down. The IDA method is a traditional method to calculate the seismic fragility analysis curve of large buildings. The calculation of this method mainly depends on the logarithmic mean and logarithmic standard deviation of structural engineering demand parameters.

Determination of MSA-Based Performance Indicators of Each Limit State
Baker [50] found that the seismic data obtained using the IDA approach can likewise be used to generate seismic fragility curves using the MSA method. Based on Figure 10 and the basic principles of MSA, the double parameter values θ and σ of each limit state performance level are obtained, as shown in Table 8.   Figure 12 shows the seismic fragility curves of the aqueduct structure calculated based on the MSA method. For limit state LS1, the exceedance probability of the structure is 48.76% at 0.11 g and close to 100% at 0.40 g based on the MSA method. The results show that the aqueduct structure selected in this paper is likely to produce small-scale damage when the IM exceeds 0.40 g. The fragility curve of the LS1 of the aqueduct structure is still very steep. In the LS2, the exceedance probability of the structure at 0.23 g is 51.96%, and at 0.51 g, the exceedance probability has reached 91.68%. It means that when the Sa(T1, 5%) reaches 0.50 g, structures such as the pier shaft and pier cap are prone to damage. The pier cap and shaft of the aqueduct need repair, likely. In the limit state LS3, when the IM reaches 1.10 g, the exceedance probability is 95.80%, showing that the earthquake will cause serious damage to individual structural parts,

Determination of MSA-Based Performance Indicators of Each Limit State
Baker [50] found that the seismic data obtained using the IDA approach can likewise be used to generate seismic fragility curves using the MSA method. Based on Figure 10 and the basic principles of MSA, the double parameter values θ and σ of each limit state performance level are obtained, as shown in Table 8.  Figure 12 shows the seismic fragility curves of the aqueduct structure calculated based on the MSA method. For limit state LS1, the exceedance probability of the structure is 48.76% at 0.11 g and close to 100% at 0.40 g based on the MSA method. The results show that the aqueduct structure selected in this paper is likely to produce small-scale damage when the IM exceeds 0.40 g. The fragility curve of the LS1 of the aqueduct structure is still very steep. In the LS2, the exceedance probability of the structure at 0.23 g is 51.96%, and at 0.51 g, the exceedance probability has reached 91.68%. It means that when the S a (T 1 , 5%) reaches 0.50 g, structures such as the pier shaft and pier cap are prone to damage. The pier cap and shaft of the aqueduct need repair, likely. In the limit state LS3, when the IM reaches 1.10 g, the exceedance probability is 95.80%, showing that the earthquake will cause serious damage to individual structural parts, affecting the safety and stability of the aqueduct structure.

Comparison of Fragility Analysis Results of Aqueduct Structure
The comparison result of the two fragility curves based on the IDA and MSA methodologies is shown in Figure 13. Probability under various earthquake intensities is shown in Table 9. It can be seen that for LS1, the exceedance probability of the structure with IM = 0.11 g, the maximum difference of two fragility curves is about 18.50%, and their tendency is roughly in line. For LS2, the exceedance probability with IM = 0.31 g obtained by the two fragility methods is 81.36% and 69.92%, respectively. The IDA-based fragility curve has a steeper upward trend, with a difference of 11.44% between the two methods, and the two curves are in good agreement. The maximum difference for LS3 is about 12.50% for IM = 0.41 g. The maximum difference between LS4 and LS5 is 17.63%. The LS1 of the two methods shows a steep trend, and the fragility curves of the LS2 (slight damage) are also steep until the fragility curves of the LS3 and LS4 gradually slow down.

Comparison of Fragility Analysis Results of Aqueduct Structure
The comparison result of the two fragility curves based on the IDA and MSA methodologies is shown in Figure 13. Probability under various earthquake intensities is shown in Table 9. It can be seen that for LS1, the exceedance probability of the structure with IM = 0.11 g, the maximum difference of two fragility curves is about 18.50%, and their tendency is roughly in line. For LS2, the exceedance probability with IM = 0.31 g obtained by the two fragility methods is 81.36% and 69.92%, respectively. The IDA-based fragility curve has a steeper upward trend, with a difference of 11.44% between the two methods, and the two curves are in good agreement. The maximum difference for LS3 is about 12.50% for IM = 0.41 g. The maximum difference between LS4 and LS5 is 17.63%. The LS1 of the two methods shows a steep trend, and the fragility curves of the LS2 (slight damage) are also steep until the fragility curves of the LS3 and LS4 gradually slow down.
In conclusion, the seismic fragility curves generated using the IDA and MSA methodologies mutually verify the aqueduct structure's seismic fragility results. The main reason for the large difference in the probability of individual damage is that MSA-based calculations are limited by the number of seismic waves exceeding the limit. In contrast, IDA-based calculations are based primarily on the mean and logarithmic standard deviation of the structural engineering demand parameter. However, their conclusions are largely in agreement, indicating that in the seismic analysis, the aqueduct structure will be difficult to maintain no damage, and it is very easy to exceed the damage index of the first limit state level. The IDA approach is very straightforward in theory in terms of calculation time, but it involves extensive calculations and complex data processing in practice. On the other hand, the MSA technique can provide the exceedance probability at any IM without complex work. In seismic calculations, IDA-based fragility analysis may fail to appropriately calculate the IM value corresponding to each limit state when the probability of the structure surpassing the damage index does not strictly increase with the IM positively. On the other hand, the MSA approach has no such flaws and can obtain multiple horizontal bands of IM-DM ground motion. The fragility curve is more precisely calculated using statistical analysis. at any IM without complex work. In seismic calculations, IDA-based fragility analysis may fail to appropriately calculate the IM value corresponding to each limit state when the probability of the structure surpassing the damage index does not strictly increase with the IM positively. On the other hand, the MSA approach has no such flaws and can obtain multiple horizontal bands of IM-DM ground motion. The fragility curve is more precisely calculated using statistical analysis.

Conclusions
The dynamic damage mechanism and damage development law of an aqueduct under the impact of a rare earthquake are investigated in this research and the weak spots of the aqueduct structure where concrete cracking may occur. Seventeen seismic waves are selected based on the factors such as magnitude, shear wave velocity, epicentral distance, and site type. We use the IDA and MSA methods to analyze the seismic fragility of the aqueduct, respectively, with Sa(T1, 5%) as the IM and the pier shaft top offset ratio under structural response to quantify the index limit, which mutually verifies the rationality of the fragility results and studies the probability of damage to the structure under different IM. The aqueduct selected in this paper is universal in the SNWDP.

Conclusions
The dynamic damage mechanism and damage development law of an aqueduct under the impact of a rare earthquake are investigated in this research and the weak spots of the aqueduct structure where concrete cracking may occur. Seventeen seismic waves are selected based on the factors such as magnitude, shear wave velocity, epicentral distance, and site type. We use the IDA and MSA methods to analyze the seismic fragility of the aqueduct, respectively, with S a (T 1 , 5%) as the IM and the pier shaft top offset ratio under structural response to quantify the index limit, which mutually verifies the rationality of the fragility results and studies the probability of damage to the structure under different IM. The aqueduct selected in this paper is universal in the SNWDP. However, the aqueduct is located in an earthquake-prone area. The research results apply to the aqueduct, bridges, and other structures similar to the structural model and conditions in this paper. The following are the main conclusions:

•
Penetrating damage is most likely to occur on both sides of the pier cap and around the pier shaft in the event of a rare earthquake, followed by the top of the aqueduct body, which will require the most concern during the earthquake. The entire aqueduct structure may collapse once the damage and cracks get severe. Therefore, in the process of seismic design, relevant damping measures should be taken to control the damage development and prevent the overall aqueduct structure from being damaged or collapsed. For easily damaged parts, appropriate structural reinforcement is required to meet the needs of earthquake resistance. • Two fragility analysis methodologies are used to investigate the aqueduct structure separately. The results reveal that the fragility curves obtained by IDA and MSA are very similar. Except for individual results, the variation between fragility curves of different LS is less than 10%, demonstrating that the fragility analysis results are rational. • Each of the two fragility methodologies has its advantages. The IDA method is simple in theory, but data processing is complicated because the logarithmic mean and logarithmic standard deviation for each seismic IM must be determined. In theory, the MSA method is more complicated, but the probability in the case of any IM may be calculated by acquiring only two parameters that correspond to the LS. In terms of calculation, it is more efficient.

•
The seismic fragility analysis of aqueduct structure is the main research content in this paper. Although the fragility analysis of different methods is realized, this paper only selects one kind of spectral acceleration and the pier shaft top offset ratio under structural as the standard. The calculation index should be added to obtain the most suitable strength index and response index in the fragility study of the aqueduct, which will also be the next work of this research.

•
The building needs need professional seismic hazard analysis, and the aqueduct structure is no exception. Combining with the research results of seismic fragility in this paper, constructing seismic hazard curves to jointly evaluate the seismic safety of the aqueduct structure is the next step of the authors.