Radial basis function based level set interpolation and evolution for deformable modelling

We present a study in level set representation and evolution using radial basis functions (RBFs) for active contour and active surface models. It builds on recent works by others who introduced RBFs into level sets for structural topology optimisation. Here, we introduce the concept into deformable models and present a new level set formulation able to handle more complex topological changes, in particular perturbation away from the evolving front. In the conventional level set technique, the initial active contour/surface is implicitly represented by a signed distance function and periodically re-initialised to maintain numerical stability. We interpolate the initial distance function using RBFs on a much coarser grid, which provides great potential in modelling in high dimensional space. Its deformation is considered as an updating of the RBF interpolants, an ordinary differential equation (ODE) problem, instead of a partial differential equation (PDE) problem, and hence it becomes much easier to solve. Re-initialisation is found no longer necessary, in contrast to conventional finite difference method (FDM) based level set approaches. The proposed level set updating scheme is efficient and does not suffer from self-flattening while evolving, hence it avoids large numerical errors. Further, more complex topological changes are readily achievable and the initial contour or surface can be placed arbitrarily in the image. These properties are extensively demonstrated on both synthetic and real 2D and 3D data. We also present a novel active contour model, implemented with this level set scheme, based on multiscale learning and fusion of image primitives from vector-valued data, e.g. colour images, without channel separation or decomposition.     

Design of distributed compliant micromechanisms with an implicit free boundary representation

In this paper, a parameterization approach is presented for structural shape and topology optimization of compliant mechanisms using a moving boundary representation. A level set model is developed to implicitly describe the structural boundary by embedding into a scalar function of higher dimension as zero level set. The compactly supported radial basis function of favorable smoothness and accuracy is used to interpolate the level set function. Thus, the temporal and spatial initial value problem is now converted into a time-separable parameterization problem. Accordingly, the more difficult shape and topology optimization of the Hamilton–Jacobi equation is then transferred into a relatively easy size optimization with the expansion coefficients as design variables. The design boundary is therefore advanced by applying the optimality criteria method to iteratively evaluate the size optimization so as to update the level set function in accordance with expansion coefficients of the interpolation. The optimization problem of the compliant mechanism is established by including both the mechanical efficiency as the objective function and the prescribed material usage as the constraint. The design sensitivity analysis is performed by utilizing the shape derivative. It is noted that the present method is not only capable of simultaneously addressing shape fidelity and topology changes with a smooth structural boundary but also able to avoid some of the unfavorable numerical issues such as the Courant–Friedrich–Levy condition, the velocity extension algorithm, and the reinitialization procedure in the conventional level set method. In particular, the present method can generate new holes inside the material domain, which makes the final design less insensitive to the initial guess. The compliant inverter is applied to demonstrate the availability of the present method in the framework of the implicit free boundary representation.     

A level-set-based topology and shape optimization method for continuum structure under geometric constraints

Recent advances in level-set-based shape and topology optimization rely on free-form implicit representations to support boundary deformations and topological changes. In practice, a continuum structure is usually designed to meet parametric shape optimization, which is formulated directly in terms of meaningful geometric design variables, but usually does not support free-form boundary and topological changes. In order to solve the disadvantage of traditional step-type structural optimization, a unified optimization method which can fulfill the structural topology, shape, and sizing optimization at the same time is presented. The unified structural optimization model is described by a parameterized level set function that applies compactly supported radial basis functions (CS-RBFs) with favorable smoothness and accuracy for interpolation. The expansion coefficients of the interpolation function are treated as the design variables, which reflect the structural performance impacts of the topology, shape, and geometric constraints. Accordingly, the original topological shape optimization problem under geometric constraint is fully transformed into a simple parameter optimization problem; in other words, the optimization contains the expansion coefficients of the interpolation function in terms of limited design variables. This parameterization transforms the difficult shape and topology optimization problems with geometric constraints into a relatively straightforward parameterized problem to which many gradient-based optimization techniques can be applied. More specifically, the extended finite element method (XFEM) is adopted to improve the accuracy of boundary resolution. At last, combined with the optimality criteria method, several numerical examples are presented to demonstrate the applicability and potential of the presented method.     

连续体结构的参数化水平集统一优化

针对传统分步式结构优化设计的不足,提出一种同时进行结构拓扑、形状和尺寸统一优化的设计方法.首先采用水平集函数描述统一的结构优化模型和几何尺寸边界,通过引入紧支径向插值基函数将结构拓扑优化变量、形状优化变量和尺寸优化变量变换为基函数的扩展系数;然后取该扩展系数为设计变量,借助一种参数的变化表达3种优化要素对结构性能的影响,将复杂的多变量优化问题变换为相对简单的参数优化问题,有利于与相对成熟的优化算法相结合提高求解效率;进一步用R函数将其融合为一个整体,构造出统一优化模型,并用最优化准则法进行求解.最后通过数值案例证明了该方法的有效性和精确性.     

A level set method for structural shape and topology optimization using radial basis functions

This paper presents an alternative level set method for shape and topology optimization of continuum structures. An implicit free boundary representation model is established by embedding structural boundary into the zero level set of a higher-dimensional level set function. An explicit parameterization scheme for the level set surface is proposed by using radial basis functions with compact support. In doing so, the originally more difficult shape and topology optimization, driven by the temporal and spatial Hamilton–Jacobi partial differential equation (PDE), is transformed into a relatively easier size optimization of the expansion coefficients of the basis functions. The design optimization is converted to an iterative numerical process that combines the parameterization with a derivation of the shape sensitivity of the design functions, so as to allow using mathematical programming algorithms to solve the level set-based design problem and avoid directly solving the Hamilton–Jacobi PDE. Furthermore, a numerically more stable and efficient volume integration scheme is proposed to implement calculations of the shape derivatives, leading to the creation of new holes which are generated initially along the boundary and then propagated to the interior of the design domain. Two widely studied examples are used to demonstrate the effectiveness of the proposed optimization method.     

Simultaneous shape and topology optimization of shell structures

In this research, Method of Moving Asymptotes (MMA) is utilized for simultaneous shape and topology optimization of shell structures. It is shown that this approach is well matched with the large number of topology and shape design variables. The currently practiced technology for optimization is to find the topology first and then to refine the shape of structure. In this paper, the design parameters of shape and topology are optimized simultaneously in one go. In order to model and control the shape of free form shells, the NURBS (Non Uniform Rational B-Spline) technology is used. The optimization problem is considered as the minimization of mean compliance with the total material volume as active constraint and taking the shape and topology parameters as design variables. The material model employed for topology optimization is assumed to be the Solid Isotropic Material with Penalization (SIMP). Since the MMA optimization method requires derivatives of the objective function and the volume constraint with respect to the design variables, a sensitivity analysis is performed. Also, for alleviation of the instabilities such as mesh dependency and checkerboarding the convolution noise cleaning technique is employed. Finally, few examples taken from literature are presented to demonstrate the performance of the method and to study the effect of the proposed concurrent approach on the optimal design in comparison to the sequential topology and shape optimization methods.     

Semi-Lagrange method for level-set-based structural topology and shape optimization

In this paper, we introduce a semi-Lagrange scheme to solve the level-set equation in structural topology optimization. The level-set formulation of the problem expresses the optimization process as a solution to a Hamilton–Jacobi partial differential equation. It allows for the use of shape sensitivity to derive a speed function for a descent solution. However, numerical stability condition in the explicit upwind scheme for discrete level-set equation severely restricts the time step, requiring a large number of time steps for a numerical solution. To improve the numerical efficiency, we propose to employ a semi-Lagrange scheme to solve level-set equation. Therefore, a much larger time step can be obtained and a much smaller number of time steps are required. Numerical experiments comparing the semi-Lagrange method with the classical explicit upwind scheme are presented for the problem of mean compliance optimization in two dimensions.     

A semi-Lagrangian level set method for structural optimization

In level set based structural optimization, semi-Lagrange method has an advantage to allow for a large time step without the limitation of Courant–Friedrichs–Lewy (CFL) condition for numerical stability. In this paper, a line search algorithm and a sensitivity modulation scheme are introduced for the semi-Lagrange method. The line search attempts to adaptively determine an appropriate time step in each iteration of optimization. With consideration of some practical characteristics of the topology optimization process, incorporating the line search into semi-Lagrange optimization method can yield fewer design iterations and thus improve the overall computational efficiency. The sensitivity modulation is inspired from the conjugate gradient method in finite-dimensions, and provides an alternative to the standard steepest descent search in level set based optimization. Two benchmark examples are presented to compare the sensitivity modulation and the steepest descent techniques with and without the line search respectively.     

Piecewise constant level set method for structural topology optimization with MBO type of projection

In this paper, we combine a Piecewise Constant Level Set (PCLS) method with a MBO scheme to solve a structural shape and topology optimization problem. The geometrical boundary of structure is represented implicitly by the discontinuities of PCLS functions. Compared with the classical level set method (LSM) for solving Hamilton–Jacobi partial differential equation (H-J PDE) we don’t need to solve H-J PDE, thus it is free of the CFL condition and the reinitialization scheme. For solving optimization problem under some constraints, Additive Operator Splitting (AOS) and Multiplicative Operator Splitting (MOS) schemes will be used. To increase the convergency speed and the efficiency of PCLS method we combine this approach with MBO scheme. Advantages and disadvantages are discussed by solving some examples of 2D structural topology optimization problems.     

基于反应扩散方程的水平集结构拓扑优化方法与实现

为实现更加先进的拓扑优化算法,研究采用反应扩散方程的水平集结构拓扑优化方法,通过理论推导给出算法中的参数选择建议.该方法允许在拓扑优化过程中生成新的孔洞,初始结构无须包含孔洞,不需要重新初始化步骤,从而可提高算法的收敛性.针对传统拓扑优化中主要采用体积约束、以柔度最小为目标和体积保留率设定存在一定主观性的问题,探究不同体积保留率下的结构应力水平的变化规律,结果显示可以依据结构最大应力水平与体积保留率的变化规律确定最优体积保留率.     

A topology optimization method based on the level set method incorporating a fictitious interface energy

This paper proposes a new topology optimization method, which can adjust the geometrical complexity of optimal configurations, using the level set method and incorporating a fictitious interface energy derived from the phase field method. First, a topology optimization problem is formulated based on the level set method, and the method of regularizing the optimization problem by introducing fictitious interface energy is explained. Next, the reaction–diffusion equation that updates the level set function is derived and an optimization algorithm is then constructed, which uses the finite element method to solve the equilibrium equations and the reaction–diffusion equation when updating the level set function. Finally, several optimum design examples are shown to confirm the validity and utility of the proposed topology optimization method.     

Topology optimization of steady Navier-Stokes flow via a piecewise constant level set method

This paper presents a piecewise constant level set method for the topology optimization of steady Navier-Stokes flow. Combining piecewise constant level set functions and artificial friction force, the optimization problem is formulated and analyzed based on a design variable. The topology sensitivities are computed by the adjoint method based on Lagrangian multipliers. In the optimization procedure, the piecewise constant level set function is updated by a new descent method, without the needing to solve the Hamilton-Jacobi equation. To achieve optimization, the piecewise constant level set method does not track the boundaries between the different materials but instead through the regional division, which can easily create small holes without topological derivatives. Furthermore, we make some attempts to avoid updating the Lagrangian multipliers and to deal with the constraints easily. The algorithm is very simple to implement, and it is possible to obtain the optimal solution by iterating a few steps. Several numerical examples for both two- and three-dimensional problems are provided, to demonstrate the validity and efficiency of the proposed method.     

Topology optimization of shell structures using adaptive inner-front (AIF) level set method

A new topology optimization using adaptive inner-front level set method is presented. In the conventional level set-based topology optimization, the optimum topology strongly depends on the initial level set due to the incapability of inner-front creation during the optimization process. In the present work, in this regard, an algorithm for inner-front creation is proposed in which the sizes, the positions, and the number of new inner-fronts during the optimization process can be globally and consistently identified. In the algorithm, the criterion of inner-front creation for compliance minimization problems of linear elastic structures is chosen as the strain energy density along with volumetric constraint. To facilitate the inner-front creation process, the inner-front creation map is constructed and used to define new level set function. In the implementation of inner-front creation algorithm, to suppress the numerical oscillation of solutions due to the sharp edges in the level set function, domain regularization is carried out by solving the edge smoothing partial differential equation (smoothing PDE). To update the level set function during the optimization process, the least-squares finite element method (LSFEM) is adopted. Through the LSFEM, a symmetric positive definite system matrix is constructed, and non-diffused and non-oscillatory solution for the hyperbolic PDE such as level set equation can be obtained. As applications, three-dimensional topology optimization of shell structures is treated. From the numerical examples, it is shown that the present method brings in much needed flexibility in topologies during the level set-based topology optimization process.     

A new method based on adaptive volume constraint and stress penalty for stress-constrained topology optimization

This paper focuses on the stress-constrained topology optimization of minimizing the structural volume and compliance. A new method based on adaptive volume constraint and stress penalty is proposed. According to this method, the stress-constrained volume and compliance minimization topology optimization problem is transformed into two simple and related problems: a stress-penalty-based compliance minimization problem and a volume-decision problem. In the former problem, stress penalty is conducted and used to control the local stress level of the structure. To solve this problem, the parametric level set method with the compactly supported radial basis functions is adopted. Meanwhile, an adaptive adjusting scheme of the stress penalty factor is used to improve the control of the local stress level. To solve the volume-decision problem, a combination scheme of the interval search and local search is proposed. Numerical examples are used to test the proposed method. Results show the lightweight design, which meets the stress constraint and whose compliance is simultaneously optimized, can be obtained by the proposed method.     

A study on X-FEM in continuum structural optimization using a level set model

In this paper, we implement the extended finite element method (X-FEM) combined with the level set method to solve structural shape and topology optimization problems. Numerical comparisons with the conventional finite element method in a fixed grid show that the X-FEM leads to more accurate results without increasing the mesh density and the degrees of freedom. Furthermore, the mesh in X-FEM is independent of the physical boundary of the design, so there is no need for remeshing during the optimization process. Numerical examples of mean compliance minimization in 2D are studied in regard to efficiency, convergence and accuracy. The results suggest that combining the X-FEM for structural analysis with the level set based boundary representation is a promising approach for continuum structural optimization.     

Stress isolation through topology optimization

This paper introduces a problem of stress isolation in structural design and presents an approach to the problem through topology optimization. We model the stress isolation problem as a topology optimization problem with multiple stress constraints in different regions. The shape equilibrium constraint approach is employed to effectively control the local stress constraints. The level set based structural optimization is implemented with the extended finite element method (X-FEM) for providing an adequately accurate stress analysis. Numerical examples of stress isolation design in two dimensions are investigated as a benchmark test of the proposed method. The results, from the force transmittance point of view, suggest that the guard “grooves” obtained can change the force path to successfully realize the stress isolation in the structure.     

Adaptive volume constraint algorithm for stress limit-based topology optimization

Homogenization or density-based topology optimization methods work by distributing a fixed amount of material to the most effective areas of the design domain so as to create an optimal structural configuration that meets the minimum compliance criteria. These topology optimization methods generally cannot control the maximum stress levels of the structure; therefore, the smoothened optimum structure is not guaranteed to be ready for immediate use. This can be because it is either unsafe if the maximum stress at this structure exceeds the strength limit, or over designed if the maximum stress is far below the stress limit. Difficult and complex shape optimization must then be done to obtain a minimum-weight structure that meets the maximum stress constraint. This paper proposes an adaptive volume constraint (AVC) algorithm, a heuristic approach, in place of traditional topology optimization methods so that the maximum stress in the optimal structural configuration will be below the predefined stress limit and the smoothened structure will be directly applicable. In order to test the applicability and robustness of the AVC algorithm, topology optimization using both a traditional fixed volume constraint and an AVC are tested on a number of configuration design problems. To further illustrate the usefulness of the AVC algorithm, shape optimizations at the maximum stress constraint are also conducted on the smooth structural models by both optimization approaches on an identical problem set.     

Simultaneous optimization of the material properties and the topology of functionally graded structures

A level set based method is proposed for the simultaneous optimization of the material properties and the topology of functionally graded structures. The objective of the present study is to determine the optimal material properties (via the material volume fractions) and the structural topology to maximize the performance of the structure in a given application. In the proposed method, the volume fraction and the structural boundary are considered as the design variables, with the former being discretized as a scalar field and the latter being implicitly represented by the level set method. To perform simultaneous optimization, the two design variables are integrated into a common objective functional. Sensitivity analysis is conducted to obtain the descent directions. The optimization process is then expressed as the solution to a coupled Hamilton-Jacobi equation and diffusion partial differential equation. Numerical results are provided for the problem of mean compliance optimization in two dimensions.     

Stress-related topology optimization via level set approach

Considering stress-related objective or constraint functions in structural topology optimization problems is very important from both theoretical and application perspectives. It has been known, however, that stress-related topology optimization problem is challenging since several difficulties must be overcome in order to solve it effectively. Traditionally, SIMP (Solid Isotropic Material with Penalization) method was often employed to tackle it. Although some remarkable achievements have been made with this computational framework, there are still some issues requiring further explorations. In the present work, stress-related topology optimization problems are investigated via a level set-based approach, which is a different topology optimization framework from SIMP. Numerical examples show that under appropriate problem formulations, level set approach is a promising tool for stress-related topology optimization problems.