def test_deviatoric_Mises(self):
"""Ensure that Mises equivalent stress depends only on deviatoric part."""
x = np.random.rand(self.n, 3, 3)
full = mechanics.equivalent_stress_Mises(x)
dev = mechanics.equivalent_stress_Mises(tensor.deviatoric(x))
assert np.allclose(full, dev)
def test_Mises(self):
"""Ensure that equivalent stress is 3/2 of equivalent strain."""
x = np.random.rand(self.n, 3, 3)
assert np.allclose(
mechanics.equivalent_stress_Mises(x) /
mechanics.equivalent_strain_Mises(x), 1.5)
def test_add_Mises_stress(self, default):
default.add_stress_Cauchy('P', 'F')
default.add_equivalent_Mises('sigma')
in_memory = mechanics.equivalent_stress_Mises(default.place('sigma'))
in_file = default.place('sigma_vM')
assert np.allclose(in_memory, in_file)