An efficient approach for reliability-based topology optimization

This article presents an efficient approach for reliability-based topology optimization (RBTO) in which the computational effort involved in solving the RBTO problem is equivalent to that of solving a deterministic topology optimization (DTO) problem. The methodology presented is built upon the bidirectional evolutionary structural optimization (BESO) method used for solving the deterministic optimization problem. The proposed method is suitable for linear elastic problems with independent and normally distributed loads, subjected to deflection and reliability constraints. The linear relationship between the deflection and stiffness matrices along with the principle of superposition are exploited to handle reliability constraints to develop an efficient algorithm for solving RBTO problems. Four example problems with various random variables and single or multiple applied loads are presented to demonstrate the applicability of the proposed approach in solving RBTO problems. The major contribution of this article comes from the improved efficiency of the proposed algorithm when measured in terms of the computational effort involved in the finite element analysis runs required to compute the optimum solution. For the examples presented with a single applied load, it is shown that the CPU time required in computing the optimum solution for the RBTO problem is 15–30% less than the time required to solve the DTO problems. The improved computational efficiency allows for incorporation of reliability considerations in topology optimization without an increase in the computational time needed to solve the DTO problem.     

Automated structural optimization system for integrated topology and shape optimization

Abstract

This paper combines previously developed techniques for image‐preprocessing and characteristic image‐interpreting together with a newly proposed automated shape‐optimization modeling technique into an integrated topology‐optimization and shape‐optimization system. As a result, structure designers are provided with an efficient and reliable automated structural optimization system (ASOS). The automated shape‐optimization modeling technique, the key technique in ASOS, uses hole‐expanding strategy, interference analysis, and hole shape‐adjusting strategy to automatically define the design variables and side constraints needed for shape optimization. This technique not only eliminates the need to manually define design variables and side constraints for shape optimization, but during the process of shape optimization also prevents interference between the interior holes and the exterior boundary. The ASOS is tested in three different structural configuration design examples.     

Multimaterial multijoint topology optimization

In this paper, a methodology that solves multimaterial topology optimization problems while also optimizing the quantity and type of joints between dissimilar materials is proposed. Multimaterial topology optimization has become a popular design optimization technique since the enhanced design freedom typically leads to superior solutions; however, the conventional assumption that all elements are perfectly fused together as a single piece limits the usefulness of the approach since the mutual dependency between optimal multimaterial geometry and optimal joint design is not properly accounted for. The proposed methodology uses an effective decomposition approach to both determine the optimal topology of a structure using multiple materials and the optimal joint design using multiple joint types. By decomposing the problem into two smaller subproblems, gradient‐based optimization techniques can be used and large models that cannot be solved with nongradient approaches can be solved. Moreover, since the joining interfaces are interpreted directly from multimaterial topology optimization results, the shape of the joining interfaces and the quantity of joints connecting dissimilar materials do not need to be defined a priori. Three numerical examples, which demonstrate how the methodology optimizes the geometry of a multimaterial structure for both compliance and cost of joining, are presented.     

结构拓扑优化中敏度过滤技术研究

为了抑制连续体结构拓扑优化结果中的棋盘格和灰度单元问题,借鉴粒子群优化算法中粒子状态的更新方法,提出一种改进的敏度更新技术.以结构的柔度最小为优化目标,构建了基于固体各项同性微惩罚结构的结构拓扑优化模型,根据结构的力学响应分析,采用优化准则法进行设计变量更新,进行载荷作用下二维连续体结构的拓扑优化设计,得到了材料在设计域内的最优分布.通过与已有敏度过滤技术的对比分析,验证了文中方法的正确性和有效性.     

Bi-directional evolutionary topology optimization for designing a neutrally buoyant underwater glider

An underwater glider is a type of autonomous profiling instrument platform used for gathering data to explore the ocean. Having a neutrally buoyant glider hull is one way to improve the glider's endurance with a passive compensation for buoyancy change. This article applies the bi-directional evolutionary structural optimization (BESO) method to the optimization of an underwater glider hull, based on two materials. Firstly, the method for determining the glider's neutral buoyancy is carried out. Secondly, the optimization problem is defined and the optimization procedure is presented. In the BESO procedure, the original design area elements with low strain energy are iteratively switched from high-value materials to low-value materials until a prescribed fraction is reached. Finally, an optimal underwater glider design is generated and the result demonstrates a reasonable material distribution of the neutrally buoyant glider hull. A 26.4% buoyancy adjustment is achieved and the mass of the glider is decreased by 31%.     

Geometrical design of thermoelectric generators based on topology optimization

This paper discusses an application of the topology optimization method for the design of thermoelectric generators. The proposed methodology provides the optimized geometry in accordance with various arbitrary conditions such as the types of materials, the volume of materials, and the temperature and shape of the installation position. By considering the coupled equations of state for the thermoelectric problem, we introduce an analytical model subject to these equations, which mimics the closed circuit composed of thermoelectric materials, electrodes, and a resistor. The total electric power applied to the resistor and the conversion efficiency are formulated as objective functions to be optimized. The proposed optimization method for thermoelectric generators is implemented as a geometrical optimization method using the solid isotropic material with penalization method used in topology optimizations. Simple relationships are formulated between the density function of the solid isotropic material with penalization method and the physical properties of the thermoelectric material. A sensitivity analysis for the objective functions is formulated with respect to the density function and the adjoint equations required for calculating it. Depending on the sensitivity, the density function is updated using the method of moving asymptotes. Finally, numerical examples are provided to demonstrate the validity of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Filters in topology optimization

In this article, a modified (‘filtered’) version of the minimum compliance topology optimization problem is studied. The direct dependence of the material properties on its pointwise density is replaced by a regularization of the density field by the mean of a convolution operator. In this setting it is possible to establish the existence of solutions. Moreover, convergence of an approximation by means of finite elements can be obtained. This is illustrated through some numerical experiments. The ‘filtering’ technique is also shown to cope with two important numerical problems in topology optimization, checkerboards and mesh dependent designs. Copyright © 2001 John Wiley & Sons, Ltd.     

Slope constrained topology optimization

The problem of minimum compliance topology optimization of an elastic continuum is considered. A general continuous density–energy relation is assumed, including variable thickness sheet models and artificial power laws. To ensure existence of solutions, the design set is restricted by enforcing pointwise bounds on the density slopes. A finite element discretization procedure is described, and a proof of convergence of finite element solutions to exact solutions is given, as well as numerical examples obtained by a continuation/SLP (sequential linear programming) method. The convergence proof implies that checkerboard patterns and other numerical anomalies will not be present, or at least, that they can be made arbitrarily weak. © 1998 John Wiley & Sons, Ltd.     

Level set topology optimization of structural problems with interface cohesion

This paper presents a finite element topology optimization framework for the design of two‐phase structural systems considering contact and cohesion phenomena along the interface. The geometry of the material interface is described by an explicit level set method, and the structural response is predicted by the extended finite element method. In this work, the interface condition is described by a bilinear cohesive zone model on the basis of the traction‐separation constitutive relation. The non‐penetration condition in the presence of compressive interface forces is enforced by a stabilized Lagrange multiplier method. The mechanical model assumes a linear elastic isotropic material, infinitesimal strain theory, and a quasi‐static response. The optimization problem is solved by a nonlinear programming method, and the design sensitivities are computed by the adjoint method. The performance of the presented method is evaluated by 2D and 3D numerical examples. The results obtained from topology optimization reveal distinct design characteristics for the various interface phenomena considered. In addition, 3D examples demonstrate optimal geometries that cannot be fully captured by reduced dimensionality. The optimization framework presented is limited to two‐phase structural systems where the material interface is coincident in the undeformed configuration, and to structural responses that remain valid considering small strain kinematics. Copyright © 2017 John Wiley & Sons, Ltd.     

Approximate reanalysis in topology optimization

In the nested approach to structural optimization, most of the computational effort is invested in the solution of finite element analysis equations. In this study, the integration of an approximate reanalysis procedure into the framework of topology optimization of continuum structures is investigated. The nested optimization problem is reformulated to accommodate the use of an approximate displacement vector and the design sensitivities are derived accordingly. It is shown that relatively rough approximations are acceptable since the errors are taken into account in the sensitivity analysis. The implementation is tested on several small and medium scale problems, including 2‐D and 3‐D minimum compliance problems and 2‐D compliant force inverter problems. Accurate results are obtained and the savings in computation time are promising. Copyright © 2008 John Wiley & Sons, Ltd.     

Simultaneous single-loop multimaterial and multijoint topology optimization

As the aerospace and automotive industries continue to strive for efficient lightweight structures, topology optimization (TO) has become an important tool in this design process. However, one ever-present criticism of TO, and especially of multimaterial (MM) optimization, is that neither method can produce structures that are practical to manufacture. Optimal joint design is one of the main requirements for manufacturability. This article proposes a new density-based methodology for performing simultaneous MMTO and multijoint TO. This algorithm can simultaneously determine the optimum selection and placement of structural materials, as well as the optimum selection and placement of joints at material interfaces. In order to achieve this, a new solid isotropic material with penalization-based interpolation scheme is proposed. A process for identifying dissimilar material interfaces based on spatial gradients is also discussed. The capabilities of the algorithm are demonstrated using four case studies. Through these case studies, the coupling between the optimal structural material design and the optimal joint design is investigated. Total joint cost is considered as both an objective and a constraint in the optimization problem statement. Using the biobjective problem statement, the tradeoff between total joint cost and structural compliance is explored. Finally, a method for enforcing tooling accessibility constraints in joint design is presented.     

A moving superimposed finite element method for structural topology optimization

Level set methods are becoming an attractive design tool in shape and topology optimization for obtaining efficient and lighter structures. In this paper, a dynamic implicit boundary‐based moving superimposed finite element method (s‐version FEM or S‐FEM) is developed for structural topology optimization using the level set methods, in which the variational interior and exterior boundaries are represented by the zero level set. Both a global mesh and an overlaying local mesh are integrated into the moving S‐FEM analysis model. A relatively coarse fixed Eulerian mesh consisting of bilinear rectangular elements is used as a global mesh. The local mesh consisting of flexible linear triangular elements is constructed to match the dynamic implicit boundary captured from nodal values of the implicit level set function. In numerical integration using the Gauss quadrature rule, the practical difficulty due to the discontinuities is overcome by the coincidence of the global and local meshes. A double mapping technique is developed to perform the numerical integration for the global and coupling matrices of the overlapped elements with two different co‐ordinate systems. An element killing strategy is presented to reduce the total number of degrees of freedom to improve the computational efficiency. A simple constraint handling approach is proposed to perform minimum compliance design with a volume constraint. A physically meaningful and numerically efficient velocity extension method is developed to avoid the complicated PDE solving procedure. The proposed moving S‐FEM is applied to structural topology optimization using the level set methods as an effective tool for the numerical analysis of the linear elasticity topology optimization problems. For the classical elasticity problems in the literature, the present S‐FEM can achieve numerical results in good agreement with those from the theoretical solutions and/or numerical results from the standard FEM. For the minimum compliance topology optimization problems in structural optimization, the present approach significantly outperforms the well‐recognized ‘ersatz material’ approach as expected in the accuracy of the strain field, numerical stability, and representation fidelity at the expense of increased computational time. It is also shown that the present approach is able to produce structures near the theoretical optimum. It is suggested that the present S‐FEM can be a promising tool for shape and topology optimization using the level set methods. Copyright © 2005 John Wiley & Sons, Ltd.     

A density field parametrization for topology optimization using Bernstein elements

A new way of describing the density field in density‐based topology optimization is introduced. The new method uses finite elements constructed from Bernstein polynomials rather than the more common Lagrange polynomials. Use of the Bernstein finite elements allows higher‐order elements to be used in the density‐field interpolation without producing unrealistic density values, ie, values lower than zero or higher than one. Results on several test problems indicate that using the higher‐order Bernstein elements produces optimal designs with sharper estimates of the optimal boundary on coarse design meshes. However, higher‐order elements are also required in the structural analysis to prevent the appearance of unrealistic material distributions. The Bernstein element density interpolation can be combined with adaptive mesh refinement to further improve design accuracy even on design domains with complex geometry.     

A fixed‐grid bidirectional evolutionary structural optimization method and its applications in tunnelling engineering

Topology optimization has exhibited an exceptional capability of improving structural design. However, several typical topology optimization algorithms are finite element (FE) based, where mesh‐dependent zigzag representation of boundaries is barely avoidable in both intermediate and final results. To tackle the problem, this paper proposes a new fixed‐grid‐based bidirectional evolutionary structural optimization method, namely FG BESO. The adoption of an FG FE framework enables a continuous boundary change in the course of topology optimization, which provides a means of dealing with not only the non‐smooth boundary of the final design but also the interpretation of intermediate densities. As a class of important practical application, it is interesting to make use of the new FG BESO method to the reinforcement design for underground tunnels. To accommodate the FG BESO technique to geological engineering applications, a nodal sensitivity is derived for a two‐phase material model comprising the artificial reinforcement and original rock. In this paper, some innovative topological designs of tunnel reinforcements are presented for minimizing the floor and sidewall heaves under different geological loading conditions. Copyright © 2007 John Wiley & Sons, Ltd.     

Smoothed finite element method for topology optimization involving incompressible materials

It is well known that the finite element method (FEM) suffers severely from the volumetric locking problem for incompressible materials in topology optimization owing to its numerical ‘overly stiff’ property. In this article, two typical smoothed FEMs with a certain softened effect, namely the node-based smoothed finite element method (NS-FEM) and the cell-based smoothed finite element method, are formulated to model the compressible and incompressible materials for topology optimization. Numerical examples have demonstrated that the NS-FEM with an ‘overly soft’ property is fairly effective in tackling the volumetric locking problem in topology optimization when both compressible and incompressible materials are involved.     

Checkerboard‐free topology optimization using non‐conforming finite elements

The objective of the present study is to show that the numerical instability characterized by checkerboard patterns can be completely controlled when non‐conforming four‐node finite elements are employed. Since the convergence of the non‐conforming finite element is independent of the Lamé parameters, the stiffness of the non‐conforming element exhibits correct limiting behaviour, which is desirable in prohibiting the unwanted formation of checkerboards in topology optimization. We employ the homogenization method to show the checkerboard‐free property of the non‐conforming element in topology optimization problems and verify it with three typical optimization examples. Copyright © 2003 John Wiley & Sons, Ltd.     

A mass constraint formulation for structural topology optimization with multiphase materials

This work is focused on the topology optimization of lightweight structures consisting of multiphase materials. Instead of adopting the common idea of using volume constraint, a new problem formulation with mass constraint is proposed. Meanwhile, recursive multiphase materials interpolation (RMMI) and uniform multiphase materials interpolation (UMMI) schemes are discussed and compared based on numerical tests and theoretical analysis. It is indicated that the nonlinearity of the mass constraint introduced by RMMI brings numerical difficulties to attain the global optimum of the optimization problem. On the contrary, the UMMI‐2 scheme makes it possible to formulate the mass constraint in a linear form with separable design variables. One such formulation favors very much the problem resolution by means of mathematical programming approaches, especially the convex programming methods. Moreover, numerical analysis indicates that fully uniform initial weighting is beneficial to seek the global optimum when UMMI‐2 scheme is used. Besides, the relationship between the volume constraint and mass constraint is theoretically revealed. The filtering technique is adapted to avoid the checkerboard pattern related to the problem with multiphase materials. Numerical examples show that the UMMI‐2 scheme with fully uniform initial weighting is reliable and efficient to deal with the structural topology optimization with multiphase materials and mass constraint. Meanwhile, the mass constraint formulation is evidently more significant than the volume constraint formulation. Copyright © 2011 John Wiley & Sons, Ltd.     

An adaptive mesh-adjustment strategy for continuum topology optimization to achieve manufacturable structural layout

This paper presents an adaptive mesh adjustment algorithm for continuum topology optimization method to describe the structural boundary using nonuniform isoparametric element. A criterion on the basis of the node movement is proposed; herein, the densities and coordinates of the nodes are defined to instruct the deformation of finite elements in subsequent optimization iterations. With such a scheme, the topology optimization can start from a regular mesh discretization then gradually yields an optimal design with clear structural boundaries. The element in the transition along the boundary is refined; on the contrary, the pure solid or void element is coarsen. The contribution of this work is to improve the resolution of the structural boundaries and decrease the percentage of transitional regions with the invariant design variable. Several 2D and 3D numerical examples indicate the effectiveness of our proposed method. Seen from the examples, the structural boundary become smoother and the intermediate densities have been reduced up to 70%. In addition, a design process based on the presented method is proposed to make the optimum solutions be fabricated conveniently and accurately by linking it with the 3D design software, ie, SolidWorks, which is also demonstrated in the numerical examples.     

A mixed FEM approach to stress‐constrained topology optimization

We present an alternative topology optimization formulation capable of handling the presence of stress constraints in a straightforward fashion. The main idea is to adopt a mixed finite‐element discretization scheme wherein not only displacements (as usual) but also stresses are the variables entering the formulation. By doing so, any stress constraint may be handled within the optimization procedure without resorting to post‐processing operation typical of displacement‐based techniques that may also cause a loss in accuracy in stress computation if no smoothing of the stress is performed. Two dual variational principles of Hellinger–Reissner type are presented in continuous and discrete form that, which included in a rather general topology optimization problem in the presence of stress constraints that is solved by the method of moving asymptotes (Int. J. Numer. Meth. Engng. 1984; 24 (3):359–373). Extensive numerical simulations are performed and ongoing extensions outlined, including the optimization of elastoplastic and incompressible media. Copyright © 2007 John Wiley & Sons, Ltd.     

The MINLP optimization approach to structural synthesis. Part I: A general view on simultaneous topology and parameter optimization

This paper presents the Mixed-Integer Non-linear Programming (MINLP) optimization approach to structural synthesis. Non-linear continuous/discrete non-convex problems of structural synthesis are proposed to be solved by means of simultaneous topology, parameter and standard dimension optimization. Part I of this three-part series of papers contains a general view of the MINLP approach to simultaneous topology and continuous parameter optimization. The MINLP optimization approach is performed through three steps. The first one includes the generation of a mechanical superstructure of different topology alternatives, the second one involves the development of an MINLP model formulation and the last one consists of a solution for the formulated MINLP problem. Some MINLP methods are also presented. A Modified OA/ER algorithm is applied to solve the MINLP problem and a simple example of a multiple cantilever beam is given to demonstrate the steps of the proposed MINLP optimization approach. As simultaneous optimization, extended to include also standard dimensions, requires additional effort, the development of suitable strategies to carry out the optimization is further discussed in Part II. The modelling of MINLP superstructures and the topology and parameter optimization of roller and sliding hydraulic steel gate structures are shown in Part III of the paper. An example of the synthesis of an already erected roller gate, i.e. the Intake Gate of Aswan II in Egypt, is presented as a comparative design research work. © 1998 John Wiley & Sons, Ltd.