Continuum Topology Optimization for Concept Design of Frame Bracing Systems

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.

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

Bi-directional evolutionary topology optimization of continuum structures subjected to inertial loads

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.

Topology optimization with meshless density variable approximations and BESO method

This paper presents an adaptive continuation method for buckling topology optimization of continuum structures using the Solid Isotropic Material with Penalization (SIMP) model. For optimization problems of minimizing structural compliance subject to constraints on material volume and buckling load factors, it has been found that the conflict between the requirements for structural stiffness and stability may have an adverse impact on the performance of existing optimization algorithms. An automatic scheme for adjusting the penalization parameter is introduced to deal with this conflict and thus achieves better designs. Based on an investigation on the effect of the penalization parameter on design evolution during the optimization process, a rule is established to determine the appropriate penalization parameter values. Using this rule, an effective scheme is developed for determining the penalization parameter values such that the buckling constraints are appropriately considered throughout the optimization process. Numerical examples are presented to illustrate the effectiveness of the proposed method.

Topology optimization of multi-material structures considering anisotropic yield strengths

SUMMARYA method for topology optimization of continuum structures based on nodal density variables and density field mapping technique is investigated. The original discrete‐valued topology optimization problem is stated as an optimization problem with continuous design variables by introducing a material density field into the design domain. With the use of the Shepard family of interpolants, this density field is mapped onto the design space defined by a finite number of nodal density variables. The employed interpolation scheme has an explicit form and satisfies range‐restricted properties that makes it applicable for physically meaningful density interpolation. Its ability to resolve more complex spatial distribution of the material density within an individual element, as compared with the conventional elementwise design variable approach, actually provides certain regularization to the topology optimization problem. Numerical examples demonstrate the validity and applicability of the proposed formulation and numerical techniques. Copyright © 2011 John Wiley & Sons, Ltd.

A method for structural topology optimization of frame structures with flexible joints is presented. A typical frame structure is a set of beams and joints assembled to carry an applied load. The problem considered in this paper is to find the stiffest frame for a given mass. By introducing design variables for beams and joints, a mass distribution for optimal structural stiffness can be found. Each beam can have several design variables connected to its cross section. One of these is an area-type design variable which is used to represent the global size of the beam. The other design variables are of length ratio type, controlling the cross section of the beam. Joints are flexible elements connecting the beams in the structure. Each joint has stiffness properties and a mass. A framework for modelling these stiffnesses is presented and design variables for joints are introduced. We prove a theorem which can be interpreted as the fact that the removal of structural elements, e.g. joints or beams, can be modelled by a small strictly positive material amount assigned to the element. This is needed for the computations of sensitivities used in the applied gradient based iterative method. Both two and three dimensional problems, as well as multiple load cases and multiple mass constraints, are treated.

In order to solve the problem that the weight of large scale theodolite's casting or steel plate welded structure is too heavy, lightweight design of large scale turntable based on topology optimization of continuum structure is investigated. Topology optimization could have better optimization efficiency because it can provide the optimized distribution of material and gain optimized structure in conceptual design phase of the structural design. In order to obtain lightweight and high stiffness structure, a truss turntable based on topology optimization of continuum structure is presented, and the validity of this lightweight structure is examined by finite element analysis. Comparing to the traditional welding steel turntable, the lightweight turntable, without changing the static displacement, can lose 386.48 kg(lightweight ratio is 26.7%), and its first-order natural frequency raises from 167.02Hz to 173.71Hz. Meanwhile, in other aspects this turntable also has better performance than traditional structure.

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.

Topology optimization of continuum structures is a recent field in structural optimization. However, an increasing research activity in this area has been developed since the statement of the very first formulations. These formulations try to obtain the most adequate material distribution that satisfies the imposed structural limitations. The existence or absence of material in each part of the domain is usually defined by using a continuum variable (the relative density) in order to avoid dealing with a discrete optimization problem. This continuum approach of the material properties present important advantages since conventionaloptimization algorithms can be used. However, numerical models must be considered in order to develop the structural analysis for intermediate values of the relative densities. In this paper, we present some improvements in a minimum weight approach of the structural topology optimization problem. The main goal of this paper is to present an improved formulation that tries to reach binary 0-1 material distributions by using a continuum approach of the design variables. Furthermore, a perimeter penalization is included in the objective function to simplify the solutions obtained. In addition, some computational aspects are considered in order to reduce the computational effort. Finally, we compare the solutions obtained by using these formulations in two application examples.

A bit-array representation method for structural topology optimization using the GA is proposed. The importance of design connectivity is further emphasized and a hierarchical violation penalty method is proposed to penalize the violated constraint functions so that the problem of representation degeneracy can be overcome and the GA search can be driven towards the combination of better structural performance, less unusable material and fewer connected objects in the design domain. An identical initialization method is also proposed to test the performance of the GA operators. With the appropriately selected GA operators, the bit-array representation GA is applied to the structural topology optimization problems of minimum weight. Numerical results demonstrate that the present GA can achieve better accuracy with less computational cost and suggest that the GA performance can be significantly improved by handling the design connectivity properly.

A method for the non-probabilistic reliability optimization on frequency of continuum structures with uncertain-but-bounded parameters is proposed. The objective function is to maximize the non-probabilistic reliability index of frequency requirement.The corresponding bi-level optimization model is built, where the constraints are applied on the material volume in the outer loop and the limit state equation in the inner loop. The non-probabilistic reliability index of frequency requirement is derived by the analytical method for the continuum structure with the uncertain elastic module and mass density. Further, the sensitivity of the non-probabilistic reliability index with respect to the design variables is analyzed. The topology optimization in the outer loop is performed by a bi-directional evolutionary structural optimization (BESO) method, where the numerical techniques and the optimization procedure of BESO method are presented. Numerical results show that the proposed BESO method is efficient, and convergent optimal solutions can be achieved for a variety of optimization problems on frequency non-probabilistic reliability of continuum structures.