• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.1
• /
• pp.55-78
• /
• 2001

Given the anisotropic Poisson equation $-{\nabla}{\cdot}{\mathcal{K}}{\nabla}p=f$, one can convert it into a system of two first order PDEs: the Darcy law for the flux $u=-{\mathcal{K}{\nabla}p$ and conservation of mass ${\nabla}{\cdot}u=f$. A very natural mixed finite volume method for this system is to seek the pressure in the nonconforming P1 space and the Darcy velocity in the lowest order Raviart-Thomas space. The equations for these variables are obtained by integrating the two first order systems over the triangular volumes. In this paper we show that such a method is really a standard finite element method with local recovery of the flux in disguise. As a consequence, we compare two approaches in analyzing finite volume methods (FVM) and shed light on the proper way of analyzing non co-volume type of FVM. Numerical results for Dirichlet and Neumann problems are included.

• Journal of the Korean Society of Industry Convergence
• /
• v.3 no.4
• /
• pp.329-336
• /
• 2000

A geometric and material nonlinear finite element analysis is developed for the layered elastomeric bearings. In this study, a mixed variational approach with separate variables is used to describe the displacement and volume change of rubber. To represent finely deformed behavior, Kirchoff stress tensors are used and converted Eulerian stress tensors to describe real physical meanings. Newton's method is utilized to solve the governing nonlinear finite element equations. Numerical test are performed in the case of compression and shear to verify the theory and to illustrate the application of this analysis. And the results of this study were compared to the results of Moore's discrete finite element analysis.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1993.10a
• /
• pp.644-649
• /
• 1993

Practical structures are subject not only to tension but also to shear and torsional loading. In this study, the mode 1 and 2 stress intensity factors of specimens were calculated by using elastic finite element mothod. The stress fields at the crack tip subjected to mixed-mode loading were also studied by usingf eleatic finite element method and were compared with theoretical results. The three-point-bending, four-point-bending, and mixed-mode-loading experiment were carried out. And, crack propagation rate da/dN and crack growth direction were examined. Also, the elastic finite element method was applied to calculate the stress intensity factors of branch crack tip and we relate the stress intenity factor range of branch crack tip(the result of FEM) to crack propagation rate(the experimental result). The .DELTA. -da/dN relation corelated with that of mode 1.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.12 no.3
• /
• pp.363-370
• /
• 1999

This study presented an modeling method for the connections in mixed structure with reinforced concrete columns and steel beam using finite element method. The contacting surfaces between concrete and steel are modeled using master-slave contact algorithm and the incompatible mode elements were used in the steel tube subjected to bending. The characteristics of mixed structure was that diaphragm was used for transferring force from beam to column. The three dimensional nonlinear analysis was performed and the analytical results compared with experimental results in order to prove modeling method.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.624-634
• /
• 1994

A numerical method which combines equal-order velocity-pressure formulation originated from SIMPLE algorithm and streamline upwinding method has been developed. To verify the proposed numerical method, we considered the lid-driven cavity flow and backward facing step flow. The trend of convergence history is stable up to the error criterion beyond which the maximum value of error is oscillatory due4 to the round-off error. In the present study, all results were obtained with the single precision calculation up to the given error criterion and it was found to be sufficient for our purpose. The present results were then compared with existing experimental results using laser doppler velocimetry and numerical results using finite difference method and mixed interpolation finite element method. It has been shown that the present method gives accurate results with less memories and execution time than the coventional finite element method.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.7 no.2
• /
• pp.95-111
• /
• 2003

BDM mixed methods are obtained for a good approximation of velocity for flow equations. In this paper, we study an implementation issue of solving the algebraic system arising from the BDM mixed finite elements. First we discuss post-processing based on the use of Lagrange multipliers to enforce interelement continuity. Furthermore, we establish an equivalence between given mixed methods and projection finite element methods developed by Chen. Finally, we present the implementation of the first order BDM on rectangular grids and show it is as simple as solving the pressure equation.

• Journal of the Korean Society of Industry Convergence
• /
• v.7 no.4
• /
• pp.355-361
• /
• 2004

The finite element method were used to determine the stress intensity factor of cracked plate. The stress method, displacement method and J Integral are most popular finte element method. ANSYS proposed another a kind of displacement method. In this paper, it was examined that the accuracy and utility of the ANSYS method could believable to determine the stress intensity factors of centered inclined crack. Generally, inclined crack has two portion of stress intensity factors, tensile mode F1 and shear mode F2. For the purpose of increasing the accuracy of stress intensity factors, examined the effect of the numbers of nodes and elements, crack tip element size and number of partition of the crack tip vicinity. It was found that the method proposed by ANSYS is useful and has high accuracy. Accuracy of calculated stress intensity factors was increased by increase of the number of nodes and elements, and at the small size of crack tip elements can get more highly accuracy.

• Proceedings of the KSME Conference
• /
• 2000.11a
• /
• pp.458-463
• /
• 2000

The mechanical structures generally have discontinuous parts such as the cracks, notches and holes owing to various reasons. In this paper, in order to analyze effectively these singularity problems using the finite element method, a mixed analysis method which an analytical solution and finite element solutions are simultaneously used is newly proposed. As the analytical solution is used in the singularity region and the finite element solutions are used in the remaining regions except this singular zone, this analysis method reasonably provides for the numerical solution of a singularity problem. Through various numerical examples, it is shown that the proposed analysis method is very convenient and gives comparatively accurate solution.

• 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.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.1
• /
• pp.11-24
• /
• 2001

We introduce an assumption about smoothing operator for mixed formulations and show that convergence of Multigrid method for the mixed finite element formulation for the Linear Elasticity. And we show that Richardson and Kaczmarz smoothing satisfy this assumption.