fix:Timoshenko beam example nu scaled for new dissipation module

GazzolaLab · May 30, 2022 · b13fd45 · b13fd45
1 parent 6ce953c
commit b13fd45
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 6 deletions.
diff --git a/examples/TimoshenkoBeamCase/convergence_timoshenko.py b/examples/TimoshenkoBeamCase/convergence_timoshenko.py
@@ -20,7 +20,7 @@ class TimoshenkoBeamSimulator(BaseSystemCollection, Constraints, Forcing):
 
 # Options
 PLOT_FIGURE = True
-SAVE_FIGURE = False
+SAVE_FIGURE = True
 SAVE_RESULTS = False
 ADD_UNSHEARABLE_ROD = False
 
@@ -38,10 +38,11 @@ def simulate_timoshenko_beam_with(
     base_length = 3.0
     base_radius = 0.25
     density = 5000
-    nu = 0.1
+    nu = 0.1 / density / (np.pi * base_radius**2)
     E = 1e6
     # For shear modulus of 1e4, nu is 99!
     poisson_ratio = 99
+    shear_modulus = E / (poisson_ratio + 1.0)
 
     shearable_rod = CosseratRod.straight_rod(
         n_elem,
@@ -53,7 +54,7 @@ def simulate_timoshenko_beam_with(
         density,
         nu,
         E,
-        poisson_ratio,
+        shear_modulus=shear_modulus,
     )
 
     timoshenko_sim.append(shearable_rod)

diff --git a/examples/TimoshenkoBeamCase/timoshenko.py b/examples/TimoshenkoBeamCase/timoshenko.py
@@ -19,9 +19,9 @@ class TimoshenkoBeamSimulator(BaseSystemCollection, Constraints, Forcing):
 
 # Options
 PLOT_FIGURE = True
-SAVE_FIGURE = False
+SAVE_FIGURE = True
 SAVE_RESULTS = False
-ADD_UNSHEARABLE_ROD = True
+ADD_UNSHEARABLE_ROD = False
 
 # setting up test params
 n_elem = 100
@@ -32,7 +32,7 @@ class TimoshenkoBeamSimulator(BaseSystemCollection, Constraints, Forcing):
 base_radius = 0.25
 base_area = np.pi * base_radius ** 2
 density = 5000
-nu = 0.1
+nu = 0.1 / density / base_area
 E = 1e6
 # For shear modulus of 1e4, nu is 99!
 poisson_ratio = 99