Level-set method for design of multi-phase elastic and thermoelastic materials

The problem of designing composite materials with desired mechanical properties is to specify the materials microstructures in terms of the topology and distribution of their constituent material phases within a unit cell of periodic microstructures. In this paper we present an approach based on a multi-phase level-set model for the geometric and material representation and for numerical solution of a least squares optimization problem. The level-set model precisely specifies the material regions and their sharp boundaries in contrast to a raster discretization of the conventional homogenization-based approaches. Combined with the classical shape derivatives, the level-set method yields a computational system of partial differential equations. In using the Eulerian computation scheme with a fixed rectilinear grid and a fixed mesh in the unit cell, the gradient descent solution of the optimization captures the interfacial boundaries naturally and performs topological changes accurately. The proposed method is illustrated with several 2D examples for the synthesis of heterogeneous microstructures of elastic and/or thermoelastic composites composed of two and three material phases.     

A level set method for shape and topology optimization of large‐displacement compliant mechanisms

A parameterization level set method is presented for structural shape and topology optimization of compliant mechanisms involving large displacements. A level set model is established mathematically as the Hamilton–Jacobi equation to capture the motion of the free boundary of a continuum structure. The structural design boundary is thus described implicitly as the zero level set of a level set scalar function of higher dimension. The radial basis function with compact support is then applied to interpolate the level set function, leading to a relaxation and separation of the temporal and spatial discretizations related to the original partial differential equation. In doing so, the more difficult shape and topology optimization problem is now fully parameterized into a relatively easier size optimization of generalized expansion coefficients. As a result, the optimization is changed into a numerical process of implementing a series of motions of the implicit level set function via an existing efficient convex programming method. With the concept of the shape derivative, the geometrical non‐linearity is included in the rigorous design sensitivity analysis to appropriately capture the large displacements of compliant mechanisms. Several numerical benchmark examples illustrate the effectiveness of the present level set method, in particular, its capability of generating new holes inside the material domain. The proposed method not only retains the favorable features of the implicit free boundary representation but also overcomes several unfavorable numerical considerations relevant to the explicit scheme, the reinitialization procedure, and the velocity extension algorithm in the conventional level set method. Copyright © 2008 John Wiley & Sons, Ltd.     

Radial basis functions and level set method for structural topology optimization

Level set methods have become an attractive design tool in shape and topology optimization for obtaining lighter and more efficient structures. In this paper, the popular radial basis functions (RBFs) in scattered data fitting and function approximation are incorporated into the conventional level set methods to construct a more efficient approach for structural topology optimization. RBF implicit modelling with multiquadric (MQ) splines is developed to define the implicit level set function with a high level of accuracy and smoothness. A RBF–level set optimization method is proposed to transform the Hamilton–Jacobi partial differential equation (PDE) into a system of ordinary differential equations (ODEs) over the entire design domain by using a collocation formulation of the method of lines. With the mathematical convenience, the original time dependent initial value problem is changed to an interpolation problem for the initial values of the generalized expansion coefficients. A physically meaningful and efficient extension velocity method is presented to avoid possible problems without reinitialization in the level set methods. The proposed method is implemented in the framework of minimum compliance design that has been extensively studied in topology optimization and its efficiency and accuracy over the conventional level set methods are highlighted. Numerical examples show the success of the present RBF–level set method in the accuracy, convergence speed and insensitivity to initial designs in topology optimization of two‐dimensional (2D) structures. It is suggested that the introduction of the radial basis functions to the level set methods can be promising in structural topology optimization. Copyright © 2005 John Wiley & Sons, Ltd.     

Level-set topology optimization for multimaterial and multifunctional mechanical metamaterials

Metamaterials are artificially engineered composites designed to have unusual properties. This article will develop a new level-set based topology optimization method for the computational design of multimaterial metamaterials with exotic thermomechanical properties. In order to generate metamaterials consisting of arrays of microstructures under periodicity, the numerical homogenization method is used to evaluate the effective properties of the microstructure, and a multiphase level-set model is used to evolve the boundaries of the multimaterial microstructure. The proposed method will produce material geometries with distinct interfaces and smoothed boundaries, which may facilitate the fabrication of the topologically optimized designs. Several numerical cases are used to demonstrate the effectiveness of the proposed method.     

Density and level set-XFEM schemes for topology optimization of 3-D structures

As the capabilities of additive manufacturing techniques increase, topology optimization provides a promising approach to design geometrically sophisticated structures. Traditional topology optimization methods aim at finding conceptual designs, but they often do not resolve sufficiently the geometry and the structural response such that the optimized designs can be directly used for manufacturing. To overcome these limitations, this paper studies the viability of the extended finite element method (XFEM) in combination with the level-set method (LSM) for topology optimization of three dimensional structures. The LSM describes the geometry by defining the nodal level set values via explicit functions of the optimization variables. The structural response is predicted by a generalized version of the XFEM. The LSM–XFEM approach is compared against results from a traditional Solid Isotropic Material with Penalization method for two-phase “solid–void” and “solid–solid” problems. The numerical results demonstrate that the LSM–XFEM approach describes crisply the geometry and predicts the structural response with acceptable accuracy even on coarse meshes.     

Level-set topology optimization for maximizing fracture resistance of brittle materials using phase-field fracture model

Fracture is one of the most common failure modes in brittle materials. It can drastically decrease material integrity and structural strength. To address this issue, we propose a level-set (LS) based topology optimization procedure to optimize the distribution of reinforced inclusions within matrix materials subject to the volume constraint for maximizing structural resistance to fracture. A phase-field fracture model is formulated herein to simulate crack initiation and propagation, in which a staggered algorithm is developed to solve such time-dependent crack propagation problems. In line with diffusive damage of the phase-field approach for fracture; topological derivatives, which provide gradient information for the topology optimization in a LS framework, are derived for fracture mechanics problems. A reaction-diffusion equation is adopted to update the LS function within a finite element framework. This avoids the reinitialization by overcoming the limitation to time step with the Courant-Friedrichs-Lewy condition. In this article, three numerical examples, namely, a L-shaped section, a rectangular slab with predefined cracks, and an all-ceramic onlay dental bridge (namely, fixed partial denture), are presented to demonstrate the effectiveness of the proposed LS based topology optimization for enhancing fracture resistance of multimaterial composite structures in a phase-field fracture context.     

Bidirectional evolutionary structural optimization for nonlinear structures under dynamic loads

This study aims to develop efficient numerical optimization methods for finding the optimal topology of nonlinear structures under dynamic loads. The numerical models are developed using the bidirectional evolutionary structural optimization method for stiffness maximization problems with mass constraints. The mathematical formulation of topology optimization approach is developed based on the element virtual strain energy as the design variable and minimization of compliance as the objective function. The suitability of the proposed method for topology optimization of nonlinear structures is demonstrated through a series of two- and three-dimensional benchmark designs. Several issues relating to the nonlinear structures subjected to dynamic loads such as material, geometric, and contact nonlinearities are addressed in the examples. It is shown that the proposed approach generates more reliable designs for nonlinear structures.     

Topology optimization of structures under multiple loading cases with a new compliance–volume product

This article proposes a new topology optimization method for the design of structures under multiple loading cases. The design is formulated as a multi-objective optimization problem by minimizing a new compliance–volume product, which optimizes the overall stiffness and volume simultaneously to avoid the empirical decision on design constraints and obtain an even lower structural volume. A normalized exponential weighted criterion (NEWC) method is included in the multi-objective optimization problem for the capture of the entire Pareto frontier. A weight evaluation method, in terms of the fuzzy multiple-attribute group decision-making (FMAGDM) theory, is incorporated in the problem to evaluate the weights of the objectives and guarantee the optimal design in an acceptable level. The solid isotropic material with penalty (SIMP) method is used to represent the dependence of elemental densities on material properties. Three typical numerical examples are employed to show the effectiveness of the proposed method.     

Simultaneous topology and fastener layout optimization of assemblies considering joint failure

This article provides a method for the simultaneous topology optimization of parts and their corresponding joint locations in an assembly. Therein, the joint locations are not discrete and predefined, but continuously movable. The underlying coupling equations allow for connecting dissimilar meshes and avoid the need for remeshing when joint locations change. The presented method models the force transfer at a joint location not only by using single spring elements but accounts for the size and type of the joints. When considering riveted or bolted joints, the local part geometry at the joint location consists of holes that are surrounded by material. For spot welds, the joint locations are filled with material and may be smaller than for bolts. The presented method incorporates these material and clearance zones into the simultaneously running topology optimization of the parts. Furthermore, failure of joints may be taken into account at the optimization stage, yielding assemblies connected in a fail-safe manner.     

Piecewise constant level set method for structural topology optimization

In this paper, a piecewise constant level set (PCLS) method is implemented to solve a structural shape and topology optimization problem. In the classical level set method, the geometrical boundary of the structure under optimization is represented by the zero level set of a continuous level set function, e.g. the signed distance function. Instead, in the PCLS approach the boundary is described by discontinuities of PCLS functions. The PCLS method is related to the phase‐field methods, and the topology optimization problem is defined as a minimization problem with piecewise constant constraints, without the need of solving the Hamilton–Jacobi equation. The result is not moving the boundaries during the iterative procedure. Thus, it offers some advantages in treating geometries, eliminating the reinitialization and naturally nucleating holes when needed. In the paper, the PCLS method is implemented with the additive operator splitting numerical scheme, and several numerical and procedural issues of the implementation are discussed. Examples of 2D structural topology optimization problem of minimum compliance design are presented, illustrating the effectiveness of the proposed method. Copyright © 2008 John Wiley & Sons, Ltd.     

Simultaneous partial topology and size optimization of a wing structure using ant colony and gradient based methods

This article presents a methodology and process for a combined wing configuration partial topology and structure size optimization. It is aimed at achieving a minimum structural weight by optimizing the structure layout and structural component size simultaneously. This design optimization process contains two types of design variables and hence was divided into two sub-problems. One is structure layout topology to obtain an optimal number and location of spars with discrete integer design variables. Another is component size optimization with continuous design variables in the structure FE model. A multi city-layer ant colony optimization (MCLACO) method is proposed and applied to the topology sub-problem. A gradient based optimization method (GBOM) built in the MSC.NASTRAN SOL-200 module was employed in the component size optimization sub-problem. For each selected layout of the wing structure, a size optimization process is performed to obtain the optimum result and feedback to the layout topology process. The numerical example shows that the proposed MCLACO method and a combination with the GBOM are effective for solving such a wing structure optimization problem. The results also indicate that significant structural weight saving can be achieved.     

Microstructural topology optimization for minimizing the sound pressure level of structural–acoustic coupled systems with multi-scale random parameters

The main aim of this article is to present a robust microstructural topology optimization methodology for structural–acoustic coupled systems with multi-scale random parameters. During the microstructural topology optimization, both the uncertainty at the macro-scale, which comes from the physical parameters of the acoustic medium or the external load, and the uncertainty existing in the constituent material properties of the microstructure at the micro-scale are considered as random parameters. A homogenization-based probabilistic finite element method (HPFEM) is first developed for quantifying the structural–acoustic system with multi-scale random parameters. The use of the HPFEM transforms the problem of microstructural topology optimization with multi-scale random parameters to an augmented deterministic microstructural topology optimization problem. This provides a computationally cheap alternative to Monte Carlo-based optimization algorithms. A numerical example of a hexahedral box is given to demonstrate the efficiency of the proposed method.     

Simultaneous optimal design of structural topology, actuator locations and control parameters for a plate structure

 Simultaneous optimization with respect to the structural topology, actuator locations and control parameters of an actively controlled plate structure is investigated in this paper. The system consists of a clamped-free plate, a H ₂ controller and four surface-bonded piezoelectric actuators utilized for suppressing the bending and torsional vibrations induced by external disturbances. The plate is represented by a rectangular design domain which is discretized by a regular finite element mesh and for each element the parameter indicating the presence or absence of material is used as a topology design variable. Due to the unavailability of large-scale 0–1 optimization algorithms, the binary variables of the original topology design problem are relaxed so that they can take all values between 0 and 1. The popular techniques in the topology optimization area including penalization, filtering and perimeter restriction are also used to suppress numerical problems such as intermediate thickness, checkerboards, and mesh dependence. Moreover, since it is not efficient to treat the structural and control design variables equally within the same framework, a nested solving approach is adopted in which the controller syntheses are considered as sub processes included in the main optimization process dealing with the structural topology and actuator locations. The structural and actuator variables are solved in the main optimization by the method of moving asymptotes, while the control parameters are designed in the sub optimization processes by solving the Ricatti equations. Numerical examples show that the approach used in this paper can produce systems with clear structural topology and high control performance. Received 16 November 2001 / Accepted 26 February 2002     

Robust topology optimization for cellular composites with hybrid uncertainties

This paper will develop a new robust topology optimization method for the concurrent design of cellular composites with an array of identical microstructures subject to random‐interval hybrid uncertainties. A concurrent topology optimization framework is formulated to optimize both the composite macrostructure and the material microstructure. The robust objective function is defined based on the interval mean and interval variance of the corresponding objective function. A new uncertain propagation approach, termed as a hybrid univariate dimension reduction method, is proposed to estimate the interval mean and variance. The sensitivity information of the robust objective function can be obtained after the uncertainty analysis. Several numerical examples are used to validate the effectiveness of the proposed robust topology optimization method.     

Geometry and topology optimization of sheet metal profiles by using a branch‐and‐bound framework

In this paper well established procedures from partial differential equation (PDE)‐constrained and discrete optimization are combined in a new way to find an optimal design of a multi‐chambered profile. Given a starting profile design, a load case and corresponding design constraints (e.g. sheet thickness, chamber sizes), the aim is to find an optimal subdivision into a predefined number of chambers with optimal shape subject to structural stiffness. In the presented optimization scheme a branch‐and‐bound tree is generated with one additional chamber in each level. Before adding the next chamber, the geometry of the profile is optimized. Then a relaxation of a topology optimization problem is solved. Based on this relaxation, a best fitting feasible topology subject to manufacturability conditions is determined using a new mixed integer method employing shortest paths. To improve the running time, the finite element simulations for the geometry optimization and topology relaxation are performed with different levels of accuracy. Finally, numerical experiments are presented including different starting geometries, load scenarios and mesh sizes.     

连续体结构的变密度拓扑优化方法研究

为了实现使连续体结构的体积约束和柔顺度最小的拓扑优化及解决采用经典变密度法引起的结构优化结果存在如灰度单元、棋盘格等数值不稳定问题，提出了一种新的拓扑优化方法。首先，采用改进的固体各向同性材料惩罚法作为材料插值方案，建立结构拓扑优化模型；其次，通过引入基于高斯权重函数的敏度过滤法和设计新灰度单元抑制算子来解决数值不稳定问题；最后，借助优化准则法求解优化模型。通过算例分析可知：新策略可以改进拓扑优化方法；新的拓扑优化方法具有收敛速度较快、能更好地获取柔顺度小且拓扑构型好的优化结构和抑制灰度单元产生等优势。研究结果为其他连续体结构的拓扑优化研究提供了新思路。     

Topology optimization for stationary fluid–structure interaction problems using a new monolithic formulation

This paper outlines a new procedure for topology optimization in the steady‐state fluid–structure interaction (FSI) problem. A review of current topology optimization methods highlights the difficulties in alternating between the two distinct sets of governing equations for fluid and structure dynamics (hereafter, the fluid and structural equations, respectively) and in imposing coupling boundary conditions between the separated fluid and solid domains. To overcome these difficulties, we propose an alternative monolithic procedure employing a unified domain rather than separated domains, which is not computationally efficient. In the proposed analysis procedure, the spatial differential operator of the fluid and structural equations for a deformed configuration is transformed into that for an undeformed configuration with the help of the deformation gradient tensor. For the coupling boundary conditions, the divergence of the pressure and the Darcy damping force are inserted into the solid and fluid equations, respectively. The proposed method is validated in several benchmark analysis problems. Topology optimization in the FSI problem is then made possible by interpolating Young's modulus, the fluid pressure of the modified solid equation, and the inverse permeability from the damping force with respect to the design variables. Copyright © 2009 John Wiley & Sons, Ltd.     

Structural shape and topology optimization of cast parts using level set method

This paper presents a level set‐based shape and topology optimization method for conceptual design of cast parts. In order to be successfully manufactured by the casting process, the geometry of cast parts should satisfy certain moldability conditions, which poses additional constraints in the shape and topology optimization of cast parts. Instead of using the originally point‐wise constraint statement, we propose a casting constraint in the form of domain integration over a narrowband near the material boundaries. This constraint is expressed in terms of the gradient of the level set function defining the structural shape and topology. Its explicit and analytical form facilitates the sensitivity analysis and numerical implementation. As compared with the standard implementation of the level set method based on the steepest descent algorithm, the proposed method uses velocity field design variables and combines the level set method with the gradient‐based mathematical programming algorithm on the basis of the derived sensitivity scheme of the objective function and the constraints. This approach is able to simultaneously account for the casting constraint and the conventional material volume constraint in a convenient way. In this method, the optimization process can be started from an arbitrary initial design, without the need for an initial design satisfying the cast constraint. Numerical examples in both 2D and 3D design domain are given to demonstrate the validity and effectiveness of the proposed method. Copyright © 2017 John Wiley & Sons, Ltd.     

Robust topology optimization under load position uncertainty

The topology optimization problem of a continuum structure on the compliance minimization objective is investigated under consideration of the external load uncertainty in its application position with a nonprobabilistic approach. The load position is defined as the uncertain-but-bounded parameter and is represented by an interval variable with a nominal application point. The structural compliance due to the load position deviation is formulated with the quadratic Taylor series expansion. As a result, the objective gradient information to the topological variables can be evaluated efficiently in a quadratic expression. Based on the maximum design sensitivity value, which corresponds to the most sensitive compliance to the uncertain loading position, a single-level optimization approach is suggested by using a popular gradient-based optimality criteria method. The proposed optimization scheme is performed to gain the robust topology optimizations of three benchmark examples, and the final configuration designs are compared comprehensively with the conventional topology optimizations under the loading point fixation. It can be observed that the present method can provide remarkably different material layouts with auxiliary components to accommodate the load position disturbances. The numerical results of the representative examples also show that the structural performances of the robust topology optimizations appear less sensitive to the load position perturbations than the traditional designs.     

New conceptual design of aeroelastic wing structures by multi-objective optimization

Internal structural layouts and component sizes of aircraft wing structures have a significant impact on aircraft performance such as aeroelastic characteristics and mass. This work presents an approach to achieve simultaneous partial topology and sizing optimization of a three-dimensional wing-box structure. A multi-objective optimization problem is assigned to optimize lift effectiveness, buckling factor and mass of a structure. Design constraints include divergence and flutter speeds, buckling factor and stresses. The topology and sizing design variables for wing internal components are based on a ground element approach. The design problem is solved by multi-objective population-based incremental learning (MOPBIL). The Pareto optimum results lead to unconventional wing structures that are superior to their conventional counterparts.