• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2003.05a
• /
• pp.461-464
• /
• 2003

A new control volume finite element method is proposed for three dimensional analysis of polymer flow. Tetrahedral finite element is employed and co-located interpolation procedure for pressure and velocity is implemented. Inclusion of pressure gradient term in the velocity shape functions prevents the checkerboard pressure field from being developed. Vectorial nature of pressure gradient is considered in the velocity shape function so that velocity profile in the limit of very small Reynolds number becomes physically meaningful. The proposed method was verified through three dimensional simulation of pipe flow problem for Newtonian and power-law fluid. Calculated pressure and velocity field showed an excellent agreement with analytic solutions for pressure and velocity. Driven-cavity problem, which is reported to yield checkerboard pressure filed when conventional finite element method is applied, could be solved without yielding checkerboard pressure field when the proposed control volume finite element method was applied. The proposed method could be successfully applied to the three dimensional mold filling problem.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.2 no.2
• /
• pp.1-19
• /
• 1998

A Chebyshev pseudospectral-finite element method is proposed for two-dimensional unsteady Navier-Stokes equation. The generalized stability and the convergence are proved strictly. The numerical results show the advantages of this method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.12 no.2
• /
• pp.69-79
• /
• 2008

Multigrid methods finite element/finite volume methods and their convergence properties are reviewed in a general setting. Some early theoretical results in simple finite element methods in variational setting method are given and extension to nonnested-noninherited forms are presented. Finally, the parallel theory for nonconforming element[13] and for cell centered finite difference methods [15, 23] are discussed.

• Journal of Power System Engineering
• /
• v.21 no.6
• /
• pp.40-45
• /
• 2017

The objective of this study is the finite element-transfer stiffness coefficient method, which is the combination of the modeling technique of finite element method and the transfer technique of transfer stiffness coefficient method, is applied in the static analyses of two dimensional curved beam structures. To confirm the effectiveness of the applied method, two computational models are selected and analyzed by using finite element method, finite element-transfer stiffness coefficient method and exact solution. The computational results of the static analyses for two computational models using finite element-transfer stiffness coefficient method are equal to those using finite element method. When the element partition number of curved beam structure is increased, the computational results of the static analyses using both methods approach the exact solution. We confirmed that the finite element-transfer stiffness coefficient method is superior to finite element method when the number of the curved beam elements is increased from the viewpoints of the computational speed and the utility of computer memory.

• Journal of the Korean Geosynthetics Society
• /
• v.20 no.3
• /
• pp.11-20
• /
• 2021

In this paper, the approximately coupled method of finite element method and boundary element method to obtain efficient and accurate analysis results is proposed for a two-dimensional elasto-static problem with a geometrically abruptly changing part. As the finite element of a two-dimensional problem, three-node and four-node plane stress element is applied, and as the boundary element of a two-dimensional problem, three-node boundary element is applied. In the modeling stage, firstly, an entire analysis target object is modeled as finite elements, and then a geometrically abruptly changing part is modeled as boundary elements. The boundary element is defined using the nodes defined for modeling finite elements. In the analysis stage, finite element analysis is firstly performed on a entire analysis target object, and boundary element analysis is automatically performed afterwards. As for the boundary conditions at boundary element analysis, displacement conditions and stress conditions, which are the results of finite element analysis, are applied. As a numerical example, the analysis results for a two-dimensional elasto-static problem, a plate with a crack, are presented and investigated.

• Structural Engineering and Mechanics
• /
• v.85 no.2
• /
• pp.265-274
• /
• 2023

The spherical steel bearings (SSBs) has been gradually replaced traditional rubber bearings and extensively applied to high-speed railway (HSR) bridges in China, due to their durability and serviceability. Nevertheless, SSB is generally simplified to the ordinary constraints in the finite element model, which cannot reflect its detailed mechanical characteristics, especially its seismic performance. To provide a more precisely simulation, an innovative and simplified finite element simulation method is proposed and the combined element group is developed in ANSYS. The primary parameters were determined by means of the performance test of SSB. The finite element model of SSB applied to a single-span HSR simply supported girder bridge was established through the proposed method. The seismic performance of the SSB was further investigated. A shake table test was conducted to evaluate the accuracy of the proposed simulation method. It is found that the numerical results could have a good agreement with the experiment, namely, the proposed method is feasible and efficient.

• Journal of the Korean Society for Precision Engineering
• /
• v.23 no.6 s.183
• /
• pp.128-135
• /
• 2006

Rational dynamic modeling and analysis method f3r complex structures are studied with special attention to slide way joints. For modeling of slide way joints, a general modeling technique is used by using the influence coefficients method which is applied to the conversion of detailed finite element model to the equivalent reduced joint model. The theoretical part of this method is illustrated and the method is applied to the structure with slide way joint. In this method, the non-linearity of the contact surfaces is considered within a proper range and the boundary effect of the joint model could be eliminated. The proposed method was applied to finite element modal analysis of a clamp jointed cantilever beam and slide way joints of the vertical type lathe. The method can also be used to other kinds of joint modeling. The results of these analysis were compared with those of Yoshimura models and rigid joint models, which demonstrated the practical applicability of the proposed method.

• Journal of applied mathematics & informatics
• /
• v.32 no.5_6
• /
• pp.761-781
• /
• 2014

An exponentially accurate nonconforming spectral element method for elasticity systems with discontinuities in the coefficients and the flux across the interface is proposed in this paper. The method is least-squares spectral element method. The jump in the flux across the interface is incorporated (in appropriate Sobolev norm) in the functional to be minimized. The interface is resolved exactly using blending elements. The solution is obtained by the preconditioned conjugate gradient method. The numerical solution for different examples with discontinuous coefficients and non-homogeneous jump in the flux across the interface are presented to show the efficiency of the proposed method.

• Coupled systems mechanics
• /
• v.2 no.4
• /
• pp.389-410
• /
• 2013

Three-dimensional Lagrangian fluid finite element is applied to seismic response analysis of an oil storage tank with a floating roof. The fluid element utilized in the present analysis is formulated based on the displacement finite element method considering only volumetric elasticity and its element stiffness matrix is derived by using one-point integration method in order to avoid volumetric locking. The method usually adds a rotational penalty stiffness to satisfy the irrotational condition for fluid motion and modifies element mass matrices through the projected mass method to suppress spurious hourglass-mode appeared in compensation for one-point integration. In the fluid element utilized in the present paper, a small hourglass stiffness is employed. The fluid and structure domains for the objective oil storage tank are modeled by eight-node solid elements and four-node shell elements, respectively, and the transient response of the floating roof structure or the free surface are evaluated by implicit direct time integration method. The results of seismic response analyses are compared with those by other method and the validation of the present analysis using three-dimensional Lagrangian fluid finite elements is shown.

• Journal of applied mathematics & informatics
• /
• v.36 no.3_4
• /
• pp.257-270
• /
• 2018

An extrapolated Crank-Nicolson characteristic finite element method is introduced for approximate solutions of nonlinear Sobolev equations with a convection term. And we obtain the higher order of convergence for approximate solutions in the temporal and the spatial directions with respect to $L^2$ norm.