回転軸に装着した工具又はワークを回転させて加工を行う工作機械において、前記回転軸の回転を制御するための回転軸制御方法であって、前記回転軸の初期回転速度及びその初期回転速度からの変更パターンをそれぞれ設定するパラメータ設定ステップと、前記初期回転速度で前記回転軸を回転させて加工を行い、加工中に発生する振動を検出する振動検出ステップと、検出された前記振動を解析して強制びびり振動か否かを判別する振動判別ステップと、前記振動が強制びびり振動であれば、前記初期回転速度から前記変更パターンに基づいて前記回転速度を変更する回転速度変更ステップと、前記回転速度の変更後に前記振動の検出及び解析を行い、前記振動が強制びびり振動であれば、前記強制びびり振動の振動量が最大となったときの振動周波数に基づいて共振周波数を特定する共振周波数特定ステップと、特定された前記共振周波数を回避する前記回転速度を算出して当該回転速度に変更する回転速度算出/変更ステップとを実行することを特徴とする工作機械の回転軸制御方法。前記パラメータ設定ステップでは、前記回転速度の上下限をそれぞれ設定し、前記共振周波数特定ステップでは、前記回転速度の上下限内での異なるハーモニクス線上で前記振動量が最大の振動周波数となる前記回転速度があるか否かを算出し、当該回転速度がある場合には前記ハーモニクス線上での前記振動周波数に基づいて共振周波数を特定し、前記回転速度算出/変更ステップでは、各ハーモニクス線上で特定された複数の前記共振周波数を回避する前記回転速度を算出して当該回転速度に変更することを特徴とする請求項1に記載の工作機械の回転軸制御方法。回転軸に装着した工具又はワークを回転させて加工を行う工作機械において、前記回転軸の回転を制御するための回転軸制御装置であって、前記回転軸の初期回転速度及びその初期回転速度からの変更パターンをそれぞれ設定するパラメータ設定手段と、前記回転軸を回転させて加工を行った際に発生する振動を検出する振動検出手段と、検出された前記振動を解析して強制びびり振動か否かを判別する振動判別手段と、前記振動が強制びびり振動であれば、前記初期回転速度から前記変更パターンに基づいて前記回転速度を変更する回転速度変更手段と、前記回転速度の変更後に検出及び解析された前記振動が強制びびり振動であれば、前記強制びびり振動の振動量が最大となったときの振動周波数に基づいて共振周波数を特定する共振周波数特定手段と、特定された前記共振周波数を回避する前記回転速度を算出して当該回転速度に変更する回転速度算出/変更手段とを備えることを特徴とする工作機械の回転軸制御装置。前記パラメータ設定手段は、前記回転速度の上下限をそれぞれ設定し、前記共振周波数特定手段は、前記回転速度の上下限内での異なるハーモニクス線上で前記振動量が最大の振動周波数となる前記回転速度があるか否かを算出し、当該回転速度がある場合には前記ハーモニクス線上での前記振動周波数に基づいて共振周波数を特定し、前記回転速度算出/変更手段は、各ハーモニクス線上で特定された複数の前記共振周波数を回避する前記回転速度を算出して当該回転速度に変更することを特徴とする請求項3に記載の工作機械の回転軸制御装置。しかし、特許文献1の方法では、シミュレーション精度が低いと実物との固有振動数の誤差が大きくなる上、加工によっては静止状態と加工状態とでは異なる特性を持つ場合があるため、事前に測定を行った場合にも誤った加工条件を提案してしまう危険がある。また、特許文献2、3の方法では、演算精度が低いと正確な分析が行えず、びびり振動の種類を正しく判別できずに誤った処置を行ってしまう危険がある一方、高い演算精度を得ようとすると、処理時間が長くなったり、演算負荷が大きくなったりして、十分なびびり振動の低減効果が得られない可能性がある。そこで、本発明は上記問題に鑑みなされたものであって、実際の加工状態から情報を得ることで精度の高い分析を行うとともに、若干の分析誤差を含んだ場合においても、びびり振動の種類を判別可能にし、発生したびびり振動に対して適切な対処を行い、びびり振動の低減が期待できる回転軸制御方法及び装置を提供することを目的としたものである。上記目的を達成するために、請求項1に記載の発明は、回転軸に装着した工具又はワークを回転させて加工を行う工作機械において、前記回転軸の回転を制御するための回転軸制御方法であって、前記回転軸の初期回転速度及びその初期回転速度からの変更パターンをそれぞれ設定するパラメータ設定ステップと、前記初期回転速度で前記回転軸を回転させて加工を行い、加工中に発生する振動を検出する振動検出ステップと、検出された前記振動を解析して強制びびり振動か否かを判別する振動判別ステップと、前記振動が強制びびり振動であれば、前記初期回転速度から前記変更パターンに基づいて前記回転速度を変更する回転速度変更ステップと、前記回転速度の変更後に前記振動の検出及び解析を行い、前記振動が強制びびり振動であれば、前記強制びびり振動の振動量が最大となったときの振動周波数に基づいて共振周波数を特定する共振周波数特定ステップと、特定された前記共振周波数を回避する前記回転速度を算出して当該回転速度に変更する回転速度算出/変更ステップとを実行することを特徴とするものである。請求項2に記載の発明は、請求項1の構成において、前記パラメータ設定ステップでは、前記回転速度の上下限をそれぞれ設定し、前記共振周波数特定ステップでは、前記回転速度の上下限内での異なるハーモニクス線上で前記振動量が最大の振動周波数となる前記回転速度があるか否かを算出し、当該回転速度がある場合には前記ハーモニクス線上での前記振動周波数に基づいて共振周波数を特定し、前記回転速度算出/変更ステップでは、各ハーモニクス線上で特定された複数の前記共振周波数を回避する前記回転速度を算出して当該回転速度に変更することを特徴とするものである。上記目的を達成するために、請求項3に記載の発明は、回転軸に装着した工具又はワークを回転させて加工を行う工作機械において、前記回転軸の回転を制御するための回転軸制御装置であって、前記回転軸の初期回転速度及びその初期回転速度からの変更パターンをそれぞれ設定するパラメータ設定手段と、前記回転軸を回転させて加工を行った際に発生する振動を検出する振動検出手段と、検出された前記振動を解析して強制びびり振動か否かを判別する振動判別手段と、前記振動が強制びびり振動であれば、前記初期回転速度から前記変更パターンに基づいて前記回転速度を変更する回転速度変更手段と、前記回転速度の変更後に検出及び解析された前記振動が強制びびり振動であれば、前記強制びびり振動の振動量が最大となったときの振動周波数に基づいて共振周波数を特定する共振周波数特定手段と、特定された前記共振周波数を回避する前記回転速度を算出して当該回転速度に変更する回転速度算出/変更手段とを備えることを特徴とするものである。請求項4に記載の発明は、請求項3の構成において、前記パラメータ設定手段は、前記回転速度の上下限をそれぞれ設定し、前記共振周波数特定手段は、前記回転速度の上下限内での異なるハーモニクス線上で前記振動量が最大の振動周波数となる前記回転速度があるか否かを算出し、当該回転速度がある場合には前記ハーモニクス線上での前記振動周波数に基づいて共振周波数を特定し、前記回転速度算出/変更手段は、各ハーモニクス線上で特定された複数の前記共振周波数を回避する前記回転速度を算出して当該回転速度に変更することを特徴とするものである。請求項1及び3に記載の発明によれば、実際の加工状態から効率良く複数の情報を得るため、計測や分析に誤差を含んだ場合においても、迅速且つ高精度にびびり振動の種類の判別を行い、適切な対処を行うことができる。よって、加工振動の低減が期待できる。請求項2及び4に記載の発明によれば、上記効果に加えて、強制びびり振動が発生する可能性が高い共振周波数を確実に回避して最適な回転速度を選択することができる。以下、本発明の実施の形態を図面に基づいて説明する。図1は、工作機械の一例である立形マシニングセンタの概略構成図で、立形マシニングセンタ1は、上方に設けた主軸頭2に回転軸としての主軸3を設け、その主軸3に取り付けた工具4によって、下方の加工テーブル5上にセットされたワーク6を加工する周知の構成で、図示しない自動工具交換装置によって工具4を自動交換可能となっている。また、立形マシニングセンタ1は、回転軸制御装置としての主軸制御装置10を具備している。この主軸制御装置10は、主軸頭2に設けられて振動を検出する振動検出手段としての振動センサここでは加速度ピックアップ7と、振動センサ7による検出値に基づいて主軸3の回転を制御する制御装置11とを備えている。まず振動センサ7は、主軸3の回転に伴って生じる時間領域の振動時間軸上の振動を検出するもので、互いに直交するX軸、Y軸、Z軸方向での時間領域の振動情報を検出可能な状態で主軸頭2に取り付けられている。制御装置11は、主軸3の回転速度の上下限と回転速度の変更パターンここでは回転速度の変更範囲及び変更回数、主軸制御に係る閾値や工具刃数を入力するパラメータ設定手段としての入力装置12と、振動センサ7により検出された振動加速度を基にフーリエ解析を行い、びびり振動が発生しているかの判断やびびり振動の種類の判別、回転速度の算出等を行う演算装置13と、主軸3の回転速度と振動センサ7から送信された測定結果及び演算装置13における演算結果とを関連付けて記憶する記憶装置14と、演算装置13からの指令を受けて主軸3の回転を制御するNC装置15とを備えている。ここでは演算装置13が振動判別手段及び共振周波数特定手段として機能し、演算装置13及びNC装置15が回転速度変更手段及び回転速度算出/変更手段として機能することになる。以上の如く構成された主軸制御装置10は、加工中にびびり振動が発生した場合において、強制びびり振動の可能性があると判断された場合に、振動量が閾値を下回るか、複数の回転速度においても強制びびり振動と考えられる挙動を示すかを監視しながら、得られた測定結果から安定な回転速度を算出し、その回転速度へ変更することで加工中に発生するびびり振動を低減する。以下、この主軸3の制御を図2のフローチャートに基づいて説明する。まず、S1で、加工開始前に入力装置12により、演算装置13で用いる振動量の閾値と、工具刃数と、回転速度の上下限と、回転速度の変更範囲例えば初期回転速度の±数の範囲内で小側から大側へ所定速度ずつ段階的に変更する等及び変更回数をそれぞれ設定するパラメータ設定ステップ。S2で加工を開始し、加工中は、演算装置13が振動センサ7で測定した振動加速度をフーリエ解析して、振動量ここでは最大加速度を算出し、振動量と閾値とを比較するS3、振動検出ステップ。振動量が閾値より小さければ制御を終了する。一方、S3の判別で振動量が閾値以上であれば、びびり振動が発生したと判断して、S4で、演算装置13は発生したびびり振動が強制びびり振動か否か再生びびり振動かを判定する振動判別ステップ。再生びびり振動と判定した場合、S5で最適回転速度へ変更してS3へ戻る。一方、S4で強制びびり振動と判定した場合、演算装置13は、S6でそのときの回転速度と振動量と振動周波数とを記憶装置14に記憶し、S7で、設定した変更パターンに従ってNC装置15を介して回転速度を変更して、S8で再び振動量と閾値との比較を行う。ここで振動量が閾値より小さくなれば制御を終了し、振動量が閾値以上であれば、S9で回転速度の変更回数が設定した回数に達したか否かを判別する。ここで設定回数に達していなければ、S6に戻って再び回転速度等を記憶してS7で回転速度の変更を行い、S8で振動量と閾値との比較を行う。このS6~S8の処理は、振動量が閾値を下回らない場合はS9の判別で変更回数が設定回数に達するまで繰り返される回転速度変更ステップ。そして、振動量が閾値を下回らないまま回転速度の変更が設定回数に達したらS9でYES、演算装置13は、S10で、S6で記憶されたデータの中で最も振動量の大きかったときの振動周波数をfmaxとして算出し、S11では、設定した回転速度の上下限の範囲内において、異なるハーモニクス線上でfmaxと成り得る回転速度が存在するか否かを判別する。ここで他のハーモニクス線上でも当該回転速度が存在する場合には、S3に戻って当該ハーモニクス線上でも回転速度を変更しながら振動量と閾値との比較を行い、fmaxを算出する。そして、S11の判別で、異なるハーモニクス線上でfmaxとなり得る回転速度が存在しないと確認されると、演算装置13は、S12で、それまでの処理で記憶した測定データに対して、S4と同様に各ハーモニクス線上での各回転速度で発生したびびり振動が強制びびり振動か再生びびり振動かを判定する。ここで、強制びびり振動と再生びびり振動が混在していると判定された場合、S13で、記憶された測定データの中で振動量が最も小さい回転速度へ変更し、制御を終了する。一方、S12の判別で強制びびり振動のみが発生していると判定された場合、S14で、各ハーモニクス線上で算出されたfmaxの値同士を比較して、一定の周波数差以内のfmaxを同一の共振モードとして分別し、分別された各共振モードにてfmaxの値の平均を共振周波数fC1、fC2として特定する共振周波数特定ステップ。次に、S15で、それらの共振周波数を回避するような入力となる回転速度Saを算出し、S16でその回転速度Saに変更して制御を終了する回転速度算出/変更ステップ。ここで上記手順において、S4及びS12での強制びびり振動か否かの判定では、例えば特許文献2に含まれる以下の式1~4の位相情報を用いて、同特許文献と同様の方法で判定することができる。具体的には、式3で得られた位相情報が0に近い値、例えば0とした定数1より大きく、且つ0.1とした定数2より小さければ、強制びびり振動が発生していると判断でき、それ以外の場合には再生びびり振動が発生していると判断できる。この場合、S5の最適回転速度への変更では、例えば式4の計算式を用いて算出された回転速度を用いることができる。ここで用いるk1値は、例えばk1値=k値1によって算出すればよい。この回転速度に変更することで、再生びびり振動を最も抑制できる。k’値=60×びびり振動数/工具刃数×主軸回転速度1k値=k’値の整数部2位相情報=k’値-k値3最適回転速度=60×びびり振動数/工具刃数×k1値4次に、S11の判別について、図3のキャンベル線図をもとに説明する。まず、入力となる各ハーモニクス線fkは、以下の式5で示されるSの一次関数として算出できる。fk=S×Z×k/60k:整数値5ここで、Sは回転速度、Zは工具刃数である。図3のS31のように、この値が共振周波数帯f1と一致すると強制びびり振動が発生する。ここで、図3のように上限回転速度Smaxと下限回転速度Sminが設定されていた場合、その範囲内において、図3のS32のように、k=3のハーモニクス線上でも共振周波数帯f1と一致して強制びびり振動が発生すると判断されるため、異なるハーモニクス線上でfmaxと成り得る回転速度が存在すると判断される。そして、異なるハーモニクス線上で試行した回転速度においてf2での振動量が大きかった場合には、f2の共振周波数帯と一致する回転速度においても上記f1と同様の動作を行う。さらに異なる共振周波数帯で振動量が大きかった場合にも同様である。また、S15での回転速度の算出、S16での変更について、図3のキャンベル線図をもとに説明する。S14までの動作に基づき、例えばf1及びf2の共振周波数帯が確認されたとする。このとき、式5から、どちらの共振周波数帯とも一致しない回転速度例えば図3の回転速度Saが算出されるため、当該回転速度Saに変更することで強制びびり振動の発生しない安定な加工を行うことができる。以下、図4をもとに制御事例を示す。まず、S1で、工具刃数Zが4枚の工具、初期回転速度1400min-1、上限回転速度2500min-1、下限回転速度1000min-1で各パラメータをそれぞれ設定し、S2で加工を開始した場合、S3で187Hzで振動量が閾値を超えて最も大きくなった。S4で、式3を用いて強制びびり振動と判定されたため、S6~S9において1400min-1から±10となる1260min-1から1540min-1で段階的に回転速度を変更して加工を行ったところ、1500min-1時に200Hzで最も振動が大きくなったS10。そのため、200Hzに共振周波数が存在すると判断でき、その±10の領域が共振周波数帯、すなわち180Hz~220Hzにおいて強制びびり振動が発生する領域であると算出される。次にS11では、以下の手順で判別を行う。まず、式5に各値を代入して変形すると、以下の式6になる。S=fk×15/k6ここで、fkに上記180Hz~220Hzを代入すると、以下の式7になり、その回転速度Sでは強制びびり振動が発生すると予測される。S=2700/k~3300/k7式6において、200Hz、k=3でS=1000となり、下限回転速度1000min-1内に含まれる範囲を有するため、S11の判別がYESとなる。よって、S3に戻って続くS6~S9では、1000min-1~1100min-1で回転速度を変更して加工を行う。ここで、例えば1000min-1~1100min-1において、600Hz~660Hzの振動が発生し、S10において1000min-1時に600Hzで最も振動が大きくなったとする。そのため、600Hzに共振周波数が存在すると判断でき、その±10の領域を共振周波数帯、すなわち540Hz~660Hzにおいて強制びびり振動が発生する領域であると算出される。よって、上述の手順と同様に計算すると、以下の式8になり、その回転速度Sでは強制びびり振動が発生することになる。S=8100/k~9900/k8式6において、600Hz、k=4~9において、上限回転速度2500min-1、下限回転速度1000min-1に含まれる範囲を有するため、S11の判別がYESとなる。よって、S3に戻ってこれらの回転速度においても同様の加工を行う。そして、このように各ハーモニクス線上でそれぞれ回転速度を変更した加工を行ったとき、S12で依然として各ハーモニクス線上で強制びびり振動が発生していたら、S14において、200Hzの共振周波数帯、及び600Hzの共振周波数帯をそれぞれ特定し、S15では、入力が各共振周波数帯に一致する回転速度を算出してそれ以外の回転速度Saを算出することで強制びびり振動を回避することができる。今回の場合、例えば2000min-1がそれに当たるため、S16ではその回転速度へ変更することで強制びびり振動が回避できる。なお、回転速度変更中に閾値を切った場合には、その回転速度で制御を終了する。このように、上記形態の主軸制御装置10によれば、実際の加工状態から効率良く複数の情報を得るため、計測や分析に誤差を含んだ場合においても、迅速且つ高精度にびびり振動の種類の判別を行い、適切な対処を行うことができる。よって、加工振動の低減が期待できる。特にここでは、異なるハーモニクス線上でも振動量が最大の振動周波数となる回転速度があれば、当該ハーモニクス線上でも共振周波数を特定して、各ハーモニクス線上で特定された複数の共振周波数を回避する回転速度を算出して変更するようにしているので、強制びびり振動が発生する可能性が高い共振周波数を確実に回避して最適な回転速度を選択することができる。なお、本発明の回転軸制御装置に係る構成は、上記実施の形態に記載した態様に何ら限定されるものではなく、本発明の趣旨を逸脱しない範囲で必要に応じて適宜変更することができる。たとえば、振動の検出には加速度ピックアップの代わりに変位計やマイクを採用して、振動加速度の代わりに振動変位や振動音を用いることも可能である。また、上記制御事例では共振周波数帯を共振周波数の±10としているが、異なる割合に変えたり、変更可能にしたり、回転速度の絶対値での範囲としたりすることも可能である。さらに、上記制御事例では異なるハーモニクス線上でも複数の回転速度での加工を試みているが、制御時間短縮のために特定された共振周波数となる回転速度のみで加工を行うことも可能である。一方、上記形態では回転速度の制御を自動で実行しているが、自動制御の代わりに手動で対話式に各動作を行ってもよいことは言うまでもない。その他、工作機械としては立形マシニングセンタに限らず、主軸に装着したワークを回転させて加工を行うNC旋盤等の他の工作機械であっても本発明は適用可能である。