Skip to content

Commit

Permalink
Add new joints to Dam app and remove old ones (#982)
Browse files Browse the repository at this point in the history
  • Loading branch information
@joaquinirazabal
joaquinirazabal authored Jul 3, 2024
1 parent af19628 commit 3edbdd7
Show file tree
Hide file tree
Showing 3 changed files with 56 additions and 126 deletions.
28 changes: 10 additions & 18 deletions kratos.gid/apps/Dam/write/write.tcl
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
115 changes: 34 additions & 81 deletions kratos.gid/apps/Dam/xml/ConstitutiveLaws.xml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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" />
Expand Down Expand Up @@ -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>
Expand All @@ -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>
39 changes: 12 additions & 27 deletions kratos.gid/apps/Dam/xml/Materials.xml
View file
Original file line number Diff line number Diff line change
Expand Up @@ -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>

Expand Down

0 comments on commit 3edbdd7

Please sign in to comment.