A parameterized level set method for structural topology optimization using the approximate re-initialization scheme

Frequency- and temperature-dependent properties of viscoelastic materials have a significant effect on vibration characteristics of constrained layer damping (CLD) structures. In this article, parametrized level set method (PLSM) is extended to solve the topology optimization problem of CLD plates, in which the viscoelastic material with frequency- and temperature-dependent characteristics is employed as damping core. The single/weighted modal loss factors, which are solved by an iterative method, are defined as objective functions, and topology optimization problem is formulated based on the finite element model for CLD plates which are described by combining parametric level set (PLS) functions. The sensitivities of objective functions with respect to expansion coefficients in PLS are calculated by the use of adjoint vector method. As the optimized results with complex constant shear modulus for viscoelastic core indicates, the proposed method is effective for topology optimization of CLD structures and insensitive to initial design. Numerical examples with the consideration of frequency- and temperature-dependent characteristics for “soft” and “hard” viscoelastic core are implemented. Through comparing the optimal results with different models for viscoelastic core, the influences of frequency- and temperature-dependent characteristics for viscoelastic core on natural frequencies, modal loss factors and optimized layouts are analyzed. Similarity indices are defined for distinguishing the difference among optimized layouts and initial design induced by modal strain energy contour. The optimal design dependence on layer thicknesses is further discussed with additional mass of CLD materials unchanged.

Non-probabilistic Reliability-based Topology Optimization (NRBTO) Scheme for Continuum Structures Based on the parameterized Level-Set method and Interval Mathematics

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

Topology optimization of multi-material structures considering anisotropic yield strengths

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.

ABSTRACTThis paper extends the energy‐based homogenization method (EBHM) to evaluate the effective elastic properties of two‐dimensional arc unit cells, where the periodic boundary conditions are constructed based on the polar strain–displacement relationships, and the displacement field is solved by the polar finite element method. Combining the extended EBHM with the parametric level set method (PLSM) and the augmented Lagrangian multiplier method (ALM), it is able to design the topological configurations of arc unit cells for seeking superior effective elastic properties. In addition, the level set function is first filtered by compactly supported radial basis functions (CSRBFs) in order to eliminate small holes in the configuration. Numerical examples are performed to demonstrate the advantages of the proposed method. The periodic boundary condition performs well when the circumferential angle of the arc unit cell is in the range of [0.0002, 0.2]. Results indicate that the effective elastic properties of optimized arc unit cells override those of geometrically transformed optimal square unit cells. This superiority mainly benefits from the fact that unisymmetrical configurations are allowed in arc unit cells but not in square unit cells. Results also imply that the capability of nucleation of new holes is retained even though small holes are eliminated by filtering. It contributes to the trade‐off between performance and manufacturability.

In conventional level set methods, the slope of the level set function needs to be well controlled to maintain the numerical stability during the topology optimization process. One common solution is to regularize the level set function to be a signed distance function, which is usually achieved by periodically implementing the so called re-initialization scheme to force the level set function to gain the desired signed distance property. However, the re-initialization scheme will bring some unwanted drawbacks to the optimization process, such as zero level set drifting, time consuming etc. In addition, re-initialization is usually implemented outside the optimization loop, which will cause convergence issues. In this paper, a distance regularization functional is introduced to the structural topology optimization objective functional to ensure the signed distance property of the level set function near the structure boundaries. This functional can also keep the level set function to be constant-value at positions far away from the structural boundaries. The radial basis function (RBF) based parameterization technique together with the mathematical programming are utilized to improve the potential capability of handling multiple constraints for the topology optimization. The combination of these two techniques makes the level set based topology optimization be capable of handling complicated multi-constrained problems with higher numerical efficiency, leaving no compromise to multiple drawbacks. To demonstrate the validity of the proposed scheme, benchmark examples on minimum compliance structural optimization are employed. This type of problem is computed by the conventional level set method with the introduced distance regularization functional, the RBF based parametric level set and at last, the distance regularized RBF based parametric level set separately to demonstrate their differences.

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.

SummaryThe parametric level set approach is an extension of the conventional level set methods for topology optimization. By parameterizing the level set function, level set methods can be directly coupled with mathematical programming to achieve better numerical robustness and computational efficiency. Moreover, the parametric level set scheme can not only inherit the primary advantages of the conventional level set methods, such as clear boundary representation and the flexibility in handling topological changes, but also alleviate some undesired features from the conventional level set methods, such as the need for reinitialization. However, in the existing radial basis function–based parametric level set method, it is difficult to identify the range of the design variables. Besides, the parametric level set evolution often struggles with large fluctuations during the optimization process. Those issues cause difficulties both in numerical stability and in material property mapping. In this paper, a cardinal basis function is constructed based on the radial basis function partition of unity collocation method to parameterize the level set function. The benefit of using cardinal basis function is that the range of the design variables can now be clearly specified as the value of the level set function. A distance regularization energy functional is also introduced, aiming to maintain the desired signed distance property during the level set evolution. With this desired feature, the level set evolution is stabilized against large fluctuations. In addition, the material properties mapped from the level set function to the finite element model can be more accurate.

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.

Topology optimization with meshless density variable approximations and BESO method

Level-set approaches are a family of domain classification techniques that rely on defining a scalar level-set function (LSF), then carrying out the classification based on the value of the function relative to one or more thresholds. Most continuum topology optimization formulations are at heart, a classification problem of the design domain into structural materials and void. As such, level-set approaches are gaining acceptance and popularity in structural topology optimization. In conventional level set approaches, finding an optimum LSF involves solution of a Hamilton-Jacobi system of partial differential equations with a large number of degrees of freedom, which in turn, cannot be accomplished without gradients information of the objective being optimized. A new approach is proposed in this paper where design variables are defined as the values of the LSF at knot points, then a Kriging model is used sto interpolate the LSF values within the rest of the domain so that classification into material or void can be performed. Perceived advantages of the Kriging-interpolated level-set (KLS) approach include alleviating the need for gradients of objectives and constraints, while maintaining a reasonable number of design variables that is independent from the mesh size. A hybrid genetic algorithm (GA) is then used for solving the optimization problem(s). An example problem of a short cantilever is studied under various settings of the KLS parameters in order to infer the best practice recommendations for tuning the approach. Capabilities of the approach are then further demonstrated by exploring its performance on several test problems.

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.

Robust topology optimization for multi-material structures considering material uncertainties

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