MACHINING SIMULATION DEVICE AND MACHINING SIMULATION METHOD

Abstract

The present disclosure provides a machining simulation device capable of obtaining a machine-tool transfer function without requiring data collection and expert knowledge. The machining simulation device includes: a transmission characteristic generation unit that generates a transmission characteristic of a machine tool; and a simulation execution unit that simulates behavior of the machine tool by using the transmission characteristic. The transmission characteristic generation unit includes: a control information acquisition unit that acquires motor information and a control parameter of a motor control system of the machine tool; a motor characteristic calculation unit that calculates a motor characteristics based on the motor information; a control characteristic calculation unit that calculates a control characteristic of the motor control system on the basis of the control parameter; and a transmission characteristic calculation unit that calculates the transmission characteristic satisfying a predetermined requirement on the basis of the motor characteristic and the control characteristic.

Claims

1. A machining simulation device, comprising: a transfer characteristic generation unit that generates a transfer function representing a transfer characteristic of a machine tool; and a simulation execution unit that simulates behavior of the machine tool using the transfer function, wherein the transfer characteristic generation unit includes: a control information acquisition unit that acquires control information, including motor information of the machine tool and a control parameter of a motor control system of the machine tool, from a storage unit; a motor characteristic calculation unit that calculates a motor characteristic, based on the motor information; a control characteristic calculation unit that calculates a control characteristic of the motor control system, based on the control parameter; and a transfer characteristic calculation unit that calculates the transfer function satisfying a predetermined requirement, based on at least one of the motor characteristic or the control characteristic, wherein the transfer function calculated by transfer characteristic calculation unit includes the at least one of the motor characteristic or the control characteristic, wherein the predetermined requirement includes at least one of a response frequency of position control, a response frequency of velocity control, a natural frequency of the machine tool, or a resonance frequency between a driving unit and a driven unit of the machine tool and is given by a user.

2. The machining simulation device according to claim 1, wherein the motor information includes at least one of motor inertia, an inertia ratio, or a spring constant, all caused by a motor, and the control parameter includes at least one of position control proportional gain, velocity loop gain, velocity control integral gain, or velocity control proportional gain.

3. The machining simulation device according to claim 2, wherein the motor characteristic calculation unit calculates the motor characteristic using a rigid-body system model including the motor inertia and the inertia ratio, or a two-inertia system model including the motor inertia, the inertia ratio, and the spring constant.

4. The machining simulation device according to claim 2, wherein the control characteristic calculation unit calculates the control characteristic using the position control proportional gain.

5. The machining simulation device according to claim 2, wherein the control characteristic calculation unit calculates the control characteristic using the velocity loop gain, or the velocity loop gain, the velocity control integral gain, and the velocity control proportional gain.

6. The machining simulation device according to claim 1, wherein the motor characteristic includes at least one of motor inertia, an inertia ratio, or a spring constant, the control characteristic includes at least one of position control proportional gain, velocity loop gain, velocity control integral gain, or velocity control proportional gain, and the transfer characteristic calculation unit changes at least one of the motor inertia, the inertia ratio, the spring constant, the position control proportional gain, the velocity loop gain, the velocity control integral gain, or the velocity control proportional gain such that the transfer function satisfies the predetermined requirement.

7. The machining simulation device according to claim 1, wherein the transfer characteristic calculation unit adds load inertia to motor inertia, or connects the motor inertia and the load inertia with a spring element, such that the transfer function satisfies the predetermined requirement.

8. A machining simulation method causing a computer to execute processing comprising: calculating a motor characteristic, based on motor information of a machine tool; calculating a control characteristic of a motor control system, based on a control parameter of the motor control system of the machine tool; calculating a transfer function representing a transfer characteristic of the machine tool satisfying a predetermined requirement, based on at least one of the motor characteristic or the control characteristic; and simulating behavior of the machine tool using the transfer function, wherein the transfer function includes the at least one of the motor characteristic or the control characteristic, wherein the predetermined requirement includes at least one of a response frequency of position control, a response frequency of velocity control, a natural frequency of the machine tool, or a resonance frequency between a driving unit and a driven unit of the machine tool and is given by a user.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0014] FIG. 1 is a configuration diagram illustrating an example of the configuration of a machining simulation system for a machine tool according to one embodiment of the present invention.

[0015] FIG. 2 is a diagram illustrating the driving unit and the driven unit of a machine tool using a rigid body model.

[0016] FIG. 3 is a diagram illustrating the driving unit and the driven unit of a machine tool using a two-inertia system model.

[0017] FIG. 4 is a block diagram illustrating the configuration of a position control loop when the motor control system, the driving unit, and the driven unit form a position control loop.

[0018] FIG. 5 is a block diagram illustrating the configuration of a simplified position control loop.

[0019] FIG. 6 is a block diagram illustrating the configuration of a velocity control loop when the motor control system, the driving unit, and the driven unit form a velocity control loop.

[0020] FIG. 7 is a characteristic diagram schematically illustrating the frequency response gain of a transfer function.

[0021] FIG. 8 is a characteristic diagram illustrating an example of the frequency response gain of velocity PI control.

[0022] FIG. 9 is a characteristic diagram illustrating an example of the frequency response gain of velocity P control.

[0023] FIG. 10 is a diagram illustrating an example of a user interface displayed on a display screen of a personal computer.

[0024] FIG. 11 is a flowchart illustrating the operation of the machining simulation device.

PREFERRED MODE FOR CARRYING OUT THE INVENTION

[0025] Hereinafter, an embodiment of the present disclosure will be described in detail with reference to the drawings. FIG. 1 is a configuration diagram illustrating an example of the configuration of a machining simulation system for a machine tool according to one embodiment of the present invention.

[0026] As illustrated in FIG. 1, the machining simulation system 10 for a machine tool includes a machine tool control device 100, and a machining simulation device 200.

[0027] The machine tool control device 100 controls the movement of the feed shaft and the rotation of the main shaft of the machine tool, based on a machining program. The control device 100 includes a motor, which serves as the driving unit of the machine tool, a driven unit, and a motor control system that controls the motor. The control device 100 stores motor information of the machine tool and a control parameter of the motor control system in a storage unit 101, The motor information and the control parameter will be described later. The motor information of the machine tool and the control parameter may be stored in a storage unit separate from the control device 100. The storage unit may be provided within the machining simulation device 200. In a case where the storage unit is provided separately from the control device 100, the machining simulation system 10 may not include the machine tool control device 100.

[0028] The machining simulation device 200 includes a transfer characteristic generation unit 210 and a simulation execution unit 220. The transfer characteristic generation unit 210 calculates the transfer function of the machine tool using the motor information and the control parameter obtained from the storage unit 101 of the machine tool control device 100. The simulation execution unit 220 uses the calculated transfer function to simulate the control executed by the machine tool control device 100, based on the machining program, the behavior of the motor as the driving unit and the driven unit, and the feedback control of the control device 100, based on the position information of the driving unit and the driven unit (for example, the position control loop described later as illustrated in FIG. 4), and outputs the position information for each shaft as simulation results.

[0029] The transfer characteristic generation unit 210 includes a control information acquisition unit 211, a motor characteristic calculation unit 212, a control characteristic calculation unit 213, and a transfer characteristic calculation unit 214. Each configuration of the transfer characteristic generation unit 210 will be described below.

(Control Information Acquisition Unit 211)

[0030] The control information acquisition unit 211 acquires control information, including motor information and a control parameter, from the storage unit 101 of the machine tool control device 100. The motor information includes, for example, at least one of motor inertia, the inertia ratio, or the spring constant, all caused by the motor. The inertia ratio refers to the ratio of load inertia to motor inertia. Load inertia is also referred to as load inertia. Motor inertia may also include the inertia of the reducer and the ball screw. The control parameter includes, for example, at least one of the position control proportional gain K.sub.P of the position control unit included in the motor control system of the control device 100, the velocity loop gain K.sub.V, the velocity control integral gain K.sub.1 of the velocity control unit, or the velocity control proportional gain K.sub.2 of the velocity control unit. The storage unit 101 may store motor inertia as part of the specification or motor characteristic database. In this case, the control information acquisition unit 211 acquires information that allows for identifying the motor, such as the motor model, from the storage unit 101 and references the specification or motor characteristic database to acquire the motor inertia.

(Motor Characteristic Calculation Unit 212)

[0031] The motor characteristic calculation unit 212 calculates a motor characteristic, based on motor information. When the mechanical model of the motor, which serves as the driving unit of the machine tool, and the driven unit is represented by the rigid-body system model illustrated in FIG. 2, and the motor information includes motor inertia J.sub.M and inertia ratio R, the motor characteristic M.sub.1 is given by the transfer function in Equation 1 (hereinafter referred to as Equation 1). Here, s represents the Laplace transform variable.

[00001]$\begin{array}{cc}{M}_1\left(s\right)=\frac{1}{\left(1+R\right)J_{M}s}& \left[\mathrm{Equation}1\right]\end{array}$

[0032] When the mechanical model of the motor, which serves as the driving unit of the machine tool, and the driven unit is represented by the two-inertia system model illustrated in FIG. 3, and the motor information includes motor inertia J.sub.M, natural frequency .sub.O, and resonance frequency .sub.P, the motor characteristic M.sub.2 is given by the transfer function in Equation 2 (hereinafter referred to as Equation 2).

[00002]$\begin{array}{cc}{M}_2\left(s\right)=\frac{1}{J_{M}s}\frac{\left(s^2+{}_O^2\right)}{\left(s^2+{}_P^2\right)}& \left[\mathrm{Equation}2\right]\end{array}$

[0033] The natural frequency .sub.O and resonance frequency .sub.P in Equation 2 are given by Equation 3 (hereinafter referred to as Equation 3) using the motor inertia J.sub.M, the inertia ratio R, and the spring constant K.sub.S of the spring element between the motor inertia and the load inertia.

[00003]$\begin{array}{cc}{}_O=\sqrt{\frac{K_S}{{\mathrm{RJ}}_M}},{}_P=\sqrt{\frac{\left(1+R\right)K_S}{{\mathrm{RJ}}_M}}& \left[\mathrm{Equation}3\right]\end{array}$

The natural frequency .sub.O of the two-inertia system model is the natural frequency of free vibration of the driven unit when the driving unit is fixed, and may be referred to as the anti-resonance frequency in some cases. The resonance frequency .sub.P of the two-inertia system model is the frequency at which the driving unit and the driven unit vibrate in opposite phases.

[0034] The transfer function of the motor inertia and load inertia connected by a spring element is disclosed, for example, in A Study on Low Frequency Vibration Suppression Control by Two-Mass System Model for Feed Axes of NC Machine Tools, Yasusuke Iwashita et al., Journal of the Japan Society for Precision Engineering, Vol. 82, No. 8, 2016.

(Control Characteristic Calculation Unit 213)

[0035] The control characteristic calculation unit 213 calculates the control characteristic of the motor control system, based on a control parameter. In a case where the motor control system, the driving unit, and the driven unit included in the control device 100 configure a position control loop, the position control loop is represented by the block diagram illustrated in FIG. 4. As illustrated in FIG. 4, the position control loop can be represented by a position control unit 111, a velocity control unit 112, a mechanical model unit 113 consisting of the driving unit and the driven unit, and an integrator 114. In many cases, the position control loop is configured as a dual-loop system with an internal velocity control loop, as illustrated in FIG. 4. In FIG. 4, y represents the position command, y represents the actual position, r represents the velocity command, u represents the actuation quantity, and w represents the actual velocity. The actuation quantity u is, for example, a torque command; however, for a servo motor, the relationship of torque=K.sub.Tcurrent holds via the torque constant K.sub.T, thus the actuation quantity u may be considered as current.

[0036] Since the response frequency of the position control loop must be lower than the response frequency of the internal velocity control loop, the position control loop can be simplified to the block diagram configuration illustrated in FIG. 4, where the closed-loop transfer function from the velocity command r to the actual velocity w is approximated as 1. In a case where the position control loop is configured as illustrated in FIG. 5, the position control loop consists of the position control unit 111 and the integrator 114, and the motor control system corresponds to the position control unit 111. The control characteristic C.sub.P of the position control unit 111, which serves as the motor control system, is expressed as C.sub.P=K.sub.P, and the control parameter is the position control proportional gain K.sub.P.

[0037] In a case where the motor control system, the driving unit, and the driven unit included in the control device 100 configure a velocity control loop, the velocity control loop is represented by the block diagram illustrated in FIG. 6. The velocity control loop can be represented by a velocity control unit 112 and a mechanical model unit 113 consisting of the driving unit and the driven unit. The motor control system corresponds to the velocity control unit 112.

[0038] The control characteristic C.sub.V of the velocity control unit 112, which serves as the motor control system, in the case of the velocity PI control, is given by Equation 4 (hereinafter referred to as Equation 4). In Equation 4, K.sub.V is the velocity loop gain, K.sub.1 is the velocity control integral gain, and K.sub.2 is the velocity control proportional gain. The control parameters are the velocity loop gain K.sub.V, the velocity control integral gain K.sub.1, and the velocity control proportional gain K.sub.2.

[00004]$\begin{array}{cc}{C}_V\left(s\right)=K_V\frac{\left(K_1+K_{2}s\right)}{s}& \left[\mathrm{Equation}4\right]\end{array}$

In the case of the velocity P control, the control characteristic C.sub.V of the velocity control unit is obtained by setting K.sub.1=0 and K.sub.2=1 in Equation 4, resulting in C.sub.V=K.sub.V. The control parameter is the velocity loop gain K.sub.V.

(Transfer Characteristic Calculation Unit 214)

[0039] The transfer characteristic calculation unit 214 calculates the transfer characteristic of the machine tool by changing at least one of the motor characteristic calculated by the motor characteristic calculation unit 212 or the control characteristic calculated by the control characteristic calculation unit 213 to satisfy a predetermined requirement. The predetermined requirement includes at least one of the response frequency of position control, the response frequency of velocity control, the natural frequency of the machine tool, or the resonance frequency between the driving unit and the driven unit of the machine tool. The predetermined requirement will be described later. The value of the predetermined requirement may be provided by the user, or may be acquired by the transfer characteristic calculation unit 214 from the storage unit 101 of the control device 100, a storage unit provided outside the control device 100, or a storage unit provided within the machining simulation device.

[0040] First, the transfer characteristic of the machine tool will be described. The closed-loop transfer function G.sub.O is expressed as G.sub.C=G.sub.O/(1+G.sub.O), where G.sub.O is the open-loop transfer function.

[0041] The transfer characteristic calculation unit 214 calculates the closed-loop transfer function G.sub.PC(s) for the position control loop and the closed-loop transfer function G.sub.VC(s) for the velocity control loop, as follows. The closed-loop transfer function G.sub.PC(s) for the position control loop, where the open-loop transfer function G.sub.PC(s) of the position control loop is G.sub.PC(s)=K.sub.P/s as illustrated in FIG. 5, is given by Equation 5 (hereinafter referred to as Equation 5).

[00005]$\begin{array}{cc}{G}_{\mathrm{PC}}\left(s\right)=\frac{G_{\mathrm{PO}}\left(s\right)}{\left(1+G_{\mathrm{PO}}\left(s\right)\right)}=\frac{K_P}{\left(s+K_P\right)}& \left[\mathrm{Equation}5\right]\end{array}$

[0042] The closed-loop transfer function G.sub.VC(s) for the velocity control loop, in the case of the velocity PI control, is given by Equation 7 (hereinafter referred to as Equation 7), where the open-loop transfer function G.sub.VC(s) is given by Equation 6 (hereinafter referred to as Equation 6). In Equations 6 and 7, M represents the motor characteristic, where M=M.sub.1 when the mechanical model is represented by the rigid-body system model, and M=M.sub.2 when the mechanical model is represented by the two-inertia system model.

[00006]$\begin{array}{cc}{G}_{\mathrm{VO}}\left(s\right)=C_V\left(s\right)M=K_V\frac{\left(K_1+K_{2}s\right)}{s}M& \left[\mathrm{Equation}6\right]\end{array}$$\begin{array}{cc}{G}_{\mathrm{VC}}\left(s\right)=\frac{K_V\frac{K_1+K_{2}s}{s}M}{1+K_V\frac{K_1+K_{2}s}{s}M}& \left[\mathrm{Equation}7\right]\end{array}$

[0043] In the case of the velocity P control, the closed-loop transfer function G.sub.VC(s) for the velocity control loop is obtained by setting K.sub.1=0 and K.sub.2=1 in Equation 7, resulting in Equation 8 (hereinafter referred to as Equation 8).

[00007]$\begin{array}{cc}{G}_{\mathrm{VC}}\left(s\right)=\frac{K_{V}M}{1+K_{V}M}& \left[\mathrm{Equation}8\right]\end{array}$

[0044] Next, the predetermined requirement will be described.

(Case where Predetermined Requirement is Response Frequency of Position Control)

[0045] When the transfer function is a rational function of s, and the command a is a sine wave with a frequency , the control variable b will also become a sine wave with the frequency . In this case, the amplitude ratio between the command a and the control variable b is referred to as the frequency response gain at the frequency . The frequency response gain of the closed-loop transfer function G.sub.C(s) can be calculated by replacing s with j in Equation 9 (hereinafter referred to as Equation 9). The unit is decibels (dB).

[00008]$\begin{array}{cc}20\log .Math.G_C\left(j\right).Math.& \left[\mathrm{Equation}9\right]\end{array}$

[0046] In general, the frequency response gain of a closed-loop transfer function is near 0 dB when the frequency response gain of the open-loop transfer function is 0 dB or higher. The case where the frequency response gain of the closed-loop transfer function is near 0 dB (amplitude ratio of 1) means that the control variable b follows the command a. Therefore, the frequency at which the frequency response gain of the closed-loop transfer function intersects 0 dB to 3 dB (amplitude ratio of 1/square root of 2) or becomes local maximum is referred to as the response frequency of the control. The frequency at which the frequency response gain of the open-loop transfer function intersects 0 dB is referred to as the response frequency of the control.

[0047] By replacing s with j in the transfer function G.sub.PC(s) of Equation 5, Equation 10 (hereinafter referred to as Equation 10) is obtained, and the response frequency =K.sub.P becomes the response frequency of position control where 20 log|G(j)|=3 dB.

[00009]$\begin{array}{cc}.Math.G_{\mathrm{PC}}\left(j\right).Math.=\frac{K_P}{\sqrt{{}^2+K_P^2}}& \left[\mathrm{Equation}10\right]\end{array}$

(Case where Predetermined Requirement is Response Frequency of Velocity Control)

[0048] When the control characteristic C is the velocity P control and the motor characteristic M is based on the two-inertia system model, the open-loop transfer function of the velocity control loop is given by Equation 11 (hereinafter referred to as Equation 11).

[00010]$\begin{array}{cc}{G}_{\mathrm{VO}}\left(s\right)=\frac{K_V}{J_{M}s}\frac{\left(s^2+{}_O^2\right)}{\left(s^2+{}_P^2\right)}& \left[\mathrm{Equation}11\right]\end{array}$

[0049] The frequency response gain of the transfer function in Equation 11 is schematically illustrated by the gain plot in FIG. 7. In FIG. 7, the horizontal axis is on a logarithmic scale of frequency. The point where the frequency response gain intersects 0 dB is the response frequency, the point where the frequency response gain reaches the local minimum value is the natural frequency, and the point where the frequency response gain reaches the local maximum value is the resonance frequency.

[0050] As illustrated in FIG. 7, in general, the relationship between the frequencies is as follows: resonance frequency>natural frequency>response frequency of velocity control. This relationship also holds for the velocity PI control. In a case where the frequency is no greater than the response frequency of velocity control, the relationship is as follows: resonance frequency .sub.P>natural frequency .sub.O>frequency (.sub.P>.sub.O>>), resulting in .sub.P.sup.2>.sub.O.sup.2, thus, the term s.sup.2 in Equation 2 may be ignored. By ignoring the term s.sup.2, Equation 2 becomes equivalent to Equation 1. Therefore, when considering the response frequency of velocity control, Equation 1 can be used, representing the rigid-body system model as the motor characteristic M.

(Case of Velocity PI Control)

[0051] Considering the velocity PI control as the control characteristic C and using the rigid-body system model as the motor characteristic M, Equation 7 is expressed as Equation 12 (hereinafter referred to as Equation 12),

[00011]$\begin{array}{cc}{G}_{\mathrm{VC}}\left(s\right)=\frac{\left(K_V\frac{K_1+K_{2}s}{s}\frac{1}{\left(1+R\right)J_{M}s}\right)}{\left(1+K_V\frac{K_1+K_{2}s}{s}\frac{1}{\left(1+R\right)J_{M}s}\right)}=\frac{K_V{K}_{2}s+K_V{K}_1}{\left(1+R\right)J_M{s}^2+K_V{K}_{2}s+K_V{K}_1}& \left[\mathrm{Equation}12\right]\end{array}$

[0052] By replacing s with j in Equation 12, Equation 12 is expressed as Equation 13 (hereinafter referred to as Equation 13).

[00012]$\begin{array}{cc}{G}_{\mathrm{VC}}\left(j\right)=\frac{K_V{K}_{2}j+K_V{K}_1}{K_V{K}_{2}j+K_V{K}_1-\left(1+R\right)J_M{}^2}& \left[\mathrm{Equation}13\right]\end{array}$

[0053] The frequency response gain is given by Equation 14 (hereinafter referred to as Equation 14) using Equation 13.

[00013]$\begin{array}{cc}20{\log }_{10}\frac{.Math.K_V{K}_{2}j+K_V{K}_1.Math.}{.Math.K_V{K}_{2}j+K_V{K}_1-\left(1+R\right)J_M{}^2.Math.}=10{\log }_{10}\frac{{\left(K_V{K}_2\right)}^2+{\left(K_V{K}_1\right)}^2}{{\left(K_V{K}_2\right)}^2+{\left\{K_V{K}_1-\left(1+R\right)J_M{}^2\right\}}^2}& \left[\mathrm{Equation}14\right]\end{array}$

[0054] In Equation 14, since the frequency response gain remains 0 dB or higher up to the values of Equation 15 (hereinafter referred to as Equation 15) at which the absolute values of the denominator and numerator are equal, the frequency of Equation 15 is used as the response frequency of the velocity PI control.

[00014]$\begin{array}{cc}=\sqrt{\frac{2K_V{K}_1}{\left(1+R\right)J_M}}& \left[\mathrm{Equation}15\right]\end{array}$

[0055] FIG. 8 is a characteristic diagram illustrating an example of the frequency response gain for the velocity PI control.

(Case of Velocity P Control)

[0056] Considering the velocity P control as the control characteristic C and using the rigid-body system model as the motor characteristic M, the frequency response gain is obtained by setting K.sub.1=0 and K.sub.2=1 in Equation 14, resulting in Equation 16 (hereinafter referred to as Equation 16).

[00015]$\begin{array}{cc}20{\log }_{10}\frac{.Math.K_{V}j.Math.}{.Math.K_{V}j-\left(1+R\right)J_M{}^2.Math.}=20{\log }_{10}\frac{K_V}{.Math.K_V+j\left(1+R\right)J_M.Math.}& \left[\mathrm{Equation}16\right]\end{array}$

[0057] In Equation 16, the value of in Equation 17 (hereinafter referred to as Equation 17), where 20 log|G(j)|=3 dB, is the response frequency of the velocity P control,

[00016]$\begin{array}{cc}=\frac{K_V}{\left(1+R\right)J_M}& \left[\mathrm{Equation}17\right]\end{array}$

[0058] FIG. 9 is a characteristic diagram illustrating an example of the frequency response gain for the velocity P control.

(Case where Predetermined Requirement is Natural Frequency or Resonance Frequency of Machine Tool)

[0059] As mentioned earlier, when the motor, which serves as the driving unit of the machine tool and the driven unit, are represented by the two-inertia system model illustrated in FIG. 3, and the motor information includes the motor inertia J.sub.M, the natural frequency , and the resonance frequency .sub.P, the motor characteristic M.sub.2 is given by Equation 2. In Equation 2, the natural frequency .sub.O and the resonance frequency op are given by Equation 3 using the motor inertia J.sub.M, the inertia ratio R, and the spring constant K.sub.S of the spring element between the motor inertia and the load inertia. The natural frequency .sub.O and the resonance frequency .sub.P represent the natural frequency and the resonance frequency of the machine tool.

[0060] The following describes the method by which the transfer characteristic calculation unit 214 calculates the transfer characteristic of the machine tool to satisfy the predetermined requirement. The transfer characteristic generation unit 210 included in the machining simulation device 200 operates as application software on a personal computer, and includes the user interface illustrated in FIG. 10.

[0061] FIG. 10 is a diagram illustrating an example of a user interface displayed on a display screen of a personal computer. As illustrated in FIG. 10, the display screen displays the control parameters such as the position control proportional gain K.sub.P, the velocity loop gain K.sub.V, the velocity control integral gain K.sub.1, and the velocity control proportional gain K.sub.2, as well as the motor information such as the motor inertia J.sub.M, the inertia ratio R, and the spring constant K.sub.S, which are read from the control information acquired by the control information acquisition unit 211. Specifically, as illustrated in FIG. 10, the position control proportional gain K.sub.P is 30, the velocity loop gain K.sub.V is 1, the velocity control integral gain K.sub.1 is 1089, the velocity control proportional gain K.sub.2 is 10.164, the motor inertia J.sub.M is 0.022 [kgm.sup.2], the inertia ratio R is 1.2, and the spring constant K.sub.S is 6878 [Nm], which are displayed on the screen. These control parameters and motor information are changed based on the predetermined requirements input by the user. On the left side of the screen, a region is displayed where the user can select the axis for which they wish to generate the model from the plurality of axes of the machine tool.

[0062] The user sets the predetermined requirement via an input form (not illustrated). For example, in FIG. 10, the result of setting the response frequency of the velocity control; target 50 Hz is displayed, with a dashed line on the graph indicating 50 Hz.

[0063] When the user presses the Apply button as illustrated in FIG. 10, the transfer characteristic calculation unit 214 changes the control characteristic or motor characteristic in accordance with the predetermined requirement and provides the generated transfer characteristic to the simulation execution unit 220. The details of the operation for changing the control characteristic or motor characteristic in accordance with the predetermined requirement will be described later.

[0064] The display screen displays the frequency response of the open-loop transfer function obtained from the results of reading the control information. When the user presses the Apply button and the transfer characteristic calculation unit 214 changes the control characteristic or motor characteristic in accordance with the predetermined requirement and generates the transfer function, the frequency response is updated based on the generated transfer characteristic. Instead of updating the frequency response, the frequency response based on the generated transfer characteristic may be overlaid. Furthermore, the frequency response may also be based on the transfer function G.sub.C(s), which takes feedback into consideration.

[0065] The timing at which the transfer characteristic calculation unit 214 changes the control characteristic or motor characteristic or alters the display of the frequency response does not have to be when the Apply button is pressed, and may be the timing at which the user changes the predetermined requirement via the input form.

[0066] The numerical values illustrated in FIG. 10 may be rounded to an appropriate number of decimal places; however, the number of decimal places does not hold special significance and is not particularly limited. In FIG. 10, although the unit of motor inertia is represented as inertia [kgm.sup.2] assuming a rotary motor, the unit of inertia may be mass [kg] in the case of using a linear motor. Similarly, although the unit of [Nm] is used for the spring constant assuming a torsional spring, the unit may be [N/m] in the case of using a tension-compression spring.

[0067] The natural frequency .sub.O and the resonance frequency .sub.P in Equation 3, the frequency in Equation 10, the frequency in Equation 15, and the frequency in Equation 16 are angular frequencies [rad/s]. On the other hand, although the unit of the frequencies illustrated in the screen display in FIG. 10 and the unit of the frequencies illustrated as the predetermined requirements are Hz, the unit may be rad/s.

[0068] The following describes an example of the operation where the transfer characteristic calculation unit 214 changes the control characteristic or motor characteristic to calculate the transfer characteristic of the machine tool to satisfy the predetermined requirement.

(1) Example of Changing the Position Control Proportional Gain K.SUB.P .to Satisfy the Predetermined Requirement Response Frequency of Position Control

[0069] In a case where the user specifies 10 Hz as the response frequency of position control, the angular frequency for 10 Hz is 102(rad/s)=62.832 (rad/s). As described using Equation 10, since the position control proportional gain K.sub.P equals the response frequency of position control, the position control proportional gain K.sub.P is set to K.sub.P=102=62.832. To satisfy the predetermined requirement response frequency of position control, the position control proportional gain K.sub.P, which is the control characteristic C.sub.P, is changed from 30 (as illustrated in FIG. 10) to 62.832.

(2) Example of Changing the Velocity Loop Gain K.SUB.V .to Satisfy the Predetermined Requirement Response Frequency of Velocity Control (in the Case of Velocity PI Control)

[0070] The left-hand side of Equation 15 is changed to fc[Hz], resulting in Equation 18 (hereinafter referred to as Equation 18).

[00017]$\begin{array}{cc}{f}_C=\frac{1}{2}\sqrt{\frac{2K_V{K}_1}{\left(1+R\right)J_M}}& \left[\mathrm{Equation}18\right]\end{array}$

[0071] In Equation 18, the motor inertia J.sub.M remains unchanged. Generally, in PI control, the velocity control integral gain K.sub.1 and the velocity control proportional gain K.sub.2 are adjusted in balance; however, here, the velocity control integral gain K.sub.1 remains unchanged. The user may choose whether to change either the velocity loop gain K.sub.V or the inertia ratio R via a detailed settings option (not illustrated in FIG. 10). Alternatively, the detailed settings may be configured to prioritize one over the other within the allowable range, and in a case where the predetermined requirement cannot be satisfied within this range, the other may be changed. The following describes an example where only the velocity loop gain K.sub.V is adjusted to satisfy the predetermined requirement of a 50 Hz response frequency for velocity control. By solving Equation 18 for the velocity loop gain K.sub.V, the velocity loop gain K.sub.V is given by Equation 19 (hereinafter referred to as Equation 19).

[00018]$\begin{array}{cc}{K}_V=\frac{2{}^2{f}_C^2\left(1+R\right)J_M}{K_1}=\frac{2{}^2{50}^{2}2.20.022}{1089}=2.193& \left[\mathrm{Equation}19\right]\end{array}$

[0072] In this manner, in order to satisfy the predetermined requirement response frequency of velocity control, the velocity loop gain K.sub.V, which is part of the control characteristic C.sub.V, is changed from 1 to 2.191, and the control characteristic C.sub.V is adjusted accordingly.

(3) Example of Changing the Velocity Loop Gain K.SUB.V .to Satisfy the Predetermined Requirement Response Frequency of Velocity Control (in the Case of Velocity P Control)

[0073] The left-hand side of Equation 17 is changed to fc [Hz], resulting in Equation 20 (hereinafter referred to as Equation 20).

[00019]$\begin{array}{cc}{f}_C=\frac{K_V}{2\left(1+R\right)J_M}& \left[\mathrm{Equation}20\right]\end{array}$

[0074] In Equation 20, the motor inertia J.sub.M remains unchanged. This example discusses the case of satisfying the predetermined requirement response frequency of velocity control (50 Hz) by changing the inertia ratio R. By solving Equation 20 for the inertia ratio R, the inertia ratio R is given by Equation 21 (hereinafter referred to as Equation 21).

[00020]$\begin{array}{cc}R=\frac{K_V}{2J_M{f}_C}-1=\frac{1}{20.02250}=-0.855& \left[\mathrm{Equation}21\right]\end{array}$

[0075] The inertia ratio R can take a value of 0 in the absence of load inertia but cannot take negative numbers. Therefore, assuming R=0 and solving Equation 20 for K.sub.V, the velocity loop gain K.sub.V is given by Equation 22 (hereinafter referred to as Equation 22).

[00021]$\begin{array}{cc}{K}_V=2f_C{J}_M=2500.022=6.912& \left[\mathrm{Equation}22\right]\end{array}$

[0076] In this manner, in order to satisfy the predetermined requirement response frequency of velocity control, the inertia ratio R of the motor characteristic M.sub.1 is changed from 1.2 to 0, and the velocity loop gain K.sub.V of the control characteristic C.sub.V is changed from 1 to 6.912, thereby adjusting both the motor characteristic M.sub.1 and the control characteristic C.sub.V.

(4) Example of Changing the Spring Constant K.sub.S to Satisfy the Predetermined Requirement Natural Frequency of the Machine Tool By solving Equation 3 for the spring constant K.sub.S, the natural frequency .sub.O is given by Equation 23 (hereinafter referred to as Equation 23).

[00022]$\begin{array}{cc}{K}_S={\mathrm{RJ}}_M{}_O^2& \left[\mathrm{Equation}23\right]\end{array}$

[0077] In a case where the predetermined requirement is a natural frequency of 45 Hz, the spring constant K.sub.S represented in Equation 24 (hereinafter referred to as Equation 24) is given by Equation 23.

[00023]$\begin{array}{cc}{K}_S={\mathrm{RJ}}_M{}_O^2=1.20.022{\left(245\right)}^2=2111\left[\mathrm{Nm}\right]& \left[\mathrm{Equation}24\right]\end{array}$

[0078] In this manner, in order to satisfy the predetermined requirement Natural Frequency of the Machine Tool, the spring constant K.sub.S of the motor characteristic is changed from 6878 [Nm] (as illustrated in FIG. 10) to 2111 [Nm], and the motor characteristic M.sub.2 is changed accordingly.

(5) Example of Changing the Spring Constant K.SUB.S .to Satisfy the Predetermined Requirement Resonance Frequency Between the Driving Unit and the Driven Unit of the Machine Tool

[0079] By solving Equation 3 for the spring constant K.sub.S, the resonance frequency op is given by Equation 25 (hereinafter referred to as Equation 25).

[00024]$\begin{array}{cc}{K}_S=J_M{}_P^2\frac{R}{1+R}& \left[\mathrm{Equation}25\right]\end{array}$

[0080] In a case where the predetermined requirement is a resonance frequency of 80 Hz, the spring constant K.sub.S represented in Equation 26 (hereinafter referred to as Equation 26) is given by Equation 25.

[00025]$\begin{array}{cc}{K}_S=J_M{}_P^2\frac{R}{1+R}=0.022{\left(280\right)}^2\frac{1.2}{2.2}=3032\left[\mathrm{Nm}\right]& \left[\mathrm{Equation}26\right]\end{array}$

[0081] In this manner, in order to satisfy the predetermined requirement Resonance Frequency of the Machine Tool, the spring constant K.sub.S of the motor characteristic is changed from 6878 [Nm] (as illustrated in FIG. 10) to 3032 [Nm], and the motor characteristic M.sub.2 is changed accordingly.

[0082] In the above examples (4) and (5), the spring constant K.sub.S was changed; however, the inertia ratio R may be changed instead. The user may choose whether to change either the spring constant K.sub.S or the inertia ratio R via a detailed settings option (not illustrated in FIG. 10). Alternatively, the detailed settings may be configured to prioritize one over the other within the allowable range, and in a case where the predetermined requirement cannot be satisfied within this range, the other may be changed.

[0083] The transfer characteristic calculation unit 214 can change at least one of the motor inertia, the inertia ratio, the spring constant, the position control proportional gain, the velocity loop gain, the velocity control integral gain, or the velocity control proportional gain such that the transfer characteristic satisfies the predetermined requirement. In a case where a rigid-body system model considering the load inertia as a motor characteristic is not used (i.e., R=0), the transfer characteristic calculation unit 214 can use a rigid-body system model that adds load inertia to motor inertia such that the transfer characteristic satisfies the predetermined requirement. In a case where a rigid-body system model is used as a motor characteristic, the transfer characteristic calculation unit 214 can use a two-inertia system model that connects the motor inertia and the load inertia with a spring element such that the transfer characteristic satisfies the predetermined requirement.

[0084] In the described embodiment, the control characteristic and the motor characteristic can be refined by considering more elements, leading to higher precision in the simulation results. For example, in academic literature, dampers (for damping) are considered when deriving the two-inertia system model first. However, even if the simulation execution unit considers dampers or nonlinearity, dampers or nonlinearity may be ignored in calculation for changing a control characteristic or obtaining a motor characteristic in the transfer characteristic calculation unit. This is because dampers and nonlinearity are less sensitive to frequency characteristics such as response frequency, natural frequency, and resonance frequency compared to control gains, inertia ratios, and spring constants.

[0085] Each configuration of the machining simulation device 200 has been described above. Next, the machining simulation method will be described. In the following description, the machining simulation method is described as being executed using the machining simulation device 200; however, the method may also be executed on devices other than the machining simulation device 200.

[0086] FIG. 11 is a flowchart illustrating the operation of the machining simulation device.

[0087] In Step S1, the control information acquisition unit 211 acquires control information that includes motor information and a control parameter. The motor information includes at least one of the motor inertia, the inertia ratio, or the spring constant, all caused by the motor. The control parameter includes, for example, at least one of the position control proportional gain K.sub.P of the position control unit, the velocity loop gain K.sub.V, the velocity control integral gain K.sub.1 of the velocity control unit, or the velocity control proportional gain K.sub.2 of the velocity control unit, which are included in the motor control system of the control unit 100.

[0088] In Step S2, the motor characteristic calculation unit 212 calculates the motor characteristic, based on the motor information.

[0089] In Step S3, the control characteristic calculation unit 213 calculates the control characteristic of the motor control system, based on the control parameter. Step S3 may be executed before Step S2, or may be executed in parallel with Step S2.

[0090] In Step S4, the transfer characteristic calculation unit 214 determines whether the predetermined requirement have been input. If the requirement has been input, the processing proceeds to Step S5; if the requirement has not been input, the processing proceeds to Step S6.

[0091] In Step S5, the transfer characteristic calculation unit 214 changes at least one of the motor characteristic calculated by the motor characteristic calculation unit 212 or the control characteristic calculated by the control characteristic calculation unit 213 so as to satisfy the predetermined requirement thus input.

[0092] In Step S6, if at least one of the motor characteristic or the control characteristic has been changed in Step S5, the transfer characteristic calculation unit 214 calculates the transfer characteristic of the machine tool, based on at least one of the motor characteristic or control characteristic thus changed. If the predetermined requirement has not been input in Step S4, the transfer characteristic calculation unit 214 calculates the transfer characteristic of the machine tool; based on at least one of the motor characteristic calculated by the motor characteristic calculation unit 212 or the control characteristic calculated by the control characteristic calculation unit 213.

[0093] In Step S7, the simulation execution unit 220 uses the calculated transfer function to simulate the control of the machine tool control device 100 based on the machining program, the behavior of the motor as the driving unit and the driven unit, and the feedback control of the control unit 100 based on the position information of the driving unit and the driven unit, and outputs the position information of each shaft as the simulation result.

[0094] The components of the machining simulation device described in the embodiment above can be implemented by hardware, software, or a combination thereof. Implementation by software means implementation by a computer that reads and executes a program. In order to implement the components of the machining simulation device by software or a combination thereof, the machining simulation device includes a processor such as a CPU (Central Processing Unit). The processor functions as the execution unit. The machining simulation device may include a plurality of processors operating in parallel. The machining simulation device also includes auxiliary storage devices such as HDDs (Hard Disk Drives) that store various programs such as application software or an OS (Operating System), and main storage devices such as RAM (Random Access Memory), which store the program and data temporarily required on the program during the execution of the functions and operations of the machining simulation device as described in FIGS. 1 and 6. The machining simulation device may include a plurality of main storage devices.

[0095] With the machining simulation device, the processor reads the application software or OS from the auxiliary storage device, loads the read application software or OS into the main storage device, and executes computation based on the application software or OS. Various hardware components of the machining simulation device are controlled based on the computation results. In this manner, the functional blocks of the present embodiment are implemented.

[0096] The components of the machining simulation device may also be implemented by hardware such as electronic circuits. When configuring the machining simulation device by hardware, some or all of the functions of the components of the machining simulation device may be implemented by ICs (Integrated Circuits) such as ASICs (Application Specific Integrated Circuits), gate arrays, FPGAs (Field Programmable Gate Arrays), or CPLDs (Complex Programmable Logic Devices).

[0097] Programs can be stored and supplied to a computer using various types of non-transitory computer-readable media. Non-transitory computer-readable media include various types of tangible storage media. Examples of non-transitory computer-readable media include magnetic storage media (such as hard disk drives), magneto-optical storage media (such as magneto-optical disks), CD-ROMs (Read Only Memory), CD-Rs, CD-RWs, and semiconductor memory (such as mask ROMS, PROMS (Programmable ROMs), EPROMS (Erasable PROMs), flash ROMs, and RAM (random access memory)). Programs may also be supplied to the computer using various types of transitory computer-readable media.

[0098] The effect of the machining simulation device and the machining simulation method described in the above embodiment is that the transfer function of the machine tool can be achieved without requiring the collection of data during test operation of the machine tool and without requiring specialized knowledge for calculating the transfer characteristic.

[0099] Although the present disclosure has been described above, the present disclosure is not limited to the individual embodiments and modifications described above. Various additions, replacements, modifications, partial deletions, etc., can be made without departing from the scope of the disclosure, as derived from the content of the claims and their equivalents. The embodiments and modifications may also be implemented in combination. For example, the order of operations and the sequence of processing described as examples in the above embodiments, and are not limited thereto.

[0100] The following additional notes are further disclosed regarding the embodiment above:

(Additional Note 1)

[0101] A machining simulation device (200) includes: a transfer characteristic generation unit (210) that generates a transfer characteristic of a machine tool; and a simulation execution unit (220) that simulates behavior of the machine tool using the transfer characteristic, in which the transfer characteristic generation unit includes: a control information acquisition unit (211) that acquires control information including motor information of the machine tool and a control parameter of the motor control system of the machine tool from a storage unit; a motor characteristic calculation unit (212) that calculates a motor characteristic, based on the motor information; a control characteristic calculation unit (213) that calculates a control characteristic of the motor control system, based on the control parameter; and a transfer characteristic calculation unit (214) that calculates the transfer characteristic satisfying a predetermined requirement, based on the motor characteristic and the control characteristic, in which the predetermined requirement includes at least one of the response frequency of position control, the response frequency of velocity control, the natural frequency of the machine tool, or the resonance frequency between the driving unit and the driven unit of the machine tool.

(Additional Note 2)

[0102] The machining simulation device as described in additional note 1, in which the motor information includes at least one of motor inertia, an inertia ratio, or a spring constant, all caused by the motor; and the control parameter includes at least one of position control proportional gain, velocity loop gain, velocity control integral gain, or velocity control proportional gain.

(Additional Note 3)

[0103] The machining simulation device as described in additional note 2, in which the motor characteristic calculation unit (212) calculates the motor characteristic using a rigid-body system model including the motor inertia and the inertia ratio, or a two-inertia system model including the motor inertia, the inertia ratio, and the spring constant.

(Additional Note 4)

[0104] The machining simulation device as described in additional note 2, in which the control characteristic calculation unit (213) calculates the control characteristic using the position control proportional gain.

(Additional Note 5)

[0105] The machining simulation device as described in additional note 2, in which the control characteristic calculation unit (213) calculates the control characteristic using the velocity loop gain, or the velocity loop gain, the velocity control integral gain, and the velocity control proportional gain.

(Additional Note 6)

[0106] The machining simulation device as described in additional note 1, in which the motor characteristic includes at least one of motor inertia, an inertia ratio, or a spring constant; the control characteristic includes at least one of position control proportional gain, velocity loop gain, velocity control integral gain, or velocity control proportional gain; and the transfer characteristic calculation unit (214) changes at least one of the motor inertia, the inertia ratio, the spring constant, the position control proportional gain, the velocity loop gain, the velocity control integral gain, or the velocity control proportional gain such that the transfer characteristic satisfies the predetermined requirement.

(Additional Note 7)

[0107] The machining simulation device as described in additional note 1, in which the transfer characteristic calculation unit (214) adds load inertia to motor inertia, or connects the motor inertia and the load inertia with a spring element, such that the transfer characteristic satisfies the predetermined requirement.

(Additional Note 8)

[0108] A machining simulation method causing a computer to execute processing including: processing of calculating a motor characteristic, based on motor information of a machine tool; processing of calculating a control characteristic of a motor control system, based on a control parameter of the motor control system of the machine tool; processing of calculating a transfer characteristic satisfying a predetermined requirement, based on the motor characteristic and the control characteristic; and processing of simulating behavior of the machine tool using the transfer characteristic, in which the predetermined requirement includes at least one of a response frequency of position control, a response frequency of velocity control, a natural frequency of the machine tool, or a resonance frequency between a driving unit and a driven unit of the machine tool.

EXPLANATION OF REFERENCE NUMERALS

[0109] 10: machining simulation system [0110] 100: machine tool control device [0111] 101: storage unit [0112] 200: machining simulation device [0113] 210: transfer characteristic generation unit [0114] 211: control information acquisition unit [0115] 212: motor characteristic calculation unit [0116] 213; control characteristic calculation unit [0117] 214: transfer characteristic calculation unit [0118] 220: simulation execution unit