粘弾性体の力学的挙動を理解・予測するためには、適切な力学モデルの知識が不可欠です。
「ばね」と「ダッシュポット」を組み合わせることで、実際の材料が示す複雑な応力緩和やクリープ挙動を数学的に表現することができます。
本記事では、粘弾性モデルの代表格であるMaxwellモデル・Voigtモデルの原理・構成方程式・応答特性から、これらを組み合わせた高度なモデルまで、体系的にわかりやすく解説します。
材料力学・高分子物理・レオロジーを学ぶすべての方にとって、基礎から応用まで役立つ内容となっています。
粘弾性モデルの基本とは？Maxwell・Voigtモデルの結論
それではまず、粘弾性モデルの基本的な考え方と、MaxwellモデルおよびVoigtモデルの本質的な違いについて解説していきます。
粘弾性モデルとは、理想的な弾性体（ばね）と理想的な粘性体（ダッシュポット）を組み合わせることで、実際の材料の力学応答を模擬する力学モデルのことです。
Maxwellモデルはばねとダッシュポットを「直列」に接続したモデルで、応力緩和の表現に優れています。一方、Voigtモデル（Kelvin-Voigtモデルとも）は「並列」接続であり、クリープ挙動の表現に適しています。
この直列・並列という接続方法の違いが、モデルが示す力学的応答の性格を根本的に変えます。
どちらのモデルも現実の材料挙動を完全に再現するものではありませんが、現象の本質的な理解と基礎的な解析には非常に有用です。
ばね要素の力学的意味
ばね要素は理想的な弾性体を表します。
フックの法則に従い、応力σと変形ひずみεは比例関係にあります。
ばね（弾性体）の構成方程式：σ ＝ E・ε
（E：ヤング率（弾性係数）、σ：応力、ε：ひずみ）
特徴：力を加えた瞬間に変形し、力を除いた瞬間に完全に回復する
ばね要素では、エネルギーが完全に保存されるため損失は一切生じません。
これは、粘弾性モデルにおける弾性成分（貯蔵弾性率E’）に対応するものです。
ダッシュポット要素の力学的意味
ダッシュポットはシリンダーとピストンで構成される機械要素で、理想的な粘性体（ニュートン流体）を表します。
ダッシュポット（粘性体）の構成方程式：σ ＝ η・(dε/dt)
（η：粘度、dε/dt：ひずみ速度）
特徴：力を加えると一定速度で変形し続け、力を除いても回復しない
ダッシュポット要素では与えられたエネルギーはすべて熱として散逸します。
これは粘弾性モデルにおける粘性成分（損失弾性率E”）に対応しています。
Maxwellモデルの構成方程式と応力緩和
Maxwellモデルは、ばねとダッシュポットを直列接続したモデルです。
直列接続では、各要素にかかる応力は等しく、全体のひずみは各要素のひずみの和になります。
Maxwellモデルの構成方程式：dε/dt ＝ (1/E)・(dσ/dt) ＋ σ/η
一定ひずみ（ε＝const）を保った場合の応力緩和：
σ(t) ＝ σ₀・exp(－t/τ)
（τ＝η/E：緩和時間）
Maxwellモデルは指数関数的に応力が緩和する挙動を表現できる点が大きな特徴です。
一方で、一定応力（クリープ）条件では変形が際限なく増大し続けるという欠点があり、固体的な材料のクリープ挙動には適さないとされています。
Voigtモデルの特徴とクリープ応答の解析
続いては、Voigtモデルの詳細とクリープ応答への適用を確認していきます。
Voigtモデル（Kelvin-Voigtモデル）は、ばねとダッシュポットを並列に接続したモデルであり、固体的な材料のクリープ挙動の表現に優れています。
Voigtモデルの構成方程式とクリープ応答
並列接続では、各要素のひずみは等しく、全体の応力は各要素の応力の和になります。
Voigtモデルの構成方程式：σ ＝ E・ε ＋ η・(dε/dt)
一定応力（σ＝const）を保った場合のクリープ応答：
ε(t) ＝ (σ₀/E)・
（τ＝η/E：遅延時間）
Voigtモデルではひずみが指数関数的に増加し、最終的にはばねの変形量（σ₀/E）に収束します。
この「有限の変形量への収束」という挙動は、多くの固体材料のクリープに見られる実際の挙動と定性的に一致しています。
ただし、Voigtモデルでは応力をゼロにした瞬間に変形が回復するため、永久変形（塑性変形）は表現できません。
MaxwellモデルとVoigtモデルの比較
| 特性 | Maxwellモデル（直列） | Voigtモデル（並列） |
|---|---|---|
| 接続方法 | 直列 | 並列 |
| 応力緩和 | 指数関数的に緩和（良好） | 一定応力を維持（不適） |
| クリープ応答 | 無限に変形し続ける（不適） | 有限量に収束（良好） |
| 適した材料 | 液体的粘弾性体（高分子溶液） | 固体的粘弾性体（架橋ゴムなど） |
| 永久変形の表現 | 可能 | 不可能 |
このように、どちらのモデルも一方の現象（応力緩和またはクリープ）の表現には優れていますが、他方には限界があります。
より高精度な材料挙動のシミュレーションには、これらを組み合わせた複合モデルが必要になります。
Voigtモデルの動的応答とDMAとの対応
Voigtモデルに正弦波変位を入力した場合の動的応答を考えると、以下の関係が得られます。
ε(t) ＝ ε₀ sin(ωt) として
σ(t) ＝ E・ε₀ sin(ωt) ＋ η・ω・ε₀ cos(ωt)
したがって：貯蔵弾性率 E’ ＝ E、損失弾性率 E” ＝ η・ω
tanδ ＝ η・ω／E
この結果は、Voigtモデルにおいてはtanδが周波数ωに比例することを示しており、周波数が高いほど粘性的に振る舞うことを意味します。
実際の高分子材料でもこれに近い傾向が観察されることがあり、モデルの妥当性の一端を確認できます。
高度な粘弾性モデル：標準線形固体と一般化Maxwell・Voigtモデル
続いては、実際の材料挙動をより精度よく表現するための高度な粘弾性モデルを確認していきます。
現実の材料は単一の緩和時間だけでは表現しきれない複雑な挙動を示すことが多く、より精緻なモデルが必要です。
標準線形固体（Zenerモデル）
標準線形固体（SLS：Standard Linear Solid）はZenerモデルとも呼ばれ、Maxwellモデルとばね要素を並列に接続した3要素モデルです。
このモデルは、応力緩和とクリープの両方をある程度適切に表現できるため、単純な2要素モデルよりも広い応用範囲を持ちます。
ゴム材料や生体組織など、固体的な長期安定性と緩和現象の両方を示す材料のモデル化に適しています。
構成方程式はMaxwellモデルやVoigtモデルより複雑になりますが、得られる応答の現実適合性は大幅に向上します。
一般化Maxwellモデル（Pronyシリーズ）
実際の高分子材料には、単一の緩和時間ではなく広い分布を持つ緩和時間スペクトルが存在します。
この現実を表現するために用いられるのが、複数のMaxwell要素を並列に接続した一般化Maxwellモデル（Prony級数表現）です。
一般化Maxwellモデルの緩和弾性率：
E(t) ＝ E∞ ＋ Σᵢ Eᵢ・exp(－t/τᵢ)
（E∞：平衡弾性率、Eᵢ：各Maxwell要素の弾性率、τᵢ：各緩和時間）
このモデルは有限要素解析（FEA）ソフトウェアでも広く採用されており、DMA測定データからPronyパラメータを同定することで、実際の製品解析に応用されています。
一般化Voigtモデル（遅延時間スペクトル）
同様に、複数のVoigt要素を直列に接続した一般化Voigtモデルも存在します。
こちらはコンプライアンス（ひずみを応力で除した値）の緩和スペクトルを表現するのに適しており、クリープコンプライアンスのフィッティングに使用されます。
一般化MaxwellモデルとVoigtモデルは数学的に等価な表現であり、緩和弾性率とクリープコンプライアンスの間にはLaplaceインバース変換によって相互変換が可能です。
どちらの表現を用いるかは、解析の目的や扱いやすさによって選択するとよいでしょう。
粘弾性モデルの実際の応用と限界
続いては、粘弾性モデルが実務・研究においてどのように活用されているか、また現実の材料挙動との乖離（限界）についても確認していきます。
粘弾性モデルは理論的な枠組みを提供しますが、あくまで現実の材料を近似するものであり、適用範囲を正しく理解した上で活用することが重要です。
有限要素解析（FEA）への組み込み
自動車部品・電子機器・航空宇宙機器の設計では、有限要素解析（FEA）を用いた構造解析が標準的な手法です。
粘弾性材料を含む部品の解析では、一般化MaxwellモデルのPronyパラメータをFEAソフトウェアに入力することで、応力緩和・クリープ・動的応答を高精度にシミュレーションできます。
DMA測定データからPronyパラメータを最適化・同定するプロセスは、実験と計算を橋渡しする重要な工程となります。
ABAQUS・ANSYS・COMSOL MultiphysicsなどのFEAソフトウェアでは、粘弾性材料モデルが標準で組み込まれており、実務設計への適用が容易になっています。
生体材料と非線形粘弾性
骨・軟骨・腱・皮膚・血管壁などの生体組織は、粘弾性挙動を示しますが、変形量が大きい領域では線形粘弾性の仮定が成立しない非線形粘弾性を示します。
この場合、MaxwellモデルやVoigtモデルといった線形粘弾性モデルでは対応できず、QLV（準線形粘弾性）モデルやFung弾性モデルなど、非線形性を取り込んだより複雑なモデルが必要になります。
バイオメカニクス分野では、これらの高度なモデルを用いた組織挙動の解析が活発に研究されています。
粘弾性モデルの限界と将来展望
ばねとダッシュポットを組み合わせた古典的粘弾性モデルには、いくつかの根本的な限界があります。
まず、温度依存性を直接組み込む枠組みがないため、温度変化を伴う問題では追加の熱－力学連成モデルが必要です。
また、塑性変形・損傷・経時変化（エイジング）などの不可逆現象は、線形粘弾性モデルでは原理的に表現できません。
近年では、分数階微積分を用いたフラクショナル粘弾性モデルや、機械学習を活用したデータ駆動型の材料モデルなど、古典的モデルの限界を超えた新しいアプローチも盛んに研究されています。
材料モデリングの世界は現在も急速に進化を続けており、今後のさらなる発展が期待される分野といえるでしょう。
まとめ
本記事では、粘弾性モデルの基本から応用まで、Maxwellモデル・Voigtモデルを中心に体系的に解説しました。
ばね（弾性体）とダッシュポット（粘性体）の組み合わせ方によって、材料の応力緩和やクリープ挙動を数学的に表現できることが、粘弾性力学モデルの本質です。
MaxwellモデルとVoigtモデルはそれぞれ長所と短所を持ち、より精密な解析には標準線形固体や一般化Maxwell/Voigtモデル（Prony級数）が用いられます。
これらのモデルは、有限要素解析・生体材料研究・制振設計など幅広い実用分野で活躍しており、材料科学・機械工学・バイオメカニクスを横断する基礎的な知識体系です。
粘弾性モデルの理解を深めることで、材料評価・設計・解析のすべての段階でより的確な判断が可能になるでしょう。
