人工波生成
[!abstract] 概述 人工波生成（人造地震动）是地震工程中根据目标反应谱或功率谱，用数值方法合成地震动加速度时程的技术。本文系统推导三角级数法（应用最广泛的人工波生成方法）的完整理论体系，涵盖从随机过程基础到迭代拟合算法的全链路，共涉及 6 个核心公式，配 3 道算例。此外简要介绍比例法和自回归法。
1. 问题引入¶
1.1 为什么需要人工波？¶
在结构抗震分析中，时程分析法需要输入地震加速度时程曲线。然而：
- 强震记录是稀有自然事件，满足同一集系的样本数量不足（来源：earthquake-engineering-hu-yuxian↗ [p.211]）
- 现有记录未必匹配具体场地的震级、距离、场地条件
- 蒙特卡罗法统计分析需要大批满足同一统计特征的地震动取样
- 振动台试验和数值计算需要可控的输入波形
工程背景：GB 50011-2010 §5.1.2 要求特别不规则建筑、甲类建筑和超限高层采用时程分析法补充计算，应选用实际强震记录和人工模拟的加速度时程曲线，实际强震记录不应少于总数的 2/3（来源：gb50011-2010-seismic-design-code↗）。
1.2 人工波需要满足的条件¶
根据 earthquake-engineering-hu-yuxian↗ [p.211] 和 ground-motion-random-process↗，一条合格的人工波应满足地震动三要素：
| 要素 | 含义 | 控制方式 |
|---|---|---|
| 振幅 | 最大加速度 \(a_{\max}\) | 按设防烈度确定（来源：design-ground-motion-parameters↗） |
| 频谱 | 反应谱或功率谱形状 | 拟合目标设计反应谱（来源：seismic-response-spectrum↗） |
| 持时 | 强震动持续时间 \(T_d\) | 与震级 \(M\) 相关（来源：ground-motion-random-process↗ [p.41-42]） |
1.3 本文公式清单¶
| 式号 | 名称 | 用途 |
|---|---|---|
| F1 | 三角级数合成公式 | 生成初始人工波 |
| F2 | 功率谱-反应谱关系（Kaul公式） | 从目标反应谱导出功率谱 |
| F3 | 迭代修正公式 | 修正傅里叶幅值以拟合目标谱 |
| F4 | 三段型强度包络函数 | 赋予非平稳性 |
| F5 | 强震动持时公式 | 确定持时参数 |
| F6 | 比例法调幅公式 | 简单调幅方法 |
2. 前置知识梳理¶
2.1 地震动随机过程（Ground Motion Random Process）¶
强震地面运动加速度记录的集合可看作随机过程。每个记录是一次实现（样本函数），由于不同震源、震级、传播途径而具有不确定性 [p.35]。
核心特征（来源：ground-motion-random-process↗ [p.35-42]）： - 一维概率分布服从正态分布（高斯分布），均值 \(m_x = 0\) - 平稳化过程：\(a(t) = f(t) \cdot y(t)\)，其中 \(f(t)\) 为确定性包络函数，\(y(t)\) 为平稳随机过程 [p.41] - 非平稳性包括强度非平稳性（振幅随时间变化）和频率特性非平稳性；平稳化表示只考虑了前者 [p.42]
2.2 功率谱密度（Power Spectral Density）¶
功率谱密度 \(S_{xx}(\omega)\) 描述随机过程在频域中的能量分布（来源：power-spectral-density↗）。
- 与自相关函数 \(R_{xx}(\tau)\) 互为傅里叶变换对（Wiener-Khinchin 定理）[p.43]
- \(\int_0^\infty S_{xx}(\omega) d\omega = \sigma_x^2\)（均方值 = 总能量）
- 白噪声：功率谱密度在全频率范围内为常数 \(S_0\)（来源：white-noise-process↗）
2.3 傅里叶谱（Fourier Spectrum）¶
加速度时程 \(a(t)\) 的傅里叶变换：
- 傅里叶振幅谱 \(|A(\omega)|\) 表示频率 \(\omega\) 分量的振幅（来源：fourier-spectrum↗ [p.68]）
- 傅里叶相位谱 \(\phi(\omega)\) 表示相位
- 无阻尼反应谱和傅里叶谱都是地面运动频率特性的表达形式，是地面运动特性而非结构特性 [p.69]
2.4 地震反应谱（Seismic Response Spectrum）¶
单质点系在给定地震加速度作用下的最大反应随自振周期变化的曲线（来源：seismic-response-spectrum↗ [p.62]）：
- 阻尼 \(\zeta = 0\) 时谱值最大，较小阻尼就能使峰点削平很多 [p.70]
- 反应谱主要决定于地震动中最强烈的一段，不能很好地反映持续时间 [p.69]
2.5 概念关系图¶
白噪声过程 ──→ 功率谱密度 S(ω)
│
▼
Kaul公式(F2) 建立 S(ω) ↔ Sₐ(T) 关系
│
▼
三角级数法(F1) + 包络函数(F4)
│
▼
初始人工波 a₀(t)
│
▼
迭代修正(F3) 拟合目标反应谱
│
▼
最终人工波 a(t)
3. 符号表¶
| 符号 | 含义 | 单位 | 首次出现于 |
|---|---|---|---|
| \(a(t)\) | 人工地震动加速度时程 | m/s² (cm/s², gal) | 式(F1) |
| \(f(t)\) | 强度包络函数 | 无量纲 | 式(F1) |
| \(\omega_k\) | 第 \(k\) 个傅里叶分量的圆频率 | rad/s | 式(F1) |
| \(A_k\) | 第 \(k\) 个傅里叶分量的振幅 | m/s² | 式(F1) |
| \(\phi_k\) | 第 \(k\) 个傅里叶分量的相位 | rad | 式(F1) |
| \(S(\omega)\) | 功率谱密度 | m²/s³ | 式(F1) |
| \(S_a(\omega, \zeta)\) | 目标加速度反应谱 | m/s² | 式(F2) |
| \(\zeta\) | 阻尼比 | 无量纲 | 式(F2) |
| \(T_d\) | 强震动持时 | s | 式(F2) |
| \(p\) | 反应不超过反应谱值的概率 | 无量纲 | 式(F2) |
| \(\Delta\omega\) | 频率增量 | rad/s | 式(F1) |
| \(N\) | 三角级数项数 | 无量纲 | 式(F1) |
| \(a_{\max}\) | 目标最大加速度 | m/s² | 式(F6) |
| \(T_p\) | 卓越周期 | s | 式(F6) |
| \(M\) | 地震震级 | 无量纲 | 式(F5) |
| \(\varepsilon\) | 迭代收敛容许误差 | m/s² | 式(F3) |
4. 核心推导过程¶
4.1 方法概述¶
人造地震动有三种数值方法（来源：earthquake-engineering-hu-yuxian↗ [p.211]）：
| 方法 | 原理 | 特点 |
|---|---|---|
| 三角级数法 | 傅里叶级数叠加 + 迭代拟合反应谱 | 应用最普遍，可用 FFT 加速 |
| 随机脉冲法 | 地震动看作随机到达的速度脉冲叠加 | 提出最早 |
| 自回归法（ARMA） | 白噪声经滤波器输出 | 七十年代后提出 |
下面以三角级数法为核心进行完整推导。
4.2 Step 1：平稳随机过程的三角级数表示¶
假设 H1：地震动的平稳部分可用有限项三角级数逼近（来源：earthquake-engineering-hu-yuxian↗ [p.211]）
根据随机过程理论，一个平稳随机过程可以展开为三角级数（通用数学知识）。对于加速度时程，取其实数形式：
参数说明
| 符号 | 含义 | 单位 | 典型取值/查表依据 | 来源 |
|---|---|---|---|---|
| \(a(t)\) | 人工地震动加速度时程 | m/s² (gal) | 输出量 | — |
| \(f(t)\) | 强度包络函数 | 无量纲 | 三段型，见式(F4) | ground-motion-random-process↗ [p.41] |
| \(A_k\) | 第 \(k\) 个傅里叶分量振幅 | m/s² | 由式(F2)功率谱导出 | earthquake-engineering-hu-yuxian↗ [p.211] |
| \(\omega_k\) | 第 \(k\) 个分量圆频率 | rad/s | \(\omega_k = k \cdot \Delta\omega\) | — |
| \(\phi_k\) | 第 \(k\) 个分量相位 | rad | \([0, 2\pi]\) 均匀分布 | earthquake-engineering-hu-yuxian↗ [p.211] |
| \(N_1, N_2\) | 频率范围下/上限序号 | 无量纲 | \(N_1 \Delta\omega < 2\pi/T_M\), \(N_2 \Delta\omega \gg 2\pi/T_1\) | — |
| \(\Delta\omega\) | 频率增量 | rad/s | 需足够小以保证收敛 | — |
为什么相位取均匀分布？ 一般人造地震动方法均取相位谱为均匀分布，即在 \(0\) 到 \(2\pi\) 之间期望值相等。日本学者大崎顺彦（Ohsaki, 1979）注意到相位差谱的形状与强度包线 \(f(t)\) 的形状十分相似，考虑相位谱的方法值得研究（来源：earthquake-engineering-hu-yuxian↗ [p.211]）。
物理直觉：式(F1)将地震动分解为 \(N\) 个不同频率的简谐波叠加，每个简谐波有独立的随机相位 \(\phi_k\)，从而产生随机性；包络函数 \(f(t)\) 赋予强度随时间先升后降的非平稳特性。
量纲校验： $$ [a(t)] = [f(t)] \cdot [A_k] = \text{无量纲} \times \text{m/s}^2 = \text{m/s}^2 \quad \checkmark $$
4.3 Step 2：振幅 \(A_k\) 与功率谱的关系¶
振幅 \(A_k\) 由功率谱密度 \(S(\omega_k)\) 确定（来源：earthquake-engineering-hu-yuxian↗ [p.211]）：
参数说明
| 符号 | 含义 | 单位 | 典型取值/查表依据 | 来源 |
|---|---|---|---|---|
| \(S(\omega_k)\) | 功率谱密度 | m²/s³ | 由目标反应谱经式(F2)导出 | power-spectral-density↗ |
| \(\Delta\omega\) | 频率增量 | rad/s | \(\Delta\omega = 2\pi / T\)，\(T\) 为总时长 | — |
量纲校验： $$ [A_k] = \sqrt{[S] \cdot [\Delta\omega]} = \sqrt{\text{m}^2/\text{s}^3 \cdot \text{rad/s}} = \sqrt{\text{m}^2/\text{s}^4} = \text{m/s}^2 \quad \checkmark $$
4.4 Step 3：功率谱与目标反应谱的关系（Kaul 公式）¶
假设 H2：反应谱控制点处的峰值近似服从跨越率的极值分布（来源：earthquake-engineering-hu-yuxian↗ [p.211]）
功率谱 \(S(\omega)\) 与目标加速度反应谱 \(S_a(\omega, \zeta)\) 之间的近似关系为 Kaul 公式：
参数说明
| 符号 | 含义 | 单位 | 典型取值/查表依据 | 来源 |
|---|---|---|---|---|
| \(S(\omega)\) | 功率谱密度 | m²/s³ | 待求量 | power-spectral-density↗ |
| \(S_a(\omega, \zeta)\) | 目标加速度反应谱 | m/s² | 按规范设计反应谱取值 | seismic-response-spectrum↗ [p.62] |
| \(\zeta\) | 阻尼比 | 无量纲 | 一般取 0.05 | — |
| \(T_d\) | 强震动持时 | s | 由式(F5)确定 | ground-motion-random-process↗ [p.41] |
| \(p\) | 反应不超过反应谱值的概率 | 无量纲 | 一般取 \(p \geq 0.85\) | earthquake-engineering-hu-yuxian↗ [p.211] |
| \(\omega\) | 圆频率 | rad/s | \(\omega = 2\pi/T\) | — |
推导说明：式(F2)建立了频域（功率谱）与时域最大值（反应谱）之间的桥梁。其物理意义是：对于一个窄频带平稳随机过程，单自由度振子的最大反应服从瑞雷分布的极值近似，由此可以反推产生该最大反应所需的功率谱强度。
为什么反应谱和功率谱"大体等价"？ 反应谱和功率谱都是强震加速度频率特性的表达形式（来源：fourier-spectrum↗ [p.69]）。功率谱描述频域能量分布，反应谱描述各频率振子的最大反应。两者通过式(F2)近似关联，但不应同时要求——规定了反应谱就不必再规定功率谱（来源：earthquake-engineering-hu-yuxian↗ [p.211]）。
量纲校验： $$ \text{右侧分子：} [S_a^2] = (\text{m/s}^2)^2 = \text{m}^2/\text{s}^4 $$ $$ \text{右侧分母：} \left[\frac{2\pi}{\omega T_d}\right] = \frac{\text{rad}}{\text{rad/s} \cdot \text{s}} = \text{无量纲}，\quad [\ln(\cdot)] = \text{无量纲} $$ $$ \text{因此 } [S] = \frac{\text{m}^2/\text{s}^4}{\text{无量纲}} = \text{m}^2/\text{s}^4 $$
⚠ 注意：此处量纲为 \(\text{m}^2/\text{s}^4\)，与严格功率谱密度 \(\text{m}^2/\text{s}^3\) 差一个 \(\text{s}^{-1}\)。这是因为式(F2)是近似公式，其中隐含了以频率 \(f\)（Hz）而非圆频率 \(\omega\)（rad/s）为自变量的换算关系。在实际应用中，该公式作为工程近似使用，各软件实现时已做内部一致性处理（来源：earthquake-engineering-hu-yuxian↗ [p.211]，标注为近似关系）。
4.5 Step 4：迭代修正算法¶
假设 H3：通过迭代修正傅里叶幅值，可以使计算反应谱逐步逼近目标反应谱（来源：earthquake-engineering-hu-yuxian↗ [p.211]）
算法流程：
Step 4.1 — 初始化：由式(F2)从目标反应谱 \(S_a(\omega)\) 计算 \(S(\omega_k)\)，再由式(F1a)计算 \(A_k^{(0)}\)，取 \(\phi_k\) 为 \([0, 2\pi]\) 均匀分布随机数，代入式(F1)得到初始人工波 \(a_0(t)\)。
Step 4.2 — 计算 \(a_0(t)\) 的反应谱 \(S_{a_0}(\omega_i)\)（\(i = 1, 2, \ldots, M\)，\(M\) 为控制点数，一般 \(M \approx 40 \sim 60\)）。
Step 4.3 — 对比计算谱与目标谱，修正傅里叶幅值（来源：earthquake-engineering-hu-yuxian↗ [p.211]）：
其中 \(\omega_i\) 为与 \(\omega_k\) 对应的反应谱控制坐标。由于控制点数 \(M\) 少于三角级数项数 \(N\)，对每个控制坐标 \(T_i\) 需修改 \(\Delta T_i\) 区间内的全部三角级数项。修改原则可以是： - 每 \(\Delta T_i\) 段内为常比例数 - 或在 \(T_{i-1}, T_i, T_{i+1}\) 三者之间线性内插
Step 4.4 — 重复 Step 4.2–4.3，直到反应谱在控制点处的最大误差满足：
参数说明
| 符号 | 含义 | 单位 | 典型取值/查表依据 | 来源 |
|---|---|---|---|---|
| \(A_k^{(n)}\) | 第 \(n\) 次迭代的第 \(k\) 个傅里叶振幅 | m/s² | 迭代量 | — |
| \(S_a(\omega_i)\) | 目标反应谱在 \(\omega_i\) 处的值 | m/s² | 给定 | — |
| \(S_{a_n}(\omega_i)\) | 第 \(n\) 次迭代波的反应谱值 | m/s² | 计算量 | — |
| \(\varepsilon\) | 收敛容许误差 | m/s² | \(\varepsilon = 0.05 S_{a,\max}\)（目标谱最大值的 5%） | earthquake-engineering-hu-yuxian↗ [p.211] |
| \(M\) | 反应谱控制点数 | 无量纲 | \(40 \sim 60\)，可达 100 | — |
量纲校验（式 F3）： $$ \frac{[S_a(\omega_i)]}{[S_{a_n}(\omega_i)]} = \frac{\text{m/s}^2}{\text{m/s}^2} = \text{无量纲} $$ $$ [A_k^{(n+1)}] = [A_k^{(n)}] \times \text{无量纲} = \text{m/s}^2 \quad \checkmark $$
为什么长周期段收敛困难？ 在长周期段（\(T \geq 2\) s），反应谱值通常已经很小，相对误差 \(\varepsilon / S_a\) 变大，不易满足精度要求。若非对长周期地震动有特殊要求，可容许较大误差（来源：earthquake-engineering-hu-yuxian↗ [p.211]）。
4.6 Step 5：强度包络函数（非平稳性）¶
假设 H4：地震动的非平稳性可用确定性包络函数 \(f(t)\) 与平稳过程相乘表示（来源：ground-motion-random-process↗ [p.41]）
平稳化过程（来源：ground-motion-random-process↗ [p.41]）：
其中 \(y(t)\) 为平稳随机过程，\(f(t)\) 为确定性包络函数。
三段型包络函数（工程中最常用）：
参数说明
| 符号 | 含义 | 单位 | 典型取值/查表依据 | 来源 |
|---|---|---|---|---|
| \(t_1\) | 上升段结束时刻 | s | 与震级、距离相关 | ground-motion-random-process↗ [p.41] |
| \(t_2\) | 平稳段结束时刻 | s | \(t_2 = t_1 + T_{\text{strong}}\) | — |
| \(t_3\) | 总持时 | s | \(t_3 = t_2 + T_{\text{decay}}\) | — |
| \(c\) | 衰减系数 | 1/s | 经验值，与震级相关 | — |
| \(T_{\text{strong}}\) | 强震动平稳段持时 | s | 由式(F5)确定 | — |
量纲校验（式 F4）： - 上升段：\((t/t_1)^2\) = 无量纲 ✓ - 平稳段：\(1\) = 无量纲 ✓ - 下降段：\(e^{-c(t-t_2)}\)，指数 \([c \cdot t]\) = \(\text{s}^{-1} \times \text{s}\) = 无量纲 ✓
4.7 Step 6：强震动持时¶
强震动持时 \(T_d\) 与震级 \(M\) 的经验关系（来源：ground-motion-random-process↗ [p.41-42]）：
参数说明
| 符号 | 含义 | 单位 | 典型取值/查表依据 | 来源 |
|---|---|---|---|---|
| \(T_d\) | 强震动持时 | s | 经验公式 | ground-motion-random-process↗ [p.41-42] |
| \(M\) | 地震震级 | 无量纲 | 按场地地震危险性确定 | — |
⚠ 注意：式(F5)为经验公式，各系数含隐含单位，仅当 \(T_d\) 按秒、\(M\) 按里氏震级代入时结果才正确。
算例 1 — 持时计算：设 \(M = 7.0\)，求强震动持时。
结果解读：7 级地震的强震动持时约 25 秒，这与实际强震记录的持时量级一致（如 1940 El Centro 记录持时约 30 s）。（构造算例，参数为典型取值）
5. 比例法（Scale Method）¶
比例法是最简单的人工波生成方法（来源：earthquake-engineering-hu-yuxian↗ [p.211]）。
5.1 方法原理¶
选择一条地质、地震条件尽可能符合要求的真实地震动记录 \(a(t)\)，通过调整时间和加速度坐标使其满足目标参数。
参数说明
| 符号 | 含义 | 单位 | 典型取值/查表依据 | 来源 |
|---|---|---|---|---|
| \(a'(t')\) | 调整后的人工波 | m/s² | 输出量 | — |
| \(a_{\max}^0\) | 目标最大加速度 | m/s² | 按设防烈度确定 | design-ground-motion-parameters↗ |
| \(a_{\max}\) | 原始记录最大加速度 | m/s² | 已知 | — |
| \(T^0\) | 目标卓越周期 | s | 按场地条件确定 | — |
| \(T\) | 原始记录卓越周期 | s | 已知 | — |
5.2 局限性¶
- 只能调整两个比例常数，只能满足最大加速度和卓越周期两个要求
- 不能满足反应谱或频谱的其他要求
- 除 \(T^0/T = T_d^0/T_d\) 外，调整后持时不等于目标持时
- 优点：利用了真实地震动记录，可能继承实际地震动的某些未知特性
量纲校验（式 F6）： $$ [a'(t')] = \frac{\text{m/s}^2}{\text{m/s}^2} \times \text{m/s}^2 = \text{m/s}^2 \quad \checkmark $$
6. 工程应用与综合算例¶
6.1 规范要求（GB 50011-2010 §5.1.2）¶
| 要求 | 具体规定 | 来源 |
|---|---|---|
| 地震波数量 | 不少于 7 组取平均值，不少于 3 组取包络值 | elastic-time-history-analysis↗ |
| 实际记录比例 | 实际强震记录 ≥ 总数的 2/3 | time-history-analysis↗ |
| 底部剪力 | 每条波底部剪力 ≥ 反应谱法的 65% | elastic-time-history-analysis↗ |
| 平均值 | 平均值 ≥ 反应谱法的 80% | elastic-time-history-analysis↗ |
| 持时 | ≥ 5 倍基本周期且 ≥ 15 s | seismic-wave-selection↗ |
6.2 三角级数项数选择¶
根据 earthquake-engineering-hu-yuxian↗ [p.211] 的讨论：
- 项数不宜过少，否则迭代收敛慢、长周期精度差
- 至少采用几百项三角级数（\(N \geq 200\)，推荐 \(N = 2^L\) 以便使用 FFT）
- 项数增多不致使计算时间成比例增加（主要计算量在反应谱计算）
6.3 综合算例¶
算例 2 — 三角级数法生成人工波（完整流程）：
已知条件（7 度设防，Ⅱ类场地，第二组）： - 目标峰值加速度：\(a_{\max} = 35 \, \text{gal} = 0.35 \, \text{m/s}^2\)（多遇地震） - 设计反应谱：按 GB 50011-2010 表 5.1.4-2 - \(\alpha_{\max} = 0.08\)（7 度多遇） - \(T_g = 0.40 \, \text{s}\)（Ⅱ类场地，第二组） - 阻尼比 \(\zeta = 0.05\) - 震级估计：\(M = 7.0\) - 收敛容许误差：\(\varepsilon = 0.05 \times S_{a,\max}\)
Step 1 — 确定持时
由式(F5)： $$ T_d = 10^{(7.0-2.5)/3.23} = 10^{1.393} = 24.7 \, \text{s} $$
考虑上升段和下降段，取总时长 \(T_{\text{total}} = 40 \, \text{s}\)。
Step 2 — 确定包络参数
取典型三段型（构造算例，参数为典型取值）： - \(t_1 = 5 \, \text{s}\)（上升段） - \(t_2 = 5 + 24.7 \approx 30 \, \text{s}\)（平稳段结束） - \(t_3 = 40 \, \text{s}\)（总时长） - \(c = 0.3 \, \text{s}^{-1}\)（衰减系数）
Step 3 — 频率离散化
- 时间步长：\(\Delta t = 0.02 \, \text{s}\)
- 采样点数：\(N = 2048 = 2^{11}\)（满足 FFT 要求）
- 总时长：\(T = N \cdot \Delta t = 40.96 \, \text{s}\)
- 频率增量：\(\Delta\omega = 2\pi / T = 0.1533 \, \text{rad/s}\)
- 频率范围：\(\omega_1 = 0.1533 \, \text{rad/s}\)（\(T_1 = 41 \, \text{s}\)），\(\omega_{2048} = 313.7 \, \text{rad/s}\)（\(T_{2048} = 0.02 \, \text{s}\)）
Step 4 — 计算目标反应谱
按 GB 50011-2010 设计反应谱公式（\(\zeta = 0.05\)）：
其中： $$ \alpha(T) = \begin{cases} 0.08 + (0.45 - 0.08) \times T / 0.1 & T < 0.1 \, \text{s} \ 0.45 & 0.1 \leq T < 0.4 \, \text{s} \ 0.45 \times (0.4/T)^{0.9} & 0.4 \leq T < 2.0 \, \text{s} \ 0.45 \times (0.4/2.0)^{0.9} \times (2.0/T)^{1.1} \times 0.6 & T \geq 2.0 \, \text{s（近似）} \end{cases} $$
Step 5 — 由式(F2)计算功率谱
取 \(p = 0.90\)，\(T_d = 24.7 \, \text{s}\)，以 \(T = 1.0 \, \text{s}\)（\(\omega = 2\pi = 6.283 \, \text{rad/s}\)）为例：
- \(\alpha(1.0) = 0.45 \times (0.4/1.0)^{0.9} = 0.45 \times 0.435 = 0.196\)
- \(S_a(1.0) = 0.35 \times 0.196 \times 9.8 = 0.672 \, \text{m/s}^2\)
- 分母：\(-\frac{2\pi}{6.283 \times 24.7} \ln\left[-\frac{\pi}{6.283 \times 24.7} \ln 0.90\right]\)
- \(= -\frac{6.283}{155.2} \ln\left[-\frac{3.1416}{155.2} \times (-0.1054)\right]\)
- \(= -0.0405 \ln[0.002138]\)
- \(= -0.0405 \times (-6.149)\)
- \(= 0.249\)
- \(S(6.283) = 0.672^2 / 0.249 = 0.4516 / 0.249 = 1.814 \times 10^{-2} \, \text{m}^2/\text{s}^4\)
Step 6 — 由式(F1a)计算振幅
Step 7 — 生成初始波并迭代
- 生成均匀分布随机相位 \(\phi_k \in [0, 2\pi]\)
- 代入式(F1)得 \(a_0(t)\)
- 计算 \(S_{a_0}(T_i)\)，与目标 \(S_a(T_i)\) 比较
- 用式(F3)修正 \(A_k\)
- 重复至满足 \(|S_{a_n} - S_a| \leq \varepsilon\)
典型收敛情况：一般需要 5~15 次迭代即可达到 5% 精度（来源：earthquake-engineering-hu-yuxian↗ [p.211]）。
结果解读：最终得到的 \(a(t)\) 即为满足目标设计反应谱的人工合成地震动加速度时程，可直接用于结构时程分析。
量纲校验（全程）： - \(S_a\) 单位：m/s² ✓ - \(S(\omega)\) 单位：m²/s⁴（近似公式，见 4.4 节说明）✓ - \(A_k\) 单位：m/s² ✓ - \(a(t)\) 单位：m/s² ✓
6.4 自回归移动平均（ARMA）法¶
白噪声经滤波器得到过滤噪声（来源：earthquake-engineering-hu-yuxian↗ [p.211]）：
其中 \(\omega_g\) 和 \(\zeta_g\) 为场地滤波器的频率和阻尼比，\(n(t)\) 为白噪声输入。这等价于 Kanai-Tajimi 模型（金井清模型），可产生具有单峰点频谱的过滤噪声地震动。
对应的离散 ARMA(m,l) 模型（来源：earthquake-engineering-hu-yuxian↗ [p.211]）：
其中 \(z_r\) 为白噪声输入，\(\alpha_i\) 和 \(\beta_i\) 为滤波器系数。
7. 知识拓展¶
7.1 三维地震动¶
沿主轴的地震动协相关函数为零，可先制造一个竖向和两个水平向地震动，再用坐标轴旋转变为结构物轴向（来源：earthquake-engineering-hu-yuxian↗ [p.211]）。
7.2 相位谱的影响¶
大崎顺彦（Ohsaki, 1979）注意到相位差谱形状与强度包线 \(f(t)\) 形状相似。考虑相位谱的人造地震动方法是值得研究的课题（来源：earthquake-engineering-hu-yuxian↗ [p.211]）。
7.3 频率非平稳性¶
平稳化过程 \(a(t) = f(t) \cdot y(t)\) 只考虑了强度非平稳性，未考虑频率特性随时间的变化。更精细的方法应考虑频率非平稳性（来源：ground-motion-random-process↗ [p.42]）。
7.4 相关知识库页面¶
- seismic-wave-selection↗ — 地震波选波方法
- earthquake-ground-motion-selection↗ — 弹塑性分析选波原则
- power-spectral-density↗ — 功率谱密度详细理论
- fourier-spectrum↗ — 傅里叶谱与反应谱的关系
- ground-motion-random-process↗ — 地震动随机过程描述
- stochastic-ground-motion-model↗ — 随机地面运动模型
- elastic-time-history-analysis↗ — 弹性时程分析规范要求
来源与置信度¶
| 来源 | 置信度 | 贡献内容 |
|---|---|---|
| earthquake-engineering-hu-yuxian↗ [p.211] | medium | 三角级数法核心公式(F1)(F2)(F3)、迭代算法、比例法(F6)、ARMA 法、相位谱讨论 |
| aseismic-engineering-shenjumin↗ [p.35-42] | medium | 随机过程描述、平稳化过程、包络函数、持时公式(F5)、功率谱密度 |
| ground-motion-random-process↗ [p.35-42] | medium | 地震动随机过程特性、正态分布、遍历性、三段型包络 |
| seismic-response-spectrum↗ [p.62-70] | medium | 反应谱定义与性质 |
| fourier-spectrum↗ [p.67-69] | medium | 傅里叶谱与反应谱的关系 |
| power-spectral-density↗ | low | 功率谱密度通用概念（3 个来源） |
| gb50011-2010-seismic-design-code↗ | high | 时程分析规范要求、设计地震动参数 |
| elastic-time-history-analysis↗ | medium | 弹性时程分析数量与精度要求 |
| seismic-wave-selection↗ | low | 选波方法与分类（天然波 T1~T7、人工波 R1~R7） |
Confidence Notes¶
- 式(F1)三角级数法：confidence = medium，知识库中两本教材（胡聿贤、沈聚敏）均有明确描述，公式来源清晰 [p.211]
- 式(F2)Kaul 公式：confidence = medium，OCR 转换的公式细节需人工复核，但公式的存在性和物理意义在两本教材中均有确认
- 式(F3)迭代修正：confidence = medium，迭代原理明确，具体修正方法（常比例/线性内插）来自胡聿贤 [p.211]
- 式(F4)三段型包络：confidence = medium，沈聚敏教材有明确描述 [p.41-42]
- 式(F5)持时公式：confidence = low，仅一个来源（沈聚敏 [p.41-42]），且为经验公式，不同研究者给出不同形式
- 式(F6)比例法：confidence = medium，胡聿贤 [p.211] 有明确描述
- 算例参数：设计反应谱参数来自 GB 50011-2010（高置信），包络参数和持时公式参数为构造的典型值
- 量纲待核：式(F2)的量纲严格来说为 m²/s⁴ 而非 m²/s³，可能是 OCR 转换误差或原书近似公式的固有特征，建议核对原始文献
- OCR 质量：胡聿贤《地震工程学》经 MinerU 从扫描型 PDF 转换，公式和表格的 OCR 精度需人工复核
⚠ 知识库覆盖评估：知识库中关于人工波生成的直接信息来自 2 个来源（胡聿贤、沈聚敏），覆盖了核心方法论。但缺少以下内容的详细资料： - 具体的 Kanai-Tajimi/Clough-Penzien 场地谱模型公式 - 现代选波软件（如 SeismoSelect、PEER）的算法细节 - GB 50011-2010 设计反应谱的完整公式表
若需更深入的专题（如场地相关反应谱、谱匹配法），建议 ingest 相关教材或规范。