Research Article  Open Access
Seismic Response Analysis of Continuous Multispan Bridges with Partial Isolation
Abstract
Partially isolated bridges are a particular class of bridges in which isolation bearings are placed only between the piers top and the deck whereas seismic stoppers restrain the transverse motion of the deck at the abutments. This paper proposes an analytical formulation for the seismic analysis of these bridges, modelled as beams with intermediate viscoelastic restraints whose properties describe the pierisolator behaviour. Different techniques are developed for solving the seismic problem. The first technique employs the complex mode superposition method and provides an exact benchmark solution to the problem at hand. The two other simplified techniques are based on an approximation of the displacement field and are useful for preliminary assessment and design purposes. A realistic bridge is considered as case study and its seismic response under a set of ground motion records is analyzed. First, the complex mode superposition method is applied to study the characteristic features of the dynamic and seismic response of the system. A parametric analysis is carried out to evaluate the influence of support stiffness and damping on the seismic performance. Then, a comparison is made between the exact solution and the approximate solutions in order to evaluate the accuracy and suitability of the simplified analysis techniques for evaluating the seismic response of partially isolated bridges.
1. Introduction
Partially restrained seismically isolated bridges (PRSI) are a particular class of bridges isolated at intermediate supports and transversally restrained at abutments. This type of partial isolation is quite common in many bridges all over the world (see, e.g., [1–4]). The growing interest on the dynamic behaviour of bridges with partial restrain is demonstrated by numerous recent experimental works [3, 4] and numerical studies [5–11] discussing the advantages and drawbacks with respect to full isolation.
An analytical model commonly employed for the analysis of the transverse behaviour of PRSI bridges consists in a continuous simply supported beam resting on discrete intermediate supports with viscoelastic behaviour representing the pierbearing systems (see, e.g., [6, 8, 9]). The damping is promoted by two different mechanisms: the bearings, usually characterized by high dissipation capacity, and the deck, characterized by a lower but widespread dissipation capacity. The strongly inhomogeneous distribution of the dissipation properties along the bridge results in nonclassical damping and this makes the rigorous analysis of the analytical PRSI bridge model very demanding.
In general, the exact solution of the seismic problem for a nonclassically damped system requires resorting to the direct integration of the equations of motion or to modal analysis [12] based on complex vibration modes. The first analysis approach is conceptually simple, though computationally costly when largescale systems are analyzed, whereas the second approach is computationally efficient but often not appealing for practical engineering applications since it involves complexvalued functions and it is difficult to implement in commercial finite element codes. For this reason, many studies have been devoted to the definition of approximate techniques of analysis and to the assessment of their accuracy [13–17]. With reference to the classical case of fully isolated bridges, the studies of Hwang et al. [18] and Lee et al. [19] have shown that acceptable estimates of the modal damping ratios are obtained by considering the real undamped modes and the diagonal terms of the modal damping matrix. Franchin et al. [20] have analyzed some realistic bridge models by showing that this decoupling approximation gives also quite accurate estimates of the response to a seismic input, as compared to the rigorous response estimates obtained through complex modal analysis. With reference to the specific case of PRSI bridges, characterized by deformation shapes and dynamic properties very different from those of fully isolated bridges, the work of Tubaldi and Dall’Asta [9] has addressed the issue of nonclassical damping within the context of the freevibration response. The authors have observed that nonclassical damping influences differently the various response parameters relevant for the performance assessment (i.e., the transverse displacement shape is less affected than the bending moment demand by the damping nonproportionality). However, the effects of the decoupling approximation on the evaluation of the seismic response of the proposed PRSI bridge model have not been investigated yet. This issue becomes of particular relevance in consequence of the numerous studies on PRSI bridges that completely disregard nonclassical damping [6–8]. Thus, a closer examination is still required to ensure whether the use of proportionally damped models provides acceptable estimates of the seismic response of these systems.
Another approximation often introduced in the analysis of PRSI bridges concerns the transverse deformation shape. In this regard, many studies are based on the assumption of a prefixed sinusoidal vibration shape [6, 8, 10]. This assumption permits deriving analytically the properties of a generalized SDOF system equivalent to the bridge and estimating the system response by expressing the seismic demand in terms of a response spectrum reduced to account for the system composite damping ratio. In Tubaldi and Dall’Asta [8, 9], it is shown that the sinusoidal approximation of the transverse displacements may be accurate for PRSI bridges if the following conditions are met: (a) the superstructure stiffness is significantly higher than the pier stiffness, (b) the variations of mass and stiffness of the deck and of the supports are not significant, (c) the span number is high, and (d) the displacement field is dominated by the first vibration mode. However, even if the displacements are well described by a sinusoidal shape, other response parameters of interest for the performance assessment such as the transverse bending moments may not exhibit a sinusoidal shape. Thus, further investigations are required to estimate the error arising due to this approximation.
The aim of this study is to develop an analytical formulation of the seismic problem of PRSI bridges, modelled as nonclassically damped continuous systems, and a rigorous solution technique based on the complex mode superposition (CMS) method [12, 21–24]. The application of this method requires the derivation of the expression of the modal orthogonality conditions and of the impulsive response specific to the problem at hand. It permits describing the seismic response in terms of superposition of the complex vibration modes, which is particularly useful for this case in consequence of the relevant contribution of higher modes of vibration to the response of PRSI bridges [8, 10]. Moreover, the proposed technique permits testing the accuracy of two simplified analysis techniques introduced in this study and commonly employed for the PRSI bridge analysis. The first approximate technique describes the transverse motion through a series expansion in terms of the classic modes of vibration, obtained by neglecting the damping of the intermediate restraints, whereas the second technique is based on a series expansion in terms of sinusoidal functions, corresponding to the vibration modes of the system without the intermediate restraints. The introduction of these approximations of the displacement field results in a coupling of the equation of motions projected in the space of the approximating functions, which is neglected in the solution to simplify the response assessment.
A realistic case study is considered and its response to a set of ground motion records is examined first through the CMS method with these two objectives: to unveil the characteristic features of the performance of PRSI bridges and to assess the influence of higher vibration modes and of the intermediate support stiffness and damping on the response of the resisting components. Successively, the seismic response estimates according to the CMS method are compared with the corresponding estimates obtained by applying the proposed simplified techniques, in order to evaluate their accuracy and reliability.
2. Dynamic Behavior of PRSI Bridges
The PRSI bridge model (Figure 1) consists of beam pinned at the abutments and resting on discrete viscoelastic supports representing the pierbearing systems.
Let be the space of functions with square integrable second derivatives in the spatial interval , the space of transverse displacement functions satisfying the kinematic boundary conditions (e.g., ), and the motion, defined in the time interval considered , belonging to the space of continuous functions and known at the initial instant together with its time derivative (initial conditions). The differential dynamic problem can be derived from the D’Alembert principle [25] and expressed in the following form:where denotes a virtual displacement consistent with the geometric restrains, denotes the number of intermediate supports, prime denotes differentiation with respect to and dot differentiation with respect to time .
The piecewise continuous functions , , and denote the mass per unit length, the transverse bending stiffness per unit length, and the deck distributed damping constant. The constants and are the stiffness and damping constant of the viscoelastic support located at the th position , while denotes the ground acceleration.
The local form of the problem is obtained by integrating by parts (1) and can be formally written aswhere , , and denote, respectively, the mass, damping, and stiffness operator. They are expressed as ( is the Dirac’s delta function)
3. Eigenvalue Problem for PRSI Bridges
The freevibrations problem of the beam is obtained by posing in (2). The corresponding differential boundary problem is then reduced to an eigenvalue problem solvable by expressing the transverse displacement as the product of a spatial function and a timedependent function :After substituting (3) and (4) into (2) for , the following transcendental equation is obtained:Equation (5) is satisfied by an infinite number of eigenvalues and eigenvectors that occur in complex conjugate pairs [26]. The th eigenvalue contains information about the system vibration frequency and damping, while the th eigenvector is the th vibration shape. The solution of the eigenvalue problem for constant deck properties is briefly recalled in the appendix. The orthogonality conditions for the complex modes are [9]
4. Seismic Response of PRSI Bridges Based on CMS Method
4.1. Series Expansion of the Response
In the CMS method, the displacement of the beam is expanded as a series of the complex vibration modes aswhere = th complex modal shape and = th complex generalized coordinate. It is noteworthy that, in practical applications, the series is truncated at the term .
Substituting (8) into (1) written for one obtainsUpon substitution of (7) into (9) the following expression is obtained
4.2. Response to Impulsive Loading
For , the generalized coordinate assumes the form [22]with and .
After substituting into (10), one obtainsIt can be noted that the righthand side of (12) vanishes for by virtue of (6). Thus, the problem can be diagonalized and the th decoupled equation reads as follows:where , , and .
By assuming that the system is at rest for , the following equation holds at :Thus, since and , one can finally express asThe corresponding expression of the complex modal impulse response function is and the sum of the contribution to the complex modal impulse response function of the th mode and of its complex conjugate yields a real function , which may be expressed aswhere , , , and denotes the impulse response function of a SDOF system with natural frequency , damping ratio , and damped frequency , whose expression is
4.3. CMS Method for Seismic Response Assessment
By taking advantage of the derived closedform expression of the impulse response function, the seismic input is expressed as a sum of Delta Dirac functions as follows:The seismic displacement response is then expressed in terms of superposition of modal impulse responses:where and denote the response of the oscillator with natural frequency and damping ratio , subjected to the seismic input .
It is noteworthy that in the case of , corresponding to intermediate supports with no damping, the system becomes classically damped and one obtains , , , , , and , where = th is the real mode participation factor. Thus (19) reduces to the wellknown expression:The expression of the other quantities that can be of interest for the seismic performance assessment of PRSI bridges can be derived by differentiating (19). In particular, the abutment reactions are obtained as and , the th pier reaction is obtained as for , and the transverse bending moments are obtained as .
5. Simplified Methods for Seismic Response Assessment
In this paragraph, two simplified approaches often employed for the seismic analysis of structural systems are introduced and discussed. These approaches are both based on the assumed modes method [27, 28] and entail using in (8) a complete set of approximating realvalued functions (denoted as for ) instead of the exact complex vibration modes for describing the displacement field. The two analysis techniques differ for the approximating function employed. The use of these functions leads to a system of coupled Galerkin equations [27], and the generic th equation reads as follows:whereThe coupling between the generalized responses and may be due to nonzero values of the “nondiagonal” damping and stiffness terms and , for . In order to evaluate the extent of coupling due to damping, the index is introduced, whose definition is [13]This index assumes high values for intermediate supports with high dissipation capacity, while in the case of deck and supports with homogeneous properties it assumes lower and lower values for increasing number of supports, since the behaviour tends to that of a beam on continuous viscoelastic restraints [9].
A second coupling index is defined for the stiffness term, whose expression isIt is noteworthy that assumes high values in the case of intermediate supports with a relatively high stiffness compared to the deck stiffness. Conversely, it assumes lower and lower values for homogenous deck and support properties and increasing number of spans, since the behaviour tends to that of a beam on continuous elastic restraints, for which .
An approximation often introduced for practical purposes [13, 16] is to disregard the offdiagonal coupling terms to obtain a set of uncoupled equations. The resulting th equation describes the motion of a SDOF system with mass , stiffness , and damping constant . Thus, traditional analysis tools available for the seismic analysis of SDOF systems can be employed to compute efficiently the generalized displacements and the response , under the given ground motion excitation.
5.1. RMS Method for Seismic Response Assessment
This method, referred to as real modes superposition (RMS) method, is based on a series expansion of the displacement field in terms of the real modes of vibration for of the undamped (or classically damped) structure. The calculation of these classic modes of vibration involves solving an eigenvalue problem less computationally demanding than that required for computing the complex modes.
It is noteworthy that the real modes of vibration are retrieved from the PRSI bridge model of Figure 1 by neglecting the damping of the intermediate supports and that the orthogonality conditions for these modes can be obtained from (6) and (7) by posing and , for , leading toThe same expressions are obtained also in the case of mass proportional damping for the deck, that is, for where is a proportionality constant.
The orthogonality conditions of (25) result in nonzero values of and zero values of . These terms are discarded to obtain a diagonal problem.
5.2. FTS Method for Seismic Response Assessment
In the second method, referred to as Fourier terms superposition (FTS) method, the Fourier sineonly series terms for are employed to describe the motion. This method is employed in [8, 10] for the design of PRSI bridges and is characterized by a very reduced computational cost, since it avoids recourse to any eigenvalue analysis. It is noteworthy that the terms of this series correspond to the vibration modes of a simply supported beam, which coincides with the PRSI bridge model without any intermediate support. The orthogonality conditions for these terms, derived from (6) and (7) by posing, respectively, , for , areThe application of these orthogonality conditions leads to a set of coupled Galerkin equations due to the nonzero terms and for . These out of diagonal terms are discarded to obtain a diagonal problem.
6. Case Study
In this section, a realistic PRSI bridge [8] is considered as case study. First, the CMS method is applied to the bridge model by analyzing its modal properties and the related response to an impulsive input. Successively the seismic response is studied and a parametric analysis is carried out to investigate the influence of the intermediate support stiffness and damping on the seismic performance of the bridge components. Finally, the accuracy of the simplified analysis techniques is evaluated by comparing the results of the seismic analyses obtained with the CMS method and the RMS and FTS methods.
6.1. Case Study and Seismic Input Description
The PRSI bridge (Figure 2), whose properties are taken from [8], consists of a fourspan continuous steelconcrete superstructure (span lengths m and m, for a total length m) and of three isolated reinforced concrete piers. The value of the deck transverse stiffness is kNm^{2}, which is an average over the values of that slightly vary along the bridge length due to the variation of the web and flange thickness. The deck mass per unit length is equal to ton/m. The circular frequency corresponding to the first mode of vibration of the superstructure vibrating alone with no intermediate supports is rad/s. The deck damping constant is such that the first mode damping factor of the deck vibrating alone with no intermediate supports is equal to . The combined pier and isolator properties are described by Kelvin models, whose stiffness and damping constant are, respectively, kN/m and kN/m for the central support and kN/m and kN/m for the other supports. These support properties are the result of the design of isolation bearings ensuring the same shear demand at the base of the piers under the prefixed earthquake input.
The seismic input is described by a set of seven records compatible with the Eurocode 81 design spectrum, corresponding to a site with a peak ground acceleration (PGA) of where is the soil factor, assumed equal to 1.15 (ground type C), and is the gravity acceleration. They have been selected from the European strong motion database [29] and fulfil the requirements of Eurocode 8 [30]. Figure 3(a) shows the pseudoacceleration spectrum of the records, the mean spectrum, and the code spectrum, whereas Figure 3(b) shows the displacement response spectrum of the records and the corresponding mean spectrum.
(a)
(b)
6.2. Modal Properties and Impulsive Response
Table 1 reports the first 5 eigenvalues of the system, the corresponding vibration periods , and damping ratios . These modal properties are determined by solving the eigenvalue problem corresponding to (5). It is noteworthy that the even modes of vibration are characterized by an antisymmetric shape and a participating factor equal to zero, and thus they do not affect the seismic response of the considered configurations (uniform support excitation is assumed).

In general, the effect of the intermediate restraints on the dynamic behaviour of the PRSI bridge is to increase the vibration frequency and damping ratio with respect to the case of the deck vibrating alone without any restraint. In fact, the fundamental vibration frequency shifts from the value rad/s to rad/s, corresponding to a first mode vibration period of about 2.40 s. The first mode damping ratio increases from to . It should be observed that the value of is significantly lower than the value of the damping ratio of the internal and external intermediate viscoelastic supports at the same vibration frequency , equal to, respectively, 0.13 and 0.12. This is the result of the low dissipation capacity of the deck and of the dual load path behaviour of PRSI bridge, whose stiffness and damping capacity are the result of the contribution of both the deck and the intermediate supports [8, 9]. The damping ratio of the higher modes is in general very low and tends to decrease significantly with the increasing mode order.
Figure 4 shows the response of the midspan transverse displacement (Figure 4(a)) and of the abutment reactions (Figure 4(b)), for a unit impulse ground motion . The analytical exact expression of the modal impulse response is reported in (16). The different response functions plotted in Figure 4 are obtained by considering the contribution of the first mode only, the contribution of the first and third modes, and the contribution of the first, third, and fifth modes. Modes higher than the 5th have a negligible influence on the response. In fact, the mass participation factors of the 1st, 3rd, and 5th modes, evaluated through the RMS method, are 81%, 9.2%, and 3.25%. While the midspan displacement response is dominated by the first mode only, higher modes strongly affect the abutment reactions. The impulse response in terms of the transverse bending moments, not reported here due to space constraints, is also influenced by the higher modes contribution, but at a less extent than the abutment reactions.
(a)
(b)
6.3. Seismic Response
In this section, the characteristic aspects of the seismic performance of PRSI bridges are highlighted by starting from the seismic analysis of the case study and by evaluating successively the response variations due to changes of the most important characteristic parameters of the bridge, that is, the decktosupport stiffness ratio and the support damping. The bridge seismic performance is evaluated by monitoring the following response parameters: the midspan transverse displacement, the transverse bending moments, and the abutment shear. The CMS method is applied to compute the maximum overtime values of these response parameters by averaging the results obtained for the seven different records considered. Only the contribution of the first three symmetric vibration modes is considered due to the negligible influence of the other modes.
Figure 5 reports the mean of the envelopes of the transverse displacements and bending moments obtained for the different records. The displacement field coincides with a sinusoidal shape. This can be explained by observing that the isolated piers have a very low stiffness compared to the deck stiffness and, thus, their restraint action on the deck is negligible. Furthermore, the displacement response is in general dominated by the first mode of vibration. The transverse bending moments exhibit a different shape because they are significantly influenced by the higher modes of vibration. The support restraint action on the deck moments is almost negligible. It is noteworthy that differently from the case of fully isolated bridges, the bending moment demand in the deck may be very significant and attain the critical value corresponding to deck yielding, a condition that should be avoided according to EC82 [30].
(a)
(b)
In order to evaluate how the support stiffness and damping influence the seismic response of the bridge components, an extensive parametric study is carried out by considering a set of bridge models with the same superstructure, but with intermediate supports having different properties. Following [9], the intermediate supports properties can be synthetically described by the following nondimensional parameters:The parameter denotes the relative support to deck stiffness. Low values of correspond to a stiff deck relative to the springs while high values correspond to a more flexible deck relative to the springs. Limit case corresponds to the simply supported beam with no intermediate restraints. The parameter describes the energy dissipation of the supports and it is equal to the ratio between the energy dissipated by the dampers and the maximum strain energy in the springs, for a uniform transverse harmonic motion of the deck with frequency . These two parameters assume the values and for the bridge considered previously.
In the following parametric study, the values of these parameters are varied by keeping the same distribution of stiffness and of damping of the original bridge, that is, by using the same scale factor for all the values of and the same scaling factor for all the values of . In particular, the values of are assumed to vary in the range between 0 (no intermediate supports) and 2 (stiff supports relative to the superstructure), whereas the values of are assumed to vary in the range between 0 (intermediate supports with no dissipation capacity) and 0.3 (very high dissipation capacity).
Figure 6(a) shows the variation with of the periods of the first three modes of vibration that participate in the seismic response, that is, modes 1, 3, and 5, normalized with respect to the values observed for . It can be seen that only the fundamental vibration period is affected significantly by the support stiffness, whereas the higher modes of vibration exhibit a constant value of the vibration period for the different values. Figure 6(b) provides some information on the seismic response and it shows the influence of on the average maximum transverse displacements in correspondence with the outer intermediate supports and at the centre of the deck, denoted, respectively, as and . These response quantities do not vary significantly with . This occurs because the displacement demand is controlled by the first mode of vibration, whose period for the different values falls within the region corresponding to the flat displacement response spectrum.
(a)
(b)
Figure 7(a) shows the variation of the bending moment demand in correspondence with the piers, , and at the centre of the deck, . These response parameters follow a trend similar to that of the displacement demand; that is, they do not vary significantly with . Figure 7(b) shows the variation with of the sum of the pier reaction forces () and the sum of the abutment reactions (). In the same figure, the contribution of the first mode only to these quantities (, ) is also reported. It can be observed that the total base shear of the system increases for increasing values of , consistently with the increase of spectral acceleration at the fundamental vibration period. Furthermore, while the values of tend to increase for increasing values of , the values of remain almost constant. This can be explained by observing that the abutment reactions are significantly affected by the contribution of higher modes. Thus, in the case of intermediate support with no dissipation capacity, increasing the intermediate support stiffness does not reduce the shear forces transmitted to the abutments.
(a)
(b)
In the following, the joint influence of support stiffness and damping on the bridge seismic response is analysed. Figure 8(a) shows the variation of the transverse displacement versus , for different values of . The values of are normalized, for each value of , by dividing them by the value corresponding to . As expected, the displacement demand decreases by increasing the dissipation capacity of the supports. Furthermore, the differences between the displacement demands for the different values reduce for tending to zero.
(a)
(b)
Figure 8(b) shows the total pier () and abutment () shear demand versus , for different values of . These values are normalized, for a given value of , by dividing, respectively, by the values and corresponding to . It can be noted that both the pier and abutment reactions decrease by increasing . Furthermore, differently from the case corresponding to , the increase of the intermediate support stiffness results in a reduction of the abutment reactions. However, the maximum decrease of the values of the abutment reactions, achieved for high levels, is inferior to the maximum decrease of the values of pier reactions. This is a consequence of the combined effects of (a) the contribution of higher modes, which is significant only for the abutment reactions and negligible for the pier reactions, and (b) the reduced efficiency of the intermediate dampers in damping the higher modes of vibration [9].
Finally, Figure 9 shows the total base shear versus , for different values of . By increasing , the total base shear increases for low support damping values () while for very high values () it first decreases and then it slightly increases. Thus, for values of higher than 0.05, the total base shear is minimized when is equal to about 0.5. This result could be useful for the preliminary design of the isolator properties ensuring the minimization of the total base shear as performance objective.
6.4. Accuracy of Simplified Analysis Techniques
In this section, the accuracy of the simplified analysis techniques is evaluated by comparing the exact estimates of the response obtained by applying the CMS method with the corresponding estimates obtained by employing the RMS and FTS methods. In the application of all these methods, the first three symmetric terms of the set of approximating functions are considered and the ground motion records are those already employed in the previous section.
Figure 10 reports the results of the comparison for the case study described at the beginning of this section. It can be observed that the RMS and FTS method provide very close estimates of the mean transverse displacements (Figure 10(a)), whose shape is sinusoidal. This result was expected, since vibration modes higher than the first mode have a negligible influence on the displacements, as already pointed out in Figure 4. With reference to the transverse bending moment envelope (Figure 10(b)), the mean shape according to RMS and FTS agrees well with the exact shape, the highest difference being in correspondence with the intermediate restraints. The influence of the third mode of vibration on the bending moments shape is well described by both simplified techniques.
(a)
(b)
Table 2 reports and compares the peak values of the midspan displacement , the transverse abutment reaction , the midspan transverse moments , and the support reactions and according to the three analysis techniques, for the 7 records employed to describe the recordtorecord variability effects. These records are characterized by different characteristics in terms of duration and frequency content (Figure 3). Thus, the values of the response parameters of interest exhibit significant variation from record to record. The level of accuracy of the simplified analysis techniques is in general very high and it is only slightly influenced by the record variability. The highest relative error, observed for the estimates of the bending moment demand, is about 0.84% for the RMS method and 3.6% for the FTS method. The coupling indexes for the damping matrix are , , and for the RMS method and , , and for the FTS method. The coupling indexes for the stiffness matrix are , , and for the FTS method. Thus, the extent of coupling in both the stiffness and damping terms is very low and this explains why the approximate techniques provide accurate response estimates.

Table 3 reports the mean value and coefficient of variation of the monitored response parameters. With reference to this application, it can be concluded that both the simplified analysis techniques provide very accurate estimates of the response statistics despite the approximations involved in their derivation. The highest relative error between the results of the simplified analysis techniques and those of the CMS method is observed for the estimates of the mean bending moment demand and it is about 0.01% for the RMS method and 0.43% for the FTS method.

The accuracy of the approximate techniques is also evaluated by considering another set of characteristic parameters, corresponding to and , and representing intermediate supports with high stiffness and dissipation capacity.
The displacement shape of the bridge, reported in Figure 11(a), is sinusoidal, despite the high stiffness and damping capacity of the supports. However, the bending moment shape, reported in Figure 11(b), is significantly influenced by the higher modes of vibration and by the restrain action of the supports. The RMS method provides very accurate estimates of this shape, whereas the FTS method is not able to describe the support restraint effect.
(a)
(b)
Table 4 reports the peak values of the response parameters of interest for each of the 7 natural records considered, according to the different analysis techniques. The highest relative error, observed for the estimates of the bending moment demand, is less than 3.6% for the RMS method and 9.2% for the FTS method.

The mean and coefficient of variation of the response parameters of interest are reported in Table 5. The highest relative error, with respect to the CMS estimate, is observed for the mean bending moment demand and it is less than 0.79% for the RMS method and 1.12% for the FTS method. The coupling indexes for the damping matrix are , , and for the RMS method and , , and for the FTS method. The coupling indexes for the stiffness matrix are , , and . In general, these coefficients assume higher values compared to the values observed in the other bridge (Tables 2 and 4). In any case, the dispersion of the response due to the uncertainties in the values of the parameters describing the pier and the bearings might be more important than the error due to the use of simplified analysis methods.

7. Conclusions
This paper proposes a formulation and different analysis technique for the seismic assessment of partially restrained seismically isolated (PRSI) bridges. The formulation is based on a description of the bridges by means of a simply supported continuous beam resting on intermediate viscoelastic supports, whose properties are calibrated to represent the pier/bearing systems. The first analysis technique developed in this paper is based on the complex modes superposition (CMS) method and provides an exact solution to the seismic problem by accounting for nonclassical damping. This technique can be employed to evaluate the influence of the complex vibration modes on the response of the bridge components and also provides a reference solution against which simplified analysis approaches can be tested. The other two proposed techniques are based on the assumed modes method and on a simplified description of the damping properties. The real mode superposition (RMS) method employs a series expansion of the displacement field in terms of the classic modes of vibration, whereas the Fourier terms superposition (FTS) method employs a series expansion based on the Fourier sineonly series terms.
A realistic PRSI bridge model is considered and its seismic response under a set of ground motion records is analyzed by employing the CMS method. This method reveals that the bridge components are differently affected by higher modes of vibration. In particular, while the displacement demand is mainly influenced by the first mode, other parameters such as the bending moments and abutment reactions are significantly affected by the contribution of higher modes. A parametric study is then carried out to evaluate the influence of intermediate support stiffness and damping on the seismic response of the main components of the PRSI bridge. These support properties are controlled by two nondimensional parameters whose values are varied in the parametric analysis. Based on the results of the study, it can be concluded that in the case of intermediate support with low damping capacity, the abutment reactions do not decrease significantly by increasing the support stiffness; increasing the support damping results in a decrease in the abutment reactions which is less significant than the decrease in the pier reactions; for moderate to high dissipation capacity of the intermediate supports, there exists an optimal value of the damping that minimizes the total shear transmitted to the foundations.
Finally, the accuracy of the simplified analyses techniques is evaluated by comparing the results of the seismic analyses carried out by employing the CMS method with the corresponding results obtained by employing the two simplified techniques. Based on the results presented in this paper, it is concluded that, for the class of bridges analyzed, the extent of nonclassical damping is usually limited and satisfactory estimates of the seismic response are obtained by considering the simplified analysis approaches.
Appendix
The analytical solution of the eigenvalue problem is obtained under the assumption that , , and are constant along the beam length and equal, respectively, to , , and . The beam is divided into a set of segments, each bounded by two consecutive restraints (external or intermediate). The function describing the motion of the th segment of length is decomposed into the product of a spatial function and of a timedependent function . The function must satisfy the following equation:where .
The solution to (A.1) can be expressed aswith , , , and to be determined based on the boundary conditions at the external supports and the continuity conditions at the intermediate restraints. This involves the calculation of higher order derivatives up to the third order.
In total, a set of conditions is required to determine the vibration shape along the whole beam. At the first span the conditions apply while at the last span the support conditions are . The boundary conditions at the intermediate spring locations areBy substituting (A.2) into the boundary (supports and continuity) conditions, a system of homogeneous equations in the constants is obtained. Since the system is homogeneous, the determinant of coefficients must be equal to zero for the existence of a nontrivial solution. This procedure yields the following frequency equation in the unknown :In the general case of nonzero damping, the solution of the equation must be sought in the complex domain. It is noteworthy that, since the system is continuous, an infinite set of eigenvalues are obtained which satisfy (A.4). However, only selected values of are significant, because they correspond to the first vibration modes which are usually characterized by the highest participation factors. It is also noteworthy that the eigensolutions appear in complex pairs. Thus, if satisfies (A.4), then also its complex conjugate satisfies (A.4). The eigenvectors corresponding to these eigenvalues are complex conjugate, too. Finally, it is observed that the undamped circular vibration frequency and damping ratio corresponding to the th mode can be obtained by recalling that and . For zero damping, and the vibration modes are real.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The reported research was supported by the Italian Department of Civil Protection within the ReluisDPC Projects 2014. The authors gratefully acknowledge this financial support.
References
 M. T. A. Chaudhary, M. Abé, and Y. Fujino, “Performance evaluation of baseisolated Yamaagé bridge with high damping rubber bearings using recorded seismic data,” Engineering Structures, vol. 23, no. 8, pp. 902–910, 2001. View at: Publisher Site  Google Scholar
 G. C. Lee, Y. Kitane, and I. G. Buckle, “Literature review of the observed performance of seismically isolated bridges,” Research Progress and Accomplishments: 2000–2001, MCEER, State University of New York, Buffalo, NY, USA, 2001. View at: Google Scholar
 R. L. Boroschek, M. O. Moroni, and M. Sarrazin, “Dynamic characteristics of a long span seismic isolated bridge,” Engineering Structures, vol. 25, no. 12, pp. 1479–1490, 2003. View at: Publisher Site  Google Scholar
 J. Shen, M. H. Tsai, K. C. Chang, and G. C. Lee, “Performance of a seismically isolated bridge under nearfault earthquake ground motions,” Journal of Structural Engineering, vol. 130, no. 6, pp. 861–868, 2004. View at: Publisher Site  Google Scholar
 K. L. Ryan and W. Hu, “Effectiveness of partial isolation of bridges for improving column performance,” in Proceedings of the Structures Congress, pp. 1–10, 2009. View at: Publisher Site  Google Scholar
 M.H. Tsai, “Transverse earthquake response analysis of a seismically isolated regular bridge with partial restraint,” Engineering Structures, vol. 30, no. 2, pp. 393–403, 2008. View at: Publisher Site  Google Scholar
 N. Makris, G. Kampas, and D. Angelopoulou, “The eigenvalues of isolated bridges with transverse restraints at the end abutments,” Earthquake Engineering and Structural Dynamics, vol. 39, no. 8, pp. 869–886, 2010. View at: Publisher Site  Google Scholar
 E. Tubaldi and A. Dall'Asta, “A design method for seismically isolated bridges with abutment restraint,” Engineering Structures, vol. 33, no. 3, pp. 786–795, 2011. View at: Publisher Site  Google Scholar
 E. Tubaldi and A. Dall'Asta, “Transverse free vibrations of continuous bridges with abutment restraint,” Earthquake Engineering and Structural Dynamics, vol. 41, no. 9, pp. 1319–1340, 2012. View at: Publisher Site  Google Scholar
 G. Della Corte, R. de Risi, and L. di Sarno, “Approximate method for transverse response analysis of partially isolated bridges,” Journal of Bridge Engineering, vol. 18, no. 11, pp. 1121–1130, 2013. View at: Publisher Site  Google Scholar
 M. N. Fardis and G. Tsionis, “Eigenvalues and modes of distributedmass symmetric multispan bridges with restrained ends for seismic response analysis,” Engineering Structures, vol. 51, pp. 141–149, 2013. View at: Publisher Site  Google Scholar
 A. S. Veletsos and C. E. Ventura, “Modal analysis of nonclassically damped linear systems,” Earthquake Engineering & Structural Dynamics, vol. 14, no. 2, pp. 217–243, 1986. View at: Publisher Site  Google Scholar
 A. M. Claret and F. VenancioFilho, “Modal superposition pseudoforce method for dynamic analysis of structural systems with nonproportional damping,” Earthquake Engineering and Structural Dynamics, vol. 20, no. 4, pp. 303–315, 1991. View at: Publisher Site  Google Scholar
 F. VenancioFilho, Y. K. Wang, F. B. Lin, A. M. Claret, and W. G. Ferreira, “Dynamic analysis of nonproportional damping structural systems time and frequency domain methods,” in Proceedings of the 16th International Conference on Structural Mechanics in Reactor Technology (SMIRT '01), Washington, DC, USA, 2001. View at: Google Scholar
 S. Kim, “On the evaluation of coupling effect in nonclassically damped linear systems,” Journal of Mechanical Science and Technology, vol. 9, no. 3, pp. 336–343, 1995. View at: Publisher Site  Google Scholar
 G. Prater Jr. and R. Singh, “Quantification of the extent of nonproportional viscous damping in discrete vibratory systems,” Journal of Sound and Vibration, vol. 104, no. 1, pp. 109–125, 1986. View at: Publisher Site  Google Scholar
 M. Morzfeld, F. Ma, and N. Ajavakom, “On the decoupling approximation in damped linear systems,” Journal of Vibration and Control, vol. 14, no. 12, pp. 1869–1884, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 J. S. Hwang, K. C. Chang, and M. H. Tsai, “Composite damping ratio of seismically isolated regular bridges,” Engineering Structures, vol. 19, no. 1, pp. 55–62, 1997. View at: Publisher Site  Google Scholar
 S. C. Lee, M. Q. Feng, S.J. Kwon, and S.H. Hong, “Equivalent modal damping of shortspan bridges subjected to strong motion,” Journal of Bridge Engineering, vol. 16, no. 2, pp. 316–323, 2011. View at: Publisher Site  Google Scholar
 P. Franchin, G. Monti, and P. E. Pinto, “On the accuracy of simplified methods for the analysis of isolated bridges,” Earthquake Engineering and Structural Dynamics, vol. 30, no. 3, pp. 380–382, 2001. View at: Google Scholar
 K. A. Foss, “Coordinates which uncouple the equations of motion of damped linear dynamic systems,” Journal of Applied Mechanics, vol. 25, no. 3, pp. 361–364, 1958. View at: Google Scholar
 G. Oliveto, A. Santini, and E. Tripodi, “Complex modal analysis of a flexural vibrating beam with viscous end conditions,” Journal of Sound and Vibration, vol. 200, no. 3, pp. 327–345, 1997. View at: Publisher Site  Google Scholar
 M. Gürgöze and H. Erol, “Dynamic response of a viscously damped cantilever with a viscous end condition,” Journal of Sound and Vibration, vol. 298, no. 12, pp. 132–153, 2006. View at: Publisher Site  Google Scholar
 E. Tubaldi, “Dynamic behavior of adjacent buildings connected by linear viscous/viscoelastic dampers,” Structural Control and Health Monitoring, 2015. View at: Publisher Site  Google Scholar
 C. Truesdell and R. A. Toupin, “The classical field theories,” in Handbuch der Physik, Band III/l, Springer, Berlin, Germany, 1960. View at: Google Scholar  MathSciNet
 G. Prater Jr. and R. Singh, “Eigenproblem formulation, solution and interpretation for nonproportionally damped continuous beams,” Journal of Sound and Vibration, vol. 143, no. 1, pp. 125–142, 1990. View at: Publisher Site  Google Scholar  MathSciNet
 M. N. Hamdan and B. A. Jubran, “Free and forced vibrations of a restrained uniform beam carrying an intermediate lumped mass and a rotary inertia,” Journal of Sound and Vibration, vol. 150, no. 2, pp. 203–216, 1991. View at: Publisher Site  Google Scholar
 P. A. Hassanpour, E. Esmailzadeh, W. L. Cleghorn, and J. K. Mills, “Generalized orthogonality condition for beams with intermediate lumped masses subjected to axial force,” Journal of Vibration and Control, vol. 16, no. 5, pp. 665–683, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 N. Ambraseys, P. Smith, R. Bernardi, D. Rinaldis, F. Cotton, and C. BergeThierry, Dissemination of European StrongMotion Data, CDROM Collection, European Council, Environment and Climate Research Programme, 2000.
 European Committee for Standardization (ECS), “Eurocode 8—design of structures for earthquake resistance,” EN 1998, European Committee for Standardization, Brussels, Belgium, 2005. View at: Google Scholar
Copyright
Copyright © 2015 E. Tubaldi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.