摩耗シミュレーション
摩耗の理論基礎
摩耗シミュレーションとは
先生、FEMで摩耗（すり減り）をシミュレーションできるんですか？
できる。Archardの摩耗則に基づいて接触面の材料除去を計算し、形状を更新する。
$V$ は摩耗体積、$s$ はすべり距離、$K$ は摩耗係数、$F_n$ は法線力、$H$ は硬さ。
すべり距離に比例して摩耗する。シンプルなモデルですね。
Archardの式は最もシンプルだが、多くの摩耗問題で実用的な精度が得られる。
FEMでの実装
摩耗シミュレーションの手順（反復的）：
1. 接触解析 — 接触圧とすべり量を計算
2. 摩耗量の計算 — Archardの式で摩耗深さ
3. メッシュの更新 — 摩耗した分だけ節点を移動
4. 再度接触解析 — 更新された形状で
5. 収束するまで反復
メッシュを毎回更新するのは大変そうですね。
AbaqusのUMESHMOTION（適応メッシュ）やLS-DYNAの*MAT_WEARで自動化可能。
まとめ
要点：
- Archardの摩耗則 — $dV/ds = K F_n / H$
- 接触→摩耗量→メッシュ更新の反復 — 摩耗の時間発展
- AbaqusのUMESHMOTION — 適応メッシュで摩耗面を更新
- ベアリング、ギア、ブレーキパッド — 主な適用
Archard摩耗則1953年
摩耗シミュレーションの基礎となるArchard摩耗則は、J.F. ArachardとW. Hirstが1953年に発表した論文に由来する。摩耗量V = k × W × s / H（k:摩耗係数、W:荷重、s:滑り距離、H:硬さ）という単純な比例則だが、60年以上たった現在でも工業的摩耗予測の標準モデルとして使われ続けている。同論文はWear誌の創刊（1957年）に先立ち、摩耗科学（tribology）という学術分野の確立に大きく貢献した。
摩耗の数値計算手法
摩耗のFEM実装
```
*ADAPTIVE MESH, TYPE=CONTACT SURFACE
wear_surface
*ADAPTIVE MESH CONTROLS, MESHING FREQUENCY=1
*UMESHMOTION $ ユーザーサブルーチンで摩耗量を計算
```
```
*MAT_WEAR
K, H $ Archard係数と硬さ
```
LS-DYNAには専用の摩耗材料モデルがあるんですね。
*MAT_WEARは接触面の要素を自動的に薄くして摩耗を表現。リメッシュなしで摩耗が進行。
まとめ
形状更新の数値スキーム
摩耗FEMの数値的な核心は「摩耗量→節点移動→メッシュ更新→再接触」の反復処理にある。Podra & Andersson（1997年）が提案したEulerianアプローチでは、1サイクル分の摩耗量をArchard則で計算し、境界節点を法線方向に後退させた後、スムーシングアルゴリズムで要素歪みを修正する。この手法は現在もABAQUSのUMESHMOTIONユーザーサブルーチンを使った実装の参照実装として参照されている。
摩耗の実務適用
摩耗の実務
実務チェックリスト
切削工具の摩耗予測
金属切削の切削工具寿命予測にArchard摩耗則を組み込んだFEMが2000年代から普及している。サンドビックが2010年代に公開した超硬合金エンドミル解析では、DEFORMコードを用いて切削速度200m/min・切込み1mm条件下での工具フランク摩耗幅VB=0.3mmに達するまでの加工長を、実験値と比較して±18%の誤差で予測した。この結果は工具コーティング最適化の設計判断に活用されている。
摩耗の摩耗のソフトウェア比較
摩耗のツール
選定ガイド
DEFORM摩耗モジュールの歴史
金属加工専用FEMコードDEFORM（Scientific Forming Technologies Corporation）は1989年の設立から切削・鍛造解析に特化してきたが、摩耗モジュールの商用提供はDEFORM V6.1（2000年頃）から始まった。同バージョンでArchard則とUsui則（切削用）の2種類が実装され、工具寿命の予測精度が業界標準となった。現行DEFORM V13では機械学習で摩耗係数を自動補正するアダプティブ機能が追加されている。
摩耗の先端研究
摩耗の先端研究
データ駆動摩耗則
2019年以降、分子動力学シミュレーションとPINNs（Physics-Informed Neural Networks）を組み合わせてナノスケールの摩耗メカニズムをマクロ摩耗則にアップスケーリングする研究が活発化している。IFToMM World Congressの2023年発表では、Fe-Fe接触のMDシミュレーション10万ケースから学習したPINNsが、ArChardモデルが再現できなかった低荷重域の非線形摩耗挙動を予測誤差15%以内で捉えることに成功している。
摩耗のトラブル対応
摩耗のトラブル
摩耗係数のばらつきと失敗
摩耗シミュレーションで最大の不確定要因は摩耗係数kの値だ。同じ鋼材同士でも潤滑状態・面粗度・温度によりkが10⁻⁶から10⁻³まで3桁変動する。2015年に国内自動車部品メーカーが経験した失敗では、文献値のk=5×10⁻⁵をそのまま使ったところ予測摩耗量が実測の7倍になり、設計判断を誤った。根本対策はPin-on-Diskトライボメーター試験で使用条件に対応したkを実測することだ。
なった
詳しく
報告