WO2020208893A1 - 数値制御装置および学習装置 - Google Patents
数値制御装置および学習装置 Download PDFInfo
- Publication number
- WO2020208893A1 WO2020208893A1 PCT/JP2020/002108 JP2020002108W WO2020208893A1 WO 2020208893 A1 WO2020208893 A1 WO 2020208893A1 JP 2020002108 W JP2020002108 W JP 2020002108W WO 2020208893 A1 WO2020208893 A1 WO 2020208893A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- unit
- phase difference
- timing signal
- spindle
- angle
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H1/00—Measuring characteristics of vibrations in solids by using direct conduction to the detector
- G01H1/003—Measuring characteristics of vibrations in solids by using direct conduction to the detector of rotating machines
Definitions
- the present disclosure relates to a numerical control device for determining the occurrence of chatter vibration in a machine tool, a method for determining the occurrence of chatter vibration, and a learning device.
- Machine tools that removes the object to be machined into a desired shape by changing the relative position between the object to be machined and the tool.
- Machine tools such as milling machines and lathes mount tools or objects to be machined on the spindle and rotate the spindle to perform machining.
- vibration called “chatter vibration” may occur.
- chatter vibration occurs, the accuracy of the finished surface deteriorates and the tool may be damaged.
- vibration is measured twice at predetermined time intervals during processing, and the second measurement result has a larger amplitude in the frequency domain and a phase in the frequency domain as compared with the first measurement result.
- a method for determining that chatter vibration has occurred when the values are different is disclosed.
- the technique described in Patent Document 1 has a problem that it takes a long time to determine whether or not chatter vibration occurs during processing.
- frequency analysis such as FFT (Fast Fourier Transform) has a restriction that the frequency resolution ⁇ F is inversely proportional to the data acquisition time T due to its calculation algorithm. Therefore, in order to obtain a frequency resolution ⁇ F sufficient for determining chatter vibration, the data acquisition time becomes long. Further, in the technique described in Patent Document 1, since it is necessary to perform the measurement twice at a predetermined time, the time for measuring the vibration during processing becomes long.
- the present disclosure has been made in view of the above, and an object of the present disclosure is to obtain a numerical control device capable of shortening the time required for determining chatter vibration.
- the numerical control device of the present disclosure includes a drive command unit that gives an operation command to the spindle and the feed shaft of the machine tool, and a tool or an object to be machined attached to the machine machine.
- a sensor signal processing unit that generates multiple types of state quantities based on the sensor signal that detects vibration, and a timing signal generation unit that periodically generates and outputs a timing signal at a timing synchronized with the rotation of the spindle. From the phase difference calculation unit that generates a non-dimensional quantity indicating a plurality of types of state quantities in the state space for each timing signal and calculates the phase difference information indicating the difference between the phases of the non-dimensional quantity, and the phase difference information.
- It has an information observation unit that creates a data set for inference and an inference unit that has a trained inference model that performs machine learning to determine the presence or absence of chatter vibration by inputting the data set for inference. It is characterized by including a vibration determination unit that determines the occurrence of chatter vibration by inputting a data set for use.
- the figure for demonstrating the filter characteristic of the filter processing part shown in FIG. The figure which shows the relationship between the timing signal output by the timing signal generation part shown in FIG. 1 and the angle command of a spindle.
- the figure which shows the dimensionless quantity calculated by the phase difference calculation part shown in FIG. The figure which shows an example of the transition of the phase difference calculated by the phase difference calculation unit shown in FIG. 1 when chatter vibration occurs.
- the figure which shows the phase difference of the dimensionless quantity obtained from the experimental result shown in FIG. It is an experimental result when chatter vibration does not occur, and is a diagram showing the waveform of acceleration generated during machining.
- the figure which shows the phase difference of the dimensionless quantity obtained from the experimental result shown in FIG. The figure for demonstrating the 1st example of the process of the vibration determination part shown in FIG.
- FIG. 1 The figure which shows the structural example for realizing the 3rd example of the processing of the vibration determination part shown in FIG. Diagram showing a configuration example of an inference model using a neural network
- FIG. 2 The figure which shows the functional structure of the learning apparatus which learns the inference model used by the vibration determination part shown in FIG.
- FIG. 1 is a diagram showing a functional configuration of the numerical control device 1 according to the first embodiment.
- the numerical control device 1 numerically controls the machine tool 100 by giving an operation command 208 to the machine tool 100.
- the machine tool 100 has a spindle 101 and a feed shaft 102 equipped with a motor driven by an operation command 208.
- the object to be machined is installed on the spindle 101, and the tool is installed on the feed shaft 102.
- the machine tool 100 outputs operation information 209 including at least the positions, speeds, and motor currents of the spindle 101 and the feed shaft 102 to the numerical control device 1.
- the number of blades of the tool is one.
- a sensor 103 is attached to the machine tool 100.
- the sensor 103 outputs a sensor signal 200 that detects the vibration of the tool or the object to be machined to the numerical control device 1.
- the sensor 103 is attached to the structure of the spindle 101 or the feed shaft 102 of the machine tool 100.
- the position where the sensor 103 is attached may be any position where vibration of the tool or the object to be machined can be detected, and is preferably near the point where the tool and the object to be machined come into contact with each other.
- the type of the sensor 103 may be any sensor capable of detecting the vibration of the tool or the object to be machined, and the sensor 103 is, for example, a displacement sensor, a speed sensor, an acceleration sensor, an angular velocity sensor, or the like.
- the senor 103 may be a force sensor that detects a cutting reaction force or a microphone that detects a cutting sound during machining.
- the sensor 103 instead of the sensor 103, at least one of the positions, speeds, and motor currents of the spindle 101 and the feed shaft 102 included in the operation information is used to detect the vibration of the tool or the object to be machined during machining. Then, information indicating the detected vibration may be output as the sensor signal 200.
- the numerical control device 1 includes a filter processing unit 10, a sensor signal processing unit 11, a phase difference calculation unit 12, a vibration determination unit 13, a command value correction unit 14, a drive command unit 15, and a timing signal generation unit 16. And have.
- the drive command unit 15 receives the operation information 209 from the machine tool 100, outputs the operation command 208 to the machine tool 100 based on the operation information 209, and outputs the spindle operation command 206 to the timing signal generation unit 16. Further, the drive command unit 15 corrects the operation command 208 based on the correction signal 205 output by the command value correction unit 14 described later.
- the operation command 208 includes at least one of an angle command and a speed command of the spindle 101 of the machine tool 100, and at least one of a position command and a speed command of the feed shaft 102. Further, the spindle operation command 206 includes at least one of an angle command and a speed command of the spindle 101.
- the filter processing unit 10 filters the sensor signal 200 output by the machine tool 100 to generate a chatter vibration component signal 201, and outputs the generated chatter vibration component signal 201 to the sensor signal processing unit 11.
- the filter processing unit 10 removes the cutting component from the sensor signal 200.
- FIG. 2 is a diagram for explaining the filter characteristics of the filter processing unit 10 shown in FIG.
- the horizontal axis of FIG. 2 indicates the frequency, and the vertical axis indicates the amplitude.
- FIG. 2 shows a chatter vibration component, a cutting component, and a filter pass band in the frequency domain.
- the cutting component is mainly a vibration component having a frequency obtained by multiplying the spindle rotation speed by the number of blades of the tool, and a harmonic component thereof.
- the filter processing unit 10 removes a band containing each cutting component with a predetermined bandwidth.
- the predetermined bandwidth can be about several Hz near the cutting component.
- the filter processing unit 10 can be designed by using a known bandpass filter, notch filter, combination of a plurality of bandpass filters, a combination of a plurality of notch filters, a comb-shaped filter, and the like. By using the pass band filter shown in FIG. 2, the filter processing unit 10 can remove the cutting component from the sensor signal 200 to generate the chatter vibration component signal 201.
- the sensor signal processing unit 11 generates a plurality of types of state quantities based on the sensor signal 200, and outputs a state quantity signal 202 indicating the generated plurality of types of state quantities to the phase difference calculation unit 12. Specifically, the sensor signal processing unit 11 generates a state quantity signal 202 based on the chatter vibration component signal 201 output by the filter processing unit 10.
- the state quantity signal 202 includes a first state quantity and a second state quantity.
- the first state quantity is the chatter vibration component signal 201
- the second state quantity is the state quantity obtained by time-differentiating the chatter vibration component signal 201 once.
- the sensor signal processing unit 11 simultaneously converts the chatter vibration component signal 201 and the chatter vibration component signal 201 into a state quantity signal 202 as a time-series signal.
- the difference in the number of times of time differentiation between the first state quantity and the second state quantity may be any combination as long as it is an odd number.
- the first state quantity is the state quantity obtained by time-differentiating the chatter vibration component signal 201 P times
- the second state quantity is the state quantity obtained by time-differentiating Q times.
- P and Q are integers.
- the difference between P and Q may be an odd number.
- the dimensions of the first state quantity and the second state quantity may be the dimensions obtained by time-integrating the chatter vibration component signal 201.
- the first state quantity may be the acceleration and the second state quantity may be the velocity.
- the first state quantity may be a jerk which is a differential value of acceleration
- the second state quantity may be a position.
- the sensor signal processing unit 11 has, for example, two or more types of state quantities of the sensor signal 200, the state quantity obtained by time-differentiating the sensor signal 200, and the state quantity obtained by time-integrating the sensor signal. Can be generated.
- the timing signal generation unit 16 determines that the angle command of the spindle 101 has passed a predetermined angle based on the spindle operation command 206 output by the drive command unit 15, the timing signal generation unit 16 outputs the timing signal 207 to the phase difference calculation unit 12. To do. The details of the method of generating the timing signal 207 will be described.
- FIG. 3 is a diagram showing the relationship between the timing signal 207 output by the timing signal generation unit 16 shown in FIG. 1 and the angle command of the spindle 101.
- the angle command of the spindle 101 is included in the spindle operation command 206 output by the drive command unit 15.
- the angle command is a signal that takes a value of 0 to 360 degrees and returns to 0 degrees when it reaches 360 degrees.
- the timing signal generation unit 16 outputs a timing signal 207 each time the angle command of the spindle 101 passes a predetermined angle ⁇ 1. As a result, the timing signal 207 is periodically output at a timing synchronized with the rotation of the spindle 101 of the machine tool 100.
- the angle ⁇ 1 may be any angle as long as it is determined at one position during one rotation of the spindle 101, and can be, for example, the angle ⁇ 1 of the spindle 101 at the time of orientation of the spindle 101.
- the spindle orientation is a reference angle at which the spindle 101 is stopped.
- the timing signal generation unit 16 can also generate the timing signal 207 by using the speed command of the spindle 101 included in the spindle operation command 206.
- the timing signal generation unit 16 outputs the timing signal 207 for each time T1 calculated from the speed command S (rpm) of the spindle 101 using the following mathematical formula (1) with the initial time t0 as the reference time. To do.
- the timing signal generation unit 16 may use the angle command and the speed command of the spindle 101 together.
- the timing signal generation unit 16 uses the operation command of the spindle 101 to generate the timing signal 207, but instead of the operation command of the spindle 101 or in addition to the operation command of the spindle 101.
- the timing signal 207 may be generated by using the operation command of the feed shaft 102.
- the phase difference calculation unit 12 receives the timing signal 207 output by the timing signal generation unit 16 and the state quantity signal 202 output by the sensor signal processing unit 11. Based on the timing signal 207 and the state quantity signal 202, the phase difference calculation unit 12 generates a dimensionless quantity indicating a plurality of types of state quantities in the state space for each timing signal 207, and the dimensionless quantity between phases. The phase difference information 203 indicating the difference between the above is calculated. The phase difference calculation unit 12 outputs the generated phase difference information 203 to the vibration determination unit 13.
- the phase difference calculation unit 12 normalizes each of the plurality of state quantities included in the state quantity signal 202.
- the phase difference calculation unit 12 normalizes each state quantity by dividing each state quantity by a predetermined maximum value of each state quantity.
- the maximum value used here may be the maximum value of each state quantity obtained in the prior machining experiment, or may be the maximum value of each state quantity obtained in the prior simulation.
- phase difference calculation unit 12 receives the timing signal 207 at time t1 and time t2. At this time, t2 is a value larger than t1.
- the phase difference calculation unit 12 normalizes the state quantity signal 202 at time t1, and sets the normalized value as the first dimensionless quantity N1.
- the phase difference calculation unit 12 normalizes the state quantity signal 202 in the same manner as at time t1, and sets the normalized value as the second dimensionless quantity N2.
- FIG. 4 is a diagram showing a dimensionless quantity calculated by the phase difference calculation unit 12 shown in FIG. FIG. 4 shows that the first dimensionless quantity N1 and the second dimensionless quantity N2 calculated by the phase difference calculation unit 12 are the dimensions of the first state quantity of the state quantity signal 202 and the second state quantity. It is represented in the state space consisting of the dimensions of.
- the phase difference calculation unit 12 sets an angle ⁇ formed by the first dimensionless quantity N1 which is the previous dimensionless quantity, the origin of the state space, and the second dimensionless quantity N2 which is the current dimensionless quantity. It is calculated as the phase difference between the first dimensionless quantity N1 and the second dimensionless quantity N2.
- the phase difference calculation unit 12 performs the above calculation for each timing signal 207, and calculates the phase difference between the latest dimensionless quantity and the previous dimensionless quantity.
- the phase difference calculation unit 12 outputs the phase difference information 203 indicating the calculated phase difference to the vibration determination unit 13.
- FIG. 5 is a diagram showing an example of the transition of the phase difference calculated by the phase difference calculation unit 12 shown in FIG. 1 when chatter vibration is generated.
- the dimensionless quantities N10 to N13 are dimensionless quantities calculated at each time t10 to t13 when the timing signal 207 is received. At this time, t13> t12> t11> t10.
- the angle ⁇ 10 indicates the angle formed by the dimensionless quantity N10, the origin of the state space, and the dimensionless quantity N11
- the angle ⁇ 11 indicates the angle formed by the dimensionless quantity N11, the origin of the state space, and the dimensionless quantity N12.
- the angle ⁇ 12 indicates the angle formed by the dimensionless quantity N12, the origin of the state space, and the dimensionless quantity N13.
- FIG. 6 is a diagram showing an example of the transition of the phase difference calculated by the phase difference calculation unit 12 shown in FIG. 1 when chatter vibration does not occur.
- the dimensionless quantities N20 to N23 are dimensionless quantities calculated at each time t20 to t23 when the timing signal 207 is received. At this time, t23> t22> t21> t20.
- the angle ⁇ 20 indicates the angle formed by the dimensionless quantity N20, the origin of the state space, and the dimensionless quantity N21
- the angle ⁇ 21 indicates the angle formed by the dimensionless quantity N21, the origin of the state space, and the dimensionless quantity N22.
- the angle ⁇ 22 indicates the angle formed by the dimensionless quantity N22, the origin of the state space, and the dimensionless quantity N23.
- the phase difference calculation unit 12 calculates a constant value angle ⁇ , that is, a phase difference when chatter vibration occurs, and an irregular value angle ⁇ , that is, when chatter vibration does not occur. The phase difference will be calculated.
- Machining conditions when chatter vibration occurs are spindle speed 2000 (RPM), feed per blade 0.1 (mm), radial depth cut 1.0 (mm), and workpiece protrusion amount 200 (mm) from the chuck. Is.
- FIG. 7 is an experimental result when chatter vibration occurs, and is a diagram showing a waveform of acceleration generated during machining.
- FIG. 8 is a diagram showing a phase difference of a dimensionless quantity obtained from the experimental results shown in FIG. 7.
- the filter processing unit 10, the sensor signal processing unit 11, the timing signal generation unit 16, and the phase difference calculation unit 12 process the acceleration data output by the acceleration sensor as the sensor signal 200, so that the dimensionless unit shown in FIG. 8 is absent.
- the dimensionless quantity and the angle ⁇ were calculated.
- the sensor signal processing unit 11 calculates and outputs a state quantity signal 202 having the first state quantity as acceleration and the second state quantity as velocity.
- the velocity which is the second state quantity, was obtained by integrating the acceleration.
- the timing signal generation unit 16 outputs a timing signal 207 each time the angle command of the spindle 101 passes the spindle angle at the time of spindle orientation.
- the phase difference calculation unit 12 normalized the state quantity signal 202 and calculated the dimensionless quantities N101 to N106 at each time t101 to t106 when each timing signal 207 was received.
- FIG. 8 shows dimensionless quantities N101 to N106 in the state space.
- the angle ⁇ 101 indicates the angle formed by the dimensionless quantity N101, the origin of the state space, and the dimensionless quantity N102
- the angle ⁇ 102 indicates the angle formed by the dimensionless quantity N102, the origin of the state space, and the dimensionless quantity N103.
- the angle ⁇ 103 indicates the angle formed by the dimensionless quantity N103, the origin of the state space, and the dimensionless quantity N104.
- the angle ⁇ 104 indicates the angle formed by the dimensionless quantity N104, the origin of the state space, and the dimensionless quantity N105
- the angle ⁇ 105 indicates the angle formed by the dimensionless quantity N105, the origin of the state space, and the dimensionless quantity N106.
- the angles ⁇ 101 to ⁇ 105 are constant. From the above experimental results, it is experimentally shown that when chatter vibration occurs, the trajectory of the dimensionless quantity in the state space becomes circular, and the angles ⁇ 101 to ⁇ 105 at each time take a constant value. It was.
- the machining conditions when chatter vibration does not occur are: spindle speed 2000 (RPM), feed 0.1 (mm) per blade, radial depth cut 1.0 (mm), and workpiece protrusion amount 100 from the chuck. (Mm).
- FIG. 9 is an experimental result when chatter vibration does not occur, and is a diagram showing a waveform of acceleration generated during machining.
- FIG. 10 is a diagram showing a phase difference of a dimensionless quantity obtained from the experimental results shown in FIG. Also in this experiment, the filter processing unit 10, the sensor signal processing unit 11, the timing signal generation unit 16, and the phase difference calculation unit 12 process the acceleration data output by the acceleration sensor as the sensor signal 200, and are shown in FIG. The dimensionless quantity and the angle ⁇ were calculated.
- the sensor signal processing unit 11 calculates and outputs a state quantity signal 202 having the first state quantity as acceleration and the second state quantity as velocity.
- the velocity which is the second state quantity, was obtained by integrating the acceleration.
- the timing signal generation unit 16 outputs a timing signal 207 each time the angle command of the spindle 101 passes the spindle angle at the time of spindle orientation.
- the phase difference calculation unit 12 normalized the state quantity signal 202 and calculated the dimensionless quantities N201 to N206 at each time t201 to t206 when each timing signal 207 was received.
- FIG. 10 shows dimensionless quantities N201 to N206 in the state space.
- the angle ⁇ 201 indicates the angle formed by the dimensionless quantity N201, the origin of the state space, and the dimensionless quantity N202
- the angle ⁇ 202 indicates the angle formed by the dimensionless quantity N202, the origin of the state space, and the dimensionless quantity N203.
- the angle ⁇ 203 indicates the angle formed by the dimensionless quantity N203, the origin of the state space, and the dimensionless quantity N204.
- the angle ⁇ 204 indicates the angle formed by the dimensionless quantity N204, the origin of the state space, and the dimensionless quantity N205
- the angle ⁇ 205 indicates the angle formed by the dimensionless quantity N205, the origin of the state space, and the dimensionless quantity N206. From the above experimental results, it is experimentally shown that when chatter vibration does not occur, the trajectory of the dimensionless quantity in the state space is irregular, and the angles ⁇ 201 to ⁇ 205 at each time take indefinite values. It was.
- the vibration determination unit 13 receives an angle ⁇ , which is a phase difference, from the phase difference calculation unit 12, and determines chatter vibration using the procedure shown below.
- the vibration determination unit 13 determines whether or not chatter vibration has occurred based on the phase difference information 203 calculated by the phase difference calculation unit 12.
- the vibration determination unit 13 outputs the vibration determination information 204 to the command value correction unit 14. The details of the vibration determination unit 13 will be described.
- the vibration determination unit 13 records the phase difference information 203 output by the phase difference calculation unit 12. This recording is repeated every time the phase difference information 203 is received, and the angle ⁇ included in the phase difference information 203 is recorded as time series data.
- the vibration determination unit 13 calculates the difference of the recorded angle ⁇ , and if the difference of the angle ⁇ is less than a predetermined threshold value, determines that chatter vibration has occurred, and uses the vibration determination information 204 as a command value. Output to the correction unit 14.
- the threshold value used by the vibration determination unit 13 will be described.
- the phase difference calculation unit 12 outputs the angle ⁇ in the range of ⁇ 180 degrees.
- the angle ⁇ takes an indefinite value in the range of ⁇ 180 degrees, so the difference in angle ⁇ takes a value of ⁇ 360 degrees.
- chatter vibration occurs, the angle ⁇ takes a constant value, so that the difference between the angles ⁇ becomes zero.
- the actual angle ⁇ includes an error due to a measurement error or the like. Therefore, if the threshold value is set to a value sufficiently larger than the error of the angle ⁇ , the presence or absence of chatter vibration can be determined. For example, when the error of the angle ⁇ is ⁇ 10 degrees, the threshold value may be a value larger than 20 degrees, for example, 30 degrees.
- the vibration determination unit 13 calculates the standard deviation of the angle ⁇ included in the plurality of phase difference information 203 repeatedly received, and when the standard deviation of the angle ⁇ is less than a predetermined threshold value, chatter vibration occurs.
- FIG. 11 is a diagram for explaining a first example of the processing of the vibration determination unit 13 shown in FIG.
- the vibration determination unit 13 calculates the standard deviation using the angle ⁇ acquired between the time t3 and the time t4, and determines the presence or absence of chatter vibration. In the first example shown in FIG. 11, the standard deviation of the angle ⁇ exceeds the threshold of the standard deviation. Therefore, the vibration determination unit 13 can determine that chatter vibration has not occurred. Further, FIG.
- the vibration determination unit 13 can determine that chatter vibration has occurred.
- the threshold value of the standard deviation can be a value larger than 10 degrees, for example, 20 degrees.
- the vibration determination unit 13 may determine the presence or absence of chatter vibration by using the dispersion of the angle ⁇ instead of the standard deviation of the angle ⁇ included in the plurality of phase difference information 203 repeatedly received. It can.
- the variance threshold can be determined in the same manner as in the case of the standard deviation.
- the vibration determination unit 13 may determine chatter vibration by using the difference between the maximum value and the minimum value of the angle ⁇ . In this case, the vibration determination unit 13 calculates the difference between the maximum value and the minimum value of the angle ⁇ received during the period of the predetermined length, and the calculated difference between the maximum value and the minimum value is predetermined. If it is smaller than the threshold value, it can be determined that chatter vibration is occurring.
- the vibration determination unit 13 may determine the chatter vibration by using an inference model in which machine learning has been performed in advance so as to determine the presence or absence of chatter vibration.
- FIG. 22 is a diagram showing a configuration example for realizing a third example of the processing of the vibration determination unit 13 shown in FIG.
- the vibration determination unit 13 has an information observation unit 301 and an inference unit 302.
- the information observation unit 301 observes the angle ⁇ included in the phase difference information 203 for a predetermined number of samplings as time series data, and generates an inference data set 303.
- the inference unit 302 inputs the inference data set 303 generated by the information observation unit 301 into an inference model that has been machine-learned in advance to output the presence or absence of chatter vibration, thereby determining vibration indicating the presence or absence of chatter vibration.
- Information 204 is output.
- the inference unit 302 may use an inference model of any algorithm.
- an inference model using a neural network will be described.
- FIG. 23 is a diagram showing a configuration example of an inference model using a neural network.
- the neural network has an input layer x1, x2, ..., Xn having n neurons and an intermediate layer y1, y2, ..., Ym having m neurons, and 1 It is composed of an output layer z1 having an individual neuron.
- FIG. 23 shows an example in which the intermediate layer is one layer, two or more intermediate layers may be provided.
- the input layers x1, x2, ..., Xn are connected to the intermediate layers y1, y2, ..., Ym, and the intermediate layers y1, y2, ..., Ym are connected to the output layer z1.
- the connection between the input layer and the intermediate layer shown in FIG. 23 is an example, and each input layer x1, x2, ..., Xn is connected to any of the intermediate layers y1, y2, ..., Ym. May be good.
- FIG. 24 is a diagram showing a functional configuration of a learning device 400 that learns an inference model used by the vibration determination unit 13 shown in FIG. 22.
- the learning device 400 has a data acquisition unit 401 and a learning processing unit 402.
- the data acquisition unit 401 acquires a learning data set 304 that associates the phase difference information 203 acquired during actual processing with the chatter vibration presence / absence information 305 indicating the presence / absence of chatter vibration.
- the presence or absence of chatter vibration can be expressed using different numerical values depending on whether or not chatter vibration has occurred. For example, a numerical value such as "1" when chatter vibration occurs and "0" when it does not occur may be used.
- the chatter vibration presence / absence information 305 included in the learning data set 304 for example, the result determined by evaluating the machined surface after machining can be used.
- the learning processing unit 402 performs learning using the phase difference information 203 included in the learning data set 304 as input data.
- the learning processing unit 402 performs so-called supervised learning so that the output of the inference model matches the numerical value of the chatter vibration presence / absence information 3-5.
- the learning processing unit 402 inputs the phase difference information 203 to the input layer and adjusts the weights A1 to Aa and B1 to Bb so that the value output from the output layer approaches the numerical value indicating the presence or absence of chatter vibration. Then, when the phase difference information 203 is input, an inference model that outputs an output value corresponding to the presence or absence of chatter vibration is learned.
- the learning processing unit 402 can use the error back propagation method as a supervised learning method used when learning an inference model. In addition, in order to improve the generalization performance of the inference model, the learning processing unit 402 uses methods such as "dropout” that randomly excludes neurons during learning and “early stopping” that monitors errors and stops learning early. You may use it.
- the learning processing unit 402 outputs the inference model for which learning has been completed as a learned inference model.
- the learned inference model output by the learning processing unit 402 can be used in the inference unit 302 shown in FIG.
- the command value correction unit 14 When the command value correction unit 14 receives the vibration determination information 204 from the vibration determination unit 13, it outputs a correction signal 205 instructing the drive command unit 15 to change the operation command 208.
- the correction signal 205 is a signal that changes the speed command of the spindle 101 by a predetermined ratio corresponding to the current spindle speed.
- the correction signal 205 may be a signal for changing the rotation speed of the spindle 101 so that the rotation speed of the spindle 101 matches an integral fraction of the resonance frequency of the mechanical structure calculated in advance by finite element analysis. ..
- the resonance frequency of the mechanical structure a value identified by a preliminary experiment may be used instead of the finite element analysis.
- the correction signal 205 may be a signal that changes the speed command location of the feed shaft 102 by a predetermined ratio. Specifically, the correction signal 205 can be a signal that reduces the speed command value of the feed shaft 102 so that the motor current of the feed shaft 102 becomes smaller than a predetermined reference value.
- time series data of a plurality of types of state quantities are generated based on the sensor signal 200 that detects the vibration of the tool or the machining object attached to the machine tool 100.
- dimensionless quantities indicating a plurality of types of state quantities are generated for each timing signal 207, and phase difference information 203 indicating an angle ⁇ which is a difference between the phases of the dimensionless quantities is calculated.
- phase difference information 203 indicating an angle ⁇ which is a difference between the phases of the dimensionless quantities is calculated.
- the occurrence of chatter vibration in the machine tool 100 is determined.
- the numerical control device 1 does not use frequency analysis typified by FFT for determining chatter vibration, it is possible to obtain sufficient frequency resolution for determining chatter vibration without being restricted by data acquisition time. Therefore, the time required for determining the chatter vibration, that is, the time from the occurrence of the chatter vibration to the determination of the occurrence of the chatter vibration can be shortened. Therefore, the time from the occurrence of the chatter vibration to the suppression of the chatter vibration is shortened, and the chatter vibration can be suppressed at high speed.
- the numerical control device 1 determines the chatter vibration based on the phase difference information 203 which is a dimensionless quantity.
- a conventional method for example, there is a method of determining chatter vibration based on the amplitude of chatter vibration.
- the amplitude of chatter vibration is used, the magnitude of the amplitude differs depending on the machining conditions, so it is necessary to set the amplitude threshold value according to the machining conditions.
- the numerical control device 1 that determines the chatter vibration based on the phase difference information 203 which is a dimensionless quantity the chatter vibration can be determined without depending on the machining conditions, and the threshold value depends on the machining conditions. There is no need to provide.
- the machine tool 100 has a configuration in which an object to be machined is installed on the spindle 101, but it is a case where a tool is installed on the spindle 101, for example, as represented by a milling machine and a turning machine. Can produce the same effect.
- one sensor 103 is used, but the sensors 103 may be installed at a plurality of locations of the machine tool 100. In this case, it is possible to perform the process described in the first embodiment on all the installed sensors 103 and determine the occurrence of chatter vibration. When the sensors 103 are installed at a plurality of locations, even if chatter vibration occurs at the plurality of locations during processing, the chatter vibration can be determined and suppressed.
- the sensor 103 installed in the machine tool 100 may detect vibration using two or more types of sensors 103 and output a sensor signal 200 including two or more types of state quantities to the filter processing unit 10.
- the sensor signal processing unit 11 can set the two types of sensor outputs as the first state quantity and the second state quantity, respectively. By directly measuring the second state quantity with the sensor 103, it is possible to determine the chatter vibration without being affected by the quantization error associated with the differentiation or integration.
- the sensor signal processing unit 11 extracts two types of state quantities, but three or more types of state quantities may be extracted.
- Embodiment 2 In the first embodiment, an example in which the number of blades of the tool is one has been described. In the second embodiment, a configuration for determining the occurrence of chatter vibration when the number of blades of the tool is two or more will be described.
- FIG. 13 is a diagram showing a functional configuration of the numerical control device 2 according to the second embodiment.
- the numerical control device 2 includes a filter processing unit 10, a sensor signal processing unit 11, a phase difference calculation unit 12, a vibration determination unit 13, a command value correction unit 14, a drive command unit 15, and a timing signal generation unit 16. It has -1 and a tool information recording unit 17.
- the numerical control device 2 has a tool information recording unit 17 for outputting tool information 210 in addition to the configuration of the numerical control device 1 according to the first embodiment, and replaces the timing signal generation unit 16 of the numerical control device 1. , Has a timing signal generation unit 16-1.
- a tool information recording unit 17 for outputting tool information 210 in addition to the configuration of the numerical control device 1 according to the first embodiment, and replaces the timing signal generation unit 16 of the numerical control device 1.
- the tool information recording unit 17 records the tool information 210, which is information about the tool installed in the machine tool 100, and outputs the recorded tool information 210 to the timing signal generation unit 16-1.
- the tool information 210 includes at least information indicating the number of blades of the tool. Further, the tool information 210 may include information indicating a tool shape such as a tool type such as an end mill and a cutting tool, a tool length, and a tool diameter.
- the tool information recording unit 17 outputs tool information 210 indicating that the number of blades of the tool is one. Further, even when a turning tool is attached to the machine tool 100, when a plurality of tools cut at the same time, the number of blades is the same as the number of tools. For example, when the machine tool 100 has a lower tool post and an upper tool post, and machining is performed with a turning tool attached to both tool stands one by one, the tool information recording unit 17 has two blades. The tool information 210 indicating that is output.
- the timing signal generation unit 16-1 determines that the spindle angle has passed a predetermined angle based on the spindle operation command 206 output by the drive command unit 15 and the tool information 210 output by the tool information recording unit 17. , The timing signal 207 is output to the phase difference calculation unit 12 at the determined timing.
- ⁇ is a natural number of 2 or more.
- FIG. 14 is a diagram showing a first example of the relationship between the timing signal 207 output by the timing signal generation unit 16-1 shown in FIG. 13 and the angle command of the spindle 101.
- the angle command is a signal that takes a value between 0 degrees and 360 degrees and returns to 0 degrees when it reaches 360 degrees.
- the timing signal 207 is output each time the angle command of the spindle 101 passes through a plurality of predetermined angles ⁇ 1 and angle ⁇ 2. That is, the timing signal 207 is output twice, which is the same number as the number of blades, while the spindle 101 makes one rotation.
- the angle ⁇ 1 may be any angle as long as it is set to 1 during one rotation of the spindle 101.
- the angle ⁇ 1 can be the spindle angle at the time of spindle orientation.
- FIG. 15 is a diagram showing a second example of the relationship between the timing signal 207 output by the timing signal generation unit 16-1 shown in FIG. 13 and the angle command of the spindle 101.
- the timing signal 207 is output each time the angle command of the spindle 101 passes through a predetermined angle ⁇ ⁇ .
- ⁇ is a natural number of 2 or more and ⁇ or less.
- the timing signal 207 is output ⁇ times while the spindle 101 makes one rotation.
- ⁇ 1 may be any angle as long as it is set to 1 during one rotation of the spindle 101.
- the timing signal generation unit 16-1 outputs a timing signal 207 each time the spindle 101 rotates by an angle obtained by dividing 360 degrees by the number of blades ⁇ .
- the timing signal generation unit 16-1 may generate the timing signal 207 by using the speed command of the spindle 101 included in the spindle operation command 206. In this case, the timing signal generation unit 16-1 generates the timing signal 207 for each time T2 calculated from the speed command S (rpm) of the spindle 101 using the following mathematical formula (3).
- the timing signal 207 is generated based on the tool information 210, even when the number of blades of the tool attached to the machine tool 100 is a plurality, the implementation is performed. Similar to the first embodiment, the time required for determining the chatter vibration, that is, the time from the occurrence of the chatter vibration to the determination of the occurrence of the chatter vibration can be shortened.
- the numerical control device 2 can determine and suppress chatter vibration by the same processing as the numerical control device 1.
- FIG. 16 is a diagram showing a functional configuration of the numerical control device 3 according to the third embodiment.
- the numerical control device 3 includes a filter processing unit 10, a sensor signal processing unit 11, a phase difference calculation unit 12, a vibration determination unit 13, a command value correction unit 14, a drive command unit 15, and a timing signal generation unit 16. It has -2.
- the numerical control device 3 has a timing signal generation unit 16-2 in place of the timing signal generation unit 16 of the numerical control device 1.
- the parts different from the numerical control device 1 will be mainly described.
- the machine tool 100 has an angle information sensor 104 attached to the spindle 101 in addition to the configuration described in the first embodiment.
- the angle information sensor 104 is a sensor that measures the spindle angle, and for example, an encoder, a potentiometer sensor, or the like can be used.
- the angle information 211 is a signal including at least the measured value of the spindle angle measured by the angle information sensor 104.
- the angle information sensor 104 outputs the generated angle information 211 to the numerical control device 3.
- the timing signal generation unit 16-2 generates the timing signal 207 based on the measured value of the spindle angle included in the angle information 211.
- the method in which the timing signal generation unit 16-2 generates the timing signal 207 is the same as that in the first embodiment, but the measured value of the spindle angle is used instead of the angle command of the spindle 101.
- the angle information sensor 104 may be a sensor that detects the measured value of the spindle speed.
- the angle information sensor 104 is an angular velocity sensor.
- the timing signal generation unit 16-2 can generate the timing signal 207 based on the measured value of the spindle speed instead of the angle command of the spindle 101.
- the timing signal generation unit 16-2 uses the initial time t0 as a reference time to generate a timing signal 207 for each time T3 calculated from the measured value S ref (rpm) of the spindle speed using the following mathematical formula (4). Output.
- the timing signal 207 is generated based on the measured value of the spindle angle. Therefore, the timing signal 207 can be generated without being affected by the control tracking error that occurs between the angle command of the spindle 101 and the actual angle. Therefore, the numerical control device 3 is the time required for determining the chatter vibration, that is, even when a control tracking error occurs between the angle command of the spindle 101 and the actual angle, as in the first embodiment. It is possible to shorten the time from the occurrence of chatter vibration to the determination of the occurrence of chatter vibration.
- timing signal generation unit 16-2 describes the spindle angle command and the measured value of the spindle angle. May be used together, or the spindle speed command and the measured value of the spindle speed may be used together.
- FIG. 17 is a diagram showing a functional configuration of the numerical control device 4 according to the fourth embodiment.
- the numerical control device 4 includes a filter processing unit 10, a sensor signal processing unit 11, a phase difference calculation unit 12, a vibration determination unit 13, a command value correction unit 14, a drive command unit 15, and a timing signal generation unit 16. It has a -3 and a tool information recording unit 17. That is, the numerical control device 4 has a tool information recording unit 17 in addition to the configuration of the numerical control device 3 described in the third embodiment, and instead of the timing signal generation unit 16-2, the timing signal generation unit 16- Has 3.
- the parts different from the numerical control device 2 will be mainly described.
- the timing signal generation unit 16-3 generates a timing signal 207 based on the angle information 211 output by the angle information sensor 104 and the tool information 210 output by the tool information recording unit 17, and generates the generated timing signal 207. Output to the phase difference calculation unit 12.
- the timing signal generation unit 16-3 can generate the timing signal 207 in the same manner as the timing signal generation unit 16-1 of the second embodiment. At this time, the timing signal generation unit 16-3 generates the timing signal 207 by using the measured value of the spindle angle instead of the angle command of the spindle 101.
- the timing signal 207 is transmitted without being affected by the control tracking error that occurs between the angle command of the spindle 101 and the actual angle, as in the third embodiment.
- the time required for determining the chatter vibration that is, the chatter vibration is generated as in the first embodiment. It is possible to shorten the time from the occurrence to the determination of the occurrence of chatter vibration.
- FIG. 18 is a diagram showing a functional configuration of the numerical control device 5 according to the fifth embodiment.
- the numerical control device 5 includes a filter processing unit 10, a sensor signal processing unit 11, a phase difference calculation unit 12, a vibration determination unit 13, a command value correction unit 14, a drive command unit 15, and a timing signal generation unit 16. It has -4 and. Further, the numerical control device 5 can acquire the machining program information 212 including the spindle rotation speed command 213 from the machining program 18 for controlling the machine tool 100.
- the parts different from the numerical control device 1 will be mainly described.
- the spindle speed command 213 is output to the timing signal generation unit 16-4, and the machining program information 212 is output to the drive command unit 15.
- the drive command unit 15 can output an operation command 208 to the machine tool 100 according to the machining program information 212.
- the timing signal generation unit 16-4 generates the timing signal 207 based on the spindle speed command 213 extracted from the machining program information 212.
- the timing signal generation unit 16-4 outputs the timing signal 207 for each time T4 calculated by using the spindle rotation speed command 213 and the mathematical formula (5) shown below with the initial time t0 as the reference time.
- the S cmd in the mathematical formula (5) indicates the spindle rotation speed indicated by the spindle rotation speed command 213.
- the timing signal 207 can be generated based on the spindle speed command 213 included in the machining program 18. Even when the timing signal 207 is generated by using the machining program 18, the time required for determining the chatter vibration, that is, the occurrence of the chatter vibration after the chatter vibration is generated is determined as in the first embodiment. It is possible to shorten the time required to do so.
- FIG. 19 is a diagram showing a functional configuration of the numerical control device 6 according to the sixth embodiment.
- the numerical control device 6 includes a filter processing unit 10, a sensor signal processing unit 11, a phase difference calculation unit 12, a vibration determination unit 13, a command value correction unit 14, a drive command unit 15, and a timing signal generation unit 16. It has -5 and a tool information recording unit 17.
- the numerical control device 6 has a tool information recording unit 17 in addition to the configuration of the numerical control device 5, and has a timing signal generation unit 16-5 in place of the timing signal generation unit 16-4.
- the timing signal generation unit 16-5 generates a timing signal 207 based on the tool information 210 and the spindle speed command 213, and outputs the generated timing signal 207 to the phase difference calculation unit 12.
- the timing signal generation unit 16-5 uses the initial time t0 as a reference time, and uses the following mathematical formula (6), the number of blades ⁇ , and the spindle rotation speed S cmd to calculate the timing signal 207 for each time T5. Is output.
- Each functional unit of the numerical control devices 1 to 6 is realized by a processing circuit.
- These processing circuits may be realized by dedicated hardware, or may be control circuits using a CPU (Central Processing Unit).
- FIG. 20 is a diagram showing dedicated hardware for realizing the functions of the numerical control devices 1 to 6 according to the first to sixth embodiments.
- the processing circuit 90 is a single circuit, a composite circuit, a programmed processor, a parallel programmed processor, an ASIC (Application Specific Integrated Circuit), an FPGA (Field Programmable Gate Array), or a combination thereof.
- ASIC Application Specific Integrated Circuit
- FPGA Field Programmable Gate Array
- FIG. 21 is a diagram showing a configuration of a control circuit 91 for realizing the functions of the numerical control devices 1 to 6 according to the first to sixth embodiments.
- the control circuit 91 includes a processor 92 and a memory 93.
- the processor 92 is a CPU, and is also called a central processing unit, a processing unit, an arithmetic unit, a microprocessor, a microcomputer, a DSP (Digital Signal Processor), or the like.
- the memory 93 is, for example, a non-volatile or volatile semiconductor memory such as RAM (Random Access Memory), ROM (Read Only Memory), flash memory, EPROM (Erasable Programmable ROM), EEPROM (registered trademark) (Electrically EPROM). Magnetic disks, flexible disks, optical disks, compact disks, mini disks, DVDs (Digital Versatile Disk), etc.
- RAM Random Access Memory
- ROM Read Only Memory
- flash memory EPROM (Erasable Programmable ROM), EEPROM (registered trademark) (Electrically EPROM).
- Magnetic disks flexible disks, optical disks, compact disks, mini disks, DVDs (Digital Versatile Disk), etc.
- the control circuit 91 When the above processing circuit is realized by the control circuit 91, it is realized by the processor 92 reading and executing the program corresponding to the processing of each component stored in the memory 93.
- the memory 93 is also used as a temporary memory in each process executed by the processor 92.
- the configuration shown in the above-described embodiment shows an example, and can be combined with another known technique, can be combined with each other, and does not deviate from the gist. It is also possible to omit or change a part of the configuration.
- 1,2,3,4,5,6 Numerical control device 10 Filter processing unit, 11 Sensor signal processing unit, 12 Phase difference calculation unit, 13 Vibration determination unit, 14 Command value correction unit, 15 Drive command unit, 16, 16-1, 16-2, 16-3, 16-4, 16-5 Timing signal generator, 17 Tool information recording unit, 18 Machining program, 90 Processing circuit, 91 Control circuit, 92 Processor, 93 Memory, 100 Machining Machine, 101 spindle, 102 feed axis, 103 sensor, 104 angle information sensor, 200 sensor signal, 201 chatter vibration component signal, 202 state quantity signal, 203 phase difference information, 204 vibration judgment information, 205 correction signal, 206 spindle operation command , 207 timing signal, 208 operation command, 209 operation information, 210 tool information, 211 angle information, 212 machining program information, 213 spindle rotation command, 301 information observation unit, 302 inference unit, 303 inference data set, 304 for learning Data set, 305 chatter vibration presence / absence information, 400 learning device, 401 data acquisition
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automatic Control Of Machine Tools (AREA)
- Numerical Control (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
数値制御装置(1)は、工作機械(100)の主軸(101)および送り軸(102)に運転指令を与える駆動指令部(15)と、工作機械(100)に取り付けられる工具または加工対象物の振動を検出したセンサ信号(200)に基づいて、複数の種類の状態量を生成するセンサ信号処理部(11)と、主軸(101)の回転と同期したタイミングで周期的にタイミング信号(207)を生成して出力するタイミング信号生成部(16)と、状態空間内において複数の種類の状態量を示す無次元量を、タイミング信号(207)毎に生成し、無次元量の位相間の差分を示す位相差情報(203)を演算する位相差演算部(12)と、位相差情報(203)から推論用データセットを作成する情報観測部と、推論用データセットを入力することでびびり振動の有無を判定する機械学習を行った学習済み推論モデルを持つ推論部とを有し、推論部に推論用データセットを入力してびびり振動の発生を判定する振動判定部(13)と、を備えることを特徴とする。
Description
本開示は、工作機械におけるびびり振動の発生を判定する数値制御装置、びびり振動の発生判定方法、および学習装置に関する。
加工対象物と工具との間の相対位置を変化させることで、加工対象物を所望の形状に除去加工する工作機械がある。フライス盤および旋盤に代表される工作機械は、工具または加工対象物を主軸に装着し、主軸を回転させて加工を行う。加工中には、「びびり振動」と呼ばれる振動が発生することがある。びびり振動が発生すると、仕上げ面精度が悪化したり、工具の破損の原因となったりする。
特許文献1には、加工中に予め定められた時間間隔をおいて振動を2回測定し、2回目の測定結果が1回目と比較して、周波数領域の振幅が大きく、かつ周波数領域の位相が異なっていた場合に、びびり振動が発生したと判定する方法が開示されている。
しかしながら、上記特許文献1に記載の技術では、加工中にびびり振動の発生有無を判定するまでの時間が長くなってしまうという問題があった。具体的には、特許文献1に記載の技術では、振動の解析に周波数解析を用いる。FFT(Fast Fourier Transform)などの周波数解析は、その計算アルゴリズム上、周波数分解能ΔFがデータ取得時間Tに反比例するという制約がある。このため、びびり振動の判定に十分な周波数分解能ΔFを得るためには、データ取得時間が長くなる。また、特許文献1に記載の技術では、予め定められた時間をおいて2回の測定を行う必要があるため、加工中の振動を計測する時間も長くなってしまう。
本開示は、上記に鑑みてなされたものであって、びびり振動の判定に要する時間を短縮することが可能な数値制御装置を得ることを目的とする。
上述した課題を解決し、目的を達成するために、本開示の数値制御装置は、工作機械の主軸および送り軸に運転指令を与える駆動指令部と、工作機械に取り付けられる工具または加工対象物の振動を検出したセンサ信号に基づいて、複数の種類の状態量を生成するセンサ信号処理部と、主軸の回転と同期したタイミングで周期的にタイミング信号を生成して出力するタイミング信号生成部と、状態空間内において複数の種類の状態量を示す無次元量を、タイミング信号毎に生成し、無次元量の位相間の差分を示す位相差情報を演算する位相差演算部と、位相差情報から推論用データセットを作成する情報観測部と、推論用データセットを入力することでびびり振動の有無を判定する機械学習を行った学習済み推論モデルを持つ推論部とを有し、推論部に推論用データセットを入力することでびびり振動の発生を判定する振動判定部と、を備えることを特徴とする。
本開示によれば、びびり振動の判定に要する時間を短縮することができるという効果を奏する。
以下に、本開示の実施の形態にかかる数値制御装置、びびり振動の発生判定方法、および学習装置を図面に基づいて詳細に説明する。なお、この実施の形態によりこの開示の技術的範囲が限定されるものではない。
実施の形態1.
図1は、実施の形態1にかかる数値制御装置1の機能構成を示す図である。数値制御装置1は、工作機械100に運転指令208を与えることによって、工作機械100を数値制御する。
図1は、実施の形態1にかかる数値制御装置1の機能構成を示す図である。数値制御装置1は、工作機械100に運転指令208を与えることによって、工作機械100を数値制御する。
工作機械100は、運転指令208により駆動するモータを備えた主軸101および送り軸102を有する。本実施の形態では、主軸101には加工対象物が設置され、送り軸102には工具が設置される。また、工作機械100は、主軸101と送り軸102の位置、速度、およびモータ電流を少なくとも含む動作情報209を数値制御装置1に出力する。なお、本実施の形態では、工具の刃数は1枚とする。
工作機械100には、センサ103が取り付けられている。センサ103は、工具または加工対象物の振動を検出したセンサ信号200を数値制御装置1に出力する。センサ103は、工作機械100の主軸101または送り軸102の構造物に取り付けられる。センサ103を取り付ける位置は、工具または加工対象物の振動を検出可能な位置であればよく、工具と加工対象物とが接触する点の近傍であることが好ましい。センサ103の種類は、工具または加工対象物の振動を検出することができるセンサであればよく、センサ103は、例えば、変位センサ、速度センサ、加速度センサ、角速度センサなどである。また、センサ103は、切削反力を検出する力センサ、加工中の切削音を検出するマイクロフォンであってもよい。或いは、センサ103の代わりに動作情報に含まれる主軸101および送り軸102の位置と速度とモータ電流とのうち、少なくとも1つ以上を用いて、加工中に生じる工具または加工対象物の振動を検出し、検出した振動を示す情報をセンサ信号200として出力してもよい。
数値制御装置1は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16とを有する。
駆動指令部15は、工作機械100から動作情報209を受け取り、動作情報209に基づいて工作機械100に運転指令208を出力し、タイミング信号生成部16へ主軸運転指令206を出力する。また、駆動指令部15は、後述の指令値補正部14が出力する補正信号205に基づいて運転指令208を補正する。運転指令208は、工作機械100の主軸101の角度指令、および速度指令の少なくともいずれかと、送り軸102の位置指令および速度指令の少なくともいずれかとを含む。また、主軸運転指令206は、主軸101の角度指令および速度指令の少なくともいずれかを含む。
フィルタ処理部10は、工作機械100が出力するセンサ信号200をフィルタリングしてびびり振動成分信号201を生成し、生成したびびり振動成分信号201をセンサ信号処理部11に出力する。フィルタ処理部10は、センサ信号200から切削成分を除去する。
図2は、図1に示すフィルタ処理部10のフィルタ特性を説明するための図である。図2の横軸は周波数を示し、縦軸は振幅を示す。図2は、周波数領域において、びびり振動成分と、切削成分と、フィルタの通過帯域とを示している。切削成分は、主に、主軸回転数に工具の刃数を乗じた周波数の振動成分、および、その高調波成分である。フィルタ処理部10は、各切削成分を含む帯域を、予め定められた帯域幅で除去する。例えば、予め定められた帯域幅は、切削成分付近の数Hz程度とすることができる。フィルタ処理部10は、公知のバンドパスフィルタ、ノッチフィルタ、複数のバンドパスフィルタの組合せ、複数のノッチフィルタの組合せ、櫛形フィルタなどを用いて設計することができる。図2に示す通過帯域のフィルタを用いることで、フィルタ処理部10は、センサ信号200から切削成分を除去してびびり振動成分信号201を生成することができる。
図1の説明に戻る。センサ信号処理部11は、センサ信号200に基づいて、複数の種類の状態量を生成し、生成した複数の種類の状態量を示す状態量信号202を位相差演算部12に出力する。具体的には、センサ信号処理部11は、フィルタ処理部10が出力するびびり振動成分信号201に基づいて、状態量信号202を生成する。状態量信号202は、第1の状態量および第2の状態量を含む。例えば、第1の状態量は、びびり振動成分信号201であり、第2の状態量は、びびり振動成分信号201を1回時間微分することにより得られる状態量である。センサ信号処理部11は、びびり振動成分信号201およびびびり振動成分信号201を1回時間微分した信号を同時に時系列信号として状態量信号202とする。
なお、第1の状態量および第2の状態量の時間微分の回数の差は、奇数であればどのような組合せであってもよい。例えば、第1の状態量がびびり振動成分信号201をP回時間微分することにより得られる状態量であり、第2の状態量がQ回時間微分することにより得られる状態量であるとする。ここでPおよびQは整数とする。この場合、PとQの差分が奇数であればよい。さらに、第1の状態量と第2の状態量の次元は、びびり振動成分信号201を時間積分して得られる次元であってもよい。具体例として、第1の状態量を加速度とし、第2の状態量を速度としてもよい。或いは、第1の状態量を加速度の微分値であるジャークとし、第2の状態量を位置としてもよい。センサ信号処理部11は、例えば、センサ信号200と、センサ信号200を時間微分することにより得られる状態量と、センサ信号を時間積分することにより得られる状態量とのうち2種類以上の状態量を生成することができる。
タイミング信号生成部16は、駆動指令部15が出力する主軸運転指令206に基づき、主軸101の角度指令が予め定められた角度を通過したと判定すると、位相差演算部12へタイミング信号207を出力する。タイミング信号207の生成方法の詳細を説明する。
図3は、図1に示すタイミング信号生成部16が出力するタイミング信号207と主軸101の角度指令との関係を示す図である。主軸101の角度指令は、駆動指令部15が出力する主軸運転指令206に含まれる。図3において角度指令は、0度から360度の値をとり、360度に達すると0度に戻る信号である。タイミング信号生成部16は、主軸101の角度指令が所定の角度ψ1を通過する度にタイミング信号207を出力する。これにより、タイミング信号207は、工作機械100の主軸101の回転と同期したタイミングで周期的に出力されることになる。角度ψ1は主軸101が1回転する間に1箇所定められていれば、任意の角度であってよく、例えば主軸101のオリエンテーション時の主軸101の角度ψ1とすることができる。主軸オリエンテーションとは、主軸101を停止させる基準角度である。
或いは、タイミング信号生成部16は、主軸運転指令206に含まれる主軸101の速度指令を用いて、タイミング信号207を生成することもできる。この場合、タイミング信号生成部16は、初期時刻t0を基準時刻として、主軸101の速度指令S(rpm)から以下に示す数式(1)を用いて算出される時間T1ごとにタイミング信号207を出力する。また、タイミング信号生成部16は、主軸101の角度指令および速度指令を合わせて使用してもよい。なお、ここではタイミング信号生成部16は、主軸101の運転指令を用いて、タイミング信号207を生成することとしたが、主軸101の運転指令に代えて、あるいは、主軸101の運転指令に加えて、送り軸102の運転指令を用いて、タイミング信号207を生成してもよい。
位相差演算部12は、タイミング信号生成部16が出力するタイミング信号207と、センサ信号処理部11が出力する状態量信号202とを受け付ける。位相差演算部12は、タイミング信号207および状態量信号202に基づいて、状態空間内において複数の種類の状態量を示す無次元量を、タイミング信号207毎に生成し、無次元量の位相間の差分を示す位相差情報203を演算する。位相差演算部12は、生成した位相差情報203を振動判定部13に出力する。
位相差演算部12は、状態量信号202に含まれる複数の状態量のそれぞれを正規化する。位相差演算部12は、予め定められたそれぞれの状態量の最大値で、各状態量を除することでそれぞれの状態量を正規化する。ここで用いる最大値は、事前の加工実験で得られた各状態量の最大値としてもよいし、事前のシミュレーションで求められた各状態量の最大値としてもよい。
位相差演算部12が時刻t1および時刻t2でタイミング信号207を受け取ったとする。このときt2はt1よりも大きい値とする。位相差演算部12は、時刻t1における状態量信号202を正規化し、正規化後の値を第1の無次元量N1とする。位相差演算部12は、時刻t2において、時刻t1と同様に状態量信号202を正規化し、正規化後の値を第2の無次元量N2とする。
図4は、図1に示す位相差演算部12が算出する無次元量を示す図である。図4は、位相差演算部12が算出する第1の無次元量N1と、第2の無次元量N2とを、状態量信号202の第1の状態量の次元と、第2の状態量の次元とからなる状態空間内で表している。位相差演算部12は、前回の無次元量である第1の無次元量N1と、状態空間の原点と、今回の無次元量である第2の無次元量N2とがなす角度θを、第1の無次元量N1と第2の無次元量N2との位相差として演算する。位相差演算部12は、タイミング信号207毎に上記の演算を行い、最新の無次元量と、前回の無次元量との位相差を演算する。位相差演算部12は、演算した位相差を示す位相差情報203を振動判定部13に出力する。
図5は、びびり振動が発生している場合に図1に示す位相差演算部12が演算する位相差の推移の一例を示す図である。無次元量N10~N13は、タイミング信号207を受け取った各時刻t10~t13で算出された無次元量である。このとき、t13>t12>t11>t10とする。また、角度θ10は、無次元量N10と状態空間の原点と無次元量N11とがなす角度を示し、角度θ11は、無次元量N11と状態空間の原点と無次元量N12とがなす角度を示し、角度θ12は、無次元量N12と状態空間の原点と無次元量N13とがなす角度を示す。
びびり振動が発生している場合、状態空間内の無次元量N10~N13の軌跡は円状となり、各無次元量の間の角度は一定となる。したがって、タイミング信号207を受け取った各時刻における角度θは一定となる。
図6は、びびり振動が発生していない場合に図1に示す位相差演算部12が演算する位相差の推移の一例を示す図である。無次元量N20~N23は、タイミング信号207を受け取った各時刻t20~t23で算出された無次元量である。このとき、t23>t22>t21>t20とする。また、角度θ20は、無次元量N20と状態空間の原点と無次元量N21とがなす角度を示し、角度θ21は、無次元量N21と状態空間の原点と無次元量N22とがなす角度を示し、角度θ22は、無次元量N22と状態空間の原点と無次元量N23とがなす角度を示す。
びびり振動が発生していない場合、状態空間内の無次元量N20~N23の軌跡と角無次元量の間の角度は不規則になる。したがって、タイミング信号207を受け取った各時刻における角度θは一定とならない。
このため、位相差演算部12は、びびり振動が発生している場合には一定値の角度θすなわち位相差を算出し、びびり振動が発生していない場合には不規則な値の角度θすなわち位相差を算出することになる。
ここで、びびり振動の有無による角度θの差異を確認するための実験結果について説明する。本実験では、一般的なCNC(Computer Numerical Control)旋盤のチャックに、炭素鋼鋼材であるS45Cの丸棒を設置し、丸棒の外周を旋削により加工した。なお、工具は超硬工具を用い、工具に加速度センサを設置することで加工中の加速度を計測した。
びびり振動が発生したときの加工条件は、主軸回転数2000(RPM)、1刃あたり送り0.1(mm)、径方向切込み1.0(mm)、チャックからのワーク突き出し量200(mm)である。
図7は、びびり振動が発生した場合の実験結果であり、加工中に生じた加速度の波形を示す図である。図8は、図7に示す実験結果から求めた無次元量の位相差を示す図である。本実験では、フィルタ処理部10、センサ信号処理部11、タイミング信号生成部16および位相差演算部12は、加速度センサの出力する加速度データをセンサ信号200として処理することにより、図8に示す無次元量および角度θが算出された。
まず、フィルタ処理部10は、主軸回転数に基づいて、2000/60=33.3(Hz)と、その高調波から成る切削成分を加速度データから除去した。次に、センサ信号処理部11は、第1の状態量を加速度、第2の状態量を速度とする状態量信号202を演算して出力した。なお、第2の状態量である速度は、加速度を積分することによって得た。
タイミング信号生成部16は、主軸101の角度指令が主軸オリエンテーション時の主軸角度を通過するたびにタイミング信号207を出力した。位相差演算部12は、状態量信号202を正規化し、各タイミング信号207を受け取った各時刻t101~t106における無次元量N101~N106を算出した。ここで、t106>t105>t104>t103>t102>t101である。
図8には、状態空間内における無次元量N101~N106が示されている。角度θ101は、無次元量N101と状態空間の原点と無次元量N102とがなす角度を示し、角度θ102は、無次元量N102と状態空間の原点と無次元量N103とがなす角度を示し、角度θ103は、無次元量N103と状態空間の原点と無次元量N104とがなす角度を示す。角度θ104は、無次元量N104と状態空間の原点と無次元量N105とがなす角度を示し、角度θ105は、無次元量N105と状態空間の原点と無次元量N106とがなす角度を示す。角度θ101~θ105は一定である。以上の実験結果から、びびり振動が発生している場合、状態空間内の無次元量の軌跡は円状になり、各時刻における角度θ101~θ105は一定の値をとることが実験的に示された。
続いて、びびり振動が発生しないときの加工条件は、主軸回転数2000(RPM)、1刃あたり送り0.1(mm)、径方向切込み1.0(mm)、チャックからのワーク突き出し量100(mm)である。
図9は、びびり振動が発生していない場合の実験結果であり、加工中に生じた加速度の波形を示す図である。図10は、図9に示す実験結果から求めた無次元量の位相差を示す図である。本実験においても、フィルタ処理部10、センサ信号処理部11、タイミング信号生成部16および位相差演算部12は、加速度センサの出力する加速度データをセンサ信号200として処理することにより、図10に示す無次元量および角度θが算出された。
まず、フィルタ処理部10は、主軸回転数に基づいて、2000/60=33.3(Hz)と、その高調波から成る切削成分を加速度データから除去した。次に、センサ信号処理部11は、第1の状態量を加速度、第2の状態量を速度とする状態量信号202を演算して出力した。なお、第2の状態量である速度は、加速度を積分することによって得た。
タイミング信号生成部16は、主軸101の角度指令が主軸オリエンテーション時の主軸角度を通過するたびにタイミング信号207を出力した。位相差演算部12は、状態量信号202を正規化し、各タイミング信号207を受け取った各時刻t201~t206における無次元量N201~N206を算出した。ここで、t206>t205>t204>t203>t202>t201である。
図10には、状態空間内における無次元量N201~N206が示されている。角度θ201は、無次元量N201と状態空間の原点と無次元量N202とがなす角度を示し、角度θ202は、無次元量N202と状態空間の原点と無次元量N203とがなす角度を示し、角度θ203は、無次元量N203と状態空間の原点と無次元量N204とがなす角度を示す。角度θ204は、無次元量N204と状態空間の原点と無次元量N205とがなす角度を示し、角度θ205は、無次元量N205と状態空間の原点と無次元量N206とがなす角度を示す。以上の実験結果から、びびり振動が発生していない場合、状態空間内の無次元量の軌跡は不規則であり、各時刻における角度θ201~θ205は不定の値をとることが実験的に示された。
図1の説明に戻る。振動判定部13は、位相差演算部12から位相差である角度θを受け取り、以下に示す手順を用いてびびり振動の判定を行う。振動判定部13は、位相差演算部12が算出した位相差情報203に基づいて、びびり振動の発生有無を判定する。びびり振動が発生していると判定した場合、振動判定部13は、振動判定情報204を指令値補正部14へ出力する。振動判定部13の詳細について説明する。
振動判定部13は、位相差演算部12が出力した位相差情報203を記録する。この記録は、位相差情報203を受け取るたびに繰り返し行われ、位相差情報203に含まれる角度θを時系列データとして記録する。
また振動判定部13は、記録された角度θの差分を計算し、角度θの差分が予め定められた閾値を下回る場合、びびり振動が発生していると判定し、振動判定情報204を指令値補正部14に出力する。ここで振動判定部13が使用する閾値について説明する。なお、ここでは位相差演算部12から角度θが±180度の範囲で出力されることとする。びびり振動が発生していない場合、角度θは±180度の範囲で不定の値をとるため角度θの差分は、±360度の値をとる。また、びびり振動が発生している場合、角度θは一定の値をとるため、角度θの差分はゼロとなる。ただし、実際の角度θには測定誤差などに由来する誤差が含まれる。このため、閾値を、角度θの誤差よりも十分に大きい値とすれば、びびり振動の有無を判定することができる。例えば、角度θの誤差が±10度ある場合、閾値は20度よりも大きい値、例えば30度とすればよい。
或いは、振動判定部13は、繰返し受信する複数の位相差情報203に含まれる角度θの標準偏差を算出し、角度θの標準偏差が予め定められた閾値を下回る場合、びびり振動が発生していると判定する。図11は、図1に示す振動判定部13の処理の第1の例について説明するための図である。振動判定部13は、時刻t3から時刻t4の間に取得された角度θを用いて標準偏差を算出し、びびり振動の有無を判定する。図11に示す第1の例では、角度θの標準偏差は、標準偏差の閾値を上回っている。このため、振動判定部13は、びびり振動が発生していないと判定することができる。また、図12は、図1に示す振動判定部13の処理の第2の例について説明するための図である。図12に示す第2の例では、角度θの標準偏差は、標準偏差の閾値を下回っている。このため、振動判定部13は、びびり振動が発生していると判定することができる。
標準偏差の閾値の定め方について説明する。びびり振動が発生する場合、角度θは一定の値をとるため、角度θの差分はゼロとなる。ただし、実際の角度θには、測定誤差などに由来する誤差が存在する。したがって、閾値に、角度θに含まれる誤差よりも十分に大きい値を設定することで、びびり振動の有無を判定することができる。例えば、角度θの誤差が±10度である場合、標準偏差の閾値は10度よりも大きい値、例えば20度とすることができる。
また、別の例として、振動判定部13は、繰返し受信する複数の位相差情報203に含まれる角度θの標準偏差に代わって、角度θの分散を用いてびびり振動の有無を判定することができる。この場合も標準偏差の場合と同様の方法で分散の閾値を定めることができる。
さらに別の例として、振動判定部13は、角度θの最大値と最小値の差を用いて、びびり振動の判定を行ってもよい。この場合、振動判定部13は、予め定められた長さの期間中に受信した角度θの最大値および最小値の差を算出し、算出した最大値および最小値の差が、予め定められた閾値よりも小さい場合、びびり振動が発生していると判定することができる。
さらに別の例として、振動判定部13は、びびり振動の有無を判定するように予め機械学習が行われた推論モデルを用いて、びびり振動の判定を行ってもよい。図22は、図1に示す振動判定部13の処理の第3の例を実現するための構成例を示す図である。第3の例において、振動判定部13は、情報観測部301と推論部302とを有する。情報観測部301は、予め定められたサンプリング数分の、位相差情報203に含まれる角度θを時系列データとして観測し、推論用データセット303を生成する。推論部302は、びびり振動の有無を出力するように予め機械学習を行った推論モデルに、情報観測部301が生成する推論用データセット303を入力することで、びびり振動の有無を示す振動判定情報204を出力する。
推論部302は、どのようなアルゴリズムの推論モデルを用いてもよい。一例として、ニューラルネットワークを用いた推論モデルについて説明する。図23は、ニューラルネットワークを用いた推論モデルの構成例を示す図である。図23に示す例において、ニューラルネットワークは、n個のニューロンを有する入力層x1,x2,・・・,xnと、m個のニューロンを有する中間層y1,y2,・・・,ymと、1個のニューロンを有する出力層z1とから構成される。なお、図23では、中間層が1層である例を示したが、2層以上の中間層を設けてもよい。
入力層x1,x2,・・・,xnは、中間層y1,y2,・・・,ymと接続されており、中間層y1,y2,・・・,ymは、出力層z1と接続されている。なお、図23に示す入力層と中間層との接続は一例であり、各入力層x1,x2,・・・,xnは、いずれの中間層y1,y2,・・・,ymと接続されてもよい。
図23に示すような3層のニューラルネットワークであれば、複数の入力が入力層x1,x2,・・・,xnに入力されると、その値に重みA1~Aaが掛けられて中間層y1,y2,・・・,ymに入力される。中間層y1,y2,・・・,ymに入力された値は、さらに重みB1~Bbが掛けられて出力層z1に入力され、出力層z1から出力される。なお、aおよびbのそれぞれは自然数であり、図23に示す例では、b=mとなる。出力結果は、重みA1~Aa,B1~Bbの値によって変わる。
図24は、図22に示す振動判定部13が用いる推論モデルを学習する学習装置400の機能構成を示す図である。学習装置400は、データ取得部401と、学習処理部402とを有する。
データ取得部401は、実加工中に取得された位相差情報203と、びびり振動の有無を示すびびり振動有無情報305とを対応づけた学習用データセット304を取得する。びびり振動の有無は、びびり振動が発生したか否かによって異なる数値を用いて表すことができる。例えば、びびり振動が発生した場合「1」、発生しなかった場合「0」といった数値が用いられてよい。なお、学習用データセット304に含まれるびびり振動有無情報305は、例えば、加工後に加工面を評価することによって判定された結果を用いることができる。
学習処理部402は、学習用データセット304に含まれる位相差情報203を入力データとして学習を行う。学習処理部402は、推論モデルの出力がびびり振動有無情報3-5の数値と一致するようにいわゆる教師あり学習を行う。学習処理部402は、入力層に位相差情報203を入力して出力層から出力された値が、びびり振動の有無を示す数値に近づくように、重みA1~Aa,B1~Bbを調整することで、位相差情報203を入力するとびびり振動の有無に対応する出力値を出力する推論モデルを学習する。
学習処理部402は、推論モデルの学習時に用いる教師あり学習手法として誤差逆伝搬法を用いることができる。また、学習処理部402は、推論モデルの汎化性能を向上させるために、学習の際にニューロンをランダムに除外する「dropout」、誤差を監視して学習を早く打ち切る「early stopping」といった手法を用いてもよい。学習処理部402は、学習が完了した推論モデルを学習済み推論モデルとして出力する。学習処理部402が出力する学習済み推論モデルは、図22に示す推論部302において用いることができる。
なお、上記では、推論モデルにニューラルネットワークを用いる場合について説明したが、他の公知の方法、例えばRNN(Recurrent Neural Network)、LSTM(Long Short-Term Memory)、SVM(Support Vector Machine)などを用いて推論モデルを構成してもよい。
指令値補正部14は、振動判定部13から振動判定情報204を受け取った場合、駆動指令部15へ運転指令208の変更を指示する補正信号205を出力する。補正信号205は、現在の主軸速度に対応して予め定められた所定の割合だけ、主軸101の速度指令を変更する信号である。
あるいは、補正信号205は、主軸101の回転数が予め有限要素解析で算出した機械構造の共振周波数の整数分の1と一致するように、主軸101の回転数を変更する信号であってもよい。なお、機械構造の共振周波数は、有限要素解析でなく事前実験により同定した値を用いてもよい。
さらに補正信号205は、送り軸102の速度指令地を所定の割合だけ変更する信号であってもよい。具体的には補正信号205は、送り軸102のモータ電流が予め定められた基準値よりも小さくなるように、送り軸102の速度指令値を減少する信号とすることができる。
以上説明したように、実施の形態1によれば、工作機械100に取り付けられる工具または加工対象物の振動を検出したセンサ信号200に基づいて、複数の種類の状態量の時系列データが生成される。そして、状態空間内において、複数の種類の状態量を示す無次元量がタイミング信号207ごとに生成され、無次元量の位相間の差分である角度θを示す位相差情報203が演算される。この位相差情報203のばらつきに基づいて、工作機械100におけるびびり振動の発生が判定される。
数値制御装置1は、びびり振動の判定にFFTに代表される周波数解析を用いないため、データ取得時間の制約を受けずに、びびり振動の判定に十分な周波数分解能を得ることができる。このため、びびり振動の判定に要する時間、つまり、びびり振動が発生してからびびり振動の発生を判定するまでの時間を短縮することができる。したがって、びびり振動が発生してから、びびり振動を抑制するまでの時間も短縮され、高速にびびり振動を抑制することができる。
また、数値制御装置1では、無次元量である位相差情報203に基づいてびびり振動を判定する。従来の方法としては、例えば、びびり振動の振幅に基づいてびびり振動を判定する方法がある。びびり振動の振幅を用いる場合、加工条件に依存して振幅の大きさが異なるため、加工条件に応じて振幅の閾値を設けなければならない。これに対して、無次元量である位相差情報203に基づいてびびり振動を判定する数値制御装置1では、加工条件に依存せずにびびり振動を判定することができ、加工条件に応じて閾値を設ける必要がなくなる。
なお、実施の形態1にかかる工作機械100は、主軸101に加工対象物を設置する構成としたが、例えばフライス盤および旋削盤に代表されるように、主軸101に工具を設置する場合であっても同様の効果を奏することができる。
また、実施の形態1では、1個のセンサ103を用いる構成としたが、工作機械100の複数箇所にセンサ103を設置してもよい。この場合、設置したすべてのセンサ103に対して、実施の形態1で説明した処理を行い、びびり振動の発生を判定することができる。複数箇所にセンサ103を設置する場合、加工中に複数の箇所でびびり振動が発生しても、びびり振動の判定および抑制を行うことができる。
なお、工作機械100に設置するセンサ103は、2種類以上のセンサ103を用いて振動を検出し、2種類以上の状態量を含むセンサ信号200をフィルタ処理部10に出力してもよい。この場合、センサ信号処理部11は、2種類のセンサ出力をそれぞれ第1の状態量および第2の状態量とすることができる。第2の状態量をセンサ103で直接計測することで、微分または積分に伴う量子化誤差の影響を受けずに、びびり振動の判定を行うことができる。
なお、実施の形態1では、センサ信号処理部11は、2種類の状態量を抽出しているが、3種類以上の状態量を抽出してもよい。
実施の形態2.
実施の形態1では、工具の刃数が1枚の例を説明した。実施の形態2では、工具の刃数が2枚以上の場合にびびり振動の発生を判定するための構成について説明する。
実施の形態1では、工具の刃数が1枚の例を説明した。実施の形態2では、工具の刃数が2枚以上の場合にびびり振動の発生を判定するための構成について説明する。
図13は、実施の形態2にかかる数値制御装置2の機能構成を示す図である。数値制御装置2は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16-1と、工具情報記録部17とを有する。
数値制御装置2は、実施の形態1にかかる数値制御装置1の構成に加えて、工具情報210を出力する工具情報記録部17を有し、数値制御装置1のタイミング信号生成部16に代えて、タイミング信号生成部16-1を有する。以下、数値制御装置1と異なる部分について主に説明する。
工具情報記録部17は、工作機械100に設置された工具に関する情報である工具情報210を記録し、記録している工具情報210をタイミング信号生成部16-1に出力する。工具情報210は、少なくとも工具の刃数を示す情報を含む。さらに、工具情報210は、エンドミル、バイトなどの工具の種類、工具長、工具径などの工具形状を示す情報を含んでもよい。なお、工作機械100に1つの旋削工具が取り付けられている場合、工具情報記録部17は、工具の刃数は1枚であることを示す工具情報210を出力する。また、工作機械100に旋削工具が取り付けられている場合であっても、複数の工具が同時に切削する場合、刃数は工具の数と同数になる。例えば工作機械100が下刃物台と上刃物台とを有し、両方の刃物台に1枚ずつ取り付けられた旋削工具で加工を行う場合、工具情報記録部17は、刃数が2枚であることを示す工具情報210を出力する。
タイミング信号生成部16-1は、駆動指令部15が出力する主軸運転指令206および工具情報記録部17が出力する工具情報210に基づき、主軸角度が予め定められた角度を通過したことを判定し、判定したタイミングでタイミング信号207を位相差演算部12へ出力する。
タイミング信号生成部16-1の処理の詳細を説明する。以下では刃数α枚の工具を、主軸101に取り付ける場合について説明する。αは、2以上の自然数とする。
まず、α=2の場合について説明する。図14は、図13に示すタイミング信号生成部16-1が出力するタイミング信号207と主軸101の角度指令との間の関係の第1の例を示す図である。角度指令は、0度から360度の間の値をとり、360度に達すると0度に戻る信号である。タイミング信号207は、主軸101の角度指令が予め定められた複数の角度である角度ψ1および角度ψ2を通過するたびに出力される。つまり、タイミング信号207は、主軸101が1回転する間に、刃数と同数である2回出力される。角度ψ1は、主軸101が1回転する間に1か所定められていれば任意の角度であってよい。例えば、角度ψ1は、主軸オリエンテーション時の主軸角度とすることができる。さらに、角度ψ1と角度ψ2との間の差は、360度を工具刃数の2で除した値となり、ψ2-ψ1=180度となる。
続いて、αが任意の自然数である場合について説明する。図15は、図13に示すタイミング信号生成部16-1が出力するタイミング信号207と主軸101の角度指令との間の関係の第2の例を示す図である。タイミング信号207は、主軸101の角度指令が予め定められた角度ψβを通過するたびに出力される。βは、2以上α以下の自然数である。タイミング信号207は、主軸101が1回転する間にα回出力される。ψ1は、主軸101が1回転する間に1か所定められていれば任意の角度であってよい。ψ1とψβ―1との間には、以下に示す数式(2)に示す制約がある。言い換えると、タイミング信号生成部16-1は、360度を刃数αで除した角度だけ主軸101が回転する毎にタイミング信号207を出力する。
なお、タイミング信号生成部16-1は、主軸運転指令206に含まれる主軸101の速度指令を用いて、タイミング信号207を生成してもよい。この場合、タイミング信号生成部16-1は、主軸101の速度指令S(rpm)から以下の数式(3)を用いて算出される時間T2ごとにタイミング信号207を生成する。
以上説明したように、実施の形態2によれば、工具情報210に基づいてタイミング信号207が生成されるため、工作機械100に取り付けられる工具の刃数が複数である場合であっても、実施の形態1と同様に、びびり振動の判定に要する時間、つまり、びびり振動が発生してからびびり振動の発生を判定するまでの時間を短縮することができる。
なお、刃数が1である場合、数値制御装置2は、数値制御装置1と同様の処理によってびびり振動の判定および抑制を行うことができる。
実施の形態3.
図16は、実施の形態3にかかる数値制御装置3の機能構成を示す図である。数値制御装置3は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16-2とを有する。数値制御装置3は、数値制御装置1のタイミング信号生成部16に代えて、タイミング信号生成部16-2を有する。以下、数値制御装置1と異なる部分について主に説明する。
図16は、実施の形態3にかかる数値制御装置3の機能構成を示す図である。数値制御装置3は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16-2とを有する。数値制御装置3は、数値制御装置1のタイミング信号生成部16に代えて、タイミング信号生成部16-2を有する。以下、数値制御装置1と異なる部分について主に説明する。
工作機械100は、実施の形態1で説明した構成に加えて、主軸101に取り付けられた角度情報センサ104を有する。角度情報センサ104は、主軸角度を測定するセンサであり、例えばエンコーダ、ポテンショセンサなどを用いることができる。角度情報211は、角度情報センサ104により測定された主軸角度の実測値を少なくとも含む信号である。角度情報センサ104は、生成した角度情報211を数値制御装置3に出力する。
タイミング信号生成部16-2は、角度情報211が含む主軸角度の実測値に基づき、タイミング信号207を生成する。タイミング信号生成部16-2がタイミング信号207を生成する方法は、実施の形態1と同様であるが、主軸101の角度指令に代わり、主軸角度の実測値が用いられる。
なお、角度情報センサ104は、主軸速度の実測値を検出するセンサであってもよい。この場合、角度情報センサ104は、角速度センサである。タイミング信号生成部16-2は、主軸101の角度指令に代わり、主軸速度の実測値に基づいてタイミング信号207を生成することができる。タイミング信号生成部16-2は、初期時刻t0を基準時刻として、主軸速度の実測値Sref(rpm)から、以下に示す数式(4)を用いて算出される時間T3ごとにタイミング信号207を出力する。
以上説明したように、実施の形態3によれば、タイミング信号207は、主軸角度の実測値に基づいて生成される。このため、主軸101の角度指令と実際の角度の間で生じる制御の追従誤差の影響を受けることなく、タイミング信号207を生成することができる。したがって、数値制御装置3は、主軸101の角度指令と実際の角度の間に制御の追従誤差が発生する場合であっても、実施の形態1と同様に、びびり振動の判定に要する時間、つまり、びびり振動が発生してからびびり振動の発生を判定するまでの時間を短縮することができる。
なお、実施の形態3では、主軸角度または主軸速度の実測値を用いて、タイミング信号207を生成する例について説明したが、タイミング信号生成部16-2は、主軸角度指令と主軸角度の実測値とを併せて用いてもよいし、主軸速度指令と主軸速度の実測値とを併せて用いてもよい。
実施の形態4.
図17は、実施の形態4にかかる数値制御装置4の機能構成を示す図である。数値制御装置4は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16-3と、工具情報記録部17とを有する。つまり、数値制御装置4は、実施の形態3で説明した数値制御装置3の構成に加えて、工具情報記録部17を有し、タイミング信号生成部16-2の代わりにタイミング信号生成部16-3を有する。以下、数値制御装置2と異なる部分について主に説明する。
図17は、実施の形態4にかかる数値制御装置4の機能構成を示す図である。数値制御装置4は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16-3と、工具情報記録部17とを有する。つまり、数値制御装置4は、実施の形態3で説明した数値制御装置3の構成に加えて、工具情報記録部17を有し、タイミング信号生成部16-2の代わりにタイミング信号生成部16-3を有する。以下、数値制御装置2と異なる部分について主に説明する。
タイミング信号生成部16-3は、角度情報センサ104が出力する角度情報211と、工具情報記録部17が出力する工具情報210とに基づいて、タイミング信号207を生成し、生成したタイミング信号207を位相差演算部12に出力する。タイミング信号生成部16-3は、実施の形態2のタイミング信号生成部16-1と同様の方法でタイミング信号207を生成することができる。このとき、タイミング信号生成部16-3は、主軸101の角度指令に代わって、主軸角度の実測値を用いて、タイミング信号207を生成する。
以上説明したように、実施の形態4によれば、実施の形態3と同様に、主軸101の角度指令と実際の角度の間で生じる制御の追従誤差の影響を受けることなく、タイミング信号207を生成することができるという効果に加えて、実施の形態2と同様に、刃数が複数である場合にも、実施の形態1と同様に、びびり振動の判定に要する時間、つまり、びびり振動が発生してからびびり振動の発生を判定するまでの時間を短縮することができる。
実施の形態5.
図18は、実施の形態5にかかる数値制御装置5の機能構成を示す図である。数値制御装置5は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16-4とを有する。また、数値制御装置5は、工作機械100を制御するための加工プログラム18から、主軸回転数指令213を含む加工プログラム情報212を取得することができる。以下、数値制御装置1と異なる部分について主に説明する。
図18は、実施の形態5にかかる数値制御装置5の機能構成を示す図である。数値制御装置5は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16-4とを有する。また、数値制御装置5は、工作機械100を制御するための加工プログラム18から、主軸回転数指令213を含む加工プログラム情報212を取得することができる。以下、数値制御装置1と異なる部分について主に説明する。
主軸回転数指令213は、タイミング信号生成部16-4に出力され、加工プログラム情報212は、駆動指令部15に出力される。駆動指令部15は、加工プログラム情報212に従って、工作機械100に運転指令208を出力することができる。
タイミング信号生成部16-4は、加工プログラム情報212から抽出された主軸回転数指令213に基づいて、タイミング信号207を生成する。タイミング信号生成部16-4は、初期時刻t0を基準時刻とし、主軸回転数指令213と以下に示す数式(5)とを用いて算出される時間T4ごとにタイミング信号207を出力する。なお、数式(5)のScmdは、主軸回転数指令213が示す主軸回転数を示す。
以上説明したように、実施の形態5によれば、加工プログラム18に含まれる主軸回転数指令213に基づいて、タイミング信号207を生成することができる。加工プログラム18を使用してタイミング信号207を生成する場合であっても、実施の形態1と同様に、びびり振動の判定に要する時間、つまり、びびり振動が発生してからびびり振動の発生を判定するまでの時間を短縮することができる。
実施の形態6.
図19は、実施の形態6にかかる数値制御装置6の機能構成を示す図である。数値制御装置6は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16-5と、工具情報記録部17とを有する。数値制御装置6は、数値制御装置5の構成に加えて、工具情報記録部17を有し、タイミング信号生成部16-4に代えてタイミング信号生成部16-5を有する。
図19は、実施の形態6にかかる数値制御装置6の機能構成を示す図である。数値制御装置6は、フィルタ処理部10と、センサ信号処理部11と、位相差演算部12と、振動判定部13と、指令値補正部14と、駆動指令部15と、タイミング信号生成部16-5と、工具情報記録部17とを有する。数値制御装置6は、数値制御装置5の構成に加えて、工具情報記録部17を有し、タイミング信号生成部16-4に代えてタイミング信号生成部16-5を有する。
タイミング信号生成部16-5は、工具情報210と、主軸回転数指令213とに基づいて、タイミング信号207を生成し、生成したタイミング信号207を位相差演算部12に出力する。
タイミング信号生成部16-5は、初期時刻t0を基準時刻として、以下に示す数式(6)と、刃数αと、主軸回転数Scmdとを用いて算出される時間T5ごとにタイミング信号207を出力する。
以上説明したように、実施の形態6によれば、実施の形態5の効果に加えて、さらに、工具の刃数が複数である場合にも、びびり振動の発生を判定することが可能になる。
続いて、実施の形態1~6にかかる数値制御装置1~6のハードウェア構成について説明する。数値制御装置1~6の各機能部は、処理回路により実現される。これらの処理回路は、専用のハードウェアにより実現されてもよいし、CPU(Central Processing Unit)を用いた制御回路であってもよい。
上記の処理回路が、専用のハードウェアにより実現される場合、これらは、図20に示す処理回路90により実現される。図20は、実施の形態1~6にかかる数値制御装置1~6の機能を実現するための専用のハードウェアを示す図である。処理回路90は、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ASIC(Application Specific Integrated Circuit)、FPGA(Field Programmable Gate Array)、またはこれらを組み合わせたものである。
上記の処理回路が、CPUを用いた制御回路で実現される場合、この制御回路は例えば図21に示す構成の制御回路91である。図21は、実施の形態1~6にかかる数値制御装置1~6の機能を実現するための制御回路91の構成を示す図である。図21に示すように、制御回路91は、プロセッサ92と、メモリ93とを備える。プロセッサ92は、CPUであり、中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、DSP(Digital Signal Processor)などとも呼ばれる。メモリ93は、例えば、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ、EPROM(Erasable Programmable ROM)、EEPROM(登録商標)(Electrically EPROM)などの不揮発性または揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、DVD(Digital Versatile Disk)などである。
上記の処理回路が制御回路91により実現される場合、プロセッサ92がメモリ93に記憶された、各構成要素の処理に対応するプログラムを読み出して実行することにより実現される。また、メモリ93は、プロセッサ92が実行する各処理における一時メモリとしても使用される。
以上の実施の形態に示した構成は、一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、実施の形態同士を組み合わせることも可能であるし、要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。
1,2,3,4,5,6 数値制御装置、10 フィルタ処理部、11 センサ信号処理部、12 位相差演算部、13 振動判定部、14 指令値補正部、15 駆動指令部、16,16-1,16-2,16-3,16-4,16-5 タイミング信号生成部、17 工具情報記録部、18 加工プログラム、90 処理回路、91 制御回路、92 プロセッサ、93 メモリ、100 工作機械、101 主軸、102 送り軸、103 センサ、104 角度情報センサ、200 センサ信号、201 びびり振動成分信号、202 状態量信号、203 位相差情報、204 振動判定情報、205 補正信号、206 主軸運転指令、207 タイミング信号、208 運転指令、209 動作情報、210 工具情報、211 角度情報、212 加工プログラム情報、213 主軸回転数指令、301 情報観測部、302 推論部、303 推論用データセット、304 学習用データセット、305 びびり振動有無情報、400 学習装置、401 データ取得部、402 学習処理部。
Claims (2)
- 工作機械の主軸および送り軸に運転指令を与える駆動指令部と、
前記工作機械に取り付けられる工具または加工対象物の振動を検出したセンサ信号に基づいて、複数の種類の状態量を生成するセンサ信号処理部と、
前記主軸の回転と同期したタイミングで周期的にタイミング信号を生成して出力するタイミング信号生成部と、
状態空間内において前記複数の種類の状態量を示す無次元量を、前記タイミング信号毎に生成し、前記無次元量の位相間の差分を示す位相差情報を演算する位相差演算部と、
前記位相差情報から推論用データセットを作成する情報観測部と、推論用データセットを入力することでびびり振動の有無を判定する機械学習を行った学習済み推論モデルを持つ推論部とを有し、前記推論部に前記推論用データセットを入力することでびびり振動の発生を判定する振動判定部と、
を備えることを特徴とする数値制御装置。 - 工作機械の主軸および送り軸に運転指令を与える駆動指令部と、前記工作機械に取り付けられる工具または加工対象物の振動を検出したセンサ信号に基づいて、複数の種類の状態量を生成するセンサ信号処理部と、前記主軸の回転と同期したタイミングで周期的にタイミング信号を生成して出力するタイミング信号生成部と、状態空間内において前記複数の種類の状態量を示す無次元量を、前記タイミング信号毎に生成し、前記無次元量の位相間の差分を示す位相差情報を演算する位相差演算部と、を備える数値制御装置から、前記位相差情報を取得し、前記位相差情報と前記工作機械におけるびびり振動の有無とを対応づけた学習用データセットを生成するデータ取得部と、
前記位相差情報を入力するとびびり振動の有無に対応する数値を出力するように推論モデルの機械学習を行う学習処理部と、
を備えることを特徴とする学習装置。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2020555927A JP6987275B2 (ja) | 2019-04-08 | 2020-01-22 | 数値制御装置および学習装置 |
CN202080025983.8A CN113646613A (zh) | 2019-04-08 | 2020-01-22 | 数控装置及学习装置 |
DE112020001795.8T DE112020001795T5 (de) | 2019-04-08 | 2020-01-22 | Numerische steuervorrichtung und lernende vorrichtung |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2019/015327 WO2020208685A1 (ja) | 2019-04-08 | 2019-04-08 | 数値制御装置およびびびり振動の発生判定方法 |
JPPCT/JP2019/015327 | 2019-04-08 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2020208893A1 true WO2020208893A1 (ja) | 2020-10-15 |
Family
ID=68532198
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2019/015327 WO2020208685A1 (ja) | 2019-04-08 | 2019-04-08 | 数値制御装置およびびびり振動の発生判定方法 |
PCT/JP2020/002108 WO2020208893A1 (ja) | 2019-04-08 | 2020-01-22 | 数値制御装置および学習装置 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2019/015327 WO2020208685A1 (ja) | 2019-04-08 | 2019-04-08 | 数値制御装置およびびびり振動の発生判定方法 |
Country Status (4)
Country | Link |
---|---|
JP (2) | JP6605185B1 (ja) |
CN (1) | CN113646613A (ja) |
DE (1) | DE112020001795T5 (ja) |
WO (2) | WO2020208685A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7179198B1 (ja) * | 2021-04-26 | 2022-11-28 | 三菱電機株式会社 | 数値制御装置、学習装置及び、びびり振動の抑制方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115315665A (zh) * | 2020-03-24 | 2022-11-08 | 三菱电机株式会社 | 数控装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020146296A1 (en) * | 2001-02-26 | 2002-10-10 | National Institute Of Standards And Technology | Method and device for avoiding chatter during machine tool operation |
JP2008137102A (ja) * | 2006-11-30 | 2008-06-19 | Matsushita Electric Works Ltd | 工作機械総合監視装置 |
JP2012056051A (ja) * | 2010-09-10 | 2012-03-22 | Makino Milling Mach Co Ltd | びびり振動検出方法及びびびり振動回避方法、並びに工作機械 |
WO2014115395A1 (ja) * | 2013-01-23 | 2014-07-31 | 株式会社日立製作所 | 切削振動抑止方法、演算制御装置、および工作機械 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1100339A (en) * | 1964-08-05 | 1968-01-24 | Rolls Royce | Gas turbine engine fuel system |
JP3217604B2 (ja) * | 1993-08-20 | 2001-10-09 | 三菱電機株式会社 | 位置決め装置 |
JPH08229772A (ja) * | 1995-03-01 | 1996-09-10 | Toyota Motor Corp | 主軸回転数の設定方法および装置 |
JP3701392B2 (ja) * | 1995-06-28 | 2005-09-28 | 旭化成ケミカルズ株式会社 | 押出方法 |
JP3835901B2 (ja) * | 1997-08-28 | 2006-10-18 | オムロン株式会社 | 状態判定装置 |
JP2005214631A (ja) * | 2004-01-27 | 2005-08-11 | Kansai Electric Power Co Inc:The | 状態監視保全装置及び状態監視保全方法 |
JP4703315B2 (ja) * | 2005-08-12 | 2011-06-15 | 国立大学法人名古屋大学 | 機械加工装置の回転数演算装置、機械加工装置のびびり振動評価装置および機械加工装置のびびり振動評価方法 |
US9046463B1 (en) * | 2006-04-21 | 2015-06-02 | University Of Washington | Method for conducting nonlinear electrochemical impedance spectroscopy |
WO2010035634A1 (ja) * | 2008-09-26 | 2010-04-01 | 本田技研工業株式会社 | 周波数成分分析装置 |
JP5742312B2 (ja) * | 2011-03-10 | 2015-07-01 | 株式会社ジェイテクト | びびり振動検出方法 |
CN102243143B (zh) * | 2011-04-20 | 2014-04-30 | 上海斐赛轴承科技有限公司 | 轴承异常声检测控制方法以及轴承振动测量仪 |
CN102314534B (zh) * | 2011-06-03 | 2012-11-07 | 东北大学 | 一种基于振动可靠性和遗传算法的齿轮齿廓修形方法 |
JP5853437B2 (ja) * | 2011-06-24 | 2016-02-09 | 株式会社ジェイテクト | びびり振動検出方法 |
JP5908342B2 (ja) * | 2012-05-17 | 2016-04-26 | オークマ株式会社 | 工作機械の加工振動抑制方法及び加工振動抑制装置 |
JP6021632B2 (ja) * | 2012-12-20 | 2016-11-09 | 三菱重工業株式会社 | 加工装置の制御装置、加工装置、加工装置の制御プログラム、加工装置の制御方法、及び加工方法 |
JP6257481B2 (ja) * | 2014-09-02 | 2018-01-10 | 三菱電機株式会社 | 数値制御装置 |
CN109270486B (zh) * | 2018-07-25 | 2020-09-22 | 中国电子科技集团公司第二十九研究所 | 一种采用两传感器旋转的频变信号二维测向方法 |
-
2019
- 2019-04-08 WO PCT/JP2019/015327 patent/WO2020208685A1/ja active Application Filing
- 2019-04-08 JP JP2019542739A patent/JP6605185B1/ja active Active
-
2020
- 2020-01-22 DE DE112020001795.8T patent/DE112020001795T5/de active Pending
- 2020-01-22 WO PCT/JP2020/002108 patent/WO2020208893A1/ja active Application Filing
- 2020-01-22 CN CN202080025983.8A patent/CN113646613A/zh active Pending
- 2020-01-22 JP JP2020555927A patent/JP6987275B2/ja active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020146296A1 (en) * | 2001-02-26 | 2002-10-10 | National Institute Of Standards And Technology | Method and device for avoiding chatter during machine tool operation |
JP2008137102A (ja) * | 2006-11-30 | 2008-06-19 | Matsushita Electric Works Ltd | 工作機械総合監視装置 |
JP2012056051A (ja) * | 2010-09-10 | 2012-03-22 | Makino Milling Mach Co Ltd | びびり振動検出方法及びびびり振動回避方法、並びに工作機械 |
WO2014115395A1 (ja) * | 2013-01-23 | 2014-07-31 | 株式会社日立製作所 | 切削振動抑止方法、演算制御装置、および工作機械 |
Non-Patent Citations (1)
Title |
---|
HINO, JUNICHI: "A study of chatter prediction in end milling process: Fuzzy neural network model with inputs of cutting conditions and sound", TRANSACTIONS OF THE JAPAN SOCIETY OF MECHANICAL ENGINEERS SERIES C, vol. 69, no. 683, 2003, pages 181 - 186 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7179198B1 (ja) * | 2021-04-26 | 2022-11-28 | 三菱電機株式会社 | 数値制御装置、学習装置及び、びびり振動の抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
JP6987275B2 (ja) | 2021-12-22 |
CN113646613A (zh) | 2021-11-12 |
DE112020001795T5 (de) | 2022-03-24 |
JPWO2020208685A1 (ja) | 2021-05-06 |
JPWO2020208893A1 (ja) | 2021-04-30 |
WO2020208685A1 (ja) | 2020-10-15 |
JP6605185B1 (ja) | 2019-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9211624B2 (en) | Vibration determination method and vibration determination device | |
US9846428B2 (en) | Controller for spindle motor | |
CN101380719B (zh) | 机床的振动抑制装置 | |
JP5160980B2 (ja) | 振動抑制方法及び装置 | |
JP4777960B2 (ja) | 振動抑制装置 | |
Eynian et al. | Chatter stability of general turning operations with process damping | |
WO2020208893A1 (ja) | 数値制御装置および学習装置 | |
JP4433422B2 (ja) | 振動抑制装置 | |
EP2497600B1 (en) | Chattering detection method for machine tool | |
CN111052015B (zh) | 数控***及电动机控制装置 | |
CN103769945B (zh) | 颤振抑制方法和机床 | |
JP6888375B2 (ja) | 振動解析システムおよび加工機 | |
CN109088570A (zh) | 电动机的控制装置 | |
JP4891150B2 (ja) | 工作機械の振動抑制装置 | |
Leal-Muñoz et al. | Accuracy of a new online method for measuring machining parameters in milling | |
JP5631779B2 (ja) | 工作機械の振動抑制方法及び装置 | |
Leal-Muñoz et al. | Identification of the actual process parameters for finishing operations in peripheral milling | |
Yoneoka et al. | Disturbance observer–based in-process detection and suppression of chatter vibration | |
JP5587707B2 (ja) | 振動抑制装置 | |
CN109702554A (zh) | 一种自适应振动塑形铣削颤振抑制方法 | |
JPH0985587A (ja) | 回転切削工具による切削状態の検出装置および検出方法 | |
JP5767931B2 (ja) | 工作機械の振動抑制方法および振動抑制装置 | |
JP7179198B1 (ja) | 数値制御装置、学習装置及び、びびり振動の抑制方法 | |
CN116140658A (zh) | 一种通过超声振动抑制车削中杆类零件自激振动的方法 | |
Schmitz et al. | Special topics in milling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ENP | Entry into the national phase |
Ref document number: 2020555927 Country of ref document: JP Kind code of ref document: A |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 20787855 Country of ref document: EP Kind code of ref document: A1 |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 20787855 Country of ref document: EP Kind code of ref document: A1 |