【０００１】
【発明の属する技術分野】
本発明は、対象物の表面、特に、生体皮膚表面の力学特性、特に硬さを評価するための評価値の算出に関する。
【０００２】
【従来の技術】
特許第３２５７５６３号には、生体の表面に一定の振幅で３０Ｈｚ以上５００Ｈｚ以下の一定の周波数で大きさが変化する力を印加して生体表面を振動させたときの表面の加速度を検出し、検出した加速度の振幅に比例する量を算出し、その逆数を生体表面の硬さの評価値とする生体表面の硬さ測定装置が開示されている。
【０００３】
また測定原理、測定器の構成は同じだが、振動発生手段の形状がベンダー型（たわみ型）であり、片持ち梁に類似した振動板の上下に圧電素子を張りつけ、一方の圧電素子を加振用、他方を加速度検出用とするとともに、振動板の支持点を２支持点とし、支持点間の距離を調整することで加振周波数を任意に設定できる歯の動揺センサ（岡 久雄、清水義和、更谷啓治、川添タカヨシ、電気情報通信学会論文誌、Ｊ８０−Ａ、１７８８−１７９３（１９９７））、インプラント動揺測定システム（岡 久雄、小野浩一、松谷義雄、龍田光弘、更谷啓治、川添タカヨシ、張 銘、施生根、王 忠義、信学技報、ＭＢＥ９９−１５７、７３−７８（２０００−０３））も提案されている。
【０００４】
【特許文献１】
特許第３２５７５６３号
【非特許文献１】
岡 久雄、清水義和、更谷啓治、川添タカヨシ、電気情報通信学会論文誌、Ｊ８０−Ａ、１７８８−１７９３（１９９７）
【非特許文献２】
岡 久雄、小野浩一、松谷義雄、龍田光弘、更谷啓治、川添タカヨシ、張 銘、施生根、王 忠義、信学技報、ＭＢＥ９９−１５７、７３−７８（２０００−０３）
【０００５】
【発明が解決しようとする課題】
前述の特許第３２５７５６３号に記載されているように、生体の各部位の中でも、例えば前頭部、胸部、肋骨上、大腿部、ふくらはぎ等の、同一被験者にあっても部位によって硬さが大きく異なる場合の硬さ評価に対しては、前述した硬さ測定装置はきわめて有用である。しかしながら、加齢や肉付きによって変化する生体皮膚表面の“硬さ〜柔らかさ”を評価する場合や、肌の乾燥などで生じた生体皮膚表面の微妙な“硬さ〜柔らかさ”の違いを評価する場合、あるいは頬のように被験者間の生体皮膚表面の硬さの違いがさほど大きくないケースでは、生体皮膚表面に振動子を軽く接触させながら測定を行わないと硬さの違いが検出しにくい。例えば、頬の皮膚表面の硬さの違いを精度良く検出するためには接触荷重を１０ｇｆ位とする必要がある。しかしながら１０ｇｆという“点”としての接触荷重で振動子を皮膚表面に接触・維持して計測を行うことはきわめて困難である。そこでデータ採取を可とする接触荷重は、例えば１０±１ｇｆとしている。データ採取を可とする接触荷重が１０±１ｇｆであることは、接触荷重に対して１０％の許容範囲をもうけていることであり、これに起因する採取データのバラツキが発生する。このことは、単にバラツキのあるデータが得られるという弊害のみならず、接触荷重が設定許容範囲内にあれば無条件で硬さ〜柔らかさの評価値としての生体機械インピーダンスのデータを採取し、その算術平均を算出して設定許容範囲の中央値で代表される接触荷重に対する硬さ〜柔らかさの評価値としているという作業を実施していることになり、設定した“点”としての接触荷重に対応する評価値が得られているかどうかの保証が十分でない。
【０００６】
また、実際問題として１０ｇｆ±１ｇｆに接触荷重を維持することは容易ではなく、測定者にかかる負担は大きい。結果として接触荷重の維持の難しさに基づく測定誤差の発生や、測定に要する時間が長くなるといった弊害が生じる。
【０００７】
さらに、前述したいずれの装置も、電気信号（駆動電圧）が機械系（振動板ないしは振動子）の振動に円滑に変換されているか否かのチェック、あるいは電気系から機械系への変換状態の変化（測定システムの定常性異常の発生）などをチェックすること、即ちシステム異常や人為的な測定ミスの発生の有無を検知し、正常と思われるデータのみを採取することができない。
【０００８】
本発明の第１の目的は、上記したような曖昧さを持つ接触荷重に対する硬さ〜柔らかさの評価値ではなく、設定した“点”としての接触荷重に対する硬さ〜柔らかさの評価値の計測を可とすることである。
【０００９】
本発明の第２の目的は、接触荷重の許容範囲をより広く（例えば、１０ｇｆ±２ｇｆ）設定できるようにして測定そのものの実施を容易にし、かつ生体皮膚の硬さ〜柔らかさをより精度高く計測することを可とすることである。
【００１０】
本発明の第３の目的は上記したシステム異常や人為的な測定ミスに対するチェック機能を付加し、自動的に正常なデータのみを採取することである。
【００１１】
【課題を解決するための手段】
本発明の硬さ測定装置は、対象物の表面に一定の振幅および所定の周波数で大きさが変化する力を印加したときの表面の加速度の代表値および接触荷重を測定する手段と、該測定手段による測定を複数回繰り返して得られる測定値から接触荷重と加速度の代表値の関係を表わす関係式を算出する手段と、該関係式を使って所定の接触荷重における加速度の代表値を算出することによって対象物の表面の硬さを評価する手段とを具備することを特徴とする。
【００１２】
前記測定手段は、印加された力の瞬時値に相当する第１の測定値とその各々に対する加速度の瞬時値に相当する第２の測定値との関係の妥当性を確認する手段を含む。
【００１３】
この妥当性確認手段は、印加された力の第１および第２の周期のそれぞれについて、複数の前記第１の測定値と、その各々に対する複数の前記第２の測定値とから得られる、２つの回帰直線がいずれも有意であり、実質的に原点を通り、かつ、両者の傾きが実質的に等しいとき、第１の測定値と第２の測定値の関係が妥当であるとする。
【００１４】
【発明の実施の形態】
特許第３２５７５６３号に記載されているように、生体表面の硬さ〜柔らかさは、生体機械インピーダンスを測定することによって評価することができる。この生体機械インピーダンスは角周波数ωの関数として複素数Ｚ（ω）として表され、生体表面に印加した力Ｆ（ω）とその応答としての生体表面の速度Ｖ（ω）から
Ｚ（ω）＝Ｆ（ω）／Ｖ（ω）
と定義され、加速度Ａ（ω）＝ｊωＶ（ω）を使うと
Ｚ（ω）＝ｊωＦ（ω）／Ａ（ω） （１）
と表される。
【００１５】
式（１）から、特定の周波数ω０ における複素機械インピーダンスＺ（ω０ ）は、
Ｚ（ω０ ）＝ｊω０ Ｆ（ω０ ）／Ａ（ω０ ）
となる。両辺の絶対値をとると
｜Ｚ（ω０ ）｜＝ω０ ｜Ｆ（ω０ ）｜／｜Ａ（ω０ ）｜
となり、ω０ および｜Ｆ（ω０ ）｜が一定であれば
｜Ｚ（ω０ ）｜＝ｋ／｜Ａ（ω０ ）｜ ； ｋ＝比例定数 （２）
式（２）は、生体表面に振幅および周波数が一定の振動をあたえた時の加速度の絶対値または実効値などの代表値を測定すれば、その周波数における生体機械インピーダンスに逆比例する量が得られること、その逆数をとれば生体機械インピーダンスに比例する量が得られることを意味している。但し、式（２）は、ｋという比例定数を含んでいるので生体機械インピーダンスそのものを表すものではない。
【００１６】
上記した生体機械インピーダンスに比例する量、すなわち生体表面の“硬さ〜柔らかさ”の指標を記号ＳＨ（Ｓｋｉｎ Ｈａｒｄｎｅｓｓ）で表すことにし
ＳＨ＝１／Ａｅ （３）
と定義する。ここでＡｅ は、測定された加速度の実効値である。加速度の実効値が大きいということは、生体表面の駆動点が動きやすいことであり、柔らかいことを意味する。
【００１７】
特許第３２５７５６３号では、式（３）のＳＨを生体表面の硬さを表すパラメータとしているが、Ａｅ を生体表面の硬さのパラメータとすることもできる。
【００１８】
また、触診で硬さ〜柔らかさを評価する場合、３段階評価（＋１，０，−１、すなわち硬い、ふつう、柔らかい）、あるいは５段階評価（＋２，＋１，０，−１，−２、すなわち硬い、やや硬い、ふつう、やや柔らかい、柔らかい）とするケースが多い。
【００１９】
式（３）で得られたＳＨの値を、例えば５段階評価値（記号ＳＨ’で表記）に変換するとしたら、ＳＨを式（４）に代入すればよい。
【００２０】
ＳＨ’＝４×（ｓｈ−ＳＨｍａｘ ）／（ ＳＨｍｉｎ−ＳＨｍａｘ ）−２ （４）
ここでＳＨｍａｘ ，ＳＨｍｉｎ は、多数のボランティアの生体表面のＳＨの測定値から得られた最大値、最小値である。またｓｈは、当該ボランティアのＳＨの実測値である。
【００２１】
図１は本発明の硬さ測定装置の一実施例に係る、皮膚の硬さ測定システムの構成を示す。
【００２２】
本システムは、プローブ１０、皮膚硬さ測定器１２、１２ビットＡ／Ｄコンバータ１４およびパーソナルコンピュータ１６から構成されている。
【００２３】
皮膚硬さ測定器１２にはＡ／Ｄコンバータ内蔵のＣＰＵ１８、一定振幅および所定の周波数の駆動電圧を発生してプローブ１０の圧電素子（後述）へ与えるための正弦波発生器２０、プローブ１０の加速度センサ（後述）からの信号から所望の帯域の信号のみを通過させてＡ／Ｄコンバータ１４へ供給するためのバンドバスフィルタ２２、及びプローブ１０の接触荷重センサ（後述）からの信号をサンプル／ホールドしてＣＰＵ１８内蔵のＡ／Ｄコンバータへ供給するサンプル／ホールド回路２４が設けられている。
【００２４】
図２はプローブ１０の詳細を示す。シム材２６の両面に駆動用圧電セラミックス２８および検出用圧電セラミックス３０が貼り付けられ、その先端に振動子３６が取り付けられる。
【００２５】
駆動用圧電素子２８に正弦波を印加して振動させた状態で振動子３６を測定対象の皮膚表面に接触させることにより、検出用圧電素子３０に加速度に比例した電圧が発生する。接触荷重は歪みゲージ３８により検出される。センサ全体は支点３２および３４の２点において支持されるが、１点において支持するようにしても良い。
【００２６】
本願発明の第１の目的である、設定した“点”としての接触荷重に対する硬さ（柔らかさ）の評価値を得るため、および第２の目的である、接触荷重の許容範囲を現実的なものとするため、図１の皮膚硬さ測定システムにおいては、プローブ１０の接触子３６を皮膚表面に軽く接触させたときの接触荷重をチェックし、それが予め設定された許容範囲内にあることを条件として接触荷重と加速度の実効値を複数回測定し、それらから図３に示すような回帰直線を算出し、これを用いて、特定の“点”としての接触荷重における加速度の実効値を算出する。図３に示した例では、接触荷重の許容範囲を１０ｇｆ±３０％として得られた複数の測定値から回帰式を求め、それを用いて接触荷重１０ｇｆにおける加速度の推定実効値６０３．１６が得られている。
【００２７】
また本発明の第３の目的である、システム異常や人為的測定ミスの発生の有無を検知するため、時々刻々得られる駆動電圧値と加速度の検出電圧値との関係に対して回帰分析を適用して妥当性をチェックする。すなわち、生体表面に振動子を接触させた場合に得られる加速度に対応する出力電圧波形の位相は、粘弾性体である生体を対象にした場合、図４に示すように、駆動電圧波形の位相に対して位相遅れが生じる。したがって駆動電圧ｘを横軸に、加速度に対応する出力電圧ｙを縦軸にプロットすると図５に示すような楕円の形をしたリサージュ図形が描ける。この場合、横軸にプロットされる駆動電圧ｘは確定変数であり、かつその振幅の大きさは一定であることから楕円の形を形成するリサージュ図形の慣性モーメントが最大となる軸の勾配の大きさは、加速度の大きさ、すなわち生体皮膚表面の硬さ〜柔らかさに関する情報を持つ。そしてシステム異常や測定ミスが発生していなければ歪みのない綺麗なリサージュ図形が得られるはずである。そこでリサージュ図形に歪みがみられるか否かをみるために図６に示すように駆動電圧ｘと加速度に対応する出力電圧ｙに関する回帰直線４０を算出する。なお、駆動電圧ｘとは正弦波発生器２０（図１）の出力をＡ／Ｄ変換器１４でディジタルに変換したものであり、出力電圧ｙとは、バンドパスフィルタ２２の出力をＡ／Ｄ変換器１４でディジタル変換したものである。
【００２８】
この場合、回帰直線に有意性がみられないとしたら、駆動電圧ｘの変化と加速度に対応する出力電圧ｙの変化に対応関係がみられないことを意味し、何らかのシステム異常や人為的な測定ミスが発生していることになる。
【００２９】
駆動電圧ｘの波形と加速度に対応する出力電圧ｙの波形の周期が同じなら、リサージュ図形のプロットから得られる回帰直線は、原点を通る。システム異常が発生しているか、あるいは人為的な測定ミスが発生している場合は、原点を通らない。原点を通っていない場合は、上記と同様、異常とする。
【００３０】
連続する２周期について回帰直線を算出した場合、同一の生体表面を測定対象にして時間的に連続計測しているわけだから、前半１周期で得られたｘとｙを用いて描いたリサージュ図形に見られる回帰直線の勾配と、後半１周期で得られたｘとｙを用いて描いたリサージュ図形に見られる回帰直線の勾配は類似しているはずである。即ち、平行性は極めて高いと予測される。この勾配の平行性の有意性の検定を行うことにより、前半で得られた皮膚の硬さの測定値と後半で得られた皮膚の硬さの測定値の同一性を判断する。
【００３１】
図７は、上記の手法に基づくパーソナルコンピュータ１６（図１）およびＣＰＵ１８における処理のフローチャートであり、図８は図７中の“妥当性チェック”の処理の詳細を示すフローチャートであり、図９は測定条件設定のための画面の一例を示す。
【００３２】
まず図９に示す画面で測定条件の設定を行なう。図９中、“サンプリング数”とは、駆動電圧波形１周期におけるサンプリング数、“測定回数”とは駆動電圧波形２周期を１回としたときの測定回数である。
【００３３】
図７において、まずパーソナルコンピュータ１６において図９に示す画面による設定が終わって、皮膚硬さ測定器１２（図１）の副電源スイッチ１９がオンにされたら（ステップ１０００）、ＣＰＵ１８がプローブ１０へ−３Ｖの電源を供給して歪みゲージ３８（図２）を通電し、正弦波発生器２０（図１）の電源をオンにする。
【００３４】
このとき、プローブ１０を測定対象に限りなく近付け測定時の向きで保持する。通電後５秒待って（ステップ１００２）、図８に詳細を示す妥当性チェックを行ない（ステップ１００４）、測定値が妥当であればそのときの接触荷重Ｘおよび加速度の実効値Ｙを算出して格納する（ステップ１００６）。接触荷重Ｘは妥当性チェックのための加速度データの採取に同期して歪みゲージ３８からサンプル／ホールド回路２４を経て採取したデータの平均値であり、加速度の実効値Ｙはそのときに得られる加速度の瞬時値の自乗平均の平方根をとることにより算出される。
【００３５】
ステップ１００４および１００６を“測定回数”で設定された回数だけ実行したら、それら複数のＸおよびＹの値に対して回帰分析を適用して回帰直線の式（回帰式）を算出する（ステップ１００８）。算出された回帰式を用いて図９の“中心荷重”において設定された荷重における加速度の実効値を算出して式（３）によりパラメータＳＨを算出する（ステップ１０１０）。式（４）によりＳＨをＳＨ’に変換して（ステップ１０１２）、結果を表示した後（ステップ１０１４）、副電源スイッチ１９をオフにする。
【００３６】
上記のステップ１００４の妥当性チェックにおいては、図８に示すように、まず、接触荷重が許容範囲内になるのを待ち（ステップ１１００）、連続２周期にわたって駆動電圧ｘおよび検出電圧ｙを採取し（ステップ１１０２）、前半１周期について回帰直線の算出を試み（ステップ１１０４）、回帰が有意でなければ採取したデータを廃棄してステップ１１００の処理へ戻る。回帰が有意であれば、回帰直線が原点を通るか、すなわち直線の切片のｙ座標の絶対値が予め定められた値以下であるかをチェックする（ステップ１１０６）。原点を通らなければステップ１１００へ戻り、原点を通れば後半１周期についても同様な処理を行なう（ステップ１１０８〜１１１０）。
【００３７】
それぞれの回帰直線に異常がなければ、両者が平行かどうか、すなわち両者の傾きの差が所定値以下であるかを調べ（ステップ１１１２）、平行でなければステップ１１１０へ戻る。平行であれば加速度検出電圧のサンプリングに同期して採取した接触荷重の平均値を算出し（ステップ１１１４）、接触荷重平均値が許容範囲内になければステップ１１１０に戻り、許容範囲内であれば妥当性チェックの処理を終了する。
【００３８】
特に歪みゲージは通電開始とともに熱が発生し、振動板への熱の伝播が続く。それに伴い測定値の継続的な増大が観測される。また、ゲージの発熱による障害の発生と比較してその影響は小さいとおもわれるが、振動板の振動に起因する発熱、振動板の金属疲労の蓄積も考えられ、測定をしていない時の振動板の振動停止は良策と考えられる。これらのことから上記のように測定時のみ、正弦波発生器２０と歪みゲージ３８の電源をＯＮとすることが望ましい。また、副電源をＯＮにして正弦波発生器に電源を供給すると所定の周波数の正弦波で振動板および振動子は駆動されるが、１次の固有モードで安定した振動状態を呈するまですなわち歪みゲージの出力電圧が一定になるまで約５秒必要とする。そのため、上記のように副電源ＯＮから５秒後に測定を開始することが望ましい。
【００３９】
本実施例の場合、本発明の第１の目的である設定した“点”としての接触荷重（例えば１０ｇｆ）に対する硬さ〜柔らかさの評価値が得られていることは、まず接触荷重Ｘにおける加速度の実効値Ｙを推定する回帰式をもとめ、得られた回帰式を用いて、所定の“点”としての接触荷重Ｘにおける加速度の実効値Ｙを算出していることから自明である。
【００４０】
本発明の第２の目的は、接触荷重の許容範囲をより広く（例えば、１０ｇｆ±２ｇｆ）設定できるようにして測定そのものの実施を容易にし、かつ生体皮膚表面の硬さ〜柔らかさをより精度高く計測することを可とすることである。
【００４１】
上記した目的が達成できたかどうかの確認を硬さの異なる２種類のシリコン材料（大きさ、直径８５ｍｍ、長さ５０ｍｍのＲＴＶ（室温硬化型）ゴム、アスカー硬度計（Ｆ型）の測定値、硬い材料＝２７、柔らかい材料＝１５）を測定試料として実施したところ表１に示すような結果が得られた。
【００４２】
【表１】
【００４３】
Ａ２ の接触荷重の許容範囲（本願実施例）は、Ａ１ の接触荷重の許容範囲の２倍になっている。単純に考えれば、測定は容易に実行できるが、データのバラツキは大きくなると予測される。しかしながら本発明の第２の目的は、接触荷重の許容範囲を広くしても、硬さ〜柔らかさをより高精度に測定できるようにすることである。加速度に対応する出力電圧を１２ビットに量子化し実効値をもとめ、ＳＮ比を算出した。得られたＡ１ ，Ａ２ の測定器のＳＮ比は、下記の如くであった。
【００４４】
結果として本願実施例の方が６．２デシベル、真数で４．２倍の利得がある。このことは、本願実施例の方法・装置で１回測定した時の測定精度は、特許第３２５７５６３号に記載した装置で４．２回測定し、その平均値を求めた時の測定精度と同じであることを意味している。すなわちこのケースでは接触荷重の許容範囲を２倍の±２ｇｆとし、接触荷重の維持に関する測定者の負担を軽くするとともに、測定の実施を容易にし、合わせて測定精度も向上させた。
【００４５】
また本発明の第３の目的であるシステム異常や測定ミスの有無を検知し、正常と思われるデータのみを自動的に採取できるシステムの構築に関しては、振動板および振動子を駆動するための駆動電圧と加速度に対応する出力電圧の整合性を統計的手法を導入し検定すること、で実現した。
【００４６】
【発明の効果】
（１）設定した“点”としての接触荷重に対する生体皮膚の硬さ〜柔らかさの評価値の計測が可能である。
【００４７】
（２）測定時の接触荷重の許容範囲を広くし、測定のしやすさ、測定者の負担の軽減をはかりつつ、同時に測定精度の向上を実現している。
【００４８】
（３）電気信号（駆動電圧）が機械系（振動板ないしは振動子）の振動に円滑に変換されているか否かのチェック、あるいは電気系から機械系への変換状態の変化（測定システムの定常性異常の発生）などをチェック、すなわちシステム異常や人為的な測定ミスの発生の有無を検知し、正常と思われるデータのみ採取できる。
【図面の簡単な説明】
【図１】本発明の一実施形態に係る皮膚の硬さ測定システムの構成を示すブロック図である。
【図２】図１のプローブ１０の詳細を示す図である。
【図３】所定の“点”としての接触荷重における加速度の算出のための回帰直線を説明するグラフである。
【図４】駆動電圧の波形と加速度検出電圧の波形の関係を示すグラフである。
【図５】駆動電圧を横軸、加速度検出電圧をたて軸とするリサージュ図形を示すグラフである。
【図６】リサージュ図形に対する回帰直線を説明するグラフである。
【図７】図１のパーソナルコンピュータ１６における処理のフローチャートである。
【図８】図７の“妥当性チェック”の詳細を示すフローチャートである。
【図９】パーソナルコンピュータ１６の設定画面を示す図である。[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to calculation of an evaluation value for evaluating mechanical properties, particularly hardness, of a surface of an object, in particular, a biological skin surface.
[0002]
[Prior art]
Japanese Patent No. 3257563 discloses detecting and detecting the acceleration of a surface of a living body when the body is vibrated by applying a force whose magnitude changes at a constant frequency of 30 Hz or more and 500 Hz or less at a constant amplitude to the surface of the living body. There is disclosed an apparatus for measuring the hardness of a living body surface, which calculates an amount proportional to the amplitude of the acceleration thus obtained and uses the reciprocal thereof as an evaluation value of the hardness of the living body surface.
[0003]
Although the measurement principle and the configuration of the measuring device are the same, the shape of the vibration generating means is a bender type (flexible type), and piezoelectric elements are attached to the top and bottom of a diaphragm similar to a cantilever, and one piezoelectric element is vibrated. And the other is used for acceleration detection, the vibration plate has two support points, and the excitation frequency can be set arbitrarily by adjusting the distance between the support points (Hisao Oka, Yoshikazu Shimizu , Keiji Saratani, Takayoshi Kawazoe, Transactions of the Institute of Electronics, Information and Communication Engineers, J80-A, 1788-1793 (1997)), implant sway measurement system (Hisao Oka, Koichi Ono, Yoshio Matsutani, Mitsuhiro Tatsuta, Keiji Saratani, Takayoshi Kawazoe , Zhang Ming, Shigene, Tadayoshi Wang, IEICE Technical Report, MBE 99-157, 73-78 (2000-03)).
[0004]
[Patent Document 1]
Patent No. 3257563 [Non-Patent Document 1]
Hisao Oka, Yoshikazu Shimizu, Keiji Saraya, Takayoshi Kawazoe, IEICE Transactions, J80-A, 1788-1793 (1997)
[Non-patent document 2]
Hisao Oka, Koichi Ono, Yoshio Matsutani, Mitsuhiro Tatsuta, Keiji Saratani, Takayoshi Kawazoe, Zhang Mei, Shione, Tadayoshi Wang, IEICE Technical Report, MBE 99-157, 73-78 (2000-03)
[0005]
[Problems to be solved by the invention]
As described in the aforementioned Japanese Patent No. 3257563, among various parts of the living body, for example, the hardness of some parts of the same subject, such as the frontal region, the chest, on the ribs, the thighs, and the calves, is high. The hardness measuring device described above is very useful for evaluating the hardness in the case of a great difference. However, when evaluating the "hardness to softness" of the living skin surface that changes with aging and flesh, and evaluating the subtle "hardness to softness" of the living skin surface caused by skin dryness In cases where the difference in hardness of the skin surface between subjects is not so large, such as in the case of cheeks, it is difficult to detect the difference in hardness unless measurement is performed while lightly contacting the vibrator with the skin surface of the body. . For example, in order to accurately detect the difference in hardness of the cheek skin surface, the contact load needs to be about 10 gf. However, it is extremely difficult to measure by contacting and maintaining the vibrator on the skin surface with a contact load of "point" of 10 gf. Therefore, the contact load at which data can be collected is, for example, 10 ± 1 gf. The fact that the contact load at which data collection is possible is 10 ± 1 gf means that a 10% tolerance is provided for the contact load, which causes variations in the collected data. This means not only the harmful effect that data with variation is obtained, but also the data of biomechanical impedance as the evaluation value of hardness to softness is unconditionally collected if the contact load is within the set allowable range, The work of calculating the arithmetic average and setting the evaluation value of hardness to softness with respect to the contact load represented by the median of the set permissible range is being carried out, and the contact load as a set “point” There is not enough guarantee that the evaluation value corresponding to is obtained.
[0006]
Further, as a practical matter, it is not easy to maintain the contact load at 10 gf ± 1 gf, and the burden on the measurer is large. As a result, there are adverse effects such as occurrence of a measurement error based on the difficulty in maintaining the contact load and an increase in the time required for the measurement.
[0007]
Further, any of the above-described devices checks whether or not the electric signal (drive voltage) is smoothly converted to the vibration of the mechanical system (diaphragm or vibrator), or checks the state of conversion from the electric system to the mechanical system. It is not possible to check for changes (occurrence of abnormalities in the measurement system), that is, to detect the occurrence of system abnormalities or artificial measurement errors, and to collect only data considered to be normal.
[0008]
A first object of the present invention is not to evaluate the hardness to softness for a contact load having the above-mentioned ambiguity, but to evaluate the hardness to softness for a contact load as a set “point”. This is to enable measurement.
[0009]
A second object of the present invention is to make it possible to set the allowable range of the contact load wider (for example, 10 gf ± 2 gf) to facilitate the measurement itself, and to more precisely measure the hardness to softness of the living skin. It is possible to measure.
[0010]
A third object of the present invention is to add a check function for the above-mentioned system abnormality or artificial measurement error, and to automatically collect only normal data.
[0011]
[Means for Solving the Problems]
The hardness measuring device of the present invention is a means for measuring a representative value of a surface acceleration and a contact load when a force that changes its magnitude at a constant amplitude and a predetermined frequency is applied to the surface of an object; Means for calculating a relational expression representing a relationship between a contact load and a representative value of acceleration from measurement values obtained by repeating the measurement by the means a plurality of times, and calculating a representative value of the acceleration at a predetermined contact load using the relational expression Means for evaluating the hardness of the surface of the object.
[0012]
The measuring means includes means for confirming the validity of the relationship between the first measured value corresponding to the instantaneous value of the applied force and the second measured value corresponding to the instantaneous value of the acceleration for each of them.
[0013]
The plausibility check means comprises, for each of the first and second periods of the applied force, a plurality of the first measurements and a plurality of the second measurements for each of the first measurements. When all the two regression lines are significant, substantially pass through the origin, and have substantially the same slope, it is assumed that the relationship between the first measurement value and the second measurement value is appropriate.
[0014]
BEST MODE FOR CARRYING OUT THE INVENTION
As described in Japanese Patent No. 3257563, hardness to softness of a living body surface can be evaluated by measuring biomechanical impedance. This biomechanical impedance is expressed as a complex number Z (ω) as a function of the angular frequency ω, and Z (ω) = F from the force F (ω) applied to the surface of the living body and the velocity V (ω) of the living body surface in response thereto. (Ω) / V (ω)
When the acceleration A (ω) = jωV (ω) is used, Z (ω) = jωF (ω) / A (ω) (1)
It is expressed as
[0015]
From equation (1), the complex mechanical impedance Z (ω 0 ) at a specific frequency ω 0 is
Z (ω 0 ) = jω 0 F (ω 0 ) / A (ω 0 )
It becomes. When the absolute values of both sides are taken, | Z (ω 0 ) | = ω 0 | F (ω 0 ) | / | A (ω 0 ) |
And if ω 0 and | F (ω 0 ) | are constant, | Z (ω 0 ) | = k / | A (ω 0 ) |; k = proportional constant (2)
Equation (2) is obtained by measuring a representative value such as an absolute value or an effective value of acceleration when a constant amplitude and frequency are applied to the surface of the living body, and obtains an amount inversely proportional to the biomechanical impedance at that frequency. In other words, taking the reciprocal thereof means that an amount proportional to the biomechanical impedance can be obtained. However, equation (2) does not represent the biomechanical impedance itself because it includes a proportional constant of k.
[0016]
The quantity proportional to the biomechanical impedance described above, that is, the index of “hardness to softness” of the living body surface is represented by the symbol SH (Skin Hardness), and SH = 1 / A e (3)
Is defined. Here, Ae is the effective value of the measured acceleration. The fact that the effective value of the acceleration is large means that the driving point on the living body surface is easy to move and is soft.
[0017]
In Japanese Patent No. 3257563, SH in Expression (3) is used as a parameter representing the hardness of the living body surface, but Ae may be used as a parameter of the hardness of the living body surface.
[0018]
When evaluating hardness to softness by palpation, a three-level evaluation (+1, 0, -1, that is, hard, ordinary, and soft) or a five-level evaluation (+2, +1, 0, -1, -2, That is, there are many cases in which it is hard, slightly hard, ordinary, slightly soft, or soft.
[0019]
If the value of SH obtained by equation (3) is converted into, for example, a five-level evaluation value (denoted by the symbol SH ′), SH may be substituted into equation (4).
[0020]
SH ′ = 4 × (sh−SH max ) / (SH min −SH max ) −2 (4)
Here, SH max and SH min are the maximum value and the minimum value obtained from the measured values of SH on the living body surfaces of many volunteers. Also, sh is the actual measured value of SH of the volunteer.
[0021]
FIG. 1 shows a configuration of a skin hardness measuring system according to an embodiment of the hardness measuring apparatus of the present invention.
[0022]
This system includes a probe 10, a skin hardness measuring device 12, a 12-bit A / D converter 14, and a personal computer 16.
[0023]
The skin hardness measuring device 12 has a CPU 18 with a built-in A / D converter, a sine wave generator 20 for generating a driving voltage having a constant amplitude and a predetermined frequency and applying the driving voltage to a piezoelectric element (described later) of the probe 10. A signal from a contact load sensor (to be described later) of the probe 10 and a band-pass filter 22 for passing only a signal of a desired band from the signal from the acceleration sensor (to be described later) to be supplied to the A / D converter 14 are sampled / A sample / hold circuit 24 for holding and supplying the data to an A / D converter built in the CPU 18 is provided.
[0024]
FIG. 2 shows details of the probe 10. The driving piezoelectric ceramics 28 and the detecting piezoelectric ceramics 30 are attached to both surfaces of the shim material 26, and the vibrator 36 is attached to the tip thereof.
[0025]
By bringing the vibrator 36 into contact with the surface of the skin to be measured in a state where a sine wave is applied to the driving piezoelectric element 28 and vibrated, a voltage proportional to the acceleration is generated in the detecting piezoelectric element 30. The contact load is detected by the strain gauge 38. The whole sensor is supported at two points of fulcrums 32 and 34, but may be supported at one point.
[0026]
The first object of the present invention is to obtain an evaluation value of hardness (softness) with respect to a set contact load as a “point”, and the second object is to set a realistic allowable range of a contact load. In the skin hardness measurement system shown in FIG. 1, the contact load when the contact 36 of the probe 10 is lightly brought into contact with the skin surface is checked, and the contact load is within a preset allowable range. The contact load and the effective value of the acceleration are measured a plurality of times under the condition, and a regression line as shown in FIG. 3 is calculated from them, and the effective value of the acceleration at the contact load as a specific “point” is calculated using the regression line. calculate. In the example shown in FIG. 3, a regression equation is obtained from a plurality of measurement values obtained by setting the allowable range of the contact load to 10 gf ± 30%, and the estimated effective value 603.16 of the acceleration at the contact load of 10 gf is obtained using the regression equation. Have been.
[0027]
In addition, in order to detect the occurrence of a system abnormality or an artificial measurement error, which is a third object of the present invention, regression analysis is applied to the relationship between the drive voltage value obtained every moment and the detected voltage value of acceleration. To check the validity. That is, when a living body which is a viscoelastic body is targeted, the phase of the output voltage waveform corresponding to the acceleration obtained when the vibrator is brought into contact with the surface of the living body is, as shown in FIG. Causes a phase delay. Therefore, when the drive voltage x is plotted on the horizontal axis and the output voltage y corresponding to the acceleration is plotted on the vertical axis, an elliptical Lissajous figure as shown in FIG. 5 can be drawn. In this case, the drive voltage x plotted on the horizontal axis is a deterministic variable, and the magnitude of the amplitude is constant, so that the magnitude of the gradient of the axis at which the moment of inertia of the Lissajous figure that forms the ellipse is maximized is large. The depth has information on the magnitude of the acceleration, that is, the hardness to the softness of the living skin surface. Unless a system error or measurement error occurs, a beautiful Lissajous figure without distortion should be obtained. Therefore, in order to determine whether or not the Lissajous figure is distorted, a regression line 40 relating to the drive voltage x and the output voltage y corresponding to the acceleration is calculated as shown in FIG. The drive voltage x is obtained by converting the output of the sine wave generator 20 (FIG. 1) into a digital signal by the A / D converter 14, and the output voltage y is obtained by converting the output of the band-pass filter 22 into an A / D signal. It is digitally converted by the converter 14.
[0028]
In this case, if the regression line has no significance, it means that there is no corresponding relationship between the change in the drive voltage x and the change in the output voltage y corresponding to the acceleration. A mistake has occurred.
[0029]
If the cycle of the waveform of the drive voltage x and the cycle of the waveform of the output voltage y corresponding to the acceleration are the same, the regression line obtained from the plot of the Lissajous figure passes through the origin. If a system error has occurred or an artificial measurement error has occurred, do not pass through the origin. If it does not pass through the origin, it is determined to be abnormal, as described above.
[0030]
When the regression line is calculated for two consecutive cycles, the same biological surface is measured continuously because of the time, so the Lissajous figure drawn using x and y obtained in the first half cycle The slope of the regression line seen should be similar to the slope of the regression line seen in the Lissajous figure drawn using x and y obtained in the last half cycle. That is, the parallelism is expected to be extremely high. By performing a test of the significance of the parallelism of the gradient, the identity of the measured value of skin hardness obtained in the first half and the measured value of skin hardness obtained in the second half is determined.
[0031]
FIG. 7 is a flowchart of processing in the personal computer 16 (FIG. 1) and the CPU 18 based on the above-described method. FIG. 8 is a flowchart showing details of the "validity check" processing in FIG. 7, and FIG. 4 shows an example of a screen for setting measurement conditions.
[0032]
First, measurement conditions are set on the screen shown in FIG. In FIG. 9, “sampling number” is the number of samplings in one cycle of the driving voltage waveform, and “measurement number” is the number of measurements in two cycles of the driving voltage waveform.
[0033]
In FIG. 7, first, when the setting on the screen shown in FIG. 9 is completed in the personal computer 16 and the sub power switch 19 of the skin hardness measuring device 12 (FIG. 1) is turned on (step 1000), the CPU 18 A power of −3 V is supplied to energize the strain gauge 38 (FIG. 2), and the power of the sine wave generator 20 (FIG. 1) is turned on.
[0034]
At this time, the probe 10 is held as close as possible to the object to be measured and in the orientation at the time of measurement. Waiting for 5 seconds after energization (step 1002), a validity check shown in detail in FIG. 8 is performed (step 1004), and if the measured values are valid, the contact load X and the effective value Y of the acceleration at that time are calculated. It is stored (step 1006). The contact load X is an average value of data collected from the strain gauge 38 via the sample / hold circuit 24 in synchronization with the collection of acceleration data for validity check, and the effective value Y of the acceleration is the acceleration obtained at that time. Is calculated by taking the square root of the root mean square of the instantaneous value of
[0035]
After executing Steps 1004 and 1006 for the number of times set in the “number of measurements”, regression analysis is applied to the plurality of X and Y values to calculate a regression line equation (regression equation) (Step 1008). . Using the calculated regression equation, the effective value of the acceleration at the load set in the “center load” of FIG. 9 is calculated, and the parameter SH is calculated by equation (3) (step 1010). SH is converted to SH ′ by equation (4) (step 1012), and after displaying the result (step 1014), the sub power switch 19 is turned off.
[0036]
In the validity check in the above step 1004, as shown in FIG. 8, first, it is waited that the contact load falls within the allowable range (step 1100), and the drive voltage x and the detection voltage y are collected over two continuous cycles. (Step 1102) Attempt to calculate a regression line for the first half cycle (Step 1104). If the regression is not significant, the collected data is discarded and the process returns to Step 1100. If the regression is significant, it is checked whether the regression line passes through the origin, that is, whether the absolute value of the y-coordinate of the intercept of the line is equal to or smaller than a predetermined value (step 1106). If it does not pass through the origin, the process returns to step 1100, and if it passes through the origin, the same processing is performed for the second half cycle (steps 1108 to 1110).
[0037]
If there is no abnormality in the respective regression lines, it is checked whether or not the two are parallel, that is, whether or not the difference between the inclinations is equal to or smaller than a predetermined value (step 1112). If not, the process returns to step 1110. If parallel, an average value of the contact loads taken in synchronization with the sampling of the acceleration detection voltage is calculated (step 1114). If the average contact load is not within the allowable range, the process returns to step 1110; End the validity check processing.
[0038]
In particular, heat is generated in the strain gauge at the start of energization, and the propagation of heat to the diaphragm continues. A continuous increase in the measured value is observed with this. In addition, the effect is considered to be small compared to the occurrence of failure due to the heat generation of the gauge.However, the heat generated by the vibration of the diaphragm and the accumulation of metal fatigue of the diaphragm are also considered, and the vibration during measurement is not considered. Stopping the vibration of the plate is considered a good measure. From these facts, it is desirable to turn on the power of the sine wave generator 20 and the strain gauge 38 only at the time of measurement as described above. When the sine wave generator is turned on with the sub power supply turned on, the diaphragm and the vibrator are driven by a sine wave of a predetermined frequency. It takes about 5 seconds until the gauge output voltage becomes constant. Therefore, it is desirable to start the measurement 5 seconds after the auxiliary power is turned on as described above.
[0039]
In the case of the present embodiment, the evaluation value of hardness to softness with respect to the contact load (for example, 10 gf) as the set “point”, which is the first object of the present invention, means that the contact load X It is obvious that the regression equation for estimating the effective value Y of the acceleration is obtained, and the effective value Y of the acceleration at the contact load X as a predetermined “point” is calculated using the obtained regression equation.
[0040]
A second object of the present invention is to make it possible to set a wider allowable range of the contact load (for example, 10 gf ± 2 gf) to facilitate the measurement itself, and to more precisely measure the hardness to the softness of the surface of the living skin. It is possible to measure high.
[0041]
Whether or not the above-mentioned object was achieved was confirmed by measuring two types of silicon materials having different hardnesses (size, diameter 85 mm, length 50 mm, RTV (room temperature curing type) rubber, Asker hardness meter (F type), When a hard material = 27 and a soft material = 15) were used as measurement samples, the results shown in Table 1 were obtained.
[0042]
[Table 1]
[0043]
Allowable range of the contact load of the A 2 (present embodiment) is twice the tolerance range of the contact load A 1. In simple terms, the measurement can be easily performed, but the variation in data is expected to be large. However, a second object of the present invention is to make it possible to measure hardness to softness with higher accuracy even when the allowable range of the contact load is widened. The output voltage corresponding to the acceleration was quantized to 12 bits, the effective value was obtained, and the SN ratio was calculated. The SN ratios of the obtained A 1 and A 2 measuring instruments were as follows.
[0044]
As a result, the embodiment of the present invention has a gain of 6.2 decibels and a gain of 4.2 times as an antilog. This means that the measurement accuracy when measuring once with the method and apparatus of the embodiment of the present application is the same as the measurement accuracy when measuring 4.2 times with the apparatus described in Japanese Patent No. 3257563 and calculating the average value. It means that That is, in this case, the allowable range of the contact load was doubled to ± 2 gf, so that the burden on the operator for maintaining the contact load was reduced, the measurement was easily performed, and the measurement accuracy was also improved.
[0045]
The third object of the present invention is to construct a system capable of automatically detecting only data considered to be normal by detecting the presence or absence of a system abnormality or a measurement error, and for driving a diaphragm and a vibrator. The consistency of the output voltage corresponding to the voltage and the acceleration was realized by introducing a statistical method and testing.
[0046]
【The invention's effect】
(1) It is possible to measure the evaluation value of the hardness to the softness of the living skin with respect to the set contact load as a “point”.
[0047]
(2) The allowable range of the contact load at the time of measurement is widened, and the measurement accuracy is improved while the ease of measurement and the burden on the measurer are reduced.
[0048]
(3) Check whether the electric signal (drive voltage) is smoothly converted to the vibration of the mechanical system (diaphragm or vibrator), or change the conversion state from the electric system to the mechanical system (steady state of the measurement system). (I.e., occurrence of sexual abnormalities), etc., that is, whether or not a system abnormality or an artificial measurement error has occurred is detected, and only data considered to be normal can be collected.
[Brief description of the drawings]
FIG. 1 is a block diagram showing a configuration of a skin hardness measurement system according to an embodiment of the present invention.
FIG. 2 is a diagram showing details of a probe 10 of FIG. 1;
FIG. 3 is a graph illustrating a regression line for calculating an acceleration at a contact load as a predetermined “point”;
FIG. 4 is a graph showing a relationship between a waveform of a driving voltage and a waveform of an acceleration detection voltage.
FIG. 5 is a graph showing a Lissajous figure with a driving voltage as a horizontal axis and an acceleration detection voltage as a vertical axis.
FIG. 6 is a graph illustrating a regression line for a Lissajous figure.
FIG. 7 is a flowchart of a process in the personal computer 16 of FIG. 1;
FIG. 8 is a flowchart showing details of “validity check” in FIG. 7;
FIG. 9 is a view showing a setting screen of the personal computer 16;