1. Introduction.
Worldwide, the design and construction of bridges is governed by regional standards and involves a series of studies prior to their design and construction that are supported by permanent research aimed at reducing the risk of collapse or premature failure; These standards include operational maintenance protocols, however, failures and even collapses frequently occur. A sample of 384 collapsed bridges worldwide in the period from 1944 to 2004 revealed that 28% collapsed due to natural causes, 21% due to design errors, 18% due to impact, 10% due to overloading, 10% due to human error, 10% from unknown causes, 2% from corrosion and 1% from vandalism (Imhof, 2004). The greatest difficulty in the inspection of bridge structures is obtaining a damage index for making maintenance, repair, or replacement decisions.
A summary of some methods to obtain mechanical response in bridges is presented below: For non-linear ranges, the dynamic plastic hinge method was proposed; It is used to obtain the modal load-deformation curve from a one-degree-of-freedom model. The structure is subjected to the acceleration of a given earthquake until the maximum response is obtained in which the plasticization of a previously proposed ball joint is reached, the system becoming a mechanism; This simplified method allows visualizing the ductility demand of the structure (E. Maldonado, et al. 1998). By means of numerical methods, an amplification factor can be obtained that considers the dynamic effects applicable to the mechanical design elements of continuous bridges, considering flexible beams and mobile loads damped by the suspension of the analyzed type trucks. motion equation of the bridge is assumed in forced vibration by the equation of motion of the test truck that has a point of application of load per axle, the beam is idealized as a series of masses concentrated in discrete intervals uniformly distributed in its total length and the effects of truck loading are computed on the front axle; This procedure makes it possible to obtain an impact factor for the mobile load, also known as the dynamic amplification factor, which results from the displacement curves generated for both static and dynamic loads.
The procedure allows to include a friction factor due to the participation of the suspension of the analyzed truck (N. Munirudrappa et al. 1999). The method frequency response functions consist of obtaining a structural characteristic matrix, which is a function of mass, damping and its stiffness. It is obtained by inverting the frequency vector functions cofactors matrix of the motion equation´s Fourier transforms and is known as frequencies response functions matrix. This matrix has only the structure properties and is independent of the excitation, so any change generated in it is due to changes in the constitution of the structure. To obtain the frequency response functions matrix of the displacement vectors and their excitation forces must be known in terms of their Fourier transforms. This is a complex method compared to the modal analysis method, since it requires a much greater number of sensors and a large numerical processing that couples the local and global effects.
The sensitivity analysis method consists of obtaining a sensitivity factor that is equivalent to the quotient between the rate of change in frequency and the rate of change in stiffness. This method helps to select the parameters with the greatest impact on the response and to rule out possible errors in the experimental instrumentation, in such a way that it allows to choose the most significant or highest impact values and to visualize the structural elements with the greatest contribution to dynamic effects. Based on this information, it is possible to determine with good precision the points where the sensors should be placed for experimental studies. A group of experimental methods that make use of numerical simulation methods, where numerical methods such as Runge-Kutta, Euler, etc. can be used to solve the differential equations of motion from the Jacobian of the quotient between the rate of change of the forces and the rate of change of the displacements whose real solutions must be negative to guarantee the stability of the method and its integration in time. In short, the finite element method is the most used to solve the motion equations, as it uses discretization and interpolation using functions that are compatible with stresses and strains, considering internal equilibrium, equilibrium at the boundary and the compatibility of deformations of the discrete elements that are solved by direct integration and assembled in linear systems of equations provided that the differential equations are ordinary (Carrión FJ et al. 1999). A model to study the effects of temperature compared to the effects of damage in the change of the natural frequency of a span bridges, using finite element models for box girders and prestressed ASHHTO type girders; The article indicates that the changes in the frequency of the superstructure of a bridge of a clearing due to thermal effects (temporary softening), are similar to the effects of damage accumulated over a long period of time (Balmes Etienne et al. , 2005).
A study was carried out by instrumenting two slab bridges and prestressed girders, both of two independent spans; data was collected by means of acceleration sensors; numerical methods were used to find the spectra of speeds and displacements. Likewise, with an adequate frequency filtering, it was possible to separate the static and dynamic effects of the obtained spectrum; The results evaluation indicates that the impact factors obtained by means of the quotient between the load dynamic effect and the load static effect, under certain mass and speed conditions of the test vehicle, exceed those established in the design codes. (Valdez J. et al. 2008). A comparative analysis to obtain the dynamic response of a reinforced concrete bridge located in Italy using a finite element model with Shell-type elements and the excitation of a typical truck at different speeds.
The theoretical analysis was carried out assuming the non-existence of cracking to ignore the contribution of reinforcing steel in estimating stiffness. The results obtained in the analysis indicated that with a single truck it was not possible to obtain the basic dynamic parameters, so a continuous convoy of trucks was used, offering results more in line with reality; the contrast was performed with the data obtained from the application in the excitation site using an electric vibrator placed in the quarter of the bridge span and placing a series of 17 acceleration sensors to obtain the response in real time. The results obtained in the analysis presented errors ranging from 1% to 46% compared with the frequencies obtained on site, the torsion frequency being the one with the lowest error and that of the fourth bending mode the one with the highest error (Veles H. et al. 2011). A methodology to determine the concrete bridges superstructure fatigue deterioration of the reinforced, using Monte Carlo techniques; The method estimates probabilistically according to the statistical data the number of load cycles and their impact on the stress level, using the most popular damage models that allow an estimate of the crack size to be given. This method is speculative and can be applied with relative ease if the operating statistics of the bridge under analysis are known, however, as it is presented, it does not include factors such as extraordinary loads, or earthquakes and corrosion effects, among others, (Crespo E. et al. 2013).
The results of the application of the environmental vibrations method (spectral ratio) applied to the instrumentation of a post-tensioned bridge that allowed them to evaluate the superstructure’s present state, the instrumentation results were compared with the results of the structural model made using commercial software. (Viviescas Al. et al. 2017).
2. Methodology.
Of the bridge infrastructure on federal highways in Mexico, 57% are built of reinforced concrete, followed by 28% built of prestressed concrete, 7% built in steel and concrete, 6% in structural steel, 1% in masonry and 1% in 3dsalb system (IMT-2014). According to the statistics of the Mexican Institute of Transportation, the largest number of bridges is built using reinforced concrete, so the present work studies two reinforced concrete bridges; one based on solid slabs and the other built on beams and slabs with diaphragms of reinforced concrete.
The present work consists of the application of a methodology that uses global parameters; includes the real present point stiffness of the structure, which is obtained from field measurements on the real structure; the value obtained is compared with the stiffness value obtained from the design parameters of the bridge that result from the stiffness invariant function; the quotient between these values corresponds to the residual capacity of the structure, allowing a damage index to be obtained.
2.1. General concepts.
The structural mechanical behavior under service loads depends on the level of efforts achieved and the number of repetitions or load cycles; Factors such as deterioration, corrosion, fatigue and / or an increase in the level of service loads, can generate permanent damage that modify the conditions of its mechanical response. Design considerations for the structural performance is based on elastic linear mechanical behavior and the structure is considered healthy if this behavior is preserved. When the accumulated damage modifies this behavior, non-linear behavior begins where the proportion between the displacements and the applied forces is no longer constant; It is under this basic principle that the proposed method offers information on the present structure state.
Figure 1 presents the structural system’s mechanical behavior under monotonic load with gradual load increase; First, a straight line can be seen beginning at the origin and reaching the coordinates
The proposed method considers that if three increasing point loads are applied, a load-displacement graph is obtained that describes a mechanical behavior like the real one; From the applied load increments, joining the coordinates of the three points, two lines are obtained that have the same slope if working in a linear range or secant to the real stress-strain curve of the structure if the work is in a non-linear range. If the load increments are small, the upper and lower areas of the real curve are remarkably like the upper and lower areas of the graphs of the lines obtained, reducing the error of the method. The areas described provide a way to measure the damage.
2.2. Damage assessment.
Figure 1 is divided into two regions, a lower one whose area corresponds to the work W that was produced during the loading process and an upper one whose area corresponds to the deformation energy U stored by the structure, which allows it to recover its original shape if the load is removed fully or partially. When the structural mechanical behavior is linear elastic, the work and the deformation energy have the same magnitude, which implies that the structure can recover its original configuration if the loads are removed; if the behavior is non-linear, the values of the work done and the stored deformation energy are different, which implies that, when the load is removed, the structure will partially recover its configuration. Disregarding the energy generated in the form of heat, if the structure has non-linear behavior, the stored deformation energy is less than the work done by the system; This consideration is applicable in reinforced concrete structures.
There are some semi-empirical expressions to evaluate the damage in reinforced concrete structures obtained from experimental tests; One of the expressions that allows obtaining a damage index considering the slope of the elastic part of the discharge curve, is defined by the following expression:
Where d is a damage index, Z is the slope of the discharge branch elastic part and
According to the proposal of the present work, the real deformation is obtained from the acceleration-time spectrum double integration that results from the impact applied to the real structure with known masses and their free fall heights, only considering the maximum amplitude of the displacement; from these values a kinetic energy vs. displacement graph can be constructed; the upper area of the graph is identified with area
2.3. Invariant elastic stiffness.
It is formed with the stiffness values below the limit of proportionality along the length of a structural element. In (3) K, corresponds to the theoretical design stiffness value, F corresponds to the applied force and δ is the resulting displacement.
If the virtual work principle is applied and considering only the contribution of the bending moment, the expression for the stiffness at a point corresponds to (4).
If the value of the force F is kept constant throughout each point of the structure, under the limit of proportionality the stiffness curve is obtained and has the character of invariant at each point.
According to figure 2, if F is applied in the coordinate X = a, it is necessary to determine the displacement δ and the stiffness K (a).
Once the displacement is known, the stiffness is obtained in the coordinate x = a; considering b = L-a, where (6) is the stiffness invariant applied for beams and (7) the invariant for solid slabs.
2.4. Measurement of actual point stiffness.
The actual point stiffness is obtained with the maximum displacement that occurs when subjecting the superstructure to impact loads applied to the center of the span; Small masses are used for excitation, which minimally modify the dynamic parameters of the structure. To calculate the impact force
Since the impact is made in a deformable medium, the magnitude of the force depends on the reaction stiffness; (9) corresponds to the kinetic energy now of impact.
From the kinetic energy-displacement graph, as previously mentioned, the upper area of the curve
The residual stiffness factor ∆k (2) is applicable to reinforced concrete structures. It corresponds to the increase in the cracking of the cross section, due to the creep of the reinforcing steel as damage accumulates which is reflected in the decrease of the compression area since the cracking as the steel receives greater deformation, grows in the tension zone. Therefore, the moment of inertia of the cross section is reduced, expressed by
With (12), the actual present stiffness
In(13),
3. Study cases.
The instrumented bridges are located on federal highway 14, Hermosillo-Moctezuma, in the state of Sonora, México; one of the bridges bears the name of "El Testarazo" figure 4, which is located at km 23 + 900. The bridge has a superstructure made of three reinforced concrete slabs, simply supported on abutments made up of reinforced concrete walls. The other instrumented bridge is called "El Gavilán" figure 5; The super structure is formed by a system of reinforced concrete beams, slabs, and diaphragms, with a section skewed at 48 ° from its transverse direction, located at Km 60 + 100 of the same roads.
3.1. Information on the case studies.
The cross sections geometric properties were obtained under the transformed section criterion, and the project data were obtained from the Type Project of Reinforced Concrete Elements, Part I, published by the now extinct SAHOP.
Bridge | bw | hf | bf | n | As(m2) | d | kd | Icrt(m4) | A(m2) |
---|---|---|---|---|---|---|---|---|---|
Testarazo | 5.2 | 0.30 | 10.0 | 8.796 | 0.035 | 0.645 | 0.1858 | 0.1108 | 2.95 |
El Gavilán | 1.08 | 0.18 | 7.5 | 8.796 | 0.0386 | 2.182 | 0.481 | 1.199 | 2.89 |
For the bridge “El Testarazo” EI = 2,595
3.2. Stiffness invariant curve.
Figures 6 and 7 show the cofactors graphs corresponding to the stiffness invariants, obtained from (4) for the “El Gavilán” bridge and from (5) for El Testarazo bridge. The
3.3. Procedure to obtain the structure real stiffness.
The field measurements were made using a structure for lifting and unloading the excitation masses. The equipment consists of a quadruped skeleton that allows the masses to be raised by means of a 2250 N capacity winch; the container holding mechanism is made up of an electromagnet with a capacity of 6 KN that enables the excitation mass to be clamped by means of a safety hook with a torsional freedom degree (fig. 8).
Sand containers were used, with plastic behavior during the impact, to avoid rebound; Three sacks of 25, 50 and 75kg capacity were manufactured, which were filled up to the test mass; the impact was achieved by raising the sacks to an average height of 1.50m, subsequently, the electromagnet electrical flow that held them was interrupted to release them and produce the impact on the bearing bridge surface; The remaining equipment consists of a low frequency acceleration sensor (0.2 Hz), with a sensitivity of 500 mV / g, placed in bridge´s mid span that allowed obtaining the response in real time. In addition, a 4-channel capture card for data reception from 0 to 25 khz was used. The card allowed the capture of analog signals produced by the sensor during measurements. The capture card was placed on a chassis with capacity for eight 11-30 V cards of 15 W, to operate from -40º to 70º C communicated to the USB port of the software carrier computer to process the analog signal where the spectra were obtained from the acceleration-time response of the structure. The analog signals captured were processed using the Labview Signal Express version 3.0 Software, license 501701A-00, which allowed the capture of the time-acceleration spectrum in a numerical matrix in ASCII code TXT format, in raw state for numerical processing. The acceleration-time spectra obtained are presented in Figures 9, 10 and the displacements were obtained by twice integrating the spectra with the Constant Average Acceleration and Linear Acceleration Methods.
3.4. Field measurements results.
Figures 9 and 10 show the acceleration-time spectra captured during field measurements.
The acceleration-time spectra of Figures 9 and 10 were subjected to a numerical treatment, which consisted, first, in the correction of the spectral baseline, and later, in the double integration with the methods already described to obtain the maximum displacements produced by the impact of the masses. Tables 3 and 4 show the results that include the test mass, the height of free fall and the displacements obtained from the numerical processing.
MASS (Kg) | HEIGHT (m) | DISPLACEMENT (m) |
---|---|---|
25 | 1.525 | -2.92927E-05 |
25 | 1.5 | -2.89101E-05 |
25 | 1.484 | -2.84981E-05 |
25 | 1.46 | -2.8037E-05 |
25 | 1.48 | -2.83313E-05 |
25 | 1.55 | -2.98126E-05 |
50 | 1.38 | -5.11297E-05 |
50 | 1.53 | -5.69176E-05 |
50 | 1.55 | -5.76141E-05 |
50 | 1.47 | -5.473E-05 |
50 | 1.548 | -5.7467E-05 |
50 | 1.56 | -5.77907E-05 |
75 | 1.51 | -8.83783E-05 |
75 | 1.49 | -8.69951E-05 |
75 | 1.5 | -8.74758E-05 |
75 | 1.51 | -8.85254E-05 |
75 | 1.49 | -8.7103E-05 |
MASS (Kg) | HEIGHT (m) | DISPLACEMENT (m) |
---|---|---|
25 | 1.54 | -0.000167776 |
25 | 1.52 | -0.000141068 |
25 | 1.55 | -0.00017 |
25 | 1.5 | -0.00013734 |
25 | 1.53 | -0.000163631 |
50 | 1.5 | -0.000260161 |
50 | 1.53 | -0.000268009 |
50 | 1.52 | -0.000266636 |
50 | 1.5 | -0.00026644 |
50 | 1.54 | -0.000273797 |
50 | 1.5 | -0.000258101 |
75 | 1.67 | -0.000491775 |
75 | 1.5 | -0.000440763 |
75 | 1.52 | -0.000448023 |
75 | 1.4 | -0.000412609 |
75 | 1.39 | -0.000408783 |
75 | 1.41 | -0.000414767 |
3.5. Statistical analysis of field data.
To know the validity of the data obtained in the measurements made on a natural scale of the study cases, we proceeded to the analysis of variance or ANOVA; The main purpose is to know with the bifactorial analysis if the displacements obtained are dependent on the mass and height in free fall used in the tests and to rule out the possibility that other factors have influenced the results. The null hypothesis
3.5.1. Analysis of variance for the data in Table 3.
Degrees of freedom | Sum of squares | Average of squares | F | Critical value of F | f | |
---|---|---|---|---|---|---|
Regression | 2 | 1.04213E-08 | 5.2107E-09 | 2866.0866 | 4.2148E-20 | 0.05146909 |
Waste | 15 | 2.72706E-11 | 1.818E-12 | |||
Total | 17 | 1.04486E-08 |
Since the distribution function f <F, the null hypothesis
Origin of variations | Sum of squares | Degrees of freedom | Average of squares | F | Probability | Critical value for F | f |
---|---|---|---|---|---|---|---|
Between groups | 20018.428 | 1 | 20018.4283 | 93.677634 | 5.0095E-11 | 4.1490974 | 0.0039 |
Within groups | 6838.2352 | 32 | 213.694853 | ||||
Total | 26856.663 | 33 |
As in the bifactorial case, from the values in Table 6, the null hypothesis is discarded, and the alternative hypothesis is accepted since f <F; it can be stated that the displacements obtained are dependent on the applied mass with a 95% confidence value.
Origin of variations | Sum of squares | Degrees of freedom | Average of squares | F | Probability | Critical value for F | f |
---|---|---|---|---|---|---|---|
Between groups | 19.17911772 | 1 | 19.1791 | 20252.1 | 2.061E-46 | 4.149097 | 0.00399 |
Within groups | 0.03030448 | 32 | 0.00095 | ||||
Total | 19.2094222 | 33 |
As can be seen in table 7, f <F which allows ruling out the null hypothesis, so the displacements obtained are dependent on the fall heights of the masses with a 95% confidence value.
3.5.2. Analysis of variances for the data in table 4.
Degrees of freedom | Sum of squares | Average of squares | F | Critical value of F | f | |
---|---|---|---|---|---|---|
Regression | 2 | 2.22636E-07 | 1.1E-07 | 322.93 | 1.93519E-12 | 0.051481683 |
Waste | 14 | 4.82596E-09 | 3.4E-10 | |||
Total | 16 | 2.27462E-07 |
Table 8 presents the results of the analysis of variance; since f <F, the null hypothesis
Origin of variations | Sum of squares | Degrees of freedom | Average of squares | F | Probability | Critical value for F | f |
---|---|---|---|---|---|---|---|
Between groups | 22518.64 | 1 | 22518.6391 | 105.37 | 1E-11 | 4.14909744 | 0.00399430 |
Within groups | 6838.235 | 32 | 213.694852 | ||||
Total | 29356.87 | 33 |
From the analysis in Table 9, the null hypothesis is discarded, and the alternative hypothesis is accepted since f <F; it can be stated that the displacements obtained are dependent on the applied mass with a 95% confidence value.
Origin of variations | Sum of squares | Degrees of freedom | Average of squares | F | Probability | Critical value for F | f |
---|---|---|---|---|---|---|---|
Between groups | 19.31294 | 1 | 19.312944 | 9230.5 | 6E-41 | 4.1490974 | 0.0039943 |
Within groups | 0.066953 | 32 | 0.0020922 | ||||
Total | 19.3799 | 33 |
As illustrated in Figure 10, f <F; the null hypothesis is discarded, so the displacements obtained are dependent on the fall heights of the masses with a confidence value of 95%.
Regarding the analysis of variance of the data obtained in the field measurements of both bridges; It can be affirmed that the results of the measured displacements depend on the masses used in the impact and their free fall heights with a confidence level of 95%, likewise, it is also affirmed that they are the factors of greatest influence with a probability greater than 99%.
3.6. Damage factor estimation.
From the the signals’ numerical processing show in Figures 8 and 9, the velocity and displacement spectra were obtained; Knowing the masses´s free fall heights, we proceeded to the construction of the kinetic energy vs. displacement graphs, using (9) to estimate the kinetic energy. Figures 11 and 12 show kinetic energy vs average displacements graphs for the study cases.
The values obtained for the reduction inertia moment by damage according to (10) and (11), are presented in table11; the areas correspond to the upper and lower surfaces of the 3 points’ enveloping.
4. Results and discussion.
To find the effective forces that defined the values of the real present stiffness obtained because of the masses’ impact, (11) was applied; Table 12 shows the results of the forces and the average displacements obtained that were used to obtain the average values of the actual superstructures point stiffness of the study cases.
BRIDGE | L (m) | E (Mpa) |
|
|
|
|
---|---|---|---|---|---|---|
El Testarazo | 12.50 | 23414 | 0.1108 | 0.1015 | 0.0000876954 | 5121.88 |
El Gavilán | 32.00 | 23414 | 1.199 | 0.9352 | 0.00043612 | 13988.70 |
Once the average values of the displacements and the real forces are obtained, (12) can be applied to obtain the real stiffness value for the elements under study. The contrast was performed with the values obtained from the EI product by the stiffness invariant cofactors illustrated in Figures 6 and 7, respectively, for each bridge; the results are presented in table 13.
PUENTE | Masa(kg) |
|
|
|
|
|
El Testarazo | 75 | 0.00876955 | 5121.88 | 58.405 | 66.432 | 12.08 |
El Gavilán | 75 | 0.043612 | 13988.70 | 32.075 | 41.136 | 22.03 |
When assessing the mechanical status of the studies cases, it is important to clarify that the damage index for healthy structures must have a null value; According to table 14, which summarizes the results obtained, the following is stated: For the bridge "El Testarazo" built with solid slabs, the damage index value is 12.08%, which indicates that the structure has 87.9% of its original capacity; It can be interpreted that, during the service life, the superstructure has lost rigidity by 12.08% compared to its original state, accumulating irreversible damage, as manifested in its pathology by excessive deflection and alternating transverse cracking at the center of the clearing. The evidence presented in fig. (10), where the three points obtained for the different masses form two lines with different slopes with decreasing behavior, showing that the state of the structure exceeded the limit of proportionality.
In the case of the “El Gavilán” bridge built with reinforced concrete beams, slabs, and diaphragms with a deviation of 48 °, the resulting damage index is 22.03%, which indicates that it has a capacity of 77.97% with respect to its design capacity, which matches cracking and excessive sag in the structure. It is important to clarify that the stiffness to vertical displacement depends on the contribution to bending and torsion (Deng Kai, 1998); Under this consideration, the calculated damage index expressed in table 13 includes the reduction in flexural and torsional stiffness, but it is not possible to distinguish which percentage corresponds to each degree of freedom.
5. Conclusions.
The use of the stiffness invariant allowed the obtention of damage index values for the study cases according to their physical conditions. The procedure is relatively simple, especially when you have the information used in its construction. The ease in the analysis to obtain the contrast data, likewise, the ease with which the real stiffness values are obtained in present time, are the main advantages of the method. The disadvantages correspond to the measurement stage in the field since it requires environmental conditions of constant temperature and the absence of wind. Another disadvantage is that for skewed bridges, a more refined analysis is required to obtain the stiffness invariant, since the reduction percentages corresponding to both flexural and torsional stiffness are not appreciated with the proposed procedure.
The method is governed by energetic principles; it can be used in reinforced and prestressed concrete structures with an acceptable approximation, since the amount of deformation energy includes the area of the region that forms the secants with the real force-displacement curve and corresponds to method error; This error is a small portion of the work done that is added to the deformation energy, so the test masses should be chosen with the smallest possible increments in order to minimize the error or, failing that, estimate the error logistically adjusting the three points and obtaining the area of the arcs to be subtracted from the deformation energy and added to the work done. The development of the adjustment by mistake is outside the scope of the present work.
Another condition in the use of the method is that the structure's own weight must be an important part of the service load, so it is recommended for use in slabs and beams bridges.