(57)【要約】 【課題】高速かつ高精度に倣い測定できる表面形状測定装置を提供する。【解決手段】倣いベクトル指令部２２０から発せられた倣いベクトル指令に基づいて駆動機構の動作状態を推定して推定動作状態量を算出する動作推定部６００と、動作推定部６００にて算出された推定動作状態量に応じて駆動センサの検出値を補正演算する補正演算部７００と、を備える。  動作推定部６００は、倣いベクトル指令が発せられてから倣いプローブの移動位置に反映されるまでの信号伝達特性であるノミナルモデルが設定されたノミナルモデル設定部３１１を有し、補正演算部７００は駆動中の駆動機構が変形することによって生じる測定誤差を補正する補正量を推定動作状態量に基づいて算出する補正量算出部７２０と、駆動センサおよび検出センサによる検出値と補正量算出部７２０にて算出された補正量とを合成して測定データとする測定データ合成部４３０と、を備える。【選択図】図２
2008089578.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも一部に円弧部分が含まれる被測定物表面に近接あるいは当接される測定子、および、前記測定子と前記被測定物表面との相対位置を検出する検出センサを有し、前記測定子と前記被測定物表面との相対位置を予め設定された基準位置に保って倣い走査する倣いプローブと、
前記被測定物表面に沿った前記倣いプローブの次の移動位置を指示する倣いベクトルを指令する倣いベクトル指令部と、
前記倣いプローブを保持するとともに三次元的に移動させる駆動軸を有し、前記倣いベクトル指令に応じて前記倣いプローブを移動させる駆動機構と、
前記駆動機構の駆動量を検出する駆動センサと、
前記倣いベクトル指令部にて指令された前記倣いベクトル指令に応じて前記駆動センサの検出値を補正演算する補正演算部と、を備え、
前記補正演算部は、駆動中の前記駆動機構が変形することによって生じる測定誤差を補正する補正量を前記倣いベクトル指令の指令値に基づいて算出する
ことを特徴とする表面形状測定装置。
【請求項２】
請求項１に記載の表面形状測定装置において、
前記補正演算部は、
前記倣いプローブに生じる加速度と前記駆動機構の変形量との関係を表す補正パラメータを記憶した補正パラメータ記憶部と、
前記駆動センサの検出値を補正する補正量を前記倣いプローブに生じる加速度に基づいて算出する補正量算出部と、
前記算出された補正量を前記駆動センサおよび前記検出センサの検出値に合成して前記被測定物表面の位置を求める測定データ合成部と、を備える
ことを特徴とする表面形状測定装置。
【請求項３】
請求項２に記載の表面形状測定装置において、
前記補正量算出部は、位相差補正量算出部を備え、
前記補正パラメータ記憶部は、位相差補正のための補正パラメータを記憶しており、
前記位相差補正のための補正パラメータ（ｂ２、ｂ１）は、
φ＝ｂ２ｆ２＋ｂ１ｆ＋φ０の関係式（φは位相差補正量、ｆは倣い測定時の回転周波数、φ０は初期位相差）と、
測定領域中の異なる位置で円ゲージの倣い測定を異なる倣い速度で行った結果と、に基づいて算出されたものである
ことを特徴とする表面形状測定装置。
【請求項４】
請求項２または請求項３に記載の表面形状測定装置において、
前記補正量算出部は、ゲイン補正量算出部を備え、
前記補正パラメータ記憶部は、ゲイン補正のための補正パラメータを記憶しており、
前記ゲイン補正のための補正パラメータｋは、
D＝ｋ・Aの関係式（ここで、ゲイン補正量をD、倣い測定時における前記倣いプローブの移動の加速度をAとする）と、
測定領域中の異なる位置で円ゲージの倣い測定を異なる倣い速度で行った結果と、に基づいて算出されたものである
ことを特徴とする表面形状測定装置。
【請求項５】
請求項２または請求項３に記載の表面形状測定装置において、
前記補正量算出部は、ゲイン補正量算出部を備え、
前記補正パラメータ記憶部は、ゲイン補正のための補正パラメータを記憶しており、
前記ゲイン補正のための補正パラメータｋは、
｜Ｇ１｜＝１−ｋ（２π・ｆ）２の関係式（ここで、｜Ｇ１｜は信号伝達特性Ｇ１のゲイン、ｆは倣い測定時の回転周波数）と、
測定領域中の異なる位置で円ゲージの倣い測定を異なる倣い速度で行った結果と、に基づいて算出されたものである
ことを特徴とする表面形状測定装置。
【請求項６】
請求項２または請求項３に記載の表面形状測定装置において、
前記補正量算出部は、ゲイン補正量算出部を備え、
前記補正パラメータ記憶部は、ゲイン補正のための補正パラメータを記憶しており、
前記ゲイン補正のための補正パラメータは、
測定領域中の異なる位置で円ゲージの倣い測定を異なる倣い速度で行った結果から導出される回転周波数とゲインとの関係から導かれる多項式の係数として表される係数である
ことを特徴とする表面形状測定装置。
【請求項７】
請求項２から請求項６のいずれかに記載の表面形状測定装置において、
前記倣いベクトル指令部から発せられた前記倣いベクトル指令に基づいて前記駆動機構の動作状態を推定して推定動作状態量を算出する動作推定部を備え、
前記動作推定部は、
前記倣いベクトル指令部により前記倣いベクトル指令が発せられてから前記倣いプローブの移動位置に反映されるまでの信号伝達特性であるノミナルモデルが設定されたノミナルモデル設定部と、
前記倣いベクトル指令部からの前記倣いベクトル指令と前記ノミナルモデル設定部に設定された前記ノミナルモデルとを用いて前記倣いプローブの位置を推定位置として求める位置推定部と、
前記位置推定部にて求められた前記倣いプローブの推定位置を２階微分して前記推定動作状態量としての加速度を算出する２階微分演算部と、を備え、
前記補正演算部は、駆動中の前記駆動機構が変形することによって生じる測定誤差を補正する補正量を前記推定動作状態量に基づいて算出する
ことを特徴とする表面形状測定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、表面形状測定装置に関する。
例えば、被測定物表面を倣い走査して被測定物の輪郭、表面粗さ、うねり等を測定する表面倣い測定装置に関する。
【背景技術】
【０００２】
被測定物表面を倣い走査して被測定物の立体的形状を測定する表面倣い測定装置が知られている。
図１９に、倣いプローブ１３０を利用した表面倣い測定装置としての測定システム１００の構成を示す。
この測定システム１００は、倣いプローブ１３０を移動させる三次元測定機１１０と、手動操作する操作部１５０と、三次元測定機１１０の動作を制御するモーションコントローラ１６０と、モーションコントローラ１６０を介して三次元測定機１１０を動作させるとともに三次元測定機１１０によって取得した測定データを処理して被測定物Ｗの寸法や形状などを求めるホストコンピュータ２００と、を備えている。
【０００３】
三次元測定機１１０は、定盤１１１と、定盤１１１に立設されて倣いプローブ１３０を三次元的に移動させる駆動機構１２０と、駆動機構１２０の駆動量を検出する駆動センサ（不図示）と、を備えて構成されている。
駆動機構１２０は、定盤１１１の両側端から定盤１１１に略垂直方向であるＺｍ軸方向に高さを有するとともに定盤１１１の側端に沿ったＹｍ軸方向へスライド可能に設けられた二本のビーム支持体１２１と、ビーム支持体１２１の上端に支持されてＸｍ軸方向に長さを有するビーム１２２と、Ｘｍ軸方向にスライド可能にビーム１２２に設けられＺｍ軸方向にガイドを有するコラム１２３と、コラム１２３内をＺ軸方向にスライド可能に設けられ下端にて倣いプローブ１３０を保持するスピンドル１２４と、を備えて構成されている。
駆動センサは、ビーム支持体１２１のＹｍ軸方向への移動を検出するＹｍ軸センサと、コラム１２３のＸｍ軸方向への移動を検出するＸｍ軸センサと、スピンドル１２４のＺｍ軸方向への移動を検出するＺｍ軸センサと、を備えている。
【０００４】
倣いプローブ１３０は、図２０に示されるように、先端に接触部（測定子）１３２を有するスタイラス１３１と、スタイラス１３１の基端を一定の範囲内でＸｐ方向、Ｙｐ方向、Ｚｐ方向にスライド可能に支持する支持部１３３と、を備えている。
支持部１３３は互いに直交する方向に移動可能なｘｐスライダ、ｙｐスライダ、ｚｐスライダを有するスライド機構（不図示）と、スライド機構の各軸方向の変位量を検出するとともに検出した変位量を出力するプローブセンサ（不図示）と、を備えている。
スタイラス１３１はスライド機構によって支持部１３３に対して一定の範囲内でスライド可能に支持されている。
なお、このような倣いプローブ１３０の構成は例えば特許文献１（特開平０５−２５６６４０号）に記載されている。
【０００５】
このような構成において、接触部１３２を基準押込量Δｒで被測定物表面Ｓに当接させた状態で倣いプローブ１３０を被測定物表面Ｓに沿って倣い移動させる。
このとき、駆動機構１２０の駆動量から倣いプローブ１３０の移動軌跡が得られる。
そして、倣いプローブ１３０の移動軌跡は、すなわち、接触部１３２の移動軌跡となるところ、接触部１３２の中心点の移動軌跡に対して所定量（ΔＱ）だけオフセットした位置に被測定物表面Ｓと接触部１３２との接触点が存在する。
そこで、駆動センサによって検出される倣いプローブ１３０の位置とプローブセンサによって検出されるスタイラス１３１の変位とを合算して接触部１３２の位置を求めた上で、接触部１３２の位置から所定のオフセット値分（ΔＱ）だけ補正することによって被測定物表面Ｓの位置が算出される。
【０００６】
ここで、表面倣い測定装置によって被測定物表面を倣い走査する際、加速度をもって駆動する部分には慣性力が作用する。
例えば、被測定物Ｗが円あるいは円弧である場合、円運動による遠心力が発生し、図２１に示されるように、駆動機構１２０（スピンドル１２４）が変形するという問題が発生する。
このような加速度による変形が生じる場合、変形が生じる分だけセンサの検出値には誤差が含まれることとなる。
例えば、遠心力が発生する場合、スピンドル１２４が外側に変形する分だけセンサによる検出値は円の内側に入ってしまい、例えば、図２２に示されるように径方向のズレが生じる。
なお、図２２中、Ｌ１はリングゲージの径であり、Ｌ２は測定データである。
このような課題は、例えば、自動車のボディを測定するような大型の三次元測定機１１０で高速倣い測定を行う必要がある場合には非常に大きな問題となって顕在化する。
【０００７】
この点、例えば、特許文献２（特開平７−３２４９２８号）では、加速度によって生じる測定誤差を補正する方法として次の構成が開示されている。
すなわち、特許文献２にあっては、測定スライダの位置と測定スライダの加速度との関数としてたわみ特性を表す補正値を予め求めておく。
例えば、測定領域の様々な位置において各加速度で半径既知のリンクゲージを測定することにより加速度とたわみ特性との関数を求めておくことができる。
そして、被測定物を測定する際には、各センサによる検出データを得ることに加えて、さらに、測定時の加速度から補正値を特定して、前記検出データを補正値で補正する。
このようにして加速度によって生じる測定誤差を補正して正しい測定値が得られる。
測定時の加速度を求める方法として、特許文献２では、測定スライダの位置の測定値を二回微分することによって得る方法（段落００３７、請求項１２）、および、加速度センサを設けて測定スライダの加速度を検出する方法（段落００４７、請求項１３）が開示されている。
【０００８】
【特許文献１】特開平０５−２５６６４０号公報
【特許文献２】特開平０７−３２４９２８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００９】
特許文献２の方法では、測定スライダの位置を二回微分することで測定時の加速度を特定することとしているが、位置の検出値の二回微分では加速度の分解能がサンプリング周波数の二乗に反比例して劣化してしまうという問題がある。
例えば、位置検出のサンプリング周波数を１０倍に高めると、求められる加速度の分解能は１／１００に劣化し、その結果、補正量の分解能も１／１００に劣化してしまう。
このように、位置の二回微分による加速度の求め方は、実用的ではなく、高速での高精度測定という要求に応えることができない。
さらには、測定スライダの位置を二回微分して測定スライダの加速度を求めているが、実際の測定にあたって加速度および変形の発生が問題となるのはスピンドル１２４の先端部あるいはプローブ部であり、測定スライダの加速度を議論の対象としていては補正の精度に原理的に限界がある。
なお、特許文献２では、加速度センサを設けることによって測定スライダの加速度を求めることも記載されているが、加速度センサの性能と設置方法については記載されておらず、現実的には困難である。
例えば、加速度センサで実測する場合、例えば、直径１００ｍｍの円を倣い速度１０ｍｍ／ｓｅｃで倣い測定する際には約５０μＧで中心方向の加速度が発生するが、５０μＧの加速度を検出可能な加速度センサを移動軸毎に一つずつ合計３個設置することは困難である。
ましてや、プローブ部の近傍にこのような加速度センサを設置することは不可能である。
【００１０】
このように、倣い測定の際に生じる変形分を正確に求めることができないため、このような変形分を補正することができなかった。
そのため、高速で測定する場合には正確に被測定物の形状を求めることができず、正確に測定するためには変形が生じない程度の倣い速度に制限されるという問題が生じていた。
特に、自動車等の大きな被測定物を高速に測定するための大型の三次元測定機１１０が必要とされているため、上記の課題を解決する手段が切望されていた。
【００１１】
本発明の目的は、従来の問題を解消し、高速かつ高精度に倣い測定できる表面形状測定装置を提供することにある。
【課題を解決するための手段】
【００１２】
本発明の表面形状測定装置は、少なくとも一部に円弧部分が含まれる被測定物表面に近接あるいは当接される測定子、および、前記測定子と前記被測定物表面との相対位置を検出する検出センサを有し、前記測定子と前記被測定物表面との相対位置を予め設定された基準位置に保って倣い走査する倣いプローブと、前記被測定物表面に沿った前記倣いプローブの次の移動位置を指示する倣いベクトルを指令する倣いベクトル指令部と、前記倣いプローブを保持するとともに三次元的に移動させる駆動軸を有し、前記倣いベクトル指令に応じて前記倣いプローブを移動させる駆動機構と、前記駆動機構の駆動量を検出する駆動センサと、前記倣いベクトル指令部にて指令された前記倣いベクトル指令に応じて前記駆動センサの検出値を補正演算する補正演算部と、を備え、前記補正演算部は、駆動中の前記駆動機構が変形することによって生じる測定誤差を補正する補正量を前記倣いベクトル指令の指令値に基づいて算出することを特徴とする。
なお、少なくとも一部に円弧部分が含まれるとは、一定の曲率を有する円弧を含むことを要せず、曲面を含むものであればよい。よって、自由曲面も本発明の測定対象に含まれる。
【００１３】
このような構成において、倣いプローブを被測定物表面に沿って移動させるための倣いベクトル指令が倣いベクトル指令部から発せられる。
すると、倣いベクトル指令に応じ、駆動機構によって倣いプローブが被測定物表面に沿って倣い移動される。
倣いプローブで被測定物表面を走査するときの駆動機構の駆動量が駆動センサで検出されるとともに、測定子と被測定物表面との相対位置が検出センサで検出される。
このような倣い測定の動作の際には加速度が生じ、加速度による力がかかった場合には駆動機構に変形が生じ、その変形分が測定誤差に含まれる恐れがある。
そこで、倣いプローブで倣い測定すると同時に、上記の変形分による誤差を補正する処理が行われる。
補正演算部により、倣いベクトル指令の指令値に基づいた補正量が算出される。
すなわち、動作時に生じる変形に伴う測定誤差分を補正する補正量が算出される。
例えば、倣いベクトル指令の指令値から導出される加速度に基づいて駆動機構の変動を見積り、補正量を算出する。
算出された補正量は、駆動センサおよび検出センサによる検出値と合成されて測定データとして出力される。
【００１４】
このような構成によれば、倣いベクトル指令の指令値に基づいて補正量が算出される。
したがって、駆動機構の駆動状態を詳細にかつ正確に求めることができ、駆動中の駆動機構に生じている変形を補正する補正量を正確に算出することができる。
従来、測定誤差の補正のために駆動機構の駆動状態を求めるにあたっては駆動センサの検出値から駆動状態を求めることが行われていた。
しかし、駆動センサの検出値では、駆動時に発生する駆動機構の変形を含んだ値となっており、正確な値とはいえず、このような誤差を含んだ値による補正量の算出では正確な補正は望めないという問題があった。
また、駆動センサの検出値に基づいて動作状態を求める場合に、例えば、位置の検出値から加速度を算出しようとすると離散化誤差が大きくなって実用的ではないという問題が生じる。
例えば、駆動センサによる位置検出のサンプリング周期を１０分の１に短くすると、動作状態量としての加速度の分解能は１００分の１に劣化してしまう。
あるいは、動作状態として加速度データを得る場合に加速度センサを設けることも考えられるが、Ｘｍ、Ｙｍ、Ｚｍの各方向に加速度センサを設けることがそもそも困難であり、加速度センサの検出能力にも限界がある。
この点、本発明では、倣いベクトル指令の指令値から駆動状態を求めることとしているので、駆動機構の動作状態を求める場合でも実測のデータによることなく、演算によって求めることができる。
よって、短い周期であっても正確に動作状態を求めることができる。
その結果、高い分解性能で補正を行うことができる。
【００１５】
本発明では、前記補正演算部は、前記倣いプローブに生じる加速度と前記駆動機構の変形量との関係を表す補正パラメータを記憶した補正パラメータ記憶部と、前記駆動センサの検出値を補正する補正量を前記倣いプローブに生じる加速度に基づいて算出する補正量算出部と、前記算出された補正量を前記駆動センサおよび前記検出センサの検出値に合成して前記被測定物表面の位置を求める測定データ合成部と、を備えることが好ましい。
【００１６】
このような構成において、補正演算に必要なパラメータは補正パラメータ記憶部に予め記憶されている。
そして、補正量算出部では、倣いベクトル指令の指令値から導出される加速度および補正パラメータに基づいて補正量を算出する。
算出された補正量が駆動センサおよび検出センサの検出値に合成されることにより正確な測定値が求められる。
【００１７】
ここで、補正量としては、駆動軸方向ごとの信号伝達特性あるいは駆動ゲインの差異によって円運動時に各軸間に生じる指令伝達の位相差を補正する位相差補正量と、加速度により駆動機構が変形することによって生じる軸ごとのゲイン誤差を補正するゲイン補正量と、がある。
そして、位相差補正量を求める構成としては次の構成が好ましい。
すなわち、本発明では、前記補正量算出部は、位相差補正量算出部を備え、前記補正パラメータ記憶部は、位相差補正のための補正パラメータを記憶しており、前記位相差補正のための補正パラメータ（ｂ２、ｂ１）は、φ＝ｂ２ｆ２＋ｂ１ｆ＋φ０の関係式（φは位相差補正量、ｆは倣い測定時の回転周波数、φ０は初期位相差）と、測定領域中の異なる位置で円ゲージの倣い測定を異なる倣い速度で行った結果と、に基づいて算出されたものであることが好ましい。
【００１８】
このように位相差補正量は倣い測定時に倣いプローブが被測定物を走査する回転周波数に相関するとして前記関係式をたてる。
そして、測定領域中の異なる位置において異なる倣い速度（回転周波数）で円ゲージを実際に測定して校正データを得る。
この校正データを上記関係式に代入してパラメータを求め、補正パラメータとして記憶する。
このように回転周波数および測定位置の違いに対応して位相ずれを補正する補正パラメータを求めることにより、位相差補正を行うことができる。
したがって、円運動時の位相誤差を補正して正確な測定データを得ることができる。
【００１９】
また、ゲイン補正量を算出する構成としては次の構成が挙げられる。
すなわち、ゲイン誤差は加速度に比例するとした関係式を用いて加速度とゲイン誤差とを関係づけるパラメータを求めてもよく、
本発明では、前記補正量算出部は、ゲイン補正量算出部を備え、前記補正パラメータ記憶部は、ゲイン補正のための補正パラメータを記憶しており、前記ゲイン補正のための補正パラメータｋは、D＝ｋ・Aの関係式（ゲイン補正量をD、倣い測定時における前記倣いプローブの移動の加速度をAとする）と、測定領域中の異なる位置で円ゲージの倣い測定を異なる倣い速度で行った結果と、に基づいて算出されたものであることが好ましい。
【００２０】
また、ゲイン誤差は倣い測定時の回転周波数に相関するとした関係式を用いて回転周波数とゲイン誤差とを関係づけるパラメータを求めてもよく、
本発明では、前記補正量算出部は、ゲイン補正量算出部を備え、前記補正パラメータ記憶部は、ゲイン補正のための補正パラメータを記憶しており、前記ゲイン補正のための補正パラメータｋは、｜Ｇ１｜＝１−ｋ（２π・ｆ）２の関係式（｜Ｇ１｜は信号伝達特性Ｇ１のゲイン、ｆは倣い測定時の回転周波数）と、測定領域中の異なる位置で円ゲージの倣い測定を異なる倣い速度で行った結果と、に基づいて算出されたものであることが好ましい。
【００２１】
さらには、ゲイン誤差と回転周波数との関係を実験的に求めてもよく、
本発明では、前記補正量算出部は、ゲイン補正量算出部を備え、前記補正パラメータ記憶部は、ゲイン補正のための補正パラメータを記憶しており、前記ゲイン補正のための補正パラメータは、測定領域中の異なる位置で円ゲージの倣い測定を異なる倣い速度で行った結果から導出される回転周波数とゲインとの関係から導かれる多項式の係数として表される係数であることが好ましい。
【００２２】
前記発明では、補正演算は倣いベクトル指令の指令値を用いるとしたが、倣いベクトル指令の指令値をそのまま用いるのみならず、倣いベクトル指令値に基づいて駆動機構および倣いプローブの動作状態を推定して、この推定値に基づいて補正演算を行ってもよい。
すなわち、本発明では、前記倣いベクトル指令部から発せられた前記倣いベクトル指令に基づいて前記駆動機構の動作状態を推定して推定動作状態量を算出する動作推定部を備え、前記動作推定部は、前記倣いベクトル指令部により前記倣いベクトル指令が発せられてから前記倣いプローブの移動位置に反映されるまでの信号伝達特性であるノミナルモデルが設定されたノミナルモデル設定部と、前記倣いベクトル指令部からの前記倣いベクトル指令と前記ノミナルモデル設定部に設定された前記ノミナルモデルとを用いて前記倣いプローブの位置を推定位置として求める位置推定部と、前記位置推定部にて求められた前記倣いプローブの推定位置を２階微分して前記推定動作状態量としての加速度を算出する２階微分演算部と、を備え、前記補正演算部は、駆動中の前記駆動機構が変形することによって生じる測定誤差を補正する補正量を前記推定動作状態量に基づいて算出することが好ましい。
【００２３】
このような構成において、倣いベクトル指令部からの倣いベクトル指令が位置推定部に入力される。
すると、位置推定部において倣いベクトル指令による位置指令Ｃと設定されたノミナルモデルＧＮとが乗算され、倣いプローブの推定位置Ｅが算出される。
算出された倣いプローブの推定位置Ｅは２階微分演算部に出力される。
そして、２階微分演算部において、倣いプローブの推定位置から推定動作状態量として加速度データが算出される。
この加速度データは、補正量算出部に出力され、駆動動作によって生じる駆動機構の変形に伴う測定誤差分を補正する補正量が算出される。
算出された補正量は、駆動センサおよび検出センサによる検出値と合成され、測定データとして出力される。
【００２４】
このような構成によれば、推定動作状態量として倣いプローブに発生している加速度を求めることができる。
実測の位置データに基づいて加速度を算出する場合にはサンプリング周期を短くするに伴って加速度データの分解能が劣化してしまうが、本発明ではノミナルモデルを用いて算出された推定位置を２階微分することで加速度を算出するので、サンプリング周期による加速度データの劣化はない。
よって、短い周期であっても加速度データの算出が可能であり、このように求められた加速度データによって補正量を細かく算出することができる。
その結果、短い周期で取得された位置データであっても動作状態に基づいて正確に補正することができ、短いピッチでの正確な形状データを取得することができる。
【００２５】
このような構成によれば、２階微分演算部にて算出した推定加速度に基づいて円を倣い測定する際の回転周波数を周波数推定部にて算出するので、円の倣い動作中に変化する回転周波数も正確に求めることができる。
例えば、回転周波数ｆの回転速度で円を倣い測定する指令が倣いベクトル指令部によって発せられても、動作の開始直後で速度が安定する前には指令された通りの速度に達していないため、指令どおりの回転周波数ｆではない。
そのため、指令されたとおりの回転周波数ｆを前提にして補正量を算出したのでは、補正を正確に行うことができないという問題が生じる。
【００２６】
この点、本発明では、指令値の回転周波数ｆをそのまま補正演算の前提とすることなく、２階微分による加速度推定を行ったうえで、推定加速度に基づいたその都度の回転周波数ｆを算出する。
その結果、円の倣い動作中に変化する回転周波数も正確に求めることができるので、円運動時の遠心力分を回転周波数に基づく補正演算で補正する場合にあっても正確な補正を行うことができる。
【００２７】
また、例えば、円の倣い測定において、指令された回転周波数に達するまでには時間がかかるので、指令された回転周波数を前提とした補正を行う場合には加速部分で取得したデータを利用することなく、倣い速度（すなわち回転周波数）が一定となった際に取得したデータのみを補正して測定データとすることで正確な測定データを得ることはできる。
しかし、この方法では、倣い速度（すなわち回転周波数）が一定になるまでの時間が無駄となり、測定の作業効率が非常に悪くなる。
【００２８】
この点、本発明では、その都度の加速度推定値から回転周波数を算出し、この算出された回転周波数に基づいて補正量を算出できるので、倣い速度（すなわち回転周波数）が一定に達するのを待つことなく、総ての取得データを適切に補正して測定データとすることができる。
その結果、測定の作業効率を向上させることができる。
【発明を実施するための最良の形態】
【００２９】
以下、本発明の実施の形態を図示するとともに図中の各要素に付した符号を参照して説明する。
【００３０】
（第１実施形態）
本発明の表面形状測定装置に係る第１実施形態について説明する。
第１実施形態として倣いプローブ１３０を用いた表面倣い測定装置としての測定システム１００を図２０に示す。
また、図１に測定システム１００の機能ブロック図を示す。
測定システム１００の概略構成は、背景技術で説明した構成に同様であり、三次元測定機１１０と、三次元測定機１１０の動作を手動操作する操作部１５０と、三次元測定機１１０の駆動制御を実行するモーションコントローラ１６０と、モーションコントローラ１６０に所定の指令を与えるとともに被測定物Ｗの形状解析等の演算処理を実行するホストコンピュータ２００と、測定条件等を入力する入力手段６１と、測定結果を出力する出力手段６２と、を備えている。
【００３１】
三次元測定機１１０は（図１９参照）、定盤１１１と、定盤１１１に立設されて倣いプローブ１３０を三次元的に移動させる駆動機構１２０と、駆動機構１２０の駆動量を検出する駆動センサ１４０と、を備えて構成されている。
駆動機構１２０は、定盤１１１の両側端から定盤１１１に略垂直方向であるＺｍ軸方向に高さを有するとともに定盤１１１の側端に沿ったＹｍ軸方向へスライド可能に設けられた二本のビーム支持体１２１と、ビーム支持体１２１の上端に支持されてＸｍ軸方向に長さを有するビーム１２２と、ビーム１２２にＸｍ軸方向にスライド可能に設けられＺｍ軸方向にガイドを有するコラム１２３と、コラム１２３内をＺ軸方向にスライド可能に設けられ下端にて倣いプローブ１３０を保持するスピンドル１２４と、を備えて構成されている。
駆動センサ１４０は、ビーム支持体１２１のＹｍ軸方向への移動を検出するＹｍ軸センサ１４１と、コラム１２３のＸｍ軸方向への移動を検出するＸｍ軸センサ１４２と、スピンドル１２４のＺｍ軸方向への移動を検出するＺｍ軸センサ１４３と、を備えている。
【００３２】
倣いプローブ１３０は（図２０参照）、先端に接触部（測定子）１３２を有するスタイラス１３１と、スタイラス１３１の基端を一定の範囲内でＸｐ方向、Ｙｐ方向、Ｚｐ方向にスライド可能に支持する支持部１３３と、と備えている。
支持部１３３は互いに直交する方向に移動可能なｘｐスライダ、ｙｐスライダ、ｚｐスライダを有するスライド機構（不図示）と、スライド機構の各軸方向の変位量を検出するとともに検出した変位量を出力するプローブセンサ１３４と、を備えている。
プローブセンサ１３４は、スタイラス１３１のＸｐ方向への移動を検出するＸｐ方向センサ１３５と、スタイラス１３１のＹｐ方向への移動を検出するＹｐ方向センサ１３６と、スタイラス１３１のＺｐ方向への移動を検出するＺｐ方向センサ１３７と、を備えている。
【００３３】
モーションコントローラ１６０は、駆動センサ１４０およびプローブセンサ１３４からの検出信号を計数するカウンタ部１６１と、ホストコンピュータ２００および操作部１５０からの指令に応じて駆動機構１２０を駆動制御する駆動制御回路１６２と、を備えている。
【００３４】
カウンタ部１６１は、駆動センサ１４０から出力されるパルス信号をカウントして駆動機構１２０の駆動量を計測する駆動カウンタ１７１と、プローブセンサ１３４から出力されるパルス信号をカウントしてスタイラス１３１のスライド量を押込量として計測するプローブカウンタ１７５と、を備えている。
駆動カウンタ１７１は、Ｙｍ軸センサ１４１からの出力を計数するＹｍ軸カウンタ１７２と、Ｘｍ軸センサ１４２からの出力を計数するＸｍ軸カウンタ１７３と、Ｚｍ軸センサ１４３からの出力を計数するＺｍ軸カウンタ１７４と、を備えている。
プローブカウンタ１７５は、Ｘｐ方向センサ１３５からの出力を計数するＸｐ方向カウンタ１７６と、Ｙｐ方向センサ１３６からの出力を計数するＹｐ方向カウンタ１７７と、Ｚｐ方向センサ１３７からの出力を計数するＺｐ方向カウンタ１７８と、を備えている。駆動カウンタ１７１によるカウント値（Ｘｍ、Ｙｍ、Ｚｍ）およびプローブカウンタ１７５によるカウント値（Ｘｐ、Ｙｐ、Ｚｐ）はそれぞれホストコンピュータ２００に出力される。
【００３５】
ホストコンピュータ２００は、入力手段６１によって設定入力される測定条件等を記憶するメモリ（記憶装置）２１０と、被測定物表面に倣う移動方向および移動速度の倣いベクトルを指令する倣いベクトル指令部２２０と、倣いベクトル指令部２２０からの指令に基づいて駆動機構１２０の動作を推定して倣いプローブ１３０に発生する加速度を算出する動作推定部６００と、動作推定部６００による動作推定に基づいて補正量を算出して測定値を補正する補正演算部７００と、被測定物Ｗの形状を解析する形状解析部５００と、演算装置および記憶装置（ＲＯＭ、ＲＡＭ）を有し所定プログラムの実行やデータ処理等を行う中央処理部（ＣＰＵ）２３０と、前記各機能部を接続するバスと、を備えている。
【００３６】
メモリ２１０は、入力手段６１から設定入力される測定条件等を記憶し、例えば、倣い走査中に駆動機構１２０の駆動量をサンプリングする間隔（サンプリングピッチ）や、接触部１３２を被測定物Ｗに対して押し込む量（基準押込量Δｒ）や、倣い走査を行う速度（倣い速度）や、被測定物Ｗの設計データなどによる輪郭データなどを記憶する。
【００３７】
倣いベクトル指令部２２０は、例えば、メモリ２１０に設定された輪郭データに基づいて被測定物Ｗを倣い走査する倣いベクトル指令を生成する。
また、倣いベクトル指令部２２０は、プローブカウンタ１７５の出力に基づいて押込量を基準押込量Δｒで所定範囲（基準位置範囲）にする押込方向の倣いベクトル指令を生成する。
倣いベクトル指令部２２０で生成された倣いベクトル指令は、駆動制御回路１６２に出力される。
【００３８】
図２は、動作推定部６００および補正演算部７００の構成を示す図である。
動作推定部６００は、加速度推定部３１０と周波数推定部６１０とを備えている。
加速度推定部３１０は、倣いベクトル指令部２２０による倣いベクトル指令が発せられてから駆動センサ１４０による位置検出が行われるまでの伝達関数であるノミナルモデルが設定されたノミナルモデル設定部３１１と、倣いベクトル指令とノミナルモデルとに基づいて倣いプローブ１３０の位置を推定する位置推定部３１４と、位置推定部３１４にて求められた推定位置のデータ値から２階微分により倣いプローブ１３０の加速度を算出する２階微分演算部３１５と、を備えている。
【００３９】
図３は、ノミナルモデル設定部の構成を示す図である。
ノミナルモデル設定部３１１には、倣いベクトル指令が発せられてから三次元測定機１１０が駆動されて指令位置に到達するまでの周波数伝達関数であるノミナルモデルが設定されている。
ノミナルモデル設定部３１１には、図３に示されるように、三次元測定機のノミナルモデルＧ１（ｓ）（駆動機構のノミナルモデル）と、倣いプローブのノミナルモデルＧ２（ｓ）とが設定されている。
倣いプローブのノミナルモデルは、倣いプローブ１３０の交換に備えて倣いプローブ１３０の種類ごとに用意されており、現に取り付けている倣いプローブ１３０に対応したノミナルモデルを選択可能となっている。
そして、倣いベクトル指令から倣いプローブ１３０の位置までの伝達関数は三次元測定機のノミナルモデルＧ１（ｓ）と倣いプローブのノミナルモデルＧ２（ｓ）との積であるＧＮ（＝Ｇ１（ｓ）×Ｇ２（ｓ））として表される。
なお、ベクトル指令、三次元測定機１１０の駆動動作およびセンサによる位置検出等はＸｍ軸、Ｙｍ軸およびＺｍ軸方向でそれぞれ行われるが、三次元測定機１１０の位置制御系（倣いベクトル指令から位置検出を行うフィードバック制御系）はＸｍ、Ｙｍ、Ｚｍ軸方向の駆動の時定数Ｔが総て同じになるように調整されていることから、Ｘｍ、Ｙｍ、Ｚｍ軸方向ともに同じノミナルモデルＧＮとする。
ノミナルモデルＧＮは設計データあるいは実験値等から適宜求められる。
例えば、ノミナルモデルＧＮが一次遅れ系であるとき、ノミナルモデルＧＮは次の式で表される。
次式においてＫはゲインを表し、ｓはラプラス演算子を表す。
【００４０】
【数１】
2008089578.tif 000003
【００４１】
位置推定部３１４は、倣いベクトル指令の位置指令Ｃ（ｃｘ、ｃｙ、ｃｚ）とノミナルモデルＧＮとにより、倣いプローブ１３０の位置を推定する。
位置推定部３１４には、倣いベクトル指令部２２０にて生成された倣いベクトル指令が刻々入力される。
そして、位置指令Ｃ（ｃｘ、ｃｙ、ｃｚ）とノミナルモデルＧＮとにより倣いプローブ１３０の推定位置Ｅ（ｅｘ、ｅｙ、ｅｚ）を次のように算出する。
すなわち、倣いプローブ１３０の推定位置ＥはノミナルモデルＧＮと位置指令Ｃとの積で表される。
【００４２】
【数２】
2008089578.tif 000004
【００４３】
推定位置を各方向（Ｘｍ軸方向、Ｙｍ軸方向、Ｚｍ軸方向）でそれぞれ分解して表すと次のようになる。
【００４４】
【数３】
2008089578.tif 000005
【００４５】
２階微分演算部３１５は、位置推定部３１４にて算出された倣いプローブ１３０の推定位置Ｅ（ｅｘ、ｅｙ、ｅｚ）を２階微分することにより、倣いプローブ１３０に発生している加速度Ａ（ａｘ、ａｙ、ａｚ）を推定動作状態量として求める。
【００４６】
【数４】
2008089578.tif 000006
【００４７】
【数５】
2008089578.tif 000007
【００４８】
周波数推定部６１０は、倣いプローブ１３０を用いて、角速度一定で円を倣い測定する際の回転周波数ｆを推定する。
ここで、２階微分演算部３１５により加速度Ａが算出されているところ、周波数推定部６１０には２階微分演算部３１５にて算出された加速度Ａ（ａｘ、ａｙ、ａｚ）が入力される。
そして、例えば、倣いベクトル指令による円運動倣い測定の回転半径をＲｓとし、２階微分演算部３１５にて算出された倣いプローブ１３０の推定加速度をＡ（ａｘ、ａｙ、ａｚ）とするとき、周波数推定部６１０は、次の式により円運動の回転周波数ｆを算出する。
【００４９】
【数６】
2008089578.tif 000008
【００５０】
この式の導出については、例えば、次のように説明される。
回転周波数ｆと角速度ωとの関係はω＝２πｆと表され、また、向心加速度をａｎと表すときａｎ＝Ｒｓω２であることから、次の式変形により導出可能である。
【００５１】
【数７】
2008089578.tif 000009
【００５２】
また、特に、ＸＹ平面における円運動である場合には、上式は次のように簡略化される。
【００５３】
【数８】
2008089578.tif 000010
【００５４】
補正演算部７００は、倣いプローブ１３０に生じた加速度と三次元測定機１１０の変形量との関係を表す補正パラメータを記憶した補正パラメータ記憶部７１０と、センサ検出値を補正する補正量Ｄを倣いプローブ１３０に生じた加速度に基づいて算出する補正量算出部７２０と、算出された補正量Ｄとセンサ検出値を合成して被測定物表面Ｓの位置を求める測定データ合成部４３０と、を備えている。
【００５５】
補正パラメータ記憶部７１０には、補正パラメータＰが設定記憶されている。
補正パラメータＰは、倣いプローブ１３０に加速度Ａが働くときに三次元測定機１１０に生じた変形量（具体的にはコラム１２３およびスピンドル１２４に生じた変形量）を加速度Ａと関係づけるパラメータである。
補正パラメータＰ（ｐｘ、ｐｙ、ｐｚ）は、各軸方向、すなわち、Ｘｍ、Ｙｍ、Ｚｍ軸方向ごとに設定されている。
なお、補正パラメータＰの設定方法については後述する変形例１、変形例２、変形例３において円運動である場合を例にして示す。
【００５６】
補正量算出部７２０は、ゲイン補正量算出部７２１と、位相差補正量算出部７２２と、を備えている。
【００５７】
ゲイン補正量算出部７２１には、動作推定部６００にて算出された加速度Ａが入力される。
そして、ゲイン補正量算出部７２１は、補正パラメータ記憶部７１０に設定されている補正パラメータＰを用いて、次のように補正量Ｄを算出する。
【００５８】
【数９】
2008089578.tif 000011
【００５９】
【数１０】
2008089578.tif 000012
【００６０】
位相差補正量算出部７２２は、駆動軸方向ごとの信号伝達特性あるいは駆動ゲインの差異によって円運動時に生じる各軸間の指令伝達の位相差を位相差補正量φとして算出する。
位相差補正量φは、回転周波数ｆの関数（例えば、φ＝ｂ２ｆ２＋ｂ１ｆ＋ｂ０）として表され、理論的あるいは実験的に求められた補正パラメータが補正パラメータ記憶部７１０に設定記憶されている。
なお、ゲイン補正量Ｄおよび位相差補正量φを算出するための補正パラメータを設定する方法については後述の変形例１、変形例２、変形例３において例を挙げて説明する。
【００６１】
測定データ合成部４３０は、駆動センサ１４０による検出値（ｘｍ、ｙｍ、ｚｍ）、補正量算出部７２０で算出された補正量Ｄ（ｄｘ、ｄｙ、ｄｚ）および位相差補正量φ、および、プローブセンサ１３４による検出値（ｘｐ、ｙｐ、ｚｐ）を合成して倣いプローブ１３０の接触部１３２の軌跡を算出する。
さらに、基準押込量Δｒを考慮して被測定物表面Ｓの位置データを求める。
測定データ合成部４３０におけるデータ合成について説明する。
説明のために、ＸＹ平面での円運動による倣い移動で被測定物表面を測定した場合を例にして説明する。
円運動の中心座標を（ｘ０、ｙ０）とし、円の半径をＲで表すとき、駆動センサ１４０による検出値（ｘｍ、ｙｍ）は、所定の位相θを用いて次のように表される。
【００６２】
【数１１】
2008089578.tif 000013
【００６３】
そして、算出されたゲイン補正量（ｄｘ、ｄｙ）と、位相差補正量φの分を補正すると、次のように表される。
なお、ｘｍ′、ｙｍ′は、補正後の駆動センサ検出値を表す。
【００６４】
【数１２】
2008089578.tif 000014
【００６５】
この補正された駆動センサ１４０の検出値（ｘｍ′、ｙｍ′）とプローブセンサ１３４の検出値とを合成して求められる接触部１３２の位置に対して基準押込量Δｒを考慮した位置に被測定物の表面位置が求められる。
算出された被測定物表面Ｓの位置データは形状解析部５００に出力される。
【００６６】
形状解析部５００は、データを合成して接触部１３２の軌跡および被測定物表面Ｓの形状データを算出する。
さらに、算出された被測定物の形状データを設計データと対比して誤差や歪みなどを求めるなどの形状解析も行う。
【００６７】
このような構成を備える第１実施形態の動作を説明する。
まず、測定に先立って測定条件を設定入力する。
測定条件としてはサンプリングピッチ、基準押込量Δｒ、被測定物の輪郭データなどが挙げられる。
また、ノミナルモデル設定部３１１に三次元測定機のノミナルモデルＧ１（ｓ）と倣いプローブのノミナルモデルＧ２（ｓ）とを設定する。
倣いプローブ１３０のノミナルモデルは複数用意されているもののなかから現に使用する倣いプローブ１３０に対応するモデルを選択する。
この状態で、測定を開始させると、まず、倣いベクトル指令部２２０で生成される倣いベクトル指令が駆動制御回路１６２に出力される。
すると、駆動制御回路１６２から駆動機構１２０に制御信号が出力されて駆動機構１２０が駆動される。
駆動機構１２０によって倣いプローブ１３０が被測定物表面Ｓに対して基準押込量Δｒまで押し込まれた状態で被測定物表面Ｓに沿って倣い移動される。
この倣い移動時に加速度が生じ、三次元測定機１１０のスピンドル１２４に変形が生じる（例えば図２１参照）。
また、倣い走査時にプローブカウンタ１７５から倣いベクトル指令部２２０に出力されるプローブカウンタ値に基づいて、押込量が基準押込量Δｒに制御される。
【００６８】
倣いプローブ１３０で被測定物表面Ｓを走査するときの駆動機構１２０の駆動量が駆動センサ１４０で検出され、スタイラス１３１の変位量がプローブセンサ１３４で検出される。
駆動センサ１４０のセンサ出力は駆動カウンタ１７１でカウントされ、プローブセンサ１３４のセンサ出力はプローブカウンタ１７５でカウントされる。
なお、このようなカウンタ部１６１でカウントされたデータは、設定されたサンプリングピッチで取得される。
カウンタ部１６１（駆動カウンタ１７１およびプローブカウンタ１７５）でカウントされて取得されたカウント値は補正演算部７００にも出力される。
【００６９】
倣いベクトル指令部２２０からの倣いベクトル指令は駆動制御回路１６２と同時に動作推定部６００にも送られる。
動作推定部６００において倣いベクトル指令は位置推定部３１４に入力される。
位置推定部３１４は、ノミナルモデル設定部３１１に設定されたノミナルモデルと倣いベクトル指令の位置指令とから周波数伝達速度を考慮した現在時点における倣いプローブ１３０の推定位置Ｅ（ｔ）を算出する。
すなわち、位置指令Ｃ（ｃｘ、ｃｙ、ｃｚ）とノミナルモデルＧＮとを用いて次のように倣いプローブ１３０の推定位置Ｅ（ｔ）（ｅｘ、ｅｙ、ｅｚ）を算出する。
【００７０】
【数１３】
2008089578.tif 000015
【００７１】
位置推定部３１４にて算出された推定位置Ｅは２階微分演算部３１５に入力される。
２階微分演算部３１５は、算出された推定位置を２階微分することにより、現在時点における倣いプローブ１３０に生じている加速度Ａ（ａｘ、ａｙ、ａｚ）を算出する。
【００７２】
【数１４】
2008089578.tif 000016
【００７３】
２階微分演算部３１５にて算出された加速度Ａは、補正演算部７００の補正量算出部７２０に出力される。
ゲイン補正量算出部７２１は、入力された加速度Ａと補正パラメータとにより、加速度によって生じる変形分を補正する補正量Ｄを軸ごとに算出する。
すなわち、加速度Ａが倣いプローブ１３０に生じている場合における駆動センサ検出値の補正量Ｄ（ｄｘ、ｄｙ、ｄｚ）は次のように算出される。
【００７４】
【数１５】
2008089578.tif 000017
【００７５】
さらに、２階微分演算部３１５にて算出された加速度Ａは、周波数推定部６１０に出力される。
周波数推定部６１０において円運動の回転周波数ｆが算出され、算出された回転周波数ｆは位相差補正量算出部７２２に出力される。
そして、位相差補正量算出部７２２により、回転周波数ｆの関数としての位相差補正量φが算出される。
算出された位相差補正量φはゲイン補正量算出部７２１で算出されたゲイン補正量Ｄ（ｄｘ、ｄｙ、ｄｚ）とともに測定データ合成部４３０に出力される。
【００７６】
測定データ合成部４３０は、駆動センサ１４０による検出値（ｘｍ、ｙｍ、ｚｍ）と、算出された補正量Ｄ（ｄｘ、ｄｙ、ｄｚ）および位相差補正量φと、プローブセンサ１３４の検出値（ｘｐ、ｙｐ、ｚｐ）と、を合成して接触部１３２の位置を求める。
さらに、接触部１３２の位置に対して所定の基準押込量Δｒを考慮した位置に被測定物表面Ｓを求める。
求められた被測定物Ｗの表面データは順次形状解析部５００に送られ、被測定物表面Ｓの形状データが求められる。
そして、被測定物Ｗの形状データを設計データとの対比により誤差や歪みなどを求めるなどの形状解析が行われる。
【００７７】
このような構成を備える第１実施形態によれば次の効果を奏する。
（１）動作推定部６００がノミナルモデル設定部３１１を有し、ノミナルモデルに基づく加速度の算出により駆動状態にある駆動機構１２０の動作状態として加速度を正確に算出することができる。
そして、この推定加速度に基づいて補正量算出部７２０にて補正量を算出するので、駆動中の駆動機構１２０に生じている変形を補正する補正量を正確に算出することができる。駆動機構１２０の動作状態を求める場合でも実測のデータによることなく、ノミナルモデルに基づいた演算によって求めることができるので、高い分解性能で補正を行うことができる。
【００７８】
（２）周波数推定部６１０にて算出した円運動時の回転周波数ｆに基づき駆動軸間の位相差を補正する。
したがって、円運動時の位相差誤差を補正して正確な測定データを得ることができる。
【００７９】
（３）２階微分演算部３１５にて算出した推定加速度に基づいて円を倣い測定する際の回転周波数を周波数推定部６１０にて算出するので、円の倣い動作中に変化する回転周波数も正確に求めることができる。
【００８０】
（４）円の倣い測定において、指令された回転周波数に達するまでには時間がかかるので、指令された回転周波数を前提とした補正を行う場合には加速部分で取得したデータを利用することなく、倣い速度（すなわち回転周波数）が一定となった際に取得したデータのみを補正して測定データとしなければならない。
例えば、図４に示すように、倣い速度（あるいは回転周波数）が一定値を示すところで１周分（３６０°）のデータを取らなければならない。
これに対し本実施形態では、その都度の加速度推定値から回転周波数ｆを算出し、この算出された回転周波数ｆに基づいて補正量を算出できるので、倣い速度（すなわち回転周波数）が一定に達するのを待つことなく、総ての取得データを適切に補正して測定データとすることができる。
例えば、図５に示すように、倣い速度（あるいは回転周波数）が変化（加速あるいは減速）している領域で取得したデータでも正確に補正することができる。
よって、全体として１周分（３６０°）のデータを取得すればよく、その結果、測定の作業効率を向上させることができる。
【００８１】
（変形例１）
次に、本発明の変形例１について説明する。
変形例１は、基本的に第１実施形態と同様であるが、補正パラメータを設定する際の補正パラメータの求め方に特徴を有する。
第１実施形態の補正演算部７００で算出する補正量としては、ゲイン補正量算出部７２１によるゲイン補正量と、位相差補正量算出部７２２による位相差補正量とがある。
【００８２】
まず、ゲイン補正量を算出するための補正パラメータの設定について説明する。
変形例１にあっては、ゲイン補正量の補正パラメータを求めるにあたり、スピンドル１２４の変形量Ｄは加速度Ａに相関するとして次の式を立てる。
【００８３】
【数１６】
2008089578.tif 000018
【００８４】
ここで、説明を単純化するために、ＸＹ面内において円運動する場合を例に説明する。
円運動時の向心加速度をａとすると、Ｘｍ軸方向の加速度ａｘは、倣い速度をＶ、半径をＲｓとして次のように表される。
【００８５】
【数１７】
2008089578.tif 000019
【００８６】
加速度が生じた場合の変形量は、測定位置にも相関することから、測定座標の２乗の項までを考慮にいれると、Ｘｍ軸方向の変形量ｄｘはＸｍ軸方向に生じている加速度ａｘとの関係で次のように表される。
【００８７】
【数１８】
2008089578.tif 000020
【００８８】
Ｙｍ軸方向の補正量ｄｙ、Ｚｍ軸方向の補正量ｄｚを含めると、次のように表される。
【００８９】
【数１９】
2008089578.tif 000021
【００９０】
ただし、変形量Ｄは、ｘ、ｙ、ｚ成分の総てに相関するため、補正の係数は次のように表される。
【００９１】
【数２０】
2008089578.tif 000022
【００９２】
このような前提のもと、実際に実験して誤差を求めたうえ、フィッティングカーブの係数として補正パラメータを求めることができる。
【００９３】
実際の例を示す。
加速度による変形量は、場所によっても異なるため、校正用データの取得にあたっては、異なる複数の位置に標準ゲージをおいて円の倣い測定を複数の倣い速度（あるいは加速度）で校正試験を行う。
例えば、図６に示されるように、Ｐ１からＰ５の複数の位置で円の倣い測定を行う。
なお、図６では説明の都合上、Ｙｍ座標は同じでＸｍ座標が異なる５点を選定しているが、測定システム１００の測定領域においてＸｍ座標、Ｙｍ座標、Ｚｍ座標のすべてをまんべんなくふることがより好ましい。
【００９４】
図７および図８は、各ポイントＰ１〜Ｐ５で円の倣い測定をした際に、倣い速度Ｖと検出された半径の減少量（半径誤差ΔＲ）との関係を示したグラフである。
ここで、円の倣い測定時には、図２２に示されるように、Ｘ軸方向とＹ軸方向とでゲインに差があるために楕円誤差が生じる。
そのため、図７の縦軸は半径誤差ΔＲのＸｍ軸方向の最大量を示し、図８の縦軸は半径誤差ΔＲのＹｍ軸方向の最大量を示す。
【００９５】
図９および図１０は、図７と図８中の倣い速度を加速度に変換して、半径の減少量（半径誤差ΔＲ）と加速度との関係を示したグラフである。
図９の縦軸は図７と同様に半径誤差ΔＲのＸｍ軸方向の最大量を示し、図１０の縦軸は図８と同様に半径誤差のＹｍ軸方向の最大量を示す。
図７および図８においては、Ｘｍ軸方向でもＹｍ軸方向でも測定半径の減少量（半径誤差ΔＲ）は倣い速度の累乗（例えば２乗）に相関していることがわかるが、図９および図１０によれば、測定半径の減少量（半径誤差ΔＲ）は加速度に比例することがわかり、上記式（Ｄ＝ｋ・Ａ）が適用できることが示されている。
【００９６】
図１１および図１２は、測定位置と測定半径の減少量（半径誤差ΔＲ）との関係を示す図である。
図１１および図１２において、縦軸は、測定半径の減少量（半径誤差ΔＲ）を加速度で割って正規化した値（ΔＲ／ａ）で示している。
そして、図１１の縦軸は半径誤差ΔＲのＸｍ軸方向の最大量（ΔＲｘ）を加速度で割った値（ΔＲｘ／ａ）であり、図１２の縦軸は半径誤差ΔＲのＹｍ軸方向の最大量（ΔＲｙ）を加速度で割った値（ΔＲｙ／ａ）である。
図１１において、Ｘｍ軸方向では測定半径の減少量（ΔＲｘ）は測定位置のＸｍ座標の一次に相関していることがわかる。
さらに、図１２においては、Ｙｍ軸方向の測定半径の減少量（ΔＲｙ）であっても、Ｘｍ軸方向における測定位置の違いに影響され、Ｘｍ座標の二次に相関していることがわかる。
【００９７】
このような実験データを収集して、Ｘｍ、Ｙｍ、Ｚｍ座標が異なる位置において、異なる倣い速度（あるいは加速度）で標準ゲージを測定したときの誤差（ｄｘ、ｄｙ、ｄｚ）を式に代入する。
そして、この式を解くことにより、ゲイン補正のパラメータが求められる。
求められたゲイン補正のパラメータは補正パラメータ記憶部７１０に設定する。
【００９８】
次に、位相差補正量φを算出するための補正パラメータの設定について説明する。
駆動軸方向ごとの信号伝達特性あるいは駆動ゲインの差異によって円運動時に各軸間で指令伝達の位相差が生じる。
この位相差を位相差補正量φとし、円運動の回転周波数ｆの関数として求める。
説明を単純化するために、ＸＹ面内において円運動する場合を例に説明する。
位相差補正量φは、円運動の回転周波数ｆの関数であり、回転周波数ｆの多項式として表される。
実際の例で示す。
図１３は、図６に示したＰ１からＰ５の複数の位置で円の倣い測定を行った結果である。
図１３においては、縦軸にＹｍ軸方向とＸｍ軸方向とで生じる位相差φをとり、横軸に円運動の回転周波数ｆをとったものである。
図１３に示されるように、Ｙｍ軸方向とＸｍ軸方向との位相差φは、円運動の回転周波数ｆに依存する。
また、影響は小さいながらも測定位置にも依存する。
そこで、位相差φを、φ＝ｂ２ｆ２＋ｂ１ｆ＋φ０として表し、フィッティングカーブから係数ｂ２、ｂ１およびφ０を求めればよい。
【００９９】
位相差φの係数ｂ２、ｂ１およびφ０は、測定位置にも影響されることから、位相差φは次のように表される。
【０１００】
【数２１】
2008089578.tif 000023
【０１０１】
このように求められた補正パラメータを用い、第１実施形態に説明したようにゲイン補正量および位相差補正量φを算出する。
そして、ゲイン補正量および位相差補正量を用いて第１実施形態で説明したように測定データの補正を行うことができる。
【０１０２】
（変形例２）
次に、本発明の変形例２について説明する。
変形例２は、基本的に変形例１と同様であるが、補正パラメータを設定する際の補正パラメータの求め方に特徴を有する。
具体的には、変形例１ではスピンドル１２４の変形量は加速度に相関するとして変形量を規定する式を立式したが、変形例２では、スピンドル１２４の変形量は回転周波数ｆに相関するとして立てた式を利用する。
【０１０３】
なお、補正量としては、ゲイン補正量と位相差補正量とがあるが、位相差補正量を算出する補正パラメータの設定にあっては変形例１と同じ方法が採用できるので説明を省略する。
【０１０４】
ゲイン補正量を算出するための補正パラメータの設定について説明する。
変形例２にあっては、ゲイン補正量は回転周波数ｆに相関するとし、駆動の信号伝達特性Ｇ１のゲイン｜Ｇ１｜と円運動の周波数ｆとの関係を次の式で表す。
【０１０５】
【数２２】
2008089578.tif 000024
【０１０６】
この式の導出としては、次のような式変形が一例として挙げられる。
半径Ｒｓの円運動を指令したときに変形量Ｄが生じている場合、信号伝達特性Ｇ１のゲインは次の式で表される。
【０１０７】
【数２３】
2008089578.tif 000025
【０１０８】
変形量Ｄは加速度Ａに比例してＤ＝ｋ・Ａであるので、さらに次のように変形される。
【０１０９】
【数２４】
2008089578.tif 000026
【０１１０】
ここで、加速度Ａは、円運動の半径Ｒｓと角速度ωとの関係から回転周波数ｆを用いた次の式で表される。
【０１１１】
【数２５】
2008089578.tif 000027
【０１１２】
したがって、回転周波数特性のゲイン｜Ｇ１｜は次のように表される。
【０１１３】
【数２６】
2008089578.tif 000028
【０１１４】
このような関係式に基づいて、ゲインと回転周波数ｆとの関係を決める補正パラメータｋを設計値あるいは実験値によって求める。
【０１１５】
ゲイン｜Ｇ１｜は、Ｘｍ軸方向分、Ｙｍ軸方向分でそれぞれ算出する。
すなわち、Ｘｍ軸方向分は｜Ｇ１Ｘ｜、Ｙｍ軸方向成分は｜Ｇ１Ｙ｜、と算出する。
実際に｜Ｇ１Ｘ｜、｜Ｇ１Ｙ｜を求めるにあたっては、変形例１と同様に、異なる位置において、異なる回転周波数ｆで円の倣い測定を行い、その結果に基づいて位置（ｘｍ、ｙｍ）と回転周波数ｆとの関数として｜Ｇ１Ｘ｜、｜Ｇ１Ｙ｜を表せばよい。
【０１１６】
さらに、このように求められたゲイン補正量（｜Ｇ１Ｘ｜、｜Ｇ１Ｙ｜）と位相差補正量φと、を考慮した測定値の補正は次のように実行する。
すなわち、駆動センサ１４０による検出値が（ｘｍ、ｙｍ）であるとき、円運動の中心座標を（ｘ０、ｙ０）とし、円の半径をＲｓで表すとき、所定の位相θを用いて駆動センサ１４０の検出値は次のように表される。
【０１１７】
【数２７】
2008089578.tif 000029
【０１１８】
そして、算出されたゲイン補正量（｜Ｇ１Ｘ｜、｜Ｇ１Ｙ｜）と、位相差φの分を補正すると、次のように表される。
なお、ｘｍ′、ｙｍ′は、補正後の駆動センサ検出値を表す。
【０１１９】
【数２８】
2008089578.tif 000030
【０１２０】
この補正された駆動センサ１４０の検出値（ｘｍ′、ｙｍ′）とプローブセンサ１３４の検出値とを合成して求められる接触部１３２の位置に対して所定の基準押込量Δｒを考慮した位置に求められる。
【０１２１】
以上説明したように変形例２の方法によってもゲイン補正量を算出するための補正パラメータを求めることができる。
さらに、ゲイン補正量および位相差補正量を用いて測定データの補正を行うことができる。
【０１２２】
（変形例３）
次に、本発明の変形例３について説明する。
変形例３は、基本的に変形例２と同様であるが、補正パラメータを設定する際の補正パラメータの求め方に特徴を有する。
具体的には、変形例２では変形量が回転周波数ｆに相関する理論式で表されるという前提のもとに補正式を求めたが、変形例３では、変形量は回転周波数ｆを用いた多項式で表されるとして補正の式を求める。
【０１２３】
すなわち、Ｘｍ軸方向のゲイン補正量｜Ｇ１Ｘ｜およびＹｍ軸方向のゲイン補正量｜Ｇ１Ｙ｜は、円の倣い測定における円運動の回転周波数ｆに関する多項式で表されるとする。
実際に実験した例を示す。
校正データの取得にあたっては、異なる複数の位置に標準ゲージをおいて複数の倣い速度で円の倣い測定を行う。
図１４は、横軸に倣い測定時の回転周波数をとり、縦軸に半径誤差ΔＲのＸｍ軸方向の最大量をとったグラフである。
図１５は、横軸に倣い測定時の回転周波数をとり、縦軸に半径誤差ΔＲのＹｍ軸方向の最大量をとったグラフである。
測定位置は変形例１と同様であり、図６に示されるＰ１からＰ５の複数の位置で円の倣い測定を行った結果である。
図１４のＸｍ軸方向、図１５のＹｍ軸方向のいずれにおいても、倣い測定の回転周波数と測定誤差（半径誤差ΔＲ）との関係を回転周波数ｆの二次の多項式でフィッティングでき、所定の係数により、Ｘｍ軸方向のゲイン誤差｜Ｇ１Ｘ｜（半径誤差ΔＲのＸｍ方向）およびＹｍ軸方向のゲイン誤差｜Ｇ１Ｙ｜（半径誤差ΔＲのＹｍ方向）は、次のように表される。
【０１２４】
【数２９】
2008089578.tif 000031
【０１２５】
さらに、変形量は測定する位置によっても変化するため測定位置の関数でもある。
図１６は、Ｘｍ軸方向のゲイン誤差｜Ｇ１Ｘ｜の係数とＸｍ座標との関係を示すグラフである。
図１７は、Ｙｍ軸方向のゲイン誤差｜Ｇ１Ｙ｜の係数とＹｍ座標との関係を示すグラフである。
すなわち、多項式の係数ψ２、ψ１、η２、η１、は、測定位置の関数であり次のように表される。
Ｘｍ軸方向ゲインを求める式に用いる係数ψ２、ψ１は、ｘ座標の関数としてそれぞれ次のように表される。
【０１２６】
【数３０】
2008089578.tif 000032
【０１２７】
Ｙｍ軸方向ゲインを求める式に用いる係数η２、η１は、ｙ座標の関数としてそれぞれ次のように表される。
【０１２８】
【数３１】
2008089578.tif 000033
【０１２９】
このような関係式に基づいて、ゲインと回転周波数との関係を表す補正パラメータγ２、γ１、η２、η１を求める。
求められたゲイン補正のパラメータは補正パラメータ記憶部７１０に設定する。
実際の測定時における測定値の補正にあっては、上記変形例２と同様に、ゲイン補正量（｜Ｇ１Ｘ｜、｜Ｇ１Ｙ｜）と位相差φとにより、次の式で補正演算を行う。
【０１３０】
【数３２】
2008089578.tif 000034
【０１３１】
このような変形例３によってもゲイン補正量を算出するための補正パラメータを求めることができる。
そして、ゲイン補正量および位相差補正量を用いて測定データの補正を行うことができる。
【０１３２】
（第２実施形態）
次に、本発明の第２実施形態について説明する。
第２実施形態の基本的構成は第１実施形態に同様であるが、第２実施形態が特徴とするところは、動作推定部６００にて算出された加速度の推定値を実際の実加速度と対比して加速度推定の正確さを判定する推定判定部８００を備えている点にある。
図１８は、第２実施形態における動作推定部、補正演算部および推定判定部の構成を示す図である。
推定判定部８００は、カウンタ部１６１にて検出された座標の測定値に基づいて実際の加速度を算出する実加速度算出部８１０と、動作推定部６００にて算出された加速度推定値を実加速度算出部８１０にて算出された実加速度に対比して差分を算出する差分算出部８２０と、差分算出部８２０にて算出された差分値を所定の閾値に対比して加速度推定値の正確さを判定する判定部８３０と、を備えている。
【０１３３】
実加速度算出部８１０には駆動カウンタ１７１からのカウント値が入力されており、実加速度算出部８１０は駆動カウンタ１７１からのカウント値を２階微分して実際の加速度を算出する。
算出した実加速度は差分算出部８２０に出力する。
【０１３４】
差分算出部８２０には、２階微分演算部３１５からの推定加速度値と実加速度算出部８１０からの実加速度とが入力されている。
差分算出部８２０は、推定加速度ａから実加速度ａＲを減算して差分値δａを算出する。
算出した差分値δａは判定部８３０に出力する。
【０１３５】
判定部８３０には、予め閾値が設定されている。
この閾値は、推定加速度ａと実加速度ａＲとの差分値δａが大きすぎる場合に加速度推定が正しく機能しておらず、推定加速度に基づく補正は不能であることを判定するための閾値である。
判定部８３０には差分算出部８２０にて算出された差分値δａが入力され、判定部８３０はこの差分値δａを閾値と対比する。
そして、差分値δａが閾値δｓよりも小さい場合には（｜δａ｜＜δｓ）、測定動作の継続を指令する。
その一方、差分値δａが閾値δｓよりも大きい場合には（｜δａ｜≧δｓ）、その結果をユーザーに対し出力手段６２を介して表示する。
あるいは、補正演算部７００における補正演算を停止させる。
【０１３６】
このような構成を備える第２実施形態によれば、推定判定部８００を備えているので、動作推定部６００による加速度の推定の正確さが判定できる。
そして、判定結果に基づいて補正演算の継続または中止を行うことができるので、測定データの誤った補正による誤った測定結果の出力を防止することができる。
【０１３７】
なお、本発明の表面形状測定装置は、上記実施形態および変形例にのみに限定されるものではなく、本発明の要旨を逸脱しない範囲において種々変更してもよいことは勿論である。
推定動作状態量としては加速度を例にして説明したが、加速度の他、倣い速度を推定動作状態量としてもよい。この場合、倣い速度に基づいて補正演算可能なように、倣い速度と補正量との関係を求めておけばよい。
上記実施形態ではノミナルモデルから位置を推定したうえで、この推定位置を２階微分して加速度を推定するとしたが、倣いベクトルの位置指令から加速度を直接に算出するノミナルモデルを設定しておいてもよい。
第２実施形態において、推定判定部８００は内部にローパスフィルタを備え、実加速度算出部８１０にて求められた実加速度から高周波ノイズをカットするようにしてもよい。これにより高周波ノイズによる判定の誤差を防止することができる。
【図面の簡単な説明】
【０１３８】
【図１】第１実施形態に係る測定システムの機能ブロック図。
【図２】第１実施形態において、動作推定部および補正演算部の構成を示す図。
【図３】ノミナルモデル設定部の構成を示す図。
【図４】倣い速度（あるいは回転周波数）が一定値を示すところで１周分（３６０°）のデータを取る場合を示す図。
【図５】倣い速度（あるいは回転周波数）が変化（加速あるいは減速）している領域でデータを取る場合を示す図。
【図６】校正データ取得のためにＰ１からＰ５の複数の位置で円の倣い測定を行う様子を示す図。
【図７】各ポイントＰ１〜Ｐ５で円の倣い測定をした際に倣い速度Ｖと検出された半径の減少量（半径誤差ΔＲ）との関係を示す図。
【図８】各ポイントＰ１〜Ｐ５で円の倣い測定をした際に倣い速度Ｖと検出された半径の減少量（半径誤差ΔＲ）との関係を示す図。
【図９】半径の減少量と加速度との関係を示す図。
【図１０】半径の減少量と加速度との関係を示す図。
【図１１】測定位置と測定半径の減少量との関係を示す図。
【図１２】測定位置と測定半径の減少量との関係を示す図。
【図１３】Ｙｍ軸方向とＸｍ軸方向とで生じる位相差φと円運動の回転周波数ｆとの関係を示す図。
【図１４】倣い測定時の回転周波数と半径誤差ΔＲとの関係を示す図。
【図１５】倣い測定時の回転周波数と半径誤差ΔＲとの関係を示す図。
【図１６】ゲイン誤差の係数と座標との関係を示す図。
【図１７】ゲイン誤差の係数と座標との関係を示す図。
【図１８】第２実施形態における動作推定部、補正演算部および推定判定部の構成を示す図。
【図１９】背景技術において倣いプローブを利用した表面倣い測定装置としての測定システムの構成を示す図。
【図２０】倣いプローブの構成および倣いプローブにて被測定物表面を測定する様子を示す図。
【図２１】倣い測定中に変形による誤差を生じている状態を示す図。
【図２２】測定誤差として径方向のズレが生じた結果を示す図。
【符号の説明】
【０１３９】
６１…入力手段、６２…出力手段、１００…測定システム、１１０…三次元測定機、１１１…定盤、１２０…駆動機構、１２１…ビーム支持体、１２２…ビーム、１２３…コラム、１２４…スピンドル、１３０…倣いプローブ、１３１…スタイラス、１３２…接触部、１３３…支持部、１３４…プローブセンサ、１３５…Ｘｐ方向センサ、１３６…Ｙｐ方向センサ、１３７…Ｚｐ方向センサ、１４０…駆動センサ、１４１…Ｙｍ軸センサ、１４２…Ｘｍ軸センサ、１４３…Ｚｍ軸センサ、１５０…操作部、１６０…モーションコントローラ、１６１…カウンタ部、１６２…駆動制御回路、１７１…駆動カウンタ、１７２…Ｙｍ軸カウンタ、１７３…Ｘｍ軸カウンタ、１７４…Ｚｍ軸カウンタ、１７５…プローブカウンタ、１７６…Ｘｐ方向カウンタ、１７７…Ｙｐ方向カウンタ、１７８…Ｚｐ方向カウンタ、２００…ホストコンピュータ、２１０…メモリ、２２０…倣いベクトル指令部、３１０…加速度推定部、３１１…ノミナルモデル設定部、３１４…位置推定部、３１５…２階微分演算部、４３０…測定データ合成部、５００…形状解析部、６００…動作推定部、６１０…周波数推定部、７００…補正演算部、７１０…補正パラメータ記憶部、７２０…補正量算出部、７２１…ゲイン補正量算出部、７２２…位相差補正量算出部、８００…推定判定部、８１０…実加速度算出部、８２０…差分算出部、８３０…判定部。
【図１】
2008089578.tif 000035
【図２】
2008089578.tif 000036
【図３】
2008089578.tif 000037
【図４】
2008089578.tif 000038
【図５】
2008089578.tif 000039
【図６】
2008089578.tif 000040
【図７】
2008089578.tif 000041
【図８】
2008089578.tif 000042
【図９】
2008089578.tif 000043
【図１０】
2008089578.tif 000044
【図１１】
2008089578.tif 000045
【図１２】
2008089578.tif 000046
【図１３】
2008089578.tif 000047
【図１４】
2008089578.tif 000048
【図１５】
2008089578.tif 000049
【図１６】
2008089578.tif 000050
【図１７】
2008089578.tif 000051
【図１８】
2008089578.tif 000052
【図１９】
2008089578.tif 000053
【図２０】
2008089578.tif 000054
【図２１】
2008089578.tif 000055
【図２２】
2008089578.tif 000056
【手続補正書】
【提出日】平成19年9月19日(2007.9.19)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項１
【補正方法】変更
【補正の内容】
【請求項１】
少なくとも一部に円弧部分が含まれる被測定物表面に近接あるいは当接される測定子を有し、前記測定子と前記被測定物表面との相対位置を予め設定された基準位置に保って倣い走査する倣いプローブと、
前記被測定物表面に沿った前記倣いプローブの次の移動位置を指示する倣いベクトルを指令する倣いベクトル指令部と、
前記倣いプローブを保持するとともに三次元的に移動させる駆動軸を有し、前記倣いベクトル指令に応じて前記倣いプローブを移動させる駆動機構と、
前記駆動機構の駆動量を検出する駆動センサと、
前記倣いベクトル指令部にて指令された前記倣いベクトル指令に応じて前記駆動センサの検出値を補正演算する補正演算部と、を備え、
前記補正演算部は、駆動中の前記駆動機構が変形することによって生じる測定誤差を補正する補正量を前記倣いベクトル指令の指令値に基づいて算出する
ことを特徴とする表面形状測定装置。
【手続補正２】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求項２
【補正方法】変更
【補正の内容】
【請求項２】
請求項１に記載の表面形状測定装置において、
前記補正演算部は、
前記倣いプローブに生じる加速度と前記駆動機構の変形量との関係を表す補正パラメータを記憶した補正パラメータ記憶部と、
前記駆動センサの検出値を補正する補正量を前記倣いプローブに生じる加速度に基づいて算出する補正量算出部と、
前記算出された補正量を前記駆動センサの検出値に合成して前記被測定物表面の位置を求める測定データ合成部と、を備える
ことを特徴とする表面形状測定装置。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００１２
【補正方法】変更
【補正の内容】
【００１２】
本発明の表面形状測定装置は、少なくとも一部に円弧部分が含まれる被測定物表面に近接あるいは当接される測定子を有し、前記測定子と前記被測定物表面との相対位置を予め設定された基準位置に保って倣い走査する倣いプローブと、前記被測定物表面に沿った前記倣いプローブの次の移動位置を指示する倣いベクトルを指令する倣いベクトル指令部と、前記倣いプローブを保持するとともに三次元的に移動させる駆動軸を有し、前記倣いベクトル指令に応じて前記倣いプローブを移動させる駆動機構と、前記駆動機構の駆動量を検出する駆動センサと、前記倣いベクトル指令部にて指令された前記倣いベクトル指令に応じて前記駆動センサの検出値を補正演算する補正演算部と、を備え、前記補正演算部は、駆動中の前記駆動機構が変形することによって生じる測定誤差を補正する補正量を前記倣いベクトル指令の指令値に基づいて算出することを特徴とする。
なお、少なくとも一部に円弧部分が含まれるとは、一定の曲率を有する円弧を含むことを要せず、曲面を含むものであればよい。よって、自由曲面も本発明の測定対象に含まれる。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００１３
【補正方法】変更
【補正の内容】
【００１３】
このような構成において、倣いプローブを被測定物表面に沿って移動させるための倣いベクトル指令が倣いベクトル指令部から発せられる。
すると、倣いベクトル指令に応じ、駆動機構によって倣いプローブが被測定物表面に沿って倣い移動される。
倣いプローブで被測定物表面を走査するときの駆動機構の駆動量が駆動センサで検出される。
このような倣い測定の動作の際には加速度が生じ、加速度による力がかかった場合には駆動機構に変形が生じ、その変形分が測定誤差に含まれる恐れがある。
そこで、倣いプローブで倣い測定すると同時に、上記の変形分による誤差を補正する処理が行われる。
補正演算部により、倣いベクトル指令の指令値に基づいた補正量が算出される。
すなわち、動作時に生じる変形に伴う測定誤差分を補正する補正量が算出される。
例えば、倣いベクトル指令の指令値から導出される加速度に基づいて駆動機構の変動を見積り、補正量を算出する。
算出された補正量は、駆動センサによる検出値と合成されて測定データとして出力される。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００１５
【補正方法】変更
【補正の内容】
【００１５】
本発明では、前記補正演算部は、前記倣いプローブに生じる加速度と前記駆動機構の変形量との関係を表す補正パラメータを記憶した補正パラメータ記憶部と、前記駆動センサの検出値を補正する補正量を前記倣いプローブに生じる加速度に基づいて算出する補正量算出部と、前記算出された補正量を前記駆動センサの検出値に合成して前記被測定物表面の位置を求める測定データ合成部と、を備えることが好ましい。
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
このような構成において、補正演算に必要なパラメータは補正パラメータ記憶部に予め記憶されている。
そして、補正量算出部では、倣いベクトル指令の指令値から導出される加速度および補正パラメータに基づいて補正量を算出する。
算出された補正量が駆動センサの検出値に合成されることにより正確な測定値が求められる。
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】００２３
【補正方法】変更
【補正の内容】
【００２３】
このような構成において、倣いベクトル指令部からの倣いベクトル指令が位置推定部に入力される。
すると、位置推定部において倣いベクトル指令による位置指令Ｃと設定されたノミナルモデルＧＮとが乗算され、倣いプローブの推定位置Ｅが算出される。
算出された倣いプローブの推定位置Ｅは２階微分演算部に出力される。
そして、２階微分演算部において、倣いプローブの推定位置から推定動作状態量として加速度データが算出される。
この加速度データは、補正量算出部に出力され、駆動動作によって生じる駆動機構の変形に伴う測定誤差分を補正する補正量が算出される。
算出された補正量は、駆動センサによる検出値と合成され、測定データとして出力される。