mises¶

Purpose:

Compute stresses and plastic strains for an elasto-plastic isotropic hardening von Mises material.

Syntax:

[es, deps, st] = mises(ptype, mp, est, st)

Description:

The mises function computes updated stresses, es, plastic strain increments deps, and state variables st for an elasto-plastic isotropic hardening von Mises material.

The input variable ptype defines the type of analysis, see also hooke. The vector mp contains the material constants:

mp = [\(\small{\, E\;\nu\;h\,}\)]

where \(E\) is the modulus of elasticity, \(\nu\) is the Poisson’s ratio, and \(h\) is the plastic modulus.

The input matrix est contains trial stresses obtained by using the elastic material matrix D in plants or a similar s-function. The input vector st contains the state parameters:

st = [\(\small{\, yi\;\sigma_y\;\varepsilon_{eff}^p\,}\)]

at the beginning of the step. The scalar \(yi\) indicates whether the material behaviour is elasto-plastic (\(yi = 1\)) or elastic (\(yi = 0\)). The current yield stress is denoted by \(\sigma_y\) and the effective plastic strain by \(\varepsilon_{eff}^p\).

The output variables es and st contain updated values obtained by integration of the constitutive equations over the actual displacement step. The increments of the plastic strains are stored in the vector deps.

If es and st contain more than one row, then every row will be treated by the command.

Note:

It is not necessary to check whether the material behaviour is elastic or elasto-plastic; this test is performed by the function. The computation is based on an Euler-Backward method, i.e., the radial return method.

Only the cases ptype = 2, 3, 4 are implemented.