[GeoMechanicsApplication] add a new oedometer test to the test_element_lab by Sabine-van-Dijk · Pull Request #14403 · KratosMultiphysics/Kratos

Sabine-van-Dijk · 2026-04-24T14:54:52Z

📝 Description

Set up a test for the oedometer experiment. + small changes to the triaxial test.

🆕 Changelog

Made a mesh and set the material and project parameters.
Implemented an assertion function by Wijtze Pieter.
Assert the displacements and normal stresses for the oedometer test.
Apply assert function also for the triaxial test.
Add a readme file & schematic.

…est (unit convention = kN)

WPK4FEM

Hi Sabine, Thank you for this nice first setup. I have a lot of small nagging comments, but am happy with the completeness of test, documentation and picture. It would be nice to also have pictures of stress development in time in the documentation, as this is the way these kind of experiments usually are assessed.

Apparently you have not experienced non-uniform deformations in the test. That is something I fear from past experience and that ( if it starts happening ) should be curable with master-slave constraints that force a uniform displacement field.

Lets await Marjan's comments before starting to change things.
Regards, Wijtze Pieter

WPK4FEM · 2026-04-29T09:35:56Z

+                    "PERMEABILITY_YY": 4.5e-30,
+                    "PERMEABILITY_XY": 0.0,
+                    "DYNAMIC_VISCOSITY": 8.9e-07,
+                    "THICKNESS": 1.0,


For plane strain and axi-symmetric models, THICKNESS should not be provided.

Removed it.
For which should it be provided?

WPK4FEM · 2026-04-29T09:38:22Z

+        "move_mesh_flag": false,
+        "reform_dofs_at_each_step": false,
+        "block_builder": true,
+        "reset_displacements": true,


In a 1 stage computation, false would be a more natural value. reset_displacements is used for setting total displacement fields back to 0.0 at the start of a later stage.

indigocoral · 2026-04-29T21:49:30Z

+            "properties_id": 1,
+            "Material": {
+                "constitutive_law": {
+                    "name": "GeoLinearElasticPlaneStrain2DLaw"


@WPK4FEM Does it make sense to have the same material properties for all these lab tests? Sabine has already kindly made the Triaxial test with the Mohr-Coulomb material, and I'm wondering if we can reuse the same material properties for this oedometer test as well? Would Mohr-Coulomb make sense, or would you prefer linear elastic? (Then we can basically have one material file in a common folder and reuse that for all types of lab tests.)

I had a brief discussion with @WPK4FEM about this and he would like to have 2 sets of tests per each lab test: one for linear elastic material and another one for Mohr-Coulomb material. For Triaxial test you already have the Mohr-Coulomb test and for this oedometer test, you have set up a linear elastic one. Perhaps after merging this, you can also add the Mohr-Coulomb model (most likely the same material file you used for the Triaxial test) to this oedometer test and also (the same) linear elastic material for the Triaxial test.

Having shared material parameters files abolutely makes sense. It would enable running all ( triaxial, oedometer, CRS etc ) tests for a material.

indigocoral · 2026-04-29T22:03:04Z

+                    "model_part_name": "PorousDomain.Top_load",
+                    "variable_name": "LINE_LOAD",
+                    "active": [false, true, false],
+                    "value": [0.0, 0.0, 0.0],


To make it stepwise load function:

Suggested change

"value": [0.0, 0.0, 0.0],

"value": [0.0, -1000.0, 0.0],

indigocoral · 2026-04-29T22:07:06Z

Very nice detailed picture! 👍🏼

indigocoral · 2026-04-29T22:16:16Z

+  End SubModelPartNodes
+End SubModelPart
+
+Begin SubModelPart Sides_L_R


Maybe a more clear name for this would be:

Suggested change

Begin SubModelPart Sides_L_R

Begin SubModelPart Sides_rollers

indigocoral

Thanks a lot, Sabine, for making this test and also producing the nice graph and nice results checking. Regarding non-uniform deformations that @WPK4FEM mentioned, yes we have not experienced non-uniform deformations in this test. Perhaps this is something to discuss later? :)
Thanks a lot Sabine!

indigocoral · 2026-04-29T22:31:14Z

+            "properties_id": 1,
+            "Material": {
+                "constitutive_law": {
+                    "name": "GeoLinearElasticPlaneStrain2DLaw"


I had a brief discussion with @WPK4FEM about this and he would like to have 2 sets of tests per each lab test: one for linear elastic material and another one for Mohr-Coulomb material. For Triaxial test you already have the Mohr-Coulomb test and for this oedometer test, you have set up a linear elastic one. Perhaps after merging this, you can also add the Mohr-Coulomb model (most likely the same material file you used for the Triaxial test) to this oedometer test and also (the same) linear elastic material for the Triaxial test.

… changes in undrained branch

Co-authored-by: Copilot <copilot@github.com>

Sabine-van-Dijk · 2026-05-01T13:22:52Z

To do:

Discuss value "BULK_MODULUS_FLUID"
Value was not used in the tests, so I removed it from the material parameter files.
For the drained triaxial test it is used so we need to make a new issue to look at that.
Discuss axisymmetric element type
For the axisymmetric element type the displacement in "the middle" (near the symmetry axis) is way more than on the outside. This needs to be constrained. --> axisymmetric loading
Apply stepwise loading
Discuss value of the permeability.
Value was not used in the tests, so I removed it from the material parameter files.

… si units (N instead of kN)

Co-authored-by: Copilot <copilot@github.com>

WPK4FEM

Hi Sabine,
Some more small remarks,
Wijtze Pieter

Sabine-van-Dijk added 8 commits April 24, 2026 14:15

make a mesh, material and project parameters file for the oedometer t…

9e54d07

…est (unit convention = kN)

insert assertion function and apply it in the triaxial test

d5e715d

add oedometer test with assertion function

db00ae7

apply assert function on displacement assertion triaxial test

a2db314

formatting of json files for the oedometer test

94b6584

remove triaxial material parameters from the json format exception list

ab63228

minor changes to triaxial test readme

b9e2305

add readme for the oedometer test incl schematic

5e2f231

Sabine-van-Dijk requested a review from a team as a code owner April 24, 2026 14:54

Sabine-van-Dijk linked an issue Apr 24, 2026 that may be closed by this pull request

[GeoMechanicsApplication] add oedometer test #14393

Closed

Sabine-van-Dijk requested review from WPK4FEM and indigocoral April 24, 2026 14:55

WPK4FEM reviewed Apr 29, 2026

View reviewed changes

indigocoral previously approved these changes Apr 29, 2026

View reviewed changes

indigocoral self-requested a review April 29, 2026 22:22

indigocoral reviewed Apr 29, 2026

View reviewed changes

indigocoral self-requested a review April 29, 2026 22:39

Sabine-van-Dijk and others added 9 commits April 30, 2026 11:41

changes in mesh, project and material files according to comments and…

772c9a4

… changes in undrained branch

add stage name drained and add functions in triaxial test

178f5a0

make equal with master

ffc54e6

Co-authored-by: Copilot <copilot@github.com>

remove old code

d501966

change to geo drainage type

a7dc6b8

change nr of checked ip's in triaxial test

dba14b7

Co-authored-by: Copilot <copilot@github.com>

make code use general helper functions

5c69788

Co-authored-by: Copilot <copilot@github.com>

change output file type

83fb795

change the poisson ratio

9cca881

Sabine-van-Dijk dismissed indigocoral’s stale review via 9cca881 May 1, 2026 13:14

Sabine-van-Dijk added 2 commits May 4, 2026 11:49

change output file name

b4baf8b

remove value for the bulk modulus fluid

2cb4988

Sabine-van-Dijk and others added 5 commits May 5, 2026 10:34

apply stepwise loading + remove unused material parameters and set to…

e4aa62d

… si units (N instead of kN)

si units in assertion and stepwise loading

c537f25

fix fixed base constraint

0ebbb37

use axisymmetric element type and loading

f34cc3b

Update readme

fa12b46

Co-authored-by: Copilot <copilot@github.com>

WPK4FEM reviewed May 5, 2026

View reviewed changes

apply review comments

4dcf3c0

WPK4FEM approved these changes May 5, 2026

View reviewed changes

Sabine-van-Dijk merged commit 7585c09 into master May 5, 2026
10 checks passed

Sabine-van-Dijk deleted the geo/add-oedometer-test branch May 5, 2026 15:10

	Begin SubModelPart Sides_L_R
	Begin SubModelPart Sides_rollers

Conversation

Sabine-van-Dijk commented Apr 24, 2026

📝 Description

🆕 Changelog

Uh oh!

WPK4FEM left a comment

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

indigocoral left a comment

Choose a reason for hiding this comment

Uh oh!

Choose a reason for hiding this comment

Uh oh!

Sabine-van-Dijk commented May 1, 2026 • edited Loading Uh oh! There was an error while loading. Please reload this page.

Uh oh!

Uh oh!

WPK4FEM left a comment

Choose a reason for hiding this comment

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

3 participants

Sabine-van-Dijk commented May 1, 2026 •

edited

Loading