An Efficient Damage IdentificationMethod for Simply Supported Beams Based on Strain Energy Information Entropy

Concrete simply supported beams arewidely used in smallandmedium-sized bridges.Hinged slabs and joints are damaged in service easily, as they are themain load-carryingmembers of bridge. In this paper, using thehigh 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 supportedbeammodel is established to verify the effectiveness of theproposed index.+enumerical 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. +e 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%.


Introduction
e 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-andgroove concrete joints to bear vehicle loads [1][2][3]. 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 [4][5][6]. 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%. erefore, 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 [9][10][11], 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. e 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. e 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. erefore, 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.

Information Entropy eory.
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 x i , i � 1, 2, . . . , n, the probability p(x i ), i � 1, 2, . . . , n, meets the following requirements: (1) e amount of information represented by each state x i is I(x i ) � −ln p(x i ), where i � 1, 2, . . . , n, 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]. e 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. us, 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.

Strain Energy Calculation.
e elastic strain energy of the concrete hinged slab element i is known from the theoretical basis of finite elements: where W si is the strain energy of the concrete hinged slab element i; n int is the number of unit integration points; σ e { } is the stress vector; ε e { } is the elastic strain vector; and vol k is the volume of unit integration point k.

Modal Strain Energy.
Taking the Euler-Bernoulli hinged slab as an example, the strain energy is where EI(x) 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 ϕ j (x i ) 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 is where k i is the extended order of the element stiffness matrix k c i in the global coordinate system.

2
Advances in Materials Science and Engineering

Relationship between Displacement Mode and Strain
Mode. e structural damage will change mechanical properties of the materials. is 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 d { } e has the following form [27]: where the superscripts e and T stand for elemental and transpose, respectively; v i and v i+1 are the nodal deflection of the ith and (i + 1)th nodes, respectively; and θ i and θ i+1 are the nodal rotation of the ith and (i + 1)th nodes, respectively. e deflection w(x i ) at local coordinate x i on the hinged slab element can be expressed as where q is the shape function vector. e strain ε(x i ) of the local coordinate x i on each hinged slab element can be expressed as where h 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., ε i and ε i+1 . erefore, the nodal strain vector ε e relationship between the nodal displacement vector d { } e and the nodal strain vector ε { } can be derived as follows: where [T] is the elemental transfer matrix: In modal strain analysis, the strain mode and displacement mode are coupled with the same transfer matrix [T]. erefore, the relationship between the jth displacement mode shape on the ith element ϕ ij and the jth strain mode shape on the ith element ϕ ε ij can be expressed as follows [28]:

Elemental MSE of Strain
Mode. Substituting (10) into (12), the elemental MSE of the strain mode can be calculated as follows: where [K i ] is defined as From (11) through (16), the elemental MSE before and after damage can be represented by the strain mode as where ϕ ε uij and ϕ ε dij are jth strain mode shape on the ith element before and after damage, respectively. Because the elemental stiffness of damaged elements is unknown, the undamaged elemental stiffness matrix [K i ] is used as an approximation instead of the damaged one in MSE dij .

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, W, as the information amount, the local probability P j i of element i at the jth hinge joint can be expressed as follows: where W j i and W j n are the deformation energies of the i and n elements at the jth hinge joint, respectively, and W m sn is the strain energy of element n on the m th concrete hinged slab.
For the jth hinge joint, a new variable λ j i 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 H j i of element i of the jth hinge joint can be expressed as Subsequently, the damage degree D j i of the i element at the jth hinge joint is calculated as where H j i,u and H j i,d are the local energy information entropy of element i of the jth hinge joint before and after damage, respectively. (16) and (17), the modal strain ratio energies before and after the damage can be, respectively, defined as new physical quantities λ ui and λ di :

Identification of Hinged Slab Damage. Based on
Substituting (22) and (23) into (2), the modal strain ratio energy entropy functions S ui and S di of element i before and after damage are, respectively, established as follows: e entropy difference can be calculated as follows: At the site of damage, the value of S i will change significantly compared to a site with no damage. erefore, we can easily identify the location and severity of damage elements.

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. e hinge joint is usually in a complicated stress state with bending moment m(x), shear g(x), and compression c(x). e 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.

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. e force diagram in Figure 1 can be simplified, as shown in Figure 2. Since m(x) and c(x) mainly act as transverse normal stresses on the hinge joints, they are equivalent to a pair of equal and opposite transverse forces na(x) and nb(x).
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. erefore, 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.

Finite Element Analysis Model Establishment.
e 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. e width and height of each hinged slab are 950 mm and 800 mm, respectively. e span is 10 m and is divided into 50 elements longitudinally. e width c and height h of the hinge joints are 75 mm and 700 mm, respectively. e concrete grade is C40. e hinge joints are composed of two springs subjected only to transverse forces and one spring that is subjected only to vertical forces. e spring stiffnesses are k abi and k gi , respectively, as shown in Figure 5. Gh.

(26)
Here, E is the elastic modulus of concrete; G is the shear modulus of concrete; l e is the length of the partition section, l e � 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 4 Advances in Materials Science and Engineering sinusoidal uniformly distributed load p(x) of a half-wave was used instead. e loaded structure is shown in Figure 6.

Analysis of Hinge Joint Damage.
rough the above theory, the strain energy information entropy index value of    Advances in Materials Science and Engineering each damage case is calculated in MATLAB programming environment.
In this paper, different damage cases were set, as shown in Table 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   Advances in Materials Science and Engineering 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%. e 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. Figures 8 and 9 show the results for the case studies 2, 3, and 4. e 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.

Multipoint Damage of Hinge Joints.
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.

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. e noise interference is inevitable in actual engineering, so it is necessary to verify the antinoise performance of the proposed method. eoretically, it is often assumed that the variation of errors follows a Gaussian distribution, but the Gaussian distribution contradicts the boundedness of actual errors. erefore, the truncated Gaussian distribution (TGD) error model is introduced [32], and the model of random variables subject to the Gaussian distribution is truncated. e simulation error is introduced as follows: where Φ i and Φ e i represent the mode shape components with and without noise error, respectively; u is the noise level. TGD(−1, 1) 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. e 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. e results of damage detection are shown in Figures 11  and 12. e damage location can be accurately identified, and the severity of damage can be evaluated with the linear   Note. "10, 30 (5% noise)" represents an addition of 5% noise error in the damage case study. 8 Advances in Materials Science and Engineering 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. e 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   Advances in Materials Science and Engineering 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.

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. is identification error is the next topic to study in future works. e 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. erefore, the practicability of the new index for this kind of structure should be further improved.
e 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). e 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
e authors declare that there are no conflicts of interest regarding the publication of this paper.