A material based model for topology optimization of thermoelastic structures

This paper presents the development of a computational model for the topology optimization problem, using a material distribution approach, of a 2-D linear-elastic solid subjected to thermal loads, with a compliance objective function and an isoperimetric constraint on volume. Defining formally the augmented Lagrangian associated with the optimization problem, the optimality conditions are derived analytically. The results of analysis are implemented in a computer code to produce numerical solutions for the optimal topology, considering the temperature distribution independent of design. The design optimization problem is solved via a sequence of linearized subproblems. The computational model developed is tested in example problems. The influence of both the temperature and the finite element model on the optimal solution obtained is analysed.     

Temperature-constrained topology optimization of thermo-mechanical coupled problems

Temperature-constrained topology optimization for thermo-mechanical coupled problems under a design-dependent temperature field considering the thermal expansion effect remains an open problem. A temperature-constrained topology optimization method is proposed for thermo-mechanical coupled problems. In this article, the temperature values at the heat sources are constrained. The numerical results reveal that the temperature constraints play an important role in topology optimization of thermo-mechanical coupled problems. The optimized structure obtained by the presented method not only has certain strength but also decreases the temperature significantly compared with structures obtained by other methods without considering temperature constraints. The proposed method is applied to the design of the cooling system of a battery package. Numerical examples verify the efficiency of the presented method.     

Multi-objective optimization of engineering systems using game theory and particle swarm optimization

This study proposes particle swarm optimization (PSO) based algorithms to solve multi-objective engineering optimization problems involving continuous, discrete and/or mixed design variables. The original PSO algorithm is modified to include dynamic maximum velocity function and bounce method to enhance the computational efficiency and solution accuracy. The algorithm uses a closest discrete approach (CDA) to solve optimization problems with discrete design variables. A modified game theory (MGT) approach, coupled with the modified PSO, is used to solve multi-objective optimization problems. A dynamic penalty function is used to handle constraints in the optimization problem. The methodologies proposed are illustrated by several engineering applications and the results obtained are compared with those reported in the literature.     

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.     

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.     

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.     

博弈分析法在建筑设计方案优化中的应用

  建筑设计是开发商投资的关键阶段。根据我国建筑设计规范与法规体系的要求,借助博弈理论,分析建筑设计质量在设计方不同态度情况下的设计方案优化问题,采用博弈中的图解分析法,说明委托方与设计方之间的博弈,提出优化建筑设计方案的对策,从而使建筑设计方案达到最优。     

Continuous approximation of material distribution for topology optimization

In this paper, we propose a checkerboard‐free topology optimization method without introducing any additional constraint parameter. This aim is accomplished by the introduction of finite element approximation for continuous material distribution in a fixed design domain. That is, the continuous distribution of microstructures, or equivalently design variables, is realized in the whole design domain in the context of the homogenization design method (HDM), by the discretization with finite element interpolations. By virtue of this continuous FE approximation of design variables, discontinuous distribution like checkerboard patterns disappear without any filtering schemes. We call this proposed method the method of continuous approximation of material distribution (CAMD) to emphasize the continuity imposed on the ‘material field’. Two representative numerical examples are presented to demonstrate the capability and the efficiency of the proposed approach against some classes of numerical instabilities. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling topology optimization problems as linear mixed 0–1 programs

This paper deals with topology optimization of discretized continuum structures. It is shown that a large class of non‐linear 0–1 topology optimization problems, including stress‐ and displacement‐constrained minimum weight problems, can equivalently be modelled as linear mixed 0–1 programs. The modelling approach is applied to some test problems which are solved to global optimality. Copyright © 2003 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.     

On a cellular division method for topology optimization

This paper concerns a comparative analysis of a novel biologically inspired method for topology optimization. The proposed methodology develops each individual topology according to a set of rules that regulate a ‘cellular division’ process. These rules are then evolved using a genetic algorithm to minimize objective functions while satisfying a set of constraints. The results reported in this work show that the methodology suits engineering design and represents an improvement over existing topology optimization methods based on evolutionary algorithms. Copyright © 2011 John Wiley & Sons, Ltd.     

多供应链间的Nash和Stackelberg博弈决策模型

以N条包含一个制造商和一个零售商的单链式供应链为研究对象,分别讨论了多供应链间Nash博弈和Stackelberg博弈问题,对N条供应链间博弈均衡解进行了分析和求解,证明了N条供应链博弈下均衡解的存在性和唯一性。在实例分析中,以两条供应链为例,比较了不同博弈框架下定价策略和系统利润的差异。结果表明,对于供应链1来说,链间Nash博弈时的零售商利润、制造商利润和供应链利润均大于链间Stackelberg博弈情形,也即无论从零售商、制造商还是整条供应链的角度,他们均乐于放弃作为领导者的先动优势;对于供应链2来说,情况则相反,也即无论从零售商、制造商还是整条供应链的角度,他们均乐于作为追随者发挥后动优势。     

Multi-material topology optimization for the transient heat conduction problem using a sequential quadratic programming algorithm

Transient heat conduction analysis involves extensive computational cost. It becomes more serious for multi-material topology optimization, in which many design variables are involved and hundreds of iterations are usually required for convergence. This article aims to provide an efficient quadratic approximation for multi-material topology optimization of transient heat conduction problems. Reciprocal-type variables, instead of relative densities, are introduced as design variables. The sequential quadratic programming approach with explicit Hessians can be utilized as the optimizer for the computationally demanding optimization problem, by setting up a sequence of quadratic programs, in which the thermal compliance and weight can be explicitly approximated by the first and second order Taylor series expansion in terms of design variables. Numerical examples show clearly that the present approach can achieve better performance in terms of computational efficiency and iteration number than the solid isotropic material with penalization method solved by the commonly used method of moving asymptotes. In addition, a more lightweight design can be achieved by using multi-phase materials for the transient heat conductive problem, which demonstrates the necessity for multi-material topology optimization.     

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.     

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.     

Game formulations for structural optimization under uncertainty

We consider structural optimization (SO) under uncertainty formulated as a mathematical game between two players –– a “designer” and “nature”. The first player wants to design a structure that performs optimally, whereas the second player tries to find the worst possible conditions to impose on the structure. Several solution concepts exist for such games, including Stackelberg and Nash equilibria and Pareto optima. Pareto optimality is shown not to be a useful solution concept. Stackelberg and Nash games are, however, both of potential interest, but these concepts are hardly ever discussed in the literature on SO under uncertainty. Based on concrete examples of topology optimization of trusses and finite element-discretized continua under worst-case load uncertainty, we therefore analyze and compare the two solution concepts. In all examples, Stackelberg equilibria exist and can be found numerically, but for some cases we demonstrate nonexistence of Nash equilibria. This motivates a view of the Stackelberg solution concept as the correct one. However, we also demonstrate that existing Nash equilibria can be found using a simple so-called decomposition algorithm, which could be of interest for other instances of SO under uncertainty, where it is difficult to find a numerically efficient Stackelberg formulation.     

Structural topology optimization of vibrating structures with specified eigenfrequencies and eigenmode shapes

In vibration optimization problems, eigenfrequencies are usually maximized in the optimization since resonance phenomena in a mechanical structure must be avoided, and maximizing eigenfrequencies can provide a high probability of dynamic stability. However, vibrating mechanical structures can provide additional useful dynamic functions or performance if desired eigenfrequencies and eigenmode shapes in the structures can be implemented. In this research, we propose a new topology optimization method for designing vibrating structures that targets desired eigenfrequencies and eigenmode shapes. Several numerical examples are presented to confirm that the method presented here can provide optimized vibrating structures applicable to the design of mechanical resonators and actuators. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Parametric topology optimization with multiresolution finite element models

We present a methodical procedure for topology optimization under uncertainty with multiresolution finite element (FE) models. We use our framework in a bifidelity setting where a coarse and a fine mesh corresponding to low- and high-resolution models are available. The inexpensive low-resolution model is used to explore the parameter space and approximate the parameterized high-resolution model and its sensitivity, where parameters are considered in both structural load and stiffness. We provide error bounds for bifidelity FE approximations and their sensitivities and conduct numerical studies to verify these theoretical estimates. We demonstrate our approach on benchmark compliance minimization problems, where we show significant reduction in computational cost for expensive problems such as topology optimization under manufacturing variability, reliability-based topology optimization, and three-dimensional topology optimization while generating almost identical designs to those obtained with a single-resolution mesh. We also compute the parametric von Mises stress for the generated designs via our bifidelity FE approximation and compare them with standard Monte Carlo simulations. The implementation of our algorithm, which extends the well-known 88-line topology optimization code in MATLAB, is provided.