Damage Identification by the Data Expansion and Substructuring Methods

Structural damage can be detected by comparing the responses before and after the damage. /e responses are transformed into curvature, strain, and stress, among others, which characterize the mechanical behavior of the structural members, and can be utilized as damage indices for damage detection. /e damage of a truss structure can rarely be detected by the displacements only at nodes. /is work investigates damage detection methods using the stress or stiffness variation rate of the truss element before and after the damage. /is paper considers three different cases according to the number of measurement locations. If the complete responses at a full set of degrees of freedom are measured, the stiffness variation rates of the elements are calculated accurately, and the damage can be explicitly detected despite external noise. If the number of measured data points is fewer than the system order, the displacements are estimated by the data expansion method, and the damage-expected regions are predicted by the stiffness variation rates. Apart from the explicitly damaged elements, the substructuring approach is adopted for closer damage detection with several measurement sensors despite external noise. It is illustrated by the examples that three cases are compared numerically. /e numerical examples compare and analyze the numerical results of the three cases.


Introduction
Structural maintenance has received significant attention in the construction industry.Existing structures have been diagnosed, repaired, reformed, and remodeled for improving the structural performance and durability.Structural health monitoring and nondestructive tests for evaluating structural performance have been developed with the advent of new measurement sensors and technologies, and techniques for structural health monitoring have been applied in the field.
e health state of a structure is evaluated by the response, internal strain, or internal stress, among others.Structural defects are detected in the process by comparing the responses before and after the structural performance deterioration.e measured responses are transformed into curvature, strain, and stress for damage detection.e measured displacements in the truss structure should be transformed into axial stresses for evaluating the existence of damage.
It is impractical to collect the responses at a full set of degrees of freedom (DOFs) because of the limited number of sensors.e number of measurement sensors is practically fewer than the system order.us, incomplete measurement data should be expanded, or the system order should be reduced.
e damage identification method has become a prominent research field.
e identification of a global structure requires a large number of sensors and unknowns; however, the substructuring method needs only several sensors for response measurement.
Blachowski et al. [1] presented a damage quantification method for the truss structure without identifying modal parameters or solving a global optimization problem.e method quantifies the damage in a given substructure with only a small subset of sensors.Park and Reich [2] provided two complementary methods for model-based structural damage detection, based on changes in the localized flexibility and invariance properties of the elemental or substructural transmission zeros.Xing and Mita [3] proposed a substructure approach for identifying localized damage on each substructure using only three sensors.Ugalde et al. [4] presented a substructuring approach to predict the health of each substructure by means of the vector autoregressive variation with the exogenous model.Hou et al. [5] proposed a substructure damage identification method using local free responses.Estimating additional unknown inputs by the algorithm, without measurements of the substructure interface DOFs, Lei et al. [6] applied the damage detection method by means of the degradation of the identified substructural element stiffness values.Lee and Eun [7] presented a modelbased substructuring method for local damage detection using the frequency response function measured at the interfaces between the adjacent substructures.Li et al. [8] proposed a substructural damage identification approach based on dynamic response reconstruction and a dynamic response sensitivity-based method.
An et al. [9] presented a stochastic damage-locating vector method as a useful tool for damage localization in steel truss bridges.An et al. [10] proposed a damage localization approach for truss structures using the curvature difference method of strain waveform fractal dimensions.Yang and Jin [11] presented a damage detection method for truss structures using incompletely measured modes, with the aid of the eigenvalue perturbation method.Kaven and Zolghadr [12] proposed a damage detection method using the differences between structure characteristics, such as natural frequencies and mode shapes of truss structures.Xu and Wu [13] developed a damage detection technique based on strain data under ambient excitation using the environmental excitation incomplete strain mode.Artar and Daloglu [14] presented a damage identification method to update the numerical model and to predict damage location and severity in the space frame using genetic algorithms.
is study considers methods for predicting damaged elements depending on complete and incomplete measurement data, including the noise effect.
e measured and estimated displacements are transformed into the internal stress of the elements, and the stiffness variation rate is used as a damage index.By dividing a full truss structure into several damage-expected substructures, this work proposes the substructuring method to estimate damaged elements.is method offers merits in reducing the model order and number of measurement locations by the static-based and dynamic-based data expansion method [15] and is insensitive to external noise.
e damage detection methods using complete and incomplete measurements are investigated in numerical examples.e validity of the proposed static-based and dynamic-based substructuring method is illustrated in examples, depending on external noise.

Substructuring Method for Damage Detection
Structural diagnosis is conducted by means of the comparison of the response, strain, stress, and stiffness in healthy and unhealthy states, or the measurement data only in the damage-expected state without the baseline data. is method is performed by a response comparison in both states.is work compares three damage detection methods, depending on the number of measurement locations to be causable in the structural performance evaluation.e three different cases are as follows: (1) damage detection using complete measurement data; (2) damage detection using the data to expand fewer measurement data than the model order; and (3) damage detection by measurement data and data expansion at the substructures to divide an entire structure.Moreover, the effect of the noise contained in the measurement data is examined.Figure 1 summarizes the flow of this study to detect damaged elements.
where  K and  u denote the n × n positive definite stiffness matrix and n × 1 displacement vector in the healthy state, respectively.Moreover, F is the n × 1 external force vector.
e displacement data at a full set of DOFs may be measured by measurement sensors or calculated by the finite element model for damage detection.
e truss structure is composed of axial members for carrying the axial load, and each member has axial displacement DOFs at both end nodes.
e measured displacements are transformed into the axial stress related to the axial stiffness, and the damage can be detected by comparing the stresses or stiffness.is complete method exhibits the limitation of requiring the same number of measurement sensors as the system order.

Incomplete Measurement Data.
It is not practical to measure the responses for the entire DOFs because of fewer measurement sensors than the system order.is indicates that the full set of displacements should be estimated by means of data expansion.Assuming that the static responses are taken as m(m < n) DOFs, constraints must be implemented to describe the static behavior of the entire system.
ese can be written as where the m × n coefficient matrix A is a Boolean matrix for defining the measurement locations and the vectors u and b are the n × 1 actual displacement vector and m × 1 measured response vector, respectively.Constraint forces for pulling the unconstrained path of the initial system into the constrained path are required.Combining equations ( 1) and (2), reference [15] derives the constrained equilibrium equation as follows: where "+" denotes the Moore-Penrose inverse.e displacements corresponding to the entire DOFs are estimated 2 Advances in Civil Engineering by using equation (3), satisfying the constraints of equation ( 2). e second term on the right-hand side of equation ( 3) represents the displacement variation deviating from the initially established trajectory.Moreover, by premultiplying the second term on the right-hand side of equation ( 3) by  K, the constraint force vector F c is obtained by e internal stresses of the truss structure can be calculated accurately by using the displacements and axial stiffness.However, the displacements for expanding the incomplete measurements rarely coincide with the accurate displacements.Moreover, the calculated internal stresses of the truss elements do not coincide with the exact values.e following example investigates the feasibility of the damage detection by comparing the internal stresses estimated by the expanded data with the initially established stresses.

Example 1.
e adequacy of the estimated responses using equation ( 3) is investigated in the numerical example of a plane truss structure, as illustrated in Figure 2. In the figure, the nodal points and members are numbered.Each node has two DOFs of the horizontal and vertical responses u and v, respectively.
e truss is composed of six nodes, nine members, and nine DOFs, except for boundary conditions.All members have the same elastic modulus of 200 GPa and cross-sectional area of 2.5 ×10 −3 m 2 .e simply supported truss has a single span.Its length is 12 m, its height is 3 m, and each bay is 4 m long.e responses in the intact state can be calculated by the finite element method under the action of the external force of 10 N in the downward direction of node 4. e elemental stiffness matrices and responses in the initial state should be saved for subsequent analysis.
e truss exhibits a 20% section loss at element ④ and 30% section loss at element ⑧. e vertical displacements v 2 , v 3 , v 4 , and v 5 , numerically simulated in the damaged state, were utilized as measurement data.Inserting the equilibrium equation in the intact state and measurement data in the damaged state into equation ( 3), the displacements for the entire DOFs are estimated.e truss structure is composed of axial members and pin joints.us, unlike the beam, the truss structure represents discontinuous responses at the element nodes.e differences between the actual and estimated responses are displayed in Figure 3. e plots illustrate that the resulting responses coincide at measured DOFs, but a discrepancy of a very small magnitude is evident at the other DOFs.e expansion method does not describe the accurate responses, except for the constraint conditions.
e utilization of inaccurately expanded data leads to incorrect analysis.Advances in Civil Engineering e stress variation rate α, defined as the stress variation with respect to the stress on an element in the damage state, can be calculated by where σ u and σ d denote the axial stresses in the undamaged and damaged states, respectively.Moreover, σ d is calculated using the cross-sectional area and axial stiffness in the intact state and measurement data in the damaged state.e stress variation rates in Figures 4(a) and 4(b) are obtained using the complete and incomplete measurement data, respectively.It is observed that accurate stress variation rates for the damaged elements can be obtained by using the complete measurement data.However, the plot using the incomplete measurement data exhibits very small variation rates, regardless of the damaged elements, and the damaged elements cannot be detected.
is expansion method exhibits the drawbacks of inaccurately estimating the expanded data and yielding inconsistent results depending on the sensor location and number.us, this work considers the analytical method for detecting the damage by dividing an entire structure into substructures in order to reduce the sensor number and detecting damage elements more reliably.

Substructuring Method.
e limited number of sensors and unclearness of the damaged locations make a more detailed diagnosis method necessary.One precise diagnosis method may be the substructuring method, whereby an entire structure is divided into several substructures, as illustrated in Figure 5, and the damage is detected by measurements on each substructure.e damage detection in the substructure is carried out by means of a similar process to those of the previous complete and incomplete methods, depending on the number of measurement locations.e measured or estimated displacement data are transformed into the stress of the truss structure.
e damage yields the variation in the element stress, and it may be utilized as a damage index for evaluating the structural performance.Assuming a constant mass before and after the damage, the stress variation rate in equation ( 5) can be applied for detecting the damage.e stress variation is related to the cross-sectional variation under a constant axial force.e axial stiffness k is also related to the variation in the cross-sectional area with the same elastic modulus and member length.us, equation ( 5) for the axial members can be expressed by where k d and k u represent the stiffness values in the damaged and undamaged states, respectively, while A d and A u are the cross-sectional areas in the damaged and undamaged states, respectively.It is observed that the stiffness variation of equation ( 6) has the same meaning as the stress variation of equation ( 5).

Dynamic Approach.
In the dynamic approach, mode shape data are utilized as measurement data as with the static displacements.e measurement data may be a full set of displacement DOFs or the data expanded from incomplete measurements.e stress variation rates in the truss elements are calculated by equation ( 5), or the stiffness variation rate is predicted by equation ( 6). e damage is detected by the stiffness variation rate plots.

Example 2.
e above methods are examined in the numerical example in order to detect the damaged elements of the plane truss structure displayed in Figure 6. e numbering of nodes and members is indicated in the figure.Each node has two DOFs of the horizontal and vertical responses u and v, respectively.e truss structure of 12 bays is simply supported and is composed of 24 nodes and 45 members.All members have the same elastic modulus of 200 GPa, cross-sectional area of 2.5 ×10 −3 m 2 , and density of 7,860 kg/m 3 , as in the previous example.e total length of the truss is 6 m, its height is 0.6 m, and each bay is 0.5 m long.e truss is subjected to an external force of 10 N in the downward direction of node 7. It has multiple damages of 20% section loss at elements 4, 14, and 42.
e numerically simulated responses and stiffness values in the intact state are saved as baseline data for comparison with the numerical results in the damaged state.e above three static methods are compared, and one dynamic method is presented in this example.Moreover, each method evaluates the effect of the noise contained in the measured data.
Degrees of freedom Displacement variation between two states (mm)  Advances in Civil Engineering 2.2.1.Complete Measurement Data. Figure 7 illustrates the numerical results calculated from the complete noise-free measurement data before and after the damage.Figure 7(a) represents the horizontal and vertical displacement variations before and after the damage.e plots display the abrupt changes at the node to move from the lower to the upper chords.e damage positions at the element represent the abrupt displacement variation owing to the stiffness deterioration.However, the continuous displacement plots at the lower and upper chords do not provide any damage information without any additional considerations, such as flexural curvature.Figure 7(b) represents the stiffness variation rates calculated by equation ( 6), respectively.It is observed that abrupt stiffness variations are located at the damaged elements, and the variation rate coincides with the damage degree of 20%.e application indicates that the damage can be clearly detected if accurate information is provided on the completely measured displacements.Advances in Civil Engineering e presence of external noise leads to a slight variation in the displacements, deviating from the accurate trajectory.
e measurement data contain errors contaminated by external noise.e i-th simulated or measurement dataset ψ i is established as where c denotes the relative magnitude of the error, ξ i is a random number variant in the range of −1, 1  , and ψ 0,i is the i-th noise-free dataset.e applicability of the proposed method is investigated in the numerical example.
Figure 8(a) represents the displacement curves containing 1% noise.In order to reduce the effect owing to external noise, the mean values of 10 repeated numerical results were utilized as the measurement data for evaluating the validity of the method.It is observed that the plots in Figure 8(a) represent the displacement variations, do not indicate any peculiarity related to the damage, and exhibit difficulties in detecting the damage.Figure 8(b) displays the stiffness variation rate plot.
e noise leads to irregular stiffness variations, unlike in the noise-free case.e damaged elements are clearly detected by the plots.e stiffness variation rates do not coincide with the actual damage degree owing to the noise.is example indicates that accurate damage detection is possible when accurate measurement data despite external noise are provided.It is demonstrated that the stiffness variation rate plots are almost the same except for the irregularities owing to noise.

Incomplete Measurement Data.
At this time, we consider damage detection using measurement data with fewer DOFs than the system order.e vertical displacements at 22 nodes, except for support nodes 1 and 13, were measured, and the mean values of 10 repeated results were recorded.By inserting the numerically simulated measurement data and equilibrium equation in the intact state into equation ( 3), the displacements for the entire DOFs are estimated.Figure 9(a) illustrates the noise-free displacement differences before and after the damage.It is observed that the curve is nearly similar to that of Figure 7(a), using the complete measurement data.e damage cannot be detected from the plots.
e estimated displacements are used to predict the internal stresses of the truss members.Figures 9(b)-9(f ) illustrate the stiffness variation rates in the lower and upper chord members, diagonal members, and vertical members before and after the damage.Overall, the plots exhibit very small variation rates, except for the lower chord members and vertical members, as indicated in Figure 9. is is the reason that the expanded method rarely describes accurate responses, and the lower chord members are affected by the external load.e slash lines in Figure 10 represent the elements with high stiffness variation rates by the members.It is indicated that the damages are located at the left half of the midspan rather than at the right half, representing higher stiffness variation rates.Moreover, higher internal stress rates are represented at the left half of the midspan than at the right half.
e highest stiffness variation rates are concentrated on the left-end support and loading point.
e abrupt stiffness change in vertical member 42 is also illustrated.It is expected that the damaged members will be positioned at the left half region of the midspan and at vertical member 42.us, the plots are merely used for estimation and do not clearly indicate the damaged elements.It can be concluded that this approach cannot detect the explicitly damaged elements.
is example also considers the noise effects included in the incomplete measurement data.e measurement data     Advances in Civil Engineering are contaminated by 1% noise at the same DOFs as before and are expanded to the full set of DOFs.e numerical results are illustrated in Figure 11.e plots are very similar to those in Figure 9, and this method is not sensitive to external noise.Damage detection should be conducted by focusing on the region representing high stiffness variation rates.us, it is desirable to perform the analysis on the damage-expected element group with several measurement sensors, which constitutes the substructuring method.

Substructuring Method.
We consider the damage detection of a substructure to be isolated from the truss structure illustrated in Figure 12.
e substructure in Figure 12 is composed of 13 elements and eight nodes.Multiple damages are located at elements 5, 17, and 29, with a 20% section loss.Measurements were taken at nine of 16 DOFs, as indicated in the figure.e number of sensors was significantly reduced, from 22 sensors in the previous test to nine in this test.e sensor locations are indicated in Figure 12. e mean values of the 10 repeated numerical experiments were taken for reducing the noise effect contained in the measurement data.
e number of measurement locations is fewer than the substructure order, and data expansion is carried out.
Figure 13 displays the stiffness variation rate plots corresponding to the noise-free and 5% noise cases.Both   Advances in Civil Engineering cases demonstrate that the damage at elements 5 and 17 of the lower and upper chord members can be explicitly detected.Moreover, the diagonal and vertical elements 28, 29, and 30 display the same damage rate of 20%, and it is believed that the damage is located at one of these.e damage of the diagonal and vertical members can be estimated by the installation and measurement of additional sensors.It is observed that the method is not sensitive to noise from the plots.e large stiffness variation rate at elements 33 and 34 is owing to the external load because these are adjacent to the loading point.

Dynamic Substructuring Approach.
e dynamic approach begins with the measured mode shape data, as with the displacements in the static approach.
e stiffness variation rate is calculated by the mode shape data. is work uses the complete or incomplete mode shape data corresponding to the first natural frequency only.e incomplete mode shape data are expanded by equation (3) with the stiffness matrix in the static approach.e damage detection on the truss structure in Figure 12 is performed by means of the substructuring method.e sensor locations are the same as in the previous example.In the noise-free case illustrated in Figure 14(a), the damage at the lower and upper chord members of 5 and 17 can be explicitly detected by the abrupt stiffness variation rate.e stiffness variation rates do not coincide with the damage rate because the analysis considers the first mode only and neglects the others.e damage at diagonal member 29 is not definitely detected because the adjacent members 28, 29, and 30 exhibit similar stiffness variations.e inadequacy of the numerical results arises from the neglect of the other modes.Similar results are indicated in the 5% noise case illustrated in Figure 14(b).It is observed that the damage at the vertical or diagonal members cannot be explicitly detected, but the damage-expected region can be predicted.e plots exhibit large variation rates at the members adjacent to the loading point.e variation rates near the loading point should be carefully investigated.
According to the static-based and dynamic-based examples, it is demonstrated that the stiffness variation rate calculated by the complete and accurate displacement data indicates the accurate damage locations, despite external noise.In the case that the number of measurement sensors is fewer than the system order, the expanded dataset provides information for predicting the damage region, including the damage-expected elements.e static-based and dynamicbased damage identification approaches by means of substructuring can detect the damaged elements more clearly, despite few measurement sensors and external noise.However, the damage of the diagonal and vertical members can be predicted by the member group.
e damage-expected region should be investigated more thoroughly by means of the installation and measurement of additional sensors.

Conclusions
is work has investigated damage detection methods depending on the number of measurement locations to be causable in the evaluation of structural performance.is study has considered three different damage detection approaches utilizing complete and incomplete measurement data in an entire structure, and incomplete measurement data in substructures were compared.Furthermore, the  Advances in Civil Engineering noise effect contained in the measurement data was investigated.e measured and expanded displacements were transformed into the internal stress of the elements, and the stiffness variation rate was used as a damage index for evaluating the damage.e complete measurement data method exhibits the drawback of requiring the same number of sensors as the system order.e incomplete measurement data method exhibits the disadvantage of unstable data expansion depending on the number and installation locations of sensors.e substructuring method offers merits in reducing the system order and number of measurement sensors and is insensitive to external noise.However, the substructuring method needs to investigate the damageexpected region more closely for accurate damage detection by means of the installation and measurement of additional sensors.

Figure 1 :
Figure 1: Flow chart of this work to detect damaged elements.

Figure 3 :
Figure 3: Difference between actual and estimated responses of the truss model.

Figure 6 :
Figure 6: Numerical model of the truss structure for damage detection.

Figure 4 :Figure 5 :
Figure 4: Stress variation rate of the truss model: (a) complete measurement data; (b) incomplete measurement data.

Figure 7 :
Figure 7: Numerical results using the complete noise-free measurement data: (a) displacement variations; (b) stiffness variation rate.In (a), the solid line indicates the horizontal displacement and the dashed line indicates the vertical displacement.

Figure 8 :
Figure 8: Numerical results using the complete measurement data contaminated by 1% noise level: (a) displacement variations; (b) stiffness variation rate.In (a), the solid line indicates the horizontal displacement and the dashed line indicates the vertical displacement.

Figure 9 :
Figure 9: Damage detection using the expanded noise-free measurement data: (a) displacement variations; (b) stiffness variation rate at lower chords; (c) stiffness variation rate at upper chords; (d) stiffness variation rate at "/"-type diagonal elements; (e) stiffness variation rate at "\"-type diagonal elements; (f ) stiffness variation rate at vertical elements.In (a), the solid line indicates the horizontal displacement and the dashed line indicates the vertical displacement.

Figure 10 :Figure 11
Figure 10: Truss members to represent high stress variation rates.

Figure 11 :
Figure 11: Damage detection using the expanded data contaminated by 1% noise: (a) displacement variations; (b) stiffness variation rate at lower chords; (c) stiffness variation rate at upper chords; (d) stiffness variation rate at diagonal elements; (e) stiffness variation rate at diagonal elements; (f ) stiffness variation rate at vertical elements.In (a), the solid line indicates the horizontal displacement and the dashed line indicates the vertical displacement.

Figure 12 :
Figure 12: A substructure isolated from the entire truss structure.

Figure 13 :
Figure 13: Damage detection by the substructuring method: (a) noise-free measurement data; (b) measurement data contaminated by 5% noise.