Framework for Flexural Rigidity Estimation in Euler-Bernoulli Beams Using Deformation Influence Lines

Structural degradation is an inevitable part of a structure’s service life. Detecting structural impairments and assessing their nature is a significant challenge. Degradations reduce structural system stiffness and subsequently affect system deformations. An appropriate structural health index that is able to capture these changes in deformation and relate them to a structural system stiffness may help engineers to adequately rate structural condition. This paper outlines a theoretical framework for the utilization of static deformation influence lines for estimating the flexural rigidity of Euler-Bernoulli beams. In the proposed technique, the relationship between the second derivative of the deformation influence line and the flexural rigidity for both statically determinate and indeterminate beam structures is presented. The proposed method provides a flexural rigidity estimate (FRE) over the entire span that is based on a single measurement location and estimates both the location and severity of impairments, regardless of the location of the measurement or the damaged zones. Noisy analytical simulations are presented with noise levels of 0%, 0.5%, 1%, 2%, 3%, and 4%; in all cases the modeled damage is quantified and localized using the FRE. A laboratory experiment is also presented that validates the theoretical framework.


Introduction
Detecting, classifying, and mitigating structural impairments are significant challenges for structural engineers, inspectors, and owners.The field of structural health monitoring (SHM) has developed to address these challenges associated with aging infrastructure [1][2][3][4][5][6][7][8][9].Rytter outlines criteria for classification of SHM algorithms into four levels based on their output information [2].At the simplest level (level-1), the SHM algorithm should provide a general indication of damage, while in levels 2-4 additional information regarding damage location, severity, and remaining service life of the structure will be provided, respectively.Algorithms utilizing appropriate structural health indices may increase the SHM classification level and aid engineers in adequately rating structural condition.Appropriate indices include the following essential features [10][11][12][13]: 1.
foundation in mechanical theory, 3. consistency in the evaluation technique, and 4.
efficacy and efficiency of the technique in practice.
Many damage scenarios reduce a structural system's stiffness; these stiffness degradations affect system deformations and lead to structural impairment [14,15].Significant advancements in state of the art of deformation estimations (e.g., camera-based systems) and their increasing applications present opportunities for monitoring frameworks that utilize deformation data streams [16][17][18].This paper presents an index and technique for locating and quantifying reductions in flexural rigidity along a beam structure, by analyzing deformation influence lines at a single location on the beam.

Influence Lines for Euler-Bernoulli Beam Evaluation
An influence line function represents the magnitude of a desired response (e.g., strain, shear, moment, deflection, rotation) at a fixed point on a structural system that is caused by a static load pattern placed anywhere on the structure.If the load pattern is limited to a single unit load, the influence line is called a Unit Influence Line (UIL).The UIL technique has been implemented widely in bridge engineering design and load rating.
A system's static response influence lines are comprehensive evaluation indices that contain many of the previously outlined features that are rich in diagnostic content [10][11][12][19][20][21].Several studies utilizing an influence line approach for bridge evaluation and monitoring are outlined in the literature.A great portion utilized the internal responses, such as strain [11,12,19,22,23] or stress influence lines [20,24].Others use global deformation influence lines [10,21,25].Changes in flexural rigidity anywhere along a span will manifest in global deformation influence lines at any other point on the beam, regardless of the determinacy of the system.
The preliminary efforts of employing deformation influence lines as a bridge condition index gathered measurements from diagnostic crawl-speed tests and derived the unit load influence line using numerical decomposition [10,21,25].Turer utilized two methods to show that a conceptual relation between flexibility and deformation influence coefficients exists [10].
Over the last decade, some studies have proposed a method to locate damage by analyzing the absolute differences between the curvature of deflection influence lines for intact and damaged beam structures [26][27][28][29][30].These studies provide influence-line-based methods to detect and localize damage that is a critical step in monitoring and evaluating a structure.
This paper presents a theoretical framework for using static global response in gravity-carrying components of beam structures for impairment detection.Utilizing the second derivative of the static deformation influence line as an impairment index, the proposed technique presents a robust method of characterizing both the location and the severity of impairments in beam structures based on deformation data streams.

Derivation for FRE
A straight beam with linear elastic material and two-dimensional linear geometric behavior is illustrated in Figure 1a.A unit force is applied on this beam at distance ξ from the left support location.The vertical deflection, y(x, ξ), (dashed lines in Figure 1a) is a function of the location of the unit load, ξ, and varies with the position along the beam.The first argument in the function y(x, ξ) represents the monitoring location along the beam.The second argument is the location of the unit load.If the first argument varies and the second is fixed, then y is a deflection function when the unit load is fixed; if the opposite is true, then y is a DIL of the fixed first argument.
The vertical deflection influence line (DIL) of an arbitrary point A located at distance x A from the left support, y(x A , ξ), is defined as the deflection of the beam at point A as a unit load located at ξ traverses the span, i.e., (0 < ξ < L).
According to Maxwell's theorem of reciprocal displacement, the deflection at point A due to a unit load acting at distance ξ (i.e., y(x A , ξ)) is equal to the beam deflection at distance ξ when the unit load is acting at point A (i.e., y(ξ, x A )). So, using Maxwell's theorem gives: If ξ varies from 0 to L, the left side of Equation ( 1) is the DIL of point A, and the right side of Equation ( 1) represents the beam's deflected shape when a unit load is applied at point A. This is illustrated in Figure 1b.

Derivation for FRE
A straight beam with linear elastic material and two-dimensional linear geometric behavior is illustrated in Figure 1a.A unit force is applied on this beam at distance from the left support location.The vertical deflection, ( , ), (dashed lines in Figure 1a) is a function of the location of the unit load, , and varies with the position along the beam.The first argument in the function ( , ) represents the monitoring location along the beam.The second argument is the location of the unit load.If the first argument varies and the second is fixed, then is a deflection function when the unit load is fixed; if the opposite is true, then is a DIL of the fixed first argument.
(a) (b)  Assuming that both sides of Equation ( 1) are twice differentiable, differentiating both sides of Equation ( 1) twice with respect to ξ yields: The left side of Equation ( 1) may be acquired practically by collecting deflection data at point A as a unit load traverses the structure and differentiating twice.
The right side of Equation ( 1) represents the deflected shape of the beam when a unit load is located at x A ; therefore, the right side of Equation ( 2) is the beam curvature when a unit load is applied at distance x A .If Euler-Bernoulli behavior is considered, then the right side of Equation ( 2) can be related to internal moment and flexural rigidity: In Equation (3), EI(ξ) is the current beam flexural rigidity at ξ and may be a variable function along the beam.Moreover, the moment m(ξ) is the internal moment when a unit load is acting at point A.
By using Maxwell's reciprocal theorem and the same procedure, it can be shown that Equation (3) is also valid for the rotation influence line (RIL).In the case of using the RIL instead of the DIL (i.e., θ(x A , ξ) instead of y(x A , ξ)), the moment m(ξ) in Equation (3) equals the internal moment when a unit moment is acting at the point A.
For simplification and generalization, the term y(x A , ξ) (or θ(x A , ξ)) is replaced with u A (ξ). Equation (3) becomes: Practically, EI(ξ) will never be equal to zero in Equation (4).So, in regions that m(ξ) is zero, u A (ξ) becomes zero and in regions that m(ξ) is not equal to zero, Equation (4) yields: Equation ( 5) is the Flexural Rigidity Estimate (FRE), which provides EI(ξ), the beam's flexural rigidity at location ξ as a function of a known internal moment m(ξ) and the second derivative of a measured DIL or RIL, u A (ξ).
By comparing the nominal flexural rigidity (EI 0 ) with the FRE calculated from Equation ( 5), the location and quantity of the impairment to the flexural rigidity can be estimated.

Alternate Derivation for FRE
The deflection of point A in Figure 1 as a function of a traversing unit load, u A (ξ), can be calculated using the flexural equation for virtual work.
In the case of linear geometry, elastic material behavior, and negligible shear deformations, the small flexural deflection u A (ξ), under the effects of a unit load at distance ξ is expressed as: In Equation (6), the moment M(x, ξ) is the beam internal moment caused by an externally applied load (i.e., the moving unit load) and m(x) is the internal moment caused by a virtual unit load applied at point A. The moment m(x) is a function of location x A , but independent of variable ξ.
The moment M(x, ξ) in Equation ( 6) is the only parameter on the right side that is a function of ξ.Both m(x) and EI(x) remain constant as a unit load traverses the beam.Hence, the first and second derivatives of u A (ξ) with respect to ξ are equal to Equations ( 7) and (8), respectively.
The moment M(x, ξ) in Equation ( 6) is the internal moment influence line at point x.According to the Muller-Breslau principle, the influence line for this internal moment is the scaled deflected shape of a beam with released moment resistance that is subject to a virtual displacement at a frictionless hinge.This modified beam with an internal hinge is depicted in Figure 2. In other words, the moment M(x, ξ) is equal to: Infrastructures 2017, 2, 23 4 of 19 In Equation ( 6), the moment ( , ) is the beam internal moment caused by an externally applied load (i.e., the moving unit load) and ( ) is the internal moment caused by a virtual unit load applied at point .The moment ( ) is a function of location , but independent of variable .
The moment ( , ) in Equation ( 6) is the only parameter on the right side that is a function of .Both ( ) and ( ) remain constant as a unit load traverses the beam.Hence, the first and second derivatives of ( ) with respect to are equal to Equations ( 7) and (8), respectively.
= ( ) = ( , ) • ( ) The moment ( , ) in Equation ( 6) is the internal moment influence line at point .According to the Muller-Breslau principle, the influence line for this internal moment is the scaled deflected shape of a beam with released moment resistance that is subject to a virtual displacement at a frictionless hinge.This modified beam with an internal hinge is depicted in Figure 2. In other words, the moment ( , ) is equal to: In Equation ( 9), . .( ) and α are the deflection curve and angular displacement at of the modified beam under the effects of a virtual displacement, as illustrated in Figure 2. The second derivative of moment with respect to is equal to: The right side of Equation ( 10) is related to the curvature of the modified beam that is subject to a small virtual displacement, as shown in Figure 2.For the case of a statically determinate beam, the moment influence line is linear, so the curvature of the deflected shape of the modified beam must be zero except for at the location of the hinge, where the curvature is infinite.Hence, the right side of Equation ( 10) is equal to: The right side of Equation ( 11) is the Dirac delta function centered at = .Substituting Equations ( 10) and ( 11) into Equation ( 8) gives: In Equation ( 9), y M.B. (ξ) and α x are the deflection curve and angular displacement at x of the modified beam under the effects of a virtual displacement, as illustrated in Figure 2. The second derivative of moment M. with respect to ξ. is equal to: The right side of Equation ( 10) is related to the curvature of the modified beam that is subject to a small virtual displacement, as shown in Figure 2.For the case of a statically determinate beam, the moment influence line is linear, so the curvature of the deflected shape of the modified beam must be zero except for at the location of the hinge, where the curvature is infinite.Hence, the right side of Equation ( 10) is equal to: The right side of Equation ( 11) is the Dirac delta function centered at ξ = x.Substituting Equations ( 10) and (11) into Equation ( 8) gives: According to the properties of the Dirac delta function, when 0 < ξ < L, Equation ( 12) yields: This is the same as Equation ( 4).Superposition is used when applying the framework to statically indeterminate structures as shown in the following section.

Calculation of Moment m(ξ)
Equation ( 5) can be used to detect the status of a beam's flexural rigidity.In practice, the influence line is extracted by solving an inverse problem.First, the desired response (e.g., vertical deflection or rotation) of a specific location is recorded.Then, by synchronizing the recorded measurement with the location of moving loads (e.g., vehicular loads) and solving an inverse problem, the influence line, u A (ξ), of a unit load can be extracted [10-12,19,21-23].
The second derivative of u A (ξ) can be estimated numerically.The moment m(ξ) in Equation ( 5) must be calculated or estimated.The moment m(ξ) will be independent of a system's flexural stiffness if the system is statically determinate; conversely, the moment m(ξ) will be dependent upon a system's flexural stiffness if the system is statically indeterminate.The methods for calculating and estimating the moment m(ξ) are as follows.

Statically Determinate Systems
In a statically determinate system, m(ξ) is not a function of beam flexural stiffness, but is a function of the location (x A ) and the type of recorded deformation.Hence, no additional information is required to calculate the m(ξ) diagram.
In the case of using a DIL: In the case of using a RIL:

Statically Indeterminate Systems
In a statically indeterminate system (e.g., continuous multi-span beam), the moment m(ξ) is a function of the system's flexural stiffness.Hence, calculation of the moment m(ξ) requires additional information.Two techniques are described to calculate and estimate the moment m(ξ) for both DIL and RIL.

Deflection Influence Line (DIL)
In the case of using a DIL, the moment m(ξ) will be the diagram of the moment that is created by a unit load that is applied at distance x A .In a linear system, the response of a statically indeterminate or continuous beam with a degree of indeterminacy (DOI, i.e., the number of interior vertical support reactions) can be calculated by the superposition of DOI + 1 number of statically determinate systems (i.e., simply supported beams).For this purpose, each vertical support (support i) is replaced by a vertical reaction R i .The moment that m(ξ) of the continuous beam is calculated by the superposition of moment of these DOI + 1 statically determinate simply supported beams: In Equation ( 16), m UL (ξ) is the moment that is caused by a unit load applied at a distance x 0 on a simply supported beam.Moment m i (ξ) is the moment that is caused by a unit load that is applied on the simply supported beam at location of the interior support i.If x i serves as the location of interior support, i, the equations for m UL and m i are as follows: Moreover, in Equation ( 16), reaction R i is the reaction of interior support i when a unit load is applied on a continuous beam at distance ξ = x A .These reactions are unknown and additional information is required.
To get this additional information, assume that a strain gage is attached to an adequate cross-section location at each support location on the continuous beam.By having knowledge about the status of the flexural rigidity on interior support locations, the internal moment at each support location can be calculated by converting the recorded strains.
Thus, when the moving unit load is located at a distance ξ = x A , the moment at each support can be measured by using these strain gages.These measured moments constitute the moment m(ξ) at each support location (i.e., m(x 1 ), m(x 2 ), . . . ,m(x DOI )).So, Or, in matrix format: As can be seen in Equation ( 20), the measured moments give DOI number of equations with DOI number of unknowns (R i ).By solving Equation (20), the unknowns R i can be solved as: Therefore, when a DIL is measured, the general diagram of moment m(ξ) of the continuous beam can constructed by substituting the resultant R i from Equation (21) in Equation ( 16).

Rotation Influence Line (RIL)
In the case of measuring the RIL, the moment m(ξ) is the moment that is caused by a concentrated unit moment at measurement location (x A ).The influence line is constructed using a unit point load; the technique proposed for using a deflection influence line is not applicable here, and the measured moments (i.e., m(x 1 ), m(x 2 ), . . . ,m(x DOI )) must be preprocessed.
The first derivative of the deflection influence line can be explained using Equation (22): The right side of Equation ( 22) can be expanded to: According to the definition of the deflection influence line, u(ξ) is the beam deformation at a location x 0 when a unit load is applied at a distance ξ.So, in Equation ( 22), the terms u(ξ+∆ξ) ∆ξ and u(ξ) ∆ξ are the beam deformations at location x 0 when a moving load with intensity equal to 1 ∆ξ is applied at distances ξ + ∆ξ and ξ, respectively.By using the principle of superposition, u in Equation ( 22) is the beam deformation under the effects of two concentrated loads that are applied in opposite directions.These concentrated loads can be replaced by a couple with intensity equal to: As ∆ξ approaches zero, the location of this couple approaches ξ.This means that u (ξ) is the beam deformation at location x 0 that is caused by a positive unit moment applied at location ξ.The first derivative of the deformation influence line under the effects of a moving unit load can be interpreted as the deformation influence line under the effects of a moving unit moment.
In a similar manner, the first derivative of the recorded moment measurement at each interior support location can be assumed as interior support moment under the effects of moving unit moment.By taking the first derivative of the recorded moment at each interior support location, the moment m(ξ) at each support location (i.e., m(x 1 ), m(x 2 ), . . . ,m(x DOI )) is determined.In this section, m i (ξ) has the same definition used in Equation (18).By gathering this information, Equation (4) still can be utilized to calculate the diagram of moment m(ξ).

Example I-Statically Determinate System
As Δ approaches zero, the location of this couple approaches .This means that ′( ) is the beam deformation at location that is caused by a positive unit moment applied at location .The first derivative of the deformation influence line under the effects of a moving unit load can be interpreted as the deformation influence line under the effects of a moving unit moment.
In a similar manner, the first derivative of the recorded moment measurement at each interior support location can be assumed as interior support moment under the effects of moving unit moment.By taking the first derivative of the recorded moment at each interior support location, the moment ( ) at each support location (i.e., ( ), ( ), … , ( )) is determined.In this section, ( ) has the same definition used in Equation (18).By gathering this information, Equation (4) still can be utilized to calculate the diagram of moment ( ).Using the arbitrary values that are listed in Table 1, a finite element (FE) model is created to solve the defined problem.By using this FE model, the DIL is calculated.In this solution, the beam is divided into 1500 elements, with flexural stiffness and two degrees of freedom at each node.Using the arbitrary values that are listed in Table 1, a finite element (FE) model is created to solve the defined problem.By using this FE model, the DIL is calculated.In this solution, the beam is divided into 1500 elements, with flexural stiffness and two degrees of freedom at each node.The resultant DIL and its first and second derivatives are plotted in Figure 4a-c, respectively.In practice, the resultant DIL is a series of points.Hence, to get the first and second derivatives, numerical methods should be utilized.To get the second derivative in this example, the second order central difference method (CDM) with three points is implemented.By substituting the calculated second derivative and Equation ( 14) into Equation ( 5), the beam's flexural rigidity can be calculated.The resultant flexural rigidity is illustrated in Figure 4d.As can be seen in Figure 4, the proposed formula is robust enough to give the exact diagram of flexural rigidity.The differences between this plot and the flexural rigidity diagram obtained from nominal conditions give the location and magnitude of damages along the beam.Moreover, the piecewise equation of moment m, the resultant moment that is caused by a unit point load applied at distance x 0 on a simply supported beam, is equal to Equation (14).

Example II-Statically Indeterminate Systems
By substituting the calculated second derivative and Equation ( 14) into Equation ( 5), the beam's flexural rigidity can be calculated.The resultant flexural rigidity is illustrated in Figure 4d.
As can be seen in Figure 4, the proposed formula is robust enough to give the exact diagram of flexural rigidity.The differences between this plot and the flexural rigidity diagram obtained from nominal conditions give the location and magnitude of damages along the beam.As can be seen in Figure 4, the proposed formula is robust enough to give the exact diagram of flexural rigidity.The differences between this plot and the flexural rigidity diagram obtained from nominal conditions give the location and magnitude of damages along the beam.

Using Deflection Influence Line (DIL)
In this section, the DIL is implemented to solve the damage detection problem in Example 2. Like Example 1, the resultant DIL and its first and second derivatives are plotted in Figure 6a-c using a FE model and the CDM.The moment, m(ξ), is a function of the system's flexural rigidity in indeterminate systems, and its calculation requires additional information.To obtain this additional information, a strain gage is attached to the cross-section at the interior support location.If the status of the flexural rigidity at this cross-section is already known, the measured strains at this cross-section can be converted to internal moment at this location.
The moment is calculated from strain that occurs when the moving unit load is located at distance ξ = x A .This moment is substituted into Equation ( 21).The resultant R 1 is then substituted into Equation ( 16) and the general function of moment m(ξ) in this two-span continuous beam is constructed.Using this procedure, the diagram of moment m(ξ) is plotted in Figure 6d.
By using Equation ( 4), the resultant calculated beam flexural rigidity is presented in Figure 6e.As can be seen from Figure 6e, the proposed formula was able to predict the flexural rigidity correctly, aside from a small region between locations x = 3.3 and x = 3.7 where the flexural rigidity curve diverges.
As stated previously, the proposed formula in Equation ( 5) is acceptable when the terms u (ξ) or m(ξ) are not equal to zero.The ratio of m(ξ)/u (ξ) is unstable around zero points of m(ξ).Due to the implementation of numerical methods in the calculation of the second derivative, the magnitudes of u are not precise enough to give the correct answers near zero values.Consequently, flexural rigidity is not precisely predicted where m(ξ) approaches zero.
into Equation ( 16) and the general function of moment ( ) in this two-span continuous beam is constructed.Using this procedure, the diagram of moment ( ) is plotted in Figure 6d.
By using Equation ( 4), the resultant calculated beam flexural rigidity is presented in Figure 6e.As can be seen from Figure 6e, the proposed formula was able to predict the flexural rigidity correctly, aside from a small region between locations = 3.3 and = 3.7 where the flexural rigidity curve diverges.As stated previously, the proposed formula in Equation ( 5) is acceptable when the terms ′′( ) or ( ) are not equal to zero.The ratio of ( ) ( ) ⁄ is unstable around zero points of ( ).Due to the implementation of numerical methods in the calculation of the second derivative, the magnitudes of ′′ are not precise enough to give the correct answers near zero values.Consequently, flexural rigidity is not precisely predicted where ( ) approaches zero.

Using Rotation Influence Line (RIL)
If the beam rotation influence line at an interior support is used, then the m(ξ) diagram will have a jump at the location of the interior support.The function m(ξ) is not defined at this point, and the location of this jump is the only point at which the corresponding flexural rigidity cannot be calculated using Equation (5).However, according to the procedure defined in Section 5.2.1, it is required that the status of beam cross-section on interior support locations be investigated in order to implement the proposed method in a statically-indeterminate system.Therefore, the only cross-section that could not be monitored by using the RIL is already diagnosed.
In other words, by using the RIL, the theoretical solution for the flexural rigidity of all cross-sections can be accurately assessed.If the central difference method is used for calculating the second derivative, then using the beam RIL at an interior support location is recommended rather than the DIL.
The resultant diagrams for Example 2 when the RIL is utilized are given in Figure 7.As can be seen in Figure 7e, the RIL was able to correctly predict the flexural rigidity of all cross-sections.
In other words, by using the RIL, the theoretical solution for the flexural rigidity of all crosssections can be accurately assessed.If the central difference method is used for calculating the second derivative, then using the beam RIL at an interior support location is recommended rather than the DIL.
The resultant diagrams for Example 2 when the RIL is utilized are given in Figure 7.As can be seen in Figure 7e, the RIL was able to correctly predict the flexural rigidity of all cross-sections.

Application
In practice, measurements can be noisy, and the finite difference method may not produce a suitable second derivative for noisy data.One approach for addressing this challenge is to use a parametric approximation of the measurements.This section introduces an algorithm to overcome the obstacles of noisy measurements and its application to a series of numerical and experimental results.

Application
In practice, measurements can be noisy, and the finite difference method may not produce a uitable second derivative for noisy data.One approach for addressing this challenge is to use a parametric approximation of the measurements.This section introduces an algorithm to overcome the obstacles of noisy measurements and its application to a series of numerical and experimental results.

Algorithm for Addressing Noisy Measurements
The approach seeks to fit a smooth parametric function to the measured noisy influence line.After curve fitting, an estimation of the second derivative of the influence line is constructed using this approximated function.The fitted function should be continuous and smooth while closely estimating the numerical values at each measurement.
A piecewise quadratic spline is considered as the parametric function.This function consists of NS number of quadratic continuous functions and is defined as Equation (25).
This piecewise function should be constrained as follows: • It will become zero at the location of each support.

•
The value of the function and its first derivative would become equal at both the left-and right-hand boundaries of the considered pieces.
The problem of fitting a smooth curve can be defined as Equation (26).
In Equation (26), the first term fits the curve to the measurements, while the second term, called the regularization term, smooths the resultant fitted curve.The value of the parameter λ makes a balance, which should be optimized [31].
This optimization problem can be simplified as a constrained linear least-squares problem, which may be solved by the robust lsqlin function implemented in the MATLAB programming software.The constraints and the objective function are defined as linear equations in the format A cts X = b cts and A obj X = b obj , respectively, where vector X is the vector of unknowns.Appendix A presents the method for calculating these matrices.Resultant matrices A obj , b obj , A cts , and b cts are input into the lsqlin function using the "Interior-Point" algorithm [32].
The regularization parameter λ should be optimized to balance fitting the measurement points with obtaining a smooth curve.The L-curve method is utilized for this optimization problem; the log-log plot of the residual of each term of the objective function for different values of parameter λ is plotted.The resultant plot will have an L-shape.The value corresponding to the corner of this curve is the optimal point.This optimal value of λ makes a balance between a smooth curve and a well-fitted curve [31].
Even after fitting a smooth curve to noisy data, some estimation error will remain in the calculation of the second derivative.The FRE is more sensitive to the second derivative calculation error in regions with small values from the unit moment curve; this occurs because the second derivative approaches small values in the denominator of Equation (5).To overcome this challenge, multiple measurements of influence line were recorded, and a weighted average of the results was used.

Application on a Numerical Model
To demonstrate the robustness of the proposed method, the approach is applied on a series of analytical models.Analytical models closely simulate the experimental test that was conducted in Section 6.3.For this purpose, a simply supported beam with a geometrical imperfection near the middle of the span is considered, as shown in Figure 8.Any normalized flexural rigidity value presented in this paper is achieved by using the constant parameter EI 0 = 1, 268, 322.8 (N-mm 2 ), representing the intact flexural rigidity defined in Figure 8.The RIL at both ends is calculated using finite element method.Equation ( 27) is used to create noisy RIL data: In Equation ( 27), the parameter is the noise level, rms(U ) is the root-mean-squared The RIL at both ends is calculated using finite element method.Equation ( 27) is used to create noisy RIL data: U noisy = U noise−free + ε•rms(U noise−free )•r (27) In Equation (27), the parameter ε is the noise level, rms(U noise−free ) is the root-mean-squared of noise-free data, and vector r is a normal random vector.Vector r has a unit standard deviation, a mean equal to zero, and the same vector length as the U noise−free vector.In this paper, calculated RILs are polluted with noise levels ε equal to 0.5%, 1%, 2%, 3%, and 4%.The 4% noise level is significantly higher than the noise level that is observed in the experiment presented in Section 6.3 and illustrates that the FRE is robust to this level.
Utilizing the proposed parametric approximation, the piecewise function is fitted to the artificial noisy measurements and the second derivative of the fitted curve is calculated.The resultant second derivatives are input into a weighted average equation: where Equation ( 28) gives the weighted average of the calculated flexural rigidity from each RIL.Weights are selected in a way that calculated values for flexural rigidity are mostly affected by the resultant values from the closer RIL.
Figure 9 presents the results of the fitted curve of the support RILs at both ends (θ A of support A is closer to damaged region, in this example).The results of FRE are shown for different noise levels.Even in higher noise levels, the method gives acceptable results regarding flexural rigidity estimation.

Application on a Real World System
A rectangular steel specimen with 1.5 mm nominal thickness and 25.4 mm nominal width is supported on the edge of an L-shape bracket in a simply supported manner, as illustrated in the Figure 10.A geometrical damage is created on this specimen by reducing the cross-section at the specified location in Figure 8.The geometry of this simply supported beam and considered nominal values for this experiment is illustrated in Figure 8.
The unit RIL was created by measuring rotations at both supports, caused by passing a hanging load of 0.5 kg weight over the beam span with the resolution of 6 mm distance (i.e., totally 81 measurements).Measurements were conducted in a static manner.
The rotation measurements are conducted by using a machine vision algorithm that is implemented in MATLAB and a robust camera calibration method that tracks the checkerboards that are attached to the ends of this beam.A camera with a resolution of 2592 × 2048 pixels, a 75 mm lens, and 15 × 16 checkerboards with 4.8 × 4.8 mm squares are utilized for this experiment setup.The RILs at the ends of this simply supported beam are calculated by tracking the movements of the checkerboards corners while the load is moved along the beam span.Zeinali et al. [33] present the details of this robust measurement method.One advantage of a camera-based displacement estimation is that the location of the load is simultaneously estimated.
A quadratic piecewise function with 82 pieces and λ opt = 0.0018 is utilized for the curve fitting.Figure 11a shows the results of RIL measurements and curve fitting to these measurements.By taking the second derivative of the resultant fitted curves and substituting into Equation (28), the flexural rigidity of the specimen is estimated, and results are plotted in Figure 11b and labeled "Estimated EI".

Application on a Real World System
A rectangular steel specimen with 1.5 mm nominal thickness and 25.4 mm nominal width is supported on the edge of an L-shape bracket in a simply supported manner, as illustrated in the Figure 10.A geometrical damage is created on this specimen by reducing the cross-section at the specified location in Figure 8.The geometry of this simply supported beam and considered nominal values for this experiment is illustrated in Figure 8.
The unit RIL was created by measuring rotations at both supports, caused by passing a hanging load of 0.5 kg weight over the beam span with the resolution of 6 mm distance (i.e., totally 81 measurements).Measurements were conducted in a static manner.
The rotation measurements are conducted by using a machine vision algorithm that is implemented in MATLAB and a robust camera calibration method that tracks the checkerboards that are attached to the ends of this beam.A camera with a resolution of 2592 × 2048 pixels, a 75 mm lens, and 15 × 16 checkerboards with 4.8 × 4.8 mm squares are utilized for this experiment setup.The RILs at the ends of this simply supported beam are calculated by tracking the movements of the checkerboards corners while the load is moved along the beam span.Zeinali et al. [33] present the details of this robust measurement method.One advantage of a camera-based displacement estimation is that the location of the load is simultaneously estimated.A quadratic piecewise function with 82 pieces and = 0.0018 is utilized for the curve fitting.
Figure 11a shows the results of RIL measurements and curve fitting to these measurements.By taking the second derivative of the resultant fitted curves and substituting into Equation ( 28), the flexural rigidity of the specimen is estimated, and results are plotted in Figure 11b and labeled "Estimated EI".A quadratic piecewise function with 82 pieces and = 0.0018 is utilized for the curve fitting.
Figure 11a shows the results of RIL measurements and curve fitting to these measurements.By taking the second derivative of the resultant fitted curves and substituting into Equation ( 28), the flexural rigidity of the specimen is estimated, and results are plotted in Figure 11b and labeled "Estimated EI".The geometry of this variable rectangular cross-section was measured using calipers with 1/100 mm accuracy at 5 mm increments.Assuming E = 200 GPa and using the measured geometry, the flexural rigidity of this beam was calculated and plotted and labeled as "Actual EI" in Figure 11b.The nominal and actual flexural rigidities shown in Figure 11b differ due to small deviations in the specimen's cross-section.Actual EI values are the reference values to which the estimated EI values from FRE results are compared.
Even in the presence of measurement noise and geometric deviations in intact sections of the specimen, the FRE method can quantify and localize the damage as a 40% reduction from 0.25 L to 0.34 L. Variations in the actual EI result from imperfections in the specimen; small variations in specimen thickness along the length can result in large departures from the nominal EI values.This result suggests that the FRE method can detect multiple reductions in EI that is caused by multiple damage locations.The geometry of this variable rectangular cross-section was measured using calipers with 1/100 mm accuracy at 5 mm increments.Assuming E = 200 GPa and using the measured geometry, the flexural rigidity of this beam was calculated and plotted and labeled as "Actual EI" in Figure 11b.The nominal and actual flexural rigidities shown in Figure 11b differ due to small deviations in the specimen's cross-section.Actual EI values are the reference values to which the estimated EI values from FRE results are compared.

Summary
Even in the presence of measurement noise and geometric deviations in intact sections of the specimen, the FRE method can quantify and localize the damage as a 40% reduction from 0.25 L to 0.34 L. Variations in the actual EI result from imperfections in the specimen; small variations

Figure 1 .
Figure 1.Initial position and deflected shape of a beam structure under the effects of a unit load located at (a) ; (b) .

Figure 1 .
Figure 1.Initial position and deflected shape of a beam structure under the effects of a unit load located at (a) ξ; (b) x A .

Figure 2 .
Figure 2. Deflected shape of the modified beam with a replaced hinge at measurement point under the effects of a virtual displacement.

Figure 2 .
Figure 2. Deflected shape of the modified beam with a replaced hinge at measurement point under the effects of a virtual displacement.

Figure 3
Figure3depicts a simply supported beam with a reduction in flexural rigidity from EI 0 (nominal rigidity) to βEI 0 (0 < β < 1) between locations L 1 and L 2 .The deflection influence line of the beam at location x 0 is available from recorded deflection data.

Figure 3
Figure 3 depicts a simply supported beam with a reduction in flexural rigidity from (nominal rigidity) to (0 < < 1) between locations and .The deflection influence line of the beam at location is available from recorded deflection data.

Figure 5
Figure 5 illustrates a statically indeterminate system similar to the one in Example I but with an additional interior support located at x 1 = 4.

Figure 5
Figure 5 illustrates a statically indeterminate system similar to the one in Example I but with an additional interior support located at = 4.

Figure 5 .
Figure 5. Considered simply supported beam in Example-II.Figure 5. Considered simply supported beam in Example-II.

Figure 5 .
Figure 5. Considered simply supported beam in Example-II.Figure 5. Considered simply supported beam in Example-II.

Figure 8 .
Figure 8. Geometry configuration and material of the numerical model.

Figure 8 .
Figure 8. Geometry configuration and material of the numerical model.

Figure 9 .
Figure 9. Results of application of the proposed method on a simply supported beam with polluted RILs.

Figure 9 .
Figure 9. Results of application of the proposed method on a simply supported beam with polluted RILs.Infrastructures 2017, 2, 23 15 of 19

Table 1 .
Value of Constant Parameters in Example-1 and Example-2.

Table 1 .
Value of Constant Parameters in Example-1 and Example-2.