**TECHNICAL PAPER**

**Dynamic behaviour of underspanned suspension road bridges under traffic loads**

**J Oliva; J M Goicolea; P Antolín M Astiz**

**ABSTRACT**

*{PRPgenerator)*to generate these profiles is presented in this paper. The structural dynamic sensitivity of underspanned suspension bridges was found to be considerable, as well as the dynamic amplification factors and deck accelerations. It was also found that vehicle speed has a relevant influence on the results. In addition, the impact of bridge deformation on vehicle vibration was addressed, and the effect on the comfort of vehicle users was shown to be negligible.

**Keywords:** underspanned suspension bridges, road roughness, vehicle– bridge dynamic interaction, deck vibrations, penalty method

**INTRODUCTION**

**Structural system**

In underspanned suspension bridges the ten-sioned cable system is located below the deck, and tendons are deflected by struts or braces that are connected with the girder and thus transmit the upward cable deviation forces that help to support it. Cables are anchored in the deck over abutments and piers, if piers exist. This system is more efficient in single-span or simply-supported structures; in continuous bridges the under-deck cable system contribution in the negative bending regions (over piers) is very small, and tendons should be placed above the girder in those parts (Ruiz-Terán & Aparicio 2007b). This paper deals only with simply-supported bridges.

Cable forces are transferred into the deck as anchorage reactions that introduce axial compression stresses, and only vertical reactions appear at bearings. Cables are prestressed to neutralise permanent loads. Therefore, tendons are tensioned, the deck is compressed and bending is reduced due to the upward deviation forces in the deck.

This typology has significant advantages. In the first place, bridge cost is reduced due to the high structural efficiency. This efficiency in the use of materials arises from the important contribution of the axial response compared with bending. Furthermore, because of this bending response reduction, it is possible to build more slender decks, so the advantages are not only economical, but also aesthetic. It is important to remark that, although deck height is reduced, the overall structural depth is increased due to the location of the cable system. Hence, substantial vertical clearance is required below the road for this type of bridge to be built.

Simply-supported underspanned suspension viaducts can span medium distances of around 80 m that would otherwise be bridged with a multi-span solution (for example a three-span bridge with lengths 24+32+24 m and a deck height of 1.5 m). Hence this solution avoids the erection of piers, with the subsequent reduction of construction time and other complications. In the case of a conventional solution with one single span, the resulting structure would be much heavier. Depending on the total weight of the bridge, an underspanned suspension structure could even be lifted into place by a crane, which would significantly reduce the disturbances produced during construction.

]]> On the other hand, their inherent lightness makes these bridges more sensitive to certain dynamic actions.*Keywords: *underspanned suspension bridges, road roughness, vehicle-bridge dynamic interaction, deck vibrations, penalty method

**Background**

First examples of bridges with the tensioned system located below the girder are found in the nineteenth century. In the 1830s G H Dufour built the Ile aux Barques Bridge in Geneva, Switzerland, with a main span of 33.5 m (Peters 1987), and J Smith erected the Micklewood Bridge (Figure 1), a 30 m single-span structure, in Scotland (Drewry 1832). In both structures the girder rests on chains that are anchored into the deck. Chains and deck are connected to each other by vertical struts that are in compression, in contrast with the vertical tensioned elements of traditional suspension bridges.

In the last quarter of the twentieth century this kind of structure started to reappear. Some authors ascribe this reappearance as having been initiated by Fritz Leonhardt's design for the Neckar Valley Bridge in the late 1970s. Leonhardt used this solution in order to avoid pier foundations in the hillside due to soil-related problems. End piers were removed in the design stage, and first and last spans are supported by cable systems located below the girder. These cables introduce vertical reactions in the deck by means of vertical elements and thus replace the eliminated piers; Figure 2 shows one of these systems. Ruiz-Terán and Aparicio (2007b) presented a wide and exhaustive state-of-the-art study of these bridges, considering the Neckar Valley viaduct as the birth of the typology (the term *under-deck cable-stayed bridges *is used by these authors), and therefore nineteenth century structures are not included in their work. Other researchers conceive these bridges as an evolution from externally prestressed bridges in which bending behaviour is enhanced by locating tendons outside cross-section bounds. These authors call these structures* bridges with highly eccentric external tendons *(Mutsuyoshi *et al *2010).

Examples of simply-supported unders-panned suspension bridges can be found all over the world, for instance the Tobu and Inachus footbridges in Japan, or the Truc de la Fare road bridge in France which spans 53 m. A very exhaustive review can be found in Ruiz-Terán and Aparicio (2007b), and several examples of footbridges are explained in Strasky (2005).

In spite of the important advantages pointed out above, underspanned suspension bridges are still unusual structures, and authorities remain reluctant to build them. This is due to limited knowledge about these viaducts. One of the first concerns is the shear capacity, as the girder depth is substantially reduced, but it has been proved experimentally that the shear resistance is even higher than in conventional girders (Mutsuyoshi *et al *2010). The structural behaviour of this kind of bridge was studied both analytically and experimentally in Witchukreangkrai *et al *(2000), Aravinthan *et al *(2001), and Ruiz-Terán and Aparicio (2007a). Some design criteria were proposed by Ruiz-Terán and Aparicio (2008) by using simply-supported 80 m span bridges.

Nowadays one of the main concerns is the lack of knowledge regarding the dynamic response of these bridges to different excitations. Ruiz-Terán and Aparicio (2009) dealt with the sudden breakage of stay cables due to a truck impacting with them, which was the principal reason why this solution for the Kirchheim overpass in Germany, designed by Jõrg Schlaich, was rejected in 1987.

Another important issue is the vibrations induced by the crossing of vehicles over the bridge. The reduction in girder depth leads to high traffic-induced accelerations. According to Ruiz-Terán and Aparicio (2010), the depth of the deck in road bridges of this kind, with short and medium spans, is governed by the serviceability limit state of vibrations. The same happens in pedestrian bridges (Tsunomoto & Ohnuma 2002). Muttoni (2002) studied the live load effects on this kind of road bridge by using static models. Ruiz-Terán and Aparicio also studied this issue in two papers (2007a; 2008) - in the 2007a study live load effects were considered by means of static models, while in the 2008 study dynamic analysis was performed with 400 kN vehicles crossing the bridge at 60 km/h (trucks were modelled as moving loads, so no dynamic interaction was taken into account).

**Research aim and limitations**

*et al*(2013a; 2013b).

The main objective of this work is to shed some light on the dynamic behaviour of these viaducts under traffic loads, as this knowledge is still limited. Important traffic-induced dynamic excitation will happen due to the extreme lightness and flexibility of this structural solution. In fact, accelerations in these decks are higher than in conventional bridges.

In this study one bridge of medium length is employed as a representative example, and only one set of support conditions is considered. Double bearings are set at both abutments, and therefore torsional rotation is prevented at those points. The use of sliding supports will not have a significant influence on the results, as vertical and torsional modes will not be affected. Another solution could be the use of only one support at one abutment. Structural torsional stiffness would decrease significantly and hence torsional frequencies of vibration would also be reduced. However, the low torsional rigidity inherent in this structural system makes this option inadvisable, and hence this solution is not considered in this work.

**VEHICLE-BRIDGE INTERACTION MODELS**

**Vehicle**

A two-axle truck is considered in this work. This vehicle model has also been used by other authors, e. g. Law and Li (2010), and is based on the H20-44 truck design loadings included in the AASHTO (1998) specifications. In this study a seat is added to the truck model in order to assess the dynamic effects on the driver.

The complete model consists of three rigid bodies that represent the box and both axles, plus one DOF mass that reproduces the driver seat. The vehicle model and its eight DOFs are depicted in Figure 3. The driver seat mass (*m _{d}*= 80 kg) is connected to the vehicle body by a vertical spring (

*k*= 10 507 N/m) and dashpot (

_{d}*c*= 876 Ν· s/m). Driver seat properties are taken from Zuo and Nayfeh (2007). Truck model mechanical properties can be found in Zhu and Law (2002).

_{d}**Structure**

An 80 m long simply-supported under-spanned suspension viaduct is considered in this work (Figure 4), as described in Ruiz-Terán and Aparicio (2009). The cable system is anchored in the deck over the two abutments and is deflected by two slightly inclined steel struts. The deck is a voided concrete girder with a depth of 1 m and a total width of 13.2 m (Figure 5). Double supports are set at both abutments, hence torsional rotation is not permitted at those points.

]]> This bridge is modelled by means of the finite element method - shells are employed to represent the deck and beams for the struts, and both kinds of elements consider shear deformation and adopt reduced integration. Cables are represented by means of truss elements. Forces and displacements induced in the bridge by the crossing vehicles are small enough to adopt linear elastic behaviour in the structure. Nonstructural masses representing other dead loads, such as pavement and safety walls, are included in the model. The damping matrix is built following the Rayleigh method by setting a damping ratio of 1% for the first mode (0.78 Hz) and also for 50.0 Hz. Figure 6 shows the first bending mode of vibration (1^{st}global mode) and the first torsion mode of vibration (3

^{rd}global mode).

**Vehicle-bridge dynamic interaction**

Vehicle-bridge vertical interaction is gathered through a node-to-surface contact at each vehicle tyre, and those contacts are implemented by means of the penalty method. This method does not introduce additional variables to the problem, and the computation is therefore faster. On the other hand, the constraint equation is only fulfilled approximately and some penetration will be unavoidable. This penetration will depend on the penalty parameter ε. When there is no contact, no forces are added and separation is reproduced in the numerical model. The whole system is set out as a fully coupled system of equations and it is solved by direct time integration with the HHT-a method (Hilber *et al *1977). This methodology has two main advantages over other fully coupled methods (Deng & Cai 2010; Yin *et al *2010; Neves *et al *2012): (1) modal superposition is not used for the bridge subsystem and hence structural non-linearities could be considered, and (2) vehicle tyres can lose contact with the deck surface. This separation capability is of interest in certain situations, as will be shown below.

As stated before, the constraint equation is fulfilled approximately when the penalty method is employed. Hence the solution is only an approximation of the correct enforcement of the constraint condition obtained with the Lagrange multiplier method (Wriggers 2002). In Figure 7(a) penetration at the right rear wheel when the vehicle crosses the bridge at 110 km/h is shown when Penalty and Lagrange multiplier methods are employed; road roughness is included. As can be seen, some penetration takes place with the penalty method. However, the contact force is the same in both cases. Figure 7(b) shows the vertical reaction under the same wheel during a short period of time, in order to facilitate visualisation. It can be concluded that the penalty method leads to correct results, although some penetration is inevitable.

**ROAD SURFACE DESCRIPTION**

Road roughness is generally the most important source of dynamic excitation in road traffic, hence a correct definition of the road surface is a key point in vehicle-bridge interaction problems. When the actual profile of a particular road stretch is not needed, but a set of profiles that are representative of a certain sort of road, stochastic definitions for the generation of synthetic profiles have to be used, as for example in Deng and Cai (2010). The fact that profiles under left and right tyres are different, but not independent, is seldom considered in this kind of simulation; therefore it is assumed that the road profile is constant across the deck width. In this paper those differences are considered by means of a procedure developed by the authors and described in Oliva *et al *(2013a). The influence of left-right dissimilarity on road vehicle-bridge interaction dynamics has been shown in Oliva *et al *(2013b).

In order to facilitate the generation and use of parallel profiles for other researchers, the authors have implemented a free-to-download program (http://w3.mecanica.upm.es/prpgenerator/index.php) that creates pairs of profiles with the mentioned procedure. This simple application named *PRPgenerator *was developed in MATLAB© and it is introduced for the first time in this paper. A brief description can be found in Appendix A.

The procedure employed in this work for considering that fact assumes hypotheses of road surface isotropy and homogeneity. In a homogeneous and isotropic road surface every straight profile has the same statistical characteristics, independent of its direction or position. Thus, parallel profiles along the road share statistical properties, but are not the same. Given that, the cross-Power Spectral Density *(G _{x}) *of the pair of parallel profiles at a certain distance is obtained from the direct Power Spectral Density or PSD

*(G).*The methodology is not explained here, for the sake of briefness, but it is detailed in Oliva

*et al*(2013a; 2013b).

The ISO-8608 (1995) PSD definition must be slightly modified in order to fulfil the isotropy admissibility conditions (Kamash & Robson 1977). It must be set constant and equal to *G(n _{a}) *at all wave numbers below a limiting wave number

*n*Thus it becomes acceptable, as its integral is bounded and the function remains monotonically non-increasing. This limit is set as

_{a}.*n*= 0.01 m

_{a}^{-1}, which is the minimum spatial frequency considered in roads (ISO-8608 1995). The PSD becomes:

where *G(n) *is the one-sided power spectral density for the spatial frequency or wave number *n *and *G(n _{0}) *is the one-sided power spectral density for the reference spatial frequency

*n*= 0.1 m

_{0}^{-1}. The value for

*G(n*is prescribed by ISO-8608 (1995) as a function of the road class.

_{0})Two parallel road profiles are generated as follows (Sayers 1998):

where *φ** _{i}*and

*θ*

*are two sets of random phase angles which are uniformly distributed from 0 to 2π. An example of two class A parallel profiles is depicted in Figure 8.*

_{t}

**NUMERICAL STUDIES**

For the truck and bridge presented in this section numerical simulations are used. Geometric nonlinearities can easily be considered in the proposed framework, but nonlinear effects in the considered scenarios have been shown to be negligible. Therefore, linear models are employed in this work because of the shorter computation time required. Six different running speeds were employed (30, 50, 60, 80, 110 and 120 km/h) in order to consider typical urban and highway vehicle speeds in different countries around the world. The road surface is assumed to be very good (class A, *G(n _{0}) *= 16·10

^{-6}m

^{3}) (ISO-8608 1995), 2000 spatial frequencies between 0.01 m

^{-1}and 10.0 m

^{-1}are used in the profile generation, and road irregularities are sampled every 2 cm in the longitudinal direction. Ten road surfaces (A01, A02, ... , A10) are generated so that the results are statistically significant. Thus, 60 different simulations had to be performed for this part of the work (6 speeds x 10 profiles). The H20-44 vehicle crosses the structure in the right lane with its right wheels at a distance of 2.55 m from the right edge of the bridge (see Figure 9). Driving on the right-hand side of the road is assumed in this work as it is the most common option in the world. The outcome of the study with left-hand traffic would be the same.

Results are obtained at every node on the bridge surface. These nodes are gathered in longitudinal lines whose transversal position and name are shown in Figure 9. The bridge surface is discretised with 14 elements in width and 81 elements in length, so that there are 15 longitudinal lines with 82 nodes per line.

]]> With regard to the Serviceability Limit State of vibration in road bridges, Eurocode EN1990 (EN1990:2002/A1+AC 2010) is very vague, and no indication is given with respect to the comfort of passengers or pedestrians. With respect to footbridges it is stated that pedestrian comfort criteria for serviceability should be defined in terms of maximum acceptable acceleration of any part of the deck. The recommended maximum value for vertical acceleration is set as 0.7 m/s^{2}. This limit is adopted in this work as a limitation to pedestrian comfort on the sidewalks of road bridges.

Maximum acceleration at every bridge surface node is obtained with the ten different road surfaces considered (A01, ... , A10). The mean value of those ten maxima is computed in order to get representative accelerations of the road class, instead of an absolute maximum that would represent a critical situation. Figure 10 shows the mean value of maximum vertical acceleration on the whole bridge surface for vehicle speeds of 110 km/h and 30 km/h. In both cases high accelerations are found right under the vehicle path due to local vibrations. At high speeds these accelerations are significantly higher than those of the rest of the surface. At low speeds the local apexes are not so markedly protruding. Clear peaks can be noticed in the two cases near the first abutment, located precisely under the truck (between lines A7 and V4); the maximum value is in line V4, because it belongs to the cross-sectional cantilever. These first abutment high values are caused by the sudden entrance of the vehicle. Vertical acceleration of the node located in line V4 at 1.0 m from the abutment, with *V *= 50 km/h and surface A01, is depicted in Figure 11 by the label *Raw; *high accelerations appear when each truck axle enters the bridge.

In Figure 12 maximum acceleration of some significant lines (V1, V4, V6 and A5) is depicted for 120 and 60 km/h. A peak at the beginning of V4 is easily noticable in both cases. For 120 km/h accelerations in nodes belonging to V4 are clearly the highest in the bridge, for 60 km/h local vibrations are less important and V4 accelerations are similar to those found in lines V1 and V6. Acceleration in the bridge longitudinal midline (A5) is always significantly lower than in the cantilevers, which fact indicates the relevance of torsion effects.

The Eurocode EN1990 limitation is generally fulfilled at 60 km/h, and is only violated near the first abutment because of the very local effects of the vehicle entrance, as explained before. At high speed (120 km/h) the maximum acceleration criterion is not satisfied in the vehicle side cantilever (V4-V6) all along the bridge length. It is remarkable that vertical acceleration along the bridge length is uniform for every vehicle speed, with the exception of the vehicle entrance region.

Two different standpoints exist regarding the acceleration effects on human comfort. The first considers that people are affected most by the largest peaks, and therefore maximum acceleration is limited as in Eurocode EN1990. According to the second school of thought it is assumed that the degree to which vibration may be noticed or tolerated is determined by some averaged effect over a period of time. International standard ISO-2631-1 (1997) specifies a method of evaluation of the effect of vibration on human beings by means of the weighted root-mean-square (RMS) acceleration which is defined as:

where *a _{w}(t) *is the weighted acceleration as a function of time and

*T*is the duration of the vibration.

Human response to vibration is a function of frequency and therefore data must be weighted in order to give greater prominence to frequencies where humans are most sensitive. Analogue transfer functions that determine the frequency weighting are specified in ISO-2631-1 (1997). Those functions are transformed into digital filters by means of bilinear transformations.

For a standing person subjected to vertical accelerations underneath their feet, as is the case with a pedestrian on a bridge, *W _{k }*weighting is applicable. Its parameters are given in ISO-2631-1 (1997). As an example, Figure 1 shows vertical weighted acceleration

*(a*in line V4 at 1.0 m from the first abutment; high-frequency peaks are eliminated by the frequency weighting. Reaction to vibration depends on many factors (frequency, duration of vibration, activity, age, etc), and therefore absolute limits cannot be established, only indications of likely human reactions. Some limits are given by the Department of Environment and Conservation of New South Wales in Australia (DEC 2006), but no specific values are provided for bridges. The highest admissible value in buildings is set as 0.08 m/s

_{w}(t))^{2}. A weighted RMS acceleration of 0.08 m/s

^{2}is also adopted in British Standard BS 6472:1992 (1992) as a value below which there exists a low probability of adverse comment from users of residential buildings during the day, with a time of exposure to vibration lower than 225 seconds (the bridge in this study is crossed by a person in approximately one minute). Building limits are not directly applicable to bridges. The same value (0.08 m/s

^{2}) is defined in the German guideline VDI 2057 (2002) as a limit above which vibration is strongly perceptible by a human being. Higher values would probably be tolerated by pedestrians on a road bridge, but it seems reasonable to adopt 0.08 m/s

^{2}as a guidance value for human comfort on viaduct decks, although presumably conservative.

As for maximum acceleration, weighted RMS acceleration is computed for every bridge surface node with the ten different road surfaces considered, and the mean value is calculated. Figure 13 shows the mean value of the weighted RMS acceleration all over the bridge surface. First the abutment peaks and high values on the vehicle trajectory disappear, and the highest values are now obtained in the bridge edges (V1 and V6), both with high and low speeds. In Figure 14, values in some significant lines (V1, V4, V6 and A5) are depicted for 120 and 60 km/h.

]]> Weighted RMS remains at acceptable levels on the whole deck surface for all vehicle speeds. It is remarkable that differences between the left and right bridge edges, which were evident when using maximum acceleration values, almost vanish if weighted RMS is employed.With respect to the Dynamic Amplification Factor (DAF) of vertical displacement, the mean value of the ten cases is depicted in Figure 15 for the bridge midline (A5) considering different vehicle speeds.

High DAF_{y} values appear and they are highly influenced by vehicle speed. DAF values as high as 1.54 are reached with *V *= 120 km/h near the first strut. These high values indicate the huge relevance of dynamic effects on this kind of structure.

Regarding vehicle vibration, Figure 16 shows vertical acceleration at the driver seat when the vehicle runs at 120 km/h with profile A07. The vertical response of the vehicle on the bridge and on a rigid road is compared under the same road surface conditions in order to assess the bridge flexibility influence. The results in Figure 16 show that bridge deflection has an effect on the driver seat.

For vertical vibrations on the seat surface, *W _{k}*weighting is also used (ISO-2631-1 1997). Figure 17 shows the mean value of weighted RMS acceleration on the driver seat considering the ten surfaces - results on the bridge and on a rigid road are compared. The increment is very small (maximum of 4%), and it can therefore be concluded that the bridge flexibility influence on driver comfort is of no significant relevance.

**Effects of pothole presence**

The eventual presence of a pothole of considerable size is considered by adding a 50 cm long and 5 cm deep defect on road surface A01, located in the midspan and affecting both vehicle sides. Tyre contact is lost in the pothole and high reactions appear when it is regained (Figure 18). The wheel-bridge separation capabilities of the interaction model are necessary for a correct simulation of this scenario. This behaviour produces very high accelerations in the bridge surface (Figure 19); even at low speeds, higher values are beyond 6.5 m/s^{2} when *V *= 50 km/h.

**Road quality influence**

Profiles of different road classes (A or Very Good, B or Good, and C or Medium) are employed in order to assess road class influence in the dynamic behaviour of the deck. Profiles of different road classes are related by a scale factor:

]]> Thus B and C surfaces can be obtained by multiplying A pairs of profiles by 2 and 4 respectively. Profile A10 is employed in this section and the vehicle speed is set as 110 km/h. Maximum vertical acceleration at lines V1 and A5 is depicted in Figure 20. As can be seen, road quality has significant influence on the deck vibration, and the increase in maximum vertical acceleration is of relevance. The Eurocode limit is not reached in the midline (A5), but is exceeded in line V1 with road level B.Road class also has influence in the weighted RMS acceleration (Figure 21) - values increase when road quality declines, and exceed 0.08 m/s^{2} manifestly in line V1 with road classes B and C. In the midline (A5) weighted RMS acceleration also increases, but values remain acceptable.

**CONCLUSIONS**

Vehicle-induced dynamics in underspanned suspension bridges were studied in this paper. Analyses were performed by means of a fully coupled vehicle-bridge dynamic interaction model. The vehicle was represented through a multibody system (MBS); the bridge was modelled with the finite element method (FEM), and interaction was gathered by means of contacts with a penalty formulation. The fully coupled system equations were solved by direct integration in time. Penetration inherent in the penalty method was shown to have no effect on the relevant results against the Lagrange multiplier method, where restriction is perfectly satisfied and no penetration takes place.

Road roughness was considered in such a way that the fact that profiles in the left and right tyres were different, but not independent, were taken into account. In order to facilitate their consideration, a parallel road profiles generation programme (*PRPgenerator*)* *was developed by the authors, and this is available on the web for free download as a standalone application.

A very good road class (A) was considered in the study (ISO-868 1995). Vertical acceleration peaks appeared directly under the vehicle path due to local vibrations. Absolute maximum values arose near the first abutment because of sudden vehicle entrance. The Eurocode EN1990:2002/A1+AC (2010) maximum acceleration criterion was generally fulfilled on the deck surface during truck runs at urban speeds (50-60 km/h). With highway speeds (110-120 km/h) this limitation was not satisfied in the vehicle proximity, but on the rest of the deck it was. Frequency weighting (ISO-2631-1 1997) reduced the high frequency peaks on acceleration histories, and the largest weighted RMS values were not found in the vehicle path but on both bridge edges. In addition, differences between left and right viaduct edges were negligible. Weighted RMS values were low enough to be tolerated on the whole deck surface for every vehicle speed.

Significant differences were found between the bridge longitudinal midline and edges. Accelerations were clearly higher on the edges due to the excitation of the deck torsion. Therefore it would be convenient to locate sidewalks along the deck longitudinal midline, if possible, or set vehicle lanes near that line in order to reduce the excitation of this rotation.

Road class has an evident influence on pedestrian comfort. Hence, good construction and maintenance are of great importance in this kind of viaduct due to its dynamic sensitivity. The presence of a significant pothole could induce very high deck accelerations. Therefore, careful conservation is important. With respect to vehicle behaviour, bridge flexibility causes additional vibration in the driver seat. However, the effect of the additional acceleration on driver comfort is very small - the maximum increase of the weighted RMS acceleration on the driver position was only 4%.

]]>

**ACKNOWLEDGEMENTS**

This work was motivated by lectures given by Dr Ana María Ruiz-Terán (Imperial College London) at the Technical University of Madrid. The authors wish to acknowledge this source of inspiration. The aid of the *Ministerio de Ciencia e Innovación *of the Spanish government through sub-programme *INNPACTO *and project *VIADINTEGRA *(Ref IPT-370000-2010-012) is greatly appreciated. The authors are grateful also for support provided by the Technical University of Madrid, Spain.

**APPENDIX A ( PRPgenerator)**

*PRPgenerator *is a simple-to-use and free application that generates pairs of parallel road profiles. It has been implemented in Matlab© and is available as stand-alone software. The program uses Power Spectral Density definition from ISO-8608 (1995) and assumes the homogeneity and isotropy of the road surface. *PRPgenerator *has been developed by the Computational Mechanics Group of the School of Civil Engineering at the Technical University of Madrid (UPM). Its theoretical background can be found in Oliva *et al *(2013a; 2013b), and the application is presented for the first time in this paper. It computes cross-Power Spectral Density and coherence function for parallel profiles and generates profiles with the total length, sample distance and number of spatial frequencies specified by the user. Different road classes can be considered and any vehicle width may be used. Figure 22 shows the main and only window of the application.

**REFERENCES**

AASHTO (American Association of State Highway and Transportation Officials) 1998. *LRFD Bridge design specifications. *Washington DC: AASHTO. [ Links ]

Aravinthan, T, Matsui, T, Hamada, Y & Shinozaki, H 2001. Structural analysis of Torisaki River Park bridge - An innovative PC bridge with large eccentric external tendons. *Proceedings, 11th Symposium on Development in Prestressed Concrete, *pp 49-54. [ Links ]

BS 6472 1992. *Guide to evaluation of human exposure to vibration in buildings. *London: British Standards Institution. [ Links ]

DEC (Department of Environment and Conservation, New South Wales) 2006. *Assessing vibration: A technical guideline. *Sydney, Australia: DEC. [ Links ]

Deng, L & Cai, C S 2010. Development of dynamic impact factor for performance evaluation of existing multi-girder concrete bridges. *Engineering Structures, *32: 21-31. [ Links ]

Drewry, C S 1832. *A memoir on suspension bridges. *London: A and R Spottiswoode. [ Links ]

EN 1990. 2002/A1+AC 2010. *Eurocode - Basis of structural design. *Brussels: European Committee for Standardization (CEN). [ Links ]

Hilber, H M, Hughes, T J R & Taylor, R L T 1977. Improved numerical dissipation for time integration algorithms in structural dynamics. *Earthquake Engineering and Structural Dynamics, *5: 283-292. [ Links ]

ISO-2631-1 1997. *Mechanical vibration and shock -Evaluation of human exposure to whole body vibration - Part 1: General requirements. *Geneva: International Organization for Standardization (ISO). [ Links ]

ISO-8608 1995. *Mechanical vibration - Road surface profiles - Reporting of measured data. *Geneva: International Organization for Standardization (ISO). [ Links ]

Kamash, K M A & Robson, J D 1977. Implications of isotropy in random surfaces. *Journal of Sound and Vibration, *54: 131-145. [ Links ]

Law, S & Li, J 2010. Updating the reliability of a concrete bridge structure based on condition assessment with uncertainties. *Engineering Structures, *32(1): 286-296. [ Links ]

Mutsuyoshi, H, Hai, N D & Kasuga, A 2010. Recent technology of prestressed concrete bridges in Japan. *Proceedings, *Joint IABSE-JSCE Conference on Advances in Bridge Engineering-II, Amin, A F M S, Okui, Y & Bhuiyan, A R (Eds), Dhaka, Bangladesh. [ Links ]

Muttoni, A 2002. Brücken mit vorgespannter stahlunterspannung. *Stahlbau, *71(8): 592-597. [ Links ]

Neves, S G M, Azevedo, A F M & Calçada, R 2012. A direct method for analyzing the vertical vehicle-structure interaction. *Engineering Structures, *34: 414-420. [ Links ]

Oliva, J, Goicolea, J M, Astiz, M A & Antolín, P 2013a. Fully three-dimensional vehicle dynamics over rough pavement. *Proceedings of the ICE - Transport, *166(3): 144-157. [ Links ]

Oliva, J, Goicolea, J M, Astiz, M A & Antolín, P 2013b. Relevance of a complete road surface description in vehicle-bridge interaction dynamics. *Engineering Structures, *56: 466-476. [ Links ]

Peters, T F 1987. *Transitions in engineering - Guillaume Henri Dufour and the early 19th century cable suspension bridges. *Basel, Switzerland: Birkhauser Verlag. [ Links ]

Ruiz-Terán, A M & Aparicio, A C 2007a. Parameters governing the response of under-deck cable-stayed bridges. *Canadian Journal of Civil Engineering, *34: 1016-1024. [ Links ]

Ruiz-Terán, A M & Aparicio, A C 2007b. Two new types of bridges: Under-deck cable-stayed bridges and combined cable-stayed bridges - The state of the art. *Canadian Journal of Civil Engineering, *34: 1003-1015. [ Links ]

Ruiz-Terán, A M & Aparicio, A C 2008. Structural behaviour and design criteria of underdeck cable-stayed bridges and combined cable-stayed bridges. Part 1: Single-span bridges. *Canadian Journal of Civil Engineering, *35: 938-950. [ Links ]

Ruiz-Terán, A M & Aparicio, A C 2009. Response of under-deck cable-stayed bridges to the accidental breakage of stay cables. *Engineering Structures, *31: 1425-1434. [ Links ]

Ruiz-Terán, A M & Aparicio, A C 2010. Developments in under-deck and combined cable-stayed bridges. *Proceedings of the Institution of Civil Engineers (UK) - Bridge Engineering, *163: 67-78. [ Links ]

Sayers, M W 1998. *Dynamic terrain inputs to predict structural integrity of ground vehicles. *Ann Arbor, MI: University of Michigan Transportation Research Institute. [ Links ]

Strasky, J 2005. *Stress ribbon and cable-supported pedestrian bridges. *London: Thomas Telford. [ Links ]

Tsunomoto, M & Ohnuma, K 2002. *Self-anchored suspended deck bridge - Pedestrian bridge of Tobu recreation resort. *National Report - Recent works on prestressed concrete structures, Osaka, Japan: Japan Prestressed Concrete Engineering Association. [ Links ]

VDI 2057 (2002). *Human exposure to mechanical vibrations - Whole-body vibration. *Dusseldorf, Germany: Association of German Engineers. [ Links ]

Witchukreangkrai, E, Mutsuyoshi, H, Aravinthan, T & Watanabe, M 2000. Analysis of the flexural behavior of externally prestressed concrete beams with large eccentricities. *Transactions of the JCI (Japan Concrete Institute), *22: 319-324. [ Links ]

Wriggers, P 2002. *Computational contact mechanics. *New York: Wiley. [ Links ]

Yin, X, Fang, Z, Cai, C S & Deng, L 2010. Non-stationary random vibration of bridges under vehicles with variable speed. *Engineering Structures, *32: 2166-2174. [ Links ]

Zhu, X Q & Law, S S 2002. Dynamic load on continuous multi-lane bridge deck from moving vehicles. *Journal of Sound and Vibration, *251: 697-716. [ Links ]

Zuo, L & Nayfeh, S A 2007. H_{2} optimal control of disturbance-delayed systems with application to vehicle suspensions. *Vehicle System Dynamics, *45(3): 233-247. [ Links ]

**Correspondence:**

Javier Oliva

Department of Mechanics and Structures

School of Civil Engineering

Technical University of Madrid

28040 Madrid

Spain

T +34 91 336 5358 ]]>
F: +34 91 336 5367

E: joliva@mecanica.upm.es

José M. Goicolea

Department of Mechanics and Structures

School of Civil Engineering

Technical University of Madrid

28040 Madrid

Spain

T +34 91 336 6761

F: +34 91 336 5367 ]]>
E: jose.goicolea@mecanica.upm.es

Pablo Antolín

Department of Mechanics and Structures

School of Civil Engineering

Technical University of Madrid

28040 Madrid

Spain

T: +34 91 336 5358

F: +34 91 336 5367

E: pablo.antolin@upm.es

Department of Mechanics and Structures

School of Civil Engineering

Technical University of Madrid

28040 Madrid

Spain

T: +34 91 336 6760

F: +34 91 336 5367

E: miguel.a.astiz@upm.es

]]>

DR JAVIER OLIVA Is a research assistant In the Department of Mechanics and Structures in the School of Civil Engineering at the Technical University of Madrid (UPM), Spain. He received his PhD from the same institute in 2011. Before joining UPM he worked as a consulting structural engineer In a private company. His research interests include structural dynamics and biomechanics. He lectures at Saint Louis University (Madrid campus), and has authored or co-authored 15 conference and scientific journal papers, and participated in several research projects with companies and administrations. Dr Oliva also collaborates with AR2V an independent firm spesia lising in the design of bridges and other structures.

PROF JOSÉ M. GOICOLEA received his PhD from the University of London in 1985. Since 1993 he has been full professor in the School of Civil Engineering at the Technical University of Madrid (UPM), Spain. His research is mainly focused on structural dynamics and cardiovascular biomechanics. He has authored more than 100 conference and scientific journal papers, and has developed around 60 research projects for companies and administrations. He is Vice-President of the Spanish Society for Numerical Methods in Engineering (SEMNI), and is also a member of the General Council ofthe International Association for Computational Mechanics (IACM).

]]>

DR PABLO ANTOLÍN is a research assistant in the Department of Mechanics and Structures in the School of Civil Engineering at the Technical University of Madrid (UPM), Spain. He received his PhD from the same institute in 2013. His research is mainly focused on structural dynamics and numerical methods, Including isogeometric analysis. During his doctorate studies he taught computational mechanics in the School of Civil Engineering at UPM. He has also participated in several research projects concerning, for example, the turbulent wind effet on tall high-speed railway bridges or membrane structures. He has authored or co-authored several conference and scientific journal papers.

PROF MIGUEL Á ASTIZ has more than 30 years' experience In major bridge design. He shares his professional activity between academia, where he currently holdsthe Chair In Bridge Design In theSchool of Civil Engineering atthe Technical University of Madrid (UPM), Spain, and professional practice with Carlos Fernández Casado S L, a prestigious global design firm. Prof Astiz has designed a multitude of major bridges, the River Suir cable-stayed bridge in Waterford, Ireland, being one of his most recent projects. He is President ofthe Spanish Scientific Association for Stuctural Concrete (ACHE).

]]>