diff --git a/kratos.gid/apps/Dam/write/write.tcl b/kratos.gid/apps/Dam/write/write.tcl
index ba287f229..a6b2f3169 100644
--- a/kratos.gid/apps/Dam/write/write.tcl
+++ b/kratos.gid/apps/Dam/write/write.tcl
@@ -69,29 +69,21 @@ proc ::Dam::write::UpdateMaterials { } {
         set constlaw [dict get $props ConstitutiveLaw]
         # Modificar la ley constitutiva
         set newconstlaw $constlaw
-        if {$constlaw eq "ElasticCohesive2DPlaneStress"} {
-            dict set matdict $mat THICKNESS  1.0000E+00
-            set newconstlaw "ElasticCohesive2DLaw"
-        }
         if {$constlaw eq "ElasticCohesive2DPlaneStrain"} {
-            dict set matdict $mat THICKNESS  1.0000E+00
+            dict set matdict $mat THICKNESS  1.0
             set newconstlaw "ElasticCohesive2DLaw"
         }
-        if {$constlaw eq "BilinearCohesive2DPlaneStress"}  {
-            dict set matdict $mat THICKNESS  1.0000E+00
-            set newconstlaw "BilinearCohesive2DLaw"
-        }
-        if {$constlaw eq "BilinearCohesive2DPlaneStrain"} {
-            dict set matdict $mat THICKNESS  1.0000E+00
-            set newconstlaw "BilinearCohesive2DLaw"
+        if {$constlaw eq "ElasticCohesive2DPlaneStress"} {
+            dict set matdict $mat THICKNESS  1.0
+            set newconstlaw "ElasticCohesive2DLaw"
         }
-        if {$constlaw eq "ExponentialCohesive2DPlaneStress"} {
-            dict set matdict $mat THICKNESS  1.0000E+00
-            set newconstlaw "ExponentialCohesive2DLaw"
+        if {$constlaw eq "IsotropicDamageCohesive2DPlaneStress"}  {
+            dict set matdict $mat THICKNESS  1.0
+            set newconstlaw "IsotropicDamageCohesive2DLaw"
         }
-        if {$constlaw eq "ExponentialCohesive2DPlaneStrain"} {
-            dict set matdict $mat THICKNESS  1.0000E+00
-            set newconstlaw "ExponentialCohesive2DLaw"
+        if {$constlaw eq "IsotropicDamageCohesive2DLawPlaneStrain"} {
+            dict set matdict $mat THICKNESS  1.0
+            set newconstlaw "IsotropicDamageCohesive2DLaw"
         }
 
         dict set matdict $mat CONSTITUTIVE_LAW_NAME $newconstlaw
diff --git a/kratos.gid/apps/Dam/xml/ConstitutiveLaws.xml b/kratos.gid/apps/Dam/xml/ConstitutiveLaws.xml
index e1490a40b..4401fb4eb 100644
--- a/kratos.gid/apps/Dam/xml/ConstitutiveLaws.xml
+++ b/kratos.gid/apps/Dam/xml/ConstitutiveLaws.xml
@@ -61,43 +61,27 @@
     <!--3D  Interface laws-->
     <CLaw n="ElasticCohesive3DLaw" pn="Joint Elastic Cohesive 3D"  Type="Interface3D" behaviour="Elastic"
     StrainSize="6"  AnalysisType="Linear,Non-Linear"
-	  help="Elastic law (no breakage)for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False"   DamageType="Local,NonLocal" >
+	  help="Elastic law (no breakage) for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False"   DamageType="Local,NonLocal" >
         <inputs>
-            <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
             <parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
-            <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e03"/>
-        </inputs>
-        <outputs>
-        </outputs>
-    </CLaw>
-
-    <CLaw n="BilinearCohesive3DLaw" pn="Joint Bilinear Cohesive 3D"  Type="Interface3D" behaviour="Elastic"
-    StrainSize="6"  AnalysisType="Linear,Non-Linear"
-	  help="Bilinear Cohesive law for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False"   DamageType="Local,NonLocal" >
-        <inputs>
-            <parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="2.6e10" />
-            <parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
-            <parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
-            <parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
-            <parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
-            <parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
-            <parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
+            <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
+            <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
         </inputs>
         <outputs>
-            <parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
         </outputs>
     </CLaw>
 
-    <CLaw n="ExponentialCohesive3DLaw" pn="Joint Exponential Cohesive 3D"  Type="Interface3D" behaviour="Elastic"
+    <CLaw n="IsotropicDamageCohesive3DLaw" pn="Joint Isotropic Damage Cohesive 3D"  Type="Interface3D" behaviour="Damage"
     StrainSize="6"  AnalysisType="Linear,Non-Linear"
-	  help="Exponential Cohesive law for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False"   DamageType="Local,NonLocal" >
+	  help="Isotropic Damage Cohesive law for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False" DamageType="Local,NonLocal" >
         <inputs>
-            <parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
-            <parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
-            <parameter n="FRACTURE_ENERGY" pn="Fracture Energy" v="90.0e06"/>
-            <parameter n="SHEAR_FRACTURE_ENERGY" pn="Shear Fracture Energy" v="350.0e06"/>
-            <parameter n="STATE_VARIABLE" pn="State Variable" v="1.0e-15"/>
-            <parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
+            <parameter n="DAMAGE_EVOLUTION_LAW" pn="Damage Evolution Law" v="1"/>
+            <parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
+            <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
+            <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
+            <parameter n="TENSILE_STRENGTH" pn="Tensile Strength" unit_magnitude="P" units="Pa" v="1e03"/>
+            <parameter n="FRACTURE_ENERGY" pn="Fracture Energy" unit_magnitude="Energy/L^2" units="J/m^2" v="1e02"/>
+            <parameter n="BETA_EQSTRAIN_SHEAR_FACTOR" pn="Beta Factor" v="1.0"/>
         </inputs>
         <outputs>
             <parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
@@ -218,9 +202,9 @@
     <CLaw n="ElasticCohesive2DPlaneStrain" pn="Simplified Joint Plane Strain"  Type="InterfaceStrain" behaviour="Elastic" StrainSize="3"   DamageType="Local,NonLocal"
           help="Simplified law for joints in Plane Strain problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False"   AnalysisType="Linear,Non-Linear">
         <inputs>
-            <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
             <parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
-            <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e03"/>
+            <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
+            <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
         </inputs>
         <outputs>
         </outputs>
@@ -229,74 +213,43 @@
     <CLaw n="ElasticCohesive2DPlaneStress" pn="Simplified Joint Plane Stress"  Type="InterfaceStress" behaviour="Elastic" StrainSize="3"  DamageType="Local,NonLocal"
           help="Simplified law for joints in Plane Stress problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False"   AnalysisType="Linear,Non-Linear">
         <inputs>
-            <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
             <parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
-            <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e03"/>
-        </inputs>
-        <outputs>
-        </outputs>
-    </CLaw>
-
-    <CLaw n="BilinearCohesive2DPlaneStrain" pn="Joint Bilinear Cohesive Plane Strain"  Type="InterfaceStrain" behaviour="Elastic" StrainSize="3"   DamageType="Local,NonLocal"
-          help="Bilinear Cohesive law for joints in Plane Strain problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False"   AnalysisType="Linear,Non-Linear">
-        <inputs>
-            <parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="2.6e10" />
-            <parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
-            <parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
-            <parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
-            <parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
-            <parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
-            <parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
-        </inputs>
-        <outputs>
-            <parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
-        </outputs>
-    </CLaw>
-
-    <CLaw n="BilinearCohesive2DPlaneStress" pn="Joint Bilinear Cohesive Plane Stress"  Type="InterfaceStress" behaviour="Elastic" StrainSize="3"  DamageType="Local,NonLocal"
-          help="Bilinear Cohesive law for joints in Plane Stress problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False"   AnalysisType="Linear,Non-Linear">
-        <inputs>
-            <parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="2.6e10" />
-            <parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
-            <parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
-            <parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
-            <parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
-            <parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
-            <parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
+            <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
+            <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
         </inputs>
         <outputs>
-            <parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
         </outputs>
     </CLaw>
 
-    <CLaw n="ExponentialCohesive2DPlaneStrain" pn="Joint Exponential Cohesive Plane Strain"  Type="InterfaceStrain" behaviour="Elastic" StrainSize="3"   DamageType="Local,NonLocal"
-          help="Exponential Cohesive law for joints in Plane Strain problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False"   AnalysisType="Linear,Non-Linear">
+    <CLaw n="IsotropicDamageCohesive2DPlaneStrain" pn="Isotropic Damage Cohesive Plane Strain"  Type="InterfaceStrain" behaviour="Elastic" StrainSize="3"   DamageType="Local,NonLocal"
+          help="Isotropic Damage Cohesive law for joints in Plane Strain problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False"   AnalysisType="Linear,Non-Linear">
         <inputs>
-            <parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
-            <parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
-            <parameter n="FRACTURE_ENERGY" pn="Fracture Energy" v="90.0e06"/>
-            <parameter n="SHEAR_FRACTURE_ENERGY" pn="Shear Fracture Energy" v="350.0e06"/>
-            <parameter n="STATE_VARIABLE" pn="State Variable" v="1.0e-15"/>
-            <parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
+            <parameter n="DAMAGE_EVOLUTION_LAW" pn="Damage Evolution Law" v="1"/>
+            <parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
+            <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
+            <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
+            <parameter n="TENSILE_STRENGTH" pn="Tensile Strength" unit_magnitude="P" units="Pa" v="1e03"/>
+            <parameter n="FRACTURE_ENERGY" pn="Fracture Energy" unit_magnitude="Energy/L^2" units="J/m^2" v="1e02"/>
+            <parameter n="BETA_EQSTRAIN_SHEAR_FACTOR" pn="Beta Factor" v="1.0"/>
         </inputs>
         <outputs>
             <parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
         </outputs>
     </CLaw>
 
-    <CLaw n="ExponentialCohesive2DPlaneStress" pn="Joint Exponential Cohesive Plane Stress"  Type="InterfaceStress" behaviour="Elastic" StrainSize="3"  DamageType="Local,NonLocal"
-          help="Exponential Cohesive law for joints in Plane Stress problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False"   AnalysisType="Linear,Non-Linear">
+    <CLaw n="IsotropicDamageCohesive2DPlaneStress" pn="Isotropic Damage Cohesive Plane Stress"  Type="InterfaceStress" behaviour="Elastic" StrainSize="3"  DamageType="Local,NonLocal"
+          help="Isotropic Damage Cohesive law for joints in Plane Stress problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False"   AnalysisType="Linear,Non-Linear">
         <inputs>
-            <parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
-            <parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
-            <parameter n="FRACTURE_ENERGY" pn="Fracture Energy" v="90.0e06"/>
-            <parameter n="SHEAR_FRACTURE_ENERGY" pn="Shear Fracture Energy" v="350.0e06"/>
-            <parameter n="STATE_VARIABLE" pn="State Variable" v="1.0e-15"/>
-            <parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
+            <parameter n="DAMAGE_EVOLUTION_LAW" pn="Damage Evolution Law" v="1"/>
+            <parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
+            <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
+            <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
+            <parameter n="TENSILE_STRENGTH" pn="Tensile Strength" unit_magnitude="P" units="Pa" v="1e03"/>
+            <parameter n="FRACTURE_ENERGY" pn="Fracture Energy" unit_magnitude="Energy/L^2" units="J/m^2" v="1e02"/>
+            <parameter n="BETA_EQSTRAIN_SHEAR_FACTOR" pn="Beta Factor" v="1.0"/>
         </inputs>
         <outputs>
             <parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
         </outputs>
     </CLaw>
-
 </ConstitutiveLaws>
diff --git a/kratos.gid/apps/Dam/xml/Materials.xml b/kratos.gid/apps/Dam/xml/Materials.xml
index 73df0d4ba..ab1d2befe 100644
--- a/kratos.gid/apps/Dam/xml/Materials.xml
+++ b/kratos.gid/apps/Dam/xml/Materials.xml
@@ -19,38 +19,23 @@
     </inputs>
   </Material>
 
-    <Material n="Simplified-Joint" MaterialType="Dam" help="">
+    <Material n="Elastic Cohesive Joint" MaterialType="Dam" help="">
     <inputs>
-        <parameter n="DENSITY" pn="Density" v="2400"/>
-        <parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="2.6e10"/>
-        <parameter n="YIELD_STRESS" pn="Yield Stress"  unit_magnitude="P" units="Pa" v="20e06"/>
-        <parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
+        <parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
+        <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
+        <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
     </inputs>
   </Material>
 
-  <Material n="Bilinear Strong-Joint" MaterialType="Dam" help="">
+  <Material n="Isotropic Damage Cohesive Joint" MaterialType="Dam" help="">
     <inputs>
-        <parameter n="DENSITY" pn="Density" v="2400"/>
-        <parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="35.0e9" />
-        <parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
-        <parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
-        <parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
-        <parameter n="YIELD_STRESS" pn="Yield Stress" v="3.5e8"/>
-        <parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
-        <parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
-    </inputs>
-  </Material>
-
-    <Material n="Bilinear Weak-Joint" MaterialType="Dam" help="">
-    <inputs>
-        <parameter n="DENSITY" pn="Density" v="2400"/>
-        <parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="1.0" />
-        <parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
-        <parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
-        <parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
-        <parameter n="YIELD_STRESS" pn="Yield Stress" v="1.0"/>
-        <parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
-        <parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
+        <parameter n="DAMAGE_EVOLUTION_LAW" pn="Damge Evolution Law" v="1"/>
+        <parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
+        <parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
+        <parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
+        <parameter n="TENSILE_STRENGTH" pn="Tensile Strength" unit_magnitude="P" units="Pa" v="1e03"/>
+        <parameter n="FRACTURE_ENERGY" pn="Fracture Energy" unit_magnitude="Energy/L^2" units="J/m^2" v="1e02"/>
+        <parameter n="BETA_EQSTRAIN_SHEAR_FACTOR" pn="Beta Equivalent Strain Shear Factor" v="1.0"/>
     </inputs>
   </Material>