Skip to content
Snippets Groups Projects
Commit f3da8c5d authored by Reza Taherdangkoo's avatar Reza Taherdangkoo
Browse files

HM_drainage 

new modifications

new modification

new

last

last

last1

last1

last2

new
parent 33d7a7ae
No related branches found
No related tags found
No related merge requests found
Expand all Hide whitespace changes
Inline Side-by-side
Showing
with 523 additions and 1 deletion
ProcessLib/HydroMechanics/Tests.cmake
+ 3
− 1
View file @ f3da8c5d
......@@ -4,7 +4,9 @@
if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE HydroMechanics/Linear/Confined_Compression/square_1e2.prj)
endif()
if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE HydroMechanics/Linear/DrainageEexcavation/HMdrainage.prj)
endif()
AddTest(
NAME HydroMechanics_HML_square_1e2_quad9_confined_compression
PATH HydroMechanics/Linear/Confined_Compression
......
Tests/Data/HydroMechanics/Linear/DrainageEexcavation/HMdrainage.prj 0 → 100644
+ 335
− 0
View file @ f3da8c5d
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<mesh>drainage_quad.vtu</mesh>
<geometry>drainage_geo.gml</geometry>
<processes>
<process>
<name>HM</name>
<type>HYDRO_MECHANICS</type>
<integration_order>4</integration_order>
<dimension>2</dimension>
<constitutive_relation>
<type>LinearElasticIsotropic</type>
<youngs_modulus>E</youngs_modulus>
<poissons_ratio>nu</poissons_ratio>
</constitutive_relation>
<process_variables>
<displacement>displacement</displacement>
<pressure>pressure</pressure>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="sigma_xx" output_name="sigma_xx"/>
<secondary_variable internal_name="sigma_yy" output_name="sigma_yy"/>
<secondary_variable internal_name="sigma_zz" output_name="sigma_zz"/>
<secondary_variable internal_name="sigma_xy" output_name="sigma_xy"/>
<secondary_variable internal_name="epsilon_xx" output_name="epsilon_xx"/>
<secondary_variable internal_name="epsilon_yy" output_name="epsilon_yy"/>
<secondary_variable internal_name="epsilon_zz" output_name="epsilon_zz"/>
<secondary_variable internal_name="epsilon_xy" output_name="epsilon_xy"/>
<secondary_variable internal_name="velocity" output_name="velocity"/>
</secondary_variables>
<specific_body_force>0 -9.81</specific_body_force>
</process>
</processes>
<media>
<medium>
<phases>
<phase>
<type>Gas</type>
<properties>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1e-3</value>
</property>
<property>
<name>density</name>
<type>Parameter</type>
<parameter_name>load-water</parameter_name>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0.35</value>
</property>
<property>
<name>density</name>
<type>Parameter</type>
<parameter_name>load-solid</parameter_name>
</property>
<property>
<name>biot_coefficient</name>
<type>Constant</type>
<value>1</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>reference_temperature</name>
<type>Constant</type>
<value>293.15</value>
</property>
<property>
<name>permeability</name>
<type>Constant</type>
<value>1e-12</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="HM">
<nonlinear_solver>basic_newton</nonlinear_solver>
<convergence_criterion>
<type>PerComponentDeltaX</type>
<norm_type>NORM2</norm_type>
<reltols>1e-12 1e-10 1e-12</reltols>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0</t_initial>
<t_end>1e5</t_end>
<timesteps>
<pair>
<repeat>10</repeat>
<delta_t>10</delta_t>
</pair>
<pair>
<repeat>50</repeat>
<delta_t>100</delta_t>
</pair>
<pair>
<repeat>50</repeat>
<delta_t>1000</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>HM_drainage</prefix>
<timesteps>
<pair>
<repeat>1</repeat>
<each_steps>100</each_steps>
</pair>
<pair>
<repeat>1000</repeat>
<each_steps>10</each_steps>
</pair>
</timesteps>
<variables>
<variable>displacement</variable>
<variable>pressure</variable>
<variable>sigma_xx</variable>
<variable>sigma_yy</variable>
<variable>sigma_zz</variable>
<variable>sigma_xy</variable>
<variable>epsilon_xx</variable>
<variable>epsilon_yy</variable>
<variable>epsilon_zz</variable>
<variable>epsilon_xy</variable>
<variable>velocity</variable>
</variables>
<suffix>_ts_{:timestep}_t_{:time}</suffix>
</output>
</time_loop>
<parameters>
<!-- Mechanics -->
<parameter>
<name>E</name>
<type>Constant</type>
<value>3e9</value>
</parameter>
<parameter>
<name>nu</name>
<type>Constant</type>
<value>.3</value>
</parameter>
<!-- Model parameters -->
<parameter>
<name>displacement0</name>
<type>Constant</type>
<values>0 0</values>
</parameter>
<parameter>
<name>zero</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>pressure_w</name>
<type>Function</type>
<expression>1000*9.81*30</expression>
</parameter>
<parameter>
<name>pressure_top</name>
<type>CurveScaled</type>
<curve>time_ramp</curve>
<parameter>pressure_w</parameter>
</parameter>
<parameter>
<name>load_min</name>
<type>Function</type>
<expression>-((2600*9.81*30*0.65)+1000*9.81*30*0.35)</expression>
</parameter>
<parameter>
<name>load_top</name>
<type>CurveScaled</type>
<curve>time_ramp</curve>
<parameter>load_min</parameter>
</parameter>
<parameter>
<name>maximum_density</name>
<type>Constant</type>
<value>2600</value>
</parameter>
<parameter>
<name>load-solid</name>
<type>CurveScaled</type>
<curve>material_ramp</curve>
<parameter>maximum_density</parameter>
</parameter>
<parameter>
<name>maximum_density_water</name>
<type>Constant</type>
<value>1000</value>
</parameter>
<parameter>
<name>load-water</name>
<type>CurveScaled</type>
<curve>material_ramp</curve>
<parameter>maximum_density_water</parameter>
</parameter>
</parameters>
<curves>
<curve>
<name>time_ramp</name>
<coords>0 100 15551</coords>
<values>1 1 0.0</values>
</curve>
<curve>
<name>material_ramp</name>
<coords>0.0 100.0 15551</coords>
<values>0.0 1.0 1.0</values>
</curve>
</curves>
<process_variables>
<process_variable>
<name>displacement</name>
<components>2</components>
<order>2</order>
<initial_condition>displacement0</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>HH_geometry</geometrical_set>
<geometry>left</geometry>
<type>Dirichlet</type>
<component>0</component>
<parameter>zero</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>HH_geometry</geometrical_set>
<geometry>left_inner</geometry>
<type>Dirichlet</type>
<component>0</component>
<parameter>zero</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>HH_geometry</geometrical_set>
<geometry>right_inner</geometry>
<type>Dirichlet</type>
<component>0</component>
<parameter>zero</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>HH_geometry</geometrical_set>
<geometry>right</geometry>
<type>Dirichlet</type>
<component>0</component>
<parameter>zero</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>HH_geometry</geometrical_set>
<geometry>bottom</geometry>
<type>Dirichlet</type>
<component>1</component>
<parameter>zero</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>HH_geometry</geometrical_set>
<geometry>top_right</geometry>
<type>Neumann</type>
<component>1</component>
<parameter>load_top</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<process_variable>
<name>pressure</name>
<components>1</components>
<order>1</order>
<initial_condition>zero</initial_condition>
<boundary_conditions>
<boundary_condition>
<geometrical_set>HH_geometry</geometrical_set>
<geometry>top_left</geometry>
<type>Dirichlet</type>
<component>0</component>
<parameter>zero</parameter>
</boundary_condition>
<boundary_condition>
<geometrical_set>HH_geometry</geometrical_set>
<geometry>top_right</geometry>
<type>Dirichlet</type>
<component>0</component>
<parameter>pressure_top</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_newton</name>
<type>Newton</type>
<max_iter>50</max_iter>
<linear_solver>general_linear_solver</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>general_linear_solver</name>
<lis>-i bicg -p jacobi -tol 1e-16 -maxiter 10000</lis>
<!--eigen>
<solver_type>BiCGSTAB</solver_type>
<precon_type>ILUT</precon_type>
<max_iteration_step>10000</max_iteration_step>
<error_tolerance>1e-16</error_tolerance>
</eigen-->
<eigen>
<solver_type>SparseLU</solver_type>
<scaling>1</scaling>
</eigen>
</linear_solver>
</linear_solvers>
<test_definition>
<vtkdiff>
<file>HM_drainage_ts_155_t_100000.000000.vtu</file>
<field>displacement</field>
<absolute_tolerance>5e-8</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
</test_definition>
</OpenGeoSysProject>
Tests/Data/HydroMechanics/Linear/DrainageEexcavation/drainage_geo.gml 0 → 100644
+ 51
− 0
View file @ f3da8c5d
<?xml version="1.0" encoding="ISO-8859-1"?>
<?xml-stylesheet type="text/xsl" href="OpenGeoSysGLI.xsl"?>
<OpenGeoSysGLI xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ogs="http://www.opengeosys.org">
<name>HH_geometry</name>
<points>
<point id="0" x="0" y="0" z="0"/>
<point id="1" x="0" y="60" z="0"/>
<point id="2" x="35" y="60" z="0"/>
<point id="3" x="35" y="20" z="0"/>
<point id="4" x="36" y="20" z="0"/>
<point id="5" x="36" y="30" z="0"/>
<point id="6" x="65" y="30" z="0"/>
<point id="7" x="65" y="0" z="0"/>
</points>
<polylines>
<polyline id="0" name="left">
<pnt>0</pnt>
<pnt>1</pnt>
</polyline>
<polyline id="1" name="left_inner">
<pnt>2</pnt>
<pnt>3</pnt>
</polyline>
<polyline id="2" name="right_inner">
<pnt>4</pnt>
<pnt>5</pnt>
</polyline>
<polyline id="3" name="right">
<pnt>6</pnt>
<pnt>7</pnt>
</polyline>
<polyline id="4" name="top_left">
<pnt>1</pnt>
<pnt>2</pnt>
</polyline>
<polyline id="5" name="top_right">
<pnt>5</pnt>
<pnt>6</pnt>
</polyline>
<polyline id="6" name="bottom_inner">
<pnt>3</pnt>
<pnt>4</pnt>
</polyline>
<polyline id="7" name="bottom">
<pnt>0</pnt>
<pnt>7</pnt>
</polyline>
</polylines>
</OpenGeoSysGLI>
Tests/Data/HydroMechanics/Linear/DrainageEexcavation/drainage_quad.vtu 0 → 100644
+ 22
− 0
View file @ f3da8c5d
This diff is collapsed.
web/content/docs/benchmarks/hydro-mechanics/HM_drainage_excavation.md 0 → 100644
+ 18
− 0
View file @ f3da8c5d
+++
author = "Reza Taherdangkoo, Thomas Nagel, Markus Barsch, Christoph Butscher"
date = "2020.11.23T14:39:39+01:00"
title = "Hydro-mechanical-drainage_excavation"
project = "HydroMechanics/Linear/DrainageEexcavation/HMdrainage.prj"
[menu]
[menu.benchmarks]
parent = "hydro-mechanics"
+++
{{< data-link >}}
## Problem description
We present the hydro-mechanical drainage of an excavation benchmark in this [this PDF](../HM_drainage_excavation.pdf).
web/content/docs/benchmarks/hydro-mechanics/HM_drainage_excavation.pdf 0 → 100644
+ 0
− 0
View file @ f3da8c5d
File added
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