More robust NortonOverstress implementation

KnutAM · Feb 13, 2024 · 7184c8a · 7184c8a
1 parent 98a2752
commit 7184c8a
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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])