• Structural Engineering and Mechanics
• /
• v.14 no.1
• /
• pp.71-84
• /
• 2002

A finite element formulation for the time-domain analysis of linear transient elastodynamic problems is presented based on the weak form obtained by applying the Galerkin's method to the integro-differential equations which contain the initial conditions implicitly and does not include the inertia terms. The weak form is extended temporally under the assumptions of the constant and linear time variations of field variables, since the time-stepping algorithms such as the Newmark method and the Wilson ${\theta}$-method are not necessary, obtaining two kinds of implicit finite element equations which are tested for numerical accuracy and convergency. Three classical examples having finite and infinite domains are solved and numerical results are compared with the other analytical and numerical solutions to show the versatility and accuracy of the presented formulation.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2001.10a
• /
• pp.58-61
• /
• 2001

An implicit finite element formulation for axisymmetric tube hydroforming is investigated. In order to describe normal anisotropy of the tube, Hill's non-quadratic yield function is employed. The frictional contact between die and tube and frictionless contact between tube and fluid are considered using the mesh-normal vector computed from finite element mesh of the tube. In order to verify the validity of the developed finite element formulation, the axisymmetric tube bulge test is simulated and simulation results are compared with experimental measurements. In the axisymmetric tube hydroforming process, an optimal hydraulic curve is pursued by performing the simulation with various internal pressures and axial forces.

• Transactions of Materials Processing
• /
• v.11 no.1
• /
• pp.75-83
• /
• 2002

Recently, the hydroforming process is widely applied to the automotive industry and rapidly spreaded to other industries. In this paper, An implicit finite element formulation for simulating axisymmetric tube hydroforming processes is performed. In order to describe normal anisotropy of the tube, Hill's non-quadratic yield function is employed. The frictional contact between die and tube and the frictionless contact between tube and fluid are considered using the mesh-normal vectors computed from the finite element mesh of the tube. The complete set of the governing relations comprising equilibrium and interfacial equations is linearized for Newton-Raphson procedure. In order to verify the validity of the developed finite element formulation, the axisymmetric tube bulge test is simulated and the simulation results are compared with experimental measurements. In a simulation of stepped circular tube hydroforming processes, an optimal hydraulic pressure curve is pursued by considering simultaneously internal pressures and axial forces.

• Transactions of Materials Processing
• /
• v.9 no.6
• /
• pp.590-597
• /
• 2000

A finite element formulation lot the simulation of axisymmetric sheet hydroforming is proposed, and an implicit program is coded. In order to describe normal anisotropy of steel sheet, Hill's non-quadratic yield function (Hill, 1979) is employed. Frictional contacts among sheet surface, rigid tool surface, and flexible hydrostatic pressure are considered using mesh normal vectors based on finite element of the sheet. Applied hydraulic pressure is also considered as a function of forming rate and time and treated as an external loading. The complete set of the governing relations comprising equilibrium and interfacial equations is approximately linearized for Newton-Raphson algorithm. In order to verify the validity of the developed finite element formulation, the axisymmetric bulge test is simulated. Simulation results are compared with other FEM results and experimental measurements and showed good agreements. In axisymmetric hydroforming processes of a disk cover, formability changes are observed according to the hydraulic pressure curve changes.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.435-447
• /
• 2019

Energy Flow Finite Element Analysis (EFFEA) is a promising tool for predicting dynamic energetics of complicated structures at high frequencies. In this paper, the Energy Flow Finite Element (EFFE) formulation of complicated Mindlin plates was newly developed to improve the accuracy of prediction of the dynamic characteristics in the high frequency. Wave transmission analysis was performed for all waves in complicated Mindlin plates. Advanced Energy Flow Analysis System (AEFAS), an exclusive EFFEA software, was implemented using $MATLAB^{(R)}$. To verify the general power transfer relationship derived, wave transmission analysis of coupled semi-infinite Mindlin plates was performed. For numerical verification of EFFE formulation derived and EFFEA software developed, numerical analyses were performed for various cases where coupled Mindlin plates were excited by a harmonic point force. Energy flow finite element solutions for coupled Mindlin plates were compared with the energy flow solutions in the various conditions.

• Bulletin of the Korean Mathematical Society
• /
• v.51 no.1
• /
• pp.267-285
• /
• 2014

An upstream scheme based on the pseudostress-velocity mixed formulation is studied to solve convection-dominated Oseen equations. Lagrange multipliers are introduced to treat the trace-free constraint and the lowest order Raviart-Thomas finite element space on rectangular mesh is used. Error analysis for several quantities of interest is given. Particularly, first-order convergence in $L^2$ norm for the velocity is proved. Finally, numerical experiments for various cases are presented to show the efficiency of this method.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1998.06a
• /
• pp.65-73
• /
• 1998

The design analysis of axisymmetric, multi-stage deep drawing dies was performed using the rigid-viscoplastic finite element formulation. In the formulation, the axisymmetric CFS algorithm was employed. Hill's non-quadratic normal anisotropic yield criterion and isotropic hardening rule were considered. For trial initial displacements and tool contact points, the geometric force equilibrium method was adopted. In order to see the validity of the formulation, the multi-stage deep drawing processes of shell-cylinder front part of hydraulic booster were simulated. The simulation showed good agreements with measurements and PAM-STAMP.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.205-216
• /
• 1997

In order to trace the lateral post-buckling behaviors of thin-wafled space frames, a geometrically nonlinear finite element formulation is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for symmetric thin-walled cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, tangent stiffness matrices of the thinwalled space frame element with 7 degrees of freedom including the restrained warping for each node are derived by using the Hermition polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines of the frame element due to rigid body rotations and incremental member forces from pure deformations. Finite element solutions for the spatial buckling and post-buckling analysis of thin-walled space frames are presented and compared with available solutions and other researcher's results.

• Computational Structural Engineering
• /
• v.10 no.1
• /
• pp.201-211
• /
• 1997

A clearly consistent finite element formulation for geometrically non-linear analysis of space frames is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for symmetric cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, elastic and geometric stiffness matrices of the space frame element are derived by using the Hermitian polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines of the frame element due to rigid body rotations and incremental member forces from pure deformaions. Finite element solutions for the spatial buckling and post-buckling analysis of space frames are compared with available solutions and other researcher's results.

• Structural Engineering and Mechanics
• /
• v.11 no.5
• /
• pp.547-564
• /
• 2001

The numerical formulation of a two-phase interface element appropriate for porous lining system is presented. The formulation is isoparametric and can be applied both for 2-D and 3-D analysis. Biot's theory is utilized as the basis for the development of the element constitutive theory. In order to be capable of simulating the reinforcing characteristics of some geotextiles utilized as lining system, a reinforcement component has also been implemented into the formulation. By employing this specially developed interface finite element, the influence of soil consolidation on the stress distribution along the lining system of a reservoir and a landfill has been investigated.