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

• Structural Engineering and Mechanics
• /
• v.58 no.3
• /
• pp.515-532
• /
• 2016

The elastoplastic response of functionally graded material (FGM) beams resting on a nonlinear elastic foundation to an eccentric axial load is investigated by using the finite element method. The FGM is assumed to be formed from ceramic and metal phases with their volume fraction vary in the thickness direction by a power-law function. A bilinear elastoplastic behavior is assumed for the metallic phase, and the effective elastoplastic properties of the FGM are evaluated by Tamura-Tomota-Ozawa (TTO) model. Based on the classical beam theory, a nonlinear finite beam element taking the shift in the neutral axis position into account is formulated and employed in the investigation. An incremental-iterative procedure in combination with the arc-length control method is employed in computing the equilibrium paths of the beams. The validation of the formulated element is confirmed by comparing the equilibrium paths obtained by using the present element and the one available in the literature. The numerical results show that the elastoplastic post-buckling of the FGM beams is unstable, and the post-buckling strength is higher for the beams associated with a higher ceramic content. Different from homogeneous beams, yielding in the FGM beam occurs in the layer near the ceramic layer before in the layer near metal surface. A parametric study is carried out to highlight the effect of the material distribution, foundation support and eccentric ratio on the elastoplastic response of the beams.

• Advances in Computational Design
• /
• v.5 no.3
• /
• pp.305-321
• /
• 2020

In this paper is presented the solution method for three-dimensional problem of transversely isotropic body's elastoplastic deformation by the finite element method (FEM). The process of problem solution consists of: determining the effective parameters of a transversely isotropic medium; construction of the finite element mesh of the body configuration, including the determination of the local minimum value of the tape width of non-zero coefficients of equation systems by using of front method; constructing of the stiffness matrix coefficients and load vector node components of the equation for an individual finite element's state according to the theory of small elastoplastic deformations for a transversely isotropic medium; the formation of a resolving symmetric-tape system of equations by summing of all state equations coefficients summing of all finite elements; solution of the system of symmetric-tape equations systems by means of the square root method; calculation of the body's elastoplastic stress-strain state by performing the iterative process of the initial stress method. For each problem solution stage, effective computational algorithms have been developed that reduce computational operations number by modifying existing solution methods and taking into account the matrix coefficients structure. As an example it is given, the problem solution of fibrous composite straining in the form of a rectangle with a system of circular holes.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.2
• /
• pp.41-49
• /
• 1993

This paper presents a geometrically non-linear and elastoplastic F.E. formulation using a total Lagrangian approach for the two dimensional isoparametric curved beam elements. The beam element is derived by using plane stress elements. The basic element geometry is constructed using the coordinates of the nodes on the element center line and the nodal point normals. The element displacement field is described using two translations of the node on the center line and a rotation about the axes normal to the plane containing the center line of the element. The layered approach is used for the elastoplastic analysis of the plane framed structure with the arbitrary cross section. The iterative load or displacement incremental method for non-linear finite element analysis of the frame structure is used. Numerical examples are presented to demonstrate the behavior and the accuracy of the proposed beam element for geometric and elastoplastic non-linear applications. Comparisons made with present theory and other published data show that tilt' beam element products accurate results with good convergence characteristics.

• Transactions of Materials Processing
• /
• v.28 no.3
• /
• pp.115-122
• /
• 2019

In this paper, finite element analysis of an automatic five-stage precision cold forging process of a yoke, a steering part of a passenger's car, is conducted with emphasis on spring back analysis at the yoke-forming stage and its experimental verification is subsequently made. An elastoplastic finite element method with MINI-element technique employed for the analysis of the entire process is explained. There is emphasis that the thin film of material formed between the punch and die in the stage may result to some errors especially in elastoplastic finite element analysis of spring back due to frequent remeshing. The numerical robustness of the spring back analysis in regards to remeshing is hence shown first through investigation into its effect on the predicted spring back. Experimental measurement of displacement due to spring back is carried out for comparison with the predicted results, and they are in a qualitative agreement with each other.

• Structural Engineering and Mechanics
• /
• v.11 no.1
• /
• pp.49-70
• /
• 2001

This paper is concerned with a study on thermo-elastoplastic characteristics of functionally graded composite. Compared to the classical layered composites, it shows a wide range of thermo-elastoplastic characteristics according to the choice of two major parameters, the thickness-wise volume fraction of constituents and the relative thickness ratio of the graded layer. Therefore, by selecting an appropriate combination of the two parameters, one is expected to design the most suitable heat-resisting composite for a given thermal circumstance. Here, we address the parametric investigation on its characteristics together with theoretical study on thermo-elastoplasticity and numerical techniques for its finite element approximations. Through the numerical experiments, we examine the influence of two parameters on the thermo-elastoplastic characteristics.

• Structural Engineering and Mechanics
• /
• v.68 no.3
• /
• pp.299-312
• /
• 2018

In this paper, an extended finite element method is proposed to analyze both geometric and material non-linear behavior of general Functionally Graded Material (FGM) plate-shell type structures. A user defined subroutine (UMAT) is developed and implemented in Abaqus/Standard to study the elastoplastic behavior of the ceramic particle-reinforced metal-matrix FGM plates-shells. The standard quadrilateral 4-nodes shell element with three rotational and three translational degrees of freedom per node, S4, is extended in the present study, to deal with elasto-plastic analysis of geometrically non-linear FGM plate-shell structures. The elastoplastic material properties are assumed to vary smoothly through the thickness of the plate-shell type structures. The nonlinear approach is based on Mori-Tanaka model to underline micromechanics and locally determine the effective FGM properties and self-consistent method of Suquet for the homogenization of the stress-field. The elasto-plastic behavior of the ceramic/metal FGM is assumed to follow Ludwik hardening law. An incremental formulation of the elasto-plastic constitutive relation is developed to predict the tangent operator. In order to to highlight the effectiveness and the accuracy of the present finite element procedure, numerical examples of geometrically non-linear elastoplastic functionally graded plates and shells are presented. The effects of the geometrical parameters and the volume fraction index on nonlinear responses are performed.

• Transactions of Materials Processing
• /
• v.12 no.5
• /
• pp.479-486
• /
• 2003

This paper presents an elastoplastic analysis for spent nuclear fuel disposal container and its 50 cm thick bentonite buffer to predict the collapse of the container while the horizontal asymmetric sudden rock movement of 10 cm is applied on the composite structure. This sudden rock movement is anticipated by the earthquake etc. at a deep underground. Elastoplastic material model is adopted. Drucker-Prager yield criterion is used for the material yield prediction of the bentonite buffer and von-Mises yield criterion is used for the material yield prediction of the container. Analysis results show that even though very large deformations occur beyond the yield point in the bentonite buffer, the container structure still endures elastic small strains and stresses below the yield strength. Hence, the asymmetric 50 cm thick bentonite buffer can protect the container safely against the 10 cm sudden rock movement by earthquake etc.. Analysis results also show that bending deformations occur in the container structure due to the shear deformation of the bentonite buffer. The finite element analysis code, NISA, is used for the analysis.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.30 no.1
• /
• pp.1-6
• /
• 2017

This paper presents development of the nine-node co-rotational(CR) planar element applicable for elastoplastic and contact analysis. The CR formulation is one of the efficient geometrically nonlinear formulations. It is based on the assumptions of small strain and large displacement. Further, it is extended to both elastoplastic analysis and contact analysis in this paper. For accurate plastic analysis, nine-node quadrilateral element, which can provide accurate stress prediction, is chosen. Bi-linear hardening rule based on Newton- Raphson return-mapping is employed. Also, Lagrange multiplier is used in order for constraints regarding the contact analysis. The present development is validated via the time transient problems. The present results are compared with those obtained by the other existing software.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2006.05a
• /
• pp.386-389
• /
• 2006

Prestressed composite girder bridges with concrete infilled steel tube at negative flexural moment region takes the advantages of enhancing local buckling and flexural resistances resulting from the lateral confining effect of concrete due to the interactive reaction in the interface layer of steel tube and concrete. The interface behavior in concrete infilled tube of the test composite girder is analyzed by 8-node zero thickness interface finite element combined with 3-D. elastoplastic concrete constitutive model and 3-D. elastoplastic Mindlin shell element. The interface effects between infillled concrete and steel tube are investigated through the comparision of the experimental and numerical results.