JP2021058945A - Machine tool, numerical control device and vibration suppression method - Google Patents
Machine tool, numerical control device and vibration suppression method Download PDFInfo
- Publication number
- JP2021058945A JP2021058945A JP2019183150A JP2019183150A JP2021058945A JP 2021058945 A JP2021058945 A JP 2021058945A JP 2019183150 A JP2019183150 A JP 2019183150A JP 2019183150 A JP2019183150 A JP 2019183150A JP 2021058945 A JP2021058945 A JP 2021058945A
- Authority
- JP
- Japan
- Prior art keywords
- vibration
- rotation speed
- work
- time
- frequency
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000001629 suppression Effects 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 80
- 238000003754 machining Methods 0.000 claims abstract description 39
- 238000012545 processing Methods 0.000 claims description 17
- 230000008569 process Effects 0.000 claims description 15
- 230000000694 effects Effects 0.000 claims description 14
- 244000145845 chattering Species 0.000 abstract description 2
- 238000005520 cutting process Methods 0.000 description 69
- 230000006870 function Effects 0.000 description 28
- 238000006073 displacement reaction Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 8
- 230000008929 regeneration Effects 0.000 description 6
- 238000011069 regeneration method Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 4
- 239000000284 extract Substances 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 230000001172 regenerating effect Effects 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 208000000705 Rift Valley Fever Diseases 0.000 description 1
- 206010039163 Right ventricular failure Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013017 mechanical damping Methods 0.000 description 1
- 238000003801 milling Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Automatic Control Of Machine Tools (AREA)
- Numerical Control (AREA)
Abstract
Description
本発明による実施形態は、工作機械、数値制御装置および振動抑制方法に関する。 Embodiments of the present invention relate to machine tools, numerical control devices and vibration suppression methods.
工作機械は、例えば、主軸に装着された工具を用いてワークを加工する。工作機械の数値制御装置（以下、ＮＣ（Numerical Controller）装置ともいう）は、主軸に指令を出力し、該主軸の動作を制御する。例えば、工作機械を用いたボーリング切削中に工具とワークとの間にびびり振動が発生すると、ワークの加工面にびびり模様が発生し、加工面品質が劣化する。また、工具の寿命が短縮したり、工具が破損したりする場合もある。従って、びびり振動が発生した場合、速やかにびびり振動を消散させる必要がある。このびびり振動を抑制する方法の一つとして、主軸の回転速度を変化させる方法（ＳＳＶ（Spindle Speed Variation））が知られている。 The machine tool processes the work using, for example, a tool mounted on the spindle. A numerical control device for a machine tool (hereinafter, also referred to as an NC (Numerical Controller) device) outputs a command to a spindle and controls the operation of the spindle. For example, if chatter vibration occurs between a tool and a work during boring cutting using a machine tool, a chatter pattern is generated on the machined surface of the work, and the quality of the machined surface deteriorates. In addition, the life of the tool may be shortened or the tool may be damaged. Therefore, when chatter vibration occurs, it is necessary to promptly dissipate the chatter vibration. As one of the methods for suppressing this chatter vibration, a method for changing the rotation speed of the spindle (SSV (Spindle Speed Variation)) is known.
しかし、ＳＳＶであっても、びびり振動を適切に抑制することが困難な場合があった。この場合、回転速度の変動周期に応じた過渡振動（Beat vibration）が発生する可能性がある。 However, even with SSV, it may be difficult to appropriately suppress chatter vibration. In this case, transient vibration (Beat vibration) may occur according to the fluctuation period of the rotation speed.
びびり振動を抑制することができる工作機械、数値制御装置および振動抑制方法を提供することを目的とする。 It is an object of the present invention to provide a machine tool capable of suppressing chatter vibration, a numerical control device, and a vibration suppression method.
本実施形態による工作機械は、ワークまたは該ワークを加工する工具を回転させる回転軸と、回転軸の回転速度を検出する回転速度センサと、加工中の工具またはワークに発生する振動を検出する振動センサと、回転軸の回転速度を制御する数値制御装置と、を備え、数値制御装置は、回転速度および振動の振動周波数を取得するセンサ情報取得部と、振動および振動のエネルギーをベッセル関数を用いて近似し、ワークを加工する第１回転速度を中心とする回転速度の変動幅および変動周波数を演算し、変動幅および変動周波数を第１回転速度に重畳させてワークを加工する第２回転速度を演算する演算部と、を有する。 The machine tool according to the present embodiment has a rotating shaft for rotating the work or a tool for machining the work, a rotation speed sensor for detecting the rotation speed of the rotating shaft, and vibration for detecting the vibration generated in the tool or the work being machined. It includes a sensor and a numerical control device that controls the rotation speed of the rotation axis. The numerical control device uses a sensor information acquisition unit that acquires the rotation speed and vibration frequency of vibration, and a vessel function for vibration and vibration energy. The fluctuation width and fluctuation frequency of the rotation speed centered on the first rotation speed for processing the work are calculated, and the fluctuation width and fluctuation frequency are superimposed on the first rotation speed to process the work. It has a calculation unit for calculating.
以下、図面を参照して本発明に係る実施形態を説明する。本実施形態は、本発明を限定するものではない。図面は模式的または概念的なものであり、各部分の比率などは、必ずしも現実のものと同一とは限らない。明細書と図面において、既出の図面に関して前述したものと同様の要素には同一の符号を付して詳細な説明は適宜省略する。 Hereinafter, embodiments according to the present invention will be described with reference to the drawings. The present embodiment is not limited to the present invention. The drawings are schematic or conceptual, and the ratio of each part is not always the same as the actual one. In the specification and the drawings, the same elements as those described above with respect to the existing drawings are designated by the same reference numerals, and detailed description thereof will be omitted as appropriate.
（第１実施形態）
図１は、第１実施形態による工作機械１００の構成の一例を示すブロック図である。工作機械１００は、主軸１と、主軸モータ２と、エンコーダ３と、振動センサ４と、数値制御装置５と、サーボ制御部６と、を備える。
(First Embodiment)
FIG. 1 is a block diagram showing an example of the configuration of the
回転軸としての主軸１は、ワークＷを加工する工具１１が設けられ、工具１１を回転させる。図１に示す例では、工具１１は、ワークＷのボーリング切削に用いられる。従って、工具１１は、回転しながらＺ軸方向に移動し、ワークＷを切削する。この場合、ワークＷは固定されている。また、工具１１の刃数は、例えば、１である。この場合、主軸１回転ごとに、ワークＷが切削される。尚、工具１１は、複数の刃を有する多刃工具であってもよい。
The
主軸モータ２は、サーボ制御部６からサーボ指令を受けて、主軸１を回転駆動させる。
The
回転速度センサとしてのエンコーダ３は、主軸１の回転速度ｎを検出する。エンコーダ３は、例えば、主軸モータ２の軸の回転位置を検出する。この回転位置から回転速度ｎが演算される。
The
振動センサ４は、加工中の工具１１またはワークＷに発生する振動を検出する。振動センサ４は、例えば、ワークＷに接触するように設けられる。振動センサ４の位置は、加工点に近いことが好ましい。尚、振動センサ４は、ワークＷが載置または固定される台に設けられてもよい。振動センサ４は、例えば、加速度センサである。しかし、これに限られず、振動センサ４は、主軸１またはワークＷを移動させる移動軸の位置（変位）を検出するリニアスケールなどであってもよい。また、振動センサ４は、例えば、主軸１に設けられてもよい。尚、以下では、工具１１に振動が発生する場合について説明する。しかし、これに限られず、工具１１に発生する振動は、ワークＷに発生する振動であってもよい。従って、ワークＷが薄肉工作物である場合等、ワークＷが振動しやすい場合であっても、工作機械１００を用いることができる。
The vibration sensor 4 detects the vibration generated in the
数値制御装置５は、主軸１の回転速度を制御する。より詳細には、数値制御装置５は、主軸１の回転速度ｎを周期的に変動させるように主軸１を制御する。回転速度ｎを周期的に変化させることは、以下では、ＳＳＶ（Spindle Speed Variation）とも呼ばれる場合がある。ＳＳＶは、加工中の工具１１またはワークＷに発生するびびり振動を抑制する方法の一つとして知られている。
The
びびり振動は、工具１１とワークＷとの間に生じる大きな振動である。びびり振動には、強制びびり振動と自励びびり振動とがある。強制びびり振動は、例えば、振動源により工具１１またはワークＷの振動が加振されて生じるびびり振動である。また、自励びびり振動の一つとして、再生型びびり振動が知られている。再生型びびり振動は、再生効果によって生じるびびり振動である。再生効果は、例えば、主軸１回転前（多刃工具では一刃前）の切削時に生じていた工具１１の振動が加工面の起伏として残り、その加工面を現在の工具１１が振動しながら削ることにより、切り取り厚さが変動することである（図２を参照）。この切り取り厚さにより切削力が変動し、工具１１が加振される。加振された工具１１の振動は、切削によって再び加工面に転写され、再生効果として主軸１回転後（多刃工具では１刃後）の切り取り厚さを変動させる。これが繰り返されることで、或る条件において振動が成長すると、大きな不安定振動（びびり振動）が発生する。切削力は、工具１１がワークＷを切削するために必要とする力である。以下では、再生型の自励びびり振動を、びびり振動と呼ぶ。尚、びびり振動の詳細については、図２を参照して、後で説明する。また、ＳＳＶの詳細については、図３を参照して、後で説明する。
The chatter vibration is a large vibration generated between the
数値制御装置５は、センサ情報取得部５１と、演算部５２と、を有する。
The
センサ情報取得部５１は、エンコーダ３および振動センサ４から、回転速度、および、振動の振動周波数を取得する。例えば、センサ情報取得部５１は、びびり振動を検知してもよい。この場合、センサ情報取得部５１は、振動センサ４の検出値を周波数解析する。センサ情報取得部５１は、例えば、振動の強度が所定の閾値を越えた場合、びびり振動が発生したと判断する。また、センサ情報取得部５１は、閾値を越えた強度の振動の周波数を、びびり振動の振動周波数ｆｃ（Ｈｚ）として取得する。センサ情報取得部５１は、びびり振動発生時の回転速度ｎを中心回転速度ｎ０として取得する。中心回転速度ｎ０は、ＳＳＶの変動の中心速度であり、ＳＳＶ開始時の回転速度でもある。センサ情報取得部５１は、中心回転速度ｎ０および振動周波数ｆｃを演算部５２に送る。
The sensor information acquisition unit 51 acquires the rotation speed and the vibration frequency of the vibration from the
演算部５２は、振動および振動のエネルギーをベッセル関数を用いて近似する。また、演算部５２は、ワークＷを加工する中心回転速度ｎ０（第１回転速度）を中心とする回転速度の変動幅および変動周波数を演算し、変動幅および変動周波数を中心回転速度ｎ０に重畳させてワークＷを加工する第２回転速度ｎ（ｔ）を演算する。例えば、演算部５２は、回転速度、振動周波数、およびベッセル関数で表される振動のエネルギーに基づいて、中心回転速度ｎ０を中心とする回転速度ｎの変動幅を演算する。振動のエネルギーは、びびり振動の要因となるエネルギーであり、例えば、切削力（加工力）に基づく加工プロセスのエネルギーである。変動幅は、ＳＳＶのパラメータであり、例えば、中心回転速度ｎ０に対する速度変動率である。変動幅は、以下では、ＲＶＡ（Ratio of the Speed Variation Amplitude）とも呼ばれる場合がある。より詳細には、演算部５２は、ベッセル関数の極小点に基づいて、変動幅および変動周波数を演算する。これにより、振動のエネルギーが低く、びびり振動を抑制する適切なＲＶＡを選定することができる。
The
また、演算部５２は、回転速度ｎの変動周波数を演算する。変動周波数は、ＲＶＡと同様に、ＳＳＶのパラメータである。変動周波数は、例えば、中心回転速度ｎ０に対する速度変動周波数比率である。変動周波数は、以下では、ＲＶＦ（Ratio of the Speed Variation Frequency）とも呼ばれる場合がある。また、例えば、演算部５２は、演算した変動幅および変動周波数でＳＳＶを実行するように、サーボ制御指令をサーボ制御部６に送る。従って、数値制御装置５は、主軸１を第２回転速度ｎ（ｔ）で制御してワークＷを加工する。尚、ＲＶＡおよびＲＶＦの演算方法の詳細については、図４を参照して、後で説明する。
Further, the
サーボ制御部６は、演算部５２からサーボ制御指令を受けて、主軸モータ２を制御する。これにより、ＳＳＶが実行される。
The servo control unit 6 receives a servo control command from the
尚、数値制御装置５を構成するセンサ情報取得部５１と、演算部５２とは、一つのＣＰＵ（Central Processing Unit）で実現されてもよく、それぞれ個別のＣＰＵで実現されていてもよい。また、サーボ制御部６は、数値制御装置５の一部であってもよい。
The sensor information acquisition unit 51 and the
次に、図２および図３を参照して、びびり振動およびＳＳＶについて説明する。 Next, chatter vibration and SSV will be described with reference to FIGS. 2 and 3.
図２は、第１実施形態による加工プロセスのモデルの一例を示す図である。図２は、工具１１がワークＷを切削する場合の例を示す。切削速度の方向は、図２の右方向である。破線ｘ（ｔ）は、加工中の現在の工具１１の加工経路を示す。実線ｘ（ｔ−τ）は、前回の加工時の工具１１の加工経路を示す。前回の切削は、例えば、主軸１回転前（多刃工具の場合、１つの刃前）を示す。また、破線ｘ（ｔ）および実線ｘ（ｔ−τ）は、工具１１の振動（変位）を示し、この振動が加工によりワークＷの加工面に転写される。尚、ｘ（ｔ）およびｘ（ｔ−τ）の軸をそれぞれＸ１軸およびＸ２軸とし、ワークＷから離れる方向をＸ１軸およびＸ２軸の正とする。また、ｘ（ｔ）およびｘ（ｔ−τ）の原点をそれぞれＯ１およびＯ２とする。ｔは、時間を示す。τは、工具１１の刃による連続する２回の加工の遅延時間を示す。すなわち、遅延時間τは、工具１１の刃が次に加工するまでの時間である。ε（ｒａｄ）は、現在の加工の変位ｘ（ｔ）と前回の加工の変位ｘ（ｔ−τ）との位相差である。ｈ（ｔ）は、工具１１による切り取り厚さを示す。Ｆは、反作用によって工具１１に作用する切削力を示す。
FIG. 2 is a diagram showing an example of a model of a machining process according to the first embodiment. FIG. 2 shows an example in which the
また、ＷｐｒｏｃｅｓｓおよびＷｍａｃｈｉｎｅは、切削力Ｆが工具１１に対してする仕事を示す。Ｗｐｒｏｃｅｓｓは、例えば、切削力Ｆとは反対方向に工具１１が変位（振動）する場合の仕事を示す。Ｗｍａｃｈｉｎｅは、例えば、切削力Ｆと同じ方向に工具１１が変位（振動）する場合の仕事を示す。
Further, W process and W machine indicate the work that the cutting force F does to the
切り取り厚さｈ（ｔ）は、加工条件として設定されるノミナル切り取り厚さ（設定切り取り厚さ）ｈ０に、変位ｘ（ｔ）と変位ｘ（ｔ−τ）との差を加えた厚さである。図２に示す例では、前回の加工における工具１１の振動が転写された加工面を、振動する工具１１が加工するため、切り取り厚さｈ（ｔ）が時間ｔによって変化する。
The cutting thickness h (t) is the thickness obtained by adding the difference between the displacement x (t) and the displacement x (t−τ) to the nominal cutting thickness (set cutting thickness) h 0 set as the processing condition. Is. In the example shown in FIG. 2, since the vibrating
仕事Ｗｍａｃｈｉｎｅでは、切削力Ｆの方向と工具１１の変位（ｘ成分）方向が平行であるため、工具１１は、切削力Ｆから正の仕事を受ける。従って、正のエネルギーが工具１１に流入する。一方、仕事Ｗｐｒｏｃｅｓｓでは、切削力Ｆの方向と工具１１の変位（ｘ成分）方向が反対であるため、工具１１は、切削力Ｆから負の仕事を受ける。従って、負のエネルギーが工具１１に流入する。すなわち、工具１１における振動のエネルギーが消散（減少）する。工具１１は、加工中に仕事Ｗｍａｃｈｉｎｅおよび仕事Ｗｐｒｏｃｅｓｓのエネルギーを交互に受け続ける。
In the work W machine , since the direction of the cutting force F and the displacement (x component) direction of the
ε＝０、２πの場合、変位ｘ（ｔ）と変位ｘ（ｔ−τ）との位相が略一致する。この場合、仕事Ｗｐｒｏｃｅｓｓと仕事Ｗｍａｃｈｉｎｅとがほぼ等しくなる。従って、工具１１が受けるエネルギーと消散されるエネルギーとがほぼ等しい。従って、振動のエネルギーはほぼ変化しないため、振動はほとんど発達も減衰もしない。従って、切り取り厚さがほとんど変動しないため、切削力Ｆの変動による機械（工具１１）の加振が行われない。この結果、びびり振動は発生しない。
When ε = 0 and 2π, the phases of the displacement x (t) and the displacement x (t−τ) are substantially the same. In this case, the work W process and the work W machine are almost equal. Therefore, the energy received by the
０＜ε≦πの場合、仕事Ｗｐｒｏｃｅｓｓの切り取り厚さは、仕事Ｗｍａｃｈｉｎｅの切り取り厚さよりも大きい。この場合、消散されるエネルギーが、工具１１に流入するエネルギーよりも大きい。従って、切り取り厚さの変動によって工具１１が加振されるものの、エネルギー的に安定であるため、振動は減衰され発達しない。この結果、びびり振動は発生しない。また、すでにびびり振動が発生している場合、位相差εが０＜ε≦πになるとびびり振動は抑制される。
When 0 < ε≤π , the cut-out thickness of the work W process is larger than the cut-out thickness of the work W machine. In this case, the dissipated energy is larger than the energy flowing into the
π＜ε＜２πの場合、仕事Ｗｍａｃｈｉｎｅの切り取り厚さは、仕事Ｗｐｒｏｃｅｓｓの切り取り厚さよりも大きい。この場合、工具１１に流入するエネルギーが、消散されるエネルギーよりも大きい。従って、図２に示される振動１サイクル毎に工具１１にエネルギーが流入する。この振動１サイクルあたりに流入するエネルギーが、機械が元来持つ機械減衰能力によって１サイクル中に吸収できるエネルギー容量を超えると、加工の継続により工具１１における振動のエネルギーは増大し続ける。この結果、工具１１はエネルギー的に不安定になる。この増大したエネルギーを消散させるため工具１１が発振することにより、びびり振動が発生する。
When π <ε <2π, the cut-out thickness of the work W machine is larger than the cut-out thickness of the work W process. In this case, the energy flowing into the
また、切削幅（図２の紙面に垂直な方向の幅）の設定が大きいほど、π＜ε＜２πにおいて１サイクル毎に工具１１に流入するエネルギー量が増大する。そのため、切削幅が或る閾値を超えるとびびり振動が発生する。尚、びびり振動の発生を避けるために切削幅の設定を下げると、工作機械１００の生産性が下がってしまう。
Further, the larger the setting of the cutting width (the width in the direction perpendicular to the paper surface of FIG. 2), the larger the amount of energy flowing into the
切削幅を下げることなくびびり振動を抑制する方法の一つとして、ＳＳＶが知られている。ＳＳＶでは、加工中に位相差εを変化させることにより、π＜ε＜２πにおける振動のエネルギーの増大を抑制する。 SSV is known as one of the methods for suppressing chatter vibration without reducing the cutting width. In SSV, the increase in vibration energy when π <ε <2π is suppressed by changing the phase difference ε during machining.
図３（Ａ）および図３（Ｂ）は、第１実施形態によるＳＳＶの一例を示す図である。図３（Ａ）の縦軸は、回転速度ｎ（ｍｉｎ−１）を示し、横軸は、時間ｔを示す。図３（Ｂ）の縦軸は、位相差εを示し、横軸は、時間ｔを示す。ｎ０は、中心回転速度である。ε０は、ＳＳＶ開始時の初期位相差である。αｓは、中心回転速度ｎ０に対するＳＳＶの変動振幅率である。ｆｓ（Ｈｚ）は、ＳＳＶの周波数である。尚、図３（Ａ）および図３（Ｂ）の横軸は、それぞれ対応しており、ＳＳＶの１周期１／ｆｓ（ｓ）の範囲の時間が示される。また、図２において説明したように、π＜ε＜２πの領域は、びびり振動が発生する不安定領域である。また、０＜ε≦πの領域は、びびり振動が抑制される安定領域である。 3 (A) and 3 (B) are diagrams showing an example of SSV according to the first embodiment. The vertical axis of FIG. 3A shows the rotation speed n (min -1 ), and the horizontal axis shows the time t. The vertical axis of FIG. 3B shows the phase difference ε, and the horizontal axis shows the time t. n 0 is the central rotation speed. ε 0 is the initial phase difference at the start of SSV. α s is the fluctuation amplitude rate of SSV with respect to the central rotation speed n 0. f s (Hz) is the frequency of SSV. The horizontal axes of FIGS. 3 (A) and 3 (B) correspond to each other, and the time in the range of 1 cycle 1 / f s (s) of SSV is shown. Further, as described with reference to FIG. 2, the region of π <ε <2π is an unstable region in which chatter vibration occurs. Further, the region of 0 <ε ≦ π is a stable region in which chatter vibration is suppressed.
図３（Ａ）に示すように、ＳＳＶにおける主軸１の回転速度ｎ（ｔ）は、時間ｔの関数として、式１により表される。
ｎ（ｔ）＝ｎ０（１＋αｓｓｉｎ２πｆｓｔ） （式１）
尚、ＳＳＶの一例として、回転速度ｎ（ｔ）が正弦波で変動する例が示されている。また、式１に示す回転速度ｎ（ｔ）は、ｎ（ｔ）＝ｎ０（１−αｓｓｉｎ２πｆｓｔ）やｎ（ｔ）＝ｎ０（１＋αｓｃｏｓ２πｆｓｔ）、ｎ（ｔ）＝ｎ０（１−αｓｃｏｓ２πｆｓｔ）等であってもよい。また、図３（Ａ）に示すように、回転速度ｎの変動振幅は２αｓｎ０であり、回転速度ｎの変動周期は１／ｆｓである。
As shown in FIG. 3A, the rotation speed n (t) of the
n (t) = n 0 (1 + α s sin2πf s t) (Equation 1)
As an example of SSV, an example in which the rotation speed n (t) fluctuates with a sine wave is shown. The rotational speed n shown in equation 1 (t) is, n (t) = n 0 (1-α s sin2πf s t) and n (t) = n 0 ( 1 + α s cos2πf s t), n (t) = N 0 (1-α s cos 2πf s t) and the like may be used. Further, as shown in FIG. 3A, the fluctuation amplitude of the rotation speed n is 2α s n 0 , and the fluctuation period of the rotation speed n is 1 / f s .
速度変動率であるＲＶＡは、変動振幅率αｓを用いて、式２により表される。
ＲＶＡ＝αｓ （式２）
また、速度変動周波数比率であるＲＶＦは、ＳＳＶ周波数ｆｓ、単位時間あたりの工具１１の刃による加工回数である中心切れ刃通過周波数ｆｚ０を用いて、式３により表される。
ＲＶＦ＝ｆｓ／ｆｚ０ （式３）
尚、「中心」は、ＳＳＶの変動の中心を示し、ＳＳＶの開始時の値でもあることも示す。また、中心切れ刃通過周波数ｆｚ０（Ｈｚ）は、工具１１の刃数Ｚおよび中心回転速度ｎ０（ｍｉｎ−１）により表される。従って、ＲＶＦは、式４により表される。
ＲＶＦ＝６０ｆｓ／（Ｚ×ｎ０） （式４）
RVA, which is the velocity volatility, is expressed by
RVA = α s (expression 2)
The RVF, which is the speed fluctuation frequency ratio, is expressed by
RVF = f s / f z0 (Equation 3)
The "center" indicates the center of fluctuation of SSV, and also indicates the value at the start of SSV. The center cutting edge passing frequency f z0 (Hz) is represented by the number of blades Z of the
RVF = 60f s / (Z × n 0) ( Equation 4)
図３（Ｂ）に示すように、回転速度ｎの変動に対応して、位相差εも変動する。ＳＳＶの開始時においてびびり振動が発生しているため、初期位相差ε０は、π＜ε０＜２πの範囲内である。回転速度ｎが中心回転速度ｎ０から上昇すると、主軸１が１回転（多刃工具の場合、１つの刃の通過）する時間が短くなる。従って、位相差εは、小さくなり、不安定領域から安定領域に移動する。その後、回転速度ｎが下降すると、主軸１が１回転する時間が長くなる。従って、位相差εは、大きくなり、安定領域から不安定領域に移動する。このように、位相差εが安定領域と不安定領域とを繰り返し通過し、ＳＳＶの１周期全体として振動のエネルギーを減少させることができる。これにより、びびり振動を抑制することができる。
As shown in FIG. 3B, the phase difference ε also fluctuates in response to the fluctuation of the rotation speed n. Since chatter vibration is generated at the start of SSV, the initial phase difference ε 0 is within the range of π <ε 0 <2π. When the rotation speed n rises from the central rotation speed n 0 , the time for the
変動幅は、工具１１の刃による連続する２回の加工における振動の位相差εが、変動する回転速度ｎの最大値である最大回転速度ｎｍにおいて、安定領域内であることを満たす。安定領域は、振動が抑制される領域であり、位相差εが０からπまでの領域である。最大回転速度ｎｍは、図３（Ａ）に示すように、式５により表される。
ｎｍ＝ｎ０（１＋αｓ） （式５）
図３（Ｂ）に示すように、最大回転速度ｎｍ付近では、回転速度ｎおよび位相差εの時間変化が小さくなる。従って、一定の回転速度ｎで加工が行われている状態とほぼ同じであり、ＳＳＶによるびびり振動の抑制が弱くなってしまう。すなわち、演算部５２は、びびり振動が発生しやすくなる最大回転速度ｎｍにおける最大速度位相差εｍが安定領域内となるように、変動振幅率αｓを選択する。これにより、最大回転速度ｎｍ付近での振動の発達を抑え、より効果的にびびり振動を抑制することができる。
The fluctuation width satisfies that the phase difference ε of the vibration in the two consecutive machining by the blade of the
n m = n 0 (1 + α s ) (expression 5)
As shown in FIG. 3B, the time change of the rotation speed n and the phase difference ε becomes small in the vicinity of the maximum rotation speed n m. Therefore, it is almost the same as the state in which processing is performed at a constant rotation speed n, and the suppression of chatter vibration by SSV is weakened. That is, the
より詳細には、位相差εが最大回転速度ｎｍにおいて高安定領域内であることが好ましい。高安定領域は、安定領域よりも振動が抑制される領域であり、位相差εが０からπ／２までの領域である。振動のエネルギーは、位相差εが安定領域内に入ってから消散され始める。また、後で説明するように、位相差εがπ／２である場合、振動のエネルギーが最も消散される。従って、高安定領域（０＜ε＜π／２）では、π／２＜ε＜πの領域よりも、振動のエネルギーが消散されている。この結果、びびり振動をより抑制することができる。 More specifically, it is preferable that the phase difference ε is within the highly stable region at the maximum rotation speed n m. The highly stable region is a region in which vibration is suppressed more than the stable region, and is a region in which the phase difference ε is from 0 to π / 2. The vibrational energy begins to dissipate after the phase difference ε enters the stable region. Further, as will be described later, when the phase difference ε is π / 2, the vibration energy is most dissipated. Therefore, in the highly stable region (0 <ε <π / 2), the vibration energy is dissipated more than in the region of π / 2 <ε <π. As a result, chatter vibration can be further suppressed.
また、図３（Ｂ）において説明したように、ＳＳＶでは、位相差εが安定領域と不安定領域とを繰り返し通過する。従って、びびり振動の抑制が弱くなる不安定領域においてびびり振動が発生してしまう場合がある。この場合、ＳＳＶ周波数ｆｓに対応する過渡振動（beat vibration）が発生する。そこで、びびり振動を適切に抑制するためには、ＲＶＡやＲＶＦのようなＳＳＶのパラメータを適切に選定する必要がある。 Further, as described in FIG. 3B, in SSV, the phase difference ε repeatedly passes through the stable region and the unstable region. Therefore, chatter vibration may occur in an unstable region where the suppression of chatter vibration is weakened. In this case, the transient vibration corresponding to SSV frequency f s (beat vibration) occurs. Therefore, in order to appropriately suppress chatter vibration, it is necessary to appropriately select SSV parameters such as RVA and RVF.
次に、図４を参照して、ＲＶＡおよびＲＶＦの演算方法（振動抑制方法）について説明する。 Next, the calculation method (vibration suppression method) of RVA and RVF will be described with reference to FIG.
図４は、第１実施形態による工作機械１００の動作例を示すフロー図である。尚、図４のフロー図には、変動振幅率αｓ（ＲＶＡ）の演算（Ｓ１０〜Ｓ３０）と、複数の制約によるＲＶＡの抽出（Ｓ４０〜Ｓ２００）と、ＲＶＦの選択（Ｓ２２０）とが含まれる。
FIG. 4 is a flow chart showing an operation example of the
まず、センサ情報取得部５１は、エンコーダ３および振動センサ４から、中心回転速度ｎ０および振動周波数ｆｃを取得する（Ｓ１０）。例えば、センサ情報取得部５１は、ワークＷの加工中にびびり振動が発生したか否かを判定する。センサ情報取得部５１は、例えば、センサ情報（振動センサ４の検出値）をＦＦＴ（Fast Fourier Transformation）演算する。センサ情報取得部５１は、或る周波数での周波数成分の大きさが閾値を超えた場合に、びびり振動が発生したと判断する。びびり振動が発生した場合、センサ情報取得部５１は、中心回転速度ｎ０および振動周波数ｆｃを取得する。
First, the sensor information acquisition unit 51, the
次に、演算部５２は、中心振動波数ｋ０および中心振動位相差ε０を演算する（Ｓ２０）。中心振動波数ｋ０は、工具１１の刃による連続する２回の加工の間におけるびびり振動の波数を示す。中心振動波数ｋ０は、例えば、主軸１回転に対するびびり振動の波数である。中心振動位相差ε０（ｒａｄ）は、工具１１の刃による連続する２回の加工に対するびびり振動の位相差である。例えば、中心振動位相差ε０は、主軸１回転の加工距離がびびり振動の波長の整数倍からずれた場合に生じる小数部である。中心振動波数ｋ０および中心振動位相差ε０は、中心切れ刃通過周波数ｆｚ０を用いて、式６により表される。
ｆｃ／ｆｚ０＝ｋ０＋ε０／２π （式６）
Next, the
f c / f z0 = k 0 + ε 0 / 2π (Equation 6)
また、中心切れ刃通過周波数ｆｚ０（Ｈｚ）は、刃数Ｚおよび中心回転速度ｎ０（ｍｉｎ−１）により表される。従って、中心振動波数ｋ０および中心振動位相差ε０は、式７により表される。
６０ｆｃ／（Ｚ×ｎ０）＝ｋ０＋ε０／２π （式７）
式７の右辺の第１項であるｋ０は、整数の商であり、第２項であるε０／２πは、余りである。尚、中心振動位相差ε０は、図３（Ｂ）に示す初期位相差ε０に対応する。
The center cutting edge passing frequency f z0 (Hz) is represented by the number of blades Z and the center rotation speed n 0 (min -1). Therefore, the central vibration wave number k 0 and the center vibration phase difference ε 0 are expressed by Equation 7.
60f c / (Z × n 0 ) = k 0 + ε 0 / 2π (Equation 7)
The first term, k 0 , on the right-hand side of Equation 7 is the quotient of an integer, and the second term, ε 0 / 2π, is the remainder. The central vibration phase difference ε 0 corresponds to the initial phase difference ε 0 shown in FIG. 3 (B).
次に、演算部５２は、ベッセル関数Ｊ０（ｍｆ）が極小値をとるｍｆを選び、変動振幅率αｓを演算する（Ｓ３０）。振動のエネルギーは、ベッセル関数を用いて表される。従って、振動のエネルギーが安定となる変動振幅率αｓを、ベッセル関数を用いて選定することができる。例えば、演算部５２は、複数のｍｆに対応する複数の変動振幅率αｓを演算し、ＲＶＡの候補とする。尚、ベッセル関数による変動振幅率αｓの演算の詳細については、図５（Ａ）および図５（Ｂ）を参照して、後で説明する。
Next, the calculation unit 52 selects m f for which the Bessel function J 0 (m f ) has a minimum value, and calculates the fluctuation amplitude rate α s (S30). The energy of vibration is expressed using the Bessel function. Therefore, the fluctuation amplitude rate α s at which the vibration energy is stable can be selected by using the Bessel function. For example, the
次に、演算部５２は、αｓ≦αｍａｘおよび許容切削速度の範囲を満たすＲＶＡの候補を抽出する（Ｓ４０）。αｍａｘは、機械的制約による変動振幅率αｓの上限である。αｍａｘは、例えば、サーボ制御部６や主軸モータ２の性能によって予め設定される。αｍａｘは、例えば、２０％である。この場合、中心回転速度ｎ０（例えば、１０００ｍｉｎ−１）から、中心回転速度ｎ０の２０％（例えば、２００ｍｉｎ−１）まで主軸１の回転速度ｎを変動させることができる。また、許容切削速度は、ワークＷや工具１１ごとに予め設定される、回転速度ｎの許容値である。
Next, the
次に、演算部５２は、最大回転速度ｎｍにおける最大速度波数ｋｍおよび最大速度位相差εｍを演算する（Ｓ５０）。最大速度波数ｋｍおよび最大速度位相差εｍは、式７と同様に、式８により表される。
６０ｆｃ／（Ｚ×ｎｍ）＝ｋｍ＋εｍ／２π （式８）
Next, the
60f c / (Z × n m ) = k m + ε m / 2π ( Equation 8)
次に、演算部５２は、位相差εの範囲εｍｂを閉区間[０，π／２]とし、変数ｆｌａｇにゼロを代入する（Ｓ６０）。変数ｆｌａｇのゼロは、範囲εｍｂが高安定領域であることを示す。
Next, the
次に、演算部５２は、最大速度位相差εｍが範囲εｍｂ内にあることを満たすＲＶＡの候補を抽出する（Ｓ７０）。この場合、最大速度位相差εｍが高安定領域内にあることを満たすＲＶＡが抽出される。
Next, the
次に、演算部５２は、変数ｆｌａｇが１であるか否かを判定する（Ｓ８０）。変数ｆｌａｇが１ではない場合（Ｓ８０のＮＯ）、演算部５２は、ＲＶＡの候補が存在するか否かを判定する（Ｓ９０）。ＲＶＡの候補が存在しない場合（Ｓ９０のＮＯ）、演算部５２は、範囲εｍｂを閉区間[０，π]とし（Ｓ１００）、変数ｆｌａｇに１を代入する（Ｓ１１０）。変数ｆｌａｇの１は、範囲εｍｂが安定領域であることを示す。すなわち、演算部５２は、最大速度位相差εｍに関する制約を高安定領域内から安定領域内に緩和する。一方、変数ｆｌａｇが１である場合（Ｓ８０のＹＥＳ）、ステップＳ１２０が実行される。
Next, the
ＲＶＡの候補が存在する場合（Ｓ９０のＹＥＳ）、演算部５２は、ＳＳＶ効率に対するＲＶＦの上限値を演算する（Ｓ１２０）。一般に、ＳＳＶ周波数ｆｓが高くなるほど、びびり振動を強く抑制することができる。しかし、ＳＳＶ周波数ｆｓが高くなりすぎると、びびり振動の抑制が弱くなり、第１実施形態によるＳＳＶの効率が減少する場合がある。従って、所定の効率を得るため、ＳＳＶ周波数ｆｓ（ＲＶＦ）の上限値が設定される。以下では、一例として、９８％の効率の上限値ＲＶＦ９８％が選定された場合について説明する。尚、ＲＶＦの上限値の演算の詳細については、図６（Ａ）〜図７を参照して、後で説明する。
When there is a candidate for RVA (YES in S90), the
次に、演算部５２は、ＳＳＶ周波数ｆｓの上限値ｆ９８％を演算する（Ｓ１３０）。ｆ９８％は、ＲＶＦ９８％に対応するＳＳＶ周波数ｆｓの上限値である。ｆ９８％は、式４を用いて、式９により表される。
ｆ９８％＝ＲＶＦ９８％×（（Ｚ×ｎ０）／６０） （式９）
Next, the
f 98% = RVF 98% × ((Z × n 0 ) / 60) (Equation 9)
次に、演算部５２は、安定制約レベルＮに４を代入する（Ｓ１４０）。安定制約レベルＮは、ＳＳＶ周波数ｆｓの下限値の制約の強さを示す。安定制約レベルＮは、弱い順から強い順に１から４まで存在する。安定制約レベルＮが強いほど、ＳＳＶ周波数ｆｓの下限値ｆｓｌｉｍＮは大きくなり、びびり振動を強く抑制することができるＳＳＶ周波数ｆｓを選定することができる。
Next, the
次に、演算部５２は、安定制約レベルＮに対応するＳＳＶ周波数ｆｓの下限値ｆｓｌｉｍＮを演算する（Ｓ１５０）。尚、下限値ｆｓｌｉｍＮの演算の詳細については、図８（Ａ）および図８（Ｂ）を参照して、後で説明する。
Next, the
次に、演算部５２は、ｆｓｌｉｍＮ≦ｆ９８％∩ｆｓｌｉｍＮ≦ｆｂａｎｄを満たすＲＶＡの候補を抽出する（Ｓ１６０）。ｆｓｌｉｍＮ≦ｆ９８％は、下限値ｆｓｌｉｍＮが上限値ｆ９８％以下であることを示す。ｆｓｌｉｍＮ≦ｆｂａｎｄは、ｆｓｌｉｍＮが機械的制約である主軸システムの帯域幅ｆｂａｎｄ以内であることを示す。帯域幅ｆｂａｎｄは、回転速度の変動を制御することができる周波数を示し、例えば、サーボ制御部６の性能によって予め設定される。
Next, the
次に、演算部５２は、ＲＶＡの候補が存在するか否かを判定する（Ｓ１７０）。ＲＶＡの候補が存在しない場合（Ｓ１７０のＮＯ）、安定制約レベルＮが１であるか否かを判定する（Ｓ１８０）。安定制約レベルＮが１である場合（Ｓ１８０のＹＥＳ）、演算部５２は、びびり振動の他の抑制方法を推奨する（Ｓ２１０）。これは、適切なＲＶＡの候補を求めることができず、発生したびびり振動をＳＳＶで抑制することが困難なためである。安定制約レベルＮが１ではない場合（Ｓ１８０のＮＯ）、演算部５２は、変数ｆｌａｇがゼロであるか否かを判定する（Ｓ１９０）。変数ｆｌａｇがゼロである場合（Ｓ１９０のＹＥＳ）、演算部５２は、範囲εｍｂを閉区間[０，π]とし（Ｓ１００）、変数ｆｌａｇに１を代入する（Ｓ１１０）。すなわち、再びステップＳ７０〜Ｓ１７０が実行され、最大速度位相差εｍに関する制約を高安定領域内から安定領域内に緩和してＲＶＡの候補が抽出される。変数ｆｌａｇがゼロでない場合（Ｓ１９０のＮＯ）、演算部５２は、安定制約レベルＮにＮ−１を代入する（Ｓ２００）。すなわち、演算部５２は、安定制約レベルＮを下げ（下限値ｆｓｌｉｍＮの制約を緩和し）、再びステップＳ１５０〜Ｓ１７０を実行する。
Next, the
演算される下限値ｆｓｌｉｍＮの大小関係は、式１０により表される。
ｆｓｌｉｍ１＜ｆｓｌｉｍ２＜ｆｓｌｉｍ３＜ｆｓｌｉｍ４ （式１０）
安定制約レベルＮが大きくなるほど下限値ｆｓｌｉｍＮが大きくなり、びびり振動を抑制しやすくなる。従って、演算部５２は、下限値ｆｓｌｉｍＮを下げて制約を緩和しながら、条件に合うＲＶＡの候補を抽出する。
The magnitude relation of the lower limit value f slimN to be calculated is expressed by Equation 10.
f slim1 <f slim2 <f slim3 <f slim4 (Equation 10)
As the stability constraint level N becomes larger, the lower limit value f slimN becomes larger, and it becomes easier to suppress chatter vibration. Therefore, the
一方、ＲＶＡの候補が存在する場合（Ｓ１７０のＹＥＳ）、演算部５２は、最終候補からＲＶＡを選択し、ｆｓｌｉｍＮ≦ｆｓ≦ｆ９８％∩ｆｓｌｉｍＮ≦ｆｓ≦ｆｂａｎｄおよび許容主軸モータ負荷を満たすＲＶＦを選択する（Ｓ２２０）。許容主軸モータ負荷は、予め設定される。演算部５２は、選択したＲＶＡおよびＲＶＦをパラメータとするＳＳＶを実行するように、サーボ制御指令をサーボ制御部６に送る。また、びびり振動の抑制効果が現れるのを待つため、数値制御装置５は、例えば、少なくとも主軸２回転分の時間以上待つ。
On the other hand, when there is a candidate for RVA (YES in S170), the
次に、センサ情報取得部５１は、びびり振動が抑制されたか否かを判定する（Ｓ２３０）。センサ情報取得部５１は、びびり振動の検知と同様の判定をすればよい。びびり振動が抑制されない場合（Ｓ２３０のＮＯ）、演算部５２は、切削幅を下げることを提案する（Ｓ２４０）。これは、実現可能な最も効果的なＳＳＶパラメータを用いたとしても、びびり振動の強度が過大でびびり振動を抑制しきれないためである。すなわち、切削幅の設定を下げることにより、びびり振動を弱くすればよい。一方、びびり振動が抑制された場合（Ｓ２３０のＹＥＳ）、工作機械１００は、加工が終了するまで、その条件で振動の抑制を継続する。すなわち、ステップＳ２２０で選択されたＲＶＡおよびＲＶＦにより、ＳＳＶが継続して実行される。
Next, the sensor information acquisition unit 51 determines whether or not the chatter vibration is suppressed (S230). The sensor information acquisition unit 51 may make a determination similar to the detection of chatter vibration. When the chatter vibration is not suppressed (NO in S230), the
このように、工作機械１００は、図４に示すフローに従って、適切なＲＶＡおよびＲＶＦを自動で選定する。従って、ワークＷまたは工具１１にびびり振動が発生した場合、適切にびびり振動を抑制することができる。
In this way, the
（ＲＶＡの演算）
以下では、図２に示すように、加工プロセスがモデル化される。まず、振動方向の比切削抵抗をＫｃとし、切削幅をｂとし、切り取り厚さをｈ（ｔ）とし、切削力が切削面積に比例すると仮定すると、切削力Ｆｃ（ｔ）は、式１１により表される。
Ｆｃ（ｔ）＝Ｋｃｂｈ（ｔ） （式１１）
また、切り取り厚さｈ（ｔ）は、式１２により表される。
ｈ（ｔ）＝ｈ０−ｘ（ｔ）＋ｘ（ｔ−τ（ｔ）） （式１２）
或る単一の周波数ωｃ（ｒａｄ／ｓ）において、支配的な大きなびびり振動が発生すると仮定すると、工具１１の変位は、正弦波を用いたモデル化により、式１３および式１４により表される。
ｘ（ｔ）＝Ｘｃｏｓ（ωｃｔ） （式１３）
ｘ（ｔ−τ（ｔ））＝Ｘｃｏｓ（ωｃｔ−ωｃτ（ｔ）） （式１４）
尚、Ｘは、仮定されたびびり振動の振幅である。また、ωｃは、振動周波数ｆｃ（２πｆｃ）に対応する。
(RVA calculation)
In the following, the machining process is modeled as shown in FIG. First, assuming that the specific cutting resistance in the vibration direction is K c , the cutting width is b, the cutting thickness is h (t), and the cutting force is proportional to the cutting area, the cutting force F c (t) is given by the equation. It is represented by 11.
F c (t) = K c bh (t) (Equation 11)
The cut thickness h (t) is represented by the
h (t) = h 0 −x (t) + x (t−τ (t)) (Equation 12)
Assuming that a dominant large chatter vibration occurs at a single frequency ω c (rad / s), the displacement of the
x (t) = Xcos (ω c t) ( Equation 13)
x (t-τ (t)) = Xcos (ω c t-ω c τ (t)) (Equation 14)
Note that X is the assumed amplitude of the vibration. Further, ω c corresponds to the vibration frequency f c (2π f c).
また、ＳＳＶ周波数ｆｓが主軸１の切れ刃通過周波数ｆｚ（ｔ）（回転速度ｎ（ｔ））に対して十分小さい場合、遅延時間τ（ｔ）は、式１５に示すように近似される。
τ（ｔ）≒１／ｆｚ（ｔ）＝６０／（Ｚ×ｎ（ｔ）） （式１５）
すなわち、遅延時間τ（ｔ）は、切れ刃通過周波数ｆｚ（ｔ）の逆数に近似される。式１５のτ（ｔ）をさらに近似し、式１４に代入することにより、変位ｘ（ｔ−τ（ｔ））は、式１６により表される。
τ (t) ≈ 1 / f z (t) = 60 / (Z × n (t)) (Equation 15)
That is, the delay time τ (t) is approximated to the reciprocal of the cutting edge passing frequency fz (t). The displacement x (t−τ (t)) is expressed by
切削力（加工プロセス）のエネルギーＵｃは、力と速度との積である仕事の時間積分により求められる。例えば、切削力のエネルギーＵｃは、式１７により表される。
ここで、式１６に示す変位ｘ（ｔ−τ（ｔ））は、ｍｆを引数とした０次の第１種ベッセル関数Ｊ０（ｍｆ）を用いて近似することができる。式１６におけるｓｉｎ（２πｆｓｔ）の係数をｍｆとする。ベッセル関数により近似された変位ｘ（ｔ−τ（ｔ））を式１７の右辺の第３項に代入すると、振動のエネルギーＵは、近似により式１８により表される。
図５（Ａ）および図５（Ｂ）は、ベッセル関数の一例を示す図である。図５（Ａ）は、０次の第１種ベッセル関数を示すグラフである。縦軸は、０次の第１種ベッセル関数Ｊ０（ｍｆ）を示し、横軸は、引数ｍｆを示す。図５（Ｂ）は、０次の第１種ベッセル関数Ｊ０（ｍｆ）が極小となる引数ｍｆを示す図である。
5 (A) and 5 (B) are diagrams showing an example of the Bessel function. FIG. 5A is a graph showing a 0th-
図５（Ａ）に示すように、ベッセル関数Ｊ０（ｍｆ）は、引数ｍｆがゼロから大きくなるにつれて、減少してゼロを中心に振動する振る舞いを示す。ベッセル関数Ｊ０（ｍｆ）が負になる引数ｍｆを選定することにより、振動のエネルギーＵを負にすることができる。 As shown in FIG. 5 (A), the Bessel function J 0 (m f ) shows the behavior of decreasing and oscillating around zero as the argument m f increases from zero. The energy U of vibration can be made negative by selecting the argument m f in which the Bessel function J 0 (m f ) is negative.
より詳細には、引数ｍｆは、ベッセル関数Ｊ０（ｍｆ）の極小点に基づいて選定されることが好ましい。これにより、振動のエネルギーＵがより低くなる引数ｍｆを選定することができる。 More specifically, the argument m f is preferably selected based on the minimum point of the Bessel function J 0 (m f). Thereby, the argument m f at which the vibration energy U becomes lower can be selected.
また、式１９を変動振幅率αｓについて整理すると、変動振幅率αｓは、式２０により表される。
（ＲＶＦの上限値の演算）
図６（Ａ）および図６（Ｂ）は、ＳＳＶによる遅延時間τの変動の一例を示す図である。図６（Ａ）および図６（Ｂ）は、ＲＶＦがそれぞれ５％および４０％である場合を示す。縦軸は、遅延時間を示し、横軸は、時間ｔを示す。τは、実際の遅延時間を示し、τ〜は、ＳＳＶ周波数ｆｓが切れ刃通過周波数ｆｚ（ｔ）（回転速度ｎ（ｔ））に対して十分小さいと近似された遅延時間を示す。尚、「τ〜」は、「τ」の上にチルダが付された文字を示す。Ａτは、実際の遅延時間τの時間変動幅を示す。Ａτ〜は、近似された遅延時間τ〜の時間変動幅を示す。尚、「τ〜」は、「τ」上にチルダが付された文字を示す。τ０は、ＳＳＶ開始時の初期の遅延時間を示す。
(Calculation of the upper limit of RVF)
6 (A) and 6 (B) are diagrams showing an example of the fluctuation of the delay time τ due to SSV. 6 (A) and 6 (B) show the case where the RVF is 5% and 40%, respectively. The vertical axis represents the delay time, and the horizontal axis represents the time t. τ indicates the actual delay time, and τ ~ indicates the delay time approximated that the SSV frequency f s is sufficiently smaller than the cutting edge passing frequency f z (t) (rotation speed n (t)). In addition, "τ ~ " indicates a character in which a tilde is added above "τ". A τ indicates the time fluctuation range of the actual delay time τ. A τ ~ indicates the time fluctuation range of the approximated delay time τ ~. In addition, " τ ~ " indicates a character in which a tilde is added on "τ". τ 0 indicates the initial delay time at the start of SSV.
図４において説明したように、演算部５２は、変動周波数の上限値と変動周波数の下限値との間の変動周波数を決定する。また、演算部５２は、工具１１の刃による連続する２回の加工の遅延時間の実際の時間変動幅Ａτと近似した時間変動幅Ａτ〜とから求めた振動の抑制効率に基づいて、変動周波数の上限値を設定する。より詳細には、演算部５２は、遅延時間τの時間変動幅Ａτと、切れ刃通過周波数ｆｚ（ｔ）の逆数の時間変動幅Ａτ〜との比率に基づいて、回転速度ｎの変動周波数の上限値を演算する。
As described with reference to FIG. 4, the
ＳＳＶでは、回転速度ｎの変動により、遅延時間τも変動する。実際の遅延時間τは、式２１により表される。
図６（Ａ）に示すように、ＲＶＦが５％と小さい場合、実際の遅延時間τと近似された遅延時間τ〜とは、ほぼ同じである。すなわち、時間変動幅Ａτおよび時間変動幅Ａτ〜は、ほぼ同じである。しかし、図６（Ｂ）に示すように、ＲＶＦが４０％と大きい場合、時間変動幅Ａτは、時間変動幅Ａτ〜に対して小さくなる。これは、ＲＶＦが大きくなるほど、振動モデルに用いたＳＳＶ周波数ｆｓが小さいとする近似が成り立たなくなるためである。従って、ＲＶＦが大きくなるほど、実際のＳＳＶによるびびり振動の抑制は、振動モデルに基づいた演算結果に対して弱くなってしまう。すなわち、ＳＳＶの効率ηが減少する。ＳＳＶ効率ηは、式２３により表される。
η＝Ａτ／Ａτ〜 （式２３）
そこで、所定のＳＳＶ効率ηを得ることができるように、ＲＶＦの上限が設定される。時間変動幅Ａτは、実際の遅延時間τと同様に、シミュレーション演算により、演算される。
As shown in FIG. 6A, when the RVF is as small as 5%, the delay time τ ~, which is approximated to the actual delay time τ, is almost the same. That is, the time fluctuation width A τ and the time fluctuation width A τ ~ are almost the same. However, as shown in FIG. 6B, when the RVF is as large as 40%, the time fluctuation width A τ becomes smaller than the time fluctuation width A τ ~. This, RVF increases, approximating to SSV frequency f s using the vibration model is small because no longer hold. Therefore, as the RVF becomes larger, the suppression of chatter vibration by the actual SSV becomes weaker than the calculation result based on the vibration model. That is, the efficiency η of SSV decreases. The SSV efficiency η is represented by the
η = A τ / A τ ~ (Equation 23)
Therefore, the upper limit of RVF is set so that a predetermined SSV efficiency η can be obtained. The time fluctuation width A τ is calculated by a simulation calculation in the same manner as the actual delay time τ.
図７は、ＲＶＡとＲＶＦの上限値との関係の一例を示すグラフである。縦軸は、所定のＳＳＶ効率ηを得ることができるＲＶＦの上限値を示し、横軸は、ＲＶＡを示す。ＲＶＦ９８％、ＲＶＦ９７％、ＲＶＦ９５％およびＲＶＦ９０％は、それぞれ、ＳＳＶ効率ηが９８％、９７％、９５％、９０％の場合におけるＲＶＦの上限値である。実線は、シミュレーション演算により演算されたＲＶＦの上限値を示す。また、図７には、シミュレーション結果をＲＶＡ（変動振幅率αｓ）でフィッティングした結果の式が示されている。フィッティングは、一例として、４次の多項式により行われている。 FIG. 7 is a graph showing an example of the relationship between the upper limit of RVA and RVF. The vertical axis represents the upper limit value of RVF capable of obtaining a predetermined SSV efficiency η, and the horizontal axis represents RVA. RVF 98% , RVF 97% , RVF 95% and RVF 90% are the upper limit values of RVF when the SSV efficiency η is 98%, 97%, 95% and 90%, respectively. The solid line shows the upper limit value of RVF calculated by the simulation calculation. Further, FIG. 7 shows an equation of the result of fitting the simulation result by RVA (fluctuation amplitude rate α s). The fitting is performed by a fourth-order polynomial as an example.
図７に示すように、許容するＳＳＶ効率ηを緩和するほど、ＲＶＦの上限値は上昇する。また、ＲＶＦの上限値は、ＲＶＡ（変動振幅率αｓ）により変化する。演算部５２は、ＲＶＡの候補を図７に示すフィッティング式に代入することにより、ＲＶＦの上限値を演算する。
As shown in FIG. 7, as the allowable SSV efficiency η is relaxed, the upper limit value of RVF increases. Further, the upper limit value of RVF changes depending on RVA (fluctuation amplitude rate α s ). The
（ＲＶＦの下限値の演算）
図８（Ａ）および図８（Ｂ）は、第１実施形態による安定制約レベルの一例を示す図。図８（Ａ）および図８（Ｂ）に示す回転速度ｎおよび位相差εのグラフは、図３（Ａ）および図３（Ｂ）と同様である。矢印Ｎ１〜Ｎ４は、それぞれ、安定制約レベルＮの１から４までの位相差εの変化を示す。
(Calculation of the lower limit of RVF)
8 (A) and 8 (B) are diagrams showing an example of the stability constraint level according to the first embodiment. The graphs of the rotation speed n and the phase difference ε shown in FIGS. 8 (A) and 8 (B) are the same as those in FIGS. 3 (A) and 3 (B). Arrows N1 to N4 indicate changes in the phase difference ε from 1 to 4 of the stability constraint level N, respectively.
演算部５２は、工具１１の刃による連続する２回の加工における振動の位相差εの時間変化、および、演算した変動幅に基づいて、回転速度ｎの変動周波数の下限値ｆｓｌｉｍＮを演算する。より詳細には、演算部５２は、切れ刃通過周波数ｆｚの逆数である切れ刃通過時間１／ｆｚにおける位相差εの変化に基づいて、下限値ｆｓｌｉｍＮを演算する。尚、切れ刃通過時間１／ｆｚは、主軸１回転（多刃工具の場合、１つの刃の通過）の時間である。びびり振動は、上記のように、前回の加工時における工具１１の振動によって発生する。すなわち、主軸１回転毎（多刃工具の場合、１刃通過毎）にびびり振動が成長する。現在の位相差εが不安定領域内に入ってしまうと、びびり振動（過渡振動）が発生してしまう可能性がある。この場合、びびり振動が発生する前に、位相差εが不安定領域を抜ける必要がある。また、一加工前の振動がびびり振動に影響するため、次の加工（切れ刃通過時間１／ｆｚ）までに位相差εが不安定領域を抜ける必要がある。
The calculation unit 52 calculates the lower limit value f slimN of the fluctuation frequency of the rotation speed n based on the time change of the phase difference ε of the vibration in the two consecutive machining by the blade of the
位相差εが不安定領域を通り抜けるための時間は、ＲＶＡおよびＲＶＦによって変化する。例えば、小さい変動振幅率αｓが選定された場合、不安定領域における位相差εの時間に対する傾斜が小さくなる。すなわち、位相差εが不安定領域を通り抜けるための時間が長くなる。この場合、ＳＳＶ周波数ｆｓを大きくする必要がある。これは、ＳＳＶ周期１／ｆｓが短くなり、不安定領域における位相差εの時間に対する傾斜を大きくすることができるためである。
The time for the phase difference ε to pass through the unstable region varies depending on RVA and RVF. For example, when a small fluctuation amplitude factor α s is selected, the slope of the phase difference ε in the unstable region with respect to time becomes small. That is, the time for the phase difference ε to pass through the unstable region becomes longer. In this case, it is necessary to increase the SSV frequency f s. This,
このように、ＳＳＶ周波数ｆｓの下限値を設定し、位相差εが不安定領域を通り抜けるための時間を調整することができる。尚、大きい変動振幅率αｓが選定された場合、不安定領域における位相差εの時間に対する傾斜が大きくなる。すなわち、位相差εが不安定領域内を通り抜けるための時間が短くなる。この場合、位相差εが不安定領域を抜ける時間に余裕があるため、ＳＳＶ周波数ｆｓの下限値は低くなる。 Thus, it is possible to set the lower limit of the SSV frequency f s, the phase difference ε to adjust the time for passing through the unstable area. When a large fluctuation amplitude factor α s is selected, the gradient of the phase difference ε in the unstable region with respect to time becomes large. That is, the time for the phase difference ε to pass through the unstable region is shortened. In this case, since the phase difference ε can afford the exit time unstable region, the lower limit of the SSV frequency f s is lower.
図４において説明したように、演算部５２は、振動の抑制効果が高い順に変動周波数の下限値を変更しながら変動周波数を決定する。振動の抑制効果が高い順は、式１０に示す大小関係における下限値ｆｓｌｉｍＮが大きい順である。例えば、演算部５２は、ＲＶＡの候補を用いて下限値ｆｓｌｉｍＮを逐次演算する。以下では、安定制約レベルＮごとのＲＶＦの下限値について説明する。
As described with reference to FIG. 4, the
まず、安定制約レベルＮ＝４について説明する。演算部５２は、第２回転速度ｎ（ｔ）の最大回転速度ｎｍの時点から振動の位相差εが２πの時点までの時間が、単位時間当たりの最大回転速度ｎｍと工具１１の刃数Ｚとの積の逆数となる様に、変動周波数の下限値ｆｓｌｉｍ４（第１下限値）を設定する。図３（Ｂ）において説明したように、最大回転速度ｎｍ付近では、びびり振動が発生しやすい。従って、最も強い安定制約レベルＮ＝４では、最大回転速度ｎｍ直後およびその後の不安定領域を、位相差εが切れ刃通過時間１／ｆｚ以内に通り抜けることを仮定する。尚、最大速度位相差εｍは、安定領域内にあるとする。ＳＳＶ周波数ｆｓが下限値ｆｓｌｉｍ４の制約を満たす場合、びびり振動が発生しやすい領域を全て抜けるため、びびり振動をより強く抑制することができる。
First, the stability constraint level N = 4 will be described. In the calculation unit 52, the time from the time when the maximum rotation speed n m of the second rotation speed n (t) to the time when the phase difference ε of the vibration is 2π is the maximum rotation speed n m per unit time and the blade of the
図８（Ａ）の矢印Ｎ４に示すように、例えば、切れ刃通過時間１／ｆｚ（６０／（Ｚ×ｎｍ））後に回転速度ｎがｎｍからｎｍ−Δｎ２になる場合を仮定する。 As indicated by an arrow N4 in FIG. 8 (A), the example, a case where the rotation speed n after cutting edge transit time 1 / f z (60 / ( Z × n m)) is made of n m to n m -Δn 2 Assume.
最大回転速度ｎｍから６０／（Ｚ×ｎｍ）後の回転速度ｎｍ−Δｎ２は、式１を用いて、式２４により表される。
ｎｍ−Δｎ２＝ｎ０（１＋αｓｓｉｎ（２πｆｓ×６０／（Ｚ×ｎｍ）＋π／２）） （式２４）
また、回転速度ｎがｎｍ−Δｎ２である時の瞬間位相差は、２πである。従って、Δｎ２は、式８と同様に、式２５により表される。
６０ｆｃ／（Ｚ×（ｎｍ−Δｎ２））＝ｋｍ＋２π／２π （式２５）
The rotation speed n m − Δn 2 after 60 / (Z × nm ) from the maximum rotation speed n m is expressed by the formula 24 using the
n m − Δ n 2 = n 0 (1 + α s sin (2πf s × 60 / (Z × nm ) + π / 2)) (Equation 24)
Further, the instantaneous phase difference when the rotation speed n is nm − Δn 2 is 2π. Therefore, Δn 2 is represented by
60f c / (Z × (n m -Δn 2)) = k m + 2π / 2π ( Formula 25)
式２４をＳＳＶ周波数ｆｓについて整理し、式２５を代入してΔｎ２を消去することにより、安定制約レベルＮ＝４における下限値ｆｓｌｉｍ４は、式２６により表される。
次に、安定制約レベルＮ＝３について説明する。演算部５２は、最大回転速度ｎｍ直後の位相差εがπの時点から位相差εが２πの時点までの時間が、位相差εがπの時の単位時間当たりの第２回転速度ｎ（ｔ）と刃数Ｚとの積の逆数となる様に、下限値ｆｓｌｉｍ４（第１下限値）よりも振動の抑制効果が低い下限値ｆｓｌｉｍ３（第２下限値）を設定する。すなわち、位相差εが安定領域から再び不安定領域に入る際に、不安定領域を切れ刃通過時間１／ｆｚ以内に通り抜けることを仮定する。
Next, the stability constraint level N = 3 will be described. The
図８（Ａ）の矢印Ｎ３に示すように、例えば、切れ刃通過時間１／ｆｚである６０／（（Ｚ×（ｎｍ−Δｎ１））後に回転速度ｎがｎｍ−Δｎ１からｎｍ−Δｎ２になる場合を仮定する。 As indicated by an arrow N3 in FIG. 8 (A), the example, rotational speed n after a cutting edge transit time 1 / f z 60 / (( Z × (n m- Δn 1)) from n m -Δn 1 It is assumed that n m − Δn 2 is obtained.
最大回転速度ｎｍからΔｔ後の回転速度ｎｍ−Δｎ１は、式１を用いて、式２７により表される。
ｎｍ−Δｎ１＝ｎ０（１＋αｓｓｉｎ（２πｆｓΔｔ＋π／２）） （式２７）
同様に、最大回転速度ｎｍからΔｔ＋６０／（Ｚ×（ｎｍ−Δｎ１））後の回転速度ｎｍ−Δｎ２は、式１を用いて、式２８により表される。
６０ｆｃ／（Ｚ×（ｎｍ−Δｎ１））＝ｋｍ＋π／２π （式２９）
The rotation speed n m − Δn 1 after the maximum rotation speed n m to Δt is expressed by the formula 27 using the
n m − Δ n 1 = n 0 (1 + α s sin (2πf s Δt + π / 2)) (Equation 27)
Similarly, the maximum from the rotational speed n m Δt + 60 / (Z × (n m -Δn 1)) rotational speed n m -Δn 2 later, using
60f c / (Z × (n m -Δn 1)) = k m + π / 2π ( Formula 29)
式２７および式２８からΔｔを消去するようにＳＳＶ周波数ｆｓについて整理し、式２５および式２９を代入してΔｎ１およびΔｎ２を消去することにより、安定制約レベルＮ＝３における下限値ｆｓｌｉｍ３は、式３０により表される。
次に、安定制約レベルＮ＝２について説明する。演算部５２は、第２回転速度ｎ（ｔ）が中心回転速度ｎ０と一致する時点から、第２回転速度ｎ（ｔ）が上昇して位相差εがπ／２になる時点までの時間が、単位時間当たりの中心回転速度ｎ０と刃数Ｚとの積の逆数となる様に、下限値ｆｓｌｉｍ３（第２下限値）よりも振動の抑制効果が低い下限値ｆｓｌｉｍ２（第３下限値）を設定する。位相差εがπ／２である場合に、切削力が工具に対して最も大きい負の仕事をする（エネルギーが最も消散される）ことが知られている（非特許文献２）。従って、安定制約レベルＮ＝２では、十分にエネルギーを消散させるため、切れ刃通過時間１／ｆｚ以内に位相差εがπ／２に達することを仮定する。
Next, the stability constraint level N = 2 will be described. The
図８（Ａ）の矢印Ｎ２に示すように、例えば、切れ刃通過時間１／ｆｚ０（６０／（Ｚ×ｎ０））後に回転速度ｎがｎ０からｎａになる場合を仮定する。 As shown by the arrow N2 in FIG. 8 (A), for example, it is assumed that the rotation speed n changes from n 0 to n a after the cutting edge passing time 1 / f z 0 (60 / (Z × n 0)).
中心回転速度ｎ０から６０／（Ｚ×ｎ０）後の回転速度ｎａは、式１を用いて、式３１により表される。
ｎａ＝ｎ０（１＋αｓｓｉｎ（２πｆｓ×６０／（Ｚ×ｎ０））） （式３１）
また、回転速度ｎがｎａである時の瞬間位相差は、π／２である。従って、ｎａは、式７と同様に、式３２により表される。
６０ｆｃ／（Ｚ×ｎａ）＝ｋ０＋（π／２）／２π （式３２）
Rotational speed n a post central rotational speed n 0 from 60 / (Z × n 0), using
n a = n 0 (1 + α s sin (2πf s × 60 / (Z × n 0 ))) (Equation 31)
Further, the instantaneous phase difference when the rotation speed n is n a is π / 2. Thus, n a is as for formula 7, represented by formula 32.
60 f c / (Z × n a ) = k 0 + (π / 2) / 2π (Equation 32)
式３１をＳＳＶ周波数ｆｓについて整理し、式３２を代入してｎａを消去することにより、安定制約レベルＮ＝２における下限値ｆｓｌｉｍ２は、式３３により表される。
次に、安定制約レベルＮ＝１について説明する。演算部５２は、第２回転速度ｎ（ｔ）が中心回転速度ｎ０と一致する時点から、第２回転速度ｎ（ｔ）が上昇して位相差εがπになる時点までの時間が、単位時間当たりの中心回転速度ｎ０と刃数Ｚとの積の逆数となる様に、下限値ｆｓｌｉｍ２（第３下限値）よりも振動の抑制効果が低い下限値ｆｓｌｉｍ１（第４下限値）を設定する。すなわち、例えば、位相差εがＳＳＶ開始時の不安定領域を切れ刃通過時間１／ｆｚ以内に通り抜けることを仮定する。
Next, the stability constraint level N = 1 will be described. The
図８（Ａ）の矢印Ｎ１に示すように、例えば、切れ刃通過時間１／ｆｚ０（６０／（Ｚ×ｎ０））後に回転速度ｎがｎ０からｎｂになる場合を仮定する。
As shown by the arrow N1 in FIG. 8A, it is assumed that the rotation speed n changes from n 0 to n b after, for example, the cutting
中心回転速度ｎ０から６０／（Ｚ×ｎ０）後の回転速度ｎｂは、式１を用いて、式３４により表される。
ｎｂ＝ｎ０（１＋αｓｓｉｎ（２πｆｓ×６０／（Ｚ×ｎ０））） （式３４）
また、回転速度がｎｂである時の瞬間位相差は、πである。従って、ｎｂは、式７と同様に、式３５により表される。
６０ｆｃ／（Ｚ×ｎｂ）＝ｋ０＋π／２π （式３５）
The rotation speed n b after the center rotation speed n 0 to 60 / (Z × n 0 ) is expressed by the formula 34 using the
n b = n 0 (1 + α s sin (2πf s × 60 / (Z × n 0 ))) (Equation 34)
Further, the instantaneous phase difference when the rotation speed is n b is π. Therefore, n b is represented by the formula 35 as in the formula 7.
60f c / (Z × n b ) = k 0 + π / 2π (Equation 35)
式３４をＳＳＶ周波数ｆｓについて整理し、式３５を代入してｎｂを消去することにより、安定制約レベルＮ＝１における下限値ｆｓｌｉｍ１は、式３６により表される。
以上のように、第１実施形態によれば、演算部５２は、振動および振動のエネルギーをベッセル関数を用いて近似し、中心回転速度ｎ０を中心とする回転速度の変動幅および変動周波数を演算する。振動のエネルギーＵは、ＳＳＶにより、解析的に解くことが困難になる。第１実施形態では、式１８に示すように、ベッセル関数を用いて振動のエネルギーＵを単純な式として解析的に解くことができる。さらに、振動のエネルギーＵが負になり、びびり振動を抑制する適切なＲＶＡを選定することができる。この結果、びびり振動を抑制することができる。
As described above, according to the first embodiment, the
また、演算部５２は、ベッセル関数の極小点に基づいて、ＲＶＡを選定することができる。これにより、振動のエネルギーが低く、びびり振動を抑制するより適切なＲＶＡを選定することができる。
Further, the
また、図５（Ａ）および図５（Ｂ）に示すように、ベッセル関数の複数の極小点に基づいて、複数の変動振幅率αｓを選定することができる。すなわち、広い範囲の変動振幅率αｓを選定することができる。図８（Ｂ）において説明したように、ＳＳＶ周波数ｆｓの下限値ｆｓｌｉｍＮは、変動振幅率αｓの大小によって決まる。従って、広い範囲で下限値ｆｓｌｉｍＮを選定することができ、機械的制約などの他の制約を満たす下限値ｆｓｌｉｍＮを得ることができる。この結果、実際の加工条件を考慮して、より適切なＲＶＡおよびＲＶＦを選定することができる。 Further, as shown in FIGS. 5 (A) and 5 (B), a plurality of fluctuation amplitude rates α s can be selected based on a plurality of minimum points of the Bessel function. That is, a wide range of fluctuation amplitude rate α s can be selected. As described in FIG. 8B, the lower limit value f slimN of the SSV frequency f s is determined by the magnitude of the fluctuation amplitude rate α s. Therefore, it is possible to select the lower limit f SlimN in a wide range, it is possible to obtain a lower limit value f SlimN satisfying other constraints such as mechanical constraints. As a result, more appropriate RVA and RVF can be selected in consideration of the actual processing conditions.
尚、工作機械１００は、ボーリング加工に限られず、ミーリング加工や旋削加工などに用いられてもよい。旋削加工の場合、回転軸は、ワークＷを回転させる。この場合、例えば、工具を一定速度で送ることによってワークＷの切削が行われる。
The
また、「ＲＶＡの演算」において、変動振幅率αｓを演算するための式は、式２０に限られない。例えば、「ＳＳＶ周波数ｆｓが主軸１の切れ刃通過周波数ｆｚ（ｔ）（回転速度ｎ（ｔ））に対して十分小さい場合」に加えて、「変動振幅率αｓが十分に小さい場合」には、式１６に示す変位ｘ（ｔ−τ（ｔ））は、式３７により表されてもよい。
ｍｆ＝αｓ（２πｋ０＋ε０） （式３８）
従って、変動振幅率αｓは、式３８を整理することにより、式３９に示すように定義される。
m f = α s (2πk 0 + ε 0 ) (Equation 38)
Therefore, the fluctuation amplitude rate α s is defined as shown in Equation 39 by rearranging Equation 38.
本実施形態による工作機械１００、数値制御装置５および振動抑制方法の少なくとも一部は、ハードウェアで構成してもよいし、ソフトウェアで構成してもよい。ソフトウェアで構成する場合には、工作機械１００、数値制御装置５および振動抑制方法の少なくとも一部の機能を実現するプログラムをフレキシブルディスクやＣＤ−ＲＯＭ等の記録媒体に収納し、コンピュータに読み込ませて実行させてもよい。記録媒体は、磁気ディスクや光ディスク等の着脱可能なものに限定されず、ハードディスク装置やメモリなどの固定型の記録媒体でもよい。また、工作機械１００、数値制御装置５および振動抑制方法の少なくとも一部の機能を実現するプログラムを、インターネット等の通信回線（無線通信も含む）を介して頒布してもよい。さらに、同プログラムを暗号化したり、変調をかけたり、圧縮した状態で、インターネット等の有線回線や無線回線を介して、あるいは記録媒体に収納して頒布してもよい。
At least a part of the
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。 Although some embodiments of the present invention have been described, these embodiments are presented as examples and are not intended to limit the scope of the invention. These embodiments can be implemented in various other forms, and various omissions, replacements, and changes can be made without departing from the gist of the invention. These embodiments and modifications thereof are included in the scope and gist of the invention, as well as in the scope of the invention described in the claims and the equivalent scope thereof.
１００ 工作機械、１ 主軸、１１ 工具、２ 主軸モータ、３ エンコーダ、４ 振動センサ、５ 数値制御装置、５１ センサ情報取得部、５２ 演算部、Ｗ ワーク、αｓ 変動振幅率、ｆｓ ＳＳＶ周波数、ε 位相差、ｎ 回転速度、ｎ０ 中心回転速度、ｎ（ｔ） 第２回転速度、ｎｍ 最大回転速度、εｍ 最大速度位相差、ｆｃ 振動周波数、η ＳＳＶ効率 100 machine machine, 1 spindle, 11 tools, 2 spindle motor, 3 encoder, 4 vibration sensor, 5 numerical control device, 51 sensor information acquisition unit, 52 calculation unit, W work, α s fluctuation amplitude rate, f s SSV frequency, epsilon phase difference, n rotational speed, n 0 central rotation speeds, n (t) the second rotational speed, n m maximum rotational speed, epsilon m maximum velocity phase difference, f c vibration frequency, eta SSV efficiency
Claims (12)
前記回転軸の回転速度を検出する回転速度センサと、
加工中の前記工具または前記ワークに発生する振動を検出する振動センサと、
前記回転軸の前記回転速度を制御する数値制御装置と、を備え、
前記数値制御装置は、
前記回転速度および前記振動の振動周波数を取得するセンサ情報取得部と、
前記振動および前記振動のエネルギーをベッセル関数を用いて近似し、前記ワークを加工する第１回転速度を中心とする前記回転速度の変動幅および変動周波数を演算し、前記変動幅および前記変動周波数を前記第１回転速度に重畳させて前記ワークを加工する第２回転速度を演算する演算部と、を有する、工作機械。 A rotating shaft that rotates the work or the tool that processes the work,
A rotation speed sensor that detects the rotation speed of the rotation shaft, and
A vibration sensor that detects vibration generated in the tool or work during machining, and a vibration sensor.
A numerical control device for controlling the rotation speed of the rotation shaft is provided.
The numerical control device is
A sensor information acquisition unit that acquires the rotation speed and the vibration frequency of the vibration,
The vibration and the energy of the vibration are approximated by using a vessel function, the fluctuation width and the fluctuation frequency of the rotation speed centered on the first rotation speed for processing the work are calculated, and the fluctuation width and the fluctuation frequency are calculated. A machine tool having a calculation unit for calculating a second rotation speed for processing the work by superimposing it on the first rotation speed.
前記回転速度、および、加工中の前記工具または前記ワークに発生する振動の振動周波数を取得するセンサ情報取得部と、
前記振動および前記振動のエネルギーをベッセル関数を用いて近似し、前記ワークを加工する第１回転速度を中心とする前記回転速度の変動幅および変動周波数を演算し、前記変動幅および前記変動周波数を前記第１回転速度に重畳させて前記ワークを加工する第２回転速度を演算する演算部と、を備える、数値制御装置。 A numerical control device that controls the rotation speed of a work or a rotation shaft that rotates a tool for processing the work.
A sensor information acquisition unit that acquires the rotation speed and the vibration frequency of the vibration generated in the tool or the work being machined.
The vibration and the energy of the vibration are approximated by using a vessel function, the fluctuation width and the fluctuation frequency of the rotation speed centered on the first rotation speed for processing the work are calculated, and the fluctuation width and the fluctuation frequency are calculated. A numerical control device including a calculation unit that calculates a second rotation speed for processing the work by superimposing it on the first rotation speed.
前記回転速度および前記振動の振動周波数をセンサ情報取得部により取得し、
前記振動および前記振動のエネルギーをベッセル関数を用いて演算部により近似し、
前記ワークを加工する第１回転速度を中心とする前記回転速度の変動幅および変動周波数を前記演算部により演算し、
前記変動幅および前記変動周波数を前記第１回転速度に重畳させて前記ワークを加工する第２回転速度を前記演算部により演算する、ことを具備する、振動抑制方法。 A rotating shaft for rotating a work or a tool for machining the work, a rotating speed sensor for detecting the rotational speed of the rotating shaft, a vibration sensor for detecting vibration generated in the tool or the work during machining, and the rotation. A method for suppressing vibration of a machine tool including a numerical control device for controlling the rotation speed of the shaft.
The rotation speed and the vibration frequency of the vibration are acquired by the sensor information acquisition unit, and the vibration frequency is acquired.
The vibration and the energy of the vibration are approximated by the calculation unit using the Bessel function.
The calculation unit calculates the fluctuation range and frequency of the rotation speed centered on the first rotation speed for processing the work.
A vibration suppression method comprising superimposing the fluctuation width and the fluctuation frequency on the first rotation speed and calculating a second rotation speed for processing the work by the calculation unit.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019183150A JP7412695B2 (en) | 2019-10-03 | 2019-10-03 | Machine tools, numerical control devices and vibration suppression methods |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2019183150A JP7412695B2 (en) | 2019-10-03 | 2019-10-03 | Machine tools, numerical control devices and vibration suppression methods |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2021058945A true JP2021058945A (en) | 2021-04-15 |
| JP7412695B2 JP7412695B2 (en) | 2024-01-15 |
Family
ID=75381064
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2019183150A Active JP7412695B2 (en) | 2019-10-03 | 2019-10-03 | Machine tools, numerical control devices and vibration suppression methods |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7412695B2 (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009064426A (en) * | 2007-08-10 | 2009-03-26 | Jtekt Corp | Chatter simulation device and chatter simulation method |
| JP2013063497A (en) * | 2011-09-20 | 2013-04-11 | Okuma Corp | Monitoring method and monitor apparatus for rotation speed of rotary shaft in machine tool, and machine tool |
| CN104647132A (en) * | 2014-12-22 | 2015-05-27 | 华中科技大学 | Active control method of milling chatter vibration based on electric spindle of magnetic suspension bearing |
| JP2015112676A (en) * | 2013-12-11 | 2015-06-22 | 国立大学法人 東京大学 | Working method |
| JP2018153907A (en) * | 2017-03-21 | 2018-10-04 | 株式会社ジェイテクト | Grinding work simulation device and method |
-
2019
- 2019-10-03 JP JP2019183150A patent/JP7412695B2/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2009064426A (en) * | 2007-08-10 | 2009-03-26 | Jtekt Corp | Chatter simulation device and chatter simulation method |
| JP2013063497A (en) * | 2011-09-20 | 2013-04-11 | Okuma Corp | Monitoring method and monitor apparatus for rotation speed of rotary shaft in machine tool, and machine tool |
| JP2015112676A (en) * | 2013-12-11 | 2015-06-22 | 国立大学法人 東京大学 | Working method |
| CN104647132A (en) * | 2014-12-22 | 2015-05-27 | 华中科技大学 | Active control method of milling chatter vibration based on electric spindle of magnetic suspension bearing |
| JP2018153907A (en) * | 2017-03-21 | 2018-10-04 | 株式会社ジェイテクト | Grinding work simulation device and method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP7412695B2 (en) | 2024-01-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Davies et al. | The stability of low radial immersion milling | |
| US8229598B2 (en) | Vibration suppressing device for machine tool | |
| JP4434137B2 (en) | Calculation method of axial depth of cut for stability limit of self-excited vibration in cutting of inclined surface | |
| JP5258921B2 (en) | Machine tool and its processing control device | |
| JP4622873B2 (en) | NC program creation method and program | |
| Urbikain et al. | Spindle speed variation technique in turning operations: Modeling and real implementation | |
| JP2012213830A5 (en) | ||
| JP4835442B2 (en) | Calculation method of cutting end coordinates in shoulder cutting using a rotary tool | |
| JP5802062B2 (en) | Machine tool control apparatus and control method | |
| JP5272569B2 (en) | Chatter simulation apparatus and chatter simulation method | |
| Lee et al. | An integrated prediction model including the cutting process for virtual product development of machine tools | |
| JP6139392B2 (en) | Processing method | |
| Yan et al. | Regenerative chatter in a plunge grinding process with workpiece imbalance | |
| JP5628072B2 (en) | Method and apparatus for controlling rotation axis of machine tool | |
| JP5406166B2 (en) | Tool path generation device, tool path generation method and program therefor | |
| JP2015085395A (en) | Processing error amount prediction device and processing error amount prediction program | |
| JP6302794B2 (en) | Rotation speed display method | |
| JP5226484B2 (en) | Chatter vibration suppression method | |
| JP5631779B2 (en) | Vibration suppression method and apparatus for machine tool | |
| JP6336265B2 (en) | Machine tool control device, machine tool control method, and program | |
| JP2021058945A (en) | Machine tool, numerical control device and vibration suppression method | |
| JP5767931B2 (en) | Vibration suppression method and vibration suppression device for machine tool | |
| JP4995115B2 (en) | Vibration suppression method and apparatus | |
| JP2019000943A (en) | Machine tool, processing method, and processing program | |
| JP2021020260A (en) | Spindle rotation number control device of machine tool and control method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A80 | Written request to apply exceptions to lack of novelty of invention |
Free format text: JAPANESE INTERMEDIATE CODE: A80 Effective date: 20191101 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220715 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220720 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20230829 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20230831 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20231005 |
|
| 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: 20231124 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20231219 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7412695 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |