JP2013202119A - 嚥下シミュレーション装置及び方法 - Google Patents

Abstract

【課題】嚥下に関する実現象を近似的に再現することが容易で、経口摂取品の物性及び挙動に係る物理量を定量化できる嚥下シミュレーション装置を提供する。
【解決手段】本発明による嚥下シミュレーション装置１００Ａは、頭頸部器官からなる動的三次元頭頸部モデルを形成する頭頸部モデリング部１０と、頭頸部モデルにおける各頭頸部器官の運動を設定する器官運動設定部３０と、解析対象としての経口摂取品及びその物性を設定する経口摂取品物性設定部４０と、経口摂取品をモデル化した擬似経口摂取品を口腔に入力する入力部８１と、頭頸部モデルにおける各頭頸部器官の運動と擬似経口摂取品の嚥下時の挙動を、粒子法を用いて三次元空間で解析する運動解析部５０と、運動解析部５０で解析された各頭頸部器官の運動と擬似経口摂取品の嚥下時の挙動の解析結果を動画面に表示する表示部８２とを備える。
【選択図】図２

Description

本発明は嚥下シミュレーション装置及び方法に関する。詳しくは、口腔ならびに咽喉を通過する際の流体・食塊の挙動について粒子法を用いて解析を行う嚥下シミュレーション装置、動的三次元頭頸部モデル、嚥下シミュレーション方法、コンピュータ読み取り可能なプログラム、経口摂取品開発支援装置、経口摂取品開発方法、経口摂取品製造方法、食育支援装置、食育方法及び診断支援装置に関する。

嚥下動作、特に嚥下時の食品物性と頭頸部器官の運動との関係は複雑であり、現象そのものを正確に把握することは非常に困難である。しかし医療や介護の現場では、高齢者や障害者の誤嚥、誤飲を防ぐため、さまざまな試行錯誤を繰り返し、誤嚥・誤飲リスクの低減に努めている。また、昨今のこんにゃくゼリーによる窒息事故を踏まえ、食品全般において客観的な数値や指標をもって飲食品の安全性を保証することが必要となってきた。

嚥下現象の解明には、嚥下造影や筋電位測定などの生体情報を直接入手する方法と、嚥下ロボットや数値シミュレーションなどを用いて間接的に情報を入手する方法の二つがある。これらは、被験者に多大な負担がかかる、ロボットの頭頸部器官の挙動や構造を変更するのが大掛かりになるなど、多くの問題を抱えていた。

また、これまで生体内での流体・固形物など食塊の挙動に関する数値解析が行われている。流体に関しては解析対象領域内をメッシュと呼ばれる格子で区切り、その格子点や格子内部における物理量（速度、温度、圧力）を解析する格子法を用いた計算が行われている（非特許文献１参照）。また食塊を半固体としてあつかう場合は、有限要素法など、機械部品の構造解析手法を用いた計算が行われている（非特許文献２参照）。

神津ら，化学工学会秋季大会講演要旨集４１回，２００９，Ｐ０９ 水沼ら，日本機械学会年次大会講演論文集，２００５（２），８３−８４

しかし、従来の数値解析の主流である格子法では、実際の嚥下時の流体・食塊に見られる、表面の大変形や飛沫化などの現象を捉えることが難しく、実現象を再現することが困難であった。

本発明は、第１に、嚥下に関する実現象、すなわち、頭頸部器官及び飲食品、医薬品、医薬部外品を含む経口摂取品の挙動を近似的に再現することが容易な嚥下シミュレーション装置、動的三次元頭頸部モデル、嚥下シミュレーション方法、コンピュータ読み取り可能なプログラム、経口摂取品開発支援装置、経口摂取品開発方法、経口摂取品製造方法、食育支援装置、食育方法及び診断支援装置を提供することを目的とする。

本願の発明者らは、口腔器官の運動、飲食品等の物性を設定し、粒子法を用いて二次元的に口腔器官及び飲食品の挙動を解析する嚥下シミュレーション装置及び嚥下シミュレーション方法を提案した（特願２０１１−１４６７８０、特願２０１１−１４６７８１、共に未公開）。この装置及び方法では、従来と比較して、嚥下に関する実現象を近似的に再現することが可能となり、嚥下現象を可視化することが可能となった。しかしながら、二次元空間でのシミュレーションよりも、嚥下時における口腔器官及び飲食品の挙動をさらに三次元的に正確に表現したいという要請があった。また、使用プログラムの制約により経口摂取品の物性及び挙動に係る物理量を定量化し難いという問題があった。

本発明は、第２に、嚥下時における頭頸部器官及び飲食品の挙動を正確に表現でき、飲食品、医薬品、医薬部外品を含む経口摂取品の物性及び挙動に係る物理量を定量化できる嚥下シミュレーション装置、動的三次元頭頸部モデル、嚥下シミュレーション方法、コンピュータ読み取り可能なプログラム、経口摂取品開発支援装置、経口摂取品開発方法、経口摂取品製造方法、食育支援装置、食育方法及び診断支援装置を提供することを目的とする。

上記課題を解決するために、本発明の第１の態様に係る嚥下シミュレーション装置１００Ａは、例えば図２に示すように（具体的な器官については適宜図３を参照）、頭頸部器官からなる動的三次元頭頸部モデル１０ａを形成する頭頸部モデリング部１０と、動的三次元頭頸部モデル１０ａにおける各頭頸部器官の運動を設定する器官運動設定部３０と、解析対象としての飲食品、医薬品又は医薬部外品を含む経口摂取品及びその物性を設定する経口摂取品物性設定部４０と、飲食品、医薬品又は医薬部外品を含む経口摂取品をモデル化した擬似経口摂取品２０を口腔に入力する入力部８１と、動的三次元頭頸部モデル１０ａにおける各頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動を、粒子法を用いて三次元空間で解析する運動解析部５０と、運動解析部５０で解析された各頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動の解析結果を動画面に表示する表示部８２とを備える。

ここにおいて、頭頸部モデルの「動的三次元」は、三次元空間において時間の経過と共に変化することを意味するものである。なお、三次元空間はコンピュータ内に仮想三次元空間として形成されても良い。また、頭頸部器官は舌１１、喉頭１２、喉頭蓋１２ａ、気管１３、咽頭１４、口蓋１５（硬口蓋１５ａと軟口蓋１５ｂから成る）、顎１６、食道１８等で構成される（図３参照）。また、口腔１７、咽頭１４及び喉頭部１２ｃ（喉頭蓋１２ａと喉頭１２から成る）に係る各器官を含むものとする。動的三次元頭頸部モデル１０ａは上記各頭頸部器官を含んで構成される。食道１８及び気管１３は入り口部分だけ含めば良い。動的三次元頭頸部モデル１０ａは口腔器官の実際の運動に合わせて形成されるのが好ましいが、粒子数が増加すると三次元空間での解析を行なうコンピュータＰＣの負荷が増大するので、解析を簡素・容易にするため可動部分を限定しても良く、各頭頸部器官についてポリゴン（ポリゴン距離関数によって作成した壁境界で、流体などから外力を受けても移動しない壁を距離関数として空間に配位した境界）を定義することで頭頸部器官を粒子として計算しなくて済むようにしても良い。より詳細な検討を行なうためには、各頭頸部器官を粒子（剛体、弾性体、塑性体又は弾塑性体）として定義して計算することも可能である（越塚誠一，計算力学レクチャーシリーズ５，粒子法，日本計算工学会編，ｐｐ．５１−６８参照）。

動的三次元頭頸部モデル１０ａは仮想空間に形成するので、解析に支障がなければ、各頭頸部器官が三次元空間で重複し合う構造としても良い。各頭頸部器官の運動には、移動、回転、周期的運動等が含まれる。経口摂取品は液体、半固体（塑性を有するが、流動性を有さない）、固体のいずれでも良い。経口摂取品の物性には、密度、粘度、表面張力、接触角、熱容量、熱伝導率、動摩擦係数等が含まれる。ただし、これら全ての物性が設定される必要はなく、いずれかの物性が設定されれば良い。擬似経口摂取品２０の嚥下時の挙動とは、典型的には口腔から咽頭１４を経て食道１８に移動する挙動をいうが、咽頭１４や食道１８に到らず口腔１７に戻る場合や、喉頭１２や鼻腔に入ってしまう場合もあり得る。擬似経口摂取品２０の嚥下時の挙動に係る物理量には時間、位置座標、速度、圧力、温度、ずり速度、垂直応力、せん断応力等が含まれる。ただし、これら全ての物理量を用いて解析が行なわれる必要はなく、いずれかの物理量を用いて解析が行なわれれば良い。これらの物性及び物理量は定量化することが好ましいが、相対値を用いても良い。

入力部８１は、例えばマウスやキーボードを含み、動的三次元頭頸部モデル１０ａの口腔内部にマウスでカーソルをドラッグ・アンド・ドロップすることにより、擬似経口摂取品（擬似飲食品、擬似医薬品又は擬似医薬部外品を含む）２０を入力する。マウスのドラッグ・アンド・ドロップによって投入する疑似経口摂取品２０の位置が指定される。また、経口摂取品入力設定部４５（図２４参照）を設け、予め擬似経口摂取品とその投入場所と投入時間を設定しておき、自動的に擬似経口摂取品２０が投入されるようにしても良い。運動解析部５０は、粒子法を用いて解析するが、例えばＭＰＳ（Ｍｏｖｉｎｇ Ｐａｒｔｉｃｌｅ−Ｓｅｍｉ−ｉｍｐｌｉｃｉｔ）法を使用できる。表示部８２の「動画面に表示」は、典型的には液晶ディスプレイ等の動画面への表示が用いられる。また、ここでは、医師や歯科医師が健康診断や医療行為においてシミュレーションを行い解析結果を評価する場合を診断と称し、一般人が行う評価と区別することとする。すなわち、本明細書においては評価は診断を含まないものとする。動画面を表示することにより、評価者又は診断者が動画面を観察して評価又は診断するのに役立つが、自動的に評価又は診断を行う場合には、コンピュータ内に擬似画面表示部８２Ａ（図２４参照）を、記憶部８３に評価診断条件記憶部８３Ａ（図２４参照）を設け、解析結果を擬似画面表示部８２Ａの仮想動画面に動的に表示し、評価診断条件記憶部８３Ａの評価条件又は診断条件と照合してコンピュータ内の評価診断部６０（図２４参照）にて評価又は診断を行うが、このような擬似画面表示部８２Ａの仮想動画面に動的に表示する場合も「動画面に表示」に含まれるものとし、擬似画面表示部８２ＡはパーソナルコンピュータＰＣと表示部８２の両者に含まれるものとする。また、動的三次元頭頸部モデル１０ａの表示は、頭頸部器官毎に色彩を変えて半透明にして表示すると、頭頸部器官毎の動きを見易くできるので好適である。

本態様のように構成すると、動的三次元頭頸部モデル１０ａについて、頭頸部器官の運動、経口摂取品の物性を設定し、粒子法を用いて経口摂取品の挙動を解析するので、嚥下に関する実現象を近似的に再現することが容易な嚥下シミュレーション装置を提供できる。また、経口摂取品の物性及び経口摂取品の挙動に関する物理量を次元のある数値として定量的に設定して解析することにより、嚥下時における頭頸部器官や経口摂取品の挙動を正確に表現でき、嚥下時の挙動に係る物理量を時間毎に定量化できる嚥下シミュレーション装置を提供できる。

また、本発明の第２の態様に係る嚥下シミュレーション装置１００Ａは、第１の態様において、頭頸部モデリング部１０は頭頸部器官をポリゴン又は粒子（剛体、粉体、弾性体、塑性体、弾塑性体）に設定し、運動解析部５０は、擬似経口摂取品２０を粒子として取り扱う。
ここにおいて、動的三次元頭頸部モデル１０ａは、口腔、咽頭、喉頭、食道の実際の運動に合わせて形成され、時間毎に変形することが精度の高い解析を行なう上で好ましく、このためにはモデル全体を粒子として取扱うのが良いのであるが、粒子数が増加すると三次元空間での解析を行なうコンピュータＰＣの負荷が増大するので、解析を簡素・容易にするためには頭頸部器官をポリゴンに設定して解析するのが効率的である。
本態様のように構成すると、三次元解析を行うに際して、頭頸部器官をポリゴン、粒子（剛体、弾性体、塑性体、弾塑性体）に設定し、擬似経口摂取品２０を粒子（流体、剛体、弾性体、塑性体、弾塑性体）として取扱うので、頭頸部器官ならびに擬似経口摂取品２０の設定を任意に設定することで目的・状況に応じた多様な解析ができる

また、本発明の第３の態様に係る嚥下シミュレーション装置１００Ａは、第２の態様において、例えば図４に示すように（構成については図２参照）、頭頸部モデリング部１０は、舌１１を前奥方向にオトガイ舌筋１１ｆの起始である下顎骨のオトガイ部１６ａを要とするｎ個（ｎは２以上の整数）の扇形部１１ａ〜１１ｅに分割する構造とし、器官運動設定部３０は、各扇形部１１ａ〜１１ｅが半径方向に所定の位相差を持って連携して振動することにより擬似経口摂取品２０を奥方向に輸送するように進行波的波動運動を設定する。
ここにおいて、舌１１の分割数ｎは、３個以上でスムーズな進行波的波動運動を表現でき、７個を超えると計算の負荷が大きくなる。よって、４個〜６個が好適である。これによって、嚥下時の舌の進行波的波動動作を近似的に再現することが可能となった。
本態様のように構成すると、舌を前奥方向にｎ個の扇形部に分割し、上記の連携振動をさせることにより、嚥下時の舌の挙動を近似的に再現できるので、より実現象に近く効率的な解析ができる。なお、舌１１の形状は実際の運動に合わせて粒子で形成され、時間毎に変形させたものを計算毎に読み込むことも可能である。

また、本発明の第４の態様に係る嚥下シミュレーション装置１００Ａは、第３の態様において、例えば図６及び図７に示すように、頭頸部モデリング部１０は、各扇形部１１ａ〜１１ｅの半径方向の振動により各扇形部１１ａ〜１１ｅの間に隙間ができないように各扇形部１１ａ〜１１ｅの表面部分が相互に重複し合う構造とし、さらに、左右両側が半径方向に延びて中央部分に窪みを有する舌表面部１９を各扇形部１１ａ〜１１ｅの表面部分に重複させ、静止状態の舌表面部１１ｇの窪みの中で各扇形部１１ａ〜１１ｅが半径方向に振動する構造とし、器官運動設定部３０は、嚥下時に各扇形部１１ａ〜１１ｅが円周方向奥側に所定の角度回転し、嚥下最終時には舌１１の最奥の扇形部１１ｅで擬似経口摂取品２０を押し込むように設定する。
ここにおいて、所定の角度は１５°が好適であるが、これに限定されず、１０〜２０°でも良い。また、オーバーラップ部分は５ｍｍが好適であるが、これに限定されず、３〜７ｍｍでも良い。本態様のように構成すると、各扇形部１１ａ〜１１ｅの間に隙間ができず、経口摂取品２０の流路が窪みの中に形成され、経口摂取品２０が口腔の奥側にスムーズに移送されるので、嚥下時の舌１１の挙動を一層近似的に再現できる。

また、本発明の第５の態様に係る嚥下シミュレーション装置１００Ａは、第３の態様又は第４の態様において、例えば図９及び図１０に示すように、頭頸部モデリング部１０は、動的三次元頭頸部モデル１０ａを舌１１、口蓋１５、咽頭１４、喉頭部１２ｃ、気管１３と食道１８に分割し、喉頭部１２ｃを喉頭蓋１２ａと喉頭１２に分割し、嚥下開始前に食道部１８の食道入口１８ａを閉塞し、嚥下時には食道入口１８ａが拡張し、咽頭１４が収縮する構造とし、器官運動設定部３０は、嚥下時に、喉頭１２を下顎骨１６のオトガイ部１６ａ方向に移動して食道入口１８ａを開放し、咽頭１４が収縮し、次いで喉頭蓋１２ａを奥側に所定の角度回転させることにより、喉頭１２を閉塞するように設定する。
ここにおいて、所定の角度回転は１３５度が実現象に近く最も好適であるが、その近傍（例えば±１５度内）も好適である。また、「咽頭が収縮」とは、具体的には咽頭が上方向に短縮すると共にその内腔が収縮するのであるが、かかる内容を含むものである。このように構成すると、嚥下時の喉頭部１２ｃの挙動と食道１８の開閉を近似的に再現できる。なお、これらの器官の動きは器官運動設定部３０から数値を与えて強制的に動かすことのほか、これらの器官に付随した筋肉の動きと連動させて動かすことも可能である。さらに、時間毎に変化する各器官のモデルを解析毎に読み込み、よりスムーズな運動を与えることも可能である。

また、本発明の第６の態様に係る嚥下シミュレーション装置１００Ａは、第１の態様ないし第５の態様のいずれかにおいて、例えば図２２に示すように（構成については図２参照）、経口摂取品物性設定部４０は、解析対象として異なる物性の液体、半固体又は固体の複数の擬似経口摂取品２０を設定し、運動解析部５０は、複数の擬似経口摂取品２０について自由表面と複数の擬似経口摂取品２０間の境界を定め、複数の擬似経口摂取品２０の連動挙動の解析を行なう。
このように構成すると、複数の擬似経口摂取品２０について、経口摂取品の連動挙動を実現象に近いように再現することができ、連動挙動の解析に有効である。

また、本発明の第７の態様に係る嚥下シミュレーション装置１００Ｂは、第１の態様ないし第６の態様のいずれかにおいて、例えば図２４に示すように、動画面における擬似経口摂取品２０の嚥下時の挙動から、経口摂取品の食べ易さ及び／又は飲み易さを評価又は診断する評価診断部６０を備え、動画面は、運動解析部５０で解析された頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動の解析結果を模擬的に表示するために嚥下シミュレーション装置により仮想空間に形成された仮想動画面であり、評価診断部６０は、仮想動画面での擬似経口摂取品２０の挙動が予め定められた条件を満たすか否かで評価又は診断する。

ここにおいて、仮想動画面とは、パーソナルコンピュータＰＣの仮想空間に形成される仮想的な動画面をいうが、表示部８２に表示した場合の動画面と同じ内容の動画面となる。自動的に評価又は診断を行う場合には、コンピュータ内に擬似画面表示部８２Ａを、記憶部８３に評価診断条件記憶部８３Ａを設け、解析結果を擬似画面表示部８２Ａの仮想動画面に動的に表示し、評価診断条件記憶部８３Ａの評価条件又は診断条件と照合して評価又は診断を行う。なお、このような仮想動画面に動的に表示する場合も「動画面に表示」に含まれるものとする。また、予め定められた条件とは、例えば、喉頭１２に入らない、咽頭１４に詰まらない、舌１１上や咽頭１４に付着しない、口腔内導入から食道入口１８ａ通過までの時間が所定の範囲内である、口腔壁面におけるせん断応力が所定の値以下である、喉頭蓋１２ａ周辺でのせん断応力が所定の範囲内である等である。所定の範囲は、例えば、口腔１７内導入から食道入口１８ａ通過までの時間が０．８〜１．０ｓｅｃ、喉頭蓋１２ａ周辺での０．８ｓｅｃでのせん断応力が０．２５〜０．４５Ｎ／ｍ^２等である。
本態様のように構成すると、仮想動画面と評価診断条件記憶部８３Ａに記憶された評価条件又は診断条件とを照合して、経口摂取品の食べ易さ及び／又は飲み易さを自動的に評価又は診断できる。

また、本発明の第８の態様に係る嚥下シミュレーション装置１００Ｂは、第１の態様又は第７の態様において、例えば図２に示すように、擬似経口摂取品２０の嚥下時の挙動の解析結果から、経口摂取品の食べ易さ及び／又は飲み易さに係る評価診断結果を記録する評価診断結果記録部８３Ｂと、評価診断結果記録部８３Ｂに記録された評価結果又は診断結果に基づいて適正とされる経口摂取品の物性を決定する物性決定部７０とを備える。

ここにおいて、物性決定は嚥下シミュレーション装置１００Ｂが評価結果又は診断結果に基づいて自動的に決定する。人（たとえば評価者又は診断者）が決定する態様も可能であるが、この時は物性決定部７０が無くても良く、あっても使用されない又は決定結果が表示部８２に表示されて評価者又は診断者の参考に供されることとなる。

上記課題を解決するために、本発明の第９の態様に係る動的三次元頭頸部モデル１０ａは、例えば図４に示すように、頭頸部器官からなる動的三次元頭頸部モデル１０ａであって、舌１１を前奥方向にオトガイ舌筋１１ｆの起始である下顎骨のオトガイ部１６ａを要とするｎ個（ｎは２以上の整数）の扇形部１１ａ〜１１ｅに分割する構造とし、各扇形部１１ａ〜１１ｅが半径方向に所定の位相差を持って連携して振動することにより擬似経口摂取品２０を奥方向に輸送する進行波的波動運動を行なうように設定されている。
このように構成すると、嚥下時の舌１１及び咽喉の挙動を近似的に再現できる動的三次元頭頸部モデル１０ａを提供できる。

また、本発明の第１０の態様に係る動的三次元頭頸部モデル１０ａは、第９の態様において、例えば図６及び図７に示すように、各扇形部１１ａ〜１１ｅの振動により各扇形部１１ａ〜１１ｅの間に隙間ができないように各扇形部１１ａ〜１１ｅの表面部分が相互に重複し合う構造とし、さらに、左右両側が半径方向に延びて中央部分に窪みを有する舌表面部１１ｇを各扇形部１１ａ〜１１ｅの表面部分に重複させ、静止状態の舌表面部１９の窪みの中で各扇形部１１ａ〜１１ｅが円周方向に振動する構造とし、嚥下時に各扇形部１１ａ〜１１ｅが円周方向奥側に所定の角度回転し、次いで舌１１の最奥の扇形部１１ｅで擬似経口摂取品２０を押し込むように設定され、咽頭１４、喉頭部１２ｃと食道１８に分割され、喉頭部１２ｃは喉頭蓋１２ａと喉頭１２に分割し、嚥下開始前に食道１８の食道入口１８ａを閉塞し、嚥下時には食道入口１８ａが拡張し、咽頭１４が収縮する構造とし、嚥下時に、喉頭１２を下顎骨のオトガイ部１６ａ方向に移動して食道入口１８ａを開放し、咽頭１４が収縮し、次いで喉頭蓋１２ａを奥側に所定の角度回転させることにより、食道入口１８ａを閉塞するように設定されている。
このように構成すると、各扇形部１１ａ〜１１ｅの間に隙間ができず、擬似経口摂取品２０の流路が窪みの中に形成され、擬似経口摂取品２０が奥側にスムーズに移送されるので、嚥下時の舌１１、咽頭１４、喉頭１２と食道１８の挙動を近似的に再現できる動的三次元頭頸部モデルを提供できる。

上記課題を解決するために、本発明の第１１の態様に係る嚥下シミュレーション方法は、例えば図２３に示すように（装置構成については図２参照）、頭頸部器官からなる動的三次元頭頸部モデル１０ａを形成する頭頸部モデリング工程（Ｓ０１０）と、動的三次元頭頸部モデル１０ａにおける各頭頸部器官の運動を設定する器官運動設定工程（Ｓ０３０）と、解析対象としての経口摂取品、及びその物性を設定する経口摂取品物性設定工程（Ｓ０４０）と、経口摂取品をモデル化した擬似経口摂取品２０を口腔に入力する入力工程（Ｓ０５０）と、頭頸部モデル１０ａにおける各頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動を、粒子法を用いて解析する運動解析工程（Ｓ０６０）と、運動解析工程（Ｓ０６０）で三次元空間で解析された頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動の解析結果を動画面に表示する表示工程（Ｓ０７０）とを備える。

このように構成すると、動的三次元頭頸部モデル１０ａについて、頭頸部器官の運動、経口摂取品の物性を設定し、粒子法を用いて経口摂取品の挙動を解析するので、嚥下に関する実現象を近似的に再現することが容易な嚥下シミュレーション方法を提供できる。

また、本発明の第１２の態様に係る嚥下シミュレーション方法は、第１１の態様において、例えば図２５に示すように、擬似経口摂取品２０の嚥下時の挙動の解析結果から、経口摂取品の食べ易さ及び／又は飲み易さに係る評価を行う評価工程（Ｓ０８０）と、評価工程（Ｓ０８０）にて評価された評価結果に基づいて適正とされる経口摂取品の物性を決定する物性決定工程（Ｓ０９０）とを備える。
本態様のように構成すると、嚥下に関する実現象を近似的に且つ容易に再現するシミュレーションを通して、経口摂取品の食べ易さ及び／又は飲み易さが適正な経口摂取品の物性を高精度で効率的に導き出すことができる。

また、本発明の第１３の態様に係るプログラムは、第１１の態様又は第１２の態様に係る嚥下シミュレーション方法をコンピュータに実行させるためのコンピュータ読み取り可能なプログラムである。
ここにおいて、プログラムは、コンピュータ内蔵の記憶部に格納しても良く、インターネットからダウンロードしても良く、コンピュータ読み取り可能な記録媒体に格納しても良い。なお、本態様に係るコンピュータには、嚥下シミュレーション装置のように、コンピュータを含んで構成される装置のコンピュータ（例えばパーソナルコンピュータＰＣ）も含まれるものとする。

また、本発明の第１４の態様に係る経口摂取品開発支援装置２００Ａは、例えば図２６に示すように、第８の態様にかかる嚥下シミュレーション装置１００Ｂと、物性決定部７０で適正と決定された物性を有するように製造条件を適宜に設定して試作を行なった結果を記録する経口摂取品試作結果記録部８３Ｃと、経口摂取品試作結果記録部８３Ｃに記録された試作結果に基づき、経口摂取品の物性を物性決定部７０で適正と決定された物性にする製造条件を決定する製造条件決定部８４とを備える。
このように構成すると、経口摂取品開発支援装置２００Ａは、動的三次元頭頸部モデル１０ａについて、頭頸部器官の運動、経口摂取品の物性を設定し、粒子法を用いて経口摂取品の挙動を解析するので、嚥下に関する実現象を再現することが容易である。そして、かかる嚥下に関する実現象を再現することが容易な嚥下シミュレーション装置を利用して経口摂取品を評価又は診断するので、食べ易さ、飲み易さが優れた経口摂取品を高い確実性をもって開発できる。

また、本発明の第１５の態様に係る経口摂取品開発支援装置２００Ａは、第１４の態様において、例えば図２６に示すように、動画面における擬似経口摂取品２０の嚥下時の挙動から、経口摂取品の食べ易さ及び／又は飲み易さを評価又は診断する評価診断部６０を備え、動画面は、運動解析部５０で解析された各頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動の解析結果を模擬的に表示するために嚥下シミュレーション装置１００Ｂにより仮想空間に形成された仮想動画面であり、評価診断部６０は、仮想動画面での擬似経口摂取品の挙動が予め定められた条件を満たすか否かで評価又は診断する。

本態様のように構成すると、経口摂取品開発支援装置２００Ａは、経口摂取品の挙動を仮想動画面と評価診断条件記憶部８３Ａに記憶された評価条件又は診断条件とを照合して、自動的に評価又は診断し、評価結果又は診断結果を経口摂取品開発に結び付けるので、効率的に、食べ易い又は飲み易い経口摂取品を開発することができる。

また、本発明の第１６の態様に係る経口摂取品開発方法は、例えば図２７に示すように、第１２の態様にかかる嚥下シミュレーション方法を備え、経口摂取品物性設定工程（Ｓ０４０）において、経口摂取品の物性を変更して設定し、それ以降の物性決定工程（Ｓ０９０）までを繰り返し行ない、又は、経口摂取品物性設定工程（Ｓ０４０）において、経口摂取品の物性を変更して設定し、それ以降の評価工程（Ｓ０８０）までを繰り返し行ない、その後に一括して物性決定工程（Ｓ０９０）を行ない、さらに、物性決定工程（Ｓ０９０）で適正と決定された物性を有するように製造条件を適宜に設定して試作を行なう経口摂取品試作工程（Ｓ１１０）と、経口摂取品試作工程（Ｓ１１０）の結果に基づき、経口摂取品の物性を物性決定工程（Ｓ０９０）で適正と決定された物性にする製造条件を決定する製造条件決定工程（Ｓ１２０）とを備える。

このように構成すると、動的三次元頭頸部モデル１０ａについて、頭頸部器官の運動、飲経口摂取品の物性を設定し、粒子法を用いて経口摂取品の挙動を解析するので、嚥下に関する実現象を再現することが容易である。そして、かかる嚥下に関する実現象を再現することが容易な嚥下シミュレーション方法を利用して経口摂取品を評価するので、食べ易さ、飲み易さが優れた飲経口摂取品を確実に開発できる。

また、本発明の第１７の態様に係る経口摂取品開発方法は、第１６の態様において、例えば図２７に示すように、動画面は、運動解析工程（Ｓ０６０）で解析された各頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動の解析結果を模擬的に表示するために嚥下シミュレーション装置１００Ｂにより仮想空間に形成された仮想動画面であり、表示工程（Ｓ０７０）は、解析結果を仮想動画面に模擬的に表示する仮想表示工程であり、評価工程（Ｓ０６０）は、仮想表示工程で模擬的に表示された仮想動画面での擬似経口摂取品２０の挙動が予め定められた条件を満たすか否かで評価する。
本態様のように構成すると、経口摂取品開発支援装置２００Ａは、仮想動画面と評価診断条件記憶部８３Ａに記憶された評価条件又は診断条件とを照合して、経口摂取品の食べ易さ、飲み易さを自動的に評価又は診断し、評価結果又は診断結果を経口摂取品の開発に結び付けるので、効率的に、食べ易い又は飲み易い経口摂取品を開発することができる。

また、本発明の第１８の態様に係る経口摂取品製造方法は、例えば図２８に示すように、第１６の態様又は第１７の態様の経口摂取品開発方法の製造条件決定工程（Ｓ１２０）にて決定された製造条件を用いて経口摂取品を製造する。

経口摂取品開発方法では、典型的にはシミュレーションを繰り返し行なって、適正な物性を確認し、そして適正な物性にするように製造条件を決定することにより、経口摂取品を開発していくのであるが、経口摂取品製造方法においては、製造工程中のいずれかの工程（例えば原料配合工程や焼成工程）で、経口摂取品開発方法の製造条件決定工程（Ｓ１２０）において決定された製造条件を用いて製造すれば、当該態様に該当することになる。
このように構成すると、食べ易さ、飲み易さが優れた経口摂取品を確実に開発できる。

上記課題を解決するために、本発明の第１９の態様に係る食育支援装置３００Ａは、例えば図２９に示すように、第１の態様又は第７の態様の嚥下シミュレーション装置１００Ａ（又は１００Ｂ）と、擬似経口摂取品２０の嚥下時の挙動の解析結果から、経口摂取品の食べ易さ、飲み易さの評価結果又は診断結果を記録する評価診断結果記録部８３Ｂと、表示部８２で動画面に表示された擬似経口摂取品２０の嚥下時の挙動を、評価診断結果記録部８３Ｂに記録された経口摂取品の評価結果又は診断結果と関連させて解説する教示部８５とを備える。

本態様は第１の態様に係る嚥下シミュレーション装置１００Ａを応用した食育支援装置３００Ａである。ここにおいて、教示は教師が動画面又は／及び教材を用いて行っても良く、食育支援装置に教示すべきことが動画像と連動して音声出力可能に組み込まれ、食育支援装置３００Ａが音声出力を通して自動的に行っても良い。また、教示内容は評価結果又は診断結果に基づいて教示部８５で自動的に作成しても良いが、教育者が更にわかり易く編集・補充するのが好適である。なお、第１の態様に係る嚥下シミュレーション装置１００Ａに代えて第７の態様に係る嚥下シミュレーション装置１００Ｂを応用しても良い。
このように構成すると、嚥下に関する実現象を再現することが容易な嚥下シミュレーション装置を用いて嚥下現象を表示するので、飲経口摂取品の食べ易さ又は飲み易さを理解しやすく、食育に有効である。

上記課題を解決するために、本発明の第２０の態様に係る食育方法は、例えば図３０に示すように、第１１に態様に係る嚥下シミュレーション方法と、擬似経口摂取品２０の嚥下時の挙動の解析結果から、経口摂取品の食べ易さ及び／又は飲み易さに係る評価を行う評価工程（Ｓ０８０）と、表示工程（Ｓ０７０）で動画面に表示された擬似経口摂取品２０の嚥下時の挙動を、評価工程（Ｓ０８０）で評価された経口摂取品の評価結果と関連させて解説する教示工程（Ｓ０８２）とを備える。

本態様は第１９の態様に係る食育支援装置３００Ａに対応する食育方法である。このように構成すると、嚥下に関する実現象を再現することが容易な嚥下シミュレーション方法を用いて嚥下現象を表示するので、教示を受ける者が飲経口摂取品の食べ易さ又は飲み易さを理解し易く、食育に有効である。

本発明の第２１の態様に係る嚥下シミュレーション装置４００Ａは、例えば図３１に示すように、第１の態様ないし第８に態様に係る嚥下シミュレーション装置１００Ａ，１００Ｂにおいて、器官運動設定部３０で設定された器官運動パラメータについて、運動解析部５０で解析された解析結果から、被診断者の症状に合う器官運動パラメータを決定する器官運動決定部７５を備える。
このように構成すると、各頭頸部器官の機能低下した部分を発見でき、診断に役立てられる。

本発明の第２２の態様に係る診断支援装置は、第２１の態様に係る嚥下シミュレーション装置を備える。
このように構成すると、嚥下シミュレーション装置を利用して、嚥下現象を考慮した迅速な診断が可能になる。

本発明によれば、第１に、嚥下に関する実現象、すなわち、頭頸部器官及び飲経口摂取品の挙動を近似的に再現することが容易な嚥下シミュレーション装置、動的三次元頭頸部モデル、嚥下シミュレーション方法、コンピュータ読み取り可能なプログラム、経口摂取品開発支援装置、経口摂取品開発方法、経口摂取品製造方法、食育支援装置、食育方法及び診断支援装置を提供できる。
第２に、嚥下時における頭頸部器官及び経口摂取品の挙動を正確に表現でき、経口摂取品の物性及び挙動に係る物理量を定量化できる嚥下シミュレーション装置、動的三次元頭頸部モデル、嚥下シミュレーション方法、コンピュータ読み取り可能なプログラム、経口摂取品開発支援装置、経口摂取品開発方法、経口摂取品製造方法、食育支援装置、食育方法及び診断支援装置を提供できる。

格子法と粒子法の違いを説明するための図である。 実施例１における嚥下シミュレーション装置の構成例を示す図である。 実施例１における動的三次元頭頸部モデルの構成例を示す図である。 舌モデルの構成例を示す図である。 舌モデルの別の構成例（初期の構成例）を示す図である。 舌構造のオーバーラップを説明するための図である。 中央に窪んだ流路を持つ舌構造を説明するための図である。 実施例１における咽頭モデルの構成例を示す図である。 動的三次元頭頸部モデルにおける喉頭の移動と食道入口の開放を示す図である。 喉頭部モデルにおける喉頭蓋の回転運動を示す図である。 喉頭と喉頭蓋の連動運動例を示す図である。 水の場合の三次元シミュレーションの解析結果を示す図である。 水、牛乳の場合の三次元シミュレーションの解析結果を示す図である。 せん断応力と流下速度についての解析領域を示す図である。 領域Ａにおけるせん断応力の経時変化を示す図である。 領域Ｂにおけるせん断応力の経時変化を示す図である。 領域Ｂにおける流下速度の経時変化を示す図である。 残存率と誤嚥率についての解析領域を示す図である。 領域Ｃ及び領域Ｄ（合計）における残存率を示す図である。 領域Ｅにおける誤嚥率を示す図である。 非ニュートン流体の場合の三次元シミュレーションの解析結果を示す図である。 固体と流体の連成解析の例を示す図である。 実施例１における嚥下シミュレーション方法の処理フロー例を示す図である。 実施例２における嚥下シミュレーション装置の構成例を示す図である。 実施例２における嚥下シミュレーション方法の処理フロー例を示す図である。 実施例５における経口摂取品開発支援装置の構成例を示す図である。 実施例５における経口摂取品開発方法の処理フロー例を示す図である。 実施例５における経口摂取品製造方法の処理フロー例を示す図である。 実施例６における食育支援装置の構成例を示す図である。 実施例６における食育方法の処理フロー例を示す図である。 実施例７における嚥下シミュレーション装置の構成例を示す図である。 実施例７における診断支援の処理フロー例を示す図である。

以下に図面に基づき本発明の実施の形態について説明する。なお、各図において、互いに同一又は相当する部分には同一符号を付し、重複した説明は省略する。

［粒子法］
本実施の形態では、液面の大変形や飛沫などの表現が可能な解析方法として、解析対象の液体や固体を粒子として扱う粒子法を用いてシミュレーションを行なう。そこで、まず、粒子法について説明する。

図１に従来の解析手法である格子法と新しい解析手法である粒子法の違いを示す。図１（ａ）に格子法の概念図を、図１（ｂ）に粒子法の概念図を示す。格子法では解析領域を格子で分割し、格子毎に物理量を計算する。すなわち、液面の変化は格子の形に沿うものとなり、飛沫や液面が大きく変形する場合の解析は困難である。これに対して、粒子法、特にＭＰＳ（Ｍｏｖｉｎｇ Ｐａｒｔｉｃｌｅ−Ｓｅｍｉ−ｉｍｐｌｉｃｉｔ）法は１９９５年に開発された比較的新しい解析手法である（Ｋｏｓｈｉｚｕｋａ ｅｔ ａｌ，Ｃｏｍｐｕｔ．Ｆｌｕｉｄ Ｄｙｎａｍｉｃｓ Ｊ，４，２９−４６，１９９５）。流体を粒子で置き換え、粒子毎に物理量を計算する。その結果、液面の微妙な変化の解析が可能となり、飛沫や液面が大きく変形する場合の解析が可能になる。しかしながら、粒子法を用いた生体内部の流体・食塊の解析はこれまで見出されていない。そこで、発明者達は、生体内流体・食塊挙動の推定に粒子法を応用したシミュレーション装置とシミュレーション方法とを開発した。

発明者達が提案した粒子法による二次元嚥下シミュレータ（特願２０１１−１４６７８０、特願２０１１−１４６７８１）では、解析手法として粒子法を採用することによって従来の技術では困難であった食塊表面の大変形や飛沫化を再現できるようになった。
しかしながら、二次元空間でのシミュレーションよりも、嚥下時における口腔器官、咽喉器官などの頭頸部器官及び飲食品等の挙動をさらに三次元的に正確に表現したいという要請があった。また、使用プログラムの制約により経口摂取品の物性及び挙動に係る物理量を定量化し難いという問題があった。

上記問題を解決するために、発明者達は三次元嚥下シミュレータを開発した。二次元嚥下シミュレータでは単一平面上に単純化した構造での解析ができた。三次元シミュレータでは、さらに、医学的知見に基づいて口腔、咽頭、喉頭部の正確な構造や挙動を再現した。また、三次元嚥下シミュレータは、二次元嚥下シミュレータでは無次元の相対的値として設定していた経口摂取品（食塊）の物性（密度、粘度、表面張力、接触角、熱容量、熱伝導率、動摩擦係数など）を次元のある数値として設定可能とするように、また、嚥下中の食塊の物理量（時間、位置座標、速度、圧力、温度、ずり速度、垂直応力、せん断応力など）についても、次元のある数値として抽出するように改善し、定量的な評価又は診断をすることが可能となった。

［装置構成］
図２に実施例１における嚥下シミュレーション装置１００Ａの構成例を示す。実施例１では、評価者又は診断者が経口摂取品の入力を行い、動画像を見て嚥下評価又は嚥下診断を行なう例を説明する。なお、動的三次元頭頸部モデル１０ａの構成例については図３を参照されたい。

嚥下シミュレーション装置１００Ａは、頭頸部器官からなる動的三次元頭頸部モデル１０ａを形成する頭頸部モデリング部１０と、動的三次元頭頸部モデル１０ａ（以下単に頭頸部モデルともいう）における各頭頸部器官の運動を設定する器官運動設定部３０と、解析対象としての経口摂取品及びその物性を設定する経口摂取品物性設定部４０と、経口摂取品をモデル化した擬似経口摂取品（擬似飲食品、擬似医薬品又は擬似医薬部外品を含む）２０を口腔に入力する入力部８１と、頭頸部モデル１０ａにおける各頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動を、粒子法を用いて三次元空間で解析する運動解析部５０と、運動解析部５０で解析された各頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動の解析結果を動画面に表示する表示部８２と、評価者又は診断者が動画面を観察して経口摂取品について評価又は診断した評価結果又は診断結果に基づいて適正とされる経口摂取品の物性を決定する物性決定部７０と、嚥下シミュレーション装置１００Ａ及びその各部を制御して、嚥下シミュレーション装置１００Ａの諸機能を実行させる制御部９０と、動的三次元頭頸部モデル１０ａ、設定条件（器官運動、経口摂取品の物性等）、解析結果、評価結果又は診断結果を記憶する記憶部８３とで構成される。また、表示部８２は頭頸部モデル１０ａ、設定条件（器官運動、経口摂取品の物性等）、評価結果又は診断結果等も画面に表示する。このうち、頭頸部モデリング部１０、器官運動設定部３０、経口摂取品物性設定部４０、運動解析部５０、物性決定部７０、制御部９０はパーソナルコンピュータＰＣで実現でき、パーソナルコンピュータＰＣ内に設けられる。評価者又は診断者は表示部８２の動画面を観察して評価又は診断を行ない、評価結果又は診断結果を入力部８１から入力する。入力された評価結果又は診断結果は記憶部８３の評価診断結果記録部８３Ｂに記録される。ここでは、医師や歯科医師が健康診断や医療行為においてシミュレーションを行い解析結果を評価する場合を診断と称し、一般人が行う評価と区別することとする。すなわち、本明細書においては評価は診断を含まないものとする。

嚥下シミュレーションの対象として、飲食品の他に医薬品、医薬部外品を含む経口摂取品を扱うことが可能である。また、場合によっては、経口摂取しない物質（例えばビー玉など）も誤飲、窒息の検証のために嚥下シミュレーションの対象として取扱うことも可能である。経口摂取品物性設定部４０は経口摂取品が液体であれば、液量・粘度・表面張力・比重・熱伝道率・比熱等の物性を、経口摂取品が固体であれば、形状・寸法・弾性係数・引っ張り強さ・降伏点・降伏点応力・粘度のずり速度依存性・動的粘弾性・静的粘弾性・圧縮応力・破断応力・破断ひずみ・硬度・付着性・凝集性・熱伝導率・比熱等の物性を、経口摂取品が半固体（可塑性があるが、流動性はない）であれば、量・粘度・比重・降伏点・降伏点応力・粘度のずり速度依存性・動的粘弾性・静的粘弾性・圧縮応力・付着性・凝集性等の物性を設定する。

動的三次元頭頸部モデル１０ａは、口腔器官の実際の運動に合わせて形成され、時間毎に変形することが精度の高い解析を行なう上で好ましく、このためにはモデル全体を粒子（剛体、粉体、弾性体、塑性体、弾塑性体）として取扱うのが良いのであるが、粒子数が増加すると三次元空間での解析を行なうコンピュータＰＣの負荷が増大するので、解析を簡素・容易にするためには頭頸部器官をポリゴンとして設定して解析するのが効率的である。

入力部８１は、マウス、キーボード等の入力機器で構成され、投入すべき擬似経口摂取品２０を口腔内に投入する。例えばマウスでカーソルを口腔内にドラッグ・アンド・ドロップし、擬似経口摂取品２０の口腔内への投入位置は例えば口腔内で歯の近傍（例えば擬似経口摂取品の長さの２倍以内）とし、ドラッグ・アンド・ドロップ直後の時間を投入時間とする。

運動解析部５０では、頭頸部器官の運動に伴う擬似経口摂取品２０（以下、単に経口摂取品ともいう）の嚥下時の挙動を解析する。舌１１の進行波的波動運動、喉頭蓋１２ａの回転運動、喉頭１２の往復運動等により、頭頸部内部に投入された経口摂取品が動かされる。経口摂取品の動きは粒子法により解析される。経口摂取品は固体・半固体・液体のいずれでもポリゴン又は粒子として取り扱われる（計算負荷軽減のためポリゴンとして取扱うこともある）。

表示部８２は経口摂取品の挙動の解析結果を動画面に表示する。動画像の１コマを静止画像として表示することも可能であり、時間を逆に辿り巻き戻し表示することも可能である。
記憶部８３は、頭頸部モデル１０ａ、器官特性、設定条件、解析結果、評価結果又は診断結果等を記憶する。

評価又は診断は評価者又は診断者が表示部８２の動画面を見て行なう。例えば表示部８２に表示された評価表又は診断表のセルに○×、ランク、点数等を入力する。評価結果又は診断結果は評価診断結果記録部８３Ｂに記録される。経口摂取品物性設定部４０で経口摂取品の物性値を変化させて評価又は診断を行なうことにより、経口摂取品の適正な物性値を求めることができる。物性決定部７０は、評価診断結果記録部８３Ｂに記録された評価結果又は診断結果に基づいて適正とされる経口摂取品の物性を自動的に決定する。なお、人が決定する場合は実施例４に示す。物性数は単数でも複数でも良い。適正な物性は例えば、適正な範囲のマップを作成して示しても良く、適正な範囲を複数のレベルに分けて（ランクＡ〜ランクＣ等）示しても良く、複数の点を示しても良く、最適な１点を示しても良い。求めるべき物性が多種類の時には、多次元の主成分分析を用いて、適正な物性の範囲を求めても良い。

制御部９０は、嚥下シミュレーション装置１００Ａ及びその各部を制御して、嚥下シミュレーション装置１００Ａの諸機能を実行させる。制御部９０は内蔵メモリに嚥下シミュレーター（解析用ソフトウェア）を保有する。

［頭頸部モデル］
図３に本実施例に係る頭頸部モデル１０ａの構成例を示す。図３（ａ）に頭頸部モデル１０ａの斜視図を、図３（ｂ）に正面図を、図３（ｃ）に側面図を、図３（ｄ）にＡ−Ａ断面図を示す。図において、１１は舌（オトガイ舌筋１１ｆ（図４参照）を含む）、１２は喉頭、１２ａは喉頭蓋、１３は気管、１４は咽頭（１４ａは咽頭の管壁、１４ｂは咽頭の粘膜）、１５は口蓋（１５ａは硬口蓋、１５ｂは軟口蓋）、１６は顎（オトガイ部１６ａを含む）、１７は口腔、１８は食道（１８ａは食道入口、１８ｂは食道の管壁）、２０は疑似経口摂取品としての球状食塊である。
粒子法では粒子数が増加するとコンピュータの負担が増加し、解析に時間がかかる。本実施例では、経口摂取品を粒子（流体）として表現し、各口腔器官を、ポリゴン距離関数によって作成した壁境界（流体などから外力を受けても移動しない壁を距離関数として空間に配位した境界）として表現することにより、頭頸部モデル１０ａを簡素化し、コンピュータの負担を軽減した。各頭頸部器官の運動（進行波的波動運動、回転運動、往復運動等）は器官運動設定部３０で設定される。

頭頸部モデル１０ａの作成に際して、解剖で理解されている構造の知見及び、ＣＴ（コンピュータ断層撮影：Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）画像より大まかに読み取ることの出来る軟口蓋１５ｂや舌１１の形態、及び気管１３の形態から、咽頭１４と食道入口１８ａの位置を推定した。舌１１、軟口蓋１５ｂ、咽頭１４と食道入口１８ａ、喉頭蓋１２ａと喉頭１２の構造を、ＣＧ（コンピュータグラフィックス：Ｃｏｍｐｕｔｅｒ Ｇｒａｐｈｉｃｓ）用ソフトウェア（Ａｕｔｏｄｅｓｋ ３ｄｓ Ｍａｘ等）を用いてモデリングした。次に、得られた動的三次元頭頸部モデル１０ａに対して、ＶＦ（嚥下造影検査：Ｖｉｄｅｏｆｌｕｏｒｏｓｃｏｐｉｃ ｅｘａｍｉｎａｔｉｏｎ ｏｆ ｓｗａｌｌｏｗｉｎｇ）による嚥下時の造影画像（正面及び側面図）を重ね合わせて、立体構造を修正させた。さらに、空間領域の三次元的な変化の概略を、嚥下のシネＭＲＩ（シネ磁気共鳴画像：Ｃｉｎｅ Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）の動画を参照して、動きを付与した。なお、シネＭＲＩは、心拍同期による心臓撮影法の原理（同期サンプリング法）の発展形である。まず、ＭＲＩ装置へのトリガーの一定間隔での入力とトリガーに同期した嚥下運動の繰り返しによって、複数の断面のＭＲＩ動画を得る。これを立体構築し、さらに時間軸上に配置することで４次元の再構築画像が得られる。

図３の頭頸部モデル１０ａを用いて嚥下運動をシミュレーションした。シミュレーションに際して頭頸部モデル１０ａを距離関数を用いて表現した。粒子法では空間点（三次元座標）にある経口摂取品（粒子）に対して、最短距離にある器官が最も強く作用し、距離が遠くなるにつれて作用が弱くなる。距離関数とは空間点（三次元座標）に対して、頭頸部モデル１０ａまでの最短の距離で定義される関数であり、距離関数を用いると頭頸部モデル１０ａの重ね合わせを容易に行うことができる（全ての頭頸部器官からの距離の最小値が空間点に定義される）。頭頸部モデル１０ａの重ね合わせが容易であるため、分割した各頭頸部器官に対して個別に移動量を設定し、後に重ね合わせることによって、嚥下運動を再現することができる。

［舌モデル］
図４に舌モデル１０ｂの構成例を示す。図４（ａ）は舌１１の筋肉（下顎骨のオトガイ部１６ａから延びるオトガイ舌筋１１ｆ）の配置を示し、図４（ｂ）はオトガイ舌筋１１ｆに沿って扇形に分割された舌モデル１０ｂを示す。舌モデル１０ｂは頭頸部モデリング部１０で形成し、舌の運動は器官運動設定部３０で設定する。舌モデル１０ｂは、舌１１による食塊の輸送を再現する必要がある。舌モデル１０ｂは、オトガイ舌筋の起始である下顎骨のオトガイ部１６ａを扇の要とし、起始から延びるオトガイ舌筋１１ｆに沿って舌１１を扇形に分割した構造とした。分割は例えば口唇を前、咽喉を後としたときの前後方向にｎ分割（本実施例ではｎ＝５）とし、横方向には分割しなかった。分割された各部を扇形部１１ａ〜１１ｅと称することとする。舌１１の運動による経口摂取品の輸送については、各扇形部１１ａ〜１１ｅを半径方向に振動させ、前方から後方に向かって振動の位相をずらしていくことで達成される。振動は、例えば、オトガイ舌筋１１ｆを中心に振幅５ｍｍで振動数１．１Ｈｚで振動させ、各扇形部１１ａ〜１１ｅは手前側から奥側に向かって順次０．１ｓずつ振動のタイミング（位相）を遅らせるように設定した。

図５に舌モデルの別の構成例（初期の構成例）を示す。図５（ａ）に分割前の舌モデル、図５（ｂ）に分割後の舌モデルを示す。図５（ｂ）に示すように、当初は舌１１を６×６＝３６分割するモデルを検討したが、予想していたよりも滑らかな舌１１の運動を再現できず、さらに３６分割としたことによってＰＣ（Ｐｅｒｓｏｎａｌ Ｃｏｍｐｕｔｅｒ）への負荷が増大し、解析が正常に行なわれなかった。このため、分割を最小にすることを考え、図４のモデルに想到した。ただし、ＰＣの処理能力が向上した場合にはかかる３６分割などの多分割モデルも使用可能性がある。
また、今回の実施例では計算負荷を低減するために舌１１を分割したモデルを使用しているが、分割せずに、形状が変化した舌モデルを各計算時間毎に読み込み、より滑らかな進行波的波動運動を模擬することも可能である。

図６は舌構造のオーバーラップを説明するための図である。図６（ａ）にオーバーラップ前の舌構造、図６（ｂ）にオーバーラップ後の舌構造、図６（ｃ）にオーバーラップの細部を示す。オーバーラップ前に単純に舌１１の各扇形部１１ａ〜１１ｅを振動運動させると、図６（ａ）に示すように各扇形部１１ａ〜１１ｅの隙間で切断面が外に現れてしまうため、図６（ｂ）及び図６（ｃ）に示すように、各扇形部１１ａ〜１１ｅの切断面に例えば高さ５ｍｍのドーム形状を作成し、２つの構造をオーバーラップさせることで切れ目が出現しないようにした。

図７は中央に窪んだ流路を持つ舌構造を説明するための図である。図７（ａ）は平常時の舌構造に中央が窪んだ舌構造（舌表面部）１１ｇをオーバーラップさせたモデルを示す斜視図、図７（ｂ）は中央が窪んだ舌構造（舌表面部）１１ｇに対する平常時の舌構造（扇形部）１１ａ〜１１ｅの動きを説明する図である。図７（ｂ）では舌１１として示す部分は各扇型部１１ａ〜１１ｅが前後方向に重なり合っている。平常時の舌構造１１ａ〜１１ｅは、図６（ｂ）、（ｃ）の舌構造であり、中央が窪んだ舌構造（舌表面部）１１ｇは、左右両側が半径方向に延びて中央部分に窪みを有する構造である。かかる中央が窪んだ舌構造１１ｇを平常時の舌構造１１ａ〜１１ｅと重複させる。そうすると、中央が窪んだ舌構造１１ｇの左右両側が半径方向に延びて壁面を構成し、平常時の舌構造（各扇形部）１１ａ〜１１ｅの表面は、中央が窪んだ舌構造１１ｇの窪みの中で半径方向に振動する（図７（ｂ）では上下に振動する）。舌１１の輸送は舌１１の中心に食塊２０を集めて流し込むという医学的な知見から、嚥下直前の中央が窪んでいる舌構造１１ｇを加えることで、食塊（擬似経口摂取品）２０の流路を確保している。

さらに、各扇形部１１ａ〜１１ｅの半径方向の振幅を例えば１０〜１５ｍｍとした。また、嚥下時に振動運動に加えて、各扇形部１１ａ〜１１ｅに円周方向奥側に例えば１５°の回転運動（振動しながら扇の要を中心にトータル１５°回転する）を付与した。位相をずらした扇形部１１ａ〜１１ｅの振動と回転運動により、食塊の輸送が行なわれる。なお、上記の振幅、振動数、振動のタイミング、オーバーラップ高さ、回転角度は１例であり、好適な範囲で（例えば±２０％）自由に設定できる。

さらに、舌１１の最奥の扇形部１１ｅでは、擬似経口摂取品２０を押し込むようにした。扇形部１１ｅの具体的な押し込みの動きは、例えば次のようである。嚥下のタイミングまではその他の扇形部１１ａ〜１１ｄと同様に振幅１０ｍｍで振動しながら、同時に回転していく。ただし、他の扇形部１１ａ〜１１ｄは０．３秒掛けて１５度回転する運動をするのに対して、最奥の扇形部１１ｅは０．１秒で回転を停止し（５度だけ回転）、その後、０．２秒掛けて１５度逆回転する。この逆回転中は扇形の要を中心とした振動運動は、逆回転開始時を最大振幅として、徐々に回転中心側に落ち込んでいく挙動を示す。その後、扇形部１１ｅは０．２秒で１０度回転し、ほぼ嚥下前の状態に戻り、次の嚥下運動を待つ。この戻りの時に擬似経口摂取品２０が押し込まれる。
この舌モデル１０ｂを用いて、重力方向を変えて解析を行い、舌１１の運動による食塊（ここでは水物性）２０の輸送が達成されることを確認できた。

［喉頭部モデル］
食道１８の入口は平常時には閉じており、嚥下時に食道１８が開放する様子を再現しなくてはならない。喉頭部モデル１０ｃは、喉頭蓋１２ａと喉頭１２の動きにより、食道１８と喉頭１２の開放と閉塞を実現する。喉頭部モデル１０ｃは頭頸部モデリング部１０で形成し、その運動は器官運動設定部３０で設定する。

図８に喉頭部モデル１０ｃの構成例を示す。喉頭部モデル１０ｃは喉頭蓋１２ａと喉頭１２で構成される。喉頭蓋１２ａには喉頭１２に連なる穴があり、喉頭１２は管状で、喉頭蓋１２ａは喉頭１２の入り口（喉頭口）の管の周りの堤防状の隆起である。喉頭１２が喉頭蓋１２ａを搭載した状態で下顎骨オトガイ部１６ａ方向とその反対側に移動する往復運動と喉頭蓋１２ａが喉頭１２上で回転する回転運動を行なう。ここで、喉頭１２の移動と喉頭蓋１２ａの回転は同時に進行する。この間に咽頭１４は上方向に短縮すると共に内腔は収縮する。

図９に動的三次元頭頸部モデル１０ａにおける喉頭１２の移動と食道入口１８ａの開放を示す。図９（ａ）は喉頭１２移動前（平常時）の食道入口１８ａが閉塞された状態を、図９（ｂ）は喉頭１２移動後の食道入口１８ａが開放された状態を示す。また、図９（ｃ）は平常時（食道入口１８ａが閉塞状態の時）の喉頭部モデル１０ｃの断面を示す図である。図９（ａ），（ｂ）に示すように、平常時は、喉頭１２は下顎骨のオトガイ部１６ａの斜下方向にあり、嚥下時には斜上方向である下顎骨のオトガイ部１６ａの方向に移動して食道入口１８ａを開放している。図９（ａ）の矢印は喉頭部１２ｃの移動方向を示す。また、図９（ｃ）に示すように、平常時に食道入口１８ａは上食道括約筋によって閉鎖し、食道１８の前壁は気管１３の後壁に接している。

図１０に喉頭部モデル１０ｃにおける喉頭蓋１２ａの回転運動を示す。図１０（ａ）は回転前、図１０（ｂ）は回転後の状態を示す。図１０（ｂ）には喉頭１２の移動方向と喉頭蓋１２ａの回転方向を矢印で示す。回転の角度は現実的な動きと照らし合わせて例えば１３５度に設定した。嚥下時には喉頭１２が下顎骨のオトガイ部１６ａ側に移動し、食道入口１８ａが開放されると共に、喉頭蓋１２ａの回転運動により喉頭１２の入口が塞がれる。

図１１に喉頭部モデル１０ｃ（喉頭部は１２ｃ）における喉頭１２と喉頭蓋１２ａの連動運動例を示す。図１１（ａ）は０．０ｓｅｃ、（ｂ）は０．２ｓｅｃ、（ｃ）は０．４ｓｅｃ、（ｄ）は０．５ｓｅｃ、（ｅ）は０．６ｓｅｃ、（ｆ）は０．７ｓｅｃ、（ｇ）は０．８ｓｅｃ、（ｈ）は０．９ｓｅｃ、（ｉ）は１．０ｓｅｃにおける状態を示す。喉頭蓋１２ａは０．２ｓｅｃで傾き始め、０．６ｓｅｃで横倒しになり、０．７ｓｅｃで最大１３５度になり、その後逆回転し、１．０ｓｅｃで元に戻る。喉頭１２は０．６ｓｅｃで移動を開始し、０．７ｓｅｃで最大になり、０．９ｓｅｃで元に戻る。図中の矢印は舌（扇形部）１１の最も上に出ている部分を示す。なお、前述の、喉頭蓋１２ａの倒れ込みの瞬間に、舌（扇形部）の最奥部１１ｅが５ｍｍ押し込みを行なう運動を組み込んだ。

［嚥下シミュレーター］
本実施の形態に係る嚥下シミュレーション方法を実行するためのシミュレーター（解析ソフトウェア）は頭頸部器官のモデリングを行い、口腔ならびに咽喉を通過する際の流体・食塊の挙動について粒子法を用いて解析を行うものである。
シミュレーターによる解析結果から、例えば、
（ａ）経口摂取品の物性値の違いによる、嚥下・誤嚥・誤飲・窒息リスクの推定、
（ｂ）経口摂取品の物性値の違いによる、嚥下時間の推定、
（ｃ）経口摂取品の物性値の違いによる、咽喉壁（気管や咽頭の管璧）にかかる力、およびせん断応力の推定、
（ｄ）上記データと官能評価又は官能診断との相関から、経口摂取品の飲み易さ・食べ易さ、飲み難さ・食べ難さの評価又は診断を行う。
評価又は診断は、評価者又は診断者が行なう場合と、嚥下シミュレーション装置が自動的に行なう場合がある。本実施例では評価者又は診断者が行なう。

嚥下シミュレーターは三次元粒子法解析ソフトを用いて作成した。本解析ソフトは流体の物性値や時間などを数値として直接入力でき、それぞれの物理量を適宜変えることができ、定量的な解析を行えるという特徴を有する。

［解析事例１］
組み込んだ嚥下運動に様々な物性の経口摂取品（食塊）を投入することによって、嚥下時の食塊の挙動を解析する。まず、疑似経口摂取品２０として水を取り上げ、水の物性値を代入して解析を行った。
表１に解析に用いた水の物性値を示す。二次元シミュレーターによる解析とは違い、物性値に単位のある数値を代入している。解析に用いた粒子の粒子径を１ｍｍとし、約３ｍｌの食塊を投入した。なお、粒子の粒子径は運動解析部５０で任意に設定可能である。

図１２に食塊が水の場合の三次元シミュレーションの解析結果を示す。図１２（ａ１）〜（ａ３）に舌１１上の食塊の挙動を、図１２（ｂ１）〜（ｂ３）に喉頭蓋１２ａ周辺での食塊の挙動を示す。（ａ１）はｔ＝０ｓｅｃ、（ａ２）はｔ＝０．４ｓｅｃ、（ｂ１）はｔ＝０．５ｓｅｃ、（ｂ２）はｔ＝０．６ｓｅｃ、（ｂ３）はｔ＝０．７ｓｅｃ、（ａ３）はｔ＝０．８ｓｅｃにおける水の挙動を示す。これらの図より、舌１１による輸送、喉頭蓋１２ａの倒れこみによる食道１８への食塊の投入が再現できていることが判る。また、一部の粒子は舌１１上や喉頭蓋１２ａ周辺に残留し、食道１８に流下せずに喉頭１２側に流入した。これらの現象について粒子法による解析では粒子の数によって数値化することが可能であり、喉残り、誤嚥のリスクなどを定量的に評価又は診断できる可能性がある。

［解析事例２］
次に、経口摂取品として、３種類の液体、水、牛乳Ａ（通常の牛乳）、牛乳Ｂ（明治社製：おいしい牛乳（株式会社明治の登録商標です））を投入して挙動を解析した。
表２に、これら３種類の液体の物性を示す。

図１３に、経口摂取品（食塊）が水、牛乳Ａ、牛乳Ｂの場合の三次元シミュレーションの解析結果を示す。図１３（ａ１）〜（ａ３）はｔ＝０．５ｓｅｃ（同一時刻）における舌１１上の食塊の挙動であり、図１３（ｂ１）〜（ｂ３）はｔ＝０．７ｓｅｃ（同一時刻）における喉頭蓋１２ａ周辺での食塊の挙動である。また、図１３（ａ１），（ｂ１）は水、図１３（ａ２），（ｂ２）は牛乳Ａ、図１３（ａ３），（ｂ３）は牛乳Ｂの挙動である。経口摂取品の物性の違いで同一時間に挙動が異なることがシミュレーションから観察できることが判った。

［せん断応力および流下速度による評価又は診断］
次に、表２の物性の経口摂取品（食塊）について、せん断応力と流下速度の解析を行ない、定量的な解析結果を得ることとした。

図１４にせん断応力と流下速度についての解析領域を示す。領域Ａは口腔内の舌１１上の領域で、せん断応力を求めた。領域Ｂは咽頭空間内の食道入口１８ａの領域で、せん断応力と流下速度を求めた。

図１５に、領域Ａにおけるせん断応力の経時変化（０．１ｓｅｃ毎に出力）を示す。領域Ａでの平均せん断応力は、領域Ａ内にある粒子にＹ方向の法線を持つ面に対してＺ方向の力が加わった際の応力（応力テンソルのσ_ＹＺ成分）の平均値によって求めることができる。舌１１による食塊の輸送によって領域Ａ内の粒子が１５以下になった時間でせん断応力を０とした。同図より舌１１上のせん断応力のピークが２つあることが判る。ここで１つ目のせん断応力のピークは食塊投入時に食塊が口腔内に広がる際に生じたせん断応力であり、口腔内での広がりの指標となると考えられる。２つ目のせん断応力のピークは嚥下時の舌１１による食塊の輸送の際に生じる舌１１上でのせん断応力であり、のど越しの１つの指標となる値であると考えられる。

図１６に領域Ｂにおけるせん断応力の経時変化（０．１ｓｅｃ毎に出力）を、図１７に領域Ｂにおける流下速度（速度のｖ_Ｙ成分）の経時変化（０．１ｓｅｃ毎に出力）を示す。領域Ｂ内での平均せん断応力は、領域Ａ内にある粒子に、Ｚ方向の法線を持つ面に対してＹ方向の力が加わった際の応力（応力テンソルのσ_ＺＹ成分）の平均値によって求めることができる。粒子が領域内で１５以下になった時間は、せん断応力および流下速度を０とした。

図１６より水と牛乳とでせん断応力に違いがあることが判る。また、牛乳Ａと牛乳Ｂでは嚥下終了タイミング（ｔ＝０．８ｓｅｃ）で大きな違いが現れた。この原因は図１７に示したように、牛乳Ｂの方がｔ＝０．８ｓｅｃでの流下速度が大きいためと考えられる。牛乳Ｂと牛乳Ａとの物性値に大きな違いはないが、液を流し込むタイミングでの挙動（液切れ、最後の数滴の飲み込み）に違いがある可能性が示唆された。
このようにして嚥下シミュレータでは微小な物性値の違いでの食塊の喉に及ぼす影響を数値化することが可能であり、これらを官能検査などの結果と比較することで、官能に及ぼす要因を考察することが可能となる。さらにせん断応力や流速などの視点から既存製品の物性による影響をマッピングすることができる。

［残存率と誤嚥率による評価又は診断］
次に、表２の物性の経口摂取品（食塊）について、残存率と誤嚥率の解析を行ない、定量的な解析結果を得ることができた。
図１８に残存率と誤嚥率についての解析領域を示す。領域Ｃは口腔１７内の舌１１上の領域で残存率を求めた。領域Ｄは喉頭蓋１２ａ周辺の領域で残存率を求めた。領域Ｅは喉頭１２出口の領域で、誤嚥率を求めた。

図１９に領域Ｃ及び領域Ｄ（両領域の合計）における残存率を、図２０に領域Ｅにおける誤嚥率を示す。粒子法は、粒子１つ１つの位置を追跡しているため、空間に存在する粒子の量を容易に求める事ができるという利点がある。そこで、嚥下動作後の舌１１上や喉頭蓋１２ａ周辺に残存した粒子の数を計数し、喉残りの指標（残存率）として評価又は診断した。さらに、喉頭１２の出口（図１８に破線で示す）に蓋をすることで喉頭１２側に混入した粒子についても計数し、誤嚥するリスク（誤嚥率）についても評価又は診断した。ただし、粒子を数えての評価方法又は診断方法については二次元嚥下シミュレータでも可能な手法であるため、三次元の利点は投入する食塊の物性値を数値で定義できるという点である。

図１９で解るように、牛乳が水よりも口腔内に残る割合が低いという結果となった。これは主に粘度が高いことによって塊を形成しやすいためと考えられた。また図２０において、牛乳Ｂの誤嚥率が低い結果となったが、この要因については未だ特定できていない。
［非ニュートン流体の解析］

図２１に食塊が非ニュートン流体の場合の三次元シミュレーションの解析結果を示す。図２１（ａ１）〜（ａ３）はｔ＝０．５ｓｅｃ（同一時刻）における舌１１上の食塊の挙動であり、図２１（ｂ１）〜（ｂ３）はｔ＝０．７ｓｅｃ（同一時刻）における喉頭蓋１２ａ周辺での食塊の挙動である。また、図２１（ａ１），（ｂ１）は濃度（水溶液中のとろみ調整食品の濃度）Ｃ＝１％、図２１（ａ２），（ｂ２）は濃度Ｃ＝２％、図２１（ａ３），（ｂ３）は濃度Ｃ＝３％の食塊の挙動である。
誤嚥率などの測定ができるようになると、とろみ調整食品の評価又は診断を行うことが可能となる。とろみ調整食品のような非ニュートン流体のレオロジー特性を持つものについても三次元嚥下シミュレータによって解析することが可能である。ここでは、とろみ調整食品（トロメイク（株式会社明治の登録商標です）ＳＰ（明治社製））について濃度１．０、１．５、２．０、３．０、３．５％（重量％）の水溶液のレオロジー特性および表面張力を測定し、その物性値での解析を行った。

［連成解析］
図２２に固体と流体（水）の連成解析の例を示す。図２２（ａ１）〜（ａ３）はｔ＝０．４ｓｅｃ（同一時刻）における舌１１上の食塊の挙動であり、図２２（ｂ１）〜（ｂ３）はｔ＝０．５ｓｅｃ（同一時刻）における喉頭蓋１２ａ周辺での食塊の挙動である。また、図２２（ａ１），（ｂ１）は固体（８ｍｍ球）のみ、図２１（ａ２），（ｂ２）は固体（４ｍｍ球）のみ、図２２（ａ３），（ｂ３）は固体（４ｍｍ球）と水との連成時の挙動である。
二次元シミュレータと同様にシミュレータ内に固体と流体を投入し、連成解析することも可能である。ここでは８ｍｍ球、４ｍｍ球および４ｍｍ球−水の連成計算を行った。剛体球の密度は１１，００ｋｇ／ｍ^３とし、喉壁面と剛体球の摩擦係数を０とした。

図２２（ｂ１）で解るように８ｍｍの固体球は大きすぎるために舌１１から咽頭１４側への輸送ができなかった。実際の人体の場合は無理をすれば飲み込める可能性があるが、口腔１７内の構造上は８ｍｍ程度の大きさのものはスムーズには輸送できないことが予想される。また、図２２（ｂ２）より４ｍｍ球に対しては問題なく輸送され、図２２（ｂ３）で解るように水との連成計算も問題なく実施できることが確認された。固体と流体の連成解析についても、固体や流体の物性を数値として設定可能となった。

図２３に実施例１における嚥下シミュレーション方法の処理フロー例を示す。まず、頭頸部器官からなる動的三次元頭頸部モデル１０ａを形成する（Ｓ０１０：頭頸部モデリング工程）。次に、頭頸部モデル１０ａにおける各頭頸部器官の運動を設定する（Ｓ０３０：器官運動設定工程）。次に、解析対象としての経口摂取品及びその物性（物性の種類及び物性値）を設定する（Ｓ０４０：経口摂取品物性設定工程）。これらの設定内容は状況に応じて自由に選定できる。設定内容は記憶部８３に記憶される。次に、経口摂取品をモデル化した擬似経口摂取品２０を口腔に入力する（Ｓ０５０：入力工程）。入力は例えば評価者又は診断者がマウスでカーソルを口腔内にドラッグ・アンド・ドロップして行う。次に頭頸部モデル１０ａにおける各頭頸部器官の運動と擬似経口摂取品２０の嚥下に係る挙動を、粒子法を用いて三次元空間で解析する（Ｓ０６０：運動解析工程）。例えばＭＰＳ法を使用できる。次に、運動解析工程（Ｓ０６０）で得られた解析結果を表示する（Ｓ０７０：表示工程）。次に、嚥下時の擬似経口摂取品２０の挙動の解析結果から、経口摂取品の食べ易さ及び／又は飲み易さの評価又は診断を行う（Ｓ０８０：評価診断工程、一方だけのときは、評価工程又は診断工程）。評価又は診断は評価者又は診断者が表示部８２の動画面を見て行なう。例えば表示部８２に表示された評価表又は診断表のセルに○×、ランク、点数等を入力する。評価又は診断を行なったら、経口摂取品物性設定工程（Ｓ０４０）に戻り、経口摂取品の物性値を変更して設定し、それ以降の評価診断工程までを繰り返し行う。変更すべき物性値は評価者又は診断者の判断で自由に選択できる。尤も１回目で適性な物性が見出されれば、その後の設定と評価又は診断を省略しても良い。次に、評価診断工程（Ｓ０８０）にて適正とされた経口摂取品の物性を決定する（Ｓ０９０：物性決定工程）。ここでは、適正な物性の範囲を示しても良く、適正な物性をランク分けしても良く、最適値を選定しても良い。

評価項目又は診断項目は、例えば次のようなものである。
（ａ）嚥下・誤嚥・誤飲・窒息リスク（口蓋壁等に付着して剥れ難い、咽喉や食道１８を塞ぐ、喉頭１２に入る）が有るか否か、
（ｂ）嚥下時間がどの位か、閾値を超えるか、
（ｃ）咽頭壁にかかる応力、およびせん断応力がどの位か、閾値を超えるか、
（ｄ）（ａ）〜（ｃ）に基づき、別にデータ取得した官能評価又は官能診断（美味しい、爽快感等）との相関性を考慮して、総合的に飲み易さ・食べ易さ、飲み難さ・食べ難さを評価又は診断する。
官能評価又は官能診断は、経口摂取品及びその物性と対応付けられて評価診断結果記録部８３Ｂに保存されており、この経口摂取品の官能評価又は官能診断を読み出して例えば嚥下時間、せん断応力等との相関性を考慮する。

以上説明したように、本実施例によれば、頭頸部モデル１０ａについて、器官特性、頭頸部器官の運動、経口摂取品の物性を設定し、粒子法を用いて経口摂取品の挙動を解析するので、嚥下に関する実現象、すなわち、頭頸部器官及び経口摂取品の挙動を近似的に再現することが容易な嚥下シミュレーション装置及び嚥下シミュレーション方法を提供できる。また、嚥下シミュレータを用いて食塊の物性や物理量を数値化して解析することにより、嚥下時における頭頸部器官及び経口摂取品の挙動を正確に表現でき、経口摂取品の物性及び挙動に係る物理量を定量化できる嚥下シミュレーション装置及び嚥下シミュレーション方法を提供できる。

実施例１では、評価者又は診断者が経口摂取品の入力を行い、動画像を見て嚥下評価又は嚥下診断を行なう例を説明したが、実施例２では、嚥下シミュレーション装置が設定に基づいて経口摂取品等の入力を自動的に行ない、嚥下シミュレーション装置が自動的に嚥下評価又は嚥下診断を行なう例を説明する。実施例１と異なる点を主に説明する（以下の実施例についても同様に異なる点を主に説明する）。

図２４に実施例２における嚥下シミュレーション装置１００Ｂの構成例を示す。実施例１（図２参照）に比して、パーソナルコンピュータ（ＰＣ）内に、、擬似経口摂取品の嚥下時の挙動の解析結果を仮想動画面に表示する擬似画面表示部８２Ａと擬似経口摂取品の入力条件を設定する経口摂取品入力設定部４５が追加され、さらに、評価診断条件記憶部８３Ａ内に記憶された評価条件又は診断条件を参照して自動的に経口摂取品の食べ易さ及び／又は飲み易さを評価又は診断する評価診断部６０（一方だけのときは、評価部又は診断部）が追加されている。その他の構成は実施例１と同様である。

図２５に嚥下シミュレーション方法の処理フロー例を示す。実施例１（図２３参照）に比して、入力工程（Ｓ０５０）の前に経口摂取品の入力条件を設定する為の経口摂取品入力設定工程（Ｓ０４５）が追加され、表示部８３に動画面表示する表示工程（Ｓ０７０）が、擬似画面表示部８２Ａの仮想動画面に表示する工程（Ｓ０７５）に代わっている。また、評価診断工程（Ｓ０８０）では、評価診断部６０が自動的に評価又は診断する。その他の工程は実施例１と同様である。

実施例２では、経口摂取品入力設定部４５に予め経口摂取品の投入位置と投入タイミングを設定しておく（Ｓ０４５：経口摂取品入力設定工程）。擬似経口摂取品の口腔内への投入位置は例えば口腔内で歯の近傍（例えば擬似経口摂取品の長さの２倍以内）とする。次に設定条件（場所、タイミング）に従って口腔内に擬似経口摂取品２０が投入される（Ｓ０５０：経口摂取品入力工程）。また、自動評価又は自動診断する場合には、予め評価条件又は診断条件を評価診断条件記憶部８３Ａに記憶しておき、パーソナルコンピュータＰＣ内の擬似画面表示部８２Ａの仮想動画面にシミュレーションの解析結果としての頭頸部モデル１０ａと擬似経口摂取品２０の挙動を表示して、評価診断部６０にて、擬似画面表示部８２Ａの表示を評価診断条件記憶部８３Ａの評価条件又は診断条件と照合することにより、評価又は診断を行なう。

評価項目又は診断項目は、例えば次のようなものである。
（ａ）嚥下・誤嚥・誤飲・窒息リスク（口蓋壁等に付着して剥れ難い、咽喉や食道１８を塞ぐ、喉頭１２に入る）が有るか否か、
（ｂ）嚥下時間がどの位か、閾値を超えるか、
（ｃ）咽頭壁にかかる応力、およびせん断応力がどの位か、閾値を超えるか、
（ｄ）（ａ）〜（ｃ）に基づき、別にデータ取得した官能評価又は官能診断（美味しい、爽快感等）との相関性を考慮して、総合的に飲み易さ・食べ易さ、飲み難さ・食べ難さを評価又は診断する。
（ａ）〜（ｃ）と官能評価又は官能診断をそれぞれ数値化しておき、それぞれ重み係数を掛けて合計数で自動的に総合評価又は総合診断する。なお、（ｃ）及び／又は（ｄ）を省略しても良い。

その他の構成及び処理フローは実施例１と同様であり、実施例１と同様に、嚥下に関する実現象を再現することが容易な嚥下シミュレーション装置及び嚥下シミュレーション方法を提供できる。

また、入力と評価又は診断の一方を人手で、他方をコンピュータで行なう場合でも、同様に実施可能であり、同様の効果を奏する。

以上の実施例では、評価又は診断を評価者又は診断者が行なう例と嚥下シミュレーション装置が自動評価又は自動診断を行なう例を説明したが、嚥下シミュレーション装置の評価結果又は診断結果（一部又は全体について）を表示して、さらに評価者又は診断者に評価又は診断を求めるようにしても良い。この場合、処理フローについては実施例２の評価診断工程（Ｓ０８０）において、表示部８２に評価表又は診断表と共に嚥下シミュレーション装置の評価結果又は診断結果（一部又は全体について）を表示する。評価者又は診断者は嚥下シミュレーション装置の評価結果又は診断結果を参照して、評価表又は診断表に自己の評価結果又は診断結果を入力する。経口摂取品入力工程（Ｓ０５０）は人手入力でも自動入力でも良い。その他の装置構成及び処理フローは実施例２と同様であり、実施例２と同様に、嚥下に関する実現象を再現することが容易な嚥下シミュレーション装置及び嚥下シミュレーション方法を提供できる。

以上の実施例では、物性決定部７０で適正とされる経口摂取品の物性を自動的に決定する例を説明したが、本実施例では人が決定する例を説明する。本実施例の装置構成については、典型的には実施例１における図２の物性決定部７０が削除される。処理フロー例については図２３を使用できる。なお、物性決定部７０があっても良く、この場合は使用しないか、物性決定部７０での決定結果を決定者（例えば評価者又は診断者）に参考として示す。決定は適正とされる経口摂取品の物性を人が決定するが、評価結果又は診断結果に基づいて行なうことに差異はないので、思考過程で多少変更される可能性はあるが、ほぼ同様な結果が予測される。その他の装置構成及び処理フローは実施例１と同様であり、実施例１と同様に、嚥下に関する実現象を再現することが容易な嚥下シミュレーション装置及び嚥下シミュレーション方法を提供できる。実施例２及び実施例３に対しても人が決定する場合、実施例１と同様のことがいえる。

［経口摂取品の開発］
次に、物性値を変化させてシミュレーションを行ない、シミュレーション結果を経口摂取品の開発に結びつける例を説明する。

図２６に経口摂取品開発支援装置２００Ａの構成例を示す。図２の嚥下シミュレーション装置１００Ａに、経口摂取品試作結果記録部８３Ｃ、決定製造条件記録部８３Ｄと製造条件決定部８４が追加されている。経口摂取品試作結果記録部８３Ｃは、記憶部８３に属し、物性決定部７０で適正と決定された物性を有するように製造条件を適宜に設定して試作を行なった結果（物性値を含む）を記録する。製造条件決定部８４は、パーソナルコンピュータＰＣ内に設けられ、経口摂取品試作結果記録部８３Ｃに記録された試作結果に基づき、経口摂取品の物性を物性決定部７０で適正と決定された物性にする製造条件を決定する。決定された製造条件は記憶部８３の決定製造条件記録部８３Ｄに記録される。これにより、物性決定部７０で適正と決定された物性にする製造条件が明確になり、食べ易い及び／又は飲み易い物性の経口摂取品を確実に開発することができる。

図２７に経口摂取品開発方法の処理フロー例を示す。図２３の嚥下シミュレーション方法の処理フロー例の後に経口摂取品試作工程（Ｓ１１０）及び製造条件決定工程（Ｓ１２０）が追加されている。経口摂取品試作工程（Ｓ１１０）では、物性決定工程（Ｓ０９０）で適正と決定された物性を有するように製造条件を適宜に設定して試作を行なうものである。例えば、菓子について、原料の配合条件（配合比、攪拌等）、焼成条件（温度、時間、雰囲気等）、冷却条件（温度、時間、雰囲気等）、寸法等を適宜に設定して試作を行なう。試作の結果は、得られた物性を製造条件と関連付けて記憶部８３の経口摂取品試作結果記録部８３Ｃに記録される。製造条件決定工程（Ｓ１２０）では、経口摂取品試作工程（Ｓ１１０）での試作結果、すなわち経口摂取品試作結果記録部８３Ｃに記録された試作結果に基づき、経口摂取品の物性を物性決定工程（Ｓ０９０）で適正と決定された物性にする製造条件を決定する。例えば、経口摂取品試作工程（Ｓ１１０）で、原料の配合比Ａ、焼成温度Ｂのときに、適正と決定された物性になるとの結果が得られていれば、適正と決定された物性にする製造条件として、原料の配合比Ａ、焼成温度Ｂと決定する。この製造条件は、１点でも良いが、範囲を指定しても良い。また、最適値を定めても良く、適正な範囲をランク分けしても良い。また、この決定は人が行なっても良く、製造条件決定部８４で経口摂取品試作結果記録部８３Ｃに記録された試作結果に基づいて自動的に行なっても良い。製造条件決定工程（Ｓ１２０）で決定された結果は、例えば記憶部８３の決定製造条件記録部８３Ｄに記録される。なお、人が決定する場合はノートに記録しても良い。

［経口摂取品の製造方法］
図２８に経口摂取品製造方法の処理フロー例を示す。経口摂取品製造工程の中に経口摂取品調製工程（Ｓ１３０）がある。最初をＡＡ工程、最後をＺＺ工程とする。この経口摂取品調製工程（Ｓ１３０）は、図２７における経口摂取品開発方法の製造条件決定工程（Ｓ１２０）において、物性決定工程（Ｓ０９０）で適正と決定された物性にする製造条件であると決定された製造条件を用いて経口摂取品を調製するものである。製造工程の中のいずれの工程でも良い。例えば菓子の製造では、原料配合工程や焼成工程等の工程で経口摂取品を調製することができる。２以上の工程を組み合わせて調製しても良い。このとき、この経口摂取品調製工程（Ｓ１３０）における製造条件が、製造条件決定工程（Ｓ１２０）において決定されたものである場合には、当該実施例に該当することとなる。製造条件決定工程（Ｓ１２０）で決定された結果は、例えば記憶部８３の決定製造条件記録部８３Ｄに記録される。これにより、食べ易い及び／又は飲み易い物性の開発経口摂取品を確実に製造することができる。

以上により、本実施例によれば、頭頸部モデル１０ａについて、器官特性、頭頸部器官の運動、経口摂取品の物性を設定し、粒子法を用いて経口摂取品の挙動を解析するので、嚥下に関する実現象を再現することが容易な嚥下シミュレーション方法を用いて嚥下現象を解析し、食べ易い又は飲み易い経口摂取品を開発することができる。

［食育装置・食育方法］
本実施例では本発明による嚥下シミュレータを食育に適用する例を説明する。
図２９に本実施例における食育支援装置３００Ａの構成例を示す。実施例１のシミュレーション装置１００Ａ（図２参照）に比して、教示部８５が追加され、物性決定部７０が削除されている。教示部８５は、表示部８２で動画面に表示された擬似経口摂取品２０の嚥下時の挙動を、評価診断結果記録部８３Ｂに記録された経口摂取品の評価結果又は診断結果と関連させて解説する。本実施例では、教示内容は評価診断結果記録部８３Ｂに記録された経口摂取品の評価結果又は診断結果に基づくもので、解説内容が予め作成され記録されている。解説内容は評価結果又は診断結果に基づいて教示部８５で自動的に作成しても良いが、教育者が更にわかり易く編集・補充するのが好適である。また、本実施例では、物性決定部７０が削除されているが、物性決定部７０を有し、物性決定について解説しても良い。なお、教示部８５は解説内容（教示内容）を記憶すると共に、運動解析部５０で解析された各頭頸部器官の運動と擬似経口摂取品２０の嚥下時の挙動の解析結果を表示部８２の動画面に表示させ、解説内容を表示部８２のスピーカ等に音声出力する。

図３０に本実施例における食育方法の処理フロー例を示す。実施例１のシミュレーション方法（図２３参照）に比して、評価診断工程（Ｓ０８０）の後に教示工程（Ｓ０８２）が追加され、物性決定工程（Ｓ０９０）が削除されている。教示工程（Ｓ０８２）では、表示工程（Ｓ０７０）で動画面に表示された擬似経口摂取品の嚥下時の挙動を、評価診断工程（Ｓ０８０）で評価又は診断された経口摂取品の評価結果又は診断結果と関連させて教示する。本実施例では、物性決定部７０及び物性決定工程（Ｓ０９０）が削除されているが、、物性決定部７０及び物性決定工程（Ｓ０９０）を有し、物性決定について解説しても良い。これにより、嚥下に関する実現象を再現することが容易な嚥下シミュレーション方法を用いて嚥下現象を表示するので、経口摂取品の食べ易さ又は飲み易さを理解しやすく、食育に大変有効である。

［診断支援］
本実施例では本発明による嚥下シミュレータを嚥下診断の支援に適用する例を説明する。

図３１に本実施例における嚥下シミュレーション装置４００Ａの構成を示す。実施例１（図２参照）の嚥下シミュレーション装置１００Ａに器管運動決定部７５及び決定器官運動条件記録部８３Ｅが追加されている。また、器管運動設定部３０がより多く使用される。器管運動設定部３０では、嚥下運動に関連する各頭頸部器官の運動特性を設定するが、例えばオトガイ舌筋１１ｆやその他の嚥下に関連する筋肉の反応速度、収縮と弛緩のタイミング、収縮距離、弾力性（しなやかさ）等を運動パラメータとして設定し、器官運動決定部７５では、シミュレーションの結果から、すなわち、運動解析部５０で解析された解析結果から、各頭頸部器官の器官運動パラメータを決定する。例えば、舌の進行波的波動運動が遅くなれば、嚥下に到るまでの時間がかかる。喉頭蓋１２ａの反応が遅れれば、経口摂取品が喉頭１２ａを経て気管１３に入り込み、誤嚥が生じるおそれがある。これにより、患者や健康診断の被診断者の各頭頸部器官の挙動や症状に良く合う各頭頸部器官の器官運動パラメータを求めることができる。決定器官運動条件記録部８３Ｅには、器官運動決定部７５で求められた器官運動パラメータが記録される。

そして、頭頸部モデル１０ａと嚥下シミュレーションの結果から、患者や健康診断対象者について、嚥下に関して機能が低下した筋肉がないか等の診断を行ない、治療に役立てることができる。なお、診断では、経口摂取品の物性より器官機能を重視するので、物性についてはループ処理行わずに固定としても良い。かかるシミュレーション装置を嚥下診断支援装置に組み込むことも可能である。例えば、患者や健康診断対象者についての診断結果を保存した医療診断結果データベースを有する嚥下診断支援装置とすれば、医療診断結果を嚥下シミュレーション装置４００Ａの評価診断結果記録部８３Ｂに記録された評価診断結果と比較することにより各頭頸部器官の機能低下した部分を発見して、迅速な診断を進めることも可能になる。

図３２に本実施例におけるシミュレーション方法の処理フロー例を示す。実施例１（図２３参照）の嚥下シミュレーション方法に、器管運動パラメータを変化させたループが追加されている。すなわち、物性のパラメータを変えたループ処理の後に、器管運動パラメータを変化させたループ処理が行われる。なお、ここでの物性決定や器官特性決定は、嚥下に適正な物性や器官運動を見出すのではなく、患者や健康診断の被診断者の症状に合う器管運動パラメータを探すループとして使用される。順次パラメータを変えてシミュレーションを行ない、患者や健康診断の被診断者の各頭頸部器官の挙動や症状に合う物性や器官特性が見出された時点で処理フローを終了する。尤も治療後のパラメータを予測してシミュレーションを継続し、治療の改善効果を求めることもできる。

［プログラム］
また、本発明は、以上の嚥下シミュレーション方法をコンピュータに実行させるためのコンピュータ読み取り可能なプログラムとしても、また、当該プログラムを記憶する記憶媒体としても実現可能である。プログラムは嚥下シミュレーション装置の制御部９０に蓄積して使用してもよく、内蔵又は外付けの記憶装置に蓄積して使用してもよく、インターネットからダウンロードして使用しても良い。

以上、本発明の実施の形態について説明したが、実施の形態は以上の例に限られるものではなく、本発明の趣旨を逸脱しない範囲で、種々の変更を加え得ることは明白である。

例えば、以上の実施例では、舌モデル、喉頭蓋モデルについて一例を示したが、これらの例に限られず、多様なモデルが可能である。例えば、時間毎に変化する各器官の形状を解析時間毎に読み込み、よりスムーズな運動を模擬したり、分割数や振動数等のパラメータを自由に変更した解析が可能である。さらに、これらの各器官は付随する筋肉と連動させて動かすことも可能である。また、実施例５の経口摂取品の開発、製造方法、実施例６の食育装置・食育方法については、主として実施例２のシミュレーション装置と組み合わせる例を説明したが、実施例１、実施例３、実施例４のシミュレーション装置と組み合わせても良い。また、実施例７お診断支援については、主として実施例１のシミュレーション装置と組み合わせる例を説明したが、実施例２ないし実施例４のシミュレーション装置と組み合わせても良い。また、経口摂取品は２個までを例示したが、３個以上を連動させて挙動を解析することも可能であり、また、ピーナッツを包んだチョコレートのように異なる物性値を持つ固体同士（固―固）の連成解析、リキュールを包んだチョコレート（固―液）の連成解析、さらにドレッシングのような液体同士（オイルと酢）の物性値が異なる（液―液）混合液の連成解析も可能である。また、唾液が加わる場合は、例えば予め加わる時期（複数回）と量を設定しておくこと及び経口摂取品が唾液により溶解することを取り込むことにより解析可能である。また、舌モデルの分割数、オーバーラップ部分の寸法、各部の移動量や移動時期、その他のパラメータ等を適切な範囲で変更可能である。

本発明は、経口摂取品の嚥下状況の解析に利用される。

１０ 頭頸部モデリング部
１０ａ 動的三次元頭頸部モデル
１０ｂ 舌モデル
１０ｃ 喉頭部モデル
１１ 舌
１１ａ〜１１ｅ 舌の扇形部
１１ｆ オトガイ舌筋
１１ｇ 舌表面部
１２ 喉頭
１２ａ 喉頭蓋
１２ｃ 喉頭部
１３ 気管
１４ 咽頭
１４ａ 咽頭の璧
１４ｂ 咽頭の粘膜
１５ 口蓋
１５ａ 硬口蓋
１５ｂ 軟口蓋
１６ 顎
１６ａ オトガイ部
１７ 口腔
１８ 食道
１８ａ 食道入口部
２０ 擬似経口摂取品
３０ 器官運動設定部
４０ 経口摂取品物性設定部
４５ 経口摂取品入力設定部
５０ 運動解析部
６０ 評価診断部
７０ 物性決定部
７５ 器官運動決定部
８１ 入力部
８２ 表示部
８２Ａ 擬似画面表示部
８３ 記憶部
８３Ａ 評価診断条件記憶部
８３Ｂ 評価診断結果記録部
８３Ｃ 経口摂取品試作結果記録部
８３Ｄ 決定製造条件記録部
８３Ｅ 決定器官運動条件記録部
８４ 製造条件決定部
８５ 教示部
９０ 制御部
１００Ａ，１００Ｂ 嚥下シミュレーション装置
２００Ａ 経口摂取品開発支援装置
３００Ａ 食育支援装置
４００Ａ 嚥下シミュレーション装置
ＰＣ パーソナルコンピュータ

図１９で解るように、牛乳が水よりも口腔内に残る割合が低いという結果となった。これは主に粘度が高いことによって塊を形成しやすいためと考えられた。また図２０において、牛乳Ｂの誤嚥率が低い結果となったが、この要因については未だ特定できていない。

［非ニュートン流体の解析］
図２１に食塊が非ニュートン流体の場合の三次元シミュレーションの解析結果を示す。図２１（ａ１）〜（ａ３）はｔ＝０．５ｓｅｃ（同一時刻）における舌１１上の食塊の挙動であり、図２１（ｂ１）〜（ｂ３）はｔ＝０．７ｓｅｃ（同一時刻）における喉頭蓋１２ａ周辺での食塊の挙動である。また、図２１（ａ１），（ｂ１）は濃度（水溶液中のとろみ調整食品の濃度）Ｃ＝１％、図２１（ａ２），（ｂ２）は濃度Ｃ＝２％、図２１（ａ３），（ｂ３）は濃度Ｃ＝３％の食塊の挙動である。
誤嚥率などの測定ができるようになると、とろみ調整食品の評価又は診断を行うことが可能となる。とろみ調整食品のような非ニュートン流体のレオロジー特性を持つものについても三次元嚥下シミュレータによって解析することが可能である。ここでは、とろみ調整食品（トロメイク（株式会社明治の登録商標です）ＳＰ（明治社製））について濃度１．０、１．５、２．０、３．０、３．５％（重量％）の水溶液のレオロジー特性および表面張力を測定し、その物性値での解析を行った。

［連成解析］
図２２に固体と流体（水）の連成解析の例を示す。図２２（ａ１）〜（ａ３）はｔ＝０．４ｓｅｃ（同一時刻）における舌１１上の食塊の挙動であり、図２２（ｂ１）〜（ｂ３）はｔ＝０．５ｓｅｃ（同一時刻）における喉頭蓋１２ａ周辺での食塊の挙動である。また、図２２（ａ１），（ｂ１）は固体（８ｍｍ球）のみ、図２２（ａ２），（ｂ２）は固体（４ｍｍ球）のみ、図２２（ａ３），（ｂ３）は固体（４ｍｍ球）と水との連成時の挙動である。
二次元シミュレータと同様にシミュレータ内に固体と流体を投入し、連成解析することも可能である。ここでは８ｍｍ球、４ｍｍ球および４ｍｍ球−水の連成計算を行った。剛体球の密度は１１，００ｋｇ／ｍ３とし、喉壁面と剛体球の摩擦係数を０とした。

Claims (22)

1. 頭頸部器官からなる動的三次元頭頸部モデルを形成する頭頸部モデリング部と；
   前記動的三次元頭頸部モデルにおける各頭頸部器官の運動を設定する器官運動設定部と；
   解析対象としての経口摂取品及びその物性を設定する経口摂取品物性設定部と；
   前記経口摂取品をモデル化した擬似経口摂取品を口腔に入力する入力部と；
   前記動的三次元頭頸部モデルにおける各前記頭頸部器官の運動と前記擬似経口摂取品の嚥下時の挙動を、粒子法を用いて三次元空間で解析する運動解析部と；
   前記運動解析部で解析された各前記頭頸部器官の運動と前記擬似経口摂取品の前記嚥下時の挙動の解析結果を動画面に表示する表示部とを備える；
   嚥下シミュレーション装置。
2. 前記頭頸部モデリング部は、前記頭頸部器官をポリゴン又は粒子に設定し；
   前記運動解析部は、前記擬似経口摂取品をポリゴン又は粒子として取り扱う；
   請求項１に記載の嚥下シミュレーション装置。
3. 前記頭頸部モデリング部は、舌を前奥方向にオトガイ舌筋の起始である下顎骨のオトガイ部を要とするｎ個（ｎは２以上の整数）の扇形部に分割する構造とし；
   前記器官運動設定部は、各扇形部が半径方向に所定の位相差を持って連携して振動することにより前記擬似経口摂取品を奥方向に輸送するように進行波的波動運動を設定する；
   請求項２に記載の嚥下シミュレーション装置。
4. 前記頭頸部モデリング部は、前記各扇形部の振動により前記各扇形部の間に隙間ができないように前記各扇形部の表面部分が相互に重複し合う構造とし、さらに、左右両側が半径方向に延びて中央部分に窪みを有する舌表面部を前記各扇形部の表面部分に重複させ、静止状態の舌表面部の前記窪みの中で前記各扇形部が半径方向に振動する構造とし；
   前記器官運動設定部は、嚥下時に各扇形部が円周方向奥側に所定の角度回転し、嚥下最終時には前記舌の最奥の扇形部で前記擬似経口摂取品を押し込むように設定する；
   請求項３に記載の嚥下シミュレーション装置。
5. 前記頭頸部モデリング部は、前記動的三次元頭頸部モデルを咽頭、喉頭部と食道に分割し、前記喉頭部を喉頭蓋と喉頭に分割し、嚥下開始前には前記食道の食道入口を閉塞し、嚥下時には前記食道入口が拡張し、咽頭が収縮する構造とし；
   前記器官運動設定部は、嚥下時に、前記喉頭を下顎骨のオトガイ部方向に移動して前記食道入口を開放し、咽頭が収縮し、次いで前記喉頭蓋を奥側に所定の角度回転させることにより、喉頭入口を閉塞するように設定する；
   請求項３又は請求項４に記載の嚥下シミュレーション装置。
6. 前記経口摂取品物性設定部は、解析対象として異なる物性の液体、半固体又は固体の複数の擬似経口摂取品を設定し；
   前記運動解析部は、前記複数の擬似経口摂取品について自由表面と前記複数の擬似経口摂取品間の境界を定め、前記複数の擬似経口摂取品の連動挙動の解析を行なう；
   請求項１ないし請求項５のいずれか１項に記載の嚥下シミュレーション装置。
7. 前記動画面における前記擬似経口摂取品の前記嚥下時の挙動から、前記経口摂取品の食べ易さ及び／又は飲み易さを評価又は診断する評価診断部を備え；
   前記動画面は、前記運動解析部で解析された各前記頭頸部器官の運動と前記擬似経口摂取品の嚥下時の挙動の解析結果を模擬的に表示するために嚥下シミュレーション装置により仮想空間に形成された仮想動画面であり；
   前記評価診断部は、前記仮想動画面での擬似経口摂取品の挙動が予め定められた条件を満たすか否かで評価又は診断する；
   請求項１ないし請求項６のいずれか１項に記載の嚥下シミュレーション装置。
8. 前記擬似経口摂取品の前記嚥下時の挙動の解析結果から、前記経口摂取品の食べ易さ及び／又は飲み易さに係る評価結果又は診断結果を記録する評価診断結果記録部と；、
   前記評価診断結果記録部に記録された評価又は診断結果に基づいて適正とされる経口摂取品の物性を決定する物性決定部とを備える；
   請求項１又は請求項７に記載の嚥下シミュレーション装置。
9. 頭頸部器官からなる動的三次元頭頸部モデルであって；
   舌を前奥方向にオトガイ舌筋の起始である下顎骨のオトガイ部を要とするｎ個（ｎは２以上の整数）の扇形部に分割する構造とし；
   各扇形部が半径方向に所定の位相差を持って連携して振動することにより前記擬似経口摂取品を奥方向に輸送する進行波的波動運動を行なうように設定されている；
   動的三次元頭頸部モデル。
10. 前記各扇形部の振動により前記各扇形部の間に隙間ができないように前記各扇形部の表面部分が相互に重複し合う構造とし、さらに、左右両側が半径方向に延びて中央部分に窪みを有する舌表面部を前記各扇形部の表面部分に重複させ、静止状態の舌表面部の前記窪みの中で前記各扇形部が円周方向に振動する構造とし：
    嚥下時に各扇形部が円周方向奥側に所定の角度回転し、次いで前記舌の最奥の扇形部で前記擬似経口摂取品を押し込むように設定され；
    咽頭、喉頭部と食道に分割され、喉頭部は喉頭蓋と喉頭に分割され、嚥下開始前に前記食道の食道入口を閉塞し、嚥下時には食道入口が拡張し、咽頭が収縮する構造とし；
    嚥下時に、前記喉頭を下顎骨のオトガイ部方向に移動して食道入口を開放し、咽頭が収縮し、次いで前記喉頭蓋を奥側に所定の角度回転させることにより、食道入口を閉塞するように設定されている；
    請求項９に記載の動的三次元頭頸部モデル。
11. 頭頸部器官からなる動的三次元頭頸部モデルを形成する頭頸部モデリング工程と；
    前記動的三次元頭頸部モデルにおける各頭頸部器官の運動を設定する器官運動設定工程と；
    解析対象としての経口摂取品、及びその物性を設定する経口摂取品物性設定工程と；
    前記経口摂取品をモデル化した擬似経口摂取品を前記口腔に入力する入力工程と；
    前記動的三次元頭頸部モデルにおける各前記頭頸部器官の運動と前記擬似経口摂取品の嚥下時の挙動を、粒子法を用いて解析する運動解析工程と；
    前記運動解析工程で三次元空間で解析された各前記頭頸部器官の運動と前記擬似経口摂取品の前記嚥下時の挙動の解析結果を動画面に表示する表示工程とを備える；
    嚥下シミュレーション方法。
12. 前記擬似経口摂取品の前記嚥下時の挙動の解析結果から、前記経口摂取品の食べ易さ及び／又は飲み易さに係る評価を行う評価工程と；
    前記評価工程にて評価された評価結果に基づいて適正とされる経口摂取品の物性を決定する物性決定工程とを備える；
    請求項１１に記載の嚥下シミュレーション方法。
13. 請求項１１又は請求項１２に記載の嚥下シミュレーション方法をコンピュータに実行させるためのコンピュータ読み取り可能なプログラム。
14. 請求項８に記載の嚥下シミュレーション装置と；
    前記物性決定部で適正と決定された物性を有するように製造条件を適宜に設定して試作を行なった結果を記録する経口摂取品試作結果記録部と；
    前記経口摂取品試作結果記録部に記録された試作結果に基づき、前記経口摂取品の物性を前記物性決定部で適正と決定された物性にする製造条件を決定する製造条件決定部とを備える；
    経口摂取品開発支援装置。
15. 前記動画面における前記擬似経口摂取品の前記嚥下時の挙動から、前記経口摂取品の食べ易さ及び／又は飲み易さを評価又は診断する評価診断部を備え；
    前記動画面は、前記運動解析部で解析された各前記頭頸部器官の運動と前記擬似経口摂取品の嚥下時の挙動の解析結果を模擬的に表示するために嚥下シミュレーション装置により仮想空間に形成された仮想動画面であり；
    前記評価診断部は、前記仮想動画面での擬似経口摂取品の挙動が予め定められた条件を満たすか否かで評価又は診断する；
    請求項１４に記載の経口摂取品開発支援装置。
16. 請求項１２に記載の嚥下シミュレーション方法を備え；
    前記経口摂取品物性設定工程において、前記経口摂取品の物性を変更して設定し、それ以降の物性決定工程までを繰り返し行ない；
    又は、前記経口摂取品物性設定工程において、前記経口摂取品の物性を変更して設定し、それ以降の評価工程までを繰り返し行ない、その後に一括して物性決定工程を行ない；
    さらに、前記物性決定工程で適正と決定された物性を有するように製造条件を適宜に設定して試作を行なう経口摂取品試作工程と；
    前記経口摂取品試作工程の結果に基づき、前記経口摂取品の物性を前記物性決定工程で適正と決定された物性にする製造条件を決定する製造条件決定工程とを備える；
    経口摂取品開発方法。
17. 前記動画面は、前記運動解析工程で解析された各前記頭頸部器官の運動と前記擬似経口摂取品の嚥下時の挙動の解析結果を模擬的に表示するために嚥下シミュレーション装置により仮想空間に形成された仮想動画面であり；
    前記表示工程は、前記解析結果を前記仮想動画面に模擬的に表示する仮想表示工程であり；
    前記評価診断工程は、前記仮想表示工程で模擬的に表示された仮想動画面での擬似経口摂取品の挙動が予め定められた条件を満たすか否かで評価する；
    請求項１６に記載の経口摂取品開発方法。
18. 請求項１６又は請求項１７に記載の経口摂取品開発方法の製造条件決定工程にて決定された製造条件を用いて前記経口摂取品を製造する；
    経口摂取品製造方法。
19. 請求項１又は請求項７に記載の嚥下シミュレーション装置と；
    前記擬似経口摂取品の前記嚥下時の挙動の解析結果から、前記経口摂取品の食べ易さ、飲み易さの評価結果又は診断結果を記録する評価診断結果記録部と；
    前記表示部で動画面に表示された前記擬似経口摂取品の嚥下時の挙動を、前記評価診断結果記録部に記録された経口摂取品の評価結果又は診断結果と関連させて解説する教示部とを備える；
    食育支援装置。
20. 請求項１１に記載の嚥下シミュレーション方法と；
    前記擬似経口摂取品の前記嚥下時の挙動の解析結果から、前記経口摂取品の食べ易さ及び／又は飲み易さに係る評価を行う評価工程と；
    前記表示工程で動画面に表示された前記擬似経口摂取品の嚥下時の挙動を、前記評価工程で評価された経口摂取品の評価結果と関連させて解説する教示工程とを備える；
    食育方法。
21. 前記器官運動設定部で設定された器官運動パラメータについて、前記運動解析部で解析された前記解析結果から、前記被診断者の器官の挙動又は症状に合う器官運動パラメータを決定する器官運動決定部を備える；
    請求項１ないし請求項８のいずれか１項に記載の嚥下シミュレーション装置。
22. 請求項２１に記載の嚥下シミュレーション装置を備える；
    診断支援装置。

Images