JP5927904B2 - Machine tool dynamic characteristics calculation device - Google Patents

Machine tool dynamic characteristics calculation device Download PDF

Info

Publication number
JP5927904B2
JP5927904B2 JP2011283254A JP2011283254A JP5927904B2 JP 5927904 B2 JP5927904 B2 JP 5927904B2 JP 2011283254 A JP2011283254 A JP 2011283254A JP 2011283254 A JP2011283254 A JP 2011283254A JP 5927904 B2 JP5927904 B2 JP 5927904B2
Authority
JP
Japan
Prior art keywords
rotary tool
tool
cutting
displacement amount
dynamic characteristic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2011283254A
Other languages
Japanese (ja)
Other versions
JP2013132706A (en
Inventor
良彦 山田
良彦 山田
賢治 濱田
賢治 濱田
孝幸 東
孝幸 東
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JTEKT Corp
Original Assignee
JTEKT Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by JTEKT Corp filed Critical JTEKT Corp
Priority to JP2011283254A priority Critical patent/JP5927904B2/en
Publication of JP2013132706A publication Critical patent/JP2013132706A/en
Application granted granted Critical
Publication of JP5927904B2 publication Critical patent/JP5927904B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Machine Tool Sensing Apparatuses (AREA)

Description

本発明は、切削加工を行う際に用いられるエンドミル等の回転工具の動特性を算出する技術に関するものである。   The present invention relates to a technique for calculating dynamic characteristics of a rotary tool such as an end mill used when cutting.

エンドミル等の回転工具により切削加工を行う場合、回転工具の動特性を取得することは、適正な加工条件を決定する上で重要である。特許文献1には、主軸にアンバランスマスタを装着して、当該アンバランスマスタの振れ量を振れ量測定器によって検出し、主軸の動剛性を算出する技術が提案されている。しかしながら、特許文献1に示される技術では、主軸の動剛性を算出することはできるが、回転工具の動特性を取得することはできなかった。このため、従来では、作業者が回転工具をハンマーで叩いて、回転工具を加振させ、回転工具の共振周波数を測定し(ハンマリング試験)、この共振周波数から動特性を算出するのが一般的であった。   When cutting with a rotary tool such as an end mill, acquiring the dynamic characteristics of the rotary tool is important in determining appropriate machining conditions. Patent Document 1 proposes a technique in which an unbalance master is attached to a main shaft, a shake amount of the unbalance master is detected by a shake amount measuring device, and dynamic stiffness of the main shaft is calculated. However, with the technique disclosed in Patent Document 1, it is possible to calculate the dynamic rigidity of the spindle, but it is not possible to acquire the dynamic characteristics of the rotary tool. For this reason, conventionally, an operator generally hits the rotary tool with a hammer, vibrates the rotary tool, measures the resonance frequency of the rotary tool (hammering test), and calculates the dynamic characteristics from this resonance frequency. It was the target.

特開2010−274375号公報(図1)JP 2010-274375 A (FIG. 1)

しかしながら、このような方法で、動特性を算出するには、回転工具に加速度計を取り付けたうえで、回転工具をハンマーで叩く必要があることから、動特性の計測に時間がかかるという問題があった。また、作業者がハンマーで回転工具叩いた後に、復元する回転工具がハンマーに接触してしまうことがあり、この場合には測定ミスとなることから、ハンマーの叩き方に経験を要するという問題があった。更に、主軸回転中においては、ハンマリング試験により動特性を計測することができず、主軸が静止時と回転中とでは、軸受の剛性が異なることから、正確な動特性を算出することができないという問題があった。   However, in order to calculate the dynamic characteristics by such a method, it is necessary to attach an accelerometer to the rotary tool and then hit the rotary tool with a hammer. there were. In addition, after the operator hits the rotary tool with a hammer, the restored rotary tool may come into contact with the hammer. In this case, a measurement error occurs. there were. In addition, the dynamic characteristics cannot be measured by a hammering test while the main spindle is rotating, and the dynamic characteristics cannot be calculated accurately because the rigidity of the bearing differs when the main spindle is stationary and rotating. There was a problem.

本発明は、このような事情に鑑みてなされたものであり、手間をかけず正確な回転工具の動特性を取得することができる工作機械の動特性算出装置を提供することを目的とする。   The present invention has been made in view of such circumstances, and an object of the present invention is to provide a machine tool dynamic characteristic calculation apparatus that can acquire accurate dynamic characteristics of a rotary tool without taking time and effort.

(請求項1)
外周側に周方向に1以上の刃部を備える回転工具を用いて、当該回転工具を軸回りに回転しながら被加工物に対して相対移動して断続的な切削加工を行う工作機械の動特性を算出する装置であって、断続的な切削加工により前記回転工具が加振された場合に、振動している状態の前記回転工具の変位量を検出する工具変位量検出手段と、前記回転工具の前記切削抵抗を算出する切削抵抗算出手段と、前記工具変位量検出手段が検出した振動している状態の前記回転工具の前記変位量、及び、前記切削抵抗算出手段によって算出された前記回転工具の前記切削抵抗に基づいて、前記回転工具の動特性である質量係数、弾性係数、及び粘性減衰係数の少なくとも1つを算出する動特性算出手段と、を備える。
(Claim 1)
Using a rotary tool having one or more blades in the circumferential direction on the outer peripheral side, the machine tool moves intermittently by moving relative to the workpiece while rotating the rotary tool about its axis. A device for calculating characteristics, wherein when the rotary tool is vibrated by intermittent cutting, a tool displacement amount detecting means for detecting a displacement amount of the rotary tool in a vibrating state, and the rotation and cutting resistance calculating means for calculating the cutting resistance of the tool, the displacement amount of the rotary tool in a state in which the tool displacement amount detecting means is vibrating the detected and the rotation that has been calculated by the cutting resistance calculating means based on the cutting resistance of the tool, it comprises a mass coefficient is a dynamic characteristic of the rotating tool, elastic modulus, and the dynamic characteristic calculating means for calculating at least one of the viscous damping coefficient, a.

(請求項2)
前記被加工物の試し加工中において、前記工具変位量検出手段によって検出された振動している状態の前記回転工具の前記変位量に基づいて、前記動特性算出手段は、前記回転工具の前記動特性を算出する。
(Claim 2)
During trial machining of the workpiece, on the basis of the displacement amount of the rotary tool in the state that vibrations detected by said tool displacement amount detecting means, the dynamic characteristic calculating means, the movement of the rotary tool Calculate the characteristics.

(請求項3)
回転工具変位量検出手段は、撮像装置であってもよい。
(Claim 3)
The rotary tool displacement amount detection means may be an imaging device.

(請求項1)
動特性算出手段は、回転工具の変位量、及び、回転工具の切削抵抗に基づいて、回転工具の動特性である質量係数、弾性係数、及び粘性減衰係数の少なくとも1つを算出する。これにより、回転工具が取り付けられる主軸の回転中における回転工具の動特性が算出される。このため、主軸が停止している状態で動特性を算出する場合と比較して、正確に動特性を取得することが可能となる。また、工具変位量検出手段は、断続的な切削加工により加振された回転工具の変位量を検出し、当該検出された回転工具の変位量に基づいて、動特性を算出する。このため、従来のように、回転工具に加速度計を取り付ける必要もなく、ハンマーの叩き方にも経験を要さないことから、手間をかけずに動特性を取得することができる。

(Claim 1)
The dynamic characteristic calculation means calculates at least one of a mass coefficient, an elastic coefficient, and a viscous damping coefficient, which are dynamic characteristics of the rotary tool, based on the displacement amount of the rotary tool and the cutting resistance of the rotary tool. Thereby, the dynamic characteristic of the rotary tool during rotation of the spindle to which the rotary tool is attached is calculated. For this reason, it is possible to obtain the dynamic characteristics more accurately than in the case where the dynamic characteristics are calculated while the main shaft is stopped. Further, the tool displacement amount detection means detects the displacement amount of the rotating tool that is vibrated by intermittent cutting, and calculates dynamic characteristics based on the detected displacement amount of the rotating tool. For this reason, it is not necessary to attach an accelerometer to the rotary tool as in the prior art, and since no experience is required in how to strike the hammer, the dynamic characteristics can be acquired without taking time and effort.

(請求項2)
被加工物の試し加工中において、前記工具変位量検出手段によって検出された振動している状態の前記回転工具の変位量に基づいて、前記動特性算出手段は、回転工具の動特性を算出する。これにより、最終製品の品質に影響すること無く、回転工具の動特性を算出することができる。
(Claim 2)
During the trial machining of the workpiece, the dynamic characteristic calculation unit calculates the dynamic characteristic of the rotary tool based on the displacement amount of the rotating tool in a vibrating state detected by the tool displacement amount detection unit. . As a result, the dynamic characteristics of the rotary tool can be calculated without affecting the quality of the final product.

(請求項3)
撮像装置によって、切削加工中の回転工具の変位量を検出する。これにより、撮像装置によって撮像された回転工具を含む画像データから、回転工具の外縁を検出し、回転工具の変位を検出する。このため、従来のように、加速度センサ等によって回転工具の変位を検出する方法と比較して、回転工具が取り付けられる回転主軸の回転に伴う振動等のノイズの影響を受けず、より高精度に回転工具の変位を検出することができ、微少な回転工具の変位も検出することができる。また、加速度センサを取り付けることができない小径の回転工具であっても、動特性を取得することができる。
(Claim 3)
The amount of displacement of the rotary tool during cutting is detected by the imaging device. Thereby, the outer edge of the rotary tool is detected from the image data including the rotary tool imaged by the imaging device, and the displacement of the rotary tool is detected. For this reason, compared with the conventional method of detecting the displacement of the rotary tool by an acceleration sensor or the like, it is not affected by noise such as vibration accompanying rotation of the rotary spindle to which the rotary tool is attached, and more accurately. The displacement of the rotary tool can be detected, and a slight displacement of the rotary tool can also be detected. Moreover, even if it is a small diameter rotary tool which cannot attach an acceleration sensor, a dynamic characteristic can be acquired.

本実施形態における工作機械の構成を示す図である。It is a figure which shows the structure of the machine tool in this embodiment. ワーク、回転工具、及び、高速度カメラの位置関係を表した上面図である。It is a top view showing the positional relationship of a workpiece | work, a rotary tool, and a high-speed camera. 回転工具に生じる切削抵抗を示す説明図である。It is explanatory drawing which shows the cutting resistance which arises in a rotary tool. 加工誤差の発生メカニズムを説明するための図であって、回転工具に生じる切削抵抗、実切込量、及び、回転工具の回転中心の変位の経過時間に対する挙動を示す。It is a figure for demonstrating the generation | occurrence | production mechanism of a process error, Comprising: The cutting resistance which arises in a rotary tool, the actual cutting amount, and the behavior with respect to the elapsed time of the displacement of the rotation center of a rotary tool are shown. 図4の各時刻における回転工具と被加工物との位置関係を示す。The positional relationship between the rotary tool and the workpiece at each time in FIG. 4 is shown. 本実施形態の工作機械の機能ブロック図である。It is a functional block diagram of the machine tool of this embodiment. 回転工具の反切込方向Yの変位Tyの経時変化を表したグラフである。It is a graph showing the time-dependent change of the displacement Ty of the anti-cutting direction Y of a rotary tool. 図6の動特性算出部による動特性算出処理を示すフローチャートである。It is a flowchart which shows the dynamic characteristic calculation process by the dynamic characteristic calculation part of FIG. 回転工具の変移、変形速度、及び、切削速度の経時変化を表したグラフである。It is a graph showing the time-dependent change of the change of a rotary tool, a deformation speed, and a cutting speed. 図1の制御部で実行される動特性算出・加工処理を示すフローチャートである。It is a flowchart which shows the dynamic characteristic calculation and the process performed by the control part of FIG.

(1.対象工作機械の構成)
本実施形態の動特性算出装置を適用する工作機械100の構成について説明する。工作機械100は、被加工物Wを回転工具により切削加工する工作機械である。その工作機械100の一例として、図1を参照して説明する。図1に示すように、工作機械100は、ベッド1と、回転工具5を保持する回転主軸4と、ベッド1上にて被加工物Wを載置するテーブル6、回転工具5を撮像する高速度カメラ9、工作機械100の統括制御を行う制御部50とを備える。
(1. Configuration of the target machine tool)
A configuration of the machine tool 100 to which the dynamic characteristic calculation apparatus of the present embodiment is applied will be described. The machine tool 100 is a machine tool that cuts the workpiece W with a rotary tool. An example of the machine tool 100 will be described with reference to FIG. As shown in FIG. 1, the machine tool 100 captures an image of a bed 1, a rotating spindle 4 that holds a rotating tool 5, a table 6 on which a workpiece W is placed on the bed 1, and a rotating tool 5. A speed camera 9 and a control unit 50 that performs overall control of the machine tool 100 are provided.

ここで、回転工具5は、外周側に周方向に1以上の刃部5a,5bを備える。回転工具5は、例えば、ボールエンドミル、スクエアエンドミル、フライスなどを含む。回転主軸4は、図示しないアクチュエータによって、X軸、Y軸、Z軸方向に移動可能となっている。工作機械100は、回転工具5を軸回りに回転しながら、被加工物Wに対して相対移動することにより、断続的な切削加工を行う。   Here, the rotary tool 5 includes one or more blade portions 5a and 5b in the circumferential direction on the outer peripheral side. The rotary tool 5 includes, for example, a ball end mill, a square end mill, and a milling cutter. The rotation main shaft 4 can be moved in the X-axis, Y-axis, and Z-axis directions by an actuator (not shown). The machine tool 100 performs intermittent cutting by moving relative to the workpiece W while rotating the rotary tool 5 about the axis.

制御部50は、各種プログラムが実行されるCPU50a、及び、CPU50aで実行されるプログラムや各種設定値が記憶される記憶部50bを有している。記憶部50bに記憶されCPU50aで実行されるプログラムは、後述の図6に示す各算出部、処理部、制御部に相当する。記憶部50bは、後述の図6に示す各記憶部に相当する。なお、後述の図6に示す各算出部、処理部を、ASIC(Application Specific Integrated Circuit)で構成しても差し支え無い。   The control unit 50 includes a CPU 50a that executes various programs and a storage unit 50b that stores programs executed by the CPU 50a and various setting values. The program stored in the storage unit 50b and executed by the CPU 50a corresponds to each calculation unit, processing unit, and control unit shown in FIG. The storage unit 50b corresponds to each storage unit shown in FIG. Each calculation unit and processing unit shown in FIG. 6 to be described later may be configured by ASIC (Application Specific Integrated Circuit).

高速度カメラ9は、CCDやCMOS等の撮像素子と、この撮像素子に被写像を結像するレンズを有していて、回転工具5の側方に配置されている。高速度カメラ9は、回転工具5を側方から連続的に(例えば5000フレーム/秒)撮像する。高速度カメラ9で撮像された回転工具5を含む画像は、制御部50に出力される。図2に示した実施形態では高速度カメラ9は、回転工具5に対して送り方向Xに配設され、回転工具5の反切込方向Yの変形が撮像できるようになっている。なお、制御部50、及び、高速度カメラ9が、回転工具5の動特性を算出する動特性算出装置に相当する。図示しないが、工作機械100は、クーラントを供給するクーラントノズル、クーラントポンプなどを備える。   The high-speed camera 9 has an image sensor such as a CCD or CMOS, and a lens that forms an image on the image sensor, and is disposed on the side of the rotary tool 5. The high-speed camera 9 images the rotary tool 5 continuously from the side (for example, 5000 frames / second). An image including the rotary tool 5 captured by the high-speed camera 9 is output to the control unit 50. In the embodiment shown in FIG. 2, the high-speed camera 9 is arranged in the feed direction X with respect to the rotary tool 5 so that the deformation of the rotary tool 5 in the anti-cutting direction Y can be imaged. The control unit 50 and the high-speed camera 9 correspond to a dynamic characteristic calculation device that calculates the dynamic characteristic of the rotary tool 5. Although not shown, the machine tool 100 includes a coolant nozzle that supplies coolant, a coolant pump, and the like.

(2.加工中の回転工具に作用する切削抵抗及び変位の説明)
以下に、図3を参照して、加工中の回転工具5に作用する切削抵抗及び変位を説明する。被加工物Wを切削中の回転工具5には、反切込方向Yの切削抵抗Fy及び反送り方向の切削抵抗Fxが生じる。つまり、回転工具5の先端側の回転中心Cは、各方向の切削抵抗Fx,Fyの合成抵抗Fxyの方向に変位Tする(図2示)。そして、切削抵抗Fyによって、回転工具5の先端側(刃部5a,5bの部位)の回転中心Cが反切込方向Yの変位Tyすることにより、被加工物Wの加工後形状が変化する。なお、図3には、スクエアエンドミルを図示しているが、ボールエンドミルの場合も同様である。
(2. Explanation of cutting force and displacement acting on rotating tool during machining)
Below, with reference to FIG. 3, the cutting force and displacement which act on the rotary tool 5 during a process are demonstrated. The rotary tool 5 that is cutting the workpiece W has a cutting resistance Fy in the anti-cutting direction Y and a cutting resistance Fx in the counterfeed direction. That is, the rotation center C on the tip side of the rotary tool 5 is displaced T in the direction of the combined resistance Fxy of the cutting resistances Fx and Fy in each direction (shown in FIG. 2). And the post-processing shape of the workpiece W changes by the displacement Ty of the anti-cutting direction Y by the cutting center F of the rotation center C of the front end side (part of blade part 5a, 5b) of the rotary tool 5 by cutting resistance Fy. Although FIG. 3 shows a square end mill, the same applies to a ball end mill.

ここで、回転工具5に生じる切削抵抗Fyが一定であれば、回転工具5の先端側の撓み量は一定となる。しかし、回転工具5による断続的な切削加工においては、回転工具5に生じる切削抵抗Fyは逐次変化する。そのため、回転工具5の先端側の回転中心Cの変位量は、主としてY方向に逐次変化するこのときの回転工具5の先端側の回転中心Cの変位量と切削抵抗Fyとは、回転工具5の動特性に依存する。工具の動特性とは、入力された力に対する変形の挙動を示すものであり、伝達関数(コンプライアンスおよび位相遅れ)もしくはそれから算出される質量係数M、粘性減衰係数C、弾性係数(バネ乗数)K、固有振動数fd、減衰率ζなどにより表される。図2に示すように、本実施形態では、上面視した場合に、回転工具5の反切込方向Yの変位Tyが、高速度カメラ9によって撮像されて検出される。   Here, if the cutting resistance Fy generated in the rotary tool 5 is constant, the amount of deflection on the tip side of the rotary tool 5 is constant. However, in the intermittent cutting with the rotary tool 5, the cutting resistance Fy generated in the rotary tool 5 changes sequentially. Therefore, the displacement amount of the rotation center C on the distal end side of the rotary tool 5 mainly changes sequentially in the Y direction. At this time, the displacement amount of the rotation center C on the distal end side of the rotary tool 5 and the cutting resistance Fy are the rotation tool 5. Depends on the dynamic characteristics of The dynamic characteristics of the tool indicate the deformation behavior with respect to the input force. The transfer function (compliance and phase lag) or the mass coefficient M, the viscous damping coefficient C, and the elastic coefficient (spring multiplier) K calculated from the transfer function (compliance and phase lag). And the natural frequency fd, the damping rate ζ, and the like. As shown in FIG. 2, in this embodiment, when viewed from the top, the displacement Ty in the anti-cutting direction Y of the rotary tool 5 is imaged and detected by the high-speed camera 9.

次に、回転工具5を回転しかつ送りながら被加工物Wの断続的な切削加工を行う際において、回転工具5に作用する切削抵抗Fy及び回転工具5の先端側の回転中心Cの変位量Tyの経過時間tに対する挙動について、図4及び図5を参照して説明する。ここでは、反切込方向(Y方向)における切削抵抗Fy及び先端側の回転中心Cの変位量Tyについて取り上げて説明する。   Next, when performing intermittent cutting of the workpiece W while rotating and feeding the rotary tool 5, the cutting resistance Fy acting on the rotary tool 5 and the displacement amount of the rotation center C on the tip side of the rotary tool 5. The behavior of Ty with respect to the elapsed time t will be described with reference to FIGS. Here, the cutting resistance Fy in the anti-cutting direction (Y direction) and the displacement amount Ty of the rotation center C on the tip side will be described.

図4の実線と図5(a)〜(e)を参照して、切削抵抗Fyの経過時間tに対する挙動について説明する。図4の実線に示すように、切削抵抗Fyは、ゼロ付近から時刻t1にて大きな値に変化し、時刻t2に再びゼロ付近に変化している。図5(a)(b)が、それぞれ図4の時刻t1,t2に対応する。図5(a)に示すように、時刻t1は、一方の刃部5aが被加工物Wに接触開始した瞬間である。つまり、時刻t1は、一方の刃部5aにより切削加工を開始した瞬間である。一方、図5(b)に示すように、時刻t2は、一方の刃部5aによる被加工物Wの切削加工を終了した瞬間である。このように、t1〜t2の間において、一方の刃部5aが切削加工している。   The behavior of the cutting resistance Fy with respect to the elapsed time t will be described with reference to the solid line in FIG. 4 and FIGS. As shown by the solid line in FIG. 4, the cutting resistance Fy changes from near zero to a large value at time t1, and again changes to near zero at time t2. FIGS. 5A and 5B correspond to times t1 and t2 in FIG. 4, respectively. As shown in FIG. 5A, time t1 is the moment when one of the blades 5a starts to contact the workpiece W. That is, time t1 is the moment when cutting is started by one of the blade portions 5a. On the other hand, as shown in FIG. 5B, the time t2 is the moment when the cutting of the workpiece W by the one blade portion 5a is finished. Thus, one blade part 5a is cutting between t1 and t2.

その後、図4に示すように、t2〜t4の間は、切削抵抗Fyがゼロとなっている。この間は、時刻t3に対応する図5(c)に示すように、両方の刃部5a,5bが被加工物Wに接触していない。つまり、回転工具5は空転している。   Thereafter, as shown in FIG. 4, the cutting resistance Fy is zero between t2 and t4. During this time, as shown in FIG. 5C corresponding to time t3, both the blade portions 5a and 5b are not in contact with the workpiece W. That is, the rotary tool 5 is idling.

その後、図4に示すように、切削抵抗Fyが時刻t4に再び大きな値に変化し、時刻t5に再びゼロ付近に変化している。図4の時刻t4には、対応する図5(d)に示すように、他方の刃部5bが被加工物Wに接触開始している。つまり、他方の刃部5bにより切削加工を開始している。また、図4の時刻t5には、対応する図5(e)に示すように、他方の刃部5bによる切削加工を終了している。このように、t4〜t5の間において、他方の刃部5bが切削加工している。   After that, as shown in FIG. 4, the cutting resistance Fy changes again to a large value at time t4 and again changes to near zero at time t5. At time t4 in FIG. 4, the other blade portion 5b starts to contact the workpiece W as shown in the corresponding FIG. 5 (d). That is, cutting is started by the other blade portion 5b. Further, at time t5 in FIG. 4, as shown in the corresponding FIG. 5 (e), the cutting with the other blade portion 5b is finished. Thus, the other blade part 5b is cutting between t4 and t5.

ここで、図5における今回の切削領域より、t1〜t2、t4〜t5の各瞬間において、実切込量h(瞬間的な切込量)が異なることが分かる。つまり、実切込量hは、切削開始から一気に多くなり、ピークに達した後に徐々に少なくなっている。より詳細には、前回切削されていない部位と前回切削された部位との境界の前後で変化している。そして、図4の切削抵抗Fyのうち急激に大きくなっている部分に示すように、切削加工中の切削抵抗Fyは、略三角形状になっており、実切込量hに応じて変化していることが分かる。   Here, it can be seen from the current cutting region in FIG. 5 that the actual cutting amount h (instantaneous cutting amount) differs at the instants t1 to t2 and t4 to t5. That is, the actual cutting amount h increases at a stroke from the start of cutting and gradually decreases after reaching the peak. In more detail, it changes before and after the boundary between the part not cut last time and the part cut last time. Then, as shown in the portion of the cutting force Fy of FIG. 4 that is abruptly increased, the cutting resistance Fy during the cutting process has a substantially triangular shape, and changes according to the actual cutting depth h. I understand that

また、回転工具5は、時刻t1,t4において切削加工を開始するということは、換言すると、時刻t1,t4において被加工物Wに衝突するということになる。つまり、回転工具5が空転状態から切削加工を開始する瞬間には、回転工具5には、被加工物Wとの衝突による断続的な切削抵抗Fyが発生する。   In addition, when the rotary tool 5 starts cutting at times t1 and t4, in other words, it collides with the workpiece W at times t1 and t4. That is, at the moment when the rotary tool 5 starts cutting from the idling state, an intermittent cutting resistance Fy is generated in the rotary tool 5 due to a collision with the workpiece W.

つまり、回転工具5の先端側の回転中心Cは、切削加工している間の切削抵抗Fyの変動によって、少なくとも反切込方向(Y方向)への加速度を生じる。さらに、断続切削であることによって、回転工具5の先端側の回転中心Cは、切削加工している間の切削抵抗Fy(衝撃力のような力)に起因して、少なくとも反切込方向(Y方向)に振動する(加振される)。   That is, the rotation center C on the distal end side of the rotary tool 5 generates an acceleration in at least the anti-cutting direction (Y direction) due to fluctuations in the cutting resistance Fy during cutting. Furthermore, due to the intermittent cutting, the rotational center C on the tip side of the rotary tool 5 is at least in the anti-cutting direction (Y) due to the cutting resistance Fy (force such as impact force) during cutting. Direction).

従って、回転工具5の先端側の回転中心Cの変位量Tyは、図4の破線に示すように、回転工具5の固有値に応じて振動している。特に、切削加工中に変動する切削抵抗Fyが加振源となり振動を起こし、回転工具5が空転している間に、減衰している。そして、再び、切削抵抗Fyにより変位量Tyが大きくなり、繰り返す。   Therefore, the displacement amount Ty of the rotation center C on the tip side of the rotary tool 5 vibrates according to the eigenvalue of the rotary tool 5 as shown by the broken line in FIG. In particular, the cutting resistance Fy, which fluctuates during the cutting process, becomes a vibration source and causes vibration, and is damped while the rotary tool 5 is idling. Then, the displacement amount Ty is increased again by the cutting resistance Fy, and the process is repeated.

(3.工作機械の機能構成)
次に、工作機械100の機能構成についての詳細を図6〜図10を参照して説明する。工作機械100は、図6の機能ブロック図に示すように構成される。以下に、図6に示す工作機械100の機能構成について説明する。
(3. Functional configuration of machine tools)
Next, details of the functional configuration of the machine tool 100 will be described with reference to FIGS. The machine tool 100 is configured as shown in the functional block diagram of FIG. The functional configuration of the machine tool 100 shown in FIG. 6 will be described below.

切削抵抗検出センサ71は、切削加工中において回転工具5に作用する反切込方向Yの切削抵抗Fyを検出し、検出信号を切削抵抗算出部32に出力する。例えば、切削抵抗検出センサ71は、荷重センサ、変位センサ、送り軸の駆動モータの消費電力検出器、供給電流センサなどを適用できる。つまり、荷重センサにより切削抵抗Fyそのものを直接検出することもできるし、変位センサやその他により間接的に切削抵抗Fyを検出することもできる。   The cutting resistance detection sensor 71 detects the cutting resistance Fy in the anti-cutting direction Y that acts on the rotary tool 5 during cutting, and outputs a detection signal to the cutting resistance calculation unit 32. For example, as the cutting resistance detection sensor 71, a load sensor, a displacement sensor, a power consumption detector of a drive motor for a feed shaft, a supply current sensor, or the like can be applied. That is, the cutting force Fy itself can be directly detected by a load sensor, or the cutting force Fy can be indirectly detected by a displacement sensor or the like.

切削抵抗算出部32は、切削抵抗検出センサ71から出力された検出信号に基づいて、切削加工中において回転工具5に作用する反切込方向Yの切削抵抗Fyを算出する。なお、切削抵抗算出部32は、回転工具5の情報や加工条件から、シミュレーションにより、切削加工中において回転工具5に作用する反切込方向Yの切削抵抗Fyを算出しても差し支え無い。   The cutting resistance calculation unit 32 calculates the cutting resistance Fy in the anti-cutting direction Y that acts on the rotary tool 5 during cutting based on the detection signal output from the cutting resistance detection sensor 71. Note that the cutting resistance calculation unit 32 may calculate the cutting resistance Fy in the anti-cutting direction Y acting on the rotary tool 5 during cutting from the information and processing conditions of the rotary tool 5 by simulation.

工具変位量算出部35は、高速度カメラ9で撮像された回転工具5の回転中心Cの反切込方向Yの変位Tyを算出する。具体的には、工具変位量算出部35は、撮像された回転工具5を含む画像データから、その明暗差等により回転工具5の画像とそれ以外(例えば背景)の部分の画像の境界(エッジ)を検出することにより、回転工具5の外縁を検出する。そして、工具変位量算出部35は、1フレーム毎に、回転工具5の軸方向所定位置において、回転工具5の外縁が反切込方向Yに何ピクセル移動したかを検出することにより、回転工具5の回転中心Cの反切込方向Yの変位Tyを算出する。このようにして、断続的な切削加工により加振され振動している状態の回転工具5の変位量(図7示)が連続的に算出される。なお、刃部5a,5bの有無によって、回転工具5の反切込方向Yの変位Tyが影響しないように、回転工具5の外縁の検出は、回転工具5の刃部5a,5bが形成されていない基部の外周面において行われる。   The tool displacement amount calculation unit 35 calculates the displacement Ty in the anti-cutting direction Y of the rotation center C of the rotary tool 5 captured by the high speed camera 9. Specifically, the tool displacement amount calculation unit 35 determines the boundary (edge) between the image of the rotary tool 5 and the image of the other portion (for example, the background) from the image data including the captured rotary tool 5 based on the brightness difference or the like. ) To detect the outer edge of the rotary tool 5. Then, the tool displacement amount calculation unit 35 detects how many pixels the outer edge of the rotary tool 5 has moved in the anti-cutting direction Y at a predetermined position in the axial direction of the rotary tool 5 for each frame. The displacement Ty of the rotation center C in the anti-cutting direction Y is calculated. In this way, the displacement amount (shown in FIG. 7) of the rotary tool 5 in the state of being vibrated and vibrated by intermittent cutting is continuously calculated. Note that the outer edge of the rotary tool 5 is detected by forming the blade parts 5a and 5b of the rotary tool 5 so that the displacement Ty in the anti-cutting direction Y of the rotary tool 5 is not affected by the presence or absence of the blade parts 5a and 5b. Not on the outer peripheral surface of the base.

工具変位量記憶部36には、図7に示すように、断続的な切削加工により加振され振動している状態の回転工具5の反切込方向Yの変位Tyの経時変化が連続的に記憶される。図7に示すように、回転工具5の反切込方向Yの変位Tyは、回転工具15の刃部5a、刃部5bのいずれかが最初に被加工物Wに当接した後に、最大の変位Ty1となり、その後、徐々に減衰しながら一定の周期Tで振動する。   As shown in FIG. 7, the tool displacement amount storage unit 36 continuously stores temporal changes in the displacement Ty in the anti-cutting direction Y of the rotary tool 5 that is vibrated and vibrated by intermittent cutting. Is done. As shown in FIG. 7, the displacement Ty in the anti-cutting direction Y of the rotary tool 5 is the maximum displacement after any one of the blade part 5a and the blade part 5b of the rotary tool 15 first comes into contact with the workpiece W. It becomes Ty1, and then vibrates at a constant period T while gradually damped.

動特性算出部37は、工具変位量記憶部36に記憶されている回転工具5の反切込方向Yの連続的な変位Ty、及び、切削抵抗算出部32で算出された切削抵抗Fyに基づいて、回転工具5の動特性を算出する。動特性算出部37による回転工具5の動特性(固有振動数fd、角固有振動数ωd、減衰率ζ、質量係数M、粘性減衰係数C、弾性係数K)の算出について、図8に示すフローチャートを参照して以下に説明する。   The dynamic characteristic calculation unit 37 is based on the continuous displacement Ty in the anti-cutting direction Y of the rotary tool 5 stored in the tool displacement amount storage unit 36 and the cutting resistance Fy calculated by the cutting resistance calculation unit 32. The dynamic characteristics of the rotary tool 5 are calculated. FIG. 8 is a flowchart illustrating calculation of dynamic characteristics (natural frequency fd, angular natural frequency ωd, damping factor ζ, mass coefficient M, viscous damping coefficient C, and elastic coefficient K) of the rotary tool 5 by the dynamic characteristic calculation unit 37. Will be described below.

まず、動特性算出部37は、工具変位量記憶部36の記憶されている回転工具5の反切込方向Yの連続的な変位Tyに基づき、変位Tyの角固有振動数ωdを算出する(S21)。具体的には、まず、式(1)に示すように、動特性算出部37は、回転工具15の最大変位Ty1、Ty2の時刻t1、t2の間隔を検出して、変位Tyの周期T[s](図7示)を検出する。   First, the dynamic characteristic calculation unit 37 calculates the angular natural frequency ωd of the displacement Ty based on the continuous displacement Ty in the anti-cutting direction Y of the rotary tool 5 stored in the tool displacement amount storage unit 36 (S21). ). Specifically, first, as shown in the equation (1), the dynamic characteristic calculation unit 37 detects the interval between the times t1 and t2 of the maximum displacements Ty1 and Ty2 of the rotary tool 15, and the period T [ s] (shown in FIG. 7) is detected.

(数1)
T[s]=t2−t1・・・(1)
(Equation 1)
T [s] = t2-t1 (1)

次に、動特性算出部37は、式(2)に示すように、変位Tyの周期Tの逆数から、変位Tyの固有振動数fd[Hz]を算出する。   Next, the dynamic characteristic calculation unit 37 calculates the natural frequency fd [Hz] of the displacement Ty from the reciprocal of the period T of the displacement Ty, as shown in Expression (2).

(数2)
fd=1/T・・・(2)
fd:固有振動数
(Equation 2)
fd = 1 / T (2)
fd: natural frequency

次に、動特性算出部37は、式(3)に示すように、変位Tyの角固有振動数ωd[rad/s]を算出する(S21)。   Next, the dynamic characteristic calculation unit 37 calculates the angular natural frequency ωd [rad / s] of the displacement Ty as shown in Expression (3) (S21).

(数3)
ωd=2πfd・・・(3)
ωd:角固有振動数
(Equation 3)
ωd = 2πfd (3)
ωd: angular natural frequency

次に、動特性算出部37は、工具変位量記憶部36の記憶されている回転工具5の反切込方向Yの連続的な変位Tyに基づき、式(4)に示すように、回転工具15の最大変位Ty1、Ty2(図7示)から、回転工具5の振動の減衰率ζを算出する(S22)。   Next, based on the continuous displacement Ty in the anti-cutting direction Y of the rotary tool 5 stored in the tool displacement amount storage unit 36, the dynamic characteristic calculation unit 37, as shown in Expression (4), the rotary tool 15 The vibration damping rate ζ of the rotary tool 5 is calculated from the maximum displacements Ty1, Ty2 (shown in FIG. 7) (S22).

Figure 0005927904
Figure 0005927904

次に、動特性算出部37は、工具変位量記憶部36の記憶されている回転工具5の反切込方向Yの連続的な変位Tyに基づき、一方の刃部5aが被加工物Wに接触開始した瞬間の時刻t0(図9示)から一方の刃部5aによる被加工物Wの切削加工を終了した瞬間の時刻t1(図9示)までの変位Tyの積分値Ak(図9示)を算出する(S23)。   Next, the dynamic characteristic calculation unit 37 makes contact with the workpiece W with one blade portion 5a based on the continuous displacement Ty in the anti-cutting direction Y of the rotary tool 5 stored in the tool displacement amount storage unit 36. The integrated value Ak (shown in FIG. 9) of the displacement Ty from the time t0 (shown in FIG. 9) at the start time to the time t1 (shown in FIG. 9) at the moment when cutting of the workpiece W by the one blade portion 5a is completed. Is calculated (S23).

次に、動特性算出部37は、工具変位量記憶部36の記憶されている回転工具5の反切込方向Yの連続的な変位Tyを時間微分することにより、回転工具5の変形速度Vy(図9示)を算出するとともに、変形速度Vyを上記時刻t0から時刻t1まで時間積分し、回転工具5の変形速度Vyの積分値Ac(図9示)を算出する(S24)。   Next, the dynamic characteristic calculation unit 37 time-differentiates the continuous displacement Ty in the anti-cutting direction Y of the rotary tool 5 stored in the tool displacement amount storage unit 36, so that the deformation speed Vy ( 9 is calculated, and the deformation speed Vy is integrated over time from the time t0 to the time t1, and an integrated value Ac (shown in FIG. 9) of the deformation speed Vy of the rotary tool 5 is calculated (S24).

次に、動特性算出部37は、切削抵抗算出部32により算出された切削抵抗Fyを上述した時刻t0から時刻t1まで時間積分し、回転工具5の切削抵抗Fyの積分値At(図9示)を算出する(S25)。   Next, the dynamic characteristic calculation unit 37 integrates the cutting resistance Fy calculated by the cutting resistance calculation unit 32 from the time t0 to the time t1 described above, and an integrated value At (see FIG. 9) of the cutting resistance Fy of the rotary tool 5. ) Is calculated (S25).

次に、動特性算出部37は、質量係数Mを算出する(S26)。以下に、質量係数M算出手法を説明する。
被加工物Wを切削加工中の回転工具5に作用する力Fは、回転工具5に作用する反切込方向Yの切削抵抗Fy、回転工具5の弾性による力Fk、及び、回転工具5の粘性減衰による力Fcの合力であるので、下式(5)で表される。
Next, the dynamic characteristic calculation unit 37 calculates the mass coefficient M (S26). Below, the mass coefficient M calculation method is demonstrated.
The force F acting on the rotary tool 5 during cutting of the workpiece W includes the cutting resistance Fy in the anti-cutting direction Y acting on the rotary tool 5, the force Fk due to the elasticity of the rotary tool 5, and the viscosity of the rotary tool 5. Since it is the resultant force Fc due to damping, it is expressed by the following equation (5).

(数5)
F=Fy+Fk+Fc・・・(5)

F:回転工具5に作用する反切込方向Yの力
Fy:回転工具5に作用する反切込方向Yの切削抵抗
Fk:回転工具5の弾性による力
Fc:回転工具5の粘性減衰による力
(Equation 5)
F = Fy + Fk + Fc (5)

F: Force in the anti-cutting direction Y acting on the rotary tool 5
Fy: Cutting resistance in the anti-cutting direction Y acting on the rotary tool 5
Fk: force due to elasticity of the rotary tool 5
Fc: Force due to viscous damping of the rotary tool 5

ここで、Fkは下式(6)で表される。   Here, Fk is expressed by the following formula (6).

(数6)
Fk=−KTy・・・(6)

Fk:回転工具5の弾性による力
K:回転工具5の弾性係数
Ty:回転工具5の反切込方向Yの変位Ty
(Equation 6)
Fk = −KTy (6)

Fk: force due to elasticity of the rotary tool 5
K: elastic coefficient of the rotary tool 5
Ty: Displacement Ty in the anti-cutting direction Y of the rotary tool 5

また、Fcは下式(7)で表される。   Fc is represented by the following formula (7).

(数7)
Fc=−CVy・・・(7)

Fc:回転工具5の粘性減衰による力
C:回転工具5の粘性減衰係数
Vy:回転工具5の反切込方向Yの変形速度Vy
(Equation 7)
Fc = −CVy (7)

Fc: Force due to viscous damping of the rotary tool 5
C: Viscous damping coefficient of rotating tool 5
Vy: Deformation speed Vy in the anti-cutting direction Y of the rotary tool 5

次に、回転工具5の質量係数Mとし、回転工具5の加速度をaとすると、上式(5)及び下式(8)から式(9)が導かれる。   Next, when the mass coefficient M of the rotary tool 5 is assumed and the acceleration of the rotary tool 5 is a, the formula (9) is derived from the above formula (5) and the following formula (8).

(数8)
F=Ma・・・(8)

F:回転工具5に作用する反切込方向Yの力
M:回転工具5の質量係数
a:回転工具5の加速度
(Equation 8)
F = Ma (8)

F: Force in the anti-cutting direction Y acting on the rotary tool 5
M: Mass coefficient of the rotary tool 5
a: Acceleration of the rotary tool 5

(数9)
a=(Fy+Fk+Fc)/M・・・(9)
(Equation 9)
a = (Fy + Fk + Fc) / M (9)

ここで、上記時刻t0での回転工具5の変形速度をV0(図9示)、上記時刻t1での回転工具5の変形速度をV1(図9示)とすると、その間の回転工具5の速度変化量ΔVy=(V1−V0)は、回転工具5の加速度aを上記時刻t0から時刻t1まで時間積分した値であるので、上式(9)に上式(6)、(7)を代入し、上記時刻t0から時刻t1まで時間積分すると、下式(10)が導かれる。   Here, if the deformation speed of the rotary tool 5 at the time t0 is V0 (shown in FIG. 9) and the deformation speed of the rotary tool 5 at the time t1 is V1 (shown in FIG. 9), the speed of the rotary tool 5 during that time Since the change amount ΔVy = (V1−V0) is a value obtained by integrating the acceleration a of the rotary tool 5 from time t0 to time t1, the above expressions (6) and (7) are substituted into the above expression (9). Then, when time integration is performed from the time t0 to the time t1, the following expression (10) is derived.

Figure 0005927904
Figure 0005927904

上式(10)に、S25において算出した切削抵抗Fyの積分値At、S23において算出した変位Tyの積分値Ak、及び、S24において算出した変形速度Vyの積分値Acを代入すると、下式(11)が導かれる。   Substituting the integrated value At of the cutting force Fy calculated in S25, the integrated value Ak of the displacement Ty calculated in S23, and the integrated value Ac of the deformation speed Vy calculated in S24 into the above formula (10), 11) is derived.

(数11)
MΔVy=At−KAk-CAc…(11)
(Equation 11)
MΔVy = At−KAk−CAc (11)

ここで、減衰率ζが十分に小さいとすると、減衰振動の運動方程式から、角固有振動数ωdは、下式(12)で表される。   Here, if the damping rate ζ is sufficiently small, the angular natural frequency ωd is expressed by the following equation (12) from the equation of motion of the damping vibration.

(数12)
ωd=√(K/M)…(12)
ωd:角固有振動数
K:弾性係数
C:粘性減衰係数
(Equation 12)
ωd = √ (K / M) (12)
ωd: angular natural frequency
K: elastic modulus
C: Viscous damping coefficient

上式(12)を変形させると、弾性係数Kは下式(13)で表される。   When the above equation (12) is deformed, the elastic coefficient K is expressed by the following equation (13).

(数13)
K=ωdM…(13)
(Equation 13)
K = ωd 2 M (13)

一方で、減衰振動の運動方程式から、減衰率ζは、下式(14)で表される。   On the other hand, the damping rate ζ is expressed by the following equation (14) from the equation of motion of the damped vibration.

(数14)
ζ=C/2√(K・M)…(14)
ζ:減衰率
K:弾性係数
C:粘性減衰係数
(Equation 14)
ζ = C / 2√ (K · M) (14)
ζ: Decay rate
K: elastic modulus
C: Viscous damping coefficient

上式(14)を変形し、上式(12)を代入すると、粘性減衰係数Cは、下式(15)で表される。   When the above equation (14) is transformed and the above equation (12) is substituted, the viscous damping coefficient C is expressed by the following equation (15).

(数15)
C=2ζ√(K・M)=2ζωdM…(15)
(Equation 15)
C = 2ζ√ (K · M) = 2ζωdM (15)

次に、上式(11)に、上式(13)及び上式(15)を代入すると、下式(16)が導き出される。   Next, when the above equation (13) and the above equation (15) are substituted into the above equation (11), the following equation (16) is derived.

(数16)
At=MΔVy+MωdAk+2MζωdAc…(16)
(Equation 16)
At = MΔVy + Mωd 2 Ak + 2MζωdAc (16)

そして、上式(16)を変形すると、下式(17)が導きだされる。
(数17)
M=At/(ΔVy+ωdAk+2ζωdAc)…(17)
And if the above equation (16) is transformed, the following equation (17) is derived.
(Equation 17)
M = At / (ΔVy + ωd 2 Ak + 2ζωdAc) (17)

次に、上式(17)に、S21で算出した角固有振動数ωd、S22で算出した減衰率ζを代入すると、回転工具5の質量係数Mが算出される。   Next, when the angular natural frequency ωd calculated in S21 and the damping factor ζ calculated in S22 are substituted into the above equation (17), the mass coefficient M of the rotary tool 5 is calculated.

次に、動特性算出部37は、上式(13)に、S21で算出した角固有振動数ωd、及び、S26で算出した質量係数Mを代入することにより、回転工具5の弾性係数Kを算出する(S27)。   Next, the dynamic characteristic calculator 37 substitutes the elastic coefficient K of the rotary tool 5 by substituting the angular natural frequency ωd calculated in S21 and the mass coefficient M calculated in S26 into the above equation (13). Calculate (S27).

最後に、動特性算出部37は、上式(15)に、S21で算出した角固有振動数ωd、S22で算出した減衰率ζ、S26で算出した質量係数Mを代入することにより、回転工具5の粘性減衰係数Cを算出する(S28)。   Finally, the dynamic characteristic calculation unit 37 substitutes the angular natural frequency ωd calculated in S21, the damping factor ζ calculated in S22, and the mass coefficient M calculated in S26 into the above equation (15), thereby rotating the rotary tool. 5 is calculated (S28).

工具動特性記憶部38は、動特性算出部37によって算出された回転工具5の固有振動数fd、角固有振動数ωd、減衰率ζ、質量係数M、粘性減衰係数C、弾性係数(バネ乗数)Kからなる回転工具5の動特性を記憶する。   The tool dynamic characteristic storage unit 38 includes a natural frequency fd, an angular natural frequency ωd, a damping factor ζ, a mass coefficient M, a viscous damping coefficient C, an elastic coefficient (spring multiplier) calculated by the dynamic characteristic calculating unit 37. ) The dynamic characteristics of the rotary tool 5 made of K are stored.

加工条件決定処理部39は、工具動特性記憶部38に記憶されている動特性に基づいて、切削加工の加工条件を決定する。例えば、加工条件決定処理部39は、回転工具5の固有振動数fd(共振周波数)に基づいて、回転工具5の固有振動数fdを避けるように、回転工具5の回転速度S(回転主軸4の回転速度)を決定する。なお、加工条件決定処理部39は、切削加工の加工条件を決定して、NCデータそのものを変更することもできるし、切削加工中に切削加工の加工条件を補正することもできる。   The machining condition determination processing unit 39 determines the machining conditions for cutting based on the dynamic characteristics stored in the tool dynamic characteristic storage unit 38. For example, the machining condition determination processing unit 39, based on the natural frequency fd (resonance frequency) of the rotary tool 5, avoids the natural frequency fd of the rotary tool 5, so that the rotational speed S (rotational spindle 4) of the rotary tool 5 is avoided. The rotation speed). Note that the machining condition determination processing unit 39 can determine the machining conditions for the cutting process, change the NC data itself, and can correct the machining conditions for the cutting process during the cutting process.

機械制御部41は、NCデータ、及び、加工条件決定処理部39にて決定された切削加工の加工条件に基づいて、各駆動部61を制御する。駆動部61には、回転主軸4を回転させるモータ、回転主軸4をX軸、Y軸、Z軸方向に移動するアクチュエータが含まれる。   The machine control unit 41 controls each drive unit 61 based on the NC data and the machining conditions determined by the machining condition determination processing unit 39. The drive unit 61 includes a motor that rotates the rotation main shaft 4 and an actuator that moves the rotation main shaft 4 in the X-axis, Y-axis, and Z-axis directions.

(4.動特性算出・加工処理の説明)
次に、図6及び図10のフローチャートを参照して、図1や図6に示す制御部50の各算出部、加工条件決定処理部39、及び、機械制御部41が実行する動特性算出・加工処理について説明する。まず、機械制御部41は、駆動部61に制御信号を出力することにより、試し加工を実行する(S101)。この試し加工は、製品の品質に影響しない箇所の加工や仕上げ加工前の粗加工、或いは、テストピースによる切削が含まれる。次に、工具変位量算出部35は、上述した処理によって、高速度カメラ9が撮像した回転工具15を含む撮像データに基づいて、回転工具5の反切込方向Yの変位Tyを算出して、試し加工中において加振され振動している回転工具5の変位Tyを検出する(S102)。
(4. Description of dynamic characteristics calculation and processing)
Next, referring to the flowcharts of FIGS. 6 and 10, calculation of dynamic characteristics executed by each calculation unit of the control unit 50, the machining condition determination processing unit 39, and the machine control unit 41 shown in FIGS. 1 and 6 is performed. Processing will be described. First, the machine control unit 41 performs a trial machining by outputting a control signal to the drive unit 61 (S101). This trial processing includes processing of a portion that does not affect the quality of the product, rough processing before finishing processing, or cutting with a test piece. Next, the tool displacement amount calculation unit 35 calculates the displacement Ty in the anti-cutting direction Y of the rotary tool 5 based on the imaging data including the rotary tool 15 imaged by the high-speed camera 9 by the above-described processing. The displacement Ty of the rotating tool 5 that is vibrated and vibrated during the trial machining is detected (S102).

次に、切削抵抗算出部32は、上述した処理によって、試し加工中の切削抵抗Fyを算出する(S103)。次に、動特性算出部37は、上述した処理によって、S102で検出した試し加工中の回転工具5の反切込方向Yの変位Ty、及び、S103で算出した試し加工中の切削抵抗Fyに基づいて、回転工具5の固有振動数fd、角固有振動数ωd、減衰率ζ、質量係数M、粘性減衰係数C、及び、弾性係数Kからなる回転工具5の動特性を算出する(S103)。   Next, the cutting resistance calculation unit 32 calculates the cutting resistance Fy during the trial machining by the above-described processing (S103). Next, the dynamic characteristic calculation unit 37 is based on the displacement Ty in the anti-cutting direction Y of the rotary tool 5 during the trial machining detected in S102 and the cutting resistance Fy during the trial machining calculated in S103 by the above-described processing. Thus, the dynamic characteristics of the rotary tool 5 including the natural frequency fd, the angular natural frequency ωd, the damping factor ζ, the mass coefficient M, the viscous damping coefficient C, and the elastic coefficient K of the rotary tool 5 are calculated (S103).

次に、加工条件決定処理部39は、S103で算出された動特性に基づいて、加工条件を決定する(S105)。次に、機械制御部41は、NCデータ、及び、S105で決定された加工条件に基づいて、駆動部61に制御信号を出力することにより、本加工を実行し(S106)、製品を完成させる。このように、S101において、最終製品の品質に影響しない箇所の被加工物Wを切削加工することにより、S102において試し加工中の回転工具5の変位量Tyを検出し、S104において回転工具5の動特性を算出し、S105において、算出された動特性によって加工条件を決定し、S106において、加工条件が最適化された状態で、本加工を行うので、最終製品の品質に影響することなく、より精度高く切削加工を行うことができる。   Next, the machining condition determination processing unit 39 determines a machining condition based on the dynamic characteristic calculated in S103 (S105). Next, the machine control unit 41 outputs the control signal to the drive unit 61 based on the NC data and the processing conditions determined in S105, thereby executing the main processing (S106) and completing the product. . Thus, in S101, by cutting the workpiece W at a location that does not affect the quality of the final product, in S102, the displacement amount Ty of the rotary tool 5 being trial-worked is detected. The dynamic characteristics are calculated, and in S105, the processing conditions are determined based on the calculated dynamic characteristics. In S106, the main processing is performed with the processing conditions optimized, so that the quality of the final product is not affected. Cutting can be performed with higher accuracy.

上述した動特性算出装置によれば、以下のような効果を奏する。図6に示す動特性算出部37(動特性算出手段)は、工具変位量算出部35が算出した切削加工中の回転工具5の反切込方向Yの変位量Ty、切削抵抗算出部32によって算出された回転工具5の反切込方向Yの切削抵抗Fyに基づいて、回転工具5の動特性を算出する。これにより、回転主軸4が回転中の回転工具5の動特性が算出される。このため、回転主軸4が停止している状態で動特性を算出する場合と比較して、正確に動特性を算出して取得することが可能となる。また、高速度カメラ9(工具変位量検出手段)は、断続的な切削加工により加振された回転工具5の反切込方向Yの変位量Tyを検出し、当該検出された回転工具5の変位量Tyに基づいて、動特性を算出する。このため、従来のように、回転工具5に加速度計を取り付ける必要もなく、ハンマーの叩き方にも経験を要さないことから、手間をかけずに動特性を取得することができる。   According to the above-described dynamic characteristic calculation apparatus, the following effects can be obtained. The dynamic characteristic calculation unit 37 (dynamic characteristic calculation unit) illustrated in FIG. 6 is calculated by the displacement amount Ty in the anti-cutting direction Y of the rotary tool 5 during cutting calculated by the tool displacement amount calculation unit 35 and the cutting resistance calculation unit 32. Based on the cutting resistance Fy of the rotating tool 5 in the anti-cutting direction Y, the dynamic characteristic of the rotating tool 5 is calculated. As a result, the dynamic characteristics of the rotary tool 5 while the rotary spindle 4 is rotating are calculated. For this reason, it is possible to accurately calculate and acquire the dynamic characteristic as compared with the case where the dynamic characteristic is calculated in a state where the rotation main shaft 4 is stopped. The high-speed camera 9 (tool displacement amount detection means) detects the displacement amount Ty in the anti-cutting direction Y of the rotary tool 5 vibrated by intermittent cutting, and the detected displacement of the rotary tool 5 is detected. Based on the quantity Ty, the dynamic characteristic is calculated. For this reason, it is not necessary to attach an accelerometer to the rotary tool 5 as in the prior art, and since no experience is required in how to strike the hammer, the dynamic characteristics can be acquired without taking time and effort.

また、被加工物Wの試し加工中において、高速度カメラ9(工具変位量検出手段)によって検出された振動している状態の回転工具5の変位量に基づいて、動特性算出部37(動特性算出手段)は、回転工具5の動特性を算出する。これにより、最終製品の品質に影響すること無く、回転工具5の動特性を算出することができる。   Further, during the trial machining of the workpiece W, based on the displacement amount of the rotating tool 5 in the vibrating state detected by the high-speed camera 9 (tool displacement amount detection means), the dynamic characteristic calculation unit 37 (dynamic (Characteristic calculating means) calculates the dynamic characteristics of the rotary tool 5. Thereby, the dynamic characteristics of the rotary tool 5 can be calculated without affecting the quality of the final product.

また、本実施形態では、高速度カメラ9(撮像装置)によって、切削加工中の回転工具5の変位量を検出している。そして、高速度カメラ9によって撮像された回転工具5を含む画像データから、回転工具5の外縁を検出し、回転工具5の反切込方向Yの変位Tyを検出する。このため、従来のように、加速度センサ等によって回転工具5の変位を検出する方法と比較して、回転主軸4の回転に伴う振動等のノイズの影響を受けず、より高精度に回転工具5の反切込方向Yの変位を検出することができ、微少な回転工具5の変位も検出することができる。また、加速度センサを取り付けることができない小径の回転工具5であっても、回転工具5の動特性を取得することができる。   In this embodiment, the amount of displacement of the rotary tool 5 during cutting is detected by the high-speed camera 9 (imaging device). Then, the outer edge of the rotary tool 5 is detected from the image data including the rotary tool 5 imaged by the high-speed camera 9, and the displacement Ty in the anti-cutting direction Y of the rotary tool 5 is detected. For this reason, compared with the conventional method of detecting the displacement of the rotary tool 5 by an acceleration sensor or the like, the rotary tool 5 is not affected by noise such as vibration accompanying the rotation of the rotary spindle 4 and more accurately. The displacement in the anti-cutting direction Y can be detected, and a slight displacement of the rotary tool 5 can also be detected. Moreover, even if it is the small diameter rotary tool 5 which cannot attach an acceleration sensor, the dynamic characteristic of the rotary tool 5 is acquirable.

以上説明した実施形態では、被加工物Wを切削加工中の回転工具5の反切込方向Yの変位量Tyを検出する工具変位量検出手段として、高速度カメラ9を用いている。しかし、前記工具変位量検出手段として、レーザ発光素子とレーザ受光素子からなり、レーザ発光素子が切削加工中の回転工具5にレーザを照射し、回転工具5によって反射されレーザ受光素子によって受光したレーザにより、切削加工中の回転工具5の反切込方向Yの変位量Tyを検出するレーザ式変位センサであっても差し支え無い。或いは、前記工具変位量検出手段として、回転主軸4内に設けられ、回転主軸4の変位量を検出する主軸内部変位センサと、回転主軸4の変位量より切削加工中の回転工具5の反切込方向Yの変位量Tyを算出する工具変位量算出部とからなる主軸内部変位センサシステムであっても差し支え無い。   In the embodiment described above, the high-speed camera 9 is used as the tool displacement amount detection means for detecting the displacement amount Ty in the anti-cutting direction Y of the rotary tool 5 during cutting the workpiece W. However, the tool displacement amount detecting means includes a laser light emitting element and a laser light receiving element, and the laser light emitting element irradiates the rotating tool 5 being cut with laser, and is reflected by the rotating tool 5 and received by the laser light receiving element. Thus, a laser displacement sensor that detects the displacement amount Ty in the anti-cutting direction Y of the rotary tool 5 during cutting may be used. Alternatively, as the tool displacement amount detection means, a spindle internal displacement sensor provided in the rotation main shaft 4 for detecting the displacement amount of the rotation main shaft 4 and the anti-cutting of the rotary tool 5 being cut from the displacement amount of the rotation main shaft 4. There may be a spindle internal displacement sensor system including a tool displacement amount calculation unit for calculating the displacement amount Ty in the direction Y.

なお、回転工具5の共振点が複数ある場合には、回転工具5の反切込方向Yの変位Tyのうち、回転工具5の回転周期の整数倍分だけ切り出して、拡散フーリエ展開(DFT)によって回転工具5の動特性を算出することが好ましい。   In addition, when there are a plurality of resonance points of the rotary tool 5, only the integral multiple of the rotation period of the rotary tool 5 is cut out from the displacement Ty in the anti-cutting direction Y of the rotary tool 5, and the diffusion Fourier expansion (DFT) is performed. It is preferable to calculate the dynamic characteristics of the rotary tool 5.

5:回転工具、 5a,5b:刃部、 9…高速度カメラ(撮像装置、工具変位量検出手段)、 32:切削抵抗算出部(切削抵抗算出手段)、 37:動特性算出部(動特性算出手段)、 50:制御部、 100…工作機械、 W:被加工物   5: Rotary tool, 5a, 5b: Blade part, 9 ... High-speed camera (imaging device, tool displacement amount detection means), 32: Cutting resistance calculation part (cutting resistance calculation means), 37: Dynamic characteristic calculation part (dynamic characteristic) Calculation means), 50: control unit, 100 ... machine tool, W: workpiece

Claims (3)

外周側に周方向に1以上の刃部を備える回転工具を用いて、当該回転工具を軸回りに回転しながら被加工物に対して相対移動して断続的な切削加工を行う工作機械の動特性を算出する装置であって、
断続的な切削加工により前記回転工具が加振された場合に、振動している状態の前記回転工具の変位量を検出する工具変位量検出手段と、
前記回転工具の切削抵抗を算出する切削抵抗算出手段と、
前記工具変位量検出手段が検出した振動している状態の前記回転工具の前記変位量、及び、前記切削抵抗算出手段によって算出された前記回転工具の前記切削抵抗に基づいて、前記回転工具の動特性である質量係数、弾性係数、及び粘性減衰係数の少なくとも1つを算出する動特性算出手段と、
を備える工作機械の動特性算出装置。
Using a rotary tool having one or more blades in the circumferential direction on the outer peripheral side, the machine tool moves intermittently by moving relative to the workpiece while rotating the rotary tool about its axis. A device for calculating characteristics,
A tool displacement amount detecting means for detecting a displacement amount of the rotating tool in a vibrating state when the rotating tool is vibrated by intermittent cutting;
Cutting force calculating means for calculating cutting resistance of the rotary tool;
The displacement amount of the rotary tool in a state in which the tool displacement amount detecting means is vibrating the detected and, based on the cutting force of the rotary tool calculated by the cutting resistance calculating means, movement of the rotary tool Dynamic characteristic calculating means for calculating at least one of a characteristic mass coefficient, elastic coefficient, and viscous damping coefficient ;
A machine tool dynamic characteristic calculation apparatus comprising:
請求項1において、
前記被加工物の試し加工中において、前記工具変位量検出手段によって検出された振動している状態の前記回転工具の前記変位量に基づいて、前記動特性算出手段は、前記回転工具の前記動特性を算出する工作機械の動特性算出装置。
In claim 1,
During trial machining of the workpiece, on the basis of the displacement amount of the rotary tool in the state that vibrations detected by said tool displacement amount detecting means, the dynamic characteristic calculating means, the movement of the rotary tool Machine tool dynamic characteristic calculation device for calculating characteristics.
請求項1又は請求項2において、
前記工具変位量検出手段は、撮像装置である工作機械の動特性算出装置。
In claim 1 or claim 2,
The tool displacement amount detection means is a dynamic characteristic calculation device of a machine tool which is an imaging device.
JP2011283254A 2011-12-26 2011-12-26 Machine tool dynamic characteristics calculation device Active JP5927904B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2011283254A JP5927904B2 (en) 2011-12-26 2011-12-26 Machine tool dynamic characteristics calculation device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2011283254A JP5927904B2 (en) 2011-12-26 2011-12-26 Machine tool dynamic characteristics calculation device

Publications (2)

Publication Number Publication Date
JP2013132706A JP2013132706A (en) 2013-07-08
JP5927904B2 true JP5927904B2 (en) 2016-06-01

Family

ID=48909775

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2011283254A Active JP5927904B2 (en) 2011-12-26 2011-12-26 Machine tool dynamic characteristics calculation device

Country Status (1)

Country Link
JP (1) JP5927904B2 (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6330324B2 (en) * 2013-12-26 2018-05-30 株式会社ジェイテクト Machine tool dynamic characteristic calculation device and dynamic characteristic calculation method
EP3225355B1 (en) * 2014-11-26 2019-09-11 Mitsubishi Electric Corporation Numerical control apparatus
JP6578195B2 (en) * 2015-11-26 2019-09-18 Dmg森精機株式会社 Method for deriving natural frequency of cutting tool, method for creating stability limit curve, and device for deriving natural frequency of cutting tool
JP6266144B1 (en) 2017-01-27 2018-01-24 Dmg森精機株式会社 Vibration state detection device and machine tool provided with vibration state detection device
JP6549646B2 (en) 2017-07-06 2019-07-24 ファナック株式会社 Machine tool and home position correction method
DE102018005612B4 (en) 2017-07-19 2024-05-16 Fanuc Corporation Reporting of violations
CN110032794B (en) * 2019-04-10 2022-06-24 哈尔滨理工大学 Method for constructing and verifying dynamic cutting force model of milling cutter under vibration action
WO2023047437A1 (en) * 2021-09-21 2023-03-30 株式会社ジェイテクト Processing estimation device
WO2024038750A1 (en) * 2022-08-15 2024-02-22 ソニーセミコンダクタソリューションズ株式会社 Vibration monitoring device and machining system
WO2024075284A1 (en) * 2022-10-07 2024-04-11 株式会社ジェイテクト Contact dynamic stiffness calculation system, machining estimation device, and proccessing system
CN116765639B (en) * 2023-08-21 2023-11-14 山东中玮动力科技有限公司 High-precision five-axis laser processing machine tool

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001225243A (en) * 2000-02-16 2001-08-21 Toshiba Corp Machining work optimizing device
JP4434137B2 (en) * 2005-12-20 2010-03-17 株式会社日立製作所 Calculation method of axial depth of cut for stability limit of self-excited vibration in cutting of inclined surface

Also Published As

Publication number Publication date
JP2013132706A (en) 2013-07-08

Similar Documents

Publication Publication Date Title
JP5927904B2 (en) Machine tool dynamic characteristics calculation device
US11733673B2 (en) Numerical control system and motor drive controller
JP5189806B2 (en) Surface shape measuring device
JP6445070B2 (en) Machine tool control system
JP6192849B2 (en) Vibration mode measuring device
JP2007315897A (en) Measurement device, and surface property measurement method and program thereof
CN103522125B (en) Dynamic characteristic computing device and dynamic characteristic computational methods for lathe
Salehi et al. Indirect model based estimation of cutting force and tool tip vibrational behavior in milling machines by sensor fusion
JP5983112B2 (en) Machine tool dynamic characteristic calculation device and dynamic characteristic calculation method
JP5748412B2 (en) Cutting tool edge diagnosis device for machine tools
JP5983113B2 (en) Machine tool dynamic characteristic calculation device and dynamic characteristic calculation method
TWI638251B (en) Modal detection system
CN111836697B (en) Balance and vibration adjustment system for rotary tool, balance and vibration measurement device, balance and vibration adjustment method, and tool holder
JP2020189388A (en) Main spindle abnormality detection device
JP2000310527A (en) Surface-properties measuring device
JP6330324B2 (en) Machine tool dynamic characteristic calculation device and dynamic characteristic calculation method
JP7479028B2 (en) Method and apparatus for measuring changes in vibration characteristics of a machine tool
JP2007007752A (en) Displacement correcting method of spindle tool tip
Jemielniak et al. New method of dynamic cutting force coefficients determination
JP7190343B2 (en) Three-dimensional measuring device and measuring method
JP2021133473A (en) Cut product manufacturing method
JP7094509B1 (en) Information processing equipment, control programs and information processing methods
KR102191653B1 (en) System for in-process monitoring of rotating tools using phase-shifted stroboscopy
US10620164B2 (en) Diagnostic device and diagnostic method
JP2024068895A (en) Processing condition setting method and processing condition setting device

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20141124

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20150928

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20151104

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20151217

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160329

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160411

R150 Certificate of patent (=grant) or registration of utility model

Ref document number: 5927904

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150