update, new claw

Author: farrufat-cimne

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

@@ -459,8 +459,24 @@ proc DEM::write::writeMaterialsParts { } {
     variable partsProperties
     set xp1 "[spdAux::getRoute [GetAttribute conditions_un]]/condition\[@n = 'PartsCont'\]/group"
     set partsProperties $::write::mat_dict
-    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]
-    foreach group [dict keys $partsProperties] {
+    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 \
+						SHEAR_ENERGY_COEF \
+						LOOSE_MATERIAL_YOUNG_MODULUS \
+						FRACTURE_ENERGY \
+						INTERNAL_FRICTION_ANGLE]
+
+	foreach group [dict keys $partsProperties] {
         write::WriteString "Begin Properties [dict get $partsProperties $group MID]"
         foreach {prop val} [dict get $partsProperties $group] {
             if {$prop in $printable} {

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

@@ -5,7 +5,7 @@
         <inputs>
             <parameter n="CONTACT_TAU_ZERO" pn="Tangential Strength" units="MPa" unit_magnitude="P" v="25"/>
             <parameter n="CONTACT_SIGMA_MIN" pn="Nornal Tensile Strength" units="MPa" unit_magnitude="P" v="5"/>
-            <parameter n="CONTACT_INTERNAL_FRICC" pn="Internal Friction Coeff" v="1"/>
+            <parameter n="CONTACT_INTERNAL_FRICC" pn="Contact Internal Friction" v="1"/>
         </inputs>
         <outputs></outputs>
     </CLaw>
@@ -14,7 +14,7 @@
         <inputs>
             <parameter n="CONTACT_TAU_ZERO" pn="Tangential Strength" units="MPa" unit_magnitude="P" v="25"/>
             <parameter n="CONTACT_SIGMA_MIN" pn="Nornal Tensile Strength" units="MPa" unit_magnitude="P" v="5"/>
-            <parameter n="CONTACT_INTERNAL_FRICC" pn="Internal Friction Coeff" v="1"/>
+            <parameter n="CONTACT_INTERNAL_FRICC" pn="Contact Internal Friction" v="1"/>
         </inputs>
         <outputs></outputs>
     </CLaw>
@@ -49,6 +49,30 @@
         <outputs></outputs>
     </CLaw>
 
+    <CLaw n="DEM_KDEM_with_damage_parallel_bond" pn="KDEM with damage" ProductionReady="ProductionReady" help="" ImplementedInApplication="DEMApplication" App="CDEM" Dimension="3D" OutputMode="Materials">
+        <inputs>
+            <parameter n="CONTACT_TAU_ZERO" pn="Mohr-Coulomb Cohesion" units="MPa" unit_magnitude="P" v="10"/>
+            <parameter n="CONTACT_SIGMA_MIN" pn="Normal Tensile Strength" units="MPa" unit_magnitude="P" v="13"/>
+            <parameter n="CONTACT_INTERNAL_FRICC" pn="Contact Internal Friction" v="0.44"/>
+            <parameter n="SHEAR_ENERGY_COEF" pn="Shear energy" v="4.0"/>
+            <parameter n="LOOSE_MATERIAL_YOUNG_MODULUS" pn="Loose Young modulus" units="Pa" unit_magnitude="P" v="14e8"/>
+            <parameter n="FRACTURE_ENERGY" pn="Fracture energy" v="20.0"/>
+        </inputs>
+        <outputs></outputs>
+    </CLaw>
+
+    <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">
+        <inputs>
+            <parameter n="CONTACT_TAU_ZERO" pn="Mohr-Coulomb Cohesion" units="MPa" unit_magnitude="P" v="10"/>
+            <parameter n="CONTACT_SIGMA_MIN" pn="Normal Tensile Strength" units="MPa" unit_magnitude="P" v="13"/>
+            <parameter n="CONTACT_INTERNAL_FRICC" pn="Contact Internal Friction" v="0.44"/>
+            <parameter n="SHEAR_ENERGY_COEF" pn="Shear energy" v="4.0"/>
+            <parameter n="LOOSE_MATERIAL_YOUNG_MODULUS" pn="Loose Young modulus" units="Pa" unit_magnitude="P" v="14e8"/>
+            <parameter n="FRACTURE_ENERGY" pn="Fracture energy" v="20.0"/>
+        </inputs>
+        <outputs></outputs>
+    </CLaw>
+
     <!-- <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">
         <inputs>
             <parameter n="DENSITY" pn="Density" unit_magnitude="Density" units="kg/m^3" v="7850"/>