JP6880322B2 - Numerical control device - Google Patents

Numerical control device Download PDF

Info

Publication number
JP6880322B2
JP6880322B2 JP2020525069A JP2020525069A JP6880322B2 JP 6880322 B2 JP6880322 B2 JP 6880322B2 JP 2020525069 A JP2020525069 A JP 2020525069A JP 2020525069 A JP2020525069 A JP 2020525069A JP 6880322 B2 JP6880322 B2 JP 6880322B2
Authority
JP
Japan
Prior art keywords
cutting force
unit
disturbance
control
control device
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
JP2020525069A
Other languages
Japanese (ja)
Other versions
JPWO2019239594A1 (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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric 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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Publication of JPWO2019239594A1 publication Critical patent/JPWO2019239594A1/en
Application granted granted Critical
Publication of JP6880322B2 publication Critical patent/JP6880322B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/404Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by control arrangements for compensation, e.g. for backlash, overshoot, tool offset, tool wear, temperature, machine construction errors, load, inertia
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/41Servomotor, servo controller till figures
    • G05B2219/41376Tool wear, flank and crater, estimation from cutting force

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)
  • Feedback Control In General (AREA)

Description

本発明は、工作機械を制御する数値制御装置に関する。 The present invention relates to a numerical control device that controls a machine tool.

数値制御される工作機械は、アクチュエータと、アクチュエータによって駆動されるテーブル、工具などの被駆動体とを備える。アクチュエータは、回転サーボモータ、リニアサーボモータなどである。工作機械に接続される数値制御装置から出力される位置、経路、速度、トルクなどを示す指令に従って、アクチュエータが制御されることによって、被駆動体が動作する。これにより、テーブルに固定される加工対象である被加工物が工具により加工される。被加工物に対する工具の位置が、指令された経路である指令軌跡に正確に追従するように、被駆動体を駆動する制御は、軌跡制御又は輪郭運動制御と呼ばれる。軌跡制御又は輪郭運動制御は、数値制御装置によって精密に実行される。軌跡制御又は輪郭運動制御が実行されることによって、工作機械に設けられる主軸、送り軸といった複数の軸のそれぞれに繋がるサーボモータが動作して、被加工物に対する工具の位置が制御される。 A numerically controlled machine tool includes an actuator and a driven body such as a table or a tool driven by the actuator. The actuator is a rotary servo motor, a linear servo motor, or the like. The driven body operates by controlling the actuator according to a command indicating a position, a path, a speed, a torque, and the like output from a numerical control device connected to the machine tool. As a result, the workpiece to be machined, which is fixed to the table, is machined by the tool. The control for driving the driven body so that the position of the tool with respect to the workpiece accurately follows the command locus, which is the commanded path, is called locus control or contour motion control. Trajectory control or contour motion control is precisely performed by a numerical control device. By executing the trajectory control or the contour motion control, the servomotors connected to each of the plurality of axes such as the spindle and the feed axis provided in the machine tool operate to control the position of the tool with respect to the workpiece.

工作機械では、被加工物の加工時に、工具の被加工物への接触部分に切削力が生じる。工具の寿命を延ばし、また被加工物の加工精度を向上させるためには、切削力を一定に保つことが重要である。しかしながら、切削工具が例えばエンドミル、ドリルなどである場合、切削工具の軸に設けられる切削刃が軸の回転方向に対して離れて配列されるため、加工中の工具の被加工物への接触部分における切削抵抗が変動して、被加工物への加工が断続的に行われる。従って、エンドミル、ドリルなどの切削工具が用いられる場合、切削力を一定に保つことが困難であり、切削力の変動が、主軸、送り軸を駆動するサーボモータの制御系の外乱となる。切削力の変動によって、例えば送り軸に繋がるサーボモータの制御系に振動等が生じた場合、被加工物を基準とした場合の工具の相対位置も振動的となり、加工精度又は加工面品位が悪化するだけでなく、工具の摩耗が進む。特許文献1には、モータを制御する制御系の周波数応答帯域を、工具摩耗が許容範囲内となるような周波数応答帯域に設定して、被加工物を加工する技術が開示される。特許文献1に開示される技術によれば、加工中の切削抵抗の変動に起因して制御系に外乱が生じた場合でも、工具の摩耗が一定範囲内に抑えられ、加工品質が安定化すると共に加工効率が向上する。 In a machine tool, when machining a work piece, a cutting force is generated at a portion where the tool comes into contact with the work piece. In order to extend the life of the tool and improve the machining accuracy of the workpiece, it is important to keep the cutting force constant. However, when the cutting tool is, for example, an end mill or a drill, the cutting blades provided on the shaft of the cutting tool are arranged apart from each other in the direction of rotation of the shaft. The cutting resistance in the work fluctuates, and the machining to the workpiece is performed intermittently. Therefore, when a cutting tool such as an end mill or a drill is used, it is difficult to keep the cutting force constant, and the fluctuation of the cutting force causes disturbance of the control system of the servomotor that drives the spindle and the feed shaft. When the control system of the servo motor connected to the feed shaft vibrates due to fluctuations in cutting force, the relative position of the tool with respect to the workpiece also becomes vibrating, and the machining accuracy or surface quality deteriorates. Not only does it, but the tool wears out. Patent Document 1 discloses a technique for processing a workpiece by setting the frequency response band of a control system for controlling a motor to a frequency response band such that tool wear is within an allowable range. According to the technique disclosed in Patent Document 1, even if a disturbance occurs in the control system due to fluctuations in cutting resistance during machining, tool wear is suppressed within a certain range and machining quality is stabilized. At the same time, the processing efficiency is improved.

特開2013−250866号公報Japanese Unexamined Patent Publication No. 2013-250866

しかしながら、特許文献1に開示される技術では、周波数応答帯域が狭くなることによって、工具位置の指令軌跡への追従性が低下して、加工精度に悪影響を及ぼす可能性があるという課題がある。 However, the technique disclosed in Patent Document 1 has a problem that the narrowing of the frequency response band may reduce the followability of the tool position to the command locus, which may adversely affect the machining accuracy.

本発明は、上記に鑑みてなされたものであって、指令軌跡への追従性の低下を抑制することができる数値制御装置を得ることを目的とする。 The present invention has been made in view of the above, and an object of the present invention is to obtain a numerical control device capable of suppressing a decrease in followability to a command locus.

上述した課題を解決し、目的を達成するために、本発明の数値制御装置は、被加工物を加工する工具により被加工物の加工を行う機械系を制御する数値制御装置であって、位置指令に基づいて、被加工物が設けられる被駆動体の位置が位置指令に追従するように被駆動体に接続されるモータを駆動するサーボ制御部を備える。数値制御装置は、工具と被加工物と間に生じる切削力を測定又は推定して出力する切削力出力部と、被駆動体の位置と切削力とに基づき切削力ベクトルを演算する切削力ベクトル演算部と、切削力に対する機械系の動特性を同定する機械特性同定部を備える。数値制御装置は、切削力ベクトルと機械特性同定部で同定された機械系の動特性とに基づき、サーボ制御部の制御特性を変更する制御特性変更部とを備えることを特徴とする。 In order to solve the above-mentioned problems and achieve the object, the numerical control device of the present invention is a numerical control device that controls a mechanical system that processes a workpiece by a tool that processes the workpiece, and is a position. A servo control unit that drives a motor connected to the driven body so that the position of the driven body on which the workpiece is provided follows the position command based on the command is provided. The numerical control device measures or estimates and outputs the cutting force generated between the tool and the workpiece, and the cutting force that calculates the cutting force vector based on the position of the driven body and the cutting force. It is equipped with a vector calculation unit and a mechanical characteristic identification unit that identifies the dynamic characteristics of the mechanical system with respect to the cutting force. The numerical control device is characterized by including a control characteristic changing unit that changes the control characteristics of the servo control unit based on the cutting force vector and the dynamic characteristics of the mechanical system identified by the mechanical characteristic identification unit.

本発明によれば、指令軌跡への追従性の低下を抑制することができるという効果を奏する。 According to the present invention, there is an effect that a decrease in followability to a command locus can be suppressed.

本発明の実施の形態1に係る数値制御装置の構成を示す図The figure which shows the structure of the numerical control apparatus which concerns on Embodiment 1 of this invention. 図1に示すモータ及び機械系の構成を示す図The figure which shows the structure of the motor and the mechanical system shown in FIG. 図1に示す模擬応答演算部の構成を示す図The figure which shows the structure of the simulated response calculation part shown in FIG. 図1に示す模擬応答演算部の変形例を示す図The figure which shows the modification of the simulated response calculation part shown in FIG. 図1に示す機械特性同定部で算出される機械特性情報を説明するための第1図FIG. 1 for explaining the machine characteristic information calculated by the machine characteristic identification unit shown in FIG. 図1に示す機械特性同定部で算出される機械特性情報を説明するための第2図FIG. 2 for explaining the machine characteristic information calculated by the machine characteristic identification unit shown in FIG. 図1に示す制御特性変更部の動作を説明するためのフローチャートA flowchart for explaining the operation of the control characteristic changing unit shown in FIG. 切削力により、図1に示す機械端位置に外乱振動が生じている例を示す第1図FIG. 1 shows an example in which disturbance vibration is generated at the machine end position shown in FIG. 1 due to a cutting force. 切削力により、図1に示す機械端位置に外乱振動が生じている例を示す第2図FIG. 2 shows an example in which disturbance vibration is generated at the machine end position shown in FIG. 1 due to a cutting force. 本発明の実施の形態2に係る数値制御装置の構成を示す図The figure which shows the structure of the numerical control apparatus which concerns on Embodiment 2 of this invention. 本発明の実施の形態3に係る数値制御装置の構成を示す図The figure which shows the structure of the numerical control apparatus which concerns on Embodiment 3 of this invention. 図2に示すモータ及び機械系を2慣性のモデルで近似したブロック図A block diagram of the motor and mechanical system shown in FIG. 2 approximated by a two-inertia model. 実施の形態1から3に係る数値制御装置のハードウェア構成例を示す図The figure which shows the hardware configuration example of the numerical control apparatus which concerns on Embodiments 1 to 3.

以下に、本発明の実施の形態に係る数値制御装置を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。 Hereinafter, the numerical control device according to the embodiment of the present invention will be described in detail with reference to the drawings. The present invention is not limited to this embodiment.

実施の形態1.
図1は本発明の実施の形態1に係る数値制御装置の構成を示す図である。図2は図1に示すモータ及び機械系の構成を示す図である。数値制御装置100−1は、位置指令200及び機械端位置201に基づきトルク指令202を生成し、生成したトルク指令202を工作機械100に設けられるモータ7へ与えることにより、制御対象である機械系8を駆動する。数値制御装置100−1は、サーボ制御部1、切削力出力部2、切削力ベクトル演算部3、模擬応答演算部4、機械特性同定部5、制御特性変更部6及び減算部9を備える。
Embodiment 1.
FIG. 1 is a diagram showing a configuration of a numerical control device according to a first embodiment of the present invention. FIG. 2 is a diagram showing the configuration of the motor and the mechanical system shown in FIG. The numerical control device 100-1 generates a torque command 202 based on the position command 200 and the machine end position 201, and gives the generated torque command 202 to the motor 7 provided in the machine tool 100 to control the mechanical system. Drive 8. The numerical control device 100-1 includes a servo control unit 1, a cutting force output unit 2, a cutting force vector calculation unit 3, a simulated response calculation unit 4, a mechanical characteristic identification unit 5, a control characteristic change unit 6, and a subtraction unit 9.

サーボ制御部1はモータ7に接続されており、モータ7にトルク指令202を与えることにより、制御対象である機械系8を駆動する。図2に示すようにモータ7は、サーボモータ71と、モータ位置検出器72とを備える。サーボモータ71のロータ71aにはシャフト71bが設けられる。モータ位置検出器72には、ロータ71aの回転位置を検出するロータリエンコーダが用いられる。機械系8は、テーブル81、ボールねじナット82、ボールねじ83、カップリング84、機械端位置検出器85及び機械端位置検出器ヘッド86を備える。機械系8のボールねじ83は、カップリング84を介して、サーボモータ71のシャフト71bに接続される。ボールねじ83には、ボールねじナット82が嵌め合わされ、ボールねじナット82には、被駆動体であるテーブル81が固定される。テーブル81、不図示の案内機構によって、ボールねじ83がびる方向に移動可能に支持されている。機械端位置検出器85及び機械端位置検出器ヘッド86は、数値制御装置100−1の制御対象であるテーブル81の位置を検出するためのエンコーダである。機械端位置検出器85は、例えばリニアエンコーダであり、機械端位置検出器ヘッド86は、例えばリニアエンコーダヘッドである。サーボ制御部1に入力される機械端位置201は、モータ位置検出器72により検出されたロータ71aの回転位置を示す情報と、機械端位置検出器85により検出されたテーブル81の位置を示す情報との少なくとも一方を示す。 The servo control unit 1 is connected to the motor 7, and by giving the torque command 202 to the motor 7, the mechanical system 8 to be controlled is driven. As shown in FIG. 2, the motor 7 includes a servomotor 71 and a motor position detector 72. A shaft 71b is provided on the rotor 71a of the servomotor 71. As the motor position detector 72, a rotary encoder that detects the rotational position of the rotor 71a is used. The mechanical system 8 includes a table 81, a ball screw nut 82, a ball screw 83, a coupling 84, a mechanical end position detector 85, and a mechanical end position detector head 86. The ball screw 83 of the mechanical system 8 is connected to the shaft 71b of the servomotor 71 via the coupling 84. A ball screw nut 82 is fitted to the ball screw 83, and a table 81, which is a driven body, is fixed to the ball screw nut 82. Table 81, by a guide mechanism (not shown), a ball screw 83 is movably supported in the extension building direction. The machine end position detector 85 and the machine end position detector head 86 are encoders for detecting the position of the table 81, which is the control target of the numerical control device 100-1. The mechanical end position detector 85 is, for example, a linear encoder, and the mechanical end position detector head 86 is, for example, a linear encoder head. The machine end position 201 input to the servo control unit 1 is information indicating the rotation position of the rotor 71a detected by the motor position detector 72 and information indicating the position of the table 81 detected by the machine end position detector 85. At least one of them is shown.

なお、機械端位置検出器85ではテーブル81の移動距離が測定されるのに対して、モータ位置検出器72ではロータ71aの回転角度が測定される。しかし、この回転角度にモータ7の1回転当たりのテーブル移動距離であるボールねじリードを乗じて、ロータ71aの1回転の角度2π[rad]で除することで、サーボ制御部1では、ロータ71aの回転角度が、テーブル81の移動方向の長さに換算される。実施の形態1では、ロータ71aの回転角度をテーブル81の移動方向に換算した値が用いられるものとする。 The mechanical end position detector 85 measures the moving distance of the table 81, whereas the motor position detector 72 measures the rotation angle of the rotor 71a. However, by multiplying this rotation angle by the ball screw lead, which is the table movement distance per rotation of the motor 7, and dividing by the angle of one rotation of the rotor 71a, 2π [rad], the rotor 71a in the servo control unit 1 The rotation angle of is converted into the length of the table 81 in the moving direction. In the first embodiment, it is assumed that a value obtained by converting the rotation angle of the rotor 71a into the moving direction of the table 81 is used.

機械端位置201として、モータ位置検出器72により検出された位置情報のみ用いるフィードバック制御をセミクローズドループ制御と呼ぶ。機械端位置201として、モータ位置検出器72により検出された位置情報と、機械端位置検出器85により検出された位置情報との両方を用いるフィードバック制御、又は機械端位置検出器85により検出された位置情報のみ用いるフィードバック制御を、フルクローズドループ制御と呼ぶ。実施の形態1に係る数値制御装置100−1には、セミクローズドループ制御とフルクローズドループ制御との何れを適用してもよい。以下では、フルクローズドループ制御が適用される数値制御装置100−1の構成例について説明する。また以下では、フィードバック制御に用いられる位置情報を機械端位置と呼ぶ。なお、機械系8の構成は一例であり、機械系8の構成は図2に示される例に限定されず、数値制御装置100−1は、複数の機械系を制御対象とすることも可能である。 Feedback control that uses only the position information detected by the motor position detector 72 as the machine end position 201 is called semi-closed loop control. As the machine end position 201, feedback control using both the position information detected by the motor position detector 72 and the position information detected by the machine end position detector 85, or the machine end position detector 85 detected. Feedback control that uses only position information is called full closed loop control. Either semi-closed loop control or fully closed loop control may be applied to the numerical control device 100-1 according to the first embodiment. Hereinafter, a configuration example of the numerical control device 100-1 to which the full closed loop control is applied will be described. In the following, the position information used for feedback control will be referred to as the machine end position. The configuration of the mechanical system 8 is an example, and the configuration of the mechanical system 8 is not limited to the example shown in FIG. 2, and the numerical control device 100-1 can control a plurality of mechanical systems. is there.

図1に示すサーボ制御部1には位置指令200及び機械端位置201が入力される。サーボ制御部1は、機械端位置201を位置指令200に追従させるトルク指令202を生成する。サーボ制御部1は、例えば、位置指令200と機械端位置201との差である位置偏差を求め、位置偏差に対して、例えば比例制御といった位置制御処理を実行して、速度指令を算出する。さらに、速度指令に対して、例えば比例積分制御といった速度制御処理を実行することによって、トルク指令202を算出して出力する。 The position command 200 and the machine end position 201 are input to the servo control unit 1 shown in FIG. The servo control unit 1 generates a torque command 202 that causes the machine end position 201 to follow the position command 200. For example, the servo control unit 1 obtains a position deviation which is a difference between the position command 200 and the machine end position 201, executes a position control process such as proportional control for the position deviation, and calculates a speed command. Further, the torque command 202 is calculated and output by executing a speed control process such as proportional integration control for the speed command.

図1に示す切削力出力部2は、例えばテーブル81に設けられる被加工物90の加工時に、工具101の被加工物90への接触部分に生じる切削力を測定し、測定された切削力を示す情報である切削力203を出力する。 The cutting force output unit 2 shown in FIG. 1 measures, for example, the cutting force generated at the contact portion of the tool 101 with the workpiece 90 when the workpiece 90 provided on the table 81 is machined, and measures the measured cutting force. The cutting force 203, which is the information to be shown, is output.

切削力203の測定には、例えばロードセル、動力計などが用いられる。ロードセル、動力計などは、切削力203を直接測定することが可能な測定器である。なお、切削力203はモータ電流からも推定することも可能であり、モータ電流から切削力203を推定する場合、外乱オブザーバと呼ばれる手法が用いられる。なお、切削力出力部2には、切削力203を直接測定する手段の代わりに、切削力203を間接的に測定する手段を用いてもよい。 For the measurement of the cutting force 203, for example, a load cell, a power meter, or the like is used. A load cell, a power meter, or the like is a measuring instrument capable of directly measuring a cutting force 203. The cutting force 203 can also be estimated from the motor current, and when the cutting force 203 is estimated from the motor current, a method called a disturbance observer is used. The cutting force output unit 2 may use a means for indirectly measuring the cutting force 203 instead of a means for directly measuring the cutting force 203.

切削力ベクトル演算部3は、測定された切削力203と機械端位置201とに基づき、例えば、3つの切削力ベクトルを演算する。3つの切削力ベクトルの内、第1切削力ベクトルは、例えば機械のX軸、Y軸、Z軸などの駆動軸に平行な方向を基準とするベクトルである。以下では、第1切削力ベクトルを機械軸方向切削力ベクトルと呼ぶ。3つの切削力ベクトルの内、第2切削力ベクトルは、工具101の送り方向に平行な切削力ベクトルである。以下では、第2切削力ベクトルを工具送り方向切削力ベクトルと呼ぶ。3つの切削力ベクトルの内、第3切削力ベクトルは、工具101の送り方向に垂直な方向の切削力ベクトルである。以下では、第3切削力ベクトルを垂直方向切削力ベクトルと呼ぶ。 The cutting force vector calculation unit 3 calculates, for example, three cutting force vectors based on the measured cutting force 203 and the machine end position 201. Of the three cutting force vectors, the first cutting force vector is a vector based on a direction parallel to the drive axis such as the X-axis, Y-axis, and Z-axis of the machine. Hereinafter, the first cutting force vector is referred to as a machine axial cutting force vector. Of the three cutting force vectors, the second cutting force vector is a cutting force vector parallel to the feed direction of the tool 101. Hereinafter, the second cutting force vector is referred to as a tool feed direction cutting force vector. Of the three cutting force vectors, the third cutting force vector is a cutting force vector in the direction perpendicular to the feed direction of the tool 101. Hereinafter, the third cutting force vector is referred to as a vertical cutting force vector.

ここでは、動力計を用いて切削力を測定する場合の例について説明する。動力計を用いて切削力を測定する場合、一般的には、機械に設置される動力計の向きを調整することによって、機械の駆動軸、例えばX軸、Y軸、Z軸のそれぞれがびる方向と平行な方向の切削力が測定される。そのため、動力計からの出力には、機械軸方向切削力ベクトルの成分が含まれる。なお、動力計には、3つの軸の方向に対応する機械軸方向切削力ベクトルを測定できるものだけでなく、X軸、Y軸、Z軸の内、1つの軸の方向に対応する機械軸方向切削力ベクトルのみ測定できるものも存在する。そして、測定されない方向の切削力は0とみなされる。なお、5つの回転軸を有する5軸制御工作機械では、動力計の設置方向が回転軸の向きによって変化する場合があり、このような場合、回転軸の角度が、機械軸方向切削力ベクトルの計算に利用される。 Here, an example in which the cutting force is measured using a power meter will be described. When measuring cutting force using a power meter, generally, by adjusting the direction of the power meter installed in the machine, the drive shafts of the machine, such as the X-axis, Y-axis, and Z-axis, are extended. The cutting force in the direction parallel to the direction of vibration is measured. Therefore, the output from the power meter includes a component of the machine axial cutting force vector. The power meter is not only capable of measuring the mechanical axial cutting force vector corresponding to the directions of the three axes, but also the mechanical axis corresponding to the direction of one of the X-axis, Y-axis, and Z-axis. Some can measure only the directional cutting force vector. Then, the cutting force in the unmeasured direction is regarded as 0. In a 5-axis control machine having five rotation axes, the installation direction of the power meter may change depending on the direction of the rotation axis. In such a case, the angle of the rotation axis is the cutting force vector in the machine axis direction. Used for calculation.

切削力ベクトル演算部3では、機械端位置201を微分することにより、複数の工具101のそれぞれの軸の移動速度を計算できる。工作機械100が例えば3つの工具101を備える場合、切削力ベクトル演算部3は、3つの工具101のそれぞれの軸の移動速度を成分とする移動速度ベクトルvを演算する。また切削力ベクトル演算部3では、工具送り方向切削力ベクトルFtが(1)式及び(2)式により演算される。「u」は単位ベクトルである。「v」は移動速度ベクトルである。「F」は機械軸方向切削力ベクトルである。「・」は内積を示す。また切削力ベクトル演算部3では、垂直方向切削力ベクトルFnが(1)式及び(3)式により演算される。
u=v/|v|・・・(1)式
Ft=(F・u)・・・(2)式
Fn=F−Ft・・・(3)式
The cutting force vector calculation unit 3 can calculate the moving speed of each axis of the plurality of tools 101 by differentiating the machine end position 201. When the machine tool 100 includes, for example, three tools 101, the cutting force vector calculation unit 3 calculates a movement speed vector v whose component is the movement speed of each axis of the three tools 101. Further, in the cutting force vector calculation unit 3, the tool feed direction cutting force vector Ft is calculated by the equations (1) and (2). “U” is a unit vector. “V” is a moving speed vector. "F" is a machine axial cutting force vector. "・" Indicates the inner product. Further, in the cutting force vector calculation unit 3, the vertical cutting force vector Fn is calculated by the equations (1) and (3).
u = v / | v | ... (1) formula Ft = (F · u) ... (2) formula Fn = F-Ft ... (3) formula

なお、切削力ベクトル演算部3は、機械端位置201を微分して得られる移動速度の代わりに、速度制御に用いる速度フィードバック情報を用いて、工具送り方向切削力ベクトル及び垂直方向切削力ベクトルを演算するように構成してもよい。また切削力ベクトル演算部3は、機械軸方向切削力ベクトルFの大きさを演算し、演算された機械軸方向切削力ベクトルFの大きさを、切削力の大きさとする。 The cutting force vector calculation unit 3 uses the speed feedback information used for speed control instead of the moving speed obtained by differentiating the machine end position 201 to obtain the tool feed direction cutting force vector and the vertical cutting force vector. It may be configured to calculate. Further, the cutting force vector calculation unit 3 calculates the magnitude of the machine axial cutting force vector F, and sets the calculated magnitude of the machine axial cutting force vector F as the magnitude of the cutting force.

さらに、切削力ベクトル演算部3は、現時点から一定時間遡った時点までに演算された複数の軸のそれぞれに対応する機械軸方向切削力ベクトルの成分を示すデータを保持して、保持されたデータを利用して、各周波数に対する切削力のパワースペクトルを計算する。なお、周波数に対するパワースペクトルは、フーリエ変換を行うことにより演算することができる。 Further, the cutting force vector calculation unit 3 holds data indicating the components of the machine axial cutting force vector corresponding to each of the plurality of axes calculated from the present time to the time point retroactive to a certain time, and the held data. Is used to calculate the power spectrum of the cutting force for each frequency. The power spectrum with respect to frequency can be calculated by performing a Fourier transform.

切削力ベクトル演算部3は、各切削力ベクトル、切削力の大きさ、及び各機械軸方向切削力ベクトルの各成分のパワースペクトルを、切削力ベクトル情報204として、制御特性変更部6に出力する。 The cutting force vector calculation unit 3 outputs the power spectrum of each component of each cutting force vector, the magnitude of the cutting force, and each machine axial cutting force vector to the control characteristic changing unit 6 as the cutting force vector information 204. ..

図3は図1に示す模擬応答演算部の構成を示す図である。模擬応答演算部4は、位置指令200に基づいて機械端位置201の応答を模擬し、機械端位置201が追従するべき理想的な位置の計算結果であるモデル位置205を算出して出力する。模擬応答演算部4は、位置制御模擬部41及び積分演算部42を備える。位置制御模擬部41は、位置指令200と積分演算部42の出力との差に対して、サーボ制御部1で行われる位置制御処理と同じ処理を実行する。積分演算部42は、位置制御模擬部41の処理結果を積分することにより、機械端位置201の値を求めて、モデル位置205として出力する。 FIG. 3 is a diagram showing the configuration of the simulated response calculation unit shown in FIG. The simulated response calculation unit 4 simulates the response of the machine end position 201 based on the position command 200, calculates and outputs the model position 205, which is the calculation result of the ideal position that the machine end position 201 should follow. The simulated response calculation unit 4 includes a position control simulation unit 41 and an integration calculation unit 42. The position control simulating unit 41 executes the same processing as the position control processing performed by the servo control unit 1 with respect to the difference between the position command 200 and the output of the integration calculation unit 42. The integration calculation unit 42 obtains the value of the machine end position 201 by integrating the processing result of the position control simulation unit 41 and outputs it as the model position 205.

図4は図1に示す模擬応答演算部の変形例を示す図である。図4に示す模擬応答演算部4Aは、図3に示す位置制御模擬部41及び積分演算部42に加えて、速度制御模擬部43、積分演算器44及び比例定数演算部45を備える。模擬応答演算部4Aは、モデル位置205に加えて、機械系8の理想的な速度の値を示すモデル速度208と、モータトルクの理想的な出力を示す情報であるモデルトルク209とを演算する。比例定数演算部45では、モデル速度208にトルク変換のための定数を積算することにより、モデルトルク209が演算される。Kpp、Kvp、Kvi、Kt及びJのそれぞれは比例定数を示す。sは微分を表すラプラス演算子である。1/sは積分を表す。位置制御模擬部41では一般的な比例制御が行われる。速度制御模擬部43では一般的な比例積分制御が行われる。速度制御模擬部43の出力は加速度の次元になっており、速度制御模擬部43の出力は、積分演算器44に入力され、積分演算器44及び積分演算部42によって、2段積分される。これにより機械端位置が模擬される。また速度制御模擬部43の出力は、比例定数演算部45に入力される。比例定数演算部45は、速度制御模擬部43の出力に比例定数Ktを掛けることにより、速度制御模擬部43の出力をトルクに変換する。 FIG. 4 is a diagram showing a modified example of the simulated response calculation unit shown in FIG. The simulated response calculation unit 4A shown in FIG. 4 includes a speed control simulation unit 43, an integral calculation unit 44, and a proportional constant calculation unit 45 in addition to the position control simulation unit 41 and the integration calculation unit 42 shown in FIG. In addition to the model position 205, the simulated response calculation unit 4A calculates the model speed 208 indicating the ideal speed value of the mechanical system 8 and the model torque 209 which is information indicating the ideal output of the motor torque. .. The proportional constant calculation unit 45 calculates the model torque 209 by integrating the constant for torque conversion with the model speed 208. Each of Kpp, Kbp, Kvi, Kt and J shows a proportionality constant. s is a Laplace operator that represents the derivative. 1 / s represents the integral. In the position control simulating unit 41, general proportional control is performed. In the speed control simulation unit 43, general proportional integration control is performed. The output of the speed control simulation unit 43 is in the dimension of acceleration, and the output of the speed control simulation unit 43 is input to the integration calculation unit 44, and is integrated in two stages by the integration calculation unit 44 and the integration calculation unit 42. This simulates the machine end position. Further, the output of the speed control simulation unit 43 is input to the proportional constant calculation unit 45. The proportional constant calculation unit 45 converts the output of the speed control simulation unit 43 into torque by multiplying the output of the speed control simulation unit 43 by the proportional constant Kt.

減算部9は、模擬応答演算部4で演算されたモデル位置205を機械端位置201から減算することにより、切削力等の制御への外乱によって生じる位置外乱206を算出する。 The subtraction unit 9 calculates the position disturbance 206 caused by the disturbance to the control such as the cutting force by subtracting the model position 205 calculated by the simulated response calculation unit 4 from the machine end position 201.

機械特性同定部5は、位置外乱206及び切削力203を用いて、切削力203に対する工作機械100の機械特性を同定する。実施の形態1では、切削力203に対する機械特性として、切削力203から制御系における位置外乱206までの伝達関数を同定する例について説明する。なお、機械特性同定部5は、位置外乱206の代わりに、図4に示す模擬応答演算部4Aにより演算されるモデル速度208、モデルトルク209などを用いて、速度外乱、トルク外乱などを算出し、切削力203から速度外乱までの伝達関数、又は切削力203からトルク外乱までの伝達関数を同定するように構成してもよい。 The machine characteristic identification unit 5 identifies the mechanical characteristics of the machine tool 100 with respect to the cutting force 203 by using the position disturbance 206 and the cutting force 203. In the first embodiment, an example of identifying a transfer function from the cutting force 203 to the position disturbance 206 in the control system as the mechanical characteristics with respect to the cutting force 203 will be described. The mechanical characteristic identification unit 5 calculates the speed disturbance, torque disturbance, etc. by using the model speed 208, the model torque 209, etc. calculated by the simulated response calculation unit 4A shown in FIG. 4 instead of the position disturbance 206. , The transfer function from the cutting force 203 to the velocity disturbance, or the transfer function from the cutting force 203 to the torque disturbance may be identified.

機械特性同定部5は、切削力203及び位置外乱206に対してフーリエ変換などの数値処理を行い、周波数伝達関数を算出する。機械特性同定部5は、算出した周波数伝達関数を機械特性情報207として制御特性変更部6に出力する。図5は図1に示す機械特性同定部で算出される機械特性情報を説明するための第1図である。図5の横軸は対数周波数を表し、図5の縦軸はデシベル値で示されるゲインを表す。図5では、周波数伝達関数がゲイン線図により表される。ゲイン線図では、共振点が正のピークとして表れ、反共振点が負のピークとして現れる。図6は図1に示す機械特性同定部で算出される機械特性情報を説明するための第2図である。図6の横軸は対数周波数を表し、図6の縦軸は位相を表す。図6では、周波数伝達関数が位相線図により表される。なお、周波数伝達関数は、横軸を対数周波数とし、縦軸をフーリエ変換後の実部及び虚部によって表されるコクアド線図でも表すことができる。 The mechanical property identification unit 5 performs numerical processing such as Fourier transform on the cutting force 203 and the position disturbance 206 to calculate the frequency transfer function. The mechanical characteristic identification unit 5 outputs the calculated frequency transfer function as mechanical characteristic information 207 to the control characteristic changing unit 6. FIG. 5 is a diagram for explaining the machine characteristic information calculated by the machine characteristic identification unit shown in FIG. The horizontal axis of FIG. 5 represents the logarithmic frequency, and the vertical axis of FIG. 5 represents the gain represented by the decibel value. In FIG. 5, the frequency transfer function is represented by a gain diagram. In the gain diagram, the resonance point appears as a positive peak and the antiresonance point appears as a negative peak. FIG. 6 is a second diagram for explaining the machine characteristic information calculated by the machine characteristic identification unit shown in FIG. The horizontal axis of FIG. 6 represents the logarithmic frequency, and the vertical axis of FIG. 6 represents the phase. In FIG. 6, the frequency transfer function is represented by a phase diagram. The frequency transfer function can also be represented by a coquad diagram whose horizontal axis is a logarithmic frequency and whose vertical axis is represented by a real part and an imaginary part after Fourier transform.

制御系における外乱には、切削力の変動による外乱以外にも、摩擦による外乱が存在し、摩擦によっても位置外乱206が生じる。代表的な例として、円弧の象限切り替え位置において、送り軸の移動方向が反転する際に生じる象限突起がある。機械特性同定部5では、切削力に対する機械特性を同定するために、摩擦の影響を除去する必要がある。そのため、送り軸の移動方向が反転する際に生じる摩擦の影響、又は送り軸が停止状態から動き出すときに生じる摩擦の影響が大きい場合を除いて、切削力に対する機械特性を同定することが望ましい。 Disturbances in the control system include disturbances due to friction in addition to disturbances due to fluctuations in cutting force, and positional disturbance 206 is also generated by friction. As a typical example, there is a quadrant protrusion that occurs when the moving direction of the feed axis is reversed at the quadrant switching position of the arc. In the mechanical property identification unit 5, it is necessary to remove the influence of friction in order to identify the mechanical property with respect to the cutting force. Therefore, it is desirable to identify the mechanical properties for the cutting force unless the influence of friction generated when the moving direction of the feed shaft is reversed or the influence of friction generated when the feed shaft starts moving from the stopped state is large.

制御特性変更部6は、切削力ベクトル演算部3で演算された切削力ベクトル情報204と、機械特性同定部5で同定された機械特性情報207とに基づき、サーボ制御部1の周波数特性の変更を行う。工作機械100において2軸以上で行われる軌跡制御では、制御する軸の目標値応答性が異なると、指令軌跡に対して軌跡誤差が生じる。実施の形態1では、このような軌跡誤差を生じないため、サーボ制御部1に2自由度制御が採用される。2自由度制御とは、目標値への追従と外乱の抑制のために、2つの独立した補償器を利用する制御方式である。外乱抑制のための補償器の周波数応答が制御軸によって異なる場合でも軌跡誤差は生じないため、制御特性変更部6は、外乱抑制のための補償器である外乱応答補償器の周波数応答を変更することにより、外乱抑制を実現する。外乱応答補償器は、一般的にP(Proportional)制御器、PI(Proportional Integral)制御器などで構成される。すなわち、サーボ制御部1において、2自由度制御を採用することにより、軌跡誤差を生じることなく、切削力による外乱を抑制することが可能となる。 The control characteristic changing unit 6 changes the frequency characteristics of the servo control unit 1 based on the cutting force vector information 204 calculated by the cutting force vector calculation unit 3 and the mechanical characteristic information 207 identified by the mechanical characteristic identification unit 5. I do. In the locus control performed on two or more axes in the machine tool 100, if the target value responsiveness of the controlled axes is different, a locus error occurs with respect to the command locus. In the first embodiment, since such a trajectory error does not occur, the servo control unit 1 employs two-degree-of-freedom control. The two-degree-of-freedom control is a control method that uses two independent compensators to follow the target value and suppress disturbance. Since no locus error occurs even when the frequency response of the compensator for disturbance suppression differs depending on the control axis, the control characteristic changing unit 6 changes the frequency response of the disturbance response compensator, which is a compensator for disturbance suppression. As a result, disturbance suppression is realized. The disturbance response compensator is generally composed of a P (Proportional) controller, a PI (Proportional Integral) controller, and the like. That is, by adopting the two-degree-of-freedom control in the servo control unit 1, it is possible to suppress the disturbance due to the cutting force without causing a locus error.

機械特性同定部5で同定される周波数伝達関数にピークが存在し、このピークが共振点である場合、機械端位置で観測可能な範囲において、切削力に起因する振動が生じている可能性がある。フルクローズドループ制御の場合、図2に示すモータ7及び機械系8において振動が生じており、セミクローズドループ制御の場合、モータ7において振動が生じていることとなる。この場合、共振点における外乱応答補償器の応答性を高める、すなわち、共振点における外乱応答補償器の制御ゲインを上げることにより、外乱を抑制することができる。 If there is a peak in the frequency transfer function identified by the machine characteristic identification unit 5 and this peak is the resonance point, there is a possibility that vibration due to the cutting force is occurring in the range observable at the machine end position. is there. In the case of full closed loop control, vibration is generated in the motor 7 and the mechanical system 8 shown in FIG. 2, and in the case of semi-closed loop control, vibration is generated in the motor 7. In this case, the disturbance can be suppressed by increasing the responsiveness of the disturbance response compensator at the resonance point, that is, increasing the control gain of the disturbance response compensator at the resonance point.

一方、周波数伝達関数のピークが反共振点である場合、機械端位置では観測不可能な範囲において、切削力に起因する振動が発生する可能性がある。例えば、テーブル81上に設置されている被加工物90、工具101などが振動している可能性がある。この場合、反共振点における外乱応答補償器の応答性を下げる、すなわち、反共振点における外乱応答補償器の制御ゲインを下げることにより、切削力による外乱を抑制できる。 On the other hand, when the peak of the frequency transfer function is the antiresonance point, vibration due to the cutting force may occur in a range that cannot be observed at the mechanical end position. For example, there is a possibility that the workpiece 90, the tool 101, etc. installed on the table 81 are vibrating. In this case, the disturbance due to the cutting force can be suppressed by lowering the responsiveness of the disturbance response compensator at the antiresonance point, that is, lowering the control gain of the disturbance response compensator at the antiresonance point.

次に制御特性変更部6の動作について説明する。図7は図1に示す制御特性変更部の動作を説明するためのフローチャートである。まず制御特性変更部6は、機械特性同定部5で算出された周波数伝達関数において、データサンプリング周波数の1/2の周波数の範囲で、ピークを探索する(ステップS1)。例えば、サンプリング周期を1kHzとした場合、ピークを探索する範囲は、1kHzの1/2の周波数までの範囲、すなわち0〜500Hzとなる。そして、500Hzから1kHzまでの範囲ではエリアシングが起きているため、制御特性変更部6ではピークの探索が行われない。 Next, the operation of the control characteristic changing unit 6 will be described. FIG. 7 is a flowchart for explaining the operation of the control characteristic changing unit shown in FIG. First, the control characteristic changing unit 6 searches for a peak in the frequency range of 1/2 of the data sampling frequency in the frequency transfer function calculated by the mechanical characteristic identification unit 5 (step S1). For example, when the sampling period is 1 kHz, the range for searching for a peak is a range up to 1/2 of 1 kHz, that is, 0 to 500 Hz. Since aliasing occurs in the range from 500 Hz to 1 kHz, the control characteristic changing unit 6 does not search for the peak.

制御特性変更部6は、周波数伝達関数に反共振点又は共振点が存在するか否かを判定する(ステップS2)。例えば、ゲインの値が低下から増加に転じる第1ピーク点が存在するか否か、ゲインの値が増加から低下に転じる第2ピーク点が存在するか否かにより、反共振点又は共振点の有無が判断される。 The control characteristic changing unit 6 determines whether or not an anti-resonance point or a resonance point exists in the frequency transfer function (step S2). For example, depending on whether or not there is a first peak point at which the gain value changes from decreasing to increasing, and whether or not there is a second peak point where the gain value changes from increasing to decreasing, the antiresonance point or resonance point of the resonance point. The presence or absence is judged.

周波数伝達関数に反共振点又は共振点が存在しない場合(ステップS2,No)、制御特性変更部6は、ステップS1,S2の処理を繰り返す。 When the anti-resonance point or the resonance point does not exist in the frequency transfer function (steps S2 and No), the control characteristic changing unit 6 repeats the processes of steps S1 and S2.

周波数伝達関数に反共振点又は共振点が存在する場合(ステップS2,Yes)、制御特性変更部6は、周波数伝達関数のピークが共振点であるか、反共振点であるかを判定する(ステップS3)。例えば、ゲインの値が低下から増加に転じる第1ピーク点が存在する場合、制御特性変更部6は、周波数伝達関数に反共振点が存在すると判定する。ゲインの値が増加から低下に転じる第2ピーク点が存在する場合、制御特性変更部6は、周波数伝達関数に共振点が存在すると判定する。なお、走査中に複数の第1ピーク点が存在する場合、制御特性変更部6は、複数の第1ピーク点を記憶し、その中で最低のものを反共振点と判断する。また、走査中に複数の第2ピーク点が存在する場合、制御特性変更部6は、複数の第2ピーク点を記憶し、その中で最大のものを共振点と判断する。 When the frequency transfer function has an anti-resonance point or a resonance point (steps S2, Yes), the control characteristic changing unit 6 determines whether the peak of the frequency transfer function is an anti-resonance point or an anti-resonance point (step S2, Yes). Step S3). For example, when there is a first peak point at which the gain value changes from decreasing to increasing, the control characteristic changing unit 6 determines that an antiresonance point exists in the frequency transfer function. When there is a second peak point at which the gain value changes from increasing to decreasing, the control characteristic changing unit 6 determines that the resonance point exists in the frequency transfer function. When a plurality of first peak points are present during scanning, the control characteristic changing unit 6 stores the plurality of first peak points, and determines that the lowest of them is an antiresonance point. When a plurality of second peak points are present during scanning, the control characteristic changing unit 6 stores the plurality of second peak points and determines that the largest of them is the resonance point.

周波数伝達関数のピークが共振点の場合(ステップS3,Yes)、制御特性変更部6は、切削力ベクトル演算部3で算出された、ピークの周波数における切削力のパワースペクトルが閾値以上であるか閾値未満であるかを判断する(ステップS4)。 When the peak of the frequency transfer function is a resonance point (steps S3, Yes), the control characteristic changing unit 6 determines whether the power spectrum of the cutting force at the peak frequency calculated by the cutting force vector calculation unit 3 is equal to or greater than the threshold value. It is determined whether the frequency is less than the threshold value (step S4).

切削力のパワースペクトルが閾値未満の場合(ステップS4,No)、制御特性変更部6はステップS8の処理を行う。切削力のパワースペクトルが閾値以上の場合(ステップS4,Yes)、制御特性変更部6は、ピークの周波数におけるサーボ制御部1の外乱応答補償器の制御ゲインを上げる(ステップS6)。 When the power spectrum of the cutting force is less than the threshold value (steps S4 and No), the control characteristic changing unit 6 performs the process of step S8. When the power spectrum of the cutting force is equal to or greater than the threshold value (step S4, Yes), the control characteristic changing unit 6 raises the control gain of the disturbance response compensator of the servo control unit 1 at the peak frequency (step S6).

ステップS3において、周波数伝達関数のピークが反共振点の場合(ステップS3,No)、制御特性変更部6は、切削力ベクトル演算部3で算出された、ピークの周波数における切削力のパワースペクトルが閾値以上であるか閾値未満であるかを判断する(ステップS5)。 In step S3, when the peak of the frequency transfer function is an antiresonance point (steps S3, No), the control characteristic changing unit 6 has the power spectrum of the cutting force at the peak frequency calculated by the cutting force vector calculation unit 3. It is determined whether it is equal to or greater than the threshold value or less than the threshold value (step S5).

切削力のパワースペクトルが閾値未満の場合(ステップS5,No)、制御特性変更部6はステップS8の処理を行う。切削力のパワースペクトルが閾値以上の場合(ステップS5,Yes)、制御特性変更部6は、ピークの周波数におけるサーボ制御部1の外乱応答補償器の制御ゲインを下げる(ステップS7)。制御特性変更部6は、走査対象の周波数の範囲内でピークの探索が完了するまでステップS1からS7までの処理を繰り返す(ステップS8)。 When the power spectrum of the cutting force is less than the threshold value (steps S5 and No), the control characteristic changing unit 6 performs the process of step S8. When the power spectrum of the cutting force is equal to or greater than the threshold value (steps S5 and Yes), the control characteristic changing unit 6 lowers the control gain of the disturbance response compensator of the servo control unit 1 at the peak frequency (step S7). The control characteristic changing unit 6 repeats the processes from steps S1 to S7 until the search for the peak is completed within the frequency range of the scanning target (step S8).

ステップS4,5において、周波数特性の変更を行うか否かを切削力のパワースペクトルにより判定されている理由は、切削力が小さい周波数の領域では、コヒーレントが低いため、周波数伝達関数のピークが正しく計算できていない可能性が高く、不要な周波数特性の変更を防ぐためである。ステップS4,5の処理では、制御特性を変更する駆動軸に対応する機械軸方向切削力ベクトルの各成分のパワースペクトルが用いられる。 The reason why whether or not to change the frequency characteristics is determined by the power spectrum of the cutting force in steps S4 and S5 is that the peak of the frequency transfer function is correct because the coherence is low in the frequency region where the cutting force is small. This is because there is a high possibility that the calculation has not been completed and unnecessary changes in frequency characteristics are prevented. In the processing of steps S4 and 5, the power spectrum of each component of the mechanical axial cutting force vector corresponding to the drive shaft whose control characteristics are changed is used.

なお、ステップS1の処理では、前述したコクアド線図に示される周波数伝達関数を用いてもよい。またステップS4,5の処理では、切削力のパワースペクトルの代わりに切削力の大きさを用いてもよい。 In the process of step S1, the frequency transfer function shown in the above-mentioned coquad diagram may be used. Further, in the processing of steps S4 and 5, the magnitude of the cutting force may be used instead of the power spectrum of the cutting force.

ステップS6における制御系の外乱応答を高める手段には、以下に説明する第一手段を利用してもよい。第一手段としては、外乱応答補償器の制御パラメータを変更することを例示できる。サーボ制御部1の外乱応答補償器は、一般的にP制御器、PI制御器などにより構成される。例えば、PI制御の場合、比例ゲインと積分ゲインとが外乱応答補償器の制御パラメータとなり、比例ゲインと積分ゲインとが共に、高い値に設定されることにより、外乱応答補償器の制御ゲインを上げることができる。しかし、この方法では、外乱応答補償器の全周波数域のゲインが上がり、制御系が不安定となる場合がある。 As a means for enhancing the disturbance response of the control system in step S6, the first means described below may be used. As a first means, it can be exemplified to change the control parameters of the disturbance response compensator. The disturbance response compensator of the servo control unit 1 is generally composed of a P controller, a PI controller, and the like. For example, in the case of PI control, the proportional gain and the integrated gain become control parameters of the disturbance response compensator, and both the proportional gain and the integrated gain are set to high values to increase the control gain of the disturbance response compensator. be able to. However, in this method, the gain in the entire frequency range of the disturbance response compensator increases, and the control system may become unstable.

ステップS6における制御系の外乱応答を高める手段には、以下に説明する第二手段を利用してもよい。第手段としては、狭帯域のゲインを上げるフィルタ、すなわち逆ノッチフィルタを、外乱応答補償器に付加することを例示できる。逆ノッチフィルタは例えば、(4)式の伝達関数で表される。sは微分を表すラプラス演算子、ωはゲインを上げる範囲の中心の角周波数である。Qはゲインを上げる周波数帯を決定する係数、aはゲインの上げる量を決定する係数である。共振点のピークの周波数としてωを設定することにより、外乱応答補償器の外乱が生じている周波数におけるゲインを上げることができる。
(s+ωs/Q+ω)/(s+as+ω)・・・(4)式
As a means for enhancing the disturbance response of the control system in step S6, the second means described below may be used. As a second means, it can be exemplified that a filter for increasing the gain in a narrow band, that is, a reverse notch filter is added to the disturbance response compensator. The inverse notch filter is represented by, for example, the transfer function of Eq. (4). s is the Laplace operator representing the derivative, and ω is the angular frequency at the center of the range for increasing the gain. Q is a coefficient for determining the frequency band for increasing the gain, and a is a coefficient for determining the amount for increasing the gain. By setting ω as the frequency of the peak of the resonance point, the gain at the frequency at which the disturbance of the disturbance response compensator is occurring can be increased.
(S 2 + ω s / Q + ω 2 ) / (s 2 + as + ω 2 ) ... Equation (4)

さらに、ステップS6における制御系の外乱応答を高める手段には、以下に説明する第三手段を利用してもよい。第三手段としては、特定の範囲の周波数のみ通過させ、かつ、特定の範囲の周波数以外の周波数は通過させないバンドパスフィルタ、又は特定の範囲の周波数のみ通過させ、かつ、特定の範囲の周波数以外の周波数は減衰させるバンドパスフィルタを用いることを例示できる。機械端位置201をバンドパスフィルタに通過させることにより、バンドパスフィルタを通過した機械端位置201に対する外乱応答補償器の制御パラメータが変更される。例えば、PI制御の場合、比例ゲインと積分ゲインとが外乱応答補償器の制御パラメータとなり、比例ゲインと積分ゲインとが共に、高い値に設定される。バンドパスフィルタを通過しない機械端位置に対しては、制御パラメータが変更されないため、第一手段より制御系が不安定となりにくい。 Further, as a means for enhancing the disturbance response of the control system in step S6, a third means described below may be used. As a third means, a bandpass filter that passes only frequencies in a specific range and does not pass frequencies other than the specific range, or passes only frequencies in a specific range and does not pass frequencies other than the specific range. It can be exemplified by using a bandpass filter that attenuates the frequency of. By passing the machine end position 201 through the bandpass filter, the control parameters of the disturbance response compensator for the machine end position 201 that has passed through the bandpass filter are changed. For example, in the case of PI control, the proportional gain and the integrated gain are control parameters of the disturbance response compensator, and both the proportional gain and the integrated gain are set to high values. Since the control parameters are not changed for the mechanical end position that does not pass through the bandpass filter, the control system is less likely to become unstable than the first means.

ステップS7においては、以下に説明する第一手段を利用してもよい。ステップS7の第一手段としては、外乱応答補償器の制御パラメータを変更することを例示できる。サーボ制御部1の外乱応答補償器がPI制御の場合、比例ゲインと積分ゲインとに低い値が設定される。 In step S7, the first means described below may be used. As the first means of step S7, changing the control parameter of the disturbance response compensator can be exemplified. When the disturbance response compensator of the servo control unit 1 is PI controlled, low values are set for the proportional gain and the integrated gain.

またステップS7においては、以下に説明する第二手段を利用してもよい。ステップS7の第二手段としては、狭帯域のゲインを下げるノッチフィルタを外乱応答補償器に付加することを例示できる。ノッチフィルタは、例えば(5)式で表される。sは微分を表すラプラス演算子、ωはゲインを下げる範囲の中心の角周波数である。Qはゲインを下げる周波数帯を決定する係数、bはゲインの下げる量を決定する係数である。共振点のピークの周波数としてωを設定することにより、外乱応答補償器の外乱が生じている周波数におけるゲインを下げることができる。
(s+bs+ω)/(s+ωs/Q+ω)・・・(5)式
Further, in step S7, the second means described below may be used. As the second means of step S7, it can be exemplified that a notch filter for lowering the gain in a narrow band is added to the disturbance response compensator. The notch filter is represented by, for example, Eq. (5). s is the Laplace operator representing the derivative, and ω is the angular frequency at the center of the range in which the gain is lowered. Q is a coefficient for determining the frequency band for lowering the gain, and b is a coefficient for determining the amount for lowering the gain. By setting ω as the frequency of the peak of the resonance point, the gain at the frequency at which the disturbance of the disturbance response compensator is generated can be lowered.
(S 2 + bs + ω 2 ) / (s 2 + ω s / Q + ω 2 ) ... Equation (5)

さらに、ステップS7には、以下に説明する第三手段を利用してもよい。第三手段としては、特定の範囲の周波数のみ通過させ、かつ、特定の範囲の周波数以外の周波数は通過させないバンドパスフィルタ、又は特定の範囲の周波数のみ通過させ、かつ、特定の範囲の周波数以外の周波数は減衰させるバンドパスフィルタを用いることを例示できる。機械端位置201をバンドパスフィルタに通過させることにより、バンドパスフィルタを通過した機械端位置201に対する外乱応答補償器の制御パラメータが変更される。例えば、PI制御の場合、比例ゲインと積分ゲインとが外乱応答補償器の制御パラメータとなり、比例ゲインと積分ゲインとが共に、低い値に設定される。 Further, in step S7, the third means described below may be used. As a third means, a bandpass filter that passes only frequencies in a specific range and does not pass frequencies other than the specific range, or passes only frequencies in a specific range and does not pass frequencies other than the specific range. It can be exemplified by using a bandpass filter that attenuates the frequency of. By passing the machine end position 201 through the bandpass filter, the control parameters of the disturbance response compensator for the machine end position 201 that has passed through the bandpass filter are changed. For example, in the case of PI control, the proportional gain and the integrated gain are control parameters of the disturbance response compensator, and both the proportional gain and the integrated gain are set to low values.

図8は切削力により、図1に示す機械端位置に外乱振動が生じている例を示す第1図である。図8には、時々刻々に変化する機械端位置を表す工具軌跡300が示され、さらに図8には、変化する各機械端位置における切削力ベクトル310が矢印で示される。なお、図8には、ピークの周波数における切削力のパワースペクトルが閾値以上のときの切削力ベクトル310のみ示される。 FIG. 8 is FIG. 1 showing an example in which disturbance vibration is generated at the machine end position shown in FIG. 1 due to a cutting force. FIG. 8 shows a tool locus 300 representing a machine end position that changes from moment to moment, and FIG. 8 shows an arrow indicating a cutting force vector 310 at each machine end position that changes. Note that FIG. 8 shows only the cutting force vector 310 when the power spectrum of the cutting force at the peak frequency is equal to or greater than the threshold value.

図9は切削力により、図1に示す機械端位置に外乱振動が生じている例を示す第2図である。図9には、位置指令200に対する切削力による外乱振動が拡大して表示されている。図9において、工具の送り方向は矢印Dで示される。図9に示すように切削中に工具の送り方向に垂直な方向に外乱振動が生じると、工具側面で加工される被加工物の加工面の面性状が悪化する。 FIG. 9 is FIG. 2 showing an example in which disturbance vibration is generated at the machine end position shown in FIG. 1 due to a cutting force. In FIG. 9, the disturbance vibration due to the cutting force with respect to the position command 200 is enlarged and displayed. In FIG. 9, the feed direction of the tool is indicated by an arrow D. As shown in FIG. 9, if disturbance vibration occurs in a direction perpendicular to the feed direction of the tool during cutting, the surface properties of the machined surface of the workpiece to be machined on the side surface of the tool deteriorate.

実施の形態1に係る数値制御装置100−1では、このように加工面の面性状が悪化する外乱振動が生じた場合、制御特性変更部6によって、サーボ制御部1の外乱応答特性が変更されることにより、外乱振動が低減される。さらに、数値制御装置100−1では、切削力ベクトル演算部3、機械特性同定部5、制御特性変更部6の処理が一定の周期で行われる。これは、加工位置の変化と、切削加工を行うことによる被加工物、工具及び機械の状態の変化とにより、切削力を、変更を行うべき周波数の変化に対応させるためである。 In the numerical control device 100-1 according to the first embodiment, when the disturbance vibration that deteriorates the surface properties of the machined surface occurs in this way, the disturbance response characteristic of the servo control unit 1 is changed by the control characteristic changing unit 6. As a result, disturbance vibration is reduced. Further, in the numerical control device 100-1, the processing of the cutting force vector calculation unit 3, the mechanical characteristic identification unit 5, and the control characteristic change unit 6 is performed at a constant cycle. This is to make the cutting force correspond to the change of the frequency to be changed by the change of the machining position and the change of the state of the workpiece, the tool and the machine due to the cutting.

以上のように、実施の形態1の数値制御装置100−1によれば、外乱抑制に用いる伝達関数を切削加工中に逐次計算できるため、切削加工中に計算される切削力ベクトル情報及び機械特性情報に基づいて、外乱応答補償器の周波数応答を変更することができる。これにより、切削力による機械系の位置外乱を抑制することができ、加工の精度向上が向上し、また工具の摩耗を抑制することができる。 As described above, according to the numerical control device 100-1 of the first embodiment, the transfer function used for disturbance suppression can be sequentially calculated during the cutting process, so that the cutting force vector information and the mechanical characteristics calculated during the cutting process can be calculated. The frequency response of the disturbance response compensator can be modified based on the information. As a result, it is possible to suppress the position disturbance of the mechanical system due to the cutting force, improve the accuracy of machining, and suppress the wear of the tool.

また、実施の形態1の数値制御装置100−1によれば、外乱抑制に用いる伝達関数を切削加工中に逐次計算できるため、切削条件が変わる度に周波数応答帯域を予め測定する必要がない。そのため、周波数応答帯域の事前測定が不要になり、また切削加工の際に適切な周波数応答帯域を変更する手間がなくなり、数値制御装置100−1のユーザの負担を軽減することができる。 Further, according to the numerical control device 100-1 of the first embodiment, since the transfer function used for the disturbance suppression can be sequentially calculated during the cutting process, it is not necessary to measure the frequency response band in advance every time the cutting conditions change. Therefore, the pre-measurement of the frequency response band becomes unnecessary, and the trouble of changing the appropriate frequency response band at the time of cutting is eliminated, and the burden on the user of the numerical control device 100-1 can be reduced.

また、実施の形態1の数値制御装置100−1では、切削力ベクトル演算部3、機械特性同定部5及び制御特性変更部6のそれぞれの処理動作が一定の周期で行われるため、加工位置の変化と、切削加工を行うことによる被加工物、工具及び機械の状態の変化とにより、切削力を対応させて変えることができる。 Further, in the numerical control device 100-1 of the first embodiment, each processing operation of the cutting force vector calculation unit 3, the machine characteristic identification unit 5, and the control characteristic change unit 6 is performed at a constant cycle, so that the machining position is determined. The cutting force can be changed in correspondence with the change and the change in the state of the workpiece, the tool and the machine due to the cutting process.

また実施の形態1の数値制御装置100−1のサーボ制御部1に2自由度制御が採用される場合、軌跡誤差を生じることなく、切削力による外乱を抑制することが可能となる。 Further, when the servo control unit 1 of the numerical control device 100-1 of the first embodiment adopts the two-degree-of-freedom control, it is possible to suppress the disturbance due to the cutting force without causing a locus error.

また実施の形態1の数値制御装置100−1によれば、従来技術のように周波数応答帯域を狭くする必要がないため、工具位置の指令軌跡への追従性を低下させることなく切削力変動によって生じる外乱を抑制できる。 Further, according to the numerical control device 100-1 of the first embodiment, it is not necessary to narrow the frequency response band as in the prior art, so that the cutting force fluctuates without deteriorating the followability of the tool position to the command locus. The generated disturbance can be suppressed.

実施の形態2.
図10は本発明の実施の形態2に係る数値制御装置の構成を示す図である。実施の形態2に係る数値制御装置100−2は、実施の形態1に係る数値制御装置100−1が備える構成に加えて、表示部10及び条件入力部11を備える。実施の形態2において、実施の形態1と同様の構成要素については実施の形態1と同一の名称及び記号を用いて説明する。
Embodiment 2.
FIG. 10 is a diagram showing a configuration of a numerical control device according to a second embodiment of the present invention. The numerical control device 100-2 according to the second embodiment includes a display unit 10 and a condition input unit 11 in addition to the configuration provided by the numerical control device 100-1 according to the first embodiment. In the second embodiment, the same components as those in the first embodiment will be described using the same names and symbols as those in the first embodiment.

条件入力部11では、数値制御装置100−2のユーザが、制御特性変更部6における制御特性を変更するための条件と、制御特性変更部6における制御特性の変更を許可するか否かを指定することができる。制御特性を変更するための条件とは、例えば、前述した制御特性変更部6のステップS4及びステップS5の処理動作における、パワースペクトルの閾値、切削力の大きさの閾値などである。 In the condition input unit 11, the user of the numerical control device 100-2 specifies a condition for changing the control characteristic in the control characteristic changing unit 6 and whether or not to allow the control characteristic to be changed in the control characteristic changing unit 6. can do. The conditions for changing the control characteristics include, for example, a threshold value of the power spectrum and a threshold value of the magnitude of the cutting force in the processing operations of step S4 and step S5 of the control characteristic changing unit 6 described above.

また工作機械100を利用して被加工物の加工を行うためのプログラムは、一般的にGコードと呼ばれる指令コードが用いられる。Gコードには、位置決め指令(G00)、直線補間切削指令(G01)、円弧補間切削指令(G02、G03)などの指令がある。一般的に、位置決め指令(G00)は非切削時に用いられる指令である。条件入力部11では、例えば、位置決め指令の場合に制御特性変更部6による制御特性の変更を実施させ、又は切削指令の場合に制御特性変更部6による制御特性の変更を実施させるなど、制御特性変更部6による制御特性の変更の許可条件をユーザが指定できる。 A command code generally called a G code is used as a program for processing a workpiece using the machine tool 100. The G code includes commands such as a positioning command (G00), a linear interpolation cutting command (G01), and an arc interpolation cutting command (G02, G03). Generally, the positioning command (G00) is a command used during non-cutting. In the condition input unit 11, for example, the control characteristic change unit 6 changes the control characteristic in the case of a positioning command, or the control characteristic change unit 6 changes the control characteristic in the case of a cutting command. The user can specify the permission condition for changing the control characteristic by the change unit 6.

また条件入力部11では、常時、制御特性の変更を許可するか否かを変更することも可能であり、例えば、制御特性の変更が許可され、かつ、制御特性変更を行う条件を満たす場合のみ、制御特性変更部6によるサーボ制御部1の制御特性の変更が行われる。 Further, the condition input unit 11 can always change whether or not to allow the change of the control characteristic. For example, only when the change of the control characteristic is permitted and the condition for changing the control characteristic is satisfied. , The control characteristic change unit 6 changes the control characteristic of the servo control unit 1.

表示部10では、機械特性同定部5の出力である機械特性情報207と、切削力ベクトル演算部3の出力である切削力ベクトル情報204の一部又は全部のデータとが表示される。数値制御装置100−2のユーザは、表示部10に表示されるこれらの情報を把握できる。例えば、表示部10には、図7に示される工具軌跡300の単位長さ(例えば1mm)が、特定の長さで表示される。また表示部10には、図7に示される切削力ベクトル310の単位力(1N)が、特定の長さで表示される。これにより、ユーザは、切削力の大きさと向きを確認することができ、切削加工中の状況を把握できると共に、制御特性変更部6において制御特性の変更を行うべきかどうかを容易に判定することができる。 On the display unit 10, the machine characteristic information 207 which is the output of the machine characteristic identification unit 5 and a part or all the data of the cutting force vector information 204 which is the output of the cutting force vector calculation unit 3 are displayed. The user of the numerical control device 100-2 can grasp the information displayed on the display unit 10. For example, on the display unit 10, the unit length (for example, 1 mm) of the tool locus 300 shown in FIG. 7 is displayed with a specific length. Further, the unit force (1N) of the cutting force vector 310 shown in FIG. 7 is displayed on the display unit 10 with a specific length. As a result, the user can confirm the magnitude and direction of the cutting force, grasp the situation during the cutting process, and easily determine whether or not the control characteristic change unit 6 should change the control characteristic. Can be done.

また、表示部10には、図5に示されるゲイン線図と図6に示される位相線図とのそれぞれの周波数の常用対数の単位量が、単位長さで表示される。これにより、ユーザは、どの周波数で振動が生じているかを確認でき、切削加工中の状況を把握できると共に、制御特性変更部6において制御特性の変更を行うべきかどうかを容易に判定することができる。 Further, the display unit 10 displays the unit amount of the common logarithm of each frequency of the gain diagram shown in FIG. 5 and the phase diagram shown in FIG. 6 in unit length. As a result, the user can confirm at which frequency the vibration is occurring, grasp the situation during cutting, and easily determine whether or not the control characteristic change unit 6 should change the control characteristic. it can.

以上のように、実施の形態2に係る数値制御装置100−2によれば、ユーザが条件入力部11より、制御特性の変更を行う条件と指定すると共に、制御特性の変更を許可するか否かとを指定することにより、ユーザの意図しない処理を防ぐことができる。さらに、表示部10に機械特性情報、切削力ベクトル情報などを表示させることにより、ユーザが制御特性の変更を行うべきかどうかを容易に判断することができる。 As described above, according to the numerical control device 100-2 according to the second embodiment, whether or not the user specifies the condition for changing the control characteristic from the condition input unit 11 and permits the change of the control characteristic. By specifying the heel, it is possible to prevent unintended processing by the user. Further, by displaying the machine characteristic information, the cutting force vector information, and the like on the display unit 10, it is possible to easily determine whether or not the user should change the control characteristics.

実施の形態3.
図11は本発明の実施の形態3に係る数値制御装置の構成を示す図である。位置の制御帯域に含まれる周波数の振動に対してノッチフィルタを利用した場合、指令軌跡への追従性が低下するため、望ましくない。そこで、制御帯域に含まれる低周波数の振動に関しては、切削力に起因する振動を異なる方法で抑制する必要がある。そこで実施の形態3に係る数値制御装置100−3では、図1に示される制御特性変更部6の代わりに、加算部12及び外乱補正量算出部13が用いられる。実施の形態3において、実施の形態1と同様の構成要素については実施の形態1と同一の名称及び記号を用いて説明する。
Embodiment 3.
FIG. 11 is a diagram showing a configuration of a numerical control device according to a third embodiment of the present invention. When a notch filter is used for vibration of a frequency included in the position control band, it is not desirable because the followability to the command locus is lowered. Therefore, regarding the low frequency vibration included in the control band, it is necessary to suppress the vibration caused by the cutting force by different methods. Therefore, in the numerical control device 100-3 according to the third embodiment, the addition unit 12 and the disturbance correction amount calculation unit 13 are used instead of the control characteristic changing unit 6 shown in FIG. In the third embodiment, the same components as those in the first embodiment will be described using the same names and symbols as those in the first embodiment.

外乱補正量算出部13は、機械特性情報207と切削力ベクトル情報204とに基づき、切削力に起因する外乱を推定し、推定した外乱から補正量を算出する。算出された補正量は加算部12により位置指令200へ加算され、加算部12の出力、すなわち補正量が加算された位置指令200は、サーボ制御部1及び模擬応答演算部4に入力される。数値制御装置100−3は、補正量が加算された位置指令200を用いて、工作機械100を制御することにより、工具位置の指令軌跡への追従性の低下を抑制しながら、切削力の変動に起因する外乱を抑制することで、加工精度の低下を抑制することができる。 The disturbance correction amount calculation unit 13 estimates the disturbance caused by the cutting force based on the mechanical characteristic information 207 and the cutting force vector information 204, and calculates the correction amount from the estimated disturbance. The calculated correction amount is added to the position command 200 by the addition unit 12, and the output of the addition unit 12, that is, the position command 200 to which the correction amount is added is input to the servo control unit 1 and the simulated response calculation unit 4. The numerical control device 100-3 controls the machine tool 100 by using the position command 200 to which the correction amount is added, thereby suppressing the decrease in the followability of the tool position to the command trajectory and fluctuating the cutting force. By suppressing the disturbance caused by the above, it is possible to suppress a decrease in processing accuracy.

図12は図2に示すモータ及び機械系を2慣性のモデルで近似したブロック図である。図12に示されるモデルにおいて、Kは機械系8の弾性要素のバネ定数であり、Cは機械系8のダンピング要素の粘性摩擦係数であり、Jmはモータ7のイナーシャであり、Jlは機械系8のイナーシャである。Jmの位置はモータ端位置を模擬しており、Jlの位置は機械端位置を模擬している。Tmは、トルク指令に従ってサーボモータ71が出力するトルクであり、xlは機械端位置を示す変数であり、xmはモータ端位置を示す変数であり、sは微分を表すラプラス演算子である。xm、xl及びfcを、ラプラス変換したものはX、X及びFcとする。FIG. 12 is a block diagram in which the motor and the mechanical system shown in FIG. 2 are approximated by a two-inertia model. In the model shown in FIG. 12, K is the spring constant of the elastic element of the mechanical system 8, C is the viscous friction coefficient of the damping element of the mechanical system 8, Jm is the inertia of the motor 7, and Jl is the mechanical system. 8 inertias. The position of Jm simulates the position of the end of the motor, and the position of Jl simulates the position of the end of the machine. Tm is the torque output by the servomotor 71 according to the torque command, xl is a variable indicating the machine end position, xm is a variable indicating the motor end position, and s is a Laplace operator representing differentiation. xm, the xl and fc, those Laplace transform to X M, X L and Fc.

切削力により弾性要素K、ダンピング要素Cなどが変化することにより、テーブル81に位置外乱が生じる。さらに、弾性要素K、ダンピング要素Cなどにより、切削力が伝達されるためモータ位置にも位置外乱が生じる。モータ7、機械系8が2慣性モデルでモデル化される場合、切削力から位置外乱までの伝達関数は、例えば(6)式で表される。
G(s)=1/(Jls+Cs+K)・・・(6)式
Positional disturbance occurs in the table 81 due to changes in the elastic element K, damping element C, and the like due to the cutting force. Further, since the cutting force is transmitted by the elastic element K, the damping element C, and the like, the position disturbance occurs at the motor position. When the motor 7 and the mechanical system 8 are modeled by a two-inertia model, the transfer function from the cutting force to the position disturbance is expressed by, for example, Eq. (6).
G (s) = 1 / (Jls 2 + Cs + K) ... Equation (6)

外乱補正量算出部13は、機械特性同定部5で同定された伝達関数に対して、例えば最小二乗法、スペクトル解析法、部分空間法などを用いて(6)式の各係数を同定する。なお、(6)式では、機械系8が2慣性モデルに近似されているが、これよりも低次又は高次のモデルによって伝達関数を表現してもよい。外乱補正量算出部13は、切削力出力部2により測定される切削力203に対して、上記のように同定された伝達関数型のフィルタを適用することにより、切削力により生じる外乱を算出することができる。ここで計算される外乱の計算値には、高周波成分も含まれている。高周波成分は制御できないため、外乱補正量算出部13では、適当な周波数のローパスフィルタを用いて高周波成分が除去される。外乱補正量算出部13は、制御周期で逐次外乱を算出する。外乱補正量算出部13は、算出された外乱の正負の値を反転させることにより、反転後の外乱を補正量とする。補正は加算部12に入力され、加算部12では、補正が位置指令200に加算される。補正量Cは(7)式及び(8)式により算出される。Fは機械軸切削力ベクトルである。Tはローパスフィルタの時定数である。
C=−GLPF(s)G(s)F・・・(7)式
LPF(s)=1/(Ts+1)・・・(8)式
The disturbance correction amount calculation unit 13 identifies each coefficient of the equation (6) with respect to the transfer function identified by the mechanical property identification unit 5 by using, for example, the least squares method, the spectrum analysis method, the subspace method, or the like. In Eq. (6), the mechanical system 8 is approximated to a two-inertia model, but the transfer function may be expressed by a lower-order or higher-order model. The disturbance correction amount calculation unit 13 calculates the disturbance generated by the cutting force by applying the transfer function type filter identified as described above to the cutting force 203 measured by the cutting force output unit 2. be able to. The calculated value of the disturbance calculated here also includes a high frequency component. Since the high frequency component cannot be controlled, the disturbance correction amount calculation unit 13 removes the high frequency component by using a low-pass filter having an appropriate frequency. The disturbance correction amount calculation unit 13 calculates the sequential disturbance in the control cycle. The disturbance correction amount calculation unit 13 reverses the calculated positive and negative values of the disturbance, and sets the disturbance after the inversion as the correction amount. The correction amount is input to the addition unit 12, and the correction amount is added to the position command 200 in the addition unit 12. The correction amount C is calculated by the equations (7) and (8). F is a machine axis cutting force vector. T is the time constant of the low-pass filter.
C = -G LPF (s) G (s) F ... (7) formula G LPF (s) = 1 / (Ts + 1) ... (8) formula

加工面の面性状に悪影響を及ぼすのは、工具の送り方向に対して垂直な方向の振動である。一方、工具の送り方向に対して補正を行った場合、工具の送り方向に対する工具の移動速度が一定ではなくなり、加工面の面性状に悪影響を及ぼす可能性がある。そこで、外乱補正量算出部13では、工具の送り方向に垂直な方向のみ補正が行われる。この場合、(7)式の機械軸方向切削力ベクトルFの代わりに、(3)式に示される垂直方向切削力ベクトルFnを用いればよい。 It is the vibration in the direction perpendicular to the feed direction of the tool that adversely affects the surface properties of the machined surface. On the other hand, when the correction is performed for the feed direction of the tool, the moving speed of the tool with respect to the feed direction of the tool becomes not constant, which may adversely affect the surface properties of the machined surface. Therefore, the disturbance correction amount calculation unit 13 corrects only in the direction perpendicular to the feed direction of the tool. In this case, the vertical cutting force vector Fn shown in the equation (3) may be used instead of the mechanical axial cutting force vector F in the equation (7).

以上のように実施の形態3に係る数値制御装置100−3によれば、指令軌跡への追従性を低下させることなく、位置制御の帯域に含まれる低周波な外乱振動を抑制できる。これにより、工具寿命が延び、また被加工物の加工精度を向上させることができる。また数値制御装置100−3では、実施の形態1に係る数値制御装置100−1と同様に、切削力ベクトル演算部3、機械特性同定部5及び制御特性変更部6のそれぞれの処理動作が一定の周期で行われるため、加工位置の変化と、切削加工を行うことによる被加工物、工具及び機械の状態の変化とにより、切削力を対応させて変えることができる。 As described above, according to the numerical control device 100-3 according to the third embodiment, it is possible to suppress low-frequency disturbance vibration included in the position control band without deteriorating the followability to the command locus. As a result, the tool life can be extended and the machining accuracy of the workpiece can be improved. Further, in the numerical control device 100-3, the processing operations of the cutting force vector calculation unit 3, the mechanical characteristic identification unit 5, and the control characteristic change unit 6 are constant, as in the numerical control device 100-1 according to the first embodiment. Since it is performed in the cycle of, the cutting force can be changed correspondingly by the change of the machining position and the change of the state of the workpiece, the tool and the machine due to the cutting.

図13は実施の形態1から3に係る数値制御装置のハードウェア構成例を示す図である。数値制御装置100−1,100−2,100−3は、演算処理を行うCPU(Central Processing Unit)51と、CPU51がワークエリアに用いるメモリ52と、プログラム、情報などを記憶可能な記憶装置53と、ユーザからの入力を受け付ける入力装置54と、表示装置55とを備える。メモリ52とは、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリー、EPROM(Erasable Programmable Read Only Memory)、EEPROM(登録商標)(Electrically Erasable Programmable Read Only Memory)といった不揮発性又は揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、DVD(Digital Versatile Disk)が該当する。入力装置54は、キーボード、マウスが例示され、表示装置55は、モニタ、ディスプレイが例示される。なお、入力装置54と表示装置55とが一体化されて、タッチパネルなどにより実現されてもよい。 FIG. 13 is a diagram showing a hardware configuration example of the numerical control device according to the first to third embodiments. The numerical control devices 100-1, 100-2, 100-3 are a CPU (Central Processing Unit) 51 that performs arithmetic processing, a memory 52 that the CPU 51 uses for a work area, and a storage device 53 that can store programs, information, and the like. And an input device 54 for receiving input from the user, and a display device 55. The memory 52 is a RAM (Random Access Memory), a ROM (Read Only Memory), a flash memory, an EPROM (Erasable Programmable Read Online Memory), an EEPROM (registered trademark), an Electrically volatile memory, or an EEPROM (registered trademark). This includes semiconductor memory, magnetic disk, flexible disk, optical disk, compact disk, mini disk, and DVD (Digital Versaille Disk). The input device 54 is exemplified by a keyboard and a mouse, and the display device 55 is exemplified by a monitor and a display. The input device 54 and the display device 55 may be integrated and realized by a touch panel or the like.

サーボ制御部1、切削力ベクトル演算部3、模擬応答演算部4、機械特性同定部5、制御特性変更部6、減算部9、外乱補正量算出部13及び加算部12は、CPU51が記憶装置53に格納されたプログラムを実行することにより実現される。 The CPU 51 stores the servo control unit 1, the cutting force vector calculation unit 3, the simulated response calculation unit 4, the mechanical characteristic identification unit 5, the control characteristic change unit 6, the subtraction unit 9, the disturbance correction amount calculation unit 13, and the addition unit 12. It is realized by executing the program stored in 53.

以上の実施の形態に示した構成は、本発明の内容の一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、本発明の要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。 The configuration shown in the above-described embodiment shows an example of the content of the present invention, can be combined with another known technique, and is one of the configurations without departing from the gist of the present invention. It is also possible to omit or change the part.

1 サーボ制御部、2 切削力出力部、3 切削力ベクトル演算部、4,4A 模擬応答演算部、5 機械特性同定部、6 制御特性変更部、7 モータ、8 機械系、9 減算部、10 表示部、11 条件入力部、12 加算部、13 外乱補正量算出部、41 位置制御模擬部、42 積分演算部、43 速度制御模擬部、44 積分演算器、45 比例定数演算部、51 CPU、52 メモリ、53 記憶装置、54 入力装置、55 表示装置、71 サーボモータ、71a ロータ、71b シャフト、72 モータ位置検出器、81 テーブル、82 ボールねじナット、83 ボールねじ、84 カップリング、85 機械端位置検出器、86 機械端位置検出器ヘッド、90 被加工物、100 工作機械、100−1,100−2,100−3 数値制御装置、101 工具、200 位置指令、201 機械端位置、202 トルク指令、203 切削力、204 切削力ベクトル情報、205 モデル位置、206 位置外乱、207 機械特性情報、208 モデル速度、209 モデルトルク、300 工具軌跡、310 切削力ベクトル。 1 Servo control unit, 2 Cutting force output unit, 3 Cutting force vector calculation unit, 4, 4A simulated response calculation unit, 5 Mechanical characteristic identification unit, 6 Control characteristic change unit, 7 Motor, 8 Mechanical system, 9 Subtraction unit, 10 Display unit, 11 condition input unit, 12 addition unit, 13 disturbance correction amount calculation unit, 41 position control simulation unit, 42 integration calculation unit, 43 speed control simulation unit, 44 integration calculation unit, 45 proportional constant calculation unit, 51 CPU, 52 Memory, 53 Storage Device, 54 Input Device, 55 Display Device, 71 Servo Motor, 71a Rotor, 71b Shaft, 72 Motor Position Detector, 81 Table, 82 Ball Screw Nut, 83 Ball Screw, 84 Coupling, 85 Machine End Position detector, 86 Machine end position detector head, 90 Work piece, 100 Machine, 100-1, 100-2, 100-3 Numerical controller, 101 Tool, 200 Position command, 201 Machine end position, 202 Torque Command, 203 Cutting Force, 204 Cutting Force Vector Information, 205 Model Position, 206 Position Disturbance, 207 Mechanical Characteristic Information, 208 Model Speed, 209 Model Torque, 300 Tool Trajectory, 310 Cutting Force Vector.

Claims (9)

被加工物を加工する工具により前記被加工物の加工を行う機械系を制御する数値制御装置であって、
位置指令に基づいて、前記被加工物が設けられる被駆動体の位置が前記位置指令に追従するように前記被駆動体に接続されるモータを駆動するサーボ制御部と、
前記工具と前記被加工物と間に生じる切削力を測定又は推定して出力する切削力出力部と、
前記被駆動体の位置と前記切削力とに基づき切削力ベクトルを演算する切削力ベクトル演算部と、
前記切削力に対する前記機械系の動特性を同定する機械特性同定部と、
前記切削力ベクトルと前記機械特性同定部で同定された前記機械系の動特性とに基づき、前記サーボ制御部の制御特性を変更する制御特性変更部と、
を備えることを特徴とする数値制御装置。
A numerical control device that controls a mechanical system that processes the workpiece with a tool that processes the workpiece.
A servo control unit that drives a motor connected to the driven body so that the position of the driven body on which the workpiece is provided follows the position command based on the position command.
A cutting force output portion to the cutting force measured or estimated and outputs generated between the workpiece and the tool,
A cutting force vector calculation unit that calculates a cutting force vector based on the position of the driven body and the cutting force.
A mechanical characteristic identification unit that identifies the dynamic characteristics of the mechanical system with respect to the cutting force,
A control characteristic changing unit that changes the control characteristics of the servo control unit based on the cutting force vector and the dynamic characteristics of the mechanical system identified by the mechanical characteristic identification unit.
A numerical control device characterized by being provided with.
前記機械特性同定部は、切削力から制御系の外乱までの周波数領域におけるゲインと位相により示される前記サーボ制御部の周波数応答特性を同定することを特徴とする請求項1に記載の数値制御装置。 The numerical control device according to claim 1, wherein the mechanical characteristic identification unit identifies the frequency response characteristic of the servo control unit indicated by a gain and a phase in a frequency region from a cutting force to a disturbance of a control system. .. 前記制御特性変更部は、前記機械特性同定部で同定された前記周波数応答特性において、前記周波数応答特性のピーク形状が互いに異なる共振点又は反共振点を探索することによって外乱要因を特定して、特定した前記外乱要因に応じて制御系の周波数特性を変更することを特徴とする請求項2に記載の数値制御装置。 The control characteristic changing unit identifies a disturbance factor by searching for a resonance point or an antiresonance point in which the peak shapes of the frequency response characteristics are different from each other in the frequency response characteristic identified by the mechanical characteristic identification unit. The numerical control device according to claim 2, wherein the frequency characteristics of the control system are changed according to the identified disturbance factor. 前記制御特性変更部は、前記周波数応答特性の変更の可否を前記切削力ベクトルにより判断することを特徴とする請求項2又は3に記載の数値制御装置。 The numerical control device according to claim 2 or 3, wherein the control characteristic changing unit determines whether or not the frequency response characteristic can be changed based on the cutting force vector. 前記制御特性変更部に対して、前記サーボ制御部の制御特性を変更する条件を入力する条件入力部を備えることを特徴とする請求項1から4の何れか一項に記載の数値制御装置。 The numerical control device according to any one of claims 1 to 4, wherein the control characteristic changing unit is provided with a condition input unit for inputting a condition for changing the control characteristic of the servo control unit. 被加工物を加工する工具により前記被加工物の加工を行う機械系を制御する数値制御装置であって、
位置指令に基づいて、前記被加工物が設けられる被駆動体の位置が前記位置指令に追従するように前記被駆動体に接続されるモータを駆動するサーボ制御部と、
前記工具と前記被加工物と間に生じる切削力を測定又は推定して出力する切削力出力部と、
前記被駆動体の位置と前記切削力とに基づき切削力ベクトルを演算する切削力ベクトル演算部と、
前記切削力に対する前記機械系の動特性を同定する機械特性同定部と、
前記切削力により生じる制御系の外乱を補正する補正量を算出する外乱補正量算出部と、
を備えることを特徴とする数値制御装置。
A numerical control device that controls a mechanical system that processes the workpiece with a tool that processes the workpiece.
A servo control unit that drives a motor connected to the driven body so that the position of the driven body on which the workpiece is provided follows the position command based on the position command.
A cutting force output portion to the cutting force measured or estimated and outputs generated between the workpiece and the tool,
A cutting force vector calculation unit that calculates a cutting force vector based on the position of the driven body and the cutting force.
A mechanical characteristic identification unit that identifies the dynamic characteristics of the mechanical system with respect to the cutting force,
A disturbance correction amount calculation unit that calculates a correction amount that corrects the disturbance of the control system caused by the cutting force, and a disturbance correction amount calculation unit.
A numerical control device characterized by being provided with.
前記機械特性同定部は、前記切削力から制御系の外乱までの周波数応答特性を同定し、
前記外乱補正量算出部は、前記切削力に、前記周波数応答特性を示す伝達関数型フィルタを適用することにより、外乱を補正するための補正量を演算することを特徴とする請求項6に記載の数値制御装置。
The mechanical characteristic identification unit identifies the frequency response characteristics from the cutting force to the disturbance of the control system.
The sixth aspect of claim 6, wherein the disturbance correction amount calculation unit calculates a correction amount for correcting disturbance by applying a transfer function type filter exhibiting the frequency response characteristic to the cutting force. Numerical control device.
前記外乱補正量算出部は、前記切削力ベクトルに基づき前記工具の送り方向に直交する方向の補正量を算出することを特徴とする請求項6又は7に記載の数値制御装置。 The numerical control device according to claim 6 or 7, wherein the disturbance correction amount calculation unit calculates a correction amount in a direction orthogonal to the feed direction of the tool based on the cutting force vector. 前記制御系の外乱は、前記位置指令、速度指令及びトルク指令の何れか1つに含まれる外乱であることを特徴とする請求項2、3、4、6、7又は8に記載の数値制御装置。 The numerical control according to claim 2, 3, 4, 6, 7 or 8, wherein the disturbance of the control system is a disturbance included in any one of the position command, the speed command and the torque command. apparatus.
JP2020525069A 2018-06-15 2018-06-15 Numerical control device Active JP6880322B2 (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2018/022975 WO2019239594A1 (en) 2018-06-15 2018-06-15 Numerical control device

Publications (2)

Publication Number Publication Date
JPWO2019239594A1 JPWO2019239594A1 (en) 2020-12-17
JP6880322B2 true JP6880322B2 (en) 2021-06-02

Family

ID=68842489

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020525069A Active JP6880322B2 (en) 2018-06-15 2018-06-15 Numerical control device

Country Status (4)

Country Link
JP (1) JP6880322B2 (en)
CN (1) CN112236729B (en)
DE (1) DE112018007736T5 (en)
WO (1) WO2019239594A1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023026411A1 (en) * 2021-08-25 2023-03-02 ファナック株式会社 State determination device and state determination method
CN115493734B (en) * 2022-09-29 2024-07-26 中国机械总院集团云南分院有限公司 Electric control system for reliability test of numerical control machine tool

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101952084B (en) * 2008-02-20 2012-08-08 兄弟工业株式会社 Backlash quantity detecting method for feeding drive device, and backlash quantity detecting device for the feeding drive device
JP5874261B2 (en) * 2011-09-14 2016-03-02 株式会社ジェイテクト Machining error calculation device, machining error calculation method, machining control device, and machining control method
JP5947113B2 (en) * 2012-06-01 2016-07-06 Dmg森精機株式会社 Workpiece processing method using machine tool
JP5739400B2 (en) * 2012-11-30 2015-06-24 ファナック株式会社 Servo control device having position correction function of driven body
JP6243260B2 (en) * 2013-12-16 2017-12-06 国立大学法人 東京大学 Spindle motor control device
CN103792888B (en) * 2014-01-26 2016-12-07 西安交通大学苏州研究院 Vibration suppression controller based on feed system positions feedback signal
JP6399866B2 (en) * 2014-09-03 2018-10-03 オークマ株式会社 Servo control device
JP6439542B2 (en) * 2015-03-30 2018-12-19 ブラザー工業株式会社 Numerical control device and control method
CN106489105B (en) * 2015-06-18 2018-06-22 三菱电机株式会社 Control parameter adjusting apparatus
JP6568035B2 (en) * 2016-10-04 2019-08-28 ファナック株式会社 Servo motor control device, servo motor control method, and servo motor control program
WO2018073873A1 (en) * 2016-10-17 2018-04-26 三菱電機株式会社 Servo control device
JP6370866B2 (en) * 2016-12-01 2018-08-08 ファナック株式会社 Servo motor control device, servo motor control method, and servo motor control program

Also Published As

Publication number Publication date
WO2019239594A1 (en) 2019-12-19
CN112236729A (en) 2021-01-15
JPWO2019239594A1 (en) 2020-12-17
CN112236729B (en) 2024-03-15
DE112018007736T5 (en) 2021-03-04

Similar Documents

Publication Publication Date Title
JP4581096B2 (en) Friction compensation method, friction compensator, and motor control device
JP5566469B2 (en) Numerical control method
JP6567205B1 (en) Machine learning device, correction parameter adjusting device, and machine learning method
JP3703664B2 (en) Backlash correction device
JP6804657B2 (en) Numerical control system and motor control device
JP6592143B2 (en) Electric motor control device
JP6880322B2 (en) Numerical control device
JP2018077778A (en) Numerical control system
US9921568B2 (en) Trajectory measuring device, numerical control device, and trajectory measuring method
JP2021057030A (en) Control system
JP6903485B2 (en) Vibration damping device and processing machine
JP5836206B2 (en) Servo control device
JP7034383B2 (en) Servo controller
WO2018077593A1 (en) Method for determining a deadzone angle of a backlash in a mechanical drive-train system, method for controlling a drive motor controller as well as drive-train system
JP6333495B1 (en) Servo control device
JP5407435B2 (en) Motor control device
JP4576530B2 (en) Servo gain calculation method, servo gain calculation program, and servo gain calculation device
JP2003157114A (en) Method and device for lost motion correction
JP2012056005A (en) Method and apparatus for setting parameter
WO2022186051A1 (en) Adjustment assistance device, control system, and adjustment assistance method
JP2007179364A (en) Apparatus for correcting amount of displacement of main shaft
WO2023153446A1 (en) Proposal device, proposal system, proposal method, and program
Chong et al. AR-CM NCTF Control for Precision Positioning Systems--Concept and Results
JP2023051318A (en) Numerical control device and control method
JP2023122980A (en) Numerical control device and control method

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200601

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200909

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: 20210406

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210430

R150 Certificate of patent or registration of utility model

Ref document number: 6880322

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250