• Journal of Korean Association for Spatial Structures
• /
• v.12 no.1
• /
• pp.77-85
• /
• 2012

The membrane structure maintains stable form by giving initial tension to ductile membrane and increasing the stiffness of exterior that is much adopted in the large span spatial structure by making its thickness thin. This kind of membrane structure has characteristic that can express free-form curve, so the selection of structural form is very important. So, this paper proposes the expression of free-form surface based on NURBS basis function and the finite element method considering geometrical nonlinearity for the deduction of large deformation result. Also, for minimizing the approximation of the surface that is derived from the form-finding result, the interface method that change finite element mesh to NURBS is proposed. So, the optimum surface of free-form membrane is derived.

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

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.9 no.3
• /
• pp.292-304
• /
• 2017

High-speed vessels are prone to the surf-riding in adverse quartering seas. The possibility of mitigating the surf-riding of the ITTC A2 fishing vessel in the design stage is investigated using the 6-DOF weakly non-linear model developed for surf-riding simulations in quartering seas. The longitudinal position of the ship's center of buoyancy (LCB) is chosen as the design parameter. The adjusting of LCB is achieved by changing frame area curves, and hull surfaces are reconstructed accordingly using the Radial Basis Function (RBF). Surf-riding motions in regular following seas for cases with different LCBs and Froude numbers are simulated using the numerical model. Results show that the surf-riding cannot be prevented by the adjusting of LCB. However, it occurs with a higher threshold speed when ship's center of buoyancy (COB) is moved towards stem compared to moving towards stern, which is mainly due to the differences on wave resistance caused by the adjusting of LCB.

• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.35-40
• /
• 2010

In the present work, LES with new variational multiscale method is conducted on the fully developed channel flow with Reynolds number, 180 based on the friction velocity and the channel half width. Incompressible Navier-Stokes equations are integrated using finite element method with the basis function of NURBS. To solve space-time equations, Newton's method with two stage predictor multicorrector algorithm is employed. The code is parallelized using MPI. The computational domain is a rectangular box of size $2{\pi}{\times}2{\times}4/3{\pi}$ in the streamwise, wall normal and spanwise direction. Mean velocity profiles and velocity fluctuations are compared with the data of DNS. The results agree well with those of DNS and other traditional LES.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.49-55
• /
• 2011

Turbulent flows with wavy wall are simulated using Residual-based Variational Multiscale Method (RB-VMS) which is proposed by Bazilves et al(2007) as new Large Eddy Simulation methodology. Incompressible Navier-Stokes equations are integrated using Isogeometric analysis which adopt the basis function as NURBS. The Reynolds number is 6760 based on the bulk velocity and averaged channel height. And the amplitude (${\alpha}/{\lambda}$) of wavy wall is 0.05. The computational domain is $2{\lambda}{\times}1.05{\lambda}{\times}{\lambda}$ in the streamwise, wall normal and spanwise direction. Mean quantities and turbulent statistics near wavy wall are compared with DNS results of Cherukat et al.(1998). The predicted results show good agreement with reference data.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.56-59
• /
• 2009

In the present work, LES with new variational multiscale method is conducted on the fully developed channel flow with Reynolds number is 180 based on the friction velocity and the channel half width. Incompressible Navier-Stokes equations are integrated using finite element method with the basis function of NURBS. To solve space-time equations, Newton's method with two stage predictor multicorretor algorithm is employed. The code is parallelized using MPI. The computational domain is a rectangular box of size $2{\pi}{\times}2{\times}4/3{\pi}$ in the streamwise, wall normal and spanwise direction. Mean velocity profiles and velocity fluctuations are compared with the data of DNS. The results agree well with those of DNS and other traditional LES.

• Architectural research
• /
• v.13 no.3
• /
• pp.41-47
• /
• 2011

A study on the free vibration analysis of elastic bar is described in this paper. In order to determine the natural frequencies of bars, a bar element is developed by using isogeometric formulation. The B-spline is introduced to represent the geometry of bar and the same geometric definition is also used to define its unknown displacement field in isogeometric formulation. Therefore, the stiffness and mass matrices are derived by the order-free B-spline basis function. The efficiency and accuracy of the present isogeometric bar elementis demonstrated by using several numerical tests. From numerical results, it is found to be that the present isogeometric element produces very accurate natural frequencies of bars. Finally, the present isogeometric solutions are provided as future reference solutions.

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

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

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