def test_Elastic(self):
K = 12.3
G = 45.6
gamma = 0.02
epsm = 0.12
Eps = np.array([[epsm, gamma, 0.0], [gamma, epsm, 0.0],
[0.0, 0.0, epsm]])
Sig = np.array([[3.0 * K * epsm, 2.0 * G * gamma, 0.0],
[2.0 * G * gamma, 3.0 * K * epsm, 0.0],
[0.0, 0.0, 3.0 * K * epsm]])
self.assertTrue(np.isclose(float(GMat.Epsd(Eps)), gamma))
mat = GMat.Elastic(K, G)
mat.setStrain(Eps)
self.assertTrue(np.allclose(mat.Stress(), Sig))
def test_Smooth(self):
K = 12.3
G = 45.6
gamma = 0.02
epsm = 0.12
Eps = np.array([[epsm, gamma, 0.0], [gamma, epsm, 0.0],
[0.0, 0.0, epsm]])
Sig = np.array([[3.0 * K * epsm, 0.0, 0.0], [0.0, 3.0 * K * epsm, 0.0],
[0.0, 0.0, 3.0 * K * epsm]])
self.assertTrue(np.isclose(float(GMat.Epsd(Eps)), gamma))
mat = GMat.Smooth(K, G, [0.01, 0.03, 0.10])
mat.setStrain(Eps)
self.assertTrue(np.allclose(mat.Stress(), Sig))
self.assertTrue(np.isclose(mat.epsp(), 0.02))
self.assertTrue(mat.currentIndex() == 1)