reijii/klippain-shaketune-telegramm
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

More precise method for damping ration calculation

Browse Source
This commit is contained in:
shinanca
2024-01-25 16:00:30 +03:00
parent eab10ce5da
commit 1a4fea3c8c
1 changed files with 4 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
K-ShakeTune/scripts/common_func.py
View File

@@ -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]) 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 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 return fr, zeta, max_power_index