HHT-α法（Hilber-Hughes-Taylor）
HHT-α法（Hilber-Hughes-Taylor）的理论基础
HHT-α法是什么
老师，HHT-α法是Newmark法的改进版吗？
是的。Hilber, Hughes, Taylor（1977）在Newmark法的基础上增加了数值衰减参数 $\alpha$。各求解器的默认积分法与设置取决于产品和应用。
Newmark法的问题点
Newmark法（$\beta=1/4, \gamma=1/2$）具有在线性问题中具有二阶精度且无条件稳定，但数值耗散为零。物理阻尼较小时，高频数值噪声可能持续存在。接触冲击或急剧的载荷变化容易产生高频噪声。
高频噪声不消失的话很麻烦啊。
HHT-α法解决了这个问题。在保持低频精度的同时，选择性地衰减高频成分。
HHT-α法的算法
修正后的运动方程：
参数关系：
当 $\alpha = 0$ 时就与Newmark法一致了呢。
$\alpha$ 的范围是 $-1/3 \leq \alpha \leq 0$。$\alpha = 0$: 无数值衰减（Newmark法）。$\alpha = -0.05$: 温和的高频衰减。$\alpha = -1/3$: 最大的高频衰减（衰减较强，需检查时间步误差）。
实践推荐: $\alpha = -0.05$ 左右。这样可以保持二阶精度，同时有效衰减高频噪声。
Abaqus
*DYNAMIC, ALPHA=-0.05
0.001, 1.0这只是分析步内的设置片段，模型、边界条件、载荷和分析步定义需要另行提供。 Abaqus官方：*DYNAMIC
Abaqus/Standard无接触时通常默认TRANSIENT FIDELITY（HHT-TF，默认α = −0.05）；有接触时通常默认MODERATE DISSIPATION（HHT-MD）。准静态应用使用另一积分法。
Ansys
TRNOPT,FULL,,,,,HHT ! 选择HHT积分
TINTP,0.05 ! GAMMA幅值衰减系数（第一个参数）Ansys Mechanical APDL用TRNOPT选择积分法，TINTP第一个参数GAMMA设置幅值衰减系数；它与理论中的Newmark γ并非相同含义。
Nastran
Nastran设置和参数含义随产品、版本和求解过程而异，请查阅所用求解器的官方文档。
总结
要点：
- Newmark法 + 数值衰减 $\alpha$ — 选择性衰减高频噪声
- 推荐 $\alpha = -0.05$ — 保持二阶精度的同时抑制高频
- Abaqus/Ansys中可用，须核对积分设置
- $\alpha = 0$ 时退化为Newmark法 — 无数值衰减
- 接触或急剧载荷变化产生高频噪声时，HHT-α法有效
原来很多工程师是在“不知不觉中使用”HHT-α法啊。
Abaqus/Standard的HHT-TF默认α=−0.05；有接触的问题默认HHT-MD，准静态应用使用BWE。应核对实际方案及时间步。
HHT-α是1977年诞生的数值衰减方案
HHT-α扩展Newmark积分以抑制高频数值振荡。α=0时回到平均加速度法；数值衰减随时间步及频率变化。
数值解法与实现
HHT-α法的数值特性
请详细讲解一下HHT-α法的数值特性。
| $\alpha$ | 数值衰减 | 精度 | 用途 |
|---|---|---|---|
| 0 | 无 | 2阶 | 纯粹的Newmark法 |
| -0.05 | 温和 | 2阶 | 推荐（标准） |
| -0.1 | 中等 | 2阶 | 噪声较多的问题 |
| -0.33 | 接近最大 | 2阶 | 仅限特殊用途 |
满足上面的β和γ关系时，表中各α均为二阶精度。收敛阶数不同于有限时间步下的振幅和相位误差；应针对具体问题检查时间步收敛。 Sandia: HHT integration theory
数值衰减增强时，即使保持二阶精度，响应误差也可能改变。
数值衰减也会影响低频。$|\alpha|$ 越大，低频响应衰减也越多。$\alpha = -0.05$ 是“衰减高频，低频几乎不受影响”的良好平衡点。
与广义α法的关系
Chung-Hulbert（1993）的广义α法进一步推广了HHT-α法。可以独立控制低频精度和高频衰减。Abaqus/Standard的APPLICATION=MODERATE DISSIPATION选择HHT-MD。
总结
α的选择会影响分析精度
Nastran设置和参数含义随产品、版本和求解过程而异，请查阅所用求解器的官方文档。
HHT-α法（Hilber-Hughes-Taylor）HHT-α法（Hilber-Hughes-Taylor）实践指南
HHT-α法的实践
实际分析应检查积分方法和时间步，并比较设置变化对所关注响应的影响。
应调整$\alpha$的场景
| 情况 | $\alpha$的调整 |
|---|---|
| 响应中包含高频噪声 | 加强至 $\alpha = -0.1$ |
| 低频精度很重要（颤振等） | 恢复为 $\alpha = 0$（Newmark法） |
| 接触冲击产生尖峰噪声 | 设为 $\alpha = -0.05 \sim -0.1$ |
| 从默认设置开始 | 验证时间步和响应 |
实践检查清单
参考资料：Abaqus/Standard与Ansys Mechanical APDL官方文档： Abaqus *DYNAMIC · Ansys TINTP · Ansys TRNOPT