Topology Optimization in Structural and Continuum Mechanics

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 deals with simultaneous topology and shape optimization of elastic contact problems. The structural optimization problem for an elastic contact problem is formulated. Shape as well as topological derivatives formulae of the cost functional are provided using material derivative and asymptotic expansion methods, respectively. These derivatives are employed to formulate necessary optimality condition for simultaneous shape and topology optimization and to calculate a descent direction in numerical algorithm. Level set based numerical algorithm for the solution of this optimization problem is proposed. Numerical examples are provided and discussed.

This work presents the framework to optimally design a cantilever for torsion mode frequency maximization. A cantilever design problem is formulated by topology and shape optimization methods. The torsion mode frequency is selected as an objective function, and the volume of the cantilever and the first bending mode frequency are constrained. Two optimization problems are defined and sequentially solved for the specific values. A new idea in this work is using a final topology obtained in the topology optimization problem as an initial shape in the shape optimization problem. The torsional mode frequency of the optimized cantilever is well improved in comparison with a nominal cantilever.

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.

According to the requirements of the aerospace industry for high strength, high stiffness, and lightweight structural parts, topology optimization has been proved to be an effective product design method. As one of the most conceptual and prospective structural optimization design methods, topology optimization intends to seek the optimal layout of materials in an allowed design region under a given load and boundary conditions. Thus, the object of study in the article is the method of topological optimization of aircraft structures. The goal of this article is to analyze the existing approaches, algorithms, as well as application of the method of topological optimization in the aerospace field in applied problems. The tasks are to describe the existing various approaches methods, features, and research directions of topological optimization as well as to study the possibility of application in the manufacturing process of composite structures. The following results were obtained. The optimization methods are briefly explained and compared, and the advantages and limitations of each approach are discussed. The various ways of simultaneous optimization of fiber orientation and structural topology were described and analyzed. The features of different methods of continuous fiber orientation optimization method were reviewed. The discrete fiber orientation optimization methods were represented. The possibility of multi-scale concurrent topological optimization was described. The combination of topology optimization and additive manufacturing was considered. Finally, the topology optimization of FRC structures which have been resolved in literature are reviewed and the potential research fields requiring more investigation are pointed out. Conclusions. In the article, a comprehensive review of the topology optimization design of FRC structures was presented. The promising way is to combine topology optimization with additive manufacturing techniques. However, these proposed methods may not suitable for other more complex problems, such as bucking stability and natural frequency. Hence, the topology optimization design of complex FRC components under complicated conditions is the main challenge in the future. This can be a new trend in the topology design of FRC structures.

This paper will give a short survey about topology optimization of structures. It is particularly focused on topological shape optimization of structures using level-set methods, including the level-set based standard methods and the level-set based alternative methods. The former often directly solve the Hamilton-Jacobi partial differential equation (H-J PDE) to obtain the boundary velocity field using Finite Differential Methods (FDM), and the later commonly employ parametric or equivalent methods to evaluate the velocity field without directly solving the H-J PDE. The unique characteristics of the level-set based topology optimization methods are discussed, and a future perspective and prospects in this research area is also included. A benchmark numerical example is used to showcase the effectiveness of the level-set based methods.

Structural design for obtaining desired eigenfrequencies by using the topology and shape optimization method

In recent years, additive manufacturing has been gaining ground among traditional manufacturing methods in many different applications. This is mainly due to its ease of prototyping, material savings and above all the highly complex geometries that can be achieved. Despite its many advantages, the use of additive manufacturing in electrical machines has been limited. The lack of maturity of the manufacturing processes and the increase in losses when applying additive approaches to some parts, mainly stators, have led little implementation in the field of electrical machines. The geometric freedom offered by additive manufacturing paves the way for the use of novel optimisation techniques, such as topology optimisation. Traditionally, topology optimisation has not been widely used because the complex geometries obtained by these algorithms are not easily produced by traditional manufacturing methods. However, thanks to additive manufacturing, topology optimisation has started to be used, mainly for mechanical applications and it has demonstrated its ability to reduce weight without compromising mechanical stability in several cases. Nevertheless, topology optimisation for physics other than mechanical has not been thoroughly studied, let alone for the simultaneous optimisation of different physics. In view of the scenario described, this thesis investigates the benefits that additive manufacturing and topology optimisation can bring to the design of electrical machines. For the additive manufacturing of electrical machines, a thorough literature review is carried out on the application of this manufacturing approach to different parts of electrical machines, active and non-active. Two main case studies are presented: in the first, the rotor, the shaft, the stator and the housing of an aerospace actuator are manufactured by L-PBF in FeCoV and the assembled actuator is tested. In the second case, a rotor is manufactured in FeSi by LP- DED. Finally, not considered an additive manufacturing case study itself in this thesis, but two electrical conductor prototypes are manufactured in CuCr1Zr and an additional geometry is manufactured in AlSi10Mg by L-PBF. With regard to topology optimisation of electrical machines, a description of the main topology optimisation methods used is given and their application to electrical machine components is presented. In a similar way to additive manufacturing, two case studies are analysed. First, a novel multiphysics - mechanical and electromagnetic - topology optimisation method is presented, in which the rotor and the stator of a permanent magnet motor are optimised simultaneously. This method is compared with two methods found in the literature, SIMP and on-off. Secondly, another topology optimisation method is described and applied to an electrical conductor model. The proposed algorithm involves a multiphysics - thermal and electromagnetic - hybrid parametric topology optimisation approach. Two reduced length geometries and one full length prototype are built via additive manufacturing. Finally, conclusions regarding additive manufacturing and topology optimisation of electrical machines are drawn from the work presented in this thesis. Future lines of work for additive manufacturing and topology optimisation of electrical machines are also presented.

Structural topology optimization is a crucial approach for achieving lightweight design. An effective topology optimization algorithm must strike a balance between the objective functions, constraints, and design variables, which essentially reflects the symmetry and tradeoff between the objective and constraints. In this study, a topology optimization method grounded in load path theory is proposed. Element bearing capacity is quantified using the element birth and death method, with an explicit formulation derived via finite element theory. The effectiveness in evaluating structural performance is assessed through comparisons with stress distributions and topology optimization density maps. In addition, a novel evaluation index for element bearing capacity is proposed as the objective function in the topology optimization model, which is validated through thin plate optimization. Subsequently, sensitivity redistribution mitigates checkerboard patterns, while mesh filtering suppresses multi-branch structures and prevents local optima. The method is applied for the lightweight design of a triangular arm, with results benchmarked against the variable density method, demonstrating the feasibility and effectiveness of the proposed method. The element bearing capacity seeks to homogenize the load distribution of each element; the technique in this study can be extended to the optimization of symmetric structures.

Stress-based topology optimization and nonlinear structural topology optimization is gaining increasing attention in order to make topology optimization more realistic. Thus, this paper extends current concepts of topology optimization to the design of structures made of nonlinear materials. An extended bi-directional evolutionary structural optimization (BESO) method for stress minimization topology optimization of material nonlinear structures is proposed in this work. BESO method based on discrete variables can effectively avoid the well-known singularity problem in density-based methods with low-density elements. The maximum von Mises stress is approximated by the p-norm global stress. The sensitivity information for designing variable updates is derived in detail by adjoint method. As for the highly nonlinear stress behavior, the updated scheme takes advantages from two filters respectively of the sensitivity and topological variables to improve convergence. Moreover, the filtered sensitivity numbers are combined with their historical sensitivity information to further stabilize the optimization process. The effectiveness of the proposed method is demonstrated by several 2D benchmark design problems.

Topology optimization is the procedure of determining the connectivity, shape and location of voids inside a given design domain. This chapter proposes a quantum-inspired evolutionary algorithm (QEA) for structural topology optimization. It reviews the formulation of the continuum structural topology optimization model and the characteristics of quantum computing. The chapter also proposes the detailed implementation of a QEA-based topology optimization. It discusses two issues: the density-based continuum structural topology optimization formulation and the characteristics of quantum computing, which may be conducive to reducing the computational expense of topology optimization. The entire theoretical procedure of quantum-inspired topology optimization, which mainly combines FEM analysis, design sensitivity solution and quantum-inspired search, is accomplished for continuum structures. As the experimental implementation of quantum computation needs initialization, coherent manipulation, control and read of the fragile quantum system, practically building quantum computers has proved extremely difficult.

The center of mass and the inertia tensor are important parameters characterizing the kinetic properties of many structural systems, e.g., the rotor, spacecraft, and aircraft. In these occasions, designs not consistent with the requirements on the center of mass and inertia tensor cannot be accepted. In this paper, mathematical models are established to carry out the topology optimization of continua considering the mass and inertia characteristics. Two numerical examples, one for minimizing the static compliance considering locations of the center of mass and the other for maximizing the fundamental eigenfrequency considering both locations of the center of mass and values of the inertia products, are presented to demonstrate that the requirements on mass and inertia characteristics could have a profound influence on the final optimized layouts in structural topology optimization problems. On the other hand, the mass and inertia characteristics can also be tuned by topology optimization methods. In the second numerical problem, the smooth minimum operator is used to approximate the fundamental eigenfrequency with the first 4 eigenfrequencies so that the possible non-differentiability of the repeated eigenfrequencies can be avoided. Due to its extreme importance in engineering practice, the work of this paper can be useful to carry out designs in accordance with the requirements on the mass and inertia characteristics under the framework of traditional density-based structural topology optimization methods.

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.

With the rapid advancement of computing technology, mathematical modeling and computer simulations have become common and standard methods for design of complex systems. To exercise these models intelligently and eliminate the burden of manual iteration, manipulating inputs and reviewing outputs, optimization strategies are applied in a simulation-based design environment. In reality, the knowledge that engineers have about the design problem is imperfect and incomplete, uncertainty always exists in simulation-based design. The conventional deterministic optimization techniques do not consider the impact of such uncertainties. If the system responses driving the design decisions are very sensitive to these uncertainties, serious problems might confront the designer. Performance of the actual system may differ greatly from expected optimum performance. Again, if the design solution is near one or more constraint boundaries, then even slight uncertainties or changes in the operating environment could produce failed, unsafe designs, and/or result in substantial performance degradation. In recent years, the concept of reliability and robust design has become very popular. Probabilistic design analysis and statistical modeling allow the identification of designs that qualify as not only feasible, but as consistently feasible in the face of uncertainty. Reliability based topology optimization (RBTO) extends the reliability notion to the area of structural topology optimization. Although reliability based design optimization (RBDO) is a well established research area, very few work has been presented in the area of RBTO. In this study, we give an overview of the various works done in the field of RBTO. Ever since Bendsoe and Kikuchi introduced the topology optimization using a homogenization method, many topology optimization methods have been developed for the both linear and nonlinear structures. Traditional topology optimization methods drive the topology of a structure to an optimum design based on the constraint on mass. Kharmanda et al. introduced reliability constraints in deterministic topology optimization problem. They proposed a heuristic strategy that aims to reduce mass while improving the reliability level of the structure without greatly increasing its weight. But the limit state function used by them was not based on failure criteria for the structure. Their formulation considered uncertainty with respect to ∗Design Automation Laboratory †Advisor geometrical dimension and applied load only. Also their reliability analysis seems to be independent of the boundary and loading condition, so their results showed similar values for the uncertain variables for different structures. Tovar et al. have developed the Hybrid Cellular Automaton (HCA) method for structural synthesis of continuum material where the state of each cell is defined by both density and strain energy. In Agarwal, a decoupled RBDO approach is employed such that the topology optimization is separate from the reliability analysis. Patel et al. showed the use of RBTO using the gradient free Hybrid Cellular Automata (HCA) method. Their formulation incorporates uncertainty with respect to material property also. They considered limit state function based on failure modes on the output displacements. This study aims at giving a comparative study on these methods. It points out the advantages and disadvantages of the methods, and their significance in the field of RBTO. As large numbers of design variables are associated with continuum topology optimization problems, RBTO methods are inherently computationally expensive because of additional required system analysis associate with RBDO. Based on the observations, this study presents an efficient method for incorporating reliability analysis in topology optimization.

This paper presents a performance index for topology and shape optimization of plate bending problems with displacement constraints. The performance index is developed based on the scaling design approach. This performance index is used in the Performance-Based Optimization (PBO) method for plates in bending to keep track of the performance history when inefficient material is gradually removed from the design and to identify optimal topologies and shapes from the optimization process. Several examples are provided to illustrate the validity and effectiveness of the proposed performance index for topology and shape optimization of bending plates with single and multiple displacement constraints under various loading conditions. The topology optimization and shape optimization are undertaken for the same plate in bending, and the results are evaluated by using the performance index. The proposed performance index is also employed to compare the efficiency of topologies and shapes produced by different optimization methods. It is demonstrated that the performance index developed is an effective indicator of material efficiency for bending plates. From the manufacturing and efficient point of view, the shape optimization technique is recommended for the optimization of plates in bending.