Fixing bug with Tendon stiffness

mbonatte · Feb 9, 2024 · c4f5510 · c4f5510
1 parent 5688bcb
commit c4f5510
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/secan/geometry.py b/secan/geometry.py
@@ -391,7 +391,8 @@ def get_normal_stress(self, strain=0):
     def get_stiffness(self, e0, k, center):
         strain = e0 + k * (center - self.center[1])
         if strain > self.strain_ULS:
-            return 0
+            return np.array(([0, 0],
+                             [0, 0]))
         strain += self.initial_strain  #This still has to be tested
         normal = self.area * self.material.get_stiff(strain)
         dist = (center-self.center[1])