(57)【特許請求の範囲】
【請求項１】
工具を回転して被削材を加工する機械を制御する数値制御装置において、
前記被削材の実加工に対応するテスト加工を行い、該テスト加工中において所定時間毎に前記工具の回転数と前記機械の外乱トルクを夫々測定する測定手段と、
前記測定手段が測定した前記工具の回転数と前記外乱トルクの情報である測定情報に基づき、前記工具にびびり振動が生じにくい前記工具の回転数である安定回転数を確定する確定手段と、
前記確定手段が確定した前記安定回転数の情報である安定回転数情報を出力する出力手段と
を備え、
前記テスト加工は、
連続的に前記工具の回転数を増加又は減少しながら、前記被削材と同一材質のテスト用被削材を加工する条件で行う加工であって、
前記確定手段は、
前記測定情報のうち、前記外乱トルクの周波数解析を行ってパワースペクトル値を解析する解析手段と、
前記解析手段が解析した前記パワースペクトル値に基づき、前記所定時間毎の最大パワースペクトル値を算出する算出手段と、
前記算出手段が算出した前記最大パワースペクトル値の時間変化の中で、傾きが負から正に変化する変化時間に対応する前記工具の回転数を、前記安定回転数として確定する回転数確定手段と
を備え、
前記出力手段は、
前記回転数確定手段が確定した前記安定回転数の情報を、前記安定回転数情報として出力すること
を特徴とする数値制御装置。
【請求項２】
前記確定手段は、
前記算出手段が算出した前記所定時間毎の前記最大パワースペクトル値の夫々に対応するピーク周波数を演算する演算手段を備え、
前記回転数確定手段は、
前記演算手段が演算した前記ピーク周波数の時間変化の中で、閾値よりも大きい変化がある前記変化時間に対応する前記回転数を、前記安定回転数として確定すること
を特徴とする請求項１に記載の数値制御装置。
【請求項３】
前記工具の一刃当たりの送り量は、前記工具の送り速度を、前記工具の刃数と前記回転数で除して得られる値であって、
前記テスト加工では、
前記工具の一刃当たりの送り量が一定となるように、前記工具の送り速度を変化させること
を特徴とする請求項１又は２に記載の数値制御装置。
【請求項４】
前記解析手段は、
前記所定時間毎の前記外乱トルクを、周波数解析処理を用いて、前記パワースペクトル値を解析すること
を特徴とする請求項１から３の何れか一つに記載の数値制御装置。
【請求項５】
前記測定手段は、
前記機械の駆動軸の位置情報とトルクに基づき前記外乱トルクを推定する外乱オブザーバを備えたこと
を特徴とする請求項１から４の何れか一つに記載の数値制御装置。
【請求項６】
前記テスト加工は、
連続的に前記工具の回転数を減少しながら、前記被削材と同一材質のテスト用被削材を加工する条件で行う加工であること
を特徴とする請求項１から５の何れか一つに記載の数値制御装置。
【請求項７】
工具を回転して被削材を加工する機械を制御する数値制御装置の制御方法において、
前記被削材の実加工に対応するテスト加工を行い、該テスト加工中において所定時間毎に前記工具の回転数と前記機械の外乱トルクを夫々測定する測定工程と、
前記測定工程で測定した前記工具の回転数と前記外乱トルクの情報である測定情報に基づき、前記工具にびびり振動が生じにくい前記工具の回転数である安定回転数を確定する確定工程と、
前記確定工程で確定した前記安定回転数の情報である安定回転数情報を出力する出力工程と
を備え、
前記テスト加工は、
連続的に前記工具の回転数を増加又は減少しながら、前記被削材と同一材質のテスト用被削材を加工する条件で行う加工であって、
前記確定工程は、
前記測定情報のうち、前記外乱トルクの周波数解析を行ってパワースペクトル値を解析する解析工程と、
前記解析工程で解析した前記パワースペクトル値に基づき、前記所定時間毎の最大パワースペクトル値を算出する算出工程と、
前記算出工程で算出した前記最大パワースペクトル値の時間変化の中で、傾きが負から正に変化する変化時間に対応する前記工具の回転数を、前記安定回転数として確定する回転数確定工程と
を備え、
前記出力工程は、
前記回転数確定工程で確定した前記安定回転数の情報を、前記安定回転数情報として出力すること
を特徴とする制御方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、数値制御装置と制御方法に関する。
【背景技術】
【０００２】
数値制御装置が制御する工作機械は、剛性の低い工具を使用した時、又は荒加工のような重切削を行う時、びびり振動を発生する。びびり振動は工具と被削材間の振動である。びびり振動は、加工精度の低下、工具寿命の短縮等に繋がるので望ましくない。特許文献１が開示する振動検出装置は、加工中の推定外乱力をデジタルフィルタ処理することで、びびり振動の発生と発生したびびり振動の種類を検出する。特許文献２が開示する回転速度表示装置は、加工中に発生したびびり振動を解析する。解析した結果、強制びびり振動であれば工具接触周波数からびびり振動が生じやすい周波数を計算する。再生びびり振動であればびびり振動の位相情報から系の固有振動数と減衰比を計算で求めて安定限界線図を作成することで、作業者に望ましい加工条件を提示する。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１２−２０６２３０号公報
【特許文献２】特開２０１２−１９６７４１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１の振動検出装置は、どの加工条件が最適となるかについてテスト加工を繰り返さないと分からないという問題点があった。特許文献２の回転速度表示装置は、現在最も強く発生している一つの周波数の振動についてのみ解析を行う。それ故、複数の振動モードが影響（モードカップリング）する再生びびり振動については、安定限界線図を正しく推定できないという問題点があった。
【０００５】
本発明の目的は、実加工時にびびり振動が生じにくい工具回転数を求めることができる数値制御装置と制御方法を提供することにある。
【課題を解決するための手段】
【０００６】
本発明の請求項１に係る数値制御装置は、工具を回転して被削材を加工する機械を制御する数値制御装置において、前記被削材の実加工に対応するテスト加工を行い、該テスト加工中において所定時間毎に前記工具の回転数と前記機械の外乱トルクを夫々測定する測定手段と、前記測定手段が測定した前記工具の回転数と前記外乱トルクの情報である測定情報に基づき、前記工具にびびり振動が生じにくい前記工具の回転数である安定回転数を確定する確定手段と、前記確定手段が確定した前記安定回転数の情報である安定回転数情報を出力する出力手段とを備え、前記テスト加工は、連続的に前記工具の回転数を増加又は減少しながら、前記被削材と同一材質のテスト用被削材を加工する条件で行う加工であって、前記確定手段は、前記測定情報のうち、前記外乱トルクの周波数解析を行ってパワースペクトル値を解析する解析手段と、前記解析手段が解析した前記パワースペクトル値に基づき、前記所定時間毎の最大パワースペクトル値を算出する算出手段と、前記算出手段が算出した前記最大パワースペクトル値の時間変化の中で、傾きが負から正に変化する変化時間に対応する前記工具の回転数を、前記安定回転数として確定する回転数確定手段とを備え、前記出力手段は、前記回転数確定手段が確定した前記安定回転数の情報を、前記安定回転数情報として出力することを特徴とする。数値制御装置は実加工に対応するテスト加工を一回行うことで、実加工でびびり振動が生じにくい安定回転数を確定できる。作業者は出力手段が出力した安定回転数を参照し、他条件を加味した上で最適な工具回転数を決定できる。それ故、数値制御装置は実加工において被削材を良好に加工できる。
【０００７】
請求項２に係る数値制御装置は、請求項１に記載の発明の構成に加え、前記確定手段は、前記算出手段が算出した前記所定時間毎の前記最大パワースペクトル値の夫々に対応するピーク周波数を演算する演算手段を備え、前記回転数確定手段は、前記演算手段が演算した前記ピーク周波数の時間変化の中で、閾値よりも大きい変化がある時間に対応する前記回転数を、前記安定回転数として確定することを特徴とする。数値制御装置はピーク周波数の時間変化の中で閾値よりも大きい変化を確定し易いので、安定回転数を容易に確定できる。
【０００８】
請求項３に係る数値制御装置は、請求項１又は２に記載の発明の構成に加え、前記工具の一刃当たりの送り量は、前記工具の送り速度を、前記工具の刃数と前記回転数で除して得られる値であって、前記テスト加工では、前記工具の一刃当たりの送り量が一定となるように、前記工具の送り速度を変化させることを特徴とする。テスト加工では、工具の回転数が変化しても一刃当たりの送り量が一定であることから工具に働く切削力の大きさは一定である。それ故、数値制御装置はテスト加工中の最大パワースペクトル値を、工具の回転数に依らずに所定時間毎に比較できる。
【０００９】
請求項４に係る数値制御装置は、請求項１から３の何れか一つに記載の発明の構成に加え、前記解析手段は、前記所定時間毎の前記外乱トルクを、周波数解析手段を用いて、前記パワースペクトル値を解析することを特徴とする。それ故、数値制御装置は外乱トルクのパワースペクトル値を迅速且つ良好に解析できる。
【００１０】
請求項５に係る数値制御装置は、請求項１から４の何れか一つに記載の発明の構成に加え、前記測定手段は、前記機械の駆動軸の位置情報とトルクに基づき前記外乱トルクを推定する外乱オブザーバを備えたことを特徴とする。それ故、数値制御装置は機械に生じる外乱トルクを容易に推定できる。数値制御装置は外部センサを必要としないので、外部センサにかかる費用を削減できる。
【００１１】
請求項６に係る数値制御装置は、請求項１から５の何れか一つに記載の発明の構成に加え、前記テスト加工は、連続的に前記工具の回転数を減少しながら、前記被削材と同一材質のテスト用被削材を加工する条件で行う加工であることを特徴とする。それ故、数値制御装置は、テスト加工中に測定した外乱トルクの周波数解析の結果において、パワースペクトル値のピークを明瞭に確定できるので、安定回転数を明確に確定できる。
【００１２】
請求項７に係る制御方法は、工具を回転して被削材を加工する機械を制御する数値制御装置の制御方法において、前記被削材の実加工に対応するテスト加工を行い、該テスト加工中において所定時間毎に前記工具の回転数と前記機械の外乱トルクを夫々測定する測定工程と、前記測定工程で測定した前記工具の回転数と前記外乱トルクの情報である測定情報に基づき、前記工具にびびり振動が生じにくい前記工具の回転数である安定回転数を確定する確定工程と、前記確定工程で確定した前記安定回転数の情報である安定回転数情報を出力する出力工程とを備え、前記テスト加工は、連続的に前記工具の回転数を増加又は減少しながら、前記被削材と同一材質のテスト用被削材を加工する条件で行う加工であって、前記確定工程は、前記測定情報のうち、前記所定時間毎の前記外乱トルクの周波数解析を行ってパワースペクトル値を解析する解析工程と、前記解析工程で解析した前記パワースペクトル値に基づき、前記所定時間毎の最大パワースペクトル値を算出する算出工程と、前記算出工程で算出した前記最大パワースペクトル値の時間変化の中で、傾きが負から正に変化する変化時間に対応する前記工具の回転数を、前記安定回転数として確定する回転数確定工程とを備え、前記出力工程は、前記回転数確定工程で確定した前記安定回転数の情報を、前記安定回転数情報として出力することを特徴とする。それ故、数値制御装置は上記制御方法を行うことで、請求項１に記載の効果を得ることができる。
【図面の簡単な説明】
【００１３】
【図１】工作機械１の斜視図。
【図２】数値制御装置３０と工作機械１の電気的構成を示すブロック図。
【図３】ＲＡＭ３３の各種記憶領域を示す概念図。
【図４】テスト加工の切削条件の一例を示す図。
【図５】テスト加工の切削方法を側方から見た図。
【図６】テスト加工の切削方法を上方から見た図。
【図７】安定回転数検出処理の流れ図。
【図８】測定処理の流れ図。
【図９】確定処理の流れ図。
【図１０】外乱オブザーバを適用した外乱トルクの推定手順を示す図。
【図１１】ＲＡＭに記憶した測定情報３３３１の一例を示す図。
【図１２】テスト加工における工具回転数、パワースペクトル値、ピーク周波数、ピーク周波数のパワースペクトル値の夫々の変化を示す図表。
【図１３】表示部２５に表示した安定回転数検出結果表の図。
【図１４】３０００ｒｐｍ／ｓの変化量で工具回転数を増減した時のＳＤＦＴ解析結果を示す図表。加工した
【図１５】５０００ｒｐｍ／ｓの変化量で工具回転数を増減した時のＳＤＦＴ解析結果を示す図表。
【図１６】７５００ｒｐｍ／ｓの変化量で工具回転数を増減した時のＳＤＦＴ解析結果を示す図表。
【発明を実施するための形態】
【００１４】
本発明の実施形態を図面を参照して説明する。以下説明は、図中に矢印で示す上下、左右、前後を使用する。工作機械１の左右方向、前後方向、上下方向は、夫々、Ｘ軸方向、Ｙ軸方向、Ｚ軸方向である。図１に示す工作機械１は主軸９に装着した工具４を高速回転し、作業台１０上に保持した被削材３に切削加工を施す機械である。数値制御装置３０（図２参照）は工作機械１の動作を制御する。
【００１５】
図１を参照し、工作機械１の構造を説明する。工作機械１は、基台２、コラム５、主軸ヘッド７、主軸９、作業台１０、工具交換装置２０、制御箱６、操作パネル（図示略）等を備える。基台２は金属製の略直方体状の土台である。コラム５は基台２上部後方に立設する。主軸ヘッド７はコラム５前面に沿ってＺ軸方向に移動可能に設ける。主軸ヘッド７は内部に主軸９を回転可能に支持する。主軸９は主軸ヘッド７下部に装着穴（図示略）を有する。主軸９は該装着穴に工具４を装着し、主軸モータ５２（図２参照）の駆動で回転する。主軸モータ５２は主軸ヘッド７に設ける。主軸ヘッド７はコラム５前面に設けたＺ軸移動機構（図示略）でＺ軸方向に移動する。数値制御装置３０はＺ軸モータ５１の駆動を制御することで、主軸ヘッド７をＺ軸方向に移動制御する。
【００１６】
作業台１０は基台２上部中央に設け、Ｘ軸モータ５３（図２参照）、Ｙ軸モータ５４（図２参照）、ガイド機構（図示略）等により、Ｘ軸方向とＹ軸方向に移動可能である。数値制御装置３０はＸ軸モータ５３とＹ軸モータ５４の夫々の駆動を制御することで、作業台１０をＸ軸方向とＹ軸方向に移動制御する。
【００１７】
工具交換装置２０は主軸ヘッド７の前側に設け、円盤型の工具マガジン２１を備える。工具マガジン２１は外周に複数の工具（図示略）を放射状に保持し、工具交換指令が指示する工具を工具交換位置に割り出す。工具交換指令はＮＣプログラムで指令する。工具交換位置は工具マガジン２１の最下部位置である。工具交換装置２０は主軸９に装着した工具４と工具交換位置にある工具とを入れ替え交換する。
【００１８】
制御箱６は数値制御装置３０(図２参照）を格納する。数値制御装置３０は、Ｚ軸モータ５１、主軸モータ５２、Ｘ軸モータ５３、Ｙ軸モータ５４を夫々制御し、作業台１０上に保持した被削材３と主軸９に装着した工具４を相対移動することで各種加工を被削材３に施す。各種加工とは、例えばドリル、タップ等を用いた穴空け加工、エンドミル、フライス等を用いた側面加工等である。操作パネル（図示略）は、例えば工作機械１を覆うカバー（図示略）の外壁に設ける。操作パネルは入力部２４と表示部２５(図２参照）を備える。入力部２４は各種情報、操作指示等の入力を受け付ける。表示部２５は各種画面を表示する。
【００１９】
図２を参照し、数値制御装置３０と工作機械１の電気的構成を説明する。数値制御装置３０は、ＣＰＵ３１、ＲＯＭ３２、ＲＡＭ３３、記憶装置３４、入出力部３５、駆動回路５１Ａ〜５５Ａ等を備える。ＣＰＵ３１は数値制御装置３０を統括制御する。ＲＯＭ３２は、メインプログラム、安定回転数検出プログラムを含む各種プログラム等を記憶する。メインプログラムは、メイン処理を実行するものである。メイン処理は、ＮＣプログラムを１ブロックずつ読み込んで各種動作を実行する。ＮＣプログラムは各種制御指令を含む複数のブロックで構成し、工作機械１の軸移動、工具交換等を含む各種動作をブロック単位で制御するものである。安定回転数検出プログラムは、後述する安定回転数検出処理（図７参照）を実行するものである。
【００２０】
ＲＡＭ３３は後述する各種記憶領域（図３参照）を備え、各種情報を記憶する。記憶装置３４は不揮発性であり、ＮＣプログラムと各種データを記憶する。ＣＰＵ３１は作業者が入力部２４で入力したＮＣプログラムに加え、外部入力で読み込んだＮＣプログラム等を記憶装置３４に記憶できる。
【００２１】
駆動回路５１ＡはＺ軸モータ５１とエンコーダ５１Ｂに接続する。駆動回路５２Ａは主軸モータ５２とエンコーダ５２Ｂに接続する。駆動回路５３ＡはＸ軸モータ５３とエンコーダ５３Ｂに接続する。駆動回路５４ＡはＹ軸モータ５４とエンコーダ５４Ｂに接続する。駆動回路５５Ａは工具マガジン２１を駆動するマガジンモータ５５とエンコーダ５５Ｂに接続する。駆動回路５１Ａ〜５５ＡはＣＰＵ３１から指令を受け、対応する各モータ５１〜５５に駆動電流を夫々出力する。駆動回路５１Ａ〜５５Ａはエンコーダ５１Ｂ〜５５Ｂからフィードバック信号を受け、位置と速度のフィードバック制御を行う。入出力部３５は入力部２４と表示部２５に夫々接続する。
【００２２】
図３を参照し、ＲＡＭ３３の記憶領域を説明する。ＲＡＭ３３はバッファ領域３３１、パラメータ記憶領域３３２、測定情報記憶領域３３３等を備える。バッファ領域３３１は、後述する安定回転数検出処理を実行する時の作業領域である。パラメータ記憶領域３３２は後述する各種パラメータを記憶する。測定情報記憶領域３３３は、後述する測定情報３３３１（図１１参照）を記憶する。
【００２３】
テスト加工を説明する。数値制御装置３０は、後述する安定回転数検出処理の中で、実加工を想定した一回のテスト加工を実施し、該テスト加工中に生じるびびり振動を周波数解析することで、びびり振動が生じにくい工具４の回転数（以下、安定回転数と呼ぶ）を確定する。実加工とは、工作機械１が被削材３に実施する加工を意味する。テスト加工を実施する為、作業者は、安定回転数検出処理を実行する前の準備として、テスト加工の切削条件を設定する必要がある。
【００２４】
図４〜図６を参照し、テスト加工の切削条件を説明する。図４に示すテスト加工の切削条件は、φ１０（ｍｍ）のスクエアエンドミルを用いて被削材３の側面を切削する実加工を想定したものである。図５，図６は、図４に示す切削条件で実施するテスト加工の切削方法を示す。図４に示す如く、作業者はテスト加工の切削条件について、例えば、一刃当たりの送り量（ｍｍ）、半径方向切込量（ｍｍ）、切削幅（ｍｍ）、切削油の有無、工具種類、突出量（ｍｍ）、刃数、被削材を夫々設定する。突出量は、主軸９の先端から工具４先端までの長さである（図５参照）。切削幅はテスト加工用の被削材３００を工具４で切削する部分の幅であり、工具４の軸方向に平行な長さである（図５参照）。刃数は工具４の刃数である。
【００２５】
例えば、一刃あたりの送り量は０．０８７５（ｍｍ）、半径方向切込量は０．５（ｍｍ）、切削幅は１５（ｍｍ）、切削油は無、工具はφ１０（ｍｍ）のスクエアエンドミル、突出量は６０（ｍｍ）、刃数は２、被削材３００はアルミニウム（Ａ２０１７）である。テスト加工の切削条件は、実加工時の切削条件に近いものにするのが好ましい。図５に示す如く、作業者は、テスト加工の切削条件(図４参照)に従い、被削材３００と工具４を用意する。作業者は、被削材３００を作業台１０上に冶具（図示略）を用いて固定し、工具４を主軸９に装着する。
【００２６】
図７〜図９を参照し、ＣＰＵ３１が実行する安定回転数検出処理を説明する。テスト加工の準備が完了した時、作業者は操作パネルの入力部２４で安定回転数検出処理の実行開始の指示を入力する。ＣＰＵ１１は入力部２４で該指示の入力を受け付けた時、ＲＯＭ３２から安定回転数検出プログラムを読み込んで本処理を実行する。図７に示す如く、先ずＣＰＵ１１は測定処理を実行する（Ｓ１）。
【００２７】
図８を参照し、測定処理を説明する。ＣＰＵ３１はパラメータ入力処理を実行する（Ｓ１１）。パラメータ入力処理では、ＣＰＵ３１は操作パネルの入力部２４で作業者が入力する五つのパラメータを受け付け、ＲＡＭ３３のパラメータ記憶領域３３２（図３参照）に記憶する。五つのパラメータは以下の通りである。
・開始回転数：テスト加工開始時の工具４の回転数（ｒｐｍ）
・終了回転数：テスト加工終了時の工具４の回転数（ｒｐｍ）
・一刃当たりの送り量：工具４の一刃あたりの工具４の送り量（ｍｍ）
・回転数変更速度：工具４の回転数の減少速度（ｒｐｍ／ｓ）
・刃数：工具４の刃数
本実施形態では、開始回転数は９０００（ｒｐｍ）、終了回転数は５４００（ｒｐｍ）、回転数変更速度は１０８０（ｒｐｍ／ｓ）、刃数は２とする。パラメータ入力処理は入力部２４で作業者によるパラメータ入力の確定指示を受け付けた時に完了する。
【００２８】
パラメータ入力処理の完了後、ＣＰＵ３１はテスト加工を開始する為、主軸９の回転数を開始回転数に設定する（Ｓ１２）。主軸９と共に工具４は９０００（ｒｐｍ）で回転する。ＣＰＵ３１は送り速度を以下の数１で演算する（Ｓ１３）。
（数１）
送り速度（ｍｍ／ｍｉｎ）＝現在の主軸の回転数（ｒｐｍ）×一刃あたりの送り量（ｍｍ）×工具の刃数
一刃当たりの送り量と工具の刃数は、ＲＡＭ３３のパラメータ記憶領域３３２に記憶した値を用いる。後述の如く、工具４の回転数は時間と共に減少するが、一刃当たりの送り量と工具の刃数は予め設定した値であるので、送り速度が変化する。即ち、一刃当たりの送り量は一定であるので、工具４における工具に働く切削力の大きさは一定である。
【００２９】
ＣＰＵ３１はＳ１３の処理で演算した送り速度に従い、被削材３００を移動する（Ｓ１４）。本実施形態の工作機械１は、主軸９がＺ軸方向に移動し、作業台１０がＸ軸方向とＹ軸方向に移動する機械である。それ故、図５，図６に示す如く、ＣＰＵ３１は、回転する工具４に対し、作業台１０をＸ軸方向に演算した送り速度で移動することで、被削材３００の側面を工具４で切削する。例えば、主軸がＸ軸方向、Ｙ軸方向、Ｚ軸方向に移動可能である工作機械の場合、ＣＰＵ３１は被削材３００に対して、回転する工具４をＸ軸方向に移動すればよい。本実施形態のテスト加工では、被削材３００をＸ軸方向に移動して行うが、Ｙ軸方向でもよい。
【００３０】
ＣＰＵ３１は所定周期毎に外乱トルクを測定し、現在の時間、現在の回転数と共に、測定情報３３３１（図１１参照）としてＲＡＭ３３の測定情報記憶領域３３３（図３参照）に記憶する（Ｓ１５）。なお、本実施形態の外乱トルクを測定する所定周期は０．１２５ｍｓ（周波数：８０００Ｈｚ）であるが、これに限定しない。外乱トルクの測定について、ＣＰＵ３１は主軸モータ５２の電流モニタ値と主軸角速度に「外乱オブザーバ」を適用する。
【００３１】
外乱オブザーバを用いた外乱トルクの測定原理を説明する。工作機械１の主軸９における運動方程式は，モータトルクＴｍと負荷トルクＴｌ（切削トルク、摩擦トルクの合計）を考慮して以下の数２のように表すことができる。
【数２】
0006311635.tif 000002
ω［ｒａｄ／ｓ］は主軸角速度である。Ｉａ［Ａ］は電流モニタ値である。Ｊ［ｋｇ／ｍ２］は主軸慣性モーメントである。Ｋｔ［Ｎｍ／Ａ］は主軸モータ５２のトルク定数である。
【００３２】
慣性モーメントとトルク定数は、機械的な構成やトルクリップルなどにより、夫々ΔＪ、ΔＫｔのばらつきがある。ばらつきは、通常、切削負荷と比べて微少であり無視できる。故に外乱トルクＴｄｉｓは切削トルクＴｃｕｔと摩擦トルクＴｆｒｉｃの合計と定義される。それ故、以下の数３のように、電流モニタ値と主軸角速度とから外乱トルクＴｄｉｓを推定できる。
【数３】
0006311635.tif 000003
【００３３】
ＣＰＵ３１は、Ｓ１５の処理の中で、電流モニタ値をトルク情報として主軸モータ５２から取得する。ＣＰＵ３１は取得した電流モニタ値とトルク定数によりトルク（Ｔｍ）を算出する。更にＣＰＵ３１は主軸モータ５２に設けたエンコーダ５２Ｂ（図２参照）より駆動軸の位置情報として主軸角速度（ω）を取得する。
【００３４】
主軸角速度を微分すると高周波域におけるノイズは大きくなる。図１０に示す如く、ＣＰＵ３１は図８に示す測定処理のＳ１５においてローパスフィルタ（ＬＰＦ）を用いて高周波ノイズを遮断して外乱トルクを推定する。推定可能な外乱トルクの周波数はＬＰＦに依存する。ＣＰＵ３１は外乱トルクを推定する為に、ＬＰＦの遮断周波数をびびり振動周波数よりも高く設定する。それ故、ＣＰＵ３１は切削加工中の外乱トルクを推定できる。
【００３５】
図８の測定処理に戻り、ＣＰＵ３１は、ＲＡＭ３３のパラメータ記憶領域３３２（図３参照）から回転数変更速度を読み込み、該回転数変更速度に従い回転数を減少する（Ｓ１６）。ＣＰＵ３１は現在の回転数が終了回転数に到達したか判断する（Ｓ１７）。終了回転数に到達していない時（Ｓ１７：ＮＯ）、ＣＰＵ３１はＳ１３に戻り、減少した回転数で送り速度を演算する。ＣＰＵ３１は演算した送り速度に変更し、回転する工具４に対して作業台１０を引き続き移動する（Ｓ１４）。ＣＰＵ３１は所定周期毎に外乱トルクを測定し続け、現在の時間と回転数と共に、ＲＡＭ３３の測定情報記憶領域３３３に記憶する（Ｓ１５）。
【００３６】
テスト加工中、ＣＰＵ３１は回転数変更速度に従い回転数を減少し（Ｓ１６）、現在の回転数が終了回転数に到達するまで（Ｓ１７：ＮＯ）、ＣＰＵ３１はＳ１３〜Ｓ１６の処理を繰り返す。現在の回転数が終了回転数に到達した時（Ｓ１７：ＹＥＳ）、ＣＰＵ３１は工具４から被削材３００を退避し、主軸９の回転を停止してテスト加工を終了する（Ｓ１８）。ＣＰＵ３１は測定処理を終了し、図３のＳ２に処理を進め、確定処理を実行する。
【００３７】
図９を参照し、確定処理を説明する。ＲＡＭ３３に記憶した測定情報３３３１において、ＣＰＵ３１は先頭行から一定時間分（例えば０．１秒分）の測定情報を纏めて取得し、ＲＡＭ３３のバッファ領域３３１（図３参照）に読み込む（Ｓ２１）。ＣＰＵ３１はバッファ領域３３１に読み込んだ一定時間分の測定情報について、外乱トルクのＳＤＦＴ（スライディング離散フーリエ変換）による周波数解析を実行する（Ｓ２２）。ＳＤＦＴとは、離散フーリエ変換による信号の周波数解析を行う際に、解析区間をサンプル時間ずつずらすことで時間方向の分解能を高めるための手法である。ＳＤＦＴ解析では、例えば図１２（ｂ）に示す如く、テスト加工中に発生したびびり振動の周波数（Ｈｚ）とパワースペクトル密度（（Ｎｍ）２／Ｈｚ）を解析する。
【００３８】
ＣＰＵ３１は、ＳＤＦＴ解析結果のうち最大となる最大パワースペクトル値を算出する（Ｓ２３）。ＣＰＵ３１は算出した最大パワースペクトル値に対応する周波数であるピーク周波数を確定する（Ｓ２４）。ＣＰＵ３１は最大パワースペクトル値とピーク周波数をＲＡＭ３３に記憶する。ＣＰＵ３１は、今回のピーク周波数と１ステップ前のピーク周波数との差分の絶対値が閾値よりも大きいか判断する（Ｓ２５）。１ステップ前のピーク周波数とは、前回確定したピーク周波数であり、ＲＡＭ３３に記憶した一つ前のピーク周波数を意味する。閾値は記憶装置３４に予め記憶し、作業者が自由に設定可能である。差分の絶対値が閾値以下の時（Ｓ２５：ＮＯ）、びびり振動の周波数は連続しているとみなし、ＣＰＵ３１はＲＡＭ３３に記憶する測定情報３３３１を全て処理したか判断する（Ｓ２８）。
【００３９】
未処理の測定情報が有る時（Ｓ２８：ＮＯ）、ＣＰＵ３１はＲＡＭ３３のバッファ領域３３１に読み込んだ測定情報のうち最も時間の古い０．１２５（ｍｓ）分の測定情報を削除し（Ｓ２９）、測定情報記憶領域３３３に記憶した測定情報３３３１から次の０．１２５（ｍｓ）分の測定情報をバッファ領域３３１に読み込む（Ｓ３０）。ＣＰＵ３１はＳ２２に戻り、新しい時間分を加えた測定情報について、外乱トルクのＳＤＦＴ解析を実行し、上記処理を繰り返す（Ｓ２２〜Ｓ２４）。
【００４０】
Ｓ２５において差分の絶対値が閾値より大きい時（Ｓ２５：ＹＥＳ）、ピーク周波数は不連続であるとみなすことができる。ピーク周波数の不連続な部分では、びびり振動が小さい可能性がある。それ故、ＣＰＵ３１はバッファ領域３３１に読み込んだ測定情報において、回転数の平均値を算出する（Ｓ２６）。ＣＰＵ３１は算出した平均値を安定回転数として、記憶装置３４に記憶する（Ｓ２７）。
【００４１】
図１２を参照し、ピーク周波数の不連続な部分における最大パワースペクトル値の時間変化を検証する。図１２（ａ）はテスト加工中の工具４の回転数の変化を示す。加工開始後、工具４の回転数は９０００（ｒｐｍ）から５４００（ｒｐｍ）まで、１０８０（ｒｐｍ／ｓ）の回転数変更速度に従い減少している。
【００４２】
図１２（ｂ）はテスト加工中に測定した外乱トルクのＳＤＦＴ解析結果を示す。テスト加工開始後、２０００、２２００、２３００（Ｈｚ）付近に、主に三つのパワースペクトルが出現する。びびり振動は夫々の周波数で生じている。回転数が減少を開始するに伴い、三つのパワースペクトルの周波数は徐々に低下するが、テスト加工開始後から２．２秒後、３．１秒後、３．９秒後、４．４秒後に、最初の周波数付近まで瞬間的に戻り、再度徐々に低下することを繰り返している。
【００４３】
図１２（ｃ）は、図１２（ｂ）のＳＤＦＴ解析結果において、所定時間（例えば０．１２５ｍｓ）毎に最大パワースペクトル値を算出し、算出した最大パワースペクトル値の周波数であるピーク周波数を所定時間毎にプロットしたものである。ピーク周波数は、図１２（ｂ）のパワースペクトルの変化に対応するように低下するが、テスト加工開始から２．２秒後、３．１秒後、３．９秒後、４．４秒後において不連続な部分（図中点線で囲った部分を参照）を生じている。不連続な部分とは、ピーク周波数の傾きが負から正に瞬間的に切り換わる部分であって閾値より大きく変化する部分である。
【００４４】
図１２（ｄ）は、図１２（ｃ）に示すピーク周波数の時間変化を、パワースペクトル値の時間変化に置き換えたものである。パワースペクトル値は、ピーク周波数が低下する部分に対応して山なりに上昇するが、ピーク周波数の不連続な部分においては低下している（図中点線で囲った部分参照）。パワースペクトル値が低い部分では、びびり振動は小さい。それ故、本実施形態は、ピーク周波数の時間変化の中で不連続な部分が生じる変化時間に対応する工具４の回転数を、びびり振動が生じにくい安定回転数として確定できる。
【００４５】
図９に戻り、ＣＰＵ３１はＲＡＭ３３に記憶する測定情報を全て処理した時（Ｓ２８：ＹＥＳ）、記憶装置３４は、安定回転数情報として、例えば８２００（ｒｐｍ）、７０５０（ｒｐｍ）、５９００（ｒｐｍ）を記憶する。ＣＰＵ３１は確定処理を終了し、図７のＳ３に処理を進める。
【００４６】
ＣＰＵ３１は表示処理を実行する（Ｓ３）。表示処理では、ＣＰＵ３１は確定処理で記憶装置３４に記憶した安定回転数情報を表示部２５に出力する。例えば、図１３に示す如く、ＣＰＵ３１は、表示部２５に、安定回転数検出結果表を表示する。安定回転数検出結果表は、安定回転数の候補として、例えば、８２００（ｒｐｍ）、７０５０（ｒｐｍ）、５９００（ｒｐｍ）を夫々表示する。作業者は表示部２５に表示した安定回転数検出結果表と、他の条件を勘案した上で、びびり振動が生じにくい最良の回転数を選択できる。また、数値制御装置３０は、実加工でびびり振動が発生した時に、記憶装置３４に記憶した安定回転数情報に基づき、自動的にびびり振動が生じにくい工具回転数に変更するようにしてもよい。
【００４７】
図１４〜図１６を参照し、テスト加工において工具４の回転数を増加させた場合と減少させた場合で、安定回転数の確定のし易さを比較した比較実験を説明する。上記実施形態は、図８に示す測定処理のＳ１５，Ｓ１６において、工具４の回転数を回転数変更速度に従い減少しながら外乱トルクを測定する。本実施形態は、工具４の回転数を増加しながら外乱トルクを測定してもよいが、安定回転数をより明確に確定する為には、回転数を減少しながら外乱トルクを測定するのが好ましい。そこで、上記のテスト加工において、工具４の回転数を増加させた場合と減少させた場合とで、ＳＤＦＴ解析結果におけるパワースペクトル値のピークの明瞭性について比較実験を行った。加工条件は以下の通りである。
・軸方向切込量（切削幅）：１５ｍｍ
・半径方向切込量：０．５ｍｍ
上記の加工条件で、工具４の回転数に対し、３０００ｒｐｍ／ｓ、５０００ｒｐｍ／ｓ、７５００ｒｐｍ／ｓの変化量（加速度に相当）を夫々与え、ＳＤＦＴ解析を行った。比  
【００４８】
実験結果を説明する。図１４は３０００ｒｐｍ／ｓの変化量を与えた場合のＳＤＦＴ解析結果、図１５は７５００ｒｐｍ／ｓの変化量を与えた場合のＳＤＦＴ解析結果、図１６は７５００ｒｐｍ／ｓの変化量を与えた場合のＳＤＦＴ解析結果である。
【００４９】
図１４に示す如く、３０００ｒｐｍ／ｓの変化量で工具４の回転数を増減した場合、５４００ｒｐｍから９０００ｒｐｍに到達するまでの時間、及び９０００ｒｐｍから５４００ｒｐｍに到達するまでの時間は１．２秒である。ＳＤＦＴ解析結果を見ると、工具回転数が増加する時よりも減少する時の方が、パワースペクトル値のピークを明瞭に確認できた。
【００５０】
図１５に示す如く、５０００ｒｐｍ／ｓの変化量で工具４の回転数を増減した場合、５４００ｒｐｍから９０００ｒｐｍに到達するまでの時間、及び９０００ｒｐｍから５４００ｒｐｍに到達するまでの時間は０．７２秒である。ＳＤＦＴ解析結果を見ると、図１４と同様に、工具回転数が増加する時よりも減少する時の方が、パワースペクトル値のピークを明瞭に確認できた。
【００５１】
図１６に示す如く、７５００ｒｐｍ／ｓの変化量で工具４の回転数を増減した場合、５４００ｒｐｍから９０００ｒｐｍに到達するまでの時間、及び９０００ｒｐｍから５４００ｒｐｍに到達するまでの時間は０．４８秒である。ＳＤＦＴ解析結果を見ると、図１４，図１５と同様に、工具回転数が増加する時よりも減少する時の方が、パワースペクトル値のピークは強く出ているが、一部のピークが消失して確認できなかった（図中点線で囲った部分参照）。原因は回転数変更速度が速過ぎてびびり振動が十分発達するだけの時間がなかったことが考えられる。
【００５２】
上記結果により、テスト加工では回転数を減少しながら外乱トルクを測定する方が、パワースペクトル値のピークを明瞭に確定できるので、安定回転数をより明確に確定できることが実証された。また、回転数を減少する際の回転数変更速度は、少なくとも７５００ｒｐｍ／ｓよりも低い速度に設定した方がよいことが実証された。
【００５３】
上記説明にて、図１に示す工作機械１が本発明の機械に相当し、被削材３が本発明の被削材に相当し、被削材３００が本発明のテスト用被削材に相当する。図７のＳ１の処理を実行するＣＰＵ３１が本発明の測定手段に相当し、Ｓ２の処理を実行するＣＰＵ３１が本発明の確定手段に相当し、Ｓ３の処理を実行するＣＰＵ３１が本発明の出力手段に相当し、図９のＳ２２の処理を実行するＣＰＵ３１が本発明の解析手段に相当し、Ｓ２３の処理を実行するＣＰＵ３１が本発明の算出手段に相当し、Ｓ２５〜Ｓ２７の処理を実行するＣＰＵ３１が本発明の回転数確定手段に相当し、Ｓ２４の処理を実行するＣＰＵ３１が本発明の演算手段に相当する。図７のＳ１の工程が本発明の測定工程に相当し、Ｓ２の工程が本発明の確定工程に相当し、Ｓ３の工程が本発明の出力工程に相当し、図９のＳ２２の工程が本発明の解析工程に相当し、Ｓ２３の工程が本発明の工程に相当し、Ｓ２５〜Ｓ２７の工程が本発明の回転数確定工程に相当する。
【００５４】
以上説明の如く、本実施形態の数値制御装置３０は、工具４を回転して被削材３を加工する工作機械１を制御する。数値制御装置３０のＣＰＵ３１は安定回転数検出処理を実行することで、びびり振動が起きにくい工具４の回転数を検出できる。安定回転数検出処理は、測定処理、確定処理、表示処理を備える。測定処理では、ＣＰＵ３１は、被削材３の実加工に対応するテスト加工を行い、該テスト加工中において所定時間毎に工具４の回転数と工作機械１に生じる外乱トルクを夫々測定する。テスト加工は、連続的に工具４の回転数を減少しながら、被削材３と同一材質のテスト用の被削材３００を加工する条件で行う。確定処理では、ＣＰＵ３１は、測定処理で測定した工具４の回転数と外乱トルクの情報である測定情報３３３１に基づき、安定回転数を確定する。出力処理では、ＣＰＵ３１は、確定処理で確定した安定回転数の情報である安定回転数情報を表示部２５に出力する。
【００５５】
具体的には、確定処理にて、ＣＰＵ３１は測定情報３３３１のうち、所定時間毎の外乱トルクの周波数解析を行ってパワースペクトル値を解析する。ＣＰＵ３１は解析したパワースペクトル値に基づき、所定時間毎の最大パワースペクトル値を算出する。ＣＰＵ３１は最大パワースペクトル値の夫々に対応するピーク周波数を演算し、演算したピーク周波数の時間変化の中で、閾値よりも大きい変化がある変化時間に対応する回転数を安定回転数として確定する。
【００５６】
それ故、数値制御装置３０は実加工に対応するテスト加工を一回行うことで、実加工でびびり振動が生じにくい安定回転数を確定できる。作業者は表示部２５に表示した安定回転数情報を参照し、他条件を加味した上で最適な工具回転数を決定できる。それ故、数値制御装置３０は実加工において被削材３を良好に加工できる。また、数値制御装置３０はピーク周波数の時間変化の中で閾値よりも大きい変化を確定し易いので、安定回転数を容易に確定できる。
【００５７】
また本実施形態では特に、工具４の一刃当たりの送り量は、工具４の送り速度を、工具４の刃数と回転数で除して得られる値である。テスト加工では工具４の一刃当たりの送り量が一定となるように、数値制御装置３０は工具４の送り速度を変化させる。テスト加工では、工具４の回転数が変化しても一刃当たりの送り量が一定であることから工具に働く切削力の大きさは一定である。それ故、数値制御装置３０はテスト加工中の最大パワースペクトル値を、工具４の回転数に依らずに所定時間毎に比較できる。
【００５８】
また本実施形態では特に、数値制御装置３０は、所定時間毎の外乱トルクを、高速フーリエ変換処理を用いて、パワースペクトル値を解析する。それ故、数値制御装置３０は外乱トルクのパワースペクトル値を迅速且つ良好に解析できる。
【００５９】
また本実施形態では特に、数値制御装置３０は、工作機械１の主軸モータ５２の位置情報とトルクに基づき外乱トルクを推定する外乱オブザーバを備える。それ故、数値制御装置３０は工作機械１に生じる外乱トルクを容易に推定できる。数値制御装置３０は外部センサを必要としないので、外部センサにかかる費用を削減できる。
【００６０】
また本実施形態では特に、テスト加工は連続的に工具４の回転数を減少しながら、被削材３と同一材質のテスト用の被削材３００を加工する条件で行う。それ故、数値制御装置３０は、テスト加工中に測定した外乱トルクのＳＤＦＴ解析結果において、パワースペクトル値のピークを明瞭に確定できるので、安定回転数を明確に確定できる。
【００６１】
なお、本発明は上記実施形態に限らず、各種の変形が可能なことはいうまでもない。上記実施形態は、工作機械１の動作を制御する数値制御装置であるが、びびり振動の検出専用の振動検出装置であってもよく、所謂パーソナルコンピュータ等の汎用型の装置であってもよい。
【００６２】
また上記実施形態では、主軸モータ５２の位置情報とトルク情報を用いて説明したが、例えば、Ｘ軸モータ５３、Ｙ軸モータ５４、Ｚ軸モータ５１の情報を用いてもよいし、主軸モータ５２、Ｘ軸モータ５３、Ｙ軸モータ５４、Ｚ軸モータ５１の内から複数の情報を元に総合的に判断してもよい。Ｘ軸、Ｙ軸、Ｚ軸は、被削材３に対して工具８を直交する３軸（Ｘ軸方向、Ｙ軸方向、Ｚ軸方向）に移動させる制御軸である。Ｘ軸モータ５３、Ｙ軸モータ５４、Ｚ軸モータ５１は、これら制御軸の駆動モータである。
【００６３】
また上記実施形態では、縦型の工作機械１を例示して説明したが、本発明は主軸が略水平方向に延びる横型の工作機械であってもよい。
【００６４】
また上記実施形態では、図９の確定処理のＳ２５において、今回のピーク周波数と１ステップ前のピーク周波数との差分の絶対値が閾値よりも大きいか判断するが、例えば、図１２（ｄ）に示す如く、所定時間毎のピーク周波数のパワースペクトル値を算出し、該算出したパワースペクトル値の時間変化を求めた上で、以下の二つの条件を何れも満たすか判断するようにしてもよい。
・（ピーク周波数のパワースペクトル値）−（１ステップ前のピーク周波数のパワースペクトル値）の符号が正であること。
・（１ステップ前のパワースペクトル値）−（２ステップ前のパワースペクトル値）の符号が負又は０であること。
上記二つの条件を何れも満たす時、図１２（ｄ）に示す如く、ピーク周波数のパワースペクトル値の時間変化の傾きは負から正に変わる時である（図中点線で囲った部分参照）。それ故、ＣＰＵ３１は処理をＳ２６，Ｓ２７に進め、その時間に対応する工具４の回転数を安定回転数として確定できる。
【００６５】
また上記実施形態では、図９の確定処理のＳ２４，Ｓ２５において、ＣＰＵ３１は最大パワースペクトル値の夫々に対応するピーク周波数を演算し、演算したピーク周波数の時間変化の中で、閾値よりも大きい変化がある変化時間に対応する回転数を安定回転数として確定するが、例えば、最大パワースペクトル値の時間変化の中で、傾きが負から正に変化する変化時間に対応する工具４の回転数を安定回転数として確定してもよい。
【００６６】
また上記実施形態では、図７のＳ３において、ＣＰＵ３１は安定回転数情報を表示部２５に表示するが、例えば、入出力部３５を介して外部装置（図示略）に出力するようにしてもよい。
【００６７】
また上記実施形態では、測定情報３３３１の外乱トルクについてＳＤＦＴ解析を行うが、例えば、ＦＦＴ解析、ウェーブレット解析等の方法で周波数解析を行ってもよい。
【符号の説明】
【００６８】
１  工作機械
３  被削材
４  工具
３０  数値制御装置
３１  ＣＰＵ
３００  被削材
【図１】
0006311635.tif 000004
【図２】
0006311635.tif 000005
【図３】
0006311635.tif 000006
【図４】
0006311635.tif 000007
【図５】
0006311635.tif 000008
【図６】
0006311635.tif 000009
【図７】
0006311635.tif 000010
【図８】
0006311635.tif 000011
【図９】
0006311635.tif 000012
【図１０】
0006311635.tif 000013
【図１１】
0006311635.tif 000014
【図１２】
0006311635.tif 000015
【図１３】
0006311635.tif 000016
【図１４】
0006311635.tif 000017
【図１５】
0006311635.tif 000018
【図１６】
0006311635.tif 000019