From 1a4fea3c8c5572e0400652848dc27bf56c05aa22 Mon Sep 17 00:00:00 2001 From: shinanca Date: Thu, 25 Jan 2024 16:00:30 +0300 Subject: [PATCH] More precise method for damping ration calculation --- K-ShakeTune/scripts/common_func.py | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/K-ShakeTune/scripts/common_func.py b/K-ShakeTune/scripts/common_func.py index 7c74109..203054b 100644 --- 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