• Structural Engineering and Mechanics
• /
• v.3 no.2
• /
• pp.185-192
• /
• 1995

A numerical algorithm for plane stress and shell elasto-plasticity is presented in this paper. The proposed strain decomposition (SD) algorithm is an elastic predictor/plastic corrector algorithm, and in the context of operator splitting, is a return mapping algorithm. However, it differs significantly from other return mapping algorithms in that only the necessary response functions are used without invoking their gradients, and the stress increment is updated only at the end of the time step. This makes the proposed SD algorithm more suitable for materials with complex yield surfaces and will guard against error accumulation during the time step. Comparative analyses of structural systems using the proposed strain decomposition (SD) algorithm and the iterative radial return (IRR) algorithm are presented. The results demonstrate the accuracy and usefulness of the proposed algorithm.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.2
• /
• pp.197-208
• /
• 2000

Choi et al./sup 1)/ presented the total Lagrangian formulation based upon the degenerated shell element. Geometrically correct formulation is developed by updating the direction of normal vectors and taking into account the second order rotation terms in the incremental displacement field. Assumed strain concept is adopted in order to overcome the shear locking phenomena and to eliminate the spurious zero energy mode. In this paper, for the ultimate strength analysis of stiffened shell structures considering effects of residual stresses, the return mapping algorithm based on the consistent elasto-plastic tangent modulus is applied to anisotropic shell structures. In addition, the load/displacement incremental scheme is adopted for non-linear F.E. analysis. Based on such methodology, the computer program is developed and numerical examples to demonstrate the accuracy and the effectiveness of the proposed shell element are presented and compared with the results in literatures.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.6
• /
• pp.1455-1464
• /
• 1994

A new and efficient algorithm for the integration of the thermal-elasto-plastic constitutive equation is proposed. While it falls into the category of the return mapping method, the algorithm adopts the three point approximation of plastic corrector within one time increment step. The results of its application to a von Mises-type thermal-elasto-plastic model with combined hardening and temperature-dependent material properties show that the accurate iso-error maps are obtained for both angular and radial errors. The accuracy achieved is because the predicted stress increment in a single step calculation follows the exact value closely not only at the end of the step but also through the whole path. Also, the comparison of the computational time for the new and other algorithms shows that the new one is very efficient.

• Structural Engineering and Mechanics
• /
• v.75 no.1
• /
• pp.101-108
• /
• 2020

The present study intends to analyze damage in thin-walled steel cylinders undergoing constant internal pressure and thermal cycles through use of anisotropic continuum damage mechanics (CDM) model coupled with nonlinear kinematic hardening rule of Chaboche. Materials damage in each direction was defined based on plastic strain and its direction. Stress and strain distribution over wall-thickness was described based on the CDM model and the return mapping algorithm was employed based on the consistency condition. Plastic zone expansion across the wall thickness of cylinders was noticeably affected with change in internal pressure and temperature gradients. Expansion of plastic zone over wall-thickness at inner and outer surfaces and their boundaries demarking elastic and plastic regions was attributed to the magnitude of damage induced over thermomechanical cycles on the thin-walled samples tested at various pressure stresses.

• Transactions of Materials Processing
• /
• v.7 no.3
• /
• pp.233-245
• /
• 1998

A new approach has been proposed for the incremental analysis of the nonsteady state large deformation of planar anisotropic elastic-plastic sheet forming. A mathematical brief review of a constitutive law for the incremental deformation theory has been presented from flow theory using the minimum plastic work path for elastic-plastic material. Since the material embedded coordinate system(Lagrangian quantity) is used in the proposed theory the stress integration procedure is completely objective. A new return mapping algorithm has been also developed from the general midpoint rule so as to achieve numerically large strain increment by successive control of yield function residuals. Some numerical tests for the return mapping algorithm were performed using Barlat's six component anisotropic stress potential. Performance of the proposed algorithm was shown to be good and stable for a large strain increment, For planar anisotropic sheet forming updating algorithm of planar anisotropic axes has been newly proposed. In order to show the effectiveness and validity of the present formulation earing simulation for a cylindrical cup drawing and front fender stamping analysis are performed. From the results it has been shown that the present formulation can provide a good basis for analysis for analysis of elastic-plastic sheet metal forming processes.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.9 no.6 s.46
• /
• pp.59-65
• /
• 2005

In this paper a constitutive model for large cracks in concrete and other brittle materials subject to dynamic and cyclic leading is presented. The suggested model is based on the plastic-damage model for cyclic leading. A numerical formulation based on the three-step return-mapping algorithm for the proposed large crack model is also present. The numerical examples show that the present algorithm works appropriately under dynamic leading and should be used in large crack problems to prevent excessive tensive plastic strain development causing unrealistic results.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.5 no.3
• /
• pp.57-64
• /
• 2001

For the dynamic nonlinear analysis of stiffened plate and shell structures, total Lagrangian formulation is presented based upon the degenerated shell element considering finite rotation effects. Assumed strain concept is adopted in order to overcome shear locking phenomena and to eliminate spurious zero energy mode. In the elasto-plastic analysis, the return mapping algorithm based on the consistent elasto-plastic tangent modulus is applied to collapse analysis of shell structures. Newmark integration method is used for dynamic nonlinear analysis of shell structures under dynamic forces.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.483-488
• /
• 2002

In this paper numerical algorithm for the continuum large crack model is proposed based on the return-mapping formulation. The numerical test results show that the present algorithm works appropriately under cyclic loading. It is also shown that in continuum damage models a large crack model to prevent excessive tensile plastic strain should be used to have realistic cyclic loading simulation results.

• Computational Structural Engineering
• /
• v.10 no.3
• /
• pp.241-250
• /
• 1997

The strain localization of concrete is a phenomenon such that the deformation of concrete is localized in finite region along with softening behavior. The objective of this paper is to develop a plasticity and damage algorithm for the finite element analysis of the strain-localization in concrete. In this paper, concrete member under strain localization is modeled with localized zone and non-localized zone. For modeling of the localized zone in concrete under strain localization, a general Drucker-Prager failure criterion by which the nonlinear strain softening behavior of concrete after peak-stress can be considered is introduced in a thermodynamic formulation of the classical plasticity model. The return-mapping algorithm is used for the integration of the elasto-plastic rate equation and the consistent tangent modulus is also derived. For the modeling of non-localized zone in concrete under strain localization, a consistent nonlinear elastic-damage algorithm is developed by modifying the free energy in thermodynamics. Using finite element program implemented with the developed algorithm, strain localization behaviors for concrete specimens under compression are simulated.

• Journal of the Korean Geotechnical Society
• /
• v.22 no.9
• /
• pp.45-59
• /
• 2006

The successful performance of any numerical geotechnical simulation depends on the accuracy and efficiency of the numerical implementation of constitutive model used to simulate the stress-strain (constitutive) response of the soil. The corner stone of the numerical implementation of constitutive models is the numerical integration of the incremental form of soil-plasticity constitutive equations over a discrete sequence of time steps. In this paper a well known two-surface soil plasticity model is implemented using a generalized implicit return mapping algorithm to arbitrary convex yield surfaces referred to as the Closest-Point-Projection method (CPPM). The two-surface model describes the nonlinear behavior of coarse-grained materials by incorporating a bounding surface concept together with isotropic and kinematic hardening as well as fabric formulation to account for the effect of fabric formation on the unloading response. In the course of investigating the performance of the CPPM integration method, it is proven that the algorithm is an accurate, robust, and efficient integration technique useful in finite element contexts. It is also shown that the algorithm produces a consistent tangent operator $\frac{d\sigma}{d\varepsilon}$ during the iterative process with quadratic convergence rate of the global iteration process.