[論文レビュー] Hessian-Enhanced Alternating Frequency/Time method for Computing Resonance Backbone Curves

Wei Wen, Weikai Qi+1 more|arXiv (Cornell University)|Feb 6, 2026

Bladed Disk Vibration Dynamics被引用数 0

ひとこと要約

本論文は、Lagrange乗数形式と拡張AFTを用いて解析 Hessian を導出し、Harmonic Balance 内で共振バックボーン曲線を計算する Hessian強化最適化フレームワークを提案する。2自由度・ビーム・ブリンクスモデルで実証。

ABSTRACT

Computing resonance and anti-resonance backbone curves in complex nonlinear mechanical systems is of high engineering relevance but remains computationally challenging, especially for large finite-element (FE) models. Existing manifold-based approaches often rely on polynomial parameterizations, limiting their effectiveness for general smooth, non-polynomial nonlinearities. To overcome these limitations, we develop a direct optimization framework that employs a Lagrange multiplier formulation to determine the resonance backbone curve on the response surface constrained by the harmonic balance governing equations. Crucially, solving this formulation efficiently requires second-order sensitivity information. Therefore, the primary innovation of this work is the derivation of a analytical Hessian Tensor for generic $C^2$-continuous nonlinear elements. This is achieved by combining an extended Alternating Frequency/Time (AFT) method for computing second-order derivatives with local-coordinate tensor transformations. By integrating this analytical Hessian into the solver, the proposed framework ensures robust convergence and significantly reduces runtime, making it practical for large-scale models where numerical differentiation is computationally prohibitive. The method is validated on three benchmarks of increasing complexity: a two-degree-of-freedom (2-DOF) system with cubic nonlinearity, a beam with cubic stiffness or hyperbolic tangent (tanh) friction nonlinearities, and an industrial-scale finite element model of a compressor bladed disk (blisk) with a friction ring damper. Results demonstrate that the proposed framework accurately and efficiently computes both resonance and anti-resonance backbone curves, providing a robust frequency-domain tool for structures with non-polynomial nonlinearities.

研究の動機と目的

• 複雑な非線形機械系における共振および反共振バックボーン曲線の堅牢な計算を動機づける。
• Harmonic Balance governing equationsを制約とする直接最適化フレームワークを開発する。
• 拡張AFTアプローチを用いて解析Hessianテンソルを導出・組み込み、2次微分を高速かつ正確に行う。
• 多項式ではない非線形性を含む大規模有限要素モデルへの適用性を示す。
• 非滑らかで非多項式の非線形性を持つ系の共振バックボーンを識別するためのスケーラブルで効率的な方法を提供する。

提案手法

• 周波数領域のHB方程式を制約とするバックボーン曲線計算を制約付き最適化問題として定式化する。
• 応答振幅を目的とするラグランジュ関数とラグランジュ乗数を用いて定常条件を得る。
• 非線形要素の解析Hessianテンソルを導出・組み立てるために拡張AFT法を適用する。
• フーリエ-ガレルゴン射影による離散化と、拡張AFTワークフローを通じて二次微分を計算する。
• 局所-全体座標変換による元素HessianをグローバルHessianテンソルへ組み立てる。
• 3つのベンチマーク（2-DOF3乗オシレーター、立方またはtan h 摩擦非線形性を持つカンチレータビーム、摩擦ダンパを有する大規模FEブリンク）でフレームワークを検証する。

実験結果

リサーチクエスチョン

• RQ1 非多項式非線形性を持つ系の共振バックボーン曲線を、Harmonic Balanceに制約された直接最適化フレームワークで正確に追跡できるか。
• RQ2 拡張AFTによる解析Hessianの組み込みは、大規模FEモデルの収束のロバスト性を改善し、実行時間を短縮するか。
• RQ3 低次元と高忠実度の非線形モデル（2-DOF, ビーム, ブリンク）の共振/バックボーン予測は、非滑らか/摩擦非線形下で一貫性があるか。
• RQ4 大規模系の二次情報を効率的に計算・組み立てて、実用的なバックボーン曲線追跡を可能にするにはどうするか。

主な発見

• Hessian強化フレームワークは共振および反共振バックボーン曲線を正確に計算する。
• 拡張AFTによる解析Hessianテンソル導出は解法のロバスト性を向上させ、収束を高速化する。
• 2-DOF、ビーム、産業規模のブリンクを横断して非多項式非線性性への適用性を検証。
• 拡張AFTによる時系列非線性評価は、純粋な周波数領域アプローチと比べて計算コストを削減する。
• 複素値縮約ベースは大規模非線形ブレードディスク系に適したモデル次元削減を可能にする。
• 本手法は変位依存・速度依存の非線形性を扱い、高忠実度FEモデルでも有効である。