建立与设计加速度反应谱等价的地震动功率谱python代码

好的，以下是一个简单的 Python 代码，用于将加速度反应谱转换为地震动功率谱。代码基于 NumPy 和 Matplotlib 库，需要先安装这些库。

import numpy as np
import matplotlib.pyplot as plt

# 加速度反应谱
Sa = np.array([1.5, 2.3, 3.8, 5.5, 7.8, 9.8, 10.9, 10.3, 8.5, 6.5, 5.0])
# 对应周期
T = np.array([0.01, 0.02, 0.05, 0.10, 0.20, 0.30, 0.40, 0.50, 0.80, 1.00, 2.00])

# 将 Sa 转换为加速度谱加权地震动功率谱
Sd = Sa * (2 * np.pi / T) ** 2

# 画图
plt.loglog(T, Sd)
plt.xlabel('Period (s)')
plt.ylabel('Power Spectral Density (m^2/s^3)')
plt.show()

这段代码将 `Sa` 数组作为输入，其中包含加速度反应谱值， `T` 数组包含对应的周期。代码将 `Sa` 转换为加速度谱加权地震动功率谱，并将结果绘制为对数-对数图形。