(57)【要約】【課題】物体からその表面に流れる流体により発生する騒音の位置を高精度で検出することができるようにする。【解決手段】音源探査装置は回転羽根車１０を回転駆動させて気流を発生させた状態のもとで回転羽根車１０から発生する騒音源の位置を検出する。回転羽根車１０に向けて流れる気流には粒子が供給され、回転羽根車１０の回転方向の位置はフォトセンサ１４により検出される。フォトセンサ１４からの信号により回転羽根車１０の被測定部位Ａにレーザ光を照射すると同時に被測定部位Ａの対をなす画像を時間差を持ってＣＣＤカメラ１７により撮影し、被測定部位Ａからの音圧をマイクロフォン１８により測定する。撮影された画像に基づいて粒子画像流速測定法により流体速度ベクトルを演算し、このデータと音圧データとに基づいて音圧と流体速度の相互相関計算により音源分布マップを演算する。【選択図】図１
2008064692.tif 000002
【特許請求の範囲】
【請求項１】
回転羽根車を回転駆動させて気流を発生させた状態のもとで前記回転羽根車から発生する騒音源の位置を検出する音源探査方法であって、
前記回転羽根車の周辺に粒子発生機から可視化粒子を供給して前記回転羽根車の周辺を可視化粒子で充満させ、
前記回転羽根車の位置を検出する位置検出手段からの信号に基づいて前記回転羽根車の被測定部位に光源からの光を照射するとともに光の照射に同期させて前記被測定部位の少なくとも２つで一対をなす画像を時間差を持って撮影手段により複数回撮影し、
前記撮影手段によるそれぞれの撮影時における前記被測定部位からの音圧をマイクロフォンにより測定し、
前記撮影手段により撮影された画像に基づいて粒子画像流速測定法などにより流体速度ベクトルを演算し、
前記マイクロフォンにより測定された音圧と前記流体速度ベクトルとのデータに基づいて音圧と気流速度の相互相関計算により音源分布マップを演算し、
前記音源分布マップを音源分布マップ出力部に出力することを特徴とする音源探査方法。
【請求項２】
被測定物に流体を流した状態のもとで前記被測定物から発生する騒音源の位置を検出する音源探査方法であって、
前記被測定物に向けて流される流体中に可視化粒子を供給して前記被測定物を前記可視化粒子で充満させ、
前記被測定物に光源から照射される被測定部位の少なくとも２つで一対をなす画像を撮影手段により複数回撮影し、
前記撮影手段によるそれぞれの撮影時における前記被測定部位からの音圧をマイクロフォンにより測定し、
前記撮影手段により撮り込まれた画像に基づいて粒子画像流速測定法などにより流体速度ベクトルを演算し、
前記マイクロフォンにより測定された音圧と前記流体速度ベクトルとのデータに基づいて音圧と気流速度の相互相関計算により音源分布マップを演算し、
前記音源分布マップを音源分布マップ出力部に出力することを特徴とする音源探査方法。
【請求項３】
回転羽根車を回転駆動させて気流を発生させた状態のもとで前記回転羽根車から発生する騒音源の位置を検出する音源探査装置であって、
前記回転羽根車の周辺に可視化粒子を供給して前記回転羽根車の周辺を可視化粒子で充満させる粒子発生機と、
前記回転羽根車の位置を検出する位置検出手段と、
前記位置検出手段からの信号に基づいて前記回転羽根車の被測定部位に光を照射する光源と、
光の照射に同期させて前記被測定部位の少なくとも２つで１対をなす画像を時間差を持って複数回撮影する撮影手段と、
前記撮影手段によるそれぞれの撮影時における前記被測定部位からの音圧を測定するマイクロフォンと、
前記撮影手段により撮影された画像に基づいて粒子画像流速測定法などにより流体速度ベクトルを演算する流体速度ベクトル演算部と、
前記マイクロフォンにより測定された音圧と前記流体速度ベクトルとのデータに基づいて音圧と流体速度の相互相関計算により音源分布マップを演算するマップ演算部と、
前記音源分布マップを出力する音源分布マップ出力部とを有することを特徴とする音源探査装置。
【請求項４】
被測定物に流体を流した状態のもとで前記被測定物から発生する騒音源の位置を検出する音源探査装置であって、
前記被測定物に向けて流される流体中に可視化粒子を供給して前記被測定物を前記可視化粒子で充満させる粒子発生機と、
前記被測定物の特定の被測定部位に光を照射する光源と、
光の照射に同期させて前記被測定部位の少なくとも２つで一対をなす画像を時間差を持って複数回撮影する撮影手段と、
前記撮影手段によるそれぞれの撮影時における前記被測定部位からの音圧を測定するマイクロフォンと、
前記撮影手段により撮影された画像に基づいて粒子画像流速測定法などにより流体速度ベクトルを演算する流体速度ベクトル演算部と、
前記マイクロフォンにより測定された音圧データと前記流体速度ベクトルとに基づいて音圧と流体速度の相互相関計算により音源分布マップを演算するマップ演算部と、
前記音源分布マップを音源分布マップ出力部に出力する音源分布マップ出力部とを有することを特徴とする音源探査装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、回転送風機や圧縮機における回転羽根車のように空気等の流体が流れる物体に起因した騒音発生源の位置を検出する音源探査技術に関する。
【背景技術】
【０００２】
機械的エネルギーを空気等の流体のエネルギーに変換する流体圧機器としては、送風機、圧縮機および真空ポンプなどがあり、送風機には吐出圧力と吸込圧力との比である圧力比が1.1未満のファンと、圧力比が1.1以上２未満のブロワとがあり、圧縮機は圧力比を2.2以上に空気やガスを加圧する流体圧機器である。これらの流体圧機器はモータにより回転駆動されてモータからの機械的エネルギーを流体のエネルギーに変換するために、ファンブレードなどの回転羽根車つまり回転翼を有している。
【０００３】
流体圧機器の作動時には回転羽根車から騒音が発生することになるので、静粛性に優れた流体圧機器を設計するためには、回転羽根車のどの位置から騒音が発生しているかを検出することが必要となっている。回転羽根車などの物体の表面または表面から離間した周辺の空気による騒音源の位置を検出するために音響ホログラフィ法や音響インテンシティ法などが用いられているが、これらの音響ホログラフィ法および音響インテンシティ法は、被測定物が高速で移動や回転しているような場合には、騒音源を正確に検出することは困難であった。
【０００４】
特許文献１には、騒音発生源から２つのマイクロフォンに到達する時間差から騒音発生源までの距離を求めることにより音源を探査するようにした技術が記載されている。
【特許文献１】特開２００２−３３３３６７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、この技術では風洞を用いて風速分布がない状態のもとで測定する必要があり、二次元翼の解析に限られて、回転羽根車のように羽根車の径方向内周部と外周部とで風速分布に差のある物体の解析を行うことが困難である。また、従来の回転翼の解析方法には、Ｘ−Ｙの２軸の測定面にマイクロフォンを多数並べることにより二次元座標の音圧を測定する方法があるが、この方法ではマイクロフォンの組立体が大型になったり、また、無指向性マイクロフォンを用いた場合は周辺の音も検知してしまい、音源の探査精度に限度がある。更に、マイクロフォンのみによる音源探査方法では、測定できるのはあくまでも音圧であり、例えば音源位置で起きている流れ現象を捉えることはできない欠点もある。
【０００６】
本発明の目的は、物体からその表面または表面から離間した周辺の空間に流れる流体により発生する騒音の位置を高精度で検出することができるようにすることにある。
【０００７】
本発明の目的は、回転羽根車からその表面または表面から離間した周辺の空間に流れる流体により発生する騒音の位置と流れ構造を高精度で検出することができるようにすることにある。
【課題を解決するための手段】
【０００８】
本発明の音源探査方法は、回転羽根車を回転駆動させて気流を発生させた状態のもとで前記回転羽根車から発生する騒音源の位置を検出する音源探査方法であって、前記回転羽根車の周辺に粒子発生機から可視化粒子を供給して前記回転羽根車の周辺を可視化粒子で充満させ、前記回転羽根車の位置を検出する位置検出手段からの信号に基づいて前記回転羽根車の被測定部位に光源からの光を照射するとともに光の照射に同期させて前記被測定部位の少なくとも２つで一対をなす画像を時間差を持って撮影手段により複数回撮影し、前記撮影手段によるそれぞれの撮影時における前記被測定部位からの音圧をマイクロフォンにより測定し、前記撮影手段により撮影された画像に基づいて粒子画像流速測定法などにより流体速度ベクトルを演算し、前記マイクロフォンにより測定された音圧と前記流体速度ベクトルとのデータに基づいて音圧と気流速度の相互相関計算により音源分布マップを演算し、前記音源分布マップを音源分布マップ出力部に出力することを特徴とする。
【０００９】
本発明の音源探査方法は、被測定物に流体を流した状態のもとで前記被測定物から発生する騒音源の位置を検出する音源探査方法であって、前記被測定物に向けて流される流体中に可視化粒子を供給して前記被測定物を前記可視化粒子で充満させ、前記被測定物に光源から照射される被測定部位の少なくとも２つで一対をなす画像を撮影手段により複数回撮影し、前記撮影手段によるそれぞれの撮影時における前記被測定部位からの音圧をマイクロフォンにより測定し、前記撮影手段により撮り込まれた画像に基づいて粒子画像流速測定法などにより流体速度ベクトルを演算し、前記マイクロフォンにより測定された音圧と前記流体速度ベクトルとのデータに基づいて音圧と気流速度の相互相関計算により音源分布マップを演算し、前記音源分布マップを音源分布マップ出力部に出力することを特徴とする。
【００１０】
本発明の音源探査装置は、回転羽根車を回転駆動させて気流を発生させた状態のもとで前記回転羽根車から発生する騒音源の位置を検出する音源探査装置であって、前記回転羽根車の周辺に可視化粒子を供給して前記回転羽根車の周辺を可視化粒子で充満させる粒子発生機と、前記回転羽根車の位置を検出する位置検出手段と、前記位置検出手段からの信号に基づいて前記回転羽根車の被測定部位に光を照射する光源と、光の照射に同期させて前記被測定部位の少なくとも２つで１対をなす画像を時間差を持って複数回撮影する撮影手段と、前記撮影手段によるそれぞれの撮影時における前記被測定部位からの音圧を測定するマイクロフォンと、前記撮影手段により撮影された画像に基づいて粒子画像流速測定法などにより流体速度ベクトルを演算する流体速度ベクトル演算部と、前記マイクロフォンにより測定された音圧と前記流体速度ベクトルとのデータに基づいて音圧と流体速度の相互相関計算により音源分布マップを演算するマップ演算部と、前記音源分布マップを出力する音源分布マップ出力部とを有することを特徴とする。
【００１１】
本発明の音源探査装置は、被測定物に流体を流した状態のもとで前記被測定物から発生する騒音源の位置を検出する音源探査装置であって、前記被測定物に向けて流される流体中に可視化粒子を供給して前記被測定物を前記可視化粒子で充満させる粒子発生機と、前記被測定物の特定の被測定部位に光を照射する光源と、光の照射に同期させて前記被測定部位の少なくとも２つで一対をなす画像を時間差を持って複数回撮影する撮影手段と、前記撮影手段によるそれぞれの撮影時における前記被測定部位からの音圧を測定するマイクロフォンと、前記撮影手段により撮影された画像に基づいて粒子画像流速測定法などにより流体速度ベクトルを演算する流体速度ベクトル演算部と、前記マイクロフォンにより測定された音圧データと前記流体速度ベクトルとに基づいて音圧と流体速度の相互相関計算により音源分布マップを演算するマップ演算部と、前記音源分布マップを音源分布マップ出力部に出力する音源分布マップ出力部とを有することを特徴とする。
【発明の効果】
【００１２】
被測定部位を撮影した画像に基づいてＰＩＶ法つまり粒子画像流速測定法により流体速度ベクトルを演算し、撮影時の被測定部位からの音圧のデータと流体速度ベクトルのデータとに基づいて相互相関計算により音源分布マップを演算するようにしたので、音源分布マップにより被測定部位における騒音源の位置を高い精度で検出することができるとともに騒音源における流体の流れ現象や騒音発生メカニズムを判定することができる。
【００１３】
被測定物としては、その周囲に流体が流れるものであれば、固定物体とすることもでき、ファンブレードなどの回転羽根車とすることもできる。回転羽根車を被測定物とするときには、所定の回転毎に被測定部位の画像を撮影するとともに撮影時の被測定部位からの音圧を測定することにより、経時的変化による音源分布マップにより騒音源の位置を高い精度で検出することができる。
【００１４】
水などの液体中に配置された被測定物に沿って流れる液体に起因した騒音源の位置を探査する場合には、画像撮影時の音圧を測定するためには、水中マイクロフォンなどを使用することになる。
【発明を実施するための最良の形態】
【００１５】
図１は本発明の一実施の形態である音源探査装置を示す概略斜視図であり、図２は図１の音源探査装置の制御回路を示すブロック図である。
【００１６】
被測定物は送風機に使用されるファンブレードつまり回転羽根車１０であり、この回転羽根車１０は測定室内に配置されたサーボモータ１１の主軸１２に取り付けられて回転駆動されるようになっている。回転羽根車１０はサーボモータ１１により矢印Ｒで示す方向に回転駆動され、矢印Ｂで示す方向に空気流を発生させる。回転羽根車１０の回転方向の位置を検出するために、回転羽根車１０には反射フィルム１３が貼り付けられ、この反射フィルム１３に向けて光を照射する発光素子と、反射フィルム１３からの反射光を受光する受光素子とが設けられたフォトセンサ１４が回転羽根車１０に対向して配置されており、フォトセンサ１４の受光素子が反射フィルム１３からの反射光を受光することにより回転羽根車１０における各部位の回転方向の位置が検出される。
【００１７】
測定室内には可視化粒子発生機としての煙発生機１５から煙が供給されるようになっており、この煙発生機１５により生成された煙の微小粒子Ｐは、平均粒子径が１μｍの煙粒子であって、矢印Ｓで示すように回転羽根車１０に向かう気流に供給され、回転羽根車１０の周辺や近傍でほぼ均一に満遍なく満たされる。
【００１８】
回転羽根車１０の被測定部位Ａに向けてレーザ光を照射するために被測定部位Ａに対向させてレーザ発光装置１６が光源として配置されている。このレーザ発光装置１６としては、イットリウムアルミニウムガーネットを母体とし、Ｎｄイオンを活性触媒とするＹＡＧレーザが用いられており、レーザ発光装置１６は２つのレーザ光発生ヘッド部を有し、相互に例えば５μＳ（秒）の時間差を持ってレーザ光を被測定部位Ａに向けて照射する。レーザ発光装置１６からレーザ光が照射されたときにおける被測定部位Ａの画像を撮り込むために被測定部位Ａに対向して撮影手段としてのＣＣＤカメラ１７が配置されている。
【００１９】
被測定部位Ａからの音圧を検出するために被測定部位Ａに対向してマイクロフォン１８が配置されている。このマイクロフォン１８は風圧の影響を受けないように、回転羽根車１０の負圧側に回転羽根車１０から２５〜５０ｍｍの距離に配置されている。
【００２０】
マイクロフォン１８は音源と思われる被測定部位Ａの位置を定めるためにも使用することができる。被測定部位Ａの位置を定めるには、３軸方向に移動するテーブルにマイクロフォン１８を装着し、回転羽根車１０を例えば２０００ｒｐｍで回転させてマイクロフォン１８からの信号により音圧を計算する。音圧Ｐａ（パスカル）は、以下の式(1)により算出される。
【００２１】
【数１】
2008064692.tif 000003
ただし、Ｎは測定点数であり、Ｘｉはマイク電圧値、ｘバーはマイク電圧平均値であり、Ｃは電圧値から音圧パスカルへの変換係数である。
【００２２】
マイクロフォン１８により多数の部位からの音圧を測定することによって、音圧が最も大きい部位を大まかに被測定領域つまり被測定部位Ａとして設定することができる。
【００２３】
図２に示すように、フォトセンサ１４からの出力信号はパルスジェネレータ１９に送られるようになっており、回転羽根車１０の回転に伴ってフォトセンサ１４の前方を反射フィルム１３が通過すると、フォトセンサ１４からの検出信号がパルスジェネレータ１９に送られる。パルスジェネレータ１９からはレーザ発光装置１６とＣＣＤカメラ１７とにトリガー信号が送られるとともに、マイクロフォン１８からのアナログ信号をデジタル信号に変換するためのＡ／Ｄコンバータ２０にもトリガー信号が送られる。被測定部位Ａがレーザ発光装置１６に対向する位置となったときには、トリガー信号に基づいて、レーザ発光装置１６に設けられた２つのレーザ光発生ヘッド部から５μＳの時間差を持ってレーザ光が被測定部位Ａに向けて照射され、この時間差における２つの画像を一対とした被測定部位Ａの画像がＣＣＤカメラ１７に撮り込まれる。
【００２４】
サーボモータ１１、レーザ発光装置１６、ＣＣＤカメラ１７およびマイクロフォン１８にはそれぞれ図示しない駆動制御部から電力が供給されるようになっており、パルスジェネレータ１９からの信号によりレーザ発光装置１６とＣＣＤカメラ１７とＡ／Ｄコンバータ２０の作動が制御される。サーボモータ１１には回転数を制御するためのパルス信号を発生するエンコーダが設けられており、ここからの信号がコントローラ２１に送られてサーボモータ１１の回転数は精密に制御される。コントローラ２１は、制御信号を演算するマイクロプロセッサＣＰＵと、制御プログラム、演算式、データ等が格納されるメモリであるＲＯＭと、一時的にデータを格納するメモリであるＲＡＭとを有するとともに、信号の入出力ポート等を有している。
【００２５】
コントローラ２１にはＣＣＤカメラ１７による撮影信号とＡ／Ｄコンバータ２０を介してデジタル信号化されたマイクロフォン１８による測定信号とが送られるようになっており、それぞれのデータがコントローラ２１内のデータ記録部に格納されるようになっている。
【００２６】
図３はフォトセンサ１４からの信号に基づくＣＣＤカメラ１７とマイクロフォン１８の作動タイミングを示すタイムチャートである。図３に示すように、フォトセンサ１４が反射フィルム１３を検出してから遅れ時間Ｔ１が経過したらパルスジェネレータ１９は第１回目のパルス信号を出力し、第１回目のパルス信号が出力した後に５μｍが経過したら第２回目のパルス信号を出力する。これらのパルス信号に基づいてレーザ発光装置１６に組み込まれた２つのレーザ発生ヘッド部からは５μｍの時間差でレーザ光が照射され、ＣＣＤカメラ１７にはこの時間差における２つの画像Ａ１，Ａ２を取り込み、それぞれの画像データはコントローラ２１内のＲＡＭ等のメモリに格納される。
【００２７】
フォトセンサ１４がオンすると、Ａ／Ｄコンバータ２０に作動信号が送られて所定のＴ２時間にわたりマイクロフォン１８により測定されたアナログ信号の音圧データをＡ／Ｄコンバータ２０でデジタル信号に変換し、コントローラ２１のメモリに格納する。Ｔ２時間は例えば６２６μＳに設定されており、カメラ撮影時における被測定部位Ａからの音圧が測定される。ただし、この時間Ｔ２は５〜２０００μＳの間の任意の時間に設定することができる。マイクロフォン１８は回転羽根車１０に対して上述したように所定の距離だけ離れて配置されているので、回転羽根車１０からマイクロフォン１８に騒音が到達するまでにはＴ３時間がかかり、ＣＣＤカメラ１７により測定されたときの画像から発生した騒音はＴ３時間後の音圧に対応する。マイクロフォン１８を被測定部位Ａから３０ｍｍ離すと、この遅れ時間Ｔ３は８７μＳとなる。
【００２８】
回転羽根車１０はサーボモータ１１により回転駆動されており、フォトセンサ１４からの信号により被測定部位Ａがレーザ光照射位置となる毎に対をなす画像を取り込むとともにマイクロフォン１８により音圧を測定する。例えば、回転羽根車１０を２０００ｒｐｍの速度で回転させた場合には、フォトセンサ１４からの信号により被測定部位Ａがレーザ光照射位置となる毎に対をなす画像を撮り込んでから、約０．２Ｓ経過直後のフォトセンサ１４からの信号により被測定部位Ａがレーザ光照射位置となる毎に対をなす画像を撮り込む動作を繰り返し、複数回、例えば合計１５００個の画像を撮影してそのデータをメモリに格納する。それぞれの画像撮影時における１５００個の音圧データも同様にしてメモリに格納される。
【００２９】
メモリに格納された対をなす２つの画像Ａ１，Ａ２のデータに基づいて、ＰＩＶ(Particle Image Velocimetry)つまり粒子画像流速測定法により被測定部位Ａにおける流体速度ベクトルがコントローラ２１における流体速度ベクトル演算部により演算される。流体速度ベクトルは対をなす２画像からｕ（軸方向）とｖ（半径方向）の二次元方向のベクトルについて演算される。
【００３０】
この粒子画像流速測定法は、空気や水などの流体についての公知の速度計測法であり、流体と同じ速度で運動する微小粒子つまり可視化粒子を流れの中に混入し、その粒子の速度を画像計測法の手法で測定することにより流体自体の速度と方向の分布つまり流体速度ベクトルを求めることができる。粒子画像流速測定法には、二次元ＰＩＶ法と三次元ＰＩＶ法とがあり、二次元ＰＩＶ法には、時間差撮影された２枚の粒子画像から測定点を中心とする画像を切り出して、粒子の平均移動量を流れの速度とする相互相関ＰＩＶ法等がある。
【００３１】
図１に示す音源探査装置においては、２つの画像Ａ１，Ａ２のデータに基づいてＰＩＶ法により二次元の流体速度ベクトルを演算するようにしているが、ＣＣＤカメラ１７を２台用いることにより円周方向を含めた三次元の流体速度ベクトルを演算することができる。これが三次元ＰＩＶ法である。
【００３２】
上述のように１５００個の画像に基づいてメモリ内に格納された流体速度ベクトルのデータとそれぞれの撮影時の音圧のデータとに基づいてコントローラ２１のマップ演算部において、音圧と流体速度の相互相関計算を行って音源分布マップを演算する。軸方向の相互相関計算（Ｒｕｐ）は以下の式(2)により行われ、半径方向の相互相関計算（Ｒｖｐ）は以下の式(3)により行われる。
【００３３】
【数２】
2008064692.tif 000004
【００３４】
【数３】
2008064692.tif 000005
式(2)(3)において、ｕｉは任意軸方向速度、ｕバーは軸方向の平均速度、ｐｉは任意の音圧、ｐバーは音圧の平均、Ｖｔは翼先端速度、ｖｉは任意半径方向速度、ｖバーは半径方向速度の平均をそれぞれ示し、ｕ′は軸方向速度変動、ｖ′は半径方向の速度変動、ｐ′は音圧変動、Ｖｔは翼先端周速度をそれぞれ示す。ｐrmsは音圧の直流成分を除いた実行値を示し、式(1)により算出される。
【００３５】
このようにして得られた音源分布マップは、被測定部位Ａにおいて流体速度つまり気流速度と音圧との相関性の高低分布を示す画像であり、コントローラ２１から音源分布マップ出力部２２に取り出される。音源分布マップ出力部２２としては、音源分布マップを画面に表示するディスプレイ、音源分布マップを印刷するプリンタ、または音源分布マップのデータを格納するＣＤＲ等の記憶媒体があり、音源分布マップに示された流体流速と音圧との相関性の高低分布に基づいて作業者が音源の位置を判定することができる。
【００３６】
このような音源分布マップを回転翼つまり回転羽根車１０の複数の部位について求めると、回転羽根車１０についての騒音発生状態と発生メカニズムを高精度に判定することができる。
【００３７】
上述した音源探査装置によって音源を探査する手順について説明すると、暗室などのように外部から光が入り込まないようにした測定室において、回転羽根車１０の音源と推測される複数の位置に対向させてマイクロフォン１８を用いて回転羽根車１０から発生する騒音を測定する。そのときには、フォトセンサ１４からの基準位置信号により所定時間、例えば６２６μＳにわたりサンプリングしてコントローラ２１において音圧を計算する。このようにして最も音圧レベルの高い領域を調べてその部位を測定すべき特定の被測定部位Ａと設定する。音圧は上記式(1)により算出される。
【００３８】
次いで、被測定部位Ａに向けてレーザ光が到達するようにレーザ発光装置１６を配置するとともに被測定部位Ａに向けてＣＣＤカメラ１７を配置し、回転羽根車１０を回転させるとともに煙発生機１５から回転羽根車１０により生成される気流に向けて可視化微粒子である煙を供給して測定室内に煙を充満させる。この状態のもとで、レーザ発光装置１６に設けられた２つのレーザ光発光ヘッド部から５μＳの時間差でレーザ光を被測定部位Ａに向けて照射し、その画像をＣＣＤカメラ１７により撮影し、その撮影に同期させて撮影時における被測定部位Ａからの音圧をマイクロフォン１８により測定する。
【００３９】
回転羽根車１０は例えば２０００ｒｐｍの速度で回転し、約０．２Ｓ毎に被測定部位Ａの画像をＣＣＤカメラ１７により撮影するとともに撮影に同期させて撮影時の音圧を測定する。約０．２Ｓ毎にレーザ光を照射するようにしたのは、一度の照射から次の照射までのレーザ発光装置１６の充電時間に０．２Ｓ程度が必要のためである。このようにして、例えば１５００回の画像撮影とそれぞれの画像撮影時の被測定部位Ａからの音圧の測定とを行った。
【００４０】
ＣＣＤカメラ１７により撮り込まれた一対の画像に基づいてコントローラ２１の流体速度ベクトル演算部において粒子画像流速測定法により流体速度ベクトルｕ，ｖを演算するとともに、コントローラ２１の音圧演算部においてマイクロフォン１８により測定された騒音ノイズに基づいて音圧を演算する。これらの演算を全て撮り込まれた画像と測定された音圧とについて行う。次いで、これらの流体速度ベクトルのデータと音圧のデータとに基づいて上述した相互相関計算により音源分布マップを作成する。
【００４１】
図４（Ａ）は粒子画像流速測定法により得られた速度分布の一例を示す速度分布図であり、図４（Ｂ）は流体速度ベクトルのデータと音圧のデータとに基づいて相互相関計算により得られた音源分布マップの一例を示すマップ図である。図４（Ａ）において、濃淡の濃い部分は流体の速度が高いことを示している。この速度分布図には流体の速度をも矢印で表示することができるが、図４（Ａ）においては矢印の標記は省略されている。図４（Ｂ）において、濃度の濃い部分は、流体速度と音圧との相関性が高いことを示しており、相関性の高さから、濃度が濃い領域は音源であることを示す。音源の判定は、音源分布マップをディスプレイに表示したり、プリンタにより印字することにより行うことができる。このように、マイクロフォン１８により測定される音圧と、ＰＩＶ法による微粒子の流れを測定した物体表面に沿う流体の速度ベクトルとに基づいてこれらの相互相関計算により音源の位置を判定することができるので、騒音源の位置を高い精度で検出することができる。しかも、音源の流れ現象や発生メカニズムを判定することができる。
【００４２】
コントローラ２１のマイクロプロセッサは、上述した流体速度ベクトル演算部、音圧演算部およびマップ演算部等の機能部を有しており、それぞれはコントローラ２１に送られたデータに基づいてＲＯＭに格納された演算式により演算される。図２に示したコントローラ２１は、送り込まれた流体速度ベクトルのデータと音圧のデータとを格納する記録部と、それぞれのデータに基づいて演算する機能を有しているが、サーボモータ１１、レーザ発光装置１６、ＣＣＤカメラ１７およびマイクロフォン１８に対する駆動信号をも供給する機能をコントローラ２１に備えるようにしても良い。
【００４３】
本発明は前記実施の形態に限定されるものではなく、その要旨を逸脱しない範囲で種々変更可能である。例えば、上述した実施の形態においては回転羽根車１０からの音源位置を検出するようにしているが、回転することなく、固定された状態のものを被測定物としてその表面から発生する音源の位置を検出するようにしても良い。また、図示する実施の形態においては、被測定物に流れる気流に起因した音源の位置を検出するようにしているが、液体の中に配置される物体を被測定物としてそこから発生する音源の位置を検出するようにしても良い。
【００４４】
また、式(2)(3)の音圧の項については、マイクロフォンによる音圧以外に、音響粒子速度も含めても良い。流体速度の項（u,v）については、通常の速度(ｍ/ｓ)だけでなく、乱流境界層から発生する二重極子音(壁面圧力変動)や四重極子音(気流乱れ)などとする場合も入れることができる。
【００４５】
例えば、二重極子音uIdおよび四重極子音uIqは以下の各式であらわされる。これらの式において、ρは密度、ｒは音源からの距離、ａは音速、Ｌは音源の代表長さ、ｕは流体速度を示す。
【００４６】
さらに、流体速度の項は、流体速度から求めた渦度ω、歪み速度ｅij、レイノルズ応力τ、乱流強度ｑを求めるようにしても良い。
【００４７】
【数４】
2008064692.tif 000006
【図面の簡単な説明】
【００４８】
【図１】本発明の一実施の形態である音源探査装置を示す概略斜視図である。
【図２】図１の音源探査装置の制御回路を示すブロック図である。
【図３】フォトセンサからの信号に基づくＣＣＤカメラとマイクロフォンの作動タイミングを示すタイムチャートである。
【図４】（Ａ）は粒子画像流速測定法により得られた速度分布の一例を示す速度分布図であり、（Ｂ）は流体速度ベクトルのデータと音圧のデータとに基づいて相互相関計算により得られた音源分布マップの一例を示すマップ図である。
【符号の説明】
【００４９】
１０      回転羽根車（被測定物）
１１      サーボモータ
１３      反射フィルム
１４      フォトセンサ
１５      煙発生機（粒子発生機）
１６      レーザ発光装置（光源）
１７      ＣＣＤカメラ（撮影手段）
１８      マイクロフォン
１９      パルスジェネレータ
２０      Ａ／Ｄコンバータ
２１      コントローラ
２２      音源分布マップ出力部
【図１】
2008064692.tif 000007
【図２】
2008064692.tif 000008
【図３】
2008064692.tif 000009
【図４】
2008064692.tif 000010