[GeomechanicsApplication] Geo/linear elastic strategy

applications/GeoMechanicsApplication/custom_conditions/U_Pw_normal_lysmer_absorbing_condition.cpp

@@ -32,7 +32,7 @@ Condition::Pointer UPwLysmerAbsorbingCondition<TDim, TNumNodes>::Create(IndexTyp
 
 template <unsigned int TDim, unsigned int TNumNodes>
 void UPwLysmerAbsorbingCondition<TDim, TNumNodes>::CalculateLeftHandSide(MatrixType& rLeftHandSideMatrix,
-    const ProcessInfo& rCurrentProcessInfo)
+                                                                         const ProcessInfo& rCurrentProcessInfo)
 {
     ElementMatrixType stiffness_matrix;
 
@@ -41,7 +41,6 @@ void UPwLysmerAbsorbingCondition<TDim, TNumNodes>::CalculateLeftHandSide(MatrixT
     this->AddLHS(rLeftHandSideMatrix, stiffness_matrix);
 }
 
-
 template <unsigned int TDim, unsigned int TNumNodes>
 void UPwLysmerAbsorbingCondition<TDim, TNumNodes>::CalculateLocalSystem(MatrixType& rLhsMatrix,
                                                                         VectorType& rRightHandSideVector,

applications/GeoMechanicsApplication/custom_conditions/U_Pw_normal_lysmer_absorbing_condition.hpp

@@ -85,12 +85,11 @@ class KRATOS_API(GEO_MECHANICS_APPLICATION) UPwLysmerAbsorbingCondition
     void CalculateRightHandSide(VectorType& rRightHandSideVector, const ProcessInfo& rCurrentProcessInfo) override;
 
     /**
-    * @brief Calculates LHS stiffness part of absorbing boundary
-    * @param rLeftHandSideMatrix Global left hand side matrix
-    * @param rCurrentProcessInfo Current process information
-    */
-    void CalculateLeftHandSide(MatrixType& rLeftHandSideMatrix,
-        const ProcessInfo& rCurrentProcessInfo) override;
+     * @brief Calculates LHS stiffness part of absorbing boundary
+     * @param rLeftHandSideMatrix Global left hand side matrix
+     * @param rCurrentProcessInfo Current process information
+     */
+    void CalculateLeftHandSide(MatrixType& rLeftHandSideMatrix, const ProcessInfo& rCurrentProcessInfo) override;
 
     /**
      * @brief Calculates LHS Damping part of absorbing boundary

applications/GeoMechanicsApplication/custom_python/add_custom_strategies_to_python.cpp

@@ -22,9 +22,11 @@
 #include "custom_strategies/strategies/geo_mechanics_newton_raphson_erosion_process_strategy.hpp"
 #include "custom_strategies/strategies/geo_mechanics_newton_raphson_strategy.hpp"
 #include "custom_strategies/strategies/geo_mechanics_ramm_arc_length_strategy.hpp"
+#include "custom_strategies/strategies/residualbased_newton_raphson_strategy_linear_elastic_dynamic.hpp"
 #include "solving_strategies/strategies/solving_strategy.h"
 
 // builders and solvers
+#include "custom_strategies/builder_and_solvers/residualbased_block_builder_and_solver_linear_elastic_dynamic.h"
 #include "custom_strategies/builder_and_solvers/residualbased_block_builder_and_solver_with_mass_and_damping.h"
 
 // schemes
@@ -33,6 +35,7 @@
 #include "custom_strategies/schemes/backward_euler_quasistatic_U_Pw_scheme.hpp"
 #include "custom_strategies/schemes/generalized_newmark_T_scheme.hpp"
 #include "custom_strategies/schemes/newmark_dynamic_U_Pw_scheme.hpp"
+#include "custom_strategies/schemes/incremental_newmark_linear_elastic_U_scheme.hpp"
 #include "custom_strategies/schemes/newmark_quasistatic_Pw_scheme.hpp"
 #include "custom_strategies/schemes/newmark_quasistatic_U_Pw_scheme.hpp"
 #include "custom_strategies/schemes/newmark_quasistatic_damped_U_Pw_scheme.hpp"
@@ -60,6 +63,7 @@ void AddCustomStrategiesToPython(pybind11::module& m)
     using NewmarkQuasistaticDampedUPwSchemeType =
         NewmarkQuasistaticDampedUPwScheme<SparseSpaceType, LocalSpaceType>;
     using NewmarkDynamicUPwSchemeType = NewmarkDynamicUPwScheme<SparseSpaceType, LocalSpaceType>;
+	using IncrementalNewmarkLinearElasticUSchemeType = IncrementalNewmarkLinearElasticUScheme<SparseSpaceType, LocalSpaceType>;
     using NewmarkQuasistaticPwSchemeType = NewmarkQuasistaticPwScheme<SparseSpaceType, LocalSpaceType>;
     using NewmarkQuasistaticTSchemeType = GeneralizedNewmarkTScheme<SparseSpaceType, LocalSpaceType>;
 
@@ -75,6 +79,8 @@ void AddCustomStrategiesToPython(pybind11::module& m)
         GeoMechanicsRammArcLengthStrategy<SparseSpaceType, LocalSpaceType, LinearSolverType>;
     using GeoMechanicsNewtonRaphsonErosionProcessStrategyType =
         GeoMechanicsNewtonRaphsonErosionProcessStrategy<SparseSpaceType, LocalSpaceType, LinearSolverType>;
+    using GeoMechanicNewtonRaphsonStrategyLinearElasticDynamicType =
+        GeoMechanicNewtonRaphsonStrategyLinearElasticDynamic<SparseSpaceType, LocalSpaceType, LinearSolverType>;
 
     //----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
 
@@ -89,6 +95,10 @@ void AddCustomStrategiesToPython(pybind11::module& m)
     py::class_<NewmarkDynamicUPwSchemeType, typename NewmarkDynamicUPwSchemeType::Pointer, BaseSchemeType>(
         m, "NewmarkDynamicUPwScheme", py::module_local())
         .def(py::init<double, double, double>());
+
+	py::class_<IncrementalNewmarkLinearElasticUSchemeType, typename IncrementalNewmarkLinearElasticUSchemeType::Pointer, BaseSchemeType>(
+	    m, "IncrementalNewmarkLinearElasticUScheme", py::module_local())
+	    .def(py::init<double, double>());
 
     py::class_<NewmarkQuasistaticPwSchemeType, typename NewmarkQuasistaticPwSchemeType::Pointer, BaseSchemeType>(
         m, "NewmarkQuasistaticPwScheme")
@@ -127,13 +137,26 @@ void AddCustomStrategiesToPython(pybind11::module& m)
                       BuilderAndSolverType::Pointer, Parameters&, int, bool, bool, bool>())
         .def("UpdateLoads", &GeoMechanicsRammArcLengthStrategyType::UpdateLoads);
 
+    py::class_<GeoMechanicNewtonRaphsonStrategyLinearElasticDynamicType,
+               typename GeoMechanicNewtonRaphsonStrategyLinearElasticDynamicType::Pointer, BaseSolvingStrategyType>(
+        m, "GeoMechanicNewtonRaphsonStrategyLinearElasticDynamic")
+        .def(py::init<ModelPart&, BaseSchemeType::Pointer, LinearSolverType::Pointer, ConvergenceCriteriaType::Pointer,
+                      BuilderAndSolverType::Pointer, Parameters&, int, bool, bool>());
+
     using ResidualBasedBlockBuilderAndSolverWithMassAndDampingType =
         ResidualBasedBlockBuilderAndSolverWithMassAndDamping<SparseSpaceType, LocalSpaceType, LinearSolverType>;
     py::class_<ResidualBasedBlockBuilderAndSolverWithMassAndDampingType,
                ResidualBasedBlockBuilderAndSolverWithMassAndDampingType::Pointer, BuilderAndSolverType>(
         m, "ResidualBasedBlockBuilderAndSolverWithMassAndDamping")
         .def(py::init<LinearSolverType::Pointer>())
         .def(py::init<LinearSolverType::Pointer, Parameters>());
+
+    using ResidualBasedBlockBuilderAndSolverLinearElasticDynamicType =
+        ResidualBasedBlockBuilderAndSolverLinearElasticDynamic<SparseSpaceType, LocalSpaceType, LinearSolverType>;
+    py::class_<ResidualBasedBlockBuilderAndSolverLinearElasticDynamicType,
+               ResidualBasedBlockBuilderAndSolverLinearElasticDynamicType::Pointer, BuilderAndSolverType>(
+        m, "ResidualBasedBlockBuilderAndSolverLinearElasticDynamic")
+        .def(py::init<LinearSolverType::Pointer, double, double, bool>());
 }
 
 } // Namespace Kratos::Python