From 46a2ebce58213f5f797397dd481ab9e17521906d Mon Sep 17 00:00:00 2001
From: FZill <florian.zill@ufz.de>
Date: Tue, 27 Feb 2024 15:32:25 +0100
Subject: [PATCH] [propertylib] updated tests
---
tests/test_propertylib.py | 45 ++++++++++++++++++++++++++++++++-------
1 file changed, 37 insertions(+), 8 deletions(-)
diff --git a/tests/test_propertylib.py b/tests/test_propertylib.py
index e751aa2f9..c94b2565f 100644
--- a/tests/test_propertylib.py
+++ b/tests/test_propertylib.py
@@ -5,8 +5,11 @@ import unittest
import numpy as np
from pint.facets.plain import PlainQuantity
+from ogstools.meshplotlib.examples import mesh_mechanics
from ogstools.propertylib import presets as pp
+from ogstools.propertylib.matrix import Matrix
from ogstools.propertylib.property import Scalar, u_reg
+from ogstools.propertylib.vector import Vector
qty = u_reg.Quantity
@@ -77,10 +80,10 @@ class PhysicalPropertyTest(unittest.TestCase):
def test_von_mises(self):
"""Test von_mises_stress property."""
sig_3D = np.array([3, 1, 1, 1, 1, 1]) * 1e6
- self.equality(pp.von_mises_stress, sig_3D, qty(2, "MPa"))
- self.equality(pp.von_mises_stress, [sig_3D, sig_3D], qty([2, 2], "MPa"))
+ self.equality(pp.stress.von_Mises, sig_3D, qty(2, "MPa"))
+ self.equality(pp.stress.von_Mises, [sig_3D, sig_3D], qty([2, 2], "MPa"))
sig_2D = np.array([2, 1, 1, 1]) * 1e6
- self.equality(pp.von_mises_stress, sig_2D, qty(1, "MPa"))
+ self.equality(pp.stress.von_Mises, sig_2D, qty(1, "MPa"))
def test_eff_pressure(self):
"""Test effective_pressure property."""
@@ -91,13 +94,39 @@ class PhysicalPropertyTest(unittest.TestCase):
def test_qp_ratio(self):
"""Test qp_ratio property."""
sig = np.array([4, 4, 1, 1, 1, 1]) * 1e6
- self.equality(pp.qp_ratio, sig, qty(-100, "percent"))
- self.equality(pp.qp_ratio, [sig] * 2, qty([-100] * 2, "percent"))
+ self.equality(pp.stress.qp_ratio, sig, qty(-100, "percent"))
+ self.equality(pp.stress.qp_ratio, [sig] * 2, qty([-100] * 2, "percent"))
- def test_dilatancy(self):
- """Test dilatancy property."""
+ def test_integrity_criteria(self):
+ """Test integrity criteria."""
sig = np.array([4, 1, 2, 1, 1, 1]) * 1e6
- self.assertRaises(TypeError, pp.dilatancy_alkan_eff, sig)
+ # not working for arrays (only works for meshes)
+ self.assertRaises(TypeError, pp.dilatancy_alkan, sig)
+ self.assertGreater(np.max(pp.dilatancy_alkan(mesh_mechanics)), 0)
+ self.assertGreater(np.max(pp.dilatancy_alkan_eff(mesh_mechanics)), 0)
+ self.assertGreater(np.max(pp.dilatancy_critescu(mesh_mechanics)), 0)
+ self.assertGreater(np.max(pp.dilatancy_critescu_eff(mesh_mechanics)), 0)
+ self.assertGreater(np.max(pp.fluid_pressure_crit(mesh_mechanics)), 0)
+
+ def test_tensor_attributes(self):
+ """Test that the access of tensor attributes works."""
+ # data needs to be a 2D array
+ sig = np.asarray([[4, 1, 2, 1, 1, 1]]) * 1e6
+ self.assertTrue(np.all(pp.stress.eigenvalues(sig) >= 0))
+ for i in range(3):
+ np.testing.assert_allclose(
+ pp.stress.eigenvectors[i].magnitude(sig), 1.0
+ )
+ self.assertGreater(pp.stress.det(sig), 0)
+ self.assertGreater(pp.stress.I1(sig), 0)
+ self.assertGreater(pp.stress.I2(sig), 0)
+ self.assertGreater(pp.stress.I3(sig), 0)
+ self.assertGreater(pp.stress.mean(sig), 0)
+ self.assertGreater(pp.stress.deviator.magnitude(sig), 0)
+ self.assertGreater(pp.stress.J1(sig), 0)
+ self.assertGreater(pp.stress.J2(sig), 0)
+ self.assertGreater(pp.stress.J3(sig), 0)
+ self.assertGreater(pp.stress.octahedral_shear(sig), 0)
def test_simple(self):
"""Test call functionality."""
--
GitLab