(57)【要約】【課題】この発明は、粘性のある粘性流動体に混入された複数の棒状体の分布状況を精度よく解析できる棒状体分布解析方法を提供することを目的とする。【解決手段】フレッシュコンクリートに混入された複数の鋼繊維１００の分布状況を解析する鋼繊維分布解析処理であって、ふたつの粒子要素１１同士の相対距離が一定である鋼繊維モデル１０で鋼繊維１００をモデル化する鋼繊維モデル生成工程（ステップｓ２）、フレッシュコンクリートを複数の粘性セルＣｖの集合体としてモデル化するコンクリートモデルＦｃを生成するコンクリートモデル生成工程（ステップｓ２）、コンクリートモデルＦｃの解析領域Ｓにおける粘性セルＣｖの運動を算出する流動体運動算出工程（ステップｓ３）、及び算出された運動を鋼繊維モデル１０に外力として作用させて鋼繊維モデル１０の運動を算出する鋼繊維運動算出工程（ステップｓ５）を行った【選択図】図２
2017191035.tif 000002
【特許請求の範囲】
【請求項１】
粘性のある粘性流動体に混入された複数の棒状体の分布状況を解析する棒状体分布解析方法であって、
棒状に沿って配置された複数の要素体で構成するとともに、隣り合う前記要素体同士の相対距離が一定である棒状体モデルでひとつの前記棒状体をモデル化する棒状体モデル生成工程、
前記粘性流動体を複数の粘性セルの集合体としてモデル化する粘性流動体モデルを生成する粘性流動体モデル生成工程、
前記粘性流動体モデルの解析領域における前記粘性セルの運動を算出する流動体運動算出工程、及び
該流動体運動算出工程で算出された前記運動を前記棒状体モデルに外力として作用させて前記棒状体モデルの運動を算出する棒状体モデル運動算出工程を行う
棒状体分布解析方法。
【請求項２】
複数の前記要素体の合計質量が前記棒状体の質量である
請求項１に記載の棒状体分布解析方法。
【請求項３】
前記棒状体モデルは、
複数の前記要素体と、隣り合う前記要素体同士を連結する無質量の線状連結体とで構成する剛状モデルである
請求項１又は２に記載の棒状体分布解析方法。
【請求項４】
前記流動体運動算出工程において、
有限体積法によって前記運動を前記粘性セルの運動ベクトルとして算出する
請求項１乃至３のうちいずれかに記載の棒状体分布解析方法。
【請求項５】
前記棒状体モデルの条件及び前記粘性流動体モデルの条件のうち少なくとも一方の条件設定する条件設定手段を備えた
請求項１乃至４のうちいずれかに記載の棒状体分布解析方法。
【請求項６】
前記粘性流動体に対する前記棒状体の混入量が２％体積比以下である
請求項１乃至５のうちいずれかに記載の棒状体分布解析方法。
【請求項７】
請求項１乃至６のうちいずれかに記載の棒状体分布解析方法における前記粘性流動体が硬化前のフレッシュコンクリートであるとともに、
前記棒状体が前記フレッシュコンクリートに混入する繊維材である
繊維補強コンクリートの繊維材分布解析方法。
【請求項８】
前記流動体運動算出工程において算出する前記粘性流動体の運動を、
前記フレッシュコンクリートの打ち込み速度を初期速度として算出する
請求項７に記載の繊維補強コンクリートの繊維材分布解析方法。
【請求項９】
粘性のある粘性流動体に混入された複数の棒状体の分布状況を解析する棒状体分布解析装置であって、
棒状に沿って配置された複数の要素体で構成するとともに、隣り合う前記要素体同士の相対距離が一定である棒状体モデルでひとつの前記棒状体をモデル化する棒状体モデル生成手段、
前記粘性流動体を複数の粘性セルの集合体としてモデル化する粘性流動体モデルを生成する粘性流動体モデル生成手段、
前記粘性流動体モデルの解析領域における前記粘性セルの運動を算出する流動体運動算出手段、及び
該流動体運動算出手段で算出された前記運動を前記棒状体モデルに外力として作用させて前記棒状体モデルの運動を算出する棒状体モデル運動算出手段を備えた
棒状体分布解析装置。
【請求項１０】
前記棒状体モデルの条件及び前記粘性流動体モデルの条件のうち少なくとも一方の条件設定する条件設定手段を備えた
請求項９に記載の棒状体分布解析装置。
【請求項１１】
請求項９または１０に記載の棒状体分布解析装置における前記粘性流動体が硬化前のフレッシュコンクリートであるとともに、前記棒状体が前記フレッシュコンクリートに混入する繊維材である
繊維補強コンクリートの繊維材分布解析装置。
【請求項１２】
前記流動体運動算出手段において算出する前記粘性流動体の運動は、
前記フレッシュコンクリートの打ち込み速度を初期速度として算出される
請求項１１に記載の繊維補強コンクリートの繊維材分布解析装置。
【請求項１３】
粘性のある粘性流動体に混入された複数の棒状体の分布状況を解析する棒状体分布解析装置として構成するコンピュータに、
棒状に沿って配置された複数の要素体で構成するとともに、隣り合う前記要素体同士の相対距離が一定である棒状体モデルでひとつの前記棒状体をモデル化する棒状体モデル生成工程、
前記粘性流動体を複数の粘性セルの集合体としてモデル化する粘性流動体モデルを生成する粘性流動体モデル生成工程、
前記粘性流動体モデルの解析領域における前記粘性セルの運動を算出する流動体運動算出工程、及び
該流動体運動算出工程で算出された前記運動を前記棒状体モデルに外力として作用させて前記棒状体モデルの運動を算出する棒状体モデル運動算出工程を実行させる
棒状体分布解析プログラム。
【請求項１４】
請求項１３に記載の棒状体分布解析プログラムにおける前記粘性流動体が硬化前のフレッシュコンクリートであるとともに、前記棒状体が前記フレッシュコンクリートに混入する繊維材である
繊維補強コンクリートの繊維材分布解析プログラム。
【請求項１５】
前記流動体運動算出工程において算出する前記粘性流動体の運動を、
前記フレッシュコンクリートの打ち込み速度を初期速度として算出する
請求項１４に記載の繊維補強コンクリートの繊維材分布解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
例えば、繊維補強コンクリートや繊維強化プラスチック（ＦＲＰ）などフレッシュコンクリートやプラスチックなどの粘性のある粘性流動体に混入された鋼繊維などの複数の棒状体の分布状況を解析する棒状体分布解析方法、繊維補強コンクリートの繊維材分布解析方法、棒状体分布解析装置、繊維補強コンクリートの繊維材分布解析装置、棒状体分布解析プログラム、及び繊維補強コンクリートの繊維材分布解析プログラムに関する。
【背景技術】
【０００２】
従来から、繊維補強コンクリートや繊維強化プラスチック（ＦＲＰ）などフレッシュコンクリートやプラスチックなどに鋼繊維やガラス繊維などを混入して強度を向上させた複合材料が用いられている。例えば、特許文献１には、有機合成繊維が混入された繊維補強樹脂成形物について開示され、特許文献２には、炭素繊維を混入したコンクリートについて開示されている。
【０００３】
これらのように、硬化前の粘性のある粘性流動体に、繊維などの複数の棒状体を混入して強度を向上させた複合材料は、棒状体の分布状況によっては成形品の強度に斑が生じるおそれがあるものの、分布状況を把握することは困難であり、実証実験によって分布状況を推測するしかなく、大きな手間やコストを要していた。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平５−２２０８５３号公報
【特許文献２】特開平８−１４３３５０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
そこで、この発明は、粘性のある粘性流動体に混入された複数の棒状体の分布状況を精度よく解析できる棒状体分布解析方法、繊維補強コンクリートの繊維材分布解析方法、棒状体分布解析装置、繊維補強コンクリートの繊維材分布解析装置、棒状体分布解析プログラム、及び繊維補強コンクリートの繊維材分布解析プログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
この発明は、粘性のある粘性流動体に混入された複数の棒状体の分布状況を解析する棒状体分布解析方法、棒状体分布解析装置及び棒状体分布解析プログラムであって、棒状に沿って配置された複数の要素体で構成するとともに、隣り合う前記要素体同士の相対距離が一定である棒状体モデルでひとつの前記棒状体をモデル化する棒状体モデル生成工程（手段）、前記粘性流動体を複数の粘性セルの集合体としてモデル化する粘性流動体モデルを生成する粘性流動体モデル生成工程（手段）、前記粘性流動体モデルの解析領域における前記粘性セルの運動を算出する流動体運動算出工程（手段）、及び該流動体運動算出工程（手段）で算出された前記運動を前記棒状体モデルに外力として作用させて前記棒状体モデルの運動を算出する棒状体モデル運動算出工程（手段）を行うことを特徴とする。
【０００７】
上述の複数の要素体は、２つ以上の要素体であればよい。  
上述粘性流動体は、硬化前のコンクリートであるフレッシュコンクリート、フレッシュモルタル、フレッシュセメントミルク、硬化前のプラスチックやエポキシ樹脂などの樹脂を含むものとする。
【０００８】
上記棒状体は、鋼繊維、炭素繊維、アラミド繊維などの合成繊維などを含み、いわゆる短繊維とすることができる。  
上記運動は、粘性セルや棒状体モデルの移動、速度、加速度、回転並びに向きの少なくともいずれかとすることができる。
【０００９】
この発明により、粘性のある粘性流動体に混入された複数の棒状体の分布状況を精度よく解析することができる。  
詳述すると、ひとつの棒状体を、棒状に沿って配置された複数の要素体で構成し、隣り合う前記要素体同士の相対距離が一定である棒状体モデルでモデル化するとともに、前記粘性流動体を複数の粘性セルの集合体として粘性流動体モデルでモデル化するため、複雑な挙動を示す粘性流体と棒状体を簡略化したモデルとして生成することができる。
【００１０】
また、前記粘性流動体モデルの解析領域における前記粘性セルの運動を算出するとともに、前記粘性セルの算出された前記運動を前記棒状体モデルに外力として作用させて前記棒状体モデルの運動を算出するため、粘性セルの複雑な挙動を運動として容易に算出し、さらに、粘性セルの算出した運動を棒状体モデルに外力として作用させるため、複雑な挙動を示す粘性セルの運動の影響を考慮した棒状体の運動を算出することで、棒状体の分布状態を精度よく解析することができる。  
このように、本願発明により、複雑な挙動を示す粘性セルの運動を複雑化させずに棒状体モデルに作用させることで、粘性セルの影響を考慮した棒状体モデルの運動を容易に算出することができる。
【００１１】
この発明の態様として、複数の前記要素体の合計質量を前記棒状体の質量とすることができる。  
この発明により、棒状体の質量の影響を考慮し、現実に即したより忠実な解析結果を得ることができる。殊に、複数の前記要素体が同質量である場合はさらに忠実な解析結果を得ることができる。
【００１２】
またこの発明の態様として、前記棒状体モデルを、複数の前記要素体と、隣り合う前記要素体同士を連結する無質量の線状連結体とで構成する剛状モデルとすることができる。  
この発明により、長さのある棒状体を、変形等の複雑化する要素を軽減し、解析上簡略化して適切にモデル化することができる。
【００１３】
またこの発明の態様として、前記流動体運動算出工程（手段）において、有限体積法によって前記運動を前記粘性セルの運動ベクトルとして算出することができる。  
この発明により、簡易な解析プロセスで、棒状体モデルの運動を適切に評価することができる。
【００１４】
またこの発明の態様として、前記棒状体モデルの条件及び前記粘性流動体モデルの条件のうち少なくとも一方の条件設定する条件設定手段を備えることができる。  
上述の棒状体モデルの条件は、棒状体モデルの長さ、質量、径、形状、硬度、可撓性あるいは摩擦係数など棒状体の性状に関する条件や、棒状体モデルの粘性流動体モデルに対する混入量や混入タイミング、混入する位置など棒状体の粘性流動体への混入に関する条件とすることができる。
【００１５】
上述の前記粘性流動体モデルの条件は、粘性流動体モデルの量、粘性、密度（比重）、スランプ、ワーカビリティ、フローあるいは自己充填性などの粘性流動体の性状に関する条件や粘性流動体で構成する形状や大きさ、さらには粘性流動体を打ち込む打ち込み箇所の位置、大きさ、形状あるいは高さなどの粘性流動体を打ち込む際の条件などとすることができる。
【００１６】
この発明により、粘性流動体に対して棒状体が適切に分布する条件を設定することができる。  
詳しくは、前記棒状体モデルの条件及び前記粘性流動体モデルの条件のうち少なくとも一方を条件設定手段で設定することで、設定された条件での棒状体の分布状況を精度よく解析することができるため、適切な分布状況となる条件をシミュレーションして、より良い分布状況となる条件を設定することができる。
【００１７】
またこの発明の態様として、前記粘性流動体に対する前記棒状体の混入量を２％体積比以下とすることができる。  
この発明により、混入する棒状体モデルが粘性体モデルの運動に影響しないため、上述の説明における解析においてより精度のよい解析結果を得ることができる。
【００１８】
またこの発明は、上述の棒状体分布解析方法、棒状体分布解析装置及び棒状体分布解析プログラムにおける前記粘性流動体が硬化前のフレッシュコンクリートであるとともに、前記棒状体が前記フレッシュコンクリートに混入する繊維材である繊維補強コンクリートの繊維材分布解析方法、繊維補強コンクリートの繊維材分布解析装置、及び繊維補強コンクリートの繊維材分布解析プログラムとすることを特徴とする。  
この発明により、繊維補強コンクリートにおける繊維材の分布を精度よく解析することができる。
【００１９】
この発明の態様として、前記流動体運動算出工程（手段）において算出する前記粘性流動体の運動を、前記フレッシュコンクリートの打ち込み速度を初期速度として算出することができる。  
この発明により、流動するフレッシュコンクリートの運動を適切に算出して、繊維材の分布をさらに精度よく解析することができる。
【発明の効果】
【００２０】
この発明によれば、粘性のある粘性流動体に混入された複数の棒状体の分布状況を精度よく解析できる棒状体分布解析方法、繊維補強コンクリートの繊維材分布解析方法、棒状体分布解析装置、繊維補強コンクリートの繊維材分布解析装置、棒状体分布解析プログラム、及び繊維補強コンクリートの繊維材分布解析プログラムを提供することができる。
【図面の簡単な説明】
【００２１】
【図１】繊維材分布解析システムの概要説明図。
【図２】繊維材分布解析方法のフロー図。
【図３】繊維運動解析方法のフロー図。
【図４】繊維導入解析方法のフロー図。
【図５】解析条件入力画面の概略図。
【図６】条件設定部の説明図。
【図７】鋼繊維のモデル化の説明図。
【図８】領域モデルにおける流れ場ベクトルの概略図。
【図９】領域モデルにおける流れ場ベクトルの概略図。
【図１０】領域モデルにおける流れ場ベクトルの概略図。
【図１１】繊維要素に流れ場ベクトルを作用させることについての説明図。
【図１２】領域モデルにおける繊維分布状況の概略出力図。
【図１３】領域モデルにおける繊維分布状況の概略出力図。
【図１４】領域モデルにおける繊維分布状況の概略出力図。
【発明を実施するための形態】
【００２２】
この発明の一実施形態を以下図面とともに説明する。  
図１は繊維材分布解析システム１の概要説明図を示している。詳しくは、図１（ａ）は繊維材分布解析システム１の概略図を示し、図１（ｂ）は繊維材分布解析システム１を構成する解析用コンピュータ２（以下において解析ＰＣ２という）の概略ブロック図を示している。
【００２３】
また、図２は繊維材分布解析方法のフロー図を示し、図３は、繊維材分布解析方法におけるステップｓ５に対応する繊維運動解析方法のフロー図を示し、図４は、繊維運動解析方法におけるステップｔ４に対応する繊維導入解析方法のフロー図を示している。
【００２４】
また、図５は表示画面３００の概略図を示し、図６は条件設定部３１０の説明図を示している。具体的には、図６（ａ）は打ち込み条件を入力する際の条件設定部３１０を示し、図６（ｂ）はフレッシュコンクリートの性状を入力する際の条件設定部３１０を示し、図６（ｃ）は鋼繊維の性状を入力する際の条件設定部３１０を示している。なお、図５における表示画面３００の条件設定部３１０は解析領域Ｓのサイズを設定する画面を示している。
【００２５】
図７は鋼繊維１００のモデル化の説明図を示し、詳しくは、図７（ａ）は、鋼繊維１００の拡大斜視図を示し、図７（ｂ）は鋼繊維１００をモデル化した鋼繊維モデル１０の拡大斜視図を示している。
【００２６】
図８乃至図１０は解析領域ＳにおけるコンクリートモデルＦｃの概略図を示し、図１１は鋼繊維モデル１０の繊維要素１１に流れ場ベクトルＶｒを作用させることについての説明図を示し、図１２乃至図１４は解析領域Ｓにおける鋼繊維モデル１０の分布状況の概略出力図を示している。
【００２７】
なお、図１１（ａ）はコンクリートモデルＦｃに鋼繊維モデル１０が導入された状態の初期状態の概略斜視図を示し、図１１（ｂ）はコンクリートモデルＦｃにおいて、流れ場ベクトルＶｒが作用して移動した鋼繊維モデル１０の概略斜視図を示している。
【００２８】
また、図１１（ｃ）は初期状態の鋼繊維モデル１０及びコンクリートモデルＦｃの平面図を示し、図１１（ｄ）は同状態の正面図を示し、図１１（ｅ）は同状態の右側面図を示している。図１１（ｆ）は流れ場ベクトルＶｒが作用して移動した鋼繊維モデル１０の平面図を示し、図１１（ｇ）は同状態の正面図を示し、図１１（ｈ）は同状態の右側面図を示している。  
図１２乃至図１４は解析領域Ｓにおける鋼繊維モデル１０の分布状況のモデル表示部３２０への概略出力図を示している。
【００２９】
フレッシュコンクリートに鋼繊維１００を混入した繊維補強コンクリートにおける鋼繊維１００に分布状況を解析する鋼繊維分布解析装置である繊維材分布解析システム１は、鋼繊維分布解析プログラムを実行するコンピュータであり、解析ＰＣ２、液晶モニタ又はＣＲＴディスプレイ等で構成して各種情報を表示する表示部３、マウス４ｂやキーボード４ａ等の操作入力部４で構成している。
【００３０】
解析ＰＣ２は、ＣＰＵ、ＲＯＭ、及びＲＡＭで構成する制御部２１とともに、内部に備えた記憶部２２で構成し、記憶部２２に格納した、プログラムに従って各種処理を実行する装置である。
【００３１】
なお、制御部２１は、流れ場ベクトルＶｒを算出する流れ場ベクトル算出部２１１と、鋼繊維１００の運動を算出する繊維運動算出部２１２などの算出手段を有する演算処理手段である。  
記憶部２２には、繊維材分布解析プログラムなどの各種プログラムや、各種装置を制御する制御プログラムに加え、各種データ等のステップｓ１で入力された条件入力データを記憶する設定条件記憶部２２１、ステップｓ３で算出された流れ場ベクトルデータを記憶する流れ場ベクトル記憶部２２２、及びステップｓ５で算出された鋼繊維モデル１０の運動を記憶する繊維運動記憶部２２３などを記憶している。
【００３２】
表示部３に表示される表示画面３００は、図５に示すように、解析条件を入力する条件設定部３１０と、コンクリートモデルＦｃを表示するモデル表示部３２０とで構成している。  
条件設定部３１０は、設定対象選択部３１１と、条件項目入力部３１２とで構成され、設定対象選択部３１１で選択され設定対象に応じた条件項目が条件項目入力部３１２に表示され、表示された条件項目を入力するように構成されている。
【００３３】
具体的には、設定対象選択部３１１は、解析領域Ｓのサイズを設定するための解析モデル設定ボタン３１１ａ、フレッシュコンクリートを打ち込む条件を設定する打ち込み条件ボタン３１１ｂ、フレッシュコンクリートの性状を設定するコンクリートボタン３１１ｃ、及び鋼繊維１００の性状を設定する鋼繊維ボタン３１１ｄで構成されている。
【００３４】
解析モデル設定ボタン３１１ａが押下された場合、条件項目入力部３１２には、図５に示すように、解析領域Ｓのサイズを示す“長さ”、“幅”及び“高さ”がそれぞれ入力されるように構成され、打ち込み条件ボタン３１１ｂが押下された場合、図６（ａ）に示すように、解析領域Ｓへのフレッシュコンクリートの打ち込み条件として、打ち込み口の中心位置を設定する“打ち込み口中心位置Ｘ座標”、“打ち込み口中心位置Ｙ座標”及び“打ち込み口中心位置Ｚ座標”並びに、打ち込み口のサイズを設定する“打ち込み口長さ”及び“打ち込み口幅”、さらには打ち込み条件となる“打ち込み速度”、“打ち込み時間”及び解析の際の時間ステップを設定する“時間ステップ”がそれぞれ入力されるように構成されている。
【００３５】
コンクリートボタン３１１ｃが押下された場合、図６（ｂ）に示すように、フレッシュコンクリートの性状として、密度及び動粘性性係数がそれぞれ入力されるように構成され、鋼繊維ボタン３１１ｄが押下された場合、図６（ｃ）に示すように、鋼繊維１００の性状として、密度、質量、長さ、直径及び長さ／直径比がそれぞれ入力されるように構成されている。  
コンクリートモデルＦｃを表示するモデル表示部３２０は、条件設定部３１０で設定された条件に応じたコンクリートモデルＦｃを表示するように構成されている。
【００３６】
このように構成した繊維材分布解析システム１を用い、記憶部２２に記憶した繊維材分布解析プログラムを実行して行う繊維材分布解析処理について、次に説明する。  
まず、繊維材分布解析処理を行うにあたって、実行する繊維材分布解析プログラムの入力画面で解析条件を入力する（ステップｓ１）。
【００３７】
具体的には、表示部３に表示される表示画面３００の条件設定部３１０に解析条件を入力して設定する。  
なお、上述の入力条件により、鋼繊維補強コンクリートを構成するフレッシュコンクリーに対して鋼繊維１００の混入量は２％体積比以下となるように設定されている。
【００３８】
続いて、条件入力された繊維材分布解析プログラムは、入力条件に基づいて解析領域Ｓ、コンクリートモデルＦｃ及び鋼繊維モデル１０を生成する（ステップｓ２）。  
解析領域Ｓは、粘性流動体であるフレッシュコンクリートを打設して製造するコンクリート加工体のサイズであり、本実施形態では、図８に示すように、正面視略正方形断面で奥行き方向Ｄに長い直方体形状としている。なお、奥行き方向Ｄに長い解析領域Ｓにおける奥側に打ち込み口ｘを設けている。
【００３９】
そして、コンクリートモデルＦｃは、解析領域Ｓに対して、入力条件に応じて、打設したフレッシュコンクリートを、打設開始から打設完了までを所定時間ステップごとに各時間ごとにモデル化している。また、コンクリートモデルＦｃは、粘性流動体であるフレッシュコンクリートを複数の粘性セルＣｖの集合体としてモデル化している。
【００４０】
鋼繊維１００をモデル化した鋼繊維モデル１０は、図７に示すように、２つの粒子要素１１と、粒子要素１１同士を連結する連結体１２とで構成している。  
粒子要素１１は、鋼繊維１００の質量の半分の質量を有する要素であり、図７（ｂ）では、球状に図示しているが、大きさ及び形状の概念を有さないモデルである。
【００４１】
連結体１２は、ふたつの粒子要素１１同士の相対間距離Ｌ（以下において要素間距離Ｌという）を、鋼繊維１００の長さ（繊維長Ｌ０）に合わせて拘束するものであり、図７（ｂ）では棒状に図示しているが、太さや形状の概念を有さず、無質量のモデルである。なお、連結体１２は、ふたつの粒子要素１１同士の相対間距離が規制できれば、連結体１２自体をモデル化せず、相対間距離が拘束されたふたつの粒子要素１１のみで鋼繊維モデル１０をモデル化してもよい。
【００４２】
このように構成された鋼繊維モデル１０は、ふたつの粒子要素１１の相対方向は変化できるものの、連結体１２によってその相対間距離を一定に保つように構成されている。つまり、鋼繊維モデル１０は、同じ質量を有するふたつの粒子要素１１と、ふたつの粒子要素１１の相対距離を拘束する連結体１２とで剛状モデルに構成されている。
【００４３】
鋼繊維モデル１０及びコンクリートモデルＦｃのモデル化が完了すると、図７乃至図９に示すように、打ち込み口ｘから解析領域Ｓに対して解析条件で設定したフレッシュコンクリートの打ち込み条件での打設状況を解析にて再現してモデル表示部３２０に表示し、時間ステップごとにコンクリートモデルＦｃの流れ場ベクトルＶｒを算出するとともに（ステップｓ３）、算出された流れ場ベクトルＶｒを記憶部２２に記憶する（ステップｓ４）。なお、図７乃至図９では、時間ステップごとに刻々と変化する打設状況における一部のみを図示している。
【００４４】
ここでいう流れ場ベクトルＶｒは、粘性流動体であるフレッシュコンクリートを複数の粘性セルＣｖの集合体として捉え、各粘性セルＣｖの運動ベクトル、つまり、粘性セルＣｖの移動速度及び移動方向を表している。
【００４５】
なお、このとき、鋼繊維１００は混入されていないフレッシュコンクリート単体での打設状況を再現することとなるが、鋼繊維１００の混入量をフレッシュコンクリートに対する２％体積比以下となるように設定しているため、流れ場ベクトルＶｒは、鋼繊維１００が混入された場合と変わりなく算出することができる。
【００４６】
具体的には、打ち込み速度及び打ち込み高さなどの入力条件に応じたコンクリートモデルＦｃに対して、解析領域Ｓ内の質量や運動量が保存されるように離散化を行う、有限体積法（finite volume method、以下においてＦＶＭ法という）によって、各粘性セルＣｖの流れ場ベクトルＶｒを時間ステップごとに算出する。
【００４７】
なお、フレッシュコンクリートをモデル化したコンクリートモデルＦｃの流動は、非ニュートン粘性流体（Non-Newtonian viscous fluid）でBingham流体モデルより、以下の式のように表現することができ、これらの式より、上述のＦＶＭ法によって、各粘性セルＣｖの流れ場ベクトルＶｒを時間ステップごとに算出する。
【００４８】
【数１】
2017191035.tif 000003
この数式において、t：時間，ρ：流体（コンクリート粒子）の密度，U：流体の速度ベクトル，p：流体の圧力，SM：外部の運動源による運動量（external momentum sources），μ：流体の動的粘性係数を示している。
【００４９】
次に、ＦＶＭ法によって算出された、時間ステップごとの各粘性セルＣｖの流れ場ベクトルＶｒを、鋼繊維モデル１０に対して外力として作用させて、鋼繊維モデル１０の運動を算出する（ステップｓ５）。
【００５０】
具体的には、図１１（ａ），（ｃ）乃至（ｅ）に示すように、複数の粘性セルＣｖで構成される一部のコンクリートモデルＦｃにおける鋼繊維モデル１０の粒子要素１１に対して、それぞれの粒子要素１１が位置する粘性セルＣｖの流れ場ベクトルＶｒが各粒子要素１１に外力として作用することとなる。
【００５１】
その結果、流れ場ベクトルＶｒが作用した粒子要素１１は、図１１（ｂ），（ｆ）乃至（ｈ）に示すように、各粒子要素１１は、連結体１２に要素間距離Ｌが拘束された状態で、移動することとなる。
【００５２】
なお、本解析では、図３に示すように、まずは、時刻Ｔｎにおける鋼繊維モデル１０を、ランダム関数を用いてランダムに生成する（ステップｔ１）。ランダム関数については一般的なランダム関数を用いることができるため、詳細な説明は省略する。
【００５３】
鋼繊維モデル１０のランダム生成方法は、図４に示すように、解析領域Ｓ内における点Ａの生成（ステップｕ１）、及び単位ベクトルＲの生成（ステップｕ２）を行い、生成された点Ａと単位ベクトルＲ及び連結体１２に対応する繊維長Ｌに基づいて点Ｂを生成する（ステップｕ３）。
【００５４】
生成された点Ｂが解析領域Ｓ内でない場合（ステップｕ４：Ｎｏ）はステップｕ２に戻り、解析領域Ｓ内である場合（ステップｕ４：Ｙｅｓ）は鋼繊維モデル１０の生成を完了し、生成された鋼繊維モデル１０を用いる。これを、時刻Ｔｎにおける鋼繊維モデル１０の本数分実行し、所定数の鋼繊維モデル１０をランダムに生成する。
【００５５】
次に、記憶部２２に記憶された時刻Ｔｎにおける流れ場ベクトルＶｒを読込み、ステップｔ１で生成した鋼繊維モデル１０に対して、読み込んだ流れ場ベクトルＶｒを外力として作用させて、鋼繊維モデル１０の運動を算出する（ステップｔ３）。
【００５６】
詳しくは、鋼繊維モデル１０の両端の粒子要素１１の受けた力と運動量との関係式は下記のとおりとなる。  
【００５７】
【数２】
2017191035.tif 000004
【００５８】
粒子の受けた流体抗力FDは次式によってあらわすことができる。  
【００５９】
【数３】
2017191035.tif 000005
u：流体（コンクリート）の速度，ν：粒子の速度，
τp：粒子の速度応答の時刻暦（velocity  response  time）で次式による。  
【００６０】
【数４】
2017191035.tif 000006
ここで、ρp：粒子の密度  dp：粒子の直径  CD：流体抗力係数  RCr：粒子の相対レイノルズ数である。  
【００６１】
【数５】
2017191035.tif 000007
【００６２】
また、A（Sp)、B（Sp)、C（Sp)、D（Sp)は下記の経験式により求められる。
【００６３】
【数６】
2017191035.tif 000008
【００６４】
このようにして算出された粒子要素１１の運動に基づき、連結体１２による拘束条件で粒子要素１１の運動を修正する（ステップｔ４）。  
詳しくは、上述の運動方程式に基づき、T＝Tn時刻における両粒子要素１１に作用した受けた力と加速度を算出するとともに、Tn+1=Tn＋ΔTの位置と速度を算出するとともに、両粒子要素１１の両粒子間の中心位置を算出する。
【００６５】
そして、両粒子要素１１の中心位置及び両粒子からなる方向を固定とし、要素間距離Ｌが鋼繊維１００の繊維長Ｌ０と等しくなるように調整する。すなわち、
要素間距離Ｌ＞繊維長Ｌ０の場合は、両粒子要素１１を中心位置に向かって（要素間距離Ｌ−繊維長Ｌ０）／２に応じた距離を移動させる。  
逆に、要素間距離Ｌ＜繊維長Ｌ０の場合は、両粒子要素１１を中心位置に離れて（繊維長Ｌ０―要素間距離Ｌ）／２に応じた距離を移動させる。
【００６６】
このようにして、両粒子要素１１の新しい位置をTn+1時刻の空間位置として、図１２乃至図１４に示すように、コンクリートモデルＦｃに対して各鋼繊維モデル１０が表示されたモデルを表示画面３００のモデル表示部３２０に表示するとともに記憶部２２に記憶する（ステップｔ５）。
【００６７】
なお、図１２乃至図１４は、時間ステップごとに刻々と変化する打設状況における一部のみを図示した図７乃至図９に対応して、鋼繊維モデル１０が表示されたコンクリートモデルＦｃを図示している。
【００６８】
これを、最終ステップまで繰り返すが（ステップｔ６：Ｎｏ）、算出された空間位置を、次のTn+2=Tn+1＋ΔTにおける粒子要素１１の初期条件として、上述の算定を最終ステップまで繰り返す（ステップｔ６：Ｙｅｓ）。繊維両端粒子の位置と速度を算出し、最終の時間ステップまで続く。
【００６９】
このようにして、最終ステップまで繰り返した粒子要素１１の位置情報、つまり算出結果を記憶部２２に記憶するとともに、上述のコンクリートモデルＦｃにおける制御部２１０に表示して、本鋼繊維分布解析を終了する。
【００７０】
このように、粘性のあるフレッシュコンクリートに混入された複数の鋼繊維１００の分布状況を解析する鋼繊維分布解析処理、鋼繊維分布解析装置及び鋼繊維分布解析プログラムは、棒状に沿って配置されたふたつの粒子要素１１で構成するとともに、隣り合う粒子要素１１同士の相対距離が一定である鋼繊維モデル１０でひとつの鋼繊維１００をモデル化する鋼繊維モデル生成工程（ステップｓ２）、フレッシュコンクリートを複数の粘性セルＣｖの集合体としてモデル化するコンクリートモデルＦｃを生成するコンクリートモデル生成工程（ステップｓ２）、コンクリートモデルＦｃの解析領域における粘性セルＣｖの運動を算出する流動体運動算出工程（ステップｓ３）、及び流動体運動算出工程（ステップｓ３）で算出された運動を鋼繊維モデル１０に外力として作用させて鋼繊維モデル１０の運動を算出する鋼繊維運動算出工程（ステップｓ５）を行うため、粘性のあるフレッシュコンクリートに混入された複数の鋼繊維１００の分布状況を精度よく解析することができる。
【００７１】
詳述すると、ひとつの鋼繊維１００を、ふたつの粒子要素１１で構成するとともに、隣り合う粒子要素１１同士の相対距離が一定である鋼繊維モデル１０でモデル化するとともに、フレッシュコンクリートを複数の粘性セルＣｖの集合体としてコンクリートモデルＦｃでモデル化するため、複雑な挙動を示す粘性流体と鋼繊維１００を簡略化したモデルとして生成することができる。
【００７２】
また、コンクリートモデルＦｃの解析領域における粘性セルＣｖの運動（流れ場ベクトルＶｒ）を算出するとともに、粘性セルＣｖの算出された運動（流れ場ベクトルＶｒ）を外力として作用させて鋼繊維モデル１０の運動を算出するため、粘性セルＣｖの複雑な挙動を運動（流れ場ベクトルＶｒ）として容易に算出し、さらに、粘性セルＣｖの算出した運動（流れ場ベクトルＶｒ）を鋼繊維モデル１０に外力として作用させるため、複雑な挙動を示す粘性セルＣｖの運動（流れ場ベクトルＶｒ）の影響を考慮した鋼繊維モデル１０の運動を算出することで、鋼繊維１００の分布状態を精度よく解析することができる。
【００７３】
このように、本願発明により、複雑な挙動を示す粘性セルＣｖの運動（流れ場ベクトルＶｒ）を複雑化させずに鋼繊維モデル１０に作用させることで、粘性セルＣｖの影響を考慮した鋼繊維モデル１０の運動を容易に算出することができる。
【００７４】
また、ふたつの粒子要素１１の合計質量を鋼繊維１００の質量とすることにより、鋼繊維１００の質量の影響を考慮し、現実に即したより忠実な解析結果を得ることができる。殊に、ふたつの粒子要素１１が同質量である場合はさらに忠実な解析結果を得ることができる。
【００７５】
また、鋼繊維モデル１０を、ふたつの粒子要素１１と、隣り合う粒子要素１１同士を連結する無質量の連結体１２とで構成する剛状モデルとすることにより、長さのある鋼繊維１００を、変形等の複雑化する要素を軽減し、解析上簡略化して適切にモデル化することができる。
【００７６】
また、流動体運動算出工程（ステップｓ３）において、有限体積法（ＦＶＭ法）によって運動を粘性セルＣｖの流れ場ベクトルＶｒとして算出することにより、簡易な解析プロセスで、鋼繊維モデル１０の運動を適切に評価することができる。
【００７７】
また、鋼繊維モデル１０の条件及びコンクリートモデルＦｃの条件を設定する条件設定部３１０を備えているため、粘性流動体に対して鋼繊維１００が適切に分布する条件を設定することができる。  
詳しくは、鋼繊維モデル１０の条件及びコンクリートモデルＦｃの条件を条件設定部３１０で設定することで、設定された条件での鋼繊維１００の分布状況を精度よく解析することができるため、適切な分布状況となる条件をシミュレーションして、より良い分布状況となる条件を設定することができる。
【００７８】
また、フレッシュコンクリートに対する鋼繊維１００の混入量を２％体積比以下としているため、混入する鋼繊維モデル１０が粘性体モデルの運動に影響せず、上述の説明における解析においてより精度のよい解析結果を得ることができる。
【００７９】
この発明の構成と、上述の実施例との対応において、この発明の粘性流動体はフレッシュコンクリートに対応し、
以下同様に、
棒状体は鋼繊維１００に対応し、
要素体は粒子要素１１に対応し、
棒状体モデルは鋼繊維モデル１０に対応し、
棒状体モデル生成工程はステップｓ２に対応し、
粘性流動体モデルはコンクリートモデルＦｃに対応し、
粘性流動体モデル生成工程はステップｓ２に対応し、
流動体運動算出工程はステップｓ３に対応し、
棒状体モデル運動算出工程はステップｓ５に対応し、
棒状体分布解析方法は繊維材分布解析処理に対応し、
線状連結体は連結体１２に対応し、
条件設定手段は、条件設定部３１０に対応し、
運動ベクトルは流れ場ベクトルＶｒに対応し、
繊維材は鋼繊維１００に対応し、
繊維補強コンクリートの繊維材分布解析方法は繊維材分布解析処理に対応し、
棒状体分布解析装置は繊維材分布解析システム１に対応し、
棒状体モデル生成手段、及び粘性流動体モデル生成手段はステップｓ２を実行する制御部２１に対応し、
流動体運動算出手段は流れ場ベクトル算出部２１１に対応し、
棒状体モデル運動算出手段は繊維運動算出部２１２に対応し、
コンピュータは繊維材分布解析システム１に対応し、
棒状体分布解析プログラムは、繊維材分布解析プログラムに対応する、
この発明は、上述の実施形態の構成のみに限定されるものではなく、多くの実施の形態を得ることができる。
【００８０】
例えば、上述の説明では、鋼繊維１００のモデルとして、ふたつの粒子要素１１を連結体１２で連結して剛状モデルとして鋼繊維モデル１０を構成したが、３つ以上の粒子要素１１を連結体１２で連結して鋼繊維モデル１０を構成してもよい。繊維補強コンクリートに混入する繊維として鋼繊維１００をモデル化して鋼繊維モデル１０を用いたが、鋼繊維１００の他、炭素繊維、アラミド繊維などの合成繊維などをモデル化して用いてもよく、さらには、屈曲可能な繊維をモデル化するために連結体１２に接点を設けて、屈曲可能な鋼繊維モデル１０としてモデル化してもよい。
【００８１】
また、繊維補強コンクリートの鋼繊維の配向を解析するためにフレッシュコンクリートをモデル化したコンクリートモデルＦｃを用いて解析したが、フレッシュコンクリートのみならず、フレッシュモルタル、フレッシュセメントミルク、硬化前のプラスチックやエポキシ樹脂などの樹脂をモデル化して解析してもよい。
【００８２】
また、鋼繊維モデル１０の条件として、鋼繊維モデル１０の長さ、質量のみならず、径、形状、硬度、可撓性あるいは摩擦係数など鋼繊維１００の性状に関する条件や、鋼繊維モデル１０のコンクリートモデルＦｃに対する混入タイミングや混入する位置などを条件設定部
３１０で設定できるように構成してもよく、コンクリートモデルＦｃの条件として、コンクリートモデルＦｃのサイズのみならず、粘性、密度（比重）、スランプ、ワーカビリティ、フローあるいは自己充填性などのフレッシュコンクリートの性状に関する条件や、フレッシュコンクリートを打ち込む打ち込み箇所の位置のみならず、大きさ、形状あるいは高さなどのフレッシュコンクリートを打ち込む際の条件などを条件設定部３１０で設定できるように構成することで、鋼繊維モデル１０の分布状況を高精度で解析することができる。
【００８３】
さらに、上述の説明では、フレッシュコンクリートに鋼繊維１００を混入し、直体状の構造物を構築するために直方体状の解析領域Ｓを用いたが、床版状の構造物を構築数ために、版状の解析領域Ｓを用いて解析を行ってもよい。
【００８４】
また、上述の説明では、有限体積法（ＦＶＭ法）によって粘性セルＣｖの流れ場ベクトルＶｒを算出したが、有限要素法（ＦＥＭ法）などの別の解析方法で粘性セルＣｖの流れ場ベクトルＶｒを算出してもよい。
【符号の説明】
【００８５】
１…繊維材分布解析システム
１０…鋼繊維モデル
１１…粒子要素
１２…連結体
２１…制御部
１００…鋼繊維
３１０…条件設定部
Ｃｖ…粘性セル
Ｆｃ…コンクリートモデル
Ｓ…解析領域
Ｖｒ…流れ場ベクトル
【図１】
2017191035.tif 000009
【図２】
2017191035.tif 000010
【図３】
2017191035.tif 000011
【図４】
2017191035.tif 000012
【図５】
2017191035.tif 000013
【図６】
2017191035.tif 000014
【図７】
2017191035.tif 000015
【図８】
2017191035.tif 000016
【図９】
2017191035.tif 000017
【図１０】
2017191035.tif 000018
【図１１】
2017191035.tif 000019
【図１２】
2017191035.tif 000020
【図１３】
2017191035.tif 000021
【図１４】
2017191035.tif 000022