• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 1997.04a
• /
• pp.128-135
• /
• 1997

The indirect boundary integral equation is formulated to analyze the behavior of a cavity in a multilayered half-plane subjected to earthquake waves. This formulation uses the fundamental solutions that are numerically calculated by the generalized transmission and reflection coefficient method. The free surface of the cavity without external excitation influences the behavior of the half-plane. Consequently this analysis adds the consideration of scattering-field into the analysis and the total motion field of the cavity is decomposed into the free-field and scattering-field motions. The free-field motion is obtained from the modification of the transmission and reflection coefficient method. The scattering-field motion is calculated is calculated by the indirect boundary value problem which has the ficticious boundaries and sources. In this study, P wave, SV wave, SH wave, and Rayleigh wave are analyzed respectively.

• Advances in Computational Design
• /
• v.9 no.1
• /
• pp.73-80
• /
• 2024

The aim of this paper is to remove the rigid body motion in the interior boundary value problem (BVP) of plane elasticity by solving the interior and exterior BVPs simultaneously. First, we formulate the interior and exterior BVPs simultaneously. The tractions applied on the contour in two problems are the same. After adding and subtracting the two boundary integral equations (BIEs), we will obtain a couple of BIEs. In the coupled BIEs, the properties of relevant integral operators are modified, and those integral operators are generally invertible. Finally, a unique solution for boundary displacement of interior region can be obtained.

• Journal of computational fluids engineering
• /
• v.9 no.1
• /
• pp.18-24
• /
• 2004

The unified simulation for the multiphase flow by predictor-corrector scheme based on CIP method is introduced. In this algorithm, the interface between different phases is identified by a density function and tracked by solving an advection equation. Solid body motion is modeled by the translation and angular motion. The mathematical formulation and numerical results are also described. To verify the efficiency, accuracy and capability of proposed algorithm, two dimensional incompressible cavity flow, the motion of a floating ball into water and a single rising bubble by buoyancy force are numerically simulated by the present scheme. As results, it is confirmed that the present scheme gives an efficient, stable and reasonable solution in the multiphase flow problem.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.66-72
• /
• 1999

A numerical simulation method has been under development for solving turbulent flows around a ship model in maneuvering motion using the Reynolds Averaged Navier-Stokes equations. The method used second-order finite differences, collocated grids, pressure-Poisson equation and four-stage Runge-Kutta scheme as key components of the solution method. A modified Baldwin-Lomax model is used for the turbulence closure. This paper presents a preliminary result of the computational study on turbulent flows past a ship model in drift motion. Calculations are carried out for a Series 60 $C_B=0.6$ ship model, for which detailed experimental data are available. The results of the present calculations are compared with the experimental data for hydrodynamic forces acting on the model as well as velocity distributions at longitudinal sections. Only fair agreements has been achieved. The computational results show the complex asymmetrical shear flow patterns including three-dimensional separations followed by formation of bilge vortices both in bow and stern regions.

• Geomechanics and Engineering
• /
• v.2 no.4
• /
• pp.229-252
• /
• 2010

In this paper, Soil-Structure Interaction (SSI) effect is investigated using a new and integrated approach. Faster solution of time dependant differential equation of motion is achieved using numerical representation of wavelet theory while dynamic Infinite Elements (IFE) concept is utilized to effectively model the unbounded soil domain. Combination of the wavelet theory with IFE concept lead to a robust, efficient and integrated technique for the solution of complex problems. A direct method for soil-structure interaction analysis in a two dimensional medium is also presented in time domain using the frequency dependent transformation matrix. This matrix which represents the far field region is constructed by assembling stiffness matrices of the frequency dependant infinite elements. It maps the problem into the time domain where the equations of motion are to be solved. Accuracy of results obtained in this study is compared to those obtained by other SSI analysis techniques. It is shown that the solution procedure discussed in this paper is reliable, efficient and less time consuming as compared to other existing concepts and procedures.

• The Journal of Korea Robotics Society
• /
• v.17 no.3
• /
• pp.288-295
• /
• 2022

This study shows a system identification method of a balancing loading device for a stair climbing delivery robot. The balancing loading device is designed as a 2-axes gimbal structure and is interpreted as two independent pendulum structures for simplifying. The loading device's properties such as mass, moment of inertia, and position of the center of gravity are changeable for luggage. The system identification process of the loading device is required, and the controller should be optimized for the system in real-time. In this study, the system identification method is based on least squares method to estimate the unknown parameters of the loading device's dynamic equation. It estimates the unknown parameters by calculating them that minimize the error function between the real system's motion and the estimated system's motion. This study improves the accuracy of parameter estimation using a null space solution. The null space solution can produce the correct parameters by adjusting the parameter's relative sizes. The proposed system identification method is verified by the simulation to determine how close the estimated unknown parameters are to the real parameters.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.4
• /
• pp.9-17
• /
• 1996

This paper develops a sequential updating method of the Euler parameter generalized coordinates for the machine kinematics and dynamics, The Newton's method is slightly modified so as to utilize the Jacobian matrix with respect to the virtual rotation instead of this with repect to the Euler parameters. An intermediate variable is introduced and the modified Newton's method solves for the variable first. Relational equation of the intermediate variable is then solved for the Euler parameters. The solution process is carried out efficiently by symoblic inversion of the relational equation of the intermediate variable and the iteration equation of the Euler parameter normalization constraint. The proposed method is applied to a kinematic and dynamic analysis with the Generalized Coordinate Partitioning method. Covergence analysis is performed to guarantee the local convergence of the proposed method. To demonstrate the validity and practicalism of the proposed method, kinematic analysis of a motion base system and dynamic analysis of a vehicle are carried out.

• Journal of the Korean Society of Industry Convergence
• /
• v.4 no.1
• /
• pp.87-94
• /
• 2001

In multibody dynamics, differential and algebraic equations which can satisfy both equation of motion and kinematic constraint equation should be solved. To solve this equation, coordinate partitioning method and constraint stabilization method are commonly used. The coordinate partitioning method divides the coordinate into independent and dependent coordinates. The most typical coordinate partitioning method arc LU decomposition, QR decomposition, projection method and SVD(sigular value decomposition).The objective of this research is to find a efficient coordinate partitioning method in flexible multibody systems and a hybrid decomposition algorithm which employs both LU and projection methods is proposed. The accuracy of the solution algorithm is checked with a slider-crank mechanism.

• Proceedings of the Korean Institute of Navigation and Port Research Conference
• /
• 2006.10a
• /
• pp.134-143
• /
• 2006

Out of all types of motions the critical motions leading to capsize is roll. The dynamic amplification in case of roll motion may be large for ships as roll natural frequency generally falls within the frequency range of wave energy spectrum typical used for estimation of motion spectrum. Roll motion is highly non-linear in nature. Den are various representations of non-linear damping and restoring available in literature. In this paper an uncoupled non-linear roll equations with three representation of damping and cubic restoring term is solved using a perturbation technique. Damping moment representations are linear plus quadratic velocity damping, angle dependant damping and linear plus cubic velocity dependant damping. Numerical value of linear damping coefficient is almost same for all types but non-linear damping is different. Linear and non-linear damping coefficients are obtained form free roll decay tests. External rolling moment is assumed as deterministic with sinusoidal form. Maximum roll amplitude of non-linear roll equation with various representations of damping is calculated using analytical procedure and compared with experimental results, which are obtained form forced tests in regular waves by varying frequency with three wave heights. Experiments indicate influence of non-linearity at resonance frequency. Both experiment and analytical results indicates increase in maximum roll amplitude with wave slope at resonance. Analytical results are compared with experiment results which indicate maximum roll amplitude analytically obtained with angle dependent and cubic velocity damping are equal and difference from experiments with these damping are less compared to non-linear equation with quadratic velocity damping.

• Structural Engineering and Mechanics
• /
• v.73 no.1
• /
• pp.83-95
• /
• 2020

In this work, a mathematical model of beam-column system carrying a double eccentric end mass system is investigated, and solved analytically based on the exact solution analysis. The model considers the case in which the double eccentric end mass is a rigid storage tank containing fluid. Both Timoshenko and Bernoulli-Euler beam bending theories are considered. Equation of motion, general solution and boundary conditions for the present system model are developed and presented in dimensional and non-dimensional format. Several important non-dimensional design parameters are introduced. Symbolic and/or explicit formulae of the frequency and mode shape equations are formulated. To the authors knowledge, the present reduced closed form symbolic and explicit frequency equations have not appeared in literature. For different applications, the results are validated using commercial finite element package, namely ANSYS. The beam-column system investigated in this paper is significant for many engineering applications, especially, in mechanical and structural systems.