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$）は線形問題では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$ 程度。これで2次精度を維持しつつ高周波ノイズを効果的に減衰。
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（第1引数）Ansys Mechanical APDLのTRNOPTで積分法を選び、TINTPの第1引数GAMMAで振幅減衰係数を設定します。これは理論式のNewmark γと同一名称でも同じ意味とは限りません。
Nastran
Nastran系の設定方法とパラメータの意味は製品・バージョン・解析手順によって異なります。利用中のソルバーの公式マニュアルを確認してください。
まとめ
要点：
- Newmark法 + 数値減衰 $\alpha$ — 高周波ノイズの選択的減衰
- $\alpha = -0.05$ が推奨 — 2次精度を維持しつつ高周波を抑制
- Abaqusでは用途に応じた時間積分法を選択
- $\alpha = 0$ でNewmark法に退化 — 数値減衰なし
- 接触や急な荷重変化で高周波ノイズが出たらHHT-αが有効
多くのエンジニアがHHT-α法を「知らずに使っている」んですね。
Abaqus/StandardのHHT-TFではα=−0.05が既定ですが、接触を伴う問題ではHHT-MD、準静的用途ではBWEが選択されます。設定と時間刻みの影響を確認してください。
HHT-αは1977年生まれの数値減衰スキーム
HHT-αはNewmark法を拡張し、高周波の数値振動を抑える方法です。α=0では平均加速度法に戻ります。減衰の影響は時間刻みや周波数に依存します。
HHT-α法（Hilber-Hughes-Taylor）の数値計算手法
HHT-α法の数値特性
HHT-α法の数値特性を詳しく教えてください。
| $\alpha$ | 数値減衰 | 精度 | 用途 |
|---|---|---|---|
| 0 | なし | 2次 | 純粋なNewmark法 |
| -0.05 | 穏やか | 2次 | 推奨（標準） |
| -0.1 | 中程度 | 2次 | ノイズが多い問題 |
| -0.33 | 最大に近い | 2次 | 特殊用途のみ |
上のβ・γの関係を満たす場合、表のαはいずれも2次精度です。次数と有限の時間刻みで生じる振幅・位相誤差は別で、解析ごとに時間刻みを変えて収束を確認してください。 Sandia: HHT integration theory
数値減衰を強めると、2次精度のままでも応答の誤差は変わるのですね。
数値減衰は低周波にも影響する。$|\alpha|$ が大きいほど低周波の応答も減衰される。$\alpha = -0.05$ が「高周波は減衰、低周波はほぼ影響なし」の良いバランス。
一般化α法との関係
Chung-Hulbert（1993）の一般化α法はHHT-α法をさらに一般化。低周波の精度と高周波の減衰を独立に制御できる。Abaqus/StandardのAPPLICATION=MODERATE DISSIPATIONはHHT-MDを選択します。
まとめ
αの選び方で解析精度が変わる
Nastran系の設定方法とパラメータの意味は製品・バージョン・解析手順によって異なります。利用中のソルバーの公式マニュアルを確認してください。
HHT-α法（Hilber-Hughes-Taylor）の実務適用
HHT-α法の実務
実務では積分法と刻み幅を確認し、関心のある応答に対する結果の感度を比較する。
$\alpha$を調整すべき場面
| 状況 | $\alpha$の調整 |
|---|---|
| 高周波ノイズが応答に含まれる | $\alpha = -0.1$ に強化 |
| 低周波の精度が重要（フラッター等） | $\alpha = 0$（Newmark法）に戻す |
| 接触の衝撃でスパイクノイズ | $\alpha = -0.05 \sim -0.1$ |
| 既定設定で開始 | 時間刻みと結果を検証 |
実務チェックリスト
「デフォルトで使い、問題が出たら調整」が実務アプローチですね。
高周波ノイズと必要な応答の両方を評価し、時間刻みと数値減衰の影響を確認する。
制震建物の地震応答でHHT-αが活躍
東京スカイツリー（2012年竣工）の設計では制震システム（心柱制震）の動的解析にHHT-αベースの陰解法時刻歴解析が採用されたと報告されている。地震波はEl Centro 1940やJMA神戸1995を用い、α＝−0.1で50秒間・時間刻み0.01sの解析を実施。高周波ノイズを除去しながら長周期成分（固有周期8s超）の応答を正確に捉えることがHHT採用の主な理由とされる。
HHT-α法（Hilber-Hughes-Taylor）のソフトウェア比較
HHT-α法のツール
NastranのデフォルトはNewmark法（$\alpha=0$）なんですか。
Nastran系の設定方法とパラメータの意味は製品・バージョン・解析手順によって異なります。利用中のソルバーの公式マニュアルを確認してください。
選定ガイド
陰解法過渡解析ソルバーの比較
時間刻み、関心周波数とエネルギーの変化を確認し、設定を変えた比較解析で結果の感度を評価します。
HHT-α法（Hilber-Hughes-Taylor）の先端研究
エネルギー保存型積分法
エネルギー-運動量保存法（Simo-Tarnow, 1992）は非線形問題でのエネルギーバランスを厳密に保存する時間積分法。HHT-αは線形問題では優秀だが、非線形の大変形で数値的なエネルギードリフトが起きることがある。
適応的$\alpha$
局所誤差に基づいて$\alpha$を時間ステップごとに自動調整する適応的HHT-α法。ノイズが多い段階では$|\alpha|$を大きく、安定した段階では小さくする。
まとめ
一般化αはHHTをさらに発展させた
時間刻み、関心周波数とエネルギーの変化を確認し、設定を変えた比較解析で結果の感度を評価します。
HHT-α法（Hilber-Hughes-Taylor）のトラブル対応
高周波ノイズが消えない
$|\alpha|$を大きくする（-0.05 → -0.1）。それでも消えない場合は$\Delta t$を小さくする。
低周波の応答が減衰しすぎる
$|\alpha|$が大きすぎる。$\alpha$を0に近づける（-0.1 → -0.05 → 0）。
エネルギーが保存されない（非線形）
大変形の非線形問題でHHT-αの数値減衰がエネルギーを過度に散逸することがある。$\alpha = 0$（Newmark法、数値的な減衰なし）に切り替えるか、$\Delta t$を小さくする。
まとめ
過渡応答の初期振動はαの過小が原因
HHT-α解析で解析開始直後に高周波振動が発生し収まらない場合、αが0に近すぎて数値減衰が働いていない可能性が高い。初期条件（変位・速度）の不連続性も高周波を励起するため、初期ステップを小さく設定（Δt=Δt_stable×0.1程度）して緩やかに立ち上げるか、*CONTROLS,PARAMETERS内で初期条件のスムージングを有効にする方法が有効である。
関連トピック
なった
詳しく
報告
参考資料：Abaqus/StandardとAnsys Mechanical APDLの公式マニュアル： Abaqus *DYNAMIC · Ansys TINTP · Ansys TRNOPT