Topology design of large displacement compliant mechanisms with multiple materials and multiple output ports

Topology optimization of compliant mechanisms is presented in this paper wherein the layout design problem is addressed in its original binary or discrete (0-1) form. Design variables are modeled as discrete variables and allowed to assume values pertaining only to their void (0) or solid (1) states. Due to this discrete nature, a genetic algorithm is employed as an optimization routine. Using the barrier assignment approach, the search algorithm is extended to use with multiple materials. The layout design of compliant mechanisms is performed wherein displacements at multiple points (ports) in the design region are maximized along the respective prescribed directions. With multiple output ports and multiple materials, additional freedom in motion and force transduction can be achieved with compliant mechanisms. Geometrically large deformation analysis is employed to compute the displacement-based multiple objectives that are extremized using Nondominated Sorting in Genetic Algorithms (or NSGA). With genetic algorithms, buckling or snap through like issues with nonconvergent solutions in the population when computing nonlinear deformations can be implicitly circumvented.     

Application of spectral level set methodology in topology optimization

In this paper, two benchmark problems in structural boundary design are solved using the spectral level set methodology, which is a new approach to topology optimization of interfaces. This methodology is an extension of the level set methods, in which the interface is represented as the zero level set of a function. According to the proposed formulation, the Fourier coefficients of that function are the design variables describing the interface during the topology optimization. An advantage of the spectral level set methodology, in the case of a sufficiently regular interface, is to admit an upper bound error which is asymptotically smaller than the one for nonadaptive spacial discretizations of the level set function. Other advantages include the nucleation of holes in the interior of the interface and the avoidance of checkerboard-like designs. The theoretical framework of the methodology is presented and estimates on its convergence rate are discussed. The numerical applications consist in the design of short and long cantilevers subject to a vertical concentrated load opposite to the fixed end. The goal is to maximize the structural stiffness subject to a solid volume constraint. The zero level set of the level set function defines the structural boundary.     

Topology optimization of hinge-free compliant mechanisms with multiple outputs using level set method

A method for topology optimization of hinge-free compliant mechanisms with multiple outputs using level set method is presented in this paper. The focus of this paper is on how to prevent generating the flexible hinges during the process of topology optimization of compliant mechanisms. In the proposed method, two types of mean compliances are introduced and built in the proposed multi-objective function for topology optimization of hinge-free compliant mechanisms with multiple outputs, therefore, the spring model widely used for topology optimization of compliant mechanisms is no longer needed. Some numerical examples are presented to illustrate the validity 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.     

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 new level set method for inhomogeneous image segmentation

Intensity inhomogeneity often appears in medical images, such as X-ray tomography and magnetic resonance (MR) images, due to technical limitations or artifacts introduced by the object being imaged. It is difficult to segment such images by traditional level set based segmentation models. In this paper, we propose a new level set method integrating local and global intensity information adaptively to segment inhomogeneous images. The local image information is associated with the intensity difference between the average of local intensity distribution and the original image, which can significantly increase the contrast between foreground and background. Thus, the images with intensity inhomogeneity can be efficiently segmented. What is more, to avoid the re-initialization of the level set function and shorten the computational time, a simple and fast level set evolution formulation is used in the numerical implementation. Experimental results on synthetic images as well as real medical images are shown in the paper to demonstrate the efficiency and robustness of 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.     

Shape design of pin-jointed multistable compliant mechanisms using snapthrough behavior

A general approach is presented for generating pin-jointed multistable compliant mechanisms using snapthrough behavior. An optimization problem is formulated for minimizing the total structural volume under constraints on the displacements at the specified nodes, stiffnesses at initial and final states, and load factors to lead to snapthrough behavior. The design variables are cross-sectional areas and the nodal coordinates. It is shown in the numerical examples that several mechanisms can be naturally found as a result of optimization starting from randomly selected initial solutions. It is also shown that no local bifurcation point exists along the equilibrium path, and the obtained mechanism is not sensitive to initial imperfections.     

A level set method for oil slick segmentation in SAR images

Despite much effort and significant progress in recent years, image segmentation remains a challenging problem in image processing, especially for the low contrast, noisy synthetic aperture radar (SAR) images. This paper explores the segmentation of oil slicks using a partial differential equation (PDE)‐based level set method, which represents the slick surface as an implicit propagation interface. Starting from an initial estimation with priori information, the level set method creates a set of speed functions to detect the position of the propagation interface. Specifically, the image intensity gradient and the curvature flow are utilized together to determine the speed and direction of the propagation. This allows the front interface to propagate naturally with topological changes, significant protrusions and narrow regions, giving rise to stable and smooth boundaries that discriminate oil slicks from the surrounding water. As the speckles are removed concurrently while the front interface propagates, the pre‐filtering of noise is saved. The proposed method has been illustrated by experiments on oil slick segmentation using the ERS‐2 SAR images. Its advantages over the traditional image segmentation approaches have also been demonstrated.     

Piecewise constant level set method for topology optimization of unilateral contact problems

The paper deals with the structural optimization of the elastic body in unilateral contact with a rigid foundation using the level set approach. A piecewise constant level set method is used to represent the evolution of interfaces rather than the standard method. The piecewise constant level set function takes distinct constant values in each subdomain of a whole design domain. Using a two-phase approximation the original structural optimization problem is reformulated as an equivalent constrained optimization problem in terms of the piecewise constant level set function. Necessary optimality condition is formulated. Finite difference and finite element methods are applied as the approximation methods. Numerical examples are provided and discussed.     

Topology optimization of structures with integral compliant mechanisms for mid-frequency response

The purpose of this paper is to present a method for developing new truss-like sandwich structures that exhibit desirable mid-frequency vibratory characteristics. Specifically, a genetic algorithm optimization routine is used to determine candidate small scale structural topologies, i.e. unit cells, that may be used in the design of larger scale periodic sandwich structures. This multi-scale procedure is demonstrated starting with several unit cell topology optimization examples. From these examples a specific optimal unit cell is selected for further investigation and integration into a periodic sandwich beam. Computational results indicate that the proposed optimization approach is effective when used to design new structures for reduced mid-frequency vibratory response.     

Topological shape optimization of geometrically nonlinear structures using level set method

Using the level set method, a topological shape optimization method is developed for geometrically nonlinear structures in total Lagrangian formulation. The structural boundaries are implicitly represented by the level set function, obtainable from “Hamilton-Jacobi type” equation with “up-wind scheme,” embedded into a fixed initial domain. The method minimizes the compliance through the variations of implicit boundary, satisfying an allowable volume requirement. The required velocity field to solve the Hamilton-Jacobi equation is determined by the descent direction of Lagrangian derived from an optimality condition. Since the homogeneous material property and implicit boundary are utilized, the convergence difficulty is significantly relieved.     

Accurate adaptive level set method and sharpening technique for three dimensional deforming interfaces

In this paper, we demonstrate improved accuracy of the level set method for resolving deforming interfaces by proposing two key elements: (1) accurate level set solutions on adapted Cartesian grids by judiciously choosing interpolation polynomials in regions of different grid levels and (2) enhanced re-initialization by an interface sharpening procedure. The level set equation is solved using a fifth order WENO scheme or a second order central differencing scheme depending on availability of uniform stencils at each grid point. Grid adaptation criteria are determined so that the Hamiltonian functions at nodes adjacent to interfaces are always calculated by the fifth order WENO scheme. This selective usage between the fifth order WENO and second order central differencing schemes is confirmed to give more accurate results compared to those in literature for standard test problems. In order to further improve accuracy especially near thin filaments, we suggest an artificial sharpening method, which is in a similar form with the conventional re-initialization method but utilizes sign of curvature instead of sign of the level set function. Consequently, volume loss due to numerical dissipation on thin filaments is remarkably reduced for the test problems.     

A level-set based variational method for design and optimization of heterogeneous objects

A heterogeneous object is referred to as a solid object made of different constituent materials. The object is of a finite collection of regions of a set of prescribed material classes of continuously varying material properties. These properties have a discontinuous change across the interface of the material regions. In this paper, we propose a level-set based variational approach for the design of this class of heterogeneous objects. Central to the approach is a variational framework for a well-posed formulation of the design problem. In particular, we adapt the Mumford-Shah model which specifies that any point of the object belongs to either of two types: inside a material region of a well-defined gradient or on the boundary edges and surfaces of discontinuities. Furthermore, the set of discontinuities is represented implicitly, using a multi-phase level set model. This level-set based variational approach yields a computational system of coupled geometric evolution and diffusion partial differential equations. Promising features of the proposed method include strong regularity in the problem formulation and inherent capabilities of geometric and material modeling, yielding a common framework for optimization of the heterogeneous objects that incorporates dimension, shape, topology, and material properties. The proposed method is illustrated with several 2D examples of optimal design of multi-material structures and materials.     

多目标轮廓Mumford-Shah水平集提取

目标轮廓的快速检测进而提取其几何形状,在图形图像处理中有着重要的作用.提出了一种多目标轮廓的水平集提取方法,对基于Mumford-Shah模型的C-V方法从两方面进行了改进:增加梯度矢量场和曲线法方向的融合作为边界吸引场,生成可以驱动主动轮廓向边缘进化的双向几何变形流,保留原图像分布信息作为区域进化能,解决未考虑局部几何信息造成的区域能量捕捉信息不全,或边缘梯度场和演化曲线法线方向正交时无法实现拓扑结构变化的缺陷;对水平集函数进行修正,使得它在收敛过程中能自动进行调整,确保其满足符号距离函数的要求,扩大初始化前迭代搜索区域,减少初始化次数,提高收敛效率;最后给出所提方法的数字化求解方案.实验表明该方法可行且具有较好的鲁棒性.     

Active contours with selective local or global segmentation: A new formulation and level set method

A novel region-based active contour model (ACM) is proposed in this paper. It is implemented with a special processing named Selective Binary and Gaussian Filtering RegularizedLevel Set(SBGFRLS) method, which first selectively penalizes the level set function to be binary, and then uses a Gaussian smoothing kernel to regularize it. The advantages of our method are as follows. First, a new region-based signed pressure force (SPF) function is proposed, which can efficiently stop the contours at weak or blurred edges. Second, the exterior and interior boundaries can be automatically detected with the initial contour being anywhere in the image. Third, the proposed ACM with SBGFRLS has the property of selective local or global segmentation. It can segment not only the desired object but also the other objects. Fourth, the level set function can be easily initialized with a binary function, which is more efficient to construct than the widely used signed distance function (SDF). The computational cost for traditional re-initialization can also be reduced. Finally, the proposed algorithm can be efficiently implemented by the simple finite difference scheme. Experiments on synthetic and real images demonstrate the advantages of the proposed method over geodesic active contours (GAC) and Chan–Vese (C–V) active contours in terms of both efficiency and accuracy.     

Compliant mechanism design using multi-objective topology optimization scheme of continuum structures

Topology optimization problems for compliant mechanisms using a density interpolation scheme, the rational approximation of material properties (RAMP) method, and a globally convergent version of the method of moving asymptotes (GCMMA) are primarily discussed. First, a new multi-objective formulation is proposed for topology optimization of compliant mechanisms, in which the maximization of mutual energy (flexibility) and the minimization of mean compliance (stiffness) are considered simultaneously. The formulation of one-node connected hinges, as well as checkerboards and mesh-dependency, is typically encountered in the design of compliant mechanisms. A new hybrid-filtering scheme is proposed to solve numerical instabilities, which can not only eliminate checkerboards and mesh-dependency efficiently, but also prevent one-node connected hinges from occurring in the resulting mechanisms to some extent. Several numerical applications are performed to demonstrate the validity of the methods presented in this paper.     

Structure topology optimization: fully coupled level set method via FEMLAB

This paper presents a procedure which can easily implement the 2D compliance minimization structure topology optimization by the level set method using the FEMLAB package. Instead of a finite difference solver for the level set equation, as is usually the case, a finite element solver for the reaction–diffusion equation is used to evolve the material boundaries. All of the optimization procedures are implemented in a user-friendly manner. A FEMLAB code can be downloaded from the homepage www.imtek.de/simulation and is free for educational purposes.     

A level set method based on the Bayesian risk for medical image segmentation

This paper proposes an alternative criterion derived from the Bayesian risk classification error for image segmentation. The proposed model introduces a region-based force determined through the difference of the posterior image densities for the different classes, a term based on the prior probability derived from Kullback-Leibler information number, and a regularity term adopted to avoid the generation of excessively irregular and small segmented regions. Compared with other level set methods, the proposed approach relies on the optimum decision of pixel classification and the estimates of prior probabilities; thus the approach has more reliability in theory and practice. Experiments show that the proposed approach is able to extract the complicated shapes of targets and robust for various types of medical images. Moreover, the algorithm can be easily extendable for multiphase segmentation.