Level-set topology optimization for robust design of structures under hybrid uncertainties

This paper will develop a new robust topology optimization (RTO) method based on level sets for structures subject to hybrid uncertainties, with a more efficient Karhunen-Loève hyperbolic Polynomial Chaos–Chebyshev Interval method to conduct the hybrid uncertain analysis. The loadings and material properties are considered hybrid uncertainties in structures. The parameters with sufficient information are regarded as random fields, while the parameters without sufficient information are treated as intervals. The Karhunen-Loève expansion is applied to discretize random fields into a finite number of random variables, and then, the original hybrid uncertainty analysis is transformed into a new process with random and interval parameters, to which the hyperbolic Polynomial Chaos–Chebyshev Interval is employed for the uncertainty analysis. RTO is formulated to minimize a weighted sum of the mean and standard variance of the structural objective function under the worst-case scenario. Several numerical examples are employed to demonstrate the effectiveness of the proposed RTO, and Monte Carlo simulation is used to validate the numerical accuracy of our proposed method.     

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.     

Robust topology optimization of phononic crystals with random field uncertainty

The uncertain spatial variation of material properties can remarkably affect the band gap characteristics of phononic crystals (PnCs). It is necessary to consider this issue when designing and manufacturing PnC materials/structures. This paper investigates a robust topology optimization method for designing the microstructures of PnCs by considering random‐field material properties. Herein, the spatial distribution of the material properties is first represented by a random field and then discretized into uncorrelated stochastic variables with the expansion optimal linear estimation method; stochastic band gap analysis is then conducted with polynomial chaos expansion. Furthermore, a robust topology optimization formulation of PnCs is proposed on the basis of the relative elemental density, where a weighted objective function handles the compromise of the mean value and standard deviation of the PnC band gap. The band gap response is analyzed, employing the finite element method for each sample of polynomial chaos expansion. In this context, the sensitivities of the stochastic band gap behaviors to the design variables are also derived. Numerical examples demonstrate that the proposed method can generate meaningful optimal topologies of PnCs with a relatively large width and less sensitive band gap. Additionally, the effects of the weight factors in the objective function and the variation coefficient of material properties are discussed.     

A new level set method for topology optimization of distributed compliant mechanisms

A new level set method for topology optimization of distributed compliant mechanism is presented in this study. By taking two types of mean compliance into consideration, several new objective functions are developed and built in the conventional level set method to avoid generating the de facto hinges in the created mechanisms. Aimed at eliminating the costly reinitialization procedure during the evolution of the level set function, an accelerated level set evolution algorithm is developed by adding an extra energy function, which can force the level set function to close to a signed distance function during the evolution. Two widely studied numerical examples in topology optimization of compliant mechanisms are studied to demonstrate the effectiveness of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Robust topology optimization of continuum structures with loading uncertainty using a perturbation method

An approach for the robust topology optimization (RTO) of continuum structures with loading uncertainty is investigated. The loading uncertainties are quantified using the second order Taylor series expansion of uncertain loading magnitudes and directions, and then the response statistic mean and standard deviation of compliance are calculated using the uncertain perturbation propagation method. A robust design Lagrange function considering the compliance objective and finite element constraints is developed, and a sensitivity analysis is performed to calculate the Lagrange coefficients. The Lagrange objective function is optimized using the modified solid isotropic material with penalization (SIMP) algorithm; thus, the optimum material distribution under loading uncertainty is acquired. The proposed methodology is used for the RTO of two examples, revealing its efficiency under both concentrated and distributed uncertain loadings. The accuracy of the results is verified by comparison with similar cases found in the literature where a different modelling approach was used.     

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.     

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

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

基于损失函数的随机设施布局优化分析

现有的随机设施布局问题的优化方法或是优化总物料处理成本的均值以获得最鲁棒的布局,或是优化其方差以获得最稳定的布局,但两者均只优化了总物料处理成本随机特征中的一个方面.通过引入田口损失函数的概念,建立了新的随机设施布局优化模型,对总物料处理成本的均值和方差同时进行优化求解,相比于现有的优化方法,新方法得到的布局形式能够使系统同时获得较好的鲁棒性和稳定性.     

拓扑相关荷载作用下结构拓扑优化的水平集方法

发展了一种利用水平集演化技术求解拓扑相关荷载作用下结构拓扑优化问题的数值方法。通过引入水平集函数,我们以隐含的方式对结构的拓扑和形状作了描述,从而把拓扑优化问题转化为了寻求最优水平集函数的数学规划问题。利用基于连续体概念的灵敏度分析技术,构造了用于驱动水平集演化的速度场。由于结构的边界可以用零水平集加以描述,因此利用适当的数学变换,我们可以方便地处理施加在结构上的拓扑相关荷载,这样就避免了以往算法中繁复的边界提取工作以及为了处理拓扑相关荷载所采取的特殊技巧。文末的数值算例表明了提出的优化方法在处理此类问题时所具有的独到的优越性。     

基于混合PSO算法的桁架动力响应优化

摘 要：本文针对以结构动力响应为约束,最小重量为目标的桁架拓扑优化问题,提出了一种将微粒群算法和优化准则法结合的混合PSO算法。利用优化准则法的迭代关系找出群体中适应度最好的微粒,将其作为特殊微粒,其他微粒的寻优采用PSO的基本进化规则,位移响应约束利用特殊微粒的灵敏度信息近似计算。算例的计算结果表明,混合PSO算法适用于受简谐荷载以及脉冲荷载作用桁架结构的拓扑优化。混合PSO的计算效率比PSO算法高,其优化效果比优化准则法好。     

A surrogate assisted parallel multiobjective evolutionary algorithm for robust engineering design

A number of multi-objective evolutionary algorithms have been proposed in recent years and many of them have been used to solve engineering design optimization problems. However, designs need to be robust for real-life implementation, i.e. performance should not degrade substantially under expected variations in the variable values or operating conditions. Solutions of constrained robust design optimization problems should not be too close to the constraint boundaries so that they remain feasible under expected variations. A robust design optimization problem is far more computationally expensive than a design optimization problem as neighbourhood assessments of every solution are required to compute the performance variance and to ensure neighbourhood feasibility. A framework for robust design optimization using a surrogate model for neighbourhood assessments is introduced in this article. The robust design optimization problem is modelled as a multi-objective optimization problem with the aim of simultaneously maximizing performance and minimizing performance variance. A modified constraint-handling scheme is implemented to deal with neighbourhood feasibility. A radial basis function (RBF) network is used as a surrogate model and the accuracy of this model is maintained via periodic retraining. In addition to using surrogates to reduce computational time, the algorithm has been implemented on multiple processors using a master–slave topology. The preliminary results of two constrained robust design optimization problems indicate that substantial savings in the actual number of function evaluations are possible while maintaining an acceptable level of solution quality.     

考虑混合约束的柔顺机构拓扑优化设计

为了满足制造工艺和静强度要求,提出一种综合考虑最小尺寸控制和应力约束的柔顺机构混合约束拓扑优化设计方法。采用改进的固体各向同性材料插值模型描述材料分布,利用多相映射方法同时控制实相和空相材料结构的最小尺寸,采用最大近似函数P范数求解机构的最大应力,以机构的输出位移最大化作为目标函数,综合考虑最小特征尺寸控制和应力约束建立柔顺机构混合约束拓扑优化数学模型,利用移动渐近算法求解柔顺机构混合约束拓扑优化问题。数值算例结果表明,混合约束拓扑优化获得的柔顺机构能够同时满足最小尺寸制造约束和静强度要求,机构的von Mises等效应力分布更加均匀。     

考虑混合约束的柔顺机构拓扑优化设计

为了满足制造工艺和静强度要求,提出一种综合考虑最小尺寸控制和应力约束的柔顺机构混合约束拓扑优化设计方法。采用改进的固体各向同性材料插值模型描述材料分布,利用多相映射方法同时控制实相和空相材料结构的最小尺寸,采用最大近似函数P范数求解机构的最大应力,以机构的输出位移最大化作为目标函数,综合考虑最小特征尺寸控制和应力约束建立柔顺机构混合约束拓扑优化数学模型,利用移动渐近算法求解柔顺机构混合约束拓扑优化问题。数值算例结果表明,混合约束拓扑优化获得的柔顺机构能够同时满足最小尺寸制造约束和静强度要求,机构的von Mises等效应力分布更加均匀。     

Robust concurrent topology optimization of multiscale structure under single or multiple uncertain load cases

Concurrent topology optimization of macrostructure and material microstructure has attracted significant interest in recent years. However, most of the existing works assumed deterministic load conditions, thus the obtained design might have poor performance in practice when uncertainties exist. Therefore, it is necessary to take uncertainty into account in structural design. This article proposes an efficient method for robust concurrent topology optimization of multiscale structure under single or multiple load cases. The weighted sum of the mean and standard deviation of the structural compliance is minimized and constraints are imposed to both the volume fractions of macrostructure and microstructure. The effective properties of the microstructure are calculated via the homogenization method. An efficient sensitivity analysis method is proposed based on the superposition principle and orthogonal similarity transformation of real symmetric matrices. To further reduce the computational cost, an efficient decoupled sensitivity analysis method for microscale design variables is proposed. The bidirectional evolutionary structural optimization method is employed to obtain black and white designs for both macrostructure and microstructure. Several two-dimensional and three-dimensional numerical examples are presented to demonstrate the effectiveness of the proposed approach and the effects of load uncertainty on the optimal design of both macrostructure and microstructure.     

Solution to the topology optimization problem using a time-evolution equation

This article describes a numerical solution to the topology optimization problem using a time-evolution equation. The design variables of the topology optimization problem are defined as a mathematical scalar function in a given design domain. The scalar function is projected to the normalized density function. The adjoint variable method is used to determine the gradient defined as the ratio of the variation of the objective function or constraint function to the variation of the design variable. The variation of design variables is obtained using the solution of the time-evolution equation in which the source term and Neumann boundary condition are given as a negative gradient. The distribution of design variables yielding an optimal solution is obtained by time integration of the solution of the time-evolution equation. By solving the topology optimization problem using the proposed method, it is shown that the objective function decreases when the constraints are satisfied. Furthermore, we apply the proposed method to the thermal resistance minimization problem under the total volume constraint and the mean compliance minimization problem under the total volume constraint.     

基于变频率区间约束的结构材料优化设计

针对频率约束的结构材料优化问题,基于结构拓扑优化思想,提出变频率区间约束的结构材料优化方法。借鉴均匀化及ICM(独立、连续、映射)方法,以微观单元拓扑变量倒数为设计变量,导出宏观单元等效质量矩阵及导数,进而获得频率一阶近似展开式。结合变频率区间约束思想,获得以结构质量为目标函数、频率为约束条件的连续体微结构拓扑优化近似模型;采用对偶方法求解。通过算例验证该方法的有效性及可行性,表明考虑质量矩阵变化影响所得优化结果更合理。     

Topology optimization of fluid problems using genetic algorithm assisted by the Kriging model

A non‐gradient‐based approach for topology optimization using a genetic algorithm is proposed in this paper. The genetic algorithm used in this paper is assisted by the Kriging surrogate model to reduce computational cost required for function evaluation. To validate the non‐gradient‐based topology optimization method in flow problems, this research focuses on two single‐objective optimization problems, where the objective functions are to minimize pressure loss and to maximize heat transfer of flow channels, and one multi‐objective optimization problem, which combines earlier two single‐objective optimization problems. The shape of flow channels is represented by the level set function. The pressure loss and the heat transfer performance of the channels are evaluated by the Building‐Cube Method code, which is a Cartesian‐mesh CFD solver. The proposed method resulted in an agreement with previous study in the single‐objective problems in its topology and achieved global exploration of non‐dominated solutions in the multi‐objective problems. © 2016 The Authors International Journal for Numerical Methods in Engineering Published by John Wiley & Sons Ltd     

Robust topology optimization under multiple independent uncertainties of loading positions

The topology optimization problem of a continuum structure is further investigated under the independent position uncertainties of multiple external loads, which are now described with an interval vector of uncertain-but-bounded variables. In this study, the structural compliance is formulated with the quadratic Taylor series expansion of multiple loading positions. As a result, the objective gradient information to the topological variables can be evaluated efficiently upon an explicit quadratic expression as the loads deviate from their ideal application points. Based on the minimum (largest absolute) value of design sensitivities, which corresponds to the most sensitive compliance to the load position variations, a two-level optimization algorithm within the non-probabilistic approach is developed upon a gradient-based optimization method. The proposed framework is then performed to achieve the robust optimal configurations of four benchmark examples, and the final designs are compared comprehensively with the traditional topology optimizations under the loading point fixation. It will be observed that the present methodology can provide a remarkably different structural layout with the auxiliary components in the design domain to counteract the load position uncertainties. The numerical results also show that the present robust topology optimization can effectively prevent the structural performance from a noticeable deterioration than the deterministic optimization in the presence of load position disturbances.     

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.     

A new algorithm for the robust optimization of rotor-bearing systems

This article presents a new algorithm for the robust optimization of rotor-bearing systems. The goal of the optimization problem is to find the values of a set of parameters for which the natural frequencies of the system are as far away as possible from the rotational speeds of the machine. To accomplish this, the penalization proposed by Ritto, Lopez, Sampaio, and Souza de Cursi in 2011 is employed. Since the rotor-bearing system is subject to uncertainties, such a penalization is modelled as a random variable. The robust optimization is performed by minimizing the expected value and variance of the penalization, resulting in a multi-objective optimization problem (MOP). The objective function of this MOP is known to be non-convex and it is shown that its resulting Pareto front (PF) is also non-convex. Thus, a new algorithm is proposed for solving the MOP: the normal boundary intersection (NBI) is employed to discretize the PF handling its non-convexity, and a global optimization algorithm based on a restart procedure and local searches are employed to minimize the NBI subproblems tackling the non-convexity of the objective function. A numerical analysis section shows the advantage of using the proposed algorithm over the weighted-sum (WS) and NSGA-II approaches. In comparison with the WS, the proposed approach obtains a much more even and useful set of Pareto points. Compared with the NSGA-II approach, the proposed algorithm provides a better approximation of the PF requiring much lower computational cost.