• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.33 no.3
• /
• pp.93-100
• /
• 2021

Since the IB (Immersed Boundary) method, which can perform coupling analysis with objects and fluids having an impermeable boundary of arbitrary shape on a fixed grid system, has been developed, the IB method in various CFD models is increasing. The representative IB methods are the directing-forcing method and the ghost cell method. The directing-forcing type method numerically satisfies the boundary condition from the fluid force calculated at the boundary surface of the structure, and the ghost-cell type method is a computational method that satisfies the boundary condition through interpolation by placing a virtual cell inside the obstacle. These IB methods have a disadvantage in that the computational algorithm is complex. In this study, the simplified immersed boundary (SIB) method enables the analysis of temporary structures on a fixed grid system and is easy to expand to three proposed dimensions. The SIB method proposed in this study is based on a one-field model for immiscible two-phase fluid that assumes that the density function of each phase moves with the center of local mass. In addition, the volume-weighted average method using the density function of the solid was applied to handle moving solid structures, and the CIP method was applied to the advection calculation to prevent numerical diffusion. To examine the analysis performance of the proposed SIB method, a numerical simulation was performed on an object falling to the free water surface. The numerical analysis result reproduced the object falling to the free water surface well.

• Structural Engineering and Mechanics
• /
• v.55 no.5
• /
• pp.1055-1086
• /
• 2015

The dynamic response of two-dimensional unbounded domain on the rigid bedrock in the time domain is numerically obtained. It is realized by the modified scaled boundary finite element method (SBFEM) in which the original scaling center is replaced by a scaling line. The formulation bases on expanding dynamic stiffness by using the continued fraction approach. The solution converges rapidly over the whole time range along with the order of the continued fraction increases. In addition, the method is suitable for large scale systems. The numerical method is employed which is a combination of the time domain SBFEM for far field and the finite element method used for near field. By using the continued fraction solution and introducing auxiliary variables, the equation of motion of unbounded domain is built. Applying the spectral shifting technique, the virtual modes of motion equation are eliminated. Standard procedure in structural dynamic is directly applicable for time domain problem. Since the coefficient matrixes of equation are banded and symmetric, the equation can be solved efficiently by using the direct time domain integration method. Numerical examples demonstrate the increased robustness, accuracy and superiority of the proposed method. The suitability of proposed method for time domain simulations of complex systems is also demonstrated.

• Computers and Concrete
• /
• v.4 no.4
• /
• pp.285-298
• /
• 2007

Corrosion of Reinforcement (rebar) is nondestructively estimated by the half-cell potential measurement. As is the case with other nondestructive testings (NDT), understanding of the underlying principles should be clarified in order to obtain meaningful results. Therefore, the measurement of potentials in concrete is analytically investigated. The effect of internal defects on the potentials measured is clarified numerically by the boundary element method (BEM). Thus, a simplified inversion by BEM is applied to convert the potentials on concrete surface to those on rebars, taking into account the concrete resistivity. Because the potentials measured on concrete surface are so sensitive to moisture content, concrete resistivity and surface condition, an inverse procedure to convert the potentials on concrete surface into those on rebars is developed on the basis of BEM. It is found that ASTM criterion is practically applicable to estimate corrosion from the potential values converted. In experiments, an applicability of the procedure is examined by accelerated corrosion tests of reinforced concrete (RC) slabs. For practical use, the procedure is developed where results of IBEM are visualized by VRML (Virtual Reality modeling Language) in three-dimensional space.

• Proceedings of the Korean Geotechical Society Conference
• /
• 2005.03a
• /
• pp.25-32
• /
• 2005

Modeling techniques of piled pier were reviewed and the influences of pile cap's boundary conditions were analyzed in this study. Among various modeling techniques, equivalent cantilever method seems relatively simple for modeling pile groups and it has some problems to determine the virtual fixed points. Through the analyses, it was found that the method of nonlinear p-y model with soil springs was more appropriate than equivalent cantilever method. The method modeling a pile group using stiffness matrix seems useful for practical design, which can represent the nonlinear three-dimensional behavior of a piled pier. In this study, the stiffness matrix of a pile group could be estimated efficiently and precisely using three-dimensional nonlinear analysis programs of pile groups (FBPier 3.0, YSGroup).

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.2
• /
• pp.221-228
• /
• 1989

The objective of this paper is to develop a numerical technique for analyzing crack driving forces in multi-phase materials. The analysis was based on finite element method coupled with a virtual crack extension technique which is known as the most efficient tool in computational fracture mechanics analysis. The modified J-integral method, proposed by Miyamoto and Kikuchi for the analysis of dual-phase material was carried out by subtracting the J-values for contours surrounding each phase boundary from the J-values for overall contour. It was shown that the proposed numerical procedure, based on the modified J-integral coupled with a virtual crack extension technique, can be used as an effective numerical tool for determining crack driving forces in multi-phase materials.

• Journal of Ship and Ocean Technology
• /
• v.10 no.2
• /
• pp.24-36
• /
• 2006

Ship hullform is usually designed with a curve network, and smooth hullform surfaces are supposed to be generated by filling in (or interpolating) the curve network with appropriate surface patches. Tensor-product surfaces such as B-spline and $B\'{e}zier$ patches are typical representations to this interpolating problem. However, they have difficulties in representing the surfaces of irregular topological type which are frequently appeared in the fore- and after-body of ship hullform curve network. In this paper, we proposed a method that can automatically generate discrete $G^1$ continuous B-spline surfaces interpolating given curve network of ship hullform. This method consists of three steps. In the first step, given curve network is reorganized to be of two types: boundary curves and reference curves of surface patches. Especially, the boundary curves are specified for their surface patches to be rectangular or triangular topological type that can be represented with tensor-product (or degenerate) B-spline surface patches. In the second step, surface fitting points and cross boundary derivatives are estimated by constructing virtual iso-parametric curves at discrete parameters. In the last step, discrete $G^1$ continuous B-spline surfaces are generated by surface fitting algorithm. Finally, several examples of resulting smooth hullform surfaces generated from the curve network data of actual ship hullform are included to demonstrate the quality of the proposed method.

• Journal of the Korean Society for Precision Engineering
• /
• v.23 no.7 s.184
• /
• pp.187-193
• /
• 2006

This paper introduces and illustrates the results of a new method for offsetting the triangular mesh generated from multiple surfaces. The meshes generated from each surface are separated each other and normal directions are different. The face normal vectors are flipped to upward and the lower faces covered by upper faces are deleted. The virtual normal vectors are introduced and used to of feet boundary. It was shown that new method is better than previous methods in offsetting the triangular meshes generated from multiple surfaces. The introduced offset method was applied for 3-axis tool path generation system and tested by NC machining.

• Journal of Korea Multimedia Society
• /
• v.22 no.2
• /
• pp.311-321
• /
• 2019

Virtual reality applications in Head-Mounted Displays require high frame rates and low latency rendering techniques. Ray tracing offers many benefits, such as high-quality image generation, but has not been utilized due to lower performance than rasterization. But that can obtain good result combined with gaze-tracking technology and human visual system's operation principle. In this paper, we propose a method to optimize the foveated sampling map and to maintain the visual quality through the fast voronoi nearest interpolation. The proposed method further reduces the computational cost that has been improved by the previous foveated sampling. It also smoothes the voronoi boundary using adaptive sibson interpolation, which was not possible in real-time. As a result, the proposed method can render real-time high-quality images with low visual difference.

• Geomechanics and Engineering
• /
• v.35 no.5
• /
• pp.465-473
• /
• 2023

The higher order shear deformable model and an exact analytical method is used for analytical bending analysis of a cylindrical shell subjected to mechanical loads, in this work. The shell is modelled using sinusoidal bivariate shear strain theory, and the static governing equations are derived using changes in virtual work. The eigenvalue-eigenvector method is used to exactly solve the governing equations for a constrained cylindrical shell The proposed kinematic relation decomposes the radial displacement into bending, shearing and stretching functions. The main advantage of the method presented in this work is the study of the effect of clamping constraints on the local stresses at the ends. Stress, strain, and deformation analysis of shells through thickness and length.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.14 no.11
• /
• pp.5338-5343
• /
• 2013

In a haptic system, a virtual wall is modeled as a virtual spring. The larger the stiffness of the virtual spring is, the more improved the reality of the virtual wall is, but the more unstable the haptic system becomes. This paper shows how to increase the stiffness of the virtual spring while the stability of the haptic system is guaranteed and shows the effects of a mass (Md) and a damper (Bd) of a haptic device on the stability when first-order hold method is applied and a virtual wall is modeled as a virtual spring (Kw). The simulation results show the boundary of the virtual spring is proportional to the square root of the mass (Md) and the damper (Bd) while maintaining the stability. The relation among the virtual spring (Kw), the mass (Md) and the damper (Bd) of the haptic device, and sampling time (T) is inferred as $K_w{\leq}{1.611M_d}^{0.50}{B_d}^{0.50}T^{-1.51}$, by using the simulation results. The maximum available stiffness of the virtual spring in first-order hold method is larger than in zero-order hold method. So the reality of the virtual wall can be improved.