• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.577-582
• /
• 2007

A variational formulation for plane elasticity problems is derived based on an isogeometric approach. The isogeometric analysis is an emerging methodology such that the basis functions in analysis domain arc generated directly from NURBS (Non-Uniform Rational B-Splines) geometry. Thus. the solution space can be represented in terms of the same functions to represent the geometry. The coefficients of basis functions or the control variables play the role of degrees-of-freedom. Furthermore, due to h-. p-, and k-refinement schemes, the high order geometric features can be described exactly and easily without tedious re-meshing process. The isogeometric sensitivity analysis method enables us to analyze arbitrarily shaped structures without re-meshing. Also, it provides a precise construction method of finite element model to exactly represent geometry using B-spline base functions in CAD geometric modeling. To obtain precise shape sensitivity, the normal and curvature of boundary should be taken into account in the shape sensitivity expressions. However, in conventional finite element methods, the normal information is inaccurate and the curvature is generally missing due to the use of linear interpolation functions. A continuum-based adjoint sensitivity analysis method using the isogeometric approach is derived for the plane elasticity problems. The conventional shape optimization using the finite element method has some difficulties in the parameterization of boundary. In isogeometric analysis, however, the geometric properties arc already embedded in the B-spline shape functions and control points. The perturbation of control points in isogeometric analysis automatically results in shape changes. Using the conventional finite clement method, the inter-element continuity of the design space is not guaranteed so that the normal vector and curvature arc not accurate enough. On tile other hand, in isogeometric analysis, these values arc continuous over the whole design space so that accurate shape sensitivity can be obtained. Through 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.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.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.1
• /
• pp.65-84
• /
• 2012

We present an application of the variational multiscale methodology to the computation of concentric annular pipe flow. Isogeometric analysis is utilized for higher order approximation of the solution using Non-Uniform Rational B-Splines (NURBS) functions. The ability of NURBS to exactly represent curved geometries makes NURBS-based isogeometric analysis attractive for the application to the flow through the curved channel.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.20 no.3
• /
• pp.339-345
• /
• 2007

In this paper, a variational formulation for plane elasticity problems is derived based on an isogeometric approach. The isogeometric analysis is an emerging methodology such that the basis functions for response analysis are generated directly from NURBS (Non-Uniform Rational B-Splines) geometry. Furthermore, the solution space for the response analysis can be represented in terms of the same functions to represent the geometry, which enables to provide a precise construction method of finite element model to exactly represent geometry using B-spline base functions in CAD geometric modeling and analyze arbitrarily shaped structures without re-meshing. In this paper, a continuum-based adjoint sensitivity analysis method using the isogeometric approach is extensively derived for the plane elasticity problems. The conventional shape optimization using the finite element method has some difficulties in the parameterization of geometry In the 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. Through some numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

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

• Structural Engineering and Mechanics
• /
• v.48 no.1
• /
• pp.125-150
• /
• 2013

This paper intends to present an application of isogeometric analysis in crack problems. An isogeometric formula is developed based on NURBS basis functions - enriched and adopted via X-FEM enrichment functions. The proposed method which is represented by the combination of the two above-mentioned methods, first by using NURBS functions models the geometry exactly and then by defining level set function on domain, identifies available discontinuity in elements. Additional DOFs are allocated to elements containing the crack and X-FEM enrichment functions enrich approximate solution. Moreover, a subelement refinement technique is used to improve the accuracy of integration by the Gauss quadrature rule. Finally, several numerical examples are illustrated to demonstrate the effectiveness, robustness and accuracy of the proposed method during calculation of crack parameters.

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

• Steel and Composite Structures
• /
• v.21 no.6
• /
• pp.1389-1410
• /
• 2016

This paper proposes a hybrid of topological derivative-based level set method (LSM) and isogeometric analysis (IGA) for structural topology optimization. In topology optimization a significant drawback of the conventional LSM is that it cannot create new holes in the design domain. In this study, the topological derivative approach is used to create new holes in appropriate places of the design domain, and alleviate the strong dependency of the optimal topology on the initial design. Furthermore, the values of the gradient vector in Hamilton-Jacobi equation in the conventional LSM are replaced with a Delta function. In the topology optimization procedure IGA based on Non-Uniform Rational B-Spline (NURBS) functions is utilized to overcome the drawbacks in the conventional finite element method (FEM) based topology optimization approaches. Several numerical examples are provided to confirm the computational efficiency and robustness of the proposed method in comparison with derivative-based LSM and FEM.

• Structural Engineering and Mechanics
• /
• v.76 no.4
• /
• pp.435-449
• /
• 2020

This work presents the formulation of the isogeometric collocation method to solve the strong form equation of a unified and integrated approach of Reissner Mindlin plate theory (UI-RM). In this plate theory model, the total displacement is expressed in terms of bending and shear displacements. Rotations, curvatures, and shear strains are represented as the first, the second, and the third derivatives of the bending displacement, respectively. The proposed formulation is free from shear locking in the Kirchhoff limit and is equally applicable to thin and thick plates. The displacement field is approximated using the B-splines functions, and the strong form equation of the fourth-order is solved using the collocation approach. The convergence properties and accuracy are demonstrated with square plate problems of thin and thick plates with different boundary conditions. Two approaches are used for convergence tests, e.g., increasing the polynomial degree (NELT = 1×1 with p = 4, 5, 6, 7) and increasing the number of element (NELT = 1×1, 2×2, 3×3, 4×4 with p = 4) with the number of control variable (NCV) is used as a comparable equivalent variable. Compared with DKMQ element of a 64×64 mesh as the reference for all L/h, the problem analysis with isogeometric collocation on UI-RM plate theory exhibits satisfying results.