Prony系列粘弾性モデル
Prony系列粘弾性の理論基礎
Prony系列粘弾性
先生、Prony系列は粘弾性の標準モデルですよね。
Prony系列は一般化Maxwell模型の離散化。緩和弾性率を指数関数の和で表現：
$G_i$ は各Maxwell要素のせん断弾性率、$\tau_i$ は緩和時間。ゴム、ポリマー、接着剤の時間依存挙動を記述。
まとめ
Prony級数の歴史的起源
Prony（プロニー）級数の名はフランスの数学者Gaspard de Prony（1755〜1839）に由来するが、粘弾性のリラクゼーション関数への適用は20世紀のレオロジー分野で確立した。E(t)=E∞+Σ Eᵢ exp(-t/τᵢ)という形で、Maxwell要素の並列モデル（Generalized Maxwell）と数学的に等価であり、任意の線形粘弾性挙動を近似できる。
Prony系列粘弾性の数値計算手法
Prony系列のFEM設定
```
*VISCOELASTIC, TIME=PRONY
g1, k1, tau1
g2, k2, tau2
```
または:
```
*VISCOELASTIC, FREQUENCY=PRONY
g1, k1, tau1
```
TIME=時間領域、FREQUENCY=周波数領域。
まとめ
緩和曲線からのパラメータ同定
Pronyパラメータの同定にはDMA（動的機械分析）または応力緩和試験を用いる。周波数掃引で得たE'(ω)・E''(ω)を時間-温度換算則（WLF式：log aT=−C₁(T−Tref)/(C₂+T−Tref)）でマスターカーブに変換し、Pronyseries.pyなどのツールで離散Prony係数を最小二乗フィットする。最低8〜12τᵢ点が精度確保に必要だ。
Prony系列粘弾性の実務適用
実務チェックリスト
自動車防振ゴムへの応用
エンジンマウントゴム（EPDM系）の振動減衰解析ではProny級数粘弾性が標準。周波数10〜1000Hzの損失係数tanδを目標値（典型0.1〜0.3）に合わせてτᵢを調整する。Abaqusの周波数応答解析（*STEADY STATE DYNAMICS）とProny入力を組み合わせ、NVH（騒音・振動・乗り心地）性能を仮想評価するワークフローが日系自動車メーカーで広く採用されている。
Prony系列粘弾性のソフトウェア比較
ツール
ソルバー別Prony入力形式
Prony級数の入力形式はソルバーにより異なる。Abaqusは正規化相対弾性率gᵢ=Eᵢ/E₀と緩和時間τᵢのペアで入力、LS-DYNAのMAT_076はE∞と各Eᵢを絶対値で入力、MSC Marcは弾性率比と時定数、ANSYSはShift function込みのProny係数テーブルを用いる。相互変換ミスが異なるソルバー間での解析比較で頻発するため、正規化の有無を必ず確認すること。
選定で最も重要な3つの問い
- 「何を解くか」：Prony系列粘弾性モデルに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
Prony系列粘弾性の先端研究
先端
超弾性との結合モデル
ゴムのような大変形粘弾性にはProny+Mooney-Rivlinを組み合わせた「粘超弾性（Viscoelastic Hyperelastic）」が使われる。Abaqusでは*HYPERELASTICと*VISCOELASTICを同一*MATERIALブロックに記述することで自動的に積算定式化が適用される。脳組織の手術シミュレーション（シュレーダー 2011）でも同手法が適用され、変形予測誤差2mm以内が報告されている。
Prony系列粘弾性のトラブル対応
トラブル
時間増分とProny精度
Prony粘弾性の時間積分では、最小緩和時間τminに対して時間増分Δt<τmin/10が目安。Δtが大きすぎると短い緩和成分を積分しきれず貯蔵弾性率E'が過大評価される。Abaqusでは*VISCOステップで最大増分数とDTMAX（最大Δt）を設定することが重要で、設定忘れで発散や不正確な応力緩和結果になるトラブルが頻発する。
関連トピック
なった
詳しく
報告