• Proceedings of the Korean Society of Tribologists and Lubrication Engineers Conference
• /
• 2004.11a
• /
• pp.143-148
• /
• 2004

The numerical analysis of the hard disk drive slider is presented. The pressure distribution was calculated using the finite element method. The generalized Reynolds equation was applied in order to include the gas rarefaction effect. The balance of the air bearing force and preload force was considered. The characteristics of the small vibrations near the equilibrium were studied using the perturbation method. Triangular mesh with variable element size was employed to model the two-rail slider. The flying height, pitching angle, rolling angle, stiffness and damping of the two-rail slider were calculated for radial position changing from the inner radius to the outer radius and for a wide range of the slider crown values. It was found that the flying height, pitching angle and rolling angle were increased with radial position while the stiffness and damping coefficients were decreased. The higher values of crown resulted in increased flying height, pitching angle and damping and decreased stiffness.

• Structural Engineering and Mechanics
• /
• v.38 no.4
• /
• pp.529-547
• /
• 2011

This paper presents the global-local finite cover method (GL-FCM) that is capable of analyzing structures involving local heterogeneities and propagating cracks. The suggested method is composed of two techniques. One of them is the FCM, which is one of the PU-based generalized finite element methods, for the analysis of local cohesive crack growth. The mechanical behavior evaluated in local heterogeneous structures by the FCM is transferred to the overall (global) structure by the so-called mortar method. The other is a method of mesh superposition for hierarchical modeling, which enables us to evaluate the average stiffness by the analysis of local heterogeneous structures not subjected to crack propagation. Several numerical experiments are conducted to validate the accuracy of the proposed method. The capability and applicability of the proposed method is demonstrated in an illustrative numerical example, in which we predict the mechanical deterioration of a reinforced concrete (RC) structure, whose local regions are subjected to propagating cracks induced by reinforcement corrosion.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.43-48
• /
• 1998

The development of unsteady three-dimensional incompressible viscous solver based on unsteady physical curvilinear coordinate system is presented. A 12-point finite analytic scheme based on local uniform grid spacing is extended for nonuniform grid spacing. The formulation of a condition is suggested to avoid the oscillation of the series summations produced by the application of the method of separation of variables. SIMPLER and pressure Poisson equation techniques are used for solving a velocity-pressure coupled problem. The matrix is solved using the Generalized Minimal RESidual (GMRES) method to enhance the convergence rate of unsteady flow solver and the Kinematic boundary condition of a free surface flow. It is demonstrated that the numerical solutions of these equations are less mesh sensitive.

• Interaction and multiscale mechanics
• /
• v.3 no.2
• /
• pp.111-121
• /
• 2010

Dynamic analysis of nano/micro bio-sensors based on a multiscale atomistic/continuum theory is introduced. We use a generalized atomistic finite element method (GAFEM) to analyze a bio-sensor which has $3{\times}N_a{\times}N_p$ degrees of freedom, where $N_p$ is the number of representative unit cells and $N_a$ is the number of atoms per unit cell. The stiffness matrix is derived from interatomic potential between pairs of atoms. This work contains two studies: (1) the resonance analysis of nano bio-sensors with different amount of target analyte and (2) the dependence of resonance frequency on finite element mesh. We also examine the Courant-Friedrichs-Lewy (CFL) condition based on the highest resonance frequency. The CFL condition is the criterion for the time step used in the dynamic analysis by GAFEM. Our studies can be utilized to predict the performance of micro/nano bio-sensors from atomistic perspective.

• Bulletin of the Korean Chemical Society
• /
• v.12 no.2
• /
• pp.228-238
• /
• 1991

Generalized multichannel quantum defect theory [C. H. Greene et al. Phys. Rev., A26, 2441 (1982)] is implemented to the vibrational predissociation of triatomic van der Waals molecules. As this is the first one of such an application, the dependences of the quantum defect parameters on energy and radius are examined carefully. Calculation shows that, in the physically important region, quantum defect parameters remain smoothly varying functions of energy for this system as in atomic applications, thus allowing us very coarse energy mesh calculations for the photodissociation spectra. The choice of adiabatic or diabatic potentials as reference potentials for the calculation of quantum defect parameters as done by Mies and Julienne [J. Chem. Phys., 80, 2526 (1984)] can not be used for this system. Physically motivated reference potentials that may be generally applicable to all kinds of systems are utilized instead. In principle, implementation can be done to any other predissociation processes with the same method.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.44 no.1
• /
• pp.49-62
• /
• 2024

Reckless development of the underground by rapid urbanization causes inspection delay on replacement of existing structure and installation new facilities. However, frequent accidents occur due to deviation in construction design planned by inaccurate location information of underground structure. Meanwhile, the electrical resistivity survey, knowns as non-destructive method, is based on the difference in the electric potential of electrodes to measure the electrical resistance of ground. This method is significantly advanced with multi-electrode and deep learning for analyzing strata. However, there is no study to quantitatively assess change in electrical resistance according to geometric conditions of structures. This study evaluates changes in electrical resistance through geometric parameters of electrodes and structure. Firstly, electrical resistance numerical module is developed using generalized mesh occurring minimal errors between theoretical and numerical resistance values. Then, changes in resistances are quantitatively compared on geometric parameters including burial depth, diameter of structure, and distance electrode and structure under steady current condition. The results show that higher electrical resistance is measured for shallow depth, larger size, and proximity to the electrode. Additionally, electric potential and current density distributions are analyzed to discuss the measured electrical resistance around the terminal electrode and structure.

• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.11
• /
• pp.150-157
• /
• 2003

In this paper, the large deformation of hollow silicon rubber gasket is treated. The frictional contact occurs between groove and the outer part of hollow gasket, and the frictional self-contact exists in the inner parts of hollow gasket. The silicon rubber has the nonlinear elastic behavior and its material property is approximately incompressible. Hence, the stress analysis requires an existence of a strain energy function, which is usually defined in terms of invariants or stretch ratio such as generalized Mooney-Rivlin and Ogden model. Considering large compressive deformation and friction, Mooney-Rivlin 3rd model and Coulomb's friction model are assumed. The numerical analysis is obtained by the commercial finite element program MARC. But, due to large deformation, the elements degenerate in the inner parts of hollow gasket. This means that the analysis of subsequent increments is carried out with a very poor mesh. In order to continue the analysis with a sufficient accuracy, it is necessary to use new finite element modeling by remesh. Experiments are also performed to show the validity of present method. As a conclusion, numerical results by this research have good agreements with experiments.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.23 no.7 s.166
• /
• pp.1205-1214
• /
• 1999

A new quadrilateral plane stress element is developed for efficient and accurate analysis of thermal stress problems. It is convenient to use the same mesh and the same shape functions for thermal analysis and stress analysis. But, because of the inconsistency between deformation related strain field and thermal strain field, oscillatory responses and considerable errors in stresses are resulted in. To avoid undesired oscillations, strain approximation is enhanced by supplementing several assumed strain terms based on the variational principle. Thermal deformation is incorporated into the generalized mixed variational principle for displacement, strain and stress fields, and basic equations for the modified enhanced assumed strain method are derived. For the stress approximation of bilinear elements, the $5{\beta}$ version of Pian and Sumihara is adopted. The numerical results for several problems show that the present element behaves well and reduces oscillatory responses. it also results in almost the same magnitude of error as compared with the quadratic element.

• Nuclear Engineering and Technology
• /
• v.52 no.2
• /
• pp.223-229
• /
• 2020

The present work proposes a solution to the static Boltzmann transport equation approximated by the simplified P₃ (SP₃) on angular, and the analytic coarse mesh finite difference (ACMFD) for spatial variables. Multi-group SP₃-ACMFD equations in 3D rectangular geometry are solved using the GMRES solution technique. As the core time dependent analysis necessitates the solution of an eigenvalue problem for an initial condition, this work is hence devoted to development and verification of the proposed static SP₃-ACMFD solver. A 3D multi-group static diffusion solver is also developed as a byproduct of this work to assess the improvement achieved using the SP₃ technique. Static results are then compared against transport benchmarks to assess the proximity of SP₃-ACMFD solutions to their full transport peers. Results prove that the approach can be considered as an acceptable interim approximation with outputs superior to the diffusion method, close to the transport results, and with the computational costs less than the full transport approach. The work would be further generalized to time dependent solutions in Part II.

• Journal of Theoretical and Applied Mechanics
• /
• v.2 no.1
• /
• pp.31-50
• /
• 1996

In this paper, a finite element analysis based on the local approach concept to fracture in the continuum damage mechanics is performed to analyze ductile fracture in two dimensional quasi-static state. First an isotropic damage model based on the generalized concept of effective stress is proposed for structural materials in the context of large deformation. In this model, the stiffness degradation is taken as a measure of damage and so, the fracture phenomenon can be explained as the critical deterioration of stiffness at a material point. The modified Riks' continuation technique is used to solve incremental iterative equations. Crack propagation is achieved by removing critically damaged elements. The mesh size sensitivity analysis and the simulation of the well known shearing mode failure in plane strain state are carried out to verify the present formulation. As numerical examples, an edge cracked plate and the specimen with a circular hole under plane stress are taken. Load-displacement curves and successively fractured shapes are shown. From the results, it can be concluded that the proposed model based on the local approach concept in the continuum damage mechanics may be stated as a reasonable tool to explain ductile fracture initiation and crack propagation.