Deformations of the colder regions lithosphere mainly occur in frictional regime. In geodynamic models, plastic deformations are often highly localized (shear bands) and used as proxies for faults. However, capturing generation evolution shear bands models is troublesome. Indeed, mesh dependency lack convergence affect, to some extent, results models. Here we extend most common plasticity imple...