Efficient implicit constraint handling approaches for constrained optimization problems

Many real-world optimization problems, particularly engineering ones, involve constraints that make finding a feasible solution challenging. Numerous researchers have investigated this challenge for constrained single- and multi-objective optimization problems. In particular, this work extends the boundary update (BU) method proposed by Gandomi and Deb (Comput. Methods Appl. Mech. Eng. 363:112917, 2020) for the constrained optimization problem. BU is an implicit constraint handling technique that aims to cut the infeasible search space over iterations to find the feasible region faster. In doing so, the search space is twisted, which can make the optimization problem more challenging. In response, two switching mechanisms are implemented that transform the landscape along with the variables to the original problem when the feasible region is found. To achieve this objective, two thresholds, representing distinct switching methods, are taken into account. In the first approach, the optimization process transitions to a state without utilizing the BU approach when constraint violations reach zero. In the second method, the optimization process shifts to a BU method-free optimization phase when there is no further change observed in the objective space. To validate, benchmarks and engineering problems are considered to be solved with well-known evolutionary single- and multi-objective optimization algorithms. Herein, the proposed method is benchmarked using with and without BU approaches over the whole search process. The results show that the proposed method can significantly boost the solutions in both convergence speed and finding better solutions for constrained optimization problems.

For example, several existing CHTs in the literature are based on feasibility and infeasibility regions [4][5][6] , priority assignment 7 , and tournament selection and a selection operator 7,8 .Some studies also provide surveys of CHTs [8][9][10][11] .In contrast, implicit techniques do not necessitate the explicit form of constraints.They handle constraints inherently during the search for the optimal solution, avoiding the need for additional objective function evaluations.Implicit CHT simplifies the problem formulation, but it may face challenges in handling complex or nonlinear constraints effectively.
For instance, Diwekar and Rubin 12 proposed an implicit constraint handling technique for the ASPEN chemical process, which is based on a mixed integer mathematical programming approach consisting of master and subproblems.Since most of the relations in the above-mentioned problems are implicit 12 , proposed a method for partitioning the variables to significantly decrease computational time.
Raghavan et al. 13 developed an implicit constraint handling technique for optimization based on the proper orthogonal decomposition (POD) of shapes, further producing a bi-level reparameterization approach for structural geometries.Uemura et al. 14 proposed a real-coded genetic algorithm for implicitly constrained black-box optimization, in which the method employs the weighted mean of the best individuals in a population to find the optimal solution.
Mirabel and Lamiraux 15 proposed a method to handle all constraints explicitly and implicitly that deals with manipulation planning, where the constraints are solved explicitly as much as possible.Implicit handling is alternatively employed with few variables in case all constraints cannot be handled explicitly.Nomura et al. 16 introduced a natural evolution strategy called DX-NES-IC for implicit handling of constrained black-box optimization, which demonstrated better performance than DX-NES, xNES, and CMA-ES.
The current study extends the work of Gandomi and Deb 6 by extending and applying boundary update (BU) to single-objective and multi-objective optimization problems (MOOPs).As novel strategies, two switching mechanisms are suggested that transform the landscape and variables to the original problem when the feasible region is found.Then, the optimization process is continued without the BU method.
In the work by Gandomi and Deb 6 , the authors proposed a novel approach, boundary update (BU), as an implicit CHT that updates variable bounds by directly using the constraints and then applying them to several single-objective optimization problems.Since, the BU method is an implicit CHT, it should be coupled with an explicit CHT, that is feasibility rules in our study 17 .As for the strategy without BU, only an explicit CHT (feasibility rules) is applied to solve the problem.The BU method is an implicit constraint handling approach that aims to cut the infeasible search space over iterations to find the feasible region faster.Although the BU method directs the search operators to the feasible space and reaches the first feasible solution reasonably fast, it twists the search space, making the optimization problem more challenging.The current study tries to tackle the above-mentioned problem by proposing two switching approaches.
The current study augments existing research by introducing two novel switching mechanisms.In the initial method, called Hybrid-cvtol, the BU approach is employed until constraint violations reach zero, ensuring zero violations across the entire population.Subsequently, the algorithm transitions to the optimization process without utilizing the BU approach.The second switching mechanism, called Hybrid-ftol, in this study involves employing the BU method during the initial phase until the objective space remains unchanged for a specified number of generations.At that point, the optimization problem shifts to a state without utilizing the BU approach.
The remainder of the paper is organized as follows.Section "Proposed approach" presents the proposed approach.Section "Research methodology" discusses the research methodology.Section "Numerical examples" provides numerical examples.Finally, Sects."Discussion", "Limitation", and "Conclusion" give the discussion, limitation, and conclusion.

Proposed approach
A constrained optimization problem can be formulated in Eqs.(1-3) as follows: where F(x) is the objective vector that consists of several objectives (t is the number of objective functions); n and m are the number of inequality and equality constraints, respectively; and x is the decision variable bounded by the lower bound (LB) and upper bound (UB) vectors.These equations yield several Pareto optimal solutions rather than producing a single solution 18 .In the case of single-objective optimization, a single solution can be found as a solution to F(x) = f 1 (x) .In the optimization process, the boundaries of decision variables, represented by lower bounds (LB) and upper bounds (UB) vectors, undergo dynamic changes as the algorithm progresses through iterations.This dynamic nature enables the optimization process to adapt and explore a continually evolving solution space, illustrating the inherent variability and adaptability of decision variable boundaries throughout the iterative stages.Hence, the iterative nature of optimization involves the dynamic evolution of the i-th decision variable boundaries throughout each iteration.
The proposed method uses the constraints to narrow down variable space and then forces the algorithm to focus its search in the feasible region by limiting the viable search space for the variable(s).In the BU method, (1) Maximize(Minimize)F(x) = (f 1 (x), . . ., f t (x)) www.nature.com/scientificreports/ the boundaries change iteratively and are updated during the optimization procedure.Mathematically, it could be written as follows 6 : where l i,j and u i,j are dynamic lower bound and upper bound for the ith decision variable, respectively.Updating the bounds involves the following scenarios: Else: where ( lb u i , ub u i ) are updated boundaries.BU begins by selecting a repairing variable that can handle the most significant number of constraints without overlapping other repairing variables.As such, if there is more than one candidate, the one that handles the greatest number of constraints is selected.Also, if there is still another candidate for variable selection, one variable is chosen randomly 6 .
If a repairing variable handles the first k i constraints, then: If another repairing variable (x r ) is defined: In the BU method, a repairing variable could be substituted with a generalized semi-independent variable 6 and then rewritten with lower and upper bounds as: where p 1 , . . ., p h are semi-independent variables; and mx indicates the set of selected repairing variables.
After selecting repairing variables, the search operator will be applied to the problem, the boundaries of non-repairing variables will be checked, and the mx-vector will be updated.During the solution procedure, the boundaries of repairing variables will be updated, and then the semi-independent variables ( p i , i = 1, . . ., h) are remapped to the actual variables using updated boundaries.In the end, those constraints that were not involved in the repairing variable boundary, along with fitness values, will be evaluated using actual variables.Algorithm 1 illustrates the BU method in a constrained optimization problem.

While the criteria are not met
Applying the search algorithm Applying the BU method and updating the mx-vector Check the boundaries of non-repairing variables Update the boundaries of repairing variables Remap p variable to the original variables using updated boundaries ( , ) Evaluate the constraint violations Evaluating fitness function End while Algorithm 1. Implementation of the BU method in a constrained optimization problem.As an illustrative example, the following multi-objective optimization problem (MOOP) is considered to explain BU: Either x 1 orx 3 can be selected as a repairing variable of the single constraint.All the variables are in the range of − 2 to 2. From the variable selection strategy 6 , as explained earlier, x 3 is selected as a repairing variable.The inequality constraint is solved with respect to x 3 as the lower and upper bound functions.The repairing variable ( x 3 ) is substituted with the mapped variable, p , in the range of 0-1.Therefore, the lower and upper bounds of the variables are {− 2, − 2, 0} and {2, 2, 1}, respectively.When the search algorithm is applied to the problem, the repairing variable is substituted with the mapped variable, and the boundary of the repairing variable is mapped for the rest of the search cycle 6 .The original variable can be calculated by the following formula to compute the objectives and constraints: where ub u 3 andlb u 3 are the lower and upper bounds for the repairing variable x 3 , respectively.The updated bounds can be written as follows: The repairing variable is remapped to the actual boundary when the fitness values are evaluated.Figure 1a-d compare Pareto-optimal fronts found by the BU method coupled with a CHT (here, feasibility rules) and by the method without the BU method (only the feasibility rules approach is considered) for a population size of 100 for different generations.Figure 1 shows that the BU method produces more non-dominated solutions than the approach without BU.However, the number of non-dominated solutions is not always a good indicator for comparison.In the following sections, other indicators are explained in detail to compare the algorithms.
Although the BU method directs the search operator to the feasible space and reaches the feasible region in a reasonable time, it twists the search space, which can make the optimization problem more challenging and, hence, does not always bear good results and could potentially result in premature convergence.Consider the following example: In case of selecting variable x 2 as a repairing variable, the design space for the original variables x 1 and x 2 (Fig. 2a) is converted to the new design space (Fig. 2b).Twisting landscape may result in challenges for optimization problems, such as premature convergence.Some of these issues are presented in the following examples.www.nature.com/scientificreports/Therefore, two distinct hybrid methods are introduced, each incorporating its own unique switching mechanism.In the initial approach, the BU method is employed until all constraint violations are reduced to zero, ensuring the absence of violations across the entire population.Subsequently, the algorithm transitions to the optimization process, excluding the use of the BU approach.In the second switching mechanism, the BU method is utilized in the initial phase until the objective space remains unchanged for a specified number of generations.Following this, the optimization problem transitions to a state where the BU approach is no longer employed.In both approaches, the optimization process is halted once the criteria and threshold tolerance are met.In the subsequent phase, the optimization problem is solved without the BU method, using the final population from the previous phase as the initial population for the new optimization problem. Figure 3 illustrates that the handling of boundaries and constraints occurs after the application of search operators.
During the first phase, the proposed methods simultaneously check constraints and examine variable boundaries, making the optimization process an iterative search procedure for handling constraints and boundaries concurrently.If the criteria are met, the second phase of the optimization process is conducted without utilizing the BU approach, constituting a sequential iterative search procedure where constraint handling occurs after boundary checking.The implementation of the proposed hybrid method is detailed in Algorithm 2.

Research methodology
To evaluate the proposed method, Python and Pymoo libraries 19 were used, which include well-stabilized algorithms.The specific algorithm parameters are presented in Table 1.Furthermore, because the implemented algorithms act stochastically, each of these algorithms were ran 31 times with random initial point(s).
Several examples, including benchmarks and real-world problems, were solved using three methods (without the BU method, with the BU method, and hybrid method).The presented examples are of different types: constrained by one, two, and three objective functions.For the problem with one objective function, the proposed method was implemented using evolutionary strategy (ES), differential evolution (DE), stochastic ranking Table 1.Specific parameter settings for the algorithms used in this study.www.nature.com/scientificreports/evolutionary strategy (SRES) 20 , and biased random key genetic algorithm (BRKGA) 21 .For two objective functions, the BU method was applied using NSGA-II 17,22 .For three and five objective functions, some other EAs were used.As such, NSGA-III 23 , U-NSGA-III 24 , AGE-MOEA 25 , NSDE 26 , and NSDER 27 were applied to optimization problems.The following metrics were evaluated to assess the accuracy of the obtained results.
• Number of non-dominated solutions: The BU method was coupled with an explicit constraint handling technique and compared to the approach without BU (only feasibility rules were applied to the optimization problem).Since the BU method aims to reduce the search space, the number of non-dominated solutions was determined for comparison.• Constraint violation (CV): Compared to unconstrained optimization problems, constrained optimization problems are more challenging since a large proportion of infeasible regions appears in the search space (meaning the hit ratio is low), especially for highly constrained problems where these regions lead to some difficulties, such as feasibility-hardness and convergence-hardness. Therefore, even finding one feasible solution is a significant achievement.To test the effectiveness of the proposed method, the first feasible solutions found by the methods were compared and the whole population was tracked against generations.• Performance indicators: Generational Distance (GD) 28 , Generational Distance Plus (GD+) 29 , Inverted Gen- erational Distance (IGD) 30 , Inverted Generational Distance Plus (IGD+) 29 , and Hypervolume (HV) are the indicators that were analyzed.• Running metric 31 : It is possible to trace the difference in the objective space followed by generations.
Moreover, the population was fixed at 100, and the maximum number of experiments was set to 11.

Numerical examples
The proposed approach was applied to several single-objective and multi-objective optimization examples, including benchmarks and real-world problems.The presented examples include problems constrained by one, two, and three objective functions.The performance of the three strategies was evaluated using an example that has linear constraints.Following this numerical example, the proposed approach was evaluated to solve a surrogate model for a car-side impact design problem.This design problem consists of two versions of optimization, i.e., single-and multi-objective optimizations, and is an example of a black-box optimization problem.The speed reducer problem, a single-objective constrained engineering optimization problem, was also considered.For the multi-objective optimization problem, two benchmarks, namely OSY and BNH, and three real-world constrained optimization problems, including welded beam design, Cantilevered Beam design, and multi-objective car-side design problems were solved via the three strategies.
Several metaheuristic methods have been developed and implemented to reduce the total runtime of the optimization problems 2,32 , but still produce results that are almost as accurate and precise as conventional solving methods.Also, real-world engineering problems involve some type of optimization that is often constrained, most of which are considered MOOPs.No single solution exists for a MOOP; instead, different solutions generate trade-offs for various objectives.Furthermore, MOOPs arise naturally in most fields, and solving them has been a challenging problem for researchers 32,33 .EA methods have been identified as more effective in tackling the challenges that arise from MOOPs, for which the form of the Pareto-optimal front (discontinuity, nonconvexity, etc.) is not important 34,35 .Moreover, most multiobjective evolutionary algorithms (MOEAs) use the dominance concept 36,37 .In this work, for the single-objective optimization problems, Genetic Algorithm (GA), Differential Evolution (DE) 38 , Evolutionary Strategy (ES) 39 , Stochastic Ranking Evolutionary Strategy (SRES) 20 , and Biased Ranking Key Genetic Algorithm (BRKGA) 40 were considered.NSGA-II 41 is applied to bi-objective and multi-objective optimization problems with three objective functions and to multi-objective optimization problems with more than three objective functions, which are called many-objective optimization problems.NSGA-III 23 , U-NSGA-III 24 , and Adaptive Geometry Estimation-based Multi-objective Evolutionary Algorithm (AGE-MOEA) 25 were further implemented for optimizing the problems.For all problems, the BU method was coupled with feasibility rules proposed by Deb 41 .Table 2 presents the studied examples in this work to assess the effectiveness of the proposed method in comparison with the other methods.

Single optimization problems
In this section, three single-optimization problems, namely G1, car-side design problem, and speed reducer, are considered to be solved using the different strategies discussed above.

G1 problem
The G1 problem is a famous mathematical constrained optimization problem that includes 13 decision variables, a quadratic objective function, and nine linear constraints.The problem is an excellent example of a highly linearly constrained problem, with a feasibility ratio of 0.111% 6 (see the Appendix).In other words, it is a highly challenging to find a feasible solution for this constrained problem.The three different strategies were applied to this problem.According to the variable selection strategy, x 10 , x 11 , and x 12 were selected as repairing variables for handling all nine constraints.The global optimum for this problem is -15, and the optimization results by using GA, DE, ES, SRES, and BRKGA with the BU method, without the BU method, and hybrid methods are presented in Tables 3, 4, 5 and 6 and Fig. 4. The Friedman rank test at a significant level of 5% was used to detect differences in treatments across test attempts (Table 4).Also, the Post-hoc Nemenyi Friedman test was performed to identify exactly which groups have different means (Tables 5, 6).
The comprehensive analysis of algorithm performance on the G01 optimization problem, characterized by 13 continuous decision variables, nine linear constraints, and a quadratic objective function, unveils intricate relationships between problem features and algorithmic behavior.Each algorithm's response to the problem's unique attributes provides valuable insights into their adaptability and efficacy.
The hybrid methods, consistently showcasing superior performance across diverse algorithms, reflect an inherent versatility that aligns with the multi-faceted nature of the G01 problem.The distinct advantage of the hybridization approach becomes particularly evident in the handling of constraints.The BU method, serving Table 3. Statistical results of different methods on G1 problem (The best performing method (objective function values) is marked in bold).a "None" means no BU method or hybrid is used and only feasibility rules as an explicit constraint handling method is worked alone.www.nature.com/scientificreports/as an implicit constraint-handling technique, demonstrates a keen ability to navigate the infeasible search space effectively.This characteristic proves invaluable in problems like G01, where constraints play a crucial role in shaping the feasible region.The hybrid methods, by strategically integrating the BU method, adeptly leverage this constraint-handling mechanism, enabling faster convergence to the feasible area.

Algorithm
Considering the explicit constraints and the quadratic objective function in G01, it becomes apparent that the algorithmic adaptation to these features significantly impacts performance.Genetic algorithms (GA), when combined with the BU method (as hybrid methods), exhibit a noteworthy convergence pattern (as illustrated in Fig. 5).The initial phase, focused on constraints handling, ensures the algorithm converges to the feasible area efficiently.Subsequently, the transition to the optimization process without BU capitalizes on the algorithm's newfound knowledge of the feasible region, resulting in a streamlined path to the final solution.This adaptive approach to the problem's structure showcases how algorithmic behavior is intricately linked to specific features, allowing for a dynamic and context-aware optimization process.
The inadequacy of the SRES algorithm in this context sheds light on the sensitivity of algorithmic performance to the problem's characteristics.The SRES algorithm, designed for stochastic optimization, might face challenges when confronted with the structured nature of the G01 problem.The lack of significant improvement suggests that algorithmic suitability is not universal but rather contingent on the problem's intricacies.This observation underscores the importance of aligning algorithmic choices with the specific features and constraints of the optimization problem at hand.Furthermore, the reduction in standard deviation observed in the hybrid methods implies a more stable and reliable convergence behavior.This heightened stability is particularly valuable in the presence of complex problem features, such as the quadratic objective function, where rapid fluctuations in convergence patterns might hinder progress.The hybrid methods' ability to maintain a consistent trajectory towards the optimal solution underscores their resilience in handling the intricacies introduced by specific problem characteristics.
In summary, the algorithmic performance on the G01 problem intricately relates to its unique features.The hybrid methods showcase adaptability and versatility in addressing constraints, demonstrating a keen understanding of the problem structure.The nuanced convergence patterns, strategic switching mechanisms, and sensitivity to problem-specific attributes provide valuable insights into the dynamic interplay between algorithmic behavior and the intricacies of the optimization problem.

Car-side problem
The car-side impact design problem, a notorious challenge in engineering, stands out as a complex and timeconsuming optimization task, as depicted in Fig. 6.This problem, adaptable to both single-objective and multiobjective formulations, adds layers of intricacy due to its quadratic objective function and explicit models.Specifically, when addressing the single-objective version with the goal of minimizing car weight, the constraints and objective function take on quadratic forms, amplifying the intricacies involved.The thicknesses of the inner Table 5. Summary of the p-value of the posthoc_nemenyi_friedman over all runs.Posthoc_nemenyi_ friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.Significant values are in bold.www.nature.com/scientificreports/floor side ( X 3 ) are strategically selected as decision variables to solve the explicit models, and their bounds are dynamically adjusted to accommodate the interconnected variables, adding an additional layer of complexity.The results of this formidable optimization challenge are presented in Tables 7, 8, 9 and 10 and Figure 7, offering a quantitative assessment of algorithmic performance.Statistical tests (Tables 8, 9, 10) and visualizations in Figure 7 highlight the superiority of the hybrid methods, particularly when combined with the genetic algorithm (GA).This dominance is established through the comparative analysis, showcasing that the hybrid approach consistently outperforms both the individual methods and, notably, requires fewer function evaluations (FE) to attain the final solution.
Significantly, the examination of the results reveals nuanced differences between the differential evolution (DE) and evolutionary strategy (ES) methods, with the hybrid method showcasing slightly superior outcomes.The incorporation of the BU method, along with the hybrid methods, demonstrates a tangible advantage, resulting in enhanced solution quality.The visual representation of the switching point for GA in Figure 8 provides further clarity on the optimization process.The algorithm strategically transitions from the BU-assisted phase to the standard optimization phase, showcasing an adaptive mechanism designed to optimize the convergence www.nature.com/scientificreports/path.However, the study identifies limitations within the context of the car-side impact design problem.While GA, BU, and hybrid methods exhibit notable improvements, the performance of biased random key genetic algorithm (BRKGA) and stochastic ranking evolution strategy (SRES) fails to show enhancement when incorporated into the hybrid framework.This suggests that BRKGA and SRES may not be suitable for addressing the intricacies of this particular problem, underscoring the need for algorithmic adaptability in the face of diverse engineering challenges.
In conclusion, the car-side impact design problem proves to be a formidable optimization challenge, and the hybrid methods, particularly in conjunction with GA, emerge as the most effective approach.The adaptability of the hybrid methods, strategic switching mechanisms, and the incorporation of the BU method collectively contribute to enhanced solution quality and computational efficiency.Nonetheless, the study sheds light on the context-specific nature of algorithmic suitability, emphasizing the importance of tailoring optimization strategies to the unique characteristics of complex engineering problems.

Speed reducer
The examination of the speed reducer design problem, featuring one objective function and seven continuous decision variables, has yielded multifaceted insights into the interplay between algorithmic approaches and problem-specific characteristics.With the overarching goal of minimizing the total weight of the speed reducer, the complexity of this task is compounded by a combination of linear and nonlinear constraints.Strategically leveraging three repairing variables-specifically, the face width and shaft diameters ( x 1 , x 6 , and x 7 )-proves instrumental in managing these constraints effectively.
The results, meticulously presented in Tables 11, 12, 13 and 14 and Fig. 9, reveal a consistent trend: the hybrid approaches consistently outperform standalone methods, yielding the best objective function values across all algorithms considered.This dominance suggests that the hybridization process contributes to algorithmic robustness, offering a versatile strategy to navigate the intricate optimization landscape presented by the speed reducer design problem.
Notably, the algorithms (GA, DE, and BRKGA) exhibit improved performance when integrated into hybrid frameworks, signifying a synergistic effect where the strengths of individual algorithms complement one another.The statistical tests conducted (Tables 12, 13, 14) affirm the significance of these improvements, providing a robust validation of the observed enhancements and instilling confidence in the findings.A closer examination of GA within the hybrid methods reveals its particular suitability for this problem, consistently outperforming standalone and BU-assisted methods.This positions GA as a robust choice for addressing the complexities inherent in the optimization landscape of the speed reducer design problem.
The incorporation of the BU method consistently leads to a substantial reduction in the number of function evaluations required to reach the final solution.This reduction in computational cost is a valuable contribution,  Despite the positive outcomes, the study highlights that algorithms like BRKGA and SRES do not exhibit improved performance within the hybrid framework, suggesting their unsuitability for this particular problem.This observation underscores the importance of algorithmic adaptability and tailoring strategies to specific problem characteristics.
In conclusion, the insights gained from the speed reducer design problem analysis underscore the effectiveness of hybrid approaches in enhancing algorithmic performance.The adaptability of GA, efficiency gains with the BU method, and the nuanced trade-offs between algorithms offer valuable guidance for tackling real-world engineering optimization challenges.These findings contribute not only to the understanding of algorithmic behavior in this specific context but also pave the way for broader applications in engineering design optimization.

Multi-objective optimization problems
The following subsections provide two benchmarks and two real-world optimization problems to assess the effectiveness of the proposed hybrid method for multi-objective optimization problems.

Osyczka and Kundu (OSY)
Osyczka and Kundu 44 introduced the OSY test problem, a widely recognized benchmark featuring two nonlinear objective functions and six linear and nonlinear constraints.As highlighted earlier, the BU method is a central focus of this work, particularly in its ability to directly handle constraints and its subsequent integration with other explicit constraint-handling techniques.The BU approach strategically reduces the feasible search space by altering one or more variables, compelling the algorithm to concentrate on exploring regions that adhere to the defined constraints.For the OSY problem, a specific repairing variable selection strategy is employed, targeting variables x 1 , x 4 , and x 6 identified as effective for handling all constraints inherent in the problem.
Following the application of the BU approach, NSGA-II is employed to navigate the complex landscape of the OSY problem.NSGA-II, known for identifying viable Pareto fronts, initially achieves an entirely feasible population after 1500 evaluations.In contrast, the coupling of the BU method with the algorithm accelerates this process, achieving a fully feasible population after only 300 evaluations.Subsequently, the optimization process seamlessly transitions to routine optimization without the continued use of the BU method.
The evaluation of performance indicators, presented in Tables 15, 16, 17 and 18 and Fig. 10, offers a comprehensive understanding of the algorithmic performance.It becomes evident that the hybrid methods, particularly when incorporating the BU approach, consistently outperform other methods.Statistical analyses, including the Friedman rank and posthoc tests, validate this observation by demonstrating a significant improvement in all performance indicators (excluding IGD+) when the hybrid approach is implemented.
Examining Fig. 10 through the lens of a violin plot further illustrates the superiority of the hybrid methods.The best values for all performance indicators, with the exception of IGD, show marked improvement when the BU and hybrid methods are employed.This enhancement in performance underscores the effectiveness of the www.nature.com/scientificreports/BU method in combination with the hybrid approach, providing valuable insights into the algorithmic strategies that lead to improved convergence and solution quality for the OSY problem.In summary, the OSY problem, coupled with the BU method and hybrid approach, serves as a compelling case study.The specific variable selection strategy, the seamless transition between BU-assisted and routine optimization, and the consistent improvement in performance indicators highlight the adaptability and effectiveness of the proposed approach.These findings contribute not only to the optimization of the OSY problem but also offer broader insights into the potential of the BU method and hybridization in addressing complex, constrained optimization scenarios.

Bin and Korn Test Problem (BNH)
BNH is a test problem proposed by Binh and Korn 5 , characterized by two conflicting objectives and two constraints, introduces a challenging landscape for optimization.The nonlinear nature of the objectives and constraints adds an additional layer of complexity, requiring sophisticated algorithmic strategies to navigate the  In addressing this nonlinear bi-objective optimization problem, the algorithmic focus will be on developing robust strategies that balance the optimization of both objectives while adhering to the imposed constraints.The inherent trade-offs between conflicting objectives necessitate a careful exploration of the Pareto front to identify solutions that represent optimal compromises.This problem involves two continuous variables, x 1 and x 2 , which can be strategically chosen as repairing variables.The findings presented in Tables 19, 20, 21, 22, 23, 24 and 25 and Fig. 11 illuminate a substantial disparity between the hybrid method and alternative approaches.Particularly noteworthy is the superior performance of the hybrid method, especially when objective function tolerance is considered, showcasing noteworthy enhancements over other methods across various cases.The results unequivocally demonstrate the hybrid approaches' consistent outperformance in tackling this bi-objective optimization problem characterized by nonlinear constraints and objective functions.Furthermore, the performance indicator values associated with the hybrid approaches are significantly smaller compared to their counterparts, underscoring the effectiveness of the proposed methodology in achieving superior convergence and solution quality.These insights not only contribute to the optimization of the BNH problem but also offer valuable considerations for addressing the nuances of nonlinear bi-objective optimization challenges across diverse applications.

Welded beam design problem
The welded beam fabrication is a well-known engineering optimization benchmark problem 48 that examines the trade-off between the strength and cost of a beam.The problem includes four constraints (Eqs.32-35) and four design variables, which are: • x 1 = h, the thickness of the welds  where t max = 13600, s max = 30000, P = 6000, P c = 64746.022 3 ) .The two objective functions are to minimize the fabrication cost of the beam and minimize the deflection of the end of the beam: The four variables of this problem include: (1) thickness of the weld, (2) welded joint's length, (3) beam's width, and (4) beam's thickness.Using the three strategies, NSGA-II was applied to this problem, and fitness values were evaluated.
The exploration of sampling and crossover operators in EAs stands as a critical endeavor, influencing the convergence speed and efficacy of these algorithms, particularly in the context of noisy optimization.In the context of the Welded Beam Design problem, where the impact of noise is significant, the experimentation involved (33)   www.nature.com/scientificreports/diverse operators, including real random sampling, Latin hypercube sampling, and two distinct crossovers: real-SBX crossover and uniform crossover.Four performance indicators, namely GD, GD+ , IGD, and IGD+, were meticulously evaluated against these operators to glean insights into their comparative effectiveness.The compelling findings, detailed in Tables 26, 27, 28, 29, 30, 31 and 32 and illustrated in Fig. 12, underscore the superiority of the hybrid approach over the two alternative methods.
Examining the Pareto optimal solutions generated by the hybrid methods across different generations reveals a consistent outperformance, particularly noteworthy when considering constraint violations as tolerance with real random sampling.This suggests that the hybrid approach, when coupled with appropriate operators, significantly improves the optimization process for the Welded Beam Design problem.Notably, the hybrid method's effectiveness is evident with Latin hypercube sampling, showcasing its compatibility and superior performance within this optimization context.Furthermore, the hybrid method incorporating objective function tolerance using uniform crossover emerges as a robust option, outperforming optimization without the BU approach.
In contrast to initial expectations regarding the BU approach as a constraint handling technique, the results indicate that it may not be the optimal choice for this specific nonlinear bi-objective model.However, the hybrid approaches, through strategic integration of various operators, demonstrate a substantial enhancement in performance.This nuanced observation highlights the importance of selecting appropriate operators within hybrid approaches, suggesting that the effectiveness of the BU method can be context-dependent.
The statistical analyses, including the Friedman rank test and the Posthoc_Nemenyi_Friedman, contribute to the robustness of the findings.These tests not only affirm the significant differences between the hybrid approach www.nature.com/scientificreports/and the two other methods but also reinforce the superiority of the proposed hybrid method for solving the Welded Beam Design problem.The identification of the switching point for the second phase of the optimization process at evaluation 300 further exemplifies the efficiency of the hybrid approach compared to the without BU approach, where the feasible region is reached at evaluation 1600 (Fig. 13).
In conclusion, the investigation into sampling and crossover operators, coupled with the nuanced insights from performance indicators and statistical analyses, positions the hybrid approach as the preferred strategy for solving the Welded Beam Design problem.The adaptability of the hybrid method, evident in its performance across various operators, reinforces its efficacy in addressing the challenges posed by the specific characteristics of this nonlinear bi-objective optimization scenario.30.Summary of the p-value of the posthoc_nemenyi_friedman test evaluated over all runs (sampling-Continue). a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.

31.
Summary of p value of the posthoc_nemenyi_friedmantest evaluated over all runs (crossover).a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.

Cantilevered stepped beam design problem
The cantilevered stepped design problem is an example of large-scale size problem.In the cantilevered beam design problem, the stepped cantilever beam must be able to carry a prescribed end load 46 .This problem is originally a single-objective optimization problem considering minimizing the beam volume; however, this study has extended the original problem and added one more objective function, which minimize end deflection subject to various engineering design variables (see Appendix).The beam supports the given load, P, at a fixed Table 32.Summary of p value of the posthoc_nemenyi_friedmantest evaluated over all runs (crossover-Continue). a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.

Hybrid-ftol vs. Without BU Hybrid-ftol vs. Hybridcvtol
Hybrid-ftol vs. BU  www.nature.com/scientificreports/distance L from the support.Besides, designers of the beam can vary the width ( b i ) and height ( h i ) of each sec- tion (Fig. 14).The problem is highly constraint and large-dimensional; however, for this example five segments (N = 5) including 10 dimensions (10 constraints) are considered.In the pursuit of solving both the single-and multi-objective versions of this intricate problem, repairing variables x 2 , x 4 , x 6, x 8 , and x 10 have been identified.The application of the BU method to this constrained nonlinear optimization problem is showcased, providing a valuable demonstration of the BU method's utility.Figure 15a-d visually presents the optimal solutions discovered by NSGA-II employing different methods across various generations.

Hybrid-ftol vs. Without BU Hybrid-ftol vs. Hybridcvtol p-value h a p-value h p-value h p-value h p-value h p-value
To assess the performance of the proposed approaches, the hypervolume (HV) indicator is utilized, offering a comprehensive measure of the quality of solutions (Tables 33, 34).The results, detailed in Tables 35, 36, 37, 38 and 39 and visualized in Fig. 16, unequivocally demonstrate the efficacy of the hybrid approaches in enhancing solution quality.This improvement is particularly noteworthy in the case of Latin hypercube sampling and simulated binary crossover (SBX), where the hybrid approaches exhibit a significant enhancement in solution quality.
This example provides valuable insights into the application of the BU method in addressing complex, largescale, constrained optimization problems.The extension of the original single-objective problem to a multiobjective setting showcases the flexibility of the proposed methodology.The utilization of hypervolume as a performance indicator reaffirms the substantial improvements achieved by the hybrid approaches.Overall, this study contributes not only to the optimization of the Cantilevered Stepped Design problem but also sheds light on the broader potential of the BU method in tackling diverse engineering optimization challenges.

Multi-objective car-side impact problem formulation
In scenarios where the black-box constraints of an optimization problem exhibit highly complex behavior, the direct evaluation of these constraints can be prohibitively time-consuming.In such instances, the use of surrogates becomes a valuable strategy to mitigate the complexity and expedite the optimization process.Surrogates,   www.nature.com/scientificreports/efficiency challenges, particularly in the context of a three-objective optimization problem, the hybrid approach proves adept at overcoming these difficulties.Statistical analyses of the results reveal that, in some cases and with different algorithms, the BU method may not be as efficient, yet the hybrid approach consistently outperforms it, signifying the added benefits of the hybridization strategy.
Comparative statistical tests emphasize the significant advantage of the hybrid approaches over the BU method.Notably, although some statistical tests indicate no significant difference between the hybrid and BU methods, the hybrid method consistently achieves higher hypervolume (HV) values.This suggests that, even when statistical significance is not apparent, the hybrid approach provides a qualitative improvement in solution quality.
In conclusion, the exploration of the car side impact problem demonstrates the efficacy of hybrid approaches in addressing complex, multi-objective optimization challenges with intricate black-box constraints.The incorporation of surrogates and the application of the BU method within a hybrid framework contribute to improved diversity and solution quality.These findings not only advance the optimization of the car side impact problem but also offer valuable insights into the broader applicability of hybrid evolutionary algorithms in addressing real-world, complex optimization scenarios.

Discussion
The BU method is an implicit constraint-handling technique coupled with an explicit constraint-handling technique that aims to cut the infeasible search space over iterations to find the feasible region faster.Since optimization is an iterative process whose boundaries can be changed during the process, the i-th decision variable boundaries in each iteration exhibit a dynamic nature.By applying this philosophy to boost the optimization process in the BU method, the bounds of selected variable (s) are repaired (repairing variables) to satisfy the constraints.As a result, new solutions are created within the updated boundary of the latest search space.The above-mentioned process helps the algorithm search for a solution within the feasible region instead of the whole area.Although the BU method directs the search operator to the feasible space and reaches the feasible region in a reasonable time, it twists the search space, making the optimization problem more challenging.As such, the results are not always accurate, and the algorithm may result in premature convergence.Therefore, in this paper, two "hybrid" switching mechanisms are proposed.In the proposed approach, the BU method is applied to the problem after initializing the algorithm, and in this phase, a threshold tolerance is set.The BU method Table 40.Results of different methods on the multi-objective car-side impact problem (HV values-sampling)."None" means no BU method or hybrid is used and only feasibility rules as an explicit constraint handling method is worked alone.Significant values are in bold.• By boosting the optimization process, the hybrid method finds the first and whole populations much faster than the classic BU and without BU methods for most problems (i.e., the single-or multi-objective optimization problems of this paper).Moreover, regarding the car side impact design problem, the hybrid method obtained good diversity, whereas without the BU method searched for solutions in a specific feasible space.• For some case studies, although the hybrid method resulted in the same (or close-to-final) solution as the two other methods, it requires fewer FEs to converge to a solution.• Compared to unconstrained optimization problems, constrained optimization problems are more challenging since many infeasible regions appear in the search space (meaning the hit ratio is low).This makes solving a constrained problem more challenging, especially highly constrained problems that can lead to difficulties such as convergence-, diversity-, and feasibility-hardness. Therefore, finding just one feasible solution is a significant achievement.The primary goal of the BU method is to minimize the search space by narrowing down the variable range and eliminating infeasible options.To compare the effectiveness of the proposed approach, the first feasible solution found by different methods is compared to one another, and the whole population is tracked against generations.The hybrid method uses this philosophy to switch the optimization process and avoid reaching premature convergence.On the other hand, once the whole feasible area is found by the BU method, the final solution (global) could be achieved with fewer FEs in a reasonable time.• Regarding the multi-objective car side impact design problem, the hybrid approach works significantly bet- ter than without the BU method when incorporated in the employed multi-objective algorithms.It is worth mentioning that although there is no significant difference in the final HV values between the hybrid and BU methods, as mentioned earlier, the hybrid method could achieve the final solution in fewer FEs.• In scenarios where explicit formulation of constraints is challenging, the use of surrogate models, also known as metamodels, proves to be a valuable solution.Surrogate models have demonstrated effectiveness in situations where deriving explicit formulations for constraints is impractical.These models serve as approximations of intricate simulations, providing an efficient means to evaluate constraints during optimization without the computational overhead of running resource-intensive simulations.For instance, consider optimization problems associated with finite element models, such as those involving nodal displacements and bar stresses in truss optimization.In these cases, surrogate models can be strategically trained to emulate the intricate behavior observed in simulation results.This approach becomes particularly advantageous when the direct formulation of these constraints proves to be impractical due to their complexity or reliance on computationally expensive procedures.
In EAs, sampling and crossover are critical factors that can strongly affect convergence speed and are often used in noisy optimization problems where reducing the negative impact of noise is crucial.Moreover, sampling is a popular strategy for dealing with noise.Therefore, experiments with different sampling operators, including real random sampling, Latin hypercube sampling, and two different real-SBX and uniform crossovers, were conducted to reduce the noise impact.

Limitations
Although this study provides valuable insights into optimization problems, there are certain limitations to this work.One significant limitation of this work lies in the relatively small size of variables in the problems.It should be noted that this study aims to contribute by addressing a specific aspect of the optimization problem.In this study, the benchmark libraries have not been used and some specific problems with solvable constraints with respect to at least one variable have been used.While the limited number of variables in our optimization problems may restrict the immediate applicability of our findings to large-scale size problems, it is clear the insights gained from this research serve as a valuable foundation for further research in the field.
Another limitation of this study is related to the experimental section, where the algorithms compared with the proposed algorithm are not the most recent.It is important to note that the aim of this work is not focused on proposing and developing state-of-the-art algorithms; rather, the advancement of constraint handling techniques Table 46.Summary of the p-value of the posthoc_nemenyi_friedman test evaluated over all runs (Crossover).a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.www.nature.com/scientificreports/and the development of the BU approach are being concentrated on.Consequently, well-known algorithms such as NSGA-II in the analysis have been included.However, as a future study, it would be valuable to consider incorporating state-of-the-art algorithms to provide a more comprehensive evaluation of our proposed approach in comparison to the latest advancements in the field.Moreover, while the current selection of repairing variables is designed to be applicable even in scenarios with multiple disjoint feasible regions, the evaluation of constraints can pose challenges in situations where the behavior is non-trivial.In such cases, alternative approaches to repairing variable selection, involving bookkeeping or linking of these disjoint feasible ranges, may be necessary.This intriguing aspect has been highlighted and it is proposed as a subject for future investigation, intending to delve deeper into the nuanced evaluation of constraints in conditions assessment.

Conclusion
This study presents an extension of the work proposed by 6 for solving constrained optimization problems efficiently.The proposed methods use an implicit constraint handling approach called boundary updating (BU), which focuses on directly handling constraints and is coupled with other explicit constraint-handling techniques (here, feasibility rules).
In this study, two augmented versions of the BU method were applied to single-objective, multi-objective, and many-objective constrained optimization problems, including benchmarks and real-world problems.To accomplish this goal, two thresholds are considered, each representing different switching method.In the initial method, the optimization process shifts to a state that does not involve the BU approach once constraint violations reach zero.In the second approach, the optimization process enters a phase without BU methods when there is no longer any observed change in the objective space.For single-objective constrained optimization problems, popular algorithms like GA, DE, ES, SRES, and BRKGA were employed, and NSGA-II, a renowned algorithm, was used for bi-objective problems.For the multi-objective constrained optimization problems, NSGA-III, UNSGA-III, and AGEMOEA were applied.The results show that the hybrid approach performs significantly better than with-and without the BU methods, and can find better solutions with fewer iterations.
As a future study, it is suggested that the proposed BU method be coupled with other constraint-handling techniques, for example different types of penalty methods with different objective functions.Also, additional test problems in the literature, such as many-objective and dynamic functions, could be considered for further validating the proposed hybrid approach.After BU generalization, it can also be applied to benchmark libraries such as CEC, NEVERGRAD, BBOB, GNBG, and DIRECTGOLib.Although BU is not defined for any specific optimization algorithms, implementing BU to a wider range of optimization algorithms can be investigated in another future study.Moreover, as a prospective avenue for future research, exploring the applicability of the proposed approach to constrained problems featuring multiple disconnected regions is a promising direction.Many constrained optimization scenarios exhibit this characteristic, with only certain regions potentially close to the true optimum.Investigating the effectiveness and adaptability of our approach in addressing such complex problem structures holds significant potential for enhancing the robustness and versatility of the proposed methodology.This important consideration will be a key focus in our forthcoming studies.

4 Figure 1 .
Figure 1.Pareto-optimal fronts obtained with and without the BU method.

Figure 2 .
Figure 2. Contour plot of the search space.

Figure 3 .
Figure 3. Optimization process of the proposed method.

Figure 4 .
Figure 4. Violin plot of implementation of using different algorithms constraint handling methods (G01 problem).

Figure 5 .
Figure 5. Constraint violation implemented by GA on the G01 test problem.

Figure 6 .
Figure 6.Schematic of car side impact design problem 42 .

Figure 9 .
Figure 9. Violin plot of implementation of using different algorithms constraint handling methods (speed reducer problem).

Figure 10 .
Figure 10.Violin plot of implementation of using different indicators (OSY problem).
www.nature.com/scientificreports/• x 2 = l, the length of the welds • x 3 = t, the height of the beam • x 4 = b, the width of the beam The constraints are:

Figure 11 .
Figure 11.Violin plot of implementation of using different indicators (BNH problem).

Figure 13 .
Figure 13.Constraint violation evaluated by NSGA-II (convergence of the whole population to feasible solutions).

Figure 14 .
Figure 14.Schematic of the stepped beam design problem.

End while Algorithm 2.
Implementation of the proposed method (hybrid methods).

Table 2 .
Test problems considered in this study.

Table 6 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs (continue).Significant values are in bold.

Table 7 .
Statistical results of different methods on car side design impact problem.a "None" means no BU method or hybrid is used and only feasibility rules as an explicit constraint handling method is worked alone.Significant values are in bold.

Table 8 .
Summary of p value of the Friedman rank test over all runs.

Table 9 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs.Significant values are in bold.

Table 10 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs (continue).Significant values are in bold.

Table 11 .
Statistical results of different methods on speed reducer problem.a "None" means no BU method or hybrid is used and only feasibility rules as an explicit constraint handling method is worked alone.Significant values are in bold.

Table 12 .
Summary of the p-value of the Friedman rank test over all runs.

Table 13 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs.Significant values are in bold.

Table 14 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs (continue).

Table 15 .
Statistical results of different methods on the OSY problem.a "None" means no BU method or hybrid is used and only feasibility rules as an explicit constraint handling method is worked alone.Significant values are in bold.

Table 16 .
Summary of the p-value of the Friedman rank test over all runs.

Table 17 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs.

Table 19 .
Statistical results of different methods on the BNH problem.Significant values are in bold.

Table 20 .
Summary of the p-value of the Friedman rank test over all runs.

Table 21 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs.
a Posthoc_nemenyi_ friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.

Table 22 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs-IGD value comparison.

Table 23 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs-IGD+value comparison.

Table 24 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs-GD value comparison.

Table 25 .
Summary of the p-value of the posthoc_nemenyi_friedman over all runs-GD+value comparison.

Table 26 .
Statistical results of different methods on the welded beam design problem (sampling).a "None" means no BU method or hybrid is used and only feasibility rules as an explicit constraint handling method is worked alone.Significant values are in bold.

Table 27 .
Statistical results of different methods on the welded beam design problem (crossover).Significant values are in bold.

Table 28 .
Summary of the p-value of the Friedman rank test over all runs.

Table 29 .
Summary of the p-value of the posthoc_nemenyi_friedman test evaluated over all runs (sampling).
a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.

Table 35 .
Summary of the p-value of the Friedman rank test over all runs.

Table 36 .
Summary of the p-value of the posthoc_nemenyi_friedman test evaluated over all runs (sampling).a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.

Table 37 .
Summary of the p-value of the posthoc_nemenyi_friedman test evaluated over all runs (sampling-Continue). a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.

Table 38 .
Summary of p value of the posthoc_nemenyi_friedmantest evaluated over all runs (crossover).a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.
, such as crossover and sampling, better results by choosing the appropriate operator setting may be achieved.•Regarding the multi-objective optimization problems, the hybrid method could find more non-dominated solutions with better distribution compared to with and without the BU methods.The BNH test problems indicated that the hybrid method produced promising results for IGD and HV values compared to with BU and without BU approaches.• Although the BU and hybrid methods outperformed the other algorithms in single-and multi-objective Vol:.(1234567890) Scientific Reports | (2024) 14:4816 | https://doi.org/10.1038/s41598-024-54841-zwww.nature.com/scientificreports/settingsoptimization, they are still case-dependent, and there are EAs that are not suitable in some cases.

Table 43 .
Summary of the p-value of the posthoc_nemenyi_friedman test evaluated over all runs (sampling).a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.

Table 44 .
Summary of the p-value of the posthoc_nemenyi_friedman test evaluated over all runs (sampling).a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.

Table 45 .
Summary of the p-value of the posthoc_nemenyi_friedman test evaluated over all runs (Crossover).a Posthoc_nemenyi_friedman: +, ~, and -presents the first method performs statistically significantly better, equal, and worse than the second approach.