Shape Optimization Design for a Centrifuge Structure with Multi Topological Conﬁgurations Based on the B-Spline FCM and GCMMA

: The geotechnical centrifuge applied in various geotechnical engineering ﬁelds provides physical data for investigating mechanisms of deformation and failure and for validating analytical and numerical methods by simulating and studying the geotechnical problems. The basket, as one of the important components used to place the inspection model of centrifugal test, is designed to withstand complex loads. This paper presents an optimization design method for the basket based on the weighted B-Spline Finite Cell Method (FCM) and the globally-convergent method of moving asymptotes (GCMMA). In order to obtain a superior design solution, four topological conﬁgurations, i.e., original single web, porous dual web, open deep groove dual web, and connected closed dual web, are investigated and optimized. The mass is selected as the optimization objective, while key shape parameters and stress are regarded as design variables and the constraint, respectively. By optimization, the ﬁnal masses of the four conﬁgurations are reduced greatly compared with the initial conﬁgurations, where the greatest weight loss, in case 4, is 10.6%. This indicates that the weighted B-Spline FCM and GCMMA can be well applied for shape optimization of structure in engineering design. In contrast to the ﬁnal single web adopted in the traditional basket design in case 1, the ﬁnal conﬁguration in case 4, i.e., connected closed dual web, has the least mass. The ﬁnal mass is reduced by 133.38 kg when the centrifuge strength requirement is met. Therefore, the ﬁnal conﬁguration in case 4, where the maximum von-Mises stress is 398.72MPa and mass is 781.82 kg, is superior to the three other conﬁgurations. this research. This paper presents an e ﬀ ective and systematic structure design and shape optimization for the lifting lug with multi topological conﬁgurations based on the weighted B-Spline FCM and GCMMA. It provides a complete shape optimization process for the lifting lug of the geotechnical centrifuge, which addresses the problem of the light-weight design of the lifting lug while satisfying the requirement of the geotechnical centrifuge strength. Four topological conﬁgurations, i.e., original single web, porous dual web, open deep groove dual web, and connected closed dual web, are investigated and optimized. conclusions be


Introduction
The centrifuge simulating composite environment loads with acceleration is designed to meet the requirements of geotechnical, medical, and bio-industrial research, etc. [1]. In recent decades, along with the development of military and social demands, the centrifuge has attracted continuous research interest. An advanced centrifuge which is capable of high rotational speed and with a more precise structure is seriously needed. Due to the high rotational speed, both the stress level significantly affecting the reliability, and the weight involving the performance of the centrifuge, are sensitive parameters. Several cases involving centrifuge failure have been reported. A centrifuge was destroyed due to excessive mechanical stress in a Beckman L2-65B ultracentrifuge at Cornell University [2]. A failure analysis of the cause of the burst of the centrifuge rotor was reported in Ref. [3]. A failed centrifuge caused by the basket employed for starch production is investigated in Ref. [4]. The geotechnical centrifuge, a large-scale centrifuge, provides physical data for investigating mechanisms of deformation and failure and for validating analytical and numerical methods [5] by simulating mechanisms of deformation and failure and for validating analytical and numerical methods [5] by simulating and studying geotechnical problems. It is currently applied in various geotechnical engineering fields, such as foundation systems, earth structures, offshore systems, earthquakerelated problems, geo-environment studies, etc. [6]. With the continually increasing demand for the geotechnical centrifuge experiments in industrial engineering, the need for a centrifuge facility with advanced capabilities is urgent. One of the most important components of a geotechnical centrifuge is the basket. It is connected to the rotating arm via the pin shaft, which is used to place the inspection model. As a variety of topological configurations of this component are possible, the specific installation and basic shape will be described in detail later (as shown in Figure 1). The basket is designed to withstand complex loads, such as the tensile stress from the arm caused by the rotation, vibration torque, gravity, centrifugal force, and so on. With the high rotational speed, loads make the basket stretch and change in size, which may cause it to fail. In addition, the weight of the basket has a significant influence on the stress of the mounting hole of the arm connected with it. Several failure cases of the centrifuge have been caused by the basket [4,7]; therefore, one of the major tasks of the geotechnical centrifuge design is make the basket as light as possible while satisfying the requirement of its strength. 1. drive unit, 2. speed reducer, 3. coupling, 4. transmission support, 5. arm, 6. pin shaft one, 7. basket  Shape optimization is important in seeking to improve mechanical properties and reduce weight [8][9][10][11][12][13]. In engineering practice, the shape optimization design of the basket is mostly based on the traditional finite element method (FEM), which makes the optimization design process timeconsuming and complicated. There are two reasons for this; one is that the traditional FEM uses conforming mesh discretization and the perturbation method of structural optimization and shape Shape optimization is important in seeking to improve mechanical properties and reduce weight [8][9][10][11][12][13]. In engineering practice, the shape optimization design of the basket is mostly based on the traditional finite element method (FEM), which makes the optimization design process time-consuming and complicated. There are two reasons for this; one is that the traditional FEM uses conforming mesh discretization and the perturbation method of structural optimization and shape sensitivity analysis [14][15][16]. During the optimization process, the mesh should be updated based on shape changes. Unfortunately, conforming the mesh is a huge challenge, especially for complex structures. Statistics show that [17] the conversion of geometric models to finite element models represents 80% of structural analysis work, so the meshing and repartitioning process is very time consuming, especially for the shape optimization of complex structures. Even with some advanced mesh generation methods, such as mesh mapping, Laplacian deformation, adaptive mesh generation, and so on [18], in some cases, mesh deformation caused by the finite element approximation error [19] is still inevitable. The other reason is the difficulty of solving sensitivity in a gradient-based optimization process. Therefore, achieving an optimal basket shape in a geotechnical centrifuge is a significant engineering challenge.
In this work, motivated by the easier remeshing and sensitivity analysis of the Finite Cell Method (FCM) [20,21] compared to traditional FEM, the shape optimization of the basket is investigated and achieved by weighted B-Spline FCM [22][23][24], where the critical component, i.e., the lifting lug, is put into the computational domain divided with the fixed grid, the geometric model is represented by the level set function [25] and R-functions [26], the boundary cells are refined by the Quadtree or Octree method, and a homogeneous Dirichlet boundary condition is applied by the weighted function. Based on the weighted B-Spline FCM, an effective optimization model is built. The structure mass is selected as the optimization objective, and key shape parameters and stress are regarded as design variables and s constraint, respectively. Additionally, the highly-efficient optimization design flow for the lifting lug based on the gradient is developed, and a gradient-based optimization algorithm named the "globally convergent method of moving asymptotes" (GCMMA) [27] is implemented. With four cases containing different topological configurations, the superior design solution of the lifting lug is chosen. The main contributions are as follows.
(1) Extend the B-Spline FCM for structure shape optimization design into engineering applications; (2) Develop a complete shape optimization process of the basket of a geotechnical centrifuge based on the competitive weighted B-Spline FCM and the GCMMA, which settles the light-weighting design of the basket with multi topological configurations whilst also meeting the requirement of the geotechnical centrifuge strength.
The rest of the paper is organized as follows. In Section 2, based on the overall structure of the geotechnical centrifuge, the key component of the basket, i.e., the lifting lug, is investigated. Firstly, the parameterized geometric representation of the lifting lug is established by the level set and Rfunctions. Then, the geometric structure is divided by the fixed grid, the homogeneous Dirichlet boundary is applied exactly by the weighted function, and the force of the lifting lug is analyzed. Finally, for the four topological configurations, the stress analysis is carried out by weighted B-Spline FCM. The stress distribution and mass of each configuration are obtained. Section 3 carries out the shape optimization of the four configurations by the GCMMA. The optimization results are analyzed in Section 4. Some conclusions are given in Section 5.

The Implicit Level Set Function
Different from the classical finite element method, the element meshing is transferred to the grid identifying in the fixed grid. The implicit level set function [25] of a physical domain is usually employed to determine which grids are inside, outside, or cut by the boundary. Suppose a physical domain with a level set function Φ(x), the position of an arbitrary point x with respect to the physical domain can be determined by Equation (1). where Ω, ∂Ω, and D, denote the physical domain, the boundary of the physical domain, and the computational domain, respectively. Taking an arbitrary structure for example, the level set function of the structure is depicted in Figure 2. According to Equation (1), the property of the point x can be identified. where Ω , Ω ∂ , and D , denote the physical domain, the boundary of the physical domain, and the computational domain, respectively. Taking an arbitrary structure for example, the level set function of the structure is depicted in Figure 2. According to Equation (1), the property of the point x can be identified. For any complicated physical domain consisting of basic geometric elements, the level set function and R-functions [26] can be formulated by: in which where 1 Φ and 2 Φ denote the level set functions of the basic geometric elements.

The B-Spline Finite Cell Method
The finite cell method [20,21] is an extension to the core principle and concept using high-order basis functions to approximate the unknown physical domain. Without losing generality, this paper focuses on 3D linear elasticity problems. The original physical domain p W with boundary p W ¶ is embedded into a cubic domain W , as shown in Figure 3. The strong form of the boundary value problem for embedded domain W can be expressed as For any complicated physical domain consisting of basic geometric elements, the level set function and R-functions [26] can be formulated by: where Φ 1 and Φ 2 denote the level set functions of the basic geometric elements.

The B-Spline Finite Cell Method
The finite cell method [20,21] is an extension to the core principle and concept using high-order basis functions to approximate the unknown physical domain. Without losing generality, this paper focuses on 3D linear elasticity problems. The original physical domain Ω p with boundary ∂Ω p is embedded into a cubic domain Ω, as shown in Figure 3. where Ω , Ω ∂ , and D , denote the physical domain, the boundary of the physical domain, and the computational domain, respectively. Taking an arbitrary structure for example, the level set function of the structure is depicted in Figure 2. According to Equation (1), the property of the point x can be identified. For any complicated physical domain consisting of basic geometric elements, the level set function and R-functions [26] can be formulated by: in which where 1 Φ and 2 Φ denote the level set functions of the basic geometric elements.

The B-Spline Finite Cell Method
The finite cell method [20,21] is an extension to the core principle and concept using high-order basis functions to approximate the unknown physical domain. Without losing generality, this paper focuses on 3D linear elasticity problems. The original physical domain p  with boundary p   is embedded into a cubic domain  , as shown in Figure 3. The strong form of the boundary value problem for embedded domain  can be expressed as The strong form of the boundary value problem for embedded domain Ω can be expressed as where σ denotes the stress tensor, u is the displacement field, t expresses the traction vector defining on Γ N , g is the prescribed displacement vector defining Dirichelt boundary condition (DBC) on Γ D , f denotes the body force on Ω p , and n is the unit outward normal vector on Γ N . The scalar factor β can be defined as The weak form of Equation (4) can be stated as with the bilinear form being and the load function being in which H 1 (Ω) expresses the Sobolev space of order 1. v denotes the admissible test displacement in kinematics. ε and D are the strain tensor and the constitutive elasticity tensor of Ω p , respectively. Similar to the finite element method, the embedded domain Ω is discretized by a general Cartesian grid. As illustrated in Figure 4a, cells are simply hexahedrons, and can be classified into physical cells, fictitious cells, and boundary cells, according to the implicit level set function of Ω p .
Appl. Sci. 2020, 10, 620 where s denotes the stress tensor, u is the displacement field, t expresses the traction vector defining on N G , g is the prescribed displacement vector defining Dirichelt boundary condition (DBC) on D G , f denotes the body force on p W , and n is the unit outward normal vector on N G . The scalar factor b can be defined as The weak form of Equation (4) can be stated as with the bilinear form being and the load function being  In the B-spline FCM [22][23][24], the B-spline basis functions are used as the interpolation functions over the hexahedral cells. The standard parametric domain (see Figure 4b) is mapped to the discretization of the embedded domain. For 3D problem mapping, it can be easily constructed by taking the tensor product of corresponding univariate B-spline basis functions.
In the B-spline FCM [22][23][24], the B-spline basis functions are used as the interpolation functions over the hexahedral cells. The standard parametric domain (see Figure 4b) is mapped to the discretization of the embedded domain. For 3D problem mapping, it can be easily constructed by taking the tensor product of corresponding univariate B-spline basis functions.
, , , , ,..., k r g g g + + = Y , respectively. n , m , and k p , q , and r denote the polynomial degree, in this paper 2 p q r = = = . The following recurrence relation (10) is included by the above expression.
where v is spanned by the B-spline basis function M . A linear combination of shape functions is considered concerning the numerical approximation h u . Let us assume an arbitrary point P = x, y, z T inside Ω; this can be interpolated by the following mapping relation: where P s = x s , y s , z s T refers to the sth control point, as illustrated in Figure 4c.
The trial functions u h and the test functions v h are defined over R h , which is a finite-dimensional subspace. The discrete form of the weak form Equation (6) can be stated as where v is spanned by the B-spline basis function M. A linear combination of shape functions is considered concerning the numerical approximation u h .
Appl. Sci. 2020, 10, 620 in which U denotes the unknown displacement vector of control points. Based on Ritz-Galerkin approach, inserting Equation (13) into Equation (12), the finite cell formulation, which is identical to the FEM in form, is obtained as where K and F are the global stiffness matrix and the global load vector that can be assembled in the traditional fashion from all cell stiffness matrices K c and load vectors F c , respectively. K c and F c are acquired by (17) in which Ω c illustrates one cell domain. Γ Nc is the Neumann boundary related to one cell. B is the strain-displacement matrix, written as

The Weighted Dirichlet Boundary Conditions
Comparing the FCM with the classical FEM, the difficulty of generating meshes for complex structures shifts to the problem of implementing boundary conditions. In particular, when the boundary ∂Ω p does not conform to the cells, this problem arises. At present, the weighted method is usually employed to implement the Dirichlet boundary condition.
In this paper, the weighted B-spline is used to penalize locally the interpolation displacement field defined by Equation (13). For homogeneous Dirichlet boundary conditions, the prescribed displacement g applied on Γ D is equal to zero. With the weighted function w(x), an alternative Equation (13) can be stated as Therefore, the geometry equation can be addressed as Appl. Sci. 2020, 10, 620 and Thus, Equations (16) and (17) are modified as

Geometric Representation Based on Level-Set Function
The basket, as one of the most important components in a geotechnical centrifuge, is used to contain the inspection model. Figure 1a demonstrates the structure of a typical geotechnical centrifuge, where the basket consists of the lifting lug, mounting base plate, and a pin shaft, and is at the farthest end of the rotation radius of the geotechnical centrifuge. As shown in Figure 1b, the shape of lifting lug may be trapezoidal, but its contour line is not necessarily a simple straight line. Based on the different strength requirements, it can be designed as circular arcs, quadratic curves, and so on. The material of the lifting lug is usually 34CrNi1Mo; thus, the upper bound of the von-Mises stress with safety factor 2 is 400 MPa, and Young's modulus and Poisson's ratio are 207.5 GPa and 0.3, respectively.
In order to carry out the shape optimization of the lifting lug based on the weighted B-Spline FCM, the level set function [25], as a competitive parametric modeling method, is adopted in this paper. Based on the statement described in Section 2.1, the level set function Φ of the lifting lug established by Boolean operations is derived as follows.
in which Φ i (i = 1, 2, . . . , 12) denotes the level set function of Ω i in Figure 6, and the design variables are included by some of the level set function Φ i (i = 1, 2, . . . , 12). With the changing of design variables, the value of Φ changes, which leads to a change in the cells properties in the fixed grid. Therefore, shape optimization can be carried out.
FCM, the level set function [25], as a competitive parametric modeling method, is adopted in this paper. Based on the statement described in Section 2.1, the level set function Φ of the lifting lug established by Boolean operations is derived as follows.
in which ( 1, 2,...,12) i Φ i = denotes the level set function of i Ω in Figure 6, and the design variables are included by some of the level set function ( 1, 2,...,12) With the changing of design variables, the value of Φ changes, which leads to a change in the cells properties in the fixed grid. Therefore, shape optimization can be carried out. Figure 6. The lifting lug model constructed by Boolean operations.

Grid Identification
Based on the zero level set function of the lifting lug, the grids are classified into physical cells, boundary cells, and fictitious cells, as described in Figure 7. In order to improve the accuracy of the analysis, the boundary cells are refined by the Octree technique.

Topological Configurations and Related Design Variables
In order to obtain a superior shape solution, four topological configurations forthe lifting lug are taken into account in this paper, as illustrated in Figure 8.

Grid Identification
Based on the zero level set function of the lifting lug, the grids are classified into physical cells, boundary cells, and fictitious cells, as described in Figure 7. In order to improve the accuracy of the analysis, the boundary cells are refined by the Octree technique.
FCM, the level set function [25], as a competitive parametric modeling method, is adopted in this paper. Based on the statement described in Section 2.1, the level set function Φ of the lifting lug established by Boolean operations is derived as follows.
in which ( 1, 2,...,12) i Φ i = denotes the level set function of i Ω in Figure 6, and the design variables are included by some of the level set function ( 1, 2,...,12) With the changing of design variables, the value of Φ changes, which leads to a change in the cells properties in the fixed grid. Therefore, shape optimization can be carried out. Figure 6. The lifting lug model constructed by Boolean operations.

Grid Identification
Based on the zero level set function of the lifting lug, the grids are classified into physical cells, boundary cells, and fictitious cells, as described in Figure 7. In order to improve the accuracy of the analysis, the boundary cells are refined by the Octree technique.

Topological Configurations and Related Design Variables
In order to obtain a superior shape solution, four topological configurations forthe lifting lug are taken into account in this paper, as illustrated in Figure 8.

Topological Configurations and Related Design Variables
In order to obtain a superior shape solution, four topological configurations forthe lifting lug are taken into account in this paper, as illustrated in Figure 8. Case 1. Original single web. The single web is widely used in the lifting lug structure, supposing that the outer contour has eight design variables.
Case 2. Porous dual web. In order to reduce the weight of the lifting lug, a topological opening is carried out inside the web, where the shape of the holes is an ellipse, so each hole contains two design variables. Therefore, the structure contains a total of 17 design variables.
Case 3. Open deep groove dual web. On the basis of the initial single web, the lifting lug is grooved to achieve the purpose of weight reduction. The number of design variables is 12.

Mechanical Analysis
The lifting lug used to connect the rotating arm and the mounting base plate is taken as the research object, and the loads are analyzed in combination with the working environment. Generally, the load type mainly includes two parts, i.e., the centrifugal force caused by its self-mass, and the tensile stress generated by the centrifugal force of other connected parts. The centrifugal force C F is represented by Equation (30).  Case 1. Original single web. The single web is widely used in the lifting lug structure, supposing that the outer contour has eight design variables.
Case 2. Porous dual web. In order to reduce the weight of the lifting lug, a topological opening is carried out inside the web, where the shape of the holes is an ellipse, so each hole contains two design variables. Therefore, the structure contains a total of 17 design variables.
Case 3. Open deep groove dual web. On the basis of the initial single web, the lifting lug is grooved to achieve the purpose of weight reduction. The number of design variables is 12.

Mechanical Analysis
The lifting lug used to connect the rotating arm and the mounting base plate is taken as the research object, and the loads are analyzed in combination with the working environment. Generally, the load type mainly includes two parts, i.e., the centrifugal force caused by its self-mass, and the tensile stress generated by the centrifugal force of other connected parts. The centrifugal force F C is represented by Equation (30).
where m is the mass of the lifting lug, r is the radius of rotation of the lifting lug, and ω represents the rotary angular velocity. With the high-speed rotation of the centrifuge, the stress σ a generated by the centrifugal load of the mounting base plate and the inspection model of centrifugal test at the joint and the centrifugal force F C act simultaneously, and the combination of the two causes significant strength problems for the lifting lug. Assuming that the load caused by the centrifugal loads of the mounting base plate and the inspection model of centrifugal test are evenly distributed over the edge of the lifting lug connection, the radial stress of the load is defined by where F represents the centrifugal load of the mounting base plate and the inspection model of centrifugal test, R is the rotary radius of the lifting lug edge, i.e., the distance from the position where the mounting base plate is attached to the lifting lug to the center axis of the centrifuge, L represents the width of the edge of the lifting lug. Based on the loads and boundary conditions stated above, the weighted B-spline FCM is employed to compute the von-Mises stress distribution of the lifting lug in fixed grid. Figure 9 shows the initial von-Mises stress distribution of the four cases. A comparison of the maximum von-Mises stress and mass among the four cases is demonstrated, as described in Figure 10. From these figures, it can be noted that the maximum von-Mises stresses of the four cases appear at the junction of the lifting lug and the mounting base plate, wherein both of the maximum von-Mises stress and mass in the case 1 are the largest. The four cases with these topological configurations all satisfy the stress condition with a safety factor greater than or equal to 2. Thus, in order to obtain a superior solution of the lifting lug based on the four topological configurations, shape optimization is carried out for the four cases.
Appl. Sci. 2020, 10, 620 11 of 19 where the mounting base plate is attached to the lifting lug to the center axis of the centrifuge, L represents the width of the edge of the lifting lug. Based on the loads and boundary conditions stated above, the weighted B-spline FCM is employed to compute the von-Mises stress distribution of the lifting lug in fixed grid. Figure 9 shows the initial von-Mises stress distribution of the four cases. A comparison of the maximum von-Mises stress and mass among the four cases is demonstrated, as described in Figure 10. From these figures, it can be noted that the maximum von-Mises stresses of the four cases appear at the junction of the lifting lug and the mounting base plate, wherein both of the maximum von-Mises stress and mass in the case 1 are the largest. The four cases with these topological configurations all satisfy the stress condition with a safety factor greater than or equal to 2. Thus, in order to obtain a superior solution of the lifting lug based on the four topological configurations, shape optimization is carried out for the four cases.

Shape Optimization Model
In order to solve the problem of designing a light-weight basket which meets the requirement of the geotechnical centrifuge in terms of strength, shape optimization is implemented. The optimization model is constructed by where θ stands for the vector of design variables characterizing geometric shapes of the lifting lug. The objective function is to minimize the mass of the physical domain Ω . In the equilibrium equation = KU F , K is the stiffness matrix, U denotes the displacement vector, and F represents the load vector. The constraints concern the maximum von-Mises stress in the lifting lug being limited by the allowable value σ . q is the number of nodes of the finite cells.

The Flowchart of Shape Optimization
The whole optimization solution is based on the weighted B-Spline FCM [23,24] and the GCMMA [27,28]. The weighted B-Spline FCM analysis process is as follows: including geometric modeling based on level set function and R-functions, determining the properties of all the cells in fixed grid, refining the boundary cells with the Quadtree or Octree technique, imposing the homogeneous Dirichlet boundary condition exactly with the weighted function, and calculating the shape sensitivity. The optimization solution is implemented by the GCMMA algorithm integrated into Boss-Quattro optimization platform [29]. A flowchart of the shape optimization based on the weighted B-Spline FCM and GCMMA is presented in Figure 11.

Shape Optimization Model
In order to solve the problem of designing a light-weight basket which meets the requirement of the geotechnical centrifuge in terms of strength, shape optimization is implemented. The optimization model is constructed by where θ stands for the vector of design variables characterizing geometric shapes of the lifting lug. The objective function is to minimize the mass of the physical domain Ω. In the equilibrium equation KU = F, K is the stiffness matrix, U denotes the displacement vector, and F represents the load vector. The constraints concern the maximum von-Mises stress in the lifting lug being limited by the allowable value σ. q is the number of nodes of the finite cells. l d and l y are numbers of design variables d i and y j , respectively. y i and y i denote the lower and upper bounds of the ith design variable y i , respectively. d j and d j denote the lower and upper bounds of the jth design variable d j , respectively.

The Flowchart of Shape Optimization
The whole optimization solution is based on the weighted B-Spline FCM [23,24] and the GCMMA [27,28]. The weighted B-Spline FCM analysis process is as follows: including geometric modeling based on level set function and R-functions, determining the properties of all the cells in fixed grid, refining the boundary cells with the Quadtree or Octree technique, imposing the homogeneous Dirichlet boundary condition exactly with the weighted function, and calculating the shape sensitivity. The optimization solution is implemented by the GCMMA algorithm integrated into Boss-Quattro optimization platform [29]. A flowchart of the shape optimization based on the weighted B-Spline FCM and GCMMA is presented in Figure 11.

. Sensitivity Analysis
In shape optimization based on Equation (32), sensitivity analysis, i.e., stress sensitivity an ass sensitivity with respect to design variables, is essential for gradient-based optimizatio orithms.
For the level set function ( , ) Φ x d of a physical domain Ω , the Dirac delta function ˆ( , δ x d notes the directional derivative of the Heaviside function H in the normal direction, i.e., Then, the mass sensitivity is calculated by Figure 11. Flowchart of shape optimization based on the weighted B-Spline FCM and GCMMA.

Sensitivity Analysis
In shape optimization based on Equation (32), sensitivity analysis, i.e., stress sensitivity and mass sensitivity with respect to design variables, is essential for gradient-based optimization algorithms.
For the level set function Φ(x, d) of a physical domain Ω, the Dirac delta functionδ(x, d) denotes the directional derivative of the Heaviside function H in the normal direction, i.e., Then, the mass sensitivity is calculated by in which Γ represents the boundary of the physical domain, i.e., Γ =∂Ω. Based on the stress-strain relationship σ = Dε = DBU, the stress sensitivity is estimated by where D represents the elastic modulus, U denotes the displacement field, and B is the strain-displacement matrix corrected by the weighted function.
For the four topological configurations of the lifting lug are considered in Section 2.3. Since the design variables are all on the free boundary, B is independent of the design variables. Thus, the stress sensitivity is derived as According to the equilibrium equation KU = F, the stress sensitivity is deduced as and

Optimization Results
Based on the weighted B-Spline FCM performed in Matlab script and GCMMA implemented in the Boss-Quattro optimization platform, the shape optimization results for the four cases are obtained. The final design solutions of the four cases provide us with the minimum mass that removes material to the maximum extent. Meanwhile, the allowable von-Mises stresses are all still satisfied. The final design solutions of the four cases of the lifting lug and the von-Mises stress distributions are depicted in Figure 12. The corresponding iteration histories of the von-Mises stress and the mass of the four cases are shown in Figure 13. The corresponding iteration histories of the von-Mises stress and the mass of the four cases are shown in Figure 13.   Table 1.  In Table 1 and Figures 12-14, it can be seen that in contrast to the initial configurations, although the maximum von-Mises stresses of the four cases all increased, they still meet the strength requirements of the lifting lug, i.e., max 400MPa σ ≤ . The mass reduction by shape optimization in the four cases is very obvious. Below, the four cases are discussed in detail. Case 1: Compared with the initial configuration, the maximum stress change is very small; however, by the optimization design for the outer contour of the lifting lug, the final mass was   Table 1.    Table 1.  In Table 1 and Figures 12-14, it can be seen that in contrast to the initial configurations, although the maximum von-Mises stresses of the four cases all increased, they still meet the strength requirements of the lifting lug, i.e., max 400MPa σ ≤ . The mass reduction by shape optimization in the four cases is very obvious. Below, the four cases are discussed in detail. Case 1: Compared with the initial configuration, the maximum stress change is very small; however, by the optimization design for the outer contour of the lifting lug, the final mass was  In Table 1 and Figures 12-14, it can be seen that in contrast to the initial configurations, although the maximum von-Mises stresses of the four cases all increased, they still meet the strength requirements of the lifting lug, i.e., σ max ≤ 400MPa. The mass reduction by shape optimization in the four cases is very obvious. Below, the four cases are discussed in detail.
Case 1: Compared with the initial configuration, the maximum stress change is very small; however, by the optimization design for the outer contour of the lifting lug, the final mass was reduced by 72.31 kg, i.e., 7.3%, which is still the largest among the four cases. This indicates that although the original single web configuration can achieve the goal of weight reduction by shape optimization on the premise of meeting the strength requirements, the effect of weight reduction is limited.
Case 2: With the maximum stress increasing by 9 MPa based on the shape optimization, the mass decreases by 90.1 kg from 909.86 kg to 819.76 kg, i.e., by 10%. Compared with the optimized single web, the weight loss of the porous dual web is about 96 kg, which shows the configuration of case 2 is better than that of case 1.
Case 3: The maximum stress in the optimized configuration does not change much, but as seen in Figure 12, the stress distribution changes greatly. The maximum stress in the initial configuration occurs at the connection hole with the mounting base plate, while that in the final configuration occurs in multiple places. Additionally, the final mass declines by 94.09 kg based on shape optimization, i.e., by 10%. In contrast to the optimized single web in case 1, the final open deep groove dual web in case 3 reduces 87.69 kg, which is less than that of case 2.
Case 4: Compared with the initial configuration, the maximum stress and stress distribution in the optimized configuration do not change substantially, but the final mass drops by 92.75 kg, i.e., by 10.6%. Compared with the optimized single web, the weight loss of the connected closed dual web is 133.38 kg, which indicates that the final configuration of case 4 is the best among the four cases.
The above results indicate that while satisfying the stress requirement, the final mass of case 4 is the least. Therefore, considering the matching, cost, and other factors, the optimized shape configuration of the case 4 is superior to those of the other three cases.

Conclusions
The lifting lug, as one key structure of the geotechnical centrifuge, is the focus of this research. This paper presents an effective and systematic structure design and shape optimization for the lifting lug with multi topological configurations based on the weighted B-Spline FCM and GCMMA. It provides a complete shape optimization process for the lifting lug of the geotechnical centrifuge, which addresses the problem of the light-weight design of the lifting lug while satisfying the requirement of the geotechnical centrifuge strength. Four topological configurations, i.e., original single web, porous dual web, open deep groove dual web, and connected closed dual web, are investigated and optimized. Two main conclusions can be drawn.
First, through shape optimization design based on the weighted B-Spline FCM and GCMMA, the final mass of the four configurations is greatly reduced compared with the initial configurations, where the greatest weight loss was observed in case 4 is 10.6%. This indicates that the weighted B-Spline FCM and GCMMA can be applied for shape optimization of the complex structure in the engineering design.
Second, in contrast to the final single web adopted in the traditional basket design in case 1, the final configuration in case 4, i.e., connected closed dual web, has the least mass. The final mass was reduced by 133.38 kg with the centrifuge strength requirement being satisfied. Therefore, the final configuration in case 4 is superior to other three configurations.
In this work, the design procedure extending the weighted B-Spline FCM and GCMMA into shape optimization is developed for the rotational structure concerning the important role it plays in engineering applications. However, our mechanical analysis was restricted to the static loads; vibrating loads, also as a common load condition during the operation of centrifuges, were not tested. Further investigation on this problem will make the design of the centrifuge within the FCM and GCMMA framework more comprehensive.

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