On the vertical coupling effect of ballasted tracks in multi–span simply–supported railway bridges under operating conditions

March 2, 2022 Structure and Infrastructure Engineering SIE˙paper To appear in Structure and Infrastructure Engineering Vol. 00, No. 00, Month 20XX, 1–31 Original Paper On the vertical coupling effect of ballasted tracks in multi–span simply–supported railway bridges under operating conditions Emma Molinera,∗, Maria Dolores Mart´ınez-Rodrigoa, Pedro Galv´ınb,c, Josep Chord`a-Monson´ısa,b, Antonio Romerob aUniversitat Jaume I, Department of Mechanical Engineering and Construction, Avda. Sos Baynat s/n, ES-12071 Castell´on, Spain bEscuela T´ecnica Superior de Ingenier´ıa, Universidad de Sevilla, Camino de los Descubrimientos, ES-41092 Sevilla, Spain cENGREEN, Laboratory of Engineering for Energy and Environmental Sustainability, Universidad de Sevilla, Camino de los Descubrimientos, ES-41092 Sevilla, Spain (Received 00 Month 20XX; final version received 00 Month 20XX) Abstract This contribution investigates the vertical coupling exerted by ballasted tracks on the vertical response of bridges composed either of (i) several successive simply-supported spans with weak coupling between them due to the continuous track; or (ii) adjacent single-track decks conforming a double-track bridge, in which interaction effects are induced due to the transverse continuity of the ballast layer. To this end, 2D and 3D track-bridge interaction Finite Element models are implemented, which consider a three-layer discrete and explicit idealization of the track. The 2D track-bridge interaction model is used to perform sensitivity analyses on the track parameters, which have revealed that the ballast shear mechanisms along the track may significantly affect the train-induced vibrations under resonant conditions. Then, the influence of the ballast coupling on the response of twin adjacent decks is investigated with a 3D track-bridge interaction model. To this end, this model is updated based on the results of an experimental campaign performed on a real bridge composed of two SS spans and two single track twin adjacent decks. The numerical-experimental comparison shows an evident dynamic vertical coupling between the bridge decks and reveals the importance of including the ballasted track in the modelling process of these structures. Keywords: Multispan railway bridges, ballasted track, track-structure interaction, span coupling, experimental measurements, resonance, vertical acceleration. 1. Introduction The sustained development of transportation systems and, especially, railways for passenger and freight transport has set a milestone in most developed countries in recent decades. The progressive increase in the trains operational speeds and in the transportation capacity constitutes a challenge for infrastructures, which need to accomplish strict requirements to ensure the admissible levels of quality, safety and reliability. In particular, the Spanish railway network currently has a total of 135 km of track that runs over more than 6000 bridges. The dynamic effects on such structures, directly related with the vehicle velocity, have become an issue of interest and concern for scientists and engineers, especially since the advent of High-Speed (HS) (Fr´yba (2008)). An excessive level of vertical acceleration at the deck platform can lead to passenger discomfort, fatigue problems in the long term, to an increase of the maintenance costs of the lines, misalignment of the rails as a result ∗Corresponding author. Email addresses: molinere@uji,es, mro[email protected], p[email protected], jc[email protected], [email protected] 1 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper of premature deconsolidation of the ballast layer and even to the loss of contact between wheel and rail with the resultant increase of the risk of derailment. Consequently, the Serviceability Limit States for Traffic Safety and, in particular, the vertical acceleration at the deck has become one of the most restrictive requirements in the design of new bridges (CEN/TC250 (2005)). Bridges composed of simply–supported (SS) spans with short to medium span lengths are especially critical in this regard, which may become particularly relevant at resonance (Hoorpah (2008); Zacher and Baeb¸ler (2008)). Aside from the design of new structures, it is essential to realistically assess the performance of existing bridges when facing new traffic requirements. In this context, there is a need of (i) accurate and computationally efficient numerical models to predict the response, (ii) experimental data to provide up-to-date information on the condition of the structures, and (iii) appropriate analysis methodologies as essential tools for the decision–making process. However, the dynamic response of a bridge under the circulation of a train is complex and it is affected by several factors. The most obvious ones are the bridge geometrical and mechanical properties, the scheme of train axles, and the speed of circulation, being also the most certain (Rocha et al. (2012)). Additionally, there are other factors that are much more uncertain to determine that significantly affect the response of the bridge such as structural damping (Rocha et al. (2014)) and various interaction mechanisms which modelling is not trivial and is currently under investigation, being the most relevant vehicle– structure (Jahangiri and Zakeri (2017); Zakeri et al. (2014)), track–structure (Ticona et al. (2020); Zhai et al. (2019)) and soil–structure (Galv´ın et al. (2021); Zangeneh et al. (2019)) interaction. In engineering consultancies, simplified numerical models that generally disregard these mechanisms are most often used, given the uncertainties and the considerable computational cost involved in their modelling. This work contributes to the investigation on the vertical coupling effect exerted by the ballasted track on railway bridges composed of SS spans in two situations: (i) in SS multi–span viaducts, where a weak coupling between consecutive spans occurs through the track (see Figure 1(a)-(b)); and (ii) in double–track bridges composed of two adjacent structurally independent decks, which vertical motion is coupled through the ballast, quite common in countries like Germany (Rauert et al. (2010)) and Spain (see Figure 1(c)-(d)). The contribution of the ballasted track to the dynamic response of the bridge remains a matter of recent investigation (Bornet et al. (2015); Melo et al. (2020); Rebelo et al. (2008)). A few authors have concluded from field measurements that the presence of the rail may influence the boundary conditions at the end sections where the track continuity introduces a moment of resistance (Rebelo et al. (2008); Rigueiro et al. (2006a,b)). According to some authors, this is one of the reasons for the discrepancies found between calculated and measured modal parameters in short or medium span bridges (Bornet et al. (2015)). A few references can be found in the literature devoted to the coupling effect of the ballast over adjacent single-track decks in double-track bridges, and its influence on the dynamic response. Among them, Rebelo et al. (2008) presented the results of an experimental campaign on a number of single–span bridges of this typology. The authors pointed out that the coupling effect between twin slabs is particularly evident in skewed configurations, and attributed the high level of damping identified to this effect. They also pointed out the inherently non linear character of the interaction, resulting the stiffness associated to the ballast coupling higher for lower amplitudes of vibration. Rigueiro et al. (2010) compared the numerical predictions from a finite element (FE) plate model of two SS twin decks coupled through the ballast layer. Three discrete track models were compared. The authors concluded that for frequencies up to 10 −15 Hz the influence of the track model is minor, but for higher frequencies the track acts as a filter. These authors also expressed the need to investigate the response of such structures at resonance. Rauert et al. (2010) presented an intense experimental campaign where the load transfer mechanism between ballast-connected structures as well as the contribution of the ballast layer to the global bending stiffness of SS bridges were investigated. The authors showed that ballast–coupled bridge decks provide a considerable amount of additional stiffness, particularly during train passages. Bonif´acio et al. (2014) assessed the dynamic behaviour of a short span filler-beam railway bridge composed of two independent 2 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper Figure 1. (a)-(b) SS multi–span railway viaducts in Spanish HS lines; (c)-(d) railway bridge in Spanish conventional line composed of two adjacent structurally independent decks. decks and several SS spans subjected to the TGV-double train. They pointed out an obvious coupling between the spans in the mode shapes and an important transfer of vibrations from the loaded to the unloaded deck under operational conditions. Fewer still are the works found in the literature devoted to the so–called weak coupling induced by the track between successive spans in SS viaducts. Among them, Chellini et al. (2011) presented the experimental identification of modal parameters in two steel-concrete composite viaducts of the Torino-Milano HS line in Italy. The work reveals a clear dynamic interaction between the spans. The authors implemented a 3D FE comprehensive numerical model of a single span of each bridge and introduce artificial boundary conditions between the spans in order to reproduce the continuity of the tracks, which is found to be essential. Liu et al. (2014) evaluated the dynamic response of multi–span SS viaducts. The authors implemented and calibrated two numerical models to take into account this effect: a single-span model with modified boundary conditions and a component synthesis model (CMS) or reduced order model of the entire viaduct. A better representation of the vertical acceleration of the viaduct was obtained with the CMS model at a reasonable computational cost. Finally, Bornet et al. (2015) performed experimental tests on a SS single-span single-track steel truss railway bridge before and after the ballasted track was placed, and calibrated numerical models for both situations. They concluded that the ballast superstructure provided an additional stiffness of about 22 −27% for the three lowest modes, and showed how the stiffness of the ballast and continuity of the track had a significant effect on the second vertical bending mode due to an apparent change in the supports behaviour. Even though the effect of the ballast track on the dynamic response of railway bridges has already been a matter of study by scientists and engineers, in the opinion of the authors some important aspects have not yet been addressed in depth or with sufficient generality. Among them (i) stiffness and damping parameters in discrete track models used for similar track geometries greatly differ among publications. Moreover, no studies have been presented stating the effect and extent of these parameters variations; (ii) many authors disregard the shear transmission between the ballast mobilized masses, increasing the track-bridge stiffness and leading to energy dissipation; (iii) in multi-span simply-supported bridges the contribution of adjacent spans is generally disregarded and no studies have been found investigating under what circumstances the maximum acceleration 3 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper response in operating conditions may be higher when considering the track coupling between successive spans; (iv) the ballast coupling between adjacent twin decks and the influence of individual track parameters on the evolution of the natural frequencies and modal shapes in these types of bridges, especially for modes with an important transverse deformation of the cross-section, has not been investigated with the sufficient level of detail. In this paper, the influence of the ballasted track–bridge vertical interaction on the response of railway bridges composed of SS spans is investigated, as an extension of a preliminar study performed by the authors in Moliner et al. (2020). The original contributions of the work presented herein with respect to previous investigations are: (i) focusing on the effect of ballast shear transmission mechanisms along the track and between structurally independent parts of the bridges, i.e. consecutive spans and adjacent twin decks, and on how discrete track models reproduce these phenomena; (ii) analysing the influence and extent of the different track parameters on the bridge dynamic properties and (iii) on the bridge response at resonant, not resonant and cancellation situations; (iv) detecting under what circumstances the use of single-span models for multi-span bridges may not be on the safe side; (v) validating the previous conclusions through the comparison with experimental measurements on a real structure from a conventional railway line composed of two SS spans and two twin adjacent decks per span weakly connected through the ballast. 2. Numerical models In the range of lengths of interest (i.e. 10 −25 m) usual deck typologies in the Spanish railway system are solid or voided concrete slabs, prestressed concrete filler beams encased in a concrete pseudo-slab and precast decks composed of double-T beams (see Figure 2). This justifies the nature of the numerical models adopted in this investigation. Figure 2. Common railway bridge decks for short-to-medium lengths: (a) pre-cast double-T girder deck; (b) concrete filler beam deck; (c) solid concrete slab; (d) voided concrete slab. In the next section, the three-layer discrete track model is presented, common for the 2D and 3D FE track-bridge interaction models. In the former, the bridge deck is simulated by means of a Bernoulli-Euler (B-E) beam. This planar model is used to perform a comprehensive sensitivity analysis of the track individual parameters, and to investigate the adequacy of predicting the maximum acceleration in multi-span bridges with single-span planar models. On the other hand, in the 3D track–bridge model the deck is represented as an isotropic Kirchhoff plate, connected to B-E beams accounting for the longitudinal girders in the case of girder decks. This model is used to analyse the coupling introduced by the ballast layer along the deck borders between structurally independent parts of the bridge, i.e. adjacent twin decks and consecutive spans, and for comparison with the experimental response in the last section. 4 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper 2.1. Track model A three–layer discrete track model (see Figure 3) as the one proposed by Zhai et al. (2004, 2009) is adopted. The rail is modelled as a B-E beam resting on equidistant supports. The damping and vertical stiffness of the rail pads, Cpand Kp, mobilised ballast, Cband Kb, and subgrade, Cfand Kf, are concentrated at each sleeper location. Also, to account for the interlocking among the ballast granules, vertical spring-damper elements simulating the associated shear stiffness and damping with constants Kwand Cw, are introduced to link the relative vertical displacements between adjacent ballast masses. The parameters Kb fand Cb fstand for the vertical subgrade stiffness and damping on the bridge deck which, in this study, is set to 100·Kfand 0, respectively, assuming that the ballast rests directly on the deck slab. Finally, Er,Iyr,mrstand for the rail Young Modulus, cross-section moment of inertia with respect to the Yaxis defined in Figure 4 and linear mass. As can be derived from Figure 3, the track model used in this work only considers the response in the vertical plane, while the longitudinal track-bridge interaction is neglected in a first approach due to the high bending stiffness of railway bridge decks. This is consistent with recent publications that focus on the transverse dynamic response of the bridge (Cantero et al. (2016); Xiao et al. (2020)). In the absence of braking forces or temperature gradients, the coupling effects in the vertical direction prevail over those in the longitudinal one, and thus they may be evaluated in a decoupled manner (Wu et al. (2021)). In the 2D model described in Subsection 2.2, mass, stiffness and damping coefficients as well as the rail area and moment of inertia are multiplied by a factor of 2 (as indicated in Figure 3), as only one rail is explicitly included in the model. This idealisation admits Ahlbeck hypothesis which states that the load transmitted from a sleeper to the ballast coincides approximately with a cone distribution with an inclination defined by the ballast stress pervasion angle corresponding to the Poisson ratio (Ahlbeck et al. (1978)). For the particular expressions used to calculate Mb,Kband Kf, the reader is referred to Zhai et al. (2004). The rest of the parameters are extracted from the literature, as specified in Section 3.1. 2.2. 2D track-bridge interaction model In the planar model, the bridge is simulated by Nsp isostatic beams, where Nsp stands for the number of spans, resting on vertical elastic supports that represent the neoprene bearings. In Figure 3, ¯ Kn bi stands for the constant equivalent vertical stiffness in each external section of the deck in the i-th bridge span. Lbi,Ebi,Iybi and mbi are the span length, bridge deck Young Modulus, cross-section moment of inertia with respect to the Yaxis and linear mass of the i-th bridge span, respectively. The only interaction between successive spans is due to the track. A track length of Lr,prev = 20 m is included before and after the bridge, equivalent to more than 30 times the sleeper distance Dsl, which is considered adequate according to previous publications (Clark et al. (1982); Lou (2005)) and to a convergence test performed to this end. As for the train excitation, a constant moving load model is adopted, therefore neglecting vehiclestructure interaction (VBI) and track irregularity effects. Since the main purpose of the work is the evaluation of the vertical coupling exerted by the ballasted track in consecutive spans and adjacent decks, this is a deliberate decision as (i) it is intended to isolate the effect of track components from other interaction mechanisms that can be overlapped in resonance situations; (ii) only the bridge response is of interest; (iii) VBI effects are mainly relevant at resonance (Dom´enech et al. (2014)); (iv) reliable data for the train suspension systems included in Section 4 are not available at the time of writing. Also, previous studies have shown that the effects of track irregularities are noticeable on the vertical response of the vehicle itself and are perceptible in the soil far away from the track, being less important for the infrastructure at the bridge site in the frequency range of interest (Peixer et al. (2021); Romero et al. (2012, 2013)). The model is implemented in ANSYS R17.1. The full mass, stiffness and damping matrices are exported to MATLAB R2017b and the equations of motion are integrated in the time domain applying the Newmark-beta constant 5 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper Figure 3. 2D track-bridge interaction model. acceleration algorithm. The time step is defined as the minimum between 1/50 times the smaller period of interest and 1/20 times the travelling time of an axle load between two consecutive sleepers. 2.3. 3D track-bridge interaction model The 3D FE track-bridge interaction model is outlined in Figure 4. This model permits one to evaluate the ballast vertical shear transmission between single-track adjacent decks of the same span (from now on deck coupling) and between consecutive spans (from now on span coupling) through the shared ballast layer. The main features of the model are: •The decks are represented as isotropic Kirchhoff plates discretised into shell FE with 6 degrees of freedom (dof) per node. The element size is chosen to adequately reproduce the wavelengths of the modes with frequencies up to 30 Hz as per Eurocode (EC) CEN/TC250 (2005). •In each slab, different mass density elements are defined in order to concentrate the weight of the nonstructural elements over the corresponding contributive areas. •The longitudinal girders are represented as B-E beam elements with 6 dof per node. These nodes are connected to the plate through kinematic constraints. The distance between the plate and the beam nodes equals the existing vertical offset between the slab neutral plane and the center of mass of the girders. •Vertical interaction due to the shared ballast distributed over the twin decks and between consecutive spans is simulated as discrete spring elements with stiffness constants KwL and KwT , respectively. •The discrete track models adopted in the longitudinal (X) and transverse (Y) directions are represented in Figure 4. In the transverse direction, the ballast shear stiffness and damping constants are designated as Kwr and Cwr, respectively. •The vertical stiffness of the laminated rubber bearings is included at the centre of each bearing by means of a vertical spring with elastic constant Kn bi. •As in the 2D track-bridge model, constant moving loads are adopted for the railway excitation, therefore neglecting VBI. •To reduce computational costs, in this case the dynamic equations of motion are transformed into modal space and integrated in the time domain by the Newmark-beta linear acceleration algorithm as in Bonif´acio et al. (2014); Chellini et al. (2011); Liu et al. (2014). The time-step is defined as 1/25 times the smallest period of interest. 6 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper Cwr Kwr Figure 4. 3D track-bridge interaction model. 3. Sensitivity analysis of ballasted track parameters 3.1. Track nominal parameters and methodology In Figure 5 the values adopted by different authors of the track parameters Kp,Cp,Kband Cbare included. The high dispersion among the values is evident and constitutes one of the reasons for the sensitivity analysis presented in what follows. Regarding the ballast shear stiffness and damping, Kωand Cω, only a few references are found in the literature providing approximate values for similar track characteristics. Particularly Kω∈5.4·107,7.84 ·107N/m according to Rodrigues and Dimitrovov´a (2013), Zhang et al. (2017), Punetha et al. (2020), and Cω= 8 ·104Ns/m is most often used as in Zhai et al. (2004). First, a set of reference track parameters is defined as a basic set for subsequent studies (see Table 1). As indicated in Table 1 the majority of the parameters are calculated according to Zhai et al. (2004). This constitutes a well-established standard and it is proximate to the average of the values included in Figure 5. The rail, rail pads and sleeper properties are selected from Standards and data from the Spanish railways (CEN/TC256 (2017); Ministerio de Fomento (2010); Melis (2007); Nguyen et al. (2012)). Two single-track bridges are analysed in this section using the 2D model from Section 2.2. The first, identified in what follows as Br-12m, corresponds to Old Guadiana River Bridge, which is a real bridge from a Spanish conventional railway line described in detail in Section 4. Old Guadiana River Bridge is made up of two identical SS spans of approximately 12 m in length. In each span the deck is composed of two adjacent structurally independent single-track decks. In this section only one of the decks is analysed neglecting any transverse coupling through the ballast between the decks. The second bridge is a variant from the former, with two identical spans of 20 m (Br-20m). In both cases, the track properties are considered identical. Table 2 includes the basic properties of the beam models that represent each of the bridges. f1and Mbi represent the fundamental frequency and the total mass per span, respectively. Br-20m bridge presents a similar static bending stiffness as Br-12m (Lb1/δUIC71 = 2250). The mass of Br-20m bridge is selected from the catalogue proposed 7 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper K(MN/m) pC(MNs/m) p Kouroussis (2015) Kouroussis (2011) Kouroussis (2011) Kouroussis (2015) K(MN/m) b Zhai(1996,2004) C(MNs/m) b (a) (b) (c) (d) Figure 5. Rail pad and ballast layer stiffness and damping values admitted by different authors in the past for comparable ballast layers thicknesses and sleeper distances (Bongini and Poisson (2009); Chen et al. (2017); Esveld (2001); Jes´us et al. (2014); Kourousis et al. (2011); Kouroussis et al. (2015); Lombaert et al. (2006); Naemi et al. (2015); Nguyen et al. (2012); Rigueiro et al. (2010); Sun and Dhanasekar (2002); Zhai (1996); Zhai et al. (2004)). Notation Parameter Value Unit Source ErRail UIC60 elastic modulus 2.100 ·1011 Pa CEN/TC256 (2017) Iyr Rail UIC60 moment of inertia 3055 ·10−8m4CEN/TC256 (2017) mrRail UIC60 mass per unit of length 60.335 kg/m CEN/TC256 (2017) KpRail pad vertical stiffness 1.000 ·108N/m Nguyen et al. (2012),Melis (2007) CpRail pad damping 7.500 ·104Ns/m Zhai et al. (2004) Msl Sleeper mass 300 kg Ministerio de Fomento (2010) Dsl Sleeper distance 0.600 m Ministerio de Fomento (2010) leHalf sleeper effective suporting length 0.950 m Zhai et al. (2004) lbSleeper width 0.300 m Ministerio de Fomento (2010) αBallast stress distribution angle 35 (◦) Zhai et al. (2004) hbBallast thickness 0.300 m Ministerio de Fomento (2010) ρbBallast density 1800 kg/m3Zhai et al. (2004) MbBallast vibrating mass 317.910 kg Zhai et al. (2004) EbBalast elastic modulus 1.100 ·108Pa Zhai et al. (2004) KbBalast vertical stiffness 1.933 ·108N/m Zhai et al. (2004) CbBallast damping 5.880 ·104Ns/m Zhai et al. (2004) EfSubgrade K30 modulus 9.000 ·107Pa/m Zhai et al. (2004) KfSubgrade vertical stiffness 7.399 ·107N/m Zhai et al. (2004) CfSubgrade damping 3.115 ·104Ns/m Zhai et al. (2004) KwBallast shear stiffness 7.840 ·107N/m Zhai et al. (2004) CwBallast shear damping 8.000 ·104Ns/m Zhai et al. (2004) Table 1. Track model parameters (per rail seat). by Dom´enech et al. (2014) for single-track bridges of a similar typology. Concerning the vertical stiffness of the neoprene bearings, in both bridges this value is selected so that the ratio between the flexural stiffness of the bridges and the equivalent vertical stiffness of the bearings for dynamic loads ¯ Kn bi,dyn, is in the range shown in Equation (1), which is a common interval in this type of structures and supports according to previous works (Museros et al. (2013)). A factor of 2 is 8 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper admitted between the static and dynamic vertical stiffnesses of the bearings. κ=EbiIybiπ3 ¯ Kn bi,dynL3 bi ∈[0.05 −0.1] (1) Figure 6 shows the first six natural frequencies and mode shapes calculated with the 2D trackbridge interaction model described in Section 2.2 for the two bridges of study (Br-12m and Br-20m). The number of sleepers is not exactly the same in the two spans. This causes a slight asymmetry in the mode shapes, which is more evident in the case of the shortest bridge. In what follows a sensitivity analysis is performed with the aim of evaluating the influence of individual variations of the track parameters on the bridge harmonic response (Section 3.2) and under the passage of trains (Section 3.3). Finally, the effect of including a different number of spans in the 2D model on the bridge maximum response and on the most critical section is also investigated (Section 3.4). Parameter Unit Br-12m Br-20m Nsp [ ] 2 2 Lbi m 11.93 20.00 f1Hz 10.07 6.67 Mbi kg 106659 210407 EbiIybi MNm27087 28800 ¯ Kn bi,st MN/m 1116.5 1488.3 Table 2. Bridges Br-12m and Br-20m properties. 3.2. Harmonic response In this section a vertical harmonic force with constant amplitude F0= 210 kN is applied on the rail at mid-span of the first span of the bridges. The maximum absolute vertical displacement at the same point is determined for exciting frequencies in the range ff∈[1,600] Hz in steps of ∆ff= 0.1 Hz. This analysis is repeated considering independent variations of the track parameters with respect to the reference ones (i.e. those included in Table 1). Figure 7 shows the results of the analysis for the Br-12m bridge. No structural damping is added, apart from that of the discrete track elements. The rail is discretised into two beam elements between consecutive sleepers, and so are the bridges. From the analysis of the results, the following is observed: •Two relevant maxima predominate in the harmonic response: a narrow peak corresponding to resonance of the fundamental mode of the bridge deck close to 10 Hz (see Figure 6(a-b)) and a wider maximum in the vecinity of 160 Hz coinciding with the first modes with a high contribution of the track deformation. Smaller peaks stand out along the response curves coinciding with the bridge higher frequencies with non-zero amplitudes at mid-span. •The effect of the subgrade damping (Cf) is negligible in the full frequency range. The effect of the rail pads and ballast dampings, Cpand Cb, affects the response at high frequencies but not under 60 Hz. •The only damping parameter that significantly affects the amplitude at low frequencies is the movilised ballast shear damping Cw, leading to a decrease of the response close to the bridge fundamental frequency, as shown in Figure 7(g). •As per the track stiffness properties, the variation of the subgrade stiffness Kfis negligible, which is expected as the ballast rests directly on the bridge slab and a very high value of Kb fis admitted. The increase in the rail pads and ballast stiffnesses, Kpand Kb, leads to a slight decrease of the amplitude at high frecuencies and to an increase of the track resonant response close to 160 Hz but, again, it does not affect the bridge response at low frequencies. 9 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper 0.00 0.20 0.40 0.60 0.80 1.00 1.20 0.1 0.3 0.5 0.7 0.9 1.2 1.4 1.6 1.8 2.1 2.3 2.5 2.7 2.9 0.00 0.20 0.40 0.60 0.80 1.00 1.20 0.1 0.3 0.5 0.7 0.9 1.2 1.4 1.6 1.8 0.00 0.20 0.40 0.60 0.80 1.00 1.20 0.1 0.3 0.5 0.7 0.9 1.2 1.4 1.6 1.8 2.1 2.3 2.5 2.7 2.9 0.00 0.20 0.40 0.60 0.80 1.00 1.20 0.1 0.3 0.5 0.7 0.9 1.2 1.4 1.6 1.8 x/Lb1x/Lb1 amax amax (a) Br-12m ( )N= 2 sp ES SS ( )c Brm ( )20 = 2Nsp ES SS ( )b Br-12m ( )N= 3 sp ( )d Brm ( )20 = 3Nsp Figure 11. ¯amax vs. bridge xcoordinate for Br-12m bridge with (a) 2 spans and (b) 3 spans; ¯amax vs. bridge xcoordinate for Br-20m bridge with (c) 2 spans and (d) 3 spans. 4. Ballast vertical coupling along the deck borders In this section, the vertical coupling induced by the shared ballast along the deck borders, i.e. between adjacent decks of the same span and between consecutive spans, is investigated. The study is centered on the real structure Old Guadiana River Bridge which consists of two identical SS spans with twin adjacent single-track decks. The deck width of this bridge is common in single track decks of Spanish conventional lines and the span length is representative of short bridges which may experience important acceleration levels under operating conditions. In what follows, the bridge and the experimental campaign carried out by the authors with the purpose of characterising the structure dynamic properties is described. Then, the 3D numerical model presented in Section 2.3 is updated to reproduce the experimental modal properties. Especial attention is paid to the vertical shear transmission through the ballast between decks and spans. Finally, the experimental response of the bridge under the passage of real trains is discussed and compared with the predictions of the calibrated numerical model. 4.1. Case study: Old Guadiana River Bridge The Old Guadiana River Bridge belongs to the conventional railway line Madrid-Alc´azar de San Juan-Ja´en, in Spain. It is a double-track concrete bridge with two identical simply-supported bays with 11.93 m of span length (Figure 12). The horizontal structure is formed by two adjacent but structurally independent decks, composed of a reinforced concrete slab resting on five prestressed concrete rectangular girders (see Figure 13). Each deck accommodates a ballasted track with Iberian gauge UIC60 rails and mono-block concrete sleepers at regular distances of 0.60 m. The longitudinal girders of the decks rest on the two abutments and on a central support through laminated rubber 16 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper bearings. Figure 12. Old Guadiana River Bridge. (a) 0.71 (b) Figure 13. Old Guadiana River Bridge: (a) Cross section of the deck; (b) Lateral view. In May 2019, the authors conducted an experimental campaign with the aim of characterising the dynamic response of the structure and the surrounding soil. For a detailed description of the experimental campaign the reader is referred to Galv´ın et al. (2021). Vertical acceleration was measured under ambient and train-induced vibration at 18 points of the lower horizontal face of the pre–stressed concrete girders (points 1-18 in Figure 14). Ambient vibration data recorded for 3600 s were used for the identification of the modal parameters of the bridge using state-space models. Additionally, modal damping ratios were also identified from the free vibrations left by the train passages, using 10 s of the acquired signals. As indicated in Galv´ın et al. (2021), the technique described by Kim et al. (2005) was used, which is based in the following three steps: (i) first, a digital band-pass filter is used to isolate each of the natural frequencies of interest from the free vibration time history signals. The pass bands of the digital filter are determined from the frequency content of the recorded accelerations. Special attention is paid to avoid overlapping when the natural frequencies are proximate; (ii) once the mode shapes are identified, single-degree of freedom responses that represent the response of the structure in each mode of vibration are obtained; (iii) the damping ratios are determined by fitting an exponential decay function to the single-degree of freedom responses. 17 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper MADRID JAÉN SPAN 1 DECK 1 DECK 2 DECK 1 DECK 2 SPAN 2 Track 1 Track 2 L/4 BB L/4 L/4 L/4 L/4 L/4 L/4 L/4 A18 A13 A12 A1 x y A4 A7 A14 A2 A5 A8 A15 A17 A3 A6 A9 A16 A10 A11 Figure 14. Old Guadiana River Bridge sensors layout. Table 5 shows the damping ratios (ζexp) and the natural frequencies of the five identified modes (fexp). The mode shapes are represented in Figure 15 in a solid black trace. The lowest one in frequency order corresponds to the first longitudinal bending of each bay where the two adjacent decks vibrate in phase. The second mode is the first torsion mode of the connected decks. In the third mode, the two adjacent decks deform under independent torsion but out of phase such that they conform a combined first transverse bending mode. In the aforementioned modes, the coupling due to the ballast layer connecting the adjacent decks is very evident. The fourth and fifth modes correspond, to an in-phase torsional deformation and to the transverse bending of each independent deck, respectively. Mode 1 2 3 4 5 fexp i[Hz] 9.84 11.03 12.84 21.43 28.74 ζexp i[%] 2.8 2.6 1.5 1.5 1.0 Table 3. First five frequencies and modal damping ratios identified experimentally. -1 6 4 2 0 -2 0 10 -4 5 0 -5 -6 -10 1 Mode 1 -1 6 4 2 0 -2 0 10 -4 5 0 -5 -6 -10 1 Mode 2 -1 6 4 2 0 -2 0 10 -4 5 0 -5 -6 -10 1 Mode 3 -1 6 4 2 0 -2 0 10 -4 5 0 -5 -6 -10 1 Mode 4 -1 6 4 2 0 -2 0 10 -4 5 0 -5 -6 -10 1 Mode 5 Figure 15. Identified mode shapes (black solid line) vs. calibrated (red solid line) on the Old Guadiana River Bridge. Undeformed shape (dashed gray line). 18 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper 4.2. Model updating from experimental measurements First, the 3D track-bridge interaction model is updated in order to reproduce the natural frequencies and mode shapes identified experimentally. This procedure is performed in two steps: (i) a preliminary sensitivity analysis to identify the most influencing parameters of the model, and (ii) the updating of these parameters through an automated iterative procedure for the minimisation of an objective function by means of a genetic algorithm. Table 4 lists the main parameters, including the initial values, ranges of variation, and the final (updated) values. The choice of the initial values, lower and upper bounds is based on engineering considerations, the level of uncertainty and on reference values found in the literature, cited in the last column. In the table, Eand ρ refer to the elastic modulus and mass density of the slab and the girders. Since no information was available on the material properties of the old Guadiana River Bridge, typical values for reinforced concrete have been adopted as initial guesses and a reasonable range of variation is admitted for the elastic modulus. The parameters Iyg and Jgstand for the second moment of area and the torsional constant of the girders cross-section, respectively. These properties have been estimated from the bridge technical drawings and are not included in the updating process. Concerning the track parameters, for the sake of brevity, in Table 4 only the parameters that may differ from those previously presented in Table 1 are listed: the ballast thickness underneath the sleepers hb, the vibrating ballast mass Mband immobilised mass mb, and the ballast shear stiffness along each rail Kwand between rails Kwr. For the remaining track parameters, the reference values listed in Table 1 are used and are not included in the updating process. Notation Initial Optimisation Calibrated Unit Source value range value Slab ρs2500 - 2500 [kg/m3] Malveiro et al. (2018) Es36000 [-20,+10]% 39390 [MPa] Malveiro et al. (2018) Girders Iyg 0.011 - 0.011 [m4] Jg0.005 - 0.005 [m4] ρg2500 - 2500 [kg/m3] Malveiro et al. (2018) Eg36000 [-20,+10]% 30127 [MPa] Malveiro et al. (2018) Track hb0.32 [ 0.3, 0.45] 0.32 [m] Ministerio de Fomento (2010) Mb317.9 349.256 [kg] Zhai et al. (2004) mb584.702 581.882 [kg/m2] Kw7.84 ·107[ -30, +30] % 5.49 ·107[N/m] Zhang et al. (2017) Kwr 7.84 ·107[ 0, +30] % 1.65 ·104[N/m] Coupling KwL 1.31 ·108[ 0, +30] % 2.76·107[N /m2] parameters KwT 1.31 ·108[ 0, +30] % 1.00·105[N /m2] Supports Kn bi,st 2.23·108- 2.23·108[N/m] Kn bi,dyn 2.45·108[ 1.1, 2]×Kn bi,st 4.46·108[N/m] Manterola (2006) Table 4. Mechanical properties of the 3D track-bridge interaction FE model. As shown in Section 3.3 the ballast shear stiffness along the track Kwmay significantly affect the bridge response, but information about realistic values for this parameter is scarce. Thus, a wide range of variation is admitted centered on the value proposed by Zhai et al. (2004). Kwr represents the ballast shear transmission in the transverse direction. For this variable, the lower bound is reduced to 0 as the gauge distance is almost four times the sleeper distance. The same initial values and ranges of variation are assumed for the ballast shear stiffnesses along the borders of the deck, i.e. between consecutive spans and adjacent decks, KwT and KwL, which are unknown and one of the main focus of this research. Since these discrete springs are distributed along the deck edges, the corresponding stiffnesses are given per unit length in Table 4, with the initial value set to Kw/Dsl = 7.84 ·107/0.6=1.31 ·108N/m2. Concerning the vertical stiffness of the neoprene bearings for static loads, an initial value Kn bi,st = 2.23 ·108N/m is calculated from a pre-design, as no information about the dimensions or mechanical properties of the bearings was available. Then, the stiffness under dynamic loading Kn bi,dyn is updated by multiplying the static 19 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper value by a factor in the range [1.1,2] according to Manterola (2006). Eight parameters are chosen as possible candidates for the updating process: Es,Eg,hb,Kw,Kwr,kwL,kwT and Kn bi,dyn. In order to avoid a poorly conditioned optimisation problem, a sensitivity analysis is performed to identify the parameters that most influence the prediction of the natural frequencies and mode shapes, in terms of frequency differences between the experimental (fexp i) and paired numerical modes (fnum i) (Equation (4)) and Modal Assurance Criterion (MAC) numbers (Equation (5)). ei=fexp i−fnum i fexp i (4) MACi=(Φexp i)TΦnum i2 ((Φexp i)TΦexp i) ((Φnum i)TΦnum i)(5) In the previous equation Φexp iand Φnum iare the i-th experimental mode and the paired numerical counterpart associated with mode i, and superscript T indicates transpose. The sensitivity analysis is performed using the Latin Hypercube sampling technique. 500 samples are generated within the ranges presented in Table 4 and the Spearman linear correlation matrix between the model parameters and the frequency differences and MAC numbers is computed. Figure 16 shows the correlation matrix for the first three modes. The coefficients +1 or -1 imply a perfect monotone correlation. In the graphs, coefficients between -0.25 and +0.25 are excluded from the graphical representation. EsEgKn bi;dyn hbKwKwr KwL KwT MAC1 MAC2 MAC3 e1 e2 e3 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 16. Spearman linear correlation matrix between model parameters and modal properties As can be seen, all the selected model parameters have a clear influence on the natural frequencies and mode shapes for the ranges of variations considered, the ballast height hbbeing comparatively the least influencing parameter. It is also noticeable that the ballast coupling parameters KwL and KwT are especially relevant for an accurate calibration of the first and second modes in terms of MAC numbers. As a second step in the updating process an automated optimisation procedure is implemented in ANSYS-MATLAB for the minimisation of the objective function Fobj, defined in Equation (6), applying a genetic algorithm. Based on the Spearman correlation matrix, seven model parameters are updated (Es,Eg,Kw,Kwr,KwL,KwT and Kn bi,dyn). The ballast height hb is assumed constant and set equal to its nominal value hb= 0.32 m. The ranges of variation of the parameters are again those given in Table 4. Fobj = n X i=1 |ei|+ n X i=1 (1 −MACi) (6) 20 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper The objective function Fobj involves both natural frequencies and MAC numbers for a total number of nidentified modes. In this study the first three modes are included in the optimisation range, due to their predominant contribution to the dynamic response of the bridge under operating conditions, as will be seen in Subsection 4.4. The values of the updated parameters after several runs of the genetic algorithm to guarantee the stability of the updating procedure are included in Table 4. In each run, a population size of of 70 individuals (10 times the number of calibration parameters) and 200 generations are considered. The crossover and mutation rates are set to 0.8 and 0.02, respectively, the tournament probability to 0.7, and the scale of mutation to 0.1. The numerical frequencies, frequency differences, and MAC numbers obtained with the updated model are provided in Table 5, for the first five identified modes. Figure 15 shows, in solid red trace, the paired numerical mode shapes. The correlation with the experimental measurements is accurate, especially for the first three modes. Parameter/Mode 1 2 3 4 5 fnum i[Hz] 9.83 10.82 13.02 19.48 27.89 fexp i[Hz] 9.84 11.03 12.84 21.43 28.74 ei[%] 0 1.92 -1.45 9.09 2.94 MAC [-] 0.96 0.93 0.98 0.95 0.78 Table 5. Numerical and experimental natural frequencies, frequency differences, and MAC numbers for the paired modes. 4.3. Effect of the ballast shear coupling along the deck borders on the modal parameters As indicated by the Spearman correlation matrix (Figure 16), the introduction of the ballast shear stiffness along the longitudinal joint (i.e., between adjacent decks of the same span), KwL, is essential to reach an accurate correlation of mode 2 (first torsion mode) in terms of MAC number. A similar coupling mechanism could have been expected between consecutive spans, however it is noticeable that the difference between the updated value of the parameters KwL and KwT is significant. Clearly, the vertical coupling exerted by the interlocking ballast granules along the deck borders is softer between consecutive spans. For a better understanding of the influence of KwL and KwT on the bridge natural frequencies and mode shapes, Figure 17 shows the results of a sensitivity analysis in which the aforementioned parameters vary independently with respect to the calibrated values from Table 4, in the range 0,1·104×KwL and 0,1·104×KwT . For each variation, the frequency differences eialong with the MAC values are calculated. A logarithmic scale is used in the horizontal axis of the plots for a better visualization of the results. In Figure 17(a-b) the results of the model calibration under variations of the ballast shear interaction between adjacent decks of the same span (KwL) are shown. The effect of KwL is significant for the calibration of the second and fourth identified modes, which are antisymmetric with respect to the longitudinal axis of the bridge and thus are more influenced by structural deformability in the transverse direction. For low values of KwL the numerical model is not able to properly reproduce the second and fourth modes, and the paired numerical modes exhibit low MAC numbers and high frequency differences (10% and 40%, respectively). The increase in the shear stiffness KwL improves the calibration of these modes in terms of both frequencies and MAC numbers, and there exists an optimum value of this parameter for best calibrated MAC and numerical frequencies. The calibrated value (1.0 ×KwL =2.76·107N/m2) provides accurate numerical frequencies with reasonably good MAC numbers for these two modes. For the first and third modes, a satisfactory mode pairing is achieved that is not influenced by the coupling between adjacent slabs. 21 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper -5 0 5 10 15 20 25 30 35 40 45 1.E-06 1.E-04 1.E-02 1.E+00 1.E+02 1.E+04 Mode 1 Mode 2 Mode 3 Mode 4 ei (%)  KwL (a) 0.6 0.7 0.8 0.9 1 1.E-06 1.E-04 1.E-02 1.E+00 1.E+02 1.E+04 MACi  KwL (b) -4 -2 0 2 4 6 8 10 1.E-06 1.E-04 1.E-02 1.E+00 1.E+02 1.E+04 Mode 1 Mode 2 Mode 3 Mode 4  KwT (c) ei (%) 0.6 0.7 0.8 0.9 1 1.E-06 1.E-04 1.E-02 1.E+00 1.E+02 1.E+04  KwT (d) MACi Figure 17. Sensitivity analysis on the influence of the ballast shear coupling parameters KwL and KwT on the natural frequencies and MAC values. Figures 17(c-d) show the isolated contribution of KwT on the model calibration. It is noticeable that the increase in KwT worsens the MAC numbers for the paired mode associated to the first and second identified eigenforms after a certain threshold value, while its influence on the natural frequencies is minor. The influence of the vertical shear transfer mechanism due to the shared ballast layer between independent spans is less relevant than between adjacent decks. In any case, it is worth noting that the discrete ballast representation used in this work presents limitations and, in view of the significant dynamic coupling between the decks exhibited by this bridge, it would be of interest to continue investigating on the dynamic behaviour of this type of bridges with structurally independent parts coupled through the ballast using continuous models. 4.4. Dynamic response under operating conditions During the experimental campaign, the response of the structure was recorded under the circulation of several passenger trains. The performance of the bridge was adequate from a serviceability perspective, as the maximum acceleration did not exceed the Serviceability Limit State for traffic safety in ballasted tracks, limited to 3.5 m/s2per European Standards (CEN/TC250 (2005)). Two of the trains are included in this section, corresponding to the medium distance Renfe trains Altaria Talgo VI and S449. Both trains crossed the bridge along track 2 (see Figure 14) northbound (Manzanares-Alc´azar de San Juan). Figure 18 shows photographs of the trains and the axle scheme. Table 6 includes the arrangement of the axles and loads for each train. Train N dk[m] l1[m] l2[m] P1[kN] P2[kN] P3[kN] Altaria 7 13.14 3.44 3.3 225 70 140 S449 3 17.75 14.8 – 161 – – Table 6. Features of the Renfe medium distance Altaria Talgo VI and S449 trains. The circulation speeds were identified from the frequency associated with the bogie distance as approximately 155 km/h in both cases. According to the distance between shared axles or bogies in the passenger cars, dk, this speed is close to a third resonant speed of the fundamental mode for 22 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper 4 x P1P2P2 P3P3P3 dd Nxd 3 l173l1l2 Altaria Talgo VI 2 x P12 x P12 x P12 x P1 2.5 l12.7 Nxd 4.21 S449 2 x P1 d Figure 18. Renfe Altaria Talgo VI (left) and S449 (right) trains crossing Old Guadiana River Bridge Altaria and to a fourth resonant speed of the same mode for the S449 train (Equations (7a) and (7b), respectively). Additionally, the speed of the S449 train is also in the vicinity of the theoretical critical speed for a fifth resonance of the third identified mode (Equation (7c)). Vr,Altaria 1,3=13.14 m ·f1·3.6 3= 155.16 km/h (7a) Vr,S449 1,4=17.75 m ·f1·3.6 4= 157.19 km/h (7b) Vr,S449 3,5=17.75 m ·f3·3.6 5= 164.01 km/h (7c) The response of the bridge under the circulation of the aforementioned trains is now calculated using the numerical model updated in the previous section. The dynamic equations of motion are solved by mode superposition. Modal damping values identified during the experimental campaign under free vibration after train passages are assigned to the paired numerical modes (see Table 5). Figure 19 shows an experimental vs. numerical comparison of the vertical acceleration at sensors 2, 5, 18 and 17 under the passage of Altaria Talgo VI train in the time (a-d) and frequency (e-h) domains. The sensors are located at mid-span under the loaded deck (a, e, b, f) and under the adjacent unloaded deck (c, g, d, h) of the same span. In all the plots the experimental signal, plotted in solid black trace, is filtered applying two third-order Chebyshev filters with high-pass and low-pass frequencies of 1 Hz and 30 Hz, respectively. The numerical predictions, plotted in gray trace, include modal contributions up to 30 Hz. As can be depicted from the frequency domain plots, the vertical acceleration of the bridge is affected by several modal contributions, which can be attributable to the span length–to–width ratio that causes the bridge behaviour to differ from a beam-type structure. However, the peak amplitude associated to the fundamental mode predominates when compared to other modal responses, as expected in a resonance situation. This is especially evident at points 5 and 17, located at mid-span under the loaded and unloaded tracks, respectively. The response also presents visible peaks at low frequencies (in the range 3 Hz - 6 Hz) associated to the excitation and coinciding with the axle passing frequency and its multiples, which are more visible at the sensors installed under the loaded deck (points 2 and 5). For sensors 17 and 18, located at the unloaded deck, the response is still relevant (the maximum acceleration at sensor 17 reaches approximately 35% of the peak amplitude in the analogous loaded sensor number 5). This is clearly a consequence of the vibration transmission through the ballast between the two independent decks, which is evident 23 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper (a) 10 12 14 Time [s] -1 -0.5 0 0.5 1 Acceleration [m/s2] Altaria Talgo VI track 2 v= 155 km/h (e) 0 10 20 30 Frequency [Hz] 0 0.5 1 Acceleration [m/s2/Hz] Altaria Talgo VI track 2 v= 155 km/h 02z exp 02z num (b) 10 12 14 Time [s] -1 -0.5 0 0.5 1 Acceleration [m/s2] Altaria Talgo VI track 2 v= 155 km/h (f) 0 10 20 30 Frequency [Hz] 0 0.5 1 Acceleration [m/s2/Hz] Altaria Talgo VI track 2 v= 155 km/h 05z exp 05z num (c) 10 12 14 Time [s] -1 -0.5 0 0.5 1 Acceleration [m/s2] Altaria Talgo VI track 2 v= 155 km/h (g) 0 10 20 30 Frequency [Hz] 0 0.5 1 Acceleration [m/s2/Hz] Altaria Talgo VI track 2 v= 155 km/h 18z exp 18z num (d) 10 12 14 Time [s] -1 -0.5 0 0.5 1 Acceleration [m/s2] Altaria Talgo VI track 2 v= 155 km/h (h) 0 10 20 30 Frequency [Hz] 0 0.5 1 Acceleration [m/s2/Hz] Altaria Talgo VI track 2 v= 155 km/h 17z exp 17z num Figure 19. Old Guadiana Bridge: (a-d) time history and (e-h) frequency content of the acceleration at points 2, 5, 18 and 17 induced by Renfe Altaria train circulating along track 2. Experimental results (black line) vs. numerical predictions (gray line) both under ambient and under forced vibrations, despite the difference in amplitude between these two situations. The numerical predictions seem reasonable, especially for the sensor located at mid-span in the loaded deck (sensor 5) and in the frequency range of interest. The contribution of mode 1 prevails in all the sensors located under the loaded track (sensors 5, 4 and 6, the latter not shown in the plots) and is to some extent overestimated by the numerical model. This could be attributable to the interaction between the bridge and the vehicle suspension systems, not included in the model and more relevant in a resonance situation. The response at the same locations but measured under the unloaded track, sensors 17 and 13 (the latter not visible in Figure 19), reveals also 24 March 2, 2022 Structure and Infrastructure Engineering SIE˙paper certain overestimation of the acceleration content associated to the fundamental mode. This is more noticeable during the forced-vibration phase, according to the corresponding time-domain plots (Figure 19(d)). Conversely, the acceleration amplitude is accurately predicted in the unloaded deck during the free-vibration phase, entailing an accurate identification of structural damping for comparable amplitudes of vibration. This difference in the numerical and experimental adjustment between the forced and free vibration phases may be related not only to VBI but also to amplitudedependent damping mechanisms and other sources of nonproportional local damping, such as the ballast shear damping along the track and between the ballast-connected structures. For the sensors located at the external borders of both decks (sensors 2 and 18) the participation of several mode contributions below 30 Hz is significant, especially that of the second (first torsion) and third (first transverse bending) modes. This stresses the need to simulate the coupling exerted by the track in these types of structures. Nevertheless, the contribution of the fundamental mode on these sensors responses prevails, as a consequence of the first mode resonance induced by Altaria Talgo. 25