ufz · endJunction · Sep 7, 2016 · Sep 5, 2016 · Sep 7, 2016 · Jul 4, 2016
diff --git a/Applications/ApplicationsLib/ProjectData.cpp b/Applications/ApplicationsLib/ProjectData.cpp
@@ -32,6 +32,7 @@
 #include "ProcessLib/UncoupledProcessesTimeLoop.h"
 
 #include "ProcessLib/GroundwaterFlow/CreateGroundwaterFlowProcess.h"
+#include "ProcessLib/SmallDeformation/CreateSmallDeformationProcess.h"
 #include "ProcessLib/TES/CreateTESProcess.h"
 #include "ProcessLib/HeatConduction/CreateHeatConductionProcess.h"
 
@@ -269,6 +270,28 @@ void ProjectData::parseProcesses(BaseLib::ConfigTree const& processes_config)
             process = ProcessLib::HeatConduction::createHeatConductionProcess(
                 *_mesh_vec[0], _process_variables, _parameters, process_config);
         }
+        else if (type == "SMALL_DEFORMATION")
+        {
+            switch (process_config.getConfigParameter<int>("dimension"))
+            {
+                case 2:
+                    process = ProcessLib::SmallDeformation::
+                        createSmallDeformationProcess<2>(
+                            *_mesh_vec[0], _process_variables, _parameters,
+                            process_config);
+                    break;
+                case 3:
+                    process = ProcessLib::SmallDeformation::
+                        createSmallDeformationProcess<3>(
+                            *_mesh_vec[0], _process_variables, _parameters,
+                            process_config);
+                    break;
+                default:
+                    OGS_FATAL(
+                        "SMALL_DEFORMATION process does not support "
+                        "given dimension");
+            }
+        }
         else
         {
             OGS_FATAL("Unknown process type: %s", type.c_str());

diff --git a/Applications/CLI/Tests.cmake b/Applications/CLI/Tests.cmake
@@ -243,6 +243,52 @@ if(NOT OGS_USE_MPI)
          temperature_analytical.vtu line_60_heat_pcs_0_ts_405_t_31640625.000000.vtu Temperature_Analytical_1year temperature
     )
 
+    # Mechanics; Small deformations, linear (SDL)
+    AddTest(
+        NAME Mechanics_SDL_square_1e0_displacementBC
+        PATH Mechanics/Linear
+        EXECUTABLE ogs
+        EXECUTABLE_ARGS square_1e0.prj
+        WRAPPER time
+        TESTER vtkdiff
+        ABSTOL 1e-16 RELTOL 1e-16
+        DIFF_DATA
+        square_1e0_expected_pcs_0_ts_4_t_1.000000.vtu square_1e0_pcs_0_ts_4_t_1.000000.vtu displacement displacement
+    )
+    AddTest(
+        NAME Mechanics_SDL_square_1e2_tractionBC
+        PATH Mechanics/Linear
+        EXECUTABLE ogs
+        EXECUTABLE_ARGS square_1e2.prj
+        WRAPPER time
+        TESTER vtkdiff
+        ABSTOL 1e-16 RELTOL 1e-16
+        DIFF_DATA
+        square_1e2_expected_pcs_0_ts_4_t_1.000000.vtu square_1e2_pcs_0_ts_4_t_1.000000.vtu displacement displacement
+    )
+    AddTest(
+        NAME LARGE_Mechanics_SDL_disc_with_hole
+        PATH Mechanics/Linear
+        EXECUTABLE ogs
+        EXECUTABLE_ARGS disc_with_hole.prj
+        WRAPPER time
+        TESTER vtkdiff
+        ABSTOL 1e-16 RELTOL 1e-16
+        DIFF_DATA
+        disc_with_hole_expected_pcs_0_ts_4_t_1.000000.vtu disc_with_hole_pcs_0_ts_4_t_1.000000.vtu displacement displacement
+    )
+    AddTest(
+        NAME LARGE_Mechanics_SDL_square_1e5_tractionBC
+        PATH Mechanics/Linear
+        EXECUTABLE ogs
+        EXECUTABLE_ARGS square_1e5.prj
+        WRAPPER time
+        TESTER vtkdiff
+        ABSTOL 1e-16 RELTOL 1e-16
+        DIFF_DATA
+        square_1e5_expected_pcs_0_ts_4_t_1.000000.vtu square_1e5_pcs_0_ts_4_t_1.000000.vtu displacement displacement
+    )
+
 else()
     # MPI groundwater flow tests
     AddTest(

diff --git a/Documentation/ProjectFile/process/SMALL_DEFORMATION/c_SMALL_DEFORMATION.md b/Documentation/ProjectFile/process/SMALL_DEFORMATION/c_SMALL_DEFORMATION.md
@@ -0,0 +1 @@
+Deformation process with linear kinematics.
diff --git a/...MATION/constitutive_relation/LinearElasticIsotropic/c_LinearElasticIsotropic.md b/...MATION/constitutive_relation/LinearElasticIsotropic/c_LinearElasticIsotropic.md
@@ -0,0 +1 @@
+Isotropic linear elastic material model.
diff --git a/...LL_DEFORMATION/constitutive_relation/LinearElasticIsotropic/t_poissons_ratio.md b/...LL_DEFORMATION/constitutive_relation/LinearElasticIsotropic/t_poissons_ratio.md
@@ -0,0 +1 @@
+Poisson's ratio of a linear elastic material.
diff --git a/...LL_DEFORMATION/constitutive_relation/LinearElasticIsotropic/t_youngs_modulus.md b/...LL_DEFORMATION/constitutive_relation/LinearElasticIsotropic/t_youngs_modulus.md
@@ -0,0 +1 @@
+Young's modulus of a linear elastic material.
diff --git a/...File/process/SMALL_DEFORMATION/constitutive_relation/i_constitutive_relation.md b/...File/process/SMALL_DEFORMATION/constitutive_relation/i_constitutive_relation.md
diff --git a/.../ProjectFile/process/SMALL_DEFORMATION/process_variables/i_process_variables.md b/.../ProjectFile/process/SMALL_DEFORMATION/process_variables/i_process_variables.md
diff --git a/...n/ProjectFile/process/SMALL_DEFORMATION/process_variables/t_process_variable.md b/...n/ProjectFile/process/SMALL_DEFORMATION/process_variables/t_process_variable.md
@@ -0,0 +1 @@
+The displacement vector field.
diff --git a/Documentation/ProjectFile/process/SMALL_DEFORMATION/t_dimension.md b/Documentation/ProjectFile/process/SMALL_DEFORMATION/t_dimension.md
@@ -0,0 +1 @@
+Displacement vector dimension. The displacement dimension must be equal to the mesh dimension.
diff --git a/MaterialLib/SolidModels/CreateLinearElasticIsotropic.h b/MaterialLib/SolidModels/CreateLinearElasticIsotropic.h
@@ -0,0 +1,49 @@
+/**
+ * \copyright
+ * Copyright (c) 2012-2016, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#ifndef MATERIALLIB_SOLIDMODELS_CREATELINEARELASTICISOTROPIC_H_
+#define MATERIALLIB_SOLIDMODELS_CREATELINEARELASTICISOTROPIC_H_
+
+#include "ProcessLib/Utils/ProcessUtils.h"  // required for findParameter
+#include "LinearElasticIsotropic.h"
+
+namespace Solids
+{
+
+template <int DisplacementDim>
+std::unique_ptr<LinearElasticIsotropic<DisplacementDim>>
+createLinearElasticIsotropic(
+    std::vector<std::unique_ptr<ProcessLib::ParameterBase>> const& parameters,
+    BaseLib::ConfigTree const& config)
+{
+    config.checkConfigParameter("type", "LinearElasticIsotropic");
+    DBUG("Create LinearElasticIsotropic material");
+
+    // Youngs modulus
+    auto& youngs_modulus = ProcessLib::findParameter<double>(
+        config, "youngs_modulus", parameters, 1);
+
+    DBUG("Use '%s' as youngs_modulus parameter.", youngs_modulus.name.c_str());
+
+    // Poissons ratio
+    auto& poissons_ratio = ProcessLib::findParameter<double>(
+        config, "poissons_ratio", parameters, 1);
+
+    DBUG("Use '%s' as poissons_ratio parameter.", poissons_ratio.name.c_str());
+
+    typename LinearElasticIsotropic<DisplacementDim>::MaterialProperties mp{
+        youngs_modulus, poissons_ratio};
+
+    return std::unique_ptr<LinearElasticIsotropic<DisplacementDim>>{
+        new LinearElasticIsotropic<DisplacementDim>{mp}};
+}
+
+}  // namespace Solids
+
+#endif  // MATERIALLIB_SOLIDMODELS_CREATELINEARELASTICISOTROPIC_H_
diff --git a/MaterialLib/SolidModels/LinearElasticIsotropic.cpp b/MaterialLib/SolidModels/LinearElasticIsotropic.cpp
@@ -0,0 +1,19 @@
+/**
+ * \copyright
+ * Copyright (c) 2012-2016, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#include "LinearElasticIsotropic.h"
+
+namespace Solids
+{
+
+template class LinearElasticIsotropic<2>;
+template class LinearElasticIsotropic<3>;
+
+
+}   // namespace Solids
diff --git a/MaterialLib/SolidModels/LinearElasticIsotropic.h b/MaterialLib/SolidModels/LinearElasticIsotropic.h
@@ -0,0 +1,112 @@
+/**
+ * \copyright
+ * Copyright (c) 2012-2016, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#ifndef MATERIALLIB_SOLIDMODELS_LINEARELASTICISOTROPIC_H_
+#define MATERIALLIB_SOLIDMODELS_LINEARELASTICISOTROPIC_H_
+
+#include "MechanicsBase.h"
+
+namespace Solids
+{
+template <int DisplacementDim>
+class LinearElasticIsotropic final : public MechanicsBase<DisplacementDim>
+{
+public:
+    /// Variables specific to the material model
+    class MaterialProperties
+    {
+        using P = ProcessLib::Parameter<double>;
+        using X = ProcessLib::SpatialPosition;
+
+    public:
+        MaterialProperties(P const& youngs_modulus, P const& poissons_ratio)
+            : _youngs_modulus(youngs_modulus), _poissons_ratio(poissons_ratio)
+        {
+        }
+
+        /// Lamé's first parameter.
+        double lambda(double const t, X const& x) const
+        {
+            return _youngs_modulus(t, x)[0] * _poissons_ratio(t, x)[0] /
+                   (1 + _poissons_ratio(t, x)[0]) /
+                   (1 - 2 * _poissons_ratio(t, x)[0]);
+        }
+
+        /// Lamé's second parameter, the shear modulus.
+        double mu(double const t, X const& x) const
+        {
+            return _youngs_modulus(t, x)[0] /
+                   (2 * (1 + _poissons_ratio(t, x)[0]));
+        }
+
+    private:
+        P const& _youngs_modulus;
+        P const& _poissons_ratio;
+    };
+
+    struct MaterialStateVariables
+        : public MechanicsBase<DisplacementDim>::MaterialStateVariables
+    {
+        void pushBackState() {}
+    };
+
+    std::unique_ptr<
+        typename MechanicsBase<DisplacementDim>::MaterialStateVariables>
+    createMaterialStateVariables() override
+    {
+        return std::unique_ptr<
+            typename MechanicsBase<DisplacementDim>::MaterialStateVariables>{
+            new MaterialStateVariables};
+    }
+
+public:
+    static int const KelvinVectorSize =
+        ProcessLib::KelvinVectorDimensions<DisplacementDim>::value;
+    using KelvinVector = ProcessLib::KelvinVectorType<DisplacementDim>;
+    using KelvinMatrix = ProcessLib::KelvinMatrixType<DisplacementDim>;
+
+    explicit LinearElasticIsotropic(
+        MaterialProperties const& material_properties)
+        : _mp(material_properties)
+    {
+    }
+
+    void computeConstitutiveRelation(
+        double const t,
+        ProcessLib::SpatialPosition const& x,
+        double const /*dt*/,
+        KelvinVector const& eps_prev,
+        KelvinVector const& eps,
+        KelvinVector const& sigma_prev,
+        KelvinVector& sigma,
+        KelvinMatrix& C,
+        typename MechanicsBase<DisplacementDim>::MaterialStateVariables&
+        /*material_state_variables*/) override
+    {
+        C.setZero();
+
+        C.template topLeftCorner<3, 3>().setConstant(_mp.lambda(t, x));
+        C.noalias() += 2 * _mp.mu(t, x) * KelvinMatrix::Identity();
+
+        sigma.noalias() = sigma_prev + C * (eps - eps_prev);
+    }
+
+private:
+    MaterialProperties _mp;
+};
+
+}  // namespace Solids
+
+namespace Solids
+{
+extern template class LinearElasticIsotropic<2>;
+extern template class LinearElasticIsotropic<3>;
+}  // namespace Solids
+
+#endif  // MATERIALLIB_SOLIDMODELS_LINEARELASTICISOTROPIC_H_
diff --git a/MaterialLib/SolidModels/MechanicsBase.h b/MaterialLib/SolidModels/MechanicsBase.h
@@ -0,0 +1,110 @@
+/**
+ * \copyright
+ * Copyright (c) 2012-2016, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ */
+
+#ifndef MATERIALLIB_SOLIDMODELS_MECHANICSBASE_H_
+#define MATERIALLIB_SOLIDMODELS_MECHANICSBASE_H_
+
+#include <memory>
+
+#include "ProcessLib/Deformation/BMatrixPolicy.h"
+#include "ProcessLib/Parameter/Parameter.h"
+
+namespace MeshLib
+{
+class Element;
+}
+
+namespace Solids
+{
+/// Interface for mechanical solid material models. Provides updates of the
+/// stress for a given current state and also a tangent at that position. If the
+/// implemented material model stores an internal state, the nested
+/// MaterialStateVariables class should be used; it's only responsibility is to
+/// provide state's push back possibility.
+template <int DisplacementDim>
+struct MechanicsBase
+{
+    /// The MaterialStateVariables may store material model specific state
+    /// (other than sigma and eps), which are usually material history
+    /// dependent. The objects are stored by the user (usually in assembly per
+    /// integration point) and are created via \ref
+    /// createMaterialStateVariables().
+    struct MaterialStateVariables
+    {
+        virtual ~MaterialStateVariables() = default;
+        virtual void pushBackState() = 0;
+    };
+
+    /// Polymorphic creator for MaterialStateVariables objects specific for a
+    /// material model.
+    virtual std::unique_ptr<MaterialStateVariables>
+    createMaterialStateVariables() = 0;
+
+    using KelvinVector = ProcessLib::KelvinVectorType<DisplacementDim>;
+    using KelvinMatrix = ProcessLib::KelvinMatrixType<DisplacementDim>;
+
+    /// Dynamic size Kelvin vector and matrix wrapper for the polymorphic
+    /// constitutive relation compute function.
+    void computeConstitutiveRelation(
+        double const t,
+        ProcessLib::SpatialPosition const& x,
+        double const dt,
+        Eigen::Matrix<double, Eigen::Dynamic, 1> const& eps_prev,
+        Eigen::Matrix<double, Eigen::Dynamic, 1> const& eps,
+        Eigen::Matrix<double, Eigen::Dynamic, 1> const& sigma_prev,
+        Eigen::Matrix<double, Eigen::Dynamic, 1>& sigma,
+        Eigen::Matrix<double, Eigen::Dynamic, Eigen::Dynamic, Eigen::RowMajor>&
+            C,
+        MaterialStateVariables& material_state_variables)
+    {
+        // TODO Avoid copies of data:
+        // Using MatrixBase<Derived> not possible because template functions
+        // cannot be virtual. Maybe there is a workaround for this.  Using
+        // Map<Matrix<double, ...>> makes the interface (for the material model
+        // implementation) unnecessary difficult.
+        KelvinVector const eps_prev_{eps_prev};
+        KelvinVector const eps_{eps};
+        KelvinVector const sigma_prev_{sigma_prev};
+        KelvinVector sigma_{sigma};
+        KelvinMatrix C_{C};
+
+        computeConstitutiveRelation(t,
+                                    x,
+                                    dt,
+                                    eps_prev_,
+                                    eps_,
+                                    sigma_prev_,
+                                    sigma_,
+                                    C_,
+                                    material_state_variables);
+
+        sigma = sigma_;
+        C = C_;
+    }
+
+    /// Computation of the constitutive relation for specific material model.
+    /// This should be implemented in the derived model. Fixed Kelvin vector and
+    /// matrix size version; for dynamic size arguments there is an overloaded
+    /// wrapper function.
+    virtual void computeConstitutiveRelation(
+        double const t,
+        ProcessLib::SpatialPosition const& x,
+        double const dt,
+        KelvinVector const& eps_prev,
+        KelvinVector const& eps,
+        KelvinVector const& sigma_prev,
+        KelvinVector& sigma,
+        KelvinMatrix& C,
+        MaterialStateVariables& material_state_variables) = 0;
+
+    virtual ~MechanicsBase() = default;
+};
+}
+
+#endif  // MATERIALLIB_SOLIDMODELS_MECHANICSBASE_H_
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		Deformation process with linear kinematics.
Copy link Collaborator chleh Aug 31, 2016 • edited Loading Choose a reason for hiding this comment The reason will be displayed to describe this comment to others. Learn more. dynamics rather than kinematics? (since forces are involved.) Or maybe just say, that it's for infinitesimal deformations? Copy link Member nagelt Aug 31, 2016 Choose a reason for hiding this comment The reason will be displayed to describe this comment to others. Learn more. linear kinematics implies small deformations.
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		Displacement vector dimension. The displacement dimension must be equal to the mesh dimension.