Topology optimization of free vibrating continuum structures based on the element free Galerkin method

In this paper, the topology optimization design of the free vibrating continuum structures is formulated based on the element free Galerkin (EFG) method. Considering the relative density of nodes as design variable, and the maximization of the fundamental eigenvalue as an objective function, the mathematical formulation of the topology optimization model is developed using the solid isotropic microstructures with penalization (SIMP) interpolation scheme. The topology optimization problem is solved by the optimality criteria method. Finally, the feasibility and efficiency of the proposed method are illustrated with several 2D examples that are widely used in the topology optimization design.     

A topology optimization approach using VOF method

Topology optimization methods with continuous design variables obtained by the homogenization formula or the solid isotropic microstructure with penalty (SIMP) model are widely used in the layout of structures. In the implementation of these approaches, one must take into account several issues, e.g., irregularity of the problem, occurrence of the checkerboard pattern, and intermediate density. To suppress these phenomena, the employment of additional strategies such as the perimeter control or the filtering method will be required. In this paper, a topology optimization method which can eliminate these difficulties is developed based on the volume of fluid (VOF) method. In the method, shape design is described in terms of the VOF function. Since this function is defined by a volume fraction of material occupying each element, it can be recognized as a continuous material density in the SIMP model. Within the framework of the VOF analysis, the topology optimization procedure is reduced to a convection motion of the material density governed by a Hamilton–Jacobi equation as in the level set method. Through numerical examples, the validity of the proposed method is investigated.     

Structural topology optimization of multibody systems

Flexible multibody dynamics (FMD) has found many applications in control, analysis and design of mechanical systems. FMD together with the theory of structural optimization can be used for designing multibody systems with bodies which are lighter, but stronger. Topology optimization of static structures is an active research topic in structural mechanics. However, the extension to the dynamic case is less investigated as one has to face serious numerical difficulties. One way of extending static structural topology optimization to topology optimization of dynamic flexible multibody system with large rotational and transitional motion is investigated in this paper. The optimization can be performed simultaneously on all flexible bodies. The simulation part of optimization is based on an FEM approach together with modal reduction. The resulting nonlinear differential-algebraic systems are solved with the error controlled integrator IDA (Sundials) wrapped into Python environment by Assimulo (Andersson et al. in Math. Comput. Simul. 116(0):26–43, 2015). A modified formulation of solid isotropic material with penalization (SIMP) method is suggested to avoid numerical instabilities and convergence failures of the optimizer. Sensitivity analysis is central in structural optimization. The sensitivities are approximated to circumvent the expensive calculations. The provided examples show that the method is indeed suitable for optimizing a wide range of multibody systems. Standard SIMP method in structural topology optimization suggests stiffness penalization. To overcome the problem of instabilities and mesh distortion in the dynamic case we consider here additionally element mass penalization.     

On the influence of model reduction techniques in topology optimization of flexible multibody systems using the floating frame of reference approach

Employing the floating frame of reference formulation in the topology optimization of dynamically loaded components of flexible multibody systems seems to be a natural choice. In this formulation the deformation of flexible bodies is approximated by global shape functions, which are commonly obtained from finite element models using model reduction techniques. For topology optimization these finite element models can be parameterized using the solid isotropic material with penalization (SIMP) approach. However, little is known about the interplay of model reduction and SIMP parameterization. Also securing the model reduction quality despite major changes of the design during the optimization has not been addressed yet. Thus, using the examples of a flexible frame and a slider-crank mechanism this work discusses the proper choice of the model reduction technique in the topology optimization of flexible multibody systems.     

Topology optimization for light-trapping structure in solar cells

The limitation associated with the low optical absorption remains to be the main technical barrier that constrains the efficiency of thin–film solar cells in energy conversion. Effective design of light-trapping structure is critical to increase light absorption, which is a highly complex phenomenon governed by several competing physical processes, imposing a number of challenges to topology optimization. This paper presents a general, yet systematic approach exploiting topology optimization for designing highly efficient light-trapping structures. We first demonstrate the proposed approach using genetic algorithm (GA) based non-gradient topology optimization (NGTO), which is robust for achieving highly-efficient designs of slot-waveguide based cells with both low-permittivity and high-permittivity scattering material at single wavelength or over a broad spectrum. The optimized light-trapping structure achieves a broadband absorption efficiency of 48.1 % and more than 3-fold increase over the Yablonovitch limit. The fabrication feasibility of the optimized design is also demonstrated. Next, the gradient topology optimization (GTO) approach for designing light-trapping structure is explored based on the Solid Isotropic Material with Penalization (SIMP) method. Similar designs are obtained through both GA based NGTO and SIMP based GTO, which verifies the validity of both approaches. Insights into the application of both approaches for solving the nanophotonic design problem with optimization nonlinearity are provided.     

基于无网格局部Petrov-Galekin法的板结构拓扑优化

为将无网格法的优势集成到结构拓扑优化中,基于无网格局部Petrov-Galerkin(Meshless Local Petrov-Galerkin,MLPG)法进行板结构的拓扑优化.基于带惩罚的各向同性固体微结构(Solid Isotropic Microstructure with Penalization,SIMP...     

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.     

POD-assisted strategies for structural topology optimization

We propose a new numerical tool for structural optimization design. To cut down the computational burden typical of the Solid Isotropic Material with Penalization (SIMP) method, we apply Proper Orthogonal Decomposition on SIMP snapshots computed on a fixed grid to construct a rough structure (predictor) which becomes the input of a SIMP procedure performed on an anisotropic adapted mesh (corrector). The benefit of the proposed design tool is to deliver smooth and sharp layouts which require a contained computational effort before moving to the 3D printing production phase.     

A discrete level-set topology optimization code written in Matlab

This paper presents a compact Matlab implementation of the level-set method for topology optimization. The code can be used to minimize the compliance of a statically loaded structure. Simple code modifications to extend the code for different and multiple load cases are given. The code is inspired by a Matlab implementation of the solid isotropic material with penalization (SIMP) method for topology optimization (Sigmund, Struct Multidiscipl Optim 21:120–127, 2001). Including the finite element solver and comments, the code is 129 lines long. The code is intended for educational purposes, and in particular it could be used alongside the Matlab implementation of the SIMP method for topology optimization to demonstrate the similarities and differences between the two approaches.     

Topology optimization for hybrid additive-subtractive manufacturing

Hybrid additive-subtractive manufacturing is gaining popularity by making full use of geometry complexity produced by additive manufacturing and dimensional accuracy derived from subtractive machining. Part design for this hybrid manufacturing approach has been done by trial-and-error, and no dedicated design methodology exists for this manufacturing approach. To address this issue, this work presents a topology optimization method for hybrid additive and subtractive manufacturing. To be specific, the boundary segments of the input design domain are categorized into two types: (i) Freeform boundary segments freely evolve through the casting SIMP method, and (ii) shape preserved boundary segments suppress the freeform evolvement and are composed of machining features through a feature fitting algorithm. Given the manufacturing strategy, the topology design is produced through additive manufacturing and the shape preserved boundary segments will be processed by post-machining. This novel topology optimization algorithm is developed under a unified SIMP and level set framework. The effectiveness of the algorithm is proved through a few numerical case studies.     

Maximization of the fundamental eigenfrequency of micropolar solids through topology optimization

Aim of this work is the maximization of the fundamental eigenfrequency of 2D bodies made of micropolar (or Cosserat) materials using a topology optimization approach. A classical SIMP–like model is used to approximate the constitutive parameters of the micropolar medium. A suitable penalization is introduced for both the linear and the spin inertia of the material, to avoid the occurrence of undesired local modes. The robustness of the proposed procedure is investigated through numerical examples; the influence of the material parameters on the optimal material layouts is also discussed. The optimal layouts for Cosserat solids may differ significantly from the truss–like solutions typical of Cauchy solids, as the intrinsic flexural stiffness of the material can lead to curved beam-like material distributions. The numerical simulations show that the results are quite sensitive to the material characteristic length and the spin inertia.     

Multi-material topology optimization using ordered SIMP interpolation

In this paper an ordered multi-material SIMP (solid isotropic material with penalization) interpolation is proposed to solve multi-material topology optimization problems without introducing any new variables. Power functions with scaling and translation coefficients are introduced to interpolate the elastic modulus and the cost properties for multiple materials with respect to the normalized density variables. Besides a mass constraint, a cost constraint is also considered in compliance minimization problems. A heuristic updating scheme of the design variables is derived from the Kuhn-Tucker optimality condition (OC). Since the proposed method does not rely on additional variables to represent material selection, the computational cost is independent of the number of materials considered. The iteration scheme is designed to jump across the discontinuous point of interpolation derivatives to make stable transition from one material phase to another. Numerical examples are included to demonstrate the proposed method. Because of its conceptual simplicity, the proposed ordered multi-material SIMP interpolation can be easily embedded into any existing single material SIMP topology optimization codes.     

Sufficiency of a finite exponent in SIMP (power law) methods

A common way to perform discrete optimization in shape or topology optimization is to use a method called the artificial power law or SIMP. The focus of this paper is to show that this method gives a discrete solution under some conditions. Examples from topology optimization are included for illustrative purposes. Received December 22, 1999     

Design of absorbing material distribution for sound barrier using topology optimization

A topology optimization approach based on the boundary element method (BEM) and the optimality criteria (OC) method is proposed for the optimal design of sound absorbing material distribution within sound barrier structures. The acoustical effect of the absorbing material is simplified as the acoustical impedance boundary condition. Based on the solid isotropic material with penalization (SIMP) method, a topology optimization model is established by selecting the densities of absorbing material elements as design variables, volumes of absorbing material as constraints, and the minimization of sound pressure at reference surface as design objective. A smoothed Heaviside-like function is proposed to help the SIMP method to obtain a clear 0–1 distribution. The BEM is applied for acoustic analysis and the sensitivities with respect to design variables are obtained by the direct differentiation method. The Burton–Miller formulation is used to overcome the fictitious eigen-frequency problem for exterior boundary-value problems. A relaxed form of OC is used for solving the optimization problem to find the optimal absorbing material distribution. Numerical tests are provided to illustrate the application of the optimization procedure for 2D sound barriers. Results show that the optimal distribution of the sound absorbing material is strongly frequency dependent, and performing an optimization in a frequency band is generally needed.     

Finite prism method based topology optimization of beam cross section for buckling load maximization

The use of the finite element method (FEM) for buckling topology optimization of a beam cross section requires large numerical cost due to the discretization in the length direction of the beam. This investigation employs the finite prism method (FPM) as a tool for linear buckling analysis, reducing degrees of freedom of three-dimensional nodes of FEM to those of two-dimensional nodes with the help of harmonic basis functions in the length direction. The optimization problem is defined as the maximization problem of the lowest eigenvalue, for which a bound variable is introduced and set as the design objective to treat mode switching phenomena of multiple eigenvalues. The use of the bound formulation also helps the proposed optimization to treat beams having local plate buckling modes as the fundamental modes as well as beams having global buckling modes. The axial stress is calculated according to the distribution of material modulus which is interpolated using the SIMP approach. Optimization problems finding cross-section layouts from rectangular, L-shaped and generally-shaped design domains are solved for various beam lengths to ascertain the effectiveness of the proposed method.     

A new multi-objective programming scheme for topology optimization of compliant mechanisms

This paper presents an alternative method in implementing multi-objective optimization of compliant mechanisms in the field of continuum-type topology optimization. The method is designated as “SIMP-PP” and it achieves multi-objective topology optimization by merging what is already a mature topology optimization method—solid isotropic material with penalization (SIMP) with a variation of the robust multi-objective optimization method—physical programming (PP). By taking advantages of both sides, the combination causes minimal variation in computation algorithm and numerical scheme, yet yields improvements in the multi-objective handling capability of topology optimization. The SIMP-PP multi-objective scheme is introduced into the systematic design of compliant mechanisms. The final optimization problem is formulated mathematically using the aggregate objective function which is derived from the original individual design objectives with PP, subjected to the specified constraints. A sequential convex programming method, the method of moving asymptotes (MMA) is then utilized to process the optimization evolvement based on the design sensitivity analysis. The main findings in this study include distinct advantages of the SIMP-PP method in various aspects such as computation efficiency, adaptability in convex and non-convex multi-criteria environment, and flexibility in problem formulation. Observations are made regarding its performance and the effect of multi-objective optimization on the final topologies. In general, the proposed SIMP-PP method is an appealing multi-objective topology optimization scheme suitable for “real world” problems, and it bridges the gap between standard topological design and multi-criteria optimization. The feasibility of the proposed topology optimization method is exhibited by benchmark examples.     

基础简谐激励下的结构动响应拓扑优化

针对关于结构动响应拓扑优化问题的研究较少、有限元分析软件的拓扑优化模块无法实现的问题,采用变密度法研究连续体结构在基础简谐激励下的动响应拓扑优化.将基础简谐激励下的响应控制问题归结为结构在体积约束下目标点响应幅值最小化的优化模型;推导有阻尼结构在基础简谐激励下目标点响应幅值的灵敏度公式;采用变密度法求解该优化问题.采用多项式惩罚模型解决带惩罚的各向同性固体微结构（Solid Isotropic Microstructure with Penalization,SIMP）模型带来的附属效应现象;采用灰度过滤方法改善经典变密度法在优化过程中灰度单元收敛过慢的问题,从而减少变密度法优化的迭代步数并且使优化结果更清晰.以平面悬臂板模型为例,验证该优化方法对目标点响应幅值的优化以及灰度过滤函数对优化迭代的改善.     

A note on the theoretical convergence properties of the SIMP method

The solid isotropic material with penalization (SIMP) method is used in topology optimization to solve problems where the variables are 0 or 1. The theoretical convergence properties have not been exhaustively studied. In this paper a convergence theorem with weaker assumptions than earlier conditions is given.     

Concurrent topology optimization of multiscale structures with multiple porous materials under random field loading uncertainty

Aim of this work is the synthesis of auxetic structures using a topology optimization approach for micropolar (or Cosserat) materials. A distributed compliant mechanism design problem is formulated, adopting a SIMP–like model to approximate the constitutive parameters of 2D micropolar bodies. The robustness of the proposed approach is assessed through numerical examples concerning the optimal design of structures that can expand perpendicularly to an applied tensile stress. The influence of the material characteristic length on the optimal layouts is investigated. Depending on the inherent flexural stiffness of micropolar solids, truss–like solutions typical of Cauchy solids are replaced by curved beam–like material distributions. No homogenization technique is implemented, since the proposed design approach applies to elements made of microstructured material with prescribed properties and not to the material itself.