"回転する工具を主軸に装着して工作物を加工するときに発生するびびり振動を回避する工具回転速度選定方法であって、基準工具を装着した主軸系の固有振動数の解析値である解析基準固有振動数を計算する基準固有振動数解析工程と、前記基準工具を主軸に装着して測定した固有振動数である実測固有振動数を測定する固有振動数実測工程と、所定固有振動数について、前記解析基準固有振動数と前記実測固有振動数の差である補正値を演算する補正値演算工程と、前記基準工具を除く工具を装着した主軸系の固有振動数の解析値である解析固有振動数を計算する固有振動数解析工程と、前記解析固有振動数を前記補正値を用いて補正して予測固有振動数を演算する固有振動数補正工程と、前記予測固有振動数を用いて推奨工具回転速度を演算する、工具回転速度演算工程と、を備える工具回転速度選定方法。前記基準固有振動数解析工程前記固有振動数実測工程がコンプライアンスの周波数応答のピーク値から固有振動数を検出し、前記補正値が、固有振動数の周波数の差である周波数補正値と、コンプライアンスの差であるコンプライアンス補正値の組であり、前記所定固有振動数が所定の閾値を越えるコンプライアンスを備える固有振動数である、請求項1に記載の工具回転速度選定方法。前記基準工具として形状パターンと寸法の異なる複数の基準工具を備える、請求項1または請求項2に記載の工具回転速度選定方法。工具の動剛性が小さくて質量の小さい小型の工具の場合は工具の特性値が支配的で固有振動数が決まるため、従来技術1の補正係数による固有振動数の計算でも誤差は少ない。しかし、工具単体の動剛性が大きい場合や質量が大きい場合は、固有振動数は工具と主軸系を含む複数の質量、ばね、減衰からなる連成系の振動モデルとして決定されるので、複数の固有振動数を持つ。従来技術1では工具のみの特性値から1個の固有振動数を算出している、このため、複数の固有振動数を持つ場合は、どの固有振動数でびびり振動が発生するか不明である。さらに、正確な固有振動数の予測ができずびびり振動の防止回転速度に誤差を生じる恐れがある。本発明は上記事情に鑑みてなされたものであり、大型の工具を用いた場合にも、正確な固有振動数の推定を可能にして、びびり振動を回避する工具回転速度選定方法を提供することを目的とする。上記の課題を解決するための請求項1に係る発明の特徴は、回転する工具を主軸に装着して工作物を加工するときに発生するびびり振動を回避する工具回転速度選定方法であって、基準工具を装着した主軸系の固有振動数の解析値である解析基準固有振動数を計算する基準固有振動数解析工程と、前記基準工具を主軸に装着して測定した固有振動数である実測固有振動数を測定する固有振動数実測工程と、所定固有振動数について、前記解析基準固有振動数と前記実測固有振動数の差である補正値を演算する補正値演算工程と、前記基準工具を除く工具を装着した主軸系の固有振動数の解析値である解析固有振動数を計算する固有振動数解析工程と、前記解析固有振動数を前記補正値を用いて補正して予測固有振動数を演算する固有振動数補正工程と、前記予測固有振動数を用いて推奨工具回転速度を演算する、工具回転速度演算工程と、を備えることである。請求項2に係る発明の特徴は、請求項1に係る発明において、前記基準固有振動数解析工程前記固有振動数実測工程がコンプライアンスの周波数応答のピーク値から固有振動数を検出し、前記補正値が、固有振動数の周波数の差である周波数補正値と、コンプライアンスの差であるコンプライアンス補正値の組であり、前記所定固有振動数が所定の閾値を越えるコンプライアンスを備える固有振動数であることである。請求項3に係る発明の特徴は、請求項1または請求項2に係る発明において、前記基準工具として形状パターンと寸法の異なる複数の基準工具を備えることである。請求項1に係る発明によれば、基準工具を装着した主軸系の固有振動数の解析値と実測値の差である補正値を用いて、使用工具を装着した主軸系の固有振動数の解析値を補正するので、使用工具を装着した主軸系の正確な固有振動値が算出できる。そのため、びびり振動の発生しない推奨工具軸回転速度を正確に選定できる。請求項2に係る発明によれば、固有振動数のコンプライアンス値を用いて、複数の固有振動数についてびびり振動発生の可能性の順位を判定できる。このため、びびり振動発生の可能性の高い順番に推奨工具軸回転速度を選定できる。請求項3に係る発明によれば、形状パターンと寸法の異なる複数の基準工具を備えることで、使用工具に近い特性を持つ基準工具の補正値を用いることができ、より正確な使用工具の固有振動値が算出できる。そのため、びびり振動の発生しない推奨工具軸回転速度をより正確に選定できる。以下、本発明の実施形態を図1~図7に基づき説明する。図1は工具2を装着した主軸1の概略と、固有振動数測定装置3を用いてインパルス応答試験により固有振動数を測定する概念を示した図である。主軸1は、工具2を保持するスピンドル11と、スピンドル11を回転自在に軸受12、13を介して支持するハウジング16をそなえる。スピンドル11の軸方向の中央部にはモータロータ14が固定され、ハウジング16にモータステータ15が固定されている。インパルス応答試験は、スピンドル11を静止した状態で、センサ32を固定した工具2をハンマー31で打撃し、そのときに出力される信号を固有振動数演算装置33で演算して測定する。図2は主軸工具系のコンプライアンスの周波数応答を示す図で、横軸が周波数、縦軸がコンプライアンスである。固有振動数はコンプライアンスが図2のa、b、c、d等のピークを示す位置である。再生びびり振動は、固有振動数近辺で切屑厚みの周期的変動が発生して、切削力の変動が増大することにより発生する。このため、工具の刃先の作用周期を固有振動数に一致させて切屑厚みの変動を無くすれば防止できることが知られている。このため、固有振動数をfHz、工具刃数をe、1以上の整数をn、工具回転速度をNmin-1とすると、N=f60/e.nで算出される工具回転速度で加工するとびびり振動が発生しないことが知られている。主軸系は振動モデルとしてみた場合、工具2スピンドル11モータロータ14の質量、工具2スピンドル11軸受12軸受13のばねなどで構成されている。このため、複数の質量とばねが連成された振動系として複数の固有振動数が存在する。この固有振動数の周波数とコンプライアンスを、質量、ばね、減衰係数を定めてモーダル計算することで解析的に求めることができることが知られている。複数の固有振動数を持つ系のびびり振動は、コンプライアンスの大きな固有振動数で発生する可能性が高いが、必ずしも順位通りに決まるものではなく、コンプライアンスの上位の固有振動数の中の1つで発生することが多い。びびり振動を回避する推奨工具回転速度選定方法について、図3の工程図に基づき説明する。初めに、基準工具2をスピンドル11に装着した状態でインパルス応答試験を行い、図2の実線に示すようなコンプライアンスの周波数応答線図を測定するS1。平行して、基準工具2をスピンドル11に装着した系の解析を行い、図2の破線に示すようなコンプライアンスの周波数応答線図を演算するS2。基準工具を装着した系の実測値と解析値の対応ピークを決定する。具体的には、図2において、コンプライアンスの実測値と解析値の最大値同士を対応するピーク対aとする、ピーク対aから周波数の降順または昇順で同じ順番のピークを対応するピーク対b、c,dとする。対象とするピークは実測値において所定の閾値より大きい値の4個のピークとする。S3。基準工具を装着した系の補正値を演算する。図4にピークbの具体例を示す、ピークbの解析による周波数をMfkb、実測による周波数をMfbとすると、周波数補正値ΔfbはΔfb=Mfb-Mfkbとなる。ピークbの解析によるコンプライアンスをMCkb、実測によるコンプライアンスをMCbとすると、コンプライアンス補正値ΔCbはΔCb=MCb-MCkbとなる。同様にして夫々の対応ピーク対毎に、周波数補正値Δfa~Δfdと、コンプライアンス補正値ΔCa~ΔCdを演算するS4。使用工具をスピンドルに装着した状態の解析を行い、周波数応答線図を演算するS5。基準工具を装着した系の解析値と使用工具を装着した系の解析値の対応ピークを決定する。具体的には、図5に示すように、S3で決めた基準工具を装着した系の解析値のピークaに使用工具を装着した系の解析値の最大ピークを対応させ、ピーク対aから周波数の降順または昇順で同じ順番のピークを対応するピーク対b、c,dとするS6。使用工具を装着した系の固有振動数の周波数とコンプライアンスの予測値を演算する。具体的には、図6に示すように、使用工具を装着した系の解析のピーク値a、b、c,dの周波数とコンプライアンスの値に、夫々対応する周波数補正値Δfa~Δfdと、コンプライアンス補正値ΔCa~ΔCdを加算し、周波数予測値fha~fhdとコンプライアンス予測値Cha~Chdを演算するS7。推奨順位と推奨工具軸回転速度の決定。推奨順位はコンプライアンス予測値の大きい順とする。順位1番目の推奨工具軸回転速度N1はN1=fha60/e.nとなる。ここで、eは工具の刃数でnは1以上の整数である。S8。推奨工具軸回転速度とその推奨順位を出力するS9。以上のように、本発明によれば基準工具の解析値と実測値差から求めた補正値を用いて、使用工具の解析値を補正して求めた予測固有振動数の周波数から、びびり振動を発生しない推奨工具軸回転速度を算出するので、実測定をすることなく、正確な推奨工具軸回転速度を設定できる。さらに、複数の固有振動数のコンプライアンス順位に基づき、複数の推奨工具軸回転速度を推奨順位と共に出力するので、複数の固有振動数をもつ主軸系に対しても確実にびびり振動を発生しない推奨工具軸回転速度を算出できる。上記の説明では、1つの基準工具の補正値を用いた例について述べたが、図7に示すように複数の基準工具を準備しておき、使用工具の形状が2つの基準工具の中間に位置する場合は、2つの基準工具の補正値の平均値を用いて補正してもよい。"