• Structural Engineering and Mechanics
• /
• v.30 no.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.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.04a
• /
• 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
• /
• v.56 no.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
• /
• v.3 no.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
• /
• v.6 no.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.

• Computational Structural Engineering
• /
• v.7 no.4
• /
• pp.57-67
• /
• 1994

It is usual that underground structures are constructed within a multi-layered medium. In this paper, an efficient numerical modelling of multi-layered structural systems is studied using coupled analysis of finite elements and boundary elements. The finite elements are applied to the area in which the material nonlinearity dominates, and the boundary elements are applied to the far field where the nonlinearity is relatively weak. In the boundary element modelling of the multi-layered medium, fundamental solutions are not readily available. Thus, methods which can utilize existing Kelvin solutions are sought for the interior multi-layered domain problem. The interior domain problem which has piecewise homogeneous layers is analyzed using boundary elements with Kelvin solution, by discretizing each homogeneous subdomain and enforcing compatibility and equilibrium conditions between interfaces. Developed methodology is verified by comparing its results with those from the finite element analysis and it is concluded that coupled analysis using boundary elements and finite elements can be reasonable and efficient.

• Structural Engineering and Mechanics
• /
• v.15 no.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
• /
• v.27 no.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
• /
• v.43 no.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
• /
• 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.