• Structural Engineering and Mechanics
• /
• v.2 no.1
• /
• pp.95-112
• /
• 1994

This paper describes the formulation and application of a dynamic model for a conventional rail track subjected to arbitary loading functions that simulate wheel/rail impact forces. The rail track is idealized as a periodic elastically coupled beam system resting on a Winkler foundation. Modal parameters of the track structure are first obtained from the natural vibration characteristics of the beam system, which is discretized into a periodic assembly of a specially-constructed track element and a single beam element characterized by their exact dynamic stiffness matrices. An equivalent frequency-dependent spring coefficient representing the resilient, flexural and inertial characteristics of the rail support components is introduced to reduce the degrees of freedom of the track element. The forced vibration equations of motion of the track subjected to a series of loading functions are then formulated by using beam bending theories and are reduced to second order ordinary differential equations through the use of mode summation with non-proportional modal damping. Numerical examples for the dynamic responses of a typical track are presented, and the solutions resulting from different rail/tie beam theories are compared.

• Structural Engineering and Mechanics
• /
• v.24 no.1
• /
• pp.31-50
• /
• 2006

This paper deals with the mechanical behaviour of the thin-walled cylindrical air-spring shell (CAS) made of rubber-textile cord composite (RCC) subjected to different types of loading. An orthotropic hyperelastic constitutive model is presented which can be applied to numerical simulation for the response of biological soft tissue and of the nonlinear anisotropic hyperelastic material of the CAS used in vibroisolation of driver's seat. The parameters of strain energy function of the constitutive model are fitted to the experimental results by the nonlinear least squares method. The deformation of the inflated CAS is calculated by solving the system of five first-order ordinary differential equations with the material constitutive law and proper boundary conditions. Nonlinear hyperelastic constitutive equations of orthotropic composite material are incorporated into the finite strain analysis by finite element method (FEM). The results for the deformation analysis of the inflated CAS made of RCC are given. Numerical results of principal stretches and deformed profiles of the inflated CAS obtained by numerical deformation analysis are compared with experimental ones.

• 제어로봇시스템학회:학술대회논문집
• /
• 1991.10a
• /
• pp.376-382
• /
• 1991

In the analysis of constrained holonomic systems, the Lagange multiplier method yields a system of second-order ordinary differential equations of motion and algebraic constraint equations. Conventional holonomic or nonholonomic constraints are defined as geometric constraints in this paper. Previous works concentrate on the geometric constraints. However, if the total energy of a dynamic system can be computed from the initial energy plus the time integral of the energy input rate due to external or internal forces, then the total energy can be artificially treated as a constraint. The violation of the total energy constraint due to numerical errors can be used as information to control these errors. It is a necessary condition for accurate simulation that both geometric and energy constraints be satisfied. When geometric constraint control is combined with energy constraint control, numerical simulation of a constrained dynamic system becomes more accurate. A new convenient and effective method to implement energy constraint control in numerical simulation is developed based on the geometric interpretation of the relation between constraints in the phase space. Several combinations of energy constraint control with either Baumgarte's Constraint Violation Stabilization Method (CVSM) are also addressed.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.25 no.6
• /
• pp.384-390
• /
• 2015

The vibration analysis of rotating inward beams considering the pre-twisted is presented based on Euler-Bernoulli beam theory. The frequency equations, are calculated using hybrid deformation variable modeling along with the Rayleigh-Ritz assumed mode methods. In this study, resulting system of ordinary differential equations shows the effects of angular speed, and Young's modulus ratio. It is believed that the results will be a reference with which other researchers and commercial FE analysis program, ANSYS can compare their result.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2006.04a
• /
• pp.589-596
• /
• 2006

An improved numerical method to obtain the exact element stiffness matrix is newly proposed to perform the spatially coupled elastic and stability analyses of non-symmetric thin-walled beam-columns with two-types of elastic foundation. This method overcomes drawbacks of the previous method to evaluate the exact stiffness matrix for the spatially coupled stability analysis of thin-walled beam-column. This numerical technique is firstly accomplished via a generalized eigenproblem associated with 14 displacement parameters by transforming equilibrium equations to a set of first order simultaneous ordinary differential equations. Then exact displacement functions are constructed by combining eigensolutions and polynomial solutions corresponding to non-zero and zero eigenvalues, respectively. Consequently an exact stiffness matrix is evaluated by applying the member force-deformation relationships to these displacement functions.

• Journal of the Korean Society of Tobacco Science
• /
• v.20 no.1
• /
• pp.87-94
• /
• 1998

After we made the computer source code with mathematical model of Muramatsu et al. that was expressed by the set of simultaneous first-order ordinary differential equations in evaporation-pyrolysis zone of cigarette, we simulated the distribution profiles of temperature and density of flue-cured tobacco. Those equations were solved numerically with the Runge-Kutta-Gill algorithm assuming step size of 0.025mm by Muramatsu at at,, but in this study the advanced algorithm of Runge-Kutta 4th Order assuming step size of 0.0005mm. The initial conditions and physical parameters of Muramatsu et at. were used for solving them. The calculated values corresponded well with results of Muramatsu et al., especially the gradient of the temperature profile increased with smoldering speed and the thickness of the evaporation-pyrolysis zone decreased with increasing of smoldering speed. On the other hand, the temperature gradient decreased with increasing of the effective thermal-conductivity value and the thickness of the evaporation-pyrolysis zone increased with the effective thermal-conductivity value.

• Journal of applied mathematics & informatics
• /
• v.27 no.5_6
• /
• pp.1001-1015
• /
• 2009

In this paper, two hybrid difference schemes on the Shishkin mesh are constructed for solving a weakly coupled system of two singularly perturbed convection-diffusion second order ordinary differential equations with a small parameter multiplying the highest derivative. We prove that the schemes are almost second order convergence in the supremum norm independent of the diffusion parameter. Error bounds for the numerical solution and its derivative are established. Numerical results are provided to illustrate the theoretical results.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.06a
• /
• pp.183-188
• /
• 2005

The effect of heat conduction resistance on laminar film condensation of the pure saturated vapor in forced flow over a flat plate has been investigated as boundary layer solutions. A efficient numerical methods for water are proposed for its solution. The momentum and energy balance equations are reduced to a nonlinear system of ordinary differential equations with four parameters: the Prandtl number, Pr, Modified Jacob number, $Ja^{\ast}/Pr$, defined by an overall temperature difference, a property ratio $\sqrt{P_l{\mu}_l/P_v{\mu}_v}$ and the conjugate parameter ${\zeta}$. The similarity and simplified solutions obtained reveal the effects of the conjugate parameter.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.8
• /
• pp.742-751
• /
• 2013

The vibration analysis of a rotating cantilever beam made-up of functionally graded materials is presented based on Timoshenko beam theory. The material properties of the beams are assumed to be varied through the thickness direction following a simple power-law form. The frequency equations, which are coupled through gyroscopic coupling terms, are calculated using hybrid deformation variable modeling along with the Rayleigh-Ritz assumed mode methods. In this study, resulting system of ordinary differential equations shows the effects of power-law exponent, angular speed, length to height ratio and Young's modulus ratio. It is believed that the results will be a reference with which other researchers and commercial FE analysis program, ANSYS can compare their results.

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

The Vibration Analysis of Rotating Inward Beams Considering The Tip-Mass is presented based on Euler-Bernoulli beam theory. The frequency equations, which are coupled through gyroscopic coupling terms, are calculated using hybrid deformation variable modeling along with the Rayleigh-Ritz assumed mode methods. In this study, resulting system of ordinary differential equations shows the effects of angular speed, and Young's modulus ratio. It is believed that the results will be a reference with which other researchers and commercial FE analysis program, ANSYS can compare their results.