Shape optimization using time evolution equations

In this study, we deal with a numerical solution based on time evolution equations to solve the optimization problem for finding the shape that minimizes the objective function under given constraints. The design variables of the shape optimization problem are defined on a given original domain on which the boundary value problems of partial differential equations are defined. The variations of the domain are obtained by the time integration of the solution to derive the time evolution equations defined in the original domain. The shape gradient with respect to the domain variations are given as the Neumann boundary condition defined on the original domain boundary. When the constraints are satisfied, the decreasing property of the objective function is guaranteed by the proposed method. Furthermore, the proposed method is used to minimize the heat resistance under a total volume constraint and to solve the minimization problem of mean compliance under a total volume constraint.     

Topology optimization of thermoelastic structures using level set method

This paper describes the topology optimization of thermoelastic structures, using level set method. The objective is to minimize the mean compliance of a structure with a material volume constraint. In level set method, free boundary of a structure is considered as design variable, and it is implicitly represented via level set model. Objective function of the optimization problem is defined as a function of the shape of a structure. Sensitivity analysis based on continuum model is conducted with respect to the free boundary, which suggests the steepest descent direction. A geometric energy term is introduced to ensure smooth structural boundary. Augmented Lagrangian multiplier method is adopted to enforce volume constraint. Numerical examples are provided for 2D cases, considering design independent temperature distribution.     

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.     

随机参数平面连续体结构的频率拓扑优化设计

摘 要 研究几何和物理参数均为随机变量的平面连续体结构在结构基频约束下的拓扑优化设计问题。以结构总质量均值极小化为目标函数,以结构的形状拓扑信息为设计变量,以结构基频概率可靠性指标为约束条件,构建了随机结构拓扑优化设计数学模型。利用代数综合法,导出了随机参数结构动力响应的均值和均方差的计算表达式。采用渐进结构优化的求解策略与方法,通过两个算例验证了文中模型及求解方法的合理性和可行性。     

A projection-based method for topology optimization of structures with graded surfaces

Graded surfaces widely exist in natural structures and inspire engineers to apply functionally graded (FG) materials to cover structural surfaces for performance improvement, protection, or other special functionalities. However, how to design such structures with FG surfaces by topology optimization is a quite challenging problem due to the difficulty for determining material properties of structural surfaces with prescribed variation rule. This paper presents a novel projection-based method for topology optimization of this class of FG structures. Firstly, a projection process is proposed for ensuring the material properties of the surfaces vary with a prescribed function. A criterion of determining the values of parameters in projection process is given by a strict theoretical derivation, and then, a new interpolation function is established, which is capable of simultaneously obtaining clear substrate topologies and realizable FG surfaces. Though such structures are actually multimaterial gradient structures, only the design variables of single-material topology optimization problem are needed. In the current research, the classical compliance minimization problem with a mass constraint is considered and the robust formulation is used to control the length scale of substrates. Several 2D and 3D numerical examples illustrate the validity and applicability of the proposed method.     

考虑瞬态响应的薄板结构阻尼材料层拓扑优化

讨论了敷设阻尼材料的薄板结构考虑瞬态响应时阻尼材料层的最优布局问题。基于SIMP方法构造人工阻尼材料惩罚模型和结构拓扑优化模型,以阻尼材料的相对密度作为设计变量,在给定阻尼材料用量的条件下,最小化结构瞬态位移响应的时间积分。由于结构整体呈现非比例阻尼特性,采用逐步积分法对结构的振动方程进行求解。通过伴随变量法得到目标函数对设计变量的灵敏度表达式,在此基础上采用基于梯度的移动渐近线方法求解。数值算例验证了优化模型与算法的合理性和有效性。     

非平稳激励下薄板结构减振附加阻尼层的拓扑优化

讨论了附加阻尼层的薄板结构在非平稳随机力作用下以减振为目标的阻尼材料层的拓扑优化问题。建立了以阻尼材料的相对密度为设计变量,以结构非平稳响应位移方差最小化为目标和阻尼材料用量为约束条件的拓扑优化模型。由于结构受到非平稳随机激励作用,其随机响应可以采用时域显式法快速求解;随机响应方差对设计变量的灵敏度采用了基于伴随变量法的时域显式法进行分析,并采用优化准则法求解优化问题。数值算例验证了所提方法在非平稳随机激励作用下进行动力拓扑优化减振的可行性与有效性。     

Material nonlinear topology optimization considering the von Mises criterion through an asymptotic approach: Max strain energy and max load factor formulations

This paper addresses material nonlinear topology optimization considering the von Mises criterion by means of an asymptotic analysis using a fictitious nonlinear elastic model. In this context, we consider the topology optimization problem subjected to prescribed energy, which leads to robust convergence in nonlinear problems. Two nested formulations are considered. In the first, the objective is to maximize the strain energy of the system in equilibrium, and in the second, the objective is to maximize the load factor. In both cases, a volume constraint is imposed. The sensitivity analysis is quite effective and efficient in the sense that there is no extra adjoint equation. In addition, the nonlinear structural equilibrium problem is solved using direct minimization of the structural strain energy using Newton's method with an inexact line-search strategy. Four numerical examples demonstrate the features of the proposed material nonlinear topology optimization framework for approximating standard von Mises plasticity.     

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.     

Howard's algorithm in a phase‐field topology optimization approach

The topology optimization problem is formulated in a phase‐field approach. The solution procedure is based on the Allan–Cahn diffusion model where the conservation of volume is enforced by a global constraint. The functional defining the minimization problem is selected such that no penalization is imposed for full and void materials. Upper and lower bounds of the density function are enforced by infinite penalty at the bounds. A gradient term that introduces cost for boundaries and thereby regularizing the problem is also included in the objective functional. Conditions for stationarity of the functional are derived, and it is shown that the problem can be stated as a variational inequality or a max–min problem, both defining a double obstacle problem. The numerical examples used to demonstrate the method are solved using the FEM, whereas the double obstacle problem is solved using Howard's algorithm. Copyright © 2012 John Wiley & Sons, Ltd.     

Multiphase topology optimization with a single variable using the phase-field design method

In this study, a multimaterial topology optimization method using a single variable is proposed by combining the solid isotropic material with penalization method and the reaction-diffusion equation. Unlike ordinary multimaterial optimization, which requires several variables depending on the number of material types, this method intends to represent various materials as one variable. The proposed method combines two special functions in the sensitivity analysis of the objective function to converge the design variable into prespecified density values defined for each of the multimaterials. The composition constraint based on a normal distribution function is also introduced to estimate the distribution of each target density value in a single variable. It enables density exchange between multiple materials by increasing or decreasing the amount of a specific material. The proposed method is applied to structural and electromagnetic problems to verify its effectiveness, and its usefulness is also confirmed from the viewpoint of cost and computation time.     

基于密度分布曲线法的复合材料变刚度铺层优化

基于梯度的优化方法对复合材料层合板进行了变刚度铺层优化设计。在优化过程中需确定铺层中各单元的密度以及角度。为了使优化结果具有可制造性,优化结果需满足制造工艺约束并且铺层角度需从预定角度中选取。为了避免在优化问题中引入过多的约束并减少设计变量的数目,提出密度分布曲线法(DDCM)对层合板中各单元的密度进行参数化。根据各单元的密度以及角度设计变量并基于Bi-value Coding Parameterization(BCP)方法中的插值公式确定各单元的弹性矩阵。优化过程中以结构柔顺度作为优化目标,结构体积作为约束,优化算法采用凸规划对偶算法。对碳纤维复合材料的算例结果表明:采用DDCM可得到较理想的优化结果,并且收敛速率较快。     

A new efficient convergence criterion for reducing computational expense in topology optimization: reducible design variable method

A new efficient convergence criterion, named the reducible design variable method (RDVM), is proposed to save computational expense in topology optimization. There are two types of computational costs: one is to calculate the governing equations, and the other is to update the design variables. In conventional topology optimization, the number of design variables is usually fixed during the optimization procedure. Thus, the computational expense linearly increases with respect to the iteration number. Some design variables, however, quickly converge and some other design variables slowly converge. The idea of the proposed method is to adaptively reduce the number of design variables on the basis of the history of each design variable during optimization. Using the RDVM, those design variables that quickly converge are not considered as design variables for the next iterations. This means that the number of design variables can be reduced to save the computational costs of updating design variables. Then, the iteration will repeat until the number of design variables becomes 0. In addition, the proposed method can lead to faster convergence of the optimization procedure, which indeed is a more significant time saving. It is also revealed that the RDVM gives identical optimal solutions as those by conventional methods. We confirmed the numerical efficiency and solution effectiveness of the RDVM with respect to two types of optimization: static linear elastic minimization, and linear vibration problems with the first eigenvalue as the objective function for maximization. Copyright © 2012 John Wiley & Sons, Ltd.     

On the equivalence of optimality criterion and sequential approximate optimization methods in the classical topology layout problem

We study the ‘classical’ topology optimization problem, in which minimum compliance is sought, subject to linear constraints. Using a dual statement, we propose two separable and strictly convex subproblems for use in sequential approximate optimization (SAO) algorithms. Respectively, the subproblems use reciprocal and exponential intermediate variables in approximating the non‐linear compliance objective function. Any number of linear constraints (or linearly approximated constraints) are provided for. The relationships between the primal variables and the dual variables are found in analytical form. For the special case when only a single linear constraint on volume is present, we note that application of the ever‐popular optimality criterion (OC) method to the topology optimization problem, combined with arbitrary values for the heuristic numerical damping factor η proposed by Bendsøe, results in an updating scheme for the design variables that is identical to the application of a rudimentary dual SAO algorithm, in which the subproblems are based on exponential intermediate variables. What is more, we show that the popular choice for the damping factor η=0.5 is identical to the use of SAO with reciprocal intervening variables. Finally, computational experiments reveal that subproblems based on exponential intervening variables result in improved efficiency and accuracy, when compared to SAO subproblems based on reciprocal intermediate variables (and hence, the heuristic topology OC method hitherto used). This is attributed to the fact that a different exponent is computed for each design variable in the two‐point exponential approximation we have used, using gradient information at the previously visited point. Copyright © 2007 John Wiley & Sons, Ltd.     

Development of a simple and efficient method for robust optimization

Robust optimization problems are newly formulated and an efficient computational scheme is proposed. Both design variables and design parameters are considered as random variables about their nominal values. To ensure the robustness of objective performance, we introduce a new performance index bounding the performance together with a constraint limiting the performance variation. The constraint variations are regulated by considering the probability of feasibility. Each probability constraint is transformed into a sub‐optimization problem by the advanced first‐order second moment (AFOSM) method for computational efficiency. The proposed robust optimization method has the advantages that the mean value and the variation of the performance function are controlled simultaneously and rationally and the second‐order sensitivity information is not required even in case of gradient‐based optimization process. The suggested method is examined by solving three examples and the results are compared with those for the deterministic case and those available in the literature. Copyright © 2002 John Wiley & Sons, Ltd.     

Robust topology optimization of optical cloaks under uncertainties in wave number and angle of incident wave

This article presents a robust topology optimization method for optical cloaks under uncertainties in the wave number and angle in the incident wave. We first discuss the governing equation derived from Maxwell's equation, and extend it to the entire domain including the dielectric material and air, based on the level set-based topology optimization method. Next, a robust optimization problem is formulated as a minimization problem of the weighted sum of the scattered wave norm and its standard deviation with respect to the wave number and angle of the incident wave. The standard deviation is mathematically expressed by the Taylor series approximation and the use of the adjoint variable method. The design sensitivity of the objective functional is also derived by the adjoint variable method. An optimization algorithm is then constructed, based on the proposed formulation for robust designs of optical cloaks. Several numerical examples are finally provided to demonstrate the validity and utility of the proposed method.     

Optimality criteria-based topology optimization of a bi-material model for acoustic–structural coupled systems

This article investigates topology optimization of a bi-material model for acoustic–structural coupled systems. The design variables are volume fractions of inclusion material in a bi-material model constructed by the microstructure-based design domain method (MDDM). The design objective is the minimization of sound pressure level (SPL) in an interior acoustic medium. Sensitivities of SPL with respect to topological design variables are derived concretely by the adjoint method. A relaxed form of optimality criteria (OC) is developed for solving the acoustic–structural coupled optimization problem to find the optimum bi-material distribution. Based on OC and the adjoint method, a topology optimization method to deal with large calculations in acoustic–structural coupled problems is proposed. Numerical examples are given to illustrate the applications of topology optimization for a bi-material plate under a low single-frequency excitation and an aerospace structure under a low frequency-band excitation, and to prove the efficiency of the adjoint method and the relaxed form of OC.     

Truss layout design under nonprobabilistic reliability-based topology optimization framework with interval uncertainties

A nonprobabilistic reliability-based topology optimization (NRBTO) method for truss structures with interval uncertainties (or unknown-but-bounded uncertainties) is proposed in this paper. The cross-sectional areas of levers are defined as design variables, while the material properties and external loads are regard as interval parameters. A modified perturbation method is applied to calculate structural response bounds, which are the prerequisite to obtain structural reliability. A deviation distance between the current limit state plane and the objective limit state plane, of which the expression is explicit, is defined as the nonprobabilistic reliability index, which serves as a constraint function in the optimization model. Compared with the deterministic topology optimization problem, the proposed NRBTO formulation is still a single-loop optimization problem, as the reliability index is explicit. The sensitivity results are obtained from an analytical approach as well as a direct difference method. Eventually, the NRBTO problem is solved by a sequential quadratic programming method. Two numerical examples are used to testify the validity and effectiveness of the proposed method. The results show significant effects of uncertainties to the topology configuration of truss structures.     

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.     

概率及非概率不确定性条件下结构鲁棒设计方法

提出了在概率不确定性和非概率不确定性同时存在时的约束函数鲁棒性和目标函数鲁棒性的实现策略及结构鲁棒设计方法。将传统优化设计问题的约束条件改造成为能同时反映两类不确定性量波动变化影响的约束条件,以实现约束函数的鲁棒性;在传统优化设计问题目标函数中增加若干个关于目标函数灵敏度的新目标函数,构成一个多目标函数设计问题,以实现目标函数的鲁棒性。所提方法应用于一个10杆桁架结构设计,采用宽容排序法求解。计算结果表明,在相同的结构总质量限制条件下,目标函数鲁棒性程度随着变量不确定性程度的增加而降低;在相同的变量不确定性程度条件下,增加结构总质量能提高目标函数鲁棒性的程度。