• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.299-306
• /
• 2006

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The 'Hamilton-Jacobi (H-J)' equation and computationally robust numerical technique of 'up-wind scheme' lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H -J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes is not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2008.04a
• /
• pp.391-396
• /
• 2008

A level set based topological shape optimization using extended B-spline basis functions is developed for steady state heat conduction problems. The only inside of complicated domain is identified by the level set functions and taken into account in computation. The solution of Hamilton-Jacobi equation leads to an optimal shape according to the normal velocity field determined from the sensitivity analysis, minimizing a thermal compliance while satisfying a volume constraint. To obtain exact shape sensitivity, the precise normal and curvature of geometry need to be determined using the level set and B-spline basis functions. The nucleation of holes is possible whenever and wherever necessary during the optimization using a topological derivative concept.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.1
• /
• pp.79-87
• /
• 2013

An isogeometric topological shape optimization method is developed using the level sets and Heaviside enrichments. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set functions, which facilitates to handle complicated topological shape changes. The Heaviside enrichment improves the isogeometric analysis by adding some enrichment functions to model the internal boundaries. The proposed topological shape optimization method has several benefits: exact geometric models can be obtained using the isogeometric approach and the limitation of tensor-product patches can be overcome using the Heaviside enrichments to represent the internal voids. Even in a single patch, discontinuous displacement fields as well as smooth stress field can be obtained. Since the level sets offer the implicit moving boundary inside the domain, it is easy to represent the topological shape variations in the isogeometric analysis using Heaviside enrichments.

• The Transactions of The Korean Institute of Electrical Engineers
• /
• v.58 no.10
• /
• pp.1936-1941
• /
• 2009

This paper deals with a new methodology for topology optimization in which the topology of the design domain may change during the shape optimization process. To achieve this, the concept of the topological gradient is introduced to compute the sensitivity of an objective function when a small hole is drilled in the domain. Based on shape and topological sensitivity values, the shape and topology of the design domain may be simultaneously changed during design iterations if necessary. To verify the advantages and also to facilitate understanding of the method itself, two electrostatic design problems have been tested by using 2D finite element analysis: the first is the inverse problem of a simple dielectric model and the second is the rotor design of a MEMS actuator.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.3
• /
• pp.239-245
• /
• 2008

A level set based topological shape optimization using extended B-spline basis functions is developed for steady-state heat conduction problems. The only inside of complicated domain identified by the level set functions is taken into account in computation, so we can remove the effects of domain outside parts in heat conduction problem. The solution of Hamilton-Jacobi equation leads to an optimal shape according to the normal velocity field determined from the sensitivity analysis, minimizing a thermal compliance while satisfying a volume constraint. To obtain exact shape sensitivity, the precise normal and curvature of geometry need to be determined using the level set and B-spline basis functions. Using topological derivative concept, the nucleation of holes for topological changes can be made whenever and wherever necessary during the optimization.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.9-16
• /
• 2014

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The "Hamilton-Jacobi(H-J)" equation and computationally robust numerical technique of "up-wind scheme" lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H-J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes are not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.

• Korean Journal of Computational Design and Engineering
• /
• v.12 no.5
• /
• pp.376-381
• /
• 2007

The geometric configuration in the weight-reduced structure is very required to be started from the conceptual design with low cost, high performance and quality. In this point, a structural-topological shape concerned with conceptual design of structure is important. The method used in this paper combines three optimization techniques, where the shape and physical dimensions of the structure and material distribution are hierachically optimized, with the maximum rigidity of structure and lightweight.

• Proceedings of the KIEE Conference
• /
• 2008.10c
• /
• pp.23-25
• /
• 2008

This paper presents a topological shape optimization for electromagnetic system using a Level Set method. The optimization is progressed by updating the implicit Level Set function from the Hamilton-Jacobi equation. The up-wind scheme is used for numerical implementation of the Hamilton-Jacobi equation. In order to validate the proposed optimization, the core part of a C-core actuator is optimized by three cases using different materials; (single steel), (two steels), and (steel and magnet).

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.6
• /
• pp.559-567
• /
• 2012

A level set based topological shape optimization method for nonlinear structure considering hyper-elastic problems is developed. To relieve significant convergence difficulty in topology optimization of nonlinear structure due to inaccurate tangent stiffness which comes from material penalization of whole domain, explicit boundary for exact tangent stiffness is used by taking advantage of level set function for arbitrary boundary shape. For given arbitrary boundary which is represented by level set function, a Delaunay triangulation scheme is used for current structure discretization instead of using implicit fixed grid. The required velocity field in the actual domain to update the level set equation is determined from the descent direction of Lagrangian derived from optimality conditions. The velocity field outside the actual domain is determined through a velocity extension scheme based on the method suggested by Adalsteinsson and Sethian(1999). The topological derivatives are incorporated into the level set based framework to enable to create holes whenever and wherever necessary during the optimization.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2008.11a
• /
• pp.343-347
• /
• 2008

Design of shape for visco-elastic vibration isolation elements, which are very cost-effective and so popular in many applications is frequently based on experiences, intuitions, or trial and errors. Such traditions in shape design make it difficult for drastic changes or new concepts to come out. In this paper, both topological method and shape optimization method are combined together to find out a most desirable isolator shape efficiently by using two commercial engineering programs. ABAQUS and MATLAB. The procedure is divided into two steps. At the first step, a topology optimization method is employed to find an initial shape. where density of either 0 or 1 for finite elements is used for physical realizability. At the second step, based on the initial shape, finer tuning of the shape is done by boundary movement method. An illustration of the procedure will be presented for a mount of an air-conditioner compressor system and the effectiveness will be discussed.