# An Efficient Damage Identification Method for Simply Supported Beams Based on Strain Energy Information Entropy

Authors: Chang-Sheng Xiang; Ling-Yun Li; Yu Zhou; Cong Dang
Journal: Advances in Materials Science and Engineering
 (2020)

---

## 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%.


---

## Body

## 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 [1–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–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%. 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 [9–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. 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 variableX has n possible states xi,i=1,2,…,n, the probability pxi,i=1,2,…,n, meets the following requirements:(1)0≤pxi≤1.∑i=1npxi=1.The amount of information represented by each statexi is Ixi=−lnpxi, 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]. The information entropy can be defined as the expected value of the variable, as follows:(2)HX=EIxi=−∑i=1npxilnpxi.It can also be regarded as a probability distribution function of variableX. 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 elementi is known from the theoretical basis of finite elements:(3)Wsi=12∑k=1nintσeTεevolk,where Wsi is the strain energy of the concrete hinged slab element i; nint is the number of unit integration points; σe is the stress vector; εe is the elastic strain vector; and volk 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 is(4)SE=12∫0LEIx∂2wx∂x22dx,where EIx 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 ϕjxi of a linear hinged slab with N units and m degrees of freedom is(5)MSEj=12∫0LEIx∂2ϕjx∂x22dx.Based on finite element method basic principle, thejth mode strain energy of the element i is(6)MSEij=12ϕjTkiϕj,where ki is the extended order of the element stiffness matrix kic 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 de has the following form [27]:(7)de=viθivi+1θi+1T,where the superscripts e and T stand for elemental and transpose, respectively; vi and vi+1 are the nodal deflection of the ith and i+1th nodes, respectively; and θi and θi+1 are the nodal rotation of the ith and i+1th nodes, respectively.The deflectionwxi at local coordinate xi on the hinged slab element can be expressed as(8)wxi=qde,where q is the shape function vector.The strainεxi of the local coordinate xi on each hinged slab element can be expressed as(9)εxi=±hd2wxidx2=±hd2qdx2de,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 theith 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. Therefore, the nodal strain vector ε is(10)ε=εiεi+1T.The relationship between the nodal displacement vectorde and the nodal strain vector ε can be derived as follows:(11)de=Tε=Tεiεi+1T,where T is the elemental transfer matrix:(12)T=13hl2+4l26l+2l2−l2−6l−l3−l26l+l3l2−6l−2l3.In modal strain analysis, the strain mode and displacement mode are coupled with the same transfer matrixT. Therefore, 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]:(13)ϕij=Tϕijε.


#### 2.2.3. Elemental MSE of Strain Mode

Substituting (10) into (12), the elemental MSE of the strain mode can be calculated as follows:(14)MSEij=12ϕijεTTkiTϕijε,=12ϕijεKiϕijε,where Ki is defined as(15)Ki=TTkiT=2lEI3h22112.From (11) through (16), the elemental MSE before and after damage can be represented by the strain mode as(16)MSEuij=∑j=1mϕuijεTKiϕuijε.(17)MSEdij=∑j=1mϕdijεTKiϕdijε,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 Ki is used as an approximation instead of the damaged one in MSEdij.


## 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 variableX has n possible states xi,i=1,2,…,n, the probability pxi,i=1,2,…,n, meets the following requirements:(1)0≤pxi≤1.∑i=1npxi=1.The amount of information represented by each statexi is Ixi=−lnpxi, 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]. The information entropy can be defined as the expected value of the variable, as follows:(2)HX=EIxi=−∑i=1npxilnpxi.It can also be regarded as a probability distribution function of variableX. 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 elementi is known from the theoretical basis of finite elements:(3)Wsi=12∑k=1nintσeTεevolk,where Wsi is the strain energy of the concrete hinged slab element i; nint is the number of unit integration points; σe is the stress vector; εe is the elastic strain vector; and volk 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 is(4)SE=12∫0LEIx∂2wx∂x22dx,where EIx 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 ϕjxi of a linear hinged slab with N units and m degrees of freedom is(5)MSEj=12∫0LEIx∂2ϕjx∂x22dx.Based on finite element method basic principle, thejth mode strain energy of the element i is(6)MSEij=12ϕjTkiϕj,where ki is the extended order of the element stiffness matrix kic 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 de has the following form [27]:(7)de=viθivi+1θi+1T,where the superscripts e and T stand for elemental and transpose, respectively; vi and vi+1 are the nodal deflection of the ith and i+1th nodes, respectively; and θi and θi+1 are the nodal rotation of the ith and i+1th nodes, respectively.The deflectionwxi at local coordinate xi on the hinged slab element can be expressed as(8)wxi=qde,where q is the shape function vector.The strainεxi of the local coordinate xi on each hinged slab element can be expressed as(9)εxi=±hd2wxidx2=±hd2qdx2de,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 theith 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. Therefore, the nodal strain vector ε is(10)ε=εiεi+1T.The relationship between the nodal displacement vectorde and the nodal strain vector ε can be derived as follows:(11)de=Tε=Tεiεi+1T,where T is the elemental transfer matrix:(12)T=13hl2+4l26l+2l2−l2−6l−l3−l26l+l3l2−6l−2l3.In modal strain analysis, the strain mode and displacement mode are coupled with the same transfer matrixT. Therefore, 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]:(13)ϕij=Tϕijε.


### 2.2.3. Elemental MSE of Strain Mode

Substituting (10) into (12), the elemental MSE of the strain mode can be calculated as follows:(14)MSEij=12ϕijεTTkiTϕijε,=12ϕijεKiϕijε,where Ki is defined as(15)Ki=TTkiT=2lEI3h22112.From (11) through (16), the elemental MSE before and after damage can be represented by the strain mode as(16)MSEuij=∑j=1mϕuijεTKiϕuijε.(17)MSEdij=∑j=1mϕdijεTKiϕdijε,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 Ki is used as an approximation instead of the damaged one in MSEdij.


## 2.2.1. Modal Strain Energy

Taking the Euler–Bernoulli hinged slab as an example, the strain energy is(4)SE=12∫0LEIx∂2wx∂x22dx,where EIx 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 ϕjxi of a linear hinged slab with N units and m degrees of freedom is(5)MSEj=12∫0LEIx∂2ϕjx∂x22dx.Based on finite element method basic principle, thejth mode strain energy of the element i is(6)MSEij=12ϕjTkiϕj,where ki is the extended order of the element stiffness matrix kic 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 de has the following form [27]:(7)de=viθivi+1θi+1T,where the superscripts e and T stand for elemental and transpose, respectively; vi and vi+1 are the nodal deflection of the ith and i+1th nodes, respectively; and θi and θi+1 are the nodal rotation of the ith and i+1th nodes, respectively.The deflectionwxi at local coordinate xi on the hinged slab element can be expressed as(8)wxi=qde,where q is the shape function vector.The strainεxi of the local coordinate xi on each hinged slab element can be expressed as(9)εxi=±hd2wxidx2=±hd2qdx2de,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 theith 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. Therefore, the nodal strain vector ε is(10)ε=εiεi+1T.The relationship between the nodal displacement vectorde and the nodal strain vector ε can be derived as follows:(11)de=Tε=Tεiεi+1T,where T is the elemental transfer matrix:(12)T=13hl2+4l26l+2l2−l2−6l−l3−l26l+l3l2−6l−2l3.In modal strain analysis, the strain mode and displacement mode are coupled with the same transfer matrixT. Therefore, 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]:(13)ϕij=Tϕijε.


## 2.2.3. Elemental MSE of Strain Mode

Substituting (10) into (12), the elemental MSE of the strain mode can be calculated as follows:(14)MSEij=12ϕijεTTkiTϕijε,=12ϕijεKiϕijε,where Ki is defined as(15)Ki=TTkiT=2lEI3h22112.From (11) through (16), the elemental MSE before and after damage can be represented by the strain mode as(16)MSEuij=∑j=1mϕuijεTKiϕuijε.(17)MSEdij=∑j=1mϕdijεTKiϕdijε,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 Ki is used as an approximation instead of the damaged one in MSEdij.


## 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, W, as the information amount, the local probability Pij of element i at the jth hinge joint can be expressed as follows:(18)Pij=Wij∑m=j−ij∑n=i−1iWsnm+∑n=i−1i+1Wnj,where Wij and Wnj are the deformation energies of the i and n elements at the jth hinge joint, respectively, and Wsnm is the strain energy of element n on the mth concrete hinged slab.For thejth hinge joint, a new variable λij is introduced to stand for the weight of the i element local probability.(19)λij=Pij∑k=1nPkj∑k=1nλkj=1.During the calculation, the largern is, the more data near support elements will be discarded. By substituting (19) into (2), the local energy information entropy Hij of element i of the jth hinge joint can be expressed as(20)Hij=pij∑k=1nPkjln∑k=1nPkj−lnPij.Subsequently, the damage degreeDij of the i element at the jth hinge joint is calculated as(21)Dij=0Hi,uj≤Hi,dj1−Hi,djHi,ujHi,uj>Hi,dj,where Hi,uj and Hi,dj are the local energy information entropy of element i of the jth 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 λui and λdi:(22)λui=MSEuij∑i=1nMSEuij.(23)λui=MSEuij∑i=1nMSEuij.Substituting (22) and (23) into (2), the modal strain ratio energy entropy functions Sui and Sdi of element i before and after damage are, respectively, established as follows:(24)Sui=−∑i=1nλuilnλui,=−∑i=1nMSEuij∑i=1nMSEuijlnMSEuij∑i=1nMSEuij,Sdi=−∑i=1nλdilnλdi,=−∑i=1nMSEdij∑i=1nMSEdijlnMSEdij∑i=1nMSEdij.The entropy difference can be calculated as follows:(25)Si=Sdi−Sui.At the site of damage, the value ofSi will change significantly compared to a site with no damage. Therefore, we can easily identify the location and severity of damage elements.


## 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, W, as the information amount, the local probability Pij of element i at the jth hinge joint can be expressed as follows:(18)Pij=Wij∑m=j−ij∑n=i−1iWsnm+∑n=i−1i+1Wnj,where Wij and Wnj are the deformation energies of the i and n elements at the jth hinge joint, respectively, and Wsnm is the strain energy of element n on the mth concrete hinged slab.For thejth hinge joint, a new variable λij is introduced to stand for the weight of the i element local probability.(19)λij=Pij∑k=1nPkj∑k=1nλkj=1.During the calculation, the largern is, the more data near support elements will be discarded. By substituting (19) into (2), the local energy information entropy Hij of element i of the jth hinge joint can be expressed as(20)Hij=pij∑k=1nPkjln∑k=1nPkj−lnPij.Subsequently, the damage degreeDij of the i element at the jth hinge joint is calculated as(21)Dij=0Hi,uj≤Hi,dj1−Hi,djHi,ujHi,uj>Hi,dj,where Hi,uj and Hi,dj are the local energy information entropy of element i of the jth 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 λui and λdi:(22)λui=MSEuij∑i=1nMSEuij.(23)λui=MSEuij∑i=1nMSEuij.Substituting (22) and (23) into (2), the modal strain ratio energy entropy functions Sui and Sdi of element i before and after damage are, respectively, established as follows:(24)Sui=−∑i=1nλuilnλui,=−∑i=1nMSEuij∑i=1nMSEuijlnMSEuij∑i=1nMSEuij,Sdi=−∑i=1nλdilnλdi,=−∑i=1nMSEdij∑i=1nMSEdijlnMSEdij∑i=1nMSEdij.The entropy difference can be calculated as follows:(25)Si=Sdi−Sui.At the site of damage, the value ofSi 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 momentmx, shear gx, and compression cx. The hinge joint is a crucial component to ensure that the slabs work together, as shown in Figure 1.Figure 1
Force diagram of simply supported beam.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 Figure1 can be simplified, as shown in Figure 2. Since mx and cx mainly act as transverse normal stresses on the hinge joints, they are equivalent to a pair of equal and opposite transverse forces nax and nbx.Figure 2
Simplification of the forces at hinge joints.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 springsa and b, which only bear transverse tension or pressure, and a connecting rod system which only bears vertical shear, as shown in Figure 3.Figure 3
Hinge joint model.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 Figure3, with numbered elements, as shown in Figure 4.Figure 4
Structure segmentation and unit numbering.


### 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 widthc 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.(26)kabi=12cEhle.kgi=12αcGh.Figure 5
Finite element analysis model.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 px of a half-wave was used instead. The loaded structure is shown in Figure 6.(27)px=2piLsinπxL.Figure 6
Finite element analysis model loading.


## 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 Figure1 can be simplified, as shown in Figure 2. Since mx and cx mainly act as transverse normal stresses on the hinge joints, they are equivalent to a pair of equal and opposite transverse forces nax and nbx.Figure 2
Simplification of the forces at hinge joints.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 springsa and b, which only bear transverse tension or pressure, and a connecting rod system which only bears vertical shear, as shown in Figure 3.Figure 3
Hinge joint model.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 Figure3, with numbered elements, as shown in Figure 4.Figure 4
Structure segmentation and unit numbering.


## 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 widthc 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.(26)kabi=12cEhle.kgi=12αcGh.Figure 5
Finite element analysis model.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 px of a half-wave was used instead. The loaded structure is shown in Figure 6.(27)px=2piLsinπxL.Figure 6
Finite element analysis model loading.


## 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 Table1.Table 1
Damage cases in hinge joints.
Damage caseDamage element no.Damage severity (%)1eh26152eh2615eh381103eh2615eh38110eh391104eh26110eh38115eh1225Note. eh261 represents element 26 of the first hinge joint.


#### 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%. Figure7 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%.Figure 7
Calculation results of damage degree of hinge joints for case study 1.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

Figures8 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.Figure 8
Calculation results of damage degree of hinge joints. (a) Case study 2. (b) Case study 3.
(a)(b)Figure 9
Calculation results of damage degree of hinge joints for case study 4.When two hinge joints are damaged at the same time, the result of damage identification is shown in Figure9. 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.Table2 shows the calculated damage degree values.Table 2
Damage degree values for multipoint damage case studies.
Damage caseDamage element no.Damage degree values (%)Identification error (%)2eh2616.4422.36eh38112.0817.213eh2617.1429.97eh38112.4119.41eh39113.0623.434eh26112.8825.79eh38117.4019.45eh1227.4224.56


### 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:(28)Φie=Φi1+u⋅TGD−1,1,where Φi and Φie 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.Table 3
Damage case studies in the hinged slab.
Damage caseDamage element no.Damage severity (%)5251062530710, 3510810, 3530910, 3530 (5% noise)105, 15, 25, 3530 (10% noise)115, 15, 25, 3520, 20, 20, 20125, 15, 25, 3510, 20, 30, 40135, 15, 25, 3540, 10, 20, 10Note. “10, 30 (5% noise)” represents an addition of 5% noise error in the damage case study.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 Figure10 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.Figure 10
Damage detection results of single damage cases. (a) Case study 5. (b) Case study 6.
(a)(b)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 Figures11 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
Damage detection results of double-damage cases. (a) Case study 7. (b) Case study 8.
(a)(b)Figure 12
Antinoise verification of damage index.Figure11 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.Figure12 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 Figure13. 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.Figure 13
Damage detection results of cases 10–12.Figure 14
Damage relation curve of case 11.


## 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 Table1.Table 1
Damage cases in hinge joints.
Damage caseDamage element no.Damage severity (%)1eh26152eh2615eh381103eh2615eh38110eh391104eh26110eh38115eh1225Note. eh261 represents element 26 of the first hinge joint.


### 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%. Figure7 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%.Figure 7
Calculation results of damage degree of hinge joints for case study 1.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

Figures8 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.Figure 8
Calculation results of damage degree of hinge joints. (a) Case study 2. (b) Case study 3.
(a)(b)Figure 9
Calculation results of damage degree of hinge joints for case study 4.When two hinge joints are damaged at the same time, the result of damage identification is shown in Figure9. 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.Table2 shows the calculated damage degree values.Table 2
Damage degree values for multipoint damage case studies.
Damage caseDamage element no.Damage degree values (%)Identification error (%)2eh2616.4422.36eh38112.0817.213eh2617.1429.97eh38112.4119.41eh39113.0623.434eh26112.8825.79eh38117.4019.45eh1227.4224.56


## 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%. Figure7 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%.Figure 7
Calculation results of damage degree of hinge joints for case study 1.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

Figures8 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.Figure 8
Calculation results of damage degree of hinge joints. (a) Case study 2. (b) Case study 3.
(a)(b)Figure 9
Calculation results of damage degree of hinge joints for case study 4.When two hinge joints are damaged at the same time, the result of damage identification is shown in Figure9. 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.Table2 shows the calculated damage degree values.Table 2
Damage degree values for multipoint damage case studies.
Damage caseDamage element no.Damage degree values (%)Identification error (%)2eh2616.4422.36eh38112.0817.213eh2617.1429.97eh38112.4119.41eh39113.0623.434eh26112.8825.79eh38117.4019.45eh1227.4224.56


## 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:(28)Φie=Φi1+u⋅TGD−1,1,where Φi and Φie 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.Table 3
Damage case studies in the hinged slab.
Damage caseDamage element no.Damage severity (%)5251062530710, 3510810, 3530910, 3530 (5% noise)105, 15, 25, 3530 (10% noise)115, 15, 25, 3520, 20, 20, 20125, 15, 25, 3510, 20, 30, 40135, 15, 25, 3540, 10, 20, 10Note. “10, 30 (5% noise)” represents an addition of 5% noise error in the damage case study.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 Figure10 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.Figure 10
Damage detection results of single damage cases. (a) Case study 5. (b) Case study 6.
(a)(b)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 Figures11 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
Damage detection results of double-damage cases. (a) Case study 7. (b) Case study 8.
(a)(b)Figure 12
Antinoise verification of damage index.Figure11 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.Figure12 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 Figure13. 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.Figure 13
Damage detection results of cases 10–12.Figure 14
Damage relation curve of case 11.


## 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.



---

*Source: 9283949-2020-09-07.xml*