[論文レビュー] Extended Hamilton's principle

Jin-Kyu Kim|arXiv (Cornell University)|Apr 2, 2012

Numerical methods for differential equations参考文献 15被引用数 5

ひとこと要約

この論文は、混合変分法とレイリーの散逸関数を用いて、保存的および非保存的力学系の両方のための整合性のある初期条件を含めるために、ハミルトンの原理を拡張する。この手法により、運動方程式を導出する統一的なフレームワークが可能となり、集中パrameter系のシミュレーションにおいて安定性と精度が実証されたガラーキンに基づく数値法が開発される。

ABSTRACT

Hamilton's principle is extended to have compatible initial conditions to the strong form. To use a number of computational and theoretical benefits for dynamical systems, the mixed variational formulation is preferred in the systems other than particle systems. With this formulation and the Rayleigh's dissipation function, we could have all the pertinent initial/boundary conditions for both conservative and non-conservative dynamical system. Based upon the extension framework of Hamilton's principle, the numerical method for representative lumped parameter models is also developed through applying Galerkin's method to time domain with the discussion of its numerical properties and simulation results.

研究の動機と目的

• 保存的および非保存的系の両方のための整合性のある初期条件を満たすように、ハミルトンの原理を拡張すること。
• 混合変分法を用いて、力学系における初期条件と境界条件の取り扱いを統一すること。
• レイリーの散逸関数を変分枠組み内に組み込み、非保存的力（例：粘性ダンピング）をモデル化すること。
• 集中パrameter系の時間領域シミュレーションに適した安定的かつ高精度なガラーキンに基づく数値法を開発すること。
• 提案手法の数値的性質および代表的なモデルにおけるシミュレーション性能を実証すること。

提案手法

• 追加の状態変数を導入することで、整合性のある初期条件を実現する混合変分原理を定式化する。
• 非保存的力（例：粘性ダンピング）をモデル化するため、レイリーの散逸関数を変分枠組みに統合する。
• 時間領域におけるガラーキン法を用いて、拡張された変分原理を離散化し、常微分方程式系を導出する。
• 重み付き残差法を用いて、有限要素に類似した空間的および時間的離散化に適した弱形式を導出する。
• 変分構造を保存するように、導出されたシステム行列および時間積分スキームを導出する。
• ベンチマークとなる集中パrameterモデルを用いた数値的シミュレーションにより、手法の妥当性を検証し、安定性および収束性を分析する。

実験結果

リサーチクエスチョン

• RQ1 ハミルトンの原理は、運動方程式の強形式と初期条件の両者を整合させるために、どのように拡張可能か？
• RQ2 保存的および非保存的力学系の両方を一貫して取り扱える変分原理は何か？
• RQ3 レイリーの散逸関数は、初期条件を含む変分枠組みにどのように自然に組み込むことができるか？
• RQ4 拡張された変分原理から導出されたガラーキンに基づく時間離散化法の数値的性質は何か？
• RQ5 標準的手法と比較して、提案手法は代表的な集中パrameter系のシミュレーションにおいて、どのように性能を発揮するか？

主な発見

• 拡張された変分原理は、初期条件と運動方程式の強形式との間の整合性を的確に確保した。
• 変分枠組み内にレイリーの散逸関数を組み込むことで、非保存的力の一貫したモデル化が可能となった。
• ガラーキンに基づく時間離散化法は、線形および非線形の集中パrameter系のシミュレーションにおいて、良好な数値的安定性と精度を示した。
• 変分構造が保存されることで、減衰を伴う物理的系と整合したエネルギー減衰挙動が得られた。
• シミュレーション結果は、複雑な初期条件を有する系を含む、さまざまなテストケースにおいて、手法の有効性と頑健性を確認した。
• 本手法は、単一の変分原理から運動方程式と数値スキームを統一的に導出するフレームワークを提供する。