Abstract

Concrete simply supported beams are widely used in small- and medium-sized bridges. Hinged slabs and joints are damaged in service easily, as they are the main load-carrying members of bridge. In this paper, using the high sensitivity and strong antinoise characteristics of strain energy to damage, a strain energy information entropy index is proposed for identifying local damage in hinge joints and hinged slabs. Moreover, a new damage index of the local strain energy entropy function is constructed according to the strain energy of a hinge joint element and its surrounding elements. A simply supported beam model is established to verify the effectiveness of the proposed index. The numerical simulation results show that the proposed damage index of hinge joints can accurately locate and quantify the degree of damage of single-hinge joints and multihinge joints. The improved modal strain energy information entropy index is sensitive to single-hinge slab damage and can accurately identify damage with noise levels of 5% and 10%.

1. Introduction

The simply supported beam is widely used in small- and medium-sized span bridges due to its simple structure, convenient construction, and mature technology. It consists of a number of concrete slabs connected by tongue-and-groove concrete joints to bear vehicle loads [13]. With time, the bridge inevitably suffers structural damage under the passing of heavy loads. Most of the damage occurs in the essential parts of the bridge; specifically, the joints between hinged slabs are easily damaged [46]. In one area of China, 402 out of 429 simply supported beam bridges have cracks at the bottom of hinged slabs, with an incidence of about 93.5%, as reported in the survey data of Qin et al. [7]. Moreover, 412 bridges had hinge joints damaged by water seepage and concrete debris, with an incidence of about 96.0%. Therefore, the working state of simply supported beam bridges should be considered and evaluated for further operation and maintenance.

With the development of technology, a number of structural damage detection techniques have been developed in the past few decades. One of these, the vibration-based techniques, have drawn much attention by researchers, such as the methods based on modal frequency [8], modal shape or mode shape derivatives [911], modal flexibility [12], modal strain energy (MSE) [13], and so on. Modal strain energy performs particularly well, showing a good sensitivity to damage of structural. Shi et al. [13] were the first to use the MSE method to identify structural damage; only a few lower displacement modes were used to calculate the MSE before and after damage. Shi et al. [14] also proposed an improved MSE method that significantly reduces the modal truncation error and finite-element modeling errors from the analytical model in the computations. A two-stage structural damage detection method based on MSE and particle swarm optimization was proposed by Seyedpoor [15] to improve the accuracy of damage detection. Jahangiri [16] proposed a novel identification method for structural damage assessment using the first vibration mode data, which not only predicts the damage location properly but also estimates the damage intensity successfully. The MSE method usually requires complete information of structural mode shapes. However, complete modal shape information is hard to obtain due to the limited number of sensors and the difficulty of measuring rotational vibration. Wu [17] proposed a novel approach, where the displacement mode and rotational mode included in the definition of the elemental MSE are replaced by strain mode, effectively improving the accuracy of the MSE.

The randomness and uncertainty of actual structural damage generally bring great challenge to the damage condition assessment. Although the traditional damage identification method based on structural vibration characteristics performs well in damage location, there are still many defects in damage quantitative analysis. Therefore, the information entropy theory, which describes the uncertainty of the system, is combined with the characteristic parameters of structural vibration to measure the damage of the structure from the perspective of probability. Recently, the information entropy theory has been introduced into structural damage identification and has been applied in dam safety analysis [18], surrounding rock stability [19], mechanical bearing damage [20], and other fields. However, there are few published studies on the field of structural beam damage. Liu et al. [21, 22] introduced the information entropy theory into the dynamic damage identification of concrete structures, providing a new idea for solving problems such as insufficient sensitivity and difficulty in nonlinear analysis for damage identification. Moreover, some scholars have successfully applied the information entropy theory to the research on crack propagation mechanism of concrete structures and achieved good results [23, 24].

In this paper, by combining the strain energy theory with the theory of information entropy, an index of information entropy of strain energy is proposed and studied to identify damage at hinge joints and hinged slabs.

2. Theoretical Basis of Damage Identification Method

2.1. Information Entropy Theory

Entropy is a special parameter in thermodynamics and statistical physics; it was proposed by the famous German physicist Clausius in 1865. Entropy is a state variable, and the change in entropy relates to the dispersion of matter; furthermore, its macroscopic definition indicates the degree of uniformity of the energy distribution in the system. In 1948, C. E. Shannon proposed the concept of information entropy together with the theory of probability and statistics, which has been used as a parameter to measure the complexity of time series, and it can also be used to solve the problem of quantitative measurement [25].

If the variable X has n possible states , the probability , meets the following requirements:

The amount of information represented by each state is , where , which is known as the self-information of variable X. If there is high uncertainty in a specific state of X, then the probability of this state pattern will be small, and vice versa [26]. The information entropy can be defined as the expected value of the variable, as follows:

It can also be regarded as a probability distribution function of variable X. If the variable is orderly, the value of entropy is small. Conversely, if the variable is irregular, the value of entropy is high. Thus, the entropy function can be regarded as a parameter to measure the complexity of the time series, and it can also be used to solve the problem of quantitative measurement.

2.2. Strain Energy Calculation

The elastic strain energy of the concrete hinged slab element i is known from the theoretical basis of finite elements:where is the strain energy of the concrete hinged slab element i; is the number of unit integration points; is the stress vector; is the elastic strain vector; and is the volume of unit integration point k.

2.2.1. Modal Strain Energy

Taking the Euler–Bernoulli hinged slab as an example, the strain energy iswhere and L are the bending rigidity and length of the hinged slab, respectively. From (4), the modal strain energy of the jth modal vibration mode of a linear hinged slab with N units and m degrees of freedom is

Based on finite element method basic principle, the jth mode strain energy of the element i iswhere is the extended order of the element stiffness matrix in the global coordinate system.

2.2.2. Relationship between Displacement Mode and Strain Mode

The structural damage will change mechanical properties of the materials. This is reflected by the loss of local structural stiffness, resulting in a change of the MSE in the damaged region [13]. In this study, the hinged slab structure is modeled as a planar Euler–Bernoulli hinged slab with a constant cross section. For the ith hinged slab element, the elemental nodal displacement vector has the following form [27]:where the superscripts e and T stand for elemental and transpose, respectively; and are the nodal deflection of the and nodes, respectively; and and are the nodal rotation of the and nodes, respectively.

The deflection at local coordinate on the hinged slab element can be expressed aswhere is the shape function vector.

The strain of the local coordinate on each hinged slab element can be expressed aswhere represents the distance from the neutral axis of the hinged slab to the location of the strain gauge.

To develop a damage indicator, the displacement mode and rotation mode are replaced by the strain mode, and the nodal displacement vector of the ith hinged slab element is replaced by a nodal strain vector, including the strains on the structural surface of the ith hinged slab element, i.e., and . Therefore, the nodal strain vector is

The relationship between the nodal displacement vector and the nodal strain vector can be derived as follows:where is the elemental transfer matrix:

In modal strain analysis, the strain mode and displacement mode are coupled with the same transfer matrix . Therefore, the relationship between the displacement mode shape on the element and the strain mode shape on the ith element can be expressed as follows [28]:

2.2.3. Elemental MSE of Strain Mode

Substituting (10) into (12), the elemental MSE of the strain mode can be calculated as follows:where is defined as

From (11) through (16), the elemental before and after damage can be represented by the strain mode aswhere and are jth strain mode shape on the element before and after damage, respectively. Because the elemental stiffness of damaged elements is unknown, the undamaged elemental stiffness matrix is used as an approximation instead of the damaged one in .

3. Structural Damage Identification Index

3.1. Identification of Hinge Joint Damage

According to the basic theory of local information entropy [29, 30], when calculating the local probability of a hinge joint element i, the strain energy of the adjacent elements around i should be considered. Considering the strain energy of each element, , as the information amount, the local probability of element i at the jth hinge joint can be expressed as follows:where and are the deformation energies of the i and n elements at the jth hinge joint, respectively, and is the strain energy of element n on the concrete hinged slab.

For the hinge joint, a new variable is introduced to stand for the weight of the i element local probability.

During the calculation, the larger n is, the more data near support elements will be discarded. By substituting (19) into (2), the local energy information entropy of element i of the hinge joint can be expressed as

Subsequently, the damage degree of the i element at the jth hinge joint is calculated aswhere and are the local energy information entropy of element i of the hinge joint before and after damage, respectively.

3.2. Identification of Hinged Slab Damage

Based on (16) and (17), the modal strain ratio energies before and after the damage can be, respectively, defined as new physical quantities and :

Substituting (22) and (23) into (2), the modal strain ratio energy entropy functions and of element i before and after damage are, respectively, established as follows:

The entropy difference can be calculated as follows:

At the site of damage, the value of Si will change significantly compared to a site with no damage. Therefore, we can easily identify the location and severity of damage elements.

4. Finite Element Numerical Simulation

Due to the existence of hinge joints, the force transfer mechanism of a simply supported beam is different from that of other types of bridges. The hinge joint is usually in a complicated stress state with bending moment , shear , and compression . The hinge joint is a crucial component to ensure that the slabs work together, as shown in Figure 1.

Local strain energy information entropy and modal strain energy information entropy are used in this paper to identify the main damage in hinge joints and hinged slabs, respectively.

4.1. Mechanical Model of Hinge Joint

It is obvious that traditional hinge plate theory cannot define the actual force at the hinge joint because it only considers shear force. The force diagram in Figure 1 can be simplified, as shown in Figure 2. Since and mainly act as transverse normal stresses on the hinge joints, they are equivalent to a pair of equal and opposite transverse forces and .

Bearing shear action is the main function of hinge joints. However, due to the small size of the hinge joints, it is assumed that the shear force is evenly distributed at the hinge joints. Moreover, when the hinge joints are compressed due to the deformation of the concrete plates, the directions of the upper and lower parts of the hinge joints may be inconsistent; there are two kinds of forces: tension and compression. Therefore, the hinge joint is equivalent to two springs a and b, which only bear transverse tension or pressure, and a connecting rod system which only bears vertical shear, as shown in Figure 3.

For the entire beam, the calculation model of the hinge joints and hinged slabs can be obtained by segmenting and simulating the hinge joints of each section, as shown in Figure 3, with numbered elements, as shown in Figure 4.

4.2. Finite Element Analysis Model Establishment

The finite element analysis model was established in ANSYS, which consists of three pieces of concrete hinged slabs connected with the hinge joints 1 and 2. The width and height of each hinged slab are 950 mm and 800 mm, respectively.

The span is 10 m and is divided into 50 elements longitudinally. The width c and height h of the hinge joints are 75 mm and 700 mm, respectively. The concrete grade is C40. The hinge joints are composed of two springs subjected only to transverse forces and one spring that is subjected only to vertical forces. The spring stiffnesses are kabi and kgi, respectively, as shown in Figure 5.

Here, E is the elastic modulus of concrete; G is the shear modulus of concrete; le is the length of the partition section, le = 200; and α is the rectangular section of constant maturity, α = 1.5.

When the finite element model is loaded, a concentrated force Pi = 100 kN is applied at the center of the structure. However, to avoid the occurrence of a concentrated strain in the structure, according to the method in literature [31], a sinusoidal uniformly distributed load of a half-wave was used instead. The loaded structure is shown in Figure 6.

5. Results of Damage Identification Analysis

5.1. Analysis of Hinge Joint Damage

Through the above theory, the strain energy information entropy index value of each damage case is calculated in MATLAB programming environment.

In this paper, different damage cases were set, as shown in Table 1.

5.1.1. Single-Point Damage of Single-Hinge Joint

In the first case study, the damage is simulated by reducing the stiffness of the joint element, and the damage degree is preset at 5%. Figure 7 shows the analysis results. From the results, the damage position of the 26th hinge joint element of hinge joint 1 can be clearly seen, and the damage degree is 6.26%, with an identification error of 20.12%.

The results show that the hinge joint damage identification method developed in this paper can accurately identify the location of single-hinge joint damage, and the calculated damage degree is higher than the preset damage degree.

5.1.2. Multipoint Damage of Hinge Joints

Figures 8 and 9 show the results for the case studies 2, 3, and 4. The hinge joint damage identification method proposed in this paper is also effective for detecting the multipoint damage of hinge joints. It can be clearly seen from the results that the damage position is the same as the preset damage position. In calculation results of damage degree of case study 3, continuous damage appears in the 38th and 39th hinge joint elements, which has a significant impact on the damage identification result of the 26th hinge joint element, and the damage identification error increases.

When two hinge joints are damaged at the same time, the result of damage identification is shown in Figure 9. It can be seen that the damage index can accurately locate the damage position of the two hinge joints, but the calculated damage degree error is slightly larger.

Table 2 shows the calculated damage degree values.

5.2. Analysis of Hinged Slab Damage

Hinged slab 1 was removed to identify local damage. Similarly, the elastic modulus of the structure was reduced to simulate damage. In this section, the first three modes of the structure are evaluated.

The noise interference is inevitable in actual engineering, so it is necessary to verify the antinoise performance of the proposed method. Theoretically, it is often assumed that the variation of errors follows a Gaussian distribution, but the Gaussian distribution contradicts the boundedness of actual errors. Therefore, the truncated Gaussian distribution (TGD) error model is introduced [32], and the model of random variables subject to the Gaussian distribution is truncated. The simulation error is introduced as follows:where and represent the mode shape components with and without noise error, respectively; is the noise level. are random variables subject to TGD; the truncated upper and lower limits are taken as −1 and 1, respectively. Table 3 defines each damage case.

For case studies 5 and 6, element 25 is assumed to have damage, represented by reductions in the elasticity modulus equal to 10% and 30%, respectively.

The results of damage detection shown in Figure 10 indicate that singular damage can be identified from the strain energy data with the proposed method successfully. When 10% and 30% damage occurred to the structure, the damage index values were 0.0075 and 0.029, respectively. With this linear relationship, the severity of the damage can be evaluated by calibration.

Case studies 7 and 8 were used to study double-damage cases. For case studies 7 and 8, elements 10 and 35 are assumed to have a stiffness reduction equal to 10%, 10%, and 10%, 30%, respectively. Meanwhile, the noise of 5% and 10% was introduced into for case study 7 to verify the antinoise performance of the damage index.

The results of damage detection are shown in Figures 11 and 12. The damage location can be accurately identified, and the severity of damage can be evaluated with the linear relationship between the values of the peak index and extent of the damage.

Figure 11 shows that, for case study 7, 10% damage occurred in elements 10 and 35 with damage indexes of 0.0095 and 0.0093, respectively. Note that the difference is small. For case study 8, 10% and 30% damage occurred in elements 10 and 35, respectively. The corresponding damage indexes are 0.0086 and 0.028, and the linear relationship of the damage index approaches the true extent of the damage.

Figure 12 shows that the proposed damage index can still locate the damage position and judge the damage degree well under noise intensities of 5% and 10%.

In order to verify the effectiveness of the index in more complex damage cases, different damage cases are set in elements 5, 15, 25, and 35 of the hinged slab, respectively; the damage identification results are shown in Figure 13. For case study 11, the damage degree of the damage position was all 20%, and the damage index values were 0.0201, 0.0192, 0.0213, and 0.0189, respectively. For case study 12, the damage degree of elements increased from 10% to 40% in turn, and the corresponding index values were 0.011, 0.023, 0.032, and 0.044, which basically showed a linear increase; the damage relationship curve is shown in Figure 14. It can be seen that the damage index is sensitive to all damage locations and can well reflect the damage degree of the structure.

6. Conclusions

In this paper, the strain energy entropy index was proposed by combining the strain energy and information entropy theories. Based on the new index used in damage identification of a simply supported beam, the following conclusions can be drawn.

According to the finite element analysis and calculation, the proposed method can locate well single-point and multipoint damage of hinge joints. Moreover, the hinge joint damage identification method is less affected by other hinge joints with this method, showing strong stability. When calculating the extent of hinge joint damage, the calculated value is often higher than the set value. This identification error is the next topic to study in future works.

The modal strain energy can accurately identify the damage location of the hinged slab and effectively judge the extent of the damage by the peak value of the index, and it is effective despite the noise. However, the stress form of the simply supported beam is relatively complex. Therefore, the practicability of the new index for this kind of structure should be further improved.

The stress analysis and hinge joint simulation of a simply supported beam should be further explored. For accurate calculation of the extent of hinge joint damage, artificial intelligence methods are promising choices, such as artificial neural network (ANN). The damage identification method based on artificial intelligence well be the focus of future works.

Data Availability

All data are available within the article or can be obtained from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Key Natural Science Project of Anhui Provincial Education Department (Grant no. KJ2019A0746), Natural Science Foundation of China (Grant no. 51868045), Natural Science Foundation Youth Fund of Anhui Province (Grant no. 2008085QE247), and Startup Fund for Doctor of Anhui Jianzhu University (Grant no. 2019QDZ08).