• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.247-250
• /
• 2007

Flexible media such as the paper, the film, etc. are thin, light and very flexible. They behave in geometrically nonlinear. Any of small force makes large deformation. So we must including aerodynamic effect when its behavior is predicted. Thus, it becomes fully coupled fluid-structure interaction(FSI) problem. In FSI problems, where the fluid mesh near the structure undergoes large deformations and becomes unacceptably distorted, which drive the time step to a very small value for explicit calculations, the arbitrary Lagrangian-Eulerian(ALE) methods or rezoning are used to create a new undistorted mesh for the fluid domain, which allows the calculations to continue. In this paper, FE sheet model considering geometric nonlinearity is formulated to simulate the behavior of the flexible media. Aerodynamic force to the media by surrounding air is calculated by solving the incompressible Navier-Stokes equations. Q2Q1(Taylor-Hood) element which means biquadratic for velocity and bilinear for pressure is used for fluid domain. Q2Q1 element satisfies LBB condition and any stabilization technique is not needed. In this paper, cantilevered sheet in the viscous incompressible Navier-Stokes flow is simulated to check the mesh motion and numerical integration scheme, and then falling paper in the air is simulated and the effects of some representative parameters are investigated.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.1
• /
• pp.1-8
• /
• 2005

For the simulation of tsunami propagation an active dispersion-correction two-dimensional finite element model has been developed based on a shallow-water wave equation. This model employs an arbitrary triangular mesh and an explicit time integration scheme. However, the physical dispersion effects as included in the Boussinesq equations can be taken into account in the computation. The validity of the dispersion-correction scheme developed in this study is verified through the comparison of numerical solutions calculated using the new scheme with analytical ones considering dispersion effect of waves. As a result, the present model is shown to be considerably accurate.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2000.04a
• /
• pp.195-204
• /
• 2000

In this paper a complicated structural behavior in collision and its effect of energy translation to the collision bulkhead was examined through a methodology of the numerical simulation to obtain a ideal bow construction and a location of collision bulkhead against head on collision. In the present the bow structure is normally designed in consideration of its specific structural arrangements and internal and external loads in these area such as hydrostatic and dynamic pressure, wave impact and bottom slamming in accordance with the Classification rules, and the specific location of collision bulkhead by SOLAS requirement. By these studies the behavior of the bow collapse due to collision was synthetically evaluated for the different size of tankers and its operational speed limits, and by the result of these simulation it provides the optimal design concept for the low construction to prevent the subsequent plastic deformation onto or near to the collision bulkhead boundary and to determine the rational location of collision bulkhead.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.7
• /
• pp.1426-1437
• /
• 2002

Using the Lagrangian explicit time-integration finite element code NET3D which can treat three-dimensional high-velocity impact problems, oblique penetration processes of long rod projectile with yaw against thin plate are simulated. Through the comparison of simulation result with experimental result and other code's computational result, the adaptability and accuracy of NET3D is evaluated under the complex situation in which yaw angle and oblique angle exist simultaneously. Main research contents to be handled in this paper include the followings. First, the accuracy and efficiency estimation of NET3D code result obtained from the oblique penetration simulations of long rod projectile with yaw against thin plate. Second, the effect of increasing impact velocity. Third, the effect of initial yaw for the spaced-plate target. Residual velocities, residual lengths, angular velocities, and final deformed configurations obtained from the NET3D computations are compared with the experimental results and other code's computational results such as Eulerian code MESA and Lagrangian code EPIC. As a result of comparisons, it has been found that NET3D code is superior to EPIC code and MESA code in the prediction capability of residual velocity and residual length of penetrator. The key features obtained from the experiment can be successfully reproduced through NET3D simulations. Throughout the study, the applicability and accuracy of NET3D as a metallic armor system design tool is verified.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.21 no.10
• /
• pp.1589-1597
• /
• 1997

The high-velocity impact forging process with eccentric loading condition is analyzed using the explicit time integration finite element method. In order to consider the strain hardening, strain rate hardening and thermal softening effects, which are frequently observed in high-velocity deformation phenomena, the Johnson-Cook constitutive model is applied to model the workpiece. It is assumed that the material response of the dies is elastic in the study. As a result of the eccentric loading simulation, it is found that the increase of the eccentric ratio and the allowable tilting angle cause the decrease of the maximum forging load and the blow efficiency, and it is also found that the forging load and the blow efficiency generated in the high-velocity impact forging process with three-dimensional geometry can be obtained efficiently.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.14 no.5
• /
• pp.424-431
• /
• 2004

An implicit direct-time integration method for obtaining transient responses of general dynamic systems is described. The conventional Newmark method cannot be directly applied to state-space first-order differential equations, which contain no explicit acceleration terms. The method proposed here is the state-space Newmark method that incorporates the average velocity concept, and can be applied to an analysis of general dynamic systems that are expressed by state-space first-order differential equations. It is also readily coded into a program. Stability and accuracy analyses indicate that the method is numerically unconditionally stable like the conventional Newmark method, and has a period error of 2nd-order accuracy for small damping and 4th-order for large damping and an amplitude error of 2nd-order, regardless of damping. In addition, its utility and validity are confirmed by two application examples. The results suggest that the proposed state-space Newmark method based on average velocity be generally applied to the analysis of transient responses of general dynamic systems with a high degree of reliability with respect to stability and accuracy.

• Journal of Ship and Ocean Technology
• /
• v.5 no.3
• /
• pp.36-44
• /
• 2001

The collision energy absorbing characteristics of side structure of the LNG carriers which have the cargo containment systems of the spherical and the membrane types are compared. A failure mechanism of the double hull side structures of 130, 000 $m^3$ class LNG carriers under sideways collision event has been simulated by using the detailed finite element calculations. In ship collision analysis, the finite element method based on explicit time integration has been use[1 with much success. Finite element modeling techniques for detail description of structural members antral ship motion regarding the dynamic behavior allowed to investigate the effect of bow shape and the initial contact position on side shell of collided ship. In the numerical simulations of the ship-to-ship sideways collision, the effect of the colliding bow shapes and the change of the colliding ship draft are investigated. The critical collision energy which is absorbed by a side structure of a collided ship until the fore-end of colliding ship arrives at the boundary of the cargo tank is calculated. The critical speed of specified colliding ships which can not penetrate the boundary of the LNG cargo tank of the collided ship under collision accident if evaluated.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.14 no.5
• /
• pp.100-106
• /
• 2006

NiTi alloy known as its shape memory effect also has superelastic characteristic, which makes it possible to be elastic under large deformation. Since the tensile strength of the alloy is very high and density is low compared to carbon steel, it can be applied to lightweight structural design. In order to design structures with shape memory alloy, finite element analysis is used and a constitutive algorithm based on Aurrichio's model is added to LS-DYNA as a user subroutine. Explicit time integration and shell element formulation are used to simulate thin-walled structures. The algorithm uses Drucker-Prager type loading condition to calculate martensite volume fraction during the transformation. The implemented algorithm is verified in uni-axial loading condition and martensite phase transformation can be detected well with the algorithm. In this study, as a energy absorbing structure, thin-walled tube is modeled with finite elements and the deformation behavior is studied. Simulation results has shown that the martensite transformation was generated in loading condition. After plastic deformation reached, the load decreases linearly without reverse martensite transformation.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.9 no.5
• /
• pp.995-1001
• /
• 2005

Implicit Integration has become a standard approach to efficient cloth animation, and it guarantees the stability of the system so that large steps can be used. Therefore, it is regarded as the best method for the real-time or interactive animation of cloth. Since the implicit method was introduced for stable cloth animation, various cloth animation techniques based on the method have been proposed. It is now possible to generate the real-time animation of cloth model with thousands of mass-point in general PC environments. Although the implicit method guarantees the stability, the implementation of the implicit method is generally more difficult than that of the explicit method. Even worse, it is very difficult to parallelize the computation process of the implicit method. The cloth animation with implicit method can be formalized as a linear system solving. In this paper we propose an stable and efficient cloth animation techniques based on the implicit method. The proposed method can be easily parallelized. Self-collision is another important issue in cloth animation, we also propose an efficient self-collision avoidance techniques.