Skip to content

add function to read from a csv to test helper #14486

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
Sabine-van-Dijk wants to merge 12 commits into
base: master
Choose a base branch
from
Open

add function to read from a csv to test helper #14486

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
12 commits
Select commit Hold shift + click to select a range
79cd7a5
added function get_values_from_csv to test_helper
Sabine-van-Dijk Jun 5, 2026
e891461
adjusted code for node (one key) to be applicable on the CROW case.
Sabine-van-Dijk Jun 5, 2026
7a69a7d
Merge remote-tracking branch 'origin/master' into 14470-geomechanicsa…
Sabine-van-Dijk Jun 5, 2026
985f5bf
add csv files for the expected values of triaxial
Sabine-van-Dijk Jun 5, 2026
a6ea7e6
add a function to get the values from the csv as a vector including a…
Sabine-van-Dijk Jun 5, 2026
f4f26b4
rewrite the triaxial test to use the new functions and csv's
Sabine-van-Dijk Jun 5, 2026
57b2d07
add function to read gouped csv files e.g. per phase
Sabine-van-Dijk Jun 8, 2026
0cfb729
add csv files to test the new functions on the crs test
Sabine-van-Dijk Jun 8, 2026
9f61e42
changed the order of the stresses and strains in the triaxial test
Sabine-van-Dijk Jun 8, 2026
6e6798b
Merge remote-tracking branch 'origin/master' into 14470-geomechanicsa…
Sabine-van-Dijk Jun 8, 2026
bc35ccc
Minor changes to the function names
Sabine-van-Dijk Jun 8, 2026
a7c5f07
changed the ordering of the tensor element naming
Sabine-van-Dijk Jun 11, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
129 changes: 59 additions & 70 deletions applications/GeoMechanicsApplication/tests/test_element_lab.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13,52 +13,52 @@ class KratosGeoMechanicsLabElementTests(KratosGeoUnittest.TestCase):
"""
def test_triaxial_drained(self):
"""Regression test for the triaxial experiment on a mohr coulomb model with a constant pore water pressure."""
expected_disp = [[0.0, -0.2, 0.0], [0.0527776, -0.2, 0.0], [0.0, -0.100033, 0.0], [0.0524025, -0.0996931, 0.0], [0.0, 0.0, 0.0], [0.105197, -0.2, 0.0], [0.105114, -0.100049, 0.0], [0.0524406, 0.0, 0.0], [0.104632, 0.0, 0.0]]
expected_stress = [[-99.9808, -252.622, -99.9806, 0.193199, 0.0, 0.0], [-99.8828, -252.381, -99.8827, 0.0448209, 0, 0], [-100.045, -252.78, -100.045, 0.0314238, 0, 0], [-99.9991, -252.668, -99.9991, 0.00846584, 0, 0], [-99.967, -252.589, -99.967, -0.0103576, 0, 0], [-100.064, -252.827, -100.064, 0.137162, 0, 0]]
expected_strain = [[0.104863, -0.19973, 0.104946, 0.000440186, 0.0, 0.0], [0.104939, -0.199861, 0.105022, 0.000100528, 0, 0], [0.10444, -0.199182, 0.104773, 6.24562e-05, 0, 0], [0.1055, -0.200303, 0.104922, 3.84218e-05, 0, 0], [0.104915, -0.200106, 0.105298, -4.97694e-06, 0, 0], [0.105631, -0.200818, 0.105335, 0.000287088, 0, 0]]
self._run_triaxial_regression_test('drained', 'triaxial_test_output.post.res', expected_disp, expected_stress,
expected_strain, 4, assert_all_integration_points=True)
file_path = test_helper.get_file_path(os.path.join('test_element_lab', 'test_triaxial', 'drained'))
expected_disp = test_helper.get_values_from_csv_as_vectors(
Path(file_path) / "expected_disp.csv",
["node_id"],
["disp_x", "disp_y", "disp_z"])
expected_stress = test_helper.get_values_from_csv_as_vectors(
Path(file_path) / "expected_stress.csv",
["element_id", "ip_index"],
["stress_xx", "stress_yy", "stress_zz", "stress_xy", "stress_yz", "stress_xz"])
expected_strain = test_helper.get_values_from_csv_as_vectors(
Path(file_path) / "expected_strain.csv",
["element_id", "ip_index"],
["strain_xx", "strain_yy", "strain_zz", "strain_xy", "strain_yz", "strain_xz"])
self._run_triaxial_regression_test(file_path, 'triaxial_test_output.post.res', expected_disp, expected_stress,
expected_strain, 4)

def test_triaxial_undrained(self):
"""
Regression test for the undrained triaxial experiment.
"""
expected_disp = [[0.0, -0.2, 0.0], [0.045, -0.2, 0.0], [0.0, -0.1, 0.0], [ 0.045, -0.1, 0.0],[0.0, 0.0, 0.0], [0.09, -0.2, 0.0], [0.09, -0.1, 0.0], [0.045, 0.0, 0.0], [0.09, 0.0, 0.0]]
expected_stress = [[-100.0, -4740.0, -100.0, 0.0, 0.0, 0.0], [-100.0, -4740.0, -100.0, 0.0, 0.0, 0.0]]
expected_strain = [[0.09, -0.2, 0.09, 0.0, 0.0, 0.0], [0.09, -0.2, 0.09, 0.0, 0.0, 0.0]]
self._run_triaxial_regression_test('undrained', 'triaxial_undrained_test_output.post.res', expected_disp,
file_path = test_helper.get_file_path(os.path.join('test_element_lab', 'test_triaxial', 'undrained'))
expected_disp = test_helper.get_values_from_csv_as_vectors(
Path(file_path) / "expected_disp.csv",
["node_id"],
["disp_x", "disp_y", "disp_z"])
expected_stress = test_helper.get_values_from_csv_as_vectors(
Path(file_path) / "expected_stress.csv",
["element_id", "ip_index"],
["stress_xx", "stress_yy", "stress_zz", "stress_xy", "stress_yz", "stress_xz"])
expected_strain = test_helper.get_values_from_csv_as_vectors(
Path(file_path) / "expected_strain.csv",
["element_id", "ip_index"],
["strain_xx", "strain_yy", "strain_zz", "strain_xy", "strain_yz", "strain_xz"])
self._run_triaxial_regression_test(file_path, 'triaxial_undrained_test_output.post.res', expected_disp,
expected_stress, expected_strain, 4)

def _run_triaxial_regression_test(self, stage_name, output_file_name, expected_displacement, expected_stress,
expected_strain, precision_places, assert_all_integration_points=False):
test_name = 'test_triaxial'
file_path = test_helper.get_file_path(os.path.join('test_element_lab', test_name, stage_name))
def _run_triaxial_regression_test(self, file_path, output_file_name, expected_displacement, expected_stress, expected_strain, precision_places):
test_helper.run_kratos(file_path)

# Read the output files from the simulation for comparison.
reader = GiDOutputFileReader()
result = reader.read_output_from(os.path.join(file_path, output_file_name))
time = 1.0

# Assert the displacement in all nodes in all directions.
for node_id, expected_node_displacement in enumerate(expected_displacement, start=1):
node_displacement = reader.nodal_values_at_time("DISPLACEMENT", time, result, [node_id])[0]
self.assertVectorAlmostEqual(node_displacement, expected_node_displacement, precision_places)

if assert_all_integration_points:
self.assert_integration_point_tensors(result, "CAUCHY_STRESS_TENSOR",
self._make_integration_point_tensor_entries(expected_stress, num_elements=2, num_integration_points_per_element=3),
time ,precision_places)
self.assert_integration_point_tensors(result, "ENGINEERING_STRAIN_TENSOR",
self._make_integration_point_tensor_entries(expected_strain, num_elements=2, num_integration_points_per_element=3),
time, precision_places)
else:
self.assert_integration_point_tensors(result, "CAUCHY_STRESS_TENSOR",
self._make_integration_point_tensor_entries(expected_stress, num_elements=2, num_integration_points_per_element=1),
time, precision_places)
self.assert_integration_point_tensors(result, "ENGINEERING_STRAIN_TENSOR",
self._make_integration_point_tensor_entries(expected_strain, num_elements=2, num_integration_points_per_element=1),
time, precision_places)
node_ids = sorted(expected_displacement)

@avdg81 avdg81 Jun 12, 2026

Copy link
Copy Markdown
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Here, I'm getting a bit confused. Why would the node IDs be equal to the sorted expected displacements? I would have expected something like (since the node IDs are the keys in the map with the expected displacements, right?):

Suggested change
node_ids = sorted(expected_displacement)
node_ids = sorted(expected_displacement.keys())

I can guess why the code works, but I'm not 100% sure. Perhaps iterating over a map yields the keys? If yes, then I would prefer the more explicit expression, i.e. by adding .keys(). Thank you.

I'm also wondering whether the node IDs need to be sorted or not. It seems to me that the code would also work if the node IDs are not sorted. Or am I overlooking something?

@Sabine-van-Dijk Sabine-van-Dijk Jun 12, 2026 •
edited
Loading

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

From the python sorting techniques it works like the following example:
sorted({1: 'D', 2: 'B', 3: 'B', 4: 'E', 5: 'A'})
[1, 2, 3, 4, 5]
They do not need to be sorted necessarily, but it was just a way to get the keys. I will add the .keys() for clarity. I can also do just expected_displacement.keys() without sorted. that returns a dict object then instead of a list but it also works here. Please let me know if that is even more clear.

@Sabine-van-Dijk Sabine-van-Dijk Jun 12, 2026

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Okay maybe not a real dict object but I don't know what to call it.
Just expected_displacement.keys()
returns: dict_keys([1, 2, 3, 4, 5, 6, 7, 8, 9])

expected_disp = [expected_displacement[node_id] for node_id in node_ids]
self.assert_nodal_values_at_time(result, "DISPLACEMENT", node_ids, expected_disp, time, precision_places)
self.assert_integration_point_tensors(result, "CAUCHY_STRESS_TENSOR", expected_stress, time, precision_places)
self.assert_integration_point_tensors(result, "ENGINEERING_STRAIN_TENSOR", expected_strain, time, precision_places)

def test_oedometer_drained(self):
"""Regression test for the oedometer experiment on a linear elastic model with constant pore water pressure."""
Expand Down Expand Up @@ -122,38 +122,34 @@ def test_crs_drained(self):
"""Regression test for the CRS experiment with constant pore water pressure."""
stage_name = 'drained'
nr_of_phases = 5
nr_of_elements = 2
nr_of_integration_points_per_element = 3
total_integration_points = nr_of_elements * nr_of_integration_points_per_element

expected_strains = [
self._repeat_tensor([0.0, -0.1, 0.0, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([0.0, -0.05, 0.0, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([0.0, -0.25, 0.0, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([0.0, -0.25, 0.0, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([0.0, -0.40, 0.0, 0.0, 0.0, 0.0], total_integration_points),
]
expected_water_pressures = [[0.0] * 9 for _ in range(nr_of_phases)]

expected_stresses = [
self._repeat_tensor([-4e+05, -12e+05, -4e+05, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([-2e+05, -6e+05, -2e+05, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([-1e+06, -3e+06, -1e+06, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([-1e+06, -3e+06, -1e+06, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([-1.6e+06, -4.8e+06, -1.6e+06, 0.0, 0.0, 0.0], total_integration_points),
]
linear_elastic_file_path = test_helper.get_file_path(os.path.join('test_element_lab', 'test_crs', stage_name, 'linear_elastic'))
linear_elastic_stresses = test_helper.get_values_from_csv_grouped(
Path(linear_elastic_file_path) / "expected_stress.csv",
['phase'], ['element_id', 'ip_index'],
["stress_xx", "stress_yy", "stress_zz", "stress_yz", "stress_xz", "stress_xy"])
linear_elastic_strains = test_helper.get_values_from_csv_grouped(
Path(linear_elastic_file_path) / "expected_strain.csv",
['phase'], ['element_id', 'ip_index'],
["strain_xx", "strain_yy", "strain_zz", "strain_yz", "strain_xz", "strain_xy"])
expected_stresses = [linear_elastic_stresses[phase] for phase in range(1, nr_of_phases + 1)]
expected_strains = [linear_elastic_strains[phase] for phase in range(1, nr_of_phases + 1)]

self._run_crs_test_for_model(stage_name, "linear_elastic")
self._check_crs_results(stage_name, "linear_elastic", nr_of_phases, expected_stresses, expected_strains, expected_water_pressures, delta_stress=1.0)

# mohr coulomb expectations per phase
expected_stresses = [
self._repeat_tensor([-9.01878e+05, -2.22843e+06, -9.01878e+05, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([-5.01878e+05, -1.02843e+06, -5.01878e+05, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([-2.25697e+06, -5.56725e+06, -2.25697e+06, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([-2.25697e+06, -5.56725e+06, -2.25697e+06, 0.0, 0.0, 0.0], total_integration_points),
self._repeat_tensor([-3.61205e+06, -8.90607e+06, -3.61205e+06, 0.0, 0.0, 0.0], total_integration_points)
]
mohr_coulomb_file_path = test_helper.get_file_path(os.path.join('test_element_lab', 'test_crs', stage_name, 'mohr_coulomb'))
mohr_coulomb_stresses = test_helper.get_values_from_csv_grouped(
Path(mohr_coulomb_file_path) / "expected_stress.csv",
['phase'], ['element_id', 'ip_index'],
["stress_xx", "stress_yy", "stress_zz", "stress_xy", "stress_yz", "stress_xz"])
mohr_coulomb_strains = test_helper.get_values_from_csv_grouped(
Path(mohr_coulomb_file_path) / "expected_strain.csv",
['phase'], ['element_id', 'ip_index'],
["strain_xx", "strain_yy", "strain_zz", "strain_xy", "strain_yz", "strain_xz"])
expected_stresses = [mohr_coulomb_stresses[phase] for phase in range(1, nr_of_phases + 1)]
expected_strains = [mohr_coulomb_strains[phase] for phase in range(1, nr_of_phases + 1)]

self._run_crs_test_for_model(stage_name, "mohr_coulomb")
self._check_crs_results(stage_name, "mohr_coulomb", nr_of_phases, expected_stresses, expected_strains, expected_water_pressures, delta_stress=1.0)
Expand All @@ -171,13 +167,9 @@ def _check_crs_results(self, stage_name, model_name, nr_of_phases, expected_stre
node_ids = list(range(1, 10))
for i in range(nr_of_phases):
self.assert_integration_point_tensors(
result, "CAUCHY_STRESS_TENSOR",
self._make_integration_point_tensor_entries(expected_stresses[i], num_elements=2, num_integration_points_per_element=3),
times[i], places=places_stress, delta=delta_stress)
result, "CAUCHY_STRESS_TENSOR", expected_stresses[i], times[i], places=places_stress, delta=delta_stress)
self.assert_integration_point_tensors(
result, "ENGINEERING_STRAIN_TENSOR",
self._make_integration_point_tensor_entries(expected_strains[i], num_elements=2, num_integration_points_per_element=3),
times[i], places=places_strain)
result, "ENGINEERING_STRAIN_TENSOR", expected_strains[i], times[i], places=places_strain)
self.assert_nodal_values_at_time(result, "WATER_PRESSURE", node_ids, expected_water_pressures[i], times[i], places=places_water_pressure)

def test_triaxial_comp_6n(self):
Expand Down Expand Up @@ -249,9 +241,6 @@ def _assert_oedometer_effective_stresses(self, effective_stresses, expected_yy,
self.assertAlmostEqual(expected_yy, integration_point[1,1], places_yy)
self.assertAlmostEqual(0.0, integration_point[2,2], 3)

def _repeat_tensor(self, tensor, count):
return [list(tensor) for _ in range(count)]

def _make_integration_point_tensor_entries(self, expected_tensors, num_elements, num_integration_points_per_element):
result = {}
idx = 0
Expand Down
31 changes: 31 additions & 0 deletions ...icsApplication/tests/test_element_lab/test_crs/drained/linear_elastic/expected_strain.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,31 @@
phase,element_id,ip_index,strain_xx,strain_yy,strain_zz,strain_xy,strain_yz,strain_xz
1,1,0,0.0,-0.1,0.0,0.0,0.0,0.0
1,1,1,0.0,-0.1,0.0,0.0,0.0,0.0
1,1,2,0.0,-0.1,0.0,0.0,0.0,0.0
1,2,0,0.0,-0.1,0.0,0.0,0.0,0.0
1,2,1,0.0,-0.1,0.0,0.0,0.0,0.0
1,2,2,0.0,-0.1,0.0,0.0,0.0,0.0
2,1,0,0.0,-0.05,0.0,0.0,0.0,0.0
2,1,1,0.0,-0.05,0.0,0.0,0.0,0.0
2,1,2,0.0,-0.05,0.0,0.0,0.0,0.0
2,2,0,0.0,-0.05,0.0,0.0,0.0,0.0
2,2,1,0.0,-0.05,0.0,0.0,0.0,0.0
2,2,2,0.0,-0.05,0.0,0.0,0.0,0.0
3,1,0,0.0,-0.25,0.0,0.0,0.0,0.0
3,1,1,0.0,-0.25,0.0,0.0,0.0,0.0
3,1,2,0.0,-0.25,0.0,0.0,0.0,0.0
3,2,0,0.0,-0.25,0.0,0.0,0.0,0.0
3,2,1,0.0,-0.25,0.0,0.0,0.0,0.0
3,2,2,0.0,-0.25,0.0,0.0,0.0,0.0
4,1,0,0.0,-0.25,0.0,0.0,0.0,0.0
4,1,1,0.0,-0.25,0.0,0.0,0.0,0.0
4,1,2,0.0,-0.25,0.0,0.0,0.0,0.0
4,2,0,0.0,-0.25,0.0,0.0,0.0,0.0
4,2,1,0.0,-0.25,0.0,0.0,0.0,0.0
4,2,2,0.0,-0.25,0.0,0.0,0.0,0.0
5,1,0,0.0,-0.4,0.0,0.0,0.0,0.0
5,1,1,0.0,-0.4,0.0,0.0,0.0,0.0
5,1,2,0.0,-0.4,0.0,0.0,0.0,0.0
5,2,0,0.0,-0.4,0.0,0.0,0.0,0.0
5,2,1,0.0,-0.4,0.0,0.0,0.0,0.0
5,2,2,0.0,-0.4,0.0,0.0,0.0,0.0
31 changes: 31 additions & 0 deletions ...icsApplication/tests/test_element_lab/test_crs/drained/linear_elastic/expected_stress.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,31 @@
phase,element_id,ip_index,stress_xx,stress_yy,stress_zz,stress_xy,stress_yz,stress_xz
1,1,0,-400000,-1200000,-400000,0.0,0.0,0.0
1,1,1,-400000,-1200000,-400000,0.0,0.0,0.0
1,1,2,-400000,-1200000,-400000,0.0,0.0,0.0
1,2,0,-400000,-1200000,-400000,0.0,0.0,0.0
1,2,1,-400000,-1200000,-400000,0.0,0.0,0.0
1,2,2,-400000,-1200000,-400000,0.0,0.0,0.0
2,1,0,-200000,-600000,-200000,0.0,0.0,0.0
2,1,1,-200000,-600000,-200000,0.0,0.0,0.0
2,1,2,-200000,-600000,-200000,0.0,0.0,0.0
2,2,0,-200000,-600000,-200000,0.0,0.0,0.0
2,2,1,-200000,-600000,-200000,0.0,0.0,0.0
2,2,2,-200000,-600000,-200000,0.0,0.0,0.0
3,1,0,-1000000,-3000000,-1000000,0.0,0.0,0.0
3,1,1,-1000000,-3000000,-1000000,0.0,0.0,0.0
3,1,2,-1000000,-3000000,-1000000,0.0,0.0,0.0
3,2,0,-1000000,-3000000,-1000000,0.0,0.0,0.0
3,2,1,-1000000,-3000000,-1000000,0.0,0.0,0.0
3,2,2,-1000000,-3000000,-1000000,0.0,0.0,0.0
4,1,0,-1000000,-3000000,-1000000,0.0,0.0,0.0
4,1,1,-1000000,-3000000,-1000000,0.0,0.0,0.0
4,1,2,-1000000,-3000000,-1000000,0.0,0.0,0.0
4,2,0,-1000000,-3000000,-1000000,0.0,0.0,0.0
4,2,1,-1000000,-3000000,-1000000,0.0,0.0,0.0
4,2,2,-1000000,-3000000,-1000000,0.0,0.0,0.0
5,1,0,-1600000,-4800000,-1600000,0.0,0.0,0.0
5,1,1,-1600000,-4800000,-1600000,0.0,0.0,0.0
5,1,2,-1600000,-4800000,-1600000,0.0,0.0,0.0
5,2,0,-1600000,-4800000,-1600000,0.0,0.0,0.0
5,2,1,-1600000,-4800000,-1600000,0.0,0.0,0.0
5,2,2,-1600000,-4800000,-1600000,0.0,0.0,0.0
31 changes: 31 additions & 0 deletions ...anicsApplication/tests/test_element_lab/test_crs/drained/mohr_coulomb/expected_strain.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,31 @@
phase,element_id,ip_index,strain_xx,strain_yy,strain_zz,strain_xy,strain_yz,strain_xz
1,1,0,0.0,-0.1,0.0,0.0,0.0,0.0
1,1,1,0.0,-0.1,0.0,0.0,0.0,0.0
1,1,2,0.0,-0.1,0.0,0.0,0.0,0.0
1,2,0,0.0,-0.1,0.0,0.0,0.0,0.0
1,2,1,0.0,-0.1,0.0,0.0,0.0,0.0
1,2,2,0.0,-0.1,0.0,0.0,0.0,0.0
2,1,0,0.0,-0.05,0.0,0.0,0.0,0.0
2,1,1,0.0,-0.05,0.0,0.0,0.0,0.0
2,1,2,0.0,-0.05,0.0,0.0,0.0,0.0
2,2,0,0.0,-0.05,0.0,0.0,0.0,0.0
2,2,1,0.0,-0.05,0.0,0.0,0.0,0.0
2,2,2,0.0,-0.05,0.0,0.0,0.0,0.0
3,1,0,0.0,-0.25,0.0,0.0,0.0,0.0
3,1,1,0.0,-0.25,0.0,0.0,0.0,0.0
3,1,2,0.0,-0.25,0.0,0.0,0.0,0.0
3,2,0,0.0,-0.25,0.0,0.0,0.0,0.0
3,2,1,0.0,-0.25,0.0,0.0,0.0,0.0
3,2,2,0.0,-0.25,0.0,0.0,0.0,0.0
4,1,0,0.0,-0.25,0.0,0.0,0.0,0.0
4,1,1,0.0,-0.25,0.0,0.0,0.0,0.0
4,1,2,0.0,-0.25,0.0,0.0,0.0,0.0
4,2,0,0.0,-0.25,0.0,0.0,0.0,0.0
4,2,1,0.0,-0.25,0.0,0.0,0.0,0.0
4,2,2,0.0,-0.25,0.0,0.0,0.0,0.0
5,1,0,0.0,-0.4,0.0,0.0,0.0,0.0
5,1,1,0.0,-0.4,0.0,0.0,0.0,0.0
5,1,2,0.0,-0.4,0.0,0.0,0.0,0.0
5,2,0,0.0,-0.4,0.0,0.0,0.0,0.0
5,2,1,0.0,-0.4,0.0,0.0,0.0,0.0
5,2,2,0.0,-0.4,0.0,0.0,0.0,0.0
31 changes: 31 additions & 0 deletions ...anicsApplication/tests/test_element_lab/test_crs/drained/mohr_coulomb/expected_stress.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,31 @@
phase,element_id,ip_index,stress_xx,stress_yy,stress_zz,stress_xy,stress_yz,stress_xz
1,1,0,-901878,-2228430,-901878,0.0,0.0,0.0
1,1,1,-901878,-2228430,-901878,0.0,0.0,0.0
1,1,2,-901878,-2228430,-901878,0.0,0.0,0.0
1,2,0,-901878,-2228430,-901878,0.0,0.0,0.0
1,2,1,-901878,-2228430,-901878,0.0,0.0,0.0
1,2,2,-901878,-2228430,-901878,0.0,0.0,0.0
2,1,0,-501878,-1028430,-501878,0.0,0.0,0.0
2,1,1,-501878,-1028430,-501878,0.0,0.0,0.0
2,1,2,-501878,-1028430,-501878,0.0,0.0,0.0
2,2,0,-501878,-1028430,-501878,0.0,0.0,0.0
2,2,1,-501878,-1028430,-501878,0.0,0.0,0.0
2,2,2,-501878,-1028430,-501878,0.0,0.0,0.0
3,1,0,-2256970,-5567250,-2256970,0.0,0.0,0.0
3,1,1,-2256970,-5567250,-2256970,0.0,0.0,0.0
3,1,2,-2256970,-5567250,-2256970,0.0,0.0,0.0
3,2,0,-2256970,-5567250,-2256970,0.0,0.0,0.0
3,2,1,-2256970,-5567250,-2256970,0.0,0.0,0.0
3,2,2,-2256970,-5567250,-2256970,0.0,0.0,0.0
4,1,0,-2256970,-5567250,-2256970,0.0,0.0,0.0
4,1,1,-2256970,-5567250,-2256970,0.0,0.0,0.0
4,1,2,-2256970,-5567250,-2256970,0.0,0.0,0.0
4,2,0,-2256970,-5567250,-2256970,0.0,0.0,0.0
4,2,1,-2256970,-5567250,-2256970,0.0,0.0,0.0
4,2,2,-2256970,-5567250,-2256970,0.0,0.0,0.0
5,1,0,-3612050,-8906070,-3612050,0.0,0.0,0.0
5,1,1,-3612050,-8906070,-3612050,0.0,0.0,0.0
5,1,2,-3612050,-8906070,-3612050,0.0,0.0,0.0
5,2,0,-3612050,-8906070,-3612050,0.0,0.0,0.0
5,2,1,-3612050,-8906070,-3612050,0.0,0.0,0.0
5,2,2,-3612050,-8906070,-3612050,0.0,0.0,0.0
10 changes: 10 additions & 0 deletions ...ns/GeoMechanicsApplication/tests/test_element_lab/test_triaxial/drained/expected_disp.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,10 @@
node_id,disp_x,disp_y,disp_z
1,0.0,-0.2,0.0
2,0.0527776,-0.2,0.0
3,0.0,-0.100033,0.0
4,0.0524025,-0.0996931,0.0
5,0.0,0.0,0.0
6,0.105197,-0.2,0.0
7,0.105114,-0.100049,0.0
8,0.0524406,0.0,0.0
9,0.104632,0.0,0.0
7 changes: 7 additions & 0 deletions .../GeoMechanicsApplication/tests/test_element_lab/test_triaxial/drained/expected_strain.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
element_id,ip_index,strain_xx,strain_yy,strain_zz,strain_xy,strain_yz,strain_xz
1,0,0.104863,-0.19973,0.104946,0.000440186,0.0,0.0
1,1,0.104939,-0.199861,0.105022,0.000100528,0.0,0.0
1,2,0.10444,-0.199182,0.104773,6.24562e-05,0.0,0.0
2,0,0.1055,-0.200303,0.104922,3.84218e-05,0.0,0.0
2,1,0.104915,-0.200106,0.105298,-4.97694e-06,0.0,0.0
2,2,0.105631,-0.200818,0.105335,0.000287088,0.0,0.0
7 changes: 7 additions & 0 deletions .../GeoMechanicsApplication/tests/test_element_lab/test_triaxial/drained/expected_stress.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,7 @@
element_id,ip_index,stress_xx,stress_yy,stress_zz,stress_xy,stress_yz,stress_xz
1,0,-99.9808,-252.622,-99.9806,0.193199,0.0,0.0
1,1,-99.8828,-252.381,-99.8827,0.0448209,0.0,0.0
1,2,-100.045,-252.78,-100.045,0.0314238,0.0,0.0
2,0,-99.9991,-252.668,-99.9991,0.00846584,0.0,0.0
2,1,-99.967,-252.589,-99.967,-0.0103576,0.0,0.0
2,2,-100.064,-252.827,-100.064,0.137162,0.0,0.0
10 changes: 10 additions & 0 deletions .../GeoMechanicsApplication/tests/test_element_lab/test_triaxial/undrained/expected_disp.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,10 @@
node_id,disp_x,disp_y,disp_z
1,0.0,-0.2,0.0
2,0.045,-0.2,0.0
3,0.0,-0.1,0.0
4,0.045,-0.1,0.0
5,0.0,0.0,0.0
6,0.09,-0.2,0.0
7,0.09,-0.1,0.0
8,0.045,0.0,0.0
9,0.09,0.0,0.0
3 changes: 3 additions & 0 deletions ...eoMechanicsApplication/tests/test_element_lab/test_triaxial/undrained/expected_strain.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
element_id,ip_index,strain_xx,strain_yy,strain_zz,strain_xy,strain_yz,strain_xz
1,0,0.09,-0.2,0.09,0.0,0.0,0.0
2,0,0.09,-0.2,0.09,0.0,0.0,0.0
3 changes: 3 additions & 0 deletions ...eoMechanicsApplication/tests/test_element_lab/test_triaxial/undrained/expected_stress.csv
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,3 @@
element_id,ip_index,stress_xx,stress_yy,stress_zz,stress_xy,stress_yz,stress_xz
1,0,-100.0,-4740.0,-100.0,0.0,0.0,0.0
2,0,-100.0,-4740.0,-100.0,0.0,0.0,0.0
Loading
Loading