Merge pull request #671 from KratosMultiphysics/cdem/automatic-bonded-domain-groups

[CDEM] Added procedure to create and assign bonded domains to all volumes

Author: ipouplana

kratos.gid/apps/CDEM/examples/examples.tcl

@@ -9,6 +9,7 @@ proc CDEM::examples::Init { } {
 proc CDEM::examples::UpdateMenus { } {
     GiDMenu::InsertOption "Kratos" [list "---"] 8 PRE "" "" "" insertafter =
     GiDMenu::InsertOption "Kratos" [list "Stone block and sand" ] 8 PRE [list ::CDEM::examples::ContinuumDrop] "" "" insertafter =
+    GiDMenu::InsertOption "Kratos" [list "Bulk groups" ] 8 PRE [list ::CDEM::xml::BulkGroup] "" "" insertafter =
     GiDMenu::UpdateMenus
 }
 

kratos.gid/apps/CDEM/images/stone.png

@@ -36,6 +36,7 @@ proc ::CDEM::LoadMyFiles { } {
     variable dir
 
     uplevel #0 [list source [file join $dir xml GetFromXML.tcl]]
+    uplevel #0 [list source [file join $dir xml BulkGroup.tcl]]
     uplevel #0 [list source [file join $dir write write.tcl]]
     uplevel #0 [list source [file join $dir write writeMDPA_Parts.tcl]]
     uplevel #0 [list source [file join $dir write writeMDPA_Walls.tcl]]
@@ -63,6 +64,7 @@ proc ::CDEM::CustomToolbarItems { } {
     if {$::Model::SpatialDimension eq "2D"} {
         Kratos::ToolbarAddItem "Example" [file join $dir images drop.png] [list -np- ::CDEM::examples::ContinuumDrop] [= "Example\nRocks fall"]
     }
+    Kratos::ToolbarAddItem "Bulk grouping" [file join $dir images stone.png] [list -np- ::CDEM::xml::BulkGroup] [= "Plugin\nBulk grouping"]
 }
 
 ::CDEM::Init

kratos.gid/apps/CDEM/start.tcl

@@ -0,0 +1,58 @@
+
+proc ::CDEM::xml::BulkGroup {args} {
+    
+
+    if {![Kratos::IsModelEmpty]} {
+        set entity surfaces
+        if {$::Model::SpatialDimension eq "3D"} {
+            set entity volumes
+        }
+        set txt "Bulk grouping of all the current $entity is about to be executed"
+        set retval [tk_messageBox -default ok -icon question -message $txt -type okcancel]
+        if { $retval == "cancel" } { return }
+    }
+
+    CreateAndAssignBondedGroups
+    
+    GiD_Process 'Redraw
+    GidUtils::UpdateWindow GROUPS
+    GidUtils::UpdateWindow LAYER
+    GiD_Process 'Zoom Frame
+    
+    spdAux::RequestRefresh
+}
+
+proc ::CDEM::xml::CreateAndAssignBondedGroups { } {
+
+    # Prepare the variables
+    set entity surface
+    set condition_id "DEM-Cohesive2D"
+    if {$::Model::SpatialDimension eq "3D"} {
+        set entity volume
+        set condition_id "DEM-Cohesive"
+    }
+    set group_base Bonded_domain_groups
+    # Locate the condition
+    set DEMConditions [spdAux::getRoute "DEMConditions"]
+    set cohesive_cond "$DEMConditions/condition\[@n='$condition_id'\]"
+
+    # Clear previous assignations
+    if {[GiD_Groups exists $group_base]} {GiD_Groups delete $group_base}
+    gid_groups_conds::delete "$cohesive_cond/group"
+
+    # Get the list of entities (surfaces/volumes)
+    set entity_list [GiD_Geometry -v2 list $entity]
+    # Foreach entity, create subgroup and assign it to the condition
+    set entity_list_length [objarray length $entity_list]
+    for {set i 0} {$i < $entity_list_length} {incr i} {
+        # Get the entity id
+        set entity_id [objarray get $entity_list $i] 
+        # Create the subgroup
+        set group_id "$group_base//SG$entity_id"
+        GiD_Groups create $group_id
+        # Assign entity to subgroup
+        GiD_EntitiesGroups assign $group_id $entity $entity_id
+        # Assign the subgroup to the tree condition
+        set cohesive_group [customlib::AddConditionGroupOnXPath $cohesive_cond $group_id]
+    }
+}