• Proceedings of the KSME Conference
• /
• 2001.06c
• /
• pp.402-407
• /
• 2001

Sensitivity analysis scheme is developed in the elasto-plastic finite element method with explicit time integration using direct differentiation method. The direct differentiation is concerned with the time integration, constitutive relation, shell element with reduced integration and the contact scheme. Sensitivity analysis results are mainly examined with the highly nonlinear and quasi-static problem with the complicated contact condition. The result shows stable sensitivity especially in the sheet metal forming analysis.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.474-479
• /
• 2004

According to our previous study, it is confirmed that the Petrov-Galerkin natural element method (PGNEM) completely resolves the numerical integration inaccuracy in the conventional Bubnov-Galerkin natural element method (BG-NEM). This paper is an extension of PG-NEM to two-dimensional nonlinear dynamic problem. For the analysis, a constant average acceleration method and a linearized total Lagrangian formulation is introduced with the PG-NEM. At every time step, the grid points are updated and the shape functions are reproduced from the relocated nodal distribution. This process enables the PG-NEM to provide more accurate and robust approximations. The representative numerical experiments performed by the test Fortran program, and the numerical results confirmed that the PG-NEM effectively and accurately approximates the nonlinear dynamic problem.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.12
• /
• pp.2982-2994
• /
• 1993

Low-velocity impact responses of composite laminates are investigated using the finite element method based on various theories. In two-dimensional nonlinear analysis, a displacement field considering higher order shear deformation and large deflection of the laminate is assumed and a finite element formulation is developed using a C$^{o}$-continuous 9-node plate element. Also, three-dimensional linear analysis based on the infinitesimal strain-displacement assumptions is performed using 8-node brick elements with incompatible modes. A modified Hertzian contact law is incorporated into the finite element program to evaluate the impact force. In the time integration, the Newmark constant acceleration algorithm is used in conjuction with successive iterations within each time step. Numerical results from static analysis as well as the impact response analysis are presented including impact force histories, deflections, strains in the laminate. Impact responses according to two typical low-velocity impact conditions are compared each other.

• Journal of Mechanical Science and Technology
• /
• v.18 no.3
• /
• pp.460-470
• /
• 2004

A computational analysis of engineering problems with moving domain or/and boundary according to either Lagrangian or Eulerian approach may encounter inherent numerical difficulties, the extreme mesh distortion in the former and the material boundary indistinctness in the latter. In order to overcome such defects in classical numerical approaches, the ALE(arbitrary Lagrangian Eulerian) method is widely being adopted in which the finite element mesh moves with arbitrary velocity. This paper is concerned with the ALE finite element formulation, aiming at the dynamic response analysis of baffled fuel-storage container in yawing motion, for which the coupled time integration scheme, the remeshing and smoothing algorithm and the mesh velocity determination are addressed. Numerical simulation illustrating theoretical works is also presented.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.51-56
• /
• 2000

The Finite Element Solutions Is reported on solid-liquid phase change in porous media with natural convection including freezing. The model is based on volume averaged transport equations, while phase change is assumed to occur over a small temperature range. The FEM (Finite Element Method) algorithm used in this study is 3-step time-splitting method which requires much less execution time and computer storage the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the energy equation. For natural convection including melting and solidification the numerical results show reasonable agreement with FDM (Finite Difference Method) results.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.8 no.3
• /
• pp.151-162
• /
• 2000

This study is concerned with computerized multi-level substructuring methods and stress analysis of coil springs. The purpose of substructuring methods is to reduce computing time and capacity of computer memory by multiple level reduction of the degrees of freedom in large size problems that are modeled by three dimensional continuum finite elements. In this paper, the spring super element developed is investigated with tension, torsion, and bending of a cylindrical bar in order to verify its accuracy and efficiency for the multi-level substructuring method. And then the algorithm is applied to finite element analysis of coil springs. The result demonstrates the validity of the multi-level substructuring method and the efficiency in computing time and memory by providing good computational results in coil spring analysis.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.11 s.242
• /
• pp.1439-1444
• /
• 2005

A two-dimensional boundary element method was used for the scattering analysis of side-drilled hole(SDH). The far-field scattering amplitude was calculated for shear vertical(SV) wave, and their frequency and time-domain results were presented. The time-domain scattering amplitude showed the directly reflected wave from the SDH leading edge as well as the creeping wave. In an immersion, pulse-echo testing, two measurement models were introduced to predict the response from SDHs. The 2-D boundary element scattering amplitude was converted to the 3-D amplitude to be used in the measurement model. The receiver voltage was calculated fer SV wave incidence at 45$^{\circ}C$ on the 1 m diameter SDH, and the result was compared with experiment.

• Structural Engineering and Mechanics
• /
• v.17 no.6
• /
• pp.735-749
• /
• 2004

This work presents a direct integration scheme, based on a fourth order finite difference approach, for elastodynamics. The proposed scheme was chosen as an alternative for attenuating the errors due to the use of the central difference method, mainly when the time-step length approaches the critical time-step. In addition to eliminating the spurious numerical oscillations, the fourth order finite difference scheme keeps the advantages of the central difference method: reduced computer storage and no requirement of factorisation of the effective stiffness matrix in the step-by-step solution. A study concerning the stability of the fourth order finite difference scheme is presented. The Finite Element Method and the Boundary Element Method are employed to solve elastodynamic problems. In order to verify the accuracy of the proposed scheme, two examples are presented and discussed at the end of this work.

• Structural Engineering and Mechanics
• /
• v.47 no.1
• /
• pp.27-44
• /
• 2013

A finite element computational procedure for the accurate analysis of quasistatic thermorheological complex structures response is developed. The geometrical nonlinearity, arising from large displacements and rotations (but small strains), is accounted for by the total Lagrangian description of motion. The Schapery's nonlinear single-integral viscoelastic constitutive model is modified for a time-stress-temperature-dependent behavior. The nonlinear thermo-viscoelastic constitutive equations are incrementalized leading to a recursive relationship and thereby the resulting finite element equations necessitate data storage from the previous time step only, and not the entire deformation history. The Newton-Raphson iterative scheme is employed to obtain a converged solution for the non-linear finite element equations. The developed numerical model is verified with the previously published works and a good agreement with them is found. The applicability of the developed model is demonstrated by analyzing two examples with different thermal/mechanical loading histories.

• Journal of the Korean Society for Precision Engineering
• /
• v.14 no.9
• /
• pp.90-100
• /
• 1997

A full three-dimensional thermo-coupled rigid-viscoplastic finite element method and the currently developed microstructural evolution system which includes semi-empirical equations suggested by different research groups were used together to form an integrated system of process and micro- structure simulation of hot rolling. The distribution and time histroy of the momechanical variables such as temperature, strain, strain rate, and time during pass and between passes were obtained from the finite element analysis of multipass hot rolling processes. The distribution of metallurgical variables were calculated on the basis of instantaneous thermomechanical data. For the verification of this method the evolution of microstructure in plate rolling and shape rolling was simulated and their results were compared with the data available in the literature. Consequently, this approach makes it possible to describe the realistic evolution of microstructure by avoiding the use of erroneous average value and can be used in CAE of multipass hot rolling.