| Abstract: | 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. |
