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

• Bulletin of the Korean Mathematical Society
• /
• v.29 no.2
• /
• pp.277-283
• /
• 1992

Let GF(q) be a finite field with q elements where q=p$^{n}$ for a prime number p and a positive integer n. Consider an arbitrary function .phi. from GF(q) into GF(q). By using the Largrange's Interpolation formula for the given function .phi., .phi. can be represented by a polynomial which is congruent (mod x$^{q}$ -x) to a unique polynomial over GF(q) with the degree < q. In [3], Wells characterized all polynomial over a finite field which commute with translations. Mullen [2] generalized the characterization to linear polynomials over the finite fields, i.e., he characterized all polynomials f(x) over GF(q) for which deg(f) < q and f(bx+a)=b.f(x) + a for fixed elements a and b of GF(q) with a.neq.0. From those papers, a natural question (though difficult to answer to ask is: what are the explicit form of f(x) with zero terms\ulcorner In this paper we obtain the exact form (together with zero terms) of a polynomial f(x) over GF(q) for which satisfies deg(f) < p$^{2}$ and (1) f(x+a)=f(x)+c for the fixed nonzero elements a and c in GF(q).

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.405-410
• /
• 2007

We present a new global/local analysis with the aid of MLS(Moving Least Square)-based finite elements which can handle an arbitrary number of nodes on every element side. It give a great flexibility in constructing finite element meshes at the specified local regions without remeshing. Compared to other type global/local analysis, it does not require any superimposed mesh or need not solve the equilibrium equation twice as well as shows an excellent accuracy. To demonstrate the performance of proposed scheme, we will show several examples in relation to capturing highly local stress field.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.10a
• /
• pp.444-451
• /
• 2003

In this study, a new three-dimensional finite element analysis model of high-speed railway bridges considering train-bridge interaction, in which various improved finite elements are used for modeling structural members, is proposed. The box-type bridge deck of a railway bridge is modeled by the NFS(Nonconforming Flat Shell) elements with 6 degrees of freedom. Track structures are idealized using the beam finite elements with the offset of beam nodes and those on Winkler foundation with two parameters. And, the vehicle model devised for a high-speed train is employed, which has an articulated bogie system. By Lagrange's equations of motion, the equations of motion of a bridge-train system can be formulated. Finally, by deriving the equations of the forces acting on a bridge considering bridge-train interaction the complete system matrices of total bridge-train system can be constructed. As numerical examples of this study, 2-span PC box-girder bridge is analyzed and results are compared with experimental results.

• 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.73 no.1
• /
• pp.1-16
• /
• 2020

Finite elements based on the partition of unity (PU) approximation have powerful capabilities for p-adaptivity and solutions with high smoothness without remeshing of the domain. Recently, the PU approximation was successfully applied to the three-node shell finite element, properly eliminating transverse shear locking and showing excellent convergence properties and solution accuracy. However, the enrichment with the PU approximation results in a significant increase in the number of degrees of freedom; therefore, it requires greater computational cost, thus making it less suitable for practical engineering. To circumvent this disadvantage, we propose a new strategy to decrease the total number of degrees of freedom in the existing PU-based shell element, without loss of optimal convergence and accuracy. To alleviate the locking phenomenon, we use the method of mixed interpolation of tensorial components and perform convergence studies to show the accuracy and capability of the proposed shell element. The excellent performances of the new shell elements are illustrated in three benchmark problems.

• Structural Engineering and Mechanics
• /
• v.5 no.2
• /
• pp.137-148
• /
• 1997

In computing eigenvalues for a large finite element system it has been observed that the eigenvalue extractors produce eigenvectors that are in some sense more accurate than their corresponding eigenvalues. From this observation the paper uses a patch type technique based on the eigenvector for one mesh quality to provide an eigenvalue error indicator. Tests show this indicator to be both accurate and reliable. This technique was first observed by Stephen and Steven for an error estimation for buckling and natural frequency of beams and two dimensional in-plane and out-of-plane structures. This paper produces and error indicator for the more difficult problem of three dimensional brick elements.

• Proceeding of EDISON Challenge
• /
• 2015.03a
• /
• pp.234-238
• /
• 2015

When we use a Finite Elements Method (FEM) to solve a linear static analysis problem, number of elements need to be sufficiently small for convergence of the solution. If we analysis a part, whose curvature is varying heavily, we face to determine how small the elements size is, because the calculated stress is increased as the elements are smaller. In this case, we need to analysis with mesh insensitive method, stress linearization. We can get a solution that is not varying with the elements size if the size is smaller than a certain level. In this paper, we evaluate a pressure vessel having geometrical discontinuities using stress linearization. First, we analysis the vessel with global model, including all part of the vessel, using large shell elements. Second, we analysis the local part of the vessel, which is the small part occurring maximum stress, using small continuum elements. Last, we evaluate the safety of the pressure vessel according to the ASME Sec. VIII Div 2.

• Structural Engineering and Mechanics
• /
• v.10 no.2
• /
• pp.93-110
• /
• 2000

An alternative interpretation of the completeness requirements for the higher order elements is presented. Apart from the familiar condition, $\sum_iN_i=1$, some additional conditions to be satisfied by the shape functions of higher order elements are identified. Elements with their geometry in the natural form, i.e., without any geometrical distortion, satisfy most of these additional conditions inherently. However, the geometrically distorted elements satisfy only fewer conditions. The practical implications of the satisfaction or non-satisfaction of these additional conditions are investigated with respect to a 3-node bar element, and 8- and 9-node quadrilateral elements. The results suggest that non-satisfaction of these additional conditions results in poorer performance of the element when the element is geometrically distorted. Based on the new interpretation of completeness requirements, a 3-node element and an 8-node rectangular element that are insensitive to mid-node distortion under a quadratic displacement field have been developed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.183-189
• /
• 2014

In this paper, a polyhedral element is presented to solve three-dimensional problems by developing shape functions based on Wachspress coordinates and moving least square approximation. A subdivision of polyhedrons into tetrahedral domains is performed for the construction of shape functions of polyhedral elements, and numerical integration of the weak form is carried out consistently over the tetrahedral domains. The weight functions for moving least square approximation are defined by solving Laplace equation with boundary values based on Wachspress coordinates on polyhedral element faces. Polyhedral elements presented in this paper have similar properties to conventional finite element regarding the continuity, the completeness, the node-element connectivity and the inter-element compatibility. Numerical examples show the effectiveness of the present method for solving three-dimensional problems using polyhedral elements.