Combined SGC-Ball Interpolation Curves: Construction and IGEO-Based Shape Optimization

: With the swift advancement of the geometric modeling industry and computer technology, traditional generalized Ball curves and surfaces are challenging to achieve the geometric modeling of various complex curves and surfaces. Constructing an interpolation curve for the given discrete data points and optimizing its shape have important research value in engineering applications. This article uses an improved golden eagle optimizer to design the shape-adjustable combined generalized cubic Ball interpolation curves with ideal shape. Firstly, the combined generalized cubic Ball interpolation curves are constructed, which have global and local shape parameters. Secondly, an improved golden eagle optimizer is presented by integrating Lévy ﬂ ight, sine cosine algorithm, and di ﬀ erential evolution into the original golden eagle optimizer; the three mechanisms work together to increase the precision and convergence rate of the original golden eagle optimizer. Finally, in view of the criterion of minimizing curve energy, the shape optimization models of combined generalized cubic Ball interpolation curves that meet the C 1 and C 2 smooth continuity are instituted. The improved golden eagle optimizer is employed to deal with the shape optimization models, and the combined generalized cubic Ball interpolation curves with minimum energy are a tt ained. The superiority and competitiveness of improved golden eagle optimizer in solving the optimization models are veri ﬁ ed through three representative numerical experiments.


Introduction
Computer-aided geometric design (CAGD) was first proposed by R. Bernhill and R. Riesenfeld in 1974 [1], which is an emerging, borderline interdisciplinary field involving mathematics, computer science, industrial design, and manufacturing. CAGD has not only been widely applied in the three major manufacturing fields of aviation, shipbuilding, and automotive but also, with the continuous improvement and development of CAGD theory, it has been widely applied in mechanical design, computer vision, bioengineering, animation production, three-dimensional medical imaging, military simulation bioengineering, terrain, and other related fields, with very broad application prospects [2][3][4][5]. The emergence and development of CAGD have led geometry from the traditional era to the digital information age, with significant cross-generational significance.
The representation and approximation of curve and surface modeling is an important research topic in CAGD, mainly focusing on the geometric shapes of various products. Interpolation is an important method for approximating discrete functions, and constructing an interpolation curve to interpolate the given discrete data points has important research value. Therefore, many scholars have discussed this issue, and polynomial interpolation is the most common function interpolation method. Lagrange interpolation polynomial and Newton interpolation polynomial are classic polynomial interpolation expressions. Hermite, piecewise, and spline interpolation are among the other wellknown interpolation techniques. The most popular interpolation method among them is spline interpolation, which is segmented and can make the interpolation curve smoother. As a special representation of spline curves, Bézier curves have many excellent properties, such as simplicity, flexibility, and ease of implementation, which have been widely used as interpolation curves [6][7][8][9][10][11][12][13]. Among them, Jaklič studied the Lagrangian interpolation problem of cubic Bézier curves in rational spaces [6]. Mao proposed the subdivision scheme based on normal cubic Bézier curve interpolation [7], which effectively solves the problem of cubic Bézier curve interpolation. Harada applied the Bézier curve to data interpolation, and the resulting interpolation curve was able to fit and approximate the data well [10]. Wahab combined the concept of intuitive fuzzy sets to construct a cubic Bézier interpolation curve [12]. Through experiments, it is found that the fitting effect of Bézier curves is ideal. The Ball curve, as a generalization of the Bézier curve, not only inherits many excellent characteristics of the Bézier curve but also greatly beats the Bezier curve in terms of computational speed and efficiency, so using Ball curves as interpolation curves has important research significance.
In recent years, numerous academics have conducted in-depth research on the Ball curve. Traditional generalized Ball curves and surfaces face the limitation of inflexible shape adjustment; this is because the basis functions of generalized Ball curves and surfaces do not carry shape parameters, so control vertices and corresponding basis functions are the only foundation for determining the shape of generalized Ball curves and surfaces. Because adjusting the appearance of the curve and surface requires modifying the position of the control points, the interpolation curve and surface design are particularly inflexible and difficult to adjust to the requirements of the daily industrial design. In response to this issue, famous scholars both domestically and internationally have conducted indepth research on it and have created shape-adjustable generalized Ball curves and surfaces, successfully solving the shape modification problem of curves and surfaces. In 2021, Hu introduced multiple control parameters and proposed new generalized cubic Ball basis functions [14]. Based on this, the SGC-Ball curve was constructed, which flexibly controlled the shape of the Ball curve through multiple parameters, greatly improving the adjustability of the curve shape. Meanwhile, By combining the effects of several form parameters, Hu introduced the CG-Ball curve in 2022 [15], which can flexibly adjust the appearance of the curve via multiple shape parameters. The cubic generalized Ball surface, including multiple shape parameters, commonly known as the CCG-Ball surface, was created by Hu in 2023 and makes it simpler to build complex surface shapes [16]. Therefore, this paper constructs a class of combined SGC-Ball interpolation curves, which contain overall and local shape parameters and can better control the curve shape. Second, the shape optimization of combined SGC-Ball interpolation curves with different modelings is studied based on the energy method [17][18][19][20][21]. We develop the shape optimization models of combined SGC-Ball interpolation curves with C 1 and C 2 smooth continuity, taking the minimal energy of combined SGC-Ball interpolation curves as the objective function. The aim function of shape optimization models is non-linear, which is inconvenient to solve directly with ordinary techniques, so we use an intelligent optimization algorithm to solve them. The meta-heuristic optimization algorithms have been extensively employed to handle various complex global optimization problems [22], and they are frequently applied in engineering [23], image processing [24], medical diagnosis [25], feature selection, and other fields [26,27]. Famous meta-heuristic algorithms include PSO [28], DE [29], GSA [30], etc., and they have demonstrated good performance in dealing with issues involving global optimization. There are different meta-heuristic algorithms, according to NFL [31], such as GA [32], SA [33], GBO [34], ASO [35], ACO [36], GWO [37], HHO [38], SOA [39], CSA [40,41], AHA [42], ARO [43], SCSO [44], TLBO [45], GJO [46], and DTCSMO [47] etc.
Among them, GEO is a swarm intelligence optimization algorithm proposed in 2021, which has been improved by many scholars due to its competitiveness. It is worth mentioning that an improved golden eagle optimizer (IGEO) was proposed, which boosted the precision and convergence rate of the algorithm [48]. The established shape optimization models of combined SGC-Ball interpolation curves will be solved based on IGEO in this study. The following are the main contributions of this paper:

•
Based on the SGC-Ball basis functions, a class of shape-adjustable combined SGC-Ball interpolation curves is constructed, which have global and local shape parameters; • IGEO is proposed by integrating LF, SCA, and DE into the original GEO; • The shape optimization models of combined SGC-Ball interpolation curves that meet the C 1 and C 2 smooth continuity are instituted based on the criterion of the curve's energy minimization; • IGEO is used to deal with the instituted shape optimization models, and the combined SGC-Ball interpolation curves with minimum energy are achieved.
The remainder of this paper is structured as follows: Section 2 gives the related literature review of Ball curves and the GEO algorithm. Introducing the constructed combined SGC-Ball interpolation curves and studying their smooth splicing conditions in Section 3. Section 4 introduces the detailed process of IGEO. In Section 5, the shape optimization models of combined SGC-Balll interpolation curves are built according to the energy method, the concrete steps of IGEO solving the models are given, and three representative numerical examples are displayed. The last part summarizes and looks forward to the work of this paper.

Related Works on the Ball Curve
Mathematician Ball [49][50][51] created the rational cubic Ball curve in 1974 and utilized it as the foundational theory for CONSURF body surface modeling in order to satisfy the design needs of complex curves and surfaces in industrial products. However, as the traditional Ball curve only has three degrees, it is difficult to satisfy the requirements of complex curve modeling design in the geometric of industrial products. Therefore, Wang proposed the Wang-Ball curve in 1987, which provided a powerful method for the evaluation of high-order curves [52]. In 1989, mathematician Said proposed Said-Ball curve to extend the traditional cubic Ball curve, which proved that the curve has many excellent geometric properties [53]. Considering that the curve is only limited to odd degrees, Hu extended it from odd degrees to even degrees in 1996 [54]. Subsequently, the dual basis functions of generalized Said-Ball curves were independently studied by Othnan [55], Xi [56], and Ding [57]. The n-times Said-Bézier and Wang-Said Ball curves were proposed by Wu in 2000 [58]; both of these curves allow for flexible adjustment of the curve's shape through the use of positional parameters. The cubic Ball curve carrying one parameter was proposed in 2008 [59], and the quartic Wang-Ball curve containing the shape parameters was defined by Wang in 2009 [60]. In addition to sharing characteristics with cubic/quartic generalized Ball curves, they also allow the user to modify the shape via the shape parameter, which increases the flexibility and applicability of creating complex curves and surfaces. The n-th Said-Ball curve with shape parameters was suggested by Xiong in 2012 [61], while the n-th Wang-Ball curve with shape parameters was proposed in 2013 [62]. In addition to having the same great qualities as the Bézier curve and surface, the generalized Ball curve and surface can additionally tackle issues like degree elevation and higher order curve evaluation more quickly and effectively [63].
In recent years, numerous academics have further studied Ball interpolation curves, which are the extensions of Ball curves [64][65][66][67][68]. To address the issue of positive data visualization, Jaafar created the C 1 piecewise rational cubic Ball interpolation curve with three shape parameters [64]. Hasan proposed the C 1 rational cubic Ball interpolation curve, the shape of the curve can be adjusted through parameters, and the condition of maintaining monotonicity was derived [65]. A positivity-preserving rational cubic Ball interpolation scheme was proposed by Jamil and Piah; the scheme had parameters that can be used either to preserve the shape of the curve or to increase its smoothness [66]. The application of the GC 1 cubic Ball interpolation curve in data convexity preserving was studied, and sufficient condition was given for the cubic Ball interpolation curve to maintain convexity [67]. Karim discussed the problem of using a rational cubic Ball interpolation curve with two parameters to maintain positivity and provided the sufficient condition for the positivity of one parameter while the other parameter is a free parameter that controls the final shape of the interpolation curve, resulting in the smoothness of C 1 [68].

Related Works on GEO
In 2021, Mohammadi proposed the golden eagle optimizer (GEO) [69], a revolutionary swarm intelligence system that copies the natural hunting action of golden eagles. The primary characteristic of GEO is the golden eagle's spiral cruising technique, which allows for seamless and clever transitions between exploration and development while attacking prey in a straight route. Due to the relatively competitive performance of GEO, numerous academics have improved it and successfully employed it to resolve a variety of optimization issues. In order to deal with the path planning issue of unmanned aerial vehicle power inspection, Pan suggested a golden eagle optimizer (GEO-DLS) with a dual learning technique based on GEO [70]. In 2021, Soheil used the multi-objective golden eagle optimizer to optimize formulas and achieve ideal results [71]. In 2022, Chandran optimized dual-channel capsules using GEO to generate adversarial networks for evaluating bone age in children from X-ray images of the hands [72]. Kumar used the gradient descent-based golden eagle optimizer (HGDGEO) to solve the efficient heterogeneous resource scheduling problem of cloud-based big data [73]. In 2023, Boriratrit proposed using a golden eagle-optimized extreme learning machine combined with adaptive meta-learning to predict incomplete data in solar irradiance, which can predict ideal data results [74]. Hu proposed a competitively improved golden eagle optimizer (IGEO) based on the original GEO introduced into Lévy flight (LF), sine cosine algorithm (SCA), and differential evolution (DE), which was suggested in 2022 to more effectively address the global optimization problem [48]. IGEO enhanced the performance of standard GEO, which improved the precision and convergence rate of the algorithm while successfully resolving global optimization problems.

Combined SGC-Ball Interpolation Curves
Given a group of discrete data points ) ,..., The expression of the SCG-Ball curve [34] is as follows: where where λ is the global shape parameter of the k-th SGC-Ball interpolation curve, are the local shape parameters of the k-th SGC-Ball interpolation curve.
Secondly, the four control vertexes of the k-th SGC-Ball interpolation curve can be determined through simple derivation and calculation according to Equation (2), resulting in: Therefore, the definition of combined SGC-Ball curves are as follows: Definition 1. Given n ordered interpolation points in the space, the combined SGC-Ball interpolation curves can be defined as: among them, in Equation (6), λ represents the global shape parameter, which controls the overall shape of the entire combined SGC-Ball interpolation curves.
(k = 1,2,..., n) are the local shape parameters, which can be employed to modify the local shape of combined SGC-Ball interpolation curves.
To ensure that the constructed combined SGC-Ball interpolation curves are meaningful, the following theorems and their proof process are given below.
are C 1 smooth and continuous.
Proof of Theorem 1. Firstly, based on the properties of the endpoints, a simple calculation can obtain the following equation Therefore, the combined SGC-Ball interpolation curves are C 1 continuous at the interpolation points ) , , in other words, the combined SGC-Ball interpolation curves are C 1 continuous as a whole. □ shape parameters of the k-th and the (k + 1)-th SGC-Ball interpolation curves meet: If the combined SGC-Ball interpolation curves are C 2 continuous at all interpolation points ) , , then the entire combined SGC-Ball interpolation curves are C 2 smooth continuous.

Proof of Theorem 2.
If the combined SGC-Ball interpolation curves are C 2 continuous, then the curves should first be C 1 continuous, as given via Theorem 1.
Secondly, if the k-th and (k + 1)-th SGC-Ball interpolation curves are C 2 continuous at interpolation point According to the properties of the endpoint, it can be concluded that According to Equations (12)-(14), it can be concluded that The combined SGC-Ball interpolation curves are C 2 continuous curve at the interpolation point k H , then Theorem 2 is proven. □

Standard Golden Eagle Optimizer
A recently developed meta heuristic technique called golden eagle optimizer (GEO) [68] simulates the golden eagles' hunting habits. The three stages of the golden eagle's hunting behavior are the attack stage, the cruise stage, and the transition step from exploration to exploitation.

Population Initialization
To create the initial population X of golden eagles, GEO adopts the random initialization method using the following formula: where r is an even number generated between [0, 1].

Attack Stage
The entire process of golden eagle attacking prey can be simulated using Equation (17): where, ] ,..., , [ represents the attack vector, → * f X represents the optimal prey position, and ] ,..., , [ is the present location of the golden eagle.

Cruise Stage
The golden eagle is presently cruising to locate its prey in the best possible spot. Equation (18) provides the hyperplane of golden eagle's cruise vector: denotes the location chosen by the golden eagle for its prey.
The following equation can be used to simulate the cruise vector ] ,..., , [ in Equation (19), k c represents the k-th element of i → C . The cruise vector leads the golden eagle to search for prey in areas beyond memory, emphasizing the exploration stage of the algorithm.

Transition from Exploration to Development
The golden eagle starts the hunt by cruising in search of prey, but as the hunt progresses, it gradually shifts from cruising to attack and attack. The following expression can be used to model this procedure: among them,

Update Position
Golden eagle's step vector in iteration t is described as follows: in which, are the euclidean norms of the attack vector and the cruise vector, respectively, and their calculation formulas are as Equation (22).
Replace the new location of the iteration t + 1 Golden Eagle using Equation (23):

Improved Golden Eagle Optimizer
Although the original GEO performed successfully in the area of global optimization, it still had certain drawbacks in the domain of handling some more challenging optimization issues, such as premature convergence and low precision. An improved golden eagle optimizer (IGEO), based on GEO, was suggested in order to make the algorithm more competitive. Incorporating LF mechanism into the initialization period can effectively enhance the initialization quality of solution and strengthen the exploration capability of the algorithm. The introduction of SCA in the cruise phase is conducive to the golden eagle's random walk, thus increasing the convergence speed of original GEO. In purpose of preventing GEO from sinking into the local optima too early, the idea of DE was incorporated into the attack stage of golden eagle, expanding the probability of population mutation and the ability to learn from the best individuals, accelerating convergence rate and improving solution precision. The following are the precise IGEO algorithm implementation steps:

Lévy Flight Strategy
A random walk approach called Lévy flight (LF) [75] can significantly improve population diversity since it complies with the Lévy distribution. The update formula is as follows: where, best x represents the optimal position, 5 . 1 = β is the Lévy random path. Introducing the LF technique during the initialization phase can enhance the quality of the initial population and increase the rate of algorithm convergence towards the optimal solution. The detailed process is as Equation (25).

Sine Cosine Algorithm
GEO determines the location of prey through information exchange between populations, which greatly limits the search scope and limits the accuracy of the solution. A popular stochastic optimization technique with straightforward concepts, robust operability, and simplicity of implementation is the Sine Cosine Algorithm (SCA) [76]. It has been extensively used to enhance numerous algorithms. The update equation is as follows: , and t i p is the position of the optimal individual in the golden eagle population.
in Equation (27), 2 = a is a constant. Using SCA to strengthen the exploration and development of GEO, boost population diversity, raise the algorithm's capacity to quit local optima, and more evenly distribute the exploration and development of algorithms.

Differential Evolution
In the exploration phase of GEO, the golden eagle did not consider learning from other individuals in the golden eagle population when searching for new prey, which has certain limitations. The first global optimization approach put forth by Storn et al. is called differential evolution (DE), and it primarily uses mutation, crossover and selection procedures [29]. Its benefits include quick convergence, high precision and simple implementation. The following are the precise operations of this algorithm: (1) Mutation operation. It is customary to randomly choose three individuals from the group, scale their vectors, and then combine them with the individuals who will undergo the mutation in order to increase the population's variety through mutation. The expression is as follows: where r1, r2 and r3 are between [1,T], is called the scale factor, F0 = 0.5 is mutation rate.
(2) Cross operation. The following is the expression for mimicking mutation operations: among them, the probability that at least one member of the population will experience a mutation is known as the crossover rate , and rand j is within [1, D].
(3) Selection operation. The following are examples of the greedy selection strategy: Integrating differential evolution into the development stage of GEO allows golden eagle to learn from other random individuals, increase the mutation rate of the population, and enrich the diversity of the population. This can reduce the probability of premature convergence, greatly improve the accuracy, and enable the algorithm to converge towards the optimal solution faster and more accurately.

The Shape Optimization Models of Combined SGC-Ball Interpolation Curves
The energy measure is a common and effective method to establish the curve's optimization model. With the objective function of minimizing curve's energy, this part establishes the shape optimization models for combined SGC-Ball interpolation curves. Assuming is the energy value of the k-th segment in the combined SGC-Ball interpolation curves, it can be defined as: where, ) , , , ; ( Then Ẽ represents the energy value of the entire combined SGC-Ball interpolation curves, which can be expressed as: Taking the minimum energy of combined SGC-Ball interpolation curves as the objective function, the shape optimization model of combined SGC-Ball interpolation curves is established as follows:  (33) in Equation (33), there is Therefore, using the energy method as the measurement standard and regarding the shape parameters as the variables to be optimized, the shape optimization models of combined SGC-Ball interpolation curves that satisfie C 1 and C 2 smooth continuity can be established, and their objective functions can be recorded as: The constraints of C 1 and C 2 smooth continuity are given below: 1. Overall C 1 smooth and continuous 2. Overall C 2 smooth and continuous  (38) Due to the high nonlinearity and complexity of the objective function of the shape optimization models established by Equations (36)- (38), solving them directly using traditional optimization techniques is a great challenge, so this paper considers utilizing IGEO to deal with the shape optimization models, to obtain the interpolation curves with the optimal shape.

Steps for IGEO to Solve Combined SGC-Ball Interpolation Curves Shape Optimization Models
This section will provide a thorough explanation of how to solve the established combined SGC-Ball interpolation curves shape optimization model utilizing the proposed IGEO algorithm. These are the precise steps: Step 1. Set the relevant initial parameter parameters; Step 2. Initialization. In combination with Lévy flight to initialize population, the energy value of the combined SGC-Ball interpolation curves is considered as the fitness function value, and the energy value of every individual is calculated; Step 3. Update a p and c p by Equation (20); Step 4. Obtain the attack vector → A by Equation (17), and apply Equation (19) to achieve the cruising vector → C ; Step 5. Calculate the mutant population ) , by Equation (28) and Equation (29), and Calculate the sub-population ) , by Equation (30); Step 6. Apply Equation (26) to find the location of the golden eagle in help of SCA, denoted as SCA i x ; Step 7. If Step 9. Calculate the iteration step x Δ of golden eagle using Equations (21) and (22); Step 10. Replace the location of the golden eagle via Equation (23); Step 11. Duplicate step 3 to step 11 untill T is attained, and output the minimum energy value Ẽ and corresponding shape parameters value of the combined SGC-Ball interpolation curves.

Modeling Examples
This part gives three representative shape optimization examples of combined SGC-Ball interpolation curves. Firstly, the combined SGC-Ball curves are used to interpolate the given data points, and then the shape optimization model of the combined SGC-Ball interpolation curves is builded. Finally, IGEO is employed to resolve the shape optimization model, further confirming the outstanding of the proposed IGEO. Example 1. The coordinates of the data points ) 10 ,..., to be interpolated are given as follows:   (39) By using the combined SGC-Ball interpolation curves to interpolate the given data points ) 10 ,..., , the model graph satisfying C 1 smooth and continuous "snake" shape was obtained.
Different optimization algorithms were applied to handle the shape optimization model of combined SGC-Ball interpolation curves, and the C 1 smooth and continuous combined SGC-Ball interpolation curves with the lowest energy were obtained. Figure 1a shows the model diagram of the C 1 smooth concatenation of the combined SGC-Ball interpolation curve when 10 , , 1 , The optimization curves obtained by solving the C 1 smooth concatenated combined SGC-Ball interpolation curves using various optimization algorithms are shown in Figure 1bh. It can be observed that the constructed "snake" interpolation curves consist of 10 SGC-Ball interpolation curves, with different colored curves representing different SGC-Ball interpolation curves, as well as the shape of the interpolation curves optimized by different algorithms is different, but the difference is not significant, which is caused by different parameter values. The energy convergence changes for all algorithms are drawn in Figure 1i. It is discovered that the IGEO-obtained energy convergence curves can swiftly converge to the least energy, demonstrating that IGEO is better to other comparison methods. Meanwhile, Table 1 provides the minimum energy and related shape parameter values of each algorithm for solving C 1 smooth continuation "snake" interpolation curves. The results obtained are the average values of the results obtained by each algorithm after 30 independent runs. According to Table 1, the energy obtained by IGEO is the minimum, with a value of 69.6892, which means it has the smoothest curve shape. Compared with other algorithms, IGEO has the best performance, competitiveness, and superiority in solving the C 1 smooth concatenation shape optimization model of combined SGC-Ball interpolation curves.  to be interpolated are given as follows:  Table 2 displays the minimum energy values and corresponding parameters after 30 runs. Observing these data, it is evident that compared with other algorithms, IGEO obtained the smoothest curve, with an energy value of 90.9885. Through this numerical example, it is proven that IGEO achieves the best results in solving the established combined SGC-Ball interpolation curves, and can obtain curve shapes with ideal shapes.  to be interpolated are given as follows:  h depict the smooth curves obtained by seven optimization methods. From a visual perspective, it can be seen that the curves are optimized by IGEO is the smoothest, indicating that IGEO has certain advantage in solving the optimization model. Figure 3i describes the energy change convergence curves of all algorithms, and it is found that IGEO exhibits better convergence performance when solving shape optimization problems.  Table A1 displays the experimental results of each technique to handle the curves' shape optimization model, including optimal energy value and related parameter values, in order to more thoroughly examine the optimization impact of every technique to solve the curves shape model (see Appendix A). Compared with other competitive algorithm, IGEO achieved the minimum energy value of 358.0383, which is the combined SGC-Ball interpolation curves with the most ideal shape.

Conclusions and Future Works
This article first constructs the combined SGC-Ball interpolation curves and studies the conditions of the C 1 and C 2 smooth continuity of two adjacent SGC-Ball interpolation curves. The constructed combined SGC-Ball interpolation curves have global and local shape parameters, and have high flexibility in shape adjustment. Secondly, three improvements were integrated into the GEO: (1) The variety of the population increased after LF was introduced during the population's initiation phase. (2) SCA enhances the development capacity and search scope of original GEO. (3) Using DE to escape local optima and increase the solution's precision. The above three strategies work together on GEO, greatly boosting the effectiveness of GEO. Finally, depending on the minimum energy of the curve, the shape optimization models of combined SGC-Ball interpolation curves were established. The shape parameters in the models were optimized using the proposed IGEO to achieve the optimal interpolation curves.
The solution results of IGEO were compared with other algorithms, and three numerical examples verified the superior performance of IGEO.
At present, the theoretical system and key technologies of generalized Ball interpolation curves and surfaces with multiple shape parameters are not yet perfect. This article only studies the construction, C 1 and C 2 smooth splicing, as well as shape optimization problems of combined SGC-Ball interpolation curves, further research is needed on the segmentation and extension algorithms of this curve in future work. In addition, we will consider extending the research technique of combined SGC-Ball interpolation curves to the CQGS-Ball surfaces in [77], and utilizing intelligent algorithms in [78] to investigate the shape optimization problem of the surfaces. Data Availability Statement: All data generated or analyzed during this study were included in this published article.