• Computers and Concrete
• /
• v.8 no.1
• /
• pp.111-123
• /
• 2011

Severe element distortion problem is observed in finite element mesh while performing numerical simulations of high velocity steel projectiles penetration/perforation of concrete targets using finite element method (FEM). This problem of element distortion in Lagrangian formulation of FEM can be resolved by using element erosion methodology. Element erosion approach is applied in the finite element program by defining failure parameters as a condition for element elimination. In this study strain parameters for both compression and tension at failure are used as failure criteria. Since no direct method exists to determine these values, a calibration approach is used to establish suitable failure strain values while performing numerical simulations of ogive-nose steel projectile penetration/perforation into concrete target. A range of erosion parameters is suggested and adopted in concrete penetration/perforation tests to validate the suggested values. Good agreement between the numerical and field data is observed.

• Journal of the Korean Institute of Telematics and Electronics
• /
• v.26 no.11
• /
• pp.1644-1651
• /
• 1989

In this paper, a two-dimensional numerical analysis of HEMT's with gate length of 0.6um is performed. In this case, Control Volume Formulation method which has been used in the analysis of heat transfer and fluid flow is used as a numerical method. As a mobility model, empirical formula including the velocithy overshoot phenomena is used instead of two-piece mobility model. The results obtained from this numerical analysis(i.e., the region in which cahnnel is formed, the strength of electric field in the channel, the distribution of potential, and the distribution of electron concentration etc.)are in good agreement with the previous analytic results. And our results also show the parasitic MESFET's operation in the range of the high gate voltage.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.9
• /
• pp.1230-1243
• /
• 1997

A numerical technique for simulating incompressible viscous flow with free surface is presented. The flow field is obtained by penalty finite element formulation. In this work, a modified volume of fluid (VOF) method which is compatible with 4-node element is proposed to track the moving free surface. This scheme can be applied to irregular mesh system, and can be easily extended to three dimensional geometries. Numerical analyses were done for two benchmark examples, namely the broken dam problem and the solitary wave propagation problem. The numerical results were in close agreement with the existing data. Illustrative examples were studied to show the effectiveness of the proposed numerical scheme.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1154-1158
• /
• 2003

Finite element analysis of welding processes, which entail phase evolution, heat transfer and deformations, is considered in this paper. Attention focuses on numerical implementation of the thermo-elastic-plastic constitutive equation proposed by Leblond in consideration of the transformation plasticity. Based upon the multiplicative decomposition of deformation gradient, hyperelastic formulation is employed for efficient numerical integration, and the algorithmic consistent moduli for elastic-plastic deformations including transformation plasticity are obtained in the closed form. The convergence behavior of the present implementation is demonstrated via a couple of numerical examples.

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

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.4 no.3
• /
• pp.130-138
• /
• 1996

The transient incompressible flow behind the widely-spaced co-axial jet is numerically simulated using the random vortex method(RVM). This numerical approach is based on the Lagrangian approach for the vorticity formulation of the unsteady Navier-Stokes equations, utilizing vortex elements to account for the convection and diffusion processes. The effects of the mass flow rate of an annular air jet and a central fuel jet on the co-axial jet flow dynamics is investigated. To validate the present procedure, the numerical results are compared with the available experimental data the present procedure, the numerical results are compared with the available experimental data in terms of the centerline and off-centerline profiles of the mean axial velocity. Discrepancies between the RVM results and the measurements are discussed in detail.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.7 no.5
• /
• pp.141-153
• /
• 1999

A finite element procedure for the analysis of rubber-like hyperelastic material is developed. The volumetric incompressiblity conditions of the rubber deformation is included in the formulation by using penalty method. In this paper, the behavior of the rubber deformation is represented by hyperelastic constitutive relations based on a generalized Mooney-Rivlin model. The principle of virtual work is used to derive nonlinear finite element equation for the large displacement problem and presented in total-Lagrangian description. The finite element procedure using analytic differentiation resulted in very close solution to the result of the well known commercial packages NISAII AND ABAQUS. Numerical tests show that the results from the numerical differentiation method coincide very well with those from the analytic method and the well known commercial packages in static analysis. The convergency of rubber usingν iteration method is also discussed.

• Steel and Composite Structures
• /
• v.19 no.3
• /
• pp.679-695
• /
• 2015

A numerical solution using finite difference method to evaluate the thermal buckling of simply supported FGM plate with variable thickness is presented in this research. First, the governing differential equation of thermal stability under uniform temperature through the plate thickness is derived. Then, the governing equation has been solved using finite difference method. After validating the presented numerical method with the analytical solution, the finite difference formulation has been extended in order to include variable thickness. The accuracy of the finite difference method for variable thickness plate has been also compared with the literature where a good agreement has been found. Furthermore, a parametric study has been conducted to analyze the effect of material and geometric parameters on the thermal buckling resistance of the FGM plates. It was found that the thickness variation affects isotropic plates a bit more than FGM plates.

• Proceedings of the KIEE Conference
• /
• 1999.07a
• /
• pp.220-222
• /
• 1999

This paper presents a numerical analysis of the eddy current testing with cracked pipe using finite element method (FEM). ${\vec{A}},\;{\phi}-{\vec{A}}$ method is adopted for the formulation of 3-dimensional(3-D) FEM with the brick element. The cracks investigated here are the inner and outer surface of axial symmetry, 90 degree circular one. The algorithm of 3-D numerical analysis is employed for the axisymmetric pipe with the cracks. In order to verify the validity of 3-D numerical analysis, the results are compared with those of 2-D analysis with the same type of the model. The differential impedance is obtained by using energy method and its locus are various 8-shaped curves for each cracks. The ICCG method is used for the calculation of a matrix.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.528-531
• /
• 2004

For a simple one-dimensional ignition problem a mathematical model is described to investigate the difficulties in numerical simulations. Some computation results are obtained and comparison is made with analytical solution. Discussions are made on topics such as 1) coordinate transformation, 2) gas-phase and solid-phase analysis; (divergence form of the governing system, a finite-volume discretization, implicit time integration, upwind split flux, spatial accuracy improvement are described. Mass, reagent mass, and energy conservations are solved.), and 3) method to determine quantities on the burning surface (matching). Results obtained for small values of the non-dimensional pressure show a steady-combustion and good agreement with the analytical solution. Numerical instability appeared for larger values of the pressure, discussion on the cause of the problem is made. This effort is a part of a study of flame spread phenomena on solid propellant surface.