Skip to content
Snippets Groups Projects
Commit caec7a3d authored by Thomas Nagel's avatar Thomas Nagel Committed by Dmitri Naumov
Browse files

[T/SD] Guenther Salzer model triaxial test.

parent fb00536d
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 465 additions and 0 deletions
MaterialLib/SolidModels/MFront/CMakeLists.txt
+ 1
− 0
View file @ caec7a3d
......@@ -9,6 +9,7 @@ mfront_behaviours_check_library(
BDT
DruckerPrager
Elasticity
GuentherSalzer
Lubby2
Lubby2mod
ModCamClay_semiExplParaInit
......
ProcessLib/SmallDeformation/Tests.cmake
+ 1
− 0
View file @ caec7a3d
......@@ -76,6 +76,7 @@ if (OGS_USE_MFRONT)
OgsTest(PROJECTFILE Mechanics/ModifiedCamClay/square_1e0_shear.prj)
OgsTest(PROJECTFILE Mechanics/ModifiedCamClay/square_1e0_biax.prj)
OgsTest(PROJECTFILE Mechanics/ModifiedCamClay/model_triaxtest.prj)
OgsTest(PROJECTFILE Mechanics/GuentherSalzer/model_triaxtest.prj)
# Linear elastic, no internal state variables, no external state variables.
AddTest(
......
Tests/Data/Mechanics/GuentherSalzer/geometry_bottom.vtu 0 → 100644
+ 34
− 0
View file @ caec7a3d
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
<!--This file was created by meshio v4.3.5-->
<UnstructuredGrid>
<Piece NumberOfPoints="5" NumberOfCells="4">
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary">
AQAAAACAAAB4AAAAJgAAAA==eJxjYMAHLB3QRf61/AcCCQxxnvkgcQ0M8YxbQOH9RhjiABHwEM8=
</DataArray>
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="binary">
AQAAAACAAABAAAAAFgAAAA==eJxjYIAAJjSaGY1mQaMZoTQAApAAFA==
</DataArray>
<DataArray type="Int64" Name="offsets" format="binary">
AQAAAACAAAAgAAAAEwAAAA==eJxjYoAAFijNBqU5oDQAAWAAFQ==
</DataArray>
<DataArray type="Int64" Name="types" format="binary">
AQAAAACAAAAgAAAADgAAAA==eJxjZoAAZhw0AAEQAA0=
</DataArray>
</Cells>
<PointData>
<DataArray type="UInt64" Name="bulk_node_ids" format="binary">
AQAAAACAAAAoAAAAFgAAAA==eJxjYIAARijNAqVZoTQblAYAASgAEQ==
</DataArray>
</PointData>
<CellData>
<DataArray type="UInt64" Name="bulk_elem_ids" format="binary">
AQAAAACAAAAgAAAAEwAAAA==eJxjYIAAYSitBqUtoTQABhAAcw==
</DataArray>
</CellData>
</Piece>
</UnstructuredGrid>
</VTKFile>
Tests/Data/Mechanics/GuentherSalzer/geometry_left.vtu 0 → 100644
+ 34
− 0
View file @ caec7a3d
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
<!--This file was created by meshio v4.3.5-->
<UnstructuredGrid>
<Piece NumberOfPoints="20" NumberOfCells="19">
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary">
AQAAAACAAADgAQAAoQAAAA==eJxjYCAGRDqgi1ipnsrWWB+OIV7iNeV6QFwYhvi6/DjnSt5QDPEXUzTWL9wTjCGutOuT1KnsIAzxmPt72j5JBWKIa7PzPTWd4o8hLmi7CKjBB0P8S5EZ0KFeGOK3V4D85Y4hvv8eyP2uGOJLRT4D3emEIf47H+QeewzxJZtA9lphiPt/BZlviiH+AWyOPoa4O1i9Koa4A1hcFC4OACQZPiE=
</DataArray>
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="binary">
AQAAAACAAAAwAQAANQAAAA==eJxdxTcSgCAAADDEgiBY8P9/ZfBckiVT+ESeeeGVN068c+bCB1dufPLFNz/c+eXfAKwoAXw=
</DataArray>
<DataArray type="Int64" Name="offsets" format="binary">
AQAAAACAAACYAAAANQAAAA==eJwtxdEGgDAAAMDJTDLJZCbJTPr/T+xhdy+3hCk6efXm7N2Hi09XN1++/bh7+PXnH1O4AX0=
</DataArray>
<DataArray type="Int64" Name="types" format="binary">
AQAAAACAAACYAAAADwAAAA==eJxjZoAA5kFGAwASaAA6
</DataArray>
</Cells>
<PointData>
<DataArray type="UInt64" Name="bulk_node_ids" format="binary">
AQAAAACAAACgAAAAMwAAAA==eJwtxckNgCAAADAS10AFQcFr//V40H4awrQ4evXm3cnZh4urT19u7r79+PXn3wO2UAKV
</DataArray>
</PointData>
<CellData>
<DataArray type="UInt64" Name="bulk_elem_ids" format="binary">
AQAAAACAAACYAAAAMQAAAA==eJxVxUkCgBAAAECpiKLl/4/tYC7mMl8YXj++3d18+XR18eHs5N2bV0cvYfYDR9gArA==
</DataArray>
</CellData>
</Piece>
</UnstructuredGrid>
</VTKFile>
Tests/Data/Mechanics/GuentherSalzer/geometry_right.vtu 0 → 100644
+ 34
− 0
View file @ caec7a3d
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
<!--This file was created by meshio v4.3.5-->
<UnstructuredGrid>
<Piece NumberOfPoints="20" NumberOfCells="19">
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary">
AQAAAACAAADgAQAAqwAAAA==eJxjYAABSwcGrAAmHumALn7iU6xzJa8ohnjcD5C4KoZ4+huep6ZT9DHE49+A1JtiiHclmU65HmCFqf4myBx7DPF5hp+kTmU7YYj/+QYy3xVDPHbPyWyN9e4Y4nsaQfZ6YYhLuy/c0/bJB0O8kocX6B5/DHGTIqByqUAMce14kDuDMMSVfDTWL9wTjCEuZRkHdH8ohriQGtCZcWEY4hzCp4D+CoeLAwCUPkQ7
</DataArray>
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="binary">
AQAAAACAAAAwAQAANAAAAA==eJxdxTECQCAAAMAQRSL5/18Nme6WC2GYeeHIK2+cOPPOBxc+ufLFNzd+uPPL0/8HqwABfA==
</DataArray>
<DataArray type="Int64" Name="offsets" format="binary">
AQAAAACAAACYAAAANQAAAA==eJwtxdEGgDAAAMDJTDLJZCbJTPr/T+xhdy+3hCk6efXm7N2Hi09XN1++/bh7+PXnH1O4AX0=
</DataArray>
<DataArray type="Int64" Name="types" format="binary">
AQAAAACAAACYAAAADwAAAA==eJxjZoAA5kFGAwASaAA6
</DataArray>
</Cells>
<PointData>
<DataArray type="UInt64" Name="bulk_node_ids" format="binary">
AQAAAACAAACgAAAAMgAAAA==eJwtxbcBgCAAADBksSAIFuT/Sx1MlkzhFz178erNybuziw9XN5++fPtx9+vhD0YgARs=
</DataArray>
</PointData>
<CellData>
<DataArray type="UInt64" Name="bulk_elem_ids" format="binary">
AQAAAACAAACYAAAAMAAAAA==eJwtxUECQCAAALC+pYgUiv//x6HtsiVM0cmrN2fvPlx8urr58u3H3cOvP/92twTn
</DataArray>
</CellData>
</Piece>
</UnstructuredGrid>
</VTKFile>
Tests/Data/Mechanics/GuentherSalzer/geometry_top.vtu 0 → 100644
+ 34
− 0
View file @ caec7a3d
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
<!--This file was created by meshio v4.3.5-->
<UnstructuredGrid>
<Piece NumberOfPoints="5" NumberOfCells="4">
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary">
AQAAAACAAAB4AAAAKwAAAA==eJxjYAABSwcwxRAJpdEBQvyKAZA4YISh/nkviNTEEPc5DCIlMcQB72MI8g==
</DataArray>
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="binary">
AQAAAACAAABAAAAAFgAAAA==eJxjYIAAJjSaGY1mQaMZoTQAApAAFA==
</DataArray>
<DataArray type="Int64" Name="offsets" format="binary">
AQAAAACAAAAgAAAAEwAAAA==eJxjYoAAFijNBqU5oDQAAWAAFQ==
</DataArray>
<DataArray type="Int64" Name="types" format="binary">
AQAAAACAAAAgAAAADgAAAA==eJxjZoAAZhw0AAEQAA0=
</DataArray>
</Cells>
<PointData>
<DataArray type="UInt64" Name="bulk_node_ids" format="binary">
AQAAAACAAAAoAAAAFgAAAA==eJxjYoAAZigtCaWloLQ0lAYABagAVA==
</DataArray>
</PointData>
<CellData>
<DataArray type="UInt64" Name="bulk_elem_ids" format="binary">
AQAAAACAAAAgAAAAEwAAAA==eJzzZoAACyitCqWFoDQAEaAAuw==
</DataArray>
</CellData>
</Piece>
</UnstructuredGrid>
</VTKFile>
Tests/Data/Mechanics/GuentherSalzer/model.vtu 0 → 100644
+ 29
− 0
View file @ caec7a3d
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
<!--This file was created by meshio v4.3.5-->
<UnstructuredGrid>
<Piece NumberOfPoints="100" NumberOfCells="76">
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary">
AQAAAACAAABgCQAAEQMAAA==eJx91W1IU1EYB/CHMpS5LdPMnBJiaDPFqTVLQbwtcuBLlpoa5h1IL6KEIhbYh6VhVpJBtJBYb1MTIVLQsLChoRZrX6QIrCQrwjKRqJUaVtq9Ppg8O0fPh/Pht8O5557/8+wCrDQShJX94DK/L/lc7bw0NjLrlLdk1zJe/Ebivrhln+twFe6uUvkL7i7+kj2M8WOTyjG9Rce4aVJer2e8vkhvGd6XyK5/Le+TzPjNWJfGWWpg/M+0vP8exgvtz0q1HUbG7TXyc1MZDzLa7HWudMarlCrpPJmMb6+Qlmv2Mx5pks+ZxXhourbDZs9mXJMgSuc/wLhvuHRMMZdxLz+n9F55//1ljDQ9XsxxqR4+N8hzBOPpA/IcyPjiSAyj+y+OylTeeQDay3jnBxi38N4XILSHdz8Ah97x7hMg0lPNuX+AdUlNnLwAflbEc/IFGGlzcuoBoG9U5NQPwJ31Pzj1BvC7TM2pT4CWznhOPQNkTomc+gf4trCPe78AGBfWu/cXgLDgS/3YexX7vek+9epWdHU33Ud4gK5cS5+7yoHu10XP2f8K/YZI36v6C3qOgt6DYRb9wwjtUw9veQ4UhlroPT8NQn90nOZyNgq9LZ7mmJKEfmWe9qnXXnSzg/apQ0T3HaR1db4M/UgzrcO0avTuM7RuFZfRFUW0zp029IJdtC/qO9HvhdA+un0X/4c9Kmle0Q/Rh07SvPoH0bXtKpJX7nP04XKa19hb9EYFzcs8gd7eoiJ5qWbQQ5tpXtbV8hwhvC+heel80K1xNK+BYPT8Wfq/mheBvqGf5vVJj/7iAs3rlAHd0UbzUmeit16keV0vQK8rp3nFFqMfzqF5PalEN+ykeWXXoIcE07ymR/H7eHSbW399RG/cQfPyGUdPzaB5WSfRrdE0r7Dv6Je66Xewawp9TTLNS5hFN5XQvJxz6F9jaF75co7Sd+H0DP0Ojnmie/fSvCqU6NdqaV5/fdDD02he5/zRzQ00rwAN+sQJmlfTJvQskeal24zem0Lz6tmCvlVH8zJGoVsCaF7/AIA38IE=
</DataArray>
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="binary">
AQAAAACAAACACQAAcgEAAA==eJxd0sVWFQAARVEEREJQFJRGygCku+s9urvr//+BAedM7p3stc74lpS8rxyH8H9oH8bB0D6CA6F9FPtD+xj+C+3j+De0T+Cf0D6Jv0P7FPaF9mnsDe0z2BPaZ7E7tM9hV2ifx1+hfQE7Q/sidoT2JWwP7cvYFtpbsQz91UdcwaHQvorDoX0NR0L7Oo6G9g0cC+0FHA/tRZwI7Zs4Gdq3cCq0b+N0aN/BmdC+i7OhfQ/nQvs+zof2A1wI7Ye4GNqPcCm0H+NyaG9Bf+avKvAEV0L7Ka6G9jNcC+3nuB7aL3AjtF9iIbRfYTG0X+NmaL/BrdB+i9uh/Q53Qvs97ob2B9wL7Y+4H9qf8CC0P+NhaH/Bo9D+isehvRn9mb/6gJ/wJLRX4mlor8Kz0F6N56G9Bi9C+2e8DO21eBXa6/A6tH/Bm9D+FW9Dez3ehfZveB/av+NDaG/Ax9DeiE+h/Qc+h/af+BLam/A1tJeiP3sDGe9FPQ==
</DataArray>
<DataArray type="Int64" Name="offsets" format="binary">
AQAAAACAAABgAgAAqAAAAA==eJwtxRF0QgEAAMB6LwiCIPgQBEEQBEEQBEEQBMEgCAaDIAiCIAiCIAiCIAiCYDAIgkEQBMEgCAZBEARBEATdyUVCb1HHHHfCgZNOOe2Ms84574KLLrnsiquuue4PN9z0p7/cctsdd91z3wMPPfLYE08989wLL73yt3+89sa/3nrnvQ/+89En//vsi6+++e6Hnw6F30UcdcxxJxw46ZTTzjjrnF9UqSDG
</DataArray>
<DataArray type="Int64" Name="types" format="binary">
AQAAAACAAABgAgAAEgAAAA==eJzjZIAAzlF6lKYBDQA5fQKt
</DataArray>
</Cells>
<CellData>
<DataArray type="Int32" Name="MaterialIDs" format="binary">
AQAAAACAAAAwAQAADQAAAA==eJxjYBgFpAAAATAAAQ==
</DataArray>
</CellData>
</Piece>
</UnstructuredGrid>
</VTKFile>
Tests/Data/Mechanics/GuentherSalzer/model_triaxtest.prj 0 → 100644
+ 259
− 0
View file @ caec7a3d
<?xml version="1.0" encoding="ISO-8859-1"?>
<!-- units: kg, m, s, N -->
<OpenGeoSysProject>
<meshes>
<mesh axially_symmetric="true">model.vtu</mesh>
<mesh axially_symmetric="true">geometry_left.vtu</mesh>
<mesh axially_symmetric="true">geometry_right.vtu</mesh>
<mesh axially_symmetric="true">geometry_top.vtu</mesh>
<mesh axially_symmetric="true">geometry_bottom.vtu</mesh>
</meshes>
<processes>
<process>
<name>SD</name>
<type>SMALL_DEFORMATION</type>
<integration_order>2</integration_order>
<constitutive_relation>
<type>MFront</type>
<behaviour>GuentherSalzer</behaviour>
<material_properties>
<material_property name="YoungModulus" parameter="YoungModulus"/>
<material_property name="PoissonRatio" parameter="PoissonRatio"/>
<material_property name="PrimaryPowerLawFactor" parameter="Ap"/>
<material_property name="PrimaryPowerLawExponent" parameter="np"/>
<material_property name="SecondaryPowerLawFactor1" parameter="As1"/>
<material_property name="SecondaryPowerLawExponent1" parameter="ns1"/>
<material_property name="SecondaryActivationEnergy1" parameter="Q1"/>
<material_property name="SecondaryPowerLawFactor2" parameter="As2"/>
<material_property name="SecondaryPowerLawExponent2" parameter="ns2"/>
<material_property name="SecondaryActivationEnergy2" parameter="Q2"/>
<material_property name="HardeningExponent" parameter="mup"/>
<material_property name="InitialHardening" parameter="epsV0"/>
<material_property name="ReferenceStress" parameter="sref"/>
</material_properties>
</constitutive_relation>
<solid_density>rho_sr</solid_density>
<specific_body_force>0 0</specific_body_force>
<process_variables>
<process_variable>displacement</process_variable>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="sigma" output_name="sigma"/>
<secondary_variable internal_name="epsilon" output_name="epsilon"/>
<secondary_variable internal_name="ElasticStrain" output_name="ElasticStrain"/>
<secondary_variable internal_name="HardeningStrain" output_name="HardeningStrain"/>
</secondary_variables>
<reference_temperature>T_ref</reference_temperature>
</process>
</processes>
<time_loop>
<processes>
<process ref="SD">
<nonlinear_solver>basic_newton</nonlinear_solver>
<convergence_criterion>
<type>DeltaX</type>
<norm_type>NORM2</norm_type>
<abstol>1e-12</abstol>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>IterationNumberBasedTimeStepping</type>
<t_initial>0</t_initial>
<t_end>200</t_end>
<initial_dt>1</initial_dt>
<minimum_dt>1.e-3</minimum_dt>
<maximum_dt>10.</maximum_dt>
<number_iterations>1 5 10 15</number_iterations>
<multiplier>1.2 1.0 0.8 0.6</multiplier>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>triaxtest_output</prefix>
<timesteps>
<pair>
<repeat>10000</repeat>
<each_steps>1</each_steps>
</pair>
</timesteps>
<variables>
<variable>displacement</variable>
<variable>sigma</variable>
<variable>epsilon</variable>
</variables>
</output>
</time_loop>
<parameters>
<parameter>
<name>YoungModulus</name>
<type>Constant</type>
<value>1e5</value>
</parameter>
<parameter>
<name>PoissonRatio</name>
<type>Constant</type>
<value>0.27</value>
</parameter>
<parameter>
<name>Ap</name>
<type>Constant</type>
<value>6.5e-24</value>
</parameter>
<parameter>
<name>np</name>
<type>Constant</type>
<value>11.88</value>
</parameter>
<parameter>
<name>As1</name>
<type>Constant</type>
<value>0.001</value>
</parameter>
<parameter>
<name>ns1</name>
<type>Constant</type>
<value>2.7</value>
</parameter>
<parameter>
<name>Q1</name>
<type>Constant</type>
<value>40e3</value>
</parameter>
<parameter>
<name>As2</name>
<type>Constant</type>
<value>0.001</value>
</parameter>
<parameter>
<name>ns2</name>
<type>Constant</type>
<value>6.8</value>
</parameter>
<parameter>
<name>Q2</name>
<type>Constant</type>
<value>65e3</value>
</parameter>
<parameter>
<name>mup</name>
<type>Constant</type>
<value>2.95</value>
</parameter>
<parameter>
<name>sref</name>
<type>Constant</type>
<value>1.0</value>
</parameter>
<parameter>
<name>epsV0</name>
<type>Constant</type>
<value>1.e-4</value>
</parameter>
<parameter>
<name>T_ref</name>
<type>Constant</type>
<value>293.15</value>
</parameter>
<!-- Initial and boundary values -->
<parameter>
<name>rho_sr</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>displacement0</name>
<type>Constant</type>
<values>0 0</values>
</parameter>
<parameter>
<name>zero</name>
<type>Constant</type>
<value>0.0</value>
</parameter>
<parameter>
<name>axial_pressure</name>
<type>CurveScaled</type>
<curve>ax_loading_curve</curve>
<parameter>loading_value_top</parameter>
</parameter>
<parameter>
<name>loading_value_top</name>
<type>Constant</type>
<value>-1</value> <!--Pa-->
</parameter>
</parameters>
<curves>
<curve>
<name>ax_loading_curve</name>
<coords>0.0 0.1 60.0 60.1 200</coords>
<values>0.0 16.0 16.0 25.0 25.0</values>
</curve>
</curves>
<process_variables>
<process_variable>
<name>displacement</name>
<components>2</components>
<order>1</order>
<initial_condition>displacement0</initial_condition>
<boundary_conditions>
<!--fix left in radial direction-->
<boundary_condition>
<mesh>geometry_left</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>zero</parameter>
</boundary_condition>
<!--fix bottom in axial direction-->
<boundary_condition>
<mesh>geometry_bottom</mesh>
<type>Dirichlet</type>
<component>1</component>
<parameter>zero</parameter>
</boundary_condition>
<!--compression in axial direction -->
<boundary_condition>
<mesh>geometry_top</mesh>
<type>Neumann</type>
<component>1</component>
<parameter>axial_pressure</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_newton</name>
<type>Newton</type>
<max_iter>60</max_iter>
<linear_solver>general_linear_solver</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>general_linear_solver</name>
<eigen>
<solver_type>SparseLU</solver_type>
<scaling>true</scaling>
</eigen>
</linear_solver>
</linear_solvers>
<test_definition>
<!--primary field-->
<vtkdiff>
<file>triaxtest_output_ts_199_t_199.000000.vtu</file>
<field>displacement</field>
<absolute_tolerance>2e-14</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
<!--secondary field-->
<vtkdiff>
<file>triaxtest_output_ts_199_t_199.000000.vtu</file>
<field>sigma</field>
<absolute_tolerance>9e-8</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
</test_definition>
</OpenGeoSysProject>
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment