An asymptotically concentrated method for structural topology optimization based on the SIMLF interpolation

This article proposes an efficient parameterized level set method (PLSM) for structural topology optimization design with minimum compliance as the objective function and volume fraction as the constraint condition. In the stage of establishing the optimization model, the level set function (LSF) is interpolated using compactly-supported radial basis functions (CS-RBFs) to transform the Hamilton-Jacobi partial differential equation (PDE) into ordinary differential equations (ODEs), thereby making the evolution of the LSF more convenient and efficient, and ensuring the smoothness of the optimization result boundary. The method of moving asymptotes (MMA) is used for solving the established optimization model, and meanwhile the shape sensitivity constraint factor is added to improve computational efficiency. During the evolution of the LSF, an approximate re-initialization scheme is employed to prevent the gradient of the LSF boundary from being too large or too small, thereby improving the numerical stability and the convergence speed of structural topology optimization process. Furthermore, the proposed method is also extensible and applicable to topology optimization of multi-material structures. The feasibility and effectiveness of this method have been verified through several typical numerical examples involving topology optimization of single-material structures and multi-material structures within the framework of minimum compliance design. HIGHLIGHTS An efficient parameterized level set method using the approximate re-initialization scheme is proposed for structural topology optimization. The approximate re-initialization scheme is employed when solving the parameterized level set function via the MMA algorithm. This scheme can prevents the gradient of the level set function boundary from being too large or too small, making the level set function update more stable and accelerating the convergence speed of structural topology optimization. The effectiveness and feasibility of the proposed method are demonstrated through examples of single-material and multi-material structural topology optimization.

This thesis aims at understanding and improving topology optimization techniques focusing on density-based level-set methods and geometrical nonlinearities. Central in this work are the numerical modeling of the mechanical response of a design and the consistency of the optimization process itself. Concerning the first topic, we investigate different means to improve the robustness of density-based numerical models including geometrical nonlinearities. The conventional approach (scaling the local material properties) can result in convergence problems due to excessive deformation in low-stiffness finite elements. To avoid excessive deformation, we combine the element connectivity parameterization method (adapting the connectivity between finite elements) with level-set-based topology optimization. Furthermore, we achieve greater robustness of analysis method using a second and improved approach called element deformation scaling. This approach eliminates the need for solving internal equilibrium equations (needed for the element connectivity parameterization method) via an explicit relation between the local internal and global external displacement field. The second focus of this thesis is on the optimization process of level-set-based topology optimization, and in particular its numerical consistency. We observe that signed-distance reinitialization of the level-set function affects the shape of a design in the optimization process. To minimize this effect, we propose a discrete level-set method that is based on an approximate Heaviside function and focusing on the implementation. Furthermore, we also propose a level-set-based topology optimization method using an exact Heaviside function and mathematical programming that effectively eliminates the need for reinitialization. We demonstrate that our density-based level-set method is closely related to conventional density-based topology optimization methods, while offering the advantage of more control over the geometrical complexity. On the other hand, we confirm that the dependence of the final result on the initial design which remains one of the big challenges for level-set-based topology optimization. The potential of the proposed level-set method is shown by applying it on problems with stress constraints and geometrical nonlinearities and performing manufacturing tolerant topology optimization. Finally, this thesis offers a review of level-set methods for structural topology optimization to identify and discuss the different approaches that are available in literature. We can distinguish between level-set methods by examination of their design parameterization, sensitivities, update procedures and regularization techniques. A level-set-based design parameterization offers the advantage of a crisp distinction between subdomains. For this reason, X-FEM approaches and conforming discretizations are an interesting option to retain the crisp nature of the level-set-based description of the design. Many level-set methods are combined with density-based numerical models and are, therefore, closely related to conventional density-based topology optimization methods. In particular, recently proposed projection methods have much in common with a level-set-based design description. The results of level-set-based TO methods often rely heavily on regularization techniques that introduce inconsistencies in the optimization process. Numerical consistency does not necessarily lead to the best search direction, but is essential to find a Karush-Kuhn-Tucker point of the discretized optimization problem.

During the last 20 years, structural optimization has become one of the most important topics of engineering applications. Design optimization of structure has been an interesting area of research in the field of engineering design for its ability to short the design cycle and to enhance product quality. Significant research activity has occurred in the area of structural optimization in the last decade. Especially for topology optimization of structure, many new theoretical, algorithmic, and computational contributions have resulted by researchers and engineers. Topology optimization is a powerful tool for global and multi-scale design of macrostructures, microstructures, and the cell of prescribed composite materials. The population based evolutionary algorithms have emerged as powerful mechanism for finding optimum solutions of complex optimization problems in engineering during the last two decades. Evolutionary computation is the study of computational systems which use ideas and get inspiration from natural evolution and adaptation [1]. The thinking has wide application in various engineering fields, such as computer science, artificial intelligence, operations research. Genetic algorithm is another kind of bio-inspired optimization method and it is playing an increasingly important role in studies of complex adaptive systems. Its application ranges from adaptive agents in economic theory to the use of machine learning techniques in the design of complex devices and structures, such as aircraft turbines and integrated circuits [2]. Optimization of structures can be classified into three categories: sizing, shaping, and topology optimization. In the topology optimization, it is concerned with the structure members and connectivity between members. In general, it is easily represented by discrete variables rather than by those used for continuous optimization problems. Topology optimization is the most difficult and complex among three categories and it is special useful in developing innovative conceptual designs. Structural optimization, in particular the topology optimization, has been identified as one of the most challenging tasks in structural design. Various techniques and approaches have been established during the last two decades. Topology optimization usually referred to as layout optimization or general shape optimization [3]. It lets engineers get the optimal topology of structure or new configurations during product design phase, as they are implementing the design of the size

This paper focuses on two kinds of bi-objective topology optimization problems with uniform-stress constraints: compliance-volume minimization and local frequency response–volume minimization problems. An adaptive volume constraint (AVC) algorithm based on an improved bisection method is proposed. Using this algorithm, the bi-objective uniform-stress-constrained topology optimization problem is transformed into a single-objective topology optimization problem and a volume-decision problem. The parametric level set method based on the compactly supported radial basis functions is employed to solve the single-objective problem, in which a self-organized acceleration scheme based on shape derivative and topological sensitivity is proposed to adaptively adjust the derivative of the objective function and the step length during the optimization. To solve the volume-decision problem, an improved bisection method is proposed. Numerical examples are tested to illustrate the feasibility and effectiveness of the self-organized acceleration scheme and the AVC algorithm based on the improved bisection method. An extended application to the bi-objective stress-constrained topology optimization of a structure with stress concentration is also presented.

This paper develops a non-probabilistic reliability-based topology optimization (NRBTO) framework for continuum structures under multi-dimensional convex uncertainties. Combined with the solid isotropic material with penalization (SIMP) model and the set-theoretical convex method, the uncertainty quantification (UQ) analysis is firstly conducted to obtain mathematical approximations and boundary laws of considered displacement responses. By normalization treatment of the limit-state function, a new quantified measure of the non-probabilistic reliability is then defined and further deduced by the principle of the hyper-volume ratio. For circumventing optimization difficulties arising from large-scale design variables, the adjoint vector scheme for sensitivity analysis of the reliability index with respect to design variables are discussed as well. Numerical applications eventually illustrate the applicability and the validity of the present problem statement as well as the proposed numerical techniques.

In this paper we consider the classical minimum compliance topology optimization in linear elasticity. In our approach, the objective is to minimize a cost functional that comprises the work of applied body and traction forces obtaining simultaneously the optimal topology and the optimal localization for proportional actuators on the concerning structural domain. A proportional actuator is considered as a spring and one of the objectives is to obtain the optimal local to couple it on the structure. The problem is treated as a topology optimization problem where a given cost functional must be minimized under suitable constraints for external forces, control energy, volume restriction and a set of project variables. The topology optimization procedure here performed uses the solid isotropic material with penalization approach, based on the concept of optimizing the material distribution. The solution of the optimization problem is obtained by using a sequential linear programming and the finite element method is used to discretize the design domain. Several numerical examples are presented showing the effectiveness of the proposed approach.

Topological optimization is an innovative method to realize the lightweight design. This paper proposes a hybrid topology optimization method that combines the SIMP (solid isotropic material with penalization) method and genetic algorithm (GA), called the SIMP-GA method. In the method, SIMP is used to update the chromosomes, which can accelerate convergence. The filtering scheme in the SIMP method can filter unconnected elements to ensure the connectivity of the structure. We studied the influence of varying the filtering radius on the optimized structure. Simultaneously, in the SIMP-GA method, each element is regarded as a gene, which controls the population number to a certain extent, reduces the amount of calculation, and improves the calculation efficiency. The calculation of some typical examples proves that the SIMP-GA method can obtain a better solution than the gradient-based method. Compared with the conventional genetic algorithm and GA-BESO (Bi-directional Evolutionary Structural Optimization) method, the calculation efficiency of the proposed method is higher and similar results are obtained. The innovative topology optimization method could be an effective way for structural lightweight design.

Structural topology optimization.- On basic properties of Michell's structures.- Validation of numerical method by analytical benchmarks and verification of exact solutions by numerical methods.- Introduction to shape and topology optimization.- Homogenization method for shape and topology optimization.- Level set method for shape and topology optimization.- Compliance minimization of two-material elastic structures.- Some notes on topology optimization of vibrating continuum structures.- Topology optimization of vibrating continuum structures for maximum values of simple and multiple eigenfrequencies and frequency gaps.- On optimum shape design and periodicity of band-gape beam structures.- Topological design for minimum dynamic compliance of continuum structures subjected to forced vibration.- Topological design for minimum sound emission from bi-material structures subjected to forced vibration.- Discrete material optimization of vibrating laminated composite plates for minimum sound emission.- Optimization of diffusive transport problems.- Fluid topology optimization: stokes and Navier-stokes models.- Topology optimization of coupled multi-physics problems.- Topology optimization based on the extended finite element method.- Topology optimization of meso- and nano-scale problems.- Topology optimization under uncertainty.- A brief review of numerical methods of structural topology optimization.- The free material design in planar elasticity.

This paper presents an efficient and comprehensive MATLAB code to solve two-dimensional structural topology optimization problems, including minimum mean compliance, compliant mechanism synthesis and multi-load compliance problems. The Unsmooth Variational Topology Optimization (UNVARTOP) method, developed by the authors in a previous work, is used in the topology optimization code, based on the finite element method (FEM), to compute the sensitivity and update the topology. The paper also includes instructions to improve the bisection algorithm, modify the computation of the Lagrangian multiplier by using an Augmented Lagrangian to impose the constraint, implement heat conduction problems and extend the code to three-dimensional topology optimization problems. The code, intended for students and newcomers in topology optimization, is included as an appendix (Appendix A) and it can be downloaded from https://github.com/DanielYago together with supplementary material.

SummaryA nonprobabilistic reliability‐based topology optimization (NRBTO) method for truss structures with interval uncertainties (or unknown‐but‐bounded uncertainties) is proposed in this paper. The cross‐sectional areas of levers are defined as design variables, while the material properties and external loads are regard as interval parameters. A modified perturbation method is applied to calculate structural response bounds, which are the prerequisite to obtain structural reliability. A deviation distance between the current limit state plane and the objective limit state plane, of which the expression is explicit, is defined as the nonprobabilistic reliability index, which serves as a constraint function in the optimization model. Compared with the deterministic topology optimization problem, the proposed NRBTO formulation is still a single‐loop optimization problem, as the reliability index is explicit. The sensitivity results are obtained from an analytical approach as well as a direct difference method. Eventually, the NRBTO problem is solved by a sequential quadratic programming method. Two numerical examples are used to testify the validity and effectiveness of the proposed method. The results show significant effects of uncertainties to the topology configuration of truss structures.

A new bionic topology optimization method by combining reaction-diffusion equations describing bone adaptation process with finite element analysis is presented in this study. The major idea of the present approach is to consider the structure to be optimized as a piece of bone that obeys bone adaptation and the process of finding the optimum topology of a structure is equivalent to the bone remodeling process. Two widely used numerical examples demonstrate that the proposed approach greatly improves numerical efficiency compared with the othert well known methods for structural topology optimization in open literature. The results show that the optimal designs from the present bionic topology optimization method without use mathematical programming and numerical instability control techniques. The proposed method results in a better and faster convergence.

Bone is a dynamic living tissue that undergoes continuous adaptation of its mass and structure in response to mechanical and biological environment demands. In this paper, we firstly propose a mathematical model based on cross-type reaction diffusion equations of bone adaptation during a remodeling cycle due to mechanical stimulus. The model captures qualitatively very well the bone adaptation and cell interactions during the bone remodeling. Secondly assuming the bone structure to be a self-optimizing biological material which maximizes its own structural stiffness, bone remodeling model coupled with finite element method by using the add and remove element a new topology optimization of continuum structure is presented. Two Numerical examples demonstrate that the proposed approach greatly improves numerical efficiency, compared with the others well known methods for structural topology optimization in open literatures.

Topology optimization for enhanced dynamic fracture resistance of structures

Topology optimization of multi-material structures considering anisotropic yield strengths

Discrete ground structure topology optimization design methods have to date received considerable attention in structural engineering. An alternative class of structural topology optimization methods, which have not yet received much attention in structural engineering, but which have undergone considerable development in the past decade, are the so-called continuum formulations. In this work, a continuum structural topology optimization formulation is presented and applied to the concept design optimization of structural bracing systems that are needed to stiffen tall structures against sidesway under lateral-wind and seismic-type loading. Although demonstrated here in the context of these specific design examples, continuum structural topology optimization methods are believed to hold potential as a design tool for a wide range of civil engineering type structures. A variety of continuum topology design formulations, including static compliance minimization and eigenvalue optimization, are explored, and solution parameters are varied to show that a number of design possibilities can be realized as solutions.