(57)【要約】【課題】簡便な取付けでエスカレーターのステップまたは手摺の駆動装置を診断する。【解決手段】回転軸１に巻き付けるように取り付けられたエスカレーター診断装置５は、加速度センサ６、信号記録装置７および判定装置８により構成される。加速度センサ６は周方向加速度９、遠心力方向加速度１０、軸方向加速度１１をそれぞれ測定し、それぞれのデータを出力する。判定装置８はこのデータを基に、回転軸１または軸受２またはスプロケット３またはチェーン４の異常の有無を、予め測定した正常時のデータと比較することによって判定する。【選択図】図１
2007331860.tif 000002
【特許請求の範囲】
【請求項１】
エスカレーターのステップまたは手摺の駆動装置の回転軸にかかる周方向または遠心力方向の加速度信号を検出する加速度センサと、
前記加速度信号から重力加速度成分を除去し周方向または遠心力方向の振動信号を算出する振動信号算出手段と、
前記振動信号を正常時と点検時とで比較して差異が生じた場合には、前記駆動装置を異常と判定し、差異がない場合には、前記駆動装置を正常と判定する判定手段とを設けたことを特徴としたエスカレーター診断装置。
【請求項２】
前記周方向または遠心力方向の重力加速度信号を用いて、前記駆動装置の回転軸の回転角を算出し回転角の時系列信号を出力する回転角検出手段と、
前記回転角の時系列信号を用いて、前記周方向または遠心力方向の振動信号を任意の回転角周期で各回転角毎に加算平均し加算平均信号を出力する加算平均手段と、
前記判定手段は前記加算平均信号を正常時と点検時とで比較することにより、差異を生じた場合には前記駆動装置を異常と判定し、差異を生じない場合には正常と判定することを特徴とした請求項１記載のエスカレーター診断装置。
【請求項３】
前記回転角の時間変化から前記駆動装置の回転速度を算出する回転速度算出手段を設け、
前記判定手段は前記回転速度の時系列信号を正常時と点検時とで比較することにより、差異を生じた場合には前記駆動装置を異常と判定し、差異を生じない場合には正常と判定することを特徴とした請求項２記載のエスカレーター診断装置。
【請求項４】
前記エスカレーターのステップまたは手摺の駆動装置の前記回転軸にかかる軸方向の加速度を振動信号として出力する加速度センサを設け、
前記判定手段は軸方向の前記振動信号を正常時と点検時とで比較して差異が生じた場合には前記駆動装置を異常と判定し、差異がない場合には前記駆動装置を正常と判定することを特徴とした請求項１から請求項３のうちのいずれか１項記載のエスカレーター診断装置。
【請求項５】
前記エスカレーターのステップまたは手摺の駆動装置の前記回転軸にかかる軸方向の振動信号として出力する加速度センサと、
前記回転角の時系列信号を用いて、前記軸方向の振動信号を任意の回転角周期で各回転角毎に加算平均し加算平均信号を出力する加算平均手段と、
前記判定手段は前記軸方向の加算平均信号を正常時と点検時とで比較することにより、差異を生じた場合には前記駆動装置を異常と判定し、差異を生じない場合には正常と判定することを特徴とした請求項２または請求項３記載のエスカレーター診断装置。
【請求項６】
前記判定手段は、前記回転速度の時系列信号、周方向の振動信号、遠心力方向の振動信号、軸方向の振動信号、周方向の加算平均信号、遠心力方向の加算平均信号または軸方向の加算平均信号の少なくともいずれか一つにおける統計特徴量を求め、前記統計特徴量を比較し、正常時と点検時の差異が一定値未満ならば正常、一定値以上ならば異常発生と判定するように、予め設定した閾値を用いる構成としたことを特徴とする請求項５記載のエスカレーター診断装置。
【請求項７】
前記判定手段は、前記回転速度の時系列信号、周方向の振動信号、遠心力方向の振動信号、軸方向の振動信号、周方向の加算平均信号、遠心力方向の加算平均信号または軸方向の加算平均信号の少なくともいずれか一つにおける特定周波数成分の大きさを求め、前記特定周波数成分の大きさを比較し、正常時と点検時の差異が一定値未満ならば正常、一定値以上ならば異常発生との判定結果を出力する手段で構成したことを特徴とする請求項５記載のエスカレーター診断装置。
【請求項８】
前記エスカレーターのステップまたは手摺の駆動装置の回転軸または回転軸端のローラ上に、少なくとも前記加速度センサが取り付けられて成ることを特徴とした請求項１から請求項７のうちのいずれか１項記載のエスカレーター診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、エスカレーターのステップまたは手摺の駆動装置を診断するためのエスカレーター診断装置に関するものである。
【背景技術】
【０００２】
従来のこの種の診断装置は、通過検出器および回転検出器を駆動装置内部に固定設置し、通過検出器によりチェーンリンク間の通過時間間隔を測定すると同時に、駆動機に設けられた回転検出器により回転速度を測定し、この回転速度の測定値により速度係数を算出するとともに、この速度係数に基づき上述の通過時間間隔を補正し、上記通過時間間隔と速度係数の変化および判定値に基づき異常判定する（例えば特許文献１）。
【０００３】
【特許文献１】特許第３５７０９１４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
このような従来の診断装置は、診断を行うために通過検出器および回転検出器などをエスカレーター駆動装置内部に固定設置する必要があり、これらの検出器を取り付けるために多くの手間と時間を要するため、既設エスカレーターを短時間で診断することができないという問題点があった。
【０００５】
本発明は、このような問題点を解決するためになされたものであって、診断装置の設置を簡便にすることを可能とし、かつ、高精度に駆動装置の異常・正常を診断することを可能としたエスカレーター診断装置を提供することを目的としている。
【課題を解決するための手段】
【０００６】
本発明に関わるエスカレーター診断装置は、エスカレーターのステップまたは手摺の駆動装置の回転軸に巻き付ける、もしくは回転軸端のローラ上に取り付けるだけで、診断を可能とする。
【０００７】
そのために本発明のエスカレーター診断装置は、エスカレーターのステップまたは手摺の駆動装置の回転軸にかかる周方向または遠心力方向の加速度信号を検出する加速度センサと、上記加速度信号から重力加速度成分を除去し周方向または遠心力方向の振動信号を算出する振動信号算出手段と、上記振動信号を正常時と点検時とで比較して差異が生じた場合には、上記駆動装置を異常と判定し、差異がない場合には、上記駆動装置を正常と判定する判定手段とを設ける構成としたものである。
【発明の効果】
【０００８】
本発明によれば、単にエスカレーターのステップまたは手摺の駆動装置の回転軸に巻き付ける、もしくは回転軸端のローラ側面上に取り付けるだけで、診断作業を可能とする。
また本発明によれば、エスカレーターのステップまたは手摺の駆動装置の回転軸に、単に取り付けた加速度センサの測定データだけで、該駆動装置の各軸方向振動・回転角・回転速度を測定することにより、これらの測定結果を用いて駆動装置の異常もしくは正常を判定するようにしたので、既設エスカレーターに対して簡便に診断装置を設置することが可能となるとともに、精度の高い診断を可能とするものである。
【発明を実施するための最良の形態】
【０００９】
実施の形態１．
図１（ａ）は、本発明の実施の形態１による診断装置を示す図である。エスカレーターの駆動装置は複数の回転軸１から構成されているが、その一例を示したものである。１つ以上の軸受２で支持された回転軸１には、１つ以上のスプロケット３が取り付けられており、１つ以上のチェーン４を介して、他の回転軸とつながっている。本図は、軸受２、スプロケット３およびチェーン４の配置についての一例を示したものに過ぎず、これらのうちの一つ以上の上記構成要素がなくてもよいし、上述の構成要素が回転軸１と接触された配置となっていれば、他の配置であってもよい。
診断装置５は、回転軸１に巻き付けるように取り付けられる。また診断装置５は、図１（ｂ）に示すように加速度センサ６、信号記録装置７および判定装置８により構成される。なお、加速度センサ６は図１（ｂ）に示すように、周方向加速度９、遠心力方向加速度１０および、軸方向加速度１１をそれぞれ測定し、それぞれのデータを出力する。
【００１０】
図２は、診断装置５を示すブロック図である。加速度センサ６で測定された、周方向加速度９、遠心力方向加速度１０および、軸方向加速度１１は、信号記録装置７のデジタル変換部１２において、一定時間間隔毎にＡ／Ｄ変換（アナログデジタル変換）され、メモリ１３に記憶される。メモリ１３に格納された信号をそれぞれ、周方向信号１４、遠心力方向信号１５および、軸方向振動信号１６と呼ぶ。
【００１１】
周方向信号１４、遠心力方向信号１５および、軸方向振動信号１６の一例を図３に示す。図３より明らかなように周方向信号１４または遠心力方向信号１５は、回転軸１の振動成分に加え、回転軸１の１回転を１周期とする正弦波波形が重畳していることがわかる。この重畳された正弦波波形成分は重力加速度成分を表しており、回転軸１の回転角に応じて重力加速度方向に対する加速度センサ６の角度が変化することにより、周方向信号１４または遠心力方向信号１５に重畳する重力加速度成分もまた変化することは明らかである。
【００１２】
判定装置８の回転角検出手段１７は、周方向信号１４および遠心力方向信号１５から回転角の時系列信号１８を検出する。図３を用いて回転角の時系列信号１８の検出方法について以下、説明する。
まず回転角検出手段１７は、周方向信号１４および遠心力方向信号１５に対して、回転軸１の回転周波数よりも大きいカットオフ周波数を有するローパスフィルタを適用し、重力加速度成分に相当する周方向重力加速度信号１９および遠心力方向重力加速度信号２０を検出または算出する。
【００１３】
ここで加速度センサの傾きを、図１（ｂ）の基準位置１０１を０度として反時計周りにプラスの角度であるとすると、サンプリング時刻ｔにおける加速度センサの傾きθ（度）は、下記の式（１）〜（４）により算出される。
ただし、サンプリング時刻ｔにおける周方向重力加速度信号をＸ（ｔ）、遠心力方向重力加速度信号Ｙ（ｔ）とし、図中の９および１０の矢印の方向をプラスとする。
【００１４】
Ｘ（ｔ−１）≧０、Ｙ（ｔ）≧０のとき、
θ（ｔ）＝ｔａｎ-1（｜Ｙ（ｔ）／Ｘ（ｔ）｜）／π×１８０（度）・・・（１）
Ｘ（ｔ−１）＞０、Ｙ（ｔ）≧０のとき、
θ（ｔ）＝１８０−ｔａｎ-1（｜Ｙ（ｔ）／Ｘ（ｔ）｜）／π×１８０（度）・・（２）
Ｘ（ｔ−１）＜０、Ｙ（ｔ）＜０のとき、
θ（ｔ）＝１８０＋ｔａｎ-1（｜Ｙ（ｔ）／Ｘ（ｔ）｜）／π×１８０（度）・・（３）
Ｘ（ｔ−１）≧０、Ｙ（ｔ）＜０のとき、
θ（ｔ）＝３６０−ｔａｎ-1（｜Ｙ（ｔ）／Ｘ（ｔ）｜）／π×１８０（度）・・（４）
【００１５】
したがって、θ（ｔ）および累積回転数Ｃを用いて、サンプリング時刻ｔにおける回転角Ｒ（ｔ）は、下記の式（５）により求められる。
Ｒ（ｔ）＝θ（ｔ）＋Ｃ×３６０（度）  ・・・（５）
ただし、
Ｘ（ｔ−１）≧０かつＹ（ｔ−１）＜０かつ
Ｘ（ｔ）≧０かつＹ（ｔ）≧０のとき、Ｃ＝Ｃ＋１
Ｘ（ｔ−１）≧０かつＹ（ｔ−１）≧０かつ
Ｘ（ｔ）≧０かつＹ（ｔ）＜０のとき、Ｃ＝Ｃ−１
【００１６】
なお上記では、周方向信号１４および遠心力方向信号１５が両方とも測定される場合における回転角の時系列信号１８の検出方法について述べたが、回転角の時系列信号１８の検出は周方向信号１４または遠心力方向信号１５のいずれかのみが測定される構成、すなわち周方向加速度９または遠心力方向加速度１０のいずれかの加速度センサ６のみの構成においても可能である。
【００１７】
以下は、周方向信号１４のみが測定される構成について、図４を用いて回転角の時系列信号１８の検出方法を説明する。
まず回転角検出手段１７は、上述と同様にして、重力加速度成分に相当する周方向重力加速度信号１９を得る。次に、周方向重力加速度信号１９にヒルベルト変換を行った信号であるヒルベルト変換信号１０２を算出する。ヒルベルト変換信号１０２は、周方向重力加速度信号１９の９０度分位相を遅らせた信号となる。よって、サンプリング時刻ｔにおける周方向重力加速度信号をＸ（ｔ）、ヒルベルト変換信号をＹ（ｔ）とすると、上述の式（１）〜（５）と同様の方法で、サンプリング時刻ｔにおける回転角Ｒ（ｔ）が得られる。
【００１８】
またここで回転角検出手段１７は、周方向重力加速度信号１９を反転して微分した後、周方向重力加速度信号１９の振幅と同一となるように振幅を正規化した信号である反転微分信号１０３を算出してもよい。この場合、反転微分信号１０３は、周方向重力加速度信号１９の９０度分位相を遅らせた信号となる。よって、サンプリング時刻ｔにおける周方向重力加速度信号をＸ（ｔ）、反転微分信号をＹ（ｔ）とすると、上述の式（１）〜（５）と同様の方法で、サンプリング時刻ｔにおける回転角Ｒ（ｔ）が得られる。
【００１９】
さらにまた本発明では、周方向重力加速度信号１９のゼロクロス点１０４を用いて、回転角の時系列信号１８を求めるように構成してもよい。すなわち、サンプリング時刻ｔにおける周方向重力加速度信号をＸ（ｔ）とすると、Ｘ（ｐ−１）＜０かつＸ（ｐ）≧０なるｐをゼロクロス点１０４とする。サンプリング時刻ｔにおける回転角Ｒ（ｔ）は、サンプリング時刻ｔを挟む直近のゼロクロス点ｍ＜ｔ＜ｎなるｍとｎを用いて、下記の式（６）により求められる。
【００２０】
Ｒ（ｔ）＝（ｔ−ｍ）／（ｎ−ｍ）×３６０＋Ｃ×３６０−９０（度）  ・・・（６）
ただし、
Ｘ（ｔ−１）＜０かつＸ（ｔ）≧０のとき、Ｃ＝Ｃ＋１
【００２１】
なおここで注意すべきことは、ゼロクロス点１０４については、上記では負から正へ変化する点を用いて説明したが、正から負へ変化する点を用いてもよいし、その両方を用いる構成としても良い。
【００２２】
振動信号算出手段２１は、周方向信号１４および遠心力方向信号１５から、重力加速度成分を除いた回転軸１の振動成分であるところの、周方向振動信号２２および遠心力方向振動信号２３を算出する。ここでは図３に示すように、周方向信号１４から周方向重力加速度信号１９を差し引いた信号を周方向振動信号２２とする。同様に、ここでは遠心力方向信号１５から周方向重力加速度信号１９を差し引いた信号を遠心力方向振動信号２３とする。
【００２３】
加算平均手段２４は、回転角の時系列信号１８すなわち、回転角検出手段１７で得られた各サンプリング時刻ｔにおける回転軸１の回転角Ｒ（ｔ）を用いて、周方向振動信号２２、遠心力方向振動信号２３および軸方向振動信号１６について、任意の回転角周期で信号の加算平均を行う。図５は、任意の回転角周期＝３６０度とした場合における、周方向振動信号２２の処理例を示す。図５に示したように加算平均手段２４は、３６０度毎に周方向振動信号２２を切り出し、各回転角毎に加算平均を行う。この例では、３６０度〜７２０度、７２０度〜１０８０度、１０８０度〜１４４０度、・・・の周方向振動信号２２を切り出して、各回転角毎に加算平均することにより、周方向加算平均信号２５を求めている。なおここでは加算平均手段２４は、周方向振動信号２２をそのまま加算平均したが、周方向振動信号２２の包絡線を取った信号を加算平均してもよいし、周方向振動信号２２の包絡線の絶対値を取った信号を加算平均してもよい。
同様に加算平均手段２４は、遠心力方向振動信号２３から遠心力方向加算平均信号２６を、軸方向振動信号１６から軸方向加算平均信号２７を各々算出する。なおここでの加算平均手段２４は、任意の回転角周期＝３６０度として説明したが、いずれの回転角としてもよいことはいうまでもない。
【００２４】
回転速度算出手段２８は、回転角の時系列信号１８から、回転速度の時系列信号２９を求める。すなわち回転速度の時系列信号２９は、任意の整数ｋ（＞０）を用いて、サンプリング時刻ｔ＋ｋにおける回転角Ｒ（ｔ＋ｋ）と、サンプリング時刻ｔにおける回転角Ｒ（ｔ）を用いて、回転速度Ｖ（ｔ）を以下の式（７）により求める。なお回転速度は、回転角を微分することにより数学的に回転速度が求められることは周知の通りであり、この例示した方法のみによらず、ここでは回転角から回転速度を求める手法であればいずれであってもよいことはいうまでもない。
【００２５】
Ｖ（ｔ）＝（Ｒ（ｔ＋ｋ）−Ｒ（ｔ））／（ｋ×サンプリング間隔）  ・・・（７）
【００２６】
判定手段３０は、点検時における回転速度の時系列信号２９、周方向加算平均信号２５、遠心力方向加算平均信号２６または軸方向加算平均信号２７について、予め求めておいた正常時の上記各信号と比較して、差異があれば異常発生として判定結果３１を出力する。
差異の判定方法としては、各信号における統計特徴量を求め、正常時と点検時の差異が一定値未満ならば正常、一定値以上ならば異常発生と判定し出力するよう構成してもよい。この場合、統計特徴量としては、信号のピーク値、平均振幅、標準偏差、尖度、歪度などの一般的な特徴量を用いることができる。
さらにまた、他の差異の判定方法としては、各信号の特定周波数成分の大きさを求め、正常時と点検時の差異が一定値未満ならば正常、一定値以上ならば異常発生と判定し出力するよう構成してもよい。ここで特定周波数成分の大きさは、各信号にフーリエ変換を行って周波数スペクトルを求め、周波数スペクトルのうちの特定周波数成分の大きさを抜き出すことで求めてもよいし、各信号にバンドパスフィルタを施して特定周波数成分のみを抜き出すことで、その大きさを求めてもよい。
また判定手段３０は、加算平均手段２４を用いて任意の回転角周期で信号の加算平均を行った周方向加算平均信号２５、遠心力方向加算平均信号２６または軸方向加算平均信号２７を用いて判定を行っているが、加算平均手段２４を設けずに、周方向振動信号２２、遠心力方向振動信号２３または、軸方向振動信号１６に対して、上述の判定を行うように構成してもよい。
【００２７】
このように本発明のエスカレーター診断装置によれば、エスカレーターのステップまたは手摺の駆動装置の回転軸１に取り付けた加速度センサ６ひとつで各種のデータを取得し、回転速度の時系列信号２９、周方向振動信号２２、遠心力方向振動信号２３、軸方向振動信号１６、周方向加算平均信号２５、遠心力方向加算平均信号２６または軸方向加算平均信号２７を求めて、該駆動装置の異常を判定するようにしたので、既設エスカレーターへ診断装置を簡便に設置することが可能となるとともに、回転軸１、軸受２、スプロケット３またはチェーン４などの各部位の異常の有無を高精度に判定することが可能となる。
【００２８】
実施の形態２．
図６（ａ）は、本発明の実施の形態２によるエスカレーター診断装置を示す図である。実施の形態１と異なる点は、回転軸１上に診断装置５を巻き付けるスペースがない場合において、回転軸端のローラ２０１上に取り付けるように構成した点である。
すなわちこの実施の形態２では、実施の形態１と同様に、１つ以上の軸受２で支持された回転軸１には、１つ以上のスプロケット３が取り付けられており、１つ以上のチェーン４を介して他の回転軸とつながっている。また回転軸端には、ローラ２０１が取り付けられており、診断装置５はこのローラ２０１の側面に取り付けられる。
図６（ａ）は、軸受２およびスプロケット３およびチェーン４およびローラ２０１の配置についての一例を示したものに過ぎず、これらのうちの一つ以上の上述の構成要素がなくてもよいしまた、上述の構成要素が回転軸１と接触された配置となっていれば、他の配置であってもよいことは実施の形態１の場合と同様である。
加えてこの実施の形態２では、本発明のエスカレーター診断装置５はローラ２０１の側面に取り付けられるが、このローラ２０１はここでは回転軸１と回転中心軸が同一となるように取り付けられている。
またこの診断装置５は、図６（ｂ）に示すように、実施の形態１と同じく、加速度センサ６、信号記録装置７および判定装置８により構成される。加速度センサ６は周方向加速度９、遠心力方向加速度１０、軸方向加速度１１をそれぞれ測定し、それぞれのデータを出力する。すなわち加速度センサ６は、すでに述べた実施の形態１と同様にここでは、上記回転軸１の回転に応じて、上述の各加速度を連続した時系列変化のアナログ信号データとして出力している。
その後のデータの処理、判定方法は実施の形態１の場合と同様であるので以下は省略する。
【００２９】
このように本発明実施の形態２のエスカレーター診断装置５は、エスカレーターのステップまたは手摺の駆動装置の回転軸１に直接、巻付けるように取り付けなくとも、回転軸端のローラ２０１の側面に当該回転軸１と回転中心軸が同一となるように取り付けるよう構成しているので、この種装置の機能、性能を何ら損ねることなく、実施の形態１と同等の高精度の異常の有無判定を行うことができる。その結果、回転軸１上に診断装置５を取り付けるスペースがない場合においても設置が可能となる。
ここで、この実施の形態２ではローラ２０１の側面に本発明診断装置５を取り付けるよう構成した例を示したが、本発明診断装置５の取り付け部位は、エスカレーターのステップまたは手摺の駆動装置の回転軸１から、その回転および軸振動などの変化が比較的忠実に伝達される部位であれば、他の部位であってもよい。
【００３０】
以上のように本発明のエスカレーター診断装置は、加速度センサを用いることにより小型軽量化を計ることによって既設の設備に対して簡便に取り付け可能としかつ、当該センサにより検出されたアナログ信号の「加速度値」を連続した時系列変化のデジタル信号（データ）に変換することによってコンピュータなどによる演算処理を可能とし、この種の被診断装置の回転軸１、軸受２、スプロケット３、チェーン４またはローラ２０１などの機械系部位の異常の有無の診断作業を、所定のアルゴリズムによるデジタル信号（データ）処理によって可能とする。その結果、診断作業はコンピュータなどの演算処理による高精度、高能率、高速で上述のような機械系部位の異常の有無判定結果を得ることが可能となる。
このことは、エスカレーターなどの装置が万一の故障の際においては、その復旧作業時間を著しく短縮し、かつ、一刻も早い復旧が望まれるこの種の迅速な作業において、多大な貢献を果たす効果を奏することを意味する。
【図面の簡単な説明】
【００３１】
【図１】（ａ）は本発明の実施の形態１によるエスカレーター診断装置を示す構成図である。（ｂ）は（ａ）における診断装置の拡大断面図である。
【図２】本発明の実施の形態１によるエスカレーター診断装置を示すブロック図である。
【図３】本発明の実施の形態１による回転角検出手段および振動信号検出手段を説明する図である。
【図４】本発明の実施の形態１による回転角検出手段を説明する図である。
【図５】本発明の実施の形態１による加算平均手段を説明する図である。
【図６】（ａ）は本発明の実施の形態２によるエスカレーター診断装置を示す構成図である。（ｂ）は（ａ）における診断装置の拡大断面図である。
【符号の説明】
【００３２】
１  回転軸、２  軸受、３  スプロケット、４  チェーン、５  診断装置、６  加速度センサ、７  信号記録装置、８  判定装置、９  周方向加速度、１０  遠心力方向加速度、１１  軸方向加速度、１２  デジタル変換部、１３  メモリ、１４  周方向信号、１５  遠心力方向信号、１６  軸方向振動信号、１７  回転角検出手段、１８  回転角の時系列信号、１９  周方向重力加速度信号、２０  遠心力方向重力加速度信号、２１  振動信号算出手段、２２  周方向振動信号、２３  遠心力方向振動信号、２４  加算平均手段、２５  周方向加算平均信号、２６  遠心力方向加算平均信号、２７  軸方向加算平均信号、２８  回転速度算出手段、２９  回転速度の時系列信号、３０  判定手段、３１  判定結果、１０１  基準位置、１０２  ヒルベルト変換信号、１０３  反転微分信号、１０４  ゼロクロス点、２０１  ローラ。
【図１】
2007331860.tif 000003
【図２】
2007331860.tif 000004
【図３】
2007331860.tif 000005
【図４】
2007331860.tif 000006
【図５】
2007331860.tif 000007
【図６】
2007331860.tif 000008