Loading [MathJax]/jax/output/SVG/jax.js
Research article

A parameterized level set method for structural topology optimization based on reaction diffusion equation and fuzzy PID control algorithm

  • Received: 17 November 2021 Revised: 22 April 2022 Accepted: 05 May 2022 Published: 12 May 2022
  • We propose a parameterized level set method (PLSM) for structural topology optimization based on reaction diffusion equation (RDE) and fuzzy PID control algorithm. By using the proposed method, the structural compliance minimization problem under volume constraints is studied. In this work, the RDE is used as the evolution equation of level set function, and the topological derivative of the material domain is used as the reaction term of the RDE to drive the evolution of level set function, which has little dependence on the initial design domain, and can generate holes in the material domain; the compactly supported radial basis function (CS-RBF) is used to interpolate the level set function and modify the RDE, which can improve the computational efficiency, and keep the boundary smooth in the optimization process. Meanwhile, the fuzzy PID control algorithm is used to deal with the volume constraints, so that the convergence process of the structure volume is relatively stable. Furthermore, the proposed method is applied to 3D structural topology optimization. Several typical numerical examples are provided to demonstrate the feasibility and effectiveness of this method.

    Citation: Mingtao Cui, Min Pan, Jie Wang, Pengjie Li. A parameterized level set method for structural topology optimization based on reaction diffusion equation and fuzzy PID control algorithm[J]. Electronic Research Archive, 2022, 30(7): 2568-2599. doi: 10.3934/era.2022132

    Related Papers:

    [1] Hualong Du, Qiuyu Cui, Pengfei Liu, Xin Ma, He Wang . Reformative artificial bee colony algorithm based PID controller for radar servo system. Electronic Research Archive, 2022, 30(8): 2941-2963. doi: 10.3934/era.2022149
    [2] Wenbin Zheng, Jinjin Li, Shujiao Liao . Decision self-information based on parameterized fuzzy β neighborhood and its application in three-way multi-attribute group decision-making. Electronic Research Archive, 2022, 30(12): 4553-4573. doi: 10.3934/era.2022231
    [3] Meng-Xue Chang, Bang-Sheng Han, Xiao-Ming Fan . Global dynamics of the solution for a bistable reaction diffusion equation with nonlocal effect. Electronic Research Archive, 2021, 29(5): 3017-3030. doi: 10.3934/era.2021024
    [4] Wen Li, Jian Wang, Zhanpeng Du, Hongfeng Ma, Lijie Zhang, Libin Duan . Lightweight design method and application of MEWP bracket based on multi-level optimization. Electronic Research Archive, 2022, 30(12): 4416-4435. doi: 10.3934/era.2022224
    [5] Sida Lin, Lixia Meng, Jinlong Yuan, Changzhi Wu, An Li, Chongyang Liu, Jun Xie . Sequential adaptive switching time optimization technique for maximum hands-off control problems. Electronic Research Archive, 2024, 32(4): 2229-2250. doi: 10.3934/era.2024101
    [6] Peng Gao, Pengyu Chen . Blowup and MLUH stability of time-space fractional reaction-diffusion equations. Electronic Research Archive, 2022, 30(9): 3351-3361. doi: 10.3934/era.2022170
    [7] Mingtao Cui, Wennan Cui, Wang Li, Xiaobo Wang . A polygonal topology optimization method based on the alternating active-phase algorithm. Electronic Research Archive, 2024, 32(2): 1191-1226. doi: 10.3934/era.2024057
    [8] Ziqiang Wang, Chunyu Cen, Junying Cao . Topological optimization algorithm for mechanical-electrical coupling of periodic composite materials. Electronic Research Archive, 2023, 31(5): 2689-2707. doi: 10.3934/era.2023136
    [9] Yuhai Zhong, Huashan Feng, Hongbo Wang, Runxiao Wang, Weiwei Yu . A bionic topology optimization method with an additional displacement constraint. Electronic Research Archive, 2023, 31(2): 754-769. doi: 10.3934/era.2023037
    [10] Yineng Ouyang, Zhaotao Liang, Zhihui Ma, Lei Wang, Zhaohua Gong, Jun Xie, Kuikui Gao . A class of constrained optimal control problems arising in an immunotherapy cancer remission process. Electronic Research Archive, 2024, 32(10): 5868-5888. doi: 10.3934/era.2024271
  • We propose a parameterized level set method (PLSM) for structural topology optimization based on reaction diffusion equation (RDE) and fuzzy PID control algorithm. By using the proposed method, the structural compliance minimization problem under volume constraints is studied. In this work, the RDE is used as the evolution equation of level set function, and the topological derivative of the material domain is used as the reaction term of the RDE to drive the evolution of level set function, which has little dependence on the initial design domain, and can generate holes in the material domain; the compactly supported radial basis function (CS-RBF) is used to interpolate the level set function and modify the RDE, which can improve the computational efficiency, and keep the boundary smooth in the optimization process. Meanwhile, the fuzzy PID control algorithm is used to deal with the volume constraints, so that the convergence process of the structure volume is relatively stable. Furthermore, the proposed method is applied to 3D structural topology optimization. Several typical numerical examples are provided to demonstrate the feasibility and effectiveness of this method.



    Structural topology optimization can achieve the pursuit of saving materials or improving the performance of structures with the given design domain and the relevant boundary conditions. Hence it has already been a powerful tool for the design of engineering structure. As an important branch of structural topology optimization, continuum topology optimization has attracted more and more attention. With the development of computational mechanics, there has been significant progress made in the theory and methodology of continuum topology optimization over the past 30 years. During this period, much pioneering research has been published, advancing related theories and techniques, and spurring the use of practical topology optimization methods. Accordingly, a variety of methods and their accompanying technical schemes have been proposed to deal with continuum topology optimization and obtain effective solutions. The methods for continuum topology optimization can be commonly classified into the following categories: the homogenization method (Bendsøe and Kikuchi [29]; Bendsøe [31]), the SIMP method (Bendsøe and Sigmund [30]), the evolutionary structural optimization (ESO) method (Xie and Steven [55]), the phase field method (Bourdin and Chambolle [1]; Marino et al. [27]), the level set method (Allaire et al. [7]; Wang et al. [35]), the feature driven optimization approach, alternatively named as the moving morphable voids approach (Guo et al. [53]; Zhou et al. [56]; Zhang et al. [51]; Zhang et al. [52]; Zhu et al. [20]), and so on.

    The level set method was first proposed by Osher and Sethian [45] to simulate moving boundaries. Initially, the level set method was mainly used to simulate the evolution of interface in multiphase flow (Sussman et al. [32]; Sethian and Smereka [16]) and image processing technology (Malladi et al. [43]; Osher and Paragios [46]). Sethian and Wiegmann [15] used the immersion interface method to calculate the structural stress, and then used the level set method to change the shape of the structure. Osher and Santosa [44] employed the level set method to solve the structural topology optimization problem of the vibration system.

    Virtual time t is introduced into the traditional level set method to make it more effective for topology optimization. By calculating the derivative of level set function relative to virtual time t (van Dijk et al. [38]), the Hamilton-Jacobi partial differential equation (Peng et al. [3]) (PDE) related to virtual time t is established. The upwind method can be employed to solve PDE, thereby updating the level set equation and realizing the structural topology optimization. However, when the upwind method is used to solve PDE for topology optimization, the computational efficiency is not high because the time step is limited by the Courant-Friedrichs-Lewy (CFL) condition. Moreover, the holes in the structure cannot be generated, so the evolution of the structure boundary depends on the initial conditions.

    To improve the traditional level set method, many researches have been carried out. Xia et al. [42] proposed a semi-Lagrangian method to update the level set function to solve the problem that the time step is limited by the CFL condition when using the upwind method, and compared the method with the upwind method to verify the effectiveness of the method. Zhou and Wang [37] used the line search algorithm to determine the time step adaptively, which improves the oscillation problem of the objective function and volume constrained optimization curve when using the semi-Lagrangian method. Eschenauer et al. [8] proposed a "bubble" method, which inserts a "bubble" into the structure, and uses the positioning criteria and shape optimization method to determine the optimal position and shape of the "bubble". This "bubble" method for structural topology optimization can generate holes in the structure. Sokolowski and Zochowski [19] used topological derivative to calculate the change of objective function when adding or deleting materials in the structure, and gave the topological derivative form of plane elastic system. Burger et al. [24] combined the topological derivative with the level set method to calculate the sensitivity of the objective function when the structure topology changes.

    To further improve the level set method, Allaire et al. [6] used the adjoint method to calculate the shape derivative, and took the shape derivative as the normal velocity of the free boundary in the optimization process. At the same time, in order to control that the level set function near the boundary is neither too steep nor too flat, the operation of reinitializing the level set function periodically was introduced to ensure that the structure can converge to the best shape. Luo et al. [57] applied the optimization criterion method (Zhou and Rozvany [36]) to the parametric level set method to replace the steepest descent method to update the design variables. Allaire et al. [5] combined the topological derivative with shape sensitivity, thereby establishing the algorithm model of alternating descent of shape sensitivity and topological derivative, which deals with the boundary evolution problem and meanwhile can introduce holes in the structure. Wei and Wang [40] combined the level set method with the phase field method, and implemented the piecewise constant level set method by using the numerical scheme of the additive operator splitting in the phase field method, which describes the evolution of the structural boundary. Yamada et al. [49] proposed a new method to construct the reaction diffusion equation (RDE) by using the virtual energy term in the phase field method, and gave the derivation process of RDE. Choi et al. [18] took the design sensitivity of topology optimization as the reaction term of RDE, which can generate new holes in the structure without using the topological derivative. Otomori et al. [28] provided a MATLAB implementation of level set method based on RDE, and gave the derivation process of topological derivative.

    The method of topology description function (TDF) was proposed by de Ruiter and van Keulen [26], which utilizes a set of basis function to describe the geometric structure. Through adjusting the shape, parameters and distribution of the basis function, the structure can be described reasonably by the method of TDF. Thereafter, the TDF method has evolved into the parameterized level set method (Cecil et al. [48]; Wang and Wang [47]; Wei et al. [41]; Cui et al. [33]) (PLSM). The PLSM retains the advantage of implicit representation of evolution surface by level set function, thereby transforming topology optimization problem into parameter optimization problem, and meanwhile does not need to solve complex Hamilton-Jacobi PDE directly or carry on the re-initialization operation. Thus good topology optimization effect can be achieved by using the PLSM. Moreover, radial basis function (RBF) and level set equation are used to represent the shape of the structure implicitly in the PLSM. Therefore, the selection of the basis function determines the shape type and detail of the structure boundary described by the level set equation.

    According to the supporting size, there are two types of basic functions: compactly supported radial basis function (CS-RBF) and global supported radial basis function (GS-RBF). Wang and Wang [47] used the multi-quadric (MQ) function to establish the implicit RBF model interpolating the level set function, transformed Hamilton-Jacobi PDE into ordinary differential equation (ODE), and solved it by finite element method. Luo et al. [58] proposed an element-free Galerkin level set method, where CS-RBF is used to construct the element-free shape function and the structural boundary is updated to achieve structural optimization by solving the discrete level set function in time. Ho et al. [12] proposed a PLSM based on mobile node by using inverse multi-quadric (IMQ) function, which takes the position of RBF node as the design variable and uses the steepest descent method to push the node to a new position. Wei et al. [41] provided a MATLAB code implementation of level set method based on MQ function, and extended it for different boundary conditions and different types of RBF. Cui et al. [33] proposed a PLSM combined with CS-RBF and the method of moving asymptotes (Svanberg [22]) (MMA) to update the design variables, and introduced shape sensitivity constraint factor in the iteration process of MMA algorithm, which makes the iteration step controllable and speeds up the optimization process.

    Three-dimensional (3D) structural topology optimization (Yang et al. [54]; Sigmund and Clausen [39]) can determine the geometric configuration of engineering structure and improve structural performance for engineering application. Some methods for two-dimensional (2D) structural topology optimization have thus far been successfully extended to 3D problem, which broadens the application scope of structural topology optimization. Du and Olhoff [17] studied the topology optimization of 3D structure under design related loads, and proposed an algorithm for generating loading surface of 3D structure based on iso-parametric surface theory. Abolbashari and Keshavarzmanesh [25] studied the influence of removal rate, evolution rate and element size on the optimization of 2D and 3D continuum structures by the ESO method. Zhang et al. [14] proposed a boundary search algorithm, which can effectively identify the load surface and deal with 3D topology optimization problem under design related loads. Muñoz et al. [2] used Cartesian mesh finite element method instead of standard finite element method in topology optimization of 3D structure, and combined h-adaptive mesh optimization method in optimization algorithm to improve the accuracy of boundary definition. Liu and Tovar [21] extended the SIMP method to the 3D field, and gave 169 lines of MATLAB code including sensitivity analysis, density filtering and finite element analysis. On this basis, Ferrari and Sigmund [4] improved the above-mentioned 169 lines of MATLAB code through more effective assembly of stiffness matrix to accelerate the realization of filter function, and obtained more concise and effective 3D MATLAB code. Zegard and Paulino [50] proposed a 3D topology optimization algorithm based on ground structure, which can well deal with the structural optimization problems of concave and internal holes, and gave the corresponding MATLAB code. Wang et al. [34] combined the B-spline function and the level set function, divided the whole design domain into a group of 3D subdomain cells, and used B-spline function to update design variables recursively, so as to realize the 3D topology optimization of cellular structure with spatial changes. Recently, Li et al. [9,10,11] have developed the RDE method and extended its applications in large-scale 3D topology optimization of multi-physics problems, i.e., mean compliance problem, fluid-structure interaction problem, and thermal fluid problem, by using adaptive mesh refinement or body-fitted mesh adaption.

    When the traditional level set method is used for structural topology optimization, the upwind method is often used to solve Hamilton-Jacobi PDE, and the conventional level set function is updated according to the solution of the PDE. In the optimization iterations, the time step is required to meet the CFL condition. Thus the size of the time step is limited, which reduces the computational efficiency. Moreover, the re-initialization operation needs to be performed in each of the optimization iterations to control the gradient of the level set function near the boundary, so the stability of numerical calculation is affected. However, when RDE is used to replace Hamilton-Jacobi PDE update the level set function in structural topology optimization, it allows not only shape but also topological changes, and allows the geometrical complexity of the optimal configuration to be specified. Therefore, the RDE-based level set method has significant advantages over the conventional level set method for structural topology optimization. In this work, a PLSM based on the RDE and the fuzzy adaptive PID control algorithm is proposed to address structural topology optimization of minimizing compliance with volume constraints. The RDE is used as the evolution equation of level set function, and the regularization parameters in RDE are used to adjust and control the geometry complexity of the optimized structure. The topological derivative is not only used to calculate the sensitivity of the objective function when the hole is inserted at any position of the structure, but also used as the reaction term of the RDE to drive the evolution of the level set function. In the optimization iteration process, CS-RBF is used to interpolate the level set function to ensure the smoothness of the approximate surface. Furthermore, the fuzzy PID control algorithm is employed to address the volume constraint, so as to makes the convergence process of structure volume more stable. The proposed method is also applied to 3D structural topology optimization, so as to confirm that this method can enhance the optimization efficiency and the stability of the optimization process on the premise of ensuring the accuracy of the optimization results. The proposed method retains the merits of the PLSM, the RDE and the fuzzy adaptive PID control algorithm. Therefore, it can obtain the reasonable results with relatively high speed and good stability for structural topology optimization

    The rest of this paper is organized as follows. In Section 2, the level set method and the corresponding solution process which combines the compactly supported radial basis function and fuzzy PID control algorithm are introduced. In Section 3, the mathematical model of the topology optimization problem which takes the minimum compliance as the objective function and the volume fraction as the constraint condition is established, and the role of the initial PID control parameters is studied. In Section 4, 2D and 3D structural optimization examples are used to verify the feasibility and effectiveness of the proposed method. Finally, the conclusions are drawn in Section 5.

    The level set function uses its zero-value curve to describe the structure boundary. By solving the value of level set function at each time, the evolution equation of the structure boundary is obtained, and the shape and topology of the structure in the material domain are clearly described. By taking a 2D structure as an example, the level set function is used to represent the material distribution of the structure as shown in Figure 1.

    Figure 1.  Relationship between level set equation and material distribution in design domain.

    The relationship between level set function ϕ and material distribution in design domain is expressed by the mathematical formula as follows.

    {ϕ(x) > 0xΩϕ(x)=0xΓϕ(x) < 0xDΩ (1)

    where, D represents the design domain of the structure, Ω represents the solid part of the structure, Γ represents the boundary of the structure, and x represents element node after the finite element mesh is divided. ϕ(x) > 0 represents the material domain, ϕ(x)=0 represents the structural boundary, and ϕ(x) < 0 represents the void.

    A topology optimization problem with minimum compliance as objective function and volume fraction as constraint condition is considered. As shown in Figure 2, structural topology optimization problem generally involves prescribed design domain and boundary conditions. Design domain D usually includes material domain Ω and void region. With respect to boundary conditions, displacement constraints are usually imposed at the boundary Γu, external loads F are often applied at the boundary Γt, and body force b is applied on the overall material domain. Under the volume constraint, the mathematical expression of the objective function is written as follows.

    minΩJ(Ω)=Ωj(x)dΩs.t.G(Ω)=ΩdΩVmax0 (2)
    Figure 2.  Design domain and boundary conditions for structural topology optimization.

    where, J(Ω) is the objective function, j(x) is the function in the structural design domain, G(Ω) is the volume constraint function, Vmax is the set volume constraint value, and x represents element node after the structural design domain is divided into finite element meshes.

    To introduce the virtual energy term of phase field model into the objective function and regularize the structural topology optimization problem, it needs to be assumed that the level set function ϕ and the phase field variables in the phase field method have the same distribution properties. Based on this assumption, the upper and lower bound constraints are imposed on the level set function ϕ, which stipulates that the value of level set function must be equal to 1 or −1 in the area far enough from the structure boundary, so that the smoothing effect only works on the points near the structure boundary. Therefore, the definition of level set function expressed in Eq (1) is modified, and the revised definition of level set function is shown in Eq (3).

    {0ϕ(x)1xΩϕ(x)=0xΓ1ϕ(x)0xDΩ (3)

    By introducing the virtual energy term in the phase field method into the objective function, the topology optimization problem described by Eq (2) is transformed into the topology optimization problem shown in Eq (4).

    minϕJR(χΩ(ϕ)ϕ)=Dj(x)χϕ(ϕ)dΩ+D12τ|ϕ|2dΩs.t.G(χΩ(ϕ))=Dχϕ(ϕ)dΩVmax0 (4)

    where, D12τ|ϕ|2dΩ is the added virtual energy term, τ is the regularization parameter, which is always greater than zero, χΩ is characteristic function, which can transform the structural topology optimization problem into the material allocation problem in the fixed design domain, so as to determine the optimal distribution of the structure in the design domain.

    The Lagrange method is used to deal with the topology optimization problem with the given volume constraints, and the constrained optimization problem of Eq (4) can be transformed into an unconstrained optimization problem, which is expressed as follows.

    minϕˉJR(χΩ(ϕ)ϕ)=JR+λG=D(j(x)+λ)χϕ(ϕ)dΩλVmax+D12τ|ϕ|2dΩ (5)

    where, λ denotes the Lagrange multiplier, which is always positive and used to deal with the volume constraints, and ˉJR is the Lagrange function. The optimal topology of the structure can be obtained by solving Eq (5).

    However, it is almost impossible to find the optimal solution of Eq (5) directly. To obtain the final topology optimization results, Equation (5) is solved by using the time evolution equation. The time evolution equation is expressed as follows.

    ϕt=K0(dˉJdϕτ2ϕ)inD (6)

    where, ˉJ=D(j(x)+λ)χϕ(ϕ)dΩλVmax. Equation (6) is also called the reaction-diffusion equation (RDE).

    To solve the problem conveniently, on the basis of Eq (6), the normalization coefficient C is introduced to normalize the reaction terms in RDE. Since the reaction term dˉJdϕ is equivalent to the topological derivative dtˉJ numerically, the topological derivative dtˉJ is used for numerical calculation. The level set function is independent of the outside of the fixed design domain, so it is assumed that the boundary conditions of the level set function are Dirichlet boundary conditions on the non-design boundary and Neumann boundary conditions on other boundaries. Therefore, the RDE with boundary conditions can be adjusted to the formulation of Eq (7).

    {ϕt=(CdtˉJτ2ϕ)inDϕ=0onΓ (7)

    where, the value of K0 in Eq (6) is set as 1, referring to the literature (Otomori et al. [28]). The sign of topological derivative dtˉJ is opposite to that of dˉJdϕ, because the sign definitions of topological derivative dtˉJ and level set function ϕ are different. The solving process of topological derivative will be detailed in Section 2.4. The formula for calculating the normalization coefficient C is expressed as follows.

    C=DΩD|dtJ|dΩ (8)

    Through the above derivation, solving Eq (7) can update the level set function and realize the topology optimization of the structure. When the level set method based on RDE is used for structural optimization, the finite element shape function is used to interpolate the level set equation. In order to move the structure boundary to an ideal position, a large number of iterations are needed in the actual operation.

    To obtain better interpolation effect, a continuous and differentiable CS-RBF model is introduced into the structural topology optimization. CS-RBF is used to interpolate the level set equation instead of the finite element shape function to improve the solving process of RDE, which can not only retain the advantages of the level set method based on RDE, but also improve the computational efficiency.

    The second-order continuous CS-RBF proposed by Wendland [13] is selected as the interpolation function to parameterize the level set equation. Taking the 2D structure as an example, the mathematical expression of CS-RBF is given as follows.

    g(r)=max{0,(1r)4+}(4r+1) (9)

    The support radius r is defined as

    r=didm=((xxci)2+(yyci)2)dm (10)

    where, di represents the distance from any point (x,y) within the support radius r to the fixed reference point (xci,yci). The search radius of basis function dm represents the maximum control range of CS-RBF centered on reference point (xci,yci). The selection of dm is very important. If the value of dm is too large, the number of nodes in the interpolation radius will increase and the amount of calculation in the interpolation process will increase. If the value of dm is too small, the interpolation matrix will be singular. In other words, the value of dm can lead to mesh-dependent optimal solutions. To ensure the non-singularity of the function interpolation and simultaneously guarantee the computational efficiency, an experiential criterion in selecting a radius of support is often needed (Luo et al. [57]). In this paper, the value of dm is selected as 2 times of the finite element mesh size to simplify the calculation.

    By using CS-RBF to parameterize the level set equation, N fixed reference points are selected, and the level set equation is decoupled into a series of products of CS-RBF and related expansion coefficients as follows.

    ϕ(x,t)=Ni=1gi(x)αi(t) (11)

    where, gi(x) represents the ith set of CS-RBF related to spatial variables, and αi(t) represent the expansion coefficient corresponding to the ith CS-RBF, which is closely related to virtual time t in the iterative process. Since the CS-RBF matrix is only related to the type of RBF and the position of the reference point, gi(x) only needs to be calculated once in the whole iterative process when the reference point is fixed, and the level set equation is only affected by the expansion coefficient αi(t). Therefore, the expansion coefficient αi(t) is selected as the design variable of PLSM.

    By substituting the level set function interpolated by the CS-RBF into the RDE, the following equations can be obtained

    {ϕt=(CdtˉJτ2ϕ)ϕ(x,t)=g(x)α(t) (12)

    By solving Eq (12), the following equations can be obtained

    g(x)α(ti+1)Δt=CdtˉJ+τ2g(x)α(ti+1)+g(x)α(ti)Δt (13)

    where, Δt is the time step of the equation evolution. According to Eq (13), the updated form of the expansion coefficient αi(t) is shown as follows.

    α(ti+1)=[g(x)Δtτ2g(x)]1[ΔtCdtˉJ+g(x)α(ti)] (14)

    In Eq (3), the upper and lower limit constraints are imposed on the level set function to smooth the points near the structure boundary. Therefore, in the actual calculation, the level set function update formula after CS-RBF parameterization is shown as Eq (15).

    α(ti+1) = [g(x)Δtτ2g(x)]1[ΔtCdtˉJ+g(x)α(ti)]ϕ(ti+1) = g(x)α(ti+1)ϕ(ti+1) = max[1,min[1,ϕ(ti+1)]] (15)

    Fuzzy PID controller (Zadeh [23]) is composed of fuzzy controller and PID controller. The block diagram of the fuzzy PID controller is shown in Figure 3.

    Figure 3.  Block diagram of fuzzy PID controller.

    In Figure 3, r(t) is the target setting value of the control system, and y(t) is the measured value. Firstly, the input variables of ER(t) and EC(t) are transferred to the fuzzy controller for processing, and the output variables of Δkp,Δki,Δkd are obtained. Then, the output variables are transferred to the PID controller to modify the parameters kp0,ki0,kd0, and the modified PID parameters are used to control the system.

    In structural topology optimization, the augmented Lagrange method is usually used to deal with volume constraints. However, at the end of the iteration process, the value of Lagrange multiplier cannot continue to increase. If the volume constraint still does not meet the convergence conditions, the volume fraction cannot be effectively reduced, and the iteration process will oscillate or even fail to converge.

    Therefore, in order to obtain better convergence effect of volume constraint and solve the problem that Lagrange multiplier cannot be effectively increased, fuzzy PID control algorithm is introduced into structural topology optimization. Fuzzy PID controller can calculate the value of Lagrange multiplier according to the difference between the current volume fraction and the target volume fraction, thereby obtaining good convergence effect of volume constraint.

    Fuzzy PID controller uses Mamdani fuzzy controller, and the triangular membership function is used to fuzzify the input variables because of its uniform distribution and high sensitivity in the universe. The Max-Min rule is employed for the fuzzy control decision-making, and the gravity center method is employed for the de-blurring.

    The input variables of the fuzzy controller are: volume deviation ER between the set volume fraction value and the current volume fraction, and volume deviation change rate EC. And the calculation formulas of input variables ER and EC are expressed respectively as follows.

    ERiT=voliTVmaxECiT=ERiTERiT - 1 (16)

    where, ERiT is the volume deviation value in the iTth iteration, and ECiT is the volume deviation change rate in the iTth iteration.

    The vector, consisting of the output variables Δkp,Δki,Δkd of the fuzzy controller, is defined as the correction value of PID controller, and the calculating formula of the output variables Δkp,Δki,Δkd of the fuzzy controller is expressed as follows.

    Δpid=evalfis([ER,EC],α) (17)

    where, Δpid = [ΔkpΔkiΔkd]T, α is the fuzzy rules defined according to Table 1, evalfis is the MATLAB function, which is used to calculate the output variables.

    Table 1.  Fuzzy rules of the change amount of the proportional coefficient Δkp, the change amount of the integral coefficient Δki, and the change amount of the differential coefficient Δkd.
    ERiT ECiT
    NB NM NS Z PS PM PB
    NB PB/NB/PS PB/NB/NS PM/NM/NB PN/NM/NB PS/NM/NB Z/Z/PM Z/Z/PS
    NM PB/NB/PS PB/NB/NS PM/NM/NB PM/NM/NB PS/NS/NM Z/Z/PS NS/PS/Z
    NS PM/NM/Z PM/NM/NS PM/NS/NB PM/NS/NM PS/NS/NM NS/PM/NS NS/PM/Z
    Z PS/NM/Z PM/NM/NS PS/NS/NS PS/NS/NS Z/Z/NS NM/PM/Z NM/PM/Z
    PS PS/NM/Z PS/NS/Z PS/NS/Z Z/Z/Z NS/PS/Z NM/PM/PS NM/PB/Z
    PM PS/Z/PB Z/Z/PM Z/Z/PM NS/PS/NM NM/PS/PM NM/PB/NS NB/PB/PB
    PB Z/Z/PB Z/Z/PM Z/Z/PM NM/PS/PM NM/PM/PM NM/PB/PS NB/PB/PB

     | Show Table
    DownLoad: CSV

    In Table 1, ERiT and ECiT are the input variables of the fuzzy controller. Δkp, Δki, and Δkd are the output variables of the fuzzy controller. NS, NM, NB, Z, PS, PM, and PB indicate the seven values of triangular membership functions used to fuzzify the input and output variables following the fuzzy rules, which are the abbreviated forms of negative small, negative middle, negative big, zero, positive small, positive middle, positive big, respectively.

    Then the output variables are used to modify the PID initial parameters kp0,ki0,kd0. The calculation formula of the modified PID initial parameters is formulated as follows.

    kp=kp0+Δkpki=ki0+Δkikd=kd0+Δkd (18)

    The corrected PID parameters can be used to calculate the Lagrange multiplier and deal with the volume constraint. The calculation formula of the improved Lagrange multiplier is formulated as follows.

    λiT=kpERiT+ki(iT1ERiT)+kdECiT (19)

    where, λiT is the Lagrange multiplier at the iTth iteration. By substituting the Lagrange multiplier updating formula proposed in Eq (19) into Eq (5), the volume constraint can be treated.

    The calculation of topological derivative plays an important role in the process of topology optimization. To simplify the calculating process of topological derivatives and improve the calculating efficiency of topology optimization, the topological derivative is calculated by using the variation of the objective function. The finite difference has the advantages of relatively simple equations and convenient calculation procedures, by which the grid can be easily divided and time can be saved. Moreover, it can make use of the structural topological meshes to easily construct high-precision schemes. Therefore, the finite difference is adopted for sensitivity calculation.

    When only the fixed external load is considered, the objective function, the stiffness matrix and the global displacement vector of structure will change when the internal material is added or deleted. Therefore, according to the variation of the objective function before and after the insertion of hole, the calculating formula of objective function change ΔJ can be obtained as follows.

    ΔJ=12(U+ΔU)T(K+ΔK)(U+ΔU)12UTKU=12(UTΔKU+2ΔUTKU+2ΔUTΔKU+ΔUTKΔU+ΔUTΔKΔU)=12[UTΔKU+ΔUT(K+ΔK)ΔU] (20)

    where, ΔJ is the variation of the objective function; K is the stiffness matrix of structure; U is the global displacement vector; ΔK is the variation of the stiffness matrix; ΔU is the variation of the global displacement vector; ΔUT(K+ΔK)ΔU is the high-order infinitesimal term, which can be ignored. Therefore, the calculating formula of the variation of objective function expressed by Eq (20), namely the calculating formula of the topological derivative, can be expressed as follows.

    ΔJ=dtJ=12UTΔKU (21)

    The mathematical expression of the minimum compliance problem is formulated as follows.

    minφJ=ΓtFvdΓ+DbvχφdΩs.t.a(u,v)=l(v)vU,uUG=DχφdΩVmax  (22)

    where, J represents the overall compliance value of the structure; u is the actual displacement; U is the allowable displacement in the displacement field; v is the virtual displacement; a(u,v)=l(v) is the weak form of the elastic balance equation; a(u,v) is the internal force virtual work of the elastic body, and l(v) is the virtual work of external force, which are expressed respectively as follows.

    a(u,v)=Dε(u):E:ε(v)χϕdΩl(v)=ΓtFvdΓ+ΩbvχϕdΩ (23)

    where, E is the elastic modulus of the material, ε is the strain tensor.

    In this work, we mainly study the optimization problem in the static load equilibrium state where the virtual displacement v and the actual displacement u of the elastic equilibrium equation are equal. The elastic equilibrium equation always satisfies the constraint. Therefore, the mathematical model of Eq (23) can be transformed into a minimum compliance problem with only one volume constraint as follows.

    minφJ=ΓtFvdΓ+DbvχφdΩs.t.G=DχφdΩVmax  (24)

    Here, the Lagrange multiplier method is employed to transform Eq (24) into an unconstrained structural topology optimization problem, as shown in Eq (25).

    ˉJ=ΓtFvdΓ+DbvχϕdΩ+λ(DχϕdΩVmax) (25)

    where, λ denotes the Lagrange multiplier.

    Figure 4 shows the flow chart of the PLSM based on RDE and fuzzy PID control algorithm.

    Figure 4.  Flow chart of numerical algorithm corresponding to the proposed method.

    The implementation process of the optimization algorithm shown in Figure 4 can be also described as the following iterative steps.

    Step 1: Discretize the structure design domain into finite element meshes, define the initial state of level set function ϕ, calculate the CS-RBF matrix g(x) at the reference point, and calculate the initial values of design variables α(t).

    Step 2: Initialize fuzzy PID controller: define the input variables ER and EC, the output variables Δkp,Δki,Δkd of the fuzzy adaptive PID controller according to Eqs (16) and (17). Formulate the fuzzy rules, and set the PID initial parameters kp0,ki0,kd0.

    Step 3: Prescribe the boundary constraints and load conditions of the structure, and calculate the element stiffness matrix.

    Step 4: Start iteration, and set the number of iterations as iT=1:200.

    Step 5: Calculate the structural displacement matrix U, the structural stiffness matrix K, the value of structural compliance J, and the values of topological derivative according to Eq (21).

    Step 6: Update the value of Lagrange multiplier λiT according to Eq (19), and deal with the volume constraint.

    Step 7: Update the design variables α(t) and the level set functions ϕ according to Eq (15).

    Step 8: Set the convergence value as θ=1e3, and judge whether the optimization results satisfy the convergence condition DdΩVmax and JiTJiT1θ. If the convergence condition is not satisfied, repeat steps (5) to (8); if the convergence condition is satisfied, stop the iteration and output the optimization results, the corresponding optimization curves of objective function and volume fraction.

    In this section, to verify that using fuzzy PID control algorithm to calculate Lagrange multiplier can get good effect on processing volume constraint, topology optimization of the Michell beam structure by the PLSM based on the RDE is used for case study. Where, two methods, namely the fuzzy PID control algorithm and the augmented Lagrange method, are used for processing volume constraint, respectively.

    The design domain and boundary conditions of the Michell beam structure are shown in Figure 5. The aspect ratio of the design domain is 5:2, which is divided into 100 times 40 quadrilateral meshes. The lower left boundary is fixed, whereas the lower right boundary is simply supported. The concentrated load of F=100 is applied at the midpoint of the lower side. Elastic modulus of material is set as E0=1. To avoid the singularity of structural stiffness during calculation, elastic modulus of the void is set as Emin=1e9. Poisson's ratio is set as μ=0.3. The volume fraction constraint is set as Vmax=0.5. The regularization parameter is set as τ=2e4. The initial values of fuzzy PID controller are set as kp0=10,ki0=10,kd0=5, respectively. Figure 5 shows the history curves of Lagrange multiplier and volume fraction obtained by the two different methods during topology optimization of the Michell beam structure. Where, Figure 6(a) shows the history curves of Lagrange multiplier and volume fraction obtained by the fuzzy PID control algorithm; and Figure 6(b) shows the history curves of Lagrange multiplier and volume fraction obtained by the augmented Lagrange method.

    Figure 5.  Initial design domain and boundary condition of Michell beam structure.
    Figure 6.  History curves of Lagrange multiplier and volume fraction obtained by two methods for processing volume constraint during topology optimization of the Michell beam structure.

    From Figure 6(a), it can be seen that: at the beginning of the iteration process, the Lagrange multiplier obtained by fuzzy PID control algorithm increases greatly; at the end of the iteration process, the growth rate of Lagrange multiplier gradually decreases; when the value of Lagrange multiplier remains stable, the volume fraction is optimized to its target value and remains stable.

    The iterative format of the augmented Lagrange multiplier with the curves shown in Figure 6(b) is expressed by Eq (26).

    λ = DdtFdΩDdΩexp[p(GGmax+d)] (26)

    According to Eq (26), the value of λ is not zero at the beginning of iteration. During the iterative optimization process, the topology optimization of the structure is realized by continuously deleting the elements with small topological derivative.

    From Figure 6(b), it can be seen that: at the beginning of the iteration process, the value of Lagrange multiplier gradually increases; at the end of the iteration process, the value of Lagrange multiplier oscillates obviously, which leads to the corresponding fluctuation in volume fraction. The fluctuation in volume fraction is not conducive to the convergence of topology optimization. Moreover, the volume fraction is finally optimized to 0.46, which is lower than its target value. According to the curves of Lagrange multiplier and volume fraction in Figure 6, it can be concluded that: compared with the augmented Lagrange method, fuzzy PID control algorithm has obviously good numerical stability and high optimization convergence efficiency. In other words, using fuzzy PID control algorithm to deal with volume constraints can get relatively stable Lagrange multiplier optimization curve, and make the optimization process of volume constraints relatively stable.

    The final topology optimization results of the Michell beam structure obtained by the two different methods for processing volume constraint are shown in Figure 7. Where, Figure 7 (a) shows the final optimization results obtained by the fuzzy PID control algorithm, and Figure 7(b) shows the final optimization results by the augmented Lagrange method. From Figure 7, it can be concluded that the correct and reasonable structure can be obtained by both of the two methods for processing volume constraint.

    Figure 7.  Topology optimization results of the Michell beam structure obtained by two methods for processing volume constraint.

    In fuzzy PID control algorithm, the suitable initial PID control parameters can greatly reduce the calculating burden, accelerate the convergence speed and improve the convergence stability. In this paper, the initial PID control parameters are selected by the trial method. In other words, the initial PID control parameters are determined according to the specific optimization requirements and effect. To illustrate the influence of different values of kp0,ki0,kd0 on the optimization results, topology optimization of the cantilever beam structure by the proposed method is used as a case in this section.

    The design domain and boundary condition of the cantilever beam structure are shown in Figure 8. The aspect ratio of the design domain is 5:3, which is divided into 100 times 60 quadrilateral meshes. The left side of the structure is fully fixed, whereas the other three sides are free. The concentrated load of F=100 is applied at the midpoint of the right side. Elastic modulus of material is set as E0=1. To avoid the singularity of structural stiffness during calculation, elastic modulus of the void is set as Emin=1e9. Poisson's ratio is set as μ=0.3. The regularization parameter is set as τ=2e4. The volume fraction constraint is set as Vmax=0.5.

    Figure 8.  Initial design domain and boundary conditions of the cantilever beam structure.

    To verify the influence of initial PID control parameters on the optimization results, the initial PID control parameters are taken as (kp0 = 5,ki0 = 5,kd0 = 5), (kp0 = 5,ki0 = 5,kd0 = 15), (kp0 = 5,ki0 = 8,kd0 = 15) and (kp0 = 15,ki0 = 8,kd0 = 15), respectively. Table 2 lists the numerical results during topology optimization of the cantilever beam structure with the four groups of different initial PID control parameters for comparison.

    Table 2.  Numerical results during topology optimization of the cantilever beam structure with the four groups of different initial PID control parameters.
    Initial value of PID Iterations Objective function Volume fraction
    kp0 = 5,ki0 = 5,kd0 = 5 250 3.96 times 105 0.50
    kp0 = 5,ki0 = 5,kd0 = 15 250 3.99 times 105 0.50
    kp0 = 5,ki0 = 8,kd0 = 15 250 3.99 times 105 0.50
    kp0 = 15,ki0 = 8,kd0 = 15 250 4.00 times 105 0.50

     | Show Table
    DownLoad: CSV

    Figures 9, 10 and 11 show the history curves of Lagrange multiplier, objective function and volume fraction during topology optimization of the cantilever beam structure with the four groups of different initial PID control parameters, respectively. Where, the red curves a correspond to the initial PID control parameters kp0 = 5,ki0 = 5,kd0 = 5, the blue curves b correspond to the initial PID control parameters kp0 = 5,ki0 = 5,kd0 = 15, the green curves c correspond to the initial PID control parameters kp0 = 5,ki0 = 8,kd0 = 15 and the black curves d correspond to the initial PID control parameters kp0 = 15,ki0 = 8,kd0 = 15.

    Figure 9.  History curves of Lagrange multiplier during topology optimization of the cantilever beam structure with four groups of different initial PID control parameters.
    Figure 10.  History curves of objective function during topology optimization of the cantilever beam structure with four groups of different initial PID control parameters.
    Figure 11.  History curves of volume fraction during topology optimization of the cantilever beam structure with four groups of different initial PID control parameters.

    By comparing curve a and curve b in Figure 9, it can be found that the optimization history curve of Lagrange multiplier does not change obviously when only increasing the value of initial parameter kd0. By comparing curve b and curve c in Figure 9, it can be found that the growth rate of Lagrange multiplier is obviously accelerated when only increasing the value of initial parameter ki0. By comparing curve c and curve d in Figure 9, it can be found that the growth rate of Lagrange multiplier will also increase by a relatively small margin when only increasing the value of initial parameter kp0. Accordingly, it can be concluded that the value of Lagrange multiplier can be adjusted and controlled by changing the values of initial PID control parameters. Specifically, increasing the value of initial parameters kp0 and ki0 can accelerate the growth rate of Lagrange multiplier, and the change of initial parameter ki0 can control the growth rate of Lagrange multiplier more obviously.

    From Figures 10 and 11, it can be found that: when the initial PID parameters are taken as kp0 = 5,ki0 = 5,kd0 = 5, about 150 iterations are needed to make the optimization process converge, the final objective function value is 3.96 times 105, and the optimization history curves of objective function and volume fraction have small fluctuations, resulting in long adjustment time and poor stability of optimization process.

    By comparing curve a and curve b in Figures 10 and 11, it can be found that: when keeping the values of initial parameters kp0 and ki0 unchanged and only increasing the value of initial parameter kd0, it also takes about 150 iterations to make the optimization process converge, and the final objective function value is 3.99 times 105. It is certified that increasing the value of initial parameter kd0 can make the optimization history curves of objective function and volume fraction more smooth and improve the stability of optimization process.

    By comparing curve b and curve c in Figures 10 and 11, it can be found that when keeping the values of initial parameters kp0 and kd0 unchanged and only increasing the value of initial parameter ki0, the increasing rate of objective function and the decreasing rate of volume fraction are both obviously accelerated, it takes about 80 iterations to make the optimization process converge, the final objective function value is 3.99 times 105, and the optimization history curves of objective function and volume fraction are relatively smooth, without wild fluctuations and oscillations. It is certified that increasing the value of initial parameter ki0 can reduce the number of iterations and improve the optimization efficiency.

    By comparing curve c and curve d in Figures 10 and 11, it can be found that when keeping the values of initial parameters ki0 and kd0 unchanged and only increasing the value of initial parameter kp0, the increasing rate of objective function and the decreasing rate of volume fraction are both further accelerated, it takes about 50 iterations to make the optimization process converge, the final objective function value is 4.00 times 105, and the optimization history curves of objective function and volume fraction are smooth, implying that the optimization process is stable. It is proved that increasing the value of initial parameter kp0 can also improve the optimization efficiency and obtain satisfactory results. The maximum difference between the final objective function values obtained by the four groups of initial PID control parameters is less than 2%, which proves that changing the values of the initial PID control parameters will not affect the objective function significantly.

    Therefore, according to the specific effects of the initial PID control parameters on the optimization results, the values of the initial PID control parameters can be selected appropriately, so as to obtain stable and fast topology optimization process and accurate topology optimization results.

    In this section, several typical 2D and 3D structures are optimized as numerical examples to illustrate the effectiveness of the proposed method. Taking the minimum compliance as the objective function and the volume fraction as the constraint condition, the optimization results obtained by the proposed method are analyzed and compared with those obtained by other methods to demonstrate the stability and effectiveness of this method.

    Figure 12 shows the design domain and boundary conditions of the L-beam structure, which is divided into 3200 quadrilateral meshes. The upper boundary is fixed, whereas a vertical downward force F=100 is applied at the midpoint of the right boundary. The design domain is divided into 3200 elements. Elastic modulus of material is set as E0=1. To avoid the singularity of structural stiffness during calculation, elastic modulus of the void is set as Emin=1e9. Poisson's ratio is set as μ=0.3. The volume fraction constraint is set as Vmax=0.5. The regularization parameter is set as τ=2e4. The initial values of fuzzy PID controller are set as kp0=15,ki0=12,kd0=15, respectively.

    Figure 12.  Initial design domain and boundary conditions of the L-beam structure.

    Figure 13 shows the topology optimization result of the L-beam structure obtained by the proposed method. Figure 14 shows the history curves of objective function and volume fraction during topology optimization of the L-beam structure by the proposed method.

    Figure 13.  Topology optimization result of the L-beam structure obtained by the proposed method.
    Figure 14.  History curves of objective function and volume fraction during topology optimization of the L-beam structure by the proposed method.

    From Figure 13, the reasonable topology configuration of the L-beam structure is obtained by the proposed method, and the structure boundary is smooth and clear. From Figure 14, the L-beam structure from initial volume fraction of 0.75 to target volume fraction of 0.4, it takes about 50 iterations, and the final value of objective function is 6.40 times 105. In the optimization process, the optimization curves of objective function and volume fraction are relatively smooth, almost no fluctuations, and the optimization process is stable.

    Figure 15 shows the design domain and boundary conditions of the multi-load bridge structure. The aspect ratio of the design domain is 8:3, which is divided into 80 times 30 quadrilateral meshes. The lower left boundary is fixed, whereas the lower right boundary is simply supported. A concentrated load of F=60 is applied at the midpoint of the bottom of bridge, and two concentrated loads of F=30 are applied at the one-fourth point and the three-fourth point of the bottom of bridge, respectively. Elastic modulus of material is set as E0=1. To avoid the singularity of structural stiffness during calculation, elastic modulus of the void is set as Emin=1e9. Poisson's ratio is set as μ=0.3. The volume fraction constraint is set as Vmax=0.5. The regularization parameter is set as τ=2e4. The initial values of fuzzy PID controller are set as kp0=15,ki0=12,kd0=15, respectively.

    Figure 15.  Design domain model and boundary conditions of the multi-load bridge structure.

    To further verify the effectiveness of the proposed method, the proposed method, level set method based on RDE (Otomori et al. [28]) and PLSM based on CS-RBF (Cui et al. [33]) are respectively utilized to address topology optimization of the multi-load bridge structure, and the optimization results are analyzed and compared in the meanwhile. Figure 16 shows the topology optimization results of the multi-load bridge structure obtained by the above-mentioned three methods. Figure 17 shows the history curves of objective function during topology optimization of the multi-load bridge structure by the above-mentioned three methods. Figure 18 shows the history curves of volume fraction during topology optimization of the multi-load bridge structure by the above-mentioned three methods. In Figures 17 and 18, the red curve b represents the curve obtained by the proposed method, the blue curve c represents the curve obtained by the level set method based on RDE, and the green curve d represents the curve obtained by the PLSM based on CS-RBF. Table 3 lists the numerical results in topology optimization of the multi-load bridge by the three methods for comparison.

    Figure 16.  Topology optimization results of the multi-load bridge structure by the three methods.
    Figure 17.  History curves of objective function during topology optimization of the multi-load bridge structure by the three methods.
    Figure 18.  History curves of volume fraction during topology optimization of the multi-load bridge structure by the three methods.
    Table 3.  Numerical results for topology optimization of the multi-load bridge by the three methods.
    Optimization method Final value of compliance Final value of volume fraction
    The proposed method 1.81 times 105 0.50
    LSM based on RDE 2.02 times 105 0.50
    PLSM based on CS-RBF 1.90 times 105 0.50

     | Show Table
    DownLoad: CSV

    From the image optimization effect in Figure 16, the topological structures obtained by the three optimization methods are basically the same, and the correct and reasonable optimal structures can be obtained. However, compared with the level set method based on RDE, the structure boundary obtained by the proposed method is smoother and clearer, and compared with the PLSM based on CS-RBF has simpler topological configuration. According to the optimization curves in Figure 17, Figure 18 and the optimization results in Table 3, the proposed method needs about 50 iterations when the volume constraint converges, and the objective function value is 1.81 times 105. Compared with PLSM based on CS-RBF, the objective function value is reduced by 4.7%. Compared with the level set method based on RDE, the objective function value is reduced by 11.1%, and better structural performance is obtained. It can be found that the introduction of virtual interface energy fully relaxes the topology optimization problem, obtains the optimal structure with smooth boundary, and improves the numerical stability of the optimization problem. Therefore, by comparing the three optimization methods, it can be concluded that the proposed method can achieve satisfactory optimization results in terms of optimization effect and computational efficiency.

    Figure 19 shows the design domain and boundary conditions of the 3D cantilever beam structure. The ratio of the three dimensions of the design domain is 15:10:1, which is divided into 60 times 40 times 4 elements. The left boundary is fixed, whereas the right midpoint is applied with vertical downward load F=100. Elastic modulus of material is set as E0=1. To avoid the singularity of structural stiffness during calculation, elastic modulus of the void is set as Emin=1e9. Poisson's ratio is set as μ=0.3. The target volume fraction is set as Vmax=0.5. The regularization parameter is set as τ=2e4. The initial PID parameter are set as kp0=40,ki0=20,kd0=15, respectively.

    Figure 19.  Initial design domain and boundary conditions of the 3D cantilever beam structure.

    The topology optimization result of the 3D cantilever structure obtained by the proposed method is shown in Figure 20. The history curves of the objective function and volume fraction during topology optimization of the 3D cantilever structure by the proposed method are shown in Figure 21.

    Figure 20.  Topology optimization results of the 3D cantilever beam structure obtained by the proposed method.
    Figure 21.  History curves of objective function and volume fraction during topology optimization of the 3D cantilever beam structure by the proposed method.

    It can be seen from Figure 20, using the algorithm proposed in this paper to optimize the 3D cantilever structure, clear structure boundary can be obtained, and there is no bad numerical phenomenon such as checkerboard. From the optimization curves shown in Figure 21, the objective function and the volume fraction change smoothly in the whole optimization process, and there is almost no fluctuation and oscillation phenomenon. And the final objective function value is 9.88 times 105. Therefore, it can be concluded that although the number of 3D structure meshing is several times larger than 2D structure, the proposed algorithm can improve the optimization efficiency, ensure the stability of the optimization process, and get a clear structure configuration.

    Figure 22 shows the design domain and boundary conditions of the 3D arch bridge structure. The ratio of the three dimensions of the design domain is 15:5:1, which is divided into 60 times 20 times 4 elements. The displacement of nodes within two-unit length at the left and the right ends of the lower boundary of the structure is fully constrained, whereas the plate area with one-unit thickness on the upper surface of the structure is non-design area. The vertical downward load F=5 is applied to all nodes on the upper surface of the structure. The design domain Elastic modulus of material is set as E0=1. To avoid the singularity of structural stiffness during calculation, elastic modulus of the void is set as Emin=1e9. Poisson's ratio is set as μ=0.3. The target volume fraction is set as Vmax=0.4. The regularization parameter is set as τ=2e4. The initial PID parameter are set as kp0=80,ki0=45,kd0=40, respectively.

    Figure 22.  Initial design domain and boundary conditions of the 3D arch bridge structure.

    Figure 23 shows the topology optimization results of the 3D arch bridge structure from different angles, in which (a) is the front view and (b) is the side view. Figure 24 shows the history curves of objective function and volume fraction during topology optimization of the 3D arch bridge structure by the proposed method.

    Figure 23.  Topology optimization result of the 3D arch bridge structure obtained by the proposed method from different angles.
    Figure 24.  History curves of objective function and volume fraction during topology optimization of the 3D arch bridge structure by the proposed method.

    From Figures 23 and 24, it can be seen that the final result obtained by using the algorithm proposed in this paper to optimize the arch bridge structure is similar to the arch bridge structure in practical application, which proves that the method proposed in this paper has certain reference value in practical engineering application. Due to the existence of non- design area in the optimization process, it takes about 70 iterations to reduce the volume fraction from the initial value of 0.95 to the target value of 0.4, and the final objective function value is 1.10 times 106. The results show that using the proposed algorithm for topology optimization of 3D structure can get the correct structure configuration, and the optimization efficiency is high, the convergence speed is high, and the iterative process is stable, which has a certain practicality in engineering applications.

    In this paper, a PLSM based on the RDE and the fuzzy PID control algorithm is proposed to solve structural topology optimization problem of minimum compliance under volume constraints. Firstly, the RDE is used as the evolution equation of the level set function, meanwhile, the topological derivative of the material domain is used as the reaction term of the RDE to drive the update of level set function, which can generate holes in the material domain and ensure the stability of numerical calculation. Secondly, in the process of solving the RDE, the CS-RBF is used to interpolate the level set function to correct the diffusion term in RDE, which can improve the efficiency and accuracy of calculation because the CS-RBF owns the advantages of sparse coefficient matrix and local correlation of parameters. Thirdly, the fuzzy PID control algorithm is used to adjust the value of Lagrange multiplier so as to control the volume constraint of structural topology optimization. The comparative examples illustrate that a more stable optimization history curve of Lagrange multiplier is obtained by the fuzzy PID control algorithm, which improves the stability of the convergence process of volume constraint as well as topology optimization. Thus, in addition to inheriting the merits of the PLSM, the proposed method takes the advantages of the RDE, the CS-RBF and the fuzzy PID control algorithm. By this method, not only can clear and smooth boundary of structural topology optimization result be obtained, but also the stability and efficiency of structural topology optimization process can be improved. Furthermore, this method has been applied to 3D structural topology optimization, which broadens its application scope and confirms its effectiveness in the meanwhile.

    The authors would like to express thanks to the reviewers for their valuable suggestions for improving the manuscript. This work was supported by Natural Science Foundation of Shaanxi Province (2020JM-200) and China Scholarship Council (201506965015).

    The authors declare there is no conflict of interest.



    [1] B. Bourdin, A. Chambolle, Design-dependent loads in topology optimization, ESAIM–Contr. Optim. Ca., 9 (2003), 19–48. https://doi.org/10.1051/cocv:2002070 doi: 10.1051/cocv:2002070
    [2] D. Muñoz, J. J. Ródenas, E. Nadal, J. Albelda, 3D topology optimization with h-adaptive refinement using cartesian grids finite element method (cgFEM), In: Proceedings of the 6th International Conference on Engineering Optimization, Springer, Cham, 2018. https://doi.org/10.1007/978-3-319-97773-7_68
    [3] D. P. Peng, B. Merriman, S. Osher, H. K. Zhao, M. J. Kang, A PDE-based fast local level set method, J. Comput. Phys., 155 (1999), 410–438. https://doi.org/10.1006/jcph.1999.6345 doi: 10.1006/jcph.1999.6345
    [4] F. Ferrari, O. Sigmund, A new generation 99 line Matlab code for compliance topology optimization and its extension to 3D, Struct. Multidiscip. O., 62 (2020), 2211–2228. https://doi.org/10.1007/s00158-020-02629-w doi: 10.1007/s00158-020-02629-w
    [5] G. Allaire, F. de Gournay, F. Jouve, A. M. Toader, Structural optimization using topological and shape sensitivity via a level set method, Control Cybern., 34 (2005), 59–80.
    [6] G. Allaire, F. Jouve, A. M. Toader, Structural optimization using sensitivity analysis and a level-set method, J. Comput. Phys., 194 (2004), 363–393. https://doi.org/10.1016/j.jcp.2003.09.032 doi: 10.1016/j.jcp.2003.09.032
    [7] G. Allaire, F. Jouve, A. M. Toader, A level-set method for shape optimization, Comptes Rendus Math., 334 (2002), 1125–1130. https://doi.org/10.1016/S1631-073X(02)02412-3 doi: 10.1016/S1631-073X(02)02412-3
    [8] H. A. Eschenauer, V. V. Kobelev, A. Schumacher, Bubble method for topology and shape optimization of structures, Struct. Optimization, 8 (1994), 42–51. https://doi.org/10.1007/BF01742933 doi: 10.1007/BF01742933
    [9] H. Li, T. Kondoh, P. Jolivet, K. Furuta, T. Yamada, B. Zhu, et al., Optimum design and thermal modeling for 2D and 3D natural convection problems incorporating level set-based topology optimization with body-fitted mesh, Int. J. Numer. Meth. Eng., 123 (2022), 1954–1990. https://doi.org/10.1002/nme.6923 doi: 10.1002/nme.6923
    [10] H. Li, T. Kondoh, P. Jolivet, K. Furuta, T. Yamada, B. Zhu, et al., Three-dimensional topology optimization of a fluid-structure system using body-fitted mesh adaption based on the level-set method, Appl. Math. Model., 101 (2022), 276–308. https://doi.org/10.1016/j.apm.2021.08.021 doi: 10.1016/j.apm.2021.08.021
    [11] H. Li, T. Yamada, P. Jolivet, K. Furuta, T. Kondoh, K. Izui, et al., Full-scale 3D structural topology optimization using adaptive mesh refinement based on the level-set method, Finite Elem. Anal. Des., 194 (2021), 103561. https://doi.org/10.1016/j.finel.2021.103561 doi: 10.1016/j.finel.2021.103561
    [12] H. S. Ho, M. Y. Wang, M. D. Zhou, Parametric structural optimization with dynamic knot RBFs and partition of unity method, Struct. Multidiscip. O., 47 (2013), 353–365. https://doi.org/10.1007/s00158-012-0848-7 doi: 10.1007/s00158-012-0848-7
    [13] H. Wendland, Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree, Adv. Comput. Math., 4 (1995), 389–396. https://doi.org/10.1007/BF02123482 doi: 10.1007/BF02123482
    [14] H. Zhang, S. T. Liu, X. O. Zhang, Topology optimization of 3D structures with design-dependent loads, Acta Mech. Sin., 26 (2010), 767–775. https://doi.org/10.1007/s10409-010-0370-3 doi: 10.1007/s10409-010-0370-3
    [15] J. A. Sethian, A. Wiegmann, Structural boundary design via level set and immersed interface methods, J. Comput. Phys., 163 (2000), 489–528. https://doi.org/10.1006/jcph.2000.6581 doi: 10.1006/jcph.2000.6581
    [16] J. A. Sethian, P. Smereka, Level set methods for fluid interfaces, Annu. Rev. Fluid Mech., 35 (2003), 341–372. https://doi.org/10.1146/annurev.fluid.35.101101.161105 doi: 10.1146/annurev.fluid.35.101101.161105
    [17] J. Du, N. Olhoff, Topological optimization of continuum structures with design-dependent surface loading – Part II: algorithm and examples for 3D problems, Struct. Multidiscip. O., 27 (2004), 166–177. https://doi.org/10.1007/s00158-004-0380-5 doi: 10.1007/s00158-004-0380-5
    [18] J. S. Choi, T. Yamada, K. Izui, S. Nishiwaki, J. Yoo, Topology optimization using a reaction–diffusion equation, Comput. Method. Appl. M., 200 (2011), 2407–2420. https://doi.org/10.1016/j.cma.2011.04.013 doi: 10.1016/j.cma.2011.04.013
    [19] J. Sokolowski, A. Zochowski, Topological derivative in shape optimization, Springer, Boston, MA (2009).
    [20] J. Zhu, Y. Zhao, W. Zhang, X. Gu, T. Gao, J. Kong, et al., Bio-inspired feature-driven topology optimization for rudder structure design, Engineered Sci., 5 (2019), 46–55. https://doi.org/10.30919/es8d716 doi: 10.30919/es8d716
    [21] K. Liu, A. Tovar, An efficient 3D topology optimization code written in Matlab, Struct. Multidiscip. O., 50 (2014), 1175–1196. https://doi.org/10.1007/s00158-014-1107-x doi: 10.1007/s00158-014-1107-x
    [22] K. Svanberg, The method of moving asymptotes–a new method for structural optimization, INT J. Numer. Meth. Eng., 24 (1987), 359–373. https://doi.org/10.1002/nme.1620240207 doi: 10.1002/nme.1620240207
    [23] L. A. Zadeh, Fuzzy sets, Inf. Control., 8 (1965), 338–353. https://doi.org/10.1016/S0019-9958(65)90241-X doi: 10.1016/S0019-9958(65)90241-X
    [24] M. Burger, B. Hackl, W. Ring, Incorporating topological derivatives into level set methods, J. Comput. Phys., 194 (2004), 344–362. https://doi.org/10.1016/j.jcp.2003.09.033 doi: 10.1016/j.jcp.2003.09.033
    [25] M. H. Abolbashari, S. Keshavarzmanesh, On various aspects of application of the evolutionary structural optimization method for 2D and 3D continuum structures, Finite Elem. Anal. Des., 42 (2006), 478–491. https://doi.org/10.1016/j.finel.2005.09.004 doi: 10.1016/j.finel.2005.09.004
    [26] M. J. de Ruiter, F. van Keulen, Topology optimization using a topology description function, Struct. Multidiscip. O., 26 (2004), 406–416. https://doi.org/10.1007/s00158-003-0375-7 doi: 10.1007/s00158-003-0375-7
    [27] M. Marino, F. Auricchio, A. Reali, E. Rocca, U. Stefanelli, Mixed variational formulations for structural topology optimization based on the phase-field approach, Struct. Multidiscip. O., 64 (2021), 2627–2652. https://doi.org/10.1007/s00158-021-03017-8 doi: 10.1007/s00158-021-03017-8
    [28] M. Otomori, T. Yamada, K. Izui, S. Nishiwaki, Matlab code for a level-set based topology optimization method using a reaction diffusion equation, Struct. Multidiscip. O., 51 (2015), 1159–1172. https://doi.org/10.1007/s00158-014-1190-z doi: 10.1007/s00158-014-1190-z
    [29] M. P. Bendsøe, N. Kikuchi, Generating optimal topologies in structural design using a homogenization method, Comput. Method. Appl. M., 71 (1988), 197–224. https://doi.org/10.1016/0045-7825(88)90086-2 doi: 10.1016/0045-7825(88)90086-2
    [30] M. P. Bendsøe, O. Sigmund, Material interpolation schemes in topology optimization, Arch. Appl. Mech., 69 (1999), 635–654. https://doi.org/10.1007/s004190050248 doi: 10.1007/s004190050248
    [31] M. P. Bendsøe, Optimal shape design as a material distribution problem, Struct. Optimization, 1 (1989), 193–202. https://doi.org/10.1007/BF01650949 doi: 10.1007/BF01650949
    [32] M. Sussman, P. Smereka, S. Osher, A level set approach for computing solutions to incompressible two-phase flow, J. Comput. Phys., 114 (1994), 146–159. https://doi.org/10.1006/jcph.1994.1155 doi: 10.1006/jcph.1994.1155
    [33] M. T. Cui, C. C. Luo, G. Li, M. Pan, The parameterized level set method for structural topology optimization with shape sensitivity constraint factor, Eng. Comput., 37 (2021), 855–872. https://doi.org/10.1007/s00366-019-00860-8 doi: 10.1007/s00366-019-00860-8
    [34] M. Y. Wang, H. M. Zong, Q. P. Ma, Y. Tian, M. D. Zhou, Cellular level set in B-splines (CLIBS): A method for modeling and topology optimization of cellular structures, Comput. Method. Appl. M., 349 (2019), 378–404. https://doi.org/10.1016/j.cma.2019.02.026 doi: 10.1016/j.cma.2019.02.026
    [35] M. Y. Wang, X. Wang, D. Guo, A level set method for structural topology optimization, Comput. Method. Appl. M., 192 (2003), 227–246. https://doi.org/10.1016/S0045-7825(02)00559-5 doi: 10.1016/S0045-7825(02)00559-5
    [36] M. Zhou, G. I. N. Rozvany, The COC algorithm, Part II: Topological, geometrical and generalized shape optimization, Comput. Method. Appl. M., 89 (1991), 309–336. https://doi.org/10.1016/0045-7825(91)90046-9 doi: 10.1016/0045-7825(91)90046-9
    [37] M. Zhou, M. Y. Wang, A semi-Lagrangian level set method for structural optimization, Struct. Multidiscip. O., 46 (2012), 487–501. https://doi.org/10.1007/s00158-012-0842-0 doi: 10.1007/s00158-012-0842-0
    [38] N. P. van Dijk, K. Maute, M. Langelaar, F. van Keulen, Level-set methods for structural topology optimization: a review, Struct. Multidiscip. O., 48 (2013), 437–472. https://doi.org/10.1007/s00158-013-0912-y doi: 10.1007/s00158-013-0912-y
    [39] O. Sigmund, P. M. Clausen, Topology optimization using a mixed formulation: An alternative way to solve pressure load problems, Comput. Method. Appl. M., 196 (2007), 1874–1889. https://doi.org/10.1016/j.cma.2006.09.021 doi: 10.1016/j.cma.2006.09.021
    [40] P. Wei, M. Y. Wang, Piecewise constant level set method for structural topology optimization, Int. J. Numer. Meth. Eng., 78 (2009), 379–402. https://doi.org/10.1002/nme.2478 doi: 10.1002/nme.2478
    [41] P. Wei, Z. Y. Li, X. P. Li, M. Y. Wang, An 88-line MATLAB code for the parameterized level set method based topology optimization using radial basis functions, Struct. Multidiscip. O., 58 (2018), 831–849. https://doi.org/10.1007/s00158-018-1904-8 doi: 10.1007/s00158-018-1904-8
    [42] Q. Xia, M. Y. Wang, S. Y. Wang, S. K. Chen, Semi-Lagrange method for level-set-based structural topology and shape optimization, Struct. Multidiscip. O., 31 (2006), 419–429. https://doi.org/10.1007/s00158-005-0597-y doi: 10.1007/s00158-005-0597-y
    [43] R. Malladi, J. A. Sethian, B. C. Vemuri, Shape modeling with front propagation: a level set approach, IEEE T. Pattern. Anal., 17 (1995), 158–175. https://doi.org/10.1109/34.368173 doi: 10.1109/34.368173
    [44] S. Osher, F. Santosa, Level set methods for optimization problems involving geometry and constrains I. Frequencies of a two-density inhomogeneous drum, J. Comput. Phys., 171 (2001), 272–288. https://doi.org/10.1006/jcph.2001.6789 doi: 10.1006/jcph.2001.6789
    [45] S. Osher, J. A. Sethian, Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton–Jacobi formulations, J. Comput. Phys., 79 (1988), 12–49. https://doi.org/10.1016/0021-9991(88)90002-2 doi: 10.1016/0021-9991(88)90002-2
    [46] S. Osher, N. Paragios, Geometric level set methods in imaging, vision, and graphics, Springer, New York (2003). https://doi.org/10.1007/b97541
    [47] S. Y. Wang, M. Y. Wang, Radial basis functions and level set method for structural topology optimization, Int. J. Numer. Meth. Eng., 65 (2006), 2060–2090. https://doi.org/10.1002/nme.1536 doi: 10.1002/nme.1536
    [48] T. Cecil, J. L. Qian, S. Osher, Numerical methods for high dimensional Hamilton-Jacobi equations using radial basis functions, J. Comput. Phys., 196 (2004), 327–347. https://doi.org/10.1016/j.jcp.2003.11.010 doi: 10.1016/j.jcp.2003.11.010
    [49] T. Yamada, K. Izui, S. Nishiwaki, A. Takezawa, A topology optimization method based on the level set method incorporating a fictitious interface energy, Comput. Method. Appl. M., 199 (2010), 2876–2891. https://doi.org/10.1016/j.cma.2010.05.013 doi: 10.1016/j.cma.2010.05.013
    [50] T. Zegard, G. H. Paulino, GRAND3–Ground structure based topology optimization for arbitrary 3D domains using MATLAB, Struct. Multidiscip. O., 52 (2015), 1161–1184. https://doi.org/10.1007/s00158-015-1284-2 doi: 10.1007/s00158-015-1284-2
    [51] W. Zhang, J. Chen, X. Zhu, J. Zhou, D. Xue, X. Lei, et al., Explicit three dimensional topology optimization via Moving Morphable Void (MMV) approach, Comput. Method. Appl. M., 322 (2017), 590–614. https://doi.org/10.1016/j.cma.2017.05.002 doi: 10.1016/j.cma.2017.05.002
    [52] W. Zhang, Y. Zhou, J. Zhu, A comprehensive study of feature definitions with solids and voids for topology optimization, Comput. Method. Appl. M., 325 (2017), 289–313. https://doi.org/10.1016/j.cma.2017.07.004 doi: 10.1016/j.cma.2017.07.004
    [53] X. Guo, W. Zhang, W. Zhong, Doing topology optimization explicitly and geometrically–a new Moving Morphable Components Based Frame, J. Appl. Mech., 81 (2014), 081009. https://doi.org/10.1115/1.4027609 doi: 10.1115/1.4027609
    [54] X. Y. Yang, Y. M. Xie, G. P. Steven, Evolutionary methods for topology optimisation of continuous structures with design dependent loads, Comput. Struct., 83 (2005), 956–963. https://doi.org/10.1016/j.compstruc.2004.10.011 doi: 10.1016/j.compstruc.2004.10.011
    [55] Y. M. Xie, G. P. Steven, Evolutionary Structural Optimization, Springer, London (1997)
    [56] Y. Zhou, W. Zhang, J. Zhu, Z. Xu, Feature-driven topology optimization method with signed distance function, Comput. Method. Appl. M., 310 (2016), 1–32. https://doi.org/10.1016/j.cma.2016.06.027 doi: 10.1016/j.cma.2016.06.027
    [57] Z. Luo, M. Y. Wang, S. Y. Wang, P. Wei, A level set-based parameterization method for structural shape and topology optimization, Int. J. Numer. Meth. Eng., 76 (2008), 1–26. https://doi.org/10.1002/nme.2092 doi: 10.1002/nme.2092
    [58] Z. Luo, N. Zhang, W. Gao, H. Ma, Structural shape and topology optimization using a meshless Galerkin level set method, Int. J. Numer. Meth. Eng., 90 (2012), 369–389. https://doi.org/10.1002/nme.3325 doi: 10.1002/nme.3325
  • This article has been cited by:

    1. Yuhai Zhong, Huashan Feng, Hongbo Wang, Runxiao Wang, Weiwei Yu, A bionic topology optimization method with an additional displacement constraint, 2023, 31, 2688-1594, 754, 10.3934/era.2023037
    2. Giovanna C. Andrade, Sandra A. Santos, A level-set-based topology optimization strategy using radial basis functions and a Hilbertian velocity extension, 2022, 111, 0307904X, 108, 10.1016/j.apm.2022.06.026
    3. Fengling Jin, Ye Yang, Zhibin Xiao, Boqing Gao, Haopeng Lou, Shear wall layout optimization of multi-tower buildings based on conceptual design and extended evolutionary structural optimization method, 2023, 0305-215X, 1, 10.1080/0305215X.2022.2161535
    4. Davide Cortellessa, Nicola Ferro, Simona Perotto, Stefano Micheletti, Enhancing level set-based topology optimization with anisotropic graded meshes, 2023, 447, 00963003, 127903, 10.1016/j.amc.2023.127903
    5. An Xu, Shumin Li, Jiyang Fu, Anil Misra, Ruohong Zhao, A hybrid method for optimization of frame structures with good constructability, 2023, 276, 01410296, 115338, 10.1016/j.engstruct.2022.115338
    6. Anatoliy Martynyuk, Gani Stamov, Ivanka Stamova, Yulya Martynyuk–Chernienko, Regularization scheme for uncertain fuzzy differential equations: Analysis of solutions, 2023, 31, 2688-1594, 3832, 10.3934/era.2023195
    7. Feiteng Cheng, Qinghai Zhao, Zibin Mao, Fajie Wang, Mechanical Response of Gradient Lattice Structures Based on Topology Optimization, 2024, 26, 1438-1656, 10.1002/adem.202301887
    8. Qianfeng Cao, Ning Dai, Jiuyang Bai, Hongqing Dai, A dimensionality reduction optimization method for thin-walled parts to realize mass center control and lightweight design, 2024, 0305-215X, 1, 10.1080/0305215X.2024.2350648
    9. Martin Sotola, Pavel Marsalek, David Rybansky, Jan Kopacka, Dusan Gabriel, Ondrej Jezek, Ludek Kovar, Josef Tejc, Miloslav Pasek, Radim Halama, Michal Barnovsky, Application of Surface-Based Smoothing Methods for Topology Optimization Results, 2024, 16, 1758-8251, 10.1142/S1758825124500868
    10. Mingtao Cui, Wang Li, Guang Li, Xiaobo Wang, The asymptotic concentration approach combined with isogeometric analysis for topology optimization of two-dimensional linear elasticity structures, 2023, 31, 2688-1594, 3848, 10.3934/era.2023196
    11. Qianqian Sui, Jun Yan, Zhirui Fan, Zhihui Liu, Stress-constrained topology optimization for multi-material thermo-hyperelastic compliant mechanisms based on inverse motion, 2024, 0305-215X, 1, 10.1080/0305215X.2023.2301552
    12. Mingtao Cui, Wennan Cui, Wang Li, Xiaobo Wang, A polygonal topology optimization method based on the alternating active-phase algorithm, 2024, 32, 2688-1594, 1191, 10.3934/era.2024057
    13. Mohammed Hameed Hafeeth, Ying Liu, An efficient volume-preserving binary filter for additive manufacturing in topology optimization, 2024, 0305-215X, 1, 10.1080/0305215X.2024.2382799
    14. Xingjun Gao, Longhua Li, Yingxiong Li, Meiling Dai, Robust topology optimization of multi-material structures based on the material-field series-expansion method, 2024, 0305-215X, 1, 10.1080/0305215X.2024.2340062
    15. Xingjun Gao, Longhua Li, Jingtong Chen, Yingxiong Li, Robust topology optimization of multi-material structures with overhang angle constraints using the material field series-expansion method, 2024, 69, 23520124, 107359, 10.1016/j.istruc.2024.107359
    16. J. Li, B. Kang, Y. Sun, Y. Huang, J. Guo, Z. Gao, Dynamic Characteristic Analysis of the Shield Machine Three-Stage Planetary Reducer Considering Tooth Wear, 2024, 12, 2523-3920, 7291, 10.1007/s42417-023-01188-w
    17. Yun Ni, Jinqing Zhan, Min Liu, Topological design of continuum structures with global stress constraints considering self-weight loads, 2023, 31, 2688-1594, 4708, 10.3934/era.2023241
    18. Wangyu Liu, Guanghui Huang, Weigui Xie, An efficient cross-platform multi-material topology optimization approach occupying enhanced BESO method, 2024, 0025-6455, 10.1007/s11012-024-01916-w
    19. Hongshuo Fan, Jianli Liu, Haobo Zhang, Tao Nie, Jingui Yu, Jianzhong Yang, Zhaohui Xia, Evolutionary topology optimization for elastoplastic structures via isogeometric analysis, 2025, 0305-215X, 1, 10.1080/0305215X.2024.2443738
    20. Qiwei Hu, Shengbo Hu, Mengxia Liu, A Multi-Objective Nutcracker Optimization Algorithm Based on Cubic Chaotic Map for Numerical Association Rule Mining, 2025, 15, 2076-3417, 1611, 10.3390/app15031611
    21. Rongbo Zhou, Shuli Sun, A degree-based topological optimization method for triangular meshes, 2025, 0020-7160, 1, 10.1080/00207160.2025.2466765
    22. Fangyi Li, Yuming Yin, Rong Zeng, Dachang Zhu, Topology optimization and mechanical performance analysis of P-type triply periodic minimal surface structures, 2025, 47, 1678-5878, 10.1007/s40430-025-05439-7
    23. Tianbo Wang, Zhaoyi Luo, Minxiang Wei, Jiang Wu, Wenshuai Xue, Ranran Liu, An improved scheme based on the variable density method for structural topology optimization, 2025, 0954-4062, 10.1177/09544062251324422
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6012) PDF downloads(191) Cited by(23)

Figures and Tables

Figures(24)  /  Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog