• 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 Architectural Institute of Korea Structure & Construction
• /
• v.34 no.6
• /
• pp.11-18
• /
• 2018

Isogeometric concept is introduced to find out the optimum layout of plane structure under free vibration. Eigenvalue problem is formulated and numerically solved in order to obtain natural frequencies and mode shapes of plane structures. For the exact geometric expression of the structure, the Non-Uniform Rational B-spline Surface (NURBS) basis functions is employed and it is also used to define the material density functions. A node-wise design variables is adopted to deal with the updating of material density in topology optimization (TO). The definition of modal strain energy is employed to achieve the maximization of fundamental frequency through its minimization. The verification of the proposed TO technique is performed by a series of benchmark test for plane structures.

• Structural Engineering and Mechanics
• /
• v.88 no.1
• /
• pp.83-94
• /
• 2023

An accurate and highly efficient inverse element labelled iPCB is developed based on the inverse finite element method (iFEM) for real-time shape estimation of plane-curved structures (such as arch bridges) utilizing onboard strain data. This inverse problem, named shape sensing, is vital for the design of smart structures and structural health monitoring (SHM) procedures. The iPCB formulation is defined based on a least-squares variational principle that employs curved Timoshenko beam theory as its baseline. The accurate strain-displacement relationship considering tension-bending coupling is used to establish theoretical and measured section strains. The displacement fields of the isoparametric element iPCB are interpolated utilizing nonuniform rational B-spline (NURBS) basis functions, enabling exact geometric modelling even with a very coarse mesh density. The present formulation is completely free from membrane and shear locking. Numerical validation examples for different curved structures subjected to different loading conditions have been performed and have demonstrated the excellent prediction capability of iPCBs. The present formulation has also been shown to be practical and robust since relatively accurate predictions can be obtained even omitting the shear deformation contributions and considering polluted strain measures. The current element offers a promising tool for real-time shape estimation of plane-curved structures.

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

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

Shape design optimization for linear elasticity problem is performed using isogeometric analysis method. In many design optimization problems for real engineering models, initial raw data usually comes from CAD modeler. Then designer should convert this CAD data into finite element mesh data because conventional design optimization tools are generally based on finite element analysis. During this conversion there is some numerical error due to a geometry approximation, which causes accuracy problems in not only response analysis but also design sensitivity analysis. As a remedy of this phenomenon, the isogeometric analysis method is one of the promising approaches of shape design optimization. The main idea of isogeometric analysis is that the basis functions used in analysis is exactly same as ones which represent the geometry, and this geometrically exact model can be used shape sensitivity analysis and design optimization as well. In shape design sensitivity point of view, precise shape sensitivity is very essential for gradient-based optimization. In conventional finite element based optimization, higher order information such as normal vector and curvature term is inaccurate or even missing due to the use of linear interpolation functions. On the other hands, B-spline basis functions have sufficient continuity and their derivatives are smooth enough. Therefore normal vector and curvature terms can be exactly evaluated, which eventually yields precise optimal shapes. In this article, isogeometric analysis method is utilized for the shape design optimization. By virtue of B-spline basis function, an exact geometry can be handled without finite element meshes. Moreover, initial CAD data are used throughout the optimization process, including response analysis, shape sensitivity analysis, design parameterization and shape optimization, without subsequent communication with CAD description.

• Earthquakes and Structures
• /
• v.22 no.5
• /
• pp.527-538
• /
• 2022

In this paper the free vibration frequencies of tri-directional functionally graded materials imperfect plate is investigated for Several plate geometries with two types of porosity (even and uneven) and different type of material configuration. The effect of several parameters such as power law index and boundary conditions have been investigated. For this purpose, an efficient computational method is developed and written under Matlab environment, based on a three-dimensional modeling and the isogeometric method is used for the discretization of the structure based on NURBS (Nonuniform rational B-spline) basis functions. The results obtained by the present method are validated by the comparison with the results given by several authors in the literature.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.5
• /
• pp.287-292
• /
• 2021

We utilized the isogeometric analysis (IGA) method that uses NURBS basis functions in CAD systems, to account for the geometric exactness of a geometrically exact beam deformation, on a new type of metamaterial, twist-translation coupled structure showing a large twist angle. A two-dimensional unit cell structure was embedded in a cylindrical wall, using free-form deformation and an appropriate interpolation scheme. A parametric study on the effects of the dimensions of the cylinder and the number of cells, on the twisting angle was performed. Furthermore, the mechanism of the twist-translation coupled metamaterial was explored through numerical examples.

• Steel and Composite Structures
• /
• v.37 no.6
• /
• pp.663-678
• /
• 2020

This paper proposes a novel topology optimization method generating multiple materials for external linear plane crack structures based on the combination of IsoGeometric Analysis (IGA) and eXtended Finite Element Method (X-FEM). A so-called eXtended IsoGeometric Analysis (X-IGA) is derived for a mechanical description of a strong discontinuity state's continuous boundaries through the inherited special properties of X-FEM. In X-IGA, control points and patches play the same role with nodes and sub-domains in the finite element method. While being similar to X-FEM, enrichment functions are added to finite element approximation without any mesh generation. The geometry of structures based on basic functions of Non-Uniform Rational B-Splines (NURBS) provides accurate and reliable results. Moreover, the basis function to define the geometry becomes a systematic p-refinement to control the field approximation order without altering the geometry or its parameterization. The accuracy of analytical solutions of X-IGA for the crack problem, which is superior to a conventional X-FEM, guarantees the reliability of the optimal multi-material retrofitting against external cracks through using topology optimization. Topology optimization is applied to the minimal compliance design of two-dimensional plane linear cracked structures retrofitted by multiple distinct materials to prevent the propagation of the present crack pattern. The alternating active-phase algorithm with optimality criteria-based algorithms is employed to update design variables of element densities. Numerical results under different lengths, positions, and angles of given cracks verify the proposed method's efficiency and feasibility in using X-IGA compared to a conventional X-FEM.