• Structural Engineering and Mechanics
• /
• v.3 no.4
• /
• pp.391-410
• /
• 1995

The application of adaptive finite element method to dynamic problems is investigated. Both the kinetic and strain energy errors induced by space and time discretization were estimated in a consistent manner and controlled by the simultaneous use of the adaptive mesh generation and the automatic time stepping. Also an optimal ratio of spatial discretization error to temporal discretization error was discussed. In this study it was found that the best performance can be obtained when the specified spatial and temporal discretization errors have the same value. Numerical examples are carried out to verify the performance of the procedure.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.6 s.237
• /
• pp.703-710
• /
• 2005

A finite element discretization of the advection and redistancing equations of level set method has been studied. It has been shown that Galerkin spatial discretization combined with Crank-Nicolson temporal discretization of the advection equation of level set yields a good result and that consistent streamline upwind Petrov-Galerkin(CSUPG) discretization of the redistancing equation gives satisfactory solutions for two test problems while the solutions of streamline upwind Petrov-Galerkin(SUPG) discretization are dissipated by the numerical diffusion added for the stability of a hyperbolic system. Furthermore, it has been found that the solutions obtained by CSUPG method are comparable to those by second order ENO method.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.1 s.71
• /
• pp.39-46
• /
• 2006

This paper analyzes the continuous Galerkin method for the space-time discretization of wave equation. The method of space-time finite elements enables the simple solution than the usual finite element analysis with discretization in space only. We present a discretization technique in which finite element approximations are used in time and space simultaneously for a relatively large time period called a time slab. The weighted residual process is used to formulate a finite element method for a space-time domain. Instability is caused by a too large time step in successive time steps. A stability problem is described and some investigations for chosen types of rectangular space-time finite elements are carried out. Some numerical examples prove the efficiency of the described method under determined limitations.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.4
• /
• pp.203-214
• /
• 2017

In this paper we propose and analyze two a posteriori error estimators for the stabilized $P_1/P_1$ finite element discretization of the Stokes equations. These error estimators are computed by solving local Poisson or Stokes problems on elements of the underlying triangulation. We establish their asymptotic exactness with respect to the velocity error under certain conditions on the triangulation and the regularity of the exact solution.

• Journal of applied mathematics & informatics
• /
• v.29 no.5_6
• /
• pp.1129-1141
• /
• 2011

In this paper, a class of large time-stepping method based on the finite element discretization for the Cahn-Hilliard equation with the Neumann boundary conditions is developed. The equation is discretized by finite element method in space and semi-implicit schemes in time. For the first order fully discrete scheme, convergence property is investigated by using finite element analysis. Numerical experiment is presented, which demonstrates the effectiveness of the large time-stepping approaches.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2067-2077
• /
• 1996

In adaptive finite element analysis, r- and h-methods are generally used on the basis of a discretization error estimator. In this paper, an rh-method is proposed as a new adaptive method which can improve the adaptivity performance by using both of them. This suggested rh-method moves nodal coordinates of initially given model to adjust element discretization errors and thereafter performes the h-method tdo obtain the specified accuracy of finite element solutions. Numerical experiments for various plane problems were performed using 4-noded isoparametric quadrilateral elements. As a result, the rh-method has been shown to be an accurate and efficient adaptive analysis method to obtain as improved solution.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.2
• /
• pp.141-149
• /
• 2002

Although the finite element method has become an indispensible tool for the dynamic analysis of structures, difficulty remains to quantify the errors associated with discretization. To improve the modeling accuracy, this paper proposes a method to make a combined use of finite elements and exact dynamic elements. Exact interpolation functions for the Timoshenko beam element are derived using the exact dynamic element modeling (EDEM) and compared with interpolation functions of the finite element method (FEM). The exact interpolation functions are tested with the Laplace variable varied. A combined use of finite element method and exact interpolation functions is presented to gain more accurate mode shape functions. This paper also presents a combined use of finite elements and exact dynamic elements in design/reanalysis problems. Timoshenko flames with tapered sections are tested to demonstrate the design procedure with the proposed method. The numerical study shows that the combined use of finite element model and exact dynamic element model is very useful.

• Structural Engineering and Mechanics
• /
• v.86 no.3
• /
• pp.349-359
• /
• 2023

This paper presents a technique for determining the optimal number of elements in stochastic finite element analysis based on reliability analysis. Using the change-of-variable perturbation stochastic finite element approach, the probability density function of the dynamic responses of stochastic structures is explicitly determined. This method combines the perturbation stochastic finite element method with the change-of-variable technique into a united model. To further examine the relationships between the random fields, discretization of the random field parameters, such as the variance function and the scale of fluctuation, is also performed. Accordingly, the reliability index is calculated based on the explicit probability density function of responses with Gaussian or non-Gaussian random fields in any number of elements corresponding to the random field discretization. The numerical examples illustrate the effectiveness of the proposed method for a one-dimensional cantilever reinforced concrete column and a two-dimensional steel plate shear wall. The benefit of this method is that the probability density function of responses can be obtained explicitly without the use simulation techniques. Any type of random variable with any statistical distribution can be incorporated into the calculations, regardless of the restrictions imposed by the type of statistical distribution of random variables. Consequently, this method can be utilized as a suitable guideline for the efficient implementation of stochastic finite element analysis of structures, regardless of the statistical distribution of random variables.

• Structural Engineering and Mechanics
• /
• v.75 no.6
• /
• pp.685-699
• /
• 2020

Polygonal finite element provides a great flexibility in mesh generation of crack propagation problems where the topology of the domain changes significantly. However, the control of the discretization error in such problems is a main concern. In this paper, a polygonal-FEM is presented in modeling of crack propagation problems via an automatic adaptive mesh refinement procedure. The adaptive mesh refinement is accomplished based on the Zienkiewicz-Zhu error estimator in conjunction with a weighted SPR technique. Adaptive mesh refinement is employed in some steps for reduction of the discretization error and not for tracking the crack. In the steps that no adaptive mesh refinement is required, local modifications are applied on the mesh to prevent poor polygonal element shapes. Finally, several numerical examples are analyzed to demonstrate the efficiency, accuracy and robustness of the proposed computational algorithm in crack propagation problems.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.04a
• /
• pp.1-8
• /
• 1994

An automated procedure which allows adaptation of spatial and time discretization simultaneously in finite element analysis of linear elastodynamic problems is presented. For dynamic problems having responses dominated by high frequency modes, such as those with impact, explosive, traveling and earthquake loads high gradient stress regions change their locations from time to time. And the time step size may need to vary in order to deal with whole process ranging from transient phase to steady state phase. As the sizes of elements in space vary in different regions, the procedure also permits different time stepping. In such a way, the best performance attainable by the finite element method can be achieved. In this study, we estimate both of the kinetic energy error and stran energy error induced by spatial and time discretization in a consistent manner. Numerical examples are used to demonstrate the performance of the procedure.