• Tribology and Lubricants
• /
• v.38 no.2
• /
• pp.33-40
• /
• 2022

In this study, we present rotor dynamic analysis and operation test of a turbo expander for a hydrogen liquefaction plant. The turbo expander consists of a turbine and compressor wheel connected to a shaft supported by two hydrostatic radial and thrust bearings. In rotor dynamic analysis, the shaft is modeled as a rigid body, and the equations of motion for the shaft are solved using the unsteady Reynolds equation. Additionally, the operating test of the turbo expander has been performed in the test rig. Pressurized helium is supplied to the bearings at 8.5 bar. Furthermore, we monitor the shaft vibration and flow rate of the helium supplied to the bearings. The rotor dynamic analysis result shows that there are two critical speeds related with the rigid body mode under 40,000 rpm. At the first critical speed of 36,000 rpm, the vibration at the compressor side is maximum, whereas that of the turbine is maximum at the second critical speed of 40,000 rpm. The predicted maximum shaft vibration is 3 ㎛, whereas sub-synchronous vibration is not presented. The operation test results show that there are two critical speeds under the rated speed, and the measured vibration value agrees well with predicted value. The measured flow rate of the helium supplied to the bearing is 2.0 g/s, which also agrees well with the predicted data.

• Journal of Navigation and Port Research
• /
• v.34 no.3
• /
• pp.195-203
• /
• 2010

This paper concerns about total dynamic analysis of integrated mining system. This system consists of vertical steel pipe, intermediate buffer station, flexible pipe and self-propelled miner. The self-propelled miner and buffer are assumed as rigid-body of 6-dof. Discrete models of vertical steel pipe and flexible pipe are adopted, which are obtained by means of lumped-parameter method. The motion of mining vessel is not considered. Instead, the motion of mining vessel is taken into account in form of various boundary conditions (e.g. forced excitation in slow motion and/or fast oscillation and so on). A terramechanics model of extremely cohesive soft soil is applied to the self-propelled miner. Hinged and ball constraints are used to define the connections between sub-systems (vertical steel pipe, buffer, flexible pipe, self-propelled miner). Equations of motion of the coupled model are derived with respect to the each local coordinates system. Four Euler parameters are used to express the orientations of the sub-systems. To solve the equations of motion of the total dynamic model, an incremental-iterative formulation is employed. Newmark-${\beta}$ method is used for time-domain integration. The total dynamic responses of integrated mining system are investigated.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.3
• /
• pp.421-429
• /
• 2010

In this paper, in order to analyze the dynamic response of a floating crane when it lifts a heavy cargo, the boom of the floating crane is considered as an elastic beam. The boom is divided into elements based on finite element formulation and the floating frame of reference formulation and nodal coordinates are employed to model the boom as a flexible body. As an extension of the previous study, in order to consider spatial motion in waves, the coupled equations of motions of the 6 degree of freedom (DOF) floating crane and 6 DOF cargo are developed based on the flexible multibody system dynamics. The 3 dimensional deformation of the elastic boom is considered with 18 DOF. The dynamic simulation of the floating crane and the cargo is performed under regular wave conditions with various cargo weights. Finally, the effects of the elastic boom on lifting cargo are discussed by comparing the simulation results between the elastic boom and a rigid boom.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.10 no.3
• /
• pp.385-392
• /
• 2018

This paper presents 2D numerical modeling to calculate ship-ice interactions that occur when an icebreaker advances into ice-covered water. The numerical model calculates repeated icebreaking of an ice plate and removal of small ice floes. The icebreaking of the ice plate is calculated using a ship-ice contact detection technique and fluid-structural interaction of ice plate bending behavior. The ship-ice interactions in small ice floes are calculated using a physically based modeling with 3DOF rigid body equations. The ice plate is broken in crushing, bending, and splitting mode. The ice floes drift by wind or current and by the force induced by the ship-ice interaction. The time history of ice force and ice floe distribution when an icebreaker advances into the ice-covered water are obtained numerically. Numerical results demonstrate that the time history of ice force and distribution of ice floes (ice channel width) depend on the ice floe size, ship motion and ice drifting by wind or current. It is shown that the numerical model of ship maneuvering in realistic ice conditions is necessary to obtain precise information about the ship in ice-covered water. The proposed numerical model can be useful to provide data of a ship operating in ice-covered water.

• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.2
• /
• pp.137-144
• /
• 2006

The objective of this paper is to present trajectory guidance and control system with a dynamic inversion for a small unmanned aerial vehicle (UAV). The UAV model is expressed by fixed-mass rigid-body six-degree-of-freedom equations of motion, which include the detailed aerodynamic coefficients, the engine model and the actuator models that have lags and limits. A trajectory is generated from the given waypoints using cubic spline functions of a flight distance. The commanded values of an angle of attack, a sideslip angle, a bank angle and a thrust, are calculated from guidance forces to trace the flight trajectory. To adapt various waypoint locations, a proportional navigation is combined with the guidance system. By the decision logic, appropriate guidance law is selected. The flight control system to achieve the commands is designed using a dynamic inversion approach. For a dynamic inversion controller we use the two-timescale assumption that separates the fast dynamics, involving the angular rates of the aircraft, from the slow dynamics, which include angle of attack, sideslip angle, and bank angle. Some numerical simulations are conducted to see the performance of the proposed guidance and control system.

• Journal of the Korean Society for Aviation and Aeronautics
• /
• v.18 no.4
• /
• pp.30-37
• /
• 2010

This paper deals with attitude determination and parameter estimation problems for a lunar module. For this we first derive equations of motion for the lunar module by considering allocation locations (configurations) of reaction thruster and a reaction wheel assembly. The lunar module is assumed as a rigid body. In order to include the effect of fuel sloshing on the dynamics of the lunar module, we model it as a spherical pendulum for a simple analysis. For estimating angular rates and moment of inertia of the module, we employ an extended Kalman filter and the least mean square algorithms, respectively. Finally we construct a dynamical model for the lunar module by combining all these elements.

• The KSFM Journal of Fluid Machinery
• /
• v.18 no.6
• /
• pp.5-11
• /
• 2015

In a pressurized water reactor (PWR), control rod assembly (CRA) falls into the guide tubes of a fuel assembly due to gravity for scram. Various theoretical approaches and numerical analyses have been performed because its shape is simple and its design was completely developed several decades ago. A control rod assembly for a sodium-cooled faster reactor (SFR) which is geometrically more complicated is being actively developed in Korea nowadays. Drop time and impact velocity of a CRA are important parameters with respect to reactivity insertion time and the mechanical robustness of a CRA and a guide duct. In this paper, computational method considering simultaneously the equation of motion for rigid body and the Navier-Stokes equations for fluid is suggested and verified by comparison with theoretical analysis results. Through this valuable CFD analysis method, drop time and impact velocity of initially designed SFR CRA are evaluated before performing scram tests with it.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.14 no.2
• /
• pp.150-156
• /
• 2004

In this paper, the dynamic behavior of a very flexible cable due to moving multibody system along its length is presented. The very deformable motion of a cable is presented using absolute nodal coordinate formulation, which is based on the finite element procedures and the general continuum mechanics theory to represent the elastic forces. Formulation for the sliding joint between a very flexible beam and a rigid body is derived. In order to formulate the constraint equations of this joint, a non-generalized coordinate, which has no inertia or forces associated with this coordinate, is used. The modeling of this sliding joint is very important to many mechanical applications such as the ski lifts. cable cars, and pulley systems. A multibody system moves along an elastic cable using this sliding joint. A numerical example is shownusing the developed analysis program for flexible multibody systems that include a large deformable cable.

• Wind and Structures
• /
• v.20 no.3
• /
• pp.423-448
• /
• 2015

In this paper the unsteady fluid-structure interaction (FSI) problems with large structural displacement are solved by partitioned solution approaches in the arbitrary Lagrangian-Eulerian finite element framework. The incompressible Navier-Stokes equations are solved by the characteristic-based split (CBS) scheme. Both a rigid body and a geometrically nonlinear solid are considered as the structural models. The latter is solved by Newton-Raphson procedure. The equation governing the structural motion is advanced by Newmark-${\beta}$ method in time. The dynamic mesh is updated by using moving submesh approach that cooperates with the ortho-semi-torsional spring analogy method. A mass source term (MST) is introduced into the CBS scheme to satisfy geometric conservation law. Three partitioned coupling strategies are developed to take FSI into account, involving the explicit, implicit and semi-implicit schemes. The semi-implicit scheme is a mixture of the explicit and implicit coupling schemes due to the fluid projection splitting. In this scheme MST is renewed for interfacial elements. Fixed-point algorithm with Aitken's ${\Delta}^2$ method is carried out to couple different solvers within the implicit and semi-implicit schemes. Flow-induced vibrations of a bridge deck and a flexible cantilever behind an obstacle are analyzed to test the performance of the proposed methods. The overall numerical results agree well with the existing data, demonstrating the validity and applicability of the present approaches.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.3 s.96
• /
• pp.251-258
• /
• 2005

The free vibration of a spinning flexible disk-spindle system supported by hydro dynamic bearings (HDB) in an HDD is analyzed by FEM. The spinning flexible disk is described using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. It is discretized by annular sector element. The rotating spindle which includes the clamp, hub, permanent magnet and yoke, is modeled by Timoshenko beam including the gyroscopic effect. The flexible supporting structure with a complex shape which includes stator core, housing, base plate, sleeve and thrust pad is modeled by using a 4-node tetrahedron element with rotational degrees of freedom to satisfy the geometric compatibility. The dynamic coefficients of HDB are calculated from the HDB analysis program, which solves the perturbed Reynolds equation using FEM. Introducing the virtual nodes and the rigid link constraints defined in the center of HDB, beam elements of the shaft are connected to the solid elements of the sleeve and thrust pad through the spring and damper element. The global matrix equation obtained by assembling the finite element equations of each substructure is transformed to the state-space matrix-vector equation, and the associated eigen value problem is solved by using the restarted Arnoldi iteration method. The validity of this research is verified by comparing the numerical results of the natural frequencies with the experimental ones. Also the effect of supporting structures to the natural modes of the total HDD system is rigorously analyzed.