• Journal of Mechanical Science and Technology
• /
• v.18 no.10
• /
• pp.1747-1754
• /
• 2004

In this study, a new harmonic axisymmetric thick shell element for static and dynamic analyses is proposed. The newly proposed element considering shear strain is based on a modified Hellinger-Reissner variational principle, and introduces additional nodeless degrees for displacement field interpolation in order to enhance numerical performance. The stress parameters selected via the field-consistency concept. are very important in formulating a trouble-free hybrid-mixed elements. For computational efficiency, the stress parameters are eliminated by the stationary condition and then the nodeless degrees are condensed out by the dynamic reduction. Several numerical examples confirm that the present element shows improved efficiency and yields very accurate results for static and vibration analyses.

• Journal of Mechanical Science and Technology
• /
• v.20 no.11
• /
• pp.1950-1960
• /
• 2006

Based on Banach fixed point theorem, a method to calculate nonlinear superposition for three interacting Stokes' waves is proposed in this paper. A mathematical formulation for the nonlinear superposition in deep water and some numerical solutions were investigated. The authors carried out the numerical study with three progressive linear potentials of different wave numbers and succeeded in solving the nonlinear wave profiles of their three wave-interaction, that is, using only linear wave potentials, it was possible to realize the corresponding nonlinear interacting wave profiles through iteration of the method. The stability of the method for the three interacting Stokes' waves was analyzed. The calculation results, together with Fourier transform, revealed that the iteration made it possible to predict higher-order nonlinear frequencies for three Stokes' waves' interaction. The proposed method has a very fast convergence rate.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.4
• /
• pp.559-566
• /
• 2003

In this study, we present a new highly accurate two-dimensional curved composite beam element. The present element, which is based on the Hellinger-Reissner variational principle and classical lamination theory, employs consistent stress parameters corresponding to cubic displacement polynomials with additional nodeless degrees to resolve the numerical difficulties due to the spurious constraints. The stress parameters are eliminated and the nodeless degrees are condensed out to obtain the (9x9) element stiffness matrix. It should be noted that the stacking sequences without transverse deformation to the load plane makes a two dimensional analysis of curved composite beams practically useful . Several numerical examples confirm the superior locking-free behavior of the present higher-order laminated curved beam element.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.1
• /
• pp.28-43
• /
• 2001

Developed in these two series of paper is a complex dynamic model representing the motion of aircraft on the ground and a computer program for numerical simulation. The first part of paper presents the theoretical derivation of equations of motion of the landing gear system based on the physical principle. Developed model is 'structured' in the sense that the undercarriage system is regarded as an assembly of strut, tire, and wheel, where each component is modeled by a separate module. These modules are linked with two external modules-the aircraft and the runway characteristics-to carry out dynamic analysis and numerical simulation of the aircraft motion on the ground. Three sets of coordinate system associated with strut, wheel/tire and runway are defined, and external loads to each component and response characteristics are examined. Lagrangian formulation is used to derive the undercarriage equations of motion relative to the moving aircraft, and the resultant forces and moments from the undercarriage are transformed to aircraft body axes.

• KSTLE International Journal
• /
• v.5 no.1
• /
• pp.23-27
• /
• 2004

The object of the present work is the numerical analysis of the computer hard disk slider. The pressure between slider and disk surfaces is calculated using the Boundary Fitted Coordinate System and Divergence Formulation for the nonlinear Reynolds' equation solution. The optimization scheme is applied to search for the steady state position of the slider. The simplified method is given for the case of the fixed inclined pad. The film thickness ratios and pitching and rolling angles are considered as alternative choice of the slider's coordinates. The behavior of the objective function for the Negative Pressure slider is studied in details. Methods of conjugate directions and feasible directions are applied.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2667-2677
• /
• 1995

A new numerical algorithm using equal order linear finite element and fractional step method has been developed which is capable of analyzing unsteady fluid flow and heat transfer problems. Streamline Upwind Petrov-Galerkin (SUPG) method is used for the weighted residual formulation of the Navier-Stokes equations. It is shown that fractional step method, in which pressure term is splitted from the momentum equation, reduces computer memory and computing time. In addition, since pressure equation is derived without any approximation procedure unlike in the previously developed SIMPLE algorithm based FEM codes, the present numerical algorithm gives more accurate results than them. The present algorithm has been applied preferentially to the well known bench mark problems associated with steady flow and heat transfer, and proves to be more efficient and accurate.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2006.05a
• /
• pp.1348-1353
• /
• 2006

This paper is concerned with the comparison of component mode synthesis md substructure synthesis methods. Those two methods were developed for the dynamic analysis of complex structure. However, the difference between two methods in solving the dynamic problem has never been explained. In this paper, a structure consisting of two beams is considered for the application of two methods. The dynamic formulation shows the characteristics and differences of two methods. The procedure for choosing proper substructure modes in each method is explained in detail. In addition, the advantage of the substructure synthesis method is discussed based on the numerical example. Numerical examples show how two methods are applied to the addressed problem.

• Structural Engineering and Mechanics
• /
• v.35 no.6
• /
• pp.717-733
• /
• 2010

The dynamic stiffness matrix is formulated for an axially loaded slender double-beam element in which both beams are homogeneous, prismatic and of the same length by directly solving the governing differential equations of motion of the double-beam element. The Bernoulli-Euler beam theory is used to define the dynamic behaviors of the beams and the effects of the mass of springs and axial force are taken into account in the formulation. The dynamic stiffness method is used for calculation of the exact natural frequencies and mode shapes of the double-beam systems. Numerical results are given for a particular example of axially loaded double-beam system under a variety of boundary conditions, and the exact numerical solutions are shown for the natural frequencies and normal mode shapes. The effects of the axial force and boundary conditions are extensively discussed.

• Nuclear Engineering and Technology
• /
• v.33 no.1
• /
• pp.62-72
• /
• 2001

The one-dimensional finite element program was developed to analyze the coupled behavior of heat, moisture, and air transfer in unsaturated porous media. By using this program, the simulation results were compared with those from the laboratory infiltration tests under isothermal condition and temperature gradient condition, respectively. The discrepancy of water uptake was found in the upper region of a bentonite sample under isothermal condition between numerical simulation and laboratory experiment. This indicated that air pressure was built up in the bentonite sample which could retard the infiltration velocity of liquid. In order to consider the swelling phenomena of compacted bentonite which cause the discrepancy of the distribution of water content and temperature, swelling and shrinkage factors were incorporated into the finite element formulation. It was found that these factors could be effective to represent the moisture diffusivity and unsaturated hydraulic conductivity due to volume change of bentonite sample.

• Proceedings of the Korean Nuclear Society Conference
• /
• 1997.10a
• /
• pp.617-624
• /
• 1997

This Paper Presents an effective and simple procedure for the simulation of the motion of the solid-liquid interfacial boundary and the transient temperature field during phase change process. To accomplish this purpose, an iterative implicit solution algorithm has been developed by employing the dual reciprocity boundary element method. The dual reciprocity boundary element approach provided in this paper is much simpler than the usual boundary element method applying a reciprocity principle and an available technique for dealing with domain integral of boundary element formulation simultaneously. The effectiveness of the present analysis method have been illustrated through comparisons of the calculation results of an example with its semi-analytical or other numerical solutions where available.