Development of Tolerance-Based Performance Prediction Technology and Optimization of Actuator Design Factors of a Magnet Vertical Magnetization of AVAS

: With the increasing proliferation of electric and hydrogen vehicles, noises to recognize the driving status at low speeds are legalized, so a virtual engine sound generator is required, and slimming is required for packaging it in vehicles. This study investigates an optimization method for improving the electromagnetic force performance and slimming of the magnetic circuit for the permanent magnet structure for the vertical magnetization of the actuator for the acoustic vehicle alerting system (AVAS) of a vehicle and the probabilistic optimization of manufacturing tolerance management. To investigate the impact of the design parameters of the magnetic circuit structure on the electromagnetic force performance and slimming, we performed an independent analysis based on a single variable and investigated the characteristic variations based on multiple variables using a full factorial design and derived a performance prediction regression model using the central composite design of response surface methodology, including the curvature effect, by adding a center point to verify and consider the nonlinear characteristics. Consequently, four effective design parameters were determined to analyze the electromagnetic force performance and slimming of the vertical magnetization structure of the AVAS actuator—permanent magnet thickness, magnetic force collecting plate thickness, yoke position, and yoke thickness. We then performed statistical analysis using Monte Carlo simulation and proposed an optimization management level of 3σ with excellent process capability as the design application tolerance that can occur in the manufacturing process of each design parameter, whereby the confidence level of electromagnetic force performance and slimming improved from 99.46% to 99.73% and 97.62% to 99.73%, respectively.


Introduction
The automotive industry is shifting its focus from vehicles with internal combustion engines to electric vehicles, powered by advancements in technology and increased environmental awareness of global warming. It is predicted that the use of eco-friendly vehicles such as electric vehicles (EVs) will continue to expand in the future. One of the advantages of vehicles using electric motors is that there is no noise pollution, which is in contrast to those using internal combustion engines [1]. However, the low noise level of an EV running at a low speed can lead to a situation where pedestrians cannot be acoustically alerted to the approaching vehicle, thus increasing the risk of accidents.
Vehicles generate noise owing to driving speed, powertrain, tire-to-road contact, and aerodynamic sound. When a vehicle is driven at a high speed, the noise generated by the aerodynamic sound and tire-to-road contact dominates over other noise sources, resulting in a sufficiently high noise level that can be easily heard by pedestrians and people in other vehicles. In the case of a vehicle with an internal combustion engine, pedestrians can recognize the movement of the vehicle by the noise generated by the combustion of air and fuel inside the engine even when it is driven at a low speed [2]. However, an electric vehicle running at a low speed (<30 km/h) has a low noise level, making it difficult for pedestrians to hear [3]. Moreover, when the background noise is high, pedestrians often fail to perceive the approaching vehicle early enough to avoid the risk of an accident [4,5].
In 2008, the University of California, Riverside, tested the time for a pedestrian or cyclist to acoustically perceive a vehicle running at a low speed. Consequently, it was observed that background noise, such as the noise generated by an idling vehicle motor, makes it impossible for pedestrians to recognize and cope with the approaching vehicle [6][7][8]. The US National Highway Traffic Safety Administration announced in its 2009 statistical report that when moving at a low speed, an EV is more prone to cause accidents than an internal combustion engine vehicle [1,9]. The US and EU countries mandate the installation of an acoustic vehicle alerting system (AVAS) on an electric or hybrid vehicle so that pedestrians and visually impaired people can perceive its movement from a minimum level of noise against road traffic noise and ambient background noise [10][11][12][13]. An AVAS for hybrid or electric vehicles is a useful device for cyclists and pedestrians, especially for the blind. However, there are arguments against the installation of AVASs because they can generate additional noise pollution. In this context, many studies have been conducted to minimize the side effects of AVASs and generate an effective alarm sound.
Felix and Alexander measured and analyzed the AVAS warning sound of several EVs and confirmed the difference in noise measured from the left and right of each vehicle [14]. Wagner and Gajda assumed that the manufacturer should comply with the AVAS minimum sound pressure level on all sides of the vehicle, making the warning sound excessive in a specific direction depending on the AVAS speaker installation method [15]. To address this problem, studies are underway, and it has been verified that high-quality warning sounds can be generated by creating a sound field controllable with a loudspeaker-based system [16,17]. However, although a loudspeaker can effectively emit sound, its installation requires vehicle design modification, resulting in increased production costs and interference with other vehicle systems. In addition, since the loudspeaker cone is exposed to external environments, such as wind, dust, and moisture, the durability of the system may be a problem [18]. Mounting a standard round loudspeaker inside the engine room in consideration of exposure to the external environment poses the problem of a warning sound transmitted only in a specific direction [19].
As an alternative to a speaker-based AVAS, research is being conducted on an actuator-based system that forces the vibration of a panel. The panel actuator system does not require a vehicle structure design change because it can be simply attached to an existing panel or vehicle. Moreover, it is more durable than the conventional loudspeaker-based system because of no exposure to the external environment [18,20]. Kournoutos et al. reported that an actuator system mounted on a panel reproduces the sound performance of a loudspeaker-based system and outperforms the sound characteristics of the speakerbased system at a certain frequency [21]. Given that the panel vibration AVAS has an omnidirectionally radiated sound field, the main focus of research has been on the sound field and control of the system by structural vibration [18,20,[22][23][24][25], and little research has been dedicated to the shape and performance optimization of the actuator of the panel vibration type.
In general, electrical energy is converted into vibration (acoustic) energy in a vibration actuator or loudspeaker system, and research has been conducted to improve its performance [26][27][28][29][30]. Ko et al. studied the electromagnetic force characteristics and vibration performance of a permanent magnet according to its magnetization direction [31]. Since the basic architecture of a panel-type actuator is similar to that of a general actuator and loudspeaker system, its excitation can be improved by optimizing the electromagnetic cir-cuit. The size and shape of the actuator are factors that affect panel excitation and vibration, and their performance can be improved by optimizing the magnetic circuit considering the design boundaries.
Against this background, this study compares the electromagnetic force performance of an AVAS actuator mounted on a vehicle hood depending on the magnetic circuit structure, selects the optimal magnetic circuit structure, and optimizes the design parameters and design tolerances for the electromagnetic force performance and slimming of the selected magnetic circuit structure. As a result of design factor analysis using design of experiment (DOE), the magnetic circuit design parameters applicable to electromagnetic force performance were set as the design parameters of each magnetization structure. With the selected design parameters, we optimized the magnetic circuit using the desirability function as a multi-objective function after deriving a nonlinear regression model using response surface methodology. Subsequently, we presented an optimum tolerance design capable of securing the electromagnetic force characteristics by minimizing the production and assembly process deviations, using the Monte Carlo simulation (MCS), to ensure a 3σ quality-level process capacity.

Architecture of the Virtual Engine Sound Actuator
The architecture of the virtual engine sound actuator is a structure with a vertically vibrating interior diaphragm under the influence of the magnetic field generated by the permanent magnet and coil, as illustrated in Figure 1, whereby the permanent magnet is magnetized in the vertical direction. Regarding the materials of the actuator components, magnetic substances were used as the components affecting the formation of the magnetic field (yoke, permanent magnet, magnetic force collecting plate, and cap) and paramagnetic substances for the components that do not affect the formation of the magnetic field (spring, coil, and frame), as defined in Table 1. The coil (copper wire) was wrapped with 238 turns of wire and had a resistivity of 8 Ω. Although it is not magnetized in a magnetic field, it forms part of the magnetic circuit, which affects the formation of the magnetic field when electricity is applied to the wire.

Electric Field of the Initial Model and Magnetic Force Analysis
Electromagnetic Field of Actuator and Vibration and Response Analysis Electromagnetic field analysis used FEMM 4.2 software. The virtual engine sound actuator generates electromagnetic force by the current flowing through the coil and the magnetic field of the permanent magnet in accordance with Fleming's left-hand rule, and the electromagnetic force is expressed as Equation (1) in accordance with the Lorentz force law.
where F is the electromagnetic force, I is current, l is coil length, and B ⃗ ⃗ is magnetic field. The operating conditions of the virtual engine sound actuator are as follows: rated voltage = 12.65 V, impedance = 8 Ω, and current = 1.65 A. Table 1 outlines the properties of the virtual engine sound actuator. The basic design geometry of the magnetic circuit structure has cylindrical cyclic symmetry, and we performed an electromagnetic field analysis using a two-dimensional plane model. The electromagnetic field analysis model was constructed by applying the Robin boundary condition (mixed boundary condition), which can define both the Dirichlet boundary condition of zero tangential electric potential (V = 0) and the Neumann boundary condition that constrains the magnetic flux passing in the normal direction as the boundary condition to maintain the magnetic flux, as illustrated in Figure 2. Electromagnetic field analysis was used FEMM 4.2 software. The electromagnetic field analysis enabled the calculation of the electromagnetic force generated by the coil current and the magnetic field formed by the yoke, magnetic force, and plate, as well as the magnetic flux density generated by the coil and the permanent magnet, as illustrated in Figure 3. The electromagnetic force was calculated at 19.80 N using Equation (1), which does not satisfy the design target of 20.00 N or higher. To meet the design target, a slimming optimization is necessary.

Overview of the Magnetic Circuit Optimization Process
The electromagnetic force of the actuation was optimized using DOE, as illustrated in Figure 4. To set the design parameters, we selected the factors that may be associated with electromagnetic characteristics and analyzed their effectiveness using a full factorial design. The effective factors were subjected to factor screening, and the effectiveness and curvature were analyzed based on the response values for the electromagnetic force. The nonlinearity of the design parameters was verified using the central composite design of response surface methodology, and the regression model was used for design parameter optimization with the multi-objective function of the electromagnetic force, magnetic circuit, and cap distance.

Selection of Design Parameters for the Magnetic Circuit and Factorial Design
We set the selection criteria for actuator design parameters as geometric factors that can minimize the volume necessary for the package while maximizing the electromagnetic force. As geometric variables affecting the performance [32], we selected the magnet thickness (MT), plate thickness (PT), yoke position (YP; distance between yoke and coil), and yoke thickness (YT), as shown in Figure 5, with the size based on the outer diameter (Φ = 65 mm) fixed. The range of design of the selected parameters was set as the range of design boundaries and the direction minimizing the magnetic circuit thickness, as defined in Table 2. The minimum yoke position was limited to 2.5 mm, which is the minimum separation distance necessary to prevent interference with other components. Among the selected design parameters, the one with the greatest effect was identified using factor analysis of the factorial design methodology [33]. As outlined in Table 3, the factorial design of the selected parameters was implemented as a total of 17 combinations consisting of 16 combinations of a four-factor two-level full factorial design, based on a two-level factorial design for the maximum and minimum values, and one center point of curvature analysis, expressing the electromagnetic force obtained by electromagnetic field analysis as response values. The analysis of variance (ANOVA) of the factorial design of the selected parameters is the adjusted sum of squares (Adj SS), as expressed by Equation (2).
where T1 and T0 are the sums of the response values for the high-level and low-level factors, respectively, and N is the number of experiments. Adj MS is the mean square obtained by dividing Adj SS by the degree of freedom (DF) presented in Table 4. The effect of the design parameters can be verified via the analysis presented in Table  4, based on the importance of the response values for the electromagnetic force. The Fvalue can be obtained by dividing the Adj MS of each factor by the Adj MS of the error, whereby the higher the value, the higher the effect of the related factor on the response value. The p-value is inversely correlated with the F-value, and the significance is considered established at 0.05 or lower. Accordingly, Table 4 lists the significant probability values for the lack of fit of the response values for the electromagnetic force with regard to the magnet thickness (MT), plate thickness (PT), yoke position (YP), and yoke thickness (YT). The analysis revealed that all four design parameters were statistically significant at the 5% significance level. Figure 6 shows the results of the corresponding Pareto analysis. Its nonlinearity suggests that not only the main factors but also their interactions have significant effects.

Design Parameter Analysis Using a Response Surface Design
Given the nonlinearity of the design parameters of factorial design analysis and the significant characteristics of their interactions, we optimized the performance for volume (thickness, distance) and electromagnetic force using a response surface design conducive to two-level or higher nonlinear prediction, thereby applying the face-centered design generally used to optimize nonlinear design parameters. In the face-centered design, which is a type of response surface methodology, a total of 25 DOE combinations are arrayed, and the analysis results for the electromagnetic force, magnetic circuit thickness, and cap distance are outlined in Table 5. As shown in Table 6, the ANOVA of the face-centered design confirmed the statistical significance of the main effects of the design parameters, square terms, and interaction effects. With respect to the main and interaction effects of the design parameters, from the electromagnetic force performance relationships, the yoke position (YP) and yoke thickness (YT) were identified as parameters that are more sensitive than magnetic thickness (MT) and plate thickness (PT), as shown in Figures 7 and 8, and the interaction effects were sensitive to the correlations among all parameters.   The regression prediction model of the response value of the design parameters can derive the prediction equation from the results of the face-centered design, which can be obtained with the quadratic regression model for the independent variable K as a quadratic function of the face-centered design, using Equation (3).
where X i , X j are the independent variables (design parameters), Y is the dependent variable (resulting value), and βi, βij are coefficients, which can be obtained by least squares

Main Effects Plot for Electromagnetic force
Fitted Means [34,35]. From Equation (3) The distance between the PT and cap must have the minimum value necessary to avoid vibrational interference. Therefore, their geometric relationship can be derived from where D is the distance between magnetic flux collecting plate and cap base, YP is the vertical distance between the coil and yoke, MT is the magnet thickness, and PT is the thickness of the magnetic flux collecting plate. Moreover, the volume of the actuator is correlated with its thickness, with its diameter fixed, and the thickness can be determined using Equation (6) where T is the magnetic circuit thickness, YP is the vertical distance between the coil and yoke, and YT is the yoke thickness. The coefficient of determination (R 2 ) in Equation (4) indicates that the closer the value is to 1, the higher the degree of coincidence of the prediction model of the experiment with the actual value [36]. In this regard, Joglekar and May [37] proposed a coefficient of determination of 0.8, or higher, to establish a good fit of the prediction model. The coefficient of determination of the design parameters for the electromagnetic force was 0.9985 in the regression model, and the adjusted coefficient of determination (R 2 Adj) was 0.9974. Therefore, the model fit of the derived prediction regression model can be established according to the Joglekar and May criteria.
The electromagnetic force and magnetic circuit thickness are in a trade-off relationship, and among the objective functions of this trade-off relationship, we used a multiobjective function as an optimization method. As an algorithm for the multi-objective function, we used a desirability function, which is a method proposed by Derringer and Suich (1980) [38]. After calculating the individual desirability for each response, the composite desirability was derived. The composite desirability function D is the geometric mean of the values of the desirability function (di), as expressed by Equation (7): (D = composite desirability, di = individual desirability, and n = total number of measured responses).
As shown in the above equation, the composite desirability is calculated as the geometric mean because if one of the response variables is not satisfied, the composite desirability is zero [39].
The individual desirability di of the composite desirability function is derived using Equations (8)-(10) according to the response characteristics. Equation (8) is the desirability function when the goal is to maximize the response value (y î (x)) and Equation (9) when the goal is to minimize the response value. Equation (10) is a desirability function when the response value approaches the target value.
where Li = lower limit for response i Ui = upper limit for response i Ti = target value for response i To perform optimization with regard to electromagnetic force and magnetic circuit thickness, we set the following three optimization target functions: (Q1) Maximization of electromagnetic force on the coil (≥20 N); (Q2) A minimum distance of 2.5 mm between magnet and cap; (Q3) Minimization of the magnetic circuit thickness. Also, the geometric relationship of the optimization target functions is shown in Figure 9.
The design parameters of the magnetic circuit structure were optimized using the desirability function of multi-objective functions for prediction Equation (4) and condition Equations (5) and (6). The optimization results are presented in Figure 10. The individual desirability of the electromagnetic force was 0.92, distance between the plate and base plate was 0.70, actuator thickness was 0.53, and composite desirability was 0.76. The corresponding optimum values of the design parameters were: magnet thickness (MT) = 8.0 mm, plate thickness (PT) = 3.5 mm, yoke position (YP) = 4.8 mm, and yoke thickness (YT) = 4.5 mm. That is, slimming was achieved for YP by 0.1 mm, and the minimum cap distance was secured (2.5-3.5 mm) owing to the slimming of PT by 1.0 mm, resulting in structural stabilization. The predicted electromagnetic force (F) was found to satisfy the design target of 20.02 N.

Verification of the Response Characteristics of the Optimization Analysis Model
To test the reliability of the actuator analysis model, we fabricated an experimental model, as shown in Figure 6, and its validity was verified by comparing the vibration characteristics of the numerical analysis model and experimental model.
The actuator vibration analysis was configured as a finite element model ( Figure 11) in which the diaphragm was induced into vertical vibrations by the electromagnetic force. The excitation and boundaries of the electromagnetic force are shown in Figure 12. For the vibration analysis, we analyzed the vibration response characteristics using the mode superposition method.

Stochastic Design Optimization
The methodology for analyzing and designing electronic devices and mechanical systems can be broadly classified into deterministic and stochastic methods, depending on whether the uncertainty of design or system parameters is considered. Of these, a stochastic method considers the uncertainty of design or system parameters using a stochastic model, which has been intensively researched from the perspective of reliability engineering [40][41][42]. In general, the uncertainty regarding electronic devices and mechanical systems is expressed by the term "tolerance" with respect to changes in utilization conditions and errors in the manufacturing and assembly processes.
Thus, tolerance of uncertainty causes uncertainty in the actual product, which can exhibit unexpected or unwanted performance by the designer. Therefore, to ensure more competitive product development and production, stochastic analysis and design considering uncertainties such as tolerances should be accompanied.

Definition of Tolerance and Process Capability
In the mechanical and electronic sectors, system-specific tolerance inevitably exists. The term "tolerance" generally refers to the amount to which variation is allowed. In mechanical engineering, tolerance is defined as the permissible range of variations of a given measure, i.e., the difference between the maximum and minimum values of these variations. Tolerance also means manufacturing limits within a range in which economically feasible manufacturing costs can be determined.
The term "process capability" can be defined as the ability of a manufacturing process to attain a certain level of quality and maintain the quality level in a stable state.
The process capability index (Cp) is obtained from the value calculated from the standard deviation of the upper and lower specification limits (USL and LSL) of the tolerance and the measured value of product performance attained by the related process, as defined by Equation (11).

Current Performance Prediction of the Design Parameters and Tolerance Designs
The basic tolerance of the optimum value for the design parameters of the magnetic circuit is shown in Figure 15 and Table 7. Owing to the difficulty in estimating the level of process capability level in the development stage, the probability distribution of the de-sign parameters is defined as a triangular distribution that can assume an uncertain process capability level. For performance prediction, we performed stochastic analyses of the sensitivity of design parameters, performance distribution of electromagnetic force, and magnetic circuit thickness distribution after generating 10,000,000 random numbers for design parameters using the Markov chain Monte Carlo method.  We defined the assumption of the probability distribution of the basic tolerance design as a triangular distribution that can assume an uncertain process capability level for the four design parameters, as illustrated in Figure 16. The triangular distribution of each design parameter was assumed to be a probability distribution that averages the minimum and maximum tolerances of the design parameters and the median of the tolerances.  The stochastic performance prediction of the basic tolerance design for the four design parameters is presented in Figures 17 and 18, which plot the results of the stochastic analysis of the electromagnetic force performance distribution and magnetic circuit thickness distribution by generating random numbers within the design range for the upper and lower limits of each design parameter, respectively. As show in Figures 19 and 20, considering the variations in the design parameters of the basic design tolerance, the sensitivity of the design parameters of electromagnetic force and magnetic circuit thickness was found to be mostly attributable to the effects of yoke thickness and yoke position, and the sensitivities of the permanent magnet and magnetic force collecting plate were found to be lower. It was found that the electromagnetic force performance distribution of the basic tolerance design variation attained a performance range of 20.01-20.09 N, and a magnetic circuit thickness of 19.3-19.47 mm.

Tolerance Design Optimization for Performance Assurance for Design Parameters
Optimization of the tolerance for the design parameters can maximize the performance of a virtual engine sound actuator through production quality control. In the product design process, tolerance design is performed to control product quality and manufacturing cost once the conceptual and parameter designs are determined. Therefore, design tolerance optimization is a process for determining the best trade-off between performance and quality. By optimizing the tolerance for the design parameters, we sought to analyze the actuator thickness minimization and electromagnetic force performance distribution at the 3σ level, which is an excellent process capability.

Optimization of Design Parameters Tolerance
The tolerance of design parameters for the magnetic circuit was optimized by calculating the standard deviation at the 3σ level of design tolerance for each parameter, and stochastic analyses were performed based on the sensitivity of the design parameters, electromagnetic force performance distribution, and thickness variation by the probabilistic analysis method after generating 10,000,000 random numbers. The standard deviation of each design parameter can be obtained from the 3σ level process capability index (Cp = 1) and Equation (15).
The probability distribution of the design parameters is a normal distribution that averages the standard deviations calculated using Equations (13)- (16), and the median values of the upper and lower limits of the nominal dimension, as defined in Figure 21. The mean and standard deviation of each design parameter are listed in Table 8.  Analysis of the design parameters assumed as a normal distribution, which was performed by applying the Monte Carlo method to the regression model, revealed the electromagnetic force performance distribution ( Figure 22) and magnetic circuit thickness distribution ( Figure 23). When the tolerance of the design parameters was constrained at the 3σ level, the yoke thickness (YT) and yoke position (YP) were found to be most sensitive in terms of electromagnetic force ( Figure 24) and magnetic circuit thickness ( Figure 25), with the sensitivity of YP increasing in tandem with an increase in the process capability. The electromagnetic force performance distribution showed a range of 20.02-20.08 N, and a magnetic circuit thickness distribution of 19.39-19.46 mm. The results of the comparison of the reliability of the electromagnetic force performance distribution and the magnetic circuit thickness distribution against the basic design tolerance and reliability at the 3σ level of 99.73% are plotted in Figures 26 and 27, respectively. With the electromagnetic force ranging from 20.01-20.09 N based on the reliability of the 3σ level and the basic design tolerance of 99.46%, the similarity was confirmed with a difference of approximately 0.23% in reliability. By contrast, the magnetic circuit thickness showed a difference of approximately 2.11% in reliability, with its value ranging from 19.39-19.46 mm and basic design tolerance of 97.62%. From this, it can be inferred that the tolerance design of the magnetic circuit thickness minimization at the 3σ level had a greater effect on reliability compared to the reliability of the electromagnetic force performance. This allows the assumption that process stability can be attained by minimizing the variations in yoke thickness (YT) and yoke position (YP), resulting in tolerance optimization, as shown in Figure 28.

Conclusions
In this study, we analyzed the actuator's operational performance, i.e., the effectiveness of the design parameters for the electromagnetic force performance and electromagnetic characteristics, and investigated the methods for magnetic circuit optimization and stochastic prediction of the performance of design parameter tolerance. As the research methods, we used a full factorial design for verifying the effectiveness and screening the characteristics of the design parameters and performed an experimental verification by applying the response surface design and multi-purpose functions for optimization. We also performed statistical analyses using the Markov chain Monte Carlo method to investigate the effects of the tolerances in the production process of design parameters on performance and verified its usefulness. In conclusion, the effectiveness factors for the electromagnetic force in the vertical magnetization magnetic circuit of a permanent magnet, which is the research object of this study, were identified to be the permanent magnet thickness, magnetic force collecting plate thickness, yoke position, and yoke thickness. With the effective factors showing the second-order nonlinear curvature characteristics including interactions, a response surface design was judged to be the most suitable model to attain optimization by modeling the curved surface characteristics that can consider the nonlinearity. As a result of optimizing the multi-objective function to minimize the boundary conditions of the design parameters and the magnetic circuit thickness and maximize the electromagnetic force, the magnetic circuit thickness was estimated at 19.40 mm, and the electromagnetic force was estimated at 20.02 N. It was found that the electromagnetic force performance can reach a range of 20.01-20.09 N, and the magnetic circuit slimming can reach 19.39-19.46 mm, in terms of the design tolerance and the design production tolerance at the 3σ level based on a reliability of 99.73%. A design tolerance optimization analysis revealed that minimizing the variations in yoke thickness and yoke position is advantageous for electromagnetic force performance, and minimizing the variations of yoke thickness and yoke position in the magnetic circuit thickness is advantageous for slimming an actuator. Regarding the vertical magnetization structure, an increase in yoke thickness (YT) and yoke position (YP) was found to have a positive effect on the electromagnetic field performance, which allowed the conclusion that slimming is possible up to 19.39 mm, even considering the 3σ process conditions. Taking these results together, it may be concluded that optimization using response surface design and optimization of tolerance management based on stochastic performance prediction using the Monte Carlo method are highly viable and a useful method for the optimum design of an AVAS.
Furthermore, we plan to analyze the mass production performance and perform ongoing data verification based on the accumulated data for tolerance management optimization and to pursue research that applies the AI-based optimization approach.

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