Vibration Serviceability Assessment of a Historic Suspension Footbridge

: Experimental and numerical studies for the structural and vibration serviceability assessment of a historic suspension footbridge adopting non-invasive surveys and low-cost equipment are presented. Field surveys have been carried out to determine geometric properties, ambient vibration tests have been performed to estimate the dynamic properties, and the dynamic response of the footbridge under the action of a single crossing pedestrian has been recorded. Based on ﬁeld surveys, a 3D Finite Element model was built and was then calibrated against ambient vibration test results. The experimentally-measured maximum acceleration under the action of one crossing pedestrian is compared with the ones obtained numerically and analytically. Furthermore, vibration serviceability assessment under multi-pedestrian loading is carried out, adopting the simpliﬁed procedure recom-mended by a recent guideline. Results show that low-cost non-invasive dynamic testing is suitable to correctly identify the footbridge vertical natural frequencies and mode shapes, including higher-order ones, and to draw considerations about the state of degradation of the structure. Moreover, the level of vibration under the action of a single pedestrian can be estimated with sufﬁcient accuracy using a simpliﬁed loading model, provided that the modal damping ratio is properly tuned.


Introduction
During the 19th Century and the beginning of the 20th Century, a great number of suspended bridges were built in Europe, especially in Italy, Switzerland, and France [1]. The majority of these bridges were destroyed during the 20th Century and sometimes reconstructed, mainly because they were no longer adequate to support the new traffic loadings. Fortunately, some of them still exist in their original configuration, but very often they are closed to traffic, due to the uncertain conditions concerning their current structural properties and serviceability for the current loading scenarios. These surviving bridges are an important part of our architectural heritage, due to their elegance and construction technique. The preservation of these historic bridges implies the need to investigate their structural behavior to assess their safety and serviceability under the expected traffic conditions. Actually, the extreme lightness and flexibility of these kinds of structures, combined with very low structural damping, make them highly prone to human-induced vibrations [2]. Figure 1 shows some examples of suspended bridges in the north of Italy. The Morca bridge (Figure 1b) is the only one allowing for light vehicular traffic: it was retrofitted in 2003 and an extensive program of non-destructive tests and analytical investigations on its dynamic behavior was carried out [3]. The Ramello footbridge (Figure 2), built in 1954 and located in the countryside of La Spezia in Italy, is another example of a historic suspension its dynamic behavior was carried out [3]. The Ramello footbridge (Figure 2), built in 1954 and located in the countryside of La Spezia in Italy, is another example of a historic suspension footbridge. The footbridge was employed for pedestrian crossing and light vehicles until 2019 and then closed, like many Italian bridges that were not considered safe after the collapse of the Morandi Bridge in Genova. The aim of this study is to provide a structural and vibration serviceability assessment of the Ramello footbridge, using non-invasive surveys and low-cost equipment, that could assist public administration in the preservation of the structure. In order to avoid invasive interventions, the structural assessment is carried out in an indirect way by the modal testing of the footbridge and the comparison of the obtained modal properties with those derived from a Finite Element (FE) model built according to the nominal properties of the materials. its dynamic behavior was carried out [3]. The Ramello footbridge (Figure 2), built in 1954 and located in the countryside of La Spezia in Italy, is another example of a historic suspension footbridge. The footbridge was employed for pedestrian crossing and light vehicles until 2019 and then closed, like many Italian bridges that were not considered safe after the collapse of the Morandi Bridge in Genova. The aim of this study is to provide a structural and vibration serviceability assessment of the Ramello footbridge, using non-invasive surveys and low-cost equipment, that could assist public administration in the preservation of the structure. In order to avoid invasive interventions, the structural assessment is carried out in an indirect way by the modal testing of the footbridge and the comparison of the obtained modal properties with those derived from a Finite Element (FE) model built according to the nominal properties of the materials. The aim of this study is to provide a structural and vibration serviceability assessment of the Ramello footbridge, using non-invasive surveys and low-cost equipment, that could assist public administration in the preservation of the structure. In order to avoid invasive interventions, the structural assessment is carried out in an indirect way by the modal testing of the footbridge and the comparison of the obtained modal properties with those derived from a Finite Element (FE) model built according to the nominal properties of the materials.
Modal testing can be carried out based on controlled input that is measured and used in the identification process (Experimental Modal Analysis, EMA), or on ambient vibration tests where only the response is measured and the force is due to environmental excitation (Operational Modal Analysis, OMA) [4]. A brief review of modal testing methods for bridges can be found, e.g., in [5]. With the exception of modern non-contact methods developed mainly for laboratory tests (e.g., [6]), EMA usually involves the excitation of the structure through contact methods. It has been applied for the modal identification of footbridges using hammer or shaker excitation (e.g., [7,8]). However, OMA testing techniques have now become attractive, due to their relatively low cost, speed of implementation, and the recent improvements in recording equipment and computational methods (e.g., [9][10][11][12][13][14]). The low amplitude of vibrations in operational conditions requires very sensitive, low-noise sensors and a high-performance measurement chain [4]. OMA testing techniques are based on the assumption that the excitation is a stationary random process with approximately white noise characteristics. Possible extensions to non-stationary long-term vibration monitoring have been proposed (e.g., [15]). Modal parameters can then be extracted, adopting frequency or time-domain methods [4]: the most commonly adopted are Peak Picking (PP) (e.g., [3,9]), Frequency Domain Decomposition (FDD) ( [16], e.g., [3,13,14,17]), and Stochastic Subspace Identification (SSI) ( [18], e.g., [8,11,12,14]). The identification of the structural modal parameters can also be employed in the framework of damage detection, since variations in the structural physical properties reflect variations of the modal parameters [19,20]. In particular, mode shapes are more sensitive to damage than natural frequencies, and recent research is focused on the detection of damage using damping [19].
The experimental campaign on Ramello Footbridge includes a field survey of the footbridge geometry and element properties, ambient vibration tests, and live load tests of a single pedestrian crossing. The aim is to obtain the information necessary to build a reliable FE model, to get measurements for the dynamic identification of the footbridge, and to have a preliminary assessment of the level of vibration induced by pedestrians. Notwithstanding that the low amplitude of vibrations in operational conditions requires very sensitive, low-noise sensors, ambient vibrations are measured through low-cost equipment and the reliability of such measurement equipment is assessed. A comparison of the different techniques for modal parameters extraction is carried out in order to select the most appropriate one. Then, the FE model, built according to the field survey of the footbridge geometry, is validated against experimentally determined mode shapes and natural frequencies, and the correspondence of experimental and numerical modal properties is used as an indication of the structural health.
The paper develops through the following sections. Section 2 provides a detailed description of the footbridge geometry and element cross-sections estimated through field surveys. Section 3 describes the experimental campaign conducted in 2019 and 2021 and the subsequent identification of the footbridge dynamic properties by operational modal analysis. In Section 4, a 3D FE model of the footbridge is developed and calibrated based on the results obtained from ambient vibration tests. In Section 5, numerical simulations of the footbridge's dynamic response under single pedestrian loading are carried out and a comparison with experimental results is performed. Moreover, the vibration serviceability of the footbridge is investigated based on SÉTRA guidelines. Finally, conclusions are outlined in Section 6.

Description of the Footbridge
The Ramello suspension footbridge (Figures 2 and 3) was built in 1954 in the countryside of La Spezia in Italy. Due to the lack of drawing details, the geometric properties of the footbridge and its elements (Figures 4 and 5) were obtained from field surveys. The footbridge has a span of length L = 90 m and a width of 2.28 m. Two main suspended cables, with a sag of 7.05 m, connect two reinforced concrete pylons from one side of the river to the opposite side. Each suspension cable is made of three individual spiral strands, with a nominal diameter of 40 mm. The main cables are anchored into the ground at distances of 7.9 m and 7.2 m from the left and right pylons, respectively. The bridge deck is supported by I-shaped transversal steel beams with a step of one meter (Figure 3a), which are suspended to the main cables by means of 89 hangers. The latter are made of rolled steel and have a C-shaped cross-section. The transversal steel floor beams have I-shaped cross-sections and a length of 2.4 m. They support two longitudinal rolled I-shaped steel beams, located at 1.78 m of distance, and three longitudinal timber planks, above which transverse timber planks form the floor (Figure 3a). Both transverse and longitudinal planks have a square cross-section with a side of around 80 mm. tances of 7.9 m and 7.2 m from the left and right pylons, respectively. The bridge deck is supported by I-shaped transversal steel beams with a step of one meter (Figure 3a), which are suspended to the main cables by means of 89 hangers. The latter are made of rolled steel and have a C-shaped cross-section. The transversal steel floor beams have I-shaped cross-sections and a length of 2.4 m. They support two longitudinal rolled I-shaped steel beams, located at 1.78 m of distance, and three longitudinal timber planks, above which transverse timber planks form the floor (Figure 3a). Both transverse and longitudinal planks have a square cross-section with a side of around 80 mm.
Hangers are connected to floor beams by bolted connections (Figure 3b) and to main cables by means of steel sockets (Figure 3c), which keep the three strands in contact with each other, avoiding any relative displacement. Welded connections link the bottom flanges of the longitudinal beams to the top flanges of the transversal floor beams. Handrails are made of steel circular hollow sections ( Figure 2b) and linked to the hangers with joints that restrain the vertical relative translational movement between the handrails and the hangers while enabling horizontal movement. Surface corrosion can be observed widely over the main structural elements (Figure 3b,c).
It is worth noting that the tension force of cables was not measured during the field surveys; therefore, it is considered an unknown parameter that is estimated based on the Finite Element (FE) model of the footbridge that will be discussed in the next sections.

Experimental Campaigns and Operational Modal Analysis
This section describes the experimental campaigns carried out in 2019 and 2021 with the aim of identifying the footbridge's dynamic properties and measuring the footbridge's dynamic response under the dynamic excitation of a single pedestrian.

Ambient Vibration Testing
Four ambient vibration tests with different measurement setup arrangements were conducted on the footbridge to detect its dynamic characteristics. Figure 6a shows the positions of accelerometers in the 4 setups. Suspended footbridges are characterized by extremely low natural frequencies, so the modes at risk of human-induced vibrations are usually higher modes, whose shapes are difficult to be correctly identified unless a great number of sensors are installed. Setup 1 was implemented in 2019, while the other three measurement setups were adopted in 2021 to identify more accurately the mode shapes of the footbridge. The low-cost equipment used for the tests included a 14-channel Labjack U6 data acquisition system with 10 TE Connectivity 4030 signal conditioned MEMS DC triaxial accelerometers that are able to measure accelerations in ±6 g range with a sensitivity of 333 mV/g and a nominal 0-200 Hz bandwidth. Figure 6b shows the detail of the accelerometers' mounting: they were bolted to an aluminum profile, which was rigidly connected through tight bands at the bottom of the hanger as close as possible to the joint with the transversal beam.
The dynamic response of the footbridge was measured at 10 different locations in each measurement setup: vertical acceleration was measured at all the locations, while transverse accelerations were only at four locations (locations 1, 3, 7, and 9). In each measurement setup, the ambient vibrations were simultaneously recorded for one hour with a sampling frequency of 500 Hz. It is worth noting that the tension force of cables was not measured during the field surveys; therefore, it is considered an unknown parameter that is estimated based on the Finite Element (FE) model of the footbridge that will be discussed in the next sections.

Experimental Campaigns and Operational Modal Analysis
This section describes the experimental campaigns carried out in 2019 and 2021 with the aim of identifying the footbridge's dynamic properties and measuring the footbridge's dynamic response under the dynamic excitation of a single pedestrian.

Ambient Vibration Testing
Four ambient vibration tests with different measurement setup arrangements were conducted on the footbridge to detect its dynamic characteristics. Figure 6a shows the positions of accelerometers in the 4 setups. Suspended footbridges are characterized by extremely low natural frequencies, so the modes at risk of human-induced vibrations are usually higher modes, whose shapes are difficult to be correctly identified unless a great number of sensors are installed. Setup 1 was implemented in 2019, while the other three measurement setups were adopted in 2021 to identify more accurately the mode shapes of the footbridge. The low-cost equipment used for the tests included a 14-channel Labjack U6 data acquisition system with 10 TE Connectivity 4030 signal conditioned MEMS DC triaxial accelerometers that are able to measure accelerations in ±6 g range with a sensitivity of 333 mV/g and a nominal 0-200 Hz bandwidth. Figure 6b shows the detail of the accelerometers' mounting: they were bolted to an aluminum profile, which was rigidly connected through tight bands at the bottom of the hanger as close as possible to the joint with the transversal beam.

Modal Identification of the Footbridge
The data processing and modal identification were performed by using MACEC, a MATLAB-based toolbox for the experimental modal analysis of structures [29]. The identification of the modal parameters was carried out by adopting the SSI, PP, and FDD techniques. Figure 7, referred to as Setup 1, plots the stabilization diagram derived from the SSI technique (a), the averaged normalized power spectral density function for the identification based on the PP technique (b), and the singular values for FDD (c). In Figure 7a, the power spectral density functions of the measured accelerations are superimposed on the stabilization diagrams in order to verify that the identified poles correspond to the peaks in the acceleration spectrum. It can be deduced that most of the stable poles in Figure 7a correspond to peaks in the power spectral density function, also detected in Figure 7b,c. The dynamic response of the footbridge was measured at 10 different locations in each measurement setup: vertical acceleration was measured at all the locations, while transverse accelerations were only at four locations (locations 1, 3, 7, and 9). In each measurement setup, the ambient vibrations were simultaneously recorded for one hour with a sampling frequency of 500 Hz.

Modal Identification of the Footbridge
The data processing and modal identification were performed by using MACEC, a MATLAB-based toolbox for the experimental modal analysis of structures [29]. The identification of the modal parameters was carried out by adopting the SSI, PP, and FDD techniques. Figure 7, referred to as Setup 1, plots the stabilization diagram derived from the SSI technique (a), the averaged normalized power spectral density function for the identification based on the PP technique (b), and the singular values for FDD (c). In Figure 7a, the power spectral density functions of the measured accelerations are superimposed on the stabilization diagrams in order to verify that the identified poles correspond to the peaks in the acceleration spectrum. It can be deduced that most of the stable poles in Figure 7a correspond to peaks in the power spectral density function, also detected in Figure 7b,c.  Table 1 reports the natural frequencies identified through the three techniques, together with a synthetic description of the corresponding mode shape (L = Lateral, V = Vertical, T = Torsional; S = Symmetric, A = Asymmetric). It can be deduced that the natural frequencies identified with the PP and FDD techniques are coincident and they are generally in very good agreement with the ones identified by SSI. However, some of the vibration modes are identified only through the SSI technique. This circumstance is due to the weak ambient excitation and to the significant noise measured by the accelerometers, which does not allow for the clear identification of the peaks corresponding to some of the vibration modes through frequency-domain methods. The results in Table 1 show that the SSI technique is the one that allows the identification of the greater number of modes. Furthermore, frequency-domain methods require long time histories for the reliable identification of the modal damping ratios. Thus, in the following, results obtained by adopting SSI for all the tested setups are reported. Figure 8 shows the stabilization diagrams extracted from Setups 2-4.   Table 1 reports the natural frequencies identified through the three techniques, together with a synthetic description of the corresponding mode shape (L = Lateral, V = Vertical, T = Torsional; S = Symmetric, A = Asymmetric). It can be deduced that the natural frequencies identified with the PP and FDD techniques are coincident and they are generally in very good agreement with the ones identified by SSI. However, some of the vibration modes are identified only through the SSI technique. This circumstance is due to the weak ambient excitation and to the significant noise measured by the accelerometers, which does not allow for the clear identification of the peaks corresponding to some of the vibration modes through frequency-domain methods. The results in Table 1 show that the SSI technique is the one that allows the identification of the greater number of modes. Furthermore, frequency-domain methods require long time histories for the reliable identification of the modal damping ratios. Thus, in the following, results obtained by adopting SSI for all the tested setups are reported.  The identified dynamic characteristics of the footbridge are summarized in Table 2. The natural frequencies and damping ratios are obtained by combining and averaging the results extracted from different measurement setups. The detected modal shapes consist of one lateral bending mode, two torsional modes, and ten vertical bending modes that have symmetric or asymmetric shapes. The results show that the identified modes with frequencies greater than 1.3 Hz are quite stable and reliable since those modes were found The identified dynamic characteristics of the footbridge are summarized in Table 2. The natural frequencies and damping ratios are obtained by combining and averaging the results extracted from different measurement setups. The detected modal shapes Buildings 2022, 12, 732 9 of 21 consist of one lateral bending mode, two torsional modes, and ten vertical bending modes that have symmetric or asymmetric shapes. The results show that the identified modes with frequencies greater than 1.3 Hz are quite stable and reliable since those modes were found in most of the ambient vibration measurements. There are four vertical bending modes (modes 9-12) whose natural frequencies fall within the frequency range of dynamic loading induced by the first walking harmonic [22]. Therefore, the resonant condition could occur for the footbridge due to walking pedestrians. The identified damping ratios of the modes at resonance risk for each Setup are reported in Table 3. It is worth noting that the damping ratios identified from Setups 2 to 4 are generally higher than those from Setup 1 and quite high with respect to the characteristic damping ratios of steel footbridges (around 0.5% [30]). In the Morca suspended footbridge [3], the identified damping ratios were quite high (2.73-7.69%) as well, and this unusual result was explained as due to energy dissipation in the connections between structural elements of the deck. Very large coefficients of variation of bridge damping ratios obtained from different experiments were also reported in [31], where it was observed that operating conditions affect the modal properties and, in particular, damping ratios are very sensitive to test and analysis methods. Therefore, damping average values are not considered fully reliable and will be verified and tuned within the simulation of the footbridge dynamic response (Section 5).

Measurements of Dynamic Response Due to a Single Pedestrian Crossing
In the experimental campaign in 2021, in addition to ambient vibration tests, three single walking load tests were carried out on the footbridge. Based on the results obtained from Setup 1 in 2019, a pedestrian, with a weight of 780 N, crossed the footbridge with his step frequency synchronized by a metronome to the frequencies of 1.5, 1.75, and 2.05 Hz, respectively. These frequencies do not exactly match the natural frequencies of the footbridge, which were correctly identified after the 2021 survey, but are very close to the resonant condition.
The layout of accelerometers in Setup 4 was used to record the acceleration of the footbridge deck in the lateral (L) and vertical (V) directions. As an example, the vertical acceleration response .. q(t) and its corresponding Fourier spectrum .. q( f ) measured by accelerometer 2 (V2) during the passage of the pedestrian with a step frequency of 1.75 Hz are shown in Figure 9. The Fourier spectrum for this case shows that the dynamic response of the footbridge is dominated by the mode with the natural frequency of 1.803 Hz. Moreover, to exclude the effects of noises on the peak accelerations, the responses were lowpass filtered and peak accelerations were extracted from the filtered signal.  Table 4 reports the maximum acceleration responses obtained from these single-pedestrian walking load tests. The maximum transverse and vertical acceleration responses of 0.156 and 0.830 m s ⁄ , respectively, are associated with a step frequency of 2.05 Hz.
The vibration serviceability of the footbridge under pedestrian walking loads can be evaluated based on the SÉTRA guidelines [28]. The guideline classifies footbridges according to four comfort levels based on maximum acceleration responses in both horizontal and vertical directions ( Table 5). The vertical and lateral peak accelerations measured on the footbridge both fall within the mean comfort level. This result, obtained for a single pedestrian crossing, raises concerns about the comfort level due to the crossing of multiple pedestrians, which will be investigated in Section 5.  Table 4 reports the maximum acceleration responses obtained from these singlepedestrian walking load tests. The maximum transverse and vertical acceleration responses of 0.156 and 0.830 m/s 2 , respectively, are associated with a step frequency of 2.05 Hz.

Walking Tests
Step Frequency Accelerometer The vibration serviceability of the footbridge under pedestrian walking loads can be evaluated based on the SÉTRA guidelines [28]. The guideline classifies footbridges according to four comfort levels based on maximum acceleration responses in both horizontal and vertical directions ( Table 5). The vertical and lateral peak accelerations measured on the footbridge both fall within the mean comfort level. This result, obtained for a single pedestrian crossing, raises concerns about the comfort level due to the crossing of multiple pedestrians, which will be investigated in Section 5.

Finite Element Model of the Footbridge
The FE model of the footbridge was built with ANSYS software [32] to investigate numerically the dynamic behavior of the structure. The model was built based on the field surveys and then updated according to ambient vibration test results.
The main cables and the hangers were modeled using the 3D spar element "LINK180". The effective steel area, density, modulus of elasticity, and Poisson's ratio of each spiral strand of the main cables were taken as 945 mm 2 , 7850 kg/m 3 , 160 GPa, and 0.3, respectively. Furthermore, the longitudinal and transversal floor beams were modeled using the 3D elastic beam element "BEAM188" considering the modulus of elasticity and density as 210 GPa and 7850 kg/m 3 , respectively. The same material properties were assigned to the hangers. The timber boards and handrails were assumed as nonstructural elements and modeled with the concentrated mass element "MASS21". The amount of concentrated mass applied at the deck nodes is estimated by assuming the values of density for timber boards and handrails as 700 and 7850 kg/m 3 , respectively. Table 6 summarizes the properties of all elements employed to build the numerical model. Pylons were not modeled due to the high stiffness assumed at the pylon saddles: hence, the cables were restrained at the pylon position by means of rigid constraints that allow sliding in the longitudinal direction. The ends of the cables were restrained to the ground by pinned supports. Moreover, it was assumed that the hangers were hinged to the main cables and floor beams. Longitudinal and transversal beams were connected with fixed joints. Finally, both ends of the longitudinal beams were restrained to translations and rotations. The general view of the FE model of the footbridge is shown in Figure 10.  Figure  10.

Preliminary Static Analysis
Suspended footbridges are characterized by non-linear behavior, due to the wellknown geometric effects of the main cables [33]. Therefore, in their structural analysis, it is first necessary to determine the geometric configuration resulting from dead load and cable prestress. In the case of the Ramello footbridge, it is worth recalling that the surveyed geometry refers to the deformed configuration of the footbridge under dead load and prestress.
Since tension force in the cables has not been measured, it has been determined through a parametric study. Non-linear static analyses have been carried out for different values of the cable pretension T in the range 0-360 kN and vertical deflections q have been measured in different sections of the footbridge. Figure 11 plots the vertical deflections at the abscissas x = L/2 and x = L/3 along the deck for different values of the cable prestress. The value of prestress corresponding to the minimum deflection at both monitored positions has been retained for successive dynamic analyses. Specifically, a tension force in the main cables T = 180 kN was selected.

Preliminary Static Analysis
Suspended footbridges are characterized by non-linear behavior, due to the wellknown geometric effects of the main cables [33]. Therefore, in their structural analysis, it is first necessary to determine the geometric configuration resulting from dead load and cable prestress. In the case of the Ramello footbridge, it is worth recalling that the surveyed geometry refers to the deformed configuration of the footbridge under dead load and prestress.
Since tension force in the cables has not been measured, it has been determined through a parametric study. Non-linear static analyses have been carried out for different values of the cable pretension T in the range 0-360 kN and vertical deflections q have been measured in different sections of the footbridge. Figure 11 plots the vertical deflections at the abscissas x = L/2 and x = L/3 along the deck for different values of the cable prestress. The value of prestress corresponding to the minimum deflection at both monitored positions has been retained for successive dynamic analyses. Specifically, a tension force in the main cables T = 180 kN was selected. through a parametric study. Non-linear static analyses have been carried out for different values of the cable pretension T in the range 0-360 kN and vertical deflections q have been measured in different sections of the footbridge. Figure 11 plots the vertical deflections at the abscissas x = L/2 and x = L/3 along the deck for different values of the cable prestress. The value of prestress corresponding to the minimum deflection at both monitored positions has been retained for successive dynamic analyses. Specifically, a tension force in the main cables T = 180 kN was selected.

Modal Analysis
A modal analysis was performed on the footbridge model to extract its dynamic properties. The modal analysis was conducted after a nonlinear static analysis of the footbridge, subjected to dead loads and the pretension of the cables, in order to determine the geometric tangent stiffness matrix [2,34]. The natural frequencies of the first thirty global modes are summarized in Table 7. From a direct inspection of Tables 2 and 7, it is evident that ambient vibration tests did not allow the identification of the complete set of the lateral and torsional modes, but only of a very limited number of them. This is due to the evanescent excitation of such modes, provided only by the wind, which was very weak during the ambient vibration tests. For this reason, the validation of the numerical model, based on a comparison of modal parameters (i.e., natural frequencies and mode shapes) with the experimental estimates obtained from the SSI method, is limited to the vertical bending modes. The correlation of numerical and experimental mode shapes is investigated through the Modal Assurance Criterion (MAC), which is expressed as [35]: where φ E,i is a vector that represents the i-th mode shape extracted experimentally from field vibration tests, φ N,j is a vector that represents the j-th numerical mode shape, and T stands for transpose. Generally, MAC values greater than 0.8 mean a very good match between the two mode shapes. Furthermore, the correlation of the two modes in terms of natural frequency can be investigated by calculating the percentage frequency error based on the following expression: where f E and f N are the experimental and numerical natural frequencies, respectively. According to Table 8, there is generally a good agreement between the numerically and experimentally identified modal parameters, with MAC values higher than 0.9 for almost all the considered modes, and frequency errors generally lower than 5%. This outcome demonstrates that the number of sensors and setups was sufficient to correctly detect the vertical bending modes, including higher-order ones. Moreover, the very good match between experimental and numerical mode shapes allows us to hypothesize that the footbridge is not interested in localized damage, despite the diffused surface deterioration of the structural elements. Figure 12 plots the mode shapes of the four vertical bending modes, whose frequencies fall within the range of walking excitation, while Figure 13 compares the experimental and numerical mode shapes reported in Table 8.    Figure 13. Comparison between experimental and numerical mode shapes. Figure 13. Comparison between experimental and numerical mode shapes.

Vibration Serviceability Assessment
The dynamic response of the footbridge under single and multiple pedestrian loads is calculated in order to assess its vibration serviceability. It is assumed that the footbridge system is a linear mono-dimensional system, whose dynamics are described by the equation of motion: where q(x,t) is the displacement of the footbridge, x is the abscissa along the bridge deck and t is the time, m(x) is the structural mass per unit length, C is the damping operator, L is the stiffness operator, f p (x,t) is the external force per unit length. Under the hypothesis of classical damping, Equation (3) is usually solved by applying the principal transformation and assuming that the dynamic response is dominated by one mode of vibration: where ϕ j is the j-th mode shape of the footbridge and p j is the corresponding principal coordinate. The equation of motion of the j-th principal coordinate p j is expressed as: ..
where M j , ω j , ξ j , and F j (t) are the modal mass, circular natural frequency, modal damping ratio, and modal force of the j-th mode, respectively, and L is the span length.

Single Pedestrian Excitation
The mathematical model usually adopted to describe the vertical dynamic load due to a single pedestrian walking is defined by the following expression ( [36,37]): where c is the pedestrian speed, δ is Dirac delta function, and F p (t) is the time-varying vertical force induced by a single pedestrian, which is usually defined as a sum of Fourier harmonic components: where G is the static weight of the pedestrian, h is the order number of the harmonic, H is the total number of contributing harmonics, α h is the dynamic load factor (DLF) of the h-th harmonic, f s the step frequency (Hz), and φ h is the phase angle of the h-th harmonic. It is generally accepted and confirmed by experimental tests that the dynamic response of footbridges is mainly affected by the first walking harmonic. For normal walking speeds, the DLF of the first vertical harmonic is α 1 ∼ = 0.4 and pedestrian velocity can be calculated as c ∼ = 0.9 f s (e.g., [30]). The dynamic response of the footbridge due to a single pedestrian crossing can be obtained by numerically solving the equation of motion (5), with f p (x,t) expressed as in Equations (6) and (7).
As an alternative, the dynamic response of the footbridge can be approximately calculated based on the analytical solution proposed by Piccardo and Tubino [38], which predicts the dynamic response of the footbridge due to a single pedestrian crossing in resonant conditions. Specifically, the j-th principal coordinate p j is calculated in the nondimensional form, as follows: where the nondimensional parameters t and Ω c are defined as: In order to predict the dynamic response of the footbridge due to a single pedestrian crossing in the experimental tests described in Section 3.3, the dynamic response of the footbridge to a single pedestrian walking at 1.5, 1.75, and 2.05 Hz is estimated numerically (Equation (5)) and analytically (Equation (8)) considering as the mode of interest j the one whose frequency is the nearest to the excitation frequency, i.e., experimental modes 9, 10, and 11, respectively. In both cases, the adopted damping ratio is the average damping identified from the ambient vibration tests reported in Table 2. The obtained peak accelerations are reported in Table 9, compared with experimental results. Both numerical and analytical predictions are not very accurate, especially for step frequencies of 1.5 and 2.05 Hz, with a maximum error of about 35%. In order to obtain a more accurate prediction of the peak acceleration response, the damping ratios were modified within the range of values identified in the different setups (Table 3) in order to obtain a quite perfect match between numerical and experimental maximum accelerations. The modified damping ratios for the considered three bending modes are 2.7%, 0.78%, and 1.8%, respectively. Figure 14 plots an example of the vertical acceleration responses due to a single pedestrian crossing with a step frequency of 1.75 Hz obtained through numerical and analytical approaches with average and modified damping. From a comparison between Figure 14a,b, it can be deduced that assuming a modified damping ratio (0.78%) lower than the average value (1.732%, see Table 3), both the analytical and the numerical vertical responses increase. However, the increase in the numerical response (f s = 1.75 Hz) is larger than the increase in the analytical one (f s = 1.803 Hz). This circumstance is due to the fact that the modal force associated with the moving harmonic load (Equations (5)-(7)) can be decomposed into two harmonic components, one of which is closer to the resonance condition when f s = 1.75 Hz. Furthermore, comparing Figures 14b and 9, it can be deduced that the experimentally measured time histories and the ones obtained numerically and analytically are not in perfect agreement. The difference can be due to many factors, such as the slight variations in the walking speed and step frequency in experimental tests, that are not taken into account by the analytical force model in Equations (6) and (7). Table 9 reports the peak accelerations obtained with the numerical and analytical approach by adopting the modified damping ratios. The results show that the modified damping generally allows us to obtain a more accurate analytical prediction of the peak response with a maximum error of about 16%. The obtained results confirm that, for the present footbridge, the dynamic response to a single pedestrian is mainly dominated by a single mode and that the analytical model in Equation (8) is able to predict the peak response with sufficient accuracy, despite the fact that resonant conditions were not perfectly achieved during the experimental tests. (9) In order to predict the dynamic response of the footbridge due to a single pedestrian crossing in the experimental tests described in Section 3.3, the dynamic response of the footbridge to a single pedestrian walking at 1.5, 1.75, and 2.05 Hz is estimated numerically (Equation (5)) and analytically (Equation (8)) considering as the mode of interest j the one whose frequency is the nearest to the excitation frequency, i.e., experimental modes 9, 10, and 11, respectively. In both cases, the adopted damping ratio is the average damping identified from the ambient vibration tests reported in Table 2. The obtained peak accelerations are reported in Table 9, compared with experimental results. Both numerical and analytical predictions are not very accurate, especially for step frequencies of 1.5 and 2.05 Hz, with a maximum error of about 35%. In order to obtain a more accurate prediction of the peak acceleration response, the damping ratios were modified within the range of values identified in the different setups (Table 3) in order to obtain a quite perfect match between numerical and experimental maximum accelerations. The modified damping ratios for the considered three bending modes are 2.7%, 0.78%, and 1.8%, respectively. Figure 14 plots an example of the vertical acceleration responses due to a single pedestrian crossing with a step frequency of 1.75 Hz obtained through numerical and analytical approaches with average and modified damping. From a comparison between Figure 14a and Figure 14b, it can be deduced that assuming a modified damping ratio (0.78%) lower than the average value (1.732%, see Table 3), both the analytical and the numerical vertical responses increase. However, the increase in the numerical response (fs = 1.75 Hz) is larger than the increase in the analytical one (fs = 1.803 Hz). This circumstance is due to the fact that the modal force associated with the moving harmonic load (Equations (5)- (7)) can be decomposed into two harmonic components, one of which is closer to the resonance condition when fs = 1.75 Hz. Furthermore, comparing Figures 14b and 9, it can be deduced that the experimentally measured time histories and the ones obtained numerically and analytically are not in perfect agreement. The difference can be due to many factors, such as the slight variations in the walking speed and step frequency in experimental tests, that are not taken into account by the analytical force model in Equations (6) and (7). Table 9 reports the peak accelerations obtained with the numerical and analytical approach by adopting the modified damping ratios. The results show that the modified damping generally allows us to obtain a more accurate analytical prediction of the peak response with a maximum error of about 16%. The obtained results confirm that, for the present footbridge, the dynamic response to a single pedestrian is mainly dominated by a single mode and that the analytical model in Equation (8) is able to predict the peak response with sufficient accuracy, despite the fact that resonant conditions were not perfectly achieved during the experimental tests.

Multiple Pedestrian Excitation
The serviceability of the footbridge is assessed based on the approach proposed by the SÉTRA guidelines [28]. According to SÉTRA, footbridges are classified into four classes, from urban footbridges with heavy traffic (Class I) to seldom-used footbridges (Class IV). Despite the Ramello footbridge belonging to Class IV, for which dynamic calculations are not required, the guideline suggests considering at least Class III for extremely lively footbridges to ensure a minimum amount of risk control. For Class III footbridges, serviceability assessment should be performed under the action of a sparse crowd, characterized by a pedestrian density ρ equal to 0.5 ped/m 2 . A resonant uniformly distributed harmonic load F v (t) [N/m 2 ] is defined as follows: F v (t) = αG cos 2π f j t N eq ψ, where αG = 280 N, N eq is the equivalent number of perfectly synchronized pedestrians per square meter generating the 95th percentile of the peak acceleration response induced by random pedestrians, f j is the natural frequency of the j-th mode, B is the deck width, and ψ is a reduction factor to consider that the risk of resonance reduces if the footbridge frequency is outside the interval of 1.7-2.1 Hz for vertical vibrations. This load should be applied for each vertical mode at risk with the same sign as the one of the considered mode shape to obtain the most unfavorable effect. Moreover, the modal mass should be estimated while also taking into account the mass of pedestrians. The peak acceleration of the footbridge can be predicted with the following expression: ..
where M j,tot is the total modal mass of the footbridge and pedestrians. The peak accelerations calculated for the three modes at resonance risk are reported in Table 10. According to Table 5, all the obtained values fall in the range of unacceptable comfort. It is worth pointing out that the simplified procedure proposed by the SÉTRA guidelines often leads to an overestimation of the structural response in the vertical direction since human-structure interaction is not taken into account. However, the very high values of peak accelerations suggest the need to further investigate the footbridge dynamic behavior under more realistic loading scenarios and to evaluate the possibility of installing suitable countermeasures.

Conclusions
This paper investigated the dynamic behavior and vibration serviceability of a historic suspension footbridge based on non-invasive low-cost modal testing and a numerical model.
The comparison among different frequency-and time-domain techniques for modal parameters extraction has shown that when using low-cost sensors with low vibration levels, the time-domain SSI method allows for the extraction of a larger number of modes. In particular, the operational modal analysis allowed the identification of 11 vibration modes, including one lateral, two torsional, and eight bending mode shapes within the frequency range of 0-2.6 Hz. The comparison of the modal characteristics of the numerical model and the ones estimated experimentally showed that ambient vibration tests carried out with low-cost sensors allowed a reliable identification of the bending modes, but not of the torsional and lateral ones, which were very weakly excited by ambient actions. The accordance between the experimentally identified natural frequencies and the ones obtained from the numerical model demonstrates that the stiffness of the elements corresponds to the one evaluated assuming standard values of the elastic modulus of steel, excluding a significant degradation of material properties. Furthermore, the good accordance between the experimental and numerical mode shapes confirms that the global structural behavior of the footbridge is well-captured by the numerical model and excludes significant local damages to the structural elements that would reflect on the identified mode shapes.
The footbridge is characterized by four lowly damped bending vibration modes in the range of step frequency typical of normal walking, and thus it is very sensitive to human-induced vibrations. The comparison between experimental accelerations and the ones estimated numerically based on a moving harmonic load model of a single pedestrian confirmed the reliability of such an approximated loading model. Finally, the serviceability assessment according to the SÉTRA guidelines showed that the footbridge would have an unacceptable comfort level under the crossing of a sparse crowd. The very high values of peak accelerations suggest the need to further investigate the footbridge dynamic behavior under more realistic loading scenarios and to evaluate the possibility of installing suitable vibration countermeasures (e.g., [39][40][41]).
In summary, the presented results have evidenced the following advantages (+) and drawbacks (−) of the non-invasive low-cost technology adopted: + The reliable identification of natural frequencies and mode shapes, in conjunction with the SSI modal identification technique; + The potential to draw considerations on the structural damage on the basis of the comparison between experimental and numerical modal properties; + The accurate estimation of the acceleration level for vibration serviceability assessment; − Modal identification is limited to vertical bending modes; − Modal damping ratios identified from ambient vibrations are very dispersed. The last two issues are worthy of further investigation and the possibility to perform forced vibration tests or adopt higher-level equipment should be taken into consideration in order to also identify lateral and torsional modes and better estimate damping ratios.