JPH032251B2 - - Google Patents

Info

Publication number
JPH032251B2
JPH032251B2 JP57171646A JP17164682A JPH032251B2 JP H032251 B2 JPH032251 B2 JP H032251B2 JP 57171646 A JP57171646 A JP 57171646A JP 17164682 A JP17164682 A JP 17164682A JP H032251 B2 JPH032251 B2 JP H032251B2
Authority
JP
Japan
Prior art keywords
diffraction image
particle size
size distribution
particles
function
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
JP57171646A
Other languages
English (en)
Other versions
JPS5960341A (ja
Inventor
Suezo Nakatate
Hiroyoshi Saito
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
RIKEN
Original Assignee
RIKEN
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by RIKEN filed Critical RIKEN
Priority to JP57171646A priority Critical patent/JPS5960341A/ja
Publication of JPS5960341A publication Critical patent/JPS5960341A/ja
Publication of JPH032251B2 publication Critical patent/JPH032251B2/ja
Granted legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/02Investigating particle size or size distribution
    • G01N15/0205Investigating particle size or size distribution by optical means
    • G01N15/0211Investigating a scatter or diffraction pattern

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Description

【発明の詳細な説明】 本発明は、レーザー回折像を用いる粒度分布の
測定法に関する。粒子の大きさやその分布状態を
計測することは、粉粒体工学、細胞学、大気汚染
計測のみならず、鉄鉱、食品、薬品などの製造業
においても、品質管理、製造工程制御などにとつ
て極めて重要な問題となつている。
従来、レーザー光を用いた粒度分布の測定方法
としては、光散乱分光法や、レーザー光の回折像
から行列計算により求める方法などがあつたが、
広範囲にわたつて、連続な粒度分布を導出するの
は原理的に不可能であつた。
本発明は上記に鑑みなされたものであつて、レ
ーザー回折像を用いて連続的な粒度分布を測定す
る方法を提供することを目的とする。
以下、添付図により本発明を説明する。第1図
は本発明の測定法の実施に用いる装置の一例であ
つて、被測定粒子の回折像を得るための光学系な
らびに信号処理系を示す。
1はコヒーレントな光を発する光源たとえばレ
ーザー装置であり、2は光を平行光とする光学
系、3は粒子がランダムに分布する被測定面であ
り、4はフーリエ変換レンズであり、5は粒子の
回折像(フーリエ変換像)が得られる面である。
6は得られた回折像の光強度を電気信号に変換す
る光電変換装置たとえばダイオードアレイ、TV
カメラ等であり、7は電気信号を処理する装置た
とえばコンピユータ等である。
本発明の目的は、次の工程によつて達成され
る。
ランダムに分布している粒子にレーザー光を
投射して粒子の回折像をつくる。
粒子は球状粒子を考えているので、回転対称
の回折像ができるので、回折像面上の動径ρは
任意の位置に選ぶことができる。
半径aの球状粒子1個の作る回折像強度I0
(ρ)はフーリエ変換を用いて、次のように計
算できる。
I0(ρ)=C・a4・{J1(aρ)/aρ}2 ……(1) ここで、J1(x)は1次の第1種ベツセル関
数。cは定数である。
従つて、粒度分布N(a)をもつ粒子群によ
る回折像強度I(ρ)は次のように書ける。
I(ρ)=K∫ 0I0(ρ)N(a)da =K∫ 0a4{J1(aρ)/aρ}2N(a)da ……(2) aは粒子半径、ρは回折像面上の空間座標、
N(a)は粒度分布、kは定数である。
ここで、動径軸ρ上に光検素子を配置して、
光強度を検出すると、I(ρ)に比例した電気
信号を得ることができる(第2図)。
以下は、I(ρ)を知つて、N(a)を求める
数学的な問題となるので、個々に表れる関数の
物理的な意味はなく、単に計算機内で実現され
るデータである。
新たな関数K(ρ,t)=ρ2J1(Tρ)を考える
(第3図)。このK(ρ,t)を得られた信号I
(ρ)に乗じて、ρについて〔0,∞)で積分
(データの総和)を求めると、新たに関数F
(t)が求まる(第4図)。
F(t)=∫ 0I(ρ)K(ρ,t)dρ=
K∫ 0da∫ 0dρ・a4{J1(aρ)/aρ}2・ρ2J1(t
ρ)Na =K∫ 0a2N(a)da∫ 0
J2 1(aρ)J1(tρ)dρ……(3) 数学公式集より となるので、(3)式は F(t) =K′∫ 02−(2)2N(a)da
……(4) となる。K′は定数である。
従つて、計算機内で、ρ2J1(tρ)なる関数を
作り、データ、I(ρ)との積の総和を計算し
て、新たな関数F(t)を求める。
次に、a2=A、(t/2)2=Tとなるような
座標変換を行うと(第5図)、(4)式は となる。ただし、A<Tの時は0となる。定数
K″は無視してある。
この(5)式を見ると、関数√U(A)〔U(A)
は第6図に図示した単位ステツプ関数であり、
√U(A)は、Aの値を負の領域まで拡張し
た関数(第7図)である。〕と関数N(√)/
√との相互相関関数になつている。従つて、
(5)式の両辺のフーリエ変換を求めてみる。F
(T)のフーリエ変換を(ω)とし、√U
(A)のフーリエ変換をH(ω)、N(√)/√
Aのそれを(ω)とすると、(5)式から、 (ω)=H(ω)・(ω) ……(6) となる。定数=1としているこの式から、 (ω)=(ω)・1/H(ω) となり、逆フイルター1/H(ω)を(ω)
に乗じて、逆フーリエ変換することになる。√
AU(A)のフーリエ変換H(ω)は、 と計算でき、 1/H(ω)=H(ω)・4/π・2(iω)3
……(8) と書ける。フーリエ面である関数に(iω)3を乗
ずることは、実関数面上では3回微分すること
に対応する(第8図)。従つて、(6)式を書き直
すと、 (ω)=(ω)・{(H(ω)(iω)3}……(9
) となり、この(9)式の両辺を逆フーリエ変換する
と、 となり、F(T)とh(T)の相互相関関数を計
算することにより(第9図)、N(√)/√
が求まる。
h(T)とは、H(ω)の逆フーリエ変換√
U(T)を3回微分したものでである。従つて、 h(T)=∂3/∂T3{√U(T)} である。h(T)を計算機上で求めると第10
図のようになる。
最後に座標AをA=a2となるように変換し、
√を乗ずると、N(a)が求まる(第11
図)。
いま、自然状態にあ半径aの粒子がもつ粒度分
布N(a)を、粉粒体工学で一般に用いられてい
る次式の対数正規分布(logarithmic normal
distribution)で仮定すると、 N(a)=1/a・exp(−0.5〔2log a〕2
……(11) 得られる粒子の回折像の光強度を表わす電気信
号I(ρ)は、 I(ρ)=k∫ 0a4{J1(aρ)/aρ}2N(a)da…
…(12) (ただし、ρは回折像面上の空間座標、kは定
数)で表わされ、これは第12図に示す関係曲線
となる。
この電気信号I(ρ)にρ2J1(2√ρ)を乗じ
〔0,∞〕の範囲で積分する、新たな関数F(T)
は次式 F(T)=∫ 0√−N√/√dA
……(13) (ただし、A=a2) で表わされ、これは第13図に示す関係曲線とな
る。
第14図は、√U(T)(ここで、U(T)は
ステツプ関数)を3回微分して求めたオペレータ
h(T)を示し、 h(T)=∂3{√U(T)〕/∂T3 ……(14) で表わされる。
そして、F(T)とh(T)の相互相関関数、す
なわち第13図と第14図との相互相関関数を求
めて、粒度分布N(a)を求めると、第15図の
粒度分布曲線が得られる。図中、イは初めに仮定
した(11)式の分布曲線であり、ロが上記処理からな
る本発明により求めた分布曲線である。一部分を
除いて両分布曲線が一致しており、本発明によれ
ば連続的な粒度分布を極めて精度良く測定できる
ことが理解される。
以上は、粒度分布を対数正規分布と仮定した時
の計算結果を表現したものである。
なお、(11)式で表わされる粒度分布を例にとつて
説明したが、第1図の装置を用いる場合について
説明する。先ず、ランダムに分布している粒子に
レーザー装置1及び光学系2によりレーザー光を
投射し、フーリエ変換レンズ4を介して面5上に
得られる粒子の回折像を、光電変換装置6で光強
度を表わす電気信号I(ρ)に変換する。次いで
処理装置7を用いて、電気信号I(ρ)を上記(12),
(13),(14)式に従つて処理し、(13),(14)式で
表わされるF(T)とh(T)の相互相関関数をと
つて、粒度分布N(a)を求めればよい。
以上説明したように、本発明はレーザー光の回
折像を用いるために、粒子の直径が光の半波長程
度から、第1図の光学系2,4を変えることによ
り、数cm程度までの粒子の粒度分布を測定でき
る。
また、サンプルに特別な加工を必要としないた
め、製産工程への導入が容易にでき、オンライン
の計測ができ計算機処理の高速化により実時間の
測定が可能となる。したがつて本発明は、粒度分
布が連続的に精度良く求まることと相俟つて、製
造品の異常検出や、原料の最適パラメータの決定
など、品質管理、工程管理の迅速化、高精度化に
寄与することができる。
【図面の簡単な説明】
第1図は本発明に用いる装置の光学系と信号処
理系の一例を示す。第2図〜第11図は、数式の
展開過程における関数を示すグラフ、第12図は
本発明を説明するために用いた仮定した粒度分布
から成る粒子の回折像強度を表わすグラフ、第1
3図は第12図に対してρ2J1(2√ρ)を乗じ
積分した後の関数を示すグラフ。第14図は本発
明に用いるオペレータh(T)を示すグラフ。第
15図は第13図と第14図との相互相関関数を
求めて得られた粒度分布N(a)を示すグラフ。 図中の符号:1……レーザー装置、2……光学
系、3……被測定面、4……フーリエ変換レン
ズ、5……回折像が得られる面、6……光電変換
装置、7……電気信号変換装置。

Claims (1)

  1. 【特許請求の範囲】 1 分布している粒子にレーザー光を投射して粒
    子の回折像をつくり、この回折像の光強度の空間
    分布を表わす電気信号 I(ρ)=k∫ 0a4{J1(aρ)/aρ}2N(a)da (ここで、aは粒子半径、ρは回折像が表われ
    る空間の一次変数、N(a)は粒度分布、kは定
    数)をつくり、この電気信号を次式に従つて処理
    して粒度分布N(a)を決定することを特徴とす
    るレーザー回折像を用いる粒度分布測定法。 I(ρ)にρ2J1(2√ρ)を乗じて〔0,∞〕
    の範囲で積分し、a2=Aなる変数変換を行つて
    【式】を導出し、 オペレータ【式】(ここで U(T)はステツプ関数)とF(T)との相関関数
    を導出してN(a)を求める。
JP57171646A 1982-09-30 1982-09-30 レ−ザ−回折像を用いる粒度分布測定法 Granted JPS5960341A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP57171646A JPS5960341A (ja) 1982-09-30 1982-09-30 レ−ザ−回折像を用いる粒度分布測定法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP57171646A JPS5960341A (ja) 1982-09-30 1982-09-30 レ−ザ−回折像を用いる粒度分布測定法

Publications (2)

Publication Number Publication Date
JPS5960341A JPS5960341A (ja) 1984-04-06
JPH032251B2 true JPH032251B2 (ja) 1991-01-14

Family

ID=15927066

Family Applications (1)

Application Number Title Priority Date Filing Date
JP57171646A Granted JPS5960341A (ja) 1982-09-30 1982-09-30 レ−ザ−回折像を用いる粒度分布測定法

Country Status (1)

Country Link
JP (1) JPS5960341A (ja)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0650282B2 (ja) * 1986-01-30 1994-06-29 日機装株式会社 粉粒体粒度オンライン自動分析装置
JPH03115949A (ja) * 1989-09-29 1991-05-16 Shimadzu Corp 粒度分布測定方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4095283A (en) * 1976-07-02 1978-06-13 International Business Machines Corporation First in-first out memory array containing special bits for replacement addressing
ES484355A1 (es) * 1978-09-23 1980-05-16 Stamicarbon Un dispositivo para el estirado lateral continuo de una ban-da de material reticulado a anchura constante

Also Published As

Publication number Publication date
JPS5960341A (ja) 1984-04-06

Similar Documents

Publication Publication Date Title
EP0359681B1 (en) Characterization of particles by modulated dynamic light scattering
Gabas et al. Response of laser diffraction particle sizer to anisometric particles
JP2001337028A (ja) 粒度分布測定方法および装置
Churchill et al. Light-scattering by very dense monodispersions of latex particles
CN104697906A (zh) 一种基于近场散射的颗粒粒度测量装置及方法
WO2024091159A1 (en) A method and an assembly for measuring particle size of particles suspended in a fluid
US4828388A (en) Method of measuring concentration of substances
JPH032251B2 (ja)
Khanam et al. Measurement of two-dimensional crystal shape using digital holography
JP3266107B2 (ja) 粒子個数計測方法および粒子計測装置
JP3702591B2 (ja) 粒度分布測定データの表示装置およびレーザ回折/散乱法に基づく粒度分布測定装置の表示装置
Gucker et al. The angular variation of light scattered by single dioctyl phthalate aerosol droplets
JP3235554B2 (ja) レーザ回折・散乱式粒度分布測定装置
US10643325B2 (en) System and method for image processing and feature recognition
JP3302991B2 (ja) 粒径分布測定装置
EP0590891A1 (en) Method and apparatus for estimating a mixing proportion of different powdery contents
DE112016006272T5 (de) Entferntes erfassen einer Richtung einer lateralen Bewegung eines Objekts unter Verwendung einer auf Phasendifferenz basierenden optischen Bahndrehimpuls-Spektroskopie
JP2000002644A (ja) レーザ回折・散乱式粒度分布測定装置
Marshall et al. Sizing opaque spherical particles using classical matched filters and holographic ring detectors
Barrientos et al. Dynamic three-dimensional displacement analysis of small-scale granular flows by fringe projection and digital image correlation
JPH0327065B2 (ja)
JPS5896238A (ja) 粒子群の平均粒径及び体積濃度測定方法
JP2821200B2 (ja) 粒子測定方法および装置
JP3645759B2 (ja) 粒子径分布解析方法
JP3282580B2 (ja) レーザ回折・散乱式粒度分布測定装置