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, tau1TIMEは時間領域、FREQUENCYは周波数領域の材料定義です。Abaqus/Standardでは、異なる領域の解析に指定した粘弾性応答は無視されます。周波数領域では長時間弾性率を定義します。適用する解析手順と弾性率の定義を確認してください。 *VISCOELASTIC · *STEADY STATE DYNAMICS
まとめ
*VISCOELASTIC, TIME=PRONY + *VISCO*VISCOELASTIC, FREQUENCY=PRONY + *STEADY STATE DYNAMICS, DIRECT緩和曲線からのパラメータ同定
PronyパラメータはDMA（動的機械分析）や応力緩和試験のデータから同定します。複数温度の測定を組み合わせる場合、時間温度換算則が成立する範囲でマスターカーブを作り、対象の時間・周波数帯に合う項数と緩和時間を選んでフィットします。項数に一律の最低値はなく、測定値との残差と予測に使う範囲での再現性を確認します。
Prony系列粘弾性の実務適用
実務チェックリスト
自動車防振ゴムへの応用
エンジンマウントゴム（EPDM系）の振動減衰解析ではProny級数粘弾性が標準。周波数10〜1000Hzの損失係数tanδを目標値（典型0.1〜0.3）に合わせてτᵢを調整する。Abaqusの周波数応答解析（*STEADY STATE DYNAMICS）とProny入力を組み合わせ、NVH（騒音・振動・乗り心地）性能を仮想評価するワークフローが日系自動車メーカーで広く採用されている。
Prony系列粘弾性のソフトウェア比較
ツール
ソルバー別Prony入力形式
入力形式はソルバーと材料モデルにより異なります。Abaqusの等方性Prony入力では、各行に無次元のせん断緩和係数 gᵢ、体積緩和係数 kᵢ、緩和時間 τᵢを指定します。gᵢ はせん断緩和係数 Gᵢ/G₀ であり、一般にヤング率 Eᵢ/E₀ をそのまま代用できません。各ソルバーの対象モデルと正規化の定義を確認してください。 Abaqus公式：Viscoelastic object
選定で最も重要な3つの問い
- 「何を解くか」：Prony系列粘弾性モデルに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
Prony系列粘弾性の先端研究
先端
超弾性との結合モデル
ゴムのような大変形粘弾性にはProny+Mooney-Rivlinを組み合わせた「粘超弾性（Viscoelastic Hyperelastic）」が使われる。Abaqusでは*HYPERELASTICと*VISCOELASTICを同一*MATERIALブロックに記述することで自動的に積算定式化が適用される。脳組織の手術シミュレーション（シュレーダー 2011）でも同手法が適用され、変形予測誤差2mm以内が報告されている。
Prony系列粘弾性のトラブル対応
トラブル
時間増分とProny精度
時間増分は、荷重変化と注目する緩和時間を十分に解像できるよう選び、増分を細かくして応力履歴が収束するか比較します。τmin/10は試行の出発点の一例であり、精度や安定性を保証する一律基準ではありません。Abaqusの自動増分では *VISCO のデータ行の第4項に最大時間増分を指定できます。 Abaqus公式：*VISCO
関連トピック
なった
詳しく
報告