Elastoplastic Modeling Using The Successive Linear Approxitamtion Method (original) (raw)
Related papers
Successive linear approximation for finite elasticity
2010
A method of successive Lagrangian formulation of linear approximation for solving boundary value problems of large deformation in finite elasticity is proposed. Instead of solving the nonlinear problem, by assuming time steps small enough and the reference configuration updated at every step, we can linearize the constitutive equation and reduce it to linear boundary value problems to be solved successively with incremental boundary data. Moreover, nearly incompressible elastic body is considered as an approximation to account for the condition of incompressibility. For the proposed method, numerical computations of pure shear of a square block for Mooney-Rivlin material are considered and the results are compared with the exact solutions.
An Implicit Numerical Integration Algorithm for Bai & Wierzbicki (2007) Elasto-Plastic Model
This contribution describes an implicit algorithm for numerical integration of a recently proposed model for metal plasticity and fracture [2]. The constitutive equations of the material model critically include both the effect of pressure through the triaxiality ratio and the effect of third deviatoric stress invariant through the lode angle in the description of material. These effects are directly introduced on the hardening rule of the material. The theoretical basis of the material model is presented in the first part of the paper. Then, the necessary steps required to implement the model within an implicit quasi-static finite element environment are discussed. In particular, the stress update procedure, which is based on the so-called operator split concept resulting in the standard elastic predictor/return mapping algorithm, and the computation of tangent matrix consistent with the stress update are described. Finally, the simulation of a flat grooved specimen subjected to tension [1] is presented to illustrate the robustness and efficiency of the proposed algorithm.
Numerical integration of elasto-plastic constitutive models using the extrapolation method
Proceedings of 10th International Symposium on Numerical Models in Geomechanics (Numerical Models in Geomechanics, NUMOG X): 211-217. Rhodes, Greece, 25 - 27 April 2007, Taylor & Francis, ISBN: 978-0-415-44027-1, 2007
Stress-strain integration algorithms are an important component of finite element codes. When they are robust, accurate and fast, the performance of a finite element code significantly improves, especially when advanced elasto-plastic constitutive models are used. This paper introduces a novel algorithm for the stressstrain integration of elasto-plastic soil models with automatic control of the integration error. The proposed algorithm is based on the extrapolation method used for the solution of ordinary differential equations. In this work the algorithm has been coded for the Barcelona Basic Model (a classic elasto-plastic volumetric hardening constitutive model for unsaturated soils) but its application can be easily extended to other categories of models. The efficiency and error properties of the extrapolation algorithm are assessed by integrating stresses over strain increments of different sizes and with different error tolerances. The performance of the algorithm is compared against alternative Runge-Kutta integration schemes with control of integration error.
New Numerical Solver for Elastoplastic Problems based on the Moreau-Yosida Theorem
2006
We discuss a new solution algorithm for solving elastoplastic problems with hardening. The one time-step elastoplastic problem can be formulated as a convex minimization problem with a continuous but non-smooth functional depending on unknown displacement smoothly and on the plastic strain nonsmoothly. It is shown that the functional structure allows the application of the Moreau-Yosida Theorem known in convex analysis. It guarantees that the substitution of the non-smooth plastic-strain as a function of the linear strain which depends on the displacement only yields an already smooth functional in the displacement only. Moreover, the second derivative of such functional exists in all continuum points apart from interfaces where elastic and plastic zones intersect. This allows the efficient implementation of a Newton-like method. For easy implementation most essential Matlab c © functions are provided. Numerical experiments in two dimensions state quadratic convergence of the Newton...
Reduction of the number of material parameters by ANN approximation
Computational Mechanics, 2012
Modern industrial standards require advanced constitutive modeling to obtain satisfactory numerical results. This approach however, is causing significant increase in number of material parameters which can not be easily obtained from standard and commonly known experimental techniques. Therefore, it is desirable to introduce procedure decreasing the number of the material parameters. This reduction however, should not lead to misunderstanding the fundamental physical phenomena. This paper proposes the reduction of the number of material parameters by using ANN approximation. Recently proposed viscoplasticity formulation for anisotropic solids (metals) developed by authors is used as an illustrative example. In this model one needs to identify 28 material parameters to handle particular metal behaviour under adiabatic conditions as reported by Glema et al (
International Journal of Mechanical Sciences, 2013
This paper presents an extension of the classical elastic law. The main objective of this new law is to represent linear and non-linear behaviour for computational metal forming purposes. The extension of the model is built by means of a stress-strain relationship given by an integral equation, its kernel characterising the mentioned complex behaviour. A specific application is presented for a TRIP 700 steel. Using experimental data obtained from elastoplasticity tests, the kernel of the model is formulated by means of a specific computational curve fitting procedure. The excellent agreement between experimental data and fitted model results validates the proposed model. Furthermore, a V-bending test is simulated, springback being represented by means of three different models: elastic, linear elastic with variable elastic modulus and extended elastic models. From the results differences between extended elastic law and the classical one are around 50%.
Proceedings of the XXXVIII Iberian Latin American Congress on Computational Methods in Engineering
The implementation of high-order (spectral) approximations associated with FEM is an approach to overcome the difficulties encountered in the numerical analysis of complex problems. This paper proposes the use of the spectral finite element method, originally developed for computational fluid dynamics problems, to achieve improved solutions for these types of problems. Here, the interpolation nodes are positioned in the zeros of orthogonal polynomials (Legendre, Lobatto, or Chebychev) or equally spaced nodal bases. A comparative study between the bases in the recovery of solutions to 1D and 2D elastostatic problems are performed. Examples are evaluated, and a significant improvement is observed when the SFEM, particularly the Lobatto approach, is used in comparison to the equidistant base interpolation.