• Structural Engineering and Mechanics
• /
• 제6권6호
• /
• pp.711-725
• /
• 1998

Mixed finite element formulations for incompressible materials show pressure oscillations or pressure modes in four-node quadrilateral elements. The criterion for the stability in the pressure solution is the so-called Babu$\check{s}$ka-Brezzi stability condition, and the four-node elements based on mixed variational principles do not appear to satisfy this condition. In this study, a pressure continuity residual based on the pressure discontinuity at element edges proposed by Hughes and Franca is used to study the stabilization of pressure solutions in bilinear displacement-constant pressure four-node quadrilateral elements. Also, a solid mechanics problem is presented by which the stability of mixed elements can be studied. It is shown that the pressure solutions, although stable, are shown to exhibit sensitivity to the stabilization parameters.

• Korean Journal of Mathematics
• /
• 제13권2호
• /
• pp.127-137
• /
• 2005

The objective of numerical analysis is to devise and analyze efficient algorithms or numerical methods for equations arising in mathematical modeling for science and engineering. In this article, we present some recent topics in computational mathematics, specially in the finite element method and overview the development of the mixed finite element method in the context of second order elliptic and parabolic problems. Multiscale methods such as MsFEM, HMM, and VMsM are included.

• 동력기계공학회지
• /
• 제4권1호
• /
• pp.68-73
• /
• 2000

A mixed degree finite element solutions using hierarchical elements are investigated for convergences on a 2-D simple cases. Elements are generated block by block and each block is assigned an arbitrary solution degree. The numerical study showed that a well constructed blocks can increase the convergence and accuracy of finite element solutions. Also, it has been found that for higher order elements, the convergence trends can be deteriorated for smaller mesh sizes. A procedure for a variable fixed boundary condition has been included.

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

• 대한수학회보
• /
• 제51권6호
• /
• pp.1655-1668
• /
• 2014

In this article, we propose and analyze a new nonconforming primal mixed finite element method for the stationary Stokes equations. The approximation is based on the pseudostress-velocity formulation. The incompressibility condition is used to eliminate the pressure variable in terms of trace-free pseudostress. The pressure is then computed from a simple post-processing technique. Unique solvability and optimal convergence are proved. Numerical examples are given to illustrate the performance of the method.

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

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제15권1호
• /
• pp.57-65
• /
• 2011

In this paper, a good quality mesh generation for the finite element method is investigated for artificial hip joint simulations. In general, bad meshes with a large aspect ratio or mixed elements can give rise to excessively long computational running times and extremely high errors. Typically, hexahedral elements outperform tetrahedral elements during three-dimensional contact analysis using the finite element method. Therefore, it is essential to mesh biologic structures with hexahedral elements. Four meshing schemes for the finite element analysis of an artificial hip joint are presented and compared: (1) tetrahedral elements, (2) wedge and hexahedral elements, (3) open cubic box hexahedral elements, and (4) proposed hexahedral elements. The proposed meshing scheme is to partition a part before seeding so that we have a high quality three-dimensional mesh which consists of only hexahedral elements. The von Mises stress distributions were obtained and analyzed. We also performed mesh refinement convergence tests for all four cases.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1995년도 가을 학술발표회 논문집
• /
• pp.153-160
• /
• 1995

Mixed finite element formulations for incompressible materials show pressure oscillations or pressure modes in four-node quadrilateral elements. The criterion for the stability in the pressure solution is the so-called Babufka-Brezzi stability condition, and the four-node elements based on mixed variational principles do not appear to satisfy this condition. In this study, a pressure continuity residual based on the pressure discontinuity at element edges is used to study the stabilization of pressure solutions in bilinear displacement-constant pressure four-node quadrilateral elements. It is shown that the pressure solutions, although stable, exhibit sensitivity to the stabilization parameters.

• 대한기계학회논문집B
• /
• 제26권9호
• /
• pp.1209-1217
• /
• 2002

The use of a two dimensional hierarchical elements are investigated for the incompressible flow computation. The construction of hierarchical elements are explained by both a geometric configuration and a determination of degrees of freedom. Also a systematic treatment of essential boundary values has been developed for the degrees of freedom corresponding to higher order terms. The numerical study for the poisson problem shows that the computation with hierarchical higher order elements can increase the convergence rate and accuracy of finite element solutions in more efficient manner than the use of standard first order element. for Stokes and Cavity flow cases, a mixed version of penalty function approach has been introduced in connection with the hierarchical elements. Solutions from hierarchical elements showed better resolutions with consistent trends in both mesh shapes and the order of elements.

• Structural Engineering and Mechanics
• /
• 제41권2호
• /
• pp.157-181
• /
• 2012

The paper presents a review of the application of the newly proposed mixed finite element model for seismic simulation of different types of composite frame structures. To evaluate the performance of the element, a comparison with displacement-based and force-based models is conducted. The study revealed that the mixed model is superior to the others in terms of both speed of convergence and numerical stability, and is therefore considered the most practical approach for modeling of composite structures. In this model, the element is derived using independent force and displacement shape functions. The nonlinear response of the frame element is based on the section discretization into fibers with uniaxial material models. The interfacial behavior is modeled using an inelastic interface element. Numerical examples to clarify the advantages of the model are presented for the following structural applications: anchored reinforcing bar problems, composite steel-concrete girders with deformable shear connectors, beam on elastic foundation elements, R/C girders strengthened with FRP sheets, R/C beam-columns with bond-slip, and prestressed concrete girders. These studies confirmed that the model represents a major advancement over existing elements in simulating the inelastic behavior of composite structures.