Implement Kelvin Voigt model