BBO:WIP: example 1

examples/Bagirrov_Bartman_Ochal_2024.py

@@ -179,13 +179,13 @@ def friction_bound(u_nu: float) -> float:
     )
 
 
-def main(config: Config, methods, forces, contact):
+def main(config: Config, methods, forces, contact, prefix=""):
     """
     Entrypoint to example.
 
     To see result of simulation you need to call from python `main(Config().init())`.
     """
-    PREFIX = "DBBO"
+    PREFIX = f"BBO{prefix}"
     if config.force:
         to_simulate = [(m, f) for m in methods for f in forces]
     else:
@@ -304,7 +304,7 @@ def top(x):
         Y[i] = contact.sub2derivative_normal_direction(X[i], 0.0, 0.0)
     plt.plot(X, Y)
     plt.show()
-    main(config, methods, forces, contact)
+    main(config, methods, forces, contact, prefix=str(len(layers_limits)))
 
 
 def survey(config):
@@ -329,7 +329,7 @@ def survey(config):
         28e3 * kN,
         30e3 * kN,
     )[:]
-    for i in [6] : #range(1, 10 + 1, 2):
+    for i in [0, 1, 2, 3, 4, 5, 6, 7, 8]: #range(1, 10 + 1, 2):
         thickness = 3 * mm
         layers_num = i
         # layers_limits = np.logspace(0, thickness, layers_num + 1)