(57)【要約】【課題】専用の測定器を用いないで吸音材料の物性値を測定できるようにする。【解決手段】イノベーション計算部２３は、所定の物性値（流れ抵抗）をパラメータとして解析により得られた吸音率と位相の周波数特性と測定により得られた吸音率と位相の周波数特性とを用いてイノベーションを計算し、同定部２４は、得られた複数のイノベーションを比較して最も白色信号に近いイノベーションに応じた周波数特性を与える物性値（流れ抵抗）を吸音材料の物性値として同定する。【選択図】    図１
2004347509.tif 000002
【特許請求の範囲】
【請求項１】
所定の物性値をパラメータとし、各物性値における音響インピーダンスを解析し、該音響インピーダンスより吸音率と位相の周波数特性を算出する解析部、
実際の吸音材料の音響インピーダンスを測定し、該音響インピーダンスより吸音率と位相の周波数特性を算出する測定部、
前記測定により得られた吸音率と位相の周波数特性と前記解析により得られた吸音率と位相の周波数特性の相関を演算する相関演算部、
相関が最大となる周波数特性を与える前記物性値を吸音材料の物性値として同定する同定部、
を備えた吸音材料の物性値同定装置。
【請求項２】
所定の物性値をパラメータとし、各物性値における音響インピーダンスを解析し、該音響インピーダンスより吸音率と位相の周波数特性を算出する解析部、
実際の吸音材料の音響インピーダンスを測定し、該音響インピーダンスより吸音率と位相の周波数特性を算出する測定部、
前記解析により得られた吸音率と位相の周波数特性におけるｎ個のポイントの値をＨｎ−１，…，Ｈ２，Ｈ１，Ｈ０、前記測定により得られた吸音率と位相の周波数特性におけるｎ個のポイントの値をＪｎ−１，…，Ｊ２，Ｊ１，Ｊ０とするとき、カルマンフィルタの基本システムが形成されるようにこれらＨｋ，Ｊｋ（ｋ＝０，１，２，…ｎ−１）間の関係を線形モデル化してイノベーションを計算し、同様に、前記解析により得られた別の吸音率と位相の周波数特性におけるｎ個のポイントについてイノベーションを計算するイノベーション計算部、
得られた複数のイノベーションを比較して最も白色信号に近いイノベーションに応じた周波数特性を与える前記物性値を吸音材料の物性値として同定する同定部、
を備えた吸音材料の物性値同定装置。
【請求項３】
吸音材料の物性値を同定する物性値同定装置において、
所定の物性値をパラメータとして解析により得られた吸音率と位相の周波数特性におけるｎ個のポイントの値をＨｎ−１，…，Ｈ２，Ｈ１，Ｈ０、測定により得られた吸音率と位相の周波数特性におけるｎ個のポイントの値をＪｎ−１，…，Ｊ２，Ｊ１，Ｊ０とするとき、カルマンフィルタの基本システムが形成されるようにこれらＨｋ，Ｊｋ（ｋ＝０，１，２，…ｎ−１）間の関係を線形モデル化してイノベーションを計算し、同様に、解析により得られた別の吸音率と位相の周波数特性におけるｎ個のポイントについてイノベーションを計算するイノベーション計算部、
得られた複数のイノベーションを比較して最も白色信号に近いイノベーションに応じた周波数特性を与える前記物性値を吸音材料の物性値として同定する同定部、
を備えた吸音材料の物性値同定装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、吸音材料の物性値同定装置に係わり、特に、相関演算結果あるいはカルマンフィルタのイノベーションを用いた吸音材料の物性値同定装置に関する。
【０００２】
【従来の技術】
多孔質吸音材料は、室内の防音や吸音（音の反射や外への漏出を軽減する）に用いられる壁面材料で、代表的な材質として発泡ウレタンやグラスウールが挙げられる。また、シートやドアトリム等の自動車内装これらに類似するものが有る。この多孔質吸音材料は、連結されて外部に通じている無数の孔を内部に持ち、この穴が持つ音伝搬への抵抗作用によって音を吸収する。
物性値として、発砲ウレタンやダラスウールに代表される多孔質吸音材料の性質を定量的に評価する指標として、以下三種類の物性値が規定されている。
【０００３】
▲１▼流れ抵抗
流れ抵抗は、材料表面垂直方向に一定空気流を通したときの速度と材料表面間の圧力差の比であり、次式
Ｒ＝Δｐ／ｈｖ
により定義される。ただし、Ｒ：流れ抵抗［ｒａｙｌ／ｍ］、Δｐ：材料表面間の圧力差［Ｐａ］、ｈ：材料厚さ［ｍ］、ｖ：流体流速［ｍ／ｖ］、１  ｒａｙｌ＝１  Ｎｓ／ｍ３＝１ｋｇ／ｍ２ｓである。
▲２▼有孔率
有孔率は多孔質材料全体の体積と空隙部（連続気泡）の体積の比であり、次式
Ω＝Ｖｖ／Ｖｔ
により定義される。ただし、Ω：有孔率、Ｖｖ：空隙部の体積［ｍ３］、Ｖｔ：全体の体積［ｍ３］である。
▲３▼構造係数
構造係数は、多孔質中の流体（空気等）について、流体の通常密度と多孔質中での密度の比であり、次式
Ｋｓ＝ρｅ／ρ
により定義される。ただし、Ｋｓ：構造係数、ρｅ：媒質の多孔質中密度［ｋｇ／ｍ３］、ρ：媒質の通常密度［ｋｇ／ｍ３］である。
【０００４】
以上の物性値を同定することができれば、車室壁面に発泡ウレタンを貼付した時の車室内音場（音の強弱の分布状態）がどのようになるかを実機を製作せずに評価可能になる。又、車室内音場を正確に評価可能となれば、たとえば車室内のどの箇所にスピーカを設置すると最適な聴取環境が得られるかを判断可能となる。
従って、多孔質吸音材料の物性値を正しく同定することは車室内の音響システムを構築する上で重要である。
【０００５】
従来は、専用の計測器を用いて物性値を測定するものであった（従来技術１）。すなわち、流れ抵抗、有孔率、構造係数それぞれに専用の測定器を用意してこれら物性値を測定するものであった。たとえば、流れ抵抗測定器では、厚さが既知の多孔質材料に流体（空気）を通じて、前後の差圧と流速を測定し、定義式から算出する。このとき圧力の印可状態でＤＣ法とＡＣ法に分類される。有孔率測定器では、見かけ体積が既知の多孔材料を水中に沈め、体積増加分を測定して材料の体積を求め（この時の体積は空隙を除いた材料自身の体積である）、この二つの体積から有孔率を求める。構造係数測定器では、捻れ率の平方根が構造係数であるため、捻れ率を計測する。
又、従来技術２として入力パラメータに基づいて室型モデルを構築すると共に、所定のスピーカを用いた場合における音場解析を行って表示するものがある（特許文献１）。
【特許文献１】特開平２００２−１２３２６２号
【０００６】
【発明が解決しようとする課題】
しかし、従来技術１では、専用の測定器が必要である。このため、物性値測定のためにはその測定器を所有するか、保有機関に測定を依頼しなければならない。このため測定データ取得には多大な費用と期間を要する問題がある。又、従来技術２では吸音材料の物性値を同定することはできない。
以上から本発明の目的は、専用の測定器を用いないで吸音材料の物性値を測定できるようにすることである。
【０００７】
【課題を解決するための手段】
上記課題は、本発明の第１実施態様によれば、所定の物性値をパラメータとし、各物性値における音響インピーダンスを解析し、該音響インピーダンスより吸音率と位相の周波数特性を算出する解析部、実際の吸音材料の音響インピーダンスを測定し、該音響インピーダンスより吸音率と位相の周波数特性を算出する測定部、前記測定により得られた吸音率と位相の周波数特性と前記解析により得られた吸音率と位相の周波数特性の相関を演算する相関演算部、相関が最大となる周波数特性を与える前記物性値を吸音材料の物性値として同定する同定部、を備えた吸音材料の物性値同定装置により達成される。
【０００８】
また、上記課題は、本発明の第２実施態様によれば、所定の物性値をパラメータとして解析により得られた吸音率と位相の周波数特性におけるｎ個のポイントの値をＨｎ−１，…，Ｈ２，Ｈ１，Ｈ０、測定により得られた吸音率と位相の周波数特性におけるｎ個のポイントの値をＪｎ−１，…，Ｊ２，Ｊ１，Ｊ０とするとき、カルマンフィルタの基本システムが形成されるようにこれらＨｋ，Ｊｋ（ｋ＝０，１，２，…ｎ−１）間の関係を線形モデル化してイノベーションを計算し、同様に、解析により得られた別の吸音率と位相の周波数特性におけるｎ個のポイントについてイノベーションを計算するイノベーション計算部、得られた複数のイノベーションを比較して最も白色信号に近いイノベーションに応じた周波数特性を与える前記物性値を吸音材料の物性値として同定する同定部、を備えた吸音材料の物性値同定装置により達成される。
【０００９】
【発明の実施の形態】
（Ａ）概略説明
図１は本発明の概略説明図であり、図１（Ａ）は第１実施例、図１（Ｂ）は第２実施例の概略説明図である。
図１（Ａ）において、音響解析部１１は、所定の物性値（たとえば流れ抵抗値）をパラメータとし、各流れ抵抗における音響インピーダンスを解析し、該音響インピーダンスより吸音率と位相の周波数特性を算出する。音響インピーダンス測定部１２は、実際の多孔質吸音材料の音響インピーダンスを測定し、該音響インピーダンスより実際の吸音率と位相の周波数特性を算出する。相関演算部１３は、前記測定により得られた吸音率と位相の周波数特性と前記解析により得られた吸音率と位相の周波数特性の相関を演算し、物性値同定部１４は相関が最大となる周波数特性を与える流れ抵抗を多孔質吸音材料の流れ抵抗として同定する。
【００１０】
図１（Ｂ）において、音響解析部２１は、所定の物性値（たとえば流れ抵抗値）をパラメータとし、各流れ抵抗における音響インピーダンスを解析し、該音響インピーダンスより吸音率と位相の周波数特性を算出する。音響インピーダンス測定部２２は、実際の多孔質吸音材料の音響インピーダンスを測定し、該音響インピーダンスより実際の吸音率と位相の周波数特性を算出する。イノベーション演算部２３は、前記解析により求めた吸音率と位相の周波数特性におけるｎ個のポイントの値をＨｎ−１，…，Ｈ２，Ｈ１，Ｈ０、前記測定により得られた吸音率と位相の周波数特性におけるｎ個のポイントの値をＪｎ−１，…，Ｊ２，Ｊ１，Ｊ０とするとき、カルマンフィルタの基本システムが形成されるようにこれらＨｋ，Ｊｋ（ｋ＝０，１，２，…ｎ−１）間の関係を線形モデル化してイノベーションを計算し、同様に、解析により得られた別の吸音率と位相の周波数特性におけるｎ個のポイントについてイノベーションを計算する。物性値同定部２４は、得られた複数のイノベーションを比較して最も白色信号に近いイノベーションに応じた周波数特性を与える流れ抵抗を多孔質吸音材料の流れ抵抗として同定する。
【００１１】
（Ｂ）第１実施例
図２は第１実施例における多孔質吸音材料の物性値同定装置の構成図であり、多孔質吸音材料における３つの物性値のうち流れ抵抗をパラメータとして１０３〜１０６の範囲で可変し、残りの２つの物性値である有孔率を０．９５〜０．９６の範囲の値に固定し、又、構造係数を有孔率から算出される値に固定した場合である。なお、パラメータ範囲は文献等にて絞り込み、その範囲内で２０￣３０パターンについて解析を行うものとする。
【００１２】
解析パラメータ設定部５１は、所定パターンの物性値を音響インピーダンス解析部５２に設定する。音響インピーダンス解析部５２は市販の解析ソフト（たとえばＳＹＳＮＯＩＳＥ）を用いて周波数毎に音響インピーダンス（複素インピーダンス）Ｚ＝Ｒ＋ｊＸを求め、吸音率・位相算出部５３に入力する。吸音率・位相算出部５３は、次の（１）、（２）式
【数１】
2004347509.tif 000003
により、周波数毎に吸音率α０、位相θを演算する。これにより、所定パターンの吸音率、位相それぞれの周波数特性を求めることができる。以下、同様に、２０〜３０パターンについて周波数特性を求める。ただしρｃは媒質密度と音速の積である。
【００１３】
一方、音響インピーダンス測定部５４は、解析パラメータ設定部５１で設定した有孔率、構造係数を有し、流れ抵抗が不明な試料（多孔質吸音材料）の音響インピーダンスを周波数毎に測定する。図３は音響インピーダンス測定部５４の構成図であり、円筒状の音響管３０を用いた２マイクロホン法（米国材料試験協会  ＡＳＴＭ  Ｅ−１０５０−９８）により多孔質吸音材料の音響インピーダンスを測定する。音響管３０の一端にはスピーカ３１が設けられ、他端にはバックプレート３２が設けられ、該バックプレート３２には試料３３が取付けられている。また、音響管３０の中間部には試料面から距離Ｌｘの位置に第１のマイク３４が設けられ、距離（Ｌｘ＋Ｄｘ）の位置に第２のマイク３５が設けられている。スピーカ３１から発せられた音が試料３３で反射して生じる定在波の音圧差を二つのマイクロホン３４，３５により測定して換算式から試料面の複素音響インピーダンスｒ＋ｊｘを測定する。この複素音響インピーダンスＺ０＝ｒ＋ｊｘは、第１のマイク３４の音圧をＰ１、第２のマイク３５の音圧をＰ２、サンプル面および第１マイク、第２マイクにおける粒子速度をＵ０，Ｕ１，Ｕ２とすれば次式
【数２】
2004347509.tif 000004
により求めることができる。ただし、
【数３】
2004347509.tif 000005
ｋ：波数（角速度ωを音速ｃで除す）、Ｓ：音響管の断面積、ρ：媒質（空気）密度、ｃ：音速である。
【００１４】
図３の測定装置は全国の公的試験機関等に広く採用されて稼動しているため、速やかに、かつ、安価に測定データを取得することができる。
スピーカ３１から発する音の周波数を変えて各周波数の音響インピーダンスＺ０＝ｒ＋ｊｘが測定されれば、吸音率・位相算出部５５は、次の（７）、（８）式
【数４】
2004347509.tif 000006
により、周波数毎に吸音率α０、位相θを演算する。これにより、試料の吸音率、位相それぞれの周波数特性を求めることができる。
【００１５】
ついで、相関演算部５６は、実測により求めた吸音率、位相の周波数特性と解析により求めた各パターンの吸音率、位相の周波数特性とを用いて、各パターンについて吸音率、位相それぞれの相関を演算する。全てのパターンについて相関演算が終了すれば相関評価部５７は、相関が最も１に近いパターンを求める。たとえば、相関評価部５７は、吸音率の相関のなかで最も１に近いパターンと位相の相関のなかで最も１に近いパターンをそれぞれ求め、▲１▼２つのパターンが同じであれば該パターンの相関が最も１に近いと判定する。しかし、▲２▼パターンが異なれば、相関が１に近い方を選択する。
【００１６】
物性値同定部５８は相関が最も１に近いパターンのパラメータ値（流れ抵抗）を試料の流れ抵抗と判定して出力する。図４は以上の流れ抵抗同定までのシーケンス説明図である。
図５は実測値と３パターンの周波数特性であり、（Ａ）は吸音率の周波数特性、（Ｂ）は位相の周波数特性である。それぞれにおいて、Ａは音響インピーダンスの測定により得られた周波数特性、Ｂ〜Ｄは流れ抵抗が１０ｋ，２０ｋ，３０ｋのときの解析により得られた周波数特性である。
流れ抵抗が１０ｋ，２０ｋ，３０ｋの吸音率周波数特性の相関値及び位相周波数特性の相関値は図６に示す通りである。これより、吸音率の相関のなかで最も１に近いパターン（流れ抵抗＝２０ｋ）と位相の相関のなかで最も１に近いパターン（流れ抵抗＝２０ｋ）が一致しており、試料の流れ抵抗は２０ｋ（＝２００００）と同定できる。
【００１７】
（Ｃ）第２実施例
第１実施例では相関（相関係数）に基づいて物性値（流れ抵抗）を同定したが、実測データと解析結果の吸音率・位相−周波数特性の相関を求めると、図７に示すように流れ抵抗値に拘わらず相関が１に近い値となり、最も相似と判断可能なピーク点を見つけることが難しくなり、流れ抵抗を同定できなくなることがある。尚、図７は厚さ２０ｍｍのウレタンの流れ抵抗・相関特性であり、Ａ，Ｂ，Ｃは、流れ抵抗がそれぞれ１６ｋ，２３ｋ，４１ｋ（実測値）の特性であり、いずれの場合も相関ピーク点が同定しにくくなっている。
そこで、第２実施例はイノベーション評価法を用いて物性値（たとえば流れ抵抗）を同定する。
【００１８】
図８は第２実施例における多孔質吸音材料の物性値同定装置の構成図であり、図２の第１実施例と同一部分には同一符号を付している。
カルマンフィルタ６１は、解析により求めた吸音率あるいは位相の周波数特性及び音響インピーダンス測定により得られた吸音率あるいは位相の周波数特性を用いて状態ベクトルｘｋを算出する。
解析により求めた吸音率あるいは位相の周波数特性におけるｎ個のポイントの値をＨｎ−１，…，Ｈ２，Ｈ１，Ｈ０とし、これらｎ個の要素からなるベクトルＨを次式
Ｈ＝［Ｈｎ−１，Ｈｎ−２，Ｈ  ｎ−３，・・，Ｈ０］
で表現する。同様に音響インピーダンス測定により得られた吸音率あるいは位相の周波数特性におけるｎ個のポイントの値をＪｎ−１，…，Ｊ２，Ｊ１，Ｊ０とするとき、ベクトルＪを次式
Ｊ＝［Ｊｎ−１，Ｊｎ−２，Ｊ  ｎ−３，・・，Ｊ０］
で表現する。以上のように解析データＨ，実測データＪを表現すると、実測データとＪｋと解析データＨｋを関係づけるシステム方程式は状態ベクトルｘｋを用いて下式で表される。ただし、Ｖｋは観測雑音である。
ｘｋ＋１＝ｘｋ                          （状態方程式）      （９）
Ｊｋ＝Ｈｋｘｋ＋Ｖｋ            （観測方程式）      （１０）
【００１９】
この状態ベクトルｘｋを下式によって逐次推定する。この推定演算を行うものがカルマンフィルタである。
【数５】
2004347509.tif 000007
ただし、Ｐｋは誤差共分散、Ｋｋはカルマンゲインである。また添え字のｋ及びｋ−１は現時点と１ステップ前の過去のデータを表す。ｋ−１｜ｋ−１は１ステップ前の過去までのデータによってｋ−１時点のデータを推定したことを表している。また、＾は推定を意味する。カルマンフィルタはｘｋ｜ｋを逐次演算して出力する。一般的なカルマンフィルタではｘｋ｜ｋを求めた時点で処理が終了する。しかし、本発明ではｘｋ｜ｋを利用してイノベーションを演算する。（１１）式の右辺括弧内は実測値と直前までの推定値の差であり、これをイノベーションという。
【００２０】
【数６】
2004347509.tif 000008
カルマンフィルタによる推定が行なわれた場合、イノベーションνｋは白色性を持つ（＝自己相関係数がΔ関数となる）ため、これを評価することで実測データと解析データの関係の強さが評価可能である。
【００２１】
従って、イノベーション算出部６２は（１４）式によりイノベーション系列νｋを算出し、自己相関算出部６３は次式
【数７】
2004347509.tif 000009
により自己相関Ｒ（ｋ）を算出する。又、相関算出部６４は、自己相関Ｒ（ｋ）とΔ関数との相関Ｃｒを次式
【数８】
2004347509.tif 000010
により算出する。尚、ＤｊはΔ関数で、ｊ＝１で１、それ以外は０となる関数である。バーは平均値を意味し、Δ関数の平均値はｎが大きくなるほど小さくなる。
【００２２】
以後、上記の状態ベクトルｘｋ｜ｋの算出、イノベーション系列νｋの算出、自己相関Ｒ（ｋ）の算出、相関Ｃｒの算出を、解析により求めた全パターンの周波数特性に適用する。
相関評価部６５は、全てのパターンについて上記演算が完了すれば、相関Ｃｒが最も１に近いパターンを求める。物性値同定部６８は該パターンのパラメータ値（流れ抵抗）を試料の流れ抵抗と判定して出力する。
【００２３】
図９は発泡ウレタン、グラスウールなどの代表的な多孔質吸音材料の吸音率と位相のそれぞれの実測データと解析データ（流れ抵抗を複数パターン用意して解析）を用いて、第１実施例の相関係数及び第２実施例のイノベーション白色性を求めた。そして、第１実施例において、相関が最も１に近い時の流れ抵抗を求めて図９中に示し、第２実施例では、イノベーション白色性が最も強い時の流れ抵抗を求めて示している。
図９より、第１実施例では位相相関のピーク点が求まらない場合があるが、第２実施例ではピーク点はすべて求まっている。なお、第２実施例において、理想的には、吸音率と位相の両方においてイノベーション白色性が最も強い時の流れ抵抗値が一致することであるが、異なる場合は、イノベーション白色性が強い方の流れ抵抗を選択する。
【００２４】
（Ｄ）物性値同定の効果
（１）車室内音場解析
物性値（たとえば流れ抵抗）を正確に同定できれば、車室壁面に発泡ウレタンを貼付した時の車室内音場（音の強弱の分布状態）がどのようになるかを実機を製作せずに評価可能となる。さらに、本手法は多孔質吸音材料に限らず、他の車室内装材にも適用可能なので、今後車室内装材料全般について物性値を同定することでより正確な車室内音場の評価が可能であると考えられる。図１０は運転席における音圧―周波数特性（実測データＡと解析結果Ｂ）であり、解析結果Ｂは実測データＡに類似していることがわかる。
【００２５】
（２）スピーカ設置位置の決定
車室内音場を正確に評価することが可能となれば、たとえば、車室内のどの箇所にスピーカを設置すると最適な聴取環境が得られるかを判断可能となる。図１１はセダン車の運転席におけるスピーカ位置変更前Ａと位置変更後Ｂの音圧−周波数特性であり、実機を製作しなくても、特性Ｂのようにフラットな周波数特性となるスピーカ位置を決定することができる。
以上本発明を多孔質吸音材料の物性値を測定する場合について説明したが、本発明は多孔質吸音材料に限らず多孔質以外の吸音材料の物性値の同定にも用いることできる。
又、以上では流れ抵抗をパラメータとして可変して吸音材料の流れ抵抗を求めたが、吸音材料の他の物性値も同様に求めることができる。
【００２６】
【発明の効果】
以上本発明によれば、専用の測定器を用いないで吸音材料の物性値を測定することができ、安価に吸音材料の物性値を同定することができる。
また、本発明における垂直入射吸音率測定は公的試験機関等で広く行なわれおり、専用機器を所有していなくとも安価に測定が行なえる。また、測定方法はＪＩＳ等国際規格で規定されており、測定の信頼性は高い。したがって、物性値を精度よく同定することができる。
また、本発明によれば、吸音率測定データを元にイノベーション評価によって物性値を決定することから安価、短時間かつより正確に物性値を決定することが可能となる。
【図面の簡単な説明】
【図１】本発明の概略説明図である。
【図２】第１実施例における多孔質吸音材料の物性値同定装置の構成図である。
【図３】音響インピーダンス測定部の構成図である。
【図４】流れ抵抗同定までのシーケンス説明図である。
【図５】実測値と３パターンの周波数特性である。
【図６】流れ抵抗が１０ｋ，２０ｋ，３０ｋの吸音率周波数特性及び位相周波数特性の相関値説明図表である。
【図７】第１実施例の問題点説明図である。
【図８】第２実施例における多孔質吸音材料の物性値同定装置の構成図である。
【図９】第１実施例において相関が最も１に近い時の流れ抵抗および第２実施例においてイノベーション白色性が最も強い時の流れ抵抗を示す図表である。
【図１０】運転席における音圧―周波数特性である。
【図１１】スピーカ位置変更前Ａと位置変更後Ｂの音圧−周波数特性である。
【符号の説明】
１１音響解析部
１２音響インピーダンス測定部
１３相関演算部
１４物性値同定部
２１音響解析部
２２音響インピーダンス測定部
２３イノベーション演算部
２４物性値同定部
【図１】
2004347509.tif 000011
【図２】
2004347509.tif 000012
【図３】
2004347509.tif 000013
【図４】
2004347509.tif 000014
【図５】
2004347509.tif 000015
【図６】
2004347509.tif 000016
【図７】
2004347509.tif 000017
【図８】
2004347509.tif 000018
【図９】
2004347509.tif 000019
【図１０】
2004347509.tif 000020
【図１１】
2004347509.tif 000021