JP3302122B2 - 表情しわの模擬方法および装置 - Google Patents
表情しわの模擬方法および装置Info
- Publication number
- JP3302122B2 JP3302122B2 JP24712093A JP24712093A JP3302122B2 JP 3302122 B2 JP3302122 B2 JP 3302122B2 JP 24712093 A JP24712093 A JP 24712093A JP 24712093 A JP24712093 A JP 24712093A JP 3302122 B2 JP3302122 B2 JP 3302122B2
- Authority
- JP
- Japan
- Prior art keywords
- analysis
- skin
- subcutaneous tissue
- finite element
- buckling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Medical Treatment And Welfare Office Work (AREA)
Description
【０００１】
【産業上の利用分野】本発明は、表情しわの模擬方法お
よび装置、具体的には、実人間の顔の表情、中でも、表
情しわを模擬する方法およびそのための装置に関する。
よび装置、具体的には、実人間の顔の表情、中でも、表
情しわを模擬する方法およびそのための装置に関する。
【０００２】
【従来の技術】顔は約３０前後の表情筋の選択的収縮に
よって各種の表情を作っているが、表情機構に関する力
学的側面からの研究は不十分である。加齢に伴うスキン
−皮下組織系の物性や厚さ等の変化および喜怒哀楽の各
種の表情に伴って人間の顔に生じる表情しわは、スキン
−皮下組織系の座屈および座屈後の有限変形現象とみな
せるが、これまで、実人間を対象とした顔表情、特に、
表情しわを模擬するための技術は未開発である。
よって各種の表情を作っているが、表情機構に関する力
学的側面からの研究は不十分である。加齢に伴うスキン
−皮下組織系の物性や厚さ等の変化および喜怒哀楽の各
種の表情に伴って人間の顔に生じる表情しわは、スキン
−皮下組織系の座屈および座屈後の有限変形現象とみな
せるが、これまで、実人間を対象とした顔表情、特に、
表情しわを模擬するための技術は未開発である。
【０００３】
【発明が解決しようとする課題】本発明は、生体外(in
vitro)座屈実験による解析によらず、非破壊での生体軟
組織の構成則の同定すなわち生体内（in vivo)皮膚の力
学物性の同定のための技術を用い、実人間の顔の表情し
わを、非破壊あるいは無侵襲により、模擬することがで
きる方法および装置を提供することを目的とする。
vitro)座屈実験による解析によらず、非破壊での生体軟
組織の構成則の同定すなわち生体内（in vivo)皮膚の力
学物性の同定のための技術を用い、実人間の顔の表情し
わを、非破壊あるいは無侵襲により、模擬することがで
きる方法および装置を提供することを目的とする。
【０００４】
【課題を解決するための手段】本発明による表情しわの
模擬方法は、非破壊かつ無侵襲によりスキン−皮下組織
の物性および厚さを計測し、該スキン−皮下組織の物性
の計測結果と同一系の有限要素解との比較に基づいてス
キン−皮下組織系の超弾性構成則における物性値を決定
し、決定された物性値を用いて有限要素法によりスキン
−皮下組織系の座屈および後座屈解析を行って座屈後挙
動を求めることにより、しわの生成および成長解析を行
うように構成される。
模擬方法は、非破壊かつ無侵襲によりスキン−皮下組織
の物性および厚さを計測し、該スキン−皮下組織の物性
の計測結果と同一系の有限要素解との比較に基づいてス
キン−皮下組織系の超弾性構成則における物性値を決定
し、決定された物性値を用いて有限要素法によりスキン
−皮下組織系の座屈および後座屈解析を行って座屈後挙
動を求めることにより、しわの生成および成長解析を行
うように構成される。
【０００５】そして、本発明による表情しわの模擬装置
は、三次元顔形状計測手段およびスキン−皮下組織の物
性および厚さ計測手段とを含む入力部と、その入力部か
らの物性データと同じ条件のスキン−皮下組織系の有限
要素解との比較により、そのスキン−皮下組織系の超弾
性構成則における物性値を決定する解析データ生成機能
と、その物性値を用いて有限要素法によりスキン−皮下
組織系の座屈および後座屈解析を行って座屈後挙動を求
める三次元有限要素法によるしわ発生およびしわ成長の
解析機能とを含む解析部と、少なくとも上記解析部によ
るしわ解析結果の画像表示を行うディスプレイ手段を含
む出力部とを備えている。
は、三次元顔形状計測手段およびスキン−皮下組織の物
性および厚さ計測手段とを含む入力部と、その入力部か
らの物性データと同じ条件のスキン−皮下組織系の有限
要素解との比較により、そのスキン−皮下組織系の超弾
性構成則における物性値を決定する解析データ生成機能
と、その物性値を用いて有限要素法によりスキン−皮下
組織系の座屈および後座屈解析を行って座屈後挙動を求
める三次元有限要素法によるしわ発生およびしわ成長の
解析機能とを含む解析部と、少なくとも上記解析部によ
るしわ解析結果の画像表示を行うディスプレイ手段を含
む出力部とを備えている。
【０００６】
【作用】この構成により、生体軟組織系の非圧縮・超弾
性体の有限変形問題の後座屈現象を有限要素法を用いて
解析することができ、したがって、非破壊かつ無侵襲で
得たデータに基づいて、実人間を対象とした顔表情、特
に、表情しわを模擬することができる。
性体の有限変形問題の後座屈現象を有限要素法を用いて
解析することができ、したがって、非破壊かつ無侵襲で
得たデータに基づいて、実人間を対象とした顔表情、特
に、表情しわを模擬することができる。
【０００７】
【実施例】以下、本発明による表情しわの模擬方法およ
び装置を、実人間の額部の表情しわを解析し、模擬する
実施例について、説明する。図１は、本発明による表情
しわの模擬方法および装置を実現するシステムの概略構
成を示すブロック図であり、１は必要な物性値や寸法値
を取り込むための入力部、２は有限要素法ＦＥＭ（fini
te element method)を用いて表情しわの生成を解析する
解析部、３は得られた解析結果により表情しわデータを
ディスプレイあるいは出力する出力部である。
び装置を、実人間の額部の表情しわを解析し、模擬する
実施例について、説明する。図１は、本発明による表情
しわの模擬方法および装置を実現するシステムの概略構
成を示すブロック図であり、１は必要な物性値や寸法値
を取り込むための入力部、２は有限要素法ＦＥＭ（fini
te element method)を用いて表情しわの生成を解析する
解析部、３は得られた解析結果により表情しわデータを
ディスプレイあるいは出力する出力部である。
【０００８】入力部１においては、例えば、 三次元の顔形状計測データ、 額部中の特定部分の或る表情時における移動を示す情
報すなわち変形情報データ、 額部の厚さデータ、 額部の印加荷重と変形との関係を示す特性曲線すなわ
ち構成則同定用データ等の必要なデータを取り込む。な
お、上記の他に、スキンの初期張力、筋収縮特性、初期
しわのピッチや深さ、顔の曲面効果等々の表情しわの生
成に必要な各種の情報を用いることができる。
報すなわち変形情報データ、 額部の厚さデータ、 額部の印加荷重と変形との関係を示す特性曲線すなわ
ち構成則同定用データ等の必要なデータを取り込む。な
お、上記の他に、スキンの初期張力、筋収縮特性、初期
しわのピッチや深さ、顔の曲面効果等々の表情しわの生
成に必要な各種の情報を用いることができる。
【０００９】そのため、三次元顔形状計測手段１１にお
いて、例えば撮像手段を用いて、顔の全体、部分の各種
方向からの立体形状、額部中の特定部分について、或る
表情をした時の動き等を表す表情変化情報等の画像デー
タを収集する。なお、この場合に、ズーミングを利用す
ることもできる。画像解析手段１２において、収集され
た画像データを用いて画像データ解析により立体画数値
データが生成されて、上記の三次元の顔形状計測デー
タおよび変形情報データが収集される。
いて、例えば撮像手段を用いて、顔の全体、部分の各種
方向からの立体形状、額部中の特定部分について、或る
表情をした時の動き等を表す表情変化情報等の画像デー
タを収集する。なお、この場合に、ズーミングを利用す
ることもできる。画像解析手段１２において、収集され
た画像データを用いて画像データ解析により立体画数値
データが生成されて、上記の三次元の顔形状計測デー
タおよび変形情報データが収集される。
【００１０】また、スキン−皮下組織の物性、厚さ計測
手段１３において、額部の寸法および力学物性データの
基礎となるデータを収集し、かつ、必要に応じて既存の
物性データを標準物性データベース１４を利用しつつ、
計測値解析手段１５により、得られた計測データを解
析、編集して、厚さデータおよび構成則同定のため
の基礎データが収集される。
手段１３において、額部の寸法および力学物性データの
基礎となるデータを収集し、かつ、必要に応じて既存の
物性データを標準物性データベース１４を利用しつつ、
計測値解析手段１５により、得られた計測データを解
析、編集して、厚さデータおよび構成則同定のため
の基礎データが収集される。
【００１１】ここで、額部の寸法、特に厚さデータの収
集に当たっては、例えば超音波測定技術を用いて、非破
壊かつ無侵襲により、額部の皮下構造すなわち皮膚−皮
下組織−筋肉−骨からなる構造の各肉厚情報を測定する
ことができる。図２は、スキン−皮下組織系の超弾性構
成則の同定のための基礎データを測定する態様の一例を
示す説明図であり、測定の対象とする額部について、例
えば、皮膚表面上に直径４０ｍｍの測定領域を定め、測
定領域に対応する円形の開口を有するスポーツ用粘着テ
ープを貼り付けて、それ以外の部分の変形を拘束する。
この状態において、同図（Ａ）に示されているような筋
収縮（上下）方向および左右方向の二軸引張試験、同図
（Ｂ）に示されているようなねじり試験等の構成則の同
定に必要な一連の載荷試験を行って必要な基礎データを
測定する。
集に当たっては、例えば超音波測定技術を用いて、非破
壊かつ無侵襲により、額部の皮下構造すなわち皮膚−皮
下組織−筋肉−骨からなる構造の各肉厚情報を測定する
ことができる。図２は、スキン−皮下組織系の超弾性構
成則の同定のための基礎データを測定する態様の一例を
示す説明図であり、測定の対象とする額部について、例
えば、皮膚表面上に直径４０ｍｍの測定領域を定め、測
定領域に対応する円形の開口を有するスポーツ用粘着テ
ープを貼り付けて、それ以外の部分の変形を拘束する。
この状態において、同図（Ａ）に示されているような筋
収縮（上下）方向および左右方向の二軸引張試験、同図
（Ｂ）に示されているようなねじり試験等の構成則の同
定に必要な一連の載荷試験を行って必要な基礎データを
測定する。
【００１２】以上のようにして得られた各種入力データ
を用いて、解析部２において、三次元ＦＥＭ解析による
解析データの生成機能１６により解析データを生成し、
そして、三次元ＦＥＭによるしわ発生およびしわ成長の
解析機能１７により表情しわの模擬が行われる。まず、
三次元ＦＥＭによる解析データの生成機能１６におい
て、収集された構成則基礎データから、例えば、次のよ
うな手法を用いて、スキン−皮下組織系の超弾性構成則
の同定が行われる。基礎データの収集と同じ条件におい
て、適当な皮膚物性を与えて有限要素（ＦＥＭ）解析を
行い、得られた荷重−変位曲線が測定の結果と一致する
ように、皮膚物性を変化させて、最急降下法または最小
自乗法により収束計算を行い、物性の同定を行う。
を用いて、解析部２において、三次元ＦＥＭ解析による
解析データの生成機能１６により解析データを生成し、
そして、三次元ＦＥＭによるしわ発生およびしわ成長の
解析機能１７により表情しわの模擬が行われる。まず、
三次元ＦＥＭによる解析データの生成機能１６におい
て、収集された構成則基礎データから、例えば、次のよ
うな手法を用いて、スキン−皮下組織系の超弾性構成則
の同定が行われる。基礎データの収集と同じ条件におい
て、適当な皮膚物性を与えて有限要素（ＦＥＭ）解析を
行い、得られた荷重−変位曲線が測定の結果と一致する
ように、皮膚物性を変化させて、最急降下法または最小
自乗法により収束計算を行い、物性の同定を行う。
【００１３】簡単化のため、第１ステップとしてスキン
と皮下組織の等価平均物性を求め、次いで、後述する有
限要素法プログラムを用いて生体内試験と同じ荷重およ
び変形拘束条件を与えれば、初期値として仮定した物性
値αi 、μi （ｉ＝１〜４）、ｃに対し、荷重ごとの変
位分布が計算される。そこで、試験により得られる変位
をｕi （バー）、計算によるものをｕi とし、この誤差
に対する繰り返し収束計算を最急降下法により行えば、
収束解として必要なスキン−皮下組織系の生体内物性が
定められる。
と皮下組織の等価平均物性を求め、次いで、後述する有
限要素法プログラムを用いて生体内試験と同じ荷重およ
び変形拘束条件を与えれば、初期値として仮定した物性
値αi 、μi （ｉ＝１〜４）、ｃに対し、荷重ごとの変
位分布が計算される。そこで、試験により得られる変位
をｕi （バー）、計算によるものをｕi とし、この誤差
に対する繰り返し収束計算を最急降下法により行えば、
収束解として必要なスキン−皮下組織系の生体内物性が
定められる。
【００１４】すなわち、収束解を得るためのアルゴリズ
ムとして、目標関数をｇi （α）〔但し、α（αi ，μ
i ，ｃ）〕とすれば、繰り返し収束計算のｋ＋１回目の
解α k+1 を求める関係は、
ムとして、目標関数をｇi （α）〔但し、α（αi ，μ
i ，ｃ）〕とすれば、繰り返し収束計算のｋ＋１回目の
解α k+1 を求める関係は、
【００１５】
【数１】
【００１６】但し、
【００１７】
【数２】
【００１８】なお、非線形問題のため、影響関数Ｈijは
次式により求める。
次式により求める。
【００１９】
【数３】
【００２０】ここで、後述するＦＥＭ解析における式
（５）から、αは低荷重線形域の係数α i と高荷重側の
指数べき部分μi に分けて多段階の繰り返し計算で扱
う。図３は、収束時の荷重−変位曲線の例を、測定値と
ともに、示しており、同図（Ａ）は二軸引張試験の結
果、同図（Ｂ）はねじり試験の結果を示している。ま
た、図４はその同定結果の一例を示す表である。
（５）から、αは低荷重線形域の係数α i と高荷重側の
指数べき部分μi に分けて多段階の繰り返し計算で扱
う。図３は、収束時の荷重−変位曲線の例を、測定値と
ともに、示しており、同図（Ａ）は二軸引張試験の結
果、同図（Ｂ）はねじり試験の結果を示している。ま
た、図４はその同定結果の一例を示す表である。
【００２１】このようにＦＥＭ解析法を用いて得られた
スキン−皮下組織系の超弾性構成則の同定データととも
に、入力部１において収集されている各種データを用い
て、図１における解析部２の三次元ＦＥＭによるしわ発
生およびしわ成長の解析機能１７において、後述するＦ
ＥＭ解析法によりスキン−皮下組織系の座屈および後座
屈現象を解析して表情しわの生成と成長を模擬し、同じ
く図１における出力部３において、得られたしわ解析結
果をディスプレイにより画像表示したり、あるいは、プ
リントアウトによりハードコピーする。
スキン−皮下組織系の超弾性構成則の同定データととも
に、入力部１において収集されている各種データを用い
て、図１における解析部２の三次元ＦＥＭによるしわ発
生およびしわ成長の解析機能１７において、後述するＦ
ＥＭ解析法によりスキン−皮下組織系の座屈および後座
屈現象を解析して表情しわの生成と成長を模擬し、同じ
く図１における出力部３において、得られたしわ解析結
果をディスプレイにより画像表示したり、あるいは、プ
リントアウトによりハードコピーする。
【００２２】図５は、上述したＦＥＭ解析法を用いて、
図１の解析部２における機能を実現するためのシステム
の動作フローの一例を説明する概念的なフローチャート
である。このフローチャートにおいて、ステップ１００
は図１の入力部１における計測・解析手段１１〜１５に
よる各種データの収集を示しており、ステップ１０１〜
１０７からなる構成則同定システムループが解析部２の
三次元ＦＥＭによる解析データの生成機能１６を実現
し、ステップ１１１〜１１６からなる表情変形・表情し
わ解析ループが三次元ＦＥＭによるしわ発生およびしわ
成長の解析機能１７を実現するためのものである。
図１の解析部２における機能を実現するためのシステム
の動作フローの一例を説明する概念的なフローチャート
である。このフローチャートにおいて、ステップ１００
は図１の入力部１における計測・解析手段１１〜１５に
よる各種データの収集を示しており、ステップ１０１〜
１０７からなる構成則同定システムループが解析部２の
三次元ＦＥＭによる解析データの生成機能１６を実現
し、ステップ１１１〜１１６からなる表情変形・表情し
わ解析ループが三次元ＦＥＭによるしわ発生およびしわ
成長の解析機能１７を実現するためのものである。
【００２３】以下、スキンと皮下組織の三次元ＦＥＭ解
析法について説明する。 (1) 構成方程式 座標系はスキン広がり面をｘ−ｙ面とし、ｘ軸を筋収縮
方向に合わせ、厚み方向をｚ軸と定める。生体軟組織系
の挙動は大変形であるから、グリーンの歪(Green strai
n)Ｅijと変形後の埋込座標系基準のキルヒホッフ応力(K
irchoff stress) Ｓijを用いることにすると、ｕi ：変
位、Ｘi 、ｘi ：変形前、後座標系として、グリーンの
歪は、
析法について説明する。 (1) 構成方程式 座標系はスキン広がり面をｘ−ｙ面とし、ｘ軸を筋収縮
方向に合わせ、厚み方向をｚ軸と定める。生体軟組織系
の挙動は大変形であるから、グリーンの歪(Green strai
n)Ｅijと変形後の埋込座標系基準のキルヒホッフ応力(K
irchoff stress) Ｓijを用いることにすると、ｕi ：変
位、Ｘi 、ｘi ：変形前、後座標系として、グリーンの
歪は、
【００２４】
【数４】
【００２５】歪エネルギー関数Ｗを
【００２６】
【数５】
【００２７】と表せば、キルヒホッフ応力Ｓijは一般に
【００２８】
【数６】
【００２９】ここで、試験データを参考にして、歪エネ
ルギー関数Ｗを次のように近似する。
ルギー関数Ｗを次のように近似する。
【００３０】
【数７】
【００３１】但し、
【００３２】
【数８】
【００３３】式（４）を式（３）に代入すれば、次の三
次元応力−歪関係
次元応力−歪関係
【００３４】
【数９】
【００３５】が得られる。但し、αi 、ｃ、μj は実験
定数である。 (2) 増分線形化近似 第ｎステップのグリーン歪の増分量ΔＥij n をつぎのよ
うにｎ−１ステップの既知解を用いて線形化近似する
と、式（１）は次式となる。
定数である。 (2) 増分線形化近似 第ｎステップのグリーン歪の増分量ΔＥij n をつぎのよ
うにｎ−１ステップの既知解を用いて線形化近似する
と、式（１）は次式となる。
【００３６】
【数１０】
【００３７】但し、Δｕi n は第ｎステップの変位増
分、ｕi n-1 はｎ−１ステップまでの全変位で、第ｎ増
分計算の初期値である。要素の変位関数をＮm とし、
分、ｕi n-1 はｎ−１ステップまでの全変位で、第ｎ増
分計算の初期値である。要素の変位関数をＮm とし、
【００３８】
【数１１】
【００３９】の関数を式（６）に用いてマトリックス形
で表すと、
で表すと、
【００４０】
【数１２】
【００４１】ここに、
【００４２】
【数１３】
【００４３】
【数１４】
【００４４】同様にして、増分応力−歪関係は、式
（５）を次式のようにマトリックス形で表して、
（５）を次式のようにマトリックス形で表して、
【００４５】
【数１５】
【００４６】これより、
【００４７】
【数１６】
【００４８】となる。 (3) 変分原理と増分型剛性方程式 非圧縮性の条件をラグランジュの未定係数ｈを用いて、
拘束条件とすれば、仮想仕事の原理は、
拘束条件とすれば、仮想仕事の原理は、
【００４９】
【数１７】
【００５０】但し、Ｔ ：表面力、Ｆ ：物体力、δ：
仮想量、Δ：増分、Ｉ ：グリーンの変形テンソルの行
列式である。以上により必要な関係が得られ、仮想量δ
ｕ（バー）、δｈの任意性から最終的に解くべき関係が
得られる。すなわち、
仮想量、Δ：増分、Ｉ ：グリーンの変形テンソルの行
列式である。以上により必要な関係が得られ、仮想量δ
ｕ（バー）、δｈの任意性から最終的に解くべき関係が
得られる。すなわち、
【００５１】
【数１８】
【００５２】(4) 座屈解析および後座屈解析 上記式（１０）の左辺および右辺で表される内部仕事と
外部仕事を、改めてδＵ（ｕi ）、δＷ（ｕi ）とし、
座屈発生の瞬間新たに付加される変位場をｕi * とす
る。荷重の増加なしに、座屈発生前の応力場Ｓij * が付
加変位ｕi * によりなす仕事をδＵ* （ｕi * ）とすれ
ば、座屈発生の瞬間には、 δＵ（ｕi ）＋δＷ（ｕi ）＋δＵ* （ｕi * ）＝０ （２０） 式（７）を考慮すれば、 δＵ* （ｕi * ）＝０ （２１） が成り立つ。これより、
外部仕事を、改めてδＵ（ｕi ）、δＷ（ｕi ）とし、
座屈発生の瞬間新たに付加される変位場をｕi * とす
る。荷重の増加なしに、座屈発生前の応力場Ｓij * が付
加変位ｕi * によりなす仕事をδＵ* （ｕi * ）とすれ
ば、座屈発生の瞬間には、 δＵ（ｕi ）＋δＷ（ｕi ）＋δＵ* （ｕi * ）＝０ （２０） 式（７）を考慮すれば、 δＵ* （ｕi * ）＝０ （２１） が成り立つ。これより、
【００５３】
【数１９】
【００５４】と表せる。系の特性が非線形であることか
ら、実際の計算では式（１１）の解と式（２２）を連立
させながら、座屈発生条件
ら、実際の計算では式（１１）の解と式（２２）を連立
させながら、座屈発生条件
【００５５】
【数２０】
【００５６】となる応力場Ｓij * または荷重〔ΔＰ* 〕
を求めていくこととなる。このような過程によって、式
（２３）を満足する座屈応力が定まれば、これを式（２
２）に用いて座屈モードｕi * （バー）が定められる。
系の座屈後挙動を求めるためには、得られたｕi * （バ
ー）を初期不整として式（１１）に用い、増分解析の初
期からこれを解くことにより必要な計算を行うことがで
きる。
を求めていくこととなる。このような過程によって、式
（２３）を満足する座屈応力が定まれば、これを式（２
２）に用いて座屈モードｕi * （バー）が定められる。
系の座屈後挙動を求めるためには、得られたｕi * （バ
ー）を初期不整として式（１１）に用い、増分解析の初
期からこれを解くことにより必要な計算を行うことがで
きる。
【００５７】
【発明の効果】図６は、実人間の表情しわの実験と本発
明による解析結果との比較を例示している。同中上部に
示す人間の額部において、眉吊り上げによって生じる額
横しわを対象として、スキン−皮下組織系の超音波によ
る厚さ分布の計測、生体内非破壊物性計測、無表情時の
顔の三次元形状のスリット光源法による計測、眉吊り上
げ状態で生じる額部表情変化（試験領域周辺の移動量）
の計測、および、眉吊り上げ時の表情しわの計測を行っ
た結果と、前頭筋走行方向断面内の二次元モデルとし
て、本発明により表情しわの発生、成長解析を行った結
果とが、比較して例示されている。これより明らかなよ
うに、しわの成長に対する解析と実測の結果は比較的一
致しており、本発明の解析法により生体における表情し
わを有効に推定することができる。
明による解析結果との比較を例示している。同中上部に
示す人間の額部において、眉吊り上げによって生じる額
横しわを対象として、スキン−皮下組織系の超音波によ
る厚さ分布の計測、生体内非破壊物性計測、無表情時の
顔の三次元形状のスリット光源法による計測、眉吊り上
げ状態で生じる額部表情変化（試験領域周辺の移動量）
の計測、および、眉吊り上げ時の表情しわの計測を行っ
た結果と、前頭筋走行方向断面内の二次元モデルとし
て、本発明により表情しわの発生、成長解析を行った結
果とが、比較して例示されている。これより明らかなよ
うに、しわの成長に対する解析と実測の結果は比較的一
致しており、本発明の解析法により生体における表情し
わを有効に推定することができる。
【図１】本発明による表情しわの模擬方法および装置を
実現するシステムの概略構成を示すブロック図である。
実現するシステムの概略構成を示すブロック図である。
【図２】スキン−皮下組織系の超弾性構成則の同定のた
めの基礎データを測定する態様の一例を示す説明図であ
る。
めの基礎データを測定する態様の一例を示す説明図であ
る。
【図３】スキン−皮下組織系の超弾性構成則の同定のた
めの荷重−変位曲線の例を示す特性図である。
めの荷重−変位曲線の例を示す特性図である。
【図４】スキン−皮下組織系の超弾性構成則の同定結果
の一例を示す表である。
の一例を示す表である。
【図５】本発明に従い、ＦＥＭ解析法を用いて、スキン
−皮下組織系の超弾性構成則の同定および表情しわの発
生・成長の解析機能を実現するためのシステムの動作フ
ローの一例を説明する概念的なフローチャートである。
−皮下組織系の超弾性構成則の同定および表情しわの発
生・成長の解析機能を実現するためのシステムの動作フ
ローの一例を説明する概念的なフローチャートである。
【図６】実人間の表情しわの実験と本発明による解析結
果との比較を例示する説明図である。
果との比較を例示する説明図である。
１…入力部 ２…解析部 ３…出力部 １１…三次元顔形状計測手段 １２…画像解析手段 １３…スキン−皮下組織の物性、厚さ計測手段 １４…標準物性データベース １５…計測値解析手段 １６…三次元ＦＥＭによる解析データの生成機能 １７…三次元ＦＥＭによるしわ発生およびしわ成長の解
析機能
析機能
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 平３−252775（ＪＰ，Ａ) 特開 平３−252776（ＪＰ，Ａ) (58)調査した分野(Int.Cl.7，ＤＢ名) A61B 5/107 A61B 5/00 G06T 1/00 340
Claims (2)
- 【請求項１】 非破壊かつ無侵襲によりスキン−皮下組
織の物性および厚さを計測し、該スキン−皮下組織の物
性の計測結果と同一系の有限要素解との比較に基づいて
スキン−皮下組織系の超弾性構成則における物性値を決
定し、決定された物性値を用いて有限要素法によりスキ
ン−皮下組織系の座屈および後座屈解析を行って座屈後
挙動を求めることにより、しわの生成および成長解析を
行うことを特徴とする表情しわの模擬方法。 - 【請求項２】 三次元顔形状計測手段およびスキン−皮
下組織の物性および厚さ計測手段とを含む入力部と、 上記入力部からの物性データと該スキン−皮下組織系の
有限要素解との比較により、該スキン−皮下組織系の超
弾性構成則における物性値を決定する三次元有限要素法
による解析データ生成機能と、該物性値を用いて有限要
素法によりスキン−皮下組織系の座屈および後座屈解析
を行って座屈後挙動を求める三次元有限要素法によるし
わ発生およびしわ成長の解析機能とを含む解析部と、 少なくとも上記解析部によるしわ解析結果の画像表示を
行うディスプレイ手段を含む出力部とを備えることを特
徴とする表情しわの模擬装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP24712093A JP3302122B2 (ja) | 1993-10-01 | 1993-10-01 | 表情しわの模擬方法および装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP24712093A JP3302122B2 (ja) | 1993-10-01 | 1993-10-01 | 表情しわの模擬方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH07100126A JPH07100126A (ja) | 1995-04-18 |
| JP3302122B2 true JP3302122B2 (ja) | 2002-07-15 |
Family
ID=17158735
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP24712093A Expired - Lifetime JP3302122B2 (ja) | 1993-10-01 | 1993-10-01 | 表情しわの模擬方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3302122B2 (ja) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6571003B1 (en) | 1999-06-14 | 2003-05-27 | The Procter & Gamble Company | Skin imaging and analysis systems and methods |
| JP2001306691A (ja) * | 2000-02-17 | 2001-11-02 | Shimizu Haruo | 診療情報システム |
| US6959119B2 (en) | 2000-11-03 | 2005-10-25 | Unilever Home & Personal Care Usa | Method of evaluating cosmetic products on a consumer with future predictive transformation |
| JP5156472B2 (ja) * | 2008-04-30 | 2013-03-06 | ラトックシステムエンジニアリング株式会社 | 椎体海綿骨の検査プログラムおよび検査装置 |
| JP5156516B2 (ja) * | 2008-07-16 | 2013-03-06 | 住友化学株式会社 | 材料特性の推定方法及び材料特性の推定プログラム |
| JP5234671B2 (ja) * | 2010-05-19 | 2013-07-10 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
| JP5653468B2 (ja) * | 2013-03-22 | 2015-01-14 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 画像表示装置 |
| CN103400023B (zh) * | 2013-06-28 | 2016-11-02 | 华北水利水电大学 | 软组织形变仿真方法 |
| JP7321951B2 (ja) * | 2019-03-01 | 2023-08-07 | ポーラ化成工業株式会社 | 肌解析方法、肌解析システム及び肌解析プログラム |
-
1993
- 1993-10-01 JP JP24712093A patent/JP3302122B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JPH07100126A (ja) | 1995-04-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Evans et al. | Measuring the mechanical properties of human skin in vivo using digital image correlation and finite element modelling | |
| Nikishkov et al. | An equivalent domain integral method for computing crack-tip integral parameters in non-elastic, thermo-mechanical fracture | |
| Villard et al. | Simulation of lung behaviour with finite elements: Influence of bio-mechanical parameters | |
| Grosse et al. | Techniques for modeling muscle‐induced forces in finite element models of skeletal structures | |
| Dong et al. | Stresses and stress intensities at notches:‘anomalous crack growth’revisited | |
| JP3302122B2 (ja) | 表情しわの模擬方法および装置 | |
| Van Schepdael et al. | Analytical determination of stress patterns in the periodontal ligament during orthodontic tooth movement | |
| Kao et al. | A microstructurally driven model for pulmonary artery tissue | |
| Lim et al. | New extensometer to measure in vivo uniaxial mechanical properties of human skin | |
| Jamli et al. | The sustainability of neural network applications within finite element analysis in sheet metal forming: A review | |
| JP3897477B2 (ja) | 応力−ひずみ関係シミュレート方法およびスプリングバック量予測方法 | |
| Liberson et al. | A physics based approach to the pulse wave velocity prediction in compliant arterial segments | |
| Zhang et al. | An optimized mass-spring model with shape restoration ability based on volume conservation | |
| Balokas et al. | Data-driven inverse uncertainty quantification in the transverse tensile response of carbon fiber reinforced composites | |
| Liu et al. | Real-time simulation of virtual palpation system | |
| Evans | On the implementation of a wrinkling, hyperelastic membrane model for skin and other materials | |
| Żur et al. | Large deformation of hyperelastic modified Timoshenko–Ehrenfest beams under different types of loads | |
| Jiang et al. | Extracting material parameters of silicone elastomers under biaxial tensile tests using virtual fields method and investigating the effect of missing deformation data close to specimen edges on parameter identification | |
| Wu et al. | Multiaxial and time dependent behavior of a filled rubber | |
| Wu et al. | Skin aging estimation by facial simulation | |
| Gérard et al. | Indentation for estimating the human tongue soft tissues constitutive law: application to a 3D biomechanical model | |
| He et al. | Explicit consideration of fiber recruitment in vascular constitutive formulation using beta functions | |
| Mazza et al. | Simulation of the aging face | |
| Marek et al. | Distortional hardening cyclic plasticity: Experiments and modeling | |
| JP3344835B2 (ja) | 永久表情しわの模擬方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20020319 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090426 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100426 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100426 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110426 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120426 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130426 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140426 Year of fee payment: 12 |
|
| EXPY | Cancellation because of completion of term |