2022 Fiscal Year Research-status Report
Analysis on fluid-structure-acoustic interaction of vocal fold vibration and sound generation
| Project/Area Number |
20K14648
|
|---|
| Research Institution | Toyohashi University of Technology |
|---|
Principal Investigator |
吉永 司 豊橋技術科学大学, 工学(系)研究科(研究院), 助教 (50824190)
|
|---|
| Project Period (FY) |
2020-04-01 – 2024-03-31
|
| Keywords | 声帯 / 流体構造連成解析 / ２質点モデル / 剥離 / 空力音 / ベルヌーイの定理 |
|---|
| Outline of Annual Research Achievements |
本年度は主に，２つの質点をバネとダンパーで接続した単純声帯モデルに対して，流体－構造－音響連成解析を行うことにより，提案手法の妥当性を検証した．声帯ひだの振動解析に関しては質点とバネに単純化することで計算効率を上げ，これまで種々の解析が行われてきた単純形状に対する流体解析を行うことにより，計算精度を検証しつつ，音まで十分な精度で解析できていることを確認しながら解析を進めた．気流解析の精度については，M5モデルと呼ばれる単純形状モデルに対する先行研究の表面圧力分布計測結果を用いて，剥離点まで表面圧力が十分な精度で計算できていることを確認した．その後，左右の声帯ひだの弾性率を変化させて左右非対称の振動が生まれた際に，3次元的な気流が，どの程度振動特性に影響を与えるのか調べた．この時，先行研究でよく用いられる1次元のベルヌーイの定理を用いたモデルと比較した．その結果，3次元気流解析と1次元モデルの結果はほとんど同じ傾向を示し，声帯の振動様式から発生音の声質まで，1次元気流モデルでも上気道の気流の圧縮性を考慮すれば十分な精度で解析できることがわかった．また，これらの結果は，声帯の振動様式に関して，気流の３次元性よりも声帯の生体組織をどのようにモデル化するのかが重要であることを示唆しており，今後は声帯ひだの複雑な形状や非線形的な変形様態に注力して解析を進めていく．得られた研究成果は2回の国際会議で発表し，米音響学会誌の学術論文として掲載された．
|
| Current Status of Research Progress |
Current Status of Research Progress
2: Research has progressed on the whole more than it was originally planned.
Reason
今年度目標にしていた3次元の気流解析と1次元モデルとの比較及び評価が無事終わり，論文として出版することができた．当初は気流の３次元性により，声帯の左右非対称性がより強く表れると考えていたが，左右のひだの弾性率の比が1.6程度の差では，３次元および１次元解析の結果はほとんど同じとなり，１次元解析でも十分な精度で解析できていることが示せた．この結果は，極端な気流の影響が表れない領域では，声帯組織の変形特性が声帯の振動特性を予測する上で重要であることを示唆している．そのため，今後は声帯ひだ周りの靭帯や軟骨，筋線維組織に注力してモデル化を進めていく．また，３次元解析ではスーパーコンピュータを用いて１週間程度時間がかかるのに対し，１次元解析では一般的なPCを用いて数分で計算が終わるため，流体解析を実施したことの無い研究者でも解析が活用できる．
これらの成果は，音声生成解析に関する議論が活発に行われている米音響学会及び声帯に関するシミュレーション研究が多く議論されるInternational Conference on Advances in Quantitative Laryngology, Voice and Speech Researchにて発表し，多くの研究者に認知してもらうことができた．また，共同研究先であるUCLAのZhang教授の研究室を訪問し，超弾性体を用いた声帯の非線形変形解析について，これまで実施されてきた手法を学び，来年度から進めていく方針について，より具体的なプランを立てた．
|
| Strategy for Future Research Activity |
声帯ひだに関して，超弾性体を用いた有限要素解析を実施する．この際，これまでの先行研究で用いられてきた声帯の筋線維と靭帯を分けて解析するbody-coverモデルを参考にし，有限要素モデルを構築する．振動解析に関しては，固有値解析を応用したモード法を活用し，３次元気流解析とカップリングして高速に解析できるようにする．モード法では，特定の周波数までの振動モードを考慮し，予め固有値解析結果を蓄積し，その結果と流体力の時間変化のデータを組み合わせることで過渡応答解析を行う．
また，共同研究先のZhang教授より声帯の筋線維構造を撮像した医療画像を提供してもらい，声帯の実形状を考慮した有限要素解析を実施する．そして，声の個人差まで考慮した精度で解析を行う．さらに，声帯ひだの直上にある仮声帯および喉頭蓋による狭窄を考慮して気流解析することで，それぞれの狭窄が音の発生にどの程度影響を及ぼすのか明らかにする．
|
| Causes of Carryover |
3月末日に開催された国際会議AQL2023に出席して発表するのとともに，共同研究先のUCLAを訪問する際，4月1日を跨いだ出張となったため，次年度の予算として執行する必要があった．その出張で残額を全て使用するため，その他の経費には用いない．
|
