(57)【要約】      （修正有）【課題】ＢＰＭ理論に基づくメゾの力チェーン粒子モデルの構築方法を提供する。【解決手段】物体の幾何モデルを空間で格子分割を行って、それぞれの格子中心の三次元座標を取得するステップと、得られた座標情報を使用して格子の数量に応じて相応数の剛性基礎粒子ユニットを生成し、かつ、粒子に固有属性を付与するステップと、ＢＰＭ理論によりそれぞれの剛性基礎粒子を結合し、相隣する粒子間が結合キーを形成し、かつ、結合キーに一定の力学特性を付与するステップと、構成した結合モデルに外力を印加し、力が結合キーに沿って内部へ伝達し、結合粒子モデルが外力の作用を受けたとき、力が結合キーに沿ってモデルの内部へ伝達するステップとを含む。【選択図】図７
2019106175.tif 000002
【特許請求の範囲】
【請求項１】
ＢＰＭ理論に基づくメゾの力チェーン粒子モデルの構築方法であって、
物体の幾何モデルを空間で格子分割を行って、それぞれの格子中心の三次元座標を取得するステップ１と、
得られた座標情報を使用して格子の数量に応じて相応数の剛性基礎粒子ユニットを生成し、かつ粒子に固有属性を付与するステップ２と、
ＢＰＭ理論によりそれぞれの剛性基礎粒子を結合し、相隣する粒子間が結合キーを形成し、かつ結合キーに一定の力学特性を付与し、その力学関係は、
2019106175.tif 000025
外力の作用に伴って、結合キーが時間ステップごとに発生した相応の弾性力及び力のモーメント増分は、
2019106175.tif 000026
平行キーに作用する最大引張応力及びせん断応力が、
2019106175.tif 000027
最大引張応力が引張強度を超え、即ち、
2019106175.tif 000028
であり、或いは最大せん断応力がせん断強度を超え、即ち
2019106175.tif 000029
であると、平行キーには破裂が生じて失効するステップ３と、
構成した結合モデルに外力を印加し、力が結合キーに沿って内部へ伝達し、上記公式に基いて結合粒子モデルを結び合わせて判断し、結合粒子モデルが外力の作用を受けたとき、力が結合キーに沿ってモデルの内部へ伝達し、形成した伝達経路が力チェーンと称され、力チェーンが強と弱に分けることができると共に、力が内部に伝達する過程において一定の伝達規律を表現するステップ４とを含むことを特徴とするＢＰＭ理論に基づくメゾの力チェーン粒子モデルの構築方法。
【請求項２】
前記ステップ４にて、強い／弱い力チェーンの判断方法は、
2019106175.tif 000030
力チェーン伝達過程における規律の統計方法は、
無次元キーにおける受力が
2019106175.tif 000031
であるように定義し、
統計学方法により、粒子系の内部の結合キーにおける力分布を統計分析し、その確率の計算方法が
2019106175.tif 000032
であり、
式中、Ｄｊがｆの区間[ｆｊ,ｆｊ＋１]
内の数量を示し、Ｄｆがサンプル総数を示すことを特徴とする請求項１に記載のＢＰＭ理論に基づくメゾの力チェーン粒子モデルの構築方法。
【請求項３】
前記ステップ１にて、結合キーが受けた平均力の計算方法は、
2019106175.tif 000033
であり、
式中、Ｎが粒子系の内部の結合キーの総数を示すことを特徴とする請求項１に記載のＢＰＭ理論に基づくメゾの力チェーン粒子モデルの構築方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、岩石・土壌工学の技術分野に関し、具体的には、ＢＰＭ理論に基づくメゾの力チェーン粒子モデルの構築方法に関し、特に、物質破砕過程をシミュレーションする粒子メゾ構造のモデル構築方法に関する。
【背景技術】
【０００２】
粒子材料の破砕は、複雑な変形過程であり、微小クラックの生じ、広がり及び交わりに伴って、接触理論、メゾメカニックス、統計力学などの複数の方面の知識を運用して破砕のミクロや、メゾ、マクロ変化過程を研究し、破砕の本質を掲示する必要がある。
【０００３】
メゾメカニックスの研究では、粒子材料が外部力の作用を受けるとき、力のその内部における伝達経路は、力チェーンと称される。国外の学者であるＲ．Ｈ．Ｅｗｏｌｄｔは、光弾性試験測定法により離散粒子系の内部粒子間の接触力を検出し、国内の学者である孫其誠らは、離散粒子系の内部粒子間の接触力分布特性及び接触力が形成するネットワークを研究したところ、粒子系が力を受けた条件で力チェーンのネットワークの形態進化を与える。楊涵らは、離散粒子媒体変形の局所のマクロメカニズム及びメゾメカニズムにより研究し、力チェーン構成をコリメーション効果が比較的よい力チェーンとコリメーション効果が比較的悪い力チェーンという２種に分ける。洪俊らは、離散粒子系を動力学解析することにより、離散粒子系の変形挙動と外部荷重との間の関係が得られる。劉恩龍らは、メゾメカニックスモデルを構築することにより、粒子集合体の破砕過程におけるメゾ応力−歪み構造方程式を与える。以上のメゾメカニックスの研究は、いずれも、離散粒子系について行うものであるが、実際の粒子材料は、異なる粒子と接着物質からなるので（例えば、コンクリート、鉱石など）、メゾ方面で粒子破砕についての研究が少ない。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明は、既存技術の不足を克服し、現実の必要に適応し、ＢＰＭ理論に基づくメゾの力チェーン粒子モデルの構築方法を提供することを目的とする。
【課題を解決するための手段】
【０００５】
本発明の目的を実現するために、本発明は、以下の技術案を採用する。
【０００６】
ＢＰＭ理論に基づくメゾの力チェーン粒子モデルの構築方法であって、以下のステップを含む。
【０００７】
ステップ１にて、物体の幾何モデルを空間で格子分割を行い、それぞれの格子中心の三次元座標を取得する。
【０００８】
ステップ２にて、得られた座標情報を使用して格子の数量に応じて相応数の剛性基礎粒子ユニットを生成し、かつ粒子に固有属性を付与する。
【０００９】
ステップ３にて、ＢＰＭ理論によりそれぞれの剛性基礎粒子を結合し、相隣する粒子間が結合キーを形成し、かつ結合キーに一定の力学特性を付与し、その力学関係は、
2019106175.tif 000003
【００１０】
外力の作用に伴って、結合キーが時間ステップごとに発生した相応の弾性力及び力のモーメント増分は、
2019106175.tif 000004
【００１１】
平行キーに作用する最大引張応力及びせん断応力が、
2019106175.tif 000005
【００１２】
最大引張応力が引張強度を超え、即ち、
2019106175.tif 000006
であり、或いは最大せん断応力がせん断強度を超え、即ち
2019106175.tif 000007
であると、平行キーには破裂が生じて失効する。
【００１３】
ステップ４にて、構成した結合モデルに外力を印加し、力が結合キーに沿って内部へ伝達し、上記公式に基いて結合粒子モデルを結び合わせて判断し、結合粒子モデルが外力の作用を受けたとき、力が結合キーに沿ってモデルの内部へ伝達し、形成した伝達経路が力チェーンと称され、力チェーンが強と弱に分けることができると共に、力が内部に伝達する過程において一定の伝達規律を表現する。
【００１４】
前記ステップ４にて、強い／弱い力チェーンの判断方法は、
2019106175.tif 000008
【００１５】
力チェーン伝達過程における規律の統計方法は、
無次元キーにおける受力が
2019106175.tif 000009
であるように定義する。
【００１６】
統計学方法により、粒子系の内部結合キーにおける力分布を統計分析し、その確率の計算方法が
2019106175.tif 000010
であり、
式中、Ｄｊがｆの区間[ｆｊ,ｆｊ＋１]
内の数量を示し、Ｄｆがサンプル総数を示す。
【００１７】
前記ステップ１にて、結合キーが受けた平均力の計算方法は、
2019106175.tif 000011
【００１８】
式中、Ｎが粒子系の内部の結合キーの総数を示す。
本発明は、以下の有益な効果を有する。
【００１９】
本方法により構築される粒子モデルは、外力を付勢するとき、力が物体内部の伝達経路に進化する情況及び内部力の伝達規律を体現することができるとともに、メゾで異なる付勢段階でのクラックの変化、物体破砕完了過程におけるメゾの力及びクラックの変化を示すことができるとともに、メゾでの物体破砕メカニズムの研究に対して一定の参考を提供する。
【図面の簡単な説明】
【００２０】
【図１】図１は、本発明に係るブラジリアンディスク（Ｂｒａｚｉｌｉａｎ  Ｄｉｓｃ）モデルの模式図である。
【図２】図２は、図１におけるブラジリアンディスクの幾何学格子分割の模式図である。
【図３】図３は、円盤粒子モデルである。
【図４】図４は、結合キーの空間フレームである。
【図５】図５は、タングステンの動態破砕するときの結合キーの断裂過程である。
【図６】図６は、粒子モデルの付勢過程における力チェーン伝達形態である。
【図７】図７は、異なる歪み段階の結合キーの力分布確率である。
【発明を実施するための形態】
【００２１】
以下、図面及び実施例を結び合わせて本発明をさらに説明する。
【００２２】
図１〜図７を参照して実施例を説明する。
【００２３】
タングステンのブラジリアンディスク粒子モデルの構築方法を例とし、本発明の技術案を詳しく説明し、本発明は、ＢＰＭ理論(ＢＰＭ理論−結合パラメータを定義し、一群の剛性球形粒子を結合し、相隣する２つの粒子間に結合キーを形成し、それにより、キーが粒子間の相対変位と回転を制限し、鉱石類重合体をシミュレーションすることができる。)に基づいてメゾの力チェーン粒子モデルを構築する方法を開示し、以下のステップを含む。
【００２４】
ステップ１にて、物体の幾何モデルを空間で格子分割を行い、それぞれの格子中心の三次元座標を取得する。図１に示すように、ブラジリアンディスク三次元幾何モデルを構築し、かつ図２に示すように、図１におけるブラジリアンディスクモデルに対して格子分割を行い、かつ各格子ユニットの中心点の三次元座標を抽出する。
【００２５】
ステップ２にて、得られた座標情報を使用して格子の数量に応じて相応数の剛性基礎粒子ユニットを生成し、かつ粒子に固有属性を付与する。図３に示すように、各剛性基礎粒子ユニットにタングステン属性を付加し、その基本属性は、表１に示される。
【００２６】
表１が剛性基礎粒子の固有属性である。
2019106175.tif 000012
【００２７】
ステップ３にて、ＢＰＭ理論によりそれぞれの剛性基礎粒子を結合し、相隣する粒子間が結合キーを形成し、図４は、各結合キーが円盤モデルの空間フレームを構成する構成図であり、かつ結合キーに一定の力学特性を付与し、その力学関係は、
2019106175.tif 000013
【００２８】
外力の作用に伴って、結合キーが時間ステップごとに発生した相応の弾性力及び力のモーメント増分は、
2019106175.tif 000014
【００２９】
平行キーに作用する最大引張応力及びせん断応力が、
2019106175.tif 000015
【００３０】
最大引張応力が引張強度を超え、即ち、
2019106175.tif 000016
であり、或いは最大せん断応力がせん断強度を超え、即ち
2019106175.tif 000017
であるとき、平行キーには破裂が生じて失効する。
【００３１】
本実施例における結合キーは、線弾性として表現し、結合キーの相応パラメータが表２に示される。
【００３２】
2019106175.tif 000018
【００３３】
図５は、タングステンの動態破砕するときの結合キーの断裂過程を示す。図には、ブラジリアンディスクが外部力を受けたとき、内部結合キーの空間構成ネットワークの断裂変化を体現し、それにより異なる付勢段階で、微小クラックの生じ、広がり及びマクロ断裂の進化を体現する。力の角度で力の内部空間構成における伝達経路を体現することもできる。
【００３４】
ステップ４にて、上記公式に基いて結合粒子モデルを結び合わせて判断し、結合粒子モデルが外力の作用を受けたとき、力が結合キーに沿ってモデルの内部へ伝達し、形成した伝達経路が力チェーンと称され、力チェーンが強と弱に分けることができると共に、力が内部に伝達する過程において一定の伝達規律を表現する。
【００３５】
強い／弱い力チェーンの判断方法は、
2019106175.tif 000019
【００３６】
力チェーン伝達過程における規律の統計方法は、
無次元キーにおける受力が
2019106175.tif 000020
であるように定義する。
【００３７】
統計学方法により、粒子系の内部結合キーにおける力分布を統計分析し、その確率計算方法が
2019106175.tif 000021
【００３８】
結合キーが受けた平均力の計算方法は、
2019106175.tif 000022
であり、
式中、Ｎが粒子系の内部の結合キーの総数を示す。
【００３９】
図６は、粒子モデル付勢過程における力チェーン伝達形態である。図には、ブラジリアンディスクが外部力を受けたとき、力の内部空間構成における伝達経路、力の分布特徴及び力チェーンのネットワークの進化規律を示す。
【００４０】
図７は、ブラジリアンディスクが力を受けたときに異なる歪み段階における力分布特徴を示す。図には、ブラジリアンディスクモデルは付勢過程において力の分布確率が指数関数的減衰特徴を呈することを示す。
【００４１】
図７は、力分布減衰関数は、
2019106175.tif 000023
であり、
式中、ａ、ｂは、フィッティングパラメータを示し、ｆが結合キーが受けた合力と平均合力との比を示す。
【００４２】
図７におけるデータをフィッティングすることによりフィッティングパラメータが得られ、表３に示される。
【００４３】
2019106175.tif 000024
【００４４】
以上、本発明の実施例であり、本発明の請求の範囲を限定しない。本発明の明細書及び図面に基づいて行う同等置換、或いは関連の技術分野に直接又は間接に応用されることは、本発明の請求の範囲に含まれる。
【図１】
2019106175.tif 000034
【図２】
2019106175.tif 000035
【図３】
2019106175.tif 000036
【図４】
2019106175.tif 000037
【図５】
2019106175.tif 000038
【図６】
2019106175.tif 000039
【図７】
2019106175.tif 000040
【外国語明細書】
2019106175000001.pdf