応力鎖
(force chain から転送)
出典: フリー百科事典『ウィキペディア（Wikipedia）』 (2026/01/09 13:22 UTC 版)
|
この記事には複数の問題があります。
|
|
この記事はカテゴライズされていないか、不十分です。
|
応力鎖（おうりょくさ、英語: force chain）は、粒状体や粉体材料において、粒子間接触を介して伝達される力が空間的に不均一に集中し、鎖状・網状の構造を形成したものを指す。
離散要素法（Discrete Element Method、 DEM）では、接触力と粒子応力テンソルを計算することで、応力鎖ネットワークを定量的に抽出し、材料の強度・変形・流動性などの巨視的特性との関係を評価することができる。
概要
粒状体内部では、外力が一様に分布するのではなく、接触力が一部の粒子列に集中して伝達される。このような高応力の粒子列が「応力鎖」であり、それ以外の粒子群は比較的弱い接触力で全体を支持する「弱いネットワーク」として振る舞うと理解されている。
実験的には、フォトエラスティック円盤などを用いた可視化により、光弾性しまの強度が高い準一次元的な経路として観察され、数値的には DEM による接触力解析から同様の構造が確認されている。応力鎖の長さ分布や空間配置は、せん断帯の形成、局所的な座屈・破壊、流動開始（ジャミング／アンジャミング）などと密接に関連することが報告されている。
定義と理論
粒子応力とマクロ応力
連続体力学における応力テンソル（Cauchy 応力）
Python スクリプトによる変換（liggghts-lammps-dump2vtk）
DEM のダンプを可視化するには、VTK 形式への変換が一般的である。従来、LAMMPS / LIGGGHTS 用の後処理プログラムである Pizza.py / LPP は Python 2 系で実装されており、Python 3 環境では利用しづらいという問題があった。
この制約を解消するため、粉体粒子（ふんたい つぶこ、Open DEM Japan）により、LIGGGHTS / LAMMPS の粒子ダンプ（dump.* や dump.*.gz）を VTK Legacy 形式の .vtk ファイルに一括変換する Python 3 スクリプト liggghts-lammps-dump2vtk が作成・公開されている。このスクリプトは、旧 Pizza.py（lpp.py / vtk.py / dump.py）に相当する機能を Python 3 向けに統合し、
- 単一ファイルのスクリプトとして配布されていること、
- VTK Legacy の ASCII / BINARY 出力や、多数タイムステップのチャンク分割・並列処理に対応すること、
などを特徴とする。出力される .vtk ファイルには、粒子応力や接触数といったスカラー量、粒子速度や力といったベクトル量を点データとして格納でき、ParaView や Mayavi などの可視化ソフトウェアで直接読み込める。
Windows GUI / CLI ツール dump2vtk.exe（D2V）
上記の Python 3 スクリプトをベースにした Windows 用スタンドアロンツールとして、dump2vtk.exe（通称 D2V）が配布されている[1]。開発者は前節と同じく粉体粒子であり、Python 3 版 LPP/Pizza を誰でも簡単に利用できるよう、PySide6 ベースの GUI と同等機能を備えた CLI を備える形で提供している。D2V は LIGGGHTS / LAMMPS の dump を VTK（レガシー形式、ASCII / BINARY）へ高速一括変換するツールであり、Windows 用バイナリ（インストーラ形式）として配布されている。
GUI は 3 つのタブから構成される。
- Particles – 粒子ダンプを VTK に変換するタブ。1 つの dump ファイルをドラッグ＆ドロップすると、同じプレフィクスを共有する連番ファイルを自動収集し、並列処理で一括変換する。
- Header Replace –
ITEM: ENTRIES ...行をバッチ置換するタブ。前節のヘッダー置換スクリプトと同等の機能を GUI から提供する。 - Force chain – pair/local 形式の dump からフォースネットワーク（force chain）を VTK として出力するタブ。Louvain 法によるコミュニティ検出により、力のネットワークをクラスタごとに色分けするテスト的な解析機能も備える。
特に Force chain タブは、GUI 上では「応力」や「応力鎖」の可視化として紹介されることがあるが、実際には pair/local dump に記録された接触ごとの力を基に、強い接触を抽出して線分メッシュとして描画する「力のネットワーク（force chain）」の可視化が主目的である。前節で述べた Love–Weber 型の粒子応力テンソルを通じて Cauchy 応力と対応付けることはできるものの、D2V 自体は連続体力学における応力場全体を厳密に再構成する解析コードではなく、「どの経路に大きな力が流れているか」を可視化するツールとして利用されることが多い。
CLI（CUI）では、GUI の各タブに対応する機能がサブコマンドとして提供されており、バッチ処理や自動実行に利用できる。代表的な用途としては、
- LPP/Pizza 互換の粒子 dump → VTK 変換、
ITEM: ENTRIES ...ヘッダーの一括置換、- pair/local dump からの force chain VTK 出力（Louvain 法によるコミュニティ集約量の出力を含む）、
などが挙げられる。force 系サブコマンドは、前処理済みの pair/local dump（例えば *-renamed.dump）を入力として、タイムステップごとに接触ネットワークの VTK を出力する。これにより、粒子中心を結ぶ線分として応力鎖（力のネットワーク）を可視化したり、接触力の大きさに応じて線の太さや色を変えたりすることが容易になる。
可視化と解析
VTK 形式に変換された粒子・接触データは、ParaView などの可視化ツールで以下のような解析に利用される。
- 粒子を点として描画し、粒子応力や接触数、主応力方向などのスカラー／ベクトル場を色付け・矢印表示する。
- 接触ダンプから生成された線分メッシュを利用し、接触力の大きさに応じた太さ・色で「力のネットワーク」を描画する。
- 一定しきい値以上の接触力だけを表示して「強いネットワーク（strong force network）」を抽出し、応力鎖の空間分布や長さ分布を評価する。
- Louvain 法などのグラフクラスタリングにより、応力鎖ネットワークをコミュニティ別に色分けし、局所構造とマクロ応答の関係を調べる。
このような可視化と統計解析により、接触力・応力鎖・マクロな荷重–変位応答の三者の相関を体系的に調査することが可能となる。
関連項目
脚注
- ^ D2V
外部リンク
- LIGGGHTS 公式サイト（英語） – https://www.cfdem.com/liggghts-open-source-discrete-element-method-particle-simulation-code
- Open DEM Japan – Python による LIGGGHTS ワークフローと dump2vtk 解説
- GitHub – opendemjapan/liggghts-lammps-dump2vtk
- GitHub – opendemjapan/liggghts-entries-header-replacer
- Open DEM Japan – 「Python による LIGGGHTS® ワークフロー」（dump2vtk GUI 配布ページ） – https://opendemjapan.parallel.jp/python.html
- 応力鎖のページへのリンク
