Skip to content

Commit 3edbdd7

Browse files
joaquinirazabaljoaquinirazabal
authored
Add new joints to Dam app and remove old ones (#982)
1 parent af19628 commit 3edbdd7

3 files changed

Lines changed: 56 additions & 126 deletions

File tree

‎kratos.gid/apps/Dam/write/write.tcl‎

Lines changed: 10 additions & 18 deletions
Original file line numberDiff line numberDiff line change
@@ -69,29 +69,21 @@ proc ::Dam::write::UpdateMaterials { } {
6969
set constlaw [dict get $props ConstitutiveLaw]
7070
# Modificar la ley constitutiva
7171
set newconstlaw $constlaw
72-
if {$constlaw eq "ElasticCohesive2DPlaneStress"} {
73-
dict set matdict $mat THICKNESS 1.0000E+00
74-
set newconstlaw "ElasticCohesive2DLaw"
75-
}
7672
if {$constlaw eq "ElasticCohesive2DPlaneStrain"} {
77-
dict set matdict $mat THICKNESS 1.0000E+00
73+
dict set matdict $mat THICKNESS 1.0
7874
set newconstlaw "ElasticCohesive2DLaw"
7975
}
80-
if {$constlaw eq "BilinearCohesive2DPlaneStress"} {
81-
dict set matdict $mat THICKNESS 1.0000E+00
82-
set newconstlaw "BilinearCohesive2DLaw"
83-
}
84-
if {$constlaw eq "BilinearCohesive2DPlaneStrain"} {
85-
dict set matdict $mat THICKNESS 1.0000E+00
86-
set newconstlaw "BilinearCohesive2DLaw"
76+
if {$constlaw eq "ElasticCohesive2DPlaneStress"} {
77+
dict set matdict $mat THICKNESS 1.0
78+
set newconstlaw "ElasticCohesive2DLaw"
8779
}
88-
if {$constlaw eq "ExponentialCohesive2DPlaneStress"} {
89-
dict set matdict $mat THICKNESS 1.0000E+00
90-
set newconstlaw "ExponentialCohesive2DLaw"
80+
if {$constlaw eq "IsotropicDamageCohesive2DPlaneStress"} {
81+
dict set matdict $mat THICKNESS 1.0
82+
set newconstlaw "IsotropicDamageCohesive2DLaw"
9183
}
92-
if {$constlaw eq "ExponentialCohesive2DPlaneStrain"} {
93-
dict set matdict $mat THICKNESS 1.0000E+00
94-
set newconstlaw "ExponentialCohesive2DLaw"
84+
if {$constlaw eq "IsotropicDamageCohesive2DLawPlaneStrain"} {
85+
dict set matdict $mat THICKNESS 1.0
86+
set newconstlaw "IsotropicDamageCohesive2DLaw"
9587
}
9688

9789
dict set matdict $mat CONSTITUTIVE_LAW_NAME $newconstlaw

‎kratos.gid/apps/Dam/xml/ConstitutiveLaws.xml‎

Lines changed: 34 additions & 81 deletions
Original file line numberDiff line numberDiff line change
@@ -61,43 +61,27 @@
6161
<!--3D Interface laws-->
6262
<CLaw n="ElasticCohesive3DLaw" pn="Joint Elastic Cohesive 3D" Type="Interface3D" behaviour="Elastic"
6363
StrainSize="6" AnalysisType="Linear,Non-Linear"
64-
help="Elastic law (no breakage)for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False" DamageType="Local,NonLocal" >
64+
help="Elastic law (no breakage) for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False" DamageType="Local,NonLocal" >
6565
<inputs>
66-
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
6766
<parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
68-
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e03"/>
69-
</inputs>
70-
<outputs>
71-
</outputs>
72-
</CLaw>
73-
74-
<CLaw n="BilinearCohesive3DLaw" pn="Joint Bilinear Cohesive 3D" Type="Interface3D" behaviour="Elastic"
75-
StrainSize="6" AnalysisType="Linear,Non-Linear"
76-
help="Bilinear Cohesive law for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False" DamageType="Local,NonLocal" >
77-
<inputs>
78-
<parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="2.6e10" />
79-
<parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
80-
<parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
81-
<parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
82-
<parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
83-
<parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
84-
<parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
67+
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
68+
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
8569
</inputs>
8670
<outputs>
87-
<parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
8871
</outputs>
8972
</CLaw>
9073

91-
<CLaw n="ExponentialCohesive3DLaw" pn="Joint Exponential Cohesive 3D" Type="Interface3D" behaviour="Elastic"
74+
<CLaw n="IsotropicDamageCohesive3DLaw" pn="Joint Isotropic Damage Cohesive 3D" Type="Interface3D" behaviour="Damage"
9275
StrainSize="6" AnalysisType="Linear,Non-Linear"
93-
help="Exponential Cohesive law for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False" DamageType="Local,NonLocal" >
76+
help="Isotropic Damage Cohesive law for joints in 3D problems" Dimension="3D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AllowsAnisotropy="False" DamageType="Local,NonLocal" >
9477
<inputs>
95-
<parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
96-
<parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
97-
<parameter n="FRACTURE_ENERGY" pn="Fracture Energy" v="90.0e06"/>
98-
<parameter n="SHEAR_FRACTURE_ENERGY" pn="Shear Fracture Energy" v="350.0e06"/>
99-
<parameter n="STATE_VARIABLE" pn="State Variable" v="1.0e-15"/>
100-
<parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
78+
<parameter n="DAMAGE_EVOLUTION_LAW" pn="Damage Evolution Law" v="1"/>
79+
<parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
80+
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
81+
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
82+
<parameter n="TENSILE_STRENGTH" pn="Tensile Strength" unit_magnitude="P" units="Pa" v="1e03"/>
83+
<parameter n="FRACTURE_ENERGY" pn="Fracture Energy" unit_magnitude="Energy/L^2" units="J/m^2" v="1e02"/>
84+
<parameter n="BETA_EQSTRAIN_SHEAR_FACTOR" pn="Beta Factor" v="1.0"/>
10185
</inputs>
10286
<outputs>
10387
<parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
@@ -218,9 +202,9 @@
218202
<CLaw n="ElasticCohesive2DPlaneStrain" pn="Simplified Joint Plane Strain" Type="InterfaceStrain" behaviour="Elastic" StrainSize="3" DamageType="Local,NonLocal"
219203
help="Simplified law for joints in Plane Strain problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AnalysisType="Linear,Non-Linear">
220204
<inputs>
221-
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
222205
<parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
223-
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e03"/>
206+
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
207+
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
224208
</inputs>
225209
<outputs>
226210
</outputs>
@@ -229,74 +213,43 @@
229213
<CLaw n="ElasticCohesive2DPlaneStress" pn="Simplified Joint Plane Stress" Type="InterfaceStress" behaviour="Elastic" StrainSize="3" DamageType="Local,NonLocal"
230214
help="Simplified law for joints in Plane Stress problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AnalysisType="Linear,Non-Linear">
231215
<inputs>
232-
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
233216
<parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
234-
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e03"/>
235-
</inputs>
236-
<outputs>
237-
</outputs>
238-
</CLaw>
239-
240-
<CLaw n="BilinearCohesive2DPlaneStrain" pn="Joint Bilinear Cohesive Plane Strain" Type="InterfaceStrain" behaviour="Elastic" StrainSize="3" DamageType="Local,NonLocal"
241-
help="Bilinear Cohesive law for joints in Plane Strain problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AnalysisType="Linear,Non-Linear">
242-
<inputs>
243-
<parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="2.6e10" />
244-
<parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
245-
<parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
246-
<parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
247-
<parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
248-
<parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
249-
<parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
250-
</inputs>
251-
<outputs>
252-
<parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
253-
</outputs>
254-
</CLaw>
255-
256-
<CLaw n="BilinearCohesive2DPlaneStress" pn="Joint Bilinear Cohesive Plane Stress" Type="InterfaceStress" behaviour="Elastic" StrainSize="3" DamageType="Local,NonLocal"
257-
help="Bilinear Cohesive law for joints in Plane Stress problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AnalysisType="Linear,Non-Linear">
258-
<inputs>
259-
<parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="2.6e10" />
260-
<parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
261-
<parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
262-
<parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
263-
<parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
264-
<parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
265-
<parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
217+
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
218+
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
266219
</inputs>
267220
<outputs>
268-
<parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
269221
</outputs>
270222
</CLaw>
271223

272-
<CLaw n="ExponentialCohesive2DPlaneStrain" pn="Joint Exponential Cohesive Plane Strain" Type="InterfaceStrain" behaviour="Elastic" StrainSize="3" DamageType="Local,NonLocal"
273-
help="Exponential Cohesive law for joints in Plane Strain problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AnalysisType="Linear,Non-Linear">
224+
<CLaw n="IsotropicDamageCohesive2DPlaneStrain" pn="Isotropic Damage Cohesive Plane Strain" Type="InterfaceStrain" behaviour="Elastic" StrainSize="3" DamageType="Local,NonLocal"
225+
help="Isotropic Damage Cohesive law for joints in Plane Strain problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AnalysisType="Linear,Non-Linear">
274226
<inputs>
275-
<parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
276-
<parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
277-
<parameter n="FRACTURE_ENERGY" pn="Fracture Energy" v="90.0e06"/>
278-
<parameter n="SHEAR_FRACTURE_ENERGY" pn="Shear Fracture Energy" v="350.0e06"/>
279-
<parameter n="STATE_VARIABLE" pn="State Variable" v="1.0e-15"/>
280-
<parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
227+
<parameter n="DAMAGE_EVOLUTION_LAW" pn="Damage Evolution Law" v="1"/>
228+
<parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
229+
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
230+
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
231+
<parameter n="TENSILE_STRENGTH" pn="Tensile Strength" unit_magnitude="P" units="Pa" v="1e03"/>
232+
<parameter n="FRACTURE_ENERGY" pn="Fracture Energy" unit_magnitude="Energy/L^2" units="J/m^2" v="1e02"/>
233+
<parameter n="BETA_EQSTRAIN_SHEAR_FACTOR" pn="Beta Factor" v="1.0"/>
281234
</inputs>
282235
<outputs>
283236
<parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
284237
</outputs>
285238
</CLaw>
286239

287-
<CLaw n="ExponentialCohesive2DPlaneStress" pn="Joint Exponential Cohesive Plane Stress" Type="InterfaceStress" behaviour="Elastic" StrainSize="3" DamageType="Local,NonLocal"
288-
help="Exponential Cohesive law for joints in Plane Stress problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AnalysisType="Linear,Non-Linear">
240+
<CLaw n="IsotropicDamageCohesive2DPlaneStress" pn="Isotropic Damage Cohesive Plane Stress" Type="InterfaceStress" behaviour="Elastic" StrainSize="3" DamageType="Local,NonLocal"
241+
help="Isotropic Damage Cohesive law for joints in Plane Stress problems" Dimension="2D" LargeDeformation="False" RequiresLocalAxes="False" HybridType="False" AnalysisType="Linear,Non-Linear">
289242
<inputs>
290-
<parameter n="YIELD_STRESS" pn="Yield Stress" v="0.3e06"/>
291-
<parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
292-
<parameter n="FRACTURE_ENERGY" pn="Fracture Energy" v="90.0e06"/>
293-
<parameter n="SHEAR_FRACTURE_ENERGY" pn="Shear Fracture Energy" v="350.0e06"/>
294-
<parameter n="STATE_VARIABLE" pn="State Variable" v="1.0e-15"/>
295-
<parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
243+
<parameter n="DAMAGE_EVOLUTION_LAW" pn="Damage Evolution Law" v="1"/>
244+
<parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
245+
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
246+
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
247+
<parameter n="TENSILE_STRENGTH" pn="Tensile Strength" unit_magnitude="P" units="Pa" v="1e03"/>
248+
<parameter n="FRACTURE_ENERGY" pn="Fracture Energy" unit_magnitude="Energy/L^2" units="J/m^2" v="1e02"/>
249+
<parameter n="BETA_EQSTRAIN_SHEAR_FACTOR" pn="Beta Factor" v="1.0"/>
296250
</inputs>
297251
<outputs>
298252
<parameter n="DAMAGE_VARIABLE" pn="Damage" v="False" />
299253
</outputs>
300254
</CLaw>
301-
302255
</ConstitutiveLaws>

‎kratos.gid/apps/Dam/xml/Materials.xml‎

Lines changed: 12 additions & 27 deletions
Original file line numberDiff line numberDiff line change
@@ -19,38 +19,23 @@
1919
</inputs>
2020
</Material>
2121

22-
<Material n="Simplified-Joint" MaterialType="Dam" help="">
22+
<Material n="Elastic Cohesive Joint" MaterialType="Dam" help="">
2323
<inputs>
24-
<parameter n="DENSITY" pn="Density" v="2400"/>
25-
<parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="2.6e10"/>
26-
<parameter n="YIELD_STRESS" pn="Yield Stress" unit_magnitude="P" units="Pa" v="20e06"/>
27-
<parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
24+
<parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
25+
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="0.3e06" />
26+
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
2827
</inputs>
2928
</Material>
3029

31-
<Material n="Bilinear Strong-Joint" MaterialType="Dam" help="">
30+
<Material n="Isotropic Damage Cohesive Joint" MaterialType="Dam" help="">
3231
<inputs>
33-
<parameter n="DENSITY" pn="Density" v="2400"/>
34-
<parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="35.0e9" />
35-
<parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
36-
<parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
37-
<parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
38-
<parameter n="YIELD_STRESS" pn="Yield Stress" v="3.5e8"/>
39-
<parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
40-
<parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
41-
</inputs>
42-
</Material>
43-
44-
<Material n="Bilinear Weak-Joint" MaterialType="Dam" help="">
45-
<inputs>
46-
<parameter n="DENSITY" pn="Density" v="2400"/>
47-
<parameter n="YOUNG_MODULUS" pn="Young Modulus" unit_magnitude="P" units="Pa" v="1.0" />
48-
<parameter n="POISSON_RATIO" pn="Poisson Ratio" v="0.20"/>
49-
<parameter n="INITIAL_JOINT_WIDTH" pn="Initial Joint Width" v="1e-03"/>
50-
<parameter n="CRITICAL_DISPLACEMENT" pn="Critical Displacement" v="0.05"/>
51-
<parameter n="YIELD_STRESS" pn="Yield Stress" v="1.0"/>
52-
<parameter n="DAMAGE_THRESHOLD" pn="Damage Threshold" v="0.04"/>
53-
<parameter n="FRICTION_COEFFICIENT" pn="Friction Coefficient" v="0.4"/>
32+
<parameter n="DAMAGE_EVOLUTION_LAW" pn="Damge Evolution Law" v="1"/>
33+
<parameter n="NORMAL_STIFFNESS" pn="Normal stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
34+
<parameter n="SHEAR_STIFFNESS" pn="Shear stiffness" unit_magnitude="F/L" units="N/m" v="1e06" />
35+
<parameter n="PENALTY_STIFFNESS" pn="Penalty stiffness" v="1e06"/>
36+
<parameter n="TENSILE_STRENGTH" pn="Tensile Strength" unit_magnitude="P" units="Pa" v="1e03"/>
37+
<parameter n="FRACTURE_ENERGY" pn="Fracture Energy" unit_magnitude="Energy/L^2" units="J/m^2" v="1e02"/>
38+
<parameter n="BETA_EQSTRAIN_SHEAR_FACTOR" pn="Beta Equivalent Strain Shear Factor" v="1.0"/>
5439
</inputs>
5540
</Material>
5641

0 commit comments

Comments
 (0)