バネ・ダッシュポット系 戻る
粘弾性力学

バネ・ダッシュポット系シミュレーター

Maxwell・Kelvin-Voigt・SLSモデルを選択し、クリープ・応力緩和・動的粘弾性(E'、E''、tanδ)をリアルタイムで計算・可視化。モデル模式図付きで粘弾性挙動を直感的に理解できます。

モデル・パラメータ

モデル選択
E₁ (MPa) 100
η (MPa·s) 500
E₂ (MPa) [SLS] 50
入力値 (σ₀ or ε₀) 10
緩和時間 τ (s)
瞬間弾性率 E₀
長期弾性率 E∞
ピーク tanδ
Maxwell 応力緩和
$$\sigma(t) = \sigma_0 e^{-t/\tau}, \quad \tau = \frac{\eta}{E_1}$$ KV クリープ
$$\varepsilon(t) = \frac{\sigma_0}{E_1}\left(1 - e^{-t/\tau}\right)$$

クリープ(青:ひずみ ε)または応力緩和(赤:応力 σ)の時間応答。

選択中のモデルの模式図。バネ(弾性要素)とダッシュポット(粘性要素)の接続方法。

E'(貯蔵弾性率)、E''(損失弾性率)、tanδ の角周波数ω依存性。

理解を深める会話

🧑‍🎓
「粘弾性」って結局どういう性質ですか?弾性と粘性が「混ざった」というのがよくイメージできなくて…
🎓
例で考えよう。「純粋な弾性体」はバネ——力を加えると即変形して、力を抜くと即戻る。「純粋な粘性体」は水やオイル——ゆっくり押すと抵抗が小さく、速く押すと抵抗が大きい。「粘弾性体」はその中間。たとえばシリコーンゲルは、ゆっくり押すとズルズル変形するのに、ボールにすると弾む。速さ(時間)によって挙動が変わる材料がまさに粘弾性だ。
🧑‍🎓
Maxwellモデルを「クリープ」タブで試したら変形が時間とともにどんどん増えていきますね。これって現実の材料みたいに「ある程度で止まる」挙動にならないんですか?
🎓
そうなんだ。Maxwellモデルは直列接続だから、ダッシュポット部分は際限なく変形する(液体のように流れる)。これは液体状の樹脂には合うけど、固体状の材料には合わない欠点だ。KV(Kelvin-Voigt)モデルに切り替えると、クリープが最終値で頭打ちになるのが分かる。SLSモデルは「KVの欠点(瞬時応答なし)」と「Maxwellの欠点(無限クリープ)」を両方改善した中間的なモデルなんだ。
🧑‍🎓
「動的粘弾性」タブのE'とE''って何を意味するんですか?角周波数を上げるとE'が上がるのはなぜ?
🎓
E'は弾性的な成分(エネルギーを蓄える)、E''は粘性的な成分(エネルギーを熱に変える)だ。周波数が低い(ゆっくり変形)と分子鎖が動く時間があって柔らかい。周波数が高い(速く変形)と分子鎖が追いつけず、バネのように固くなる——だからE'が上がる。ゴムハンマーで叩くと硬いのに、ゆっくり押すと柔らかいのはこれだよ。
🧑‍🎓
tanδ(タンジェントデルタ)のピークが特定の周波数にありますね。これは何を意味しますか?
🎓
tanδ = E''/E' は損失係数で、「振動エネルギーをどれだけ熱に変えるか」の指標だ。ピーク位置(ω=1/τ)が材料の「ガラス転移領域」に相当する。車のエンジンマウントやタイヤは、走行振動域(10〜100 Hz)でtanδが高い材料を選ぶことで振動を吸収する。逆に精密機器の筐体では変形させたくないからtanδが低い材料を使う。
🧑‍🎓
CAEの有限要素解析でこういう粘弾性を使うにはどうすればいいんですか?
🎓
AnsysやAbaqusでは「粘弾性材料モデル」としてProny級数展開で入力する。Prony級数は一般化Maxwellモデル(複数のMaxwell要素の並列)で、n個の緩和時間τᵢと弾性率gᵢを指定する。DMA(動的粘弾性測定装置)のデータから最小二乗フィッティングでこれらのパラメータを求め、Ansys Materialsに入力するのが実務のフローだよ。このシミュレーターのグラフが、そのフィッティングの「感覚」を掴む練習になる。

よくある質問

クリープと応力緩和は何が違いますか?
クリープは「一定応力のもとで変形が時間とともに増える」現象です(例:プラスチック棚に重荷をずっと載せる)。応力緩和は「一定変形のもとで応力が時間とともに減る」現象です(例:ボルトで締めたゴムガスケットの締め付け力が徐々に低下する)。シミュレーターで「応答タイプ」を切り替えてグラフの違いを確認してください。
緩和時間τが材料設計でなぜ重要なのですか?
τは「その材料の特徴的な時間スケール」を定めます。製品の使用時間がτより十分長い場合、材料はほぼ液体的に振る舞います(長期クリープ問題)。使用時間がτより十分短い場合は弾性的に振る舞います。地震(秒オーダー)に対する免震ゴムの設計では、τが地震継続時間と同程度になるよう制振材を設計することがあります。
実材料のτはどのくらいですか?
室温での典型値:シリコーンゴム(τ≈0.1〜10秒)、エポキシ樹脂(τ≈数時間〜数日)、粘着テープ(τ≈数秒〜数分)、コンクリート(τ≈数十年)。温度が上がるとτは大幅に短くなります(アレニウス則)。これが「高温環境でのクリープ問題」が工学的に重要な理由です。
損失係数tanδが高い材料の例は何ですか?
粘着テープの粘着剤(tanδ≈1〜5)、防振ゴム(tanδ≈0.3〜1)、アスファルト(温度依存:高温でtanδ大)、タイヤのゴム(走行域でtanδ=0.1〜0.3)が典型例です。一般に tanδが大きい材料は振動吸収に優れますが、ヒステリシス発熱が多く、タイヤでは転がり抵抗の増加(燃費悪化)につながります。
時温換算則(WLF式)とは何ですか?
高分子材料では「温度を上げることと、測定周波数を下げることは等価」という関係(時温換算則)が成り立ちます。低温・高周波での挙動は高温・低周波での挙動と等価であり、これを利用して広い周波数範囲のデータを少ない実験で取得できます(マスターカーブ)。CAEでは温度依存の粘弾性挙動を予測するために活用されます。
SLSモデルをProny級数で表すとどうなりますか?
SLSの緩和弾性率はG(t) = E₂ + E₁·exp(-t/τ)と書けます。これはProny級数の1項版G(t) = G∞ + Σ Gᵢ·exp(-t/τᵢ)に対応します。AnsysやAbaqusへ入力する際は、ΔG₁=E₁、τ₁=η/E₁、G∞=E₂に換算します。実材料では多項(n=5〜20項)のProny級数を使い、DMAデータに最小二乗フィッティングします。