• Journal of the Korean Society for Precision Engineering
• /
• v.22 no.12 s.177
• /
• pp.117-123
• /
• 2005

This paper deals with the finite element model of the monitor fur the simulation of directional free drop tests such as the backward and vertex free drop as well as the associated free-drop experiment. The model was made for an unconditional stable solution fur the explicit integration algorithm. It was found through the comparison between simulation and experimental results that the general behaviors at the time of impact were observed to well correlate with each other in terms of acceleration, displacement, contact force and stress of monitor components.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.60-67
• /
• 1998

This paper deals with the development of a variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined elements and efficiently used for construction of a refined mesh without generating distorted elements. A modified Gaussian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• The Sea:JOURNAL OF THE KOREAN SOCIETY OF OCEANOGRAPHY
• /
• v.12 no.3
• /
• pp.133-141
• /
• 2007

A finite element model is developed for the extended Boussinesq equations that is capable of simulating the dynamics of long and short waves. Galerkin weighted residual method and the introduction of auxiliary variables for 3rd spatial derivative terms in the governing equations are used for the model development. The Adams-Bashforth-Moulton Predictor Corrector scheme is used as a time integration scheme for the extended Boussinesq finite element model so that the truncation error would not produce any non-physical dispersion or dissipation. This developed model is applied to the problems of solitary wave propagation. Predicted results is compared to available analytical solutions and laboratory measurements. A good agreement is observed.

• Structural Engineering and Mechanics
• /
• v.43 no.3
• /
• pp.349-370
• /
• 2012

Variable-node finite element families, termed (4 + k + l + m + n)-node elements with an arbitrary number of nodes (k, l, m, and n) on each of their edges, are developed based on the generic point interpolation with special bases having slope discontinuities in two-dimensional domains. They retain the linear interpolation between any two neighboring nodes, and passes the standard patch test when subdomain-wise $2{\times}2$ Gauss integration is employed. Their shape functions are automatically generated on the master domain of elements although a certain number of nodes are inserted on their edges. The elements can provide a flexibility to resolve nonmatching mesh problems like mesh connection and adaptive mesh refinement. In the case of adaptive mesh refinement problem, so-called "1-irregular node rule" working as a constraint in performing mesh adaptation is relaxed by adopting the variable-node elements. Through several examples, we show the performance of the variable-node finite elements in terms of accuracy and efficiency.

• Transactions of the Society of Information Storage Systems
• /
• v.5 no.1
• /
• pp.19-35
• /
• 2009

The present work is an experimental and analytical study on a flexible disk rotating close to a rigid rotating disk in open air. In the analytical study, the air flow in the gap between the flexible disk and the rigid disk is modeled using Navier-Stokes and continuity equations while the flexible disk is modeled using the linear plate theory. The flow equations are discretized using the cell centered finite volume method (FVM) and solved numerically with semi-implicit pressure-linked equations (SIMPLE algorithm). The spatial terms in the disk equation are discretized using the finite difference method (FDM) and the time integration is performed using fourth-order Runge-Kutta method. An experimental test-rig is designed to investigate the dynamics of the flexible disk when rotating close to a co-rotating, a counter-rotating and a fixed rigid disk, which works as a stabilizer. The effects of rotational speed, initial gap height and inlet-hole radius on the flexible disk displacement and its vibration amplitude are investigated experimentally for the different types of stabilizer. Finally, the analytical and experimental results are compared.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.5
• /
• pp.388-399
• /
• 2016

A motion analysis of an underwater towed cable is a complex task due to its nonlinear nature of the problem. The major source of the nonlinearity of the underwater cable analysis is that the motion of the cable involves large rigid-body motion. This large rigid-body motion makes difficult to use standard displacement-based finite element method. In this paper, the authors apply recently developed nodal position-based finite element method which can deal with the geometric nonlinearity due to the large rigid-body motion. In order to enhance the stability of the large-scale nonlinear cable motion simulation, an efficient time-integration scheme is proposed, namely predictor/multi-corrector Newmark scheme. Three different predictors are introduced, and the best predictor in terms of stability and robustness for impulsive cable motion analysis is proposed. As a result, the nonlinear motion of underwater cable is predicted in a very efficient manner compared to the classical finite element of finite difference methods. The efficacy of the method is demonstrated with several test cases, involving static and dynamic motion of a single cable element, and also under water towed cable composed of multiple cable elements.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.10
• /
• pp.1793-1799
• /
• 2003

An optimal design of a multi-layered plate structure to endure high-velocity impact has been suggested by using size optimization after numerical simulations. The NET2D, a Lagrangian explicit time-integration finite element code for analyzing high-velocity impact, was used to find the parameters for the optimization. Three different materials such as mild steel, aluminum for a multi-layered plate structure and die steel for the pellet, were assumed. In order to consider the effects of strain rate hardening, strain hardening and thermal softening, Johnson-Cook model and Phenomenological Material Model were used as constitutive models for the simulation. It was carried out with several different gaps and thickness of layers to figure out the trend in terms of those parameters' changes under the constraint, which is against complete penetration. Also, the measuring domain has been shrunk with several elements to reduce the analyzing time. The response surface method based on the design of experiments was used as optimization algorithms. The optimized thickness of each layer in which perforation does not occur has been obtained at a constant velocity and a designated total thickness. The result is quite acceptable satisfying both the minimized deformation energy and the weight criteria. Furthermore, a conceptual idea for topology optimization was suggested for the future work.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• Journal of Mechanical Science and Technology
• /
• v.20 no.11
• /
• pp.1790-1800
• /
• 2006

Dynamic behavior of flexural-torsional coupled vibration of rotating beams using the Rayleigh-Ritz method with orthogonal polynomials as basis functions is studied. Performance of various orthogonal polynomials is compared to each other in terms of their efficiency and accuracy in determining the required natural frequencies. Orthogonal polynomials and functions studied in the present work are: Legendre, Chebyshev, integrated Legendre, modified Duncan polynomials, the special trigonometric functions used in conjunction with Hermite cubics, and beam characteristic orthogonal polynomials. A total of 5 cases of beam boundary conditions and rotation are studied for their natural frequencies. The obtained natural frequencies and mode shapes are compared to those available in various references and the results for coupled flexural-torsional vibrations are especially compared to both previously available references and with those obtained using NASTRAN finite element package. Among all the examined orthogonal functions, Legendre orthogonal polynomials are the most efficient in overall CPU time, mainly because of ease in performing the integration required for determining the stiffness and mass matrices.