Structural Assessment under Uncertain Parameters via the Interval Optimization Method Using the Slime Mold Algorithm

: Damage detection of civil and mechanical structures based on measured modal parameters using model updating schemes has received increasing attention in recent years. In this study, for uncertainty-oriented damage identiﬁcation, a non-probabilistic structural damage identiﬁcation (NSDI) technique based on an optimization algorithm and interval mathematics is proposed. In order to take into account the uncertainty quantiﬁcation, the elastic modulus is described as unknown-but-bounded interval values and the proposed new scheme determines the upper and lower bounds of the damage index. In this method, the interval bounds can provide supports for structural health diagnosis under uncertain conditions by considering the uncertainties in the variables of optimization algorithm. The model updating scheme is subsequently used to predict the interval-bound of the Elemental Stiffness Parameter (ESP). The slime mold algorithm (SMA) is used as the main algorithm for model updating. In addition, in this study, an enhanced variant of SMA (ESMA) is developed, which removes unchanged variables after a deﬁned number of iterations. The method is implemented on three well-known numerical examples in the domain of structural health monitoring under single damage and multi-damage scenarios with different degrees of uncertainty. The results show that the proposed NSDI methodology has reduced computation time, by at least 30%, in comparison with the probabilistic methods. Furthermore, ESMA has the capability to detect damaged elements with higher certainty and lower computation cost in comparison with the original SMA.


Introduction
Structural systems in civil and mechanical engineering are subjected to damage and deterioration during their service life. Damage is defined as a weakening of a structure that may result in undesired displacements, stresses, strains, or vibrations, resulting in unexpected and catastrophic consequences. Therefore, early detection of damage can improve safety and extend serviceability of infrastructures [1].
Modal parameters of a structure are usually used to detect structural damages [2][3][4][5][6][7][8][9][10], since modal parameters (modal frequencies and mode shapes) are functions of the physical parameters (mass and stiffness). As a result, the presence of damage causes changes in the structure's modal characteristics. Modal parameters can also be monitored conveniently and precisely. The experimentally obtained modal parameters are considered accurate and deterministic in conventional applications of model updating for damage detection [11][12][13][14][15][16][17][18]. However, there are always uncertainties in the recorded modal parameters in real applications, which can lead to unreliable and inaccurate structural damage assessments. As a result, uncertainties in structural damage detection must be considered [19,20].
In the domain of damage detection, research addressing uncertainty has been conducted utilizing probability analysis and fuzzy theory [21]. The randomness of the structural finite element (FE) model and the test data were employed by Xia et al. [22]. The probability distribution characteristics of the FE model were determined using the Monte Carlo technique and data perturbation. The structural damage probability was then calculated using the interval estimation approach. Despite extensive studies that have been performed, some problems have been identified by other scholars in the study of uncertain damage identification. One major difficulty is that the existing uncertain damage detection methods based on probability analysis and fuzzy theory typically require a large quantity of sample data to explain the source of uncertainty [23]. However, the probability model is commonly used to follow specific distributions in engineering practice because of insufficient experimental information and limited cost [24].
As a result, the accuracy and reliability of approaches based on probability and fuzzy theory cannot be guaranteed [25]. In interval analysis, uncertain parameters are represented as interval numbers that only require the uncertainty parameters' boundaries [26]. As a result of its convenience and feasibility, the interval analysis approach is extensively used in damage identification and other domains [27,28]. The interval optimization technique, which combines interval analysis with an intelligent algorithm, is a new way to handle interval problems [29]. Based on an auto-regressive model and interval optimization, Liu and Mao [30] presented a non-probabilistic damage detection approach. To determine interval bounds of uncertain variables, Ghiasi et al. [31] suggested a non-probabilistic Least Square Support Vector Machine (LS-SVM) approach. The results show that the proposed method can successfully identify the structural damage elements. However, the interval optimization method is still immature and suffers from the problem of high computational cost, especially when handling a large-scale structure.
In comparison with probabilistic methods, research on non-probabilistic uncertain damage identification has been relatively limited. Scholars have used various algorithms to raise the calculation precision of interval methods, but interval expansion is still a severe problem [32]. Therefore, in this study, a non-probabilistic structural damage identification (NSDI) approach based on a model updating scheme and interval mathematics is proposed in order to address the aforementioned shortcomings.
Damage estimation strategies using modal data are mainly based on model updating methods. These methods rely on a structure's parametric model and the optimization of the designated objective functions based on the differences between measured data and model predictions [33,34]. The accuracy of the FE model, the quality of the modal tests, the formulation of the optimization problem, and the capabilities of the optimization algorithm all play a role in the success of the finite element (FE) model updating process. In this paper, a new optimization algorithm is proposed.
Li et al. [35] recently developed SMA, a simple and effective continuous optimization technique for optimization problems. The performance of SMA is improved in this study by deleting unchanging variables after a certain number of iterations. This aspect increases the approach's adaptability for a larger range of practical applications while retaining the fundamental SMA's attractive features.
Based on the aforementioned discussion and in order to tackle the shortcomings of non-probabilistic based schemes, the objective of this paper is to investigate the influence Appl. Sci. 2022, 12, 1876 3 of 21 of uncertainty on damage identification of civil and mechanical structures. To consider the uncertainties in the measurement data, a novel non-probabilistic approach is introduced. Using this method, Young's modulus, which acts as an elemental stiffness parameter (ESP), is defined as a variable of the SMA algorithm. Moreover, damage indexes are defined by utilizing the natural frequencies and mode shapes of the structure. Furthermore, measurement noise of the structural responses is assumed to be coupled rather than statistically distributed using the interval analysis approach [21]. It is worth noting that, in previous work of the authors [31], LS-SVM acts as a non-probabilistic surrogate model for damage detection of structures, whereas the exact FE model is implemented in this work. Furthermore, ESMA acts as a metaheuristic search algorithm and is used as the main algorithm of the model updating process.
A 31-bar planar truss, a 120-bar dome truss, and a four-story building of Phase II of the IASC-ASCE SHM benchmark are used to demonstrate the robustness of the developed procedure. The impact of various noise levels and damage severity levels on the damage detection results are also investigated using parametric computations. The study's main contribution is to provide an effective strategy for implementing a damage detection procedure that takes uncertainty into account, based on the interval optimization method and ESMA.
The article is structured as follows. Section 2 presents the process of damage detection using a correlation index. Interval analysis method is described in Section 3. The fundamental basis of SMA and proposed enhanced version are described in Sections 4 and 5. Section 6 presents the main steps for the proposed damage detection method. In Section 7, the merits of the algorithms are assessed by solving NSDI problems. Some concluding remarks are provided in Section 8.

Multiple Damage Location Assurance Criterion (MDLAC)
Structural damage detection (SDD) techniques are generally classified into two main categories. They are the dynamic and static identification techniques, which, respectively, need dynamic and static test data. Furthermore, compared with static identification approaches, dynamic identification methods have proven to be superior [36]. Among the dynamic data, the modal analysis information of a structure such as the natural frequencies and mode shapes have been widely used for SDD [36,37]. The level of correlation between measured and estimated natural frequencies and/or mode shapes can be used to locate and evaluate structural damage. When the natural frequencies are employed to identify the damage, two parameter vectors are defined. One parameter vector consists of the ratios of the first n f natural frequency changes ∆F due to structural damage, i.e.: where F h and F d denote the natural frequency vectors of the healthy and damaged structure, respectively. Another parameter vector can be defined similarly as: where F(ESP) is a natural frequency vector that can be predicted from an analytic model and elemental stiffness parameters (ESPs) ESP = [E 1 , . . . , E i , . . . , E n ] T , which represents a damage variable vector containing the elasticity modulus of structural elements (E i , i = 1, . . . , n) of all n structural elements. There are numerous approaches to evaluate the level of correlation given a pair of parameter vectors. The multiple damage location assurance criterion (MDLAC), which is presented in the following form, is an effective technique to analyze a correlation index [38]: The MDLAC analyzes two frequency change vectors, one from the inspected structure and the other from the analytical model. The MDLAC ranges from 0 to 1. It will be at its maximum when the analytical frequency vector equals the frequency vector of the damaged structure, i.e.: In addition to using the natural frequencies to identify structural damage, the mode shape vectors have been used in several studies. [39,40]. In this case, the correlation coefficient of Equation (3) can be expressed as: where ∆ϕ is the vector of the measured mode shape changes that have a dimension equal to the product of the number of measured modes and sensor locations; and ∆ϕ is the vector of the analytical mode shape changes with the same dimension as ∆ϕ. In this paper, an objective function is defined based on a combination of MDLAC(ESP) ϕ and MDLAC(ESP) F .

Interval Analysis Method for Consideration of Uncertainties
The updated ESP is prone to uncertainty since uncertainties (noises) in the obtained vibration data are unavoidable. The uncertainties in the measured modal data are considered to be independent, normally distributed random variables with zero means and a particular covariance, as indicated in [19]. As a result, the eigenvalues and mode shapes may be written in the following [19]: where the subscript 0 represents the true values, X λi and X φi indicate relative random noises in the measured frequencies and mode shapes, respectively. The mean value of vector X is zero and the standard deviation indicates the noise level. The interval bounds can provide supports for structural health monitoring under uncertain conditions [41]. To account for epistemic and aleatory uncertainties, the main concept of interval mathematics is implemented by providing the upper and lower boundaries of input parameters to generate the upper and lower boundaries of the output parameters. The stiffness reduction factor (SRF) represents the changes in the stiffness parameter for each element of the structure, as given in: where α d is ESP value of the damaged state and α u is ESP value of the undamaged state.
The intervals of the ESPs, including the natural frequencies and mode shapes, for the undamaged and damaged state, can be formulated as follows [42]: [α] ≈ λ; φ = ESP value lower bound, (9) [α] ≈ λ; φ = ESP value upper bound, Therefore, the interval bounds for each parameter can be derived as: where the subscript c is the index of damage cases, the subscripts i, j and k are the number of modes, the number of structural nodes, and the number of segments of structures, respectively. The middle values of the parameters are denoted as: where x denotes the exact values of the modal parameters (frequencies and mode shapes) and the output parameters (ESPs). The upper and lower bars, respectively, represent the upper and lower boundaries of x.
The uncertainties are coupled with the natural frequencies (λ) and mode shapes (φ) in terms of the interval bounds. The ESPs are used as variables of the optimization algorithm, and natural frequencies and mode shapes are used as in MDLAC as damage indexes. Thus, two ESMA models, which include the lower bound and upper bound analyses, are provided. The expressions are summarized in Table 1.
The variable ω represents the modal data's level of uncertainty, with varying values for natural frequencies and mode shapes. The + and − values of the uncertainty in two separate ESMA models-ESMA 1 and ESMA 2-are used to apply the boundaries (lower and upper bounds) of the input parameters of MDLAC. α and α are the values of the variables of the optimization algorithm after its execution and represent the lower and upper bounds of the predicted ESPs of damage case c.
After determining the lower and upper bounds of the ESPs, the possibility of damage existence (PoDE) is computed, and the damage severity is determined using the damage measure index (DMI).
where α I u shows the interval bound for the undamaged ESP α uk , α uk and α I d shows the interval bound for the damaged ESP α dk , α dk . A detailed explanation of these indexes is presented in the next section.
The intersection of the intervals of the damaged and undamaged ESPs on the same axis is depicted in Figure 1, with the shaded region indicating the PoDE. The PoDE ranges from 0 to 100%, with 100% indicating a relatively high possibility of damage to that particular element and 0% indicating no damage.  Figure 2 depicts the regions of damaged and undamaged ESP on two different axes. A single rectangle depicts the probable damage region, with the failure plane for both states being equal. The shaded region reflects the ESP damage. Since the ESP damage is greater than that of the undamaged ESP, the PoDE is defined as the ratio of the area of damage region to the total area of the entire region. As a result, the quantitative measurement of the PoDE can be described as below [43]: The interval value of damaged and undamaged ESP will have large variances in reality. As a result, relying just on PoDE will not give an accurate estimation of the damage. As a result, the damage measure index (DMI) is developed [26]:

Brief Description of the Slime Mold Algorithm
Chen et al. [35] presented the slime mold algorithm (SMA) based on the oscillation mode of slime mold in nature. The suggested SMA contains a number of unique features, including a mathematical model that uses adaptive weights to simulate the process of producing positive and negative feedback in slime mold propagation waves. These characteristics are based on a bio-oscillator, which creates the best path for connecting food with high exploring ability and exploitation tendency. For complete details, please refer to the paper published by [35]. An overview of the slime mold algorithm is provided below.  Figure 2 depicts the regions of damaged and undamaged ESP on two different axes. A single rectangle depicts the probable damage region, with the failure plane for both states being equal. The shaded region reflects the ESP damage. Since the ESP damage is greater than that of the undamaged ESP, the PoDE is defined as the ratio of the area of damage region to the total area of the entire region. As a result, the quantitative measurement of the PoDE can be described as below [43]:  Figure 2 depicts the regions of damaged and undamaged ESP on two d A single rectangle depicts the probable damage region, with the failure pl states being equal. The shaded region reflects the ESP damage. Since the ES greater than that of the undamaged ESP, the PoDE is defined as the ratio o damage region to the total area of the entire region. As a result, the quantitat ment of the PoDE can be described as below [43]: The interval value of damaged and undamaged ESP will have large var ality. As a result, relying just on PoDE will not give an accurate estimation of As a result, the damage measure index (DMI) is developed [26]:

Brief Description of the Slime Mold Algorithm
Chen et al. [35] presented the slime mold algorithm (SMA) based on th mode of slime mold in nature. The suggested SMA contains a number of uni including a mathematical model that uses adaptive weights to simulate th producing positive and negative feedback in slime mold propagation waves acteristics are based on a bio-oscillator, which creates the best path for con with high exploring ability and exploitation tendency. For complete details to the paper published by [35]. An overview of the slime mold algorithm is p low. The interval value of damaged and undamaged ESP will have large variances in reality. As a result, relying just on PoDE will not give an accurate estimation of the damage. As a result, the damage measure index (DMI) is developed [26]:

Brief Description of the Slime Mold Algorithm
Chen et al. [35] presented the slime mold algorithm (SMA) based on the oscillation mode of slime mold in nature. The suggested SMA contains a number of unique features, including a mathematical model that uses adaptive weights to simulate the process of producing positive and negative feedback in slime mold propagation waves. These characteristics are based on a bio-oscillator, which creates the best path for connecting food with high exploring ability and exploitation tendency. For complete details, please refer to the paper published by [35]. An overview of the slime mold algorithm is provided below.

Approach Food
The following rule is given to represent the approaching behavior of slime mold as a mathematical equation to replicate the contraction mode: where → W is the weight of slime mold, where i ∈ 1, 2, . . . , n, S(i) represents the fitness of → X. The best fitness acquired in all iterations is denoted by the DF. The → vb formula is as follows: The formula of → W is organized as follows: where condition represents that S(i) ranks first half of the population, r denotes the random value in the interval of [0, 1], bF and wF denotes the optimal and worst fitness obtained in the current iterative process, respectively. Smell Index represents the sequence of fitness values sorted (ascends in the minimum value problem).

Wrap Food
The following describes the updating position of slime mold: where LB and UB represents the lower and upper boundaries of the searching range, rand and r denote the random value in [0, 1].

Grabble Food
As the number of iterations increases, the value of → vb oscillates randomly between [−a, a] and gradually approaches zero. The value of → vc oscillates between [−1, 1] and eventually tends to zero.
The logic of SMA is shown in Figure 3 and its pseudo-code is presented in Algorithm 1.

While ( ≤ _ )
Calculate the fitness of all slime mould; , Calculate the W by Equation (23);

Enhanced Slime Mold Algorithm
This paper proposes an enhanced version of SMA to improve its capability to handle problems with a large array of variables. For structural damage evaluation, there may be a large number of variables to consider. Damaged elements and damage extents are evaluated by an optimization procedure in an optimization-based SDD problem, such that the response of a hypothesized damaged structure equals that of an actual damaged structure. The number of elements (variables) in a real structure increases when a large-scale structure such as bridge or space structure is considered. [44]. As a result, when the optimization technique seeks to minimize the objective function, it must be able to handle a large number of variables, which may slow down the convergence speed of the algorithm. To resolve this issue, this study presents a strategy as presented in the following.
When the initial population is generated, each slime mold has a location vector that reflects its position in an n-dimensional space. The elasticity modulus of the structural components is represented by each variable in this vector. The number of variables in initial stage of the SMA algorithm is considered as the total number of elements. Then, the algorithm will converge to the exact locations and severity of damages by eliminating all intact elements at each stage. A non-zero value for a variable indicates that the − th element of the structure is damaged, whereas a zero value indicates that the − th element

Enhanced Slime Mold Algorithm
This paper proposes an enhanced version of SMA to improve its capability to handle problems with a large array of variables. For structural damage evaluation, there may be a large number of variables to consider. Damaged elements and damage extents are evaluated by an optimization procedure in an optimization-based SDD problem, such that the response of a hypothesized damaged structure equals that of an actual damaged structure. The number of elements (variables) in a real structure increases when a largescale structure such as bridge or space structure is considered. [44]. As a result, when the optimization technique seeks to minimize the objective function, it must be able to handle a large number of variables, which may slow down the convergence speed of the algorithm. To resolve this issue, this study presents a strategy as presented in the following.
When the initial population is generated, each slime mold has a location vector that reflects its position in an n-dimensional space. The elasticity modulus of the structural components is represented by each variable in this vector. The number of variables in Appl. Sci. 2022, 12, 1876 9 of 21 initial stage of the SMA algorithm is considered as the total number of elements. Then, the algorithm will converge to the exact locations and severity of damages by eliminating all intact elements at each stage. A non-zero value for a variable indicates that the i − th element of the structure is damaged, whereas a zero value indicates that the i − th element of the structure is intact. If the variables with near zero values (SRF i ≤ 0.05) remain the same for 10 iterations, it can be concluded that these variables belong to the intact elements set and therefore these near-zero variables will be eliminated in the next iteration of the SMA.
The objective function is defined as an unconstrained optimization problem, as follows: F(ESP) is the minimization problem, whereas E min and E max are the damage vector's lower and upper limits, respectively. The structure's physical behavior necessitates the use of boundaries. The damage variables are estimated using an optimization approach and by solving Equation (26).

Main Steps of the Proposed NSDI Method
The following are the main steps in the suggested NSDI approach utilizing the ESMA algorithm: (a) Create two parallel ESMA model based on Table 1 and set the initial number of variables equal to the total number of elements in structure. Equations (17) and (18).

Numerical Results of NSDI
Three structures are used as numerical examples in this work to illustrate the robustness and accuracy of the proposed NSDI approach. These structures are: 1.
four-story steel structure model.
Damage to the structure is modeled as a relative reduction in the elasticity modulus of each individual element, with the mass matrix assumed to remain constant. For the optimization process, the number of slime molds for ESMA was fixed to 20 for each run along a maximum of 100 iterations. The fit-test design in the current iteration is compared to that of 20 prior iterations as a measure of convergence for engineering optimization problems. As a result, if the difference between these two results is smaller than a certain tolerance value, it is recorded as converged. The desired value for this example is set as 10 −3 . If the converge does not occur, the algorithm will be terminated by implementing a fixed maximum number of iterations.

31-Bar Planar Truss
The standard finite element (FE) approach is used to model the 31-bar planar truss in Figure 4 without internal nodes, resulting in 25 degrees of freedom [45]. The first five vibrating modes are used for NSDI in this case. The elasticity modulus and density of the material are 70 GPa and 2770 kg/m, respectively. Two distinct damage scenarios are implemented in the structure, as shown in Table 2, and the suggested approach is evaluated in each case. material are 70 GPa and 2770 kg/m, respectively. Two distinct damage scenarios are implemented in the structure, as shown in Table 2, and the suggested approach is evaluated in each case.
The impact of noise on the accuracy of NSDI based on modal data is studied in this section. The uncertainties of the measured mode shapes of structures are typically larger than the uncertainties of the frequencies in modal testing [1]. In this study, uncertain frequencies and mode shapes are considered as the normal distributed random variables defined in Equations (6) and (7). The mean values of the relative random noises are zeros, and standard deviations indicate the noise level.   Tables 3 and 4 show PODEs and DMIs of the damage scenarios 1 and 2, respectively. The noise level is set as 2% and 5% in the frequencies and mode shapes, respectively. In scenario 1, higher PoDE values are obtained at element numbers 1 and 2 compared to the undamaged elements. Moreover, the PoDE values of elements 1 and 2 are higher than those of the undamaged elements. The DMI value of element 1 is similarly higher than that of element 2; both of these elements are correct damage locations with varying severity conditions. The same condition can be seen for scenario 2, where higher PoDE values observed at segments 11 and 25 and the DMI values are also higher at the elements with higher severity.   The impact of noise on the accuracy of NSDI based on modal data is studied in this section. The uncertainties of the measured mode shapes of structures are typically larger than the uncertainties of the frequencies in modal testing [1]. In this study, uncertain frequencies and mode shapes are considered as the normal distributed random variables defined in Equations (6) and (7). The mean values of the relative random noises X λi are zeros, and standard deviations ξ λ indicate the noise level. Tables 3 and 4 show PODEs and DMIs of the damage scenarios 1 and 2, respectively. The noise level is set as 2% and 5% in the frequencies and mode shapes, respectively. In scenario 1, higher PoDE values are obtained at element numbers 1 and 2 compared to the undamaged elements. Moreover, the PoDE values of elements 1 and 2 are higher than those of the undamaged elements. The DMI value of element 1 is similarly higher than that of element 2; both of these elements are correct damage locations with varying severity conditions. The same condition can be seen for scenario 2, where higher PoDE values observed at segments 11 and 25 and the DMI values are also higher at the elements with higher severity.  To further prove the capability of the proposed method, Tables 5 and 6 compare the PoDE values obtained by the proposed method with the PDE values for the same damage cases calculated using a probability-based method developed by Ghiasi et al. [33,46].  Ghiasi et al. [33] uses the probability density function of the existence of undamaged and damaged states to determine the probability of damage existence (PDE). The main algorithm consists of three interconnected loops based on coupling metamodeling techniques with a metaheuristic optimization algorithm for probability-based damage detection (PBDD) of structures.
These results show that the proposed technique outperforms the PBDD method by providing a more meaningful damage severity indicator than the statistical surrogate model method, which simply displays damage severity in terms of the probabilities. [21,26]. Furthermore, compared to the aforementioned surrogate methods [33,47], the proposed non-probabilistic method provides an advantage since it does not require surrogate models to predict the upper and lower boundaries of the ESP values. Surrogate models impose some uncertainties on the process, which increase the prediction error. The standard deviation of the probability density function (PDF) used to calculate the probability of damage existence is typically increased by the prediction error in the statistical surrogate model (PDE) [47]. This higher standard deviation may result in lower probability values when detecting damage in structures.
On other hand, in both scenarios, the PoDE is demonstrated to be a more accurate damage index and results in smaller errors. For example, in scenario 1, the suggested approach reveals a 0% PoDE value in element 10, which is undamaged, compared to a 7% PDE value. In scenario 2, the undamaged element 26 shows 0% damage, whereas the probability-based technique predicts 12% damage. It is also observed that for both scenarios, the proposed method provides higher PoDE values at the damaged elements compared to the PDE value. The main reason is that the proposed non-probabilistic method utilizes only two models to predict the upper and lower boundaries of the ESP values. Furthermore, it implemented the exact FE model. In comparison, the probabilistic method for the same example needs four surrogate models to calculate the probability density function of the damage [47], which leads to greater prediction errors. Furthermore, as can be seen in Tables 5 and 6, the proposed methodology reduces the computation time by at least 37%.

120-Bar Dome Truss
The second example is a 120-bar dome truss, as illustrated in Figure 5 [48]. Appl On other hand, in both scenarios, the PoDE is demonstrated to be a more accurate damage index and results in smaller errors. For example, in scenario 1, the suggested approach reveals a 0% PoDE value in element 10, which is undamaged, compared to a 7% PDE value. In scenario 2, the undamaged element 26 shows 0% damage, whereas the probability-based technique predicts 12% damage. It is also observed that for both scenarios, the proposed method provides higher PoDE values at the damaged elements compared to the PDE value. The main reason is that the proposed non-probabilistic method utilizes only two models to predict the upper and lower boundaries of the ESP values. Furthermore, it implemented the exact FE model. In comparison, the probabilistic method for the same example needs four surrogate models to calculate the probability density function of the damage [47], which leads to greater prediction errors. Furthermore, as can be seen in Tables 5 and 6, the proposed methodology reduces the computation time by at least 37%.

120-Bar Dome Truss
The second example is a 120-bar dome truss, as illustrated in Figure 5 [42] was used. As described in Table 7, two damage scenarios are considered. Figures 6 and 7a and Tables 8 and 9 show the performance of the proposed method. It is worth noting that the noise level is set as 2% and 5% in the frequencies and mode shapes, respectively. The dome's diameter and height are 31.78 m and 7 m, respectively. A seamless steel pipe with a modulus of elasticity of 30,450 ksi (210,000 MPa) and material density of 0.288 3 (7971.810 3 ) is used. The pipes have an external diameter of 0.2 m and a thickness of 0.006 m. For FE modeling and analysis, OpenSees software [42] was used. As described in Table 7, two damage scenarios are considered. Figures 6 and 7a and Tables 8 and 9 show the performance of the proposed method. It is worth noting that the noise level is set as 2% and 5% in the frequencies and mode shapes, respectively.  The results show that although the measured modal data were noise corrupted, the proposed approach provides a reliable prediction of damage occurrence. Furthermore, ESMA can effectively find the correct locations and severity of the damage. For a larger noise level, that is = 10%, the resultant PoDE values of damage scenario 2 are shown in Figure 7b and Table 10.
The numerical results show that the proposed method performs well in detecting the location and severity of the damage scenarios even with a large noise level. The PoDE values of the damaged elements are higher in NSDI approaches, whereas the PoDE values of the undamaged elements are lower. The results show that although the measured modal data were noise corrupted, the proposed approach provides a reliable prediction of damage occurrence. Furthermore, ESMA can effectively find the correct locations and severity of the damage. For a larger noise level, that is ξ λ = 10%, the resultant PoDE values of damage scenario 2 are shown in Figure 7b and Table 10.    The numerical results show that the proposed method performs well in detecting the location and severity of the damage scenarios even with a large noise level. The PoDE values of the damaged elements are higher in NSDI approaches, whereas the PoDE values of the undamaged elements are lower.

4 Story Steel Structure Model
The third example is a 3D modular structure based on Phase II of the IASC-ASCE SHM benchmark four-story building [49,50]. As illustrated in Figure 8, the model is built without side braces along the four floors. The structure has 12 degrees of freedom (DOFs) and 84 members. This structure has a base plan of 2.5×2.5 m and height of 3.6 m. The members are hot rolled grade 300 W steel with a nominal yield stress 300 MPa (42.6 kpsi). This complex model can effectively illustrate the efficiency of structural damage tracking using the proposed non-probabilistic method. As indicated in Figure 9, one damage situation is considered by replicating the damage in element 7 by reducing its Young's modulus by 25%. This complex model can effectively illustrate the efficiency of structural damage tracking using the proposed non-probabilistic method. As indicated in Figure 9, one damage situation is considered by replicating the damage in element 7 by reducing its Young's modulus by 25%.
Appl. Sci. 2022, 11, x FOR PEER REVIEW Figure 9. The damage scenario. Table 11 show PODEs, DMIs, and PDEs for the damage scenario. Further comparison between the proposed ESMA and the original SMA is shown. He simulated mode shapes and frequencies consist of 5% and 2% measurement noise tively.
The maximum numbers of function evaluations (NFEs) shown in Table 11 that determining the global optimum with ESMA is substantially less computa costly than with SMA. Based on the results displayed in Table 11, it can be determined that SMA ability to detect damaged elements with high certainty, whereas the resulting PDE and DMI for the undamaged elements is higher in comparison with ESMA. Furth  Table 11 show PODEs, DMIs, and PDEs for the damage scenario. Furthermore, a comparison between the proposed ESMA and the original SMA is shown. Herein, the simulated mode shapes and frequencies consist of 5% and 2% measurement noise, respectively. The maximum numbers of function evaluations (NFEs) shown in Table 11 indicate that determining the global optimum with ESMA is substantially less computationally costly than with SMA.
Based on the results displayed in Table 11, it can be determined that SMA has the ability to detect damaged elements with high certainty, whereas the resulting PDE, PoDE, and DMI for the undamaged elements is higher in comparison with ESMA. Furthermore, as can be seen in Table 11, the proposed methodology reduces the computation time by at least 30%.

Conclusions
Using the most recently introduced optimization algorithm (SMA), a unique approach for NSDI was proposed in this study as a one-stage procedure. There are always uncertainties in the measured modal parameters in typical applications of health monitoring approaches, which can lead to unreliable and inaccurate structural damage prediction. As a result, uncertainties are included in the developed NSDI process. An SMA method is updated to increase its performance in dealing with problems involving a large number of variables. An interval analysis is adopted for use with the ESMA to consider the uncertainties using the interval bounds of the uncertainties in the damage variables of the metaheuristic optimization model.
The numerical results demonstrate that the proposed approach is capable of accurately detecting the location and severity of various damage scenarios. Based on the results it can be concluded that the proposed interval optimization method procedure can obtain accurate PoDE and DMI and can provide robust damage identification results. In summary, the proposed non-probabilistic method is able to provide accurate damage detection results and has a small prediction cost since it does not require the utilization of surrogate models. Despite these advantages, this approach may need further improvements to be implemented in the real-time recorded response of a structure, such as adding the structural identification (SI) stage, in order to convert the acceleration response to the frequency response. Consideration of uncertainty in loadings will be pursued in our upcoming papers.
Author Contributions: R.G., the first author, was responsible for developing the methodology introduced in this work and carrying out the majority of the research project. Both analytical derivations and statistical methods were closely tested and confirmed by M.N. and W.A.A., who suggested the procedures utilized. W.A.A. provided his expertise and technical background in AI and assisted R.G. with carrying out the numerical analyses and implementation of algorithm. A.S. with his expertise in statistical methods and feature extraction worked closely with R.G. and W.A.A. throughout the project. T.W. worked closely with R.G., W.A.A., F.P. and S.-C.K. in the reviewing and the editing phase and final verification of the analytical derivations. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest.