JP4414525B2

JP4414525B2 - BEHAVIOR BEHAVIOR PREDICTION METHOD AND BEHAVIOR PREDICTION DEVICE

Info

Publication number: JP4414525B2
Application number: JP32784799A
Authority: JP
Inventors: 正俊藤中
Original assignee: Bando Chemical Industries Ltd
Current assignee: Bando Chemical Industries Ltd
Priority date: 1999-11-18
Filing date: 1999-11-18
Publication date: 2010-02-10
Anticipated expiration: 2019-11-18
Also published as: JP2001139116A

Description

【０００１】
【発明の属する技術分野】
本発明は、石炭、鉱石等の塊状ばら物や粒状の搬送物をベルトコンベア装置により搬送する際の、該搬送物の挙動を予測する挙動予測方法及び挙動予測装置に関する。
【０００２】
【従来の技術】
従来より、例えば、鉱石やコークス等を溶鉱炉の高所にある投入口に連続的に投入するために、それらばら物搬送物を急傾斜又は垂直搬送可能な桟付コンベアベルトを備えたベルトコンベア装置が使用されている（図２参照）。このような桟付コンベアベルトは、例えば図３に示すように、平ベルトからなるベースベルト１（ベルト本体）の外周の搬送面１ａに、ベルト幅方向に延びる板状の横桟２２，２２，…をベルト長さ方向の全周に亘って所定間隔を空けて立設し、この横桟２２，２２，…により図示しない搬送物を係止して、急傾斜又は垂直状態で搬送できるようにしたものである。
【０００３】
ところで、前記のようなベルトコンベア装置の受入シュートでは、上流側のコンベアベルト等から搬送物が落下して、その下方のコンベアベルト上に受け止められるようになっており、このときにコンベアベルトに加わる衝撃力がベルトコンベア装置に作用する最大の負荷となることは知られている。そして、この衝撃力は、通常、古典的な物理式や搬送物の流体力学的な近似による噴流荷重の考え方によって推定されている。
【０００４】
また、そのように推定した搬送物の衝撃荷重に基づいて、ベルトコンベア装置のレイアウト設計や使用条件等が決定されるのであるが、この決定は通常、技術者の経験に頼ってなされており、例えば、前記受入シュートにおける搬送物の投下高さを低めに設定したり、或いは搬送量を所定以下に制限したりして、搬送物の落下による衝撃荷重をできるだけ軽減するようにしている。
【０００５】
また、前記受入シュートには、通常、搬送物が飛び散らないように案内する案内板や衝撃力を緩和するための干渉板等が配置されているが、これらの案内板や干渉板による効果を事前に予測することは困難なので、これらはベルトコンベア装置の設置後に搬送物を実際に投入して見て、その結果に基づいて位置を調節するようにしている。
【０００６】
【発明が解決しようとする課題】
しかし、前記のように、ベルトコンベア装置のレイアウトや使用条件等を経験に頼って決定するようにすると、ベルトコンベア装置の搬送量や搬送速度の増大に伴い衝撃力が大きくなったときに、コンベアベルトを実際に使用してから短時間で摩耗や破損を生じることがある。
【０００７】
特に、前記桟付コンベアベルトでは、受入シュートからの落下搬送物との衝突によって横桟が破損しやすく、また、搬送のばらつきによって一部の横桟に多量の搬送物が係止されると、急傾斜又は垂直搬送状態でその横桟が荷重を支え切れなくなって変形し、多量の搬送物が落下したり、横桟そのものが破損したりする虞れがある。
【０００８】
一方、このような問題の発生を避けるために、装置レイアウトの決定やベルトの仕様選定を安全側の設計にすると、ベルトコンベア装置のコスト高、重量増等の不具合を招く。
【０００９】
また、前記桟付コンベアベルトの場合、排出シュートから排出されるべき搬送物が横桟に引っかかってしまい、排出されずに残ることがあるので、横桟の形状や搬送速度等の設定によって、そのような搬送物の残留量をできるだけ低減したいという要請もある。
【００１０】
本発明は斯かる点に鑑みてなされたものであり、その目的とするところは、ベルトコンベア装置のコンベアベルトにより搬送される搬送物挙動の予測方法に工夫を凝らし、該搬送物挙動を正確に予測し得るようにして、ベルトコンベア装置のレイアウトや使用条件、及びコンベアベルトの仕様等を事前に適切に設定可能とすることにある。
【００１１】
【課題を解決するための手段】
前記目的を達成するために、本発明の解決手段では、コンベアベルトの搬送面を仮想の壁面部材により模擬する解析モデルを用意し、その壁面部材を移動させて、コンベアベルトの移動状態を記述しながら、搬送物の挙動を予測するようにした。
【００１２】
具体的に、請求項１の発明は、ベルトコンベア装置のコンベアベルトにより搬送されるばら物搬送物の挙動を予測するベルト搬送物の挙動予測方法を前提とする。そして、前記コンベアベルトの搬送面に搬送物を係止するための横桟が設けられている場合に、少なくとも前記コンベアベルトの搬送面を、ベルト長手方向に連続的に配置した複数の仮想の壁面部材により模擬するとともに、前記横桟の少なくともベルト進行方向前側の面を仮想の壁面部材により模擬する解析モデルを用意し、この解析モデルに対して、少なくとも、コンベアベルト及び搬送物のそれぞれの材料データと幾何データとを入力するとともに、該コンベアベルト上に供給される搬送物の挙動の初期条件を設定し、その上で、前記解析モデルにおいて壁面部材を、コンベアベルトの移動状態を記述するようにベルト長手方向に移動させながら、該壁面部材と共に移動する搬送物の挙動を前記初期条件に基づいて予測する。
【００１３】
この方法では、コンベアベルトの搬送面や横桟を複数の仮想の壁面部材により模擬する解析モデルにおいて、コンベアベルト上に供給される搬送物の挙動の初期条件を設定するとともに、該コンベアベルトの移動状態を記述するように壁面部材を移動させながら、この壁面部材と共に移動する搬送物の挙動を前記初期条件に基づいて予測することで、例えば、コンベアベルト上に落下する搬送物による衝撃力の分布状態や、コンベアベルト上に載置されて搬送される搬送物の挙動、或いはコンベアベルトから排出される搬送物の挙動等を予測することができ、この予測結果に基づいて、ベルトコンベア装置のレイアウトや使用条件等を決定することにより、ベルトコンベア装置の信頼性確保と重量やコストの低減とを両立することができる。
【００１４】
また、一般的に、前記のような横桟が設けられた桟付コンベアベルトでは、コンベアベルト上への落下搬送物の衝突によって、横桟が破損しやすいので、横桟の無いものに比べて、落下搬送物の挙動を正確に把握することの重要性が大きい。しかも、搬送ばらつきが大きくなると、一部の横桟にのみ過大な荷重が作用して、該横桟が破損したり、搬送物がこぼれ落ちたりすることが懸念されるので、この点からも搬送物の挙動を正確に把握することの重要性が大きい。よって、前記の如く搬送物挙動を予測できることが、特に有効な作用効果を奏するものである。
【００１５】
好ましいのは、搬送物の挙動を個別要素法の適用によって予測することであり、こうすれば、その予測精度を極めて高いものとして、請求項１の発明の作用効果を高めることができる。
【００１６】
また、好ましいのは、３次元の解析モデルを用意し、該解析モデルにおいてコンベアベルトの移動状態がそれぞれ異なる複数の領域を設定し、この各領域毎に壁面部材の移動手順を定義する。
【００１７】
こうすることで、まず３次元の解析モデルを用いることで、搬送物の挙動を空間的に予測して、請求項１の発明の作用効果を十分に得ることができる。また、コンベアベルトの移動状態が異なる領域毎に壁面部材の移動手順を定義することで、コンベアベルトがプーリ間を並進する移動状態やプーリに巻きかけられて周回する移動状態にあるときに、その各移動状態に対応するように壁面部材を移動させて、解析モデル上でコンベアベルトの移動状態を正確に記述することができる。
【００１８】
次に、請求項２の発明は、前記請求項１の発明に係るベルト搬送物の挙動予測方法と同じ前提において、少なくとも前記コンベアベルトの搬送面を、ベルト長手方向に連続的に配置した複数の仮想の壁面部材により模擬するとともに、複数のコンベアベルトをそれぞれの配置を含めて模擬する解析モデルを用意し、この解析モデルに対して、少なくとも、コンベアベルト及び搬送物のそれぞれの材料データと幾何データとを入力するとともに、該コンベアベルト上に供給される搬送物の挙動の初期条件を設定する。そして、前記解析モデルにおいて壁面部材を、コンベアベルトの移動状態を記述するようにベルト長手方向に移動させながら、該壁面部材と共に移動する搬送物の挙動を前記初期条件に基づいて予測し、かつ前記複数のコンベアベルト間の搬送物の乗り継ぎ挙動を個別要素法の適用により予測する。
【００１９】
こうすることで、複数のコンベアベルト間の搬送物の乗り継ぎ部において、上流側コンベアベルトの搬送面から下流側コンベアベルトの搬送面上に落下する搬送物の挙動を個別要素法の適用によって正確に求めることができる。
【００２０】
次に、請求項３の発明は、ベルトコンベア装置のコンベアベルトにより搬送されるばら物搬送物の挙動を予測するベルト搬送物の挙動予測装置を前提とする。そして、少なくとも前記コンベアベルトの搬送面をベルト長手方向に連続的に配置した複数の仮想の壁面部材により模擬する解析モデルと、この解析モデルに対して、コンベアベルトの搬送面上に供給される搬送物の挙動の初期条件を設定する初期条件設定手段と、前記解析モデルにおいて壁面部材を、コンベアベルトの少なくとも幾何データに基づいてベルト長手方向に移動させることにより、該コンベアベルトの移動状態を記述するベルト移動状態記述手段と、前記解析モデルにおける搬送物の挙動を、該搬送物の少なくとも材料データ及び幾何データと、前記初期条件設定手段により設定された搬送物挙動の初期条件と、前記ベルト移動状態記述手段により記述されるコンベアベルトの移動状態とに基づいて逐次記述する搬送物挙動記述手段とを備える構成とする。
【００２１】
その上さらに、前記複数のコンベアベルトをそれぞれの配置を含めて模擬するものとし、前記初期条件設定手段は、前記複数のコンベアベルトのうちの上流側コンベアベルトの搬送面上に所定の分布状態で搬送物を生成配置する搬送物生成手段を有するものとし、前記搬送物挙動記述手段は、前記上流側コンベアベルトの搬送面から下流側コンベアベルトの搬送面上へ落下する搬送物の乗り継ぎ挙動を、前記搬送物生成手段により配置された搬送物の初期分布状態と、ベルト移動状態記述手段により記述される上流側コンベアベルトの移動状態とに基づいて、個別要素法により逐次記述する構成とする。
【００２２】
前記の構成により、コンベアベルトの搬送面を複数の仮想の壁面部材により模擬する解析モデルにおいて、コンベアベルト上に供給される搬送物の挙動の初期条件が初期条件設定手段により設定されるとともに、ベルト移動状態記述手段により該コンベアベルトの移動状態を記述するように前記壁面部材が移動され、この壁面部材の移動状態、即ちコンベアベルトの移動状態と前記搬送物挙動の初期条件とに基づいて、搬送物挙動記述手段により、搬送物の挙動が逐次記述される。これにより、コンベアライン上の搬送物の挙動をきめ細かく予測することができ、請求項１の発明と同様の作用効果が得られる。
【００２３】
また、前記搬送物生成手段により上流側コンベアベルト上に所定の分布状態で搬送物が生成配置されるとともに、該上流側コンベアベルトの移動状態がベルト移動状態記述手段により記述され、この上流側コンベアベルトの移動状態と前記搬送物の初期分布状態とに基づいて、搬送物挙動記述手段により、上流側コンベアベルトから下流側コンベアベルトの搬送面上に落下する搬送物の挙動が個別要素法の適用によって逐次記述される。このことで、まず搬送物の挙動を個別要素法の適用によって予測することで、その予測精度を極めて高いものとして、前記の発明の作用効果を高めることができる。また、下流側コンベアベルト上に落下する搬送物の挙動を搬送物挙動記述手段によって記述することで、該落下搬送物による衝撃力の分布状態を正確に求めることができ、このことによっても前記の発明の作用効果を高めることができる。
【００２４】
【発明の実施の形態】
以下、本発明の実施形態を図面に基いて説明する。
【００２５】
（コンベアラインの構成）
図２は、本発明に係るベルト搬送物の挙動予測方法を適用するコンベアラインＬ（ベルトコンベア装置）の一例を示し、このコンベアラインＬは、例えば石炭や鉱石等のばら物搬送物を略垂直上方に搬送するよう、いわゆるＳ字レイアウトとされた下流側のコンベアベルトＢ１と、このコンベアベルトＢ１に搬送物を乗り継ぎ供給する上流側のコンベアベルトＢ２とを有する。そして、前記下流側コンベアベルトＢ１は、平ベルトからなるベースベルト１（ベルト本体）を有し、このベースベルト１が、コンベア頭部（図の上端部）に設けられたドライブプーリ２とコンベア尾部（図の下端部）に設けられたテイクアッププーリ３との間に巻き掛けられるとともに、変角ディスクプーリ４，４、変角馬の背ローラ５，５，…及びリタン変角プーリ６によりそれぞれ変角されている。そして、前記ドライブプーリ２が図示しないモータ等によリ図の反時計回り方向に回転されると、コンベアベルトＢ１がドライブプーリ２により駆動され、搬送物をベースベルト１外周の搬送面１ａ（図３参照）に載置して連続的に搬送するようになっている。
【００２６】
より詳しくは、前記コンベアベルトＢ１は、コンベア尾部においてテイクアッププーリ３と変角ディスクプーリ４との間に略水平に延びるように設けられ、該テイクアッププーリ３により図に矢印ｆで示すように張力を加えられるとともに、複数（図例では５つ）のインパクトローラ８，８，…により下側から支持されていて、上方の受入シュート９（搬送物の乗り継ぎ部）から投入された搬送物をベースベルト１の搬送面１ａで受け止めるようになっている。また、コンベアベルトＢ１は変角ディスクプーリ４の周りを上方に向かって略９０°周回し、そこから略鉛直上方に延びていて、搬送面１ａに載置した搬送物を図に矢印で示すようにコンベア頭部まで連続的に搬送する。尚、符号１０，１０，…は、垂直搬送状態のコンベアベルトＢ１の内周側を支持するように上下方向に略等間隔に配置された複数（図例では６つ）のキャリアローラを示す。
【００２７】
さらに、前記コンベア頭部において、コンベアベルトＢ１は、ベルト内周側に所定のカーブを描くように配設された複数（図例では６つ）の変角馬の背ローラ５，５，…に沿って周回し、そこからドライブプーリ２まで略水平に延びて、該ドライブプーリ２の回りを図の反時計回り方向に略１８０度周回して、折り返されている。これにより、搬送物は、コンベアベルトＢ１とともにドライブプーリ２の回りを周回しながら、該ドライブプーリ２を取り囲む排出シュート１１内に落下して、排出される。一方、そのように搬送物が落下した後のコンベアベルトＢ１は、ドライブプーリ２から変角ディスクプーリ４まで略水平に延びていて、ここで、搬送面１ａに付着している搬送物がビータクリーナ１２によってふるい落とされるようになっている。
【００２８】
最後に、コンベアベルトＢ１は、変角ディスクプーリ４の周りを図の下側に向かって略９０°周回し、複数（図例では５つ）のリタンローラ１３，１３，…により内周面を支持されながら、略鉛直下方に延びて、コンベア尾部においてリタン変角プーリ６の周りを周回した後、そこから略水平に延びて、再びテイクアッププーリ３に至る。尚、図の符号１４，１４，１４は、コンベアベルトＢ１をベルト幅方向両側から支持するサイドローラを示し、また、符号１５，１６は、それぞれベースベルト１の内周面に付着した搬送物をそぎ落とすスクレーパを示している。さらに、符号１７は、搬送面１ａに残った搬送物がコンベア戻り側で落下したときに、この搬送物が飛び散るのを防ぐためのリターン案内板を示している。
【００２９】
前記コンベアベルトＢ１は、搬送物を係止して急傾斜状態でも上方に搬送可能な桟付きのものである。詳しくは図３に示すように、ベースベルト１の外周面（図の上面）における中央の所定幅の範囲が搬送面１ａとされ、この搬送面１ａのベルト幅方向両側には、それぞれ搬送物の落下を防止するための耳桟２１，２１が、ベルト長さ方向の全周に亘って設けられている。この一対の耳桟２１，２１は、それぞれベースベルト１の外周面から略垂直に突出していて、ベースベルト５が変角ディスクプーリ４，４や変角馬の背ローラ５，５，…等に沿って周回するときに伸縮できるように、蛇腹状に形成されている。また、前記左右両側の耳桟２１，２１よりもさらに外側に、それぞれ変角ディスクプーリ４，４の外周面に当接される当接面１ｂ，１ｂが形成されていて、ベースベルト１は、該当接面１ｂ，１ｂにおいて変角ディスクプーリ４，４に巻きかけられている。
【００３０】
さらに、前記ベースベルト１の搬送面１ａ上には、両側の耳桟２１，２１の間でベルト幅方向に延びる略矩形板状の横桟２２，２２，…が設けられている。この横桟２２，２２，…は、ベルト長さ方向の全周に亘って互いに所定間隔をあけて立設されていて、それぞれ、ベースベルト１の搬送面１ａから略垂直に突出する垂直壁部と、この垂直壁部の先端からベルト進行方向前側に向かって斜めに延びる折曲部とからなる。そして、横桟２２の折曲部の傾斜角度を適切に設定することで、この横桟２２によりかなり多くの搬送物を係止して、こぼれ落ちないように保持しながら、急傾斜又は垂直搬送状態であっても上方に搬送できるようになっている。
【００３１】
一方、前記上流側のコンベアベルトＢ２（以下、供給ベルトという）は平ベルトからなり、ドライブプーリ１８と図示しないドリブンプーリとに巻掛けられていて、該両プーリ間に略水平に延びている。そして、前記ドライブプーリ１８が図示しないモータ等によリ図の反時計回り方向に回転されると、そのドライブプーリ１８により供給ベルトＢ２が駆動され、この供給ベルトＢ２上の搬送物が図の左側に向かって連続的に搬送されて、コンベア終端部から受入シュート９に投入されるようになっている。尚、図の１９は、供給ベルトＢ２から受入シュート９に向かって落下する搬送物を案内する案内板である。
【００３２】
ところで、上述の如きコンベアラインＬのレイアウトや使用条件、或いはベルトＢ１，Ｂ２の仕様等は、従来までは主に技術者の経験に頼って決定されていたので、ベルトの信頼性確保と重量やコストの低減とを両立させることが難しかった。例えば、コンベアベルトＢ１は受入シュート９における落下搬送物によって大きな衝撃力重を受けることが知られているが、この衝撃力の分布状態や局所的な最大値等は正確には分からないので、ベルトの早期破損等の不具合を避けるためには、供給ベルトＢ２の高さ位置をかなり低めに制限したり、或いはコンベアベルトＢ１の仕様を過剰に安全側に設定するといった対応をせざるを得ず、結果として、コスト高や重量増といった悪影響が大きかった。
【００３３】
特に、前記のように垂直搬送が可能な桟付コンベアベルトＢ１を用いる場合、落下搬送物による横桟２２の損傷が大きいことから、この横桟２２そのものやベースベルト１との接着部分を過剰に補強することが多く、これによるベルトの重量増が大きな問題になっていた。また、供給ベルトＢ２から落下する搬送物の分量ばらつきが大きいと、この搬送物の垂直搬送状態で一部の横桟２２，２２，…に過大な荷重がかかることになり、この横桟２２が大きく変形して、搬送物がこぼれ落ちることによって、コンベアベルトＢ１が破損する虞れもあった。
【００３４】
さらに、排出シュート１１においては、上述の如く、コンベアベルトＢ１がコンベア頭部のドライブプーリ２の回りを周回する間に、搬送物が排出シュート１１内に落下するようになっているが、横桟２２，２２，…があると、この横桟２２，２２，…に搬送物が引っかかってしまい、落下せずに残留するという不具合もあった。このような不具合は、この実施形態のように途中で折れ曲った横桟２２，２２，…を設けた場合に特に顕著になり、また、ベルトの搬送速度が大きいほど残留量の増える傾向があるので、搬送物の残留量を減らすためには、横桟２２，２２，…の折曲部の傾斜角度やベルト搬送速度の設定にも十分な注意を払う必要がある。
【００３５】
（搬送物の挙動の予測）
斯かる実状に鑑み、この実施形態では、前記コンベアラインＬの設計段階において、コンベアラインＬと搬送物とを例えば図４に示すように模擬する解析モデルを用意し、この解析モデル上で、搬送物に対応する球状粒子の状態を個別要素法（Distinct Elemennt Method）によりシミュレーションするようにした。これにより、受入シュート９を介してコンベアベルトＢ１の搬送面１ａ上に投下する搬送物の挙動や、この搬送物がコンベアベルトＢ１により垂直搬送されて、排出シュート１１にて排出されるまでの挙動を高精度に予測することができ、この予測結果に基づいてコンベアラインＬのレイアウトや使用条件等を決定することで、前記のような不具合を解消することができる。
【００３６】
すなわち、この実施形態の搬送物挙動予測装置Ａは、概略図１の理論ブロック図にも示すように、前記コンベアベルトＢ１及び供給ベルトＢ２等を模擬する解析モデルＭを備えており、その解析モデルＭにおいて、上流側の供給ベルトＢ２の搬送面上に所定の分布状態で搬送物を生成配置する搬送物生成手段３０と、該供給ベルトＢ２及びコンベアベルトＢ１の移動状態をそれぞれ記述するベルト移動状態記述手段３１とが設けられ、さらに、前記供給ベルトＢ２及びコンベアベルトＢ１上でそれぞれ搬送される搬送物の挙動を、個別要素法により逐次記述する搬送物挙動記述手段３２が設けられている。
【００３７】
詳しくは、前記解析モデルＭは、前記図４に示すように、コンベアベルトＢ１の搬送面１ａや横桟２２，２２，…及び供給ベルトＢ２の搬送面を仮想の剛体壁からなる壁面部材ｗ，ｗ，…により模擬し、また、該供給ベルトＢ２からコンベアベルトＢ１へ投下される搬送物の案内板１９等を仮想の剛体壁からなる平面部材wallで模擬するとともに、搬送物そのものは所定の弾性及び減衰特性を有する球状粒子ptclにより模擬するものである。また、前記ベルト移動状態記述手段３１は、解析モデルＭにおいて供給ベルトＢ２及びコンベアベルトＢ１をそれぞれ構成する壁面部材ｗ，ｗ，…を、ベルトレイアウトデータ等に基づいてベルト長手方向に移動させることにより、前記供給ベルトＢ２及びコンベアベルトＢ１の移動状態をそれぞれ記述するものである。
【００３８】
さらに、前記搬送物挙動記述手段３２は、解析モデルＭ上で、上流側の供給ベルトＢ２の搬送面から下流側コンベアベルトＢ１の搬送面１ａ上へ落下する搬送物の乗り継ぎ挙動を、該搬送物の少なくとも材料データ及び幾何データと、前記搬送物生成手段３０により配置された搬送物の初期分布状態と、前記ベルト移動状態記述手段３１により記述される供給ベルトＢ２の移動状態とに基づいて、個別要素法により逐次記述する。また、該搬送物挙動記述手段３２は、前記のように記述した落下搬送物の挙動と、該搬送物の少なくとも材料データ及び幾何データと、前記ベルト移動状態記述手段３１により記述されるコンベアベルトＢ１の移動状態とに基づいて、コンベアベルトＢ１が受入シュート９から排出シュート１１まで移動する間の一連の搬送物の挙動を、個別要素法により逐次記述する。
【００３９】
そして、前記搬送物挙動記述手段３２からの出力結果として、受入シュート９においてコンベアベルトＢ１上に落下する搬送物の位置及び速度データを該コンベアベルトＢ１のＦＥＭ解析の初期条件として与えるとともに、受入シュート９から排出シュート１１に至るコンベアベルトＢ１の移動状態及び搬送物挙動をそれぞれディスプレー装置（画像表示装置）に出力して、表示させるようにしている。
【００４０】
ここで、個別要素法とは、多数の粒状物が集合してなる不連続体の挙動を解析するための計算手法の一種であり、この実施形態のようなばら物搬送物の挙動を従来までの連続体近似による計算手法等と比較して高精度に求めることができるものである。
【００４１】
具体的に、粒状物の要素を球状粒子ptclとした場合について、図５に基づいて説明すると、同図（ａ）に示すように、或る時刻ｔにおいて、座標位置Ａに並進速度ｖA、回転速度ωAの球状粒子ptclが存在し、隣り合う他の球状粒子（図に仮想線で示す）と点Ｐにおいて接しているとする。この場合、まず、微小時間Δｔが経過した後の球状粒子ptclの位置Ａ′及び速度ｖA′，ωA′を求めるために、そのときに球状粒子ptclに作用している外力Ｆが微小時間Δｔの経過する間一定であると仮定し、この外力Ｆ、座標位置Ａ、並進速度ｖA、及び回転速度ωAに基づいて、球状粒子ptclの一般的な運動方程式を解くことで、移動後の座標位置Ａ′、並進速度ｖA′、及び回転速度ωA′を演算する。
【００４２】
次に、前記の運動方程式の演算を解析モデルにおける全ての球状粒子について行った後、各球状粒子の位置座標に基づいて、隣接する球状粒子同士の干渉度合いを調べる。すなわち、同図（ｂ）に示すように、隣接する球状粒子間で重なりや滑りが発生しているときには、その重なり量ｄ及び滑り量Δｓに加えて、各球状粒子の弾性や摩擦係数等の材料特性を考慮し、一般的な線形ばねモデルに基づく力・変位関係式を解いて、新たな接触点ｐ″における半径方向及び接線方向の接触力Ｆn，Ｆsをそれぞれ演算する。この力・変位関係式は、例えば接触点Ｐ″における半径方向及び接線方向のそれぞれについて、減衰のある１自由度の線形ばねモデルを仮定し、フックの法則により力と変位との関係を定めたものとすればよい。
【００４３】
次に、球状粒子ptclに対し図示しないがその周囲に隣接する全ての球状粒子からそれぞれ作用する接触力Ｆn，Ｆsの合力を求め、この合力と球状粒子ptclに作用する重力とを該球状粒子ptclの新しい運動方程式における外力Ｆとして、この新たな運動方程式を解く。すなわち、図６に模式的に示すように、前記運動方程式の演算結果に基づく球状粒子等の位置座標の更新、この更新結果に基づく力・変位関係式の演算、及びその演算結果に基づく接触力の更新、この更新された接触力に基づく新たな運動方程式の演算というように、微小時間Δｔ毎に運動方程式及び力・変位関係式を交互に解くことで、球状粒子ptclの挙動を時々刻々と正確に記述することができるのである。
【００４４】
以下、個別要素法の適用によって、移動するコンベアベルトＢ１，Ｂ２とこのコンベアベルトＢ１，Ｂ２上の搬送物の挙動とを記述するシミュレーション制御の手順について、図７〜図９に基づいて具体的に説明する。
【００４５】
（全体的な手順）
まず、シミュレーション制御の全体的な手順は図７のフローチャート図に示すようになり、このメインフローにおけるスタート後のステップＳ１において、コンベアラインＬの基本的なレイアウトや球状粒子ptclの直径（幾何データ）、それらの密度、ばね定数、減衰率及び摩擦係数（材料データ）、球状粒子ptclの個数、ベルト移動速度（使用条件）等のデータを入力する。尚、前記コンベアラインＬの基本的なレイアウトとしては、例えば、供給ベルトＢ２及び搬送用コンベアベルトＢ１の寸法及び配置、該コンベアベルトＢ１上の横桟２２の間隔、並びに該横桟２２の寸法形状等である。
【００４６】
続いて、ステップＳ２において、前記ステップＳ１で入力したレイアウトデータに基づいて、コンベアベルトＢ１及び供給ベルトＢ２の基本レイアウトを生成する。即ち、例えば図９に模式的に示すコンベアレイアウト図において、レイアウト基準点Ｐ０〜Ｐ５がそれぞれ決定されて、前記図４に示すような３次元の仮想のコンベアラインＬが構成される。この仮想のコンベアラインＬにおいて、下流側のコンベアベルトＢ１については、ベースベルト１の搬送面１ａを表す２枚の壁面部材ｗ，ｗと、横桟２２の垂直壁部及び折曲部の表裏両面をそれぞれ表す４枚の壁面部材ｗ，ｗ，…との合計６枚の仮想の壁面部材ｗを１つの構成単位とし、この構成単位を排出シュート１１付近から受入シュート９付近までベルト長手方向に連続的に配置することで、コンベアベルト全体を模擬するようにしている。また、上流側の供給ベルトＢ２については、同様にその搬送面を表す壁面部材ｗ，ｗ，…をベルト長手方向に連続的に配置して、ベルト全体を模擬するようにしている。
【００４７】
続いて、ステップＳ３では、シミュレーション制御の制御サイクルのカウント値ｉが最終値ｉendになったかどうか判定し、この判定がＹＥＳであれば制御を終了する一方、判定がＮＯならばステップＳ４に進む。このステップＳ４では、供給ベルトＢ２の搬送面上に所定数の球状粒子ptclを例えばガウス分布になるように生成配置し、続くステップＳ５において供給ベルトＢ２を移動させ（供給用コンベアベルト移動）、続くステップＳ６において、コンベアベルトＢ１を移動させる（搬送用コンベアベルト移動）。この両ベルトＢ１，Ｂ２の移動状態を記述するサブルーチンについては後述する。
【００４８】
続いて、ステップＳ７において、前記供給ベルトＢ２及びコンベアベルトＢ１の移動に伴い搬送される球状粒子ptclの挙動を、上述の如き個別要素法の適用によって演算する。すなわち、微小時間Δｔの経過する間の球状粒子ptcl及び壁面部材ｗに関する運動方程式と力・変位関係式とを解いて、球状粒子ptclの位置座標や並進速度、回転速度等を更新する。続いて、ステップＳ８では、球状粒子ptcl及び壁面部材ｗの位置座標をディスプレー装置に出力し、或いは球状粒子ptcl位置座標をＦＥＭ解析装置に出力して、続くステップＳ９において、制御サイクルのカウント値ｉをインクリメントし、しかる後に前記ステップＳ３にリターンする。
【００４９】
そうして、前記ステップＳ３からステップＳ９までの制御手順を繰り返し行うことで、球状粒子ptclが供給ベルトＢ２に搬送されて水平方向に移動し、この供給ベルトＢ２からコンベアベルトＢ１に向かって落下し、その後、該コンベアベルトＢ１のベースベルト１の搬送面１ａ上に載置されて略鉛直上方に搬送され、排出シュート１１において、コンベアベルトＢ１から離れて落下するというような、一連の球状粒子ptclの位置座標や速度を時々刻々と求めることができる。つまり、コンベアラインＬにおける搬送物に対応する球状粒子ptcl群の挙動を正確に記述することができる。
【００５０】
前記図７に示すメインフローのステップＳ４は、解析モデルＭにおいて、供給ベルトＢ２上に所定の分布状態で球状粒子ptcl（搬送物）を生成配置する搬送物生成手段３０に対応している。
【００５１】
また、ステップＳ５及びステップＳ６は、それぞれ、前記解析モデルＭにおける供給ベルトＢ２及びコンベアベルトＢ１の移動状態を、該解析モデルＭに入力されたベルトＢ１，Ｂ２の幾何データ等に基づいて、逐次記述するベルト移動状態記述手段３１に対応している。
【００５２】
さらに、ステップＳ７は、前記解析モデルＭにおいて、解析モデルＭ２に入力された球状粒子ptclの少なくとも材料データ及び幾何データと、搬送物生成手段３０により供給ベルトＢ２上に配置された球状粒子ptclの初期分布状態と、ベルト移動状態記述手段３１により記述された供給ベルトＢ２の移動状態とに基づいて、該供給ベルトＢ２からコンベアベルトＢ１上に落下する球状粒子ptclの挙動を個別要素法により逐次記述するとともに、このようにして記述された球状粒子ptcl（落下搬送物）の挙動を初期条件として、コンベアベルトＢ１により受入シュート９から排出シュート１１まで搬送される球状粒子ptclの挙動を、同様に個別要素法により記述する搬送物挙動記述手段３２に対応している。
【００５３】
加えて、前記ステップＳ９は、前記ベルト移動状態記述手段３１により記述されたコンベアベルトＢ１及び供給ベルトＢ２の移動状態と、前記搬送物挙動記述手段３２により記述された球状粒子ptclの挙動とを、それぞれ、図外のディスプレー装置に表示させる画像処理手段に対応している。
【００５４】
（ベルト移動状態の記述）
次に、前記図７に示すメインフローのステップＳ５，Ｓ６において、供給ベルトＢ２及びコンベアベルトＢ１の移動状態をそれぞれ記述するサブルーチンについて、前記図４を参照しながら、図８及び図９に基づいて詳細に説明する。尚、以下の説明においては、図４及び図９における左右方向をＸ方向、図の紙面に直交する方向をｙ方向、図の上下方向をｚ方向とも呼ぶものとする。
【００５５】
まず、図８に示すフローのスタート後のステップＴ１では、仮想の壁面部材ｗの座標位置を代表する配列変数ｗｐの先頭位置へ移動する。続いて、ステップＴ２において、該配列変数ｗｐの値が空（ｗｐ null）かどうか判別し、この判別結果がＹＥＳならば制御を終了する（ＥＮＤ）一方、判別結果がＮＯならば、ステップＴ３に進み、当該壁面部材ｗが供給ベルトＢ２及びコンベアベルトＢ１のいずれを構成するものかを表すＩＤと該壁面部材ｗの現在の座標値とを読み込む。
【００５６】
続いて、ステップＴ４では、前記ステップＴ３で読み込んだＩＤに基づいて、壁面部材ｗが供給ベルトＢ２を構成するものかどうか判定し、この判定結果がＮＯならば、図示しないが、コンベアベルトＢ１の移動状態を記述する制御手順に進む一方、判定結果がＹＥＳで供給ベルトＢ２であれば、ステップＴ５に進んで供給ベルトＢ２の移動状態の記述を開始する。すなわち、まず、前記ステップＴ３で読み込んだ座標値に基づいて、壁面部材ｗの座標値がレイアウト基準点Ｐ１よりも下側（図９における下側：以下、同様とする）にあるかどうか判定する。ここで、レイアウト基準点Ｐ１は、供給ベルトＢ２を駆動するドライブプーリ１８の回転軸心に対応する位置にあり、前記の判定結果がＹＥＳならば、ステップＴ６に進んで壁面部材ｗを消去し（壁消去）、ステップＴ２にリターンする一方、判定結果がＮＯならばステップＴ７に進む。つまり、解析モデルＭ２おいて供給ベルトＢ２のドライブプーリ１８の回転軸心よりも低い位置については、該供給ベルトＢ２は記述しない。
【００５７】
一方、ステップＴ７では、今度は壁面部材ｗの座標値が前記点Ｐ１よりも左側にあるかどうか判定し、この判定結果がＮＯならばステップＴ１１に進む一方、判定結果がＹＥＳならばステップＴ８に進んで、ｗｐのｘ方向並進速度を一旦、零とし、かつ点Ｐ１の周りを図の反時計周り方向に回転するように、一定の角速度を与える（Ｐ１を中心に回転）。そして、ステップＴ９において、前記一定の角速度に基づいて壁面部材ｗの移動量を計算し、続くステップＴ１０では、前記ステップＴ３で読み込んだ座標値とステップＴ９で計算した移動量とに基づいて、配列変数ｗｐの座標値を新しい値に更新して、しかる後にステップＴ２にリターンする。つまり、点Ｐ１よりも左側の区間（イ）では、壁面部材ｗを点Ｐ１を中心として一定の角速度で回転移動させる。これにより、供給ベルトＢ２はドライブプーリ１８の周りを回転する移動状態になる。
【００５８】
また、前記ステップＴ７において、壁面部材ｗがレイアウト基準点Ｐ１よりも右側にあると判定して進んだステップＴ１１では、ｗｐにＸ方向の一定の並進速度を与え、続くステップＴ１２において壁面部材ｗの移動量を計算する。続いて、ステップＴ１３において、前記ステップＴ３で読み込んだ座標値とステップＴ１３で計算した移動量とに基づいて、壁面部材ｗの座標値を新しい値に更新し、しかる後にステップＴ２にリターンする。つまり、点Ｐ１よりも右側の区間（ロ）では、壁面部材ｗを一定の並進速度でＸ方向に並進させて、供給ベルトＢ２を一定速度の並進移動状態とする。尚、点Ｐ０よりも右側の区間に（ハ）に壁面部材ｗが１つも無くなれば、新しい壁面部材ｗを生成する。
【００５９】
つまり、供給ベルトＢ２については、該供給ベルトＢ２の先頭位置からベルト進行方向後側に向かって順番に壁面部材ｗを配置し、この壁面部材ｗをそれぞれ進行方向前側に移動させることにより、該供給ベルトＢ２の移動状態を記述する。また、供給ベルトＢ２の先頭の壁面部材ｗがドライブプーリ１８の回転軸心に対応する基準点Ｐ１よりも低い位置まで進めば、この壁面部材ｗを消去する一方、最後尾の壁面部材ｗが生成基準点Ｐ０よりも進行方向前側に移動したときには、その後ろに新しい壁面部材ｗを生成する。尚、詳しい説明は省略するが、コンベアベルトＢ１の移動状態についても同様に壁面部材ｗ，ｗ，…を用いて記述すればよい。
【００６０】
（実施形態の作用効果）
したがって、この実施形態に係るベルト搬送物の挙動予測装置Ａを使用して、上述の如く、３次元の解析モデルＭにおいて壁面部材ｗ，ｗ，…によりコンベアベルトＢ１及び供給ベルトＢ２の移動状態を記述しながら、この壁面部材ｗ，ｗ，…の移動に伴い搬送される球状粒子ptclの挙動を個別要素法の適用によってシミュレーションすることで、それら球状粒子ptclに対応するばら物搬送物の搬送中の挙動を３次元的に正確に予測することができる。
【００６１】
これにより、例えば、受入シュート９から投下された搬送物によりコンベアベルトＢ１に作用する衝撃力をその分布状態も含めてきめ細かくかつ正確に把握することができるので、この衝撃力の分布状態を初期条件として、例えばＦＥＭ解析によりコンベアベルトＢ１の応力分布状態を求め、この解析結果に基づいて、コンベアベルトＢ１の仕様や供給ベルトＢ２の高さ位置等のレイアウト、或いは最大搬送量等の使用条件等を最適に決定することで、コンベアラインＬの信頼性確保と重量やコストの低減とを高次元で両立させることができる。しかも、前記解析モデルＭにおいてはベルトＢ１，Ｂ２だけでなく、案内板１９等も模擬されているので、コンベアベルトＢ１上への落下搬送物による衝撃力を極めて正確に求めることが可能になる。
【００６２】
また、垂直搬送中にコンベアベルトＢ１の横桟２２，２２，…からこぼれ落ちる搬送物の挙動や、排出シュート１１においてコンベアベルトＢ１の横桟２２，２２，…に引っかかって残留する搬送物の挙動も正確に把握することができるので、そのようにこぼれ落ちたり残留する搬送物の量が少なくなるよう横桟２２形状の設計を変更したり、平均の搬送量を適度に制限したり、或いは、搬送量のばらつきが小さくなるよう供給ベルトＢ２の高さや案内板１９の位置を修正するというように、コンベアレイアウトや使用条件の設定をコンベアラインＬの設置前に適切に行うことができる。
【００６３】
さらにまた、この実施形態では、前記コンベアベルトＢ１及び供給ベルトＢ２の移動状態や搬送物の挙動がディスプレー装置に表示されるので、従来まではコンベアラインＬの設置後でも直接、見ることは難しかったばら物搬送物の流れを設計技術者が視覚的に捉えることができるようになり、受入シュート９から排出シュート１１に至る搬送物の流れを直感的に把握して、設計の一部修正等の適切な対策を迅速に行うことができるという優れた設計支援効果が得られる。
【００６４】
加えて、この実施形態では、コンベアラインＬをシミュレーションする際に、コンベアベルトＢ１や供給ベルトＢ２の全体を記述することはせず、ベルト搬送面の一部や案内板等の搬送物挙動に直接的に影響する部分だけを記述するようにしているので、シミュレーション制御の計算量を減らして、計算時間を短縮することができる。
【００６５】
尚、本発明は前記実施形態に限定されるものではなく、その他種々の実施形態を包含するものである。すなわち、前記実施形態では本発明のベルト搬送物挙動予測方法を桟付コンベアベルトを用いたコンベアラインＬに適用しているが、これに限らず、本発明は横桟の無いものを含めて、種々のコンベアベルトに適用することができる。
【００６６】
また、前記実施形態では、解析モデルＭとして供給ベルトＢ２を含むものを用意して、該供給ベルトＢ２からコンベアベルトＢ１に落下する搬送物の挙動も個別要素法により記述するようにしているが、これに限らず、供給ベルトＢ２についてのシミュレーションは行わずに、例えば、コンベアベルトＢ１上に落下する球状粒子ptclの空間分布と速度分布を与えるようにしてもよい。
【００６７】
【発明の効果】
以上説明したように、請求項１の発明におけるベルト搬送物の挙動予測方法によると、コンベアベルトの搬送面を仮想の壁面部材により模擬する解析モデルを用意し、その壁面部材を移動させて、コンベアベルトの移動状態を記述しながら、搬送物の挙動を予測するようにしたので、その搬送物挙動の予測結果に基づいて、ベルトコンベア装置のレイアウトや使用条件、及びコンベアベルトの仕様等を事前に適切に設定することができる。これにより、ベルトコンベア装置の信頼性確保と重量やコストの低減とを両立できる。
【００６８】
特に、桟付コンベアベルトに適用することで、前記の如く搬送物挙動を正確に予測できることが、特に有効な効果を奏する。
【００６９】
請求項２の発明によると、複数のコンベアベルト間の搬送物の乗り継ぎ挙動を個別要素法の適用によって正確に求めることができ、これにより、本発明の効果がさらに有効なものとなる。
【００７０】
請求項３の発明におけるベルト搬送物の挙動予測装置によると、個別要素法の適用により搬送物の挙動を極めて高精度に予測することができ、また、複数のコンベアベルト間の搬送物の乗り継ぎ部で上流側コンベアベルトから下流側コンベアベルト上に落下する搬送物の挙動をきめ細かくかつ正確に求めることができる。
【図面の簡単な説明】
【図１】本発明の実施形態に係るベルト搬送物の挙動予測装置Ａの構成を示す理論ブロック図である。
【図２】本発明を適用するコンベアラインの一例を示す概略構成図である。
【図３】桟付コンベアベルトの構成を拡大して示す斜視図である。
【図４】解析モデルにおける仮想のコンベアライン及び搬送物を示す斜視図である。
【図５】個別要素法における要素粒子の運動を示す模式図である。
【図６】個別要素法における演算手順の概念図である。
【図７】搬送物の挙動予測方法の全体的な手順を示すフローチャート図である。
【図８】供給ベルトの移動状態を記述する手順を示すフローチャート図である。
【図９】解析モデルにおけるコンベアラインの基本的なレイアウト設定を示す模式図である。
【符号の説明】
Ａ搬送物挙動予測装置
Ｂ１コンベアベルト（下流側コンベアベルト）
Ｂ２供給ベルト（上流側コンベアベルト）
Ｌコンベアライン（ベルトコンベア装置）
Ｍ解析モデル
１ａコンベアベルトの搬送面
２２横桟
３０搬送物生成手段（初期条件設定手段）
３１ベルト移動状態記述手段
３２搬送物挙動記述手段
ptcl 球状粒子（搬送物）
ｗ壁面部材[0001]
BACKGROUND OF THE INVENTION
  The present invention relates to a behavior predicting method and a behavior predicting device for predicting the behavior of a bulk material such as coal and ore or a granular transported material when transported by a belt conveyor device.
[0002]
[Prior art]
  Conventionally, for example, in order to continuously put ore, coke, etc. into the inlet of the blast furnace at a high place, there is a belt conveyor apparatus provided with a conveyor belt with a beam capable of steeply or vertically conveying these bulk goods. Used (see FIG. 2). For example, as shown in FIG. 3, such a conveyor belt with a cross is formed as a plate-like horizontal cross 22, 22,... Are set up at predetermined intervals over the entire circumference in the belt length direction, and a conveying object (not shown) is locked by the horizontal rails 22, 22,. Is.
[0003]
  By the way, in the receiving chute of the belt conveyor device as described above, the conveyed product falls from the upstream conveyor belt or the like and is received on the conveyor belt below it, and is added to the conveyor belt at this time. It is known that the impact force is the maximum load acting on the belt conveyor device. And this impact force is usually estimated by the concept of the jet load by the classical physical formula and the hydrodynamic approximation of the conveyed product.
[0004]
  In addition, the layout design and usage conditions of the belt conveyor device are determined based on the estimated impact load of the conveyed product, but this determination is usually made by relying on the experience of an engineer, For example, the dropped load height of the transported object in the receiving chute is set low, or the transport amount is limited to a predetermined value or less, so that the impact load due to the fall of the transported object is reduced as much as possible.
[0005]
  In addition, the receiving chute is usually provided with a guide plate for guiding the transported object so as not to scatter, an interference plate for reducing the impact force, and the like. Since these are difficult to predict, they are actually loaded after the belt conveyor device is installed, and the position is adjusted based on the result.
[0006]
[Problems to be solved by the invention]
  However, as described above, if the layout and usage conditions of the belt conveyor device are determined based on experience, when the impact force increases with an increase in the conveyance amount or conveyance speed of the belt conveyor device, the conveyor Wear and breakage may occur in a short time after the belt is actually used.
[0007]
  In particular, in the conveyor belt with crosspieces, the crosspieces are likely to be damaged by collision with the fallen conveyed item from the receiving chute, and when a large amount of conveyed items are locked to some of the horizontal crosspieces due to conveyance variations, In the inclined or vertical conveyance state, the horizontal beam may not be able to support the load and deforms, and there is a possibility that a large amount of conveyed items may fall or the horizontal beam itself may be damaged.
[0008]
  On the other hand, in order to avoid the occurrence of such problems, if the design of the device layout and the specification of the belt are designed on the safe side, problems such as high cost and increased weight of the belt conveyor device are caused.
[0009]
  Further, in the case of the conveyor belt with a cross, the conveyed product to be discharged from the discharge chute may be caught by the horizontal cross and remain without being discharged. There is also a demand to reduce the remaining amount of transported goods as much as possible.
[0010]
  The present invention has been made in view of such points, and the object of the present invention is to devise a method for predicting the behavior of a conveyed product conveyed by a conveyor belt of a belt conveyor device, and to accurately determine the behavior of the conveyed product. In order to be able to predict, the layout and usage conditions of the belt conveyor device, the specifications of the conveyor belt, and the like can be appropriately set in advance.
[0011]
[Means for Solving the Problems]
  In order to achieve the above object, in the solution means of the present invention, an analysis model for simulating the transport surface of the conveyor belt by a virtual wall surface member is prepared, and the movement state of the conveyor belt is described by moving the wall surface member. However, the behavior of the transported object was predicted.
[0012]
  Specifically, the invention of claim 1 presupposes a behavior prediction method for a belt conveyed product that predicts the behavior of a bulk conveyed product conveyed by a conveyor belt of a belt conveyor device. AndWhen a horizontal rail for locking the conveyed product is provided on the conveying surface of the conveyor belt,At least the conveying surface of the conveyor belt is simulated by a plurality of virtual wall members arranged continuously in the belt longitudinal direction.In addition, at least the front surface of the horizontal rail in the belt traveling direction is simulated by a virtual wall surface member.An analysis model is prepared, and at least material data and geometric data of the conveyor belt and the conveyed product are input to the analytical model, and initial conditions of the behavior of the conveyed product supplied on the conveyor belt are set. Based on the initial condition, the behavior of the transported object that moves with the wall surface member while moving in the longitudinal direction of the belt so as to describe the movement state of the conveyor belt is set. Predict.
[0013]
  In this method, the conveying surface of the conveyor beltOr side pierIn an analysis model that simulates a plurality of virtual wall members, the initial condition of the behavior of the conveyed product supplied on the conveyor belt is set, and the wall member is moved so as to describe the moving state of the conveyor belt By predicting the behavior of the transported object that moves together with the wall surface member based on the initial condition, for example, the distribution of the impact force due to the transported object that falls on the conveyor belt, and the transported object placed on the conveyor belt are transported. It is possible to predict the behavior of the conveyed product or the behavior of the conveyed product discharged from the conveyor belt, and by determining the layout and usage conditions of the belt conveyor device based on this prediction result, the belt conveyor It is possible to achieve both the reliability of the device and the reduction of weight and cost.
[0014]
  Moreover, in general, in the conveyor belt with a cross provided with a horizontal cross as described above, the horizontal cross is easily damaged due to the collision of the fallen transported material on the conveyor belt. It is important to accurately grasp the behavior of the fallen transported object. In addition, if there is a large variation in conveyance, there is a concern that an excessive load will act on only some of the horizontal rails, causing damage to the horizontal rails and spilling of the conveyed items. It is important to accurately grasp the behavior of Therefore, it is particularly effective to predict the behavior of the conveyed product as described above.
[0015]
  Is preferredBy predicting the behavior of transported goods by applying the individual element methodYes, this wayThe effect of the invention of claim 1 can be enhanced by setting the prediction accuracy to be extremely high.
[0016]
  Also,Is preferredA three-dimensional analysis model is prepared, and a plurality of regions in which the moving state of the conveyor belt is different are set in the analysis model, and the wall member moving procedure is defined for each region.
[0017]
  By doing so, the behavior of the invention of claim 1 can be sufficiently obtained by first predicting the behavior of the conveyed object by using a three-dimensional analysis model. In addition, by defining the movement procedure of the wall surface member for each region where the moving state of the conveyor belt is different, when the conveyor belt is in a moving state in which the conveyor belt translates between pulleys or in a moving state in which the belt is wound around the pulley, By moving the wall member to correspond to each movement state, the movement state of the conveyor belt can be accurately described on the analysis model..
[0018]
  next,Claim2InventionIs based on the same premise as the behavior prediction method for a belt conveyed product according to the first aspect of the present invention, and at least the conveying surface of the conveyor belt is simulated by a plurality of virtual wall members continuously arranged in the belt longitudinal direction. ,Simulate multiple conveyor belts including their arrangementAn analysis model is prepared, and at least material data and geometric data of the conveyor belt and the conveyed product are input to the analytical model, and initial conditions of the behavior of the conveyed product supplied on the conveyor belt are set. Set. And predicting, based on the initial conditions, the behavior of the conveyed product moving with the wall member while moving the wall member in the longitudinal direction of the belt so as to describe the movement state of the conveyor belt in the analysis model, andThe transfer behavior of the conveyed product between the plurality of conveyor belts is predicted by applying the individual element method.
[0019]
  In this way, the behavior of the transported material falling from the transport surface of the upstream conveyor belt onto the transport surface of the downstream conveyor belt at the transfer portion of the transported material between the plurality of conveyor belts can be accurately determined by applying the individual element method. To askwear.
[0020]
  Next, the claim3The invention of the present invention is based on a behavior prediction device for a belt conveyed product for predicting the behavior of a conveyed product conveyed by a conveyor belt of a belt conveyor device. An analysis model for simulating at least the transport surface of the conveyor belt by a plurality of virtual wall members arranged continuously in the belt longitudinal direction, and transport supplied to the transport surface of the conveyor belt with respect to the analysis model An initial condition setting means for setting an initial condition of the behavior of the object, and a moving state of the conveyor belt is described by moving the wall member in the analysis model in the belt longitudinal direction based on at least the geometric data of the conveyor belt. Belt movement state description means, the behavior of the conveyed object in the analysis model, at least material data and geometric data of the conveyed object, initial conditions of the conveyed object behavior set by the initial condition setting means, and the belt movement state Conveyed object behavior description described sequentially based on the moving state of the conveyor belt described by the description means A structure and means.
[0021]
  Furthermore, the plurality of conveyor belts including the respective arrangements are simulated, and the initial condition setting means has a predetermined distribution state on the transport surface of the upstream conveyor belt of the plurality of conveyor belts. The transport object generation means for generating and arranging the transport object, the transport object behavior description means, the transfer behavior of the transport object falling from the transport surface of the upstream conveyor belt onto the transport surface of the downstream conveyor belt, Based on the initial distribution state of the conveyance object arranged by the conveyance object generation unit and the movement state of the upstream conveyor belt described by the belt movement state description unit, the structure is described sequentially by the individual element method.
[0022]
  With the above-described configuration, in the analysis model that simulates the conveying surface of the conveyor belt by a plurality of virtual wall members, the initial condition of the behavior of the conveyed product supplied on the conveyor belt is set by the initial condition setting unit, and the belt The wall surface member is moved so as to describe the movement state of the conveyor belt by the movement state description means. Based on the movement state of the wall surface member, that is, based on the movement state of the conveyor belt and the initial condition of the transported object behavior. The behavior of the transported object is described sequentially by the object behavior description means.. ThisThereby, the behavior of the conveyed product on the conveyor line can be predicted in detail, and the same effect as the invention of claim 1 can be obtained.
[0023]
  Also,SaidThe conveyed product is generated and arranged in a predetermined distribution state on the upstream conveyor belt by the conveyed product generation means, and the movement state of the upstream conveyor belt is described by the belt movement state description means, and the movement of the upstream conveyor belt is described. Based on the state and the initial distribution state of the transported object, the transported object behavior description means sequentially describes the behavior of the transported object falling from the upstream conveyor belt onto the transport surface of the downstream conveyor belt by applying the individual element method. Is done. By this, by first predicting the behavior of the conveyed object by applying the individual element method, the prediction accuracy is extremely high,SaidThe effect of this invention can be enhanced. In addition, by describing the behavior of the transported object falling on the downstream conveyor belt by the transported object behavior description means, it is possible to accurately determine the distribution state of the impact force by the dropped transported object.SaidThe effects of the invention can be enhanced.
[0024]
DETAILED DESCRIPTION OF THE INVENTION
  Hereinafter, embodiments of the present invention will be described with reference to the drawings.
[0025]
  (Conveyor line configuration)
  FIG. 2 shows an example of a conveyor line L (belt conveyor device) to which the belt conveyance object behavior prediction method according to the present invention is applied. This conveyor line L is substantially perpendicular to a loose article conveyance object such as coal or ore. A downstream conveyor belt B1 having a so-called S-shaped layout and an upstream conveyor belt B2 for transferring and feeding a conveyed product to the conveyor belt B1 so as to be conveyed upward. The downstream conveyor belt B1 has a base belt 1 (belt main body) made of a flat belt, and the base belt 1 is provided with a drive pulley 2 and a conveyor tail portion provided on a conveyor head (upper end portion in the figure). It is wound around a take-up pulley 3 provided at the lower end of the figure, and the angle is changed by the angle-changing disk pulleys 4, 4, the angle-changing horse back rollers 5, 5,. Has been. When the drive pulley 2 is rotated counterclockwise as shown in the figure by a motor or the like (not shown), the conveyor belt B1 is driven by the drive pulley 2, and the conveyed product is conveyed to the conveying surface 1a (see FIG. 3) and is continuously conveyed.
[0026]
  More specifically, the conveyor belt B1 is provided to extend substantially horizontally between the take-up pulley 3 and the beveled disk pulley 4 at the conveyor tail, and as shown by the arrow f in the figure by the take-up pulley 3. While being applied with tension, it is supported from the lower side by a plurality (five in the illustrated example) of impact rollers 8, 8,... It is received by the conveying surface 1a of the base belt 1. Further, the conveyor belt B1 circulates approximately 90 ° upward around the beveled disk pulley 4 and extends substantially vertically upward therefrom, so that the conveyed object placed on the conveying surface 1a is indicated by an arrow in the figure. To the conveyor head continuously. Reference numerals 10, 10,... Denote a plurality (six in the illustrated example) of carrier rollers arranged at substantially equal intervals in the vertical direction so as to support the inner peripheral side of the conveyor belt B1 in the vertical conveyance state.
[0027]
  Further, in the conveyor head, the conveyor belt B1 extends along a plurality of (six in the illustrated example) back rollers 5, 5,... It circulates and extends substantially horizontally from there to the drive pulley 2, and then turns around the drive pulley 2 in the counterclockwise direction of the figure by approximately 180 degrees. As a result, the conveyed product falls around the drive pulley 2 together with the conveyor belt B1, falls into the discharge chute 11 surrounding the drive pulley 2, and is discharged. On the other hand, the conveyor belt B1 after the fall of the conveyed product extends substantially horizontally from the drive pulley 2 to the beveled disk pulley 4, where the conveyed product adhering to the conveying surface 1a is the beater cleaner. 12 is screened out.
[0028]
  Finally, the conveyor belt B1 rotates about 90 ° around the beveled disk pulley 4 toward the lower side of the figure, and supports the inner peripheral surface by a plurality of (five in the illustrated example) return rollers 13, 13,. However, it extends substantially vertically downward, circulates around the return variable pulley 6 at the conveyor tail, and then extends substantially horizontally from there to reach the take-up pulley 3 again. Reference numerals 14, 14, and 14 in the figure indicate side rollers that support the conveyor belt B 1 from both sides in the belt width direction, and reference numerals 15 and 16 denote transported objects that adhere to the inner peripheral surface of the base belt 1. A scraper to be scraped off is shown. Furthermore, the code | symbol 17 has shown the return guide plate for preventing this conveyed product from scattering, when the conveyed product which remained on the conveyance surface 1a falls on the conveyor return side.
[0029]
  The conveyor belt B1 is provided with a bar that can hold the conveyed product and convey it upward even in a steeply inclined state. Specifically, as shown in FIG. 3, a range of a predetermined width at the center of the outer peripheral surface (upper surface in the drawing) of the base belt 1 is a conveyance surface 1 a, and on both sides of the conveyance surface 1 a in the belt width direction, Ear bars 21 and 21 for preventing the fall are provided over the entire circumference in the belt length direction. The pair of ear rails 21 and 21 protrude substantially perpendicularly from the outer peripheral surface of the base belt 1, and the base belt 5 extends along the beveled disk pulleys 4, 4 and the beveled horse's back rollers 5, 5. It is formed in a bellows shape so that it can expand and contract when it goes around. In addition, contact surfaces 1b and 1b that are in contact with the outer peripheral surfaces of the variable angle disk pulleys 4 and 4 are formed on the outer sides of the left and right ear rails 21 and 21, respectively. The contact surfaces 1b and 1b are wound around the variable angle disk pulleys 4 and 4, respectively.
[0030]
  Further, on the conveying surface 1a of the base belt 1, there are provided substantially rectangular plate-like horizontal bars 22, 22,... Extending in the belt width direction between the ear bars 21, 21 on both sides. The horizontal rails 22, 22,... Are vertically arranged at predetermined intervals over the entire circumference in the belt length direction, and project from the conveying surface 1a of the base belt 1 substantially vertically. And a bent portion extending obliquely from the front end of the vertical wall portion toward the front side in the belt traveling direction. Then, by setting the inclination angle of the bent part of the horizontal beam 22 appropriately, a large number of conveyed objects are locked by the horizontal beam 22 and held so as not to fall out, while being steeply inclined or vertically conveyed. Even so, it can be conveyed upward.
[0031]
  On the other hand, the upstream conveyor belt B2 (hereinafter referred to as a supply belt) is a flat belt, is wound around a drive pulley 18 and a driven pulley (not shown), and extends substantially horizontally between the pulleys. When the drive pulley 18 is rotated counterclockwise in the drawing by a motor or the like (not shown), the supply belt B2 is driven by the drive pulley 18, and the conveyed product on the supply belt B2 is moved to the left side of the drawing. In this case, the sheet is continuously conveyed toward the end of the conveyor and is introduced into the receiving chute 9 from the end of the conveyor. Reference numeral 19 in the figure denotes a guide plate that guides a conveyed product falling from the supply belt B2 toward the receiving chute 9.
[0032]
  By the way, the layout and usage conditions of the conveyor line L as described above, or the specifications of the belts B1 and B2 and the like have been determined mainly depending on the experience of engineers until now. It was difficult to achieve both cost reduction. For example, although it is known that the conveyor belt B1 receives a large impact force due to the fallen conveyance in the receiving chute 9, the distribution state of the impact force and the local maximum value are not accurately known. In order to avoid problems such as premature breakage, it is necessary to limit the height position of the supply belt B2 to be rather low, or to set the specification of the conveyor belt B1 excessively on the safe side, As a result, adverse effects such as high cost and weight increase were great.
[0033]
  In particular, in the case of using the conveyor belt B1 with the beam that can be vertically conveyed as described above, the horizontal beam 22 is greatly damaged by the fallen conveyed material, and thus the horizontal beam 22 itself and the bonding portion with the base belt 1 are excessively reinforced. In many cases, the increase in the weight of the belt due to this has been a big problem. Further, if the amount of the transported object falling from the supply belt B2 varies greatly, an excessive load is applied to some of the horizontal bars 22, 22,... There is also a possibility that the conveyor belt B1 may be damaged due to large deformation and spillage of the conveyed product.
[0034]
  Further, in the discharge chute 11, as described above, while the conveyor belt B1 circulates around the drive pulley 2 on the conveyor head, the conveyed product falls into the discharge chute 11, When there is 22, 22,..., The conveyed product is caught on the horizontal rails 22, 22,. Such a problem becomes particularly prominent when the horizontal bars 22, 22,... Are bent in the middle as in this embodiment, and the residual amount tends to increase as the belt conveyance speed increases. Therefore, in order to reduce the remaining amount of the conveyed product, it is necessary to pay sufficient attention to the setting of the inclination angle of the bent portions of the horizontal rails 22, 22,.
[0035]
  (Prediction of the behavior of the conveyed product)
  In view of such a situation, in this embodiment, in the design stage of the conveyor line L, an analysis model for simulating the conveyor line L and the conveyed product is prepared as shown in FIG. 4, for example. The state of spherical particles corresponding to the object is simulated by the discrete element method. Thereby, the behavior of the conveyed product dropped onto the conveying surface 1a of the conveyor belt B1 via the receiving chute 9 and the behavior until the conveyed product is vertically conveyed by the conveyor belt B1 and discharged by the discharge chute 11. Can be predicted with high accuracy, and by determining the layout, use conditions, and the like of the conveyor line L based on the prediction result, the above-described problems can be solved.
[0036]
  That is, the conveyed product behavior prediction apparatus A of this embodiment includes an analysis model M that simulates the conveyor belt B1, the supply belt B2, and the like as shown in the theoretical block diagram of FIG. M, a belt movement state for describing the movement state of the conveyance object generating means 30 for generating and arranging the conveyance object in a predetermined distribution state on the conveyance surface of the upstream supply belt B2, and the movement states of the supply belt B2 and the conveyor belt B1. Description means 31 is further provided, and further, conveyed object behavior description means 32 is provided which sequentially describes the behavior of the conveyed objects respectively conveyed on the supply belt B2 and the conveyor belt B1 by the individual element method.
[0037]
  Specifically, as shown in FIG. 4, the analysis model M includes a wall surface member w made of a virtual rigid wall on the conveyance surface 1a of the conveyor belt B1, the cross rails 22, 22,... And the conveyance surface of the supply belt B2. .., and the conveyance object guide plate 19 dropped from the supply belt B2 onto the conveyor belt B1 is simulated by a planar member wall made of a virtual rigid wall, and the conveyance object itself has a predetermined elasticity. And a spherical particle ptcl having a damping characteristic. Further, the belt movement state description means 31 moves the wall surface members w, w,... Constituting the supply belt B2 and the conveyor belt B1 in the analysis model M in the belt longitudinal direction based on the belt layout data and the like. The movement states of the supply belt B2 and the conveyor belt B1 are respectively described.
[0038]
  Further, the transported object behavior description means 32 indicates, on the analysis model M, the transfer behavior of the transported object falling from the transport surface of the upstream supply belt B2 onto the transport surface 1a of the downstream conveyor belt B1. Based on at least the material data and geometric data, the initial distribution state of the transported object arranged by the transported object generating means 30, and the moving state of the supply belt B2 described by the belt moving state describing means 31 Describe sequentially by the element method. Further, the transported object behavior description means 32 is a conveyor belt B1 described by the behavior of the dropped transport object described above, at least material data and geometric data of the transport object, and the belt movement state description means 31. The behavior of a series of conveyed objects while the conveyor belt B1 moves from the receiving chute 9 to the discharging chute 11 is sequentially described by the individual element method on the basis of the movement state.
[0039]
  Then, as the output result from the transport object behavior description means 32, the position and speed data of the transport object falling on the conveyor belt B1 in the receiving chute 9 are given as initial conditions for the FEM analysis of the conveyor belt B1, and the receiving chute is received. The moving state of the conveyor belt B1 from 9 to the discharge chute 11 and the behavior of the conveyed product are each output to a display device (image display device) for display.
[0040]
  Here, the individual element method is a kind of calculation method for analyzing the behavior of a discontinuous body in which a large number of granular materials are aggregated. Compared with the calculation method by the continuum approximation, etc., it can be obtained with high accuracy.
[0041]
  Specifically, the case where the element of the granular material is a spherical particle ptcl will be described with reference to FIG. 5. As shown in FIG. 5A, at a certain time t, the translational velocity vA and the rotation at the coordinate position A It is assumed that there is a spherical particle ptcl having a velocity ωA and is in contact with another adjacent spherical particle (shown by a virtual line in the figure) at a point P. In this case, first, in order to obtain the position A ′ and the velocities vA ′ and ωA ′ of the spherical particle ptcl after the minute time Δt has elapsed, the external force F acting on the spherical particle ptcl at that time is the minute time Δt. The coordinate position A after the movement is assumed by solving the general equation of motion of the spherical particle ptcl on the basis of the external force F, the coordinate position A, the translation speed vA, and the rotational speed ωA, assuming that it is constant over time. ', The translation speed vA', and the rotation speed ωA 'are calculated.
[0042]
  Next, after calculating the equation of motion for all the spherical particles in the analysis model, the degree of interference between adjacent spherical particles is examined based on the position coordinates of each spherical particle. That is, as shown in FIG. 4B, when overlap or slip occurs between adjacent spherical particles, in addition to the overlap amount d and slip amount Δs, the elasticity and friction coefficient of each spherical particle In consideration of material characteristics, a force / displacement relational expression based on a general linear spring model is solved to calculate contact forces Fn and Fs in the radial direction and tangential direction at the new contact point p ″, respectively. For example, assuming that a linear spring model with one degree of freedom with attenuation is assumed for each of the radial direction and the tangential direction at the contact point P ″, and the relationship between force and displacement is determined by Hooke's law. Good.
[0043]
  Next, although not shown in the figure for the spherical particle ptcl, the resultant force of the contact forces Fn and Fs acting from all the neighboring spherical particles is obtained, and this resultant force and the gravity acting on the spherical particle ptcl are obtained. This new equation of motion is solved as an external force F in the new equation of motion. That is, as schematically shown in FIG. 6, the position coordinates of the spherical particles and the like are updated based on the calculation result of the equation of motion, the calculation of the force / displacement relational expression based on the update result, and the contact force based on the calculation result By resolving the equation of motion and force / displacement relations every minute time Δt, the behavior of the spherical particle ptcl is changed from moment to moment, such as the calculation of a new equation of motion based on this updated contact force. It can be described accurately.
[0044]
  Hereinafter, the simulation control procedure describing the moving conveyor belts B1 and B2 and the behavior of the conveyed objects on the conveyor belts B1 and B2 by the application of the individual element method will be specifically described based on FIGS. explain.
[0045]
  (Overall procedure)
  First, the overall procedure of the simulation control is as shown in the flowchart of FIG. 7. In step S1 after the start in this main flow, the basic layout of the conveyor line L and the diameter (geometric data) of the spherical particles ptcl. , Data such as density, spring constant, damping rate and friction coefficient (material data), number of spherical particles ptcl, belt moving speed (use conditions), etc. are input. The basic layout of the conveyor line L includes, for example, the dimensions and arrangement of the supply belt B2 and the conveyor belt B1, the distance between the horizontal bars 22 on the conveyor belt B1, and the dimensions and shapes of the horizontal bars 22. Etc.
[0046]
  Subsequently, in step S2, a basic layout of the conveyor belt B1 and the supply belt B2 is generated based on the layout data input in step S1. That is, for example, in the conveyor layout diagram schematically shown in FIG. 9, layout reference points P0 to P5 are determined, and a three-dimensional virtual conveyor line L as shown in FIG. 4 is configured. In this virtual conveyor line L, for the downstream conveyor belt B1, two wall surface members w, w representing the conveying surface 1a of the base belt 1, and both the front and back surfaces of the vertical wall portion and the bent portion of the cross rail 22 are shown. A total of six virtual wall members w, each of which represents a single structural unit, and this structural unit continues from the vicinity of the discharge chute 11 to the vicinity of the receiving chute 9 in the longitudinal direction of the belt. By arranging them, the entire conveyor belt is simulated. Further, for the supply belt B2 on the upstream side, wall members w, w,... Representing the conveyance surface are similarly continuously arranged in the belt longitudinal direction so as to simulate the entire belt.
[0047]
  Subsequently, in step S3, it is determined whether or not the count value i of the control cycle of the simulation control has reached the final value iend. If this determination is YES, the control is terminated. If the determination is NO, the process proceeds to step S4. In this step S4, a predetermined number of spherical particles ptcl are generated and arranged on the conveying surface of the supply belt B2 so as to have, for example, a Gaussian distribution, and in the subsequent step S5, the supply belt B2 is moved (supply conveyor belt movement). In step S6, the conveyor belt B1 is moved (conveyor belt movement). A subroutine for describing the moving state of both belts B1 and B2 will be described later.
[0048]
  Subsequently, in step S7, the behavior of the spherical particles ptcl conveyed with the movement of the supply belt B2 and the conveyor belt B1 is calculated by applying the individual element method as described above. That is, the equation of motion and the force / displacement relational expression regarding the spherical particle ptcl and the wall surface member w during the minute time Δt are solved, and the position coordinates, translation speed, rotation speed, etc. of the spherical particle ptcl are updated. Subsequently, in step S8, the position coordinates of the spherical particles ptcl and the wall surface member w are output to the display device, or the spherical particle ptcl position coordinates are output to the FEM analyzer, and in step S9, the control cycle count value i is output. Is incremented, and then the process returns to step S3.
[0049]
  Then, by repeating the control procedure from step S3 to step S9, the spherical particles ptcl are conveyed to the supply belt B2 and moved in the horizontal direction, and dropped from the supply belt B2 toward the conveyor belt B1. Then, a series of spherical particles ptcl that are placed on the conveying surface 1a of the base belt 1 of the conveyor belt B1 and conveyed substantially vertically upward and fall off the conveyor belt B1 at the discharge chute 11 The position coordinates and speed of can be obtained from moment to moment. That is, the behavior of the spherical particle ptcl group corresponding to the conveyed product in the conveyor line L can be accurately described.
[0050]
  Step S4 of the main flow shown in FIG. 7 corresponds to the conveyed product generating means 30 that generates and arranges spherical particles ptcl (conveyed material) in a predetermined distribution state on the supply belt B2 in the analysis model M.
[0051]
  Steps S5 and S6 sequentially describe the movement states of the supply belt B2 and the conveyor belt B1 in the analysis model M based on the geometric data of the belts B1 and B2 input to the analysis model M, respectively. This corresponds to the belt movement state description means 31.
[0052]
  Further, in step S7, in the analysis model M, at least material data and geometric data of the spherical particles ptcl input to the analysis model M2, and initial spherical particles ptcl arranged on the supply belt B2 by the conveyed product generating means 30 Based on the distribution state and the movement state of the supply belt B2 described by the belt movement state description means 31, the behavior of the spherical particles ptcl falling on the conveyor belt B1 from the supply belt B2 is sequentially described by the individual element method. At the same time, the behavior of the spherical particles ptcl conveyed from the receiving chute 9 to the discharging chute 11 by the conveyor belt B1 is similarly determined by using the behavior of the spherical particles ptcl (falling conveyance object) described in this way as initial conditions. This corresponds to the transport object behavior description means 32 described by the method.
[0053]
  In addition, the step S9 includes the movement states of the conveyor belt B1 and the supply belt B2 described by the belt movement state description unit 31 and the behavior of the spherical particles ptcl described by the conveyed object behavior description unit 32. Each corresponds to image processing means displayed on a display device (not shown).
[0054]
  (Description of belt movement state)
  Next, in steps S5 and S6 of the main flow shown in FIG. 7, a subroutine for describing the movement states of the supply belt B2 and the conveyor belt B1, respectively, based on FIGS. 8 and 9 with reference to FIG. This will be described in detail. In the following description, the left-right direction in FIGS. 4 and 9 is also referred to as the X direction, the direction orthogonal to the drawing sheet as the y direction, and the up and down direction in the figure as the z direction.
[0055]
  First, in step T1 after the start of the flow shown in FIG. 8, the process moves to the head position of the array variable wp representing the coordinate position of the virtual wall member w. Subsequently, in step T2, it is determined whether the value of the array variable wp is empty (wp null). If the determination result is YES, the control is terminated (END). If the determination result is NO, the process proceeds to step T3. Then, the ID indicating whether the wall surface member w constitutes the supply belt B2 or the conveyor belt B1 and the current coordinate value of the wall surface member w are read.
[0056]
  Subsequently, in step T4, based on the ID read in step T3, it is determined whether or not the wall member w constitutes the supply belt B2. If this determination result is NO, although not shown, the conveyor belt B1 While the process proceeds to the control procedure describing the moving state, if the determination result is YES and the supply belt B2, the process proceeds to step T5 and the description of the moving state of the supply belt B2 is started. That is, first, based on the coordinate value read in step T3, it is determined whether the coordinate value of the wall surface member w is below the layout reference point P1 (lower side in FIG. 9; hereinafter the same). . Here, the layout reference point P1 is at a position corresponding to the rotational axis of the drive pulley 18 that drives the supply belt B2, and if the determination result is YES, the process proceeds to step T6 to delete the wall surface member w ( While removing the wall), the process returns to step T2, while if the determination result is NO, the process proceeds to step T7. That is, the supply belt B2 is not described for a position lower than the rotational axis of the drive pulley 18 of the supply belt B2 in the analysis model M2.
[0057]
  On the other hand, in step T7, it is determined whether or not the coordinate value of the wall surface member w is on the left side of the point P1. If the determination result is NO, the process proceeds to step T11. If the determination result is YES, the process proceeds to step T8. Then, a constant angular velocity is given so that the translational velocity of wp in the x direction is once set to zero and rotates around the point P1 in the counterclockwise direction in the figure (rotating around P1). In step T9, the amount of movement of the wall surface member w is calculated based on the constant angular velocity, and in the subsequent step T10, the arrangement based on the coordinate value read in step T3 and the amount of movement calculated in step T9. The coordinate value of the variable wp is updated to a new value, and then the process returns to step T2. That is, in the section (A) on the left side of the point P1, the wall surface member w is rotated at a constant angular velocity about the point P1. As a result, the supply belt B <b> 2 enters a moving state that rotates around the drive pulley 18.
[0058]
  Further, in step T11, which is determined by determining that the wall surface member w is on the right side of the layout reference point P1 in step T7, a constant translation speed in the X direction is given to wp, and in step T12, the wall surface member w Calculate the amount of movement. Subsequently, in step T13, the coordinate value of the wall surface member w is updated to a new value based on the coordinate value read in step T3 and the movement amount calculated in step T13, and then the process returns to step T2. That is, in the section (b) on the right side of the point P1, the wall surface member w is translated in the X direction at a constant translation speed, so that the supply belt B2 is in a translational movement state at a constant speed. If no wall member w is present in (c) in the section on the right side of the point P0, a new wall member w is generated.
[0059]
  That is, for the supply belt B2, the wall surface members w are sequentially arranged from the leading position of the supply belt B2 toward the rear side in the belt traveling direction, and the wall surface members w are respectively moved to the front side in the traveling direction. The moving state of the belt B2 is described. Further, if the leading wall member w of the supply belt B2 is advanced to a position lower than the reference point P1 corresponding to the rotational axis of the drive pulley 18, the wall member w is erased while the last wall member w is generated. When it moves to the front side in the traveling direction from the reference point P0, a new wall member w is generated behind it. In addition, although detailed description is abbreviate | omitted, what is necessary is just to describe using the wall surface members w, w, ... similarly about the movement state of the conveyor belt B1.
[0060]
  (Effect of embodiment)
  Therefore, using the behavior prediction apparatus A of the belt conveyed product according to this embodiment, as described above, the movement state of the conveyor belt B1 and the supply belt B2 by the wall surface members w, w,... While describing, the behavior of the spherical particles ptcl conveyed along with the movement of the wall surface members w, w,... Is simulated by applying the individual element method, so that the bulk material corresponding to the spherical particles ptcl is being conveyed. Can be accurately predicted three-dimensionally.
[0061]
  As a result, for example, the impact force acting on the conveyor belt B1 due to the transported object dropped from the receiving chute 9 can be grasped in a fine and accurate manner including its distribution state. For example, the stress distribution state of the conveyor belt B1 is obtained by FEM analysis, and based on the analysis result, the specifications of the conveyor belt B1, the layout of the height position of the supply belt B2, etc., or the usage conditions such as the maximum transport amount, etc. By determining optimally, the reliability of the conveyor line L and the reduction in weight and cost can be achieved at a high level. Moreover, in the analysis model M, not only the belts B1 and B2, but also the guide plate 19 and the like are simulated, so that it is possible to determine the impact force caused by the fallen conveyed material on the conveyor belt B1 very accurately.
[0062]
  In addition, the behavior of the conveyed product spilling from the horizontal rails 22, 22,... Of the conveyor belt B1 during the vertical conveyance, and the behavior of the conveyed material remaining on the discharge chute 11 by being caught by the horizontal rails 22, 22,. Can also be accurately grasped, so the design of the horizontal beam 22 can be changed so that the amount of spilled or remaining transported material can be reduced, the average transport amount can be limited appropriately, or transport The conveyor layout and use conditions can be appropriately set before the conveyor line L is installed, such as correcting the height of the supply belt B2 and the position of the guide plate 19 so as to reduce the variation in the amount.
[0063]
  Furthermore, in this embodiment, since the movement state of the conveyor belt B1 and the supply belt B2 and the behavior of the conveyed product are displayed on the display device, it has been difficult to see directly after the installation of the conveyor line L until now. The design engineer can visually grasp the flow of the conveyed product of the bulk, intuitively grasp the flow of the conveyed product from the receiving chute 9 to the discharge chute 11, and modify the design partly. An excellent design support effect is obtained in that appropriate measures can be taken quickly.
[0064]
  In addition, in this embodiment, when simulating the conveyor line L, the entire conveyor belt B1 and the supply belt B2 are not described, but directly on the behavior of a conveyed object such as a part of a belt conveying surface or a guide plate. Since only the part that affects the operation is described, the calculation amount of the simulation control can be reduced and the calculation time can be shortened.
[0065]
  In addition, this invention is not limited to the said embodiment, Other various embodiment is included. That is, in the above-described embodiment, the belt conveyance object behavior prediction method of the present invention is applied to the conveyor line L using the conveyor belt with a cross. However, the present invention is not limited thereto, and the present invention includes various methods including those without a horizontal cross. It can be applied to other conveyor belts.
[0066]
  In the embodiment, the analysis model M including the supply belt B2 is prepared, and the behavior of the conveyed product falling from the supply belt B2 to the conveyor belt B1 is described by the individual element method. For example, the spatial distribution and the velocity distribution of the spherical particles ptcl falling on the conveyor belt B1 may be given without performing the simulation of the supply belt B2.
[0067]
【The invention's effect】
  As described above, according to the method for predicting the behavior of the belt conveyed product according to the first aspect of the present invention, an analysis model for simulating the conveying surface of the conveyor belt with a virtual wall surface member is prepared, and the wall surface member is moved to Since the behavior of the conveyed product is predicted while describing the movement state of the belt, the layout and usage conditions of the belt conveyor device, the specifications of the conveyor belt, etc. are preliminarily determined based on the predicted result of the conveyed product behavior. It can be set appropriately. This makes it possible to ensure both the reliability of the belt conveyor device and the reduction in weight and cost.
[0068]
  In particularBy applying to the conveyor belt with a cross,SaidIt is particularly effective to be able to accurately predict the behavior of the conveyed product as described above.
[0069]
  Claim2According to the invention, it is possible to accurately determine the transfer behavior of the conveyed product between the plurality of conveyor belts by applying the individual element method, and thereby the effect of the present invention.ButfurtherBecome validThe
[0070]
  Claim3According to the apparatus for predicting the behavior of a belt conveyed product in the present invention,PiecesThe behavior of the conveyed product can be predicted with very high accuracy by applying the different element method, and the conveyed product that falls on the downstream conveyor belt from the upstream conveyor belt at the transfer portion of the conveyed product between multiple conveyor belts Can be determined precisely and accuratelyThe
[Brief description of the drawings]
FIG. 1 is a theoretical block diagram showing a configuration of a behavior prediction apparatus A for a belt conveyed product according to an embodiment of the present invention.
FIG. 2 is a schematic configuration diagram showing an example of a conveyor line to which the present invention is applied.
FIG. 3 is an enlarged perspective view showing a configuration of a conveyor belt with a crosspiece.
FIG. 4 is a perspective view showing a virtual conveyor line and a conveyed product in an analysis model.
FIG. 5 is a schematic diagram showing movement of element particles in the individual element method.
FIG. 6 is a conceptual diagram of a calculation procedure in an individual element method.
FIG. 7 is a flowchart showing an overall procedure of a method for predicting the behavior of a conveyed product.
FIG. 8 is a flowchart showing a procedure for describing a moving state of a supply belt.
FIG. 9 is a schematic diagram showing a basic layout setting of a conveyor line in an analysis model.
[Explanation of symbols]
A Transported object behavior prediction device
B1 conveyor belt (downstream conveyor belt)
B2 Supply belt (upstream conveyor belt)
L conveyor line (belt conveyor device)
M analysis model
1a Conveyor belt transport surface
22 side pier
30 Transported object generating means (initial condition setting means)
31 Belt movement state description means
32 Means for describing transported object behavior
ptcl Spherical particles (conveyance)
w Wall member

Claims

A method for predicting the behavior of a belt-conveyed object for predicting the behavior of a bulk-conveyed object conveyed by a conveyor belt of a belt conveyor device,
When a horizontal rail for locking the conveyed product is provided on the conveying surface of the conveyor belt,
An analysis model for simulating at least the transport surface of the conveyor belt with a plurality of virtual wall members arranged continuously in the belt longitudinal direction and simulating at least the front surface in the belt traveling direction of the horizontal rail with a virtual wall member Prepare
Input at least the material data and geometric data of the conveyor belt and the conveyed product to the analysis model, and set initial conditions for the behavior of the conveyed product supplied on the conveyor belt,
In the analysis model, the behavior of the transported object moving together with the wall surface member is predicted based on the initial condition while moving the wall surface member in the longitudinal direction of the belt so as to describe the moving state of the conveyor belt. A method for predicting the behavior of belt conveyors.

A method for predicting the behavior of a belt conveyed product for predicting the behavior of a conveyed product conveyed by a conveyor belt of a belt conveyor device,
At least the conveyor surface of the conveyor belt is simulated by a plurality of virtual wall members arranged continuously in the belt longitudinal direction, and an analysis model for simulating a plurality of conveyor belts including each arrangement is prepared,
Input at least the material data and geometric data of the conveyor belt and the conveyed product to the analysis model, and set initial conditions for the behavior of the conveyed product supplied on the conveyor belt,
In the analysis model, while the wall surface member is moved in the longitudinal direction of the belt so as to describe the moving state of the conveyor belt, the behavior of the transported object moving with the wall surface member is predicted based on the initial condition, and the plurality of A method for predicting the behavior of a belt conveyed product, wherein the connection behavior of a conveyed product between conveyor belts is predicted by applying an individual element method.

A belt conveyance object behavior prediction device for predicting the behavior of a bulk material conveyed by a conveyor belt of a belt conveyor device,
An analysis model for simulating at least a plurality of virtual wall members arranged continuously in the belt longitudinal direction on the conveying surface of the conveyor belt;
An initial condition setting means for setting an initial condition of the behavior of the conveyed product supplied on the conveying surface of the conveyor belt with respect to the analysis model;
A belt movement state description means for describing a movement state of the conveyor belt by moving a wall surface member in the analysis model in the longitudinal direction of the belt based on at least geometric data of the conveyor belt;
The behavior of the transported object in the analysis model includes at least material data and geometric data of the transported object, an initial condition of the transported object behavior set by the initial condition setting unit, and a conveyor described by the belt movement state description unit. Conveying object behavior description means sequentially describing based on the movement state of the belt,
The analysis model is a three-dimensional model that simulates a plurality of conveyor belts including their arrangements,
The initial condition setting means includes transported object generating means for generating and arranging a transported object in a predetermined distribution state on the transport surface of the upstream conveyor belt among the plurality of conveyor belts,
The transported object behavior description means indicates the transfer behavior of the transported object falling from the transport surface of the upstream conveyor belt onto the transport surface of the downstream conveyor belt, and the initial distribution state of the transported object arranged by the transported object generating means. And the moving state of the upstream conveyor belt described by the belt moving state description means are configured to sequentially describe by the individual element method.
An apparatus for predicting the behavior of a belt conveyed product.