• The Pure and Applied Mathematics
• /
• v.19 no.2
• /
• pp.147-169
• /
• 2012

We introduce a Heaviside-function formulation of the interaction between incompressible two-phase fluid and a non-deformable solid. Fluid and solid interact in two ways : fluid satises the Dirichlet boundary condition imposed by the velocity field of solid, and solid is accelerated by the surface traction exerted by fluid. The two-way couplings are formulated by the Heaviside function to the interface between solid and fluid. The cumbersome treatment of interface is taken care of by the Heaviside function, and the interaction is discretized in a simple manner. The discretization results in a stable and accurate projection method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.55-61
• /
• 2002

In the present study, a first order shear deformable Loop-subdivision triangular element which can handle transverse shear deformation of moderately thick shell and composite laminated or sandwich shells are developed. The developed element is more general than the previous one based on classical shell theory, since it includes the effect of transverse shell deformation and has standard five degrees of freedom per node. The quartic box spline function is employed as the interpolation basis function. Numerical examples for the benchmark static shell problems are analyzed to assess the performance of the developed subdivision shell element and locking trouble.

• Journal of the Korean Society of Visualization
• /
• v.18 no.3
• /
• pp.35-41
• /
• 2020

We conduct numerical simulations of the interaction of a deformable structure with two-phase compressible flow. The finite volume method (FVM) is used to simulate fluid phenomena including a shock wave, a gas bubble, and the deformation of free surface. The deformation of a floating structure is computed with the finite element method (FEM). The compressible two-phase volume of fluid (VOF) method is used for the generation and development of a cavitation bubble, and the immersed boundary method (IBM) is used to impose the effect of the structure on the fluid domain. The result of the simulation shows the generation of a shock wave, and the expansion of the bubble. Also, the deformation of the structure due to the hydrodynamic loading by the explosion is identified.

• Structural Engineering and Mechanics
• /
• v.44 no.1
• /
• pp.61-72
• /
• 2012

The paper is on contact analysis of a spherical ball with a deformable flat, considering the effect of tangent modulus on the contact parameters of a non-adhesive frictionless elastic-plastic contact. The contact analysis of this model has been carried out using analysis software Ansys and Abaqus. The contact parameters such as area of contact between two consecutive steps, volume of bulged material are evaluated from the formulated equations. The effect of the tangent modulus is considered for determining these parameters. The tangent modulus are accounted between 0.1E and 0.5E of materials E/Y value greater than 500 and less than 1750. Result shows that upto an optimal tangent modulus values the elastic core push up to the free surface in the flat. The simulation is also carried out in Abaqus and result provide evidence for the volume of bulged material in the contact region move up and flow into the free surface of the flat from the contact edge between the ball and flat. The strain energy of the whole model is varied between 20 to 40 percentage of the stipulated time for analysis.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.2 no.1
• /
• pp.23-34
• /
• 1998

A dynamic fluid-structure-soil interaction analysis method is developed to investigate the effects of translational and/or rocking motions on the seismic response of flexible rectangular liquid storage tanks founded on the deformable ground. The governing equation for the dynamics of 3-D rectangular tanks subjected to the translational and/or rocking motion is abtained by applying Rayleigh-Ritz method. The dynamic stiffness matrices of a rigid rectangular foundation resting on the surface of a stratum overlaid bedrock are calculated by hyperelement method. The seismic responses of 3-D flexible tank model founded on the deformable ground is calculated by combining the governing equation for the fluid-tank system with the dynamic stiffness matrix of th rigid surface foundation.

• Journal of International Society for Simulation Surgery
• /
• v.1 no.1
• /
• pp.37-40
• /
• 2014

Purpose For living donor liver transplantation, liver segmentation is difficult due to the variability of its shape across patients and similarity of the density of neighbor organs such as heart, stomach, kidney, and spleen. In this paper, we propose an automatic segmentation of the liver using multi-planar anatomy and deformable surface model in portal phase of abdominal contrast-enhanced CT images. Method Our method is composed of four main steps. First, the optimal liver volume is extracted by positional information of pelvis and rib and by separating lungs and heart from CT images. Second, anisotropic diffusing filtering and adaptive thresholding are used to segment the initial liver volume. Third, morphological opening and connected component labeling are applied to multiple planes for removing neighbor organs. Finally, deformable surface model and probability summation map are performed to refine a posterior liver surface and missing left robe in previous step. Results All experimental datasets were acquired on ten living donors using a SIEMENS CT system. Each image had a matrix size of $512{\times}512$ pixels with in-plane resolutions ranging from 0.54 to 0.70 mm. The slice spacing was 2.0 mm and the number of images per scan ranged from 136 to 229. For accuracy evaluation, the average symmetric surface distance (ASD) and the volume overlap error (VE) between automatic segmentation and manual segmentation by two radiologists are calculated. The ASD was $0.26{\pm}0.12mm$ for manual1 versus automatic and $0.24{\pm}0.09mm$ for manual2 versus automatic while that of inter-radiologists was $0.23{\pm}0.05mm$. The VE was $0.86{\pm}0.45%$ for manual1 versus automatic and $0.73{\pm}0.33%$ for manaual2 versus automatic while that of inter-radiologist was $0.76{\pm}0.21%$. Conclusion Our method can be used for the liver volumetry for the pre-surgery planning of living donor liver transplantation.

• 한국HCI학회:학술대회논문집
• /
• 2008.02a
• /
• pp.39-44
• /
• 2008

This paper proposes a real-time geometric correction system based on a projector to project digital images onto deformable surface. Markers use to trace lots of corresponding points would spoil the projected image when the projector projects a digital image onto the surface because they leave marks on the surface. In addition, it is difficult to build a real-time geometric correction system since bottlenecks occur through the process of the geometric correction for projecting images. In this paper, we use invisible infrared markers and a vertex shader of GPU using Cg TookKit of NVIDIA in order to eliminate disadvantage and bottlenecks in the process of markers recognition so that it is possible to project natural correction images in real-time. As a result, this system overlays an interactive virtual texture onto the real paper by using the geometric transformation. Therefore, it is possible to develop variation of AR(Augmented Reality) based on digital contents systems.

• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.12
• /
• pp.148-158
• /
• 2003

In this paper, efficient and accurate contact search algorithm is proposed for the contact problems by the finite element method. A slave node and a maser contact segment is defined using the side of a finite element on the contact surface. The specific goal is to develop techniques of reducing the nonsmoothness of the contact interactions arising from the finite element discretization of the contact surface. Contact detection is accomplished by monitoring the territory of the slave nodes throughout the calculation for possible penetration of a master surface. To establish the validity of the proposed algorithm, some different process and geometries examples were simulated. Efforts are focused on the error rate that is based on the penetrated area through the simulations fur large deformation with contact surface between deformable bodies. A proposed algorithm offers improvements in contact detection from the simulation results.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.450-455
• /
• 2004

There is few research about contact problem for a rigid surface with an arbitrary shape in SPH. The variational equation based on the virtual work principle is derived and its solution is obtained by the penalty method. It is proposed a new method that can determine the parameters for a penetration and a penetration rate used in the penalty method. The reproducing condition is adopted to correct the deficiency of kernel on the boundary. In order to calculate a penetration of particles, after checking boundary particles for deformable body boundary normal vectors were determined on the rigid surface. Numerical simulations for models which have rigid surface with an arbitrary shape were conducted to validate the proposed method in 2D. The results of those analysis represent that the contact algorithm proposed in this study works properly.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.1 s.232
• /
• pp.30-37
• /
• 2005

There is few research about contact problem for a rigid surface with an arbitrary shape in SPH. The variational equation based on the virtual work principle is derived and its solution is obtained by the penalty method. It is proposed a new method that can determine the parameters for a penetration and a penetration rate used in the penalty method. The reproducing condition is adopted to correct the deficiency of kernel on the boundary. In order to calculate a penetration of particles, after checking boundary particles for deformable body, boundary normal vectors were determined on the rigid surface. Numerical simulations for models which have rigid surface with an arbitrary shape were conducted to validate the proposed method in 2D Cartesian and cylindrical coordinate. The results of those analysis represent that the contact algorithm proposed in this study works properly.