A topology optimization algorithm based on topological derivative and level-set function for designing phononic crystals

2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Rolando Yera ◽  
Luisina Forzani ◽  
Carlos Gustavo Méndez ◽  
Alfredo E. Huespe

PurposeThis work presents a topology optimization methodology for designing microarchitectures of phononic crystals. The objective is to get microstructures having, as a consequence of wave propagation phenomena in these media, bandgaps between two specified bands. An additional target is to enlarge the range of frequencies of these bandgaps.Design/methodology/approachThe resulting optimization problem is solved employing an augmented Lagrangian technique based on the proximal point methods. The main primal variable of the Lagrangian function is the characteristic function determining the spatial geometrical arrangement of different phases within the unit cell of the phononic crystal. This characteristic function is defined in terms of a level-set function. Descent directions of the Lagrangian function are evaluated by using the topological derivatives of the eigenvalues obtained through the dispersion relation of the phononic crystal.FindingsThe description of the optimization algorithm is emphasized, and its intrinsic properties to attain adequate phononic crystal topologies are discussed. Particular attention is addressed to validate the analytical expressions of the topological derivative. Application examples for several cases are presented, and the numerical performance of the optimization algorithm for attaining the corresponding solutions is discussed.Originality/valueThe original contribution results in the description and numerical assessment of a topology optimization algorithm using the joint concepts of the level-set function and topological derivative to design phononic crystals.

Author(s):  
Yoshifumi Okamoto ◽  
Hiroshi Masuda ◽  
Yutaro Kanda ◽  
Reona Hoshino ◽  
Shinji Wakao

PurposeThe purpose of this paper is the improvement of topology optimization. The scope of the paper is focused on the speedup of optimization. Design/methodology/approachTo achieve the speedup, the method of moving asymptotes (MMA) with constrained condition of level set function is applied instead of solving the Hamilton–Jacobi equation. FindingsThe acceleration of convergence of objective function is drastically improved by the implementation of MMA. Originality/valueNormally, the level set method is solved through the Hamilton–Jacobi equation. However, the possibility of introducing mathematical programming is clear by the constrained condition. Furthermore, the proposed method is suitable for efficiently solving the topology optimization problem in the magnetic field system.


2021 ◽  
Vol ahead-of-print (ahead-of-print) ◽  
Author(s):  
Meisam Takalloozadeh ◽  
Gil Ho Yoon

Purpose Body forces are always applied to structures in the form of the weight of materials. In some cases, they can be neglected in comparison with other applied forces. Nevertheless, there is a wide range of structures in civil and mechanical engineering in which weight or other types of body forces are the main portions of the applied loads. The optimal topology of these structures is investigated in this study. Design/methodology/approach Topology optimization plays an increasingly important role in structural design. In this study, the topological derivative under body forces is used in a level-set-based topology optimization method. Instability during the optimization process is addressed, and a heuristic solution is proposed to overcome the challenge. Moreover, body forces in combination with thermal loading are investigated in this study. Findings Body forces are design-dependent loads that usually add complexity to the optimization process. Some problems have already been addressed in density-based topology optimization methods. In the present study, the body forces in a topological level-set approach are investigated. This paper finds that the used topological derivative is a flat field that causes some instabilities in the optimization process. The main novelty of this study is a technique used to overcome this challenge by using a weighted combination. Originality/value There is a lack of studies on level-set approaches that account for design-dependent body forces and the proposed method helps to understand the challenges posed in such methods. A powerful level-set-based approach is used for this purpose. Several examples are provided to illustrate the efficiency of this method. Moreover, the results show the effect of body forces and thermal loading on the optimal layout of the structures.


Author(s):  
Long Jiang ◽  
Shikui Chen ◽  
Xiangmin Jiao

The parametric level set method is an extension of the conventional level set methods for topology optimization. By parameterizing the level set function, conventional levels let methods can be easily coupled with mathematical programming to achieve better numerical robustness and computational efficiency. Furthermore, the parametric level set scheme not only can inherit the original advantages of the conventional level set methods, such as clear boundary representation and high topological changes handling flexibility but also can alleviate some un-preferred features from the conventional level set methods, such as needing re-initialization. However, in the RBF-based parametric level set method, it was difficult to determine the range of the design variables. Moreover, with the mathematically driven optimization process, the level set function often results in significant fluctuations during the optimization process. This brings difficulties in both numerical stability control and material property interpolation. In this paper, an RBF partition of unity collocation method is implemented to create a new type of kernel function termed as the Cardinal Basis Function (CBF), which employed as the kernel function to parameterize the level set function. The advantage of using the CBF is that the range of the design variable, which was the weight factor in conventional RBF, can be explicitly specified. Additionally, a distance regularization energy functional is introduced to maintain a desired distance regularized level set function evolution. With this desired distance regularization feature, the level set evolution is stabilized against significant fluctuations. Besides, the material property interpolation from the level set function to the finite element model can be more accurate.


2011 ◽  
Vol 308-310 ◽  
pp. 1046-1049 ◽  
Author(s):  
Yu Wang ◽  
Zhen Luo

This paper proposes a meshless Galerkin level set method for structural shape and topology optimization of continua. To taking advantage of the implicit free boundary representation scheme, structural design boundary is represented through the introduction of a scalar level set function as its zero level set, to flexibly handle complex shape fidelity and topology changes by maintaining concise and smooth interface. Compactly supported radial basis functions (CSRBFs) are used to parameterize the level set function and also to construct the shape functions for mesh free function approximation. The meshless Galerkin global weak formulation is employed to implement the discretization of the state equations. This provides a pathway to simplify two numerical procedures involved in most conventional level set methods in propagating the discrete level set functions and in approximating the discrete equations, by unifying the two different stages at two sets of grids just in terms of one set of scattered nodes. The proposed level set method has the capability of describing the implicit moving boundaries without remeshing for discontinuities. The motion of the free boundary is just a question of advancing the discrete level set function by finding the design variables of the size optimization in time. One benchmark example is used to demonstrate the effectiveness of the proposed method. The numerical results showcase that this method has the ability to simplify numerical procedures and to avoid numerical difficulties happened in most conventional level set methods. It is straightforward to apply the present method to more advanced shape and topology optimization problems.


2019 ◽  
Vol 30 (7) ◽  
pp. 3653-3684
Author(s):  
LanHao Zhao ◽  
Kailong Mu ◽  
Jia Mao ◽  
Khuc Hongvan ◽  
Dawei Peng

Purpose Moving interface problems exist commonly in nature and industry, and the main difficulty is to represent the interface. The purpose of this paper is to capture the accurate interface, a novel three-dimensional one-layer particle level set (OPLS) method is presented by introducing Lagrangian particles to reconstruct the seriously distorted level set function. Design/methodology/approach First, the interface is captured by the level set method. Then, the interface is corrected with only one-layer particles advected with the flow to ensure that the level set function value of the particle is equal to 0. When interfaces are merged, all particles in merged regions are deleted, while the added particles near the generated interface are used to determine the interface as the interface is separated. Findings The OPLS method is validated with well-known benchmark examples, such as the long-term advection of a sphere, the rotation of a three-dimensional slotted disk and sphere, single vortex in a box, sphere merging and separation, deformation of a sphere. The simulation results indicate that the proposed method is found to be highly reliable and accurate. Originality/value This method exhibits excellent conservation of the area bounded by the interface. The extraordinary performance is also shown in dealing with complex interface topological changes.


2018 ◽  
Vol 140 (7) ◽  
Author(s):  
Benliang Zhu ◽  
Rixin Wang ◽  
Hai Li ◽  
Xianmin Zhang

In level-set-based topology optimization methods, the spatial gradients of the level set field need to be controlled to avoid excessive flatness or steepness at the structural interfaces. One of the most commonly utilized methods is to generalize the traditional Hamilton−Jacobi equation by adding a diffusion term to control the level set function to remain close to a signed distance function near the structural boundaries. This study proposed a new diffusion term and built it into the Hamilton-Jacobi equation. This diffusion term serves two main purposes: (I) maintaining the level set function close to a signed distance function near the structural boundaries, thus avoiding periodic re-initialization, and (II) making the diffusive rate function to be a bounded function so that a relatively large time-step can be used to speed up the evolution of the level set function. A two-phase optimization algorithm is proposed to ensure the stability of the optimization process. The validity of the proposed method is numerically examined on several benchmark design problems in structural topology optimization.


Sign in / Sign up

Export Citation Format

Share Document