An implicit time integration scheme for inelastic constitutive equations with internal state variables. (English) Zbl 0584.73108

In a broad class of inelastic constitutive models for the deformation of metals the inelastic strain rates are functions of the current state of stress and internal state variables only. All known models are in some regions of application mathematically stiff and therefore difficult to integrate. The unconditionally stable implicit Euler rule is used for integration. It leads to a system of highly nonlinear algebraic equations which have to be solved by an iterative process. The general Newton- Raphson method, which converges under very broad conditions, requires repeated solution of the finite element system and is infeasible for large inelastic problems. But for the inelastic strains and internal state variables the Jacobian can be computed analytically and therefore the NRI can be used. For the stresses the Jacobian cannot be computed analytically and therefore the accelerated Jacobi iteration is used. A new method for computing the relaxation parameter is introduced which increases the rate of convergence significantly. The new algorithm is applied on Hart’s model. A comparison with prior computations using an approximation is made.


74S30 Other numerical methods in solid mechanics (MSC2010)
65D30 Numerical integration
74C15 Large-strain, rate-independent theories of plasticity (including nonlinear plasticity)
74C20 Large-strain, rate-dependent theories of plasticity
74-04 Software, source code, etc. for problems pertaining to mechanics of deformable solids


Full Text: DOI


[1] Hart, Trans. ASME. J. Engng. Mat. Technol. 98 pp 193– (1976)
[2] ’A physically based internal variable model for rate-dependent plasticity’, in Inelastic Behaviour of Pressure Vessels and Piping Components ( and , Eds.), American Society of Mechanical Engineers, New York, 1978, pp. 15-28.
[3] Miller, Trans. ASME, J. Engng. Mat. Technol. 98 pp 97– (1976)
[4] Miller, Trans. ASME, J. Engng. Mat. Technol. 98 pp 106– (1976)
[5] ’A candidate creep-recovery model for 2-1/4Cr-1Mo steel and its experimental implementation’, Report ORNL-TM-5110 Oak Ridge National Laboratory (1975).
[6] ’Numerical integration of some new unified plasticity-creep formulations’, Proc. 4th SMIRT Conference, paper M6/6, San Francisco (1977).
[7] Kumar, Int. Journ. Mech. Sci. 19 pp 713– (1977)
[8] The Finite Element Method, 3rd edn, McGraw-Hill, London, 1977.
[9] Boundary Element Methods in Creep and Fracture, Applied Science, London, 1982. · Zbl 0534.73070
[10] ’Thermoviscoplastic response of cylindrical structures using a state variable theory’, in Mechanical Behaviour of Metals–Proceedings of ICM3, Cambridge, England, Vol. 2 ( and , Eds.), Pergamon Press, Oxford, 1979, pp. 233-242.
[11] Gear, Comm. ACM 14 pp 176– (1971)
[12] I, Int. j. numer. methods eng. 9 pp 109– (1975)
[13] Construction of Integration Formulas for Initial Value Problems, North-Holland, Amsterdam, 1977.
[14] Hughes, Comp. Struct. 8 pp 169– (1978)
[15] Kumar, Trans, ASME, J. Engng. Mat. Technol. 102 pp 92– (1980)
[16] Private communication.
[17] Studienarbeit (unpublished), Technische Universität Braunschweig, 1981.
[18] Finite Elements of Nonlinear Continua, McGraw-Hill, New York, 1972. · Zbl 0235.73038
[19] Report to Deutsche Forschungsgemeinschaft (unpublished), Technische Hochschule Darmstadt (1983).
[20] P, Advan. Appl. Mech. 9 pp 16– (1973)
[21] ’Some considerations on the constitutive law in thermoplasticity’, General lecture on the 24 Polish Solid Mechanics Conf., Jachranka (1983).
[22] , and , ’Deformation in type 304 austenitic stainless steel’, Report EPRI-NP1276, Electric Power Research Institute, Palo Alto, U.S.A. (1979).
[23] and : ’Boundary element formulation for large strain-large deformation problems of plasticity and viscoplasticity’, in Developments in Boundary Element Methods, Vol. 3 ( and , Eds.), Applied Science, London, 1984.
[24] Kollmann, Ing.-Archiv. 54 pp 355– (1984)
[25] Rotating Thin Viscoplastic Discs, in preparation.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.