def von_mises_strain(self):
        """
        Equivalent strain to Von Mises Stress
        """
        eps = self - 1/3 * np.trace(self) * np.identity(3)

        return np.sqrt(np.sum(eps * eps) * 2/3)
Equivalent strain to Von Mises Stress