知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2024-07-08
(45)【発行日】2024-07-17
(54)【発明の名称】数値制御装置、及び制御方法
(51)【国際特許分類】
B23Q 15/18 20060101AFI20240709BHJP
G05B 19/404 20060101ALI20240709BHJP
【ＦＩ】
B23Q15/18
G05B19/404 K
【請求項の数】
8
(21)【出願番号】P 2019219388
(22)【出願日】2019-12-04
(65)【公開番号】P2021088024
(43)【公開日】2021-06-10
【審査請求日】2022-10-21
(73)【特許権者】
【識別番号】390008235
【氏名又は名称】ファナック株式会社
(74)【代理人】
【識別番号】100106002
【弁理士】
【氏名又は名称】正林 真之
(74)【代理人】
【識別番号】100165157
【弁理士】
【氏名又は名称】芝 哲央
(74)【代理人】
【識別番号】100160794
【弁理士】
【氏名又は名称】星野 寛明
(72)【発明者】
【氏名】澤岡 浩貴
【審査官】小川 真
(56)【参考文献】
【文献】特開２０００－３２２１１５（ＪＰ，Ａ）
【文献】特開２０１５－０９３３４６（ＪＰ，Ａ）
【文献】特開２０１６－１１０４４３（ＪＰ，Ａ）
【文献】特開平０７－１２４８４９（ＪＰ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ｂ２３Ｑ １５／１８、１７／２０
Ｇ０５Ｂ １９／４０４
(57)【特許請求の範囲】
【請求項1】
指令解析手段から受信した切削指令が示す指令座標値で工作機械に切削させる数値制御装置であって、切削されたテストワークの形状を前記工作機械に機上測定させ、測定された前記テストワークの形状を示す測定データを取得する測定手段と、
前記切削指令が示す指令形状と前記測定手段により取得された前記測定データとに基づいて、切削において前記工作機械に作用する力及び速度によって生じる動的誤差を補正する動的補正パラメータを算出する動的補正パラメータ算出手段と、
算出された前記動的補正パラメータに基づいて、前記指令座標値に対して前記動的誤差を少なくとも工作機械の直角度誤差として補正する動的補正手段と、を備え、
前記指令形状、及び前記テストワークの形状が円であり、
前記測定手段は、ＸＹ平面で前記テストワークに切削された半径Ｒ ０ の穴の前記測定データを取得し、
前記動的補正パラメータ算出手段は、
前記測定データが示す前記テストワークの形状に楕円をベストフィットして、前記楕円のＸ軸、Ｙ軸それぞれの半径Ｒ ｘ ，Ｒ ｙ を取得し、
ＸＹ平面における動的補正パラメータＷ ｚｘ 、Ｗ ｚｙ を数１式で算出する、数値制御装置。
【数1】
【請求項2】
動的補正パラメータ補間手段をさらに備え、前記測定手段は、互いに異なる複数の切削負荷の各々で切削された前記テストワークの形状の前記測定データを取得し、
前記動的補正パラメータ算出手段は、前記複数の切削負荷の各々における前記指令形状と前記測定データとに基づいて前記動的補正パラメータを算出し、
前記動的補正パラメータ補間手段は、算出された複数の前記動的補正パラメータを用いて、任意の切削負荷における動的補正パラメータを補間し、
前記動的補正手段は、補間された前記任意の切削負荷における動的補正パラメータに基づいて前記切削指令の前記指令座標値を補正する、請求項１に記載の数値制御装置。
【請求項3】
前記テストワークの切削時の負荷電流を取得し、取得した前記負荷電流、前記指令形状と前記測定データとのズレ量、及び前記動的補正パラメータに基づいて少なくともＸ軸及びＹ軸における負荷電流と切削点負荷との関係を示す関係パラメータを算出する関係パラメータ算出手段と、前記関係パラメータと前記負荷電流とに基づいて前記切削点負荷を算出する切削点負荷算出手段と、をさらに備える、請求項１に記載の数値制御装置。
【請求項4】
前記負荷電流は、前記工作機械に含まれるサーボモータの負荷電流から、前記工作機械に含まれる工具を加速させた時の加速度による負荷電流、及び前記工具を軸移動させた時の軸移動による負荷電流を引いて算出される、請求項３に記載の数値制御装置。
【請求項5】
前記加速度による負荷電流、及び前記軸移動による負荷電流は、切削を伴わない前記工具の移動時の負荷電流から算出される、請求項４に記載の数値制御装置。
【請求項6】
前記切削を伴わない前記工具の移動は、前記工具を円状に時計周りに移動させる移動、及び前記工具を円状に反時計周りに移動させる移動の両方を行う、請求項５に記載の数値制御装置。
【請求項7】
コンピュータにより実行される、指令解析手段から受信した切削指令が示す指令座標値で工作機械に切削させる制御方法であって、切削されたテストワークの形状を前記工作機械に機上測定させ、測定された前記テストワークの形状を示す測定データを取得する測定ステップと、
前記切削指令が示す指令形状と取得された前記測定データとに基づいて、切削において前記工作機械に作用する力及び速度によって生じる動的誤差を補正する動的補正パラメータを算出する動的補正パラメータ算出ステップと、
算出された前記動的補正パラメータに基づいて、前記指令座標値に対して前記動的誤差を補正する動的補正ステップと、を備え、
前記指令形状、及び前記テストワークの形状が円であり、
前記測定ステップは、ＸＹ平面で前記テストワークに切削された半径Ｒ０の穴の前記測定データを取得し、
前記動的補正パラメータ算出ステップは、
前記測定データが示す前記テストワークの形状を円でベストフィットして、前記円の半径Ｒｔを取得し、
切削負荷で生じる前記工作機械に含まれる工具のたわみ量の動的補正パラメータＷｔ、及び工具径補正量ｒｔ’を数２式で算出する、制御方法。
【数2】
【請求項8】
コンピュータにより実行される、指令解析手段から受信した切削指令が示す指令座標値で工作機械に切削させる制御方法であって、切削されたテストワークの形状を前記工作機械に機上測定させ、測定された前記テストワークの形状を示す測定データを取得する測定ステップと、
前記切削指令が示す指令形状と取得された前記測定データとに基づいて、切削において前記工作機械に作用する力及び速度によって生じる動的誤差を補正する動的補正パラメータを算出する動的補正パラメータ算出ステップと、
算出された前記動的補正パラメータに基づいて、前記指令座標値に対して前記動的誤差を補正する動的補正ステップと、を備え、
前記指令形状、及び前記テストワークの形状が円であり、
前記測定ステップは、ＸＹ平面で前記テストワークに切削された半径Ｒ０の穴の前記測定データを取得し、
前記動的補正パラメータ算出ステップは、
前記測定データが示す前記テストワークの形状に楕円をベストフィットして、前記楕円のＸ軸、Ｙ軸それぞれの半径Ｒｘ，Ｒｙを取得し、
ＸＹ平面における動的補正パラメータＷｚｘ、Ｗｚｙを数３式で算出する、制御方法。
【数3】
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、数値制御装置、及び制御方法に関する。
【背景技術】
【0002】
工作機械の剛性及び熱変形、又は工具のたわみ等により、切削されたワークの形状と、設計時の形状との間に誤差が生じてしまう。そこで、レーザー干渉計、オートコリメータ、水準器等を用いて加工されたワークの形状を予め測定し、測定した誤差に基づいて補正する方法がある。
しかしながら、上述の補正は、静的誤差を補正するものであり、切削時に発生する動的誤差を補正することが困難である。なお、動的誤差とは、工作機械に作用する力及び速度によって生じる誤差のことで、例えば切削点負荷により剛性の低い個所に生じる、工作機械の直角度誤差や工具のたわみによる誤差である。
この点、工具のたわみによる圧力をセンサを用いて検出し、検出された圧力に基づいて工具のたわみ量を補正することにより、高速度の切削加工であっても高精度に加工することができる技術が知られている。例えば、特許文献１参照。
しかしながら、上述の補正は、静的誤差を補正するものであり、切削時に発生する動的誤差を補正することが困難である。なお、動的誤差とは、工作機械に作用する力及び速度によって生じる誤差のことで、例えば切削点負荷により剛性の低い個所に生じる、工作機械の直角度誤差や工具のたわみによる誤差である。
この点、工具のたわみによる圧力をセンサを用いて検出し、検出された圧力に基づいて工具のたわみ量を補正することにより、高速度の切削加工であっても高精度に加工することができる技術が知られている。例えば、特許文献１参照。
【先行技術文献】
【特許文献】
【0003】
【文献】特開平５－３１８２８３号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかしながら、工具のたわみ量を補正するために、別途センサを用意して工作機械に取り付ける必要があり、コストがかかる。
また、補正量算出の基となる切削点負荷と加工結果のずれとの関係を算出することが難しい。
また、補正量算出の基となる切削点負荷と加工結果のずれとの関係を算出することが難しい。
【0005】
そこで、センサを用いることなく、動的誤差を精度良く補正することが望まれている。
【課題を解決するための手段】
【0006】
（１）本開示の数値制御装置の一態様は、指令解析手段から受信した切削指令が示す指令座標値で工作機械に切削させる数値制御装置であって、切削されたテストワークの形状を前記工作機械に機上測定させ、測定された前記テストワークの形状を示す測定データを取得する測定手段と、前記切削指令が示す指令形状と前記測定手段により取得された前記測定データとに基づいて、切削において前記工作機械に作用する力及び速度によって生じる動的誤差を補正する動的補正パラメータを算出する動的補正パラメータ算出手段と、算出された前記動的補正パラメータに基づいて、前記指令座標値に対して前記動的誤差を補正する動的補正手段と、を備え、前記動的補正パラメータ算出手段は、前記指令形状と前記測定データとの比較から前記動的誤差のみを取得し、取得した前記動的誤差から前記動的補正パラメータを算出する。
【0007】
（２）本開示の制御方法の一態様は、コンピュータにより実行される、指令解析手段から受信した切削指令が示す指令座標値で工作機械に切削させる制御方法であって、切削されたテストワークの形状を前記工作機械に機上測定させ、測定された前記テストワークの形状を示す測定データを取得する測定ステップと、前記切削指令が示す指令形状と取得された前記測定データとに基づいて、切削において前記工作機械に作用する力及び速度によって生じる動的誤差を補正する動的補正パラメータを算出する動的補正パラメータ算出ステップと、算出された前記動的補正パラメータに基づいて、前記指令座標値に対して前記動的誤差を補正する動的補正ステップと、を備え、前記動的補正パラメータ算出ステップは、前記指令形状と前記測定データとの比較から前記動的誤差のみを取得し、取得した前記動的誤差から前記動的補正パラメータを算出する。
【発明の効果】
【0008】
一態様によれば、センサを用いることなく、動的誤差を精度良く補正することができる。
【図面の簡単な説明】
【0009】
【図1】第１実施形態に係る数値制御装置の機能的構成例を示す機能ブロック図である。
【図2】工作機械の一例を示す図である。
【図3A】直角度誤差の一例を示す図である。
【図3B】工具のたわみによる誤差の一例を示す図である。
【図4】工作機械にテストワークを切削させる一例を示す図である。
【図5】切削点負荷が円周方向にもかかる場合のズレ量の一例を示す図である。
【図6】動的補正パラメータと工作機械との関係の一例を示す図である。
【図7】直角度誤差の動的補正パラメータについて説明する一例を示す図である。
【図8】測定データに円をベストフィットする一例を示す図である。
【図9】切削点負荷の方向の一例を示す図である。
【図10】工具のたわみによる誤差の補正の一例を示す図である。
【図11】第１実施形態に係る数値制御装置の一例を示す図である。
【図12】第１実施形態に係る数値制御装置の一例を示す図である。
【図13】工具のたわみによる誤差の補正の一例を示す図である。
【図14】第２実施形態に係る数値制御装置のうち第１実施形態の数値制御装置に追加される機能的構成例を示す機能ブロック図である。
【図15A】動的補正パラメータ補間手段の補間処理を説明する一例を示す図である。
【図15B】動的補正パラメータ補間手段の補間処理を説明する一例を示す図である。
【図16】複数の動的補正パラメータをＭ次関数で補間する場合の一例を示す図である。
【図17】第３実施形態に係る数値制御装置のうち第１実施形態の数値制御装置に追加される機能的構成例を示す機能ブロック図である。
【図18A】工具を時計周りに円運動させる場合の一例を示す図である。
【図18B】工具を反時計周りに円運動させる場合の一例を示す図である。
【発明を実施するための形態】
【0010】
＜第１実施形態＞
まず、本実施形態の概略を説明する。本実施形態では、数値制御装置は、後述するテストワークを所定の形状に切削加工して、切削されたテストワークの形状を工作機械に機上測定させ、測定されたテストワークの形状を示す測定データを取得する。数値制御装置は、切削指令が示す指令形状と取得された測定データとに基づいて、切削加工に伴い工作機械に作用する力及び速度によって生じる動的誤差を補正する動的補正パラメータを算出する。数値制御装置は、算出された動的補正パラメータに基づいて指令座標値に対して動的誤差を補正する。
まず、本実施形態の概略を説明する。本実施形態では、数値制御装置は、後述するテストワークを所定の形状に切削加工して、切削されたテストワークの形状を工作機械に機上測定させ、測定されたテストワークの形状を示す測定データを取得する。数値制御装置は、切削指令が示す指令形状と取得された測定データとに基づいて、切削加工に伴い工作機械に作用する力及び速度によって生じる動的誤差を補正する動的補正パラメータを算出する。数値制御装置は、算出された動的補正パラメータに基づいて指令座標値に対して動的誤差を補正する。
【0011】
これにより、本実施形態によれば、「工具のたわみによる圧力を検出するセンサを用いることなく、動的誤差を精度良く補正する」という課題を解決することができる。
以上が本実施形態の概略である。
以上が本実施形態の概略である。
【0012】
次に、本実施形態の構成について図面を用いて詳細に説明する。
【0013】
図１は、第１実施形態に係る数値制御装置の機能的構成例を示す機能ブロック図である。なお、制御方法の観点に基づく説明は、「手段」を「ステップ」に置き換えることで説明できるため、省略する。
【0014】
数値制御装置１０は、図示しない接続インタフェースを介して、工作機械２０と互いに直接接続されてもよい。なお、数値制御装置１０、及び工作機械２０は、ＬＡＮ（Ｌｏｃａｌ Ａｒｅａ Ｎｅｔｗｏｒｋ）やインターネット等の図示しないネットワークを介して相互に接続されていてもよい。この場合、数値制御装置１０、及び工作機械２０は、かかる接続によって相互に通信を行うための図示しない通信部を備えている。
【0015】
＜工作機械２０＞
工作機械２０は、Ｘ軸、Ｙ軸、及びＺ軸方向に主軸頭が移動する公知の直交３軸の工作機械であり、数値制御装置１０からの動作指令（切削指令）に基づいて動作する。
図２は、工作機械２０の一例を示す図である。
図２に示すように、工作機械２０は、ＸＹ平面状に配置されたテーブル（定盤）２１と、テーブル２１の両端の位置に鉛直（Ｚ軸）方向に設けられた支柱２２（１）、２２（２）と、支柱２２（１）と支柱２２（２）との間に水平（Ｘ軸）方向に設けられた支柱２３とで構成される。
工作機械２０は、Ｘ軸、Ｙ軸、及びＺ軸方向に主軸頭が移動する公知の直交３軸の工作機械であり、数値制御装置１０からの動作指令（切削指令）に基づいて動作する。
図２は、工作機械２０の一例を示す図である。
図２に示すように、工作機械２０は、ＸＹ平面状に配置されたテーブル（定盤）２１と、テーブル２１の両端の位置に鉛直（Ｚ軸）方向に設けられた支柱２２（１）、２２（２）と、支柱２２（１）と支柱２２（２）との間に水平（Ｘ軸）方向に設けられた支柱２３とで構成される。
【0016】
主軸頭２４及び主軸頭２４に取り付けられた工具２５は、Ｘ軸サーボモータ３１により支柱２３に対してＸ軸方向に移動するとともに、Ｚ軸サーボモータ３３により支柱２３に対してＺ軸方向に上下動する。また、支柱２２（１）、２２（２）、及び支柱２３で構成される門は、Ｙ軸サーボモータ３２によりＹ軸方向に移動する。
【0017】
＜数値制御装置１０＞
数値制御装置１０は、当業者にとって公知の数値制御装置であり、制御情報に基づいて動作指令を生成し、生成した動作指令を工作機械２０に送信する。これにより、数値制御装置１０は、工作機械２０の動作を制御する。
数値制御装置１０は、当業者にとって公知の数値制御装置であり、制御情報に基づいて動作指令を生成し、生成した動作指令を工作機械２０に送信する。これにより、数値制御装置１０は、工作機械２０の動作を制御する。
【0018】
図１に示すように、数値制御装置１０は、記憶部１００、及び制御部２００を有する。さらに、制御部２００は、測定手段２１０、動的補正パラメータ算出手段２２０、指令解析手段２３０、静的補正手段２４０、動的補正手段２５０、補間手段２６０、Ｘ軸用加減速制御手段２７０、Ｙ軸用加減速制御手段２８０、及びＺ軸用加減速制御手段２９０を有する。
【0019】
記憶部１００は、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）やＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）等であり、静的誤差データ１１０、及び動的補正パラメータデータ１２０を記憶する。
【0020】
静的誤差データ１１０は、例えば、後述する静的補正手段２４０による工作機械２０に作用する静的誤差を補正するため、予め測定された静的誤差である。
【0021】
動的補正パラメータデータ１２０は、例えば、後述する動的補正パラメータ算出手段２２０により算出された動的補正パラメータである。
【0022】
制御部２００は、ＣＰＵ、ＲＯＭ、ＲＡＭ、ＣＭＯＳメモリ等を有し、これらはバスを介して相互に通信可能に構成される、当業者にとって公知のものである。
ＣＰＵは数値制御装置１０を全体的に制御するプロセッサである。ＣＰＵは、ＲＯＭに格納されたシステムプログラム及びアプリケーションプログラムを、バスを介して読み出し、前記システムプログラム及びアプリケーションプログラムに従って数値制御装置１０全体を制御する。これにより、図１に示すように、制御部２００が、測定手段２１０、動的補正パラメータ算出手段２２０、指令解析手段２３０、静的補正手段２４０、動的補正手段２５０、補間手段２６０、Ｘ軸用加減速制御手段２７０、Ｙ軸用加減速制御手段２８０、及びＺ軸用加減速制御手段２９０の機能を実現するように構成される。ＲＡＭには一時的な計算データや表示データ等の各種データが格納される。ＣＭＯＳメモリは図示しないバッテリでバックアップされ、数値制御装置１０の電源がオフされても記憶状態が保持される不揮発性メモリとして構成される。
ＣＰＵは数値制御装置１０を全体的に制御するプロセッサである。ＣＰＵは、ＲＯＭに格納されたシステムプログラム及びアプリケーションプログラムを、バスを介して読み出し、前記システムプログラム及びアプリケーションプログラムに従って数値制御装置１０全体を制御する。これにより、図１に示すように、制御部２００が、測定手段２１０、動的補正パラメータ算出手段２２０、指令解析手段２３０、静的補正手段２４０、動的補正手段２５０、補間手段２６０、Ｘ軸用加減速制御手段２７０、Ｙ軸用加減速制御手段２８０、及びＺ軸用加減速制御手段２９０の機能を実現するように構成される。ＲＡＭには一時的な計算データや表示データ等の各種データが格納される。ＣＭＯＳメモリは図示しないバッテリでバックアップされ、数値制御装置１０の電源がオフされても記憶状態が保持される不揮発性メモリとして構成される。
【0023】
測定手段２１０は、例えば、後述する指令解析手段２３０が解析した加工プログラムの切削指令に基づいて工作機械２０が切削したテストワーク（図示しない）の形状を、工作機械２０に含まれる非接触式プローブ（図示しない）で機上測定させる。測定手段２１０は、測定されたテストワーク（図示しない）の形状を示す測定データを工作機械２０から取得する。
なお、機上で形状を測定することで、切削時と測定時とで同じように静的誤差が生じ、結果として機上測定では静的誤差が相殺され、動的誤差のみを測定することができる。
なお、機上で形状を測定することで、切削時と測定時とで同じように静的誤差が生じ、結果として機上測定では静的誤差が相殺され、動的誤差のみを測定することができる。
【0024】
動的補正パラメータ算出手段２２０は、テストワーク（図示しない）を切削した切削指令が示す指令形状と測定手段２１０により取得された測定データとに基づいて、動的誤差を補正する動的補正パラメータを算出する。動的補正パラメータ算出手段２２０は、算出した動的補正パラメータを記憶部１００の動的補正パラメータデータ１２０に記憶する。動的補正パラメータ算出手段２２０の動作については、後述する。
【0025】
指令解析手段２３０は、加工プログラムからＸ軸、Ｙ軸、Ｚ軸の移動の指令を含むブロックを逐次読みだして解析し、解析結果に基づいて各軸の移動の指令座標値を含む切削指令を生成する。
【0026】
静的補正手段２４０は、静的誤差データ１１０から静的誤差を読み出し、読み出した静的誤差に基づいて、指令解析手段２３０により生成された切削指令の指令座標値を補正する。
【0027】
動的補正手段２５０は、動的補正パラメータ算出手段２２０により算出された動的補正パラメータを動的補正パラメータデータ１２０から読み出し、読み出した動的補正パラメータに基づいて、切削指令の指令座標値に対して動的誤差を補正する。
【0028】
補間手段２６０は、動的補正手段２５０から出力される切削指令により指令される移動指令に基づいて、指令経路上の点を補間周期で補間計算した補間データを生成する。
【0029】
Ｘ軸用加減速制御手段２７０は、補間手段２６０から出力される補間データに基づいて、加減速処理を行い補間周期毎のＸ軸の加工速度を計算し、計算した加工速度に応じたパルスを工作機械２０のＸ軸サーボモータ３１に出力する。
【0030】
Ｙ軸用加減速制御手段２８０は、補間手段２６０から出力される補間データに基づいて、加減速処理を行い補間周期毎のＹ軸の加工速度を計算し、計算した加工速度に応じたパルスを工作機械２０のＹ軸サーボモータ３２に出力する。
【0031】
Ｚ軸用加減速制御手段２９０は、補間手段２６０から出力される補間データに基づいて、加減速処理を行い補間周期毎のＺ軸の加工速度を計算し、計算した加工速度に応じたパルスを工作機械２０のＺ軸サーボモータ３３に出力する。
【0032】
次に、動的補正パラメータ算出手段２２０による動的補正パラメータの算出について説明する。なお、工作機械２０における動的誤差には、直角度誤差と、工具２５のたわみによる誤差とがある。
図３Ａは、直角度誤差の一例を示す図である。図３Ｂは、工具２５のたわみによる誤差の一例を示す図である。
図３Ａに示すように、直角度誤差は、例えば、超精密加工機等で各軸間の剛性が低い場合で、テーブル２１と支柱２２（１）、２２（２）とのつなぎ目、支柱２２（１）、２２（２）と支柱２３とのつなぎ目、及び支柱２３と主軸頭２４とのつなぎ目等で発生する。
一方、図３Ｂに示すように、工具のたわみによる誤差は、例えば、工作機械２０の剛性は高いが切削点負荷が大きい場合、主軸頭２４に取り付けられた工具２５で発生する。
以下、直角度誤差に対する動的補正パラメータの算出、及び工具２５のたわみに対する動的補正パラメータの算出それぞれについて説明する。
図３Ａは、直角度誤差の一例を示す図である。図３Ｂは、工具２５のたわみによる誤差の一例を示す図である。
図３Ａに示すように、直角度誤差は、例えば、超精密加工機等で各軸間の剛性が低い場合で、テーブル２１と支柱２２（１）、２２（２）とのつなぎ目、支柱２２（１）、２２（２）と支柱２３とのつなぎ目、及び支柱２３と主軸頭２４とのつなぎ目等で発生する。
一方、図３Ｂに示すように、工具のたわみによる誤差は、例えば、工作機械２０の剛性は高いが切削点負荷が大きい場合、主軸頭２４に取り付けられた工具２５で発生する。
以下、直角度誤差に対する動的補正パラメータの算出、及び工具２５のたわみに対する動的補正パラメータの算出それぞれについて説明する。
【0033】
＜直角度誤差に対する動的補正パラメータの算出について＞
数値制御装置１０は、図４に示すように、直角度誤差に対する動的補正パラメータを算出するために、工作機械２０に対して、Ｚ軸方向の高さを一定にしてＸＹ平面で治具４０に固定されたテストワーク５０に半径Ｒ０の穴を切削させる。なお、テストワーク５０の切削は、例えば、実際の商品のワークの切削加工時と同じ工具２５、同じ材質のワーク等の加工条件で行われるものとする。そうすることで、算出される動的補正パラメータを実際の商品のワークの切削加工時に適用することができる。
そして、数値制御装置１０は、切削された穴を接触式プローブ（図示しない）で工作機械２０に機上測定させる。数値制御装置１０の測定手段２１０は、機上測定されたテストワーク５０の形状を示す測定データを工作機械２０から取得する。
このように、機上で形状を測定することで、切削時と測定時とで同じように静的誤差が生じ、結果として機上測定では静的誤差が相殺され、動的誤差のみを測定することができる。
数値制御装置１０は、図４に示すように、直角度誤差に対する動的補正パラメータを算出するために、工作機械２０に対して、Ｚ軸方向の高さを一定にしてＸＹ平面で治具４０に固定されたテストワーク５０に半径Ｒ０の穴を切削させる。なお、テストワーク５０の切削は、例えば、実際の商品のワークの切削加工時と同じ工具２５、同じ材質のワーク等の加工条件で行われるものとする。そうすることで、算出される動的補正パラメータを実際の商品のワークの切削加工時に適用することができる。
そして、数値制御装置１０は、切削された穴を接触式プローブ（図示しない）で工作機械２０に機上測定させる。数値制御装置１０の測定手段２１０は、機上測定されたテストワーク５０の形状を示す測定データを工作機械２０から取得する。
このように、機上で形状を測定することで、切削時と測定時とで同じように静的誤差が生じ、結果として機上測定では静的誤差が相殺され、動的誤差のみを測定することができる。
【0034】
動的補正パラメータ算出手段２２０は、テストワーク５０を切削した切削指令が示す指令形状と取得された測定データとに基づいて、動的補正パラメータを算出する。
具体的には、動的補正パラメータ算出手段２２０は、図５に示すように、実線で示す測定データにベストフィットする楕円（（ｘ／Ｒｘ）２＋（ｙ／Ｒｙ）２＝１）を、例えば、最小二乗法等で算出する。図５の破線で示す円は、半径Ｒ０の穴の指令形状を示す。
具体的には、動的補正パラメータ算出手段２２０は、図５に示すように、実線で示す測定データにベストフィットする楕円（（ｘ／Ｒｘ）２＋（ｙ／Ｒｙ）２＝１）を、例えば、最小二乗法等で算出する。図５の破線で示す円は、半径Ｒ０の穴の指令形状を示す。
【0035】
なお、テストワーク５０の切削時における切削点負荷は円の直径方向だけでなく、切込み量によっては円の円周方向にもかかる。これにより、図５に示すように、ベストフィットした楕円は角度α傾く。換言すれば、切削指令が示す移動指令位置ＲｎがＲ０（ｃｏｓθ，ｓｉｎθ）とする場合、測定データが示す加工形状位置Ｒａは（Ｒｘｃｏｓ（θ－α），Ｒｙｓｉｎ（θ－α）））となる。
そして、加工形状位置Ｒａの角度θが０度及び９０度のときのズレ量δ０及びδ９０は数１式のように表される。ここで、ズレ量δ０、δ９０の向きは各々Ｘ軸平行、Ｙ軸平行としている。
数１式のズレ量δを用いることで、動的補正パラメータＷｚｘ、Ｗｚｙは数２式のように表される。
ここで、図６に示すように、Ｈｚはテーブル（定盤）２１から主軸頭２４のＸ軸までの高さを示し、Ｈｗはテーブル（定盤）２１からテストワーク５０までの高さを示す。
そして、加工形状位置Ｒａの角度θが０度及び９０度のときのズレ量δ０及びδ９０は数１式のように表される。ここで、ズレ量δ０、δ９０の向きは各々Ｘ軸平行、Ｙ軸平行としている。
【数1】
【数2】
【0036】
なお、テストワーク５０の切削時における切削点負荷が理想的に円の直径方向のみの場合、ベストフィットする楕円は傾かず、角度α＝０となる。
そうすると、切削点負荷で直角度誤差が発生するとした場合、図７に示すように、数２式の動的補正パラメータＷｚｘ、Ｗｚｙは、数３式のように表される。すなわち、直角度誤差は、各軸の直角がずれることよって発生する誤差である。
そうすると、切削点負荷で直角度誤差が発生するとした場合、図７に示すように、数２式の動的補正パラメータＷｚｘ、Ｗｚｙは、数３式のように表される。すなわち、直角度誤差は、各軸の直角がずれることよって発生する誤差である。
【数3】
【0037】
この場合、ズレ量δは、数４式のように表される。
そして、切削点負荷Ｆ（＝（Ｆｘ，Ｆｙ））の方向は直径方向となることから、Ｆ／｜Ｆ｜＝－（ｃｏｓθ，ｓｉｎθ）となり、ズレ量δは切削点負荷Ｆを用いて数５式のように表される。
Ｈｚ－Ｈｗは、加工の高さｚである。
なお、実際の商品のワーク加工時における切削点負荷の方向、すなわち、切削点負荷が円の直径方向だけでなく円周方向にもかかる場合の方向は、公知の手法（例えば、松村隆、「切削シミュレーションの現状と課題」、精密工学会誌、Ｖｏｌ．８０、Ｎｏ．９、２０１４）を用いて、商品のＣＡＤモデルと加工プランとから加工シミュレーションを行うことで算出されてもよい。
また、実際の商品のワーク加工時における切削点負荷の方向は、第３実施形態で後述するように、例えば、Ｘ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３のトルクから見積もられるようにしてもよい。あるいは、実際の商品のワーク加工時における切削点負荷Ｆの方向は、工具２５に取り付けられたセンサを用いて検出するようにしてもよい。
【数4】
【数5】
なお、実際の商品のワーク加工時における切削点負荷の方向、すなわち、切削点負荷が円の直径方向だけでなく円周方向にもかかる場合の方向は、公知の手法（例えば、松村隆、「切削シミュレーションの現状と課題」、精密工学会誌、Ｖｏｌ．８０、Ｎｏ．９、２０１４）を用いて、商品のＣＡＤモデルと加工プランとから加工シミュレーションを行うことで算出されてもよい。
また、実際の商品のワーク加工時における切削点負荷の方向は、第３実施形態で後述するように、例えば、Ｘ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３のトルクから見積もられるようにしてもよい。あるいは、実際の商品のワーク加工時における切削点負荷Ｆの方向は、工具２５に取り付けられたセンサを用いて検出するようにしてもよい。
【0038】
そして、動的補正パラメータ算出手段２２０は、数１式及び数２式を用いて、動的補正パラメータＷｚｘ、Ｗｚｙを算出し、算出した動的補正パラメータＷｚｘ、Ｗｚｙを動的補正パラメータデータ１２０に記憶する。
【0039】
その後、動的補正手段２５０は、実際の商品のワークを切削加工する場合、動的補正パラメータデータ１２０から動的補正パラメータＷｚｘ、Ｗｚｙを読み出す。動的補正手段２５０は、読み出した動的補正パラメータＷｚｘ、Ｗｚｙに基づいて実際の商品となるワークに対する切削指令の指令座標値に対して直角度誤差の動的誤差を補正することができる。
【0040】
＜工具２５のたわみによる誤差に対する動的補正パラメータの算出について＞
数値制御装置１０は、直角度誤差の場合と同様に、工具２５のたわみによる誤差に対する動的補正パラメータを算出するために、工作機械２０に対して、Ｚ軸方向の高さを一定にしてＸＹ平面で治具４０に固定されたテストワーク５０に半径Ｒ０の穴を切削させる。なお、テストワーク５０の切削は、例えば、実際の商品のワークの切削加工時と同じ工具２５、同じ材質のワーク等の加工条件で行われるものとする。そうすることで、算出される動的補正パラメータを実際の商品のワークの切削加工時に適用することができる。
そして、数値制御装置１０は、切削された穴を接触式プローブ（図示しない）で工作機械２０に機上測定させる。数値制御装置１０の測定手段２１０は、機上測定されたテストワーク５０の形状を示す測定データを工作機械２０から取得する。
数値制御装置１０は、直角度誤差の場合と同様に、工具２５のたわみによる誤差に対する動的補正パラメータを算出するために、工作機械２０に対して、Ｚ軸方向の高さを一定にしてＸＹ平面で治具４０に固定されたテストワーク５０に半径Ｒ０の穴を切削させる。なお、テストワーク５０の切削は、例えば、実際の商品のワークの切削加工時と同じ工具２５、同じ材質のワーク等の加工条件で行われるものとする。そうすることで、算出される動的補正パラメータを実際の商品のワークの切削加工時に適用することができる。
そして、数値制御装置１０は、切削された穴を接触式プローブ（図示しない）で工作機械２０に機上測定させる。数値制御装置１０の測定手段２１０は、機上測定されたテストワーク５０の形状を示す測定データを工作機械２０から取得する。
【0041】
動的補正パラメータ算出手段２２０は、テストワーク５０を切削した切削指令が示す指令形状と取得された測定データとに基づいて、動的補正パラメータを算出する。
具体的には、工具２５のたわみには異方位性がないことから、動的補正パラメータ算出手段２２０は、図８に示すように、実線で示す測定データにベストフィットする円（ｘ２＋ｙ２＝Ｒｔ ２）を、例えば、最小二乗法等で算出する。なお、図８の破線で示す円は、図５の場合と同様に、半径Ｒ０の穴の指令形状を示す。
ただし、図９に示すように、切削点負荷Ｆの方向は、工具２５の回転と当該回転に対するテストワーク５０からの反作用により、テストワーク５０の法線に対して角度β傾く。このため、指令した位置と実際の切削位置とは、図８に示すようにずれてしまう。そこで、動的補正パラメータ算出手段２２０は、公知の手法（例えば、谷口和雄、「金属切削機構の力学的解析（第３報）」、精密機械、２９巻、第３号、１９６３）を用いて、切削点負荷Ｆの方向の傾きβを算出する。
これにより、切削点負荷で発生する工具２５のたわみ量の係数（動的補正パラメータ）Ｗｔは、数６式のように表される。
Ｌｔは工具長を示す。そして、ズレ量δは、数７式のように表される。
具体的には、工具２５のたわみには異方位性がないことから、動的補正パラメータ算出手段２２０は、図８に示すように、実線で示す測定データにベストフィットする円（ｘ２＋ｙ２＝Ｒｔ ２）を、例えば、最小二乗法等で算出する。なお、図８の破線で示す円は、図５の場合と同様に、半径Ｒ０の穴の指令形状を示す。
ただし、図９に示すように、切削点負荷Ｆの方向は、工具２５の回転と当該回転に対するテストワーク５０からの反作用により、テストワーク５０の法線に対して角度β傾く。このため、指令した位置と実際の切削位置とは、図８に示すようにずれてしまう。そこで、動的補正パラメータ算出手段２２０は、公知の手法（例えば、谷口和雄、「金属切削機構の力学的解析（第３報）」、精密機械、２９巻、第３号、１９６３）を用いて、切削点負荷Ｆの方向の傾きβを算出する。
これにより、切削点負荷で発生する工具２５のたわみ量の係数（動的補正パラメータ）Ｗｔは、数６式のように表される。
【数6】
【数7】
【0042】
動的補正パラメータ算出手段２２０は、数６式を用いて、動的補正パラメータＷｔを算出し、算出した動的補正パラメータＷｔを動的補正パラメータデータ１２０に記憶する。
【0043】
その後、動的補正手段２５０は、実際の商品のワークを切削する場合、動的補正パラメータデータ１２０から動的補正パラメータＷｔを読み出す。動的補正手段２５０は、読み出した動的補正パラメータＷｔに基づいて実際の商品となるワークに対する切削指令の指令座標値に対して工具２５のたわみによる動的誤差を補正する。
これにより、図１０に示すように、実線で示す実際の位置にいる工具２５を破線で示す指令座標値の位置に移動させることができる。
これにより、図１０に示すように、実線で示す実際の位置にいる工具２５を破線で示す指令座標値の位置に移動させることができる。
【0044】
以上により、第１実施形態の数値制御装置１０は、切削させたテストワーク５０の形状を機上測定させ、測定されたテストワーク５０の形状の測定データを取得する。数値制御装置１０は、切削指令の指令形状と取得された測定データとに基づいて動的誤差を補正する動的補正パラメータを算出する。数値制御装置１０は、算出された動的補正パラメータに基づいて指令座標値に対して動的誤差を補正することができる。
これにより、数値制御装置１０は、工具のたわみによる圧力を検出するセンサを用いることなく、動的誤差を精度良く補正することができる。
また、数値制御装置１０は、動的補正パラメータを算出することにより、センサ等で測定する切削点負荷と、工具のたわみとの関係を予め調べる必要がない。
以上、第１実施形態について説明した。
これにより、数値制御装置１０は、工具のたわみによる圧力を検出するセンサを用いることなく、動的誤差を精度良く補正することができる。
また、数値制御装置１０は、動的補正パラメータを算出することにより、センサ等で測定する切削点負荷と、工具のたわみとの関係を予め調べる必要がない。
以上、第１実施形態について説明した。
【0045】
＜第１実施形態の変形例＞
上述の第１実施形態では、数値制御装置１０は、加工プログラムを解析することにより生成された切削指令の指令座標値に対して動的誤差を補正したがこれに限定されない。例えば、数値制御装置１０は、工作機械２０のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３に出力するパルスに対して、動的誤差を補正する補正パルスを加算してもよい。
上述の第１実施形態では、数値制御装置１０は、加工プログラムを解析することにより生成された切削指令の指令座標値に対して動的誤差を補正したがこれに限定されない。例えば、数値制御装置１０は、工作機械２０のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３に出力するパルスに対して、動的誤差を補正する補正パルスを加算してもよい。
【0046】
図１１は、第１実施形態に係る数値制御装置１０の一例を示す図である。
図１１に示すように、静的補正手段２４０－１及び動的補正手段２５０－１は、Ｘ軸用加減速制御手段２７０、Ｙ軸用加減速制御手段２８０、Ｚ軸用加減速制御手段２９０の後に配置される。そして、静的補正手段２４０－１及び動的補正手段２５０－１は、Ｘ軸用加減速制御手段２７０、Ｙ軸用加減速制御手段２８０、Ｚ軸用加減速制御手段２９０の各々から出力されるパルスに対して、静的誤差を補正する補正パルス及び動的誤差を補正する補正パルスを加算することで、静的誤差及び動的誤差の補正を行う。なお、静的補正手段２４０－１及び動的補正手段２５０－１は、静的誤差の補正パルス及び動的誤差の補正パルスを加算する点を除き、図１の静的補正手段２４０及び動的補正手段２５０と同様の動作を行う。
図１１に示すように、静的補正手段２４０－１及び動的補正手段２５０－１は、Ｘ軸用加減速制御手段２７０、Ｙ軸用加減速制御手段２８０、Ｚ軸用加減速制御手段２９０の後に配置される。そして、静的補正手段２４０－１及び動的補正手段２５０－１は、Ｘ軸用加減速制御手段２７０、Ｙ軸用加減速制御手段２８０、Ｚ軸用加減速制御手段２９０の各々から出力されるパルスに対して、静的誤差を補正する補正パルス及び動的誤差を補正する補正パルスを加算することで、静的誤差及び動的誤差の補正を行う。なお、静的補正手段２４０－１及び動的補正手段２５０－１は、静的誤差の補正パルス及び動的誤差の補正パルスを加算する点を除き、図１の静的補正手段２４０及び動的補正手段２５０と同様の動作を行う。
【0047】
また、数値制御装置１０は、工具２５の工具径に動的補正の補正量を組み込んでもよい。
図１２は、第１実施形態に係る数値制御装置１０の一例を示す図である。
図１２に示す工具径動的補正手段２５０－２は、工具２５のたわみによる誤差に対する動的補正を行う。この場合、工具径動的補正手段２５０－２は、数６式により算出された動的補正パラメータＷｔを用いて算出される補正量を、数８式に基づいて工具２５の工具径ｒｔに組み込んで工具径補正量ｒｔ’を算出してもよい。
なお、β’は、実際の切削時の切削点負荷Ｆの方向の傾きを示す。そして、傾きβ’がテストワーク５０の切削時の切削点負荷Ｆの方向の傾きβと同じ場合、図１３に示すように、数８式は、数９式のように表される。
なお、図１２の数値制御装置１０は、静的補正手段２４０を省略したが、静的補正手段２４０を有してもよい。
図１２は、第１実施形態に係る数値制御装置１０の一例を示す図である。
図１２に示す工具径動的補正手段２５０－２は、工具２５のたわみによる誤差に対する動的補正を行う。この場合、工具径動的補正手段２５０－２は、数６式により算出された動的補正パラメータＷｔを用いて算出される補正量を、数８式に基づいて工具２５の工具径ｒｔに組み込んで工具径補正量ｒｔ’を算出してもよい。
【数8】
【数9】
【0048】
＜第２実施形態＞
次に、第２実施形態について説明する。第１実施形態では、テストワーク５０の加工時と、実際の商品のワークの加工時とで負荷の大きさに差があった場合、その差分は考慮されない。この場合、例えば、切削負荷の大きさの差分が一定以下（数値制御装置１０が正常に機能する範囲）になるように、特開２０１６－１３７５５７号公報等の公知の手法を用いて、自動で切削条件を変更することができるが、切削条件を変更することが余儀なくされる。そこで、第２実施形態では、数値制御装置１０は、第１実施形態の機能に加えて、テストワークを互いに異なる複数の切削負荷の各々で切削して算出された複数の動的補正パラメータを用いて、任意の切削負荷における動的補正パラメータを補間し、補間した任意の切削負荷における動的補正パラメータに基づいて切削指令の指令座標値を補正する。
次に、第２実施形態について説明する。第１実施形態では、テストワーク５０の加工時と、実際の商品のワークの加工時とで負荷の大きさに差があった場合、その差分は考慮されない。この場合、例えば、切削負荷の大きさの差分が一定以下（数値制御装置１０が正常に機能する範囲）になるように、特開２０１６－１３７５５７号公報等の公知の手法を用いて、自動で切削条件を変更することができるが、切削条件を変更することが余儀なくされる。そこで、第２実施形態では、数値制御装置１０は、第１実施形態の機能に加えて、テストワークを互いに異なる複数の切削負荷の各々で切削して算出された複数の動的補正パラメータを用いて、任意の切削負荷における動的補正パラメータを補間し、補間した任意の切削負荷における動的補正パラメータに基づいて切削指令の指令座標値を補正する。
【0049】
これにより、第２実施形態の数値制御装置１０は、テストワークの加工時の切削点負荷の大きさと、実際の商品のワークの加工時の切削点負荷の大きさが異なる場合でも、適切に動的誤差を補正することができる。
以下に、第２実施形態について説明する。
以下に、第２実施形態について説明する。
【0050】
図１４は、第２実施形態に係る数値制御装置のうち第１実施形態の数値制御装置１０に追加される機能的構成例を示す機能ブロック図である。図１の数値制御装置１０の要素と同様の機能を有する要素については、同じ符号を付し、詳細な説明は省略する。
以下、切削負荷と単位時間当たりの切削体積とは、強い相関関係があることから、切削負荷に替えて、単位時間当たりの切削体積Ｖを変数にして説明する。
以下、切削負荷と単位時間当たりの切削体積とは、強い相関関係があることから、切削負荷に替えて、単位時間当たりの切削体積Ｖを変数にして説明する。
【0051】
図１４に示すように、制御部２００は、切削負荷算出手段３００、及び動的補正パラメータ補間手段３１０をさらに備える。これらの各機能部は、制御部２００が制御部２００のＲＯＭ（図示しない）に格納されたシステムプログラム及びアプリケーションプログラムを実行することにより実現される。
【0052】
なお、測定手段２１０は、例えば、指令解析手段２３０が解析した加工プログラムの切削指令に基づいて工作機械２０が互いに異なる複数の切削負荷、すなわち複数の切削体積Ｖ１からＶＮの各々で切削したテストワーク５０の各々の形状を、工作機械２０の非接触式プローブ（図示しない）で機上測定させる。Ｎは、２以上の整数である。測定手段２１０は、測定された各テストワーク５０の形状の測定データを工作機械２０から取得する。
【0053】
動的補正パラメータ算出手段２２０は、直角度誤差の場合、複数の切削体積Ｖ１からＶＮの各々でテストワーク５０を切削した時の切削指令の指令形状と、測定手段２１０により取得された測定データとに基づき、数１式及び数２式を用いて動的補正パラメータ｛Ｗｚｘ（Ｖｉ）｜１≦ｉ≦Ｎ、Ｎは２以上の整数｝、｛Ｗｚｙ（Ｖｉ）｜１≦ｉ≦Ｎ｝を算出する。
また、動的補正パラメータ算出手段２２０は、工具２５のたわみによる誤差の場合、複数の切削体積Ｖ１からＶＮの各々でテストワーク５０を切削した時の切削指令の指令形状と、測定手段２１０により取得された測定データとに基づき、数６式から動的補正パラメータ｛Ｗｔ（Ｖｉ）｜１≦ｉ≦Ｎ｝を算出する。
そして、動的補正パラメータ算出手段２２０は、各切削点負荷の直角度誤差の動的補正パラメータ｛Ｗｚｘ（Ｖｉ）｜１≦ｉ≦Ｎ｝、｛Ｗｚｙ（Ｖｉ）｜１≦ｉ≦Ｎ｝、及び工具２５のたわみによる誤差の動的補正パラメータ｛Ｗｔ（Ｖｉ）｜１≦ｉ≦Ｎ｝を、切削体積Ｖ１からＶＮの各々と対応付けして動的補正パラメータデータ１２０に記憶する。
また、動的補正パラメータ算出手段２２０は、工具２５のたわみによる誤差の場合、複数の切削体積Ｖ１からＶＮの各々でテストワーク５０を切削した時の切削指令の指令形状と、測定手段２１０により取得された測定データとに基づき、数６式から動的補正パラメータ｛Ｗｔ（Ｖｉ）｜１≦ｉ≦Ｎ｝を算出する。
そして、動的補正パラメータ算出手段２２０は、各切削点負荷の直角度誤差の動的補正パラメータ｛Ｗｚｘ（Ｖｉ）｜１≦ｉ≦Ｎ｝、｛Ｗｚｙ（Ｖｉ）｜１≦ｉ≦Ｎ｝、及び工具２５のたわみによる誤差の動的補正パラメータ｛Ｗｔ（Ｖｉ）｜１≦ｉ≦Ｎ｝を、切削体積Ｖ１からＶＮの各々と対応付けして動的補正パラメータデータ１２０に記憶する。
【0054】
切削負荷算出手段３００は、指令解析手段２３０が解析した加工プログラムの加工条件に基づいて、単位時間当たりの切削体積Ｖを算出する。なお、加工条件から切削体積Ｖを算出する方法は公知の方法を用いることができ、説明を省略する。
【0055】
動的補正パラメータ補間手段３１０は、動的補正パラメータデータ１２０から複数の動的補正パラメータ｛Ｗｚｘ（Ｖｉ）｜１≦ｉ≦Ｎ｝、｛Ｗｚｙ（Ｖｉ）｜１≦ｉ≦Ｎ｝（又は複数の動的補正パラメータ｛Ｗｔ（Ｖｉ）｜１≦ｉ≦Ｎ｝）を読み出す。動的補正パラメータ補間手段３１０は、読み出した複数の動的補正パラメータ｛Ｗｚｘ（Ｖｉ）｜１≦ｉ≦Ｎ｝、｛Ｗｚｙ（Ｖｉ）｜１≦ｉ≦Ｎ｝（又は複数の動的補正パラメータ｛Ｗｔ（Ｖｉ）｜１≦ｉ≦Ｎ｝）を用いて、切削負荷算出手段３００により算出された切削体積Ｖにおける動的補正パラメータＷｚｘ（Ｖ）、Ｗｚｙ（Ｖ）（又は動的補正パラメータＷｔ（Ｖ））を補間する。
以下、直角度誤差に対する動的補正パラメータの補間、及び工具２５のたわみに対する動的補正パラメータの補間それぞれについて説明する。
以下、直角度誤差に対する動的補正パラメータの補間、及び工具２５のたわみに対する動的補正パラメータの補間それぞれについて説明する。
【0056】
＜直角度誤差に対する動的補正パラメータの補間について＞
図１５Ａ及び図１５Ｂは、動的補正パラメータ補間手段３１０の補間処理を説明する一例を示す図である。なお、図１５Ａは、動的補正パラメータＷｚｘ（Ｖ）の場合を示し、図１５Ｂは、動的補正パラメータＷｚｙ（Ｖ）の場合を示す。また、図１５Ａ及び図１５Ｂは、Ｎ＝２の場合を示すが、Ｎが３以上の場合も同様である。
図１５Ａ及び図１５Ｂは、動的補正パラメータ補間手段３１０の補間処理を説明する一例を示す図である。なお、図１５Ａは、動的補正パラメータＷｚｘ（Ｖ）の場合を示し、図１５Ｂは、動的補正パラメータＷｚｙ（Ｖ）の場合を示す。また、図１５Ａ及び図１５Ｂは、Ｎ＝２の場合を示すが、Ｎが３以上の場合も同様である。
【0057】
動的補正パラメータ補間手段３１０は、例えば、図１５Ａに示すように、動的補正パラメータＷｚｘ（Ｖ１）、Ｗｚｘ（Ｖ２）を線形補間して、切削負荷算出手段３００により算出された切削体積Ｖにおける動的補正パラメータＷｚｘ（Ｖ）を算出する。また、動的補正パラメータ補間手段３１０は、図１５Ｂに示すように、動的補正パラメータＷｚｙ（Ｖ１）、Ｗｚｙ（Ｖ２）を線形補間して、切削負荷算出手段３００により算出された単位時間当たりの切削体積Ｖにおける動的補正パラメータＷｚｙ（Ｖ）を算出する。動的補正パラメータ補間手段３１０は、算出した動的補正パラメータＷｚｘ（Ｖ）、Ｗｚｙ（Ｖ）を動的補正手段２５０に出力する。
なお、動的補正パラメータ補間手段３１０は、Ｎが３以上の場合、図１６に示しように、Ｍ次関数でベストフィットしてもよく、機械学習で求めてもよい（Ｍは２以上の整数）。
なお、動的補正パラメータ補間手段３１０は、Ｎが３以上の場合、図１６に示しように、Ｍ次関数でベストフィットしてもよく、機械学習で求めてもよい（Ｍは２以上の整数）。
【0058】
そして、動的補正手段２５０は、動的補正パラメータ補間手段３１０により算出された動的補正パラメータＷｚｘ（Ｖ）、Ｗｚｙ（Ｖ）に基づいて実際の商品となるワークに対する切削指令の指令座標値に対して直角度誤差の動的誤差を補正する。
【0059】
＜工具２５のたわみによる誤差に対する動的補正パラメータの補間について＞
動的補正パラメータ補間手段３１０は、直角度誤差の場合と同様に、例えば、２つの切削体積Ｖ１、Ｖ２の各々における動的補正パラメータＷｔ（Ｖ１）、Ｗｔ（Ｖ２）を線形補間して、切削負荷算出手段３００により算出された切削体積Ｖにおける動的補正パラメータＷｔ（Ｖ）を算出する。動的補正パラメータ補間手段３１０は、算出した動的補正パラメータＷｔ（Ｖ）を動的補正手段２５０に出力する。
なお、動的補正パラメータ補間手段３１０は、Ｎが３以上の場合、Ｍ次関数でベストフィットしてもよく、機械学習で求めてもよい（Ｍは２以上の整数）。
動的補正パラメータ補間手段３１０は、直角度誤差の場合と同様に、例えば、２つの切削体積Ｖ１、Ｖ２の各々における動的補正パラメータＷｔ（Ｖ１）、Ｗｔ（Ｖ２）を線形補間して、切削負荷算出手段３００により算出された切削体積Ｖにおける動的補正パラメータＷｔ（Ｖ）を算出する。動的補正パラメータ補間手段３１０は、算出した動的補正パラメータＷｔ（Ｖ）を動的補正手段２５０に出力する。
なお、動的補正パラメータ補間手段３１０は、Ｎが３以上の場合、Ｍ次関数でベストフィットしてもよく、機械学習で求めてもよい（Ｍは２以上の整数）。
【0060】
そして、動的補正手段２５０は、動的補正パラメータ補間手段３１０により算出された動的補正パラメータＷｔ（Ｖ）に基づいて実際の商品となるワークに対する切削指令の指令座標値に対して工具２５のたわみによる誤差の動的誤差を補正する。
【0061】
以上により、第２実施形態の数値制御装置１０は、互いに異なる複数の切削体積（切削負荷）の各々で切削されたテストワークの形状の測定データを取得し、複数の切削体積の各々における指令形状と測定データとに基づいて、切削体積毎の動的補正パラメータを算出する。数値制御装置１０は、算出された切削体積毎の動的補正パラメータを用いて、任意の切削体積における動的補正パラメータを補間し、補間した任意の切削体積における動的補正パラメータに基づいて切削指令の指令座標値を補正する。
これにより、数値制御装置１０は、工具のたわみによる圧力を検出するセンサを用いることなく、動的誤差を精度良く補正することができる。
また、数値制御装置１０は、テストワークの加工時の切削負荷の大きさと、実際の商品のワークの加工時の切削負荷の大きさが異なる場合でも、適切に動的誤差を補正することができる。
以上、第２実施形態について説明した。
これにより、数値制御装置１０は、工具のたわみによる圧力を検出するセンサを用いることなく、動的誤差を精度良く補正することができる。
また、数値制御装置１０は、テストワークの加工時の切削負荷の大きさと、実際の商品のワークの加工時の切削負荷の大きさが異なる場合でも、適切に動的誤差を補正することができる。
以上、第２実施形態について説明した。
【0062】
＜第２実施形態の変形例＞
上述の第２実施形態では、数値制御装置１０は、図１の数値制御装置１０に図１４の構成を追加することで、加工プログラムを解析することにより生成された切削指令の指令座標値に対して動的誤差を補正したがこれに限定されない。例えば、数値制御装置１０は、図１１の数値制御装置１０に図１４の構成を追加することで、工作機械２０のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３に出力するパルスに対して、動的誤差を補正する補正パルスを加算してもよい。
あるいは、数値制御装置１０は、図１２の数値制御装置１０に図１４の構成を追加することで、工具２５の工具径に動的補正の補正量を組み込み、動的誤差を補正してもよい。
上述の第２実施形態では、数値制御装置１０は、図１の数値制御装置１０に図１４の構成を追加することで、加工プログラムを解析することにより生成された切削指令の指令座標値に対して動的誤差を補正したがこれに限定されない。例えば、数値制御装置１０は、図１１の数値制御装置１０に図１４の構成を追加することで、工作機械２０のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３に出力するパルスに対して、動的誤差を補正する補正パルスを加算してもよい。
あるいは、数値制御装置１０は、図１２の数値制御装置１０に図１４の構成を追加することで、工具２５の工具径に動的補正の補正量を組み込み、動的誤差を補正してもよい。
【0063】
＜第３実施形態＞
次に、第３実施形態について説明する。第３実施形態では、数値制御装置１０は、第１実施形態の機能に加えて、テストワークの切削時の負荷電流を取得し、取得した負荷電流、指令形状と測定データとのズレ量、及び動的補正パラメータに基づいて各軸の負荷電流と切削点負荷との関係を示す関係パラメータを算出し、関係パラメータと負荷電流とに基づいて切削点負荷を算出する。
次に、第３実施形態について説明する。第３実施形態では、数値制御装置１０は、第１実施形態の機能に加えて、テストワークの切削時の負荷電流を取得し、取得した負荷電流、指令形状と測定データとのズレ量、及び動的補正パラメータに基づいて各軸の負荷電流と切削点負荷との関係を示す関係パラメータを算出し、関係パラメータと負荷電流とに基づいて切削点負荷を算出する。
【0064】
これにより、第３実施形態の数値制御装置１０は、ハードウェアに手を加えることなく、追加のハードウェアも不要で切削点負荷を推定することができる。
以下に、第３実施形態について説明する。
以下に、第３実施形態について説明する。
【0065】
図１７は、第３実施形態に係る数値制御装置のうち第１実施形態の数値制御装置１０に追加される機能的構成例を示す機能ブロック図である。図１の数値制御装置１０の要素と同様の機能を有する要素については、同じ符号を付し、詳細な説明は省略する。
【0066】
図１７に示すように、制御部２００は、関係パラメータ算出手段３２０、及び切削点負荷算出手段３３０をさらに備える。これらの各機能部は、制御部２００が制御部２００のＲＯＭ（図示しない）に格納されたシステムプログラム及びアプリケーションプログラムを実行することにより実現される。
【0067】
関係パラメータ算出手段３２０は、テストワーク５０の切削時のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３から負荷電流を取得する。関係パラメータ算出手段３２０は、取得された負荷電流、指令形状と測定データとのズレ量δ、及び動的補正パラメータに基づいて少なくともＸ軸及びＹ軸における負荷電流と切削点負荷との関係を示す関係パラメータを算出する。
具体的には、関係パラメータ算出手段３２０は、Ｘ軸及びＹ軸における負荷電流と切削点負荷との関係を示す関係パラメータを算出するために、図１８Ａに示すように、主軸頭２４に取り付けられた工具２５による切削加工は行わずに速度一定でＸＹ平面を時計周りに円運動させる。工具２５がＺ軸方向に移動しないことから、関係パラメータ算出手段３２０は、工具２５が円運動している時のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２の負荷電流を取得する。この場合、負荷電流Ｉ１(θ)は、（Ｉ１ｘ(θ)，Ｉ１ｙ(θ)）と表される。なお、θは円運動の位相を示す。
具体的には、関係パラメータ算出手段３２０は、Ｘ軸及びＹ軸における負荷電流と切削点負荷との関係を示す関係パラメータを算出するために、図１８Ａに示すように、主軸頭２４に取り付けられた工具２５による切削加工は行わずに速度一定でＸＹ平面を時計周りに円運動させる。工具２５がＺ軸方向に移動しないことから、関係パラメータ算出手段３２０は、工具２５が円運動している時のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２の負荷電流を取得する。この場合、負荷電流Ｉ１(θ)は、（Ｉ１ｘ(θ)，Ｉ１ｙ(θ)）と表される。なお、θは円運動の位相を示す。
【0068】
同様に、関係パラメータ算出手段３２０は、図１８Ｂに示すように、主軸頭２４に取り付けられた工具２５による切削加工は行わずに、時計周りの場合と同じ速度一定でＸＹ平面を反時計周りに円運動させる。工具２５がＺ軸方向に移動しないことから、関係パラメータ算出手段３２０は、工具２５が円運動している時のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２の負荷電流を取得する。この場合、負荷電流Ｉ２(θ)は、（Ｉ２ｘ(θ)，Ｉ２ｙ(θ)）と表される。
ここで、負荷電流Ｉ１(θ)と負荷電流Ｉ２(θ)との加速度は、法線方向で向きが同じになるのに対し、負荷電流Ｉ１(θ)と負荷電流Ｉ２(θ)との速度は、接線方向で向きが逆になる。このことから、Ｉ１(θ)＋Ｉ２(θ)は、速度の影響がキャンセルされ、加速度の影響のみになる。一方、Ｉ１(θ)－Ｉ２(θ)は、加速度の影響がキャンセルされ、速度の影響のみになる。
ここで、負荷電流Ｉ１(θ)と負荷電流Ｉ２(θ)との加速度は、法線方向で向きが同じになるのに対し、負荷電流Ｉ１(θ)と負荷電流Ｉ２(θ)との速度は、接線方向で向きが逆になる。このことから、Ｉ１(θ)＋Ｉ２(θ)は、速度の影響がキャンセルされ、加速度の影響のみになる。一方、Ｉ１(θ)－Ｉ２(θ)は、加速度の影響がキャンセルされ、速度の影響のみになる。
【0069】
そこで、関係パラメータ算出手段３２０は、加速度の影響を示すＩ１(θ)＋Ｉ２(θ)にベストフィットする楕円（（ｘ／Ｒａｘ）２＋（ｙ／Ｒａｙ）２＝１）を、例えば、最小二乗法等で算出する。また、関係パラメータ算出手段３２０は、速度の影響を示すＩ１(θ)－Ｉ２(θ)にベストフィットする楕円（（ｘ／Ｒｖｘ）２＋（ｙ／Ｒｖｙ）２＝１）を、例えば、最小二乗法等で算出する。これにより、Ｉ１(θ)＋Ｉ２(θ)、及びＩ１(θ)－Ｉ２(θ)は、数１０式のように表される。
【数10】
【0070】
一方、Ｘ軸の加速度とトルクの比例係数Ｋｘ、Ｙ軸の加速度とトルクの比例係数Ｋｙ、Ｘ軸の速度とトルクの比例係数Ｌｘ、Ｙ軸の速度とトルクの比例係数Ｌｙとする場合、負荷電流Ｉ１(θ)、Ｉ２(θ)は、数１１式のように表される。なお、ａは加速度を示す。ｖは速度を示す。
そして、数１１式の負荷電流Ｉ１(θ)、Ｉ２(θ)を用いて、Ｉ１(θ)＋Ｉ２(θ)及びＩ１(θ)－Ｉ２(θ)は、数１２式のよう表される。
【数11】
【数12】
【0071】
これにより、Ｘ軸及びＹ軸の加速度とトルクの比例係数（加速度と負荷電流との関係パラメータ）Ｋｘ、Ｋｙ、及びＸ軸及びＹ軸の速度とトルクの比例係数（速度と負荷電流との関係パラメータ）Ｌｘ、Ｌｙは、数１０式の値Ｒａｘ、Ｒａｙ、Ｒｖｘ、Ｒｖｙ、ａ及びｖを用いて数１３式のように表される。
【数13】
【0072】
換言すれば、関係パラメータ算出手段３２０は、負荷電流Ｉ１(θ) 、Ｉ２(θ)、数１０式及び数１３式から加速度と負荷電流との関係パラメータＫｘ、Ｋｙ、及び速度と負荷電流との関係パラメータＬｘ、Ｌｙを算出する。そして、関係パラメータ算出手段３２０は、算出した加速度と負荷電流との関係パラメータＫｘ、Ｋｙ、及び速度と負荷電流との関係パラメータＬｘ、Ｌｙを後述する切削点負荷算出手段３３０に出力する。なお、関係パラメータ算出手段３２０は、算出した加速度と負荷電流との関係パラメータＫｘ、Ｋｙ、及び速度と負荷電流との関係パラメータＬｘ、Ｌｙを記憶部１００に記憶してもよい。
【0073】
次に、数値制御装置１０は、第１実施形態の場合と同様に、直角度誤差に対する動的補正パラメータＷｚｘ、Ｗｚｙ、又は工具２５のたわみによる誤差に対する動的補正パラメータＷｔを算出するために、工作機械２０に対して、Ｚ軸方向の高さを一定にしてＸＹ平面で治具４０に固定されたテストワーク５０に半径Ｒ０の穴を切削させる。
そして、数値制御装置１０は、切削された穴を接触式プローブ（図示しない）で工作機械２０に機上測定させる。数値制御装置１０の測定手段２１０は、機上測定されたテストワーク５０の形状を示す測定データを工作機械２０から取得する。
動的補正パラメータ算出手段２２０は、第１実施形態の場合と同様に、テストワーク５０を切削した切削指令が示す指令形状と取得された測定データとに基づいて、直角度誤差に対する動的補正パラメータＷｚｘ、Ｗｚｙ、又は工具２５のたわみによる誤差に対する動的補正パラメータＷｔを算出する。
そして、数値制御装置１０は、切削された穴を接触式プローブ（図示しない）で工作機械２０に機上測定させる。数値制御装置１０の測定手段２１０は、機上測定されたテストワーク５０の形状を示す測定データを工作機械２０から取得する。
動的補正パラメータ算出手段２２０は、第１実施形態の場合と同様に、テストワーク５０を切削した切削指令が示す指令形状と取得された測定データとに基づいて、直角度誤差に対する動的補正パラメータＷｚｘ、Ｗｚｙ、又は工具２５のたわみによる誤差に対する動的補正パラメータＷｔを算出する。
【0074】
また、関係パラメータ算出手段３２０は、テストワーク５０の切削加工時にＸ軸サーボモータ３１、及びＹ軸サーボモータ３２から負荷電流Ｉｆ（Ｘ）、Ｉｆ（Ｙ）を取得する。関係パラメータ算出手段３２０は、取得された切削加工時の負荷電流Ｉｆ（Ｘ）、Ｉｆ（Ｙ）を用いて、Ｘ軸の負荷電流と切削点負荷Ｆとの比例係数（負荷電流と切削点負荷との関係パラメータ）Ｊｘ、及びＹ軸の負荷電流と切削点負荷Ｆとの比例係数（負荷電流との切削点負荷との関係パラメータ）Ｊｙを算出する。
以下、直角度誤差における負荷電流と切削点負荷との関係パラメータの算出、及び工具２５のたわみにおける負荷電流と切削点負荷との関係パラメータの算出それぞれについて説明する。
以下、直角度誤差における負荷電流と切削点負荷との関係パラメータの算出、及び工具２５のたわみにおける負荷電流と切削点負荷との関係パラメータの算出それぞれについて説明する。
【0075】
＜直角度誤差における負荷電流と切削点負荷との関係パラメータの算出について＞
関係パラメータ算出手段３２０は、テストワーク５０を切削加工したことにより発生する直角度誤差は工作機械２０の各軸間の剛性が低いことによるため、取得された負荷電流Ｉｆ（＝（Ｉｆ（Ｘ）、Ｉｆ（Ｙ））にベストフィットする楕円（（ｘ／Ｉｆｘ）２＋（ｙ／Ｉｆｙ）２＝１）を、例えば、最小二乗法等で算出する。なお、Ｉｆｘは楕円のＸ軸方向の半径である。また、Ｉｆｙは楕円のＹ軸方向の半径である。
ここで、数５式に基づいて、ズレ量δ（＝（Ｒ０－Ｒｘ，Ｒ０－Ｒｙ））は、加工の高さ（座標値）ｚ、直角度誤差の動的補正パラメータＷｚｘ、Ｗｚｙ、及び切削点負荷Ｆで表されることから、数１４式のように表される。
なお、負荷電流と切削点負荷との関係パラメータ（比例係数）Ｊｘ、Ｊｙを求める必要があるのは、同じ負荷電流でも発生する切削点負荷ＦがＸ軸とＹ軸とで異なるためである。
関係パラメータ算出手段３２０は、テストワーク５０を切削加工したことにより発生する直角度誤差は工作機械２０の各軸間の剛性が低いことによるため、取得された負荷電流Ｉｆ（＝（Ｉｆ（Ｘ）、Ｉｆ（Ｙ））にベストフィットする楕円（（ｘ／Ｉｆｘ）２＋（ｙ／Ｉｆｙ）２＝１）を、例えば、最小二乗法等で算出する。なお、Ｉｆｘは楕円のＸ軸方向の半径である。また、Ｉｆｙは楕円のＹ軸方向の半径である。
ここで、数５式に基づいて、ズレ量δ（＝（Ｒ０－Ｒｘ，Ｒ０－Ｒｙ））は、加工の高さ（座標値）ｚ、直角度誤差の動的補正パラメータＷｚｘ、Ｗｚｙ、及び切削点負荷Ｆで表されることから、数１４式のように表される。
【数14】
【0076】
関係パラメータ算出手段３２０は、数１４式から変形した数１５式と、算出したＩｆｘ、Ｉｆｙと、動的補正パラメータデータ１２０から読み出した動的補正パラメータＷｚｘ、Ｗｚｙと、動的補正パラメータ算出手段２２０により算出されたＲ０、Ｒｘ、Ｒｙとを用い、負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙを算出することができる。
そして、関係パラメータ算出手段３２０は、算出した負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙを切削点負荷算出手段３３０に出力する。なお、関係パラメータ算出手段３２０は、算出した負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙを記憶部１００に記憶してもよい。
【数15】
【0077】
切削点負荷算出手段３３０は、実際の商品のワークを切削加工時において、関係パラメータ算出手段３２０により算出された加速度と負荷電流との関係パラメータＫｘ、Ｋｙ、及び速度と負荷電流との関係パラメータＬｘ、Ｌｙと、Ｘ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３から取得した実際の商品のワークの切削加工時の負荷電流Ｉと、数１６式とから切削による負荷電流Ｉｆを算出する。
なお、数１６式は、負荷電流Ｉが加速度による負荷電流、速度による負荷電流、及び切削による負荷電流Ｉｆの和となることから得られる。また、ａｘはＸ軸方向の加速度を示す。ａｙはＹ軸方向の加速度を示す。ｖｘはＸ軸方向の速度を示す。ｖｙはＹ軸方向の速度を示す。
そして、切削点負荷算出手段３３０は、関係パラメータ算出手段３２０により算出された負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙと、算出した負荷電流Ｉｆと、数１７式とを用いて、切削点負荷Ｆを算出する。
【数16】
そして、切削点負荷算出手段３３０は、関係パラメータ算出手段３２０により算出された負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙと、算出した負荷電流Ｉｆと、数１７式とを用いて、切削点負荷Ｆを算出する。
【数17】
【0078】
動的補正手段２５０は、動的補正パラメータデータ１２０から読み出した動的補正パラメータＷｚｘ、Ｗｚｙ、切削点負荷算出手段３３０により算出された切削点負荷Ｆ、及び数５式からズレ量δを算出する。動的補正手段２５０は、算出したズレ量δに基づいて実際の商品のワークに対する切削指令の指令座標値に対して直角度誤差の動的誤差を補正する。
【0079】
＜工具２５のたわみにおける負荷電流と切削点負荷との関係パラメータの算出について＞
工具２５のたわみには異方位性がないが、同じ切削点負荷を発生させるための負荷電流の大きさが、軸の慣性やモータ特性により異なる。そのため、切削点負荷を発生させるのに、大きな電流が必要な軸方向（慣性が大きな軸方向）では半径が大きく、小さな電流ですむ軸方向（慣性が小さな軸方向）では半径が小さくなる。
そこで、関係パラメータ算出手段３２０は、テストワーク５０を切削加工した際に取得した負荷電流Ｉｆ（＝（Ｉｆ（Ｘ）、Ｉｆ（Ｙ））にベストフィットする惰円（（ｘ／Ｉｆｔｘ）２＋（ｙ／Ｉｆｔｙ）２＝１）を、例えば、最小二乗法等で算出する。なお、Ｉｆｔｘは楕円のＸ軸方向の半径である。また、Ｉｆｔｙは楕円のＹ軸方向の半径である。
ここで、数７式に基づいて、ズレ量δ（＝（Ｒ０－Ｒｔ，Ｒ０－Ｒｔ））は、工具長Ｌｔ、動的補正パラメータＷｔ、及び切削点負荷Ｆで表されることから、数１８式のように表される。
関係パラメータ算出手段３２０は、数１８式から変形した数１９式と、算出したＩｆｔｘ、Ｉｆｔｙと、動的補正パラメータデータ１２０から読み出した動的補正パラメータＷｔと、動的補正パラメータ算出手段２２０により算出されたＲ０、Ｒｔとを用い、負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙを算出することができる。
そして、関係パラメータ算出手段３２０は、算出した負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙを切削点負荷算出手段３３０に出力する。なお、関係パラメータ算出手段３２０は、算出した負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙを記憶部１００に記憶してもよい。
工具２５のたわみには異方位性がないが、同じ切削点負荷を発生させるための負荷電流の大きさが、軸の慣性やモータ特性により異なる。そのため、切削点負荷を発生させるのに、大きな電流が必要な軸方向（慣性が大きな軸方向）では半径が大きく、小さな電流ですむ軸方向（慣性が小さな軸方向）では半径が小さくなる。
そこで、関係パラメータ算出手段３２０は、テストワーク５０を切削加工した際に取得した負荷電流Ｉｆ（＝（Ｉｆ（Ｘ）、Ｉｆ（Ｙ））にベストフィットする惰円（（ｘ／Ｉｆｔｘ）２＋（ｙ／Ｉｆｔｙ）２＝１）を、例えば、最小二乗法等で算出する。なお、Ｉｆｔｘは楕円のＸ軸方向の半径である。また、Ｉｆｔｙは楕円のＹ軸方向の半径である。
ここで、数７式に基づいて、ズレ量δ（＝（Ｒ０－Ｒｔ，Ｒ０－Ｒｔ））は、工具長Ｌｔ、動的補正パラメータＷｔ、及び切削点負荷Ｆで表されることから、数１８式のように表される。
【数18】
【数19】
【0080】
切削点負荷算出手段３３０は、実際の商品のワークを切削時において、関係パラメータ算出手段３２０により算出された加速度と負荷電流との関係パラメータＫｘ、Ｋｙ、及び速度と負荷電流との関係パラメータＬｘ、Ｌｙと、Ｘ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３から取得した実際の商品のワークの切削加工時の負荷電流Ｉと、数１６式とから切削による負荷電流Ｉｆを算出する。
そして、切削点負荷算出手段３３０は、関係パラメータ算出手段３２０により算出された負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙと、算出した負荷電流Ｉｆと、数１７式とを用いて、切削点負荷Ｆを算出する。
そして、切削点負荷算出手段３３０は、関係パラメータ算出手段３２０により算出された負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙと、算出した負荷電流Ｉｆと、数１７式とを用いて、切削点負荷Ｆを算出する。
【0081】
動的補正手段２５０は、動的補正パラメータデータ１２０から読み出した動的補正パラメータＷｔ、切削点負荷算出手段３３０により算出された切削点負荷Ｆ、及び数７式からズレ量δを算出する。動的補正手段２５０は、算出したズレ量δに基づいて実際の商品のワークに対する切削指令の指令座標値に対して工具２５のたわみによる誤差の動的誤差を補正する。
【0082】
以上により、第３実施形態の数値制御装置１０は、テストワーク５０の切削時の負荷電流を取得し、取得した負荷電流、指令形状と測定データとのズレ量δ、及び動的補正パラメータに基づいて負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙを算出する。数値制御装置１０は、負荷電流と切削点負荷との関係パラメータＪｘ、Ｊｙと実際の商品のワークの切削時の負荷電流Ｉとに基づいて切削点負荷Ｆを算出し、算出した切削点負荷Ｆ及び動的補正パラメータに基づいて切削指令が示す指令座標値に対して動的誤差を補正する。
これにより、数値制御装置１０は、センサを用いることなく、動的誤差を精度良く補正することができる。
また、数値制御装置１０は、ハードウェアに手を加えることなく、追加のハードウェアも不要で切削点負荷を推定することができる。
以上、第３実施形態について説明した。
これにより、数値制御装置１０は、センサを用いることなく、動的誤差を精度良く補正することができる。
また、数値制御装置１０は、ハードウェアに手を加えることなく、追加のハードウェアも不要で切削点負荷を推定することができる。
以上、第３実施形態について説明した。
【0083】
＜第３実施形態の変形例＞
上述の第３実施形態では、数値制御装置１０は、図１の数値制御装置１０に図１７の構成を追加することで、加工プログラムを解析することにより生成された切削指令の指令座標値に対して動的誤差を補正したがこれに限定されない。例えば、数値制御装置１０は、図１１の数値制御装置１０に図１７の構成を追加することで、工作機械２０のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３に出力するパルスに対して、動的誤差を補正する補正パルスを加算してもよい。
あるいは、数値制御装置１０は、図１２の数値制御装置１０に図１７の構成を追加することで、工具２５の工具径に動的補正の補正量を組み込み、動的誤差を補正してもよい。
上述の第３実施形態では、数値制御装置１０は、図１の数値制御装置１０に図１７の構成を追加することで、加工プログラムを解析することにより生成された切削指令の指令座標値に対して動的誤差を補正したがこれに限定されない。例えば、数値制御装置１０は、図１１の数値制御装置１０に図１７の構成を追加することで、工作機械２０のＸ軸サーボモータ３１、Ｙ軸サーボモータ３２、Ｚ軸サーボモータ３３に出力するパルスに対して、動的誤差を補正する補正パルスを加算してもよい。
あるいは、数値制御装置１０は、図１２の数値制御装置１０に図１７の構成を追加することで、工具２５の工具径に動的補正の補正量を組み込み、動的誤差を補正してもよい。
【0084】
以上、第１実施形態から第３実施形態について説明したが、数値制御装置１０は、上述の実施形態に限定されるものではなく、目的を達成できる範囲での変形、改良等を含む。
【0085】
＜変形例＞
上述の第１実施形態から第３実施形態では、工作機械２０は、直交３軸の工作機械としたが、５軸等の工作機械でもよい。
上述の第１実施形態から第３実施形態では、工作機械２０は、直交３軸の工作機械としたが、５軸等の工作機械でもよい。
【0086】
なお、第１実施形態から第３実施形態に係る数値制御装置１０に含まれる各機能は、ハードウェア、ソフトウェア又はこれらの組み合わせによりそれぞれ実現することができる。ここで、ソフトウェアによって実現されるとは、コンピュータがプログラムを読み込んで実行することにより実現されることを意味する。
また、数値制御装置１０に含まれる各構成部は、電子回路等を含むハードウェア、ソフトウェア又はこれらの組み合わせにより実現することができる。
また、数値制御装置１０に含まれる各構成部は、電子回路等を含むハードウェア、ソフトウェア又はこれらの組み合わせにより実現することができる。
【0087】
プログラムは、様々なタイプの非一時的なコンピュータ可読媒体（Ｎｏｎ－ｔｒａｎｓｉｔｏｒｙ ｃｏｍｐｕｔｅｒ ｒｅａｄａｂｌｅ ｍｅｄｉｕｍ）を用いて格納され、コンピュータに供給することができる。非一時的なコンピュータ可読媒体は、様々なタイプの実体のある記録媒体（Ｔａｎｇｉｂｌｅ ｓｔｏｒａｇｅ ｍｅｄｉｕｍ）を含む。非一時的なコンピュータ可読媒体の例は、磁気記録媒体（例えば、フレキシブルディスク、磁気テープ、ハードディスクドライブ）、光磁気記録媒体（例えば、光磁気ディスク）、ＣＤ－ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）、ＣＤ－Ｒ、ＣＤ－Ｒ／Ｗ、半導体メモリ（例えば、マスクＲＯＭ、ＰＲＯＭ（Ｐｒｏｇｒａｍｍａｂｌｅ ＲＯＭ）、ＥＰＲＯＭ（Ｅｒａｓａｂｌｅ ＰＲＯＭ）、フラッシュＲＯＭ、ＲＡＭを含む。また、プログラムは、様々なタイプの一時的なコンピュータ可読媒体（Ｔｒａｎｓｉｔｏｒｙ ｃｏｍｐｕｔｅｒ ｒｅａｄａｂｌｅ ｍｅｄｉｕｍ）によってコンピュータに供給されてもよい。一時的なコンピュータ可読媒体の例は、電気信号、光信号、及び電磁波を含む。一時的なコンピュータ可読媒体は、電線及び光ファイバ等の有線通信路、又は、無線通信路を介して、プログラムをコンピュータに供給できる。
【0088】
なお、記録媒体に記録されるプログラムを記述するステップは、その順序に沿って時系列的に行われる処理はもちろん、必ずしも時系列的に処理されなくとも、並列的あるいは個別に実行される処理をも含むものである。
【0089】
以上を換言すると、本開示の数値制御装置は、次のような構成を有する各種各様の実施形態を取ることができる。
【0090】
（１）本開示の数値制御装置１０は、指令解析手段２３０から受信した切削指令が示す指令座標値で工作機械２０に切削させる数値制御装置であって、切削されたテストワーク５０の形状を工作機械２０に機上測定させ、測定されたテストワーク５０の形状を示す測定データを取得する測定手段２１０と、切削指令が示す指令形状と測定手段２１０により取得された測定データとに基づいて、切削において工作機械２０に作用する力及び速度によって生じる動的誤差を補正する動的補正パラメータを算出する動的補正パラメータ算出手段２２０と、算出された前記動的補正パラメータに基づいて、指令座標値に対して動的誤差を補正する動的補正手段２５０と、を備え、動的補正パラメータ算出手段２２０は、指令形状と測定データとの比較から動的誤差のみを取得し、取得した動的誤差から動的補正パラメータを算出する。
この数値制御装置１０によれば、センサを用いることなく、動的誤差を精度良く補正することができる。
この数値制御装置１０によれば、センサを用いることなく、動的誤差を精度良く補正することができる。
【0091】
（２）指令形状、及びテストワーク５０の形状が円でもよい。
そうすることで、動的補正パラメータを容易に算出することができる。
そうすることで、動的補正パラメータを容易に算出することができる。
【0092】
（３）測定手段２１０は、ＸＹ平面でテストワーク５０に切削された半径Ｒ０の穴の測定データを取得し、動的補正パラメータ算出手段２２０は、測定データが示すテストワーク５０の形状に楕円をベストフィットして、前記楕円のＸ軸、Ｙ軸それぞれの半径Ｒｘ，Ｒｙを取得し、ＸＹ平面における動的補正パラメータＷｚｘ、Ｗｚｙを数２０式で算出してもよい。
ただし、Ｈｚは工作機械２０の定盤からＸ軸までの高さを示し、Ｈｗは定盤からテストワーク５０までの高さを示し、δ０及びδ９０は測定データが示す加工形状位置（Ｒｘｃｏｓ（θ＋α），Ｒｙｓｉｎ（θ＋α））の角度θが０度及び９０度のときのズレ量を示し、αはベストフィットした楕円が円周方向の負荷により傾いているときの傾き分ずらした位相を示す。
そうすることで、直角度誤差の動的補正パラメータを算出することができる。
【数20】
そうすることで、直角度誤差の動的補正パラメータを算出することができる。
【0093】
（４）測定手段２１０は、ＸＹ平面でテストワーク５０に切削された半径Ｒ０の穴の測定データを取得し、動的補正パラメータ算出手段２２０は、測定データが示すテストワーク５０の形状を円でベストフィットして、円の半径Ｒｔを取得し、切削負荷で生じる工作機械２０に含まれる工具のたわみ量の動的補正パラメータＷｔ、及び工具径補正量ｒｔ’を数２１式で算出してもよい。
ただし、βはワーク加工時のワーク法線方向と切削点負荷方向との間の角度を示し、ｒｔは補正前の工具径を示し、Ｌｔは工具長を示す。
そうすることで、工具のたわみによる誤差の動的補正パラメータを算出することができる。
【数21】
そうすることで、工具のたわみによる誤差の動的補正パラメータを算出することができる。
【0094】
（５）動的補正パラメータ補間手段３１０をさらに備え、測定手段２１０は、互いに異なる複数の切削負荷の各々で切削されたテストワーク５０の形状の測定データを取得し、動的補正パラメータ算出手段２２０は、複数の切削負荷の各々における指令形状と測定データとに基づいて動的補正パラメータを算出し、動的補正パラメータ補間手段３１０は、算出された複数の動的補正パラメータを用いて、任意の切削負荷における動的補正パラメータを補間し、動的補正手段２５０は、補間された任意の切削負荷における動的補正パラメータに基づいて切削指令の指令座標値を補正してもよい。
そうすることで、テストワーク５０の加工時の切削負荷の大きさと、実際の商品のワークの加工時の切削負荷の大きさが異なる場合でも、適切に動的誤差を補正することができる。
そうすることで、テストワーク５０の加工時の切削負荷の大きさと、実際の商品のワークの加工時の切削負荷の大きさが異なる場合でも、適切に動的誤差を補正することができる。
【0095】
（６）テストワークの切削時の負荷電流を取得し、取得した負荷電流、指令形状と測定データとのズレ量、及び動的補正パラメータに基づいて少なくともＸ軸及びＹ軸における負荷電流と切削点負荷との関係を示す関係パラメータを算出する関係パラメータ算出手段３２０と、関係パラメータと負荷電流とに基づいて切削点負荷を算出する切削点負荷算出手段３３０と、をさらに備えてもよい。
そうすることで、ハードウェアに手を加えることなく、追加のハードウェアも不要で切削点負荷を推定することができる。
そうすることで、ハードウェアに手を加えることなく、追加のハードウェアも不要で切削点負荷を推定することができる。
【0096】
（７）負荷電流は、工作機械２０に含まれるサーボモータの負荷電流から、工作機械２０に含まれる工具２５を加速させた時の加速度による負荷電流、及び工具２５を軸移動させた時の軸移動による負荷電流を引いて算出されてもよい。
そうすることで、切削による負荷電流Ｉｆを算出することができる。
そうすることで、切削による負荷電流Ｉｆを算出することができる。
【0097】
（８）加速度による負荷電流、及び軸移動による負荷電流は、切削を伴わない工具２５の移動時の負荷電流から算出されてもよい。
そうすることで、空加工時の加速度による負荷電流、及び軸移動による負荷電流を算出することができる。
そうすることで、空加工時の加速度による負荷電流、及び軸移動による負荷電流を算出することができる。
【0098】
（９）切削を伴わない工具２５の移動は、工具２５を円状に時計周りに移動させる移動、及び工具２５を円状に反時計周りに移動させる移動の両方を行ってもよい。
そうすることで、加速度による負荷電流、及び軸移動による負荷電流を算出することができる。
そうすることで、加速度による負荷電流、及び軸移動による負荷電流を算出することができる。
【0099】
（１０）本開示の制御方法は、コンピュータにより実行される、指令解析手段２３０から受信した切削指令が示す指令座標値で工作機械に切削させる制御方法であって、切削されたテストワーク５０の形状を工作機械２０に機上測定させ、測定されたテストワーク５０の形状を示す測定データを取得する測定ステップと、切削指令が示す指令形状と取得された測定データとに基づいて、切削において工作機械２０に作用する力及び速度によって生じる動的誤差を補正する動的補正パラメータを算出する動的補正パラメータ算出ステップと、算出された動的補正パラメータに基づいて、指令座標値に対して動的誤差を補正する動的補正ステップと、を備え、動的補正パラメータ算出ステップは、指令形状と測定データとの比較から動的誤差のみを取得し、取得した動的誤差から動的補正パラメータを算出する。
この制御方法によれは、（１）と同様の効果を奏することができる。
この制御方法によれは、（１）と同様の効果を奏することができる。
【符号の説明】
【0100】
１０ 数値制御装置
２０ 工作機械
２５ 工具
５０ テストワーク
１００ 記憶部
２００ 制御部
２１０ 測定手段
２２０ 動的補正パラメータ算出手段
２３０ 指令解析手段
２５０ 動的補正手段
３１０ 動的補正パラメータ補間手段
３２０ 関係パラメータ算出手段
３３０ 切削点負荷算出手段
２０ 工作機械
２５ 工具
５０ テストワーク
１００ 記憶部
２００ 制御部
２１０ 測定手段
２２０ 動的補正パラメータ算出手段
２３０ 指令解析手段
２５０ 動的補正手段
３１０ 動的補正パラメータ補間手段
３２０ 関係パラメータ算出手段
３３０ 切削点負荷算出手段
【図1】
【図2】
【図3A】
【図3B】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15A】
【図15B】
【図16】
【図17】
【図18A】
【図18B】