6.2 Perzyna Model
319
φ
n
:= |σ
n
vp | − σ y = φ
− E λ.
(6.99)
Here the trial value of the overstress function φ
has been defined as
φ
:= |σ
vp | − σ y .
(6.100)
Next for visco-plastic loading with λ > 0 the definition for the incremental
visco-plastic multiplier is regrouped to render
φ
n
− λ
η
t n = 0.
(6.101)
Thus the incremental visco-plastic multiplier λ ≥ 0 is computed in closed form
from
λ =
φ
E + η//t n ≥ 0.
(6.102)
Observe that λ degenerates to the plastic case for η → 0, likewise λ degenerates to zero in the limit of very fast processes with t
n
→ 0. Once λ is computed
all other variables may be updated. In particular, the visco-plastic stress at the end
of the time step reads
σ
n
vp = σ
vp − E λ
σ
vp
|σ
vp |
.
(6.103)
The sensitivity of σ
n
vp = σ
n with respect to
n is denoted the algorithmic tangent
E a (thus dσ = E a d) and is computed from the product rule while noting that λ
depends implicitly on
n
∂ σ
n
vp = E − E λ ∂
σ
vp
|σ
vp |
− E
σ
vp
|σ
vp |
∂ ((λ).
(6.104)
The first derivative term on the right-hand-side computes to zero since
∂
σ
vp
|σ
vp |
=
1
|σ
vp |
E −
σ
vp
|σ
vp | 2
σ
vp
|σ
vp |
E ≡ 0.
(6.105)
It shall be noted that the corresponding tangent modulus (tensor) in more than
one dimension is different from zero. The second derivative term on the right-handside computes from requiring satisfaction of ∂ [φ
n
− λ η//t
n
] = 0 for ongoing
visco-plastic flow at the end of the time step, i.e.
∂ φ
−
E +
η
t n
∂ ((λ) =
σ
vp
|σ
vp |
E −
E +
η
t n
∂ ((λ)
.
= 0.
(6.106)