Multi-Objective Optimization of Three Di ﬀ erent SMA-LRBs for Seismic Protection of a Benchmark Highway Bridge against Real and Synthetic Ground Motions

: Many researchers have taken advantage of adding shape memory alloy (SMA) wires to base isolators to control displacements and residual deformations. In the literature, di ﬀ erent arrangements of SMA wires wrapped around the rubber bearings can be found, as examples, straight, cross and double-cross arrangements. SMA wires with various conﬁgurations and radii lead to the di ﬀ erent characteristics of the isolator system and thus various shear hysteresis. Therefore, the aim of this study is to evaluate the performance of these three SMA wire’s conﬁgurations in the seismic retroﬁtting of a benchmark highway bridge by implementing them in the bridge’s existing lead rubber bearings (LRB). This system is referred to as SMA-LRB isolator. Firstly, because of the crucial inﬂuence of the wire’s radius, this parameter is determined using a multi-objective optimization algorithm (non-dominated sorting genetic algorithm (NSGA)-II). This algorithm simultaneously minimizes the deck acceleration and mid-span displacement. Secondly, the optimized SMA-LRBs are implemented in the highway bridge and nonlinear dynamic analysis is conducted. For the nonlinear response history analysis, two strong ground motion records are selected from the PEER database, by studying the site’s conditions. In addition, ten synthetic ground acceleration time histories are generated. The result illustrates that the double-cross SMA-LRB reduces the maximum and residual displacements more than two other devices; however, it causes the largest base shear force and deck acceleration. Besides, the cross-conﬁguration results in the least displacement reduction and has the least shear force and acceleration. To ﬁnd SMA-LRB with the best overall performance, a multi-objective decision-making method is utilized and the straight SMA-LRB is recognized as the most e ﬀ ective isolator.


Introduction
Bridges play a critical role in transportation networks, and any disruption in their function causes huge direct and indirect economic and life losses. Thus, it is of paramount importance to protect bridges against earthquakes or other severe loadings. In recent major earthquakes such as Northridge 1994 [1], Kobe 1995 [2], Chi-Chi 1999 [3], Wenchuan 2008 [4], Chile 2010 [5], and Great East Japan 2011 [6] a variety of damage types have occurred in highway bridges and proved the limitation of the conventional design methods. The conventional seismic design approaches are based on the dissipation of the input earthquake energy through partial inelastic deformations in the designed areas of the structural elements, and at the same time, maintaining the overall stability and integrity of A three-dimensional finite element model of this bridge was developed using ABAQUS software, in two phases. In phase I, the semi-isolated model, the deck of the bridge is fixed to the center outrigger, and eight seismic isolators are applied to separate the bridge from the end abutments [27]. In phase II, the completely isolated model, two lead rubber bearings are installed between the deck and the center outrigger in addition to the eight aforesaid isolators [28]. Many researchers have implemented their control systems in phases I and II of the benchmark bridge. For example, Tan and Agrawal [29] presented sample designs for three control systems and applied them in phase I. Madhekar et al. [30][31][32] investigated the performance of different passive and semi-active devices on the phase I bridge. Saha et al. [33][34][35] also applied passive control systems on this bridge. Regarding sample controllers in phase II, Zhang et al. [36] proposed a damper with SMA for controlling the displacements and studied the impact of temperature changes on the SMA behavior. Casciati et al. [37] suggested a new SMA device as a passive damper to reduce the maximum displacement of the base isolators. Li et al. [38] investigated the performance of a novel negative stiffness device (NSD) and a damper system in seismic control of the benchmark bridge. In this research, the completely isolated benchmark bridge in phase II is applied, which is modelled by Nagarajaiah et al. [28]. A brief review of the bridge modelling is provided in the following section.

Structural Model
To model the completely isolated highway bridge, Nagarajaiah et al. [28] took advantage of ABAQUS finite element software. The model has a total of 108 nodes, 4 rigid links, 70 beam elements, 24 springs, 27 dashpots, and 10 user-defined bearing elements. The superstructure, including the deck and the beam, is assumed to be elastic; also, the abutments and deck-ends are assumed to be rigid and have a skew angle of 33°. Nonlinearities are included in members to capture inelastic moment-curvature behavior of columns and the shear-displacement relationship of bearings. The bearings are idealized using the bi-directional bilinear plasticity model, as described by Makris and Zhang [39], and defined as a subroutine in ABAQUS. The bearings are modelled as a shear element, which means they express infinite vertical stiffness and zero torsional rigidity and bending stiffness. In addition, the mass and damping of these bearings are eliminated. The soil-structure interaction is considered at the end of abutments and embankments, and the pile foundations at both abutments are assumed to behave like an equivalent linear visco-elastic element. More details on this model, in conjunction with their dynamic properties, are described in the definition paper [27].
After developing the three-dimensional model of the bridge with 430 degrees of freedom, all element mass matrices and initial elastic element stiffness matrices obtained are summed at nodal masses to assemble global stiffness and mass matrices within MATLAB environment. The global damping matrix C is a combination of the distributed 5% inherent Rayleigh damping in the first two modes and soil radiation damping. The fundamental period of this phase II benchmark isolated bridge is 5.1 s. The dynamic equation of this model is expressed in the form of Equation (1). M (t) + C (t) + K(t)U(t) = −Mƞ (g) + bF(t) (1) A three-dimensional finite element model of this bridge was developed using ABAQUS software, in two phases. In phase I, the semi-isolated model, the deck of the bridge is fixed to the center outrigger, and eight seismic isolators are applied to separate the bridge from the end abutments [27]. In phase II, the completely isolated model, two lead rubber bearings are installed between the deck and the center outrigger in addition to the eight aforesaid isolators [28]. Many researchers have implemented their control systems in phases I and II of the benchmark bridge. For example, Tan and Agrawal [29] presented sample designs for three control systems and applied them in phase I. Madhekar et al. [30][31][32] investigated the performance of different passive and semi-active devices on the phase I bridge. Saha et al. [33][34][35] also applied passive control systems on this bridge. Regarding sample controllers in phase II, Zhang et al. [36] proposed a damper with SMA for controlling the displacements and studied the impact of temperature changes on the SMA behavior. Casciati et al. [37] suggested a new SMA device as a passive damper to reduce the maximum displacement of the base isolators. Li et al. [38] investigated the performance of a novel negative stiffness device (NSD) and a damper system in seismic control of the benchmark bridge. In this research, the completely isolated benchmark bridge in phase II is applied, which is modelled by Nagarajaiah et al. [28]. A brief review of the bridge modelling is provided in the following section.

Structural Model
To model the completely isolated highway bridge, Nagarajaiah et al. [28] took advantage of ABAQUS finite element software. The model has a total of 108 nodes, 4 rigid links, 70 beam elements, 24 springs, 27 dashpots, and 10 user-defined bearing elements. The superstructure, including the deck and the beam, is assumed to be elastic; also, the abutments and deck-ends are assumed to be rigid and have a skew angle of 33 • . Nonlinearities are included in members to capture inelastic moment-curvature behavior of columns and the shear-displacement relationship of bearings. The bearings are idealized using the bi-directional bilinear plasticity model, as described by Makris and Zhang [39], and defined as a subroutine in ABAQUS. The bearings are modelled as a shear element, which means they express infinite vertical stiffness and zero torsional rigidity and bending stiffness. In addition, the mass and damping of these bearings are eliminated. The soil-structure interaction is considered at the end of abutments and embankments, and the pile foundations at both abutments are assumed to behave like an equivalent linear visco-elastic element. More details on this model, in conjunction with their dynamic properties, are described in the definition paper [27].
After developing the three-dimensional model of the bridge with 430 degrees of freedom, all element mass matrices and initial elastic element stiffness matrices obtained are summed at nodal masses to assemble global stiffness and mass matrices within MATLAB environment. The global damping matrix C is a combination of the distributed 5% inherent Rayleigh damping in the first two modes and soil radiation damping. The fundamental period of this phase II benchmark isolated bridge is 5.1 s. The dynamic equation of this model is expressed in the form of Equation (1). M ..
In this equation, the incremental displacement, ground acceleration, the incremental controlling force, mass, and stiffness matrices are denoted by U, .. U, F, M, and K, respectively. Note that the stiffness matrix of the structure includes linear and nonlinear sub-matrices. η and b are the loading vectors for the ground acceleration and control forces, respectively. Equation (1) can be solved with the help of the Newmark integration method. MATLAB software has various tools for solving ordinary differential equations; however, they are not able to find the nonlinear structural responses. Thus, Nagarajaiah et al. [28] developed an S-Function program in the Simulink for the bridge model. Figure 2 shows the Simulink of the passive control device. Researchers can modify the contents of the 'control device' to accommodate their control designs.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 5 of 28 In this equation, the incremental displacement, ground acceleration, the incremental controlling force, mass, and stiffness matrices are denoted by U, , F, M, and K, respectively. Note that the stiffness matrix of the structure includes linear and nonlinear sub-matrices. η and b are the loading vectors for the ground acceleration and control forces, respectively. Equation (1) can be solved with the help of the Newmark integration method. MATLAB software has various tools for solving ordinary differential equations; however, they are not able to find the nonlinear structural responses. Thus, Nagarajaiah et al. [28] developed an S-Function program in the Simulink for the bridge model. Figure 2 shows the Simulink of the passive control device. Researchers can modify the contents of the 'control device' to accommodate their control designs.

Earthquake Excitations
In this study, for the nonlinear dynamic analysis of the bridge under earthquakes, the nonlinear response history analysis is carried out by selecting a set of ground motion time histories. For this purpose, researchers often select ground motion time histories from a database of past recorded earthquakes. However, recorded motions are not always available for all types of earthquakes and in all regions. Therefore, researchers often have to choose records from locations other than the project site and modify these records by scaling or spectrum matching techniques to fit their desired conditions. Sometimes, they are forced to scale recorded motions by factors as large as 10 or 20 or modify the frequency contents to achieve desired intensity or frequency characteristics [40]. Such modified motions may not accurately characterize real earthquake ground motions. To tackle this deficiency, researchers have developed techniques for the generation of synthetic ground motions that simulate realistic ground motions for specific design scenarios [40,41].
In the benchmark bridge package, six ground motions from past earthquakes were selected without scaling or spectrum matching [27]. These earthquakes were North Palm Springs (1986), TCU084 component of Chi-Chi earthquake Taiwan (1999), El Centro component of Imperial Valley earthquake (1940), Rinaldi component of Northridge earthquake (1994), Bolu component of Duzce Turkey earthquake (1999), and Nishi-Akashi component of Kobe earthquake (1995). Because of the weaknesses of selecting ground motions from inappropriate locations with unreasonable scaling or spectrum matching, in this study, strong motion records are carefully selected to match site conditions. Besides, 10 synthetic motions are generated to be used along with the real records in the nonlinear response history analysis. The approach for selecting real records of strong ground motion and generating synthetic records are presented in the next section.

Earthquake Excitations
In this study, for the nonlinear dynamic analysis of the bridge under earthquakes, the nonlinear response history analysis is carried out by selecting a set of ground motion time histories. For this purpose, researchers often select ground motion time histories from a database of past recorded earthquakes. However, recorded motions are not always available for all types of earthquakes and in all regions. Therefore, researchers often have to choose records from locations other than the project site and modify these records by scaling or spectrum matching techniques to fit their desired conditions. Sometimes, they are forced to scale recorded motions by factors as large as 10 or 20 or modify the frequency contents to achieve desired intensity or frequency characteristics [40]. Such modified motions may not accurately characterize real earthquake ground motions. To tackle this deficiency, researchers have developed techniques for the generation of synthetic ground motions that simulate realistic ground motions for specific design scenarios [40,41].
In the benchmark bridge package, six ground motions from past earthquakes were selected without scaling or spectrum matching [27]. These earthquakes were North Palm Springs (1986), TCU084 component of Chi-Chi earthquake Taiwan (1999), El Centro component of Imperial Valley earthquake (1940), Rinaldi component of Northridge earthquake (1994), Bolu component of Duzce Turkey earthquake (1999), and Nishi-Akashi component of Kobe earthquake (1995). Because of the weaknesses of selecting ground motions from inappropriate locations with unreasonable scaling or spectrum matching, in this study, strong motion records are carefully selected to match site conditions. Besides, 10 synthetic motions are generated to be used along with the real records in the nonlinear response history analysis. The approach for selecting real records of strong ground motion and generating synthetic records are presented in the next section.

Ground Motion Time Histories
For nonlinear response history analyses and the optimization algorithm in this study, instead of selecting input ground motion time histories from past earthquakes in areas other than the bridge site and modifying them significantly to match the current conditions, real records of strong motions are utilized that resemble the following conditions: 1.
Distance to the causative fault approximately the same as the distance of the bridge to the closest active seismic fault.

2.
Size of the earthquake equal to the probable maximum magnitude at the closest fault.

3.
Soil conditions and deeper geological surroundings of the recording station similar to those beneath the bridge.
The bridge is located in Orange County of Southern California, and the two closest fault zones are the Whittier-Ellsinore, which is 11.6 km to the northeast, and the Newport-Inglewood, which is 20 km to the southwest of the bridge. Both faults are right-lateral strike-slip faults [42,43]. The Whittier-Ellsinore fault has a slip rate of 2.5 to 3.0 mm/year and is estimated to be capable of a 6.0-7.2 magnitude (M w ) earthquake, while the Newport-Inglewood fault has a slip rate of approximately 0.6 mm/year and could generate a 6.0-7.4 magnitude (M w ) earthquake [44]. Therefore, in this study, it is decided to analyze the benchmark bridge when exposed to an M w = 7.2 earthquake with the causative strike-slip fault at around 11 km from the site.
According to Makris and Zhang [45], a geotechnical exploration at the location of the piers and near the end abutments was conducted before construction, down to a depth of about 35 m. The moderately stiff soil was identified, and the shear wave velocity of 200 m/s was estimated [45]. Hence, real strong motion records at the stations with similar local soil conditions were looked for. A further requirement is for the deeper geological conditions of the recording sites to resemble the geological surroundings of the benchmark bridge. The reason is that the geological characteristics up to depths of hundreds of meters or even a few kilometers may strongly affect the severity of surface seismic waves with both the shorter and longer oscillation periods [46]. Several recent seismic micro zonation and strong motion studies have shown that this influence is also regionally dependent [47,48]. Therefore, to properly capture deeper geological conditions, it is decided to use the strong motion time histories recorded in Southern California as close to the benchmark bridge as possible.
Usually, it is very difficult to find strong motion records that will match all defined criteria above, but in this study, fortunately, such records are found. Table 1 illustrates the two records (each with three orthogonal components, i.e., two horizontal and one vertical) available at PEER Ground Motion Database [49], which closely match the defined earthquake and site conditions. Table 1 includes PEER Record Sequence Numbers (RSN), station name, pulse period, significant duration, Peak Ground Acceleration (PGA), earthquake name, magnitude, fault mechanism, distance to the rupture, and shear wave velocity for the first 30 m of the stratigraphic profile. For both records, the causative earthquake was the 2010 "El Mayor-Cucapah Mexico" earthquake with M w = 7.2 and strike-slip fault mechanism. Moreover, both were recorded around 10 km from the fault, in Southern California, with the shear wave velocity parameter, Vs30, relatively close to 200 m/s. Another substantial point is that both records represent near-source ground motions that exhibit distinguishable strong acceleration and velocity pulses. In addition to the two real records, 10 synthetic accelerograms are generated to account for stochastic features of strong ground motion. Figure 3a,b show the recorded acceleration time histories described in Table 1, and the corresponding synthesized accelerograms, for both horizontal ground motion components.   Table 1 for the station "El Centro Array #12" and five corresponding spectrum-compatible and energy-compatible synthetic acceleration time histories; (b) Acceleration time histories of the two horizontal components of the real record from Table 1 for the station "Westside Elementary School" and five corresponding spectrum-compatible and energy-compatible synthetic acceleration time histories.
The synthetic ground acceleration time histories have been generated as spectrum-compatible and energy-compatible with the two selected real records, based on the method of Li et al. [50], see Equations 2 to 6, and by using the corresponding Matlab subroutines. The method of Li et al. [50] is partially based on the same procedure used by Gasparini and Vanmarcke [51], i.e., on the fact that any periodic function can be expanded into a series of sinusoidal waves:  Table 1 for the station "El Centro Array #12" and five corresponding spectrum-compatible and energy-compatible synthetic acceleration time histories; (b) Acceleration time histories of the two horizontal components of the real record from Table 1 for the station "Westside Elementary School" and five corresponding spectrum-compatible and energy-compatible synthetic acceleration time histories.
The synthetic ground acceleration time histories have been generated as spectrum-compatible and energy-compatible with the two selected real records, based on the method of Li et al. [50], see Equations (2)- (6), and by using the corresponding Matlab subroutines. The method of Li et al. [50] is partially based on the same procedure used by Gasparini and Vanmarcke [51], i.e., on the fact that any periodic function can be expanded into a series of sinusoidal waves: where A i . is the amplitude and ϕ i is the phase angle of the i-th contributing sinusoid. By fixing an array of amplitudes and generating different arrays of phase angles, one obtains different motions with the same general appearance but different details. The computer uses a random number generator to produce strings of phase angles with the uniform likelihood in the range between 0 and 2π. In order to match the target spectrum, the amplitudes A i have to be iteratively adjusted: where j is the number of the current iteration and PSA t is the target pseudo-acceleration spectrum.
To simulate the transient character of real earthquakes, the steady-state motions are multiplied by an envelope function I(t). The artificial motion z(t) then becomes: Usually, the envelope function is a deterministic function that does not change during the iteration process. In their method, Li et al. [50] adjust I(t) according to the differences between the target and simulated energy distributions. For this purpose, they proposed iterative adjustment of I(t) based on the values of Arias Intensity (AI). The Arias Intensity is the measure of the energy contained in the accelerogram and can be defined as: where t d is the total duration of strong ground motion and a g is the acceleration time history. In order to match the target AI, the envelope I(t) has to be iteratively adjusted: where j is the number of the current iteration, a g is the target pseudo-acceleration spectrum, and p is a factor controlling the convergence speed. Table 2 demonstrates the PGA and AI values of the synthetic accelerograms compatible with the real records for the stations "El Centro Array #12" and "Westside Elementary School". Furthermore, Figures 4 and 5, respectively, show pseudo-acceleration spectra and Arias intensities of the real and synthesized accelerograms. For more details on the whole procedure, the readers are referred to [50].   Table 1 for the station "El Centro Array #12" and the spectra of five corresponding synthetic acceleration time histories shown in Figure  1a, for both horizontal ground motion components. (b) Pseudo-acceleration spectrum of the real record from Table 1 for the station "Westside Elementary School" and the spectra of five corresponding synthetic acceleration time histories shown in Figure 1b, for both horizontal ground motion components.  Table 1 for the station "El Centro Array #12" and the spectra of five corresponding synthetic acceleration time histories shown in Figure 1a, for both horizontal ground motion components. (b) Pseudo-acceleration spectrum of the real record from Table 1 Table 1 for the station "El Centro Array #12" and the Arias intensities of five corresponding synthetic acceleration time histories shown in Figure 1a, for both horizontal ground motion components. (b) Arias intensity of the real record from Table 1 for the station "Westside Elementary School" and the Arias intensities of five corresponding synthetic acceleration time histories shown in Figure 1b, for both horizontal ground motion components.

The SMA-LRBs
Adding the SMA wires to base isolators as a supplementary element can mitigate the seismic response of the structures in terms of displacement and internal forces, and offers some benefits such as re-centering capability, stability, high fatigue resistance, high energy dissipation capacity, etc. For these reasons, many researchers have used SMA wires in conjunction with LRBs, which some of them are mentioned in the introduction.
In this study, the newly discovered iron-based SMA alloy, i.e., FeNiCoAlTaB with 13.5% superelastic strain is used, and the finite element approach is applied to investigate the effects of different wire configurations on the seismic performance of a benchmark bridge. To do so, some milestones should be reached, including the development of an appropriate constitutive model for the SMA-LRBs. Due to the complexity of the shear behavior of the SMA-LRBs, Hedayati Dezfuli and Alam [24] suggested using the idea of superposition to divide the smart isolation into two separate systems, namely, SMA wire and LRB to simplify the model. Then, to obtain the shear hysteresis model for the SMA-LRB, first, the model is found for each device; second, the models are added to each other to form the hysteresis behavior of the isolation system. It is worth mentioning that Hedayati Dezfuli and Alam verified the superposition method through the finite element model validated by experimental results [24].  Table 1 for the station "El Centro Array #12" and the Arias intensities of five corresponding synthetic acceleration time histories shown in Figure 1a, for both horizontal ground motion components. (b) Arias intensity of the real record from Table 1 for the station "Westside Elementary School" and the Arias intensities of five corresponding synthetic acceleration time histories shown in Figure 1b, for both horizontal ground motion components.

The SMA-LRBs
Adding the SMA wires to base isolators as a supplementary element can mitigate the seismic response of the structures in terms of displacement and internal forces, and offers some benefits such as re-centering capability, stability, high fatigue resistance, high energy dissipation capacity, etc. For these reasons, many researchers have used SMA wires in conjunction with LRBs, which some of them are mentioned in the introduction.
In this study, the newly discovered iron-based SMA alloy, i.e., FeNiCoAlTaB with 13.5% super-elastic strain is used, and the finite element approach is applied to investigate the effects of different wire configurations on the seismic performance of a benchmark bridge. To do so, some milestones should be reached, including the development of an appropriate constitutive model for the SMA-LRBs. Due to the complexity of the shear behavior of the SMA-LRBs, Hedayati Dezfuli and Alam [24] suggested using the idea of superposition to divide the smart isolation into two separate systems, namely, SMA wire and LRB to simplify the model. Then, to obtain the shear hysteresis model for the SMA-LRB, first, the model is found for each device; second, the models are added to each other to form the hysteresis behavior of the isolation system. It is worth mentioning that Hedayati Dezfuli and Alam verified the superposition method through the finite element model validated by experimental results [24].

The Hysteresis Model of the LRBs
In order to simplify the finite element analysis procedure, the hysteresis response of the LRBs can be modelled using a bilinear model [24,37,52,53]. Here, it is assumed that the behavior of LRB can be simulated by the bilinear kinematic hardening (BKH) model. Three parameters exist in this model, including yield force F y , initial stiffness (i.e., stiffness in the elastic region) K 0 and post-yield hardening ratio r (i.e., the ratio of the stiffness in the plastic region to the initial stiffness). These parameters are required for finding the shear force-displacement curve. The BKH model is available in some software. In this case, the available code in the benchmark bridge package is applied. Characteristics of benchmark bridge LRBs are shown in Table 3.

The Hysteresis Model of SMA Wires
Researchers have wrapped SMA wires around the base isolators in three different arrangements, namely, straight [18], cross [16,26], and double-cross [24,25]. These isolators are schematically shown in Figure 6. It should be noted that, in these isolators, wires pass through steel hooks, which are mounted on the top and the bottom supporting plates. It is assumed that there is a smooth contact between the steel hook and the SMA wire. In order to reduce the complexity of the finite element simulation, Hedayati Dezfuli and Alam [16] proposed to consider the applied forces on the isolator due to the SMA wires, instead of modelling the steel hooks and their contact with the wire. They calculated the force generated in SMA wires as a function of time, and then applied the effect of SMA wire on the isolator. For calculating this axial force of the wire, first, by specifying the geometrical properties of the isolator and configuration of wire, the strain in the wire at each time step can be found. Second, with the help of the stress-strain relationship of SMA, the axial stress in the wire is determined. Finally, using the amount of the axial stress and the direction of wire at each time step, the force vectors exerted from the SMA wire to the LRB isolator as an external force are computed.

The Hysteresis Model of the LRBs
In order to simplify the finite element analysis procedure, the hysteresis response of the LRBs can be modelled using a bilinear model [24,37,52,53]. Here, it is assumed that the behavior of LRB can be simulated by the bilinear kinematic hardening (BKH) model. Three parameters exist in this model, including yield force , initial stiffness (i.e., stiffness in the elastic region) , and post-yield hardening ratio r (i.e., the ratio of the stiffness in the plastic region to the initial stiffness). These parameters are required for finding the shear force-displacement curve. The BKH model is available in some software. In this case, the available code in the benchmark bridge package is applied. Characteristics of benchmark bridge LRBs are shown in Table 3.

The Hysteresis Model of SMA Wires
Researchers have wrapped SMA wires around the base isolators in three different arrangements, namely, straight [18], cross [16,26], and double-cross [24,25]. These isolators are schematically shown in Figure 6. It should be noted that, in these isolators, wires pass through steel hooks, which are mounted on the top and the bottom supporting plates. It is assumed that there is a smooth contact between the steel hook and the SMA wire. In order to reduce the complexity of the finite element simulation, Hedayati Dezfuli and Alam [16] proposed to consider the applied forces on the isolator due to the SMA wires, instead of modelling the steel hooks and their contact with the wire. They calculated the force generated in SMA wires as a function of time, and then applied the effect of SMA wire on the isolator. For calculating this axial force of the wire, first, by specifying the geometrical properties of the isolator and configuration of wire, the strain in the wire at each time step can be found. Second, with the help of the stress-strain relationship of SMA, the axial stress in the wire is determined. Finally, using the amount of the axial stress and the direction of wire at each time step, the force vectors exerted from the SMA wire to the LRB isolator as an external force are computed. As mentioned above, the strain in SMA wire at each time step depends on the geometry of the isolator and wires arrangement. In this study, the SMA wires are wrapped around the LRBs of the benchmark bridge, so the geometry of the LRB is constant. Therefore, calculations must be performed for the three different configurations of wires separately. As mentioned above, the strain in SMA wire at each time step depends on the geometry of the isolator and wires arrangement. In this study, the SMA wires are wrapped around the LRBs of the benchmark bridge, so the geometry of the LRB is constant. Therefore, calculations must be performed for the three different configurations of wires separately.

SMA-LRB with Double-Cross Wires
In double-cross arrangement, each of the two continuous SMA wires is passed through eight steel hooks welded to the top and bottom supporting plates of the LRB, as can be seen in Figure 6. Hedayati Dezfuli and Alam proposed and utilized this configuration of wires to increase the length of the SMA wire with the purpose of reducing the wire's strain. For an LRB with specific dimensions, the required SMA wire to wrap around the isolator is the longest for the double-cross arrangement. As a result, the effective strain in the SMA wires generated due to the shear deformation in the LRB decreases because wires have a longer initial length.
The strain in the SMA wire is calculated by Equation (7). Where l 0,SMA is the initial length of the SMA wire and l SMA is the length of the SMA wire after the isolator undergoes a deflection of ∆X in the X-direction. For the double-cross arrangement of wires, these parameters are determined from Equations (8) and (9). In which l, w, and h are the distances between hooks in X, Y, and Z directions, respectively. Having calculated the strain, the axial stress in SMA wires can be determined from the idealized stress-strain relationship of shape memory alloy based on Auricchio's super-elasticity model [54]. The idealized super-elastic behavior of SMA is shown in Figure 7. After that, the axial force in the wire, F, can be found from Equation (10). The resultant force exerted to the LRB in the X-direction, and the consequent hysteresis of the SMA wires is comprehensively explained in [24] and is not repeated here. In order to achieve the hysteresis curves of the SMA wire, a MATLAB function is developed by the authors based on the presented algorithm by Hedayati Dezfuli and Alam and was verified in a previous study [23]. This algorithm calculates the force applied by the SMA wire to the LRB at each time step for a specific shear displacement.

SMA-LRB with Double-Cross Wires
In double-cross arrangement, each of the two continuous SMA wires is passed through eight steel hooks welded to the top and bottom supporting plates of the LRB, as can be seen in Figure 6. Hedayati Dezfuli and Alam proposed and utilized this configuration of wires to increase the length of the SMA wire with the purpose of reducing the wire's strain. For an LRB with specific dimensions, the required SMA wire to wrap around the isolator is the longest for the double-cross arrangement. As a result, the effective strain in the SMA wires generated due to the shear deformation in the LRB decreases because wires have a longer initial length.
The strain in the SMA wire is calculated by Equation (7). Where , is the initial length of the SMA wire and is the length of the SMA wire after the isolator undergoes a deflection of ΔX in the X-direction. For the double-cross arrangement of wires, these parameters are determined from Equations (8) and (9). In which l, w, and h are the distances between hooks in X, Y, and Z directions, respectively. Having calculated the strain, the axial stress in SMA wires can be determined from the idealized stress-strain relationship of shape memory alloy based on Auricchio's super-elasticity model [54]. The idealized super-elastic behavior of SMA is shown in Figure 7. After that, the axial force in the wire, F, can be found from Equation (10). The resultant force exerted to the LRB in the Xdirection, and the consequent hysteresis of the SMA wires is comprehensively explained in [24] and is not repeated here. In order to achieve the hysteresis curves of the SMA wire, a MATLAB function is developed by the authors based on the presented algorithm by Hedayati Dezfuli and Alam and was verified in a previous study [23]. This algorithm calculates the force applied by the SMA wire to the LRB at each time step for a specific shear displacement.

SMA-LRB with Straight Wires
This arrangement of wires was proposed by Choi et al. [18]. As can be seen in Figure 6, in straight configuration, two continuous SMA wires are wound in two opposite sides of the isolator. Wires are passed through steel hooks with smooth contact. Since the wire is continuous and is not fixed at the corners of the isolator, the initial length of the SMA wire is longer, and its effective strain is lower.
In this arrangement the strain in SMA wire is also determined by Equation (7); however, the parameters for straight wires are obtained based on its configuration as in Equations (11) and (12). After that, the strain and axial force in the SMA wire is determined like in the double-cross arrangement. The proposed algorithm by Hedayati Dezfuli and Alam [24] can be modified slightly

SMA-LRB with Straight Wires
This arrangement of wires was proposed by Choi et al. [18]. As can be seen in Figure 6, in straight configuration, two continuous SMA wires are wound in two opposite sides of the isolator. Wires are passed through steel hooks with smooth contact. Since the wire is continuous and is not fixed at the corners of the isolator, the initial length of the SMA wire is longer, and its effective strain is lower.
In this arrangement the strain in SMA wire is also determined by Equation (7); however, the arameters for straight wires are obtained based on its configuration as in Equations (11) and (12). After that, the strain and axial force in the SMA wire is determined like in the double-cross arrangement. The proposed algorithm by Hedayati Dezfuli and Alam [24] can be modified slightly to account for the straight configuration. This code is written by the authors as a function in MATLAB to input ∆X at each time step and calculate the resultant force that SMA wire exerts to the LRB.

SMA-LRB with Cross Wires
In the cross configuration, two continuous SMA wires are wound around the LRB diagonally, as can be seen in Figure 6. Hedayati Dezfuli and Alam firstly suggested and used this configuration of wires [16,30]. The reason was that the wire in this figure is longer compared with the straight arrangement, and thus, the induced strain in the wire is lower. In this case, again, the wire's strain is calculated by Equation (2), but the parameters of this equation are determined by Equations (13) and (14). Similar to the previous two configurations, by modifying the Hedayati Dezfuli and Alam [24] algorithm, a new MATLAB function is obtained that can calculate the hysteresis shear response of SMA wires with cross configuration.

Efficiency of Different Wire Configurations
An essential prerequisite for the efficient performance of SMA-LRB is for the SMA wires to stay in their super-elastic range under lateral displacements. The SMA alloy that is used in this study is FeNiCoAlTaB, with 13.5% maximum super-elastic strain. As explained above, the strain in the wire is a function of wires' configuration and LRB's dimension and is independent of the wire's radius. Here, LRBs geometry is constant; therefore, the strain is only dependent on the wires' arrangement.
In order to find the range of shear strain that each configuration can remain applicable, the induced strain in wires due to a continuous range of shear strain (i.e., lateral displacement amplitude divided by total rubber height), γ, ranging between 0% and 250% calculated and plotted for all arrangements of wires in Figure 8. The wire strain is calculated using Equations (7) to (14). As can be observed, shear strain in the cross configuration is lower than in the others, but the performance of double-cross and straight wires are almost the same. At the ultimate shear strain, the induced strain in straight, cross and double-cross wire is 10.44%, 3.23%, and 11.39%, respectively, which is below the super-elastic strain limit of the iron-based SMA alloy. Thus, with this LRB dimension and this SMA alloy, wires in all three configurations can perform effectively.
It is worth mentioning that, for shear strain levels between 50% and 150%, the strain in SMA wire with cross arrangement varies between 0.13% and 1.18%, which is not high enough to create the flag-shaped hysteresis of SMA. Therefore, in lower shear strains, the cross configuration is not as effective as the other arrangements.

Force-Displacement Curves
In the next section, the radius of the SMA wire is found for each configuration through an optimization algorithm. However, the optimized radius is different for each arrangement. Thus, to observe and compare the sole effect of various arrangements on the hysteresis diagram of the SMA wires, the force-displacement curves for the three wire's configurations with a fixed wire radius of 2.5 mm, under three different lateral displacement scenarios, are plotted in Figure 9. As illustrated, for specified excitation and wire radius, the stiffness of the double-cross configuration is more than the straight one and more than the cross one. By increasing the radius of the wire in each arrangement, the stiffness could be increased. Thus, it is expected that SMA-LRBs with cross arrangement need thicker wires, and with double-cross arrangement, require the thinnest wire. Additionally, the cross wires create less flag-shaped hysteresis because of their lower levels of strain in the SMA wires (see Figure 8). The input displacements are extracted from Hedayati Dezfuli and Alam [24] and are named in the same way as E5, E6, and E7. They selected these displacement scenarios in a way to follow different functions like ramp, step, sinusoidal, and a combination of them with various peak amplitudes to assess the model in any condition. As mentioned before, the written code in this study is based on the algorithm developed by Hedayati Dezfuli and Alam [24] and was verified in a previous paper by Hosseini et al. [23].

Force-Displacement Curves
In the next section, the radius of the SMA wire is found for each configuration through an optimization algorithm. However, the optimized radius is different for each arrangement. Thus, to observe and compare the sole effect of various arrangements on the hysteresis diagram of the SMA wires, the force-displacement curves for the three wire's configurations with a fixed wire radius of 2.5 mm, under three different lateral displacement scenarios, are plotted in Figure 9. As illustrated, for specified excitation and wire radius, the stiffness of the double-cross configuration is more than the straight one and more than the cross one. By increasing the radius of the wire in each arrangement, the stiffness could be increased. Thus, it is expected that SMA-LRBs with cross arrangement need thicker wires, and with double-cross arrangement, require the thinnest wire. Additionally, the cross wires create less flag-shaped hysteresis because of their lower levels of strain in the SMA wires (see Figure 8). The input displacements are extracted from Hedayati Dezfuli and Alam [24] and are named in the same way as E5, E6, and E7. They selected these displacement scenarios in a way to follow different functions like ramp, step, sinusoidal, and a combination of them with various peak amplitudes to assess the model in any condition. As mentioned before, the written code in this study is based on the algorithm developed by Hedayati Dezfuli and Alam [24] and was verified in a previous paper by Hosseini et al. [23].

Optimizing the SMA-LRBs
In the previous section, the approach for obtaining the hysteresis models for SMA-LRBs with three different configurations of wires was explained. As discussed, the influential parameters on the hysteresis behavior of the SMA-LRBs include the geometry of the isolator, wires arrangement, and the radius of the SMA wire (see Equation (10)). This research uses the SMA wire in three

Optimizing the SMA-LRBs
In the previous section, the approach for obtaining the hysteresis models for SMA-LRBs with three different configurations of wires was explained. As discussed, the influential parameters on the hysteresis behavior of the SMA-LRBs include the geometry of the isolator, wires arrangement, and the radius of the SMA wire (see Equation (10)). This research uses the SMA wire in three arrangements, as a supplementary element to the existed LRBs of the benchmark bridge, to enhance its seismic performance without changing its isolators. Therefore, the only unknown yet significant parameter is the radius of the wire. It should be determined carefully since it has a substantial effect on the force applied to the LRB isolator. By increasing the radius of the wire, the re-centering force generated in the SMA wire and applied to the LRB increases, and, as a result, the horizontal stiffness of the isolator increases. Consequently, the residual deformation declines, but the base shear force and the deck acceleration of the bridge rise. The main goal of this research is to minimize the residual displacements of the bridge after strong earthquakes; however, a detrimental increase in the base shear and acceleration is not acceptable. Therefore, it is decided to take advantage of a multi-objective optimization algorithm to simultaneously minimize displacement and acceleration and find the optimum radius of the wire.
An evolutionary optimization algorithm named NSGA-II (non-dominated sorting genetic algorithm) is used to find the best radius of the SMA wire for each arrangement under the most severe earthquake. The NSGA-II is a multi-objective genetic algorithm that was proposed by Deb et al. [55]. It is an extension and development of NSGA, which was proposed previously by Srinivas and Deb [56]. In the structure of the NSGA-II, in addition to genetic operators (i.e., crossover and mutation), two specified operators are defined and employed. The first one is non-dominated sorting, which sorts and partitions the population into fronts (F1, F2, etc.), where F1 indicates the approximated Pareto front. The second one is crowding distance, which specifies a mechanism of ranking amongst members of a front, which are either dominated by or are dominating each other. These ranking mechanisms are utilized with genetic selection operators to generate the population of the next generation. The flowchart of NSGA-II is demonstrated in Figure 10  The crossover and mutation rates are 0.7 and 0.4, in turn. Under these defined circumstances, the algorithm is run three times for the three various configurations of wires, and the Pareto front is achieved for each case. The Pareto front indicates the points that simultaneously and independently minimize the normalized maximum mid-span displacement and the normalized maximum deck acceleration. Figure 11 depicts the optimized point for the three wire configurations. Selecting the desired radius amongst the optimized points of a Pareto front depends on the designer's perspective. In this research, the main goal is to decrease displacements; hence, the final radius is determined based on this idea. The final optimized and rounded radii for the double-cross, straight, and cross arrangements are 5 mm, 6 mm and 7 mm, respectively. It should be mentioned that, by increasing the population size, this algorithm can find more counts of optimum points at the expense of increasing the runtime.  In the NSGA-II algorithm of this project, there are two objective functions, namely, maximum mid-span displacement and maximum deck acceleration. For simplicity, these parameters are normalized with respect to the quantities of the uncontrolled bridge. The optimization should be carried out in the worst scenario of the earthquake, so the real record from the station El Centro Array #12 (RSN 8161) is used in the algorithm, as it has larger Arias Intensity and larger PGA value (Tables 1  and 2). The population includes 30 members, and the maximum number of iterations is equal to 30. The crossover and mutation rates are 0.7 and 0.4, in turn. Under these defined circumstances, the algorithm is run three times for the three various configurations of wires, and the Pareto front is achieved for each case. The Pareto front indicates the points that simultaneously and independently minimize the normalized maximum mid-span displacement and the normalized maximum deck acceleration. Figure 11 depicts the optimized point for the three wire configurations. Selecting the desired radius amongst the optimized points of a Pareto front depends on the designer's perspective. In this research, the main goal is to decrease displacements; hence, the final radius is determined based on this idea. The final optimized and rounded radii for the double-cross, straight, and cross arrangements are 5 mm, 6 mm and 7 mm, respectively. It should be mentioned that, by increasing the population size, this algorithm can find more counts of optimum points at the expense of increasing the runtime.
increasing the runtime.

Numerical Analysis
In order to assess the performance of the three optimized SMA-LRBs, the seismic responses of the highway bridge equipped with the SMA-based isolators are calculated and compared with those

Numerical Analysis
In order to assess the performance of the three optimized SMA-LRBs, the seismic responses of the highway bridge equipped with the SMA-based isolators are calculated and compared with those of the LRB isolated and the non-isolated bridges. These responses include the base shear at the piers, the deck acceleration, relative displacement at the mid-span of the bridge, and the mid-span residual deformation. For the purpose of the comparison, five distinctive sets of dynamic time history analyses must be conducted on the (1) the non-isolated bridge, (2) LRB isolated bridge, (3) SMA-LRB isolated with double-cross wires, (4) SMA-LRB isolated with straight wires and finally (5) SMA-LRB isolated with cross wires. All these models are excited under two real and 10 synthetic ground motion records. The detailed characteristics of these records are mentioned in Section 3. These records are exerted on the bridge only in the longitudinal direction, and between the two horizontal components for each record, the one that leads to more critical responses is chosen. It is worth mentioning that the responses of the non-isolated bridge are calculated by analyzing the non-isolated model (phase I) of the benchmark bridge, which was developed and released by Agrawal et al. [27].

Base Shear
The peak base shear responses of the benchmark bridge in five different conditions, i.e., non-isolated bridge, LRB isolated bridge, and SMA-LRB isolated bridge with double-cross (denoted as DC-SMA-LRB), straight (denoted as S-SMA-LRB) and cross (denoted as C-SMA-LRB) wire arrangements resulted from the nonlinear dynamic time history analysis are tabulated in Table 4. As can be seen in this Table, base shear for the non-isolated benchmark bridge is by far the largest under all ground motion records. After implementing LRB isolators, the base shear decreases significantly due to the reduction in the stiffness and transmitted load from the base to the superstructure. The largest drop in the base shear after isolating by the LRBs occurs under the real record of RSN 8606 earthquake, with 91%. However, by isolating the superstructure, its displacements increase. That is why SMA wires are employed in this study, to control the excessive displacements of the bridge. Adding SMA wires to the LRBs leads to a marginal rise in the horizontal stiffness and, as a consequence, in the base shear. According to Figure 12, it is true that SMA-LRBs increase the base shear of the benchmark bridge more than the LRBs, but the figures are still far less than those of the non-isolated bridge. Consequently, both the LRBs and SMA-LRBs are effective measures for reducing the base shear of the bridge.   Table 5 and Figure 13 show and compare the maximum mid-span acceleration in the nonisolated, LRB isolated, and three SMA-LRB isolated benchmark bridges. As reflected, both isolation systems, i.e., LRB isolators and SMA-LRBs substantially decrease the peak acceleration under all ground motion records, compared to the non-isolated state. By introducing isolators in the central pier of the benchmark bridge, its stiffness and load transferred from the base to the superstructure reduce, so the acceleration reduces too. After wrapping SMA wires around the LRB isolators, however, the stiffness and the acceleration increase slightly. This increase is the least for SMA-LRBs with cross wires and the largest for the double-cross arrangement of the wires, and still far less than the figures for the non-isolated bridge. Furthermore, as an example, the time history curves of the mid-span acceleration for the five different benchmark bridges under the synthetic record number 3 of the RSN 8161 earthquake are plotted in Figure 14. Table 5. Maximum mid-span acceleration for the non-isolated, LRB-isolated and three SMA-LRB isolated bridges.  Table 5 and Figure 13 show and compare the maximum mid-span acceleration in the non-isolated, LRB isolated, and three SMA-LRB isolated benchmark bridges. As reflected, both isolation systems, i.e., LRB isolators and SMA-LRBs substantially decrease the peak acceleration under all ground motion records, compared to the non-isolated state. By introducing isolators in the central pier of the benchmark bridge, its stiffness and load transferred from the base to the superstructure reduce, so the acceleration reduces too. After wrapping SMA wires around the LRB isolators, however, the stiffness and the acceleration increase slightly. This increase is the least for SMA-LRBs with cross wires and the largest for the double-cross arrangement of the wires, and still far less than the figures for the non-isolated bridge. Furthermore, as an example, the time history curves of the mid-span acceleration for the five different benchmark bridges under the synthetic record number 3 of the RSN 8161 earthquake are plotted in Figure 14.

Mid-Span Displacement
The benchmark highway bridge in the non-isolated phase, in which the deck is connected to the central pier, experiences the smallest superstructure displacements, and the highest base shear and deck acceleration (see Sections 6.1 and 6.2). After isolating the deck from the central piers via LRBs, the acceleration, and base shear drop, but mid-span displacement rises dramatically. This rise is evident in Figure 15 and Table 6. For example, under synthetic record number 3 of the RSN 8606 record, the maximum mid-span displacement jumps from 81 to 650 mm after isolation. The extreme

Mid-Span Displacement
The benchmark highway bridge in the non-isolated phase, in which the deck is connected to the central pier, experiences the smallest superstructure displacements, and the highest base shear and

Mid-Span Displacement
The benchmark highway bridge in the non-isolated phase, in which the deck is connected to the central pier, experiences the smallest superstructure displacements, and the highest base shear and deck acceleration (see Sections 6.1 and 6.2). After isolating the deck from the central piers via LRBs, the acceleration, and base shear drop, but mid-span displacement rises dramatically. This rise is evident in Figure 15 and Table 6. For example, under synthetic record number 3 of the RSN 8606 record, the maximum mid-span displacement jumps from 81 to 650 mm after isolation. The extreme deck displacement resulting from the isolating process can have detrimental effects, such as the collision of the adjacent parts of the superstructure, deck unseating, or plastic deformation of the bearings. To tackle this issue, SMA wires are added to the LRBs to control displacements.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 20 of 28 cross wires reduce displacements more than the straight wires, and the cross configuration has the least displacement reduction.    As illustrated in Table 6 and Figure 15, SMA wires are quite beneficial in minimizing the mid-span displacements. For instance, the maximum mid-span displacement under synthetic record number 3 of the RSN 8606 decreases from 650 to 122 mm using double-cross SMA-LRBs. According to Table 6, SMA-LRBs with various wire configurations reduce the maximum mid-span displacements in a range from 26% to 81%. Another interesting point is that SMA-LRBs with double-cross wires reduce displacements more than the straight wires, and the cross configuration has the least displacement reduction.

Residual Mid-Span Displacement
One of the main goals of implementing SMA wires in LRBs is to reduce the residual displacement of the deck after strong earthquakes and decrease the repairing cost. Residual displacement is the lateral displacement that remains in the deck at the end of the earthquake. This value can be extracted from the displacement history curve of each point of the bridge's deck. Here, the mid-span residual displacement is calculated, which also represents the residual displacement at the place of the mid-pier isolators. The mid-span residual displacement under all ground motion records, for the LRB isolated and the three SMA-LRB isolated bridges, is calculated and illustrated in Table 7. Besides, the reduction ratios of the residual displacement of the SMA-LRBs with respect to the LRBs are shown. The most significant reductions are under synthetic records number 3 of the RSN 8606 and RSN 816 records, for which the displacement history curves are demonstrated in Figures 16 and 17. As indicated, the mid-span residual displacement in the LRB isolated benchmark bridge is 13.19 cm under synthetic record number 3 RSN 8606, and by using double-cross, straight, and cross SMA wires, it drops to 0.18 cm, 0.89 cm, and 2.76 cm, respectively.

Comparing the Three SMA-LRBs
SMA wires can be wrapped around the LRB isolators in double-cross, straight, and cross arrangements to improve the seismic performance of the isolated structure. On the one hand, it is Due to the self-centering capability of the SMA wires, the residual displacement of the benchmark bridge equipped with SMA-LRBs is much smaller than that of the bridge with LRBs. Furthermore, amongst the three wires' arrangements, the optimized double-cross is the most successful device in reducing the residual displacement, from 52.16% to 98.67% under various records. This range for the straight arrangement is between 51.55% and 93.25%, and for the cross wires is from 42.51% to 82.17%. Figure 18 graphically compares the mid-span residual displacement for different isolators.

Comparing the Three SMA-LRBs
SMA wires can be wrapped around the LRB isolators in double-cross, straight, and cross arrangements to improve the seismic performance of the isolated structure. On the one hand, it is

Comparing the Three SMA-LRBs
SMA wires can be wrapped around the LRB isolators in double-cross, straight, and cross arrangements to improve the seismic performance of the isolated structure. On the one hand, it is proven that adding SMA wires to the LRBs, results in a marginal rise in the base shear and deck acceleration. As evidenced in Tables 4 and 5, this rise is the least for the bridge with cross SMA-LRBs and is the highest for the bridge with the double-cross SMA-LRBs. On the other hand, adding SMA wires reduces the maximum and residual displacements because of their re-centering and super-elasticity capabilities. According to Tables 6 and 7, amongst the three configurations, cross wires have the least reduction ratio in the maximum and residual displacements, while double-cross wires have the highest ratio. Another difference between these three isolators is the length of the used SMA wire, which can affect the cost of the project. According to Equations (8), (11), and (13), for SMA-LRBs with double-cross, straight and cross configurations, 242.71 cm, 133.12 cm, and 205.85 cm SMA wire is required to wrap around each of the LRBs, respectively.
For a specified design perspective, in order to choose the best configuration of SMA wires to add to the LRBs, the weighted sum method (WSM) can be utilized. The WSM is the best known and commonly used multi-criteria decision-making method for evaluating a number of alternatives regarding a number of decision criteria. According to the WSM, in a decision problem with M alternatives and N criteria, a performance value can be calculated for each alternative by applying Equation (15). In which, A WSM−Score i is the WSM score of the i th . alternative, w j is the weight of importance of j th criteria, and a ij is the actual value of the i th alternative regarding the j th criterion.
a ijNormalized = Min a ij a ij (16) In this multi-criteria decision-making problem, there are three alternatives, including SMA-LRB with double-cross, straight and cross wires, and five decision criteria are set as the maximum base shear, maximum deck acceleration, maximum and residual mid-span displacement and the length of the SMA wire. Moreover, the responses under the real record RSN 8161, as illustrated in Table 8, are considered here, since this record has larger PGA values. To calculate the score of each alternative, first of all, the values must be normalized. Since for all the five criteria, the lower value is desired, then normalization must be performed using Equation (16). The normalized decision matrix is shown in Table 8.   The next step is to assign weights to the criteria. These numbers are selected based on the design perspective. Here, for example, reducing displacements is more important than a minor increase in the shear force and acceleration, and the length of the wire is presumed as the least important criteria. As a result, weights are selected as 0.15, 0.15, 0.3, 0.3, and 0.1 for base shear, acceleration, mid-span displacement, residual displacement, and wire length, respectively. The final step is to substitute the normalized values in Equation (15) and calculate the performance score for each alternative. As illustrated in Table 8, the WSM score for the SMA-LRB with double-cross, straight, and cross wires is 0.92, 0.93, and 0.74, in turn. Consequently, the straight wire arrangement has the best overall performance (rank 1), just above the double-cross arrangement (rank 2), but the SMA-LRB with cross wires displays by far the weakest performance (rank 3).

Conclusions
For nearly two decades, researchers have been using SMA wires in conjunction with isolators of the bridges to control their displacements and residual deformations. Different configurations of the SMA wires around the LRB isolators have been reported, namely, double-cross, straight, and cross. In this paper, the performance of these three SMA-LRBs is compared by implementing them in the LRBs of a benchmark highway bridge. Each configuration leads to specific characteristics of the SMA-LRBs and thus specific shear hysteresis. In addition to the wire's arrangements, LRB dimensions and the radius of the wire are also influential parameters in the SMA-LRB behavior. For a specified LRB dimension, by increasing the radius of the wire, the horizontal stiffness, and consequently, the base shear and deck acceleration of the bridge rise, but the maximum and residual displacement reduce.
To get the best performance out of SMA wires, a multi-objective non-dominated sorting genetic algorithm (NSGA-II) is applied to find the radius of the SMA wire for each configuration that simultaneously minimizes the deck acceleration and mid-span displacement. After optimizing the isolators, to evaluate and compare their performance, nonlinear dynamic analysis is conducted for the five distinctive states of the benchmark bridge, i.e., non-isolated, LRB isolated, and the three SMA-LRB isolated bridges. For the nonlinear response history analysis, researchers can select ground motions from past recorded earthquakes in places other than the bridge's site and modify them to match desired conditions. However, because of the deficiencies of this practice, in this paper, strong ground motion records are selected in a way to fit the site conditions, and also ten synthetic spectrum-compatible and energy-compatible ground acceleration time histories are generated to be used in the analysis. The summarized results are as follows.

•
As for the optimized double-cross SMA-LRB, this device decreases the maximum base shear and deck acceleration of the isolated bridge from 70.29% to 87.85% and from 20.89% to 62.57%, respectively, under various records compared with the non-isolated bridge. However, LRB still leads to the largest reduction in these figures, from 77.71% to 91.07%, and from 56.85% to 82.66% under various records. Moreover, SMA-LRBs with double-cross wires reduce the maximum mid-span and residual displacements in ranges between 39.51% and 81.23% and from 52.16% to 98.67%, in turn, under different ground motions compared with the LRB isolated bridge. In brief, amongst the three SMA-LRBs, it has the largest shear force and deck acceleration, and needs the longest SMA wire, but leads to the least maximum and residual displacements. Finally, the result of the WSM method indicates that DC-SMA-LRB is the second most successful device.

•
Regarding the optimized straight SMA-LRB, it decreases the maximum base shear and deck acceleration of the isolated bridge slightly more than the double-cross SMA-LRBs, from 72.50% to 88.28% and from 26.76% to 66.90%, respectively, under different records, compared with the non-isolated benchmark bridge. Besides, SMA-LRBs with straight wires reduce the maximum and residual mid-span displacement slightly less than the DC-SMA-LRBs, from 37.06% to 80.00% and from 51.55% to 94.67%, in turn, under various excitations, compared with the LRBs. Moreover, the required length of the SMA wire is the least for this configuration. Lastly, the WSM method shows that S-SMA-LRB has the best overall performance.

•
As regards the optimized cross SMA-LRB, it results in the reduction of the base shear and deck acceleration more than the other two SMA isolators, compared with the non-isolated bridge. However, it decreases the maximum and residual displacements less than the others and uses rather long SMA wires. Eventually, by conducting WSM, it is ranked as the least successful device.
In this research, the bridge model presented by Nagarajaiah et al. [28] is utilized in which the bearings are modelled as shear elements by infinite vertical stiffness and zero torsional rigidity, and the coupling between the horizontal and vertical responses is neglected. In future studies, more up-to-date models for the bearings should be utilized [58]. Moreover, a sensitivity analysis of the optimization methodology is needed to be done. A further future research proposal is to carry out probabilistic seismic fragility analysis or risk-based assessment on the optimized SMA-LRB isolated bridge [59,60].