Skip to content

Commit d1eb391

Browse files
farrufat-cimnefarrufat-cimne
authored
Merge pull request #664 from KratosMultiphysics/dem/initial_cond
[DEM]add initial conditions 2d/3d
2 parents 222e809 + 7c8552b commit d1eb391

11 files changed

Lines changed: 226 additions & 40 deletions

‎kratos.gid/apps/CDEM/write/writeMDPA_Parts.tcl‎

Lines changed: 152 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -31,7 +31,8 @@ proc DEM::write::WriteMDPAParts { } {
3131

3232
# SubmodelParts
3333
write::writePartSubModelPart
34-
writeVelocityMeshes
34+
if {$::Model::SpatialDimension eq "2D"} { writeVelocityMeshes2D
35+
} else {writeVelocityMeshes}
3536

3637
# CustomSubmodelParts
3738
#WriteWallCustomDEMSmp not required for dem.
@@ -292,18 +293,135 @@ proc DEM::write::writeVelocityMeshes { } {
292293
}
293294

294295

296+
proc DEM::write::writeVelocityMeshes2D { } {
297+
set i 0
298+
foreach {cond group_list} [GetSpheresGroupsListInConditions] {
299+
if {$cond eq "DEM-VelocityBC2D"} {
300+
set cnd [Model::getCondition $cond]
301+
foreach group $group_list {
302+
incr i
303+
write::WriteString "Begin SubModelPart $i // GUI DEM-VelocityBC - $cond - group identifier: $group"
304+
write::WriteString " Begin SubModelPartData // DEM-VelocityBC. Group name: $group"
305+
set xp1 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = '$cond'\]/group\[@n = '$group'\]"
306+
set group_node [[customlib::GetBaseRoot] selectNodes $xp1]
307+
308+
set prescribeMotion_flag [write::getValueByNode [$group_node selectNodes "./value\[@n='PrescribeMotion_flag'\]"]]
309+
if {[write::isBooleanTrue $prescribeMotion_flag]} {
310+
311+
# Linear velocity
312+
set velocity [write::getValueByNode [$group_node selectNodes "./value\[@n='VelocityModulus'\]"]]
313+
lassign [write::getValueByNode [$group_node selectNodes "./value\[@n='DirectionVector'\]"]] velocity_X velocity_Y
314+
lassign [MathUtils::VectorNormalized [list $velocity_X $velocity_Y]] velocity_X velocity_Y
315+
lassign [MathUtils::ScalarByVectorProd $velocity [list $velocity_X $velocity_Y] ] vx vy
316+
write::WriteString " LINEAR_VELOCITY \[3\] ($vx, $vy, 0.0)"
317+
318+
# Period
319+
set periodic [write::getValueByNode [$group_node selectNodes "./value\[@n='LinearPeriodic'\]"]]
320+
if {[write::isBooleanTrue $periodic]} {
321+
set period [write::getValueByNode [$group_node selectNodes "./value\[@n='LinearPeriod'\]"]]
322+
} else {
323+
set period 0.0
324+
}
325+
write::WriteString " VELOCITY_PERIOD $period"
326+
327+
# Angular velocity
328+
set avelocity [write::getValueByNode [$group_node selectNodes "./value\[@n='AngularVelocityModulus'\]"]]
329+
write::WriteString " ANGULAR_VELOCITY \[3\] (0.0,0.0,$avelocity)"
330+
331+
# Angular center of rotation
332+
lassign [write::getValueByNode [$group_node selectNodes "./value\[@n='CenterOfRotation'\]"]] oX oY
333+
write::WriteString " ROTATION_CENTER \[3\] ($oX,$oY,0.0)"
334+
335+
# Angular Period
336+
set angular_periodic [write::getValueByNode [$group_node selectNodes "./value\[@n='AngularPeriodic'\]"]]
337+
if {[write::isBooleanTrue $angular_periodic]} {
338+
set angular_period [write::getValueByNode [$group_node selectNodes "./value\[@n='AngularPeriod'\]"]]
339+
} else {
340+
set angular_period 0.0
341+
}
342+
write::WriteString " ANGULAR_VELOCITY_PERIOD $angular_period"
343+
344+
# set intervals
345+
set LinearStartTime [write::getValueByNode [$group_node selectNodes "./value\[@n='LinearStartTime'\]"]]
346+
set LinearEndTime [write::getValueByNode [$group_node selectNodes "./value\[@n='LinearEndTime'\]"]]
347+
set AngularStartTime [write::getValueByNode [$group_node selectNodes "./value\[@n='AngularStartTime'\]"]]
348+
set AngularEndTime [write::getValueByNode [$group_node selectNodes "./value\[@n='AngularEndTime'\]"]]
349+
write::WriteString " VELOCITY_START_TIME $LinearStartTime"
350+
write::WriteString " VELOCITY_STOP_TIME $LinearEndTime"
351+
write::WriteString " ANGULAR_VELOCITY_START_TIME $AngularStartTime"
352+
write::WriteString " ANGULAR_VELOCITY_STOP_TIME $AngularEndTime"
353+
354+
# Interval. # TODO, able to assign intervals as start and end time. Both linear and angular could be different
355+
# set interval [write::getValueByNode [$group_node selectNodes "./value\[@n='Interval'\]"]]
356+
# lassign [write::getInterval $interval] ini end
357+
# if {![string is double $ini]} {
358+
# set ini [write::getValue DEMTimeParameters StartTime]
359+
# }
360+
# write::WriteString " ${cond}_START_TIME $ini"
361+
# if {![string is double $end]} {
362+
# set end [write::getValue DEMTimeParameters EndTime]
363+
# }
364+
# write::WriteString " ${cond}_STOP_TIME $end"
365+
}
366+
367+
write::WriteString " End SubModelPartData"
368+
write::WriteString " Begin SubModelPartNodes"
369+
GiD_WriteCalculationFile nodes -sorted [dict create [write::GetWriteGroupName $group] [subst "%10i\n"]]
370+
write::WriteString " End SubModelPartNodes"
371+
write::WriteString "End SubModelPart"
372+
write::WriteString ""
373+
}
374+
} elseif {$cond eq "DEM-VelocityIC2D"} {
375+
set cnd [Model::getCondition $cond]
376+
foreach group $group_list {
377+
incr i
378+
write::WriteString "Begin SubModelPart $i // GUI DEM-VelocityIC - $cond - group identifier: $group"
379+
write::WriteString " Begin SubModelPartData // DEM-VelocityIC. Group name: $group"
380+
set xp1 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = '$cond'\]/group\[@n = '$group'\]"
381+
set group_node [[customlib::GetBaseRoot] selectNodes $xp1]
382+
383+
set prescribeMotion_flag [write::getValueByNode [$group_node selectNodes "./value\[@n='PrescribeMotion_flag'\]"]]
384+
if {[write::isBooleanTrue $prescribeMotion_flag]} {
385+
386+
# Linear velocity
387+
set velocity [write::getValueByNode [$group_node selectNodes "./value\[@n='InitialVelocityModulus'\]"]]
388+
lassign [write::getValueByNode [$group_node selectNodes "./value\[@n='iDirectionVector'\]"]] velocity_X velocity_Y
389+
lassign [MathUtils::VectorNormalized [list $velocity_X $velocity_Y]] velocity_X velocity_Y
390+
lassign [MathUtils::ScalarByVectorProd $velocity [list $velocity_X $velocity_Y] ] vx vy
391+
write::WriteString " INITIAL_VELOCITY_X_VALUE $vx"
392+
write::WriteString " INITIAL_VELOCITY_Y_VALUE $vy"
393+
write::WriteString " INITIAL_VELOCITY_Z_VALUE 0.0"
394+
395+
# Angular velocity
396+
set avelocity [write::getValueByNode [$group_node selectNodes "./value\[@n='InitialAngularVelocityModulus'\]"]]
397+
write::WriteString " INITIAL_ANGULAR_VELOCITY_X_VALUE 0.0"
398+
write::WriteString " INITIAL_ANGULAR_VELOCITY_Y_VALUE 0.0"
399+
write::WriteString " INITIAL_ANGULAR_VELOCITY_Z_VALUE $avelocity"
400+
}
401+
402+
write::WriteString " End SubModelPartData"
403+
write::WriteString " Begin SubModelPartNodes"
404+
GiD_WriteCalculationFile nodes -sorted [dict create [write::GetWriteGroupName $group] [subst "%10i\n"]]
405+
write::WriteString " End SubModelPartNodes"
406+
write::WriteString "End SubModelPart"
407+
write::WriteString ""
408+
}
409+
}
410+
}
411+
}
412+
413+
295414
proc DEM::write::GetSpheresGroupsListInConditions { } {
296415
set conds_groups_dict [dict create ]
297416
set groups [list ]
298-
299417
# Get all the groups with spheres
300418
foreach group [GetSpheresGroups] {
301-
foreach surface [GiD_EntitiesGroups get $group elements] {
302-
foreach involved_group [GiD_EntitiesGroups entity_groups elements $surface] {
303-
set involved_group_id [write::GetWriteGroupName $involved_group]
304-
if {$involved_group_id ni $groups} {lappend groups $involved_group_id}
305-
}
306-
}
419+
foreach surface [GiD_EntitiesGroups get $group elements] {
420+
foreach involved_group [GiD_EntitiesGroups entity_groups elements $surface] {
421+
set involved_group_id [write::GetWriteGroupName $involved_group]
422+
if {$involved_group_id ni $groups} {lappend groups $involved_group_id}
423+
}
424+
}
307425
}
308426

309427
# Find the relations condition -> group
@@ -321,12 +439,14 @@ proc DEM::write::GetSpheresGroupsListInConditions { } {
321439

322440
proc DEM::write::GetSpheresGroups { } {
323441
set groups [list ]
324-
set xp1 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = 'DEM-VelocityBC'\]/group"
442+
if {$::Model::SpatialDimension eq "2D"} { set xp1 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = 'DEM-VelocityBC2D'\]/group"
443+
} else {set xp1 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = 'DEM-VelocityBC'\]/group"}
325444
foreach group [[customlib::GetBaseRoot] selectNodes $xp1] {
326445
set groupid [$group @n]
327446
lappend groups [write::GetWriteGroupName $groupid]
328447
}
329-
set xp2 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = 'DEM-VelocityIC'\]/group"
448+
if {$::Model::SpatialDimension eq "2D"} { set xp2 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = 'DEM-VelocityIC2D'\]/group"
449+
} else {set xp2 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = 'DEM-VelocityIC'\]/group"}
330450
foreach group [[customlib::GetBaseRoot] selectNodes $xp2] {
331451
set groupid [$group @n]
332452
lappend groups [write::GetWriteGroupName $groupid]
@@ -339,13 +459,33 @@ proc DEM::write::writeMaterialsParts { } {
339459
variable partsProperties
340460
set xp1 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = 'PartsCont'\]/group"
341461
set partsProperties $::write::mat_dict
342-
set printable [list PARTICLE_DENSITY YOUNG_MODULUS POISSON_RATIO FRICTION COEFFICIENT_OF_RESTITUTION PARTICLE_MATERIAL ROLLING_FRICTION ROLLING_FRICTION_WITH_WALLS CONTACT_SIGMA_MIN CONTACT_TAU_ZERO CONTACT_INTERNAL_FRICC ConstitutiveLaw]
343-
foreach group [dict keys $partsProperties] {
462+
set printable [list PARTICLE_DENSITY \
463+
YOUNG_MODULUS \
464+
POISSON_RATIO \
465+
FRICTION \
466+
COEFFICIENT_OF_RESTITUTION \
467+
PARTICLE_MATERIAL \
468+
ROLLING_FRICTION \
469+
ROLLING_FRICTION_WITH_WALLS \
470+
CONTACT_SIGMA_MIN \
471+
CONTACT_TAU_ZERO \
472+
CONTACT_INTERNAL_FRICC \
473+
ConstitutiveLaw \
474+
SHEAR_ENERGY_COEF \
475+
LOOSE_MATERIAL_YOUNG_MODULUS \
476+
FRACTURE_ENERGY \
477+
INTERNAL_FRICTION_ANGLE]
478+
479+
foreach group [dict keys $partsProperties] {
344480
write::WriteString "Begin Properties [dict get $partsProperties $group MID]"
345481
foreach {prop val} [dict get $partsProperties $group] {
346482
if {$prop in $printable} {
347483
if {$prop eq "ConstitutiveLaw"} {
348484
write::WriteString " DEM_CONTINUUM_CONSTITUTIVE_LAW_NAME $val"
485+
} elseif {$prop eq "FRICTION"} {
486+
set pi $MathUtils::PI
487+
set propvalue [expr {tan($val*$pi/180.0)}]
488+
write::WriteString " FRICTION $propvalue"
349489
} else {
350490
write::WriteString " $prop $val"
351491
}

‎kratos.gid/apps/CDEM/write/writeMDPA_Walls.tcl‎

Lines changed: 6 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -45,7 +45,12 @@ proc DEM::write::WriteWallProperties { } {
4545
# write::WriteString " $prop $v"
4646
# }
4747
#}
48-
write::WriteString " FRICTION [write::getValueByNode [$group selectNodes "./value\[@n='friction_coeff'\]"]]"
48+
49+
set friction_value [write::getValueByNode [$group selectNodes "./value\[@n='friction_angle'\]"]]
50+
set pi $MathUtils::PI
51+
set propvalue [expr {tan($friction_value*$pi/180.0)}]
52+
write::WriteString " FRICTION $propvalue"
53+
#write::WriteString " FRICTION [write::getValueByNode [$group selectNodes "./value\[@n='friction_coeff'\]"]]"
4954
write::WriteString " WALL_COHESION [write::getValueByNode [$group selectNodes "./value\[@n='WallCohesion'\]"]]"
5055
set compute_wear_bool [write::getValueByNode [$group selectNodes "./value\[@n='DEM_Wear'\]"]]
5156
if {[write::isBooleanTrue $compute_wear_bool]} {
@@ -430,8 +435,6 @@ proc DEM::write::writeConditionMeshes { } {
430435
write::WriteString " IDENTIFIER [write::transformGroupName $group]"
431436

432437
set material_analysis [write::getValue DEMTestMaterial Active]
433-
W $material_analysis
434-
W "material_analysis"
435438
if {$material_analysis == "true"} {
436439
set is_material_test [write::getValueByNode [$group_node selectNodes "./value\[@n='MaterialTest'\]"]]
437440
if {$is_material_test == "true"} {

‎kratos.gid/apps/CDEM/write/writeProjectParameters.tcl‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -89,7 +89,7 @@ proc CDEM::write::getParametersDict { } {
8989
dict set project_parameters_dict "FinalTime" $FinalTime
9090
# TODO: check for inconsistencies in DEMTIMEPARAMETERS UN
9191
# dict set project_parameters_dict "ControlTime" [write::getValue DEMTimeParameters ScreenInfoOutput]
92-
# dict set project_parameters_dict "NeighbourSearchFrequency" [write::getValue DEMTimeParameters NeighbourSearchFrequency]
92+
dict set project_parameters_dict "NeighbourSearchFrequency" [write::getValue DEMTimeParameters NeighbourSearchFrequency]
9393
dict set project_parameters_dict "GraphExportFreq" [write::getValue DGraphs GraphExportFreq]
9494
dict set project_parameters_dict "VelTrapGraphExportFreq" 1e-3
9595

‎kratos.gid/apps/CDEM/xml/ConstitutiveLawsC.xml‎

Lines changed: 26 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -5,7 +5,7 @@
55
<inputs>
66
<parameter n="CONTACT_TAU_ZERO" pn="Tangential Strength" units="MPa" unit_magnitude="P" v="25"/>
77
<parameter n="CONTACT_SIGMA_MIN" pn="Nornal Tensile Strength" units="MPa" unit_magnitude="P" v="5"/>
8-
<parameter n="CONTACT_INTERNAL_FRICC" pn="Internal Friction Coeff" v="1"/>
8+
<parameter n="CONTACT_INTERNAL_FRICC" pn="Contact Internal Friction" v="1"/>
99
</inputs>
1010
<outputs></outputs>
1111
</CLaw>
@@ -14,7 +14,7 @@
1414
<inputs>
1515
<parameter n="CONTACT_TAU_ZERO" pn="Tangential Strength" units="MPa" unit_magnitude="P" v="25"/>
1616
<parameter n="CONTACT_SIGMA_MIN" pn="Nornal Tensile Strength" units="MPa" unit_magnitude="P" v="5"/>
17-
<parameter n="CONTACT_INTERNAL_FRICC" pn="Internal Friction Coeff" v="1"/>
17+
<parameter n="CONTACT_INTERNAL_FRICC" pn="Contact Internal Friction" v="1"/>
1818
</inputs>
1919
<outputs></outputs>
2020
</CLaw>
@@ -49,6 +49,30 @@
4949
<outputs></outputs>
5050
</CLaw>
5151

52+
<CLaw n="DEM_KDEM_with_damage_parallel_bond" pn="KDEM with damage" ProductionReady="ProductionReady" help="" ImplementedInApplication="DEMApplication" App="CDEM" Dimension="3D" OutputMode="Materials">
53+
<inputs>
54+
<parameter n="CONTACT_TAU_ZERO" pn="Mohr-Coulomb Cohesion" units="MPa" unit_magnitude="P" v="10"/>
55+
<parameter n="CONTACT_SIGMA_MIN" pn="Normal Tensile Strength" units="MPa" unit_magnitude="P" v="13"/>
56+
<parameter n="CONTACT_INTERNAL_FRICC" pn="Contact Internal Friction" v="0.44"/>
57+
<parameter n="SHEAR_ENERGY_COEF" pn="Shear energy" v="4.0"/>
58+
<parameter n="LOOSE_MATERIAL_YOUNG_MODULUS" pn="Loose Young modulus" units="Pa" unit_magnitude="P" v="14e8"/>
59+
<parameter n="FRACTURE_ENERGY" pn="Fracture energy" v="20.0"/>
60+
</inputs>
61+
<outputs></outputs>
62+
</CLaw>
63+
64+
<CLaw n="DEM_KDEM_with_damage_parallel_bond_2D" pn="KDEM with damage" ProductionReady="ProductionReady" help="Not implemented" ImplementedInApplication="DEMApplication" App="CDEM" Dimension="2D" OutputMode="Materials">
65+
<inputs>
66+
<parameter n="CONTACT_TAU_ZERO" pn="Mohr-Coulomb Cohesion" units="MPa" unit_magnitude="P" v="10"/>
67+
<parameter n="CONTACT_SIGMA_MIN" pn="Normal Tensile Strength" units="MPa" unit_magnitude="P" v="13"/>
68+
<parameter n="CONTACT_INTERNAL_FRICC" pn="Contact Internal Friction" v="0.44"/>
69+
<parameter n="SHEAR_ENERGY_COEF" pn="Shear energy" v="4.0"/>
70+
<parameter n="LOOSE_MATERIAL_YOUNG_MODULUS" pn="Loose Young modulus" units="Pa" unit_magnitude="P" v="14e8"/>
71+
<parameter n="FRACTURE_ENERGY" pn="Fracture energy" v="20.0"/>
72+
</inputs>
73+
<outputs></outputs>
74+
</CLaw>
75+
5276
<!-- <CLaw n="LinearElasticPlaneStrain2DLaw" pn="Linear Elastic Plane Strain" ProductionReady="ProductionReady" Type="PlaneStrain" Behaviour="Elastic" StrainSize="3" ImplementedInApplication="StructuralMechanicsApplication" help="Linear elastic behaviour in 2D plane strain" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False">
5377
<inputs>
5478
<parameter n="DENSITY" pn="Density" unit_magnitude="Density" units="kg/m^3" v="7850"/>

‎kratos.gid/apps/CDEM/xml/MaterialsC.xml‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -5,8 +5,8 @@
55
<parameter n="PARTICLE_DENSITY" pn="Density" unit_magnitude="Density" units="kg/m^3" v="2500.0" type="double" help="Density of the material (not bulk density)"/>
66
<parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="1.0e5" type="double" help="The Youngs Modulus of the material (not bulk stiffness)"/>
77
<parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20" type="double" help="The Poissons Ratio of the material (not bulk Ratio)"/>
8-
<parameter n="FRICTION" pn="Coefficient of friction" v="1.0" type="double" help="Particle friction coefficient (for Coulombs friction)"/>
9-
<parameter n="COEFFICIENT_OF_RESTITUTION" pn="Coefficient of restitution" v="0.2" type="double" help="Particle friction angle (for Coulombs friction)"/>
8+
<parameter n="FRICTION" pn="Friction angle" v="30.0" type="double" units="deg" unit_magnitude="Angle" help="Particle friction angle (for Coulombs friction)"/>
9+
<parameter n="COEFFICIENT_OF_RESTITUTION" pn="Coefficient of restitution" v="0.2" type="double" help=""/>
1010
<parameter n="PARTICLE_MATERIAL" pn="Color" v="1" type="integer"/>
1111
<parameter n="ROLLING_FRICTION" pn="Rolling friction" v="0.01" type="double" help="Rolling friction"/>
1212
<parameter n="ROLLING_FRICTION_WITH_WALLS" pn="Rolling friction With Walls" v="0.01" type="double" help="Rolling friction with walls"/>

‎kratos.gid/apps/DEM/write/writeMDPA_Parts.tcl‎

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -365,6 +365,10 @@ proc DEM::write::writeMaterialsParts { } {
365365
if {$prop in $printable} {
366366
if {$prop eq "ConstitutiveLaw"} {
367367
write::WriteString " DEM_DISCONTINUUM_CONSTITUTIVE_LAW_NAME $val"
368+
} elseif {$prop eq "FRICTION"} {
369+
set pi $MathUtils::PI
370+
set propvalue [expr {tan($val*$pi/180.0)}]
371+
write::WriteString " FRICTION $propvalue"
368372
} else {
369373
write::WriteString " $prop $val"
370374
}

‎kratos.gid/apps/DEM/write/writeMDPA_Walls.tcl‎

Lines changed: 5 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -26,7 +26,6 @@ proc DEM::write::WriteMDPAWalls { } {
2626

2727

2828
proc DEM::write::WriteWallProperties { } {
29-
#set print_list [list "FRICTION" "WALL_COHESION" "COMPUTE_WEAR" "SEVERITY_OF_WEAR" "IMPACT_WEAR_SEVERITY" "BRINELL_HARDNESS" "YOUNG_MODULUS" "POISSON_RATIO"]
3029
set wall_properties [dict create ]
3130
set cnd [Model::getCondition "DEM-FEM-Wall"]
3231

@@ -45,7 +44,11 @@ proc DEM::write::WriteWallProperties { } {
4544
# write::WriteString " $prop $v"
4645
# }
4746
#}
48-
write::WriteString " FRICTION [write::getValueByNode [$group selectNodes "./value\[@n='friction_coeff'\]"]]"
47+
set friction_value [write::getValueByNode [$group selectNodes "./value\[@n='friction_angle'\]"]]
48+
set pi $MathUtils::PI
49+
set propvalue [expr {tan($friction_value*$pi/180.0)}]
50+
write::WriteString " FRICTION $propvalue"
51+
# write::WriteString " FRICTION [write::getValueByNode [$group selectNodes "./value\[@n='friction_coeff'\]"]]"
4952
write::WriteString " WALL_COHESION [write::getValueByNode [$group selectNodes "./value\[@n='WallCohesion'\]"]]"
5053
set compute_wear_bool [write::getValueByNode [$group selectNodes "./value\[@n='DEM_Wear'\]"]]
5154
if {[write::isBooleanTrue $compute_wear_bool]} {

‎kratos.gid/apps/DEM/write/writeProjectParameters.tcl‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -89,7 +89,7 @@ proc DEM::write::getParametersDict { } {
8989
dict set project_parameters_dict "FinalTime" $FinalTime
9090
# TODO: check for inconsistencies in DEMTIMEPARAMETERS UN
9191
# dict set project_parameters_dict "ControlTime" [write::getValue DEMTimeParameters ScreenInfoOutput]
92-
# dict set project_parameters_dict "NeighbourSearchFrequency" [write::getValue DEMTimeParameters NeighbourSearchFrequency]
92+
dict set project_parameters_dict "NeighbourSearchFrequency" [write::getValue DEMTimeParameters NeighbourSearchFrequency]
9393
dict set project_parameters_dict "GraphExportFreq" [write::getValue DGraphs GraphExportFreq]
9494
dict set project_parameters_dict "VelTrapGraphExportFreq" 1e-3
9595

0 commit comments

Comments
 (0)