• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.9
• /
• pp.1327-1334
• /
• 2003

Flow characteristics of the fluid laden with many particles in the two-dimensional channel are investigated using the Navier-Stokes equations coupled with the equation of motion of particles by direct numerical simulation. A four-step fractional step method with Crank-Nicolson scheme and ALE technique is used for P2P1 mixed finite element method. The motion and distribution of particles in the fluid is virtually described as a result of direct numerical simulation and the increase of viscosity is compared with theoretical equations. The effect of channel height on the relative viscosity and the tubular pinch effect are discussed.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2011.10a
• /
• pp.130-135
• /
• 2011

A study of the bending-extension-transverse shear coupled random response of the composite beams with thin-walled open sections subjected to various types of concentrated and distributed random excitations is dealt with in this paper. First of all, equations of motion of thin-walled composite H-type cross-section beams incorporating a number of nonclassical effects of transverse shear and primary and secondary warping, and anisotropy of constituent materials are derived. On the basis of derived equations of motion, analytical expressions for the displacement response of the composite beams are derived by using normal mode method combined with frequency response function method.

• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.6
• /
• pp.148-159
• /
• 1999

A numerical method is presented for the dynamic analysis of spur gears rotating with very high angular speeds. For an efficient computation each gear is assumed to consist of a rotating rigid disk and an elastic tooth having mass, and finite element formulations are used for the equations of motion of the tooth. The geometric constraint is imposed between the rigid disk and the elastic tooth to fix them, and contact condition is imposed between the meshing teeth of the gears. At each iteration of each time step the Lagrange multiplier and contact force are revised by using the constraint error vector, and then the whole equations of motion are time integrated with the given Lagrange multiplier and contact force. For the accurate solution the velocity and acceleration constraints as well as the displacement constraint are satisfied by the monotone reductions of the constraint error vectors. Computing procedures associated with the iterative schemes are explained and numerical simulations are conducted with the spur gears.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.16 no.2 s.107
• /
• pp.198-206
• /
• 2006

The dynamic responses of both the in-plane and out-of-plane vibrations are investigated for an axially moving membrane. The equations of motion are derived for the moving membrane with no-slip boundary conditions by using the extended Hamilton principle. Based on the Galerkin method, the discretized equations of motion are derived. The generalized-time integration method is applied to compute the dynamic responses for the in-plane and out-of-plane motions. From the computed results, the responses are compared between the in-plane and out-of-plane vibrations. Furthermore. the effects of velocity and acceleration on the dynamic behaviours for displacements and stresses are presented.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2014.10a
• /
• pp.518-522
• /
• 2014

This paper is concerned with the dynamic modeling of transmission line undergoing galloping vibrations. To this end, the kinetic and potential energies of a uniform wire vibrating in space are derived. The equations of motion suitable for numerical simulations are derived using the assumed mode method and Lagrange equation. The resulting equations of motion are expressed in matrix form. To cope with bundled transmission line, the spacer was modelled by a spring element. As a numerical example, a two-wire transmission line combined by spacers was considered. Natural vibration characteristics show that the in-plane vibrations of the transmission line appeared in low frequency range, which may lead to galloping.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.9 s.102
• /
• pp.1053-1059
• /
• 2005

Dynamic stability of an axially accelerating beam structure is investigated in this paper. The equations of motion of a fixed-free beam are derived using the hybrid deformation variable method and the assumed mode method. Unstable regions due to periodical acceleration are obtained by using the Floquet's theory. Stability diagrams are presented to illustrate the influence of the dimensionless acceleration, amplitude, and frequency. Also, buckling occurs when the acceleration exceeds a certain value. It is found that relatively large unstable regions exist around the first bending natural frequency, twice the first bending natural frequency, and twice the second bending natural frequency. The validity of the stability diagram is confirmed by direct numerical integration of the equations of motion.

• Structural Engineering and Mechanics
• /
• v.13 no.5
• /
• pp.569-589
• /
• 2002

In performing the dynamic analysis, the step size used in a step-by-step integration method might be much smaller than that required by the accuracy consideration in order to capture the rapid chances of dynamic loading or to eliminate the linearization errors. It was first found by Chen and Robinson that these difficulties might be overcome by integrating the equations of motion with respect to time once. A further study of this technique is conducted herein. This include the theoretical evaluation and comparison of the capability to capture the rapid changes of dynamic loading if using the constant average acceleration method and its integral form and the exploration of the superiority of the time integration to reduce the linearization error. In addition, its advantage in the solution of the impact problems or the wave propagation problems is also numerically demonstrated. It seems that this time integration technique can be applicable to all the currently available direct integration methods.

• Journal of the Optical Society of Korea
• /
• v.8 no.3
• /
• pp.137-146
• /
• 2004

A new method of obtaining the eigenmode of an anisotropic planar waveguide is studied. The planar waveguide can be composed of an arbitrary number of isotropic or uniaxially anisotropic layers, provided all the optical axes arc lying in the incidence plane. Since the equation of motion for the TE mode is not different from that for the TE mode in an isotropic planar waveguide, only the equation of motion for the TM mode is of any concern. For this kind of device structure, the Maxwell's equations can be solved for one component of the electric field and one component of the magnetic field. The resulting coupled set of equations is linear in the propagation constant and the eigenmode can be easily obtained using canned numerical routines.

• Journal of Ocean Engineering and Technology
• /
• v.13 no.2 s.32
• /
• pp.129-137
• /
• 1999

The time of the onset of double-diffusive convection in time-dependent, nonlinear concentration fields is investigated theoretically. The initially quiescent horizontal fluid layer with a uniform temperature gradient experiences a sudden concentration change from below, but its stable thermal stratification affects concentration effects in such way to invoke convective motion. The related stability analysis, including Soret effect, is conducted on the basis of the propagation theory. Under the linear stability theory the concentration penetration depth is used as a length scaling factor, and the similarity transform for the linearized perturbation equations. The newlly obtained stability equations are solved numerically. The resulting critical time to mark the onset of regular cells are obtained as a function of the thermal Rayleigh number, the solute Rayleigh number, and the Soret effect coefficient. For a certain value of the Soret effect coefficient, the stable thermal gradient promote double-diffusive convective motion.

• Proceedings of the KIEE Conference
• /
• 2000.07d
• /
• pp.2804-2807
• /
• 2000

This paper suggests a method of the forward dynamic analysis for the computer simulation on the analysis of the dynamic behavior for biped walking robot. The equations f motion of the system or the simulation are constructed by using the Method of the multibody dynamics which is powerful method for modeling of the complex biped system. For the simplicity of simulation, we consider that the sole of the contacting foot is affected by the reaction forces for tree structure system topology instead of the addition or deletion of the kinematic constraints. The ground reaction forces can be modeled using the simple spring and damper model at the three contacting points on the sole of the foot. For minimizing the errors of numerical integration, the number of equations of motion is minimized by adding the driving constraints or a controller instead of the direct driving torques.