hardening max(0,stress_strain_curve(solid.epe+solid.mises/solid.E)-solid.sigmags) "Hardening function"
deltaEeqap -T*R_const/(zm*F_const)*log((nu*alpha*solid.epe)/N0+1) "Change in anode equilibrium potential due plastic deformation"
Eeqa Eeq0a+deltaEeqae+deltaEeqap "Anode equilibrium potential including elastic and plastic deformation terms"
ic i0c*10^(-solid.mises*Vm/(6*F_const*bc)) "Cathode exchange current density including stress factor"