• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.18 no.2
• /
• pp.141-149
• /
• 2005

This paper deals with formulations of the energy equilibrium equation by an introduction of the algebraic description, quarternion, which meets conservations of system energy for the equation of motion. Then the equation is discretized to analyze the dynamits analysis of flexible multibody systems in such a way that the work done by the constrained force completely is eliminated. Meanwhile, Rodrigues parameters we used to express the finite rotation lot the proposed method. This method lot the initial essential step to a guarantee of developments of the 3D dynamical problem provides unconditionally stable conditions for the nonlinear problems through the numerical examples.

• Journal of the Korean Society of Tobacco Science
• /
• v.2 no.2
• /
• pp.44-50
• /
• 1980

The growth of flue-cured tobacco was analyzed with the mathematical treatment. The results are summarized as followings: 1. The growth curve was fitted to the quadratic polynomial equation in improved-mulching cultivation, blot to the Gompertz equation in the other cultivations. 2. The initial point of the maximum growth phase for dry weight was about 50 days after transplanting in improved-mulching cultivation, but about 40 days in the other cultivations, and the maximum growth period was for 25 days in all cultivations. 3. The growth rate of the maximum growth period in dry weight decreased in the order of improved-mulching cultivation, mulching cultivation, and non-mulching cultivation. 4. A relative growth amount in the maximum growth period was higher in later sowing. 5. The length of maximum growth was 5 days shorter in leaf area than in dry weight. The maximum growth phase was 7 days earlier in leaf area than in dry weight.

• Water for future
• /
• v.18 no.1
• /
• pp.85-94
• /
• 1985

A method is developed for the separation of the major base flow in a river hydrograph combining the numerical techniques and the empirical methods. The linearized Boussinesq equation and the storage function are used to obtain the base flow recession. The shape of base flow curve made by the recharge of the groundwater table aquifer resulting from rainfall in determined by the Singh and Stall's graphical method, and the continuous from for the curve is approximated by the multiple and polynomial regression. this procedure was successfully tested for the separation of base flow and the establishment of hydrograph in a natural watershed. It was found that the direct numerical method applied to the homogeneous linear second order ordinary differential equation system is not suited to obtain the recession curve, and the case that the loss is generated in the partially penetrating stream can not be solved by the method of this study.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.14 no.10
• /
• pp.4733-4738
• /
• 2013

This paper proposes a method to obtain the dynamic equilibrium configuration of a constrained mechanical system by using multibody dynamic analysis. Dynamic equilibrium equations with independent coordinates are derived from the time-dependent constraint equations and dynamic equations of a multibody system. The Newton-Raphson method is used to find numerical solutions for nonlinear algebraic equations that are composed of the dynamic equilibrium and constraint equations. The proposed method is applied to obtain the dynamic equilibrium configuration of a speed governor, and the results are verified on the basis of the results from conventional dynamic analysis. Furthermore, vertical displacements at equilibrium configuration, which varied with the rotational velocity of the speed governor, are calculated, and design parameter analysis of the equilibrium configuration is presented.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.3
• /
• pp.45-52
• /
• 2003

Detailed flowfield resulting from the secondary sonic gas injection into a divergent section of supersonic conical nozzle has been numerically investigated. The three-dimensional flowfield associated with the bow-shock/boundary-layer interaction inside the nozzle has been solved by Reynolds-averaged Navier-Stokes equations with an algebraic and $\kappa$-$\varepsilon$ turbulence model. The numerical results have been compared with the experimental results for the identical flow conditions, and it is shown that the comparison is satisfactory Effects of different injection pressures of the secondary jet on the shock/boundary-layer interactions and the overall flow structure inside the nozzle have been investigated. The vortex structures behind the shock interaction and wall pressure variations have also been studied.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.2
• /
• pp.86-97
• /
• 1993

A numerical method is developed for the nonlinear motion of two-dimensional wedges and axisymmetric-forced-heaving motion using Semi-Largrangian scheme under assumption of potential flows. In two-dimensional-problem Cauchy's integral theorem is applied to calculate the complex potential and its time derivative along boundary. In three-dimensional-problem Rankine ring sources are used in a Green's theorem boundary integral formulation to salve the field equation. The solution is stepped forward numerically in time by integrating the exact kinematic and dynamic free-surface boundary condition. Numerical computations are made for the entry of a wedge with a constant velocity and for the forced harmonic heaving motion from rest. The problem of the entry of wedge compared with the calculated results of Champan[4] and Kim[11]. By Fourier transform of forces in time domain, added mass coefficient, damping coefficient, second harmonic forces are obtained and compared with Yamashita's experiment[5].

• Journal of the Korean Society of Propulsion Engineers
• /
• v.8 no.4
• /
• pp.9-15
• /
• 2004

The measured thrust could be distorted because of the transient response of thrust stand during combustion of rocket motor. As a result of the distorted thrust, it is not easy for us to know the values of thrust peak and thrust duration time. Therefore, it is of great importance to compute the true thrust from the measured thrust. In this study the method to eliminate the transient response from the measured thrust using only the measured thrust was Proposed, and also experimental data were used to approve the proposed method. The result showed that the proposed method would be available to compute the true thrust.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.1
• /
• pp.1-7
• /
• 2002

In the past, safety assessment on the site of Ulchin Nuclear Power Plants against tsunamis was carried out with probable maximum earthquake magnitude and related tsunamigenic fault parameters. Recently, however, based on the seismic gap theory, some seismologists warned about earthquakes of larger magnitudes than had been expected. In this study, we revaluated tsunami risk with a finite difference model based on linear and nonlinear shallow water equations. Firstly, we simulated the\`83 tsunami and compared the calculated water surface profile with the observed wave heights. Secondly, we evaluated the rise and drop of sea water level at the site of Ulchin Nuclear Power Plant with fault parameters of the past '83, '93 tsunamis and some dangerous faults. Finally, we showed that the cooling water intake facility of Ulchin Nuclear Power Plants would be safely operated in disastrous tsunamis.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.23 no.11
• /
• pp.943-950
• /
• 2013

This paper discusses a harmonic response estimation method on the L$\acute{e}$vy plate with two opposite edges simply supported and the other two edges having free boundary conditions. Since the equation of motion of the plate is not self-adjoint, the modes are not orthogonal to each other on the domain. Noting that the L$\acute{e}$vy plate can be expressed using one term sinusoidal function that is orthogonal to other sinusoidal functions, this paper suggested the calculation method that is equivalent to finding a least square error minimization solution of the finite number of algebraic equations. Example problems subjected to a distributed area loading and a distributed line loading are defined and their solutions are provided. The solutions are compared to those of the commercial code, ANSYS. According to the verification results, it is expected that the suggested method will be useful to predict the forced response on the L$\acute{e}$vy plate with the distributed area or line loading conditions.

• Journal of the Korean Society for Precision Engineering
• /
• v.17 no.11
• /
• pp.176-184
• /
• 2000

Recently, it becomes a very important issue to consider the mechanical systems such as high-speed vehicle and railway train moving on a flexible beam structure. Using general approach proposed in the first part of this paper, it tis possible to predict planar motion of constrained mechanical system and elastic structure with various kinds of foundation supporting condition. Combined differential-algebraic equations of motion derived from both multibody dynamics theory and Finite Element Method can be analyzed numerically using generalized coordinate partitioning algorithm. To verify the validity of this approach, results from simply supported elastic beam subjected to a moving load are compared with exact solution from a reference. Finally, parameter study is conducted for a moving vehicle model on a simply supported 3-span bridge.