[Bbar#M] A 3D test

ufz · Oct 2, 2024 · de3641c · de3641c
1 parent 7483537
commit de3641c
Show file tree

Hide file tree

Showing 10 changed files with 564 additions and 0 deletions.
diff --git a/ProcessLib/SmallDeformation/Tests.cmake b/ProcessLib/SmallDeformation/Tests.cmake
@@ -303,6 +303,21 @@ AddTest(
  axisymmetric_bbar_ts_1_t_1.000000.vtu axisymmetric_bbar_ts_1_t_1.000000.vtu analytic_eps epsilon 1.e-9 5.e-7
 )
 
+AddTest(
+ NAME Mechanics_Simple3DBbar
+ PATH Mechanics/Simple3DBbar
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS simple_3d_bbar.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ RUNTIME 1
+ DIFF_DATA
+ simple_3d_bbar_ts_1_t_1.000000.vtu simple_3d_bbar_ts_1_t_1.000000.vtu analytic_displacement displacement 1.e-10 1.e-10
+ simple_3d_bbar_ts_1_t_1.000000.vtu simple_3d_bbar_ts_1_t_1.000000.vtu analytic_sigma sigma 6.e-5 5.e-7
+ simple_3d_bbar_ts_1_t_1.000000.vtu simple_3d_bbar_ts_1_t_1.000000.vtu analytic_eps epsilon 1.e-9 5.e-7
+)
+
 if(NOT OGS_USE_PETSC)
  NotebookTest(NOTEBOOKFILE Mechanics/CooksMembrane/CooksMembraneBbar.ipynb RUNTIME 1)
  NotebookTest(NOTEBOOKFILE Mechanics/Linear/SimpleMechanics.ipynb RUNTIME 5)

diff --git a/Tests/Data/Mechanics/Simple3DBbar/analytical_solution.py b/Tests/Data/Mechanics/Simple3DBbar/analytical_solution.py
@@ -0,0 +1,36 @@
+import pyvista as pv
+
+mesh = pv.read("simple_3d_bbar_ts_1_t_1.000000.vtu")
+E = 1.0e10
+nv = 0.499
+p = 1e7
+
+sigma_analytical = mesh.point_data["sigma"]
+sigma_analytical[:, 0] = 0.0
+sigma_analytical[:, 1] = 0.0
+sigma_analytical[:, 2] = -p
+sigma_analytical[:, 3] = 0.0
+sigma_analytical[:, 4] = 0.0
+sigma_analytical[:, 5] = 0.0
+mesh.point_data["analytic_sigma"] = sigma_analytical
+
+eps_analytical = mesh.point_data["epsilon"]
+eps_1 = -p / E
+eps_0 = -eps_1 * nv
+eps_analytical[:, 0] = eps_0
+eps_analytical[:, 1] = eps_0
+eps_analytical[:, 2] = eps_1
+eps_analytical[:, 3] = 0.0
+eps_analytical[:, 4] = 0.0
+eps_analytical[:, 5] = 0.0
+mesh.point_data["analytic_eps"] = eps_analytical
+
+
+u_analytical = mesh.point_data["displacement"]
+u_analytical[:, 0] = eps_0 * mesh.points[:, 0]
+u_analytical[:, 1] = eps_0 * mesh.points[:, 1]
+u_analytical[:, 2] = eps_1 * mesh.points[:, 2]
+mesh.point_data["analytic_displacement"] = u_analytical
+
+mesh.save("simple_3d_bbar_ts_1_t_1.000000.vtu")
+# mesh.plot(scalars="analytic_displacement")
diff --git a/Tests/Data/Mechanics/Simple3DBbar/bottom.vtu b/Tests/Data/Mechanics/Simple3DBbar/bottom.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="63" NumberOfCells="40">
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="62">
+ +AEAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA=
+ </DataArray>
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="binary" RangeMin="0" RangeMax="0">
+ oAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
+ </DataArray>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="binary" RangeMin="0" RangeMax="39">
+ QAEAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAAA==
+ </DataArray>
+ <DataArray type="UInt64" Name="bulk_face_ids" format="binary" RangeMin="0" RangeMax="0">
+ QAEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA==
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="1.004987562112089">
+ 6AUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACamZmZmZmpPwAAAAAAAAAAAAAAAAAAAACamZmZmZm5PwAAAAAAAAAAAAAAAAAAAAAzMzMzMzPDPwAAAAAAAAAAAAAAAAAAAACamZmZmZnJPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADQPwAAAAAAAAAAAAAAAAAAAAAzMzMzMzPTPwAAAAAAAAAAAAAAAAAAAABmZmZmZmbWPwAAAAAAAAAAAAAAAAAAAACamZmZmZnZPwAAAAAAAAAAAAAAAAAAAADNzMzMzMzcPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADgPwAAAAAAAAAAAAAAAAAAAACamZmZmZnhPwAAAAAAAAAAAAAAAAAAAAAzMzMzMzPjPwAAAAAAAAAAAAAAAAAAAADNzMzMzMzkPwAAAAAAAAAAAAAAAAAAAABmZmZmZmbmPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADoPwAAAAAAAAAAAAAAAAAAAACamZmZmZnpPwAAAAAAAAAAAAAAAAAAAAAzMzMzMzPrPwAAAAAAAAAAAAAAAAAAAADNzMzMzMzsPwAAAAAAAAAAAAAAAAAAAABmZmZmZmbuPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJqZmZmZmak/AAAAAAAAAACamZmZmZmpP5qZmZmZmak/AAAAAAAAAACamZmZmZm5P5qZmZmZmak/AAAAAAAAAAAzMzMzMzPDP5qZmZmZmak/AAAAAAAAAACamZmZmZnJP5qZmZmZmak/AAAAAAAAAAAAAAAAAADQP5qZmZmZmak/AAAAAAAAAAAzMzMzMzPTP5qZmZmZmak/AAAAAAAAAABmZmZmZmbWP5qZmZmZmak/AAAAAAAAAACamZmZmZnZP5qZmZmZmak/AAAAAAAAAADNzMzMzMzcP5qZmZmZmak/AAAAAAAAAAAAAAAAAADgP5qZmZmZmak/AAAAAAAAAACamZmZmZnhP5qZmZmZmak/AAAAAAAAAAAzMzMzMzPjP5qZmZmZmak/AAAAAAAAAADNzMzMzMzkP5qZmZmZmak/AAAAAAAAAABmZmZmZmbmP5qZmZmZmak/AAAAAAAAAAAAAAAAAADoP5qZmZmZmak/AAAAAAAAAACamZmZmZnpP5qZmZmZmak/AAAAAAAAAAAzMzMzMzPrP5qZmZmZmak/AAAAAAAAAADNzMzMzMzsP5qZmZmZmak/AAAAAAAAAABmZmZmZmbuP5qZmZmZmak/AAAAAAAAAAAAAAAAAADwP5qZmZmZmak/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmbk/AAAAAAAAAACamZmZmZmpP5qZmZmZmbk/AAAAAAAAAACamZmZmZm5P5qZmZmZmbk/AAAAAAAAAAAzMzMzMzPDP5qZmZmZmbk/AAAAAAAAAACamZmZmZnJP5qZmZmZmbk/AAAAAAAAAAAAAAAAAADQP5qZmZmZmbk/AAAAAAAAAAAzMzMzMzPTP5qZmZmZmbk/AAAAAAAAAABmZmZmZmbWP5qZmZmZmbk/AAAAAAAAAACamZmZmZnZP5qZmZmZmbk/AAAAAAAAAADNzMzMzMzcP5qZmZmZmbk/AAAAAAAAAAAAAAAAAADgP5qZmZmZmbk/AAAAAAAAAACamZmZmZnhP5qZmZmZmbk/AAAAAAAAAAAzMzMzMzPjP5qZmZmZmbk/AAAAAAAAAADNzMzMzMzkP5qZmZmZmbk/AAAAAAAAAABmZmZmZmbmP5qZmZmZmbk/AAAAAAAAAAAAAAAAAADoP5qZmZmZmbk/AAAAAAAAAACamZmZmZnpP5qZmZmZmbk/AAAAAAAAAAAzMzMzMzPrP5qZmZmZmbk/AAAAAAAAAADNzMzMzMzsP5qZmZmZmbk/AAAAAAAAAABmZmZmZmbuP5qZmZmZmbk/AAAAAAAAAAAAAAAAAADwP5qZmZmZmbk/AAAAAAAAAAA=
+ <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+ <Value index="0">
+ 0
+ </Value>
+ <Value index="1">
+ 1.0049875621
+ </Value>
+ </InformationKey>
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="62">
+ AAUAAAAAAAAAAAAAAAAAABUAAAAAAAAAFgAAAAAAAAABAAAAAAAAAAEAAAAAAAAAFgAAAAAAAAAXAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAAXAAAAAAAAABgAAAAAAAAAAwAAAAAAAAADAAAAAAAAABgAAAAAAAAAGQAAAAAAAAAEAAAAAAAAAAQAAAAAAAAAGQAAAAAAAAAaAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAaAAAAAAAAABsAAAAAAAAABgAAAAAAAAAGAAAAAAAAABsAAAAAAAAAHAAAAAAAAAAHAAAAAAAAAAcAAAAAAAAAHAAAAAAAAAAdAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAdAAAAAAAAAB4AAAAAAAAACQAAAAAAAAAJAAAAAAAAAB4AAAAAAAAAHwAAAAAAAAAKAAAAAAAAAAoAAAAAAAAAHwAAAAAAAAAgAAAAAAAAAAsAAAAAAAAACwAAAAAAAAAgAAAAAAAAACEAAAAAAAAADAAAAAAAAAAMAAAAAAAAACEAAAAAAAAAIgAAAAAAAAANAAAAAAAAAA0AAAAAAAAAIgAAAAAAAAAjAAAAAAAAAA4AAAAAAAAADgAAAAAAAAAjAAAAAAAAACQAAAAAAAAADwAAAAAAAAAPAAAAAAAAACQAAAAAAAAAJQAAAAAAAAAQAAAAAAAAABAAAAAAAAAAJQAAAAAAAAAmAAAAAAAAABEAAAAAAAAAEQAAAAAAAAAmAAAAAAAAACcAAAAAAAAAEgAAAAAAAAASAAAAAAAAACcAAAAAAAAAKAAAAAAAAAATAAAAAAAAABMAAAAAAAAAKAAAAAAAAAApAAAAAAAAABQAAAAAAAAAFQAAAAAAAAAqAAAAAAAAACsAAAAAAAAAFgAAAAAAAAAWAAAAAAAAACsAAAAAAAAALAAAAAAAAAAXAAAAAAAAABcAAAAAAAAALAAAAAAAAAAtAAAAAAAAABgAAAAAAAAAGAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAAGQAAAAAAAAAZAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAaAAAAAAAAABoAAAAAAAAALwAAAAAAAAAwAAAAAAAAABsAAAAAAAAAGwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAHAAAAAAAAAAcAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAdAAAAAAAAAB0AAAAAAAAAMgAAAAAAAAAzAAAAAAAAAB4AAAAAAAAAHgAAAAAAAAAzAAAAAAAAADQAAAAAAAAAHwAAAAAAAAAfAAAAAAAAADQAAAAAAAAANQAAAAAAAAAgAAAAAAAAACAAAAAAAAAANQAAAAAAAAA2AAAAAAAAACEAAAAAAAAAIQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAIgAAAAAAAAAiAAAAAAAAADcAAAAAAAAAOAAAAAAAAAAjAAAAAAAAACMAAAAAAAAAOAAAAAAAAAA5AAAAAAAAACQAAAAAAAAAJAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAJQAAAAAAAAAlAAAAAAAAADoAAAAAAAAAOwAAAAAAAAAmAAAAAAAAACYAAAAAAAAAOwAAAAAAAAA8AAAAAAAAACcAAAAAAAAAJwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAKAAAAAAAAAAoAAAAAAAAAD0AAAAAAAAAPgAAAAAAAAApAAAAAAAAAA==
+ </DataArray>
+ <DataArray type="Int64" Name="offsets" format="binary" RangeMin="4" RangeMax="160">
+ QAEAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAA==
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="binary" RangeMin="9" RangeMax="9">
+ KAAAAAAAAAAJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJ
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/Simple3DBbar/geometry_origin.vtu b/Tests/Data/Mechanics/Simple3DBbar/geometry_origin.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="1" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="0" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="24" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0" offset="48" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="92" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="116" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="140" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _CAAAAAAAAAAAAAAAAAAAAA==CAAAAAAAAAAAAAAAAAAAAA==GAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA=CAAAAAAAAAAAAAAAAAAAAA==CAAAAAAAAAABAAAAAAAAAA==AQAAAAAAAAAB
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Mechanics/Simple3DBbar/hex_brick.vtu b/Tests/Data/Mechanics/Simple3DBbar/hex_brick.vtu
diff --git a/Tests/Data/Mechanics/Simple3DBbar/simple_3d_bbar.prj b/Tests/Data/Mechanics/Simple3DBbar/simple_3d_bbar.prj
@@ -0,0 +1,172 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <meshes>
+ <mesh>hex_brick.vtu</mesh>
+ <mesh>x0.vtu</mesh>
+ <mesh>y0.vtu</mesh>
+ <mesh>top.vtu</mesh>
+ <mesh>bottom.vtu</mesh>
+ </meshes> 
+ <processes>
+ <process>
+ <name>SD</name>
+ <type>SMALL_DEFORMATION</type>
+ <integration_order>3</integration_order>
+ <use_b_bar>true</use_b_bar>
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <specific_body_force>0 0 0</specific_body_force>
+ <process_variables>
+ <process_variable>displacement</process_variable>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ </secondary_variables>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="SD">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltols> 1.0e-13 1.0e-13 1e-13</reltols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial> 0.0 </t_initial>
+ <t_end>1</t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>simple_3d_bbar</prefix>
+ <timesteps>
+ <pair>
+ <repeat> 1 </repeat>
+ <each_steps> 1 </each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>displacement</variable>
+ <variable>epsilon</variable>
+ <variable>sigma</variable>
+ </variables>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ </output>
+ </time_loop>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>density</name>
+ <type>Parameter</type>
+ <parameter_name>rho</parameter_name>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ </medium>
+ </media>
+ <parameters>
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1.0E+10</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.499</value>
+ </parameter>
+ <parameter>
+ <name>rho</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>zero</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>u0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>load</name>
+ <type>Constant</type>
+ <values>-10e+6</values>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>1</order>
+ <initial_condition>u0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>bottom</mesh>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>zero</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>x0</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>y0</mesh>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>zero</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>top</mesh>
+ <type>Neumann</type>
+ <component>2</component>
+ <parameter>load</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_newton</name>
+ <type>Newton</type>
+ <max_iter>10</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>
+</OpenGeoSysProject>