• 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.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.3
• /
• pp.329-335
• /
• 2011

The finite element method (FEM) has been used for various fields like mathematics and engineering. However, the FEM has a difficulty in describing the geometric shape exactly due to its property of piecewise linear discretization. Recently, however, a so-called isogeometric analysis method that uses the non-uniform rational B-spline(NURBS) basis function has been developed. The NURBS can be used to describe the geometry exactly and play a role of basis functions for the response analysis. Nevertheless, constructing the NURBS basis functions in analysis is as costly as a meshing process in the FEM. Since the isogeometric method shares geometric data with CAD, it is possible to intactly import the model data from commercial CAD tools. In this paper, we use the Rhinoceros 3D software to create CAD models and export in the form of STEP file. The information of knot vectors and control points in the NURBS is utilized in the isogeometric analysis. Through some numerical examples, the accuracy of isogeometric method is compared with that of FEM. Also, the efficiency of the isogeometric method that includes the CAD and CAE in a unified framework is verified.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.4
• /
• pp.255-262
• /
• 2013

In this paper, a shape design sensitivity analysis(DSA) method is presented for Mindlin plates using an isogeometric approach. The isogeometric method possesses desirable advantages; the representation of exact geometry and the higher order inter-element continuity, which lead to the fast convergence of solution as well as accurate sensitivity results. Unlike the finite element methods using linear shape functions, the isogeometric method considers the exact normal vector and curvature of the CAD geometry, taking advantages of higher order NURBS basis functions. A selective reduced integration(SRI) technique is incorporated to overcome the difficulty of 'shear locking' phenomenon. This simple technique is surprisingly helpful for the accuracy of the isogeometric shape sensitivity without complicated formulation. Through the numerical examples of plate bending problems, the accuracy of the proposed isogeometric analysis method is compared with that of finite element one. Also, the isogeometric shape sensitivity turns out to be very accurate when compared with finite difference sensitivity.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2009.04a
• /
• pp.408-411
• /
• 2009

본 논문에서는 아이소-지오메트릭 해석 방법을 사용하여 응력 제한 조건이 있는 형상 최적설계 문제를 다룬다. 아이소-지오메트릭 해석 방법은 해석에 사용되는 기저 함수와 기하 모델을 구성하는 함수가 일치하여 기하학적으로 정확하기 때문에 설계민감도 해석 및 형상 최적설계에 있어서 강점이 있다. 많은 최적화 문제에서 최대 강성을 확보하는 방향으로 최적화가 진행되고 있는데 이때 응력 조건을 고려하지 않는 경우가 대부분이다. 응력 제한조건이 있는 구조물에서 아이소-지오메트릭 형상 최적설계를 적용시켜 봄으로써 그 효용성을 확인하였다.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.3
• /
• pp.275-281
• /
• 2010

In this paper, the design optimization of structures with stress constraints is performed using isogeometric shape optimization method. The stress constraints have an important role in design optimization problems since stress concentration could result in structural failure. To represent exact geometry in analysis, the isogeometric analysis method uses the same basis functions as used in the CAD geometry. The geometrically exact model can be used in both stress and design sensitivity analyses so that it can yield more precise optimal design than finite element one. Through numerical examples, the isogeometric approach turns out to be effective in shape optimization problems under stress constraints.

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

Finite element analysis is to approximate a geometry model developed in computer-aided design(CAD) to a finite element model, thus the conventional shape design sensitivity analysis and optimization using the finite element method have some difficulties in the parameterization of geometry. However, isogeometric analysis is to build a geometry model and directly use the functions describing the geometry in analysis. Therefore, the geometric properties can be embedded in the NURBS basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. In this study, the isogeometric structural analysis and shape design sensitivity analysis in the generalized curvilinear coordinate(GCC) systems are discussed for the curved geometry. Representing the higher order geometric information, such as normal, tangent and curvature, yields the isogeometric approach to be the best way for generating exact GCC systems from a given CAD geometry. The developed GCC isogeometric structural analysis and shape design sensitivity analysis are verified to show better accuracy and faster convergency by comparing with the results obtained from the conventional isogeometric method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.1
• /
• pp.1-7
• /
• 2011

In this paper, based on an isogeometric approach, we have developed a shape design optimization method for plane elasticity problems subjected to design-dependent loads. The conventional shape optimization using the finite element method has some difficulties in the parameterization of geometry. In an isogeometric analysis, however, the geometric properties are already embedded in the B-spline basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. The solution space for the response analysis can be represented in terms of the same NURBS basis functions to represent the geometry, which yields a precise analysis model that exactly represents the normal and curvature depending on the applied loads. A continuum-based isogeometric adjoint sensitivity is extensively derived for the plane elasticity problems under the design-dependent loads. Through some numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.6
• /
• pp.515-522
• /
• 2017

In this paper, an isogeometric analysis for multipatch problem is investigated, in which two or more geometries are connected at the interface in a conforming or non-conforming conditions. To express higher continuity at the patch interface, two approaches such as Nitsche based method and master-slave method are formulated for the linear elasticity problem and discretized using the isogeometric approach using NURBS basis functions. A short comparison between two approaches in formulations reveals the pros and cons of them with the applicability in the isogeometric multipatch problem. In addition, a NURBS based stress recovery is adopted to express a better stress continuity through the post-processing. Numerical examples indicate the effectiveness of Nitsche method in the non-conforming patch, following the exact solution well. For the stress concentration problem with the conforming patch, introduced two methodologies show comparative results, meanwhile the NURBS based stress recovery presents an improved smooth stress contour in the whole domain including the patch interface.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.155-162
• /
• 2014

Using an isogeometric approach, a continuum-based shape design optimization method is developed for steady state power flow problems at high frequencies. In case the isogeometric method is employed to the shape design optimization, the NURBS basis functions used in CAD geometric modeling are directly utilized to embed the exact geometry into the computational framework so that the design parameterization for shape optimization is much easier than that in the finite element method and consequently provides the enhanced smoothness of design perturbations. Thus, exact geometric models can be used in both the response and the shape sensitivity analyses, where normal vector and curvature are continuous over the whole design space so that enhanced shape sensitivity can be expected. Through numerical examples, the developed isogeometric sensitivity is compared with finite difference one to provide excellent agreement. Also, it turns out that the proposed method works very well in the shape optimization problems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.5
• /
• pp.345-352
• /
• 2014

In this paper, the CAD data for the optimal shape design obtained by isogeometric shape optimization is directly used to fabricate the specimen by using 3D printer for the experimental validation. In a conventional finite element method, the geometric approximation inherent in the mesh leads to the accuracy issue in response analysis and design sensitivity analysis. Furthermore, in the finite element based shape optimization, subsequent communication with CAD description is required in the design optimization process, which results in the loss of optimal design information during the communication. Isogeometric analysis method employs the same NURBS basis functions and control points used in CAD systems, which enables to use exact geometrical properties like normal vector and curvature information in the response analysis and design sensitivity analysis procedure. Also, it vastly simplify the design modification of complex geometries without communicating with the CAD description of geometry during design optimization process. Therefore, the information of optimal design and material volume is exactly reflected to fabricate the specimen for experimental validation. Through the design optimization examples of elasticity problem, it is experimentally shown that the optimal design has higher stiffness than the initial design. Also, the experimental results match very well with the numerical results. Using a non-contact optical 3D deformation measuring system for strain distribution, it is shown that the stress concentration is significantly alleviated in the optimal design compared with the initial design.