"粉粒体の粒度分布fを求め、この粒度分布fから下記a式により当該粉粒体の充填率pを求めることを特徴とする粉粒体の充填率の算出方法。粒度分布f1r,f2r,...をもつ2種以上の粉粒体をそれぞれ体積V1:V2:...で混合した混合粉粒体について、下記b式で与えられる混合粉粒体の度数Friを用いてa式により当該混合粉粒体の充填率pを求めることを特徴とする粉粒体の充填率の算出方法。ここでPj:粒度分布fjrをもつ粉粒体の充填率Vj:粒度分布fjrをもつ粉粒体の体積fjri:粒度分布fjrをもつ粉粒体に含まれる、半径riをもつ任意のi粒子の度数rk:混合粉粒体に含まれる任意のk粒子の半径fjrk:粒度分布fjrをもつ粉粒体に含まれる、半径rkをもつ任意のk粒子の度数請求項1または2に記載の算出方法で算出された粉粒体の充填率pから、下記c式に基づき粉粒体の空隙率δを算出することを特徴とする粉粒体の空隙率の算出方法。粉粒体がセメントであることを特徴とする請求項1~3のいずれかに記載の粉粒体の充填率または空隙率の算出方法。上述したように、セメントの配合原料粉の充填率の予想は、経験則によるデータの蓄積に基づいたマニュアルにより、最終的には勘を頼りになされているが、充填率は厳密には原材料の性質粒度分布、粒径などに依存するため、原材料の性質が変わると特定の原材料にしか成り立たない経験則はもちろん通用しなくなるため、新たな原料に適した経験則が必要となる。このような経験則を、考えられるあらゆる原材料ごとにすべて得るには膨大な数の実験が必要となり、現実的にはまず不可能である。そこで新たな原材料に対応して、それに類似した原材料の既存のデータをもとに経験と勘を頼りに予想をするしかないが、このような方法では科学的な精密さにおいて限度があり、実際に日々異なる原材料で一定の充填率を精密に確保することは至難の業である。したがって本発明の目的は、セメントなどのような粉粒体の充填率や空隙率を精密且つ簡便に算出することができる算出方法を提供することにある。上記課題を解決するための本発明の要旨は以下のとおりである。[1]粉粒体の粒度分布fを求め、この粒度分布fから下記a式により当該粉粒体の充填率pを求めることを特徴とする粉粒体の充填率の算出方法。[2]上記[1]の算出方法において、粒度分布f1r,f2r,...をもつ2種以上の粉粒体をそれぞれ体積V1:V2:...で混合した混合粉粒体について、下記b式で与えられる混合粉粒体の度数Friを用いてa式により当該混合粉粒体の充填率pを求めることを特徴とする粉粒体の充填率の算出方法。ここでPj:粒度分布fjrをもつ粉粒体の充填率Vj:粒度分布fjrをもつ粉粒体の体積fjri:粒度分布fjrをもつ粉粒体に含まれる、半径riをもつ任意のi粒子の度数rk:混合粉粒体に含まれる任意のk粒子の半径fjrk:粒度分布fjrをもつ粉粒体に含まれる、半径rkをもつ任意のk粒子の度数[3]上記[1]または[2]の算出方法で算出された粉粒体の充填率pから、下記c式に基づき粉粒体の空隙率δを算出することを特徴とする粉粒体の空隙率の算出方法。[4]上記[1]~[3]のいずれかの算出方法において、粉粒体がセメントであることを特徴とする粉粒体の充填率または空隙率の算出方法。本発明によれば、粉粒体の充填率または空隙率、特にセメントなどのような2種以上の原料粉を混合した粉粒体の充填率や空隙率を、精密且つ簡便に算出することができる。粒度分布をもつ充填率の基礎理論は統計的手法を導入することにより導出することができ、これによって粒度分布関数に基づいて粒子系の充填率を計算することが可能になる。この理論を用いることで、異なる粒度分布をもつ粒子系を混合した系の充填率を計算することも可能となる。粉粒体の充填状態を解析するために、多くの基礎的な理論や計算が、球充填系という基本分野において行われてきた。一つの基本問題として、同一半径をもつ球のランダム充填密度が計算機シミュレーションによって行われてきた。本発明者は、粒度分布をもつ系の充填率を理論的に求めるために、以下のような算出方法を創案した。一般に、粉体密度の計算は、粒子形状、粒子間相互作用などに起因する困難な問題を内蔵する。従来、多くの試みがこれらの問題に関してなされてきたが、充填密度においてきわめて影響の大きい、粒度分布を厳密に取り扱った汎用的な理論研究は未だなされていない。したがって、最も重要なテーマは粒度分布関数を伴う充填密度の基礎理論であると考えられる。それに付け加えるかたちで、粒形、粒子間相互作用などの影響を、この基礎理論に基づいて考察する必要がある。稠密充填粒子系を取り扱うにあたり、簡単のため、粒子の形状は球と仮定する。球粒子系の分布関数frは、規格化されている。図1に、粉粒体の粒度分布frを模式的に示す。この粉粒体の粒度分布frにおいて、半径riの任意のi粒子の度数がfriである。なお、通常、半径riは或る粒径範囲の代表径が用いられる。ここで、riは粉粒体に含まれる任意のi粒子の半径である。粒子系の空隙率を解析するために、粉粒体を構成する任意のj粒子半径rjの周囲に任意のi粒子半径riを稠密充填した状態すなわち、j粒子、i粒子ともに球であるとし、j粒子の周囲をi粒子が接触して取り囲んでいる状態を考える。図2に、この状態を模式的に示すが、j粒子とi粒子との間の隙間は、j粒子の中心から半径rjriの仮想球の中に形成されるので、この半径rjriの仮想球を考える。そして、j粒子の表面と仮想球面とで囲まれる空隙の合計をΔVjiとおく。このΔVjiは、次の2つの要素からなる。すなわち、一つの要素は、j粒子の中心0を原点とするi粒子に接する立体角Ωにより切り取られるi粒子の下部j粒子面側につくられる空隙部分の合計ΔVjiでありこれはj粒子に接する全てのi粒子それぞれに付随するものである、もう一つの要素は、j粒子上の空隙部分i粒子がない部分が、j粒子の中心0を原点とする立体角Ωにより切り取られる空間のうちj粒子表面と仮想球面とで囲まれる空隙部分の合計ΔVjiである。したがって、j粒子の表面と仮想球面とで囲まれる空隙の合計ΔVjiは、下記2式で表される。そこで、j粒子の周囲にi粒子が接触して配位したときi粒子が最稠密数以下で配位したとき、仮想球内における空隙率δjiを一般に、下記3式のように書くことにする。Vmは仮想球の体積である。粒度分布fをもつ粒子系の場合、半径riの粒子の数niは、下記4式で与えられる。この4式の左辺はi粒子の占有率と呼ばれる。半径rjの一つのj粒子周りに半径riのi粒子が接触して最稠密充填をなしている場合すなわち、半径rjのj粒子周りに同一の半径riをもつi粒子をj粒子に接触させて最大限に稠密に配位させた場合、j粒子周りの半径rjriの仮想球内の空隙率δjiは上記3式で求めることができ、したがって、その場合におけるj粒子周りの仮想球内の充填率pjiは、下記5式で与えられる。粒子系における半径riのi粒子が占める表面積Siの比率は、粒子系の全粒子の総表面積をSとおくとSi/Sで与えられるが、これを用いると、j粒子周りの仮想球内の充填率は、下記6式で与えられる。系を構成する粒子が球の場合、SおよびSjは次式で与えられるなお、r2は式を参照。j粒子が系の中で占める体積Vjの比率は、系の全粒子の総体積をVと置くと、Vj/Vで与えられるから、平均充填率は下記7式で与えられる。系を構成する粒子が球の場合、Vjは下式で与えられるから、充填率pは、上記7式より、下記8式で表されることになる。したがって、系の空隙率δは下記9式で与えられる。上記8式を計算するためには、上記2式のΔVjiを幾何学的計算あるいはコンピューターシミュレーションによる数値計算で求めればよい。上記8式および9式を用いることにより、任意の粒度分布fをもつ球粒子系の充填率と空隙率を計算することが可能である。このため本発明では、粉粒体の粒度分布frを求め、この粒度分布frから上記8式に基づき、当該粉粒体の充填率pを求める。また、上記9式に基づき、当該粉粒体の空隙率δを求める。ここで、粉粒体の粒度分布fを求める方法としては、SEM走査型電子顕微鏡による観察、レーザー光の回折を利用した粒度分布測定法がある。種々の粒度分布をもつ複数の粒子系の混合系については、以下のようにして算出を行う。粒度分布f1r,f2r,...をもつ2種以上の粉粒体をそれぞれ体積V1:V2:...で混合した混合粉粒体の場合、この混合系の度数Friは計算の結果、下記10式で与えられる。ここでPj:粒度分布fjrをもつ粉粒体の充填率Vj:粒度分布fjrをもつ粉粒体の体積fjri:粒度分布fjrをもつ粉粒体に含まれる、半径riをもつ任意のi粒子の度数rk:混合粉粒体に含まれる任意のk粒子の半径fjrk:粒度分布fjrをもつ粉粒体に含まれる、半径rkをもつ任意のk粒子の度数pjは粒度分布fjをもつ系の充填率である。上記8式で度数friを上記10式の度数Friで置き換えることにより、混合系の充填率を計算することができる。従来では、ある1種類の原材料について、原材料の粒度分布を、よく知られている正規分布や対数正規分布、あるいはロジンラムラー分布などの統計分布函数を用いて近似することにより充填率の算出を行ってきた。しかし、このような手法には、次のような二つの根本的な問題が内包されている。その一つは、実際に測定される粒度分布が、そのような統計函数で説明できるような場合はきわめて限られること、もう一つの問題は、そのような統計函数を用いた充填率の計算は、半経験的あるいは現象論的なものであり、実際に充填率を計算するためには経験的な定数を導入しなければならず、一般に計算の基となる理論が汎用的ではないことである。さらに、その次の段階として、粒径も粒度分布も異なるさまざまな種類の原料をある割合で混合したときの充填率を求める場合となると、そのような計算を行う手法や理論は皆無であるため、ある程度の経験則と勘に頼らざるを得なった。本発明は、ある原料についての粒度分布がわかれば、精密かつ簡便に充填率を算出できるものである。すなわち、粒度分布の測定データをインプットするだけで充填率や空隙率を算出することができる。理論計算が任意の粒度分布に適用されるため、測定データではなく粒度分布を数学的な統計函数に置き換えて計算することができるため、統計函数のパラメータと充填率との相関を研究するためのツールとしても活用することができる。実際の配合粉をシミュレーションするためには数種類の原料粉の配合比を変えて充填率を計算しなければならないが、これに対応してそれぞれの原料粉の粒度分布のデータをインプットするだけで、これらを任意の割合で混合したときの充填率を算出することができるため、どのような原料を用いる場合でもそれぞれの原料粉の粒度分布のデータさえあれば充填率を最大にする混合比率を計算によって算出することができる。このように、理論により充填率を計算できるため、原料粉の配合や操作が経験や勘によらない厳密なものとなり、一定の基準を満たすセメントを安定して製造することが可能となる。本発明は、セメント産業やセラミック製造における原料の配合、その他、粉粒体を取り扱う全ての分野に適用可能である。[実施例1]半径1mmと半径5mmの2種類の球粒子を混合した粒子系の充填率について、本発明法により充填率を算出し、実測値と比較した結果を図3に示す。図3の横軸は、粒子系における半径1mmの球粒子が占める体積V1と半径5mmの球粒子が占める体積V2の体積比V1/V2である。本発明例では、上記2式のΔVjiをコンピューターシミュレーションによる数値計算で求め、8式の計算を行った。また、実測値は粒子系を液体で浸し、それに要する液体の体積を測定することで空隙率を算出し、それに基づき求めた。図中、実線で示したのが本発明法による算出値であるが、実測値をほぼ定性的定量的に説明することが明らかである。[実施例2]図4で示すような粒度分布を有する2種類の原料A,Bを混合した混合粉粒体の充填率を、本発明法により10式および8式を用いて算出した。その結果を図5に示す。図5の横軸は、混合粉粒体中での原料Aの質量割合である。本実施例でも、上記2式のΔVjiをコンピューターシミュレーションによる数値計算で求め、10式と8式の計算を行った。図5によれば、原料A,Bを混合した混合粉粒体の充填率が最大値をとるのは、混合比A:B=0.74:0.26の場合であることが判る。"