From adf61fc2a9bfd597a1e254dfc4bef1d5e0747809 Mon Sep 17 00:00:00 2001 From: apearce Date: Thu, 14 Mar 2024 14:54:02 +0000 Subject: [PATCH 1/3] =?UTF-8?q?=F0=9F=90=9B=20fix=20typos=20in=20SIF=20equ?= =?UTF-8?q?ations?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- process/cs_fatigue.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/process/cs_fatigue.py b/process/cs_fatigue.py index 2786f1e9b1..05fea0d241 100644 --- a/process/cs_fatigue.py +++ b/process/cs_fatigue.py @@ -173,7 +173,7 @@ def surface_stress_intensity_factor(hoop_stress, t, w, a, c, phi): m1 = 1.13e0 - 0.09e0 * a_c m2 = -0.54e0 + 0.89e0 / (0.2e0 + a_c) m3 = 0.5e0 - 1.0e0 / (0.65e0 + a_c) + 14.0e0 * (1 - a_c) ** 24.0e0 - g = 1.0e0 + (0.1e0 + 0.35e0 * a_c**2.0e0) * (1.0e0 - sin_phi) ** 2.0e0 + g = 1.0e0 + (0.1e0 + 0.35e0 * a_t**2.0e0) * (1.0e0 - sin_phi) ** 2.0e0 f_phi = (a_c**2.0e0 * cos_phi_2 + sin_phi**2.0e0) ** 0.25e0 p = 0.2e0 + a_c + 0.6e0 * a_t H1 = 1.0e0 - 0.34e0 * a_t - 0.11e0 * a * a / (c * t) @@ -205,7 +205,7 @@ def surface_stress_intensity_factor(hoop_stress, t, w, a, c, phi): H2 = ( 1.0e0 + (-2.11e0 + 0.77e0 * c_a) * a_t # G21 * a / t - + (0.55e0 - 0.72e0 * c_a * 0.75e0 + 0.14e0 * c_a * 1.5e0) * a_t_2 # G22 + + (0.55e0 - 0.72e0 * c_a ** 0.75e0 + 0.14e0 * c_a * 1.5e0) * a_t_2 # G22 ) # compute the unitless geometric correction From 32fb9427a6219dead7e3d537ac5f5c62db1ad426 Mon Sep 17 00:00:00 2001 From: apearce Date: Thu, 14 Mar 2024 15:24:50 +0000 Subject: [PATCH 2/3] =?UTF-8?q?=F0=9F=9A=A8=20fix=20black?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- process/cs_fatigue.py | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/process/cs_fatigue.py b/process/cs_fatigue.py index 05fea0d241..2022acad02 100644 --- a/process/cs_fatigue.py +++ b/process/cs_fatigue.py @@ -205,7 +205,8 @@ def surface_stress_intensity_factor(hoop_stress, t, w, a, c, phi): H2 = ( 1.0e0 + (-2.11e0 + 0.77e0 * c_a) * a_t # G21 * a / t - + (0.55e0 - 0.72e0 * c_a ** 0.75e0 + 0.14e0 * c_a * 1.5e0) * a_t_2 # G22 + + (0.55e0 - 0.72e0 * c_a**0.75e0 + 0.14e0 * c_a * 1.5e0) + * a_t_2 # G22 ) # compute the unitless geometric correction From 930205963ba4d035ebfeeccd7db7f82136bacba2 Mon Sep 17 00:00:00 2001 From: apearce Date: Thu, 14 Mar 2024 17:07:05 +0000 Subject: [PATCH 3/3] =?UTF-8?q?=E2=9C=85=20update=20test=5Fcs=5Ffatigue?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- tests/unit/test_cs_fatigue.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/tests/unit/test_cs_fatigue.py b/tests/unit/test_cs_fatigue.py index d4c7193d0f..0591873c27 100644 --- a/tests/unit/test_cs_fatigue.py +++ b/tests/unit/test_cs_fatigue.py @@ -46,7 +46,7 @@ class NcycleParam(NamedTuple): t_structural_vertical=0.0063104538380405924, t_structural_radial=0.0063104538380405924, t_crack_radial=0.0026699999999999996, - expected_n_cycle=985.2543993462306, + expected_n_cycle=1113.5875631615095, expected_t_crack_radial=0.0026699999999999996, ), ),