• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제5권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.

• 한국산업융합학회 논문집
• /
• 제3권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.

• 한국정밀공학회:학술대회논문집
• /
• 한국정밀공학회 1993년도 추계학술대회 논문집
• /
• 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.

• 한국전산구조공학회논문집
• /
• 제12권3호
• /
• pp.363-370
• /
• 1999

본 연구에서는 철근콘크리트기둥과 철골보로 이루어진 혼합구조 접합부의 해석에 대한 유한 유소법을 이용한 해석 모델 방법을 제시하였다. 혼합구조 접합부에서 콘크리트와 강판이 접하는 접촉면은 두 접촉면 사이를 부착과 마찰의 개념으로 표현할 수 있는 주-종속 접촉 알고리즘(master-slave contact algorithm)을 이용하여 모델링하였다. 그리고, 휨응력의 지배를 받는 강관에는 비적합 모드 요소를 사용하였다. 본 연구에서의 혼합구조 특징은 보에서 기둥으로 힘의 전달을 원활히 하기 위하여 다이아프램이 사용되었고, 이러한 혼합구조 접합부 모델링 방법에 대한 타당성을 알아보기 위하여 3차원 비선형 해석을 행하여 실험결과와 비교한 결과 잘 일치하는 결과를 얻었다.

• 대한기계학회논문집
• /
• 제18권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
• /
• 제7권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.

• 한국산업융합학회 논문집
• /
• 제7권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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2000년도 추계학술대회논문집A
• /
• 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.

• 대한수학회보
• /
• 제51권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
• /
• 제5권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.