Half-Quadratic Gain method for damping ratio estimation

that should be more precise than the Half-Power method for higher damping ratio values (above 0.05)
Frix-x · Jan 29, 2024 · 9fa07a1 · 9fa07a1
2 parents eab10ce + 1a4fea3
commit 9fa07a1
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/K-ShakeTune/scripts/common_func.py b/K-ShakeTune/scripts/common_func.py
@@ -84,7 +84,10 @@ def compute_mechanical_parameters(psd, freqs):
     freq_above_half_power = freqs[idx_above - 1] + (half_power - psd[idx_above - 1]) * (freqs[idx_above] - freqs[idx_above - 1]) / (psd[idx_above] - psd[idx_above - 1])
 
     bandwidth = freq_above_half_power - freq_below_half_power
-    zeta = bandwidth / (2 * fr)
+    bw1 = math.pow(bandwidth/fr,2)
+    bw2 = math.pow(bandwidth/fr,4)
+
+    zeta = math.sqrt(0.5-math.sqrt(1/(4+4*bw1-bw2)))
 
     return fr, zeta, max_power_index