From 7184c8a4cc2fab6a8dc0abaf9ff504df5c5093a9 Mon Sep 17 00:00:00 2001 From: Knut Andreas Meyer Date: Tue, 13 Feb 2024 19:33:57 +0100 Subject: [PATCH] More robust NortonOverstress implementation --- src/plasticity_components/OverstressFunctions.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/plasticity_components/OverstressFunctions.jl b/src/plasticity_components/OverstressFunctions.jl index 3241120..812bfef 100644 --- a/src/plasticity_components/OverstressFunctions.jl +++ b/src/plasticity_components/OverstressFunctions.jl @@ -38,7 +38,7 @@ relaxation time and overstress sensitivty. nexp::TN end -overstress_function(ratelaw::NortonOverstress, Φ, σy) = Φ<=0 ? zero(Φ) : (1/ratelaw.tstar)*(Φ/σy)^ratelaw.nexp +overstress_function(ratelaw::NortonOverstress, Φ, σy) = (1/ratelaw.tstar) * macaulay(Φ/σy)^ratelaw.nexp MMB.get_num_params(::NortonOverstress) = 2 MMB.vector2material(v::AbstractVector, ::NortonOverstress; offset=0) = NortonOverstress(v[offset+1], v[offset+2])