Optimization of the Curved Metal Damper to Improve Structural Energy Dissipation Capacity

: Structural curved metal dampers are implemented in various applications to mitigate the damages at a speciﬁc area efﬁciently. A stable and saturated hysteretic behavior for the in-plane direction is dependent on the shape of a curved-shaped damper. However, it has been experimentally shown that the hysteretic behavior in the conventional curved-shaped damper is unstable, mainly as a result of bi-directional deformations. Therefore, it is necessary to conduct shape optimization for curved dampers to enhance their hysteretic behavior and energy dissipation capability. In this study, the ﬁnite element (FE) model built in ABAQUS, is utilized to obtain optimal shape for the curved-shaped damper. The effectiveness of the model is checked by comparisons of the FE model and experimental results. The parameters for the optimization include the curved length and shape of the damper, and the improved approach is conducted by investigating the curved sections. In addition, the design parameters are represented by B-spline curves (to ensure enhanced system performance), regression analysis is implemented to derive optimization formulations considering energy dissipation, constitutive material model, and cumulative plastic strain. Results determine that the energy dissipation capacity of the curved steel damper could be improved by 32% using shape optimization techniques compared to the conventional dampers. Ultimately, the study proposes simple optimal shapes for further implementations in practical designs. Author Contributions: Conceptualization, project supervision, writing—review and editing Y.-C.K.; formal analysis, investigation, writing—original draft S.J.M., A.F. and J.-W.H.; writing— review and editing Y.-C.K., S.J.M., A.F., J.-W.H., I.M. and P.O.A.; writing—review and J.-W.H. I.M.; writing—review and editing, acquisition,


Introduction
In recent decades, the enhancement of structural resistance to seismic load inducing significant damages has been achieved by using passive energy dissipation (PED), which also relatively reduced the initial non-structural damages. The PED, operating as structural fuses, concentrates the damages in a specific part of the structures, which keeps the primary structure in an elastic state [1]. The use of the structural fuse has shown several benefits. One of the major benefits of damper usage is the component's replaceability after an event occurred; thus, the repair works are conducted effectively without affecting the structural functionality [2][3][4][5][6][7][8][9].
The flexural yielding mechanism was originally addressed by the use of metallic fuse in various applications (e.g., Added Damping and Stiffness [10]). An added damping and stiffness showed in Figure 1 consisting of X-shaped steel plates acting as bracing systems.
For this configuration, yielding is only possible across the entire device length under severe loads causing double curvature inelastic deformations [11].
Buildings 2022, 12, x FOR PEER REVIEW 2 of 13 For this configuration, yielding is only possible across the entire device length under severe loads causing double curvature inelastic deformations [11].
Metallic fuses are extensively applied to serve as working mechanisms: flexural, axial, and shear yielding [13][14][15][16][17]. The shear panel fuse is made of a steel plate shear yielding mechanism, and it is mostly used for bracing hysteretic damping devices. An example of this system is presented in Figure 1. Several previous studies related to the use of steel plate shear walls and link beams have shown that the steel shear mechanism exhibits high strength and stiffness [18][19][20][21]. This mechanism demonstrates severe out-of-plane buckling Figure 1. Configuration of (a) ADAS damper (b) TADS damper [12].
Metallic fuses are extensively applied to serve as working mechanisms: flexural, axial, and shear yielding [13][14][15][16][17]. The shear panel fuse is made of a steel plate shear yielding mechanism, and it is mostly used for bracing hysteretic damping devices. An example of this system is presented in Figure 1. Several previous studies related to the use of steel plate shear walls and link beams have shown that the steel shear mechanism exhibits high strength and stiffness [18][19][20][21]. This mechanism demonstrates severe out-of-plane buckling behavior, which is a result of thin plate geometry, resulting in a reduction of hysteresis pinching, low ductility, and shear force. Procedures could be ascertained by either having cut-outs into a thin plate or adding stiffeners to the boundary of plates [11].
The cyclic behavior of uniform and non-uniform steel stripe dampers has been studied previously [22][23][24][25]. Figure 2 represent various damper shapes, which include (a) tapered strip, (b) dumbbell-shaped strip, and (c) hourglass-shaped strip. Test results indicate that the damage was not concentrated at the ends causing crack propagations compared to the case of the uniform sample [22]. Using appropriate design geometries for steel dampers could improve the cumulative ductility by 1.13-1.75, and energy dissipation by 1.27-2.36 [22] according to previous studies. behavior, which is a result of thin plate geometry, resulting in a reduction of hysteresis pinching, low ductility, and shear force. Procedures could be ascertained by either having cut-outs into a thin plate or adding stiffeners to the boundary of plates [11]. The cyclic behavior of uniform and non-uniform steel stripe dampers has been studied previously [22][23][24][25]. Figure 2 represent various damper shapes, which include (a) tapered strip, (b) dumbbell-shaped strip, and (c) hourglass-shaped strip.Test results indicate that the damage was not concentrated at the ends causing crack propagations compared to the case of the uniform sample [22]. Using appropriate design geometries for steel dampers could improve the cumulative ductility by 1.13-1.75, and energy dissipation by 1.27-2.36 [22] according to previous studies. Iso-geometric analysis (IGA) has been identified as effective for modeling complex acoustic domain geometries and ensuring the accuracy of the model over the optimization phase [26]. For this approach, CAD geometries and the unknown variables are defined by using the spline basis functions. B-Splines and Non-Uniform Rational B-Splines (NURBs) are among the most commonly used spines [27]. It is noted that compared with the standard FE technique, the IGA method exhibits higher accuracy as a result of exact geometry evaluations, and higher basis function continuity [28,29]. In addition, computational errors could be reduced by using this approach refining the NURBs and ensuring a higher rate of solution approximations, and continuity [27,30].
Advancement in computer technology has aided numerous optimization algorithms related to structural optimization [31]. Algorithms such as the simulated annealing (SA), genetic algorithm (GA), and multi-level or multi-disciplinary optimization techniques have been developed previously for improving the low cycle fatigue behavior of the steel links [32][33][34][35][36]. In this study, the shape optimization approach using the SA method for enhancement of low cycle fatigue behavior steel shear panel dampers (SSPDs) is used for optimization purposes.
This research paper consists of two phases. The first phase is optimization analysis, and the second phase is FEM analysis. These phases were connected to each other. In the Iso-geometric analysis (IGA) has been identified as effective for modeling complex acoustic domain geometries and ensuring the accuracy of the model over the optimization phase [26]. For this approach, CAD geometries and the unknown variables are defined by using the spline basis functions. B-Splines and Non-Uniform Rational B-Splines (NURBs) are among the most commonly used spines [27]. It is noted that compared with the standard FE technique, the IGA method exhibits higher accuracy as a result of exact geometry evaluations, and higher basis function continuity [28,29]. In addition, computational errors could be reduced by using this approach refining the NURBs and ensuring a higher rate of solution approximations, and continuity [27,30].
Advancement in computer technology has aided numerous optimization algorithms related to structural optimization [31]. Algorithms such as the simulated annealing (SA), genetic algorithm (GA), and multi-level or multi-disciplinary optimization techniques have been developed previously for improving the low cycle fatigue behavior of the steel links [32][33][34][35][36]. In this study, the shape optimization approach using the SA method for enhancement of low cycle fatigue behavior steel shear panel dampers (SSPDs) is used for optimization purposes.
This research paper consists of two phases. The first phase is optimization analysis, and the second phase is FEM analysis. These phases were connected to each other. In the optimization phase, IGA analysis is carried out and the isoparametric formulation is applied in the FEM phase. To clarify this point, it was added to the text of the article. In this study, the finite element (FE) model developed in ABAQUS, is utilized to obtain optimal shape for the curved-shaped damper. The shape optimization techniques are implemented to improve the steel damper energy dissipation and the effectiveness of the model is compared with previous initial models.

B-Spline Basis Functions
By using the B-splines in parametric space, a non-uniform rational B-spline function (NURBS) was derived based on Equations (1) and (2) [27]: and for 1 ≤ p.
For which, the basis function for NURBS is given by Equation (3): where w i = i th weight. Based on control points Bi, the NURBS curve was parameterized according to Equation (4): To ensure an effective steel shear panel dampers behavior, both end lengths and the middle section must be properly monitored. This implies that both end lengths and middle sections are among the most main influential dimensions for a typical steel damper. Thus, these dimensions are represented as the main design variables for the optimization process in this study (see Figure 1). To ensure effective energy dissipation and deformation ability of the steel shear panel dampers during cyclic loading, the optimization is performed based on maximizing the dissipated energy for plastic deformations (E d ) over the maximum equivalent plastic strain (PEEQ max ) which is determined as the optimization objective function. The objective function for the steel shear panel damper shape optimization problem is defined as indicated in Equation (5) [3].

Shape Optimization Based on Whale Optimization Algorithm
The metal curved damper is known to perform effectively based on two factors: the middle section (R) and the length of the damper end section (a) which were previously shown that as major parameters affecting the effective performance (see Figure 3). In this approach, humpback whales' social behavior was used for the whale op zation algorithm (WOA) development [37]. The algorithm is proposed by [38], adopting the bubble-net hunting procedures. The WOA algorithm has been mainly for continuous function optimization problems since it was proposed because it is e principle, operation, and implementation, and it is superior to other optimization rithms in terms of solution accuracy and algorithm stability. In detail, it is mainly ap to optimization techniques in applied mathematics and civil engineering and sk structure design in construction engineering [39]. The principle of WOA is based o preference of humpback whales to hunt a krill school or small fish around the wate face. This is achieved by causing a series of bubbles along a circle. It is noted that a in this algorithm corresponds to one feasible solution and the prey refers to the op solution. The targeted prey is encircled by the humpback whales simultaneously wi shrinking mechanism and through a spiral-shaped route. Therefore, each whale sumed to select either the shrinking mechanism or the spiral-shaped approach to m position by the probability of 50% [37] which is summarized in Equation (6).
where X(t+1) denotes updated next iteration (t + 1) position and p is a random number in * ( ) is the best solution obtained so far, is a random value greater than 1 or less th to force each search agent (whale) to move, respectively, toward or far away from a ence whale, D′ is the distance of the current position to the best one so far, t stands f current iteration, b = a constant for defining the shape of the logarithmic spiral, l is dom number in [−1, 1].
A set of random feasible solutions is made at the start of the WOA algorithm initial position of the whales from the population. At each iteration, WOA makes tions between a spiral or shrinking movements, based on the specific p-value. I shrinking mechanism, the whales update their positions based on random selectio ploration phase) or the best position obtained so far (exploitation phase). Upon a sa tory performance, in terms of fulfilling the completion criterion (based on maximum ation numbers), the WOA algorithm is then terminated accordingly [40].

Finite Element Modeling Methodology
In this study, both geometric and material nonlinear analyses were carried o In this approach, humpback whales' social behavior was used for the whale optimization algorithm (WOA) development [37]. The algorithm is proposed by [38], after adopting the bubble-net hunting procedures. The WOA algorithm has been mainly used for continuous function optimization problems since it was proposed because it is easy in principle, operation, and implementation, and it is superior to other optimization algorithms in terms of solution accuracy and algorithm stability. In detail, it is mainly applied to optimization techniques in applied mathematics and civil engineering and skeletal structure design in construction engineering [39]. The principle of WOA is based on the preference of humpback whales to hunt a krill school or small fish around the water surface. This is achieved by causing a series of bubbles along a circle. It is noted that a whale in this algorithm corresponds to one feasible solution and the prey refers to the optimal solution. The targeted prey is encircled by the humpback whales simultaneously within a shrinking mechanism and through a spiral-shaped route. Therefore, each whale is assumed to select either the shrinking mechanism or the spiral-shaped approach to main its position by the probability of 50% [37] which is summarized in Equation (6).
where X (t+1) denotes updated next iteration (t + 1) position and p is a random number in [0, 1], X * (t) is the best solution obtained so far, A is a random value greater than 1 or less than −1 to force each search agent (whale) to move, respectively, toward or far away from a reference whale, D is the distance of the current position to the best one so far, t stands for the current iteration, b = a constant for defining the shape of the logarithmic spiral, l is a random number in [−1, 1]. A set of random feasible solutions is made at the start of the WOA algorithm as an initial position of the whales from the population. At each iteration, WOA makes selections between a spiral or shrinking movements, based on the specific p-value. In the shrinking mechanism, the whales update their positions based on random selection (exploration phase) or the best position obtained so far (exploitation phase). Upon a satisfactory performance, in terms of fulfilling the completion criterion (based on maximum iteration numbers), the WOA algorithm is then terminated accordingly [40].

Finite Element Modeling Methodology
In this study, both geometric and material nonlinear analyses were carried out for the optimization in ABAQUS. For this purpose, FEM updating processes utilize only a single strip (SSD). The buckling response of the plate was mobilized by the inclusion of the initial geometric imperfections in numerical models. Thus, the most conservative initial geometric imperfection is the first buckling mode shape [41]. Based on the previous study [42], eigenvalue buckling analysis of the SSD model and the first Eigen shape under a shear force was performed using a maximum amplitude of h/100 and validated in a verification study, and subsequently incorporated in the FEM analysis. It should be noted that the described material properties during the verification study were included in the FEM analysis.
The accumulated damages are assessed based on the cyclic loading with incremental amplitudes according to FEMA-461 provisions. Based on a conservative selection, the lowest shear deformation damage state (∆ 0 = 0.1%) was taken as 0.001 radians, a value that is much lower than the SSD model buckling or yielding. For each of the increments, there are two identical cycles, and displacement amplitudes were increased by 1.4 for every loading iteration as is shown in Figure 4. Other studies [23,43] have also utilized a similar scheme of loading conditions for assessing slit dampers.
, 12, x FOR PEER REVIEW 6 of 13 the initial geometric imperfections in numerical models. Thus, the most conservative initial geometric imperfection is the first buckling mode shape [41]. Based on the previous study [42], eigenvalue buckling analysis of the SSD model and the first Eigen shape under a shear force was performed using a maximum amplitude of h/100 and validated in a verification study, and subsequently incorporated in the FEM analysis. It should be noted that the described material properties during the verification study were included in the FEM analysis. The accumulated damages are assessed based on the cyclic loading with incremental amplitudes according to FEMA-461 provisions. Based on a conservative selection, the lowest shear deformation damage state (Δ0 = 0.1%) was taken as 0.001 radians, a value that is much lower than the SSD model buckling or yielding. For each of the increments, there are two identical cycles, and displacement amplitudes were increased by 1.4 for every loading iteration as is shown in Figure 4. Other studies [23,43] have also utilized a similar scheme of loading conditions for assessing slit dampers.

Verification Study
Validation of FEM methodology is performed by using the experimental data [44]. In Figure 5, the B09-56 specimen geometrical dimensions are illustrated. In addition, for the FEM analysis, the true stress and logarithmic strain material curves obtained from the coupon test are considered (Plate C [44]). The engineering yielding stress is 347 MPa, while the ultimate stress of the material is considered to be 435 MPa. At rupture displacement of 25 mm. Along with the same lines, a shell element with three nodes and the reduced integration mode (S3R) shell element was utilized in the FEM, with a maximum size of 3mm. Also, the FEM comprised of a global geometrical imperfection (linked to the SSD model's first local buckling Eigen shape and seven strips), and a one percent strip height amplitude strip (3.5 mm) [41].

Verification Study
Validation of FEM methodology is performed by using the experimental data [44]. In Figure 5, the B09-56 specimen geometrical dimensions are illustrated. In addition, for the FEM analysis, the true stress and logarithmic strain material curves obtained from the coupon test are considered (Plate C [44]). The engineering yielding stress is 347 MPa, while the ultimate stress of the material is considered to be 435 MPa. At rupture displacement of 25 mm. Along with the same lines, a shell element with three nodes and the reduced integration mode (S3R) shell element was utilized in the FEM, with a maximum size of 3 mm. Also, the FEM comprised of a global geometrical imperfection (linked to the SSD model's first local buckling Eigen shape and seven strips), and a one percent strip height amplitude strip (3.5 mm) [41].  [41] A comparison between the results of the experimental test and the FEM model afte loading indicated more than 91% agreement for predicting experimental results. [41]. It i shown that the FEM captures the mode of behavior, the ultimate buckled shape, and th buckling mode shape.

Optimization Process
The maximum stiffness and the full stress state could be used as the optimizatio objective function for the general topology optimization analysis [45]. Alternately, doubl optimization results of specified volume reduction ratio can be obtained, and compared with the optimization results for different volume reduction rates in order to obtain th final optimization scheme. Figure 6 illustrates the specific optimization process. In this study, a set of 500 FE models is considered to optimize limit state functions Various geometrical shapes are considered during the investigation of nonlinear load-de formation performance, and the pushover behavior of the structural shear links.  A comparison between the results of the experimental test and the FEM model after loading indicated more than 91% agreement for predicting experimental results. [41]. It is shown that the FEM captures the mode of behavior, the ultimate buckled shape, and the buckling mode shape.

Optimization Process
The maximum stiffness and the full stress state could be used as the optimization objective function for the general topology optimization analysis [45]. Alternately, double optimization results of specified volume reduction ratio can be obtained, and compared with the optimization results for different volume reduction rates in order to obtain the final optimization scheme. Figure 6 illustrates the specific optimization process.
Buildings 2022, 12, x FOR PEER REVIEW 7 of 1 Figure 5. B09−56 specimen geometrical dimensions (in mm) [41] A comparison between the results of the experimental test and the FEM model afte loading indicated more than 91% agreement for predicting experimental results. [41]. It i shown that the FEM captures the mode of behavior, the ultimate buckled shape, and th buckling mode shape.

Optimization Process
The maximum stiffness and the full stress state could be used as the optimization objective function for the general topology optimization analysis [45]. Alternately, doubl optimization results of specified volume reduction ratio can be obtained, and compared with the optimization results for different volume reduction rates in order to obtain th final optimization scheme. Figure 6 illustrates the specific optimization process. In this study, a set of 500 FE models is considered to optimize limit state functions Various geometrical shapes are considered during the investigation of nonlinear load-de formation performance, and the pushover behavior of the structural shear links.

Structural analysis
Finite element (ABAQUS) to evaluate objective functions and the constraints • Minimizing total strain energy (ABAQUS and Matlab)  In this study, a set of 500 FE models is considered to optimize limit state functions. Various geometrical shapes are considered during the investigation of nonlinear loaddeformation performance, and the pushover behavior of the structural shear links.

Applying design constraints
The following geometrical parameters are considered accordingly: links length (L) of 0.5 m for all the models; slenderness ratio varying between 10-50; damper taper ratio (a/R), which describes the inelastic behavior of links varying between 0.017 and 3; The R/L ratio which is between 0.2 and 0.95 indicating a major optimization parameter for the links' performance. It is noted that damper optimization is typically conducted to improve the energy dissipation capability with less concentrated plastic strain areas yielding to more steady stress distribution in larger drift values [46].
As is shown in Figure 7, the boundary condition of the steel shear panel damper at the bottom is fixed, and also restrained against the out-of-plane displacements at the top. An initial eigenvalue analysis is utilized for applying the initial imperfection displacement (L/250) based as a result of the first buckling mode. Following this process, the displacement-control load is applied on the link top edge, and cyclic results are monitored.
The yield stress for the material model is 300 MPa, and the ultimate stress is 345 MPa. A linear hardening between the ultimate point and yielding point was assumed. As previously indicated, the monotonic loading is applied at the top boundary as shown in Figure 7 for each model which is used for extraction of the pushover curve for determination of the dissipated energy and maximum plastic strains. The sensitivity of the mesh was done to achieve high precision in plastic strain and pushover results. The shear locking and hour glassing issues are monitored and prevented. Moreover, FEM accuracy was ensured by using a minimum of 8, 20, and 30 elements, at the middle section, top length, and inclined side lengths of steel shear steel dampers [3].
Buildings 2022, 12, x FOR PEER REVIEW 8 The following geometrical parameters are considered accordingly: links length ( 0.5 m for all the models; slenderness ratio varying between 10-50; damper taper (a/R), which describes the inelastic behavior of links varying between 0.017 and 3; Th ratio which is between 0.2 and 0.95 indicating a major optimization parameter fo links' performance. It is noted that damper optimization is typically conducted to imp the energy dissipation capability with less concentrated plastic strain areas yieldin more steady stress distribution in larger drift values [46].
As is shown in Figure 7, the boundary condition of the steel shear panel damp the bottom is fixed, and also restrained against the out-of-plane displacements at the An initial eigenvalue analysis is utilized for applying the initial imperfection displace (L/250) based as a result of the first buckling mode. Following this process, the disp ment-control load is applied on the link top edge, and cyclic results are monitored.
The yield stress for the material model is 300 MPa, and the ultimate stress is 345 A linear hardening between the ultimate point and yielding point was assumed. As viously indicated, the monotonic loading is applied at the top boundary as shown in Fig for each model which is used for extraction of the pushover curve for determination o dissipated energy and maximum plastic strains. The sensitivity of the mesh was do achieve high precision in plastic strain and pushover results. The shear locking and glassing issues are monitored and prevented. Moreover, FEM accuracy was ensure using a minimum of 8, 20, and 30 elements, at the middle section, top length, and inc side lengths of steel shear steel dampers [3]. In this study, the synergy of system integration of both MATLAB and ABAQU explored for the optimization procedure. First, a python script is developed by ABA and performed the analysis. The outcome of the analysis, using the optimization tions, is added to MATLAB for WOA application and variable updates resulting in a model for the following iterations. With this process, it is observed that the more the n ber of iterations performed, the higher the accuracy of the optimal function. It is n that both FE analysis and the WOA algorithm are combined based on the Python s implementation. In this study, the synergy of system integration of both MATLAB and ABAQUS is explored for the optimization procedure. First, a python script is developed by ABAQUS and performed the analysis. The outcome of the analysis, using the optimization functions, is added to MATLAB for WOA application and variable updates resulting in a new model for the following iterations. With this process, it is observed that the more the number of iterations performed, the higher the accuracy of the optimal function. It is noted that both FE analysis and the WOA algorithm are combined based on the Python script implementation. Figure 8 shows the plastic strain distribution representing the steel shear damper model at a 5% drift ratio. The FE steel shear panel damper model is compared with the optimized model based on the proposed method. The optimization functions are generated from WOA algorithms to achieve less fracture and maximize the energy dissipation capability. In addition, a previously defined loading condition was applied to the optimized models. As a result of the analysis, the distribution of the maximum plastic strain is generally uniformly distributed around 10% at the end of the damper. It is shown that WOA algorithm optimization could lead to a 51% reduction in the plastic strains at a 5% drift ratio compared to the initial models resulting in better resistance at higher drifts. Moreover, Figure 8 shows the optimal model plastic strain behavior which indicates low plastic strain values, especially around the geometrical changes. Overall, there was uniformity in the plastic strain distribution within the models in the optimal model compared to other models. This implies that the method of optimization permits more energy dissipation at a far position from the sharper angles and critical points with high strain concentration. Based on behavior investigations, this study shows four possible types of behavior. In the first group, plastic strains are concentrated at steel shear panel damper ends (the model with the highest PEEQ compared with other dampers at similar drift). The second group comprised models having plastic strains concentrated at the middle (position of shear yielding) [18,47]. In the third group, the buckling limit state occurs leading to the lesser ability for load-bearing resistance during medium to large drifts. Finally, the last group is the optimal model, where there is a simultaneous occurrence of shear and flexural yielding without attaining a brittle limit state. In this model, fracture occurrence is avoided and energy dissipation is optimum, thereby achieving the optimized limit state function.

Discussion and Results
Buildings 2022, 12, x FOR PEER REVIEW 9 of 1 optimized model based on the proposed method. The optimization functions are gener ated from WOA algorithms to achieve less fracture and maximize the energy dissipatio capability. In addition, a previously defined loading condition was applied to the opti mized models. As a result of the analysis, the distribution of the maximum plastic strai is generally uniformly distributed around 10% at the end of the damper. It is shown tha WOA algorithm optimization could lead to a 51% reduction in the plastic strains at a 5% drift ratio compared to the initial models resulting in better resistance at higher drifts Moreover, Figure 8 shows the optimal model plastic strain behavior which indicates low plastic strain values, especially around the geometrical changes. Overall, there was uni formity in the plastic strain distribution within the models in the optimal model compared to other models. This implies that the method of optimization permits more energy dissi pation at a far position from the sharper angles and critical points with high strain con centration. Based on behavior investigations, this study shows four possible types of be havior. In the first group, plastic strains are concentrated at steel shear panel damper end (the model with the highest PEEQ compared with other dampers at similar drift). Th second group comprised models having plastic strains concentrated at the middle (posi tion of shear yielding) [18,47]. In the third group, the buckling limit state occurs leadin to the lesser ability for load-bearing resistance during medium to large drifts. Finally, th last group is the optimal model, where there is a simultaneous occurrence of shear and flexural yielding without attaining a brittle limit state. In this model, fracture occurrenc is avoided and energy dissipation is optimum, thereby achieving the optimized limit stat function. Similarly, the stress distribution of selected models following a 5% drift ratio appli cation is shown in Figure 9. As a result of the analysis, the maximum stress distributio was found to be uniformly distributed at about 350 Mpa at the end of the damper, simila to the distribution of maximum plastic strain. It is shown that the stress distribution of th optimal model is steadily demonstrated through the length of the link due to better mate rials implementation and higher energy dissipation. The study also showed that the geo metrical properties of the optimal model are responsible for the mode of transition from flexure yielding to shear yielding, in accordance with the transitions equation provided by Farzampour and Eatherton [18,47]. In addition, the optimization analysis has show that the models with thicker plates and better buckling resistance have a higher capabilit for dissipating. Therefore, the optimal steel shear panel damper models exhibit lesser frac Similarly, the stress distribution of selected models following a 5% drift ratio application is shown in Figure 9. As a result of the analysis, the maximum stress distribution was found to be uniformly distributed at about 350 Mpa at the end of the damper, similar to the distribution of maximum plastic strain. It is shown that the stress distribution of the optimal model is steadily demonstrated through the length of the link due to better materials implementation and higher energy dissipation. The study also showed that the geometrical properties of the optimal model are responsible for the mode of transition from flexure yielding to shear yielding, in accordance with the transitions equation provided by Farzampour and Eatherton [18,47]. In addition, the optimization analysis has shown that the models with thicker plates and better buckling resistance have a higher capability for dissipating. Therefore, the optimal steel shear panel damper models exhibit lesser fracture and improved energy dissipation capacity compared to other dampers. Based on the normalized maximum PEEQ results (obtained from values associate with the optimal solution), behavioral comparisons are performed. Table 1 shows the re sults of the dissipated energy and optimization function. The buckling resistance model exhibit effective energy dissipation due to the occurrence of both shear and flexure yield ing limit states simultaneously allowing uniform plastic strains and stresses distribution  Based on the normalized maximum PEEQ results (obtained from values associated with the optimal solution), behavioral comparisons are performed. Table 1 shows the results of the dissipated energy and optimization function. The buckling resistance models exhibit effective energy dissipation due to the occurrence of both shear and flexure yielding limit states simultaneously allowing uniform plastic strains and stresses distribution. Following the determination of the damper optimal form, the initial and the optimized dampers are investigated for energy dissipation and hysteretic behavior. The overall hysteretic curves of the optimized and initial models are illustrated in Figure 10. From the simulation results and hysteresis curves, it is shown that interior damper optimization possesses a lower yield force and significantly higher stiffness.
Buildings 2022, 12, x FOR PEER REVIEW 11 of 13 the simulation results and hysteresis curves, it is shown that interior damper optimization possesses a lower yield force and significantly higher stiffness. Figure 10. Hysteretic behavior comparison of the optimized dampers and initial models.

Conclusions
In this study, the optimization of the curved metal dampers for improving the energy dissipation capacity was investigated. Three types of combined low yield steel plate dampers were considered for simplifying the deficiencies (high yield point, and low adjustability) of the traditional low yield point plate steel damper. The low yield point steel plates and common steel plates with proportions were proposed. The study utilized the concept of "maximum stiffness" and "full stress state," for optimizing the dampers with interior hollows, boundary hollows, and ellipse hollows based on the alternating topology optimization technique. It is shown that the obtained three forms of combined low yield point steel plate dampers satisfied the optimization objectives.
An analysis of the FEM for various optimized combined steel plate dampers was computationally performed. The study has shown that the interior optimized damper is suitable for large deformation, and the energy dissipation capacity of the boundary optimized damper is relatively higher than conventional dampers. The study also used two types of optimization models having different material ratios for numerical simulation, and the results determine that it is possible to regulate the yield stress of the combined dampers by changing the material proportions.
In addition, over 300 butterfly-shaped links were generated following FEM validations. It should be noted that energy dissipation over concentrated PEEQ values is achieved based on the proposed optimization function. Using the WOA algorithms, it is possible to maximize energy dissipation capacity and minimize the plastic strain concentration over the links. The proposed optimization criteria show that optimal values of a geometry occur when it transforms flexural yielding to the shear yielding point. The proposed optimization procedure exhibits optimal dissipation of energy, with a minimum fracture. Finally, the optimal model demonstrates the development of the plastic strain concentration points at a far position from sharp angle geometrical sharpness resulting in desirable fracture resistance. It is noted that the proposed parametric study and optimization were developed for the specific shape of a commonly used damper; therefore, for extending and applying the work to other dampers types, further investigations are required.

Conclusions
In this study, the optimization of the curved metal dampers for improving the energy dissipation capacity was investigated. Three types of combined low yield steel plate dampers were considered for simplifying the deficiencies (high yield point, and low adjustability) of the traditional low yield point plate steel damper. The low yield point steel plates and common steel plates with proportions were proposed. The study utilized the concept of "maximum stiffness" and "full stress state," for optimizing the dampers with interior hollows, boundary hollows, and ellipse hollows based on the alternating topology optimization technique. It is shown that the obtained three forms of combined low yield point steel plate dampers satisfied the optimization objectives.
An analysis of the FEM for various optimized combined steel plate dampers was computationally performed. The study has shown that the interior optimized damper is suitable for large deformation, and the energy dissipation capacity of the boundary optimized damper is relatively higher than conventional dampers. The study also used two types of optimization models having different material ratios for numerical simulation, and the results determine that it is possible to regulate the yield stress of the combined dampers by changing the material proportions.
In addition, over 300 butterfly-shaped links were generated following FEM validations. It should be noted that energy dissipation over concentrated PEEQ values is achieved based on the proposed optimization function. Using the WOA algorithms, it is possible to maximize energy dissipation capacity and minimize the plastic strain concentration over the links. The proposed optimization criteria show that optimal values of a geometry occur when it transforms flexural yielding to the shear yielding point. The proposed optimization procedure exhibits optimal dissipation of energy, with a minimum fracture. Finally, the optimal model demonstrates the development of the plastic strain concentration points at a far position from sharp angle geometrical sharpness resulting in desirable fracture resistance. It is noted that the proposed parametric study and optimization were developed for the specific shape of a commonly used damper; therefore, for extending and applying the work to other dampers types, further investigations are required.