• 한국기계연구소 소보
• /
• s.18
• /
• pp.55-66
• /
• 1988

Many direct integration methods are used for numerical analyses of dynamic motion. In these methods, the governing equations of a dynamic system are integrated successively using a step-by-step numerical integration procedure. Time derivatives in the equations are generally approximated using difference formulas involving one or more increments of the time. Time increment has closely relationship with the accuracy of the motion analysis. In this paper, a 4th order Houbolt direct integration method is derived. For a spring-mass system, the motion of the system are analyzed from the 3rd order Houbolt and the 4th order Houbolt approaches respectively. Finally the paper proposes the optimal time-increment based on the accuracy of numerical analyses.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.6
• /
• pp.1014-1019
• /
• 2003

A dynamic modeling method for beams undergoing overall rigid body motion is presented in this paper. Two special deformation variables are introduced to represent the stretching and the curvature and are approximated by the assumed mode method. Geometric constraint equations that relate the two special deformation variables and the cartesian deformation variables are incorporated into the modeling method. By using the special deformation variables, all natural as well as geometric boundary conditions can be satisfied. It is shown that the geometric nonlinear effects of stretching and curvature play important roles to accurately predict the dynamic response when overall rigid body motion is involved.

• 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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.5
• /
• pp.904-913
• /
• 2002

In this paper, the dynamic stress intensity factors of fracture mechanics are numerically computed in time domain using the FEM. For which the finite element formulations are derived applying the Galerkin method to the equations of motion in convolution integral as has been presented in the previous paper. To assure the strain fields of r$^{-1}$ 2/ singularity near the crack tip, the triangular quarter-point singular elements are imbedded in the finite element mesh discretized by the isoparametric quadratic quadrilateral elements. Two-dimensional problems of the elastodynamic fracture mechanics under the impact load are solved and compared with the existing numerical and analytical solutions, being shown that numerical results of good accuracy are obtained by the presented 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.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.05a
• /
• pp.877-882
• /
• 2005

Dynamic stability of an axially accelerating beam stucture 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 targe 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.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.5 s.176
• /
• pp.1093-1102
• /
• 2000

Dynamic behaviors are analyzed for an automatic ball balancer with double races which is a device to reduce eccentricity of rotors. Equations of motion are derived by using the polar coordinate sys tem instead of the rectangular coordinate system which is used in other previous researches. To analyze the stability around equilibrium positions, the perturbation method is used. On the other hand, the time responses are computed from the nonlinear equations of motion by using a time integration method.

• 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.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.11a
• /
• pp.563-567
• /
• 2000

This research is concerned with the dynamic modeling of solar arrays equipped with strain energy hinges(SEH). It is found from experiments that the SEH has nonlinear dynamic characteristics and complex buckling behavior, which is difficult to explain theoretically. In this paper, we use an equivalent one-dimensional nonlinear torsional spring for the SEH. Assuming that solar panels are rigid, we developed the systematic approach for the derivation of the theoretical model for the solar arrays equipped with the multitudes of the SEH. To this end, the kinematic relation of the displacement vector of each body is derived and then applied to the equations of motion. Lagrangian equations of motion are used for the derivations.