粘弾性減衰材の周波数応答
粘弾性減衰材の周波数応答の理論基礎
粘弾性材料とは
先生、粘弾性材料の減衰って構造減衰とどう違いますか？
構造減衰は $g$ が周波数に依存しないと仮定するが、粘弾性材料の減衰は周波数と温度に強く依存する。ゴム、ポリマー、接着剤、制振材が代表例。
複素弾性率
粘弾性材料の応力-ひずみ関係：
- $E'(\omega)$ — 貯蔵弾性率（剛性。エネルギーを蓄える）
- $E''(\omega)$ — 損失弾性率（減衰。エネルギーを散逸する）
- $\eta(\omega) = E''/E'$ — 損失係数
周波数によって剛性も減衰も変わるんですね。
典型的な挙動：
- 低周波 — $E'$ が低い（柔らかい）、$\eta$ が中程度
- ガラス転移領域 — $E'$ が急増、$\eta$ がピーク（最大減衰）
- 高周波 — $E'$ が高い（硬い）、$\eta$ が低下
温度-周波数等価則（WLF則）
粘弾性材料は温度が上がると低周波側にシフトし、温度が下がると高周波側にシフトする。Williams-Landel-Ferry（WLF）式でこの等価性を記述：
$a_T$ がシフトファクター。これにより1つの温度のデータから全温度の特性を推定できる。
DMA（動的粘弾性試験）で一つの温度のデータを測定すれば、WLF式で他の温度に外挿できるんですね。
マスターカーブを作成する。基準温度でのDMAデータをWLF則でシフトし、広い周波数範囲の $E'(\omega), \eta(\omega)$ を得る。これをFEMに入力する。
まとめ
要点：
- $E^*(\omega) = E'(1+i\eta)$ — 周波数依存の複素弾性率
- $E', \eta$ とも周波数と温度に依存 — ガラス転移でピーク
- WLF式で温度-周波数等価 — マスターカーブの構築
- 直接法（周波数応答）でFEMに入力 — モード法では扱いにくい
粘弾性の記憶効果はMaxwellが1867年に発見
粘弾性材料の「ひずみが応力の歴史に依存する」という記憶効果は、ジェームズ・クラーク・マクスウェルが1867年にガス粘性の研究中に発見した。後にKelvin-Voigt（1890年）が並列モデルを提案し、Maxwell直列+Kelvin並列を組み合わせた「Zener（標準線形固体）モデル」が構造減衰のFEM実装の基礎となった。3Mのダンピングシートの特性評価にも今日このモデルが使われている。
粘弾性減衰材の周波数応答の数値計算手法
FEMでの粘弾性材料
粘弾性材料をFEMでどう扱いますか？
Prony系列で $E'(\omega), E''(\omega)$ をモデル化：
$E_i, \tau_i$ がProny系列のパラメータ。$N = 5 \sim 15$ 項で広い周波数範囲をカバー。
Abaqus
```
*VISCOELASTIC, FREQUENCY=PRONY
g_1, k_1, tau_1
g_2, k_2, tau_2
...
```
Nastran
```
MAT1, 1, ...
TABLEM1, 100, ...
$ 周波数依存のE'とηをテーブルで定義
```
Nastranでは直接テーブル（周波数 vs. E', η）で入力可能。
AbaqusのProny系列はDMAデータから直接フィッティングできますか？
Abaqus/CAEの*VISCOELASTIC, TEST DATAオプションで実測の$E'(\omega), E''(\omega)$を入力し、Prony系列を自動フィッティングできる。非常に便利。
制振材（CLD）の解析
CLD（Constrained Layer Damping）の解析：
1. 基板（鋼板等）＋粘弾性層＋拘束層の積層をモデル化
2. 粘弾性層に周波数依存のProny系列を設定
3. 直接法の周波数応答解析で減衰効果を評価
4. FRFのピーク低減を確認
まとめ
時間-温度換算則でデータを1/100に圧縮
粘弾性材料の損失係数は温度と周波数が複合して変化する。Williams-Landel-Ferry（WLF）式（1955年）による時間-温度換算則を使えば、基準温度Trでの測定データだけで任意温度・周波数の特性を予測できる。車内の内装材（エチレンプロピレンゴム系）の周波数特性は-40℃〜100℃・1Hz〜10kHzの全ての条件を20℃基準のマスターカーブ1本で表現可能。
粘弾性減衰材の周波数応答の実務適用
粘弾性減衰の実務
制振材の設計は自動車のNVHで最も重要な適用。
制振材の設計フロー
1. DMA試験で粘弾性材料の$E'(\omega, T), \eta(\omega, T)$を測定
2. WLF則でマスターカーブを構築
3. Prony系列にフィッティング
4. FEMで制振材の配置と厚さを最適化
5. FRFの共振ピーク低減を確認
実務チェックリスト
制振材は「最大ひずみの場所」に貼るんですか。
制振材はせん断変形でエネルギーを散逸する。基板の曲げひずみが最大の場所（固定端付近）に制振材を貼ると効果が最大。ひずみの小さい場所（自由端付近）に貼っても効果がない。
3Mの制振シートで新幹線騒音を5dB低減
3M社の制振材「ダイナマット」シリーズはスチール板に粘弾性ポリマー（アクリル系）を積層した構造で、100〜1000Hzの曲げ波を熱エネルギーに変換する。JR東日本のE5系新幹線（2011年）の床下鉄板に3Mダイナマット相当品を貼付し、軌道騒音由来の車内固体音を約5dB低減した。厚さわずか2mmで面密度2.4kg/m²。コスト優位性から量産車両に広く普及している。
粘弾性減衰材の周波数応答のソフトウェア比較
粘弾性解析のツール
選定ガイド
日東電工とESIが連携する粘弾性材料DB
日東電工（ニチバン・テープ系最大手）はESI Group（VA One）と提携し、主要制振材（レジェトレックスシリーズ等）の動的特性データベースをVA Oneに組み込んだ。これにより設計者は材料選定と解析を同一ソフトで完結できる。3Mも自社制振材のAbaqus Material Cardを公式配布しており、材料データの入力工数を大幅に削減できる。2022年時点で主要制振材30種以上がCAEデータベース化されている。
粘弾性減衰材の周波数応答の先端研究
分数微分粘弾性モデル
Prony系列の代替として分数微分モデルが研究されている。少ないパラメータ（3〜4個）で広い周波数範囲をカバーできる。FEMへの実装は発展途上。
トポロジー最適化で制振材配置
まとめ
拘束型制振はフリー型より損失係数が5倍高い
制振材の構成は①フリー層（単体貼り付け）と②拘束層（金属板でサンドイッチ）に大別され、拘束型はせん断ひずみを集中させるため損失係数が5〜10倍高い。Nashif（1985年）の実験データでは1mm厚の拘束型制振処理が同重量のフリー型より3倍の振動低減効果を示した。自動車ダッシュパネルには拘束型が主流で、設計には3Mや日東電工のCAEデータベースをAbaqus/Ansysに直接インポートできる。
粘弾性減衰材の周波数応答のトラブル対応
Prony系列のフィッティングが悪い
温度依存性を無視した
粘弾性は温度で劇的に変化。ガラス転移温度付近で$E'$が10倍、$\eta$がピークに。使用温度全範囲で評価すべき。
モード法で粘弾性を使ってしまった
モード法は固有モード（材料特性に依存しない基底）で展開するため、周波数依存の材料を正確に扱えない。直接法に切り替える。
まとめ
夏と冬で減衰効果が3倍変わる
粘弾性材料は温度依存性が強く、同じ制振シートでも0℃（冬季）と40℃（夏季）では損失係数が2〜3倍変動する。北海道向けの鉄道車両で夏季の試験を通過した制振対策が冬季走行で効果を失った事例が2008年に報告されている。対策は使用温度範囲全体でマスターカーブを取得し、最悪温度条件で設計余裕を確認すること。ANSYSのViscoelastic Material Modelに温度テーブルを設定するのが標準手順。
関連トピック
なった
詳しく
報告