• 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 Korean Association for Spatial Structures
• /
• v.12 no.1
• /
• pp.109-119
• /
• 2012

For membrane structure, there are three main steps in design and construction, which are form finding, statistical load analysis, and cutting patterning. Unlike the first two stages, the step of cutting pattern involves the translation of a double-curved surface in 3D space into a 2D plane with minimal error. For economic reasons, the seam lines of generated cutting patterns rely greatly on the geodesic line. Generally, as searching regions of the seam line are plane elements in the step of shape analysis, the seam line is not a smooth curve, but an irregularly divided straight line. So, it is how we make an irregularly divided straight line a smooth curve that defines the quality of the pattern. Accordingly, in this paper, we analyzed interpolation schemes using spline, and apply these methods to cutting pattern generation on the curved surface. To generate the pattern, three types of spline functions were used, i.e., cubic spline function, B-spline, and least-square spline approximation, and simple model and the catenary-shaped membrane was adopted to examine the result of generation. The result of comparing the approximation curves by the number of elements and the number of extracted nodes of simple model revealed that the seam line for less number of extracted nodes with large number of elements were more efficient, and the least-square spline approximation provided smoother seam line than other methods.

• Geomechanics and Engineering
• /
• v.2 no.2
• /
• pp.141-156
• /
• 2010

In the framework of meshfree methods, a new methodology is developed based on radial point interpolation method (RPIM). This methodology is applied to a one-dimensional contaminant transport modelling in the saturated porous media. The one-dimensional form of advection-dispersion equation involving reactive contaminant is considered in the analysis. The Galerkin weak form of the governing equation is formulated using 1D meshfree shape functions constructed using thin plate spline radial basis functions. MATLAB code is developed to obtain the numerical solution. Numerical examples representing various phenomena, which occur during migration of contaminants, are presented to illustrate the applicability of the proposed method and the results are compared with those obtained from the analytical and finite element solutions. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints. In order to test the practical applicability and performance of the RPIM, three case studies of contaminant transport through the landfill liners are presented. A good agreement is obtained between the results of the RPIM and the field investigation data.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.9 no.6
• /
• pp.87-94
• /
• 2010

In this paper, a new simulator of the biped robot for planning walking patterns was showed. And this simulator(MHBiped) is able to not only visualize the plan of patterns but also verify whether a biped robot design is suitable. In addition, MHBipd can modify various kinds of walking parameters and the trajectory of biped robot. Therefore, a new biped robot can be designed easily by the this simulator before you apply to a robot. As a result, a well-balanced parameters of walking patterns watching the movement of CG and ZMP can be obtained. Walking patterns should be changed according to both the existence of obstacles and conditions of ground and it can be described by the trajectory of hip and ankles. All those trajectorys can be also obtained by the cubic spline functions and the way of modeling walking patterns. The results of simulator, the movement function of CG and ZMP, the cubic spline functions and modeling of biped robot were introduced in this paper. And the effectiveness of this simulator was confirmed by the simulations.

• Journal of the Korean Data and Information Science Society
• /
• v.20 no.1
• /
• pp.89-96
• /
• 2009

We examined the characteristics of temperature variation in Seoul between 1961 to 2008 using the spectral heteroscedastic model. The mean function in the propsed model explains the season effect using periodic functions and the overall increase using the quadratic regression spline. The variance function also had periodic functions to explain the seasonality of variance. We found that there has been annual mean temperature increase by about $1.5^{\circ}C$ for the last 48 years. The increase of annual mean temperature was mainly caused by the increase in winter, which made the amplitude decreased.

• Korean Journal of Computational Design and Engineering
• /
• v.13 no.6
• /
• pp.440-449
• /
• 2008

This paper deals with a method calculating various hull form parameters which are required in numerical analysis for ship performance such as motion, maneuverability, resistance and propulsion, etc. After the hull form is designed, before the model tests the ship's performances are evaluated by various analysis tools in which the hull form parameters are used with many kinds of forms aside from offset data. Here, The hull form parameters characterize the properties of hull form and contain positional, differential and integral information implicitly. Generally, the commercial CAD-system has not functions enough for supporting these form parameters and therefore each shipyard uses its own in-house analysis program as well as commercial analysis software. To overcome these limitations, modules for supporting these analysis programs have developed. The modules contain cubic composite spline cure using local curve fairing, intersect algorithm, Gaussian integral, and other geometric techniques needed in calculating hull form parameters. Using our analysis-supporting modules, a complex hull form can be remodeled exactly to the hull form designed by CAD-system and any hull form parameter required in various performance analyses can be calculated.

• Structural Engineering and Mechanics
• /
• v.38 no.6
• /
• pp.733-751
• /
• 2011

In the present study, a new kind of multivariable Reissner-Mindlin plate elements with two kinds of variables based on B-spline wavelet on the interval (BSWI) is constructed to solve the static and vibration problems of a square Reissner-Mindlin plate, a skew Reissner-Mindlin plate, and a Reissner-Mindlin plate on an elastic foundation. Based on generalized variational principle, finite element formulations are derived from generalized potential energy functional. The two-dimensional tensor product BSWI is employed to form the shape functions and construct multivariable BSWI elements. The multivariable wavelet finite element method proposed here can improve the solving accuracy apparently because generalized stress and strain are interpolated separately. In addition, compared with commonly used Daubechies wavelet finite element method, BSWI has explicit expression and a very good approximation property which guarantee the satisfying results. The efficiency of the proposed multivariable Reissner-Mindlin plate elements are verified through some numerical examples in the end.

• Journal of the Society of Naval Architects of Korea
• /
• v.42 no.2 s.140
• /
• pp.113-123
• /
• 2005

A higher order panel method based on B-spline representation for both the geometry and the solution is developed for the analysis of steady flow around marine propellers. The self-influence functions due to the normal dipole and the source are desingularized through the quadratic transformation, and then shown to be evaluated using conventional numerical quadrature. By selecting a proper order for numerical quadrature, the accuracy of the present method can be increased to the machine limit. The far- and near-field influences are shown to be evaluated based on the same far-field approximation, but the near-field solution requires subdividing the panels into smaller subpanels continuously, which can be effectively implemented due to the B-spline representation of the geometry. A null pressure jump Kutta condition at the trailing edge is found to be effective in stabilizing the solution process and in predicting the correct solution. Numerical experiments indicate that the present method is robust and predicts the pressure distribution on the blade surface, including very close to the tip and trailing edge regions, with far fewer panels than existing low order panel methods.

• Korean Journal of Applied Biomechanics
• /
• v.14 no.3
• /
• pp.301-311
• /
• 2004

The purpose of this study was to verify the usefulness of the Mathcad program as a tool for the studying and teaching the sport biomechanics. A projectile motion was analyzed because it is the one of the most popular motion in sports activities. A 3 dimensional CG data for the high jump bar clear phase was used to calculate the initial velocity vector of the CG. Linear regression function and other functions such as cubic spline and derivative of Mathcad were used to calculate this vector. Finally, the approach angle to the bar and peak jump height was calculated. Programming in Mathcad was relatively easy compare to traditional computer language such as Fortran and C, because of the unique documentation method of Mathcad. Additionally the 2 and 3 dimensional graph function was very easy and useful to describe the mechanical data. If the use of Mathcad program is more popular in the field of sport biomechanics, it could greatly contribute to overcome the limit of research caused by the lack of proper programming ability.

• East Asian mathematical journal
• /
• v.34 no.5
• /
• pp.661-681
• /
• 2018

In this paper, we discuss a reduced-order modeling for the Benjamin-Bona-Mahony-Burgers (BBMB) equation and its application to a distributed feedback control problem through the centroidal Voronoi tessellation (CVT). Spatial distcritization to the BBMB equation is based on the finite element method (FEM) using B-spline functions. To determine the basis elements for the approximating subspaces, we elucidate the CVT approaches to reduced-order bases with snapshots. For the purpose of comparison, a brief review of the proper orthogonal decomposition (POD) is provided and some numerical experiments implemented including full-order approximation, CVT based model, and POD based model. In the end, we apply CVT reduced-order modeling technique to a feedback control problem for the BBMB equation.