Update cs_fatigue.py

ukaea · Nov 2, 2023 · 3674544 · 3674544
1 parent 74e12bb
commit 3674544
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/process/cs_fatigue.py b/process/cs_fatigue.py
@@ -110,7 +110,7 @@ def embedded_stress_intensity_factor(hoop_stress, t, w, a, c, phi):
             f_phi = (a_c**2.0e0 * cos_phi_2 + sin_phi_2) ** 0.25e0
         else:  # elif a > c:
             c_a = c / a
-            Q = 1.0e0 + 1.464e0 * (c_a) ** 1.65e0
+            Q = 1.0e0 + 1.464e0 * c_a**1.65e0
             m1 = numpy.sqrt(c_a)
             f_phi = (c_a**2.0e0 * sin_phi_2 + cos_phi_2) ** 0.25e0
 
@@ -169,18 +169,18 @@ def surface_stress_intensity_factor(hoop_stress, t, w, a, c, phi):
             # reuse of calc
             a_c = a / c
 
-            Q = 1.0e0 + 1.464e0 * (a_c) ** 1.65e0
+            Q = 1.0e0 + 1.464e0 * a_c**1.65e0
             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
-            f_phi = (a_c**2.0e0 * cos_phi_2 + (sin_phi) ** 2.0e0) ** 0.25e0
+            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)
             H2 = (
                 1.0
                 + (-1.22e0 - 0.12e0 * a_c) * a_t  # G21 * a / t
-                + (0.55e0 - 1.05e0 * (a_c) ** 0.75e0 + 0.47e0 * a_c**1.5e0)  # G22
+                + (0.55e0 - 1.05e0 * a_c**0.75e0 + 0.47e0 * a_c**1.5e0)  # G22
                 * a_t_2
             )
         else:  # elif a > c: