(57)【要約】      （修正有）【課題】多孔質材料の変形前における材料パラメータを用いて、当該多孔質材料の変形後における材料パラメータを推定する際の推定精度を向上させる。【解決手段】材料パラメータ推定装置５は、多孔質材料の変形前における多孔度φ（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付部５１と、前記多孔質材料の変形後における多孔度φ（Ｐ）を推定する推定部５２とを備えている。【選択図】図３
2020170000.tif 000002
【特許請求の範囲】
【請求項１】
多孔質材料の変形前における多孔度φ（１）と、前記変形に伴う前記多孔質材料の体積
変化率Ｐとの入力を受け付ける受付部と、
前記多孔質材料の変形後における多孔度φ（Ｐ）を、
【数１】
2020170000.tif 000022
により推定する推定部と
を備えた材料パラメータ推定装置。
【請求項２】
前記受付部がさらに、前記多孔質材料の変形前における流れ抵抗σ（１）の入力を受け
付けるものであり、
前記推定部がさらに、前記多孔質材料の変形後における流れ抵抗σ（Ｐ）を、前記変形
が２次元的である場合は
【数２】
2020170000.tif 000023
により推定し、前記変形が３次元的である場合は
【数３】
2020170000.tif 000024
により推定するものである、請求項１に記載の材料パラメータ推定装置。
【請求項３】
多孔質材料の変形前における迷路度α∞（１）と、前記変形に伴う前記多孔質材料の体
積変化率Ｐとの入力を受け付ける受付部と、
前記多孔質材料の変形後における迷路度α∞（Ｐ）を、
【数４】
2020170000.tif 000025
により推定する推定部と
を備えた材料パラメータ推定装置。
【請求項４】
多孔質材料の変形前における多孔度φ（１）及び粘性特性長Λ（１）と、前記変形に伴
う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付部と、
前記多孔質材料の変形後における粘性特性長Λ（Ｐ）を、
【数５】
2020170000.tif 000026
により推定する推定部と
を備えた材料パラメータ推定装置。
【請求項５】
多孔質材料の変形前における多孔度φ（１）及び熱的特性長Λ’（１）と、前記変形に
伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付部と、
前記多孔質材料の変形後における熱的特性長Λ’（Ｐ）を、
【数６】
2020170000.tif 000027
により推定する推定部と
を備えた材料パラメータ推定装置。
【請求項６】
多孔質材料の変形前における多孔度φ（１）と、前記変形に伴う前記多孔質材料の体積
変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における多孔度φ（Ｐ）を、
【数７】
2020170000.tif 000028
により推定する推定ステップと
を含む材料パラメータ推定方法。
【請求項７】
前記受付ステップにおいてさらに、前記多孔質材料の変形前における流れ抵抗σ（１）
の入力が受け付けられ、
前記推定ステップにおいてさらに、前記多孔質材料の変形後における流れ抵抗σ（Ｐ）
が、前記変形が２次元的である場合は
【数８】
2020170000.tif 000029
により推定され、前記変形が３次元的である場合は
【数９】
2020170000.tif 000030
により推定される、請求項６に記載の材料パラメータ推定方法。
【請求項８】
多孔質材料の変形前における迷路度α∞（１）と、前記変形に伴う前記多孔質材料の体
積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における迷路度α∞（Ｐ）を、
【数１０】
2020170000.tif 000031
により推定する推定ステップと
を含む材料パラメータ推定方法。
【請求項９】
多孔質材料の変形前における多孔度φ（１）及び粘性特性長Λ（１）と、前記変形に伴
う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における粘性特性長Λ（Ｐ）を、
【数１１】
2020170000.tif 000032
により推定する推定ステップと
を含む材料パラメータ推定方法。
【請求項１０】
多孔質材料の変形前における多孔度φ（１）及び熱的特性長Λ’（１）と、前記変形に
伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における熱的特性長Λ’（Ｐ）を、
【数１２】
2020170000.tif 000033
により推定する推定ステップと
を含む材料パラメータ推定方法。
【請求項１１】
多孔質材料の変形前における多孔度φ（１）と、前記変形に伴う前記多孔質材料の体積
変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における多孔度φ（Ｐ）を、
【数１３】
2020170000.tif 000034
により推定する推定ステップと
をコンピュータに実行させるための材料パラメータ推定プログラム。
【請求項１２】
前記受付ステップにおいてさらに、前記多孔質材料の変形前における流れ抵抗σ（１）
の入力が受け付けられ、
前記推定ステップにおいてさらに、前記多孔質材料の変形後における流れ抵抗σ（Ｐ）
が、前記変形が２次元的である場合は
【数１４】
2020170000.tif 000035
により推定され、前記変形が３次元的である場合は
【数１５】
2020170000.tif 000036
により推定される、請求項１１に記載の材料パラメータ推定プログラム。
【請求項１３】
多孔質材料の変形前における迷路度α∞（１）と、前記変形に伴う前記多孔質材料の体
積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における迷路度α∞（Ｐ）を、
【数１６】
2020170000.tif 000037
により推定する推定ステップと
をコンピュータに実行させるための材料パラメータ推定プログラム。
【請求項１４】
多孔質材料の変形前における多孔度φ（１）及び粘性特性長Λ（１）と、前記変形に伴
う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における粘性特性長Λ（Ｐ）を、
【数１７】
2020170000.tif 000038
により推定する推定ステップと
をコンピュータに実行させるための材料パラメータ推定プログラム。
【請求項１５】
多孔質材料の変形前における多孔度φ（１）及び熱的特性長Λ’（１）と、前記変形に
伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における熱的特性長Λ’（Ｐ）を、
【数１８】
2020170000.tif 000039
により推定する推定ステップと
をコンピュータに実行させるための材料パラメータ推定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、材料パラメータ推定装置、材料パラメータ推定方法、及び材料パラメータ推
定プログラムに関する。
【背景技術】
【０００２】
自動車、ＯＡ機器、家電製品、住宅といった様々な製品の付加価値の一つとして、「静
粛性」が求められるようになってきている。この静粛性を高めるためには、吸音、遮音と
いった作用を持つ材料が用いられることが多い。このような材料は音響材料と呼ばれてい
る。音響材料は、単一材料である場合もあれば、複数の材料を層状に積み重ねた積層材料
である場合もある。
【０００３】
ここで、吸音とは、ある材料に入射した音が反射せずに当該材料に吸収される、または
透過することをいう。吸音の性能は吸音率という値で評価される。また、遮音とは、ある
材料に入射した音が当該材料を通り抜けないことをいう。遮音の性能は音響透過損失とい
う値で評価される。吸音率及び音響透過損失はいずれも周波数の関数である。吸音率と音
響透過損失とをまとめて音響性能と呼ぶ。
【０００４】
音響材料を開発する場合には、実際に試作してみて、当該試作材料の音響性能を実測し
てその測定結果が所望のものかどうかを確認し、所望の結果が得られるまで試作と実測と
を繰り返すことが一般的である。ただ、試作を複数回行うことになるため材料のコストが
かさむ。また、試作材料の音響性能を実測するためには、専用の設備が必要となる。さら
に、試作と実測とを繰り返すため、開発期間が長期化することが多い。
【０００５】
そこで、試作及び実測を繰り返すという手法に代えて、音響材料を数理モデルとして数
学的に表現し、その数理モデルにしたがって当該音響材料の音響性能を計算するという手
法がとられることもある。この手法の場合、モデル化しようとする材料の特性、例えば当
該材料の密度、厚さといった材料パラメータの値を決め、それらの値を用いて当該材料の
音響性能を計算する。数理モデルにはいくつかあり、モデル化の対象となる材料の種類に
よって、使用する数理モデルが異なる。また、数理モデルによって必要となる材料パラメ
ータの種類も異なる。もちろん、音響材料が積層材料である場合には、材料ごとに異なる
数理モデルを用いることができる。特許文献１には、数理モデルの一つであるビオモデル
が記載されている。
【０００６】
音響材料の例として、繊維系多孔質材料が挙げられる。繊維系多孔質材料の骨格となる
繊維をそれ単体で見たときに、ポリエステル等の高分子材料や綿などの繊維は非常に柔ら
かい。その一方で、ガラス繊維や岩綿などの人造鉱物繊維等で構成されている繊維は非常
に硬い。ゆえに、多孔質材料全体の弾性は骨格の弾性に少なからず影響される。
【０００７】
この多孔質材料の弾性が非常に柔らかいとみなせるときは、音波が多孔質材料に入射す
ると該材料が音波によって励振され音波と同じように振動し、骨格の振動が多孔質材料の
音響性能に与える影響は小さいと考えられる。他方、多孔質材料の弾性が硬いとみなせる
ときは、該材料に音波が入射しても骨格は容易に振動せず、その骨格の振動を無視できる
と考えられる。
【０００８】
このように多孔質材料の骨格が柔らかい場合や硬い場合、多孔質材料の弾性を無視して
モデル化することがある。特に繊維系の多孔質材料は弾性を無視しても差し支えないこと
が多い。骨格の硬い多孔質材料の弾性を無視したモデルを提唱したのがＪｏｈｎｓｏｎ氏
、Ｃｈａｍｐｏｕｘ氏、Ａｌｌａｒｄ氏に代表される研究者たちである。この弾性を無視
したモデルは、３人の名前の頭文字をとって「ＪＣＡモデル」と呼ばれている。
【０００９】
骨格が硬いために多孔質材料全体の弾性を無視したモデルは広く剛骨格モデル（rigid 
frame model）と呼ばれている。ＪＣＡモデルは剛骨格モデルに属する。また、骨格が柔
らかいために弾性を無視したモデルは柔骨格モデル（limp frame model）と呼ばれている
。
【００１０】
ＪＣＡモデルにおいては、（１）多孔度、（２）流れ抵抗、（３）迷路度、（４）粘性
特性長、（５）熱的特性長の５つの材料パラメータが用いられる。これらの材料パラメー
タは、音響材料の流体特性、すなわち、空気伝搬音に関するものである。具体的には、多
孔度は、音響材料中の空気の占める割合である。空隙の占める割合が大きいほど音響材料
の多孔度は大きくなる。次に、流れ抵抗は、音響材料中の空気の流れにくさを表す数値で
ある。音は空気の振動であるため、音響材料の流れ抵抗が大きい場合には、当該材料内で
は空気が流れにくくなり、つまり、当該材料は音を伝えにくい材料であると言える。この
流れ抵抗は、材料パラメータの中でも、非常に大きなウェイトをしめるものである。また
、迷路度は、空隙により生まれる空気の経路の形状の複雑さを表す数値である。迷路度が
大きければ、それだけ音が材料中を長く伝わるため、より吸音されることになる。粘性特
性長及び熱的特性長は、それぞれ、粘性損失及び熱交換損失の大きさを表す数値である。
【００１１】
上述のＪＣＡモデルにより多孔質材料の音響性能を推定するには、それぞれに適切な５
つの材料パラメータが必要となる。この材料パラメータは、一般的に、同一の材料であっ
ても成型加工などの圧縮等による変形が生じた場合には、全ての材料パラメータの値が変
わってしまう。したがって、音響性能を推定したい多孔質材料の状態全てに対してその都
度材料パラメータを求め直さねばならず、非常な労力を伴うこととなる。
【００１２】
そのため、多孔質材料の変形による材料パラメータの値の変化を容易に推定する手段が
望まれる。そのような手段の例が、非特許文献１及び２に記載されている。非特許文献１
では、５つ全ての材料パラメータについて、多孔質材料が圧縮される場合の推定が試みら
れている。また、非特許文献２では、流れ抵抗について検討が加えられている。
【先行技術文献】
【特許文献】
【００１３】
【特許文献１】国際公開第２０１３／０８０９９３号
【非特許文献】
【００１４】
【非特許文献１】Bernard Castagnede, Achour Aknine, Bruno Brouard, Viggo Tarnow, Effects of compression on the sound absorption of fibrous materials., Applied Acoustics, 61(2), 173-182 (2000).
【非特許文献２】Viggo Tarnow, Airflow resistivity of models of fibrous acoustic materials., J. Acoust. Soc. Am. 100 (6), 3706-3713 (1996).
【発明の概要】
【発明が解決しようとする課題】
【００１５】
非特許文献２には、流れ抵抗以外の材料パラメータについては言及されておらず、かつ
適用可能な多孔質材料の形状が限られていると考えられる。また、非特許文献１における
推定も、適用範囲が狭くかつ検証が不十分であると考えられる。
【００１６】
本発明は、かかる不都合を解消し、多孔質材料の変形前における材料パラメータを用い
て、当該多孔質材料の変形後における材料パラメータを推定する際の推定精度を向上させ
ることを目的とする。
【課題を解決するための手段】
【００１７】
上記の目的を達成するために、本発明に係る材料パラメータ推定においては、まず、多
孔質材料の変形前における多孔度φ（１）と、前記変形に伴う前記多孔質材料の体積変化
率Ｐとの入力が受け付けられる。そして、前記多孔質材料の変形後における多孔度φ（Ｐ
）が、
【数１】
2020170000.tif 000003
により推定される。
【００１８】
さらに、前記多孔質材料の変形前における流れ抵抗σ（１）の入力が受け付けられ、前
記多孔質材料の変形後における流れ抵抗σ（Ｐ）が、前記変形が２次元的である場合は
【数２】
2020170000.tif 000004
により推定され、前記変形が３次元的である場合は
【数３】
2020170000.tif 000005
により推定されるとすることもできる。
【００１９】
材料パラメータ推定の別の形態においては、まず、多孔質材料の変形前における迷路度
α∞（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐとの入力が受け付けられる
。そして、前記多孔質材料の変形後における迷路度α∞（Ｐ）が、
【数４】
2020170000.tif 000006
により推定される。
【００２０】
材料パラメータ推定の別の形態においては、まず、多孔質材料の変形前における多孔度
φ（１）及び粘性特性長Λ（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐとの
入力が受け付けられる。そして、前記多孔質材料の変形後における粘性特性長Λ（Ｐ）が
、
【数５】
2020170000.tif 000007
により推定される。
【００２１】
材料パラメータ推定の別の形態においては、まず、多孔質材料の変形前における多孔度
φ（１）及び熱的特性長Λ’（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐと
の入力が受け付けられる。そして、前記多孔質材料の変形後における熱的特性長Λ’（Ｐ
）が、
【数６】
2020170000.tif 000008
により推定される。
【発明の効果】
【００２２】
本発明によれば、多孔質材料の変形後において材料パラメータを再測定せずに、多孔質
材料の変形前における材料パラメータを用いて、当該多孔質材料の変形後における材料パ
ラメータを推定することができるとともに、その推定精度を向上させることができる。
【図面の簡単な説明】
【００２３】
【図１】（Ａ）多孔質材料の断面図である。（Ｂ）変形後における多孔質材料の断面図である。
【図２】変形後の多孔度と体積変化率との関係を表すグラフである。
【図３】材料パラメータ推定装置の機能構成例を示す説明図である。
【図４】材料パラメータ推定装置のコンピュータハードウエア構成例を示す説明図である。
【図５】変形後の多孔度と迷路度との関係を表すグラフである。
【図６】変形後の多孔度と熱的特性長との関係を表すグラフである。
【図７】変形後の多孔度と粘性特性長との関係を表すグラフである。
【図８】変形後の多孔度と流れ抵抗との関係を表すグラフである。
【発明を実施するための形態】
【００２４】
図１に、繊維系の多孔質材料の断面を示している。多孔質材料１は、骨格２と、骨格２
の間の空隙３とを有している。多孔質材料１が、成型加工等による圧縮又は膨張といった
変形をする場合、骨格２の間の空隙３が変形することが多く、骨格２自体が変形すること
は稀である。そこで、多孔質材料の変形に伴って空隙のみが変形すると仮定する。
【００２５】
変形前の多孔質材料において、骨格が占める体積をＶｆ、空隙が占める体積をＶａ、全
体の体積をＶＴとすると、以下のような関係がある。
ＶＴ＝Ｖｆ＋Ｖａ
【００２６】
そして、多孔質材料の変形後における全体の体積ＶＴ（Ｐ）は、以下のように表される
。
ＶＴ（Ｐ）：＝ＰＶＴ＝Ｖｆ＋ＰＶａ
ただし、Ｐは、変形に伴う多孔質材料の体積変化率である。Ｐ＞１のとき体積は増加し、
Ｐ＜１のとき体積は減少する。また、変形後における空隙の体積Ｖａ（Ｐ）は以下のよう
に表される。
Ｖａ（Ｐ）：＝ＰＶａ
【００２７】
以上を踏まえた本発明の実施形態を以下に説明する。ただし、本発明は、以下の実施形
態によって限定されるものではない。
【００２８】
［多孔度］
多孔質材料の変形前における多孔度φ（１）は、その定義より、以下のように表される
。
【数７】
2020170000.tif 000009
である。そして、多孔質材料の変形後における多孔度φ（Ｐ）は、以下のように表される
。
【数８】
2020170000.tif 000010
【００２９】
このように、変形前の多孔度φ（１）と変形に伴う体積変化率Ｐとによって、変形後の
多孔度φ（Ｐ）を求めることができる。また、体積変化率Ｐが未知であって、変形前の多
孔度φ（１）と変形後の多孔度φ（Ｐ）とが与えられている場合は、式（２）から、以下
のように体積変化率Ｐを求めることができる。
【数９】
2020170000.tif 000011
【００３０】
図２に、変形前の多孔度φ（１）＝０．９９における、変形後の多孔度φ（Ｐ）と体積
変化率Ｐとの関係を示している。
【００３１】
図３に、多孔度を含む材料パラメータの、変形後における値を推定する材料パラメータ
推定装置５の機能構成例を示している。材料パラメータ推定装置５は、受付部５１と推定
部５２とを備えている。受付部５１及び推定部５２の機能は後述する。
【００３２】
図４は、材料パラメータ推定装置５のコンピュータハードウエア構成例を示している。
材料パラメータ推定装置５は、ＣＰＵ５１０と、インタフェース装置５２０と、表示装置
５３０と、入力装置５４０と、ドライブ装置５５０と、補助記憶装置５６０と、メモリ装
置５７０とを備えており、これらがバス５８０により相互に接続されている。
【００３３】
材料パラメータ推定装置５の機能を実現するプログラムは、ＣＤ−ＲＯＭ等の記録媒体
５９０によって提供される。プログラムを記録した記録媒体５９０がドライブ装置５５０
にセットされると、プログラムが記録媒体５９０からドライブ装置５５０を介して補助記
憶装置５６０にインストールされる。あるいは、プログラムのインストールは必ずしも記
録媒体５９０により行う必要はなく、ネットワークを介して他のコンピュータからダウン
ロードすることもできる。補助記憶装置５６０は、インストールされたプログラムを格納
すると共に、必要なファイルやデータ等を格納する。
【００３４】
メモリ装置５７０は、プログラムの起動指示があった場合に、補助記憶装置５６０から
プログラムを読み出して格納する。ＣＰＵ５１０は、メモリ装置５７０に格納されたプロ
グラムにしたがって材料パラメータ推定装置５の機能を実現する。インタフェース装置５
２０は、ネットワークを通して他のコンピュータに接続するためのインタフェースとして
用いられる。表示装置５３０はプログラムによるＧＵＩ（Graphical User Interface）等
を表示する。入力装置５４０はキーボード及びマウス等である。
【００３５】
多孔度の推定にあたり、材料パラメータ推定装置５内の受付部５１は、変形前の多孔度
φ（１）の入力と、変形に伴う体積変化率Ｐの入力とを受け付ける。そして、推定部５２
は、式（２）により、変形後の多孔度φ（Ｐ）を推定する。
【００３６】
なお、変形前の多孔度φ（１）は、実測値及び計算値のいずれでもよい。
【００３７】
［迷路度］
発明者は、コンピュータによる数値流体解析（以下２次元ＣＦＤと表す）を事前に行っ
た。具体的には、複素変数境界要素法（ＣＶＢＥＭ）を用いた、非圧縮性非粘性流体の２
次元理想流体定常流解析を行った。その結果、変形後の迷路度α∞（Ｐ）は、変形後の多
孔度φ（Ｐ）に対する１次式で表現できることを見いだした。この１次式の傾きをａとし
、式（１）も用いると、以下のようにして、変形後の迷路度α∞（Ｐ）を求めることがで
きる。
【数１０】
2020170000.tif 000012
ただし、α∞（１）は、変形前の多孔度である。
【００３８】
図５に、４種の繊維径の多孔質材料に対する、変形後の迷路度α∞（Ｐ）と変形後の多
孔度φ（Ｐ）との関係を示している。繊維径とは、図１（Ａ）の符号Ｄとして示している
ように、円形断面を有する骨格２の直径である。
【００３９】
図５において、丸印（○）の点は、２次元ＣＦＤにより得られた、繊維径１０μｍの場
合の迷路度を表している。四角形（□）の点は、２次元ＣＦＤにより得られた、繊維径５
μｍの場合の迷路度を表している。三角形（△）の点は、２次元ＣＦＤにより得られた、
繊維径２．５μｍの場合の迷路度を表している。ひし形（◇）の点は、２次元ＣＦＤによ
り得られた、繊維径１μｍの場合の迷路度を表している。なお、同図において、丸印の点
と四角形の点と三角形の点とは、ひし形の点と重なっているため、明瞭に表されていない
。
【００４０】
同じく図５において、グラフＧ１１は、繊維径１０μｍの場合の、式（４）を用いた、
迷路度の推定値を表している。グラフＧ１２は、繊維径５μｍの場合の、式（４）を用い
た、迷路度の推定値を表している。グラフＧ１３は、繊維径２．５μｍの場合の、式（４
）を用いた、迷路度の推定値を表している。グラフＧ１４は、繊維径１μｍの場合の、式
（４）を用いた、迷路度の推定値を表している。
【００４１】
このように、変形後の迷路度α∞（Ｐ）は、変形後の多孔度φ（Ｐ）に対する１次式で
表現することができる。変形後の多孔度φ（Ｐ）が０．９以上の場合においては特に、推
定精度が高いことがわかる。
【００４２】
迷路度の推定にあたり、材料パラメータ推定装置５内の受付部５１は、変形前の迷路度
α∞（１）の入力と、変形に伴う体積変化率Ｐの入力とを受け付ける。そして、推定部５
２は、式（４）により、変形後の迷路度α∞（Ｐ）を推定する。
【００４３】
なお、変形前の迷路度α∞（１）は、実測値及び計算値のいずれでもよい。
【００４４】
［熱的特性長］
変形前における熱的特性長をΛ’（１）は、その定義より、以下のように表される。
【数１１】
2020170000.tif 000013
ただし、Ｓｆは骨格の表面積を表しているとともに、ＶＴ（１）＝ＶＴである。つまり、
熱的特性長とは、空隙の体積と骨格の表面積の比の２倍である。変形後の多孔質材料の熱
的特性長Λ’（Ｐ）は、式（２）と式（５）から、以下のように表すことができる。
【数１２】
2020170000.tif 000014
【００４５】
このように、変形前の熱的特性長Λ’（１）が既知であれば、体積変化率Ｐと変形前に
おける多孔度φ（１）から、変形後の熱的特性長Λ’（Ｐ）を求めることができる。
【００４６】
図６に、円形断面の４種の繊維径の多孔質材料における、変形後の熱的特性長Λ’（Ｐ
）と変形後の多孔度φ（Ｐ）との関係を示している。丸印（○）の点は、２次元ＣＦＤに
より得られた、繊維径１０μｍの熱的特性長を表している。四角形（□）の点は、２次元
ＣＦＤにより得られた、繊維径５μｍの熱的特性長を表している。三角形（△）の点は、
２次元ＣＦＤにより得られた、繊維径２．５μｍの熱的特性長を表している。ひし形（◇
）の点は、２次元ＣＦＤにより得られた、繊維径１μｍの熱的特性長を表している。なお
、２次元ＣＦＤにおいては、迷路度と同様に、複素変数境界要素法（ＣＶＢＥＭ）を用い
た、非圧縮性非粘性流体の２次元理想流体定常流解析を行った。
【００４７】
同じく図６において、グラフＧ２１は、繊維径１０μｍの場合の、式（６）を用いた、
熱的特性長の推定値を表している。グラフＧ２２は、繊維径５μｍの場合の、式（６）を
用いた、熱的特性長の推定値を表している。グラフＧ２３は、繊維径２．５μｍの場合の
、式（６）を用いた、熱的特性長の推定値を表している。グラフＧ２４は、繊維径１μｍ
の場合の、式（６）を用いた、熱的特性長の推定値を表している。
【００４８】
このように、４種の繊維径のいずれにおいても、式（６）により変形後の熱的特性長を
精度良く推定できることがわかる。
【００４９】
熱的特性長の推定にあたり、材料パラメータ推定装置５内の受付部５１は、変形前の熱
的特性長Λ’（１）の入力と、変形前の多孔度φ（１）の入力と、変形に伴う体積変化率
Ｐの入力とを受け付ける。そして、推定部５２は、式（６）により、変形後の熱的特性長
Λ’（Ｐ）を推定する。
【００５０】
なお、変形前の熱的特性長Λ’（１）及び多孔度φ（１）は、実測値及び計算値のいず
れでもよい。
【００５１】
［粘性特性長］
２次元ＣＦＤの結果、変形後の粘性特性長Λ（Ｐ）は、変形後の多孔度φ（Ｐ）に対し
て熱的特性長とほぼ同じ変化をすることを発明者は見いだした。そこで、式（６）を参考
にして、変形後の粘性特性長Λ（Ｐ）は以下のように表される。
【数１３】
2020170000.tif 000015
ただし、Λ（１）は、変形前の粘性特性長である。
【００５２】
このように、変形前の粘性特性長Λ（１）が既知であれば、体積変化率Ｐと変形前の多
孔度φ（１）から、変形後の粘性特性長Λ（Ｐ）を求めることができる。
【００５３】
図７に、円形断面の４種の繊維径の多孔質材料における、変形後の粘性特性長Λ（Ｐ）
と変形後の多孔度φ（Ｐ）との関係を示している。丸印（○）の点は、２次元ＣＦＤによ
り得られた、繊維径１０μｍの粘性特性長を表している。四角形（□）の点は、２次元Ｃ
ＦＤにより得られた、繊維径５μｍの粘性特性長を表している。三角形（△）の点は、２
次元ＣＦＤにより得られた、繊維径２．５μｍの粘性特性長を表している。ひし形（◇）
の点は、２次元ＣＦＤにより得られた、繊維径１μｍの粘性特性長を表している。
【００５４】
同じく図７において、グラフＧ３１は、繊維径１０μｍの場合の、式（７）を用いた、
粘性特性長の推定値を表している。グラフＧ３２は、繊維径５μｍの場合の、式（７）を
用いた、粘性特性長の推定値を表している。グラフＧ３３は、繊維径２．５μｍの場合の
、式（７）を用いた、粘性特性長の推定値を表している。グラフＧ３４は、繊維径１μｍ
の場合の、式（７）を用いた、粘性特性長の推定値を表している。
【００５５】
このように、４種の繊維径のいずれにおいても、式（７）により変形後の粘性特性長を
精度良く推定することが可能である。
【００５６】
粘性特性長の推定にあたり、材料パラメータ推定装置５内の受付部５１は、変形前の粘
性特性長Λ（１）の入力と、変形前の多孔度φ（１）の入力と、変形に伴う体積変化率Ｐ
の入力とを受け付ける。そして、推定部５２は、式（７）により、変形後の粘性特性長Λ
（Ｐ）を推定する。
【００５７】
なお、変形前の粘性特性長Λ（１）及び多孔度φ（１）は、実測値及び計算値のいずれ
でもよい。
【００５８】
［流れ抵抗］
流れ抵抗を測定するときの多孔質材料中を流れる空気の流速は極めて遅く、そのときの
レイノルズ数は１よりも十分に小さい。そのため、多孔質材料中の空気の流れを表現する
方程式として、ナビエ−ストークス（Ｎａｖｉｅｒ−Ｓｔｏｋｅｓ）方程式における慣性
項を無視したストークス近似を採用する。この近似におけるストークス方程式は、以下の
ように表される。
【数１４】
2020170000.tif 000016
【００５９】
ただし、μは空気の粘度であり、ｕは多孔質材料の表面から裏面に向かって流れる空気
の、該多孔質材料中の流速である。多孔質材料に浸入する前の与えられた速度をＵ、その
多孔質材料の多孔度をφとすると、ｕ＝Ｕ／φと表すことができる。また、ｐは圧力を表
し、ｇｒａｄは方向微分である。そのため、圧力ｐの方向微分（ｇｒａｄ  ｐ）は、多孔
質材料の表裏面における圧力差を表すと考えることができる。すなわち、以下の通りであ
る。
【数１５】
2020170000.tif 000017
ただし、Ｐ１は多孔質材料の表面における圧力を表し、Ｐ２は多孔質材料の裏面における
圧力を表し、Ｌは多孔質材料の表面から裏面までの長さである。
【００６０】
さらに、式（８）におけるΔは方向の２階微分を表す。したがって、式（８）が示すと
ころは、圧力差ｇｒａｄ  ｐは、多孔質材料中の流速ｕに比例し、長さＬの２乗に逆比例
すると解釈できる。すなわち、以下の通りである。
【数１６】
2020170000.tif 000018
【００６１】
他方、変形前の流れ抵抗σ（１）は、その定義から、式（９）も考慮すると、以下のよ
うに表される。
【数１７】
2020170000.tif 000019
【００６２】
次に、変形後の流れ抵抗を考えるのだが、変形に伴って体積がＰ倍になるということは
、２次元的な変形を考えると、正方形の多孔質材料の一辺の長さは√Ｐ倍となる。ゆえに
変形後の流れ抵抗σ（Ｐ）は、式（１０）を参照して、以下のように表される。
【数１８】
2020170000.tif 000020
【００６３】
３次元的な変形（前後左右上下から力が加わる）の場合は、立方体の多孔質材料の一辺
の長さはＰの立方根倍となる。そのため、３次元的な変形後の流れ抵抗σ（Ｐ）は、以下
のように表される。
【数１９】
2020170000.tif 000021
【００６４】
すなわち、式（２）も考えると、変形前の流れ抵抗σ（１）が既知であるならば、体積
変化率Ｐと変形前の多孔度φ（１）から、変形後の流れ抵抗σ（Ｐ）を求めることができ
る。
【００６５】
図８に、円形断面の４種の繊維径の多孔質材料における、変形後の流れ抵抗σ（Ｐ）と
変形後の多孔度φ（Ｐ）との関係を示している。丸印（○）の点は、２次元ＣＦＤにより
得られた、繊維径１０μｍの流れ抵抗を表している。四角形（□）の点は、２次元ＣＦＤ
により得られた、繊維径５μｍの流れ抵抗を表している。三角形（△）の点は、２次元Ｃ
ＦＤにより得られた、繊維径２．５μｍの流れ抵抗を表している。ひし形（◇）の点は、
２次元ＣＦＤにより得られた、繊維径１μｍの流れ抵抗を表している。なお、２次元ＣＦ
Ｄにおいては、非圧縮性粘性流体の２次元定常流解析を行った。
【００６６】
同じく図８において、グラフＧ４１は、繊維径１０μｍの場合の、式（１１）を用いた
、流れ抵抗の推定値を表している。グラフＧ４２は、繊維径５μｍの場合の、式（１１）
を用いた、流れ抵抗の推定値を表している。グラフＧ４３は、繊維径２．５μｍの場合の
、式（１１）を用いた、流れ抵抗の推定値を表している。グラフＧ４４は、繊維径１μｍ
の場合の、式（１１）を用いた、流れ抵抗の推定値を表している。
【００６７】
このように、４種の繊維径のいずれにおいても、式（１１）により変形後の流れ抵抗を
精度良く推定することができる。変形後の多孔度φ（Ｐ）が０．９以上の場合においては
特に、推定精度が高い。
【００６８】
流れ抵抗の推定にあたり、材料パラメータ推定装置５内の受付部５１は、変形前の流れ
抵抗σ（１）の入力と、変形前の多孔度φ（１）の入力と、変形に伴う体積変化率Ｐの入
力とを受け付ける。そして、推定部５２は、式（２）により、変形後の多孔度φ（Ｐ）を
推定するとともに、式（１１）により、変形後の流れ抵抗σ（Ｐ）を推定する。
【００６９】
なお、変形前の流れ抵抗σ（１）及び多孔度φ（１）は、実測値及び計算値のいずれで
もよい。
【００７０】
以上のような材料パラメータ推定の手法によれば、多孔質材料の変形後において材料パ
ラメータを再測定せずに、多孔質材料の変形前における材料パラメータを用いて、当該多
孔質材料の変形後における材料パラメータを推定することができるとともに、その推定精
度を向上させることができる。
【００７１】
［その他］
上記５つの材料パラメータは、ＪＣＡモデルだけではなく、剛骨格モデル、柔骨格モデ
ル、ビオモデルなどの他の数理モデルでも用いられる。前述した材料パラメータ推定の手
法は、５つの材料パラメータを用いる任意の数理モデルに適用可能である。また、本手法
が対象とする音響材料は、繊維系多孔質材料を含む任意の多孔質材料とすることができる
。多孔質材料中の骨格の分布状態に関わらず、また、骨格の断面形状が円形であるか否か
に関わらず、本手法を用いることができる。
【００７２】
前述した材料パラメータ推定の手法は、材料パラメータ推定装置としての側面だけでは
なく、材料パラメータ推定方法及び材料パラメータ推定プログラムとしての側面をも有し
ている。
【００７３】
前述した材料パラメータ推定装置の機能的構成及び物理的構成は、前述の態様に限られ
るものではなく、例えば、各機能や物理資源を統合して実装したり、逆に、さらに分散し
て実装したりすることも可能である。
【符号の説明】
【００７４】
１      多孔質材料
２      骨格
３      空隙
５      材料パラメータ推定装置
５１    受付部
５２    推定部
【図１】
2020170000.tif 000040
【図２】
2020170000.tif 000041
【図３】
2020170000.tif 000042
【図４】
2020170000.tif 000043
【図５】
2020170000.tif 000044
【図６】
2020170000.tif 000045
【図７】
2020170000.tif 000046
【図８】
2020170000.tif 000047
【手続補正書】
【提出日】令和2年7月8日(2020.7.8)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
多孔質材料の変形前における多孔度φ（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付部と、
前記多孔質材料の変形後における多孔度φ（Ｐ）を、
【数１】
2020170000.tif 000048
により推定する推定部と
を備えた材料パラメータ推定装置。
【請求項２】
多孔質材料の変形前における迷路度α∞（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付部と、
前記多孔質材料の変形後における迷路度α∞（Ｐ）を、
【数２】
2020170000.tif 000049
により推定する推定部と
を備えた材料パラメータ推定装置。
【請求項３】
多孔質材料の変形前における多孔度φ（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における多孔度φ（Ｐ）を、
【数３】
2020170000.tif 000050
により推定する推定ステップと
を含む材料パラメータ推定方法。
【請求項４】
多孔質材料の変形前における迷路度α∞（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における迷路度α∞（Ｐ）を、
【数４】
2020170000.tif 000051
により推定する推定ステップと
を含む材料パラメータ推定方法。
【請求項５】
多孔質材料の変形前における多孔度φ（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における多孔度φ（Ｐ）を、
【数５】
2020170000.tif 000052
により推定する推定ステップと
をコンピュータに実行させるための材料パラメータ推定プログラム。
【請求項６】
多孔質材料の変形前における迷路度α∞（１）と、前記変形に伴う前記多孔質材料の体積変化率Ｐとの入力を受け付ける受付ステップと、
前記多孔質材料の変形後における迷路度α∞（Ｐ）を、
【数６】
2020170000.tif 000053
により推定する推定ステップと
をコンピュータに実行させるための材料パラメータ推定プログラム。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００１８
【補正方法】削除
【補正の内容】
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００２０
【補正方法】削除
【補正の内容】
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００２１
【補正方法】削除
【補正の内容】