ボンデッドパーティクルモデル

Bonded particle model（ボンデッドパーティクルモデル、bonded particle model、BPM）は、粒子間にばねとダッシュポットからなる結合を導入することにより、岩石・脆性材料などの破壊挙動を表現する離散要素法の一種である。

ここでは、2 粒子間ボンドの力学的モデル化について述べる。

 基本構成 


 結合断面積と断面二次モーメント 

2つの粒子  と  の結合点に対して，接合内側半径 ，接合外側半径  を粒子  と粒子  の粒子半径  と  を用いて次のように定義する。

結合断面積  は

極断面二次モーメント  は

極断面二次モーメント（断面全体に対してねじりに抵抗する抵抗）は

主断面二次モーメント（断面をある軸まわりに曲げモーメントがかかったときの抵抗）は

 ボンド長と各種剛性定数 

ボンド長さ  は

法線バネ定数  は

接線バネ定数  は

ねじりバネ定数  は

曲げバネ定数  は

ここで , , , ,  は定数である。Potyondy は ,  とし， と  の関係はラメ定数で決めた。また， と仮定した。接合ヤング率は

と定義した。Horabik は ， で条件を振って脆性・半脆性・延性を表現できると報告している。

 有効質量と慣性モーメント 

各粒子の質量は

有効質量  は

各粒子の極慣性モーメントは

有効極慣性モーメント  は

 相対運動の分解 

粒子間ベクトルは

その長さは

法線方向の単位ベクトルは

 直線運動による相対速度 

相対速度ベクトルは

法線方向成分は

接線方向成分は

 回転運動による相対速度の補正 

回転によるせん断相対速度の補正角速度は

接線補正後の相対速度は

 曲げ・ねじり回転の相対速度 

曲げとねじりの相対角速度は

ねじり成分は

曲げ成分は

 力とモーメント 

 法線方向ばね力 

伸縮量（ひずみ）は

法線ばね力は

指数補正を加味する場合は

ここで

である。

 接線方向ばね力（せん断） 

 ねじりトルク・曲げトルク 

ねじりトルクは

曲げトルクは

 ダンピング項 

法線ダンピング力は

接線ダンピング力は

ねじりダンピングトルクは

曲げダンピングトルクは

 履歴ベクトルの回転（ベクトル表記） 

粒子間の結合軸  が毎ステップ傾くたび，前ステップの履歴  を新しい軸に射影し直す操作を「回転」と呼ぶ。ここで「'」や「''」のついたベクトルは計算用の一時的なベクトル変数である。

 定義 

履歴ベクトルを

とする（上付き  は時刻  を示す）。

 増分加算（オイラー陽） 

力とトルクの更新にはバネだけ考え，ダッシュポットは考えない。

これは

と書き換えられる。

 法線力の回転（未使用） 

法線力については毎ステップ重なり量から新しく計算し直すだけなので，過去の大きさを温存する必要がない。

 せん断力の回転 

全体  から法線成分を引くと接線成分になる。これは直交化（Gram–Schmidt 法）と呼ばれる。履歴ばねとして「前ステップまでに蓄えられたせん断バネ伸び（＝エネルギー）」をそっくり次ステップへ持ち越す必要がある。

 ねじりモーメントの回転 

ここでは  ではなく  を使う。 を使ってしまうと符号が逆になる可能性があるためである。履歴ばねとして「前ステップまでに蓄えられたねじりバネ伸び（＝エネルギー）」をそっくり次ステップへ持ち越す必要がある。

 曲げモーメントの回転 

全体  から法線成分を引くと接線成分になる。これも直交化（Gram–Schmidt 法）であり，履歴ばねとして「前ステップまでに蓄えられた曲げバネ伸び（＝エネルギー）」をそっくり次ステップへ持ち越す必要がある。

 最新増分の加算 

最新ステップで計算された増分を加えて履歴ベクトルを更新する。

 結合破壊判定 

力とトルクによる破壊判定にはバネだけ考え，ダッシュポットは考えない。， はそれぞれ臨界法線応力，臨界接線応力である。

 法線応力破壊 

破壊条件は

である。

 せん断応力破壊 

破壊条件は

である。

 粒子方程式への寄与 

接触  から粒子  に作用する合力・合トルクは

接触  から粒子  に作用する合力・合トルクは

 脚注