From 260bd524e3211b87b2c004001de9aae70633863d Mon Sep 17 00:00:00 2001 From: Piotr Bartman Date: Thu, 11 May 2023 20:04:35 +0200 Subject: [PATCH] BBO: reproduced --- examples/Bagirov_Bartman_Ochal_2023.py | 4 ++-- examples/Sofonea_Ochal_Bartman_2023.py | 6 ++++++ 2 files changed, 8 insertions(+), 2 deletions(-) diff --git a/examples/Bagirov_Bartman_Ochal_2023.py b/examples/Bagirov_Bartman_Ochal_2023.py index a4879346..d3e2ae2e 100644 --- a/examples/Bagirov_Bartman_Ochal_2023.py +++ b/examples/Bagirov_Bartman_Ochal_2023.py @@ -46,7 +46,7 @@ def potential_normal_direction(u_nu: float) -> float: if u_nu < 1 * mm: return k10 * u_nu ** 2 + k11 * u_nu if u_nu < 2 * mm: - return k20 * u_nu ** 2 + k21 * u_nu + 4.0 + return k20 * u_nu ** 2 + k21 * u_nu + 4 return 16 @staticmethod @@ -95,7 +95,7 @@ def friction_bound(u_nu: float) -> float: def main(save: bool = False): setup = StaticSetup(mesh_type="cross") - for method in ("Powell", "BFGS", "qsm")[2:]: + for method in ("Powell", "BFGS", "CG", "qsm"): for force in np.arange(25e3 * kN, 26e3 * kN + 1, 1e3 * kN) * surface: def outer_forces(x, t=None): if x[1] >= 0.0099: diff --git a/examples/Sofonea_Ochal_Bartman_2023.py b/examples/Sofonea_Ochal_Bartman_2023.py index 5f1054e3..59134f06 100644 --- a/examples/Sofonea_Ochal_Bartman_2023.py +++ b/examples/Sofonea_Ochal_Bartman_2023.py @@ -161,6 +161,12 @@ def zero_relaxation(t=None): setup.relaxation = examples[name]["relaxation"] runner = QuasistaticRelaxation(setup, solving_method="schur") + bid = runner.body.mesh.contact_boundary[:, 0] + eid = runner.body.mesh.contact_boundary[:, 1] + bx = runner.body.mesh.initial_nodes[bid][:, 0] + ex = runner.body.mesh.initial_nodes[eid][:, 0] + length = np.max(ex - bx) + print(length) states = runner.solve( n_steps=examples[name]["n_steps"],