Seismic Response of Resilient Bridges with SMA-Based Rocking ECC-Reinforced Piers

Post-earthquake investigation shows that numerous reinforced concrete (RC) bridges were demolished because of large residual displacements. Improving the self-centering capability and hence resilience of these bridges located in earthquake-prone regions is essential. In this regard, a resilient bridge system incorporating engineered cementitious composites (ECC) reinforced piers and shape memory alloy (SMA) energy dissipation components, i.e., SMA washers, is proposed to enhance its resilience when subjected to strong earthquakes. This study commences with a detailed introduction of the resilient SMA-washer-based rocking bridge system with ECC-reinforced piers. Subsequently, a constitutive model of the ECC material is implemented into OpenSees and the constitutive model is validated by test data. The working principle and constitutive model of the SMA washers are also introduced. A series of dynamic analysis on the conventional and resilient rocking bridge systems with ECC-reinforced piers under a suite of ground motions at E1 and E2 earthquake levels are conducted. The analysis results indicate that the resilient rocking bridge system with ECC-reinforced piers has superior resilience and damage control capacities over the conventional one.


Introduction
Post-earthquake field investigations on the damaged bridges revealed that many reinforced concrete (RC) bridges, although designed conforming to the ductility design philosophy, commonly experienced overly large residual displacement which is difficult to recover. For example, more than 100 RC piers were demolished because they suffered from large permanent drift ratio (i.e., 1.5%) after the 1995 Kobe earthquake [1]. Lessons drawn from these events enlighten us that only satisfying the seismic ductility demand is not enough for engineering structures because their residual deformation after earthquake still significantly jeopardizes their normal functionality [2][3][4]. To guarantee service operation of the structures after earthquakes, resilient capacity is being paid more attention in the seismic codes of many countries (e.g., US, Japan, and New Zealand) [5]. Rocking component, as a resilient structural member, has been attracting extensive experimental and numerical studies [6][7][8][9][10]. For instance, the conventional post-tensioned (PT) rocking bridges have been studied by shake table tests recently [11][12][13]. These studies revealed that these self-centering bridge systems were capable of sustaining a large drift ratio of up to 10% but only experienced small residual drift ratio (i.e., 0.5%) with non-critical damages [14]. Subsequently, a series of innovative devices had been presented to further improve the selfcentering and energy dissipation capacities of the rocking piers under extreme earthquake events [15][16][17][18]. Although the PT tendons together with various energy dissipaters can provide excellent recoverability and energy dissipation capacity to the rocking pier [19,20], the energy dissipater may be damaged and thus should be replaced after earthquakes, leading to compromised rescue efficiency. Additionally, some damage patterns such as Materials 2021, 14, 6500 2 of 21 relaxation and environmental corrosion of the PT tendons are difficult to fix. In this regard, shape memory alloy (SMA) that is characterized by super-elasticity has been recently considered for various devices (i.e., SMA tendons, bars, and springs) used in resilient bridge structures [21][22][23][24][25][26][27][28][29][30][31] as well as other types of structural systems [32][33][34][35][36][37][38][39]. In particular, a bridge system with SMA-washer based rocking pier was recently proposed to achieve self-centering functionality during earthquakes [40]. The SMA washers provided restoring force for the RC pier, which can eliminate some inherent shortcomings, such as corrosion and relaxation, induced by the PT tendon. However, the reinforcing steel embedded in the plastic hinge of the pier was still vulnerable to yield due to large bending moment during severe earthquakes. Varela and Saiidi [41] integrated SMA bars with elastomeric rubber bearing to replace the conventional plastic hinge of the RC pier. The test results indicated that, except for the bucking of the SMA bars, the RC pier experienced almost no damage even under 2.5 times the design's seismic loading. One drawback of the SMA bars is their sensitivity to manufacture imperfections and vulnerable to brittle fracture when sustaining large strains [42,43]. Apart from the SMA material, an innovative engineered cementitious composites (ECC) material was used to replace the conventional concrete in the plastic hinge of the pier to enhance its strength and ductility in both tension and compression [44]. The metal-like strain strengthening property of the ECC material makes it an attractive alternative to the conventional concrete for earthquake resilient structures.
Inspired by the existing studies on rocking piers and the application of emerging materials such as ECC and SMA, a resilient self-centering structural bridge system with SMA-washer-based ECC-reinforced rocking piers is proposed. The SMA washers allows the pier to rock in a large rotation angle under reliable control. The ECC material in the pier can significantly enhance the resilience of the rocking piers and remarkably reduces structural cracks during earthquakes. In this paper, the configuration and self-centering rocking mechanism of the resilient bridge system with ECC reinforced pier is firstly proposed. Subsequently, the constitutive model of the ECC material is introduced and validated by an experimental study on a 1/5 scaled ECC-reinforced pier. An FE model of the prototype rocking bridge with the SMA-washer-based ECC-reinforced rocking pier is established in OpenSees [45] to study the superiority of the proposed solution in alleviating seismic damage over the conventional bridges.

Configuration of the Rocking Bridge System
The configuration of the resilient bridge system with the ECC-reinforced rocking pier is presented in Figure 1. The main body of the pier and two separated pile caps (i.e., the upper pile cap and the lower pile cap) are the most critical components in the resilient bridge system. Each SMA washer set consists of a group of SMA washer springs and associated nuts, anchor rebar and shim plates. The anchor rebar was casted into the concrete of the lower pile cap and passed through the upper one via a plastic tube which was embedded in the upper pile cap. The SMA washers which were assembled in appropriate styles (e.g., in parallel, in series or in combination) passed through the anchor rebar and were tightly precompressed by the nuts. This SMA-washer based rocking solution is capable of providing stable recoverability and moderate energy dissipation capacity for the bridge system when subjected to severe earthquake events. The fixed bearings were installed on the top of the bent cap but sliding bearings are installed on the top of two abutments for the girder deformation induced by temperature variance. A comprehensive experimental study on reduced scale SMA-washer based rocking piers (with normal concrete) has been recently conducted, where the details are reported elsewhere [40]. When ECC is used for the pier, extra benefits are enabled such as enhanced tension and compression strength and ductility to sustain flexural and shear deformation during seismic excitations. In this way, the residual deformation of the resilient bridge system tends to decrease and the damage to the pier will be alleviated. For the prototype bridge, six sets of the SMA washer devices are equipped in the pile. The configuration of a typical SMA washer set is shown in Figure 1b.

Rocking Mechanism and Seismic Design Objectives of Resilient Bridge System
There are three seismic design objectives under different earthquake intensities (i.e., small earthquake, moderate earthquake and large earthquake). The first objective is that the interface between two pile caps is close under small earthquakes, which ensures that the rocking bridge system exhibits similar function to the conventional bridge with fixed base pier. The second one is that the maximum drift ratio of the girder is not more than 1.0% under moderate earthquakes (i.e., E1 level). Such a small drift angle would induce limited damage to the structure. The last one is that the maximum drift ratio of the girder is not more than 2.0% under large earthquakes (i.e., E2 level). When the bridge subjects to severe earthquakes, the resulting inertial force will cause the bridge to rock around the two base corners of the upper pile cap and meanwhile the compressed SMA washers together with the gravity of the bridge will provide restoring force for the bridge to return to its original state. The maximum rocking angle of the pier can be controlled by designing appropriate assembles (i.e., series, parallel or both of them) of the SMA washer set. The natural period of the resilient bridge is much larger than that of the conventional bridge, which makes it far away from the dominant periods of the earthquake ground motions. Thus, the resilient bridge system is expected to significantly reduce seismic damage.
The rocking mechanism of the pier and the seismic design objectives of the bridge can be interpreted by Figure 2. The total height of the pier and the pier cap is H. To satisfy the first seismic design objective, an appropriate prestressed force should be imposed on the SMA washer sets, as shown in Figure 2a. The uniform distribution of the resulting reacting force at the bottom of the upper pile cap is also shown in Figure 2a. The upper pile cap and the bottom pile cap always touch tightly under dead load of the superstructure (i.e., G) and servicing loads such as vehicle and temperature. When the horizontal force F h continuously increases, the reacting force at one base side gradually reduces to zero, as shown in Figure 2b. If F h further increases, the pier will uplift with a drift ratio of α = ∆/H, where ∆ is the horizontal displacement of the pier, as shown in Figure 2c. The horizontal displacement is commonly composed of two parts, of which one part is the deformation of the pier and another part is the contribution of rocking. If the horizontal deformation of the RC pier is too large, the bottom of the pier may yield. Therefore, the ECC material is proposed in this study as an alternative to the conventional concrete material for the pier. The ultimate drift ratio α u can be calculated by ∆ u /H, as shown in Figure 2d. When the designated ultimate drift ratio (i.e., 2.0%) reaches, the SMA washer sets at one side of the pier is fully flattened. If the drift ratio of the bridge exceeds 2.0% during earthquakes, the pier will yield. The proposed resilient rocking bridge system is capable of self-locking, which is a unique property over the conventional one.

Constitutive Models of ECC Material
In order to study the seismic response of the ECC-reinforced structures, a constitutive model [45,46], which is capable of considering the computational efficiency and accuracy, is essential for establishing finite element (FE) model, e.g., OpenSees model. The envelop curves of the constitutive model is composed of two parts, of which one is related to the tension state (i.e., O-A-F-I) and the other part is associated with the compression state (i.e., O-J-P), as shown in Figure 3. The envelop curve in tension is multilinear curves which are expressed as Equation (1), and the envelop curve in compression is also the multilinear curves as written in Equation (2). Segment A-F in the tension region is a strain hardening stage, which indicates that the ECC material is more resilient than the normal concrete. The three segments on the envelope curve within the tensile region (i.e., segments O-A-F-I) are given in where E is Young's modulus. The first microcrack will occur if the strain ε is larger than the cracking strain ε t0 at point A. The loading path further follows the segment A-F until the strain reaches the peak tensile strain ε tp at point F, where the stress simultaneously reaches the peak tensile stress σ tp . However, a soft stiffness (i.e., segment F-I) will occur if the tensile strain exceeds ε tp but is smaller than the ultimate tensile strain ε tu at point I, where the corresponding stress becomes zero. It is worth noting that the loading path will continuously move forward along the positive abscissa axis if the strain is larger than ε tu . The envelope curve consists of multilinear curves in the compressive region (i.e., segments O-J-P) which can be expressed as where σ cp and ε cp are the peak compressive stress and corresponding strain, respectively. ε cu is the ultimate compressive strain. It is noted that loading path will further move forward along the negative abscissa axis if the strain is smaller than ε cu . The loading, the unloading and the reloading rules of the ECC model in the tension region (see Figure 4a) are given by where α t is a constant that is larger than or equal to 1. It can be calibrated using the experimental data. The segment B-C-E is the initial unloading path expressed by the Equation (4). ε tm is the maximum strain reached in history on the envelop curve where the unloading is triggered. ε tl is the strain corresponding to the stress vanishing on the initial unloading path (i.e., the strain associated with the point E or H). The value ε tl = β t ·ε tm , where β t is a constant. Segment C-D is a typically partial reloading path, where the stress at point C is not zero, expressed by Equation (5), which ensures the extension of the segment C-D passing through the historically reached maximum strain point on the envelop curve (i.e., point B). The ε tr and ε tu are strains at points C and D, respectively. The unloading path starting from D is controlled by Equation (4). The subscript on the parameter ε tm denotes that ε tm should be set as ε tm or ε tr when they are used to define the initial unloading path or partially reloading path, respectively. The same specification applies to the parameters ε tl and ε tu . The unloading and reloading paths after point F (i.e., F-G) are a linear curve given by Equation (6). The loading, the unloading and the reloading paths of the ECC model in the compression region (see Figure 4b) are expressed as: The segment K-N, as presented in Figure 4b, is formulated by Equation (8). ε cm is the strain on the envelope curve where unloading is triggered (i.e., the strain at point K) and the ε cl is the strain on the initial unloading path corresponding to zero stress (i.e., the strain at point L). The value ε cl = β c ·ε cm , where β c is a constant. Segment N-M is formulated by Equation (9). ε cr and ε cu are strains of the points M and N, respectively. The unloading path from N is given by Equation (8). The parameter ε cm should be set to ε cm or ε cr when they are employed to define the initial unloading path or partially reloading path.

ECC-Reinforced Column and Numerical Verification
A 1/5 scaled ECC-reinforced column [47,48], as shown in Figure 5a, is used to verify the effectiveness of the ECC constitutive model proposed in Section 3.1. To provide lateral load on the top of the cantilever column, a rigid ECC transverse beam is monolithically casted with the cantilever base. This loading configuration is chosen to promote a flexural deformation mode in the specimen and to investigate the effect of ECC material properties on the expected plastic hinge region in particular. Longitudinal steel reinforcement was bent at a 90º angle at the bottom of the transverse beam and further extended to provide sufficient anchorage. A total of four reinforcing steels with a diameter of 10 mm are arranged at the four corners of the cross section of the pier, resulting in a reinforcement ratio of 3.  The FE model of the ECC-reinforced column is established in OpenSees (Version 2.4.1). A total of 5 displacement-based beam-column elements are used to model the ECC column. The ECC02 material model that has been developed and implemented into OpenSees is employed to capture the response of the column. The analysis results regarding the horizontal force against drift ratio together with the test results are displayed in Figure 6. It indicates that the numerical simulation results agree well with the test results. It confirms that the ECC constitutive model proposed in this study is sufficiently accurate for further dynamic analysis.

Constitutive Models of SMA Washer
The chemical composition of the SMA washer is 55.87% nickel and 44.13% titanium alloys by weight and it is supposed to exhibit super elasticity at room temperature (austenite finish temperature is 4.5 • C). The constitutive model of the SMA washer has been developed in [40], which has a flag-shaped envelop curves, including loading stage (i.e., O-A-B-E) and unloading stage (i.e., E-B-C-E-O), as shown in Figure 7, where F 2 and δ 2 represent the "yield" force and the corresponding deformation, respectively; F 3 and δ 3 are the force and deformation when the SMA washer is fully flattened. Similarly, F 4 , δ 4 , F 1 , and δ 1 represent the characteristic forces and deformations during the reverse plateau. Once the SMA washer reaches its maximum compressive deformation, the axial stiffness increases abruptly (i.e., BE). The loading stag on the envelop curve is composed of three segments, which are expressed as: where E w is the elastic stiffness of the single SMA washer; k 1 and k 2 are constants which can be determined from test. The unloading stage on the envelop curve consists of four segments, which are given as follows: The uploading and reloading stiffness away from the envelop curve is E w .
The key parameters for a single SMA washer are shown in Figure 7.
The diagram of a typical SMA washer is presented in Figure 8, where the outer and inner diameters of the SMA washer are D and d, respectively. The thickness and the height of the SMA washer are h and t, respectively. These parameters of a SMA washer can be appropriately adjusted according to the seismic design objectives of rocking bridges. Once all the parameters are designated, the maximum compressive deformation and the restoring force provided by a SMA washer can be calculated.

Validation of the Numerical Simulation Method for Capturing Rocking Behavior
A 1/4 scale RC rocking pier specimen equipped with SMA washer sets [29] is employed to validate the numerical simulation method used in this study for capturing the overall rocking behavior. The configuration of the rocking pier is shown in Figure 9. The test specimen includes one bent cap, the main body of a pier, two separated pile caps, four SMA washer sets, and a tendon. The diameter of the RC pier is 0.3 m. The strength of the normal concrete was 39.0 MPa at the test day. A total of 12 longitudinal HRB400 rebars with a diameter of 16 mm were uniformly arranged along the perimeter, resulting in a reinforcement ratio of 0.3%. The average yielding strength of the longitudinal rebar was 400.4 MPa. The diameter of the transverse reinforcement was 10 mm and the space between two adjacent stirrups was 75 mm. A total of four SMA washer sets were installed on the top of the upper pile cap and each SMA washer set composed of the loading protocol is shown in Figure 10. Other information regarding this rocking pier can be find elsewhere [29].  To validate the numerical simulation method proposed in this study for capturing the overall rocking behavior, the experimental results of the rocking pier specimen are used to make comparison with the numerical analysis results. The RC pier was divided into seven displacement-based fiber beam-column elements, of which the normal concrete fiber was assigned with Concrete02 material model and the reinforcement fiber was assigned with Steel02 material model. The connection between the pier and the bent cap or the pile cap is modeled by rigid elements. The interface between upper and lower pile caps was simulated by eight zero-length spring elements, of which the tension strength is ignored. Each SMA washer set was modeled by a zero-length spring element, of which the assigned material was the Self-centering material model. The FE model is shown in Figure 10. It is worth noting that the nonlinear material properties of the normal concrete and the reinforcement can be considered by the Concrete02 material model and Steel02 material model, respectively. The geometric nonlinearity of the RC pier can also be accounted for by the co-rotational geometric transformations [39]. The numerical simulation results and the test ones are presented in Figure 11, which shows that the analysis results match well with the test regarding the amplitude of the lateral load and the trend of the hysteretic curve. It confirms that the numerical simulation method proposed in this study can be used to capture the overall rocking behavior of the rocking bridge system. Figure 11. Comparison between numerical simulation and test results.

Seismic Responses of the Resilient Bridges with SMA-Based Rocking Piers
In order to demonstrate the superiorities of the proposed bridge system in resilience enhancement, four different bridge systems are studied. The first bridge system is a conventional bridge with a RC fixed pier. The second one is a conventional bridge with an ECC-reinforced fixed pier. The third one is a resilient bridge with a RC rocking pier with SMA washers. The last one is a resilient bridge with an ECC-reinforced rocking pier with SMA washers.

FE Model of Resilient Bridge Systems
A comprehensive FE model of the resilient continuous RC box girder bridge is established in OpenSees, as demonstrated in Figure 12. The RC girder (i.e., 20 + 20 m) of the bridge is modeled by 35 elastic beam-column elements. The rocking RC or ECC-reinforced pier is 8.0 m high which is modeled by six displacement-based beam-column fiber elements. The Concrete02 material model is used to model the properties of the normal concrete, of which the tension strength is considered. The ECC02 material model is employed to simulate the ECC material. The Steel02 material model is used for modeling the reinforcing steel. The heights of the two separated pile caps are both 2.0 m, which are modeled using rigid elements. The width of each expansion joint is 0.25 m, which is modeled by a zero-length element. Each bent cap and abutment are modeled by two rigid beam elements along the transverse direction of the bridge. Similarly, the transverse beams in the girder are simulated by two rigid beam elements. Two conventional sliding bearings are placed on each abutment and two fixed bearings are installed on the bent cap above the middle pier. The frictional coefficient of the sliding bearing is 0.02. Each bearing is simulated by a zero-length element connecting the girder and the bent caps or the abutments. The soil and structural interaction of each pile foundation has been considered by six zero-length elements [49,50].  Two additional conventional bridge systems and a resilient bridge system with a RC rocking pier are also considered to demonstrate the super resilience of the innovative bridge with an ECC-reinforced rocking pier. The main difference between the conventional bridge and the resilient bridge is that the 10 m high fixed pier in the conventional bridge is replaced by an 8 m high pier and a 2 m rocking upper pile cap in height. The FE models of the other three bridge systems can be easily adjusted from that shown in Figure 12 and are therefore not elaborated.

Earthquake Ground Motions
A suite of ground motions being compatible to the acceleration spectra (i.e., E1 and E2 levels) [51,52] are generated and each suite includes seven earthquake ground motions, as shown in Figure 13a,b, respectively. The exceedance probabilities of E1 and E2 level earthquakes in a recurrence interval of 50 years are 10% and 2.5%, respectively. The damping ratio of the acceleration spectra is 5%.

Comparison of Seismic Responses between the Conventional and Resilient Bridges
To evaluate the damage state of the RC and ECC-reinforced pier after earthquake excitation, sectional characteristic analyses are conducted prior to analysis of bridge system. The equivalent bending moments versus curvatures of the sections are shown in  The curvature ductility and the drift ratio of the pier as well as the shear deformation of the bearing are usually selected as the damage indicators for seismic performance assessment of bridges with rocking piers. Therefore, these indicators are used to make a comparison between the conventional and the resilient bridge systems under the considered earthquakes. The curvature ductility indicates the damage state of the plastic hinge region of the pier during earthquakes. Drift ratio is employed to assess the lateral deformation of the bridge system. The uplift ratio indicates the rocking amplitude of the rocking pier, which is expressed as c/(a + b) shown in Figure 15, where c is the maximum uplift distance of the pier; a is the distance between the outmost edge of the upper pile cap and the left side of the stress block; b is the width of the compressive stress block (i.e., rocking zone). The maximum seismic responses regarding the aforementioned damage indicators of the conventional (fixed base) and the resilient (rocking), bridges with the normal RC piers under E1 and E2 earthquakes are summarized in Tables 1 and 2, respectively.  The results summarized in Table 1 reveal that the average maximum values of the curvature ductility of the pier, bearing deformation, the drift ratio, and the residual drift ratio of the conventional RC bridge are all larger than those of the resilient RC bridge under E1 level earthquakes. The curvature ductility responses of the RC piers of the conventional bridge and the resilient bridge are both less than 1.0, which means that the RC piers in two bridge systems keep linear state under E1 level earthquakes. The drift ratio of the resilient bridge is 0.88%, which satisfied the principle of the seismic design objective. All the seismic responses confirm that the two bridge systems are both safe under E1 level earthquakes. The average maximum uplift ratio is 0.11%. Even though the earthquake intensity (i.e., E1 level) is not large, the unique property, such as rocking, of the resilient bridge is well exhibited.
When the earthquake intensity increases from E1 to E2 level, the average maximum curvature ductility of the conventional RC bridge sharply increases from 0.94 to 2.99. It implies that the RC pier experiences severe damage when the bridge is subjected to E2 level earthquakes although it performs linearly under E1 level earthquakes. A similar trend also can be found from the seismic performance of the resilient RC bridge system even though the increasing amplitude is not so large compared with the conventional bridge. To gain a thorough understanding about the rocking mechanism of the pier, the typical shear force against drift ratio of the bridge, and the bending moment versus curvature ductility of the section at the bottom of the RC pier under a typical ground motion (i.e., Earthquake No. 1 at E2 level) are shown in Figures 16 and 17, respectively. From these two figures, it can be concluded that the drift ratio of the conventional RC bridge (i.e., 1.71) is slightly larger than that of the resilient RC bridge (i.e., 1.63), but the maximum curvature ductility of the RC pier and the residual drift ratio of the conventional bridge (i.e., 3.24 and 0.039%) are much larger than those of the resilient bridge (i.e., 1.47 and 0.009%). It confirms that the resilient rocking bridge system can significantly alleviate the seismic damage to the RC pier because the dominate period of the rocking bridge is elongated. The seismic response of a SMA washer set is displayed in Figure 18, where the compressive deformation is 0.021 m that is less the uplift threshold of the rocking pier (i.e., 0.03 m). The hysteretic loop response as shown in Figure 18 demonstrates that the SMA washer sets dissipates a lot of energy during earthquakes. 0.009%). It confirms that the resilient rocking bridge system can significantly alleviate the seismic damage to the RC pier because the dominate period of the rocking bridge is elongated. The seismic response of a SMA washer set is displayed in Figure 18, where the compressive deformation is 0.021 m that is less the uplift threshold of the rocking pier (i.e., 0.03 m). The hysteretic loop response as shown in Figure 18 demonstrates that the SMA washer sets dissipates a lot of energy during earthquakes.      The maximum and average maximum seismic responses regarding the aforementioned damage indicators of the conventional and the resilient bridges with the ECC-reinforced piers under E1 and E2 level earthquakes are summarized in Tables 3  and 4, respectively. The results presented in Table 3 show that all the average maximum responses of the resilient bridge except for the curvature ductility of the ECC-reinforced pier are smaller than those of the conventional bridge under E1 earthquakes. Even though the maximum average curvature ductility of the resilient bridge is larger than that of the conventional one, it is still less than 1.0. In other words, the ECC-reinforced pier in the resilient bridge always experiences elastic state under E1 earthquakes.
When the earthquake intensity increased from E1 to E2 level, all the average maximum responses sharply enhanced except the curvature ductility of the ECC-reinforced pier of the resilient bridge increases slightly (i.e., from 0.74 to 0.98). The average maximum responses of the bearing deformation of the conventional and the resilient ECC-reinforced bridges are similar under E2 level earthquakes, but the average maximum curvature ductility of the ECC-reinforced pier and the residual drift ratio of the resilient bridge (i.e., 0.98 and 0.011%) are much smaller than those of the conventional bridge (i.e., 1.80 and 0.073%). The ECC-reinforced pier in the resilient bridge system can always keep in linear performance under E2 earthquakes but it may yield in the conventional bridge system under the same level earthquakes. The average maximum drift ratio of the conventional bridge is larger in comparison with the resilient bridge because the ECC-reinforced pier in the conventional bridge yields but the ECC-reinforced pier in resilient bridge keeps a linear behavior. To take a close look at the seismic response, the typical time history responses of the lateral seismic force versus drift ratio of the ECC-reinforced conventional and the ECC-reinforced resilient bridges under a typical earthquake (i.e., Earthquake No.1 at E2 level) are illustrated in Figure 19. It can be confirmed that although the two bridges experienced similar lateral displacement, the rocking effect contributes to the most part of the total lateral displacement but the fixed pier mainly depended on the flexural deformation of the pier itself. The seismic responses of the section at the bottom of the plastic hinge region in two bridge systems are presented in Figure 20, which implies that the fixed base ECC-reinforced pier suffered severe damage but the rocking pier stayed elastic. The seismic response of one SMA washer set is displayed in Figure 21, where the self-locking effect is triggered when the compressive deformation reached 0.03 m. Once the maximum drift ratio of the resilient bridge reaches 2.0%, the self-locking effect will act and the ECC-reinforced pier will yield simultaneously. In other words, if the drift ratio of the resilient bridge exceeds 2.0%, the incremental lateral displacement will completely depend on the yielding deformation of the ECC-reinforced pier.
ECC-reinforced resilient bridges under a typical earthquake (i.e., Earthquake No.1 at E2 level) are illustrated in Figure 19. It can be confirmed that although the two bridges experienced similar lateral displacement, the rocking effect contributes to the most part of the total lateral displacement but the fixed pier mainly depended on the flexural deformation of the pier itself. The seismic responses of the section at the bottom of the plastic hinge region in two bridge systems are presented in Figure 20, which implies that the fixed base ECC-reinforced pier suffered severe damage but the rocking pier stayed elastic. The seismic response of one SMA washer set is displayed in Figure 21, where the self-locking effect is triggered when the compressive deformation reached 0.03 m. Once the maximum drift ratio of the resilient bridge reaches 2.0%, the self-locking effect will act and the ECCreinforced pier will yield simultaneously. In other words, if the drift ratio of the resilient bridge exceeds 2.0%, the incremental lateral displacement will completely depend on the yielding deformation of the ECC-reinforced pier.

Comparison of Seismic Responses between the RC and ECC-Reinforced Resilient Bridges
The aforementioned comparisons demonstrates that the ECC-reinforced bridge systems are more resilient than the conventional bridge systems. The following part will make a comparison between the two resilient bridge systems with the RC rocking pier and the ECC-reinforced rocking pier. The average maximum drift ratio of the resilient bridge with RC rocking bridge under E2 level earthquakes is 1.49 but the value of the ECC-reinforced resilient rocking bridge is 1.70. The reason is that the yielding strength of the RC pier is smaller than the ECC-reinforced pier so that it cannot sustain large rocking amplitude. It can reconfirm by the response of the average maximum curvature ductility of two bridges. For instance, the average maximum curvature ductility of the resilient bridge with RC rocking bridge under E2 level earthquake is 1.50, whereas the value of the ECC-reinforced resilient rocking bridge is 0.98. A case is selected for further investigation: the drift ratio versus lateral seismic force responses of two resilient bridges subjected to a typical earthquake (i.e., Earthquake No. 1 at E2 level) is shown in Figure 22. From Figure 23, it can be recognized that the maximum drift ratio of the resilient bridge with the ECC-reinforced pier is 2.07, but the corresponding value is only 1.63. The curvature ductility versus bending moment responses of two resilient bridges under a typical earthquake (i.e., Earthquake No. 1 at E2 level) is shown in Figure 23. The maximum curvature ductility of the RC pier is 1.47 but the counterpart of the ECC-reinforced pier is 1.07. The earthquake-induced damage of the ECC-reinforced pier is so small (i.e., 1.07), which needs no repair but the RC pier should be retrofitted after earthquake.

Conclusions
This study presented a resilient self-centering bridge system with the ECC-reinforced rocking pier. An experimental study on the ECC-reinforced column was employed to demonstrate the accuracy of the proposed constitutive model of the ECC material implemented in OpenSees. The advantages of the resilient ECC-reinforced rocking bridge system, such as low damage and negligible residual deformation, were verified by conducting a series of nonlinear seismic analyses on both resilient and conventional bridges with RC or ECC-reinforced piers. Several comments and conclusions are summarized as follows: 1.
All the conventional and the resilient rocking bridge systems with the RC pier or the ECC-reinforced pier can satisfy the seismic design objectives under E1 level earthquakes; 2.
The proposed resilient rocking bridge system with the RC pier or the ECC-reinforced pier shows superior seismic performance over the conventional bridge in terms of the response such as the curvature ductility of the pier, bearing deformation, drift ratio, and residual drift ratio of the bridge under E2 level earthquakes; 3.
The resilient rocking bridge with the ECC-reinforced pier can achieve superior damage retrofitting capacity than the resilient rocking bridge with RC pier. The average maximum curvature ductility of the ECC-reinforced pier was only 0.98 but the counterpart of the RC pier was 1.5 which indicates the RC experienced yielding state during E2 level earthquakes; 4.
Due to the protection mechanism (i.e., self-locking effect) of the SMA washer spring device against overload and the super resilient property of the ECC material, the resilient rocking bridge system shows the excellent damage control capacity; 5.
The SMA washer spring device cannot only provide restoring force for the resilient bridge system but also can dissipate moderate earthquake energy input.