Merge pull request #12798 from KratosMultiphysics/servol_loading_with…

…_periodic_boundary [DEMApplication] [Fast] Servo-control algorithm with periodic boundaries
KratosMultiphysics · Oct 30, 2024 · ebfd53a · ebfd53a
2 parents 5289b98 + 1b0b1af
commit ebfd53a
Show file tree

Hide file tree

Showing 17 changed files with 515 additions and 28 deletions.
diff --git a/applications/DEMApplication/DEM_application.cpp b/applications/DEMApplication/DEM_application.cpp
@@ -101,6 +101,7 @@ KRATOS_CREATE_VARIABLE(int, ROTATION_OPTION)
 KRATOS_CREATE_VARIABLE(int, VIRTUAL_MASS_OPTION)
 KRATOS_CREATE_VARIABLE(int, SEARCH_CONTROL)
 KRATOS_CREATE_VARIABLE(bool, IS_TIME_TO_PRINT)
+KRATOS_CREATE_VARIABLE(bool, IS_TIME_TO_UPDATE_CONTACT_ELEMENT)
 KRATOS_CREATE_VARIABLE(double, COORDINATION_NUMBER)
 KRATOS_CREATE_VARIABLE(double, CONTINUUM_SEARCH_RADIUS_AMPLIFICATION_FACTOR)
 KRATOS_CREATE_VARIABLE(double, MAX_AMPLIFICATION_RATIO_OF_THE_SEARCH_RADIUS)
@@ -153,6 +154,7 @@ KRATOS_CREATE_VARIABLE(bool, DELTA_OPTION)
 KRATOS_CREATE_VARIABLE(double, SKIN_SPHERE)
 KRATOS_CREATE_VARIABLE(int, PROPERTIES_ID)
 KRATOS_CREATE_VARIABLE(int, CONTACT_MESH_OPTION)
+KRATOS_CREATE_VARIABLE(int, BOUNDING_BOX_SERVO_LOADING_OPTION)
 KRATOS_CREATE_VARIABLE(double, MAX_NUMBER_OF_INTACT_BONDS_TO_CONSIDER_A_SPHERE_BROKEN)
 //KRATOS_CREATE_VARIABLE(int, FAILURE_CRITERION_OPTION)
 KRATOS_CREATE_VARIABLE(int, CONCRETE_TEST_OPTION)
@@ -612,6 +614,7 @@ void KratosDEMApplication::Register() {
     KRATOS_REGISTER_VARIABLE(VIRTUAL_MASS_OPTION)
     KRATOS_REGISTER_VARIABLE(SEARCH_CONTROL)
     KRATOS_REGISTER_VARIABLE(IS_TIME_TO_PRINT)
+    KRATOS_REGISTER_VARIABLE(IS_TIME_TO_UPDATE_CONTACT_ELEMENT)
     KRATOS_REGISTER_VARIABLE(COORDINATION_NUMBER)
     KRATOS_REGISTER_VARIABLE(CONTINUUM_SEARCH_RADIUS_AMPLIFICATION_FACTOR)
     KRATOS_REGISTER_VARIABLE(MAX_AMPLIFICATION_RATIO_OF_THE_SEARCH_RADIUS)
@@ -660,6 +663,7 @@ void KratosDEMApplication::Register() {
     KRATOS_REGISTER_VARIABLE(SKIN_SPHERE)
     KRATOS_REGISTER_VARIABLE(PROPERTIES_ID)
     KRATOS_REGISTER_VARIABLE(CONTACT_MESH_OPTION)
+    KRATOS_REGISTER_VARIABLE(BOUNDING_BOX_SERVO_LOADING_OPTION)
     KRATOS_REGISTER_VARIABLE(MAX_NUMBER_OF_INTACT_BONDS_TO_CONSIDER_A_SPHERE_BROKEN)
     //KRATOS_REGISTER_VARIABLE(FAILURE_CRITERION_OPTION)
     KRATOS_REGISTER_VARIABLE(CONCRETE_TEST_OPTION)

diff --git a/applications/DEMApplication/DEM_application_variables.h b/applications/DEMApplication/DEM_application_variables.h
@@ -68,6 +68,7 @@ namespace Kratos
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, int, VIRTUAL_MASS_OPTION)
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, int, SEARCH_CONTROL)
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, bool, IS_TIME_TO_PRINT)
+  KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, bool, IS_TIME_TO_UPDATE_CONTACT_ELEMENT)
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, double, COORDINATION_NUMBER)
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, double, CONTINUUM_SEARCH_RADIUS_AMPLIFICATION_FACTOR)
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, double, MAX_AMPLIFICATION_RATIO_OF_THE_SEARCH_RADIUS)
@@ -125,6 +126,7 @@ namespace Kratos
 
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, int, PROPERTIES_ID)
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, int, CONTACT_MESH_OPTION)
+  KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, int, BOUNDING_BOX_SERVO_LOADING_OPTION)
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, double, MAX_NUMBER_OF_INTACT_BONDS_TO_CONSIDER_A_SPHERE_BROKEN)
   //KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, int, FAILURE_CRITERION_OPTION)
   KRATOS_DEFINE_APPLICATION_VARIABLE(DEM_APPLICATION, int, CONCRETE_TEST_OPTION)

diff --git a/applications/DEMApplication/custom_elements/spheric_particle.cpp b/applications/DEMApplication/custom_elements/spheric_particle.cpp
@@ -894,6 +894,11 @@ void SphericParticle::ComputeBallToBallContactForceAndMoment(SphericParticle::Pa
                 if ((i < (int)mNeighbourElements.size()) && this->Id() < neighbour_iterator_id) {
                     CalculateOnContactElements(i, LocalContactForce, GlobalContactForce);
                 }
+            } else if (r_process_info[IS_TIME_TO_UPDATE_CONTACT_ELEMENT] && r_process_info[CONTACT_MESH_OPTION] == 1) {
+                unsigned int neighbour_iterator_id = data_buffer.mpOtherParticle->Id();
+                if ((i < (int)mNeighbourElements.size()) && this->Id() < neighbour_iterator_id) {
+                    CalculateOnContactElements(i, LocalContactForce, GlobalContactForce);
+                }
             }
 
             if (r_process_info[ENERGY_CALCULATION_OPTION]){

diff --git a/applications/DEMApplication/custom_python/DEM_python_application.cpp b/applications/DEMApplication/custom_python/DEM_python_application.cpp
@@ -71,6 +71,7 @@ PYBIND11_MODULE(KratosDEMApplication,m)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, AUTOMATIC_SKIN_COMPUTATION)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, SKIN_FACTOR_RADIUS)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, IS_TIME_TO_PRINT)
+    KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, IS_TIME_TO_UPDATE_CONTACT_ELEMENT)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, COORDINATION_NUMBER)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, CONTINUUM_SEARCH_RADIUS_AMPLIFICATION_FACTOR)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, MAX_AMPLIFICATION_RATIO_OF_THE_SEARCH_RADIUS)
@@ -125,6 +126,7 @@ PYBIND11_MODULE(KratosDEMApplication,m)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, COHESIVE_GROUP)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, PROPERTIES_ID)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, CONTACT_MESH_OPTION)
+    KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, BOUNDING_BOX_SERVO_LOADING_OPTION)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, MAX_NUMBER_OF_INTACT_BONDS_TO_CONSIDER_A_SPHERE_BROKEN)
     //KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, FAILURE_CRITERION_OPTION)
     KRATOS_REGISTER_IN_PYTHON_VARIABLE(m, CONCRETE_TEST_OPTION)

diff --git a/applications/DEMApplication/custom_strategies/strategies/explicit_solver_strategy.cpp b/applications/DEMApplication/custom_strategies/strategies/explicit_solver_strategy.cpp
@@ -459,6 +459,9 @@ namespace Kratos {
         if (is_time_to_print_results && r_process_info[CONTACT_MESH_OPTION] == 1) {
             CreateContactElements();
             InitializeContactElements();
+        }else if (r_process_info[IS_TIME_TO_UPDATE_CONTACT_ELEMENT] && r_process_info[CONTACT_MESH_OPTION] == 1) {
+            CreateContactElements();
+            InitializeContactElements();
         }
 
         //RebuildPropertiesProxyPointers(mListOfSphericParticles);