Investigation of Stochastic Seismic Response and Index Correlation of an Arch Dam Using Endurance Time Analysis Method

,e seismic dynamic response and correlation between various indicators of an arch dam under a series of stochastic earthquake excitation are analysed. Seismic response assessment of a concrete arch dam is conducted through nonlinear endurance time analysis (ETA) of a typical arch dam.,e 289m high Baihetan arch dam in Southwest China (as a case study) is subjected to a set of 10 three-component endurance time acceleration histories (ETAHs) synthesised on the basis of the design site acceleration response spectrum to increasing peak ground acceleration (PGA). ,e effects of dam-reservoir-foundation interaction, radiation damping of infinite rock foundation, and nonlinearities of dam material and joint contact elements are considered. ,e indices to failure modes of damage propagation through the dam body and the deformation and joint opening beyond the limit states of the arch dam are identified. Dynamic responses, including deformation, joint opening, and damage distribution, are statistically analysed. ,e dispersion and correlation of indicators are analysed by employing the coefficient of variation and coefficient of correlation, respectively. ,e deterministic relationships between contraction joint opening, displacement, and damage volume ratio are obtained by using a multivariate fitting method. ,e results show that those failure indices have similar exponential trends, and seismic intensity parameters exert a significant influence on seismic response and damage under increasing PGA.,e results of correlation and fragility analysis demonstrate that damage volume ratio and the sum of joint opening are consistent, which can provide a strong scientific basis for predicting earthquake damage.


Introduction
e seismic safety evaluation and performance analysis are the primary considerations when assessing the safety margin of a structure during an earthquake [1,2], and it also becomes one of the most significant concerns among hydraulic structural engineers worldwide. e dynamic response of a structure has uncertainty deriving from the randomness of ground motions, which will have an influence on the seismic capacity of the structure [3]. In particular, the damage response can reflect the characteristics of the damage distribution over time, resulting in eventual global instability; however, for an arch dam, challenges lie in the complexity of the dam-reservoir-foundation system and the extensive distributions of certain parameters so that specific failure modes are unclear in previous studies. is warrants further research to understand the severity of damage and realise fully the behaviours manifested in an arch dam subjected to earthquake loading.
Currently, several researchers investigated the input approach to an arch dam-foundation system [4] and explored the damage thereto by defining different dynamic responses [5] when subjected to strong ground motions. Various seismic evaluation indices were presented to estimate structural seismic performance accounting for structural nonlinearity, such as displacement, damage, and contraction joint opening indices [6]. Xu et al. [7,8] performed the dynamic analysis of an arch dam to analyse the stress and contraction joint opening considering different input models. Hariri-Ardebili and Kianoush [9,10] developed a seismic safety evaluation for arch dams and concrete gravity dams to simulate crack propagation behaviour. Ren et al. [11] proposed the use of a mixed finite element method for Shuanghe arch dam to simulate crack propagation therein and verify actual measured data with reasonable precision. Wang et al. [12,13] developed a comprehensive procedure with which to analyse the seismic nonlinear behaviour of Pacoima and Dagangshan arch dams based on an arch dam-water-foundation rock system. Pan et al. [14] described a nonlinear damage system for featuring damage-crack behaviour under earthquake loading using both a massless foundation input model and a viscousspring boundary model. Other studies (Karalar and Cavusli [15][16][17][18][19], Kartal et al. [20,21], and Demirtas et al. [22]) revealed the effects of ground motions on the seismic response and the damage assessment of structures. Most current studies only consider the dynamic response and damage extension under seismic ground motion, but research into the correlation of indicators remains sparse. In addition, the mechanism underpinning the dynamic behaviour of an arch dam is rarely analysed under strong ground motion and studies seldom consider both material nonlinearity and the verification of ground motion input to guarantee accuracy of the results.
Moreover, it is worth noting that many researchers paid more attention to seismic performance analysis methods considering the hazards of structural failure due to the complexity and bulk of the structure. In terms of structural analysis methods, incremental dynamic analysis (IDA) [23,24], pushover analysis [25], and multiple stripe analysis [26] are widely used. Alembagheri et al. [27,28] adopted the IDA method to analyse the damage propagation behaviour of both a concrete gravity dam and an arch dam and chose different indices to compare dynamic response requirements and structural capabilities with increasing earthquake intensity. Wang et al. [29] analysed a dynamic damage model of the Dagangshan arch dam-reservoir-foundation system to illustrate the damage characteristic and subsequently determined the seismic fragility curves under three limit states by using the IDA method. Pang et al. [30] evaluated the seismic resistance of high concrete face rockfill dams using the IDA method regarding different damage indices; however, it should be noted that these methods mainly established the relationship between dynamic responses and earthquakes at different intensities with high accuracy, but the computational burden and solution time were onerous. In the present research, this burden is a key issue considered while developing a method which can satisfy the demands imposed by analysis of different intensities and response spectra. Estekanchi et al. [31,32] modelled a new seismic performance procedure called endurance time analysis (ETA) that can express different intensities and spectra to measure structural health. Valamanesh et al. [33] compared the seismic behaviour, dynamic responses, and crack tracing of Folsom and Koyna dams in the endurance time acceleration history (ETAH) input with real earthquakes to a reasonable degree of approximation. Hariri-Ardebili et al. [34,35] used an ETA method for steel moment frames and concrete arch dam to assess structural seismic power effects in comparison to other methods to verify the effectiveness and precision of the ETA method. Other researchers (Hariri-Ardebili et al. [36,37], Riahi et al. [38,39], and Estekanchi et al. [40]) also made many probabilistic damage analysis of structures by using the ETA method. e key to the ETA method is primarily to generate a set of ETAHs, which are selected as ground motion inputs indicating different levels of intensity. Due to the fact that time history in the ETA method corresponds to several different sets of the peak ground acceleration and response spectra, the amplitude modulation results can be obtained by a single endurance time acceleration history, which considerably saves computing time compared with the IDA method. Although the IDA method can more accurately describe the seismic performance of structures, the dramatically high computational effort leads to the limitation of the use of this method. In addition, the results of reference [39] show that the ETA method is a good tool to appraise IDA results with reasonably good accuracy using a few number of timehistory analyses. Furthermore, the scatter of results by the ETA method is considerably less than real earthquakes used in IDA analysis especially in the state of dynamic instability. erefore, it can be seen that the ETA method is an effective method, with reasonable precision, and can thus be used as the input acceleration history of the base of dam in the present work.
In engineering practice, seismic acceleration histories are synthesised on the basis of the design site acceleration response spectrum. Seismic safety for dams can be derived on the basis of these artificial accelerograms inputs; however, the discrete nature of the dynamic results obtained by this method is key to the reliability of the method. In this work, a seismic performance ETA method suggests that different times correspond to the different acceleration response spectra via a scale factor, and the numerical optimisation method is adopted to generate ETAHs. e seismic fragility study of the Baihetan arch dam-reservoir-foundation system is presented using sets of ETAHs. e objective is the investigation of a case study that models the radiation damping of infinite rock foundation, contraction joint opening, and damage cracking of the concrete in the dam, to determine its seismic performance under different PGAs. e increasing intensity in earthquake ground motions is taken into account using the ETA method. e dam is modelled using finite elements along with its reservoir and foundation in 3D space. e plastic-damage model proposed by Lee and Fenves [41] is utilised for nonlinear analyses and determination of damage propagation in the dam body.
ree components of ETAHs corresponding to the site response spectrum for increasing PGAs are simultaneously input. e dynamic responses of displacement, damage cracking, contraction joint opening, and so on for different failure modes are explored considering the randomness of ETAHs and the uncertainty of these responses of the dam, and these are analysed on the basis of their coefficients of variation and correlation. e relationship between each index is elucidated. 2 Advances in Civil Engineering

Seismic Response Analysis
e objective of seismic response analysis of dams is to determine failure mode, to propose fragility indices, and to analyse probability distributions when dams are subjected to a given seismic intensity measure.

ETA Method.
e PGA and acceleration response spectrum are two important factors for seismic intensity measure. ETA is a time history-based dynamic pushover procedure in which special predesigned intensifying acceleration functions are used for analysis and estimation of structural responses at various performance levels by using just a single analysis. us, an endurance time acceleration history (ETAH) generated by an unconstrained optimisation method can consider these two factors (PGA and acceleration response spectrum) simultaneously for different seismic intensity measures.
In the ETA method, an ETAH is used to replace a series of seismic histories with different acceleration peaks. As shown in Figure 1, the ETA method uses different interval times to satisfy the linear growth in design acceleration response spectrum with the factor t/t target . e relationship among the target acceleration response spectrum S aT (T, t), the target displacement response spectrum S uT (T, t), and the codified design acceleration response spectrum S aC (T) is where t target is the target duration of ETAH, T is the period of free vibration, S aC (T) is the design maximum acceleration spectrum, and S aT (T, t) and S uT (T, t) are the target acceleration response spectrum and the displacement spectrum of interval times of ETAH from 0 s to t. e ETAH is optimised by the ETA method, using an unconstrained optimisation formula in the time domain as follows: where α is the weighted parameter. Here, the displacement spectrum is not considered so α � 0. Figure 1 shows the incremental relationship of acceleration response spectrum to ETAH in the ETA method.

Seismic Dynamic Indices.
In an ETA analysis, a series of strong ground motions are applied to the arch dam to assess the resistance of the structure. After applying increasing motion, the dam will change from a linear elastic state of stress to elastoplastic and ultimate collapse failure. e appropriate indices should be chosen to predict and monitor the crack propagation over its life-span in advance in case of sudden failure of the dam. In this way, the degree of damage and dynamic characteristic of the structure can be controlled within a certain range to ensure an adequate safety margin, maintain good performance, and achieve an efficient equilibrium state.
In general, arch dams mainly exhibit three failure modes: damage propagation through the dam body and deformation and joint opening beyond the limit states ( Figure 2). us, this work mainly focuses on these three failure modes in terms of seismic fragility analysis. Furthermore, distinctive dynamic responses including damage volume ratio, contraction joint opening, and radial and vertical deformation are used as the appropriate indices for assessing the performance of an arch dam. e stream displacement at the crest centre is also used to evaluate the deformation, and the maximum and the sum of contraction joint opening are used to evaluate the joint opening. e damage volume ratios used to evaluate the damage from cracking are defined as where d t is the tensile damage at each individual element and v e is the volume of the element of the dam.

Probability
Distributions. e fragility curves are assumed to take the form of a two-parameter lognormal distribution function [42][43][44]: where P LS (y) is the probability that a ground motion intensity y will result in the limit state, Φ is the standard normal distribution, and m and β are the mean and logarithmic standard deviation of y, respectively. m and β are determined by the least squares regression.

Seismic Analysis
Step. is investigation presents the seismic analysis of the Baihetan arch dam, a 289 m high double curvature arch dam in Southwest China. Figure 3 summarises the step in the analysis: modelling of the damwater-foundation system, synthesising ETAHs to the design site acceleration response spectrum and PGA for increasing seismic intensity measure levels, nonlinear seismic analyses, construction of fragility indicators and fragility curves, and analysing dispersion and correlation between indicators.

Constitutive Damage Model for Concrete and Nonlinear
Contraction Joint. Based on the plastic-damage theory Advances in Civil Engineering 3 proposed by Lee and Fenves [41], for concrete, the stressstrain relationship under uniaxial tensile load is as shown in Figure 4.
A contact boundary model is used to simulate the contraction joint [14,[45][46][47].    surface is associated with a unique anchor point M′ on the master surface. e relationship between the slave node M and its corresponding anchor point M′ is used to analyse the interaction between the slave and master surfaces. e tangential relative displacement s is assumed to be zero. e tangential displacements between the two surfaces of the element are not allowed due to the shear key action at the contraction joints in the dam. e overlap of node M normal to the master surface is denoted by l. Figure 5(b) shows the normal pressure-overlap constitutive relation, which is used to prevent the slave node from overlapping with the master surface.

Viscous-Spring Artificial Boundary Condition Input.
To simulate the radiation damping of the infinite mass rock, the dynamic artificial boundaries (such as transmission boundary, viscous boundary, and viscous-spring boundary) are used at the surface of the sufficient range foundation. Compared with the viscous boundary (also called quiet boundary) and the free-field input model generally used in software FLAC3D, the viscous-spring artificial boundary model, which not only can dissipate extroverted waves on the boundary but also simulates an elastic restoration effect of the faraway foundation, is generally employed in software ABAQUS. In addition, a typical numerical example is provided to illustrate the accuracy of the viscous-spring artificial boundary model.

Viscous-Spring Artificial Boundary Condition Input
Method.
e viscous-spring artificial boundary condition [8] is applied as a 3D space seismic input model. e seismic input method corresponding to the viscous-spring boundary is determined by several researchers, and the accuracy and applicability thereof are verified. With the help of the spring and the damping element in the finite element analysis software, springs and dampers are arranged in parallel in the normal and tangential directions of the artificial boundary nodes. e parameters of springs and dashpots of node m on Advances in Civil Engineering the viscous-spring artificial boundary are obtained by following equations for solving seismic input problems.
In the normal direction, In the tangential direction, where n and s refer to the normal and tangential directions of the boundary plane, K is the elastic stiffness of the spring, C is the viscous damping, G is the Lamé constant, c p is the propagation velocity of a P-wave, c s is the S-wave velocity, ρ is the mass density, and r is the distance from the wave source, which takes the approximate value of the perpendicular distance from the centre of the structure to the nodes on the artificial boundary. e recommended values of parameters α n and α s are 1.33 and 0.67, respectively, according to Liu et al. [48].
Here, the explicit finite element method combined with time-domain integral analysis of a viscous-spring boundary is used. It is assumed that the displacement-time histories of plane p-wave and s-waves vertically incident at the bottom boundary are u P (t) and u S (t), respectively, and the equivalent force exerted on each artificial boundary node m is calculated analytically according to the following formula and schematic diagram: where u (x m , y m , z m , t) and σ (x m , y m , z m , t) denote the displacement and stress fields at the node m in the artificial boundary model and σ 0 (x m , y m , z m , t) is the stress field caused by the free wave field in the semi-infinite space.

Verification of Viscous-Spring Artificial Boundary
Condition Input Method. MATLAB ™ is used to realise a model incorporating viscous-spring boundary theory and the corresponding equivalent force in the direction at each artificial boundary node m. e rationality and correctness of the implementation of this viscous-spring boundary and the input method of ground motion herein are verified. It is worth noting that the appropriate time steps Δt combined with finite element mesh size Δx should be taken to ensure the stability and accuracy of the solution in numerical analysis. With regard to the near-field fluctuation problem, the use of the following formula is proposed [49]: where λ T is the minimum wavelength. In general, the time step is as follows: where c is the wave velocity of seismic propagation in the medium.
In this example, the 3D free-field model employing SV wave motion is used to verify the accuracy of the wave simulation process. To meet the requirements of finite element accuracy as mentioned above, a cubic solid element model with a side length of 180 m and an element size of 6 m is adopted. e time integral step Δt � 0.005 s for the seismic wave is used in this validation. e aforementioned model can meet the requirement of calculation accuracy. For the free-field model material, Young's modulus is 335 MPa, Poisson's ratio is 0.3, and the density is 1800 kg/m 3 . e six observed points of bottom and top surface are selected to extract the time history of displacement, as shown in Figure 6(a). e time history of displacement P (τ) imposed on the bottom model is expressed as follows [50]: is the Heaviside ladder function, and the total duration of the input pulse T 0 is 1 s. e time history of the input Dirac displacement pulse function is shown in Figure 6(b). e time history of displacement responses of observed points is depicted in Figure 7, and the numerical and analytical solutions of free-field objects are also compared. It can be seen that the comparison shows a good consistency between numerical and analytical solutions. e maximum displacement of top surface is almost twice that of the input displacement history, and the time delay effect can be properly reflected. erefore, it can be concluded from the above that the compiled program can be further applied to the analysis of an arch dam.

Baihetan Finite Element Model and Seismic Input.
e Baihetan arch dam (289 m in height and 715 m along the crest) is analysed: as shown in Figure 8, the maximum thicknesses at the crest and base of the dam crown cantilever are 14 m and 63 m, and the corresponding water level varies from 545 m at base to 834 m at the crest of dam above sea level, respectively. e Baihetan dam-reservoirfoundation system is composed of an arch dam, a rock foundation with radiation damping. e far boundary of the foundation is at a distance from the dam of approximately 1.5 times the height of dam in all directions. e nonlinearities of the dam material and joint contact are also considered. e normal depth of reservoir water is 825 m above sea level. Finite element models, consisting of 25,776 hexagonal solid elements, are employed in the simulation of dam-foundation system. e hexagonal eight-node element is mainly used, while the prism element is used in the irregular part at the portions of the model. e arch dam body and foundation are discretized into 4,656 solid elements and 21,120 solid elements, respectively. To simulate the damage process, the average mesh size of the 3D dam body is approximately 0.6 m. Moreover, to express the contact mechanism and simulate the possible opening and slippage of the dam, the 30 contraction joints are incorporated in the model, thus dividing the dam into 31 blocks ( Figure 9). For each joint, the initial contact distance c � 0.005 mm and the pressure value p 0 � 50 MPa at the zero opening are selected based on Wang's article [46].
For the concrete in the dam, the density is 2400 kg/m 3 , Young's modulus is 36 GPa, and Poisson's ratio is 0.167. For the rock foundation material, the density, the elastic modulus, and Poisson's ratio are 2800 kg/m 3 , 26 GPa, and 0.24, respectively. e tensile strength f t of the concrete used in the analysis is 3.6 MPa. e fracture energy G F is assumed to be 200 N/m based on the parameters in [10,51]. e limiting tensile strain ε f is set to 0.0003, and the characteristic length l c is 0.6 m. e damping of the concrete dam system is assumed to be of Rayleigh type with a damping ratio of 5%. e modal frequencies of the Baihetan model are listed in Table 1. According to the research of Hariri-Ardebili and Kianoush [9], the modal responses of arch dam-reservoirfoundation system are determined, which suggested the first and tenth natural vibration modes as the effective range. en, the Rayleigh vibration parameters can be calculated from the first and tenth natural frequencies and damping ratio. e main static load on the high arch dam is the selfweight of the dam body and the hydrostatic pressure. e sliding at the dam-foundation rock contact interface is not considered in the numerical analysis. e endurance time acceleration histories are generated by the ETA method. e viscous-spring artificial boundary is applied on the outermost boundary of the rock foundation. In the dynamic calculations, the interaction of the dam-reservoir, namely, the action of the hydrodynamic pressure due to earthquakes, is simulated by the Westergaard added mass technique [52] at the upstream surface nodes. e added mass of node l on the upstream of dam during an earthquake is given as follows [53]: where ρ is the water density, H 0 is the water depth of the reservoir, h is the depth above node, and A l is the tributary area of node l with n being the normal direction of the curved surfaces containing node l. e diagonal terms of added mass of node l are then obtained as follows: m lx � m ln cos(n, x), m ly � m ln cos(n, y), m lz � m ln cos(n, z), (14) where m li is the added mass component in i-indirection (i � x, y, z) and cos (n, i) is the direction cosines of the normal direction n of the curved surfaces.
Ten sets of ETAHs are synthesised based on the Baihetan arch dam site acceleration response spectrum: these have linearly increasing PGAs, as shown in Figure 10  Advances in Civil Engineering 7 three-directional ETAHs, in which the X, Y, and Z-directions represent the cross-stream, stream, and vertical directions, respectively. According to the current hydraulic standard, the peak acceleration in the vertical direction is two-thirds of that in the horizontal direction. Figure 10 shows that the acceleration response spectrum for 0-5 s, 0-10 s, 0-15 s, and 0-20 s follows a monotonic linear relationship which establishes that between the interval times of the ETAH and PGAs. In this model, the acceleration response spectrum for 0 to 20 s corresponds with a PGA of 1.2 g, while the acceleration response spectrum from 0 to 5 s corresponds to 0.3 g, and other acceleration response spectra follow the same trend with a uniform scale factor. e dynamic response of Baihetan arch dam-reservoir-foundation system under ten sets of ETAHs is analysed by the commercial finite element software ABAQUS. Figure 11 shows the damage profiles of the arch dam under four typical PGAs for two groups of ETAHs. When the PGA is 0.3 g, 0.6 g, 0.9 g, and 1.2 g, the damage distribution characteristics on the downstream face, upstream face, and crown cantilever are also shown in Figure 11, respectively. Considering that the compressive strength is significantly higher than the  Advances in Civil Engineering tensile strength, only the tensile damage to concrete is considered [12]. e damaget in the graphics represents tensile damage of the arch dam, which is in correspondence with d t in Figure 4. e results show that the damage to the arch dam increases gradually with increasing seismic intensity measure.

Numerical Analysis of the Arch Dam
e damage first appears at the bottom of dam on its upstream face and, subsequently, develops in the middle of the dam on its downstream face with increasing PGA. e main region showing damage lies on the arch base surface. Moreover, the upper region of the downstream face, with its greater curvature, is more damage prone. In particular, as seen in Figure 11, even if the seismic intensity measure is the same, the damage distributions differ slightly because of the randomness of the ETAH; however, the overall damage distributions are similar.

Damage Indices Analysis.
To quantify the degree of damage to the arch dam, two damage indices are calculated using equations (4) and (5), respectively. Figure 12 shows the maximum, minimum, mean, and mean ± standard deviation (SD) of the two damage indices. Figure 12 shows that the two damage indices have a cumulative effect and follow similar exponential evolutionary trends, suggesting that these damage indices share a certain equivalence. e data show that damage indicator D g,1 is more sensitive to seismic intensity measure than D g,2 . us, D g,1 is the better indicator, and the SDs are approximately 50% and 25% of the means of the two damage indices at 10 s and 20 s, respectively. is suggests that the discreteness of these damage indices will rapidly decline at a large seismic intensity measure.   Advances in Civil Engineering 9

Deformation
to t. e deformation indices can be defined by the following formula: where Ω is maximum absolute operator and f(t) denotes the time history of the considered response. Figure 13 shows the distribution of maximum stream deformation along the crest and the crown cantilever caused by 10 sets of ETAHs. e maximum, minimum, mean, and mean ± SD are also plotted. As shown in Figure 13, the stream deformation response of the dam increases from the base to the crest as well as from the abutments to the centre. e peak values generally occur at the crest centre; therefore, the stream deformation at the crest centre can be used as an important fragility deformation indicator to different PGAs. Figure 13 also shows that the difference of stream deformation (the maximum, minimum, mean, and mean ± SD) in the downstream (+) and upstream (−) direction increases from a PGA of 0.3 g to 0.6 g but less so from 0.6 g to 1.2 g.
is suggests that the stream deformation is sensitive to initial damage cracking propagation in the dam.

Deformation Indices Analysis.
e stream deformation at the crest centre is used as a fragility deformation indicator. Figure 14 shows the stream deformation at the crest centre (the maximum, minimum, mean, and mean ± SD in the downstream and upstream directions at different PGAs). From 0 s to 20 s, the stream deformation at the crest centre in the downstream (+) and upstream (−) directions increases in a quasilinear manner and the SD is approximately 60% of the means of the stream deformation at the crest centre from 10 s to 20 s. is suggests that the discreteness of these stream deformation indices is relatively stable under a large seismic intensity measure. Figure 15 shows the statistical distribution of the maximum of joint opening, which includes the mean, maximum, minimum, and mean ± SD values for PGAs of 0.3 g, 0.6 g, 0.9 g, and 1.      respectively. At a PGA of 0.6 g, as the PGA increased, the third joint near the left abutment and 18 th joint approximately at the crest centre both undergo large opening displacement. e mean and maximum openings increase nearly five-fold, reaching approximately 4 mm and 20 mm, respectively. At a PGA of 0.9 g, the third joint near the left abutment, 28 th joint near the right abutment, and 18 th joint at approximately the crest centre all undergo large opening displacements. e mean and maximum openings increase nearly four-fold, reaching approximately 15 mm and 62 mm, respectively. At a PGA of 1.2 g, the distributions are similar to those at 0.9 g. e mean and maximum openings still increase slightly, reaching approximately 24 mm and 62 mm, respectively, suggesting that when the PGA reaches a certain value, such as 0.9 g in this case, the joint opening will slowly increase and no obvious change in the distribution can be seen. In general, the peak joint openings occur at the crest centre and near the abutment. Figure 16 depicts two joint opening indices the maximum of the joint opening (MJO) and the sum of the joint opening (SJO) that include the mean, maximum, minimum, and mean ± SD values. Figure 16 shows that the two joint opening indices follow similar exponential evolutionary trends. is suggests that these joint opening indices also share a certain equivalence. e scatter is more significant than that of the deformation indices and damage indices, as shown in Figure 16. Furthermore, the scatter of the MJO is more significant than that of the SJO, suggesting that the SJO is a better joint opening indicator: the SDs are approximately 90% and 30% of the means for the MJOs at 10 s and 20 s but approximately 75% and 30% of the means for the SJO at 10 s and 20 s. is suggests that the discreteness of these joint opening indices will rapidly decline at large PGA.

Dispersion Analysis.
To explore the relationship between the strong motion intensity parameters and the seismic resistance of this arch dam, ten sets of dynamic response curves corresponding to each PGA obtained by ETA numerical analysis are determined. e damage volume ratios (D g,1 and D g,2 ), joint opening (SJO and MJO), and maximum deformation (MRD and MVD) are also selected in this study. MRD and MVD responses represent the maximum radial deformation and vertical deformation of the dam crest, respectively. Subsequently, the response curves are used in a spline interpolation method at each intensity measure. On the basis of ten groups of dynamic responses analysed by ETA, the median and mean values are calculated to represent the central tendency of the response indicators.   Advances in Civil Engineering e horizontal axis shows the damage indices in terms of the dynamic response, and the vertical axis denotes the PGA of the stream component which represents the intensity of ground motion. Meanwhile, the cumulative frequency curves (16% and 84% fractiles) are also depicted in Figure 17.
It is inferred from Figure 17 that the mean value of the aforementioned indices has an increasing trend and most curves between the 16% and 84% fractiles show excellent consistency between each index. ree limit states of the arch dam can be determined by dynamic response indices. e tendency of median curve of various dynamic responses including damage volume ratio, joint opening, and deformation is analysed to determine three limit states. According to the mean joint opening value, there is little cracking up to 0.3 g, indicating that the joint is not reflected in the nonlinear analysis nor is any such nonlinear resistance mobilised to maintain the stability of the dam. Meanwhile, slight damage occurs at the heel of the arch dam ( Figure 11). Contraction joints are known to offer little resistance and will open and close repeatedly under strong ground motion. e joint opening can be used to describe the failure mode of an arch dam; therefore, when the intensity reaches 0.3 g, there is no obvious opening of contraction joints, and the limit state can be defined as "minor" grade. en, the slope and growth rate observed in the figure also show a significant effect of various indices, and the yield point also can be determined therefrom. e damage index, joints opening, and deformation increase significantly beyond approximately 0.7 g, especially D g,1 and SJO. e downstream surface of the dam suffers a small range of damage. us, the mean value of indicators corresponding to the intensity of 0.6 g is regarded as "moderate" grade. An intensity of 1.2 g causes the  most damage and is regarded as "serious" grade. To perform an accurate assessment of this arch dam, three limit states are proposed based on dynamic results and damage distribution. e coefficient of variation of indices can be found from the variance and mean, as follows: where CV is the coefficient of variation of indices, X is response index, D (·) is the variance of the indices, and E (·) is the mean value thereof.
e dispersion curves of ETA responses ( Figure 17) at different intensities are plotted in Figure 18. e degree of dispersion at each intensity measure can thus be deduced: a low dispersion can reflect the rationality of the optimal indices. e dispersion of the deformation index data is less than that of the damage index and less than that related to joint opening. In the state of minor, moderate, and serious grade damage, response results exhibit relatively little dispersion

Correlation Analysis between Different Indicators.
To reveal the correlation between different indices, the correlation coefficient is calculated using the following formula: where X and Y are response indices, ρ XY is the correlation coefficient, and Cov (X, Y) is the covariance of X and Y. Figure 19 shows the variation in correlation coefficient with ground motion intensity. e correlation coefficient of the identical category indices is strong and stable at each intensity level, such as D g,1 and D g,2 , SJO and MJO, and MRD and MVD indices, suggesting similar dynamic characteristics among these indices. e correlation between D g,1 and joint opening is greater than that with D g,2 , indicating that the sensitivity of D g,1 is generally good for contraction joints, and this can be used to evaluate the overall damage to the structure. For the comparison of deformation, vertical deformation has a better correlation with other indices than radial deformation. Compared with other indices, there is little difference between the sum of joint opening and the maximum of the joint opening, illustrating that both the overall evaluation index and local value can reflect the extension of transverse joints within the dam body. e damage index D g,1 correlates with joint opening indices, while D g,2 has a relatively good correlation with vertical deformation contrary to contraction joints indices at each level of intensity measure. In conclusion, damage and joint opening indices can be selected as effective indicators by comparing the coefficients thereof.
In practice, the measured response of the arch dam is mainly seen in terms of the deformation index (such as displacement and contraction joint opening), a damage index that can help the engineer to evaluate the overall degree of damage to the structure, but the damage index is difficult to detect in practice, so how to predict the damage index from the deformation index is a key problem that needs to be solved.
To demonstrate the previous research and establish the relationship between response indicators, the linear statistical curves and factors are modelled. e correlation coefficient between various indices is summarised in Table 2. As depicted, D g,1 has a remarkable effect on the SJO index in contraction joint opening, followed by that of the MVD index. e relationships between damage index and SJO index and MVD index are shown in Figure 20. e composite index using SJO and MVD indices is used to predict the damage index through multivariate fitting. In particular, the trend line and R 2 value in Figure 20 show that damage index D g,1 has a significant correlation with joint opening under a sequence of increasing intensity levels. In addition, R 2 for the composite index is 0.926; this is greater than that of each individual index (0.906 and 0.624) including the SJO and MVD indices, indicating that the comprehensive model containing SJO and MVD indices has a relatively good correlation with D g,1 . e statistical formulae linking SJO, MVD, and the composite index are given in equations (18)- (20), respectively. In this way, the overall damage characteristics of the structure can be deduced by specific physical quantities such as the joint opening and maximum vertical deformation:

18
Advances in Civil Engineering

Fragility Analysis
ree limit states are defined by Table 3, and the descriptions of three limit states are given based on response results of nonlinear dynamic analysis. According to the description, the corresponding phenomena and the value of dynamic response threshold are determined for the usage of fragility analysis. Regarding the mean value as an index and dividing the cases into three failure grades according to the intensity of ground motions, limit state 1-minor (<0.3 g), limit state 2moderate (0.3-0.6 g), and limit state 3-serious (0.6-1.2 g) are defined. e mean damage indices under three limit states are listed in Table 4. In conjunction with the samples of fragility indices and the selection of limit states value, the failure probability of indices is obtained at different intensities by assuming these fragility indices obey a lognormal normal distribution (equation (6)). On the basis of the investigation above, two indicators with a strong correlation are compared. Figure 21 describes the evolution of fragility indices such as damage indices and joint opening indices including D g,1 and SJO. As shown, the probability curves of these fragility indices exhibit remarkable consistency because of the strong correlation between the two indicators.
In conclusion, the dynamic results show that the dam firstly suffers initial damage to its base surface near the I Minor e contraction joint may be opened, but with a small degree of opening that will have no additional impact on the performance of dam II Moderate e nonlinear behaviour of arch dam is shown up in the form of contraction opening and cracking; nonpenetrating cracks appear on the arch dam, and the length of cracks is small; damage may occur downstream face in the centre of dam, which has a certain influence on the seismic performance of the dam body, but it will reach the normal state of use in short time after slight repair III Serious Penetrating cracks appear in the abutment and centre of the arch dam, and the length of cracks is large; the arch dam is collapsed, which has a serious impact on the seismic performance of dam body  upstream face under a history with an intensity of 0.3 g and then suffers damage in the middle of its downstream surface at 0.6 g. us, the indices D g,2 and deformation are very sensitive to the initial damage on the base surface near the upstream face of the dam. e fragility indices D g,1 , MJO, and SJO are very sensitive to the damage in the middle of the downstream surface, and the fragility indices D g,2 and deformation have some intrinsic relationship as do the fragility indices D g,1 , MJO, and SJO. In the practice of dam monitoring, deformation monitoring is the main method recommended.
us, these fragility indices SJO and deformation can reveal the failure mode of such an arch dam.

Conclusions
e analysis of seismic dynamic response and performance for the Baihetan arch dam-reservoir-foundation system is made by determining failure mode, proposing fragility indices, and analysing probability distributions when the dam is subjected to ETAHs for increasing seismic intensity measure. Several conclusions and recommendations are proposed: (1) Cracking firstly occurs on the base surface near the upstream face and then in the middle of the downstream surface with increasing seismic intensity measure. e damage indices D g,1 and D g,2 follow similar exponential evolutionary trends with increasing seismic intensity measure. e damage indicator D g,1 is more sensitive to seismic intensity than D g,2 . e discreteness of these damage indices will rapidly decline at large seismic intensity measure.
(2) e peak joint openings firstly occur at the crest centre and then near the abutment with increasing seismic intensity measure. e joint opening indices MJO and SJO follow similar exponential evolutionary trends with increasing seismic intensity measure. e scatter of the MJO is more significant than that of the SJO, suggesting that SJO is a better joint opening indicator. (3) e fragility indices D g,2 and deformation are very sensitive to the initial cracking on the base surface near the upstream face of the dam. e fragility indices D g,1 , MJO, and SJO are very sensitive to the cracking in the middle of the downstream surface. In the practice of dam monitoring, joint opening and deformation monitoring are the main indices recommended because these fragility indices can reveal the types of damage to such an arch dam. (4) e coefficients of variation and correlation are used to analyse each index and show that the damage index D g,1 has a good correlation with the SJO index. ree limit states are proposed according to the damage distribution and yield points seen in ETA response curves. In addition, the linear equivalent relationship between the two indices is established, and the derived fragility curves also confirmed the consistency of the two indices, which is of great help in engineering practice. (5) On the basis of a multivariate regression analysis, the relationship between damage index Dg,1 and joint opening SJO index, deformation MVD index, is determined, and the coefficient R 2 is greater than that of any single index, which provides a foundation for structural damage prediction and seismic performance evaluation; however, the effects of other factors such as earthquake wave duration on the seismic fragility of an arch dam are still not considered. More research is thus required in the future.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.