A velocity field level set method for shape and topology optimization

In this paper, we propose a new implementation of the level set shape and topology optimization, the velocity field level set method. Therein, the normal velocity field is constructed with specified basis functions and velocity design variables defined on a given set of points that are independent of the finite element mesh. A general mathematical programming algorithm can be employed to find the optimal normal velocities on the basis of the sensitivity analysis. As compared with conventional level set methods, mapping the variational boundary shape optimization problem into a finite‐dimensional design space and the use of a general optimizer makes it more efficient and straightforward to handle multiple constraints and additional design variables. Moreover, the level set function is updated by the Hamilton‐Jacobi equation using the normal velocity field; thus, the inherent merits of the implicit representation is retained. Therefore, this method combines the merits of both the general mathematical programming and conventional level set methods. Integrated topology optimization of structures with embedded components of designable geometries is considered to show the capability of this method to deal with general design variables. Several numerical examples in 2D or 3D design domains illustrate the robustness and efficiency of the method using different basis functions.     

基于遗传算法的离散型结构拓扑优化设计

采用遗传算法求解包括桁架结构和框架结构的离散型结构拓扑优化问题。在遗传算法的基础上,通过引入拓扑变量并修改被删除杆件的材料弹性模量,提出了一个受多工况荷载作用,能同时考虑应力、稳定及位移等约束的离散型结构拓扑优化问题统一数学模型。该模型不但能同时适用于桁架结构和框架结构等离散型结构拓扑优化问题,而且还能解决奇异最优解问题。结合上述统一数学模型和遗传算法,给出了求解离散型结构拓扑优化问题的优化方法。算例结果表明,采用该文提出的拓扑优化方法可有效、方便地对桁架结构、框架结构等离散型结构进行拓扑优化设计。     

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.     

Topology optimization of continuum structures with displacement constraints based on meshless method

In this paper, the element free Galerkin method (EFG) is applied to carry out the topology optimization of continuum structures with displacement constraints. In the EFG method, the matrices in the discretized system equations are assembled based on the quadrature points. In the sense, the relative density at Gauss quadrature point is employed as design variable. Considering the minimization of weight as an objective function, the mathematical formulation of the topology optimization subjected to displacement constraints is developed using the solid isotropic microstructures with penalization interpolation scheme. Moreover, the approximate explicit function expression between topological variables and displacement constraints are derived. Sensitivity of the objective function is derived based on the adjoint method. Three numerical examples are used to demonstrate the feasibility and effectiveness of the proposed method.     

A nodal variable method of structural topology optimization based on Shepard interpolant

A method for topology optimization of continuum structures based on nodal density variables and density field mapping technique is investigated. The original discrete‐valued topology optimization problem is stated as an optimization problem with continuous design variables by introducing a material density field into the design domain. With the use of the Shepard family of interpolants, this density field is mapped onto the design space defined by a finite number of nodal density variables. The employed interpolation scheme has an explicit form and satisfies range‐restricted properties that makes it applicable for physically meaningful density interpolation. Its ability to resolve more complex spatial distribution of the material density within an individual element, as compared with the conventional elementwise design variable approach, actually provides certain regularization to the topology optimization problem. Numerical examples demonstrate the validity and applicability of the proposed formulation and numerical techniques. Copyright © 2011 John Wiley & Sons, Ltd.     

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

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

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.     

Optimization of truss topology using tabu search

A design procedure for integrating topological considerations in the framework of structural optimization is presented. The proposed approach is capable of considering multiple load conditions, stress, displacement and local/global buckling constraints, and multiple objective functions in the problem formulation. Further, since the proposed method permits members to be added to or deleted from an existing topology and the topology is not defined by member areas, the difficulty of not being able to reach singular optima is also avoided. These objectives are accomplished using a discrete optimization procedure which uses 0–1 topological variables to optimize alternate designs. Since the topological variables are discrete in nature and the member cross-sections are assumed to be continuous, the topological optimization problem has mixed discrete-continuous variables. This non-linear programming problem is solved using a memory-based combinatorial optimization technique known as tabu search. Numerical results obtained using tabu search for single and multiobjective topological optimization of truss structures are presented. To model the multiple objective functions in the problem formulation, a cooperative game theoretic approach is used. The results indicate that the optimum topologies obtained using tabu search compare favourably, and in some instances, outperform the results obtained using the ground–structure approach. However, this improvement occurs at the expense of a significant increase in computational burden owing to the fact that the proposed approach necessitates that the geometry of each trial topology be optimized.     

Structural shape and topology optimization based on level-set modelling and the element-propagating method

This article introduces the element-propagating method to structural shape and topology optimization. Structural optimization based on the conventional level-set method needs to solve several partial differential equations. By the insertion and deletion of basic material elements around the geometric boundary, the element-propagating method can avoid solving the partial differential equations and realize the dynamic updating of the material region. This approach also places no restrictions on the signed distance function and the Courant–Friedrichs–Lewy condition for numerical stability. At the same time, in order to suppress the dependence on the design initialization for the 2D structural optimization problem, the strain energy density is taken as a criterion to generate new holes in the material region. The coupled algorithm of the element-propagating method and the method for generating new holes makes the structural optimization more robust. Numerical examples demonstrate that the proposed approach greatly improves numerical efficiency, compared with the conventional level-set method for structural topology optimization.     

Continuum structural topological optimizations with stress constraints based on an active constraint technique

Stress‐related problems have not been given the same attention as the minimum compliance topological optimization problem in the literature. Continuum structural topological optimization with stress constraints is of wide engineering application prospect, in which there still are many problems to solve, such as the stress concentration, an equivalent approximate optimization model and etc. A new and effective topological optimization method of continuum structures with the stress constraints and the objective function being the structural volume has been presented in this paper. To solve the stress concentration issue, an approximate stress gradient evaluation for any element is introduced, and a total aggregation normalized stress gradient constraint is constructed for the optimized structure under the r?th load case. To obtain stable convergent series solutions and enhance the control on the stress level, two p‐norm global stress constraint functions with different indexes are adopted, and some weighting p‐norm global stress constraint functions are introduced for any load case. And an equivalent topological optimization model with reduced stress constraints is constructed,being incorporated with the rational approximation for material properties, an active constraint technique, a trust region scheme, and an effective local stress approach like the qp approach to resolve the stress singularity phenomenon. Hence, a set of stress quadratic explicit approximations are constructed, based on stress sensitivities and the method of moving asymptotes. A set of algorithm for the one level optimization problem with artificial variables and many possible non‐active design variables is proposed by adopting an inequality constrained nonlinear programming method with simple trust regions, based on the primal‐dual theory, in which the non‐smooth expressions of the design variable solutions are reformulated as smoothing functions of the Lagrange multipliers by using a novel smoothing function. Finally, a two‐level optimization design scheme with active constraint technique, i.e. varied constraint limits, is proposed to deal with the aggregation constraints that always are of loose constraint (non active constraint) features in the conventional structural optimization method. A novel structural topological optimization method with stress constraints and its algorithm are formed, and examples are provided to demonstrate that the proposed method is feasible and very effective. © 2016 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons Ltd.     

基于物质点描述的双向渐进式拓扑优化方法

:为了克服连续体结构拓扑优化中的数值不稳定现象,定义了表征物质点及其领域有无的物质点拓扑变量,提出基于物质点描述的双向渐进式拓扑优化方法.基于过滤法构造拓扑变量场的插值函数,从而在拓扑优化模型中自然消除了棋盘格现象.为适用于不同单元类型和网格离散形式等,重新定义了灵敏度密度.通过二维数值算例对理论方法进行验证.结果表明:方法在连续体结构拓扑优化设计中具有可行性和有效性.     

OPTIMUM SHAPE DESIGN AND POSITIONING OF FEATURES USING THE BOUNDARY INTEGRAL EQUATION METHOD

A design optimization procedure is developed using the boundary integral equation (BIE) method for linear elastostatic two-dimensional domains. Optimal shape design problems are treated where design variables are geometric parameters such as the positions and sizing dimensions of entire features on a component or structure. A fully analytical approach is adopted for the design sensitivity analysis where the BIE is implicitly differentiated. The ability to evaluate response sensitivity derivatives with respect to design variables such as feature positions is achieved through the definition of appropriate design velocity fields for these variables. How the advantages of the BIE method are amplified when extended to sensitivity analysis for this category of shape design problems is also highlighted. A mathematical programming approach with the penalty function method is used for solving the overall optimization problem. The procedure is applied to three example problems to demonstrate the optimum positioning of holes and optimization of radial dimensions of circular arcs on structures.     

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.     

A hybrid cellular automaton–bi-directional evolutionary optimization algorithm for topological optimization of crashworthiness

This article proposes a new algorithm for topological optimization under dynamic loading which combines cellular automata with bi-directional evolutionary structural optimization (BESO). The local rules of cellular automata are used to update the design variables, which avoids the difficulty of obtaining gradient information under nonlinear collision conditions. The intermediate-density design problem of hybrid cellular automata is solved using the BESO concept of 0–1 binary discrete variables. Some improvement strategies are also proposed for the hybrid algorithm to solve certain problems in nonlinear topological optimization, e.g. numerical oscillation. Some typical examples of crashworthiness problems are provided to illustrate the efficiency of the proposed method and its ability to find the final optimal solution. Finally, numerical results obtained using the proposed algorithms are compared with reference examples taken from the literature. The results show that the hybrid method is computationally efficient and stable.     

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.     

约束阻尼结构的改进准则法拓扑减振动力学优化

旨在为减振设计提供理论基础,研究约束阻尼结构拓扑动力学优化。以阻尼材料用量、振动特征方程、模态频率为约束,以多模态损耗因子倒数的加权和最小为目标,建立了约束阻尼结构拓扑优化模型,引入MAC因子控制结构的振型跃阶。在引入质量阵惩罚因子基础上推导出优化目标灵敏度。考虑到优化目标函数的非凸性,采用常规准则法(OC)寻优可能会使拓扑变量出现负值或陷入局部优化,故引入数学规划移动渐近技术对OC法进行改进,从而将全体拓扑变量纳入改进算法的优化迭代全过程。编程实现了约束阻尼板改进OC法拓扑动力学优化并对改进法性能进行了仿真。结果显示,改进算法可得到更合理的约束阻尼层构形,可使结构取得更佳减振效果。研究表明,改进算法迭代稳定性更好、寻优效率更高、更具全域最优性。     

On a marching level-set method for extended discontinuous Galerkin methods for incompressible two-phase flows: Application to two-dimensional settings

In this work a solver for two-dimensional, instationary two-phase flows on the basis of the extended discontinuous Galerkin (extended DG/XDG) method is presented. The XDG method adapts the approximation space conformal to the position of the interface. This allows a subcell accurate representation of the incompressible Navier-Stokes equations in their sharp interface formulation. The interface is described as the zero set of a signed-distance level-set function and discretized by a standard DG method. For the interface, resp. level-set, evolution an extension velocity field is used and a two-staged algorithm is presented for its construction on a narrow-band. On the cut-cells a monolithic elliptic extension velocity method is adapted and a fast-marching procedure on the neighboring cells. The spatial discretization is based on a symmetric interior penalty method and for the temporal discretization a moving interface approach is adapted. A cell agglomeration technique is utilized for handling small cut-cells and topology changes during the interface motion. The method is validated against a wide range of typical two-phase surface tension driven flow phenomena in a 2D setting including capillary waves, an oscillating droplet and the rising bubble benchmark.     

基于Matlab的两级星型齿轮传动的优化设计

建立了两级内外啮合星型齿轮传动优化设计的数学模型,分析了系统的传动方案和约束条件,以中心距最小为目标函数,确定了主要优化设计参数。采用了以齿数、模数为离散变量和变位系数,齿数比等为连续变量的优化设计方法。过多的等式约束造成最优化时可行域减少,甚至找不到可行域,严重影响了优化设计的运行结果。鉴于此,发展了分级优化的方法,以齿数、模数为离散变量,变位系数、齿数比等为连续变量,运用Matlab中非线性有约束的多元函数fmincon求解,经过优化计算,得到主要优化设计参数。