• Structural Engineering and Mechanics
• /
• 제30권1호
• /
• pp.119-129
• /
• 2008

A small strain and elastoplastic formulation of Polygonal Element Method (PEM) is developed for efficient analysis of elastoplastic solids. In this work, the polygonal elements are constructed based on traditional triangular finite meshes. The construction method of polygonal mesh can directly utilize the sophisticated triangularization algorithm and reduce the difficulty in generating polygonal elements. The Wachspress rational finite element basis function is used to construct the approximations of polygonal elements. The incremental variational form and a von Mises type model are used for non-linear elastoplastic analysis. Several small strain elastoplastic numerical examples are presented to verify the advantages and the accuracy of the numerical formulation.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1997년도 봄 학술발표회 논문집
• /
• pp.3-10
• /
• 1997

Meshless particle method is a numerical technique which does not use the concept of element. This method can easily handle special engineering problems which cause difficulty in the use of finite element method, however it has a drawback that essential boundary condition is not satisfied. In this paper, several studies for satisfying essential boundary conditions and enhancing the accuracy of solutions are discussed. Particular emphasis is placed on a new numerical technique in which finite elements are used on the boundaries to satisfy the essential boundary conditions and meshless particle method is used in the interior domain. For coupling of the two methods interface elements are introduced into the zone between the subdomains using meshless particle method and finite element method. The shape functions and the approximated displacement functions of the interface element are derived with the ramp function based on the shape function of finite elements. The whole numerical procedures are formulated by Galerkin method. Several numerical examples for enhancing the accuracy of solution in the meshless particle method and a new coupling method are presented.

• Structural Engineering and Mechanics
• /
• 제56권2호
• /
• pp.275-291
• /
• 2015

Adherence between reinforcement and the surrounding concrete is usually ignored in finite element analysis (FEA) of reinforced concrete (RC) members. However, load transition between the reinforcement and surrounding concrete effects RC members' behavior a great deal. In this study, the effects of bond-slip on the FEA of RC members are examined. In the analyses, three types of bond-slip modeling methods (perfect bond, contact elements and spring elements) and three types of reinforcement modeling methods (smeared, one dimensional line and three dimensional solid elements) were used. Bond-slip behavior between the reinforcement and surrounding concrete was simulated with cohesive zone materials (CZM) for the first time. The bond-slip relationship was identified experimentally using a beam bending test as suggested by RILEM. The results obtained from FEA were compared with the results of four RC beams that were tested experimentally. Results showed that, in FE analyses, because of the perfect bond occurrence between the reinforcement and surrounding concrete, unrealistic strains occurred in the longitudinal reinforcement. This situation greatly affected the load deflection relationship because the longitudinal reinforcements dominated the failure mode. In addition to the spring elements, the combination of a bonded contact option with CZM also gave closer results to the experimental models. However, modeling of the bond-slip relationship with a contact element was quite difficult and time consuming. Therefore bond-slip modeling is more suitable with spring elements.

• Journal of Theoretical and Applied Mechanics
• /
• 제3권1호
• /
• pp.66-80
• /
• 2002

The bilinear formulation proposed earlier by Peters and Izadpanah to develop finite elements in time to solve undamped linear systems, Is extended (and found to be readily amenable) to develop time finite elements to obtain transient responses of both linear and nonlinear, and damped and undamped systems. The formulation Is used in the h-, p- and hp-versions. The resulting linear and nonlinear algebraic equations are differentiated to obtain the first- and second-order sensitivities of the transient response with respect to various system parameters. The present developments were tested on a series of linear and nonlinear examples and were found to yield, when compared with results obtained using other methods, excellent results for both the transient response and Its sensitivity to system parameters. Mostly. the results were obtained using the Legendre polynomials as basis functions, though. in some cases other orthogonal polynomials namely. the Hermite. the Chebyshev, and integrated Legendre polynomials were also employed (but to no great advantage). A key advantage of the time finite element method, and the one often overlooked in its past applications, is the ease In which the sensitivity of the transient response with respect to various system parameters can be obtained. The results of sensitivity analysis can be used for approximate schemes for efficient solution of design optimization problems. Also. the results can be applied to gradient-based parameter identification schemes.

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

Advection-dominated transport problems possess difficulties in the design of numerical methods for solving them. Because of the hyperbolic nature of advective transport, many characteristic numerical methods have been developed such as the classical characteristic method, the Eulerian-Lagrangian method, the transport diffusion method, the modified method of characteristics, the operator splitting method, the Eulerian-Lagrangian localized adjoint method, the characteristic mixed method, and the Eulerian-Lagrangian mixed discontinuous method. In this paper relationships among these characteristic methods are examined. In particular, we show that these sometimes diverse methods can be given a unified formulation. This paper focuses on characteristic finite element methods. Similar examination can be presented for characteristic finite difference methods.

• 전산구조공학
• /
• 제7권4호
• /
• pp.57-67
• /
• 1994

지하구조물의 주위지반은 일반적으로 퇴적층의 형성 또는 지각의 변동에 의해 다층구조를 가지게 되므로, 구조물 및 주위지반의 거동을 정확히 예측하기 위해서는 해석에 다층구조의 영향을 반영해야 한다. 본 연구에서는 다층으로 구성된 지하구조계를 대상으로 하여 구조물과 그 주변에는 비선형 유한요소를 사용하고, 비선형성이 상대적으로 미약한 주변 다층지반에는 선형 경계요소를 사용하여 재료의 비선형성과 비균질성을 고려한 효율적인 조합해석방법을 개발하고자 한다. 반무한영역에 설정되는 다층구조계를 경계요소로 해석할 경우 그 기본해가 제한되어 있으므로, 본 연구에서는 기존의 무한기본해를 이용하는 방법을 사용하였다. 무한기본해를 이용하는 내부영역문제의 경우 각각의 균질한 층을 부영역(subdomain)으로 분할하고 계방정식을 구성한 뒤에 접합면에 대하여 평형조건과 적합조건을 만족시켜 주는 방법을 사용하여 비균질성을 고려한다. 부영역으로 층을 분할한 내부영역문제의 경계요소해석 결과는 선형 유한요소해석 결과와 비교하여 검증하였고, 검증된 경계요소 프로그램을 비선형 유한요소 프로그램과 조합하였다. 조합해석 결과, 굴착부 주변의 응력집 중부에는 비선형 유한요소를 사용하고, 비선형의 영향이 미소한 주변의 다층지반에 대해서는 부영역에 의한 선형 경계요소를 사용하는 조합해석방법이 합리적이고 효율적임을 알 수 있었다.

• Structural Engineering and Mechanics
• /
• 제15권5호
• /
• pp.579-608
• /
• 2003

This paper deals with adaptive finite element analysis of linearly elastic structures using different error estimators based on flux projection (or best guess stress values) and residual methods. Presentations are given on a typical h-type adaptive analysis, a mesh refinement scheme and the coupling of adaptive finite element analysis with automatic mesh generation. Details about different error estimators are provided and their performance, reliability and convergence are studied using six node quadratic triangular elements. Several examples are presented to demonstrate the reliability of different error estimators.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제27권4호
• /
• pp.245-271
• /
• 2023

In this paper we discuss some recovery of H(div)-conforming flux approximations from the equilibrated fluxes of Ainsworth and Oden for quadratic finite element methods of second-order elliptic problems. Combined with the hypercircle method of Prager and Synge, these flux approximations lead to a posteriori error estimators which provide guaranteed upper bounds on the numerical error. Furthermore, we prove some superconvergence results for the flux approximations and asymptotic exactness for the error estimator under proper conditions on the triangulation and the exact solution. The results extend those of the previous paper for linear finite element methods.

• Structural Engineering and Mechanics
• /
• 제43권4호
• /
• pp.411-437
• /
• 2012

In this paper, we propose a continuum mechanics based 3-D beam finite element with cross-sectional discretization allowing for warping displacements. The beam element is directly derived from the assemblage of 3-D solid elements, and this approach results in inherently advanced modeling capabilities of the beam element. In the beam formulation, warping is fully coupled with bending, shearing, and stretching. Consequently, the proposed beam elements can consider free and constrained warping conditions, eccentricities, curved geometries, varying sections, as well as arbitrary cross-sections (including thin/thick-walled, open/closed, and single/multi-cell cross-sections). We then study the modeling and predictive capabilities of the beam elements in twisting beam problems according to geometries, boundary conditions, and cross-sectional meshes. The results are compared with reference solutions obtained by analytical methods and solid and shell finite element models. Excellent modeling capabilities and solution accuracy of the proposed beam element are observed.

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