P Richard

Sophie Mongabure
email: philippe.mongabure@cea.fr

F Capdevielle

Ragueneau

C Lejouad

B Richard

P Mongabure

S Capdevielle
email: sophie.capdevielle@univ-grenoble-alpes.fr

F Ragueneau
email: frederic.ragueneau@ens-paris-saclay.fr

Experimental study of corroded RC beams: dissipation and equivalent viscous damping ratio identification

Keywords: Steel corrosion, shaking table, reinforced concrete, dissipation, damping

scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

Corrosion of steel reinforcement is one of the most widespread pathologies that leads to a loss of structural performance of reinforced concrete (RC) components. This statement is particularly true for RC structures located in marine areas with exposure to chloride ions and humidity.

At the material scale, the corrosion consequences are numerous: reduction of the steel reinforcement cross-section and ductility [START_REF] Almusallam | Effect of degree of corrosion on the properties of reinforcing steel bars[END_REF], concrete cover cracking [START_REF] Andrade | Cover cracking as a function of bar corrosion: Part I-Experimental test[END_REF] and bond strength variation between steel and concrete [START_REF] Almusallam | Effect of reinforcement corrosion on bond strength[END_REF].

Corrosion degradation leads to structural modifications at the member scale. In order to assess these effects, several experimental campaigns have been conducted. These experimental programs considered specimens of various types and dimensions, from small scale RC beams [START_REF] Mancini | Local reinforcing bar damage in r.c. members due to accelerated corrosion and loading[END_REF] to large scale RC columns [START_REF] Meda | Experimental evaluation of the corrosion influence on the cyclic behaviour of RC columns[END_REF] [6] [START_REF] Malumbela | Behaviour of RC beams corroded under sustained service loads[END_REF] [START_REF] Rinaldi | Experimental evaluation of the flexural behavior of corroded P/C beams[END_REF].

The experimental campaigns carried out so far have investigated the effect of steel reinforcement corrosion on the quasi-static behavior of RC structures using either monotonous [START_REF] Zhu | Corrosion of the reinforcement and its influence on the residual structural performance of a 26-year-old corroded RC beam[END_REF] [START_REF] Malumbela | Behaviour of RC beams corroded under sustained service loads[END_REF] or cyclic loadings [START_REF] Guo | Experimental investigation on the cyclic performance of reinforced concrete piers with chloride-induced corrosion in marine environment[END_REF] [START_REF] Ma | Behavior of corrosion damaged circular reinforced concrete columns under cyclic loading[END_REF]. Based on the experimental results, a decrease of the static properties, namely the bearing capacity, the ductility and the stiffness with an increasing corrosion ratio has been observed.

Several coastal regions are located in seismic zones such as in Italy [START_REF] Casolo | Comparative seismic vulnerability analysis on ten masonry towers in the coastal Po Valley in Italy[END_REF]. For this reason, the study of the dynamic behavior of corroded RC structures is of crucial importance. The energy dissipation capacity is an important indicator when studying the seismic behavior of RC structures. It quantifies the ability of a structure to accommodate the seismic demand and therefore, to dissipate energy without compromising the structural stability, especially regarding safety. The equivalent viscous damping ratio (EVDR), which characterizes the dissipation due to the non-linear (hysteretic) behavior, provides knowledge about the seismic response of RC elements. In common practice, the EVDR is deduced from the quasi-static cyclic response [START_REF] Guo | Experimental investigation on the cyclic performance of reinforced concrete piers with chloride-induced corrosion in marine environment[END_REF] [11] using Jacobsen's area method [START_REF] Jacobsen | Damping in composite structures[END_REF].

Few experimental campaigns using dynamic loadings have been performed. Among them, an experimental campaign was carried out by [START_REF] Zou | Experimental investigation on the dynamic properties of RC structures affected by the reinforcement corrosion[END_REF] on corroded RC frames in order to assess the influence of corrosion on the modal properties of RC structures. Using hammer shock tests, an increase in the modal damping ratio is observed for an increasing corrosion degree. However, at a critical corrosion degree, this quantity of interest tends to decrease. Contradictory results were obtained in another study [START_REF] Razak | The effect of corrosion on the natural frequency and modal damping of reinforced concrete beams[END_REF], in which the authors observed a decrease of the damping ratio for the corroded specimens. Based upon this observation, there is no consensus in the literature on the effect of corrosion on the damping capability of RC structures and additional studies should be performed to clarify this issue. In [START_REF] Yuan | Shaking table tests of coastal bridge piers with different levels of corrosion damage caused by chloride penetration[END_REF], shaking table tests were performed on corroded bridge piers. An increase in the damping ratio have been observed as the corrosion gravity gets worst. From this brief literature review, it is indicated that there are no universal agreements on the effect of reinforcement corrosion on the damping capacity of RC structures. Hence, more experimental investigations are required to address the inconsistency between these different results. Some studies on RC structures have identified these important parameters relying on both the quasistatic cyclic and the dynamic response [START_REF] Zou | Experimental investigation on the dynamic properties of RC structures affected by the reinforcement corrosion[END_REF] keeping the same experimental setup. The non-parametric identification method is [START_REF] Heitz | Dissipations in reinforced concrete components: Static and dynamic experimental identification strategy[END_REF] used in this type of studies to evaluate the dissipation capacity. According to the authors' knowledge, no similar experimental study has been performed on corroded RC structures.

Experimental campaigns using cyclic loadings are a significant means to predict the dynamic response of structures. However the inertial effect is not taken into account in such experiments. Furthermore, if viscous effects are introduced, they cannot be considered in such type of tests. For these reasons, an experimental campaign has been carried out on 18 large-scale corroded RC beams, in addition to two non-corroded RC beams. The corroded beams are subjected to a cyclic four-point bending loading as well as to dynamic loadings using a shaking table and using the same experimental setup. This study aims to:

(1) provide numerical models with some input data regarding the dissipation capacity of corroded RC components; (2) compare the trends of specific structural properties (eigenfrequency, damping ratio) identified from quasi-static and dynamic responses; (3) compare different identification methods to estimate damping based on quasi-static and dynamic experimental results. This paper will first give an overview of the experimental campaign carried out. Then, the experimental results in terms of dissipation and damping capacity, obtained from quasi-static as well as dynamic testing, are underlined. Finally, the concordance between the results of dynamic and quasi-static testing is evaluated.

Experimental campaign

Scope

The experimental campaign DYSBAC, a French acronym for "Dynamic behavior of corroded RC structures", is performed by means of the AZALEE shaking table and the strong floor, which are parts of the TAMARIS experimental facility operated by the French Alternative Energies and Atomic Energy Commission (CEA) located in Saclay, France. The main objective of this experimental campaign is to study the influence of corrosion on the dynamic behavior of RC beams. In particular, the energy dissipation and the damping capacities are evaluated based on the results acquired from both quasi-static and dynamic tests. Thereafter, the results are cross-checked in order to assess the consistency between the two types of identifications.

Specimen details

In order to be representative of real RC structures and overcome any scale effect, the choice of largescale RC beams has been made. Each beam is 4. The 20 specimens are cast with a C25/30 concrete class according to the Eurocodes 2 [START_REF] Mosley | Reinforced concrete design: to Eurocode 2[END_REF]. The considered concrete is formulated with a high water cement ratio of 0.59. This leads to a porous and low strength concrete, representative of concrete in aged RC structures.

Accelerated corrosion of the specimens

Since using natural corrosion process to obtain corroded specimens is highly time-demanding [START_REF] Zhu | Corrosion of the reinforcement and its influence on the residual structural performance of a 26-year-old corroded RC beam[END_REF], some technics are used in order to reduce the required time exposure [21] [22] [23] [START_REF] Vidal | Analyzing crack width to predict corrosion in reinforced concrete[END_REF] [25] by accelerating the corrosion phenomenon. The differences between each accelerated corrosion technique and the natural corrosion process have been widely discussed in the literature [23] [26]. The imposed current technique is the most commonly used when studying the structural effect of steel reinforcement corrosion. However, many studies have shown that the choice of the current density value has a direct consequence on the similitude with natural corrosion. Especially, the crack patterns or the chemical nature of rust products is affected. However, most of the existing studies suggest not to exceed 100 μA.cm -2 to have similar structural effects to natural corrosion [2] [27].

The imposed current technique consists in applying an electrical current from a Direct Current (DC) power supply between the cathode, which is a stainless steel grid, and the anode, which is the reinforcement inside the RC specimen [START_REF] El Maaddawy | Effectiveness of impressed current technique to simulate corrosion of steel reinforcement in concrete[END_REF]. The whole specimen is immersed in an electrolytic solution containing chlorides in order to ensure electrical conduction and to be representative of corrosion by chlorides.

One of the objectives of this study is to separately assess the effect of corrosion of the different reinforcement parts (longitudinal and transversal reinforcements) on the mechanical response of the specimens. For this reason, three beam configurations are considered:

• 𝐶 1 for the longitudinal reinforcement corrosion;

• 𝐶 2 for the stirrups corrosion;

• 𝐶 3 for the full reinforcement corrosion.

In order to reach the corrosion targets for the three configurations, different parts of reinforcement were electrically insulated and different cathode settings were adopted, depending on the configuration.

As recommended in [START_REF] Caré | Influence of impressed current on the initiation of damage in reinforced mortar due to corrosion of embedded steel[END_REF], the current density was limited to 100 𝜇𝐴. 𝑐𝑚 -2 in order to obtain rust products as similar as possible to those formed in natural conditions. Despite the fact that the chemical nature of rust product is not strictly equivalent to the one observed in case of natural corrosion, the mechanical consequences of rust formation are the same [START_REF] Ožbolt | Accelerated Corrosion of Steel Reinforcement in Concrete: Experimental Tests and Numerical 3D FE Analysis[END_REF]. Three corrosion rates, expressed in terms of mass losses are targeted: • 5% , because many studies demonstrated a degradation of the bond between steel and concrete for a corrosion rate between 1.5% [START_REF] Coccia | Influence of corrosion on the bond strength of steel rebars in concrete[END_REF] and 5% [START_REF] Auyeung | Bond behavior of corroded reinforcement bars[END_REF] [3]; • 10% , rate from which civil engineering maintenance operations start [START_REF]Projet APPLET : Durée de vie des ouvrages en béton : Approches prédictives performantielles et probabilistes[END_REF]; • 15% , demonstrated in some studies to be the threshold from which a change of the steel failure mechanism is observed [1] [33].

All the RC beams are immersed in a 3.5% NaCl solution. The exposure time is estimated for each type of beam and each corrosion rate using Faraday's law (Equation 1) with 𝛼 = 1.3 [START_REF] Meda | Experimental evaluation of the corrosion influence on the cyclic behaviour of RC columns[END_REF]. The required exposure duration varies from 31 days to 141 days, depending on the targeted corrosion rate and on the beam configuration.

∆𝑡 =

𝛼 .∆𝑤 .𝑧.𝐹 𝑀.𝐼

In Equation (1), Δ𝑤 is the mass of steel consumed due to corrosion (kg.m -2 ), 𝐼 is the current density (A.m -2 ), 𝛥𝑡 is the exposure time (s), 𝐹 is the Faraday constant 96 500 (C.mol -1 ), 𝑧 is the ionic charge (2 for Fe), 𝑀 is the atomic weight of steel (g.mol -1 ), 𝛼 is a coefficient usually taken between 1 and 2 to take into account the duration of chloride ingress into concrete before reaching the rebar.

Since the actual corrosion degree is difficult to measure experimentally on the tested beams due to their large dimensions as well as the high reinforcement density, 9 dedicated RC specimens were cast. The same concrete formulation and steel reinforcement as the real tested beams are used. The specimens were corroded at 5%, 10% and 15% using the same setup previously described. The rebars were removed from the concrete and mechanically cleaned afterwards. Table 1 summarizes the mean weighing results. Since the chloride RC corrosion is not uniform, the mass loss ratio is not a sufficient indicator to characterize the corrosion state. For this reason, using the same bars dedicated to the mass loss measurement, the distribution of the steel cross section diameter has been measured. The measurements have been realized by the means of a profilometer bench with a running transverse laser beam along the bar.

As long as the bending behavior is studied in this experimental campaign, only the longitudinal rebars are inspected. Table 2 summarizes statistical indicators coming from the measured diameter distributions. For the sake of clarity, each beam is labeled hereafter with respect to its corrosion configuration (𝐶 1 , 𝐶 2 , or 𝐶 3 ) followed by the targeted corrosion rate (5%, 10% or 15%). To be more specific, the label Cn_R refers to the configuration n and to the corrosion rate R. For example, the label C1_15 corresponds to a longitudinal corrosion beam corroded at 15%.

Experimental setup and test program

Experimental testing consists of a quasi-static as well as a dynamic characterization. 9 corroded specimens (3 corrosion rates and 3 corrosion configurations) and a unique reference beam are tested by means of the strong floor. The dynamic testing of the 9 other corroded beams (3 corrosion rates and 3 corrosion configurations) in addition to another virgin beam is performed on the AZALEE shaking table.

The setup is similar to the one used for the IDEFIX campaign [START_REF] Heitz | Dissipations in reinforced concrete components: Static and dynamic experimental identification strategy[END_REF]. Some modifications have been made to adapt the experimental setup for accommodating high range displacements and rotations.

The beams are excited for both dynamic and quasi-static tests along their weakest flexural axis; the boundary conditions are the followings:

• rotating supports allowing the rotation at the beam extremities around the vertical axis. These supports are designed to ensure a maximum rotation angle of ±24° at the beam boundaries. This choice is consistent with the expected maximal displacement imposed at the loading points of ±400 𝑚𝑚; • two additional masses, of 94 𝑘𝑔 weight each, fixed at the intermediate supports. These masses ensure a first natural frequency of the beams in the range of 12 and 13 𝐻𝑧. The total mass of each intermediate support including the two additional masses is around 310 kg; • two air-cushion systems to bear the beam weight and to drastically reduce the friction between the beam and the shaking table's or strong floor's upper plate. These measures are taken to prevent the beam from cracking under the addition of its self-weight and the lumped masses weight.

Quasi-static testing

Regarding the quasi-static tests, the corroded and non-corroded beams are subjected to a classical fourpoint reverse cyclic bending test on the TAMARIS strong floor. The loading is applied by means of a long-stroke actuator. This actuator is linked with a reinforced metal beam able, through swivels at its ends, to distribute the loading on two points of the DYSBAC beam. The actuator of 100 𝐾𝑁 capacity has a maximum displacement of ± 500 𝑚𝑚 and a maximum velocity of 1.7 𝑚. 𝑠 -1 . A general view on the experimental setup is shown in Figure 1. The choice of four-point rather than three-point bending test is justified by the desire to have a distributed instead of a localized damage. This choice does not necessarily constrain the beam failure to occur at midspan. The applied loading is composed of blocks of three identical cycles during which the displacement is controlled, with an increasing amplitude between two consecutive blocks (Figure 2). The displacement speed is controlled at 0.5 mm/s. Each cycle involves 4 phases: loading in one direction, unloading, loading in the other direction and unloading. Unloading is performed up to recovering a zero force measured with the actuator load cell. The goal of this loading protocol is to have 3 identical cycles to stabilize the new damage level associated with the current block before moving on to the next one. The amplitude range varies from 0.4 mm up to 200 mm. Hammer shock tests are also performed between the blocks in order to get the evolution of the modal properties with damage. Two impact positions are considered: at midspan and near an intermediate support. The modal properties presented in this paper are the mean values of the modal properties determined with each hammer test using the measured acceleration at the different sensors positions. In order to decrease the intrinsic discrepancy related to these tests, they have been performed by the same operator all along the experimental campaign.

Dynamic testing

The dynamic tests are performed on the AZALEE shaking table. It is a 6 × 6 𝑚² shaking table able to reproduce seismic signals up to 3 𝑔 depending on the payload. The table is controlled on the 6 degrees of freedom (3 rotations, 3 translations) (Figure 3). The dynamic loading consists in a synthetic signal able to excite only the first natural mode of the beam. This bandlimited signal has been chosen in order to have a quasi-constant power spectral density over the entire frequency domain of interest [2Hz, 13 Hz] (according to the displacement limits of the shaking table). This choice has been made in order to take into account the frequency shift due to the damage growth during the test. 5 peak acceleration levels are used: 0.125 𝑔, 0.5 𝑔, 0.8 𝑔, 1.25 𝑔 and 2 𝑔. A modal characterization of each beam is performed before testing and between two consecutive testing sequences using a white noise (WN) signal (𝑃𝐺𝐴 = 0.1 𝑔).

Measurements

In order to fully characterize the mechanical response of the specimens during the tests, different types of sensors are used:

• 5 displacement wire sensors;

• 2 high precision displacement wire sensors;

• 2 six-axis load cells at the beam ends;

• 7 three-axis accelerometers;

• Actuator load cell.

In addition, the digital image correlation (DIC) technique is used. It consists of black and white strips painted all along the upper surface of the beam. The whole beam is filmed from above using a stereovision rig. In this way, the deformed shape of the beam during the tests can be post processed after testing.

The damping capability identification

The equivalent damping is an interesting feature to study when assessing the seismic behavior of RC structures. It measures the capability of a structure to dissipate the input energy. This parameter depends on the structure's ductility demand as well as on the plastic hinges position [START_REF] Priestley | Displacement-based seismic assessment of reinforced concrete buildings[END_REF]. The EVDR ( 𝜉 𝑒𝑞 ) is calculated as the sum of an elastic damping 𝜉 𝑒𝑙𝑎𝑠 and an hysteretic one 𝜉 ℎ𝑦𝑠𝑡 [START_REF] Smyrou | Modelling of elastic damping in nonlinear time-history analyses of cantilever RC walls[END_REF] (see Equation 2).

𝜉 𝑒𝑞 = 𝜉 𝑒𝑙𝑎𝑠 + 𝜉 ℎ𝑦𝑠𝑡 (2)
𝜉 𝑒𝑙𝑎𝑠 accounts for the damping when the RC element remains in the elastic domain. It is commonly taken equal to 2% [START_REF] Richard | SMART 2008: Overview, synthesis and lessons learned from the International Benchmark[END_REF]. 𝜉 ℎ𝑦𝑠𝑡 accounts for the damping induced by hysteretic phenomena.

In a predictive numerical context, the EVDR is calibrated in order to dissipate the right amount of energy by a viscous force field, acting proportionally to the velocity field. However, a more classical approach consists of assigning a damping ratio value to each eigenmode of the structure.

In this part, several identification methods are used to estimate the specimens damping. These methods can be divided into two categories:

-methods using the structure response subjected to low level loads (hammer shock excitation or low level white noise signal). The identification of the modal properties using low level loads is of primary importance when dealing with dynamic monitoring of existing structures. In this study, the covariance-driven stochastic subspace identification (SSI) method [START_REF] Li | Modal identification of civil structures via covariance-driven stochastic subspace method[END_REF] is used for both hammer shock and white noise tests. The main objective is to compare the trends revealed by both types of testing; -methods based on the structure response subjected to high level loads (quasi-static or dynamic).

This type of identification methods is useful for numerical model calibration. For quasi-static loads, Jacobsen's area method [START_REF] Jacobsen | Damping in composite structures[END_REF] is used to evaluate dissipated and stored energy and therefore the equivalent viscous damping ratio. For dynamic loads, the non-parametric identification method [START_REF] Heitz | Dissipations in reinforced concrete components: Static and dynamic experimental identification strategy[END_REF] applied on the structure displacement response can provide a valuable information about the evolution of modal properties (natural frequency and damping ratio) along time depending on the damage sustained by the structure along the loading.

Since the specimen exhibits a behavior that mostly includes bending and only the first eigenmode of the beam is excited during the dynamic testing, the focus is made on 1 st natural mode exclusively. However, the acquired data allow a similar analysis for the 2 nd eigenmode which can be performed if needed.

It is to be noted that the C3_15% beam is excluded from the results discussion. This is due to an incident that occurred during the corrosion process.

Covariance-driven stochastic subspace identification method

The experimental modal analysis is a powerful tool for describing, understanding and modeling the dynamic behavior of a structure. It consists of determining the natural frequencies, mode shapes and modal damping ratios based on the measured experimental data.

The covariance-driven stochastic subspace identification (SSI) method is an effective instrument for the experimental modal identification driven by output-only records [START_REF] Pappa | Some experiences with the eigensystem realization algorithm[END_REF] [39] [START_REF] Liu | Modal parameter estimation using the state space method[END_REF]. It is of a great interest when the input information is unknown (e.g., traffic load or wind load). Thus, the excitation is assumed to be a realization of a stochastic process. In [START_REF] Li | Modal identification of civil structures via covariance-driven stochastic subspace method[END_REF], a detailed overview on the covariance-driven SSI method basic theory is given.

Quasi-static testing condition

Thanks to the hammer shock tests performed during the quasi-static campaign after each loading block, the modal damping ratio can be identified. The objective of such characterization is to assess the evolution of the damping capability as a function of the damage level experienced by the structure. In addition, the influence of the corrosion ratio on the damping capability is investigated. The covariancedriven SSI method is used (see section 3.1) to identify the corresponding modal damping ratio for each RC beam subjected to hammer shock tests.

Table 3 summarizes the 1 st natural frequency as well as the 1 st damping ratio of each specimen prior to any load testing. Regarding the effect of corrosion on the 1 st natural frequency, no particular trend could be noted. It is generally known that natural frequencies have low sensitivity to damage [41] [42].

On the contrary, an increase in the 1 st damping ratio is basically observed as the corrosion level gets higher, especially for the 𝐶 1 and 𝐶 3 configuration beams. This observation is consistent with a previous study reported in the literature [START_REF] Zou | Experimental investigation on the dynamic properties of RC structures affected by the reinforcement corrosion[END_REF].

However, another previous study dealing with the effect of general corrosion on the modal properties of RC elements, has shown an opposite trend of the 1 st damping ratio [START_REF] Razak | The effect of corrosion on the natural frequency and modal damping of reinforced concrete beams[END_REF]. The inconsistency between the two studies may be explained by the fact that the corrosion degradation type is different, leading to non-similar structural effects. In the present study, the localized corrosion induces a local loss of bond strength, which is not the case for generalized corrosion where the degradation is more global. In generalized corrosion, the steel cross-section loss is uniform and thus remains small even if the total corrosion ratio is high. In contrast, the localized corrosion is characterized by a sharp local reduction in the rebar cross-section. The coupling between these two effects of localized corrosion, which are the localized bond strength loss and the localized corrosion pits, is studied in [START_REF] Castel | Mechanical behaviour of corroded reinforced concrete beams-Part 2: Bond and notch effects[END_REF]. Figure 4 depicts the relationship between the 1 st modal damping ratio and the ductility level for each beam configuration and corrosion rate. To quantify the ductility level, the following index is introduced:

𝜇 = max|𝑢(𝑡)| 𝑢 𝑦 (3)
where 𝑚𝑎𝑥|𝑢(𝑡)| is the maximum of the absolute value of the relative displacement measured at midspan with respect to the supports during the loading block prior to the hammer shock test, and 𝑢 𝑦 is the midspan displacement corresponding to the first yield occurrence in the steel rebars for the concerned beam.

The procedure adopted in this work for the evaluation of 𝑢 𝑦 is the Park's method [START_REF] Park | Evaluation of ductility of structures and structural assemblages from laboratory testing[END_REF]. First, it consists of computing both secant and post-yield stiffness, and second determining their intersection points. The yield displacement is computed as the average of the two values obtained at the intersecting points. The main advantage of this method is the possibility to take into account the asymmetries observed in the two loading directions.

For the 𝐶 1 configuration, the corroded beams exhibit a smaller damping capability compared to the reference beam, as shown in Figure 4. a. In particular, when the steel yielding occurs, which corresponds to a ductility level equal to 1, the 1 st damping ratio decreases with an increasing corrosion degree. This trend can be explained by the fact that the bonding between steel and concrete in enhanced for lightly corroded specimens [START_REF] Almusallam | Effect of reinforcement corrosion on bond strength[END_REF].

In Figure 4. b, it can be noticed that 𝐶 2 corroded beams curves give similar trend regarding the 1 st modal damping ratio evolution, regardless of the degree of corrosion. Moreover, a drop-off in the damping capability is observed at each ductility level for each corroded beam compared to the virgin one.

It can be noticed that the stirrups corrosion induces the highest damping decrease compared to the other corrosion configurations. Furthermore, the similarity between the values of the 1 st modal damping ratio for each 𝐶 2 beam at each ductility level confirms that for low corrosion ratios the bond between steel and concrete is enhanced [3] [45].

Regarding the 𝐶 3 corroded beams, a decrease in the 1 st modal damping ratio is observed for the C3_5%, C3_10% at each ductility level after yielding. It is important to note that for the 𝐶 3 corroded beams, the corrosion process using the imposed current induces a heterogeneous corrosion state at the longitudinal reinforcement and stirrups level. This is because only one DC power supply is used. Therefore, the distribution of the electrical current depends on the resistivity of the different reinforcement parts. Based on the values displayed in Figure 4. c, it could be conjectured that for C3_10% beam the mainly corroded part is the longitudinal rebars, whereas it is stirrups for C3_5% beam.

For all the tested beams, a change in the shape of the curves is observed near the yielding point (ductility = 1). The 1 st damping ratio increases first, then slightly decreases in the vicinity of the yielding point till the loading ends. This observation is in good agreement with the fact that the modal damping ratio increases with an increasing crack depth, as demonstrated in [START_REF] Panteliou | Damping factor as an indicator of crack severity[END_REF].

Dynamic testing condition

During the dynamic campaign, the WN applied on each tested beam after each testing sequence made possible the identification of the modal properties using the covariance-driven SSI method introduced in section 3.1. The objectives are the same as those targeted in section 3.1.1.

The 1 st natural frequency as well as 1 st natural damping ratio identified with white noise modal testing of each beam before being damaged by the dynamic loading is presented Table 4. Based on these results, an increase in the modal damping ratio with an increasing corrosion degree is observed. These results are consistent with those issued from hammer shock tests in section 3.1.1 and with a recent study reported in the literature [START_REF] Yuan | Shaking table tests of coastal bridge piers with different levels of corrosion damage caused by chloride penetration[END_REF]. As far as the 1 st frequency is concerned, no specific effect of corrosion could be observed. These results are consistent with those coming from hammer shock testing in section 3.1.1. The evolutions of the 1 st damping ratios of each beam configuration and each corrosion rate, as a function of the ductility level are presented in Figure 5. The ductility level is evaluated in the same way as for quasi-static testing according to Equation 3. The maximum displacement involved in Equation 3is measured during the dynamic test performed before the concerned WN test. The computation of the yielding displacement was inspired from the Park's method [START_REF] Park | Evaluation of ductility of structures and structural assemblages from laboratory testing[END_REF] applied on the capacity curves issued from dynamic testing.

After the beams dynamic loading, an increase of the damping ratio with an increasing damage or ductility level is noticed [START_REF] Yuan | Shaking table tests of coastal bridge piers with different levels of corrosion damage caused by chloride penetration[END_REF]. Regarding the effect of the corrosion degree, the damping capability is reduced with an increasing corrosion degree especially for the 𝐶 2 beams. These observations are fully consistent with the quasi-static testing results.

However, the results of the current are not consistent with a study carried out on corroded RC bridge piers [START_REF] Yuan | Shaking table tests of coastal bridge piers with different levels of corrosion damage caused by chloride penetration[END_REF]. In the latter, an increase in the damping ratio at each ductility level with the increase of the corrosion degree is observed. The incoherence between the two studies may be due to the difference between the tested specimen (beam/column and corrosion along the beam/partial corrosion of the column) as well as the targeted corrosion degrees which are much higher (6 ~ 23% for longitudinal rebars and 10 ~ 46% for stirrups) [START_REF] Yuan | Shaking table tests of coastal bridge piers with different levels of corrosion damage caused by chloride penetration[END_REF].

Discussion

In this part, a comparison between the results exposed in sections 3.1.1 and 3.1.2 is carried out. This comparison is relevant since the same experimental setup was used for both quasi-static and dynamic testing and the modal damping ratio was computed based on the same identification method.

Prior to load testing, we can observe an increase in the modal damping ratio with an increasing corrosion degree. This conclusion is drawn by both quasi-static [START_REF] Zou | Experimental investigation on the dynamic properties of RC structures affected by the reinforcement corrosion[END_REF] and dynamic testing [START_REF] Yuan | Shaking table tests of coastal bridge piers with different levels of corrosion damage caused by chloride penetration[END_REF].

When the beams are loaded, both types of testing confirm the increase of the damping capability with an increasing damage level seen in [START_REF] Panteliou | Damping factor as an indicator of crack severity[END_REF].

Regarding the values of the 1 st modal damping ratios measured on the basis of quasi-static and dynamic tests, a slight difference is observed for a given beam at a given ductility level. This discrepancy may be explained by the fact that the injected energy during a hammer shock is not similar to the one injected during a WN signal.

The difference between the two levels of energy is both quantitative and qualitative. Not the same amount of energy is injected in the system during hammer shock or WN tests. In addition, the hammer shock testing imposes an impact force on the beam, whereas the WN signal leads to a distributed force that acts on the entire beam.

Jacobsen's areas method

3.2.1. Method Jacobsen [START_REF] Jacobsen | Damping in composite structures[END_REF] suggested a method to estimate the EVDR. This method is applicable for a single degree of freedom (SDOF) oscillator excited by mono-harmonic force. The EVDR is calculated according to Equation 4.

𝜉 𝑒𝑞 = 1 4.𝜋 . 𝜔 0 𝜔 . 𝐸 𝑑 𝐸 𝑠 (4)
where 𝐸 𝑑 is the dissipated energy during a cycle, 𝐸 𝑠 is the stored elastic potential energy, 𝜔 is the excitation frequency and 𝜔 0 is the eigenfrequency of the oscillator.

In a displacement-force diagram, the dissipated energy 𝐸 𝑑 is calculated as the area enclosed in the loop. Whereas, the stored energy 𝐸 𝑠 is computed by different methods based on the particularity of the tested elements [START_REF] Kumar | Cyclic response of sand using stress controlled cyclic triaxial tests[END_REF] [48] [10] [START_REF] Ma | Behavior of corrosion damaged circular reinforced concrete columns under cyclic loading[END_REF].

For both quasi-static and dynamic loading, asymmetries can be observed. For this reason, specific methods [START_REF] Kumar | Cyclic response of sand using stress controlled cyclic triaxial tests[END_REF] [48] derived from the former Jacobsen's area method are used to take into account these asymmetries in the computation of dissipated and stored energy. Since the EVDR depends on the excitation frequency 𝜔 and the eigenfrequency 𝜔 0 , the use of Jacobsen's area method in quasi-static testing needs to justified. For this reason, the deformed shape of DYSBAC specimen subjected to 4 point-bending test as well as the 1 st mode shape are calculated relying on the models shown in Figure 7. As depicted in Figure 8, the beam deformed shape under four point bending test and the beam 1 st mode shape are similar. The local error between the two shapes (see Equation 5) does not exceed 3%.

𝜀 𝑒𝑟𝑟 (𝑥) = 𝑢(𝑥)-𝜙(𝑥) 𝜙(𝑥) (5)
where 𝑥 is the position along the beam, 𝑢 is deformed shape and 𝜙 is the 1 st modeshape.

The normalized experimentally measured displacement by the DIC technique is also shown in Figure 8. a, for the reference beam at 0.8 mm imposed displacement. The differences noted between the theoretical and the experimental deformed shapes is due to the quality of the measurement itself (see section 2.4.3), the concrete cracking even in early stages of loading as well as the backlashes presence in the experimental assembly. Since the calculated local error between the two shapes is less than 3%, we can assume that the beam subjected to the four point bending test deforms as a SDOF oscillator excited by a mono-harmonic force with 𝜔 = 𝜔 0 .

Dissipated energy

Based on quasi-static tests, the dissipated energy is evaluated for each beam and for each loading block.

Based on the beam response to a given loading block, the area of the third cycle, where the damage is assumed to be stabilized, is computed. It should be noted that only loading blocks with an amplitude higher or equal to 4 mm have been considered in this part. This is because the measurement uncertainty is high for low imposed displacements. This may be explained by the precision of the measurement system itself as well as the unavoidable presence of backlashes in the experimental setup.

Figure 9 illustrates the evolution of the dissipated energy as a function of the ductility level for each tested beam. For all the tested beams, the dissipated energy increases with an increasing ductility level. This is due to the increase of damage level due to the concrete cracking and thus the dissipation capability. This observation is in a good agreement with the results reported in previous studies [START_REF] Panteliou | Damping factor as an indicator of crack severity[END_REF].

It is also observed that an increasing corrosion degree leads to a decrease in the dissipation capability, as shown in Figure 9. These results are consistent with the previous studies [10] [11] that showed that the corroded RC elements are characterized by a low hysteretic response and a reduced dissipation capability.

Regarding the difference between each beam configuration dissipation capability, no major differences could be observed between two different beam configurations corroded at the same degree.

Stored energy

The elastic stored energy is evaluated according to the method exposed in section 3.2.1 for each tested beam at each loading level. Figure 10 shows the stored energy evolution as a function of ductility level for each tested beam.

Based on the curves obtained, no major differences could be observed between the different tested beams as far as the elastic stored energy is concerned. These results are not consistent with those observed in [START_REF] Guo | Experimental investigation on the cyclic performance of reinforced concrete piers with chloride-induced corrosion in marine environment[END_REF]. The experimental study carried out on corroded bridge piers showed a decrease in the stored elastic energy as the severity of corrosion gets worse.

The discrepancy between the results of the current study and the conclusions drawn by [START_REF] Guo | Experimental investigation on the cyclic performance of reinforced concrete piers with chloride-induced corrosion in marine environment[END_REF] may be explained by the difference of the considered specimens in addition to the boundary and loading conditions. In [START_REF] Guo | Experimental investigation on the cyclic performance of reinforced concrete piers with chloride-induced corrosion in marine environment[END_REF], the specimens are cantilever columns fixed to the building rigid foundation, loaded by a constant vertical force and subjected to a cyclic horizontal loads. These settings affect the storing energy capability of the specimen as well as the appearance of plastic hinges [START_REF] Crambuer | Contribution à l'identification de l'amortissement: approches expérimentales et numériques[END_REF].

EVDR

The EVDR has been evaluated according to the Jacobsen's area method for each tested beam at each loading level. The curves describing the evolution the EVDR as a function of the ductility level are shown in Figure 11.

A decrease in the EVDR for the corroded specimens in comparison with the reference beam is observed. As exposed in section 3.1.1, the steel/concrete bonding enhancement in the case of slightly corroded structures can explain this conclusion. The obtained results are also in a good agreement with the conclusions drawn by [START_REF] Ma | Behavior of corrosion damaged circular reinforced concrete columns under cyclic loading[END_REF].

Non-parametric identification

In this part, the beam responses to high level dynamic loads are investigated. The non-parametric identification method is used to determine the evolution of natural frequency and the damping ratio over the testing time. This identification strategy is suggested in [START_REF] Heitz | Dissipations in reinforced concrete components: Static and dynamic experimental identification strategy[END_REF] and was applied on the response of RC beams subjected to an intermediate and high level WN signal on the shaking table .   The main advantage of this method is that no surrogate model has to be assumed in order to identify the modal properties evolution. The displacement field coming from DIC technique (see section 2.4.3) is projected on the analytical modal basis computed by the beam theory. The modal basis can be truncated or enriched to represent the nonlinearities occurring during the dynamic tests.

Time is divided into time windows. For each time window and each eigenmode i, the i th SDOF oscillator response, subjected to the i th modal force, is computed for a set of parameters (i th eigenfrequency, i th modal damping ratio), according to the Newmark method [START_REF] Newmark | A method of computation for structural dynamics[END_REF]. Subsequently, the resulting displacement is compared to the modal experimental displacement. An optimization algorithm is used to determine the modal properties (i th eigenfrequency, i th modal damping ratio) for each time-window and the modal displacement for each considered eigenmode. This algorithm is based on the minimization of an error function between the projected experimental displacement and the i th computed displacement.

The modal displacements are combined to reconstruct an approximate displacement field afterwards. An illustration of the algorithm output results is presented in Figure 12. This example corresponds to the reference beam subjected to an acceleration level of 0.5𝑔.

The choice of the time window size has a considerable effect on the results quality [START_REF] Heitz | Dissipations in reinforced concrete components: Static and dynamic experimental identification strategy[END_REF]. Therefore, rather than choosing a constant time window size [START_REF] Heitz | Dissipations in reinforced concrete components: Static and dynamic experimental identification strategy[END_REF], in the current study a variable window time size is considered. The size of each time window depends on the previous time window resulting eigenfrequency. It should include two periods of the previously computed SDOF oscillator (see Equation 7).

{ 𝐿 0 = 0.3 𝑠 𝐿 𝑖 = 2 𝑓 𝑖-1 (7)
where 𝐿 0 is the initial time window size, 𝐿 𝑖 is the size of the current time window and 𝑓 𝑖-1 the 1 st eigenfrequency of the computed SDOF oscillator in the previous time window.

Since the final result, which is the approximate displacement field, is the combination of the different approximate modal displacements, the size of each time window is determined using the 1 st natural frequency. The so-computed size of the successive time windows is then kept for the other eigenmodes.

In the present study, a good match between the numerical and experimental displacement field could be found by considering only the first two eigenmodes.

From Figure 12.a, in a first phase a decrease in the 1 st natural frequency along with time is observed. This result can be explained by the cracks development and the damage growth due to the seismic loading. A second phase where the natural frequency increases over time follows. This is indicative of the crack closure and the gain in stiffness due to the decrease of the excitation level at the end of each testing sequence. The shape of the curves is found for all the tested beams at each acceleration level. The evolution of the 1 st modal damping ratio over time has been determined using the non-parametric identification method for each tested beam. Four accelerations levels have been considered for all the beams: 0.125𝑔 , 0.5𝑔 , 0.8𝑔 and 1.25𝑔.

Figure 13 shows the evolution of the mean 1 st modal damping ratio at each considered acceleration level as a function of the level. From the second acceleration level (2 nd point in each curve), it can be observed that, the general trend is an increase of the modal damping ratio as the damage level increases for all the tested specimens. This finding is similar to the results mentioned in part 3.1. The high mean damping ratio for 0.125𝑔 acceleration level is explained by the limited quality of the displacement field measurement. At this acceleration level, the measured displacement is in the range of 0 -1 𝑚𝑚, which made difficult the displacement field acquisition with a good precision.

Regarding the effect of the reinforcement corrosion on the modal damping ratio, no particular trend could be discerned using this specific identification method. It is probably due to the choice made of the average value as a representative parameter of the damping ratio time evolution. Given the fact that the mean values smooth out the local variability, which is predominant in our case of study. However, it should be noticed that the advantage of this method is the evaluation of the modal damping during a seismic loading. These identified values are different from those identified from a hammer shock or a WN This is mainly because the injected energy in the system during a seismic sequence is higher than the injected energy during the modal analysis testing.

Conclusions

In this paper, an extensive experimental campaign aiming at assessing the influence of corrosion on the dissipation capability of RC components has been detailed. 18 corroded beams and 2 additional reference beams have been subjected to quasi-static and dynamic testing using the same experimental setup.

The main results can be summed up as follows:

• when the specimens are unloaded, an increase in the modal damping ratio with an increasing corrosion degree is observed. The same conclusion can be drawn based on both quasi-static and dynamic testing; • when the specimens are loaded, either by quasi-static or dynamic loading, the modal damping ratio measured on corroded specimens decreases as the corrosion severity becomes worse; • the modal damping ratio increases with an increasing damage level for all the tested beams.

This conclusion is confirmed by quasi-static and dynamic and quasi-static tests; • the modal damping ratios have been computed using the same identification method on the results issued from both hammer shock and WN tests. A slight difference in the damping ratio values is noticed. This is mainly due to the difference between the two excitations (local/distributed); • for the quasi-static tests, the dissipated energy and the equivalent viscous damping ratio decrease with an increasing corrosion degree; • for the dynamic tests, the non-parametric identification method was applied. The effect of the corrosion on the average value of the modal damping ratio evolution during a test sequence could not be discerned. However, it should be noticed that this method made possible the identification of the modal damping ratio evolution during the seismic sequence. The obtained values are different from the modal damping ratios coming from the modal analysis tests; • the covariance-driven SSI method gives consistent results for hammer shock and white noise tests. Jacobsen's area method applied on quasi-static cyclic responses, gives the same trends as revealed by the hammer shock and white noise characterization. The nonparametric identification method reveals the high variability of the damping ratio of the tested specimens during a given seismic sequence. Hence, this identification method seems to be the most suitable to identify a non-linear model.

As a major outlook of the present study, a simplified numerical model being able to take into account the different effects of corrosion at the material scale is under development. The main objective of this simplified model is to reproduce the different structural consequences observed by this experimental campaign [START_REF] Lejouad | A numerical model to predict the effect of corrosion on the dynamic behavior of reinforced concrete beams[END_REF] .

Figure 1 :

 1 Figure 1: General view of the experimental setup for the quasi-static tests.

Figure 2 :

 2 Figure 2: Applied loading in quasi-static testing.

Figure 3 :

 3 Figure 3: General view of the experimental setup for dynamic tests.

Figure 4 :

 4 Evolution of first damping ratio obtained from hammer shock tests: (a) 𝐶 1 configuration, (b) 𝐶 2 configuration, (c) 𝐶 3 configuration.

Figure 5 :

 5 Evolution of first damping ratio obtained from dynamic tests: (a) 𝐶 1 configuration, (b) 𝐶 2 configuration, (c) 𝐶 3 configuration.

Figure 6

 6 Figure 6 is a representation of the dissipated and stored energy involved in the calculation of EVDR in this paper. The plotted displacement-force diagram corresponds to the response of the reference beam subjected to a cycle of 5 𝑐𝑚 amplitude.

Figure 6 :

 6 Figure 6: Jacobsen's areas method applied to the virgin beam response subjected to a 5 cm amplitude cycle.

Figure 7 :

 7 DYSBAC beam modeling: (a) Quasi-static, (b) Dynamic.

Figure 8 :

 8 Comparison between the computed deformed shape under four point bending test and 1 st modeshape : (a) Normalized displacement, (b) Local error

Figure 9 :

 9 Evolution of the dissipated energy: (a) 𝐶 1 configuration, (b) 𝐶 2 configuration, (c) 𝐶 3 configuration.

Figure 11 :

 11 Evolution of the EVDR: (a) 𝐶 1 configuration, (b) 𝐶 2 configuration, (c) 𝐶 3 configuration. (a) (b) (c) Figure 10: Evolution of the stored energy: (a) 𝐶 1 configuration, (b) 𝐶 2 configuration, (c) 𝐶 3 configuration.

Figure 12 :

 12 Example of the algorithm output results for the reference beam at 0.5g acceleration level: (a) the 1 st eigenfrequency evolution, (b) the 1 st modal damping ratio evolution, (c) the 1 st modal displacement evolution, (d) the midspan displacement evolution.

Figure 13 :

 13 Evolution of the 1 st modal damping ratio: (a) 𝐶 1 configuration, (b) 𝐶 2 configuration, (c) 𝐶 3 configuration.

  [START_REF] Meda | Experimental evaluation of the corrosion influence on the cyclic behaviour of RC columns[END_REF] 𝑚 long with a cross-section of 20 × 40 𝑐𝑚 2 .

The reinforcement pattern is designed according to the European standards Eurocodes 2 and 8

[START_REF] Mosley | Reinforced concrete design: to Eurocode 2[END_REF]

[START_REF] Fardis | Seismic design, assessment and retrofitting of concrete buildings: based on EN-Eurocode[END_REF]

. The steel reinforcement consists of ribbed B500A steel bars according to the French steel classification NF A 35-080-1

[START_REF] Afnor | Steel for the reinforcement of concreteweldable reinforcing steel -General[END_REF]

. 4 bars of 12 𝑚𝑚 diameter constitute the longitudinal reinforcements and the transversal reinforcement is formed of 8 𝑚𝑚 diameter stirrups with a spacing of 10 𝑐𝑚. The concrete cover thickness is 3 𝑐𝑚.

Table 1 :

 1 Mean weighting results.

	Target corrosion rate (%) Measured corrosion rate (%)
	5	3,1
	10	6,3
	15	9,4

Table 2 :

 2 Statistical indicators of the measured diameter distributions.

	Targeted mass loss	Minimum	Maximum	The average	Standard
	(%)	diameter (mm)	diameter (mm)	diameter (mm)	deviation (mm)
	0	11.693	13.555	12.096	0.144
	5	10.901	12.598	12.083	0.167
	10	9.820	12.983	11.997	0.300
	15	9.656	12.766	11.797	0.295

Table 3 :

 3 Natural damping ratio of the specimens before loading based on hammer shock modal testing.

	Configuration	Reference		𝐶 1			𝐶 2		𝐶 3
	Corrosion ratio (%)	0	5	10	15	5	10	15	5	10
	1 st natural frequency (Hz)	12.31	12.33 12.72 12.45 12.35 12.40 12.38 12.12 12.42
	1 st modal damping	2.5	2.7	3.4	3.4	2.7	2.6	2.7	2.5	4.2
	ratio (%)

Table 4 :

 4 Natural damping ratio of the specimens before loading based on white noise modal testing.

	Configuration	Reference		𝐶 1			𝐶 2		𝐶 3
	Corrosion ratio (%)	0	5	10	15	5	10	15	5	10
	1 st natural frequency (Hz)	12.21	12.34 12.69 12.01 12.30 11.92 12.03 12.32 12.151
	1 st modal damping	1.4	1.8	2.1	2.1	1.4	2	3.1	2.8	2.9
	ratio (%)

Acknowledgements

The French Institute for Radioprotection and Nuclear Safety (IRSN) and the Nuclear Energy Division of the French Sustainable Energies and Atomic Energy Commission (CEA/DEN) are kindly thanked for financial supports.