The first part of the mathematical model proposed is based on a saturation model, where β is defined as a positive material constant. The rate dependency is described through the flow stressσf, which is supposed to be given by the following expression