We propose a mathematical model for relaxation modulus and its numerical solution. The model formula is extended from sigmoidal function considering nonlinear strain hardening. Its physical meaning ...