(57)【要約】【課題】  流れ抵抗σの測定のみによって、広い周波数範囲（１００〜５ｋＨｚ程度）において、且つ、密度の大小にかかわらずで、その音響特性を高い精度で予測することが可能な音響特性予測方法と音響特性予測装置を提供すること。【解決手段】  任意の材料の流れ抵抗σ（Ｐａ・ｓ／ｍ２）を測定し、上記測定した流れ抵抗σ（Ｐａ・ｓ／ｍ２）と常用対数を用いた算出法によって上記任意の材料の特性インピーダンスＺｃと伝搬定数γを算出し、上記特性インピーダンスＺｃと伝搬定数γを使用して上記任意の材料の音響特性を予測するようにしたもの。【選択図】  図１
2007057515.tif 000002
【特許請求の範囲】
【請求項１】
任意の材料の流れ抵抗σ（Ｐａ・ｓ／ｍ２）を測定し、
上記測定した流れ抵抗σ（Ｐａ・ｓ／ｍ２）と常用対数を用いた算出法によって上記任意の材料の特性インピーダンスＺｃと伝搬定数γを算出し、
上記特性インピーダンスＺｃと伝搬定数γを使用して上記任意の材料の音響特性を予測するようにしたことを特徴とする音響特性予測方法。
【請求項２】
請求項１記載の音響特性予測方法において、
上記特性インピーダンスＺｃと伝搬定数γは次の式（ＩＩ）、（ＩＩＩ）によって定義されるものであり、又、該式（ＩＩ）、（ＩＩＩ）における変数Ｒ、Ｘ、α、βは常用対数を用いた算出法によって算出されるものであることを特徴とする音響特性予測方法。
ＺＣ＝ρ０・ｃ０（Ｒ＋ｊＸ）―――（ＩＩ）
γ＝α＋ｊＢ              ―――（ＩＩＩ）
但し、
ρ０    ：空気の密度〔ｋｇ／ｍ３〕
ｃ０    ：音速〔ｍ／ｓ〕
Ｒ    ：特性インピーダンスＺｃの実数部  
Ｘ    ：特性インピーダンスのＺｃの虚数部
α    ：伝搬定数の実数部〔ｎｅｐｅｒｓ／ｍ〕
β    ：伝搬定数の虚数部〔ｒａｄｉａｎ／ｓ〕
【請求項３】
請求項２記載の音響特性予測方法において、
上記変数Ｒ、Ｘ、α、βは次の式（ＶＩＩＩ）〜（ＸＩ）によって定義されるものであることを特徴とする音響特性予測方法。
2007057515.tif 000014
但し、
f      ：周波数(Ｈz)
ω    ：角周波数（radian／s）
ａ    ：定数
ｂ    ：定数  
ｃ    ：定数
ｃ０   ：音速（ｍ／ｓ）
ｄ    ：定数
ｐ    ：定数
ｑ    ：定数
ｒ    ：定数
ｓ    ：定数
【請求項４】
請求項１〜請求項３の何れかに記載の音響特性予測方法において、
上記特性インピーダンスＺｃと伝搬定数γと上記材料の厚さｂと室温ｔとから所望の音響特性を算出するようにしたことを特徴とする音響特性予測方法。
【請求項５】
請求項４記載の音響特性予測方法において、
上記所望の音響特性は周波数毎の音響インピーダンスＺ１と垂直入射吸音率α1であり、夫々次に示す式（ＸＩＩ）、（ＸＩＩＩ）によって算出するようにしたことを特徴とする音響特性予測方法。
2007057515.tif 000015
但し、
Ｚ２    ：材料背面の比音響インピーダンス
b       ：材料の厚み
Ｚｃ    ：特性インピーダンス
γ      ：伝搬定数
2007057515.tif 000016
但し、
Ｚ０      ：空気の特性インピーダンス
【請求項６】
請求項１〜請求項５の何れかに記載の音響特性予測方法において、
上記任意の材料は吸音材料であることを特徴とする音響特性予測方法。
【請求項７】
任意の材料に関して測定された流れ抵抗σ（Ｐａ・ｓ／ｍ２）と常用対数を用いた算出法によって上記任意の材料の特性インピーダンスＺｃと伝搬定数γを算出する特性インピーダンス・伝搬定数算出手段と、
上記特性インピーダンス・伝搬定数算出手段によって算出された特性インピーダンスＺｃと伝搬定数γを使用して上記任意の材料の音響特性を算出する音響特性算出手段と、
を具備したことを特徴とする音響特性予測装置。
【請求項８】
請求項７記載の音響特性予測装置において、
上記特性インピーダンスＺｃと伝搬定数γは次の式（ＩＩ）、（ＩＩＩ）によって定義されるものであり、又、該式（ＩＩ）、（ＩＩＩ）における変数Ｒ、Ｘ、α、βは常用対数を用いた算出法によって算出されるものであることを特徴とする音響特性予測装置
ＺＣ＝ρ０・ｃ０（Ｒ＋ｊＸ）―――（ＩＩ）
γ＝α＋ｊＢ              ―――（ＩＩＩ）
但し、
ρ０    ：空気の密度〔ｋｇ／ｍ３〕
ｃ０    ：音速〔ｍ／ｓ〕
Ｒ    ：特性インピーダンスＺｃの実数部  
Ｘ    ：特性インピーダンスのＺｃの虚数部
α    ：伝搬定数の実数部〔ｎｅｐｅｒｓ／ｍ〕
β    ：伝搬定数の虚数部〔ｒａｄｉａｎ／ｓ〕
【請求項９】
請求項８記載の音響特性予測装置において、
上記変数Ｒ、Ｘ、α、βは次の式（ＶＩＩＩ）〜（ＸＩ）によって定義されるものであることを特徴とする音響特性予測装置。
2007057515.tif 000017
但し、
f      ：周波数(Ｈz)
ω    ：角周波数（radian／s）
ａ    ：定数
ｂ    ：定数  
ｃ    ：定数
ｃ０   ：音速（ｍ／ｓ）
ｄ    ：定数
ｐ    ：定数
ｑ    ：定数
ｒ    ：定数
ｓ    ：定数
【請求項１０】
請求項８又は請求項９記載の音響特性予測装置において、
上記特性インピーダンスＺｃと伝搬定数γと上記材料の厚さｂと室温ｔとから所望の音響特性を算出するようにしたことを特徴とする音響特性予測装置。
【請求項１１】
請求項１０記載の音響特性予測装置において、
上記所望の音響特性は周波数毎の音響インピーダンスＺ１と垂直入射吸音率α1であり、夫々次に示す式（ＸＩＩ）、（ＸＩＩＩ）によって算出するようにしたことを特徴とする音響特性予測装置。
2007057515.tif 000018
但し、
Ｚ２    ：材料背面の比音響インピーダンス
b       ：材料の厚み
Ｚｃ    ：特性インピーダンス
γ      ：伝搬定数
2007057515.tif 000019
但し、
Ｚ０      ：空気の特性インピーダンス
【請求項１２】
請求項８〜請求項１１の何れかに記載の音響特性予測装置において、
上記任意の材料は吸音材料であることを特徴とする音響特性予測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、吸音材料等の任意の材料の音響特性、例えば、周波数毎の比音響インピーダンスＺ１や垂直入射吸音率α１を、上記材料の流れ抵抗σ（Ｐａ・ｓ／ｍ２）を使用して算出する音響特性予測方法と音響特性予測装置に係り、特に、その予測精度を向上させることができるように工夫したものに関する。
【背景技術】
【０００２】
例えば、吸音材料の性能を示す吸音率は、音の入射角度によって「垂直入射吸音率」、「斜入射吸音率」、「ランダム入射吸音率」の３種類に分類される。そのうち「垂直入射吸音率」の測定は、音響インピーダンス管による測定方法（「ＪＩＳ Ａ １４０５」や「ＩＳＯ １０５３４」）によって行われる。この場合、測定装置が比較的小規模で、且つ、試料も少量で済むため、多孔質吸音材料等の品質管理や研究開発に応用されている。
但し、正確な測定を行うためには、吸音材料のカッティングや音響インピーダンス管内での吸音材料の支持の方法等において注意すべき点が多い。
【０００３】
ところで、音は空気振動の伝達であるから、空気が通り難いところでは伝わり難くなる。上記流れ抵抗σはその流れ難さを表す評価量であり吸音材料の音響特性を把握するための最も重要な物性値の一つである。
【０００４】
又、多孔質の吸音材料の吸音率が流れ抵抗σと関係が深いことは以前から知られていた。上記流れ抵抗σの測定は「垂直入射吸音率」の測定に比べて測定原理が簡単であるが、反面、微風速や微差圧の測定が難しいという欠点もあり、これまで研究開発や品質管理の現場でそれほど利用されてこなかった。これに対して、最近では、流量計や差圧センサの性能が上がり、コンピュータを使った計測システムを構築することにより、比較的簡単な装置によって精度の高い流れ抵抗σを測定することができるようになっている。
【０００５】
そして、流れ抵抗σという物性値から吸音材料の「垂直入射吸音率」を精度良く予測できれば、材料開発においてどのような工夫をすれば吸音性能を向上させることが可能かについてより明確に確認することができ、それによって、開発コストを大幅に縮減することができるものである。
【０００６】
又、流れ抵抗σのみから「垂直入射吸音率」を予測するときに必要となる式は過去に２つ発表されている。一つは１９７０年にイギリスの「Ｄｅｌａｎｙ and Ｂａｚｌｅｙ」が発表した式であり、他の一つは１９９０年に「幹」が発表した式である。これらの式を使用することによりグラスウール等の繊維系多孔質材料の吸音率を比較的容易に、且つ、高い精度で予測することができる。
【０００７】
ところで、流れ抵抗σは単位空気流れに対する単位厚さ当たりの圧力低下の形で表される物性値であり、「ＩＳＯ９０５３」により規定されている。そして、流れ抵抗σの測定原理は次の通りである。
まず、被測定試料を取り付けた管内に空気を流す。この時、試料が空気の流れを妨げるので、試料の前後に圧力差が発生する。その時、空気の流れの速さと圧力差を測定することによって流れ抵抗σは求められる。すなわち、流れ抵抗σ（Ｐａ・ｓ／ｍ２）は次の式（Ｉ）によって定義される。
σ＝△ｐ／ｄ・ｖ          ―――（Ｉ）    
但し、
△ｐ  ：材料前後に発生する差圧（Ｐａ）  
ｄ    ：材料厚さ（ｍ）
ｖ    ：材料中を流れる空気の流速（ｍ／ｓ）
である。
【０００８】
又、「垂直入射吸音率」は、吸音材料の特性インピーダンスＺｃ、伝搬定数γの他、吸音材料の厚さｂと室温ｔから容易に算出できる。特性インピーダンスＺｃと伝搬定数γは吸音材料の素材固有の音響特性を示す複素数であり、次の式（ＩＩ）、（ＩＩＩ）によって表される。
ＺＣ＝ρ０ｃ０（Ｒ＋ｊＸ）  ―――（ＩＩ）
γ＝α＋ｊＢ              ―――（ＩＩＩ）
但し、
ρ０    ：空気の密度〔ｋｇ／ｍ３〕
ｃ０    ：音速〔ｍ／ｓ〕
Ｒ    ：特性インピーダンスＺｃの実数部  
Ｘ    ：特性インピーダンスのＺｃの虚数部
α    ：伝搬定数の実数部〔ｎｅｐｅｒｓ／ｍ〕
β    ：伝搬定数の虚数部〔ｒａｄｉａｎ／ｓ〕
である。
【０００９】
そして、流れ抵抗σという物性値から、特性インピーダンスＺｃと伝搬定数γを精度良く予測できれば、流れ抵抗σから「垂直入射吸音率」を高い精度で予測できることになる。流れ抵抗σから特性インピーダンスＺｃと伝搬定数γを予測する式が、「Ｄｅｌａｎｙ  and  Ｂａｚｌｅｙ」と「幹」の２式である。
【００１０】
「Ｄｅｌａｎｙ and   Ｂａｚｌｅｙ」と「幹」の予測式では、式（ＩＩ）、（ＩＩＩ）の変数Ｒ、Ｘ、α、βは夫々（ｆ／σ）の関数として式（ＩＶ）、（Ｖ）、（ＶＩ）、（ＶＩＩ）の形で表される。又、係数と次数の値はそれぞれ「表１」、「表２」に示す通りである。「表１」が「Ｄｅｌａｎｙ and   Ｂａｚｌｅｙ」の場合であり、「表2」が「幹」の場合である。
【００１１】
2007057515.tif 000003
但し、
f      ：周波数(Ｈz)
ω    ：角周波数（radian／s）
ａ    ：定数
ｂ    ：定数  
ｃ    ：定数
ｃ０   ：音速（ｍ／ｓ）
ｄ    ：定数
ｐ    ：定数
ｑ    ：定数
ｒ    ：定数
ｓ    ：定数
【００１２】
【表１】
2007057515.tif 000004
【表２】
2007057515.tif 000005
【００１３】
尚、上記「Ｄｅｌａｎｙ and   Ｂａｚｌｅｙ」と「幹」の予測式に関しては、非特許文献１、非特許文献2に開示されている。
【００１４】
【非特許文献１】Ｍ．Ｅ．Ｄｅｌａｎｄｙ ｅｔ  ａｌ．,  Ａｐｐｌ．Ａｃｏｕｓｔ. ３，  １０５−１１６  （１９７０）
【非特許文献２】Ｍｉｋｉ．, Ｊ．  Ａｃｏｕｓｔ．  Ｓｏｃ．  Ｊｐｎ. (E), １１．１，  １９−２４  （１９９０）
【発明の開示】
【発明が解決しようとする課題】
【００１５】
上記従来の構成によると次のような問題があった。
すなわち、「Ｄｅｌａｎｙ  and Ｂａｚｌｅｙ」と「幹」の式から、特性インピーダンスＺｃ、伝搬定数γ、「垂直入射吸音率」等の予測を行ってみると、広い周波数範囲での予測を網羅できないという問題があった。又、密度の低い吸音材料や密度が高い吸音材料に対して予測誤算が大きいという問題点があった。このため吸音材料の正確な吸音率を知るためには、「垂直入射吸音率」の測定が必要不可欠であった。
これに対しては、精度の高い予測方法があれば「垂直入射吸音率」の測定を行うことなく流れ抵抗σのみの測定から吸音材料の吸音率等必要な音響特性を 確認することができるわけであるが、従来そのような予測方法は実現されていなかった。
【００１６】
本発明はこのような点に基づいてなされたものでその目的とするところは、流れ抵抗σの測定のみによって、広い周波数範囲（１００〜５ｋＨｚ程度）において、且つ、密度の大小にかかわらずで、その音響特性を高い精度で予測することが可能な音響特性予測方法と音響特性予測装置を提供することにある。
【課題を解決するための手段】
【００１７】
上記目的を達成するべく本願発明の請求項１による音響特性予測方法は、任意の材料の流れ抵抗σ（Ｐａ・ｓ／ｍ２）を測定し、上記測定した流れ抵抗σ（Ｐａ・ｓ／ｍ２）と常用対数を用いた算出法によって上記任意の材料の特性インピーダンスＺｃと伝搬定数γを算出し、上記特性インピーダンスＺｃと伝搬定数γを使用して上記任意の材料の音響特性を予測するようにしたことを特徴とするものである。
又、請求項２による音響特性予測方法は、請求項１記載の音響特性予測方法において、上記特性インピーダンスＺｃと伝搬定数γは次の式（ＩＩ）、（ＩＩＩ）によって定義されるものであり、又、該式（ＩＩ）、（ＩＩＩ）における変数Ｒ、Ｘ、α、βは常用対数を用いた算出法によって算出されるものであることを特徴とするものである。
ＺＣ＝ρ０・ｃ０（Ｒ＋ｊＸ）―――（ＩＩ）
γ＝α＋ｊＢ              ―――（ＩＩＩ）
又、請求項３による音響特性予測方法は、請求項２記載の音響特性予測方法において、上記変数Ｒ、Ｘ、α、βは次の式（ＶＩＩＩ）〜（ＸＩ）によって定義されるものであることを特徴とするものである。
2007057515.tif 000006
又、請求項４による音響特性予測方法は、請求項１〜請求項３の何れかに記載の音響特性予測方法において、上記特性インピーダンスＺｃと伝搬定数γと上記材料の厚さｂと室温ｔとから所望の音響特性を算出するようにしたことを特徴とするものである。
又、請求項５による音響特性予測方法は、請求項４記載の音響特性予測方法において、上記所望の音響特性は周波数毎の音響インピーダンスＺ１と垂直入射吸音率α1であり、夫々次に示す式（ＸＩＩ）、（ＸＩＩＩ）によって算出するようにしたことを特徴とするものである。
2007057515.tif 000007
又、請求項６による音響特性予測方法は、請求項１〜請求項５の何れかに記載の音響特性予測方法において、上記任意の材料は吸音材料であることを特徴とするものである。
又、請求項７による音響特性予測方法は、任意の材料に関して測定された流れ抵抗σ（Ｐａ・ｓ／ｍ２）と常用対数を用いた算出法によって上記任意の材料の特性インピーダンスＺｃと伝搬定数γを算出する特性インピーダンス・伝搬定数算出手段と、上記特性インピーダンス・伝搬定数算出手段によって算出された特性インピーダンスＺｃと伝搬定数γを使用して上記任意の材料の音響特性を算出する音響特性算出手段と、を具備したことを特徴とするものである。
又、請求項８による音響特性予測方法は、請求項７記載の音響特性予測装置において、上記特性インピーダンスＺｃと伝搬定数γは次の式（ＩＩ）、（ＩＩＩ）によって定義されるものであり、又、該式（ＩＩ）、（ＩＩＩ）における変数Ｒ、Ｘ、α、βは常用対数を用いた算出法によって算出されるものであることを特徴とするものである。
ＺＣ＝ρ０・ｃ０（Ｒ＋ｊＸ）―――（ＩＩ）
γ＝α＋ｊＢ              ―――（ＩＩＩ）
又、請求項９による音響特性予測方法は、請求項８記載の音響特性予測装置において、上記変数Ｒ、Ｘ、α、βは次の式（ＶＩＩＩ）〜（ＸＩ）によって定義されるものであることを特徴とするものである。
2007057515.tif 000008
又、請求項１０による音響特性予測方法は、請求項８又は請求項９記載の音響特性予測装置において、上記特性インピーダンスＺｃと伝搬定数γと上記材料の厚さｂと室温ｔとから所望の音響特性を算出するようにしたことを特徴とするものである。
又、請求項１１による音響特性予測方法は、請求項１０記載の音響特性予測装置において、上記所望の音響特性は周波数毎の音響インピーダンスＺ１と垂直入射吸音率α1であり、夫々次に示す式（ＸＩＩ）、（ＸＩＩＩ）によって算出するようにしたことを特徴とするものである。
2007057515.tif 000009
又、請求項１２による音響特性予測方法は、請求項８〜請求項１１の何れかに記載の音響特性予測装置において、上記任意の材料は吸音材料であることを特徴とするものである。
尚、各式における変数、定数は特許請求の範囲において定義した通りのものである。
【発明の効果】
【００１８】
以上本発明による音響特性予測方法と音響特性予測装置によると、流れ抵抗σの測定のみによって、吸音材料の音響特性を高い精度で予測することができるようになった。これは、常用対数を使用した算出方法を採用したことによるものである。
一つの物性値である流れ抵抗σのみから材料の吸音率等の音響特性を予測できるようになったので、例えば、吸音材料の開発に要するコストを大幅に低減させることが可能となり、又、開発に要する時間を大幅に短縮させることができるものである。
【発明を実施するための最良の形態】
【００１９】
以下、本発明の一実施の形態を図１乃至図８を参照して説明する。図１は本実施の形態による音響特性予測方法の概略を説明するための図である。図１に示すように、まず、流れ抵抗測定部１において「ＩＳＯ９０５３」に準じた流れ抵抗σの測定が実行される。上記流れ抵抗σは、従来例の説明で挙げた式（Ｉ）によって規定されるものである。
又、特性インピーダンスＺｃと伝搬定数γの予測部２で、測定された流れ抵抗σの値を使って、周波数毎の特性インピーダンスＺｃと伝搬定数γの予測がコンピュータ等の演算手段によって実行される。この特性インピーダンスＺｃと伝搬定数γの予測部２で使う式が、従来例の説明で挙げた式（ＩＩ）、（ＩＩＩ）と示す式（ＶＩＩＩ）〜（ＸＩ）である。
【００２０】
2007057515.tif 000010
但し、
f      ：周波数(Ｈz)
ω    ：角周波数（radian／s）
ａ    ：定数
ｂ    ：定数  
ｃ    ：定数
ｃ０   ：音速（ｍ／ｓ）
ｄ    ：定数
ｐ    ：定数
ｑ    ：定数
ｒ    ：定数
ｓ    ：定数
【００２１】
又、上記式（ＶＩＩＩ）〜（ＸＩ）において、吸音材料がグラスウールの場合には、上記各定数は次の表３に示すようなものとなる。
【００２２】
【表３】
2007057515.tif 000011
【００２３】
上記式（ＶＩＩＩ）〜（ＸＩ）に示すように、この実施の形態では、変数Ｒ、Ｘ、α、βを算出するに際して常用対数を使用している。これによって、広い周波数範囲（１００〜５ｋＨｚ程度）において、且つ、密度の大小にかかわらずで、吸音材料の音響特性を高い精度で予測することが可能になる。
【００２４】
そして、既に述べたように、従来例の説明で挙げている式（ＩＩ）、式（ＩＩＩ）と上記式（ＶＩＩＩ）〜式（ＸＩ）を使用して必要な演算を行なう。それによって、特性インピーダンスＺｃと伝搬定数γを算出する。それら特性インピーダンスＺｃと伝搬定数γの各変数Ｒ、Ｘ、α、βを図２〜図５に示す。
まず、図２では、変数Ｒと（ｆ／σ）との関係を示していて、従来方式である「Ｄ＆Ｂモデル」の場合を破線で示し、同じく従来方式である「Ｍｉｋｉも出る」を一点鎖線で示し、本実施の形態の場合、すなわち、「Ｎｅｗモデル」を実線で示している。同様に、図３では、変数Ｘと（ｆ／σ）との関係を示していて、従来方式である「Ｄ＆Ｂモデル」の場合を破線で示し、同じく従来方式である「Ｍｉｋｉも出る」を一点鎖線で示し、本実施の形態の場合、すなわち、「Ｎｅｗモデル」を実線で示している。同様に、図４では、変数αと（ｆ／σ）との関係を示していて、従来方式である「Ｄ＆Ｂモデル」の場合を破線で示し、同じく従来方式である「Ｍｉｋｉも出る」を一点鎖線で示し、本実施の形態の場合、すなわち、「Ｎｅｗモデル」を実線で示している。同様に、図５では、変数βと（ｆ／σ）との関係を示していて、従来方式である「Ｄ＆Ｂモデル」の場合を破線で示し、同じく従来方式である「Ｍｉｋｉも出る」を一点鎖線で示し、本実施の形態の場合、すなわち、「Ｎｅｗモデル」を実線で示している。
これら図２〜図５から明らかなように、従来例の場合には各変数Ｒ、Ｘ、α、βの変化が直線的なものとなっているが、本実施の形態の場合には各定数Ｒ、Ｘ、α、βの変化が曲線的なものとなっている。それによって、広い周波数範囲（１００〜５ｋＨｚ程度）において、且つ、密度の大小にかかわらずで、吸音材料の音響特性を高い精度で予測することが可能になる。
尚、図１における情報処理はコンピュータ等の演算手段を用いて行う。
【００２５】
次に、上記得られた特性インピーダンスＺｃ、伝搬定数γと、吸音材料の厚さｂ及び室温ｔから吸音材料の吸音率を求める。その工程を図６のブロック図に示す。図６に示すように、予測された特性インピーダンスＺｃと伝搬定数γに、吸音材料の厚さｂと室温ｔの入力部３から吸音材料の厚さｂと室温ｔの２つのデータを追加して、比音響インピーダンスＺ１の算出部４、垂直入射吸音率α１の算出部６の処理を順に行う。
尚、図６の一連の情報処理はコンピュータ等の演算手段を用いて行う。
【００２６】
上記比音響インピーダンスＺ１の算出部４は次に示す式（ＸＩＩ）に示すような演算を実行するものであり、又、上記垂直入射吸音率α１の算出部６は次に示す式（ＸＩＩＩ）に示すような演算を実行するものである。
2007057515.tif 000012
但し、
Ｚ２    ：材料背面の比音響インピーダンス
b       ：材料の厚み
Ｚｃ    ：特性インピーダンス
γ      ：伝搬定数    
2007057515.tif 000013
但し、
Ｚ０      ：空気の特性インピーダンス
である。
尚、上記空気の特性インピーダンスＺ０の値は室温ｔによって変化するものである。
【００２７】
以上のように、流れ抵抗σ、吸音材料の厚さｂ及び室温ｔから、比音響インピーダンスＺ１と垂直入射吸収率α１を予測することができる。
実際に流れ抵抗σからグラスウールの垂直入射吸収率α１を予測して、実測値と比較した結果を図７、図８に示す。
図７は、流れ抵抗σが６６００Ｐａ・ｓ／ｍ２，厚さ５３．５ｍｍの試料について室温２４．１℃での結果である。図７中「Ｄ＆Ｂモデル」と「Ｍｉｋｉモデル」は従来方式による予測結果を示す線図であり、「Ｄ＆Ｂモデル」を破線で示すと共に「Ｍｉｋｉモデル」を一点鎖線で示している。又、「Ｍｅａｓｕｒｅｄ」は実測値を示す線図であり二点鎖線で示している。又、「Ｎｅｗ  モデル」が本実施の形態の場合の予測結果を示しているものであり、実線で示している。
それによると、１０００Ｈｚの誤差を見ると過去の２つの予測では７％であるのに対して新しい予測では４％である。
図８は、流れ抵抗σが５３１００Ｐａ・ｓ／ｍ２，厚さ２６．５ｍｍの試料について室温２４．９℃での結果である。図８中「Ｄ＆Ｂモデル」と「Ｍｉｋｉモデル」は従来方式による予測結果を示す線図であり、「Ｄ＆Ｂモデル」を破線で示すと共に「Ｍｉｋｉモデル」を一点鎖線で示している。又、「Ｍｅａｓｕｒｅｄ」は実測値を示す線図であり二点鎖線で示している。又、「Ｎｅｗ  モデル」が本実施の形態の場合の予測結果を示しているものであり、実線で示している。
それによって、同様に１０００Ｈｚの誤差を見ると過去の予測では１０％であるのに対して新しい予測では２％である。
つまり、新しい方法によって、流れ抵抗σの低い材料から高い材料において、過去の２つの方法よりも高い精度の予測が可能となっていることがわかる。
【００２８】
尚、本発明は前記一実施の形態に限定されるものではない。
まず、前記一実施の形態において示した式（ＶＩＩＩ）〜（ＸＩ）は一例であって、常用対数を使用した他の合理的な式の使用も考えられる。
又、対象物としては、吸音材料以外にも様々なものが想定される。
又、音響特性としても比音響インピーダンスＺ１と垂直入射吸収率α１に限定されるものではなく、その他様々な特性が想定される。
【産業上の利用可能性】
【００２９】
本発明は、吸音材料の音響特性、例えば、周波数毎の音響インピーダンスと複素音圧反射係数と垂直入射吸音率を流れ抵抗σ（Ｐａ・ｓ／ｍ２）を使用して算出する音響特性予測方法と音響特性予測装置に係り、特に、その予測精度を向上させることができるように工夫したものに関し、優れた吸音材料を開発する上で有効である。
【図面の簡単な説明】
【００３０】
【図１】本発明の一実施の形態を示す図で、音響特性予測方法と音響特性予測装置を説明するための図である。
【図２】本発明の一実施の形態を示す図で、変数Ｒと周波数ｆを流れ抵抗σ（Ｐａ・ｓ／ｍ２）で除した値との関係を示す図で、従来例を細線で示し本実施の形態の場合を太線で示した特性図である。
【図３】本発明の一実施の形態を示す図で、変数Ｘと周波数ｆを流れ抵抗σ（Ｐａ・ｓ／ｍ２）で除した値との関係を示す図で、従来例を細線で示し本実施の形態の場合を太線で示した特性図である。
【図４】本発明の一実施の形態を示す図で、変数αと周波数ｆを流れ抵抗σ（Ｐａ・ｓ／ｍ２）で除した値との関係を示す図で、従来例を細線で示し本実施の形態の場合を太線で示した特性図である。
【図５】本発明の一実施の形態を示す図で、変数βと周波数ｆを流れ抵抗σ（Ｐａ・ｓ／ｍ２）で除した値との関係を示す図で、従来例を細線で示し本実施の形態の場合を太線で示した特性図である。
【図６】本発明の一実施の形態を示す図で、流れ抵抗σ（Ｐａ・ｓ／ｍ２）と吸音材料の厚さｂと室温ｔから吸音率を求める手順を示すブロック図である。
【図７】本発明の一実施の形態を示す図で、本実施の形態による吸音率の周波数特性の予測値と従来例による予測値と実測値とを比較して示す特性図である。
【図８】本発明の一実施の形態を示す図で、本実施の形態による吸音率の周波数特性の予測値と従来例による予測値と実測値とを比較して示す特性図である。
【符号の説明】
【００３１】
１  流れ抵抗測定部
２  特性インピーダンスと伝播定数の予測部
３  材料厚さと室温の入力部
４  比音響インピーダンスの算出部
６  垂直入射吸音率の算出部
【図１】
2007057515.tif 000020
【図２】
2007057515.tif 000021
【図３】
2007057515.tif 000022
【図４】
2007057515.tif 000023
【図５】
2007057515.tif 000024
【図６】
2007057515.tif 000025
【図７】
2007057515.tif 000026
【図８】
2007057515.tif 000027