[go: up one dir, main page]

WO2013031353A1 - 加工異常検知方法および加工装置 - Google Patents

加工異常検知方法および加工装置 Download PDF

Info

Publication number
WO2013031353A1
WO2013031353A1 PCT/JP2012/066102 JP2012066102W WO2013031353A1 WO 2013031353 A1 WO2013031353 A1 WO 2013031353A1 JP 2012066102 W JP2012066102 W JP 2012066102W WO 2013031353 A1 WO2013031353 A1 WO 2013031353A1
Authority
WO
WIPO (PCT)
Prior art keywords
cutting
threshold value
cutting force
amount
harmonic
Prior art date
Application number
PCT/JP2012/066102
Other languages
English (en)
French (fr)
Inventor
中須 信昭
英明 小野塚
Original Assignee
株式会社日立製作所
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社日立製作所 filed Critical 株式会社日立製作所
Priority to US14/126,198 priority Critical patent/US20140288882A1/en
Priority to JP2013531140A priority patent/JP5740475B2/ja
Publication of WO2013031353A1 publication Critical patent/WO2013031353A1/ja

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/406Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by monitoring or safety
    • G05B19/4065Monitoring tool breakage, life or condition
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/37Measurements
    • G05B2219/37242Tool signature, compare pattern with detected signal
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/37Measurements
    • G05B2219/37355Cutting, milling, machining force
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/50Machine tool, machine tool null till machine tool work handling
    • G05B2219/50203Tool, monitor condition tool

Definitions

  • the present invention relates to a method and a processing apparatus for monitoring a processing state during machining and detecting an abnormality.
  • Machining is a general processing method used for various metal processing, and a cutting blade attached to a rotary tool is cut into a work material and processed into various shapes by removing the material.
  • the removal amount increases, so that the cutting efficiency, the feed speed, and the tool rotation speed are increased to increase the efficiency.
  • Patent Document 1 Japanese Patent Application Laid-Open No. Hei 10 (1999) discloses a method in which a change pattern of the motor drive current value is grasped in advance by experiments and simulations, and a threshold value is set for each machining path from this change pattern. 5-337790).
  • the method of setting a threshold value in advance for each machining pass is applicable only when the cutting depth in a single machining pass is constant, and is applicable when the cutting load changes and the machining load changes. Can not. Further, in the processing of a complicated three-dimensional shape, it is necessary to divide a large number of short processing paths, and it is difficult to set a threshold value for each processing path.
  • An object of the present invention is to provide a method capable of dynamically determining a cutting force abnormality detection threshold value even in a machining path in which the cutting amount changes from moment to moment.
  • the present application includes a plurality of means for solving the above problems.
  • the signal generated by the rotary cutting is measured, and the cutting force component including the fundamental wave and the high frequency is extracted from the measured signal.
  • Calculates the threshold value for abnormality detection based on the ratio of the fundamental wave and high frequency of the cutting force component calculates the cutting force based on the cutting force component, and compares the cutting force with the threshold value to determine machining abnormality. To do.
  • the cutting force abnormality detection threshold value can be dynamically determined according to the change of the cutting amount, the setting accuracy of the abnormality detection threshold value is improved and the machining accuracy is improved. be able to.
  • FIG. 2 shows an apparatus configuration of a general machining apparatus used in this embodiment.
  • a three-axis control machining apparatus will be described as an example, but the number of control axes and the apparatus configuration are not limited thereto.
  • the machining apparatus 100 holds a casing 101, a machining tool 104, a spindle 103 that holds and rotates the machining tool 104, a spindle stage 102 that moves the spindle 103 in the axial direction, a workpiece 105, and a workpiece.
  • a table 106 to be moved and a controller 107 for controlling the machining apparatus 100 are included.
  • the MPU (not shown) in the controller, by executing software, performs a frequency conversion unit, a cutting force component extraction unit, a cutting force calculation unit, an abnormality determination unit, a cutting amount calculation unit, an abnormality detection threshold, which will be described later.
  • the memory (not shown) functions as a calculation unit, and includes a machining condition storage unit, a cutting amount conversion coefficient storage unit, and a threshold conversion coefficient storage unit.
  • the machining apparatus 100 is configured to machine the shape of the work material 105 by rotating the work tool 104 to cut into the work material 105 and removing the work material 105.
  • the machining tool 104 vibrates the machining tool 104, the casing 101, and the like, causing problems such as a reduction in surface roughness of the machining surface and breakage of the machining tool 104.
  • Fig. 1 shows the processing flow of the machining abnormality determination method.
  • cutting state quantity measurement (S1) is implemented, and frequency conversion (S2) of the measured signal is implemented.
  • the cutting force component extraction (S3) is performed, and the cutting amount calculation (S4) is performed from the extracted signal.
  • the cutting force component extracted by the cutting force component extraction (S3) is subjected to inverse frequency conversion to calculate the cutting force (S6).
  • abnormality determination (S7) is performed in which an abnormal state is determined by comparing the cutting force calculated in cutting force calculation (S6) with the threshold calculated in abnormality detection threshold calculation (S5).
  • a cutting state quantity is measured using a sensor (not shown).
  • a sensor output such as a force sensor signal, a driving current value of a spindle motor, an acceleration sensor signal, an acoustic signal, and acoustic emission can be generally used.
  • the force sensor can be installed by being built in the table 106 or the spindle stage 102 or disposed so as to be sandwiched between the work material 105 and the table 106. Since the drive current value of the spindle motor becomes a value proportional to the force for rotating the machining tool 104, the machining load can be measured.
  • the acceleration sensor and the acoustic emission are mainly attached to the casing 101, the spindle stage 102, and the table 106, and measure the vibration of the apparatus.
  • the acoustic signal is for collecting sounds generated by the vibration of the apparatus with a microphone or the like.
  • the processing tool 104 has a structure in which a tip 121 having a cutting edge is attached to a rotating shaft 122.
  • the processing tool 104 rotates about the rotation center C, and the chip 121 is cut into the work material 105 and processed.
  • 3 and 4 show an example in which two chips 121 are attached, the number of chips may be different depending on the tool.
  • the axial direction used for signal analysis is the axial cutting direction (perpendicular to the paper surface), the feed direction of the machining tool 104, and the radial cutting direction perpendicular to them.
  • the tool feed direction X is a substantially constant direction and the moving average line 32 of the locus 31 drawn by the rotation center of the rotating shaft 122 is substantially a straight line, the tool feed direction X Can be fixed. Also, as shown in FIG.
  • the measurement signal may be coordinate-transformed so that the tangential direction is Fx and the perpendicular direction is Fy.
  • abnormality detection it is not always necessary to determine abnormality in three directions, and it is sufficient to make determination using a signal component Fy in a representative direction, for example, a radial cutting direction. Or you may determine by the signal component of the direction where the fluctuation
  • the direction in which the variation in the amount of cutting state appears significantly is determined by the attachment angle of the tip 121, the tool movement direction, and the like.
  • the frequency conversion unit in the controller 107 performs frequency conversion on the measured cutting state quantity measurement value.
  • the frequency conversion method general techniques such as discrete Fourier transform and fast Fourier transform can be used.
  • the cutting force component extraction (S3) the cutting force component extraction unit in the controller 107 extracts a frequency component related to the cutting force. Taking the force sensor output as an example, the measured signal includes a cutting force generated when the workpiece is removed and a vibration force generated by tool vibration or the like.
  • the vibration force frequency determined by the natural frequency of the machining tool 104 can be separated. That is, in the cutting force component extraction (S3), the rotational speed of the machining tool is calculated based on the rotational speed of the spindle motor, and the frequency corresponding to the value obtained by multiplying the rotational speed by the number of blades is used as the fundamental wave. Then, the fundamental wave and the frequency near the integer multiple are extracted from the measured signal as a cutting force component.
  • the cutting amount calculation unit in the controller 107 calculates the cutting amount in the radial direction. This will be described with reference to FIGS. 5A, 5B, 5C, 6A, 6B, and 6C.
  • 5A to 5C show an example in which the diameter cutting amount h is small, and the diameter cutting amount h is about the same as the radius of the machining tool 104.
  • FIG. 5A to 5C show an example in which the diameter cutting amount h is small, and the diameter cutting amount h is about the same as the radius of the machining tool 104.
  • FIG. 5B shows an example of the cutting force signal when the tool is rotated at a rotational speed of 3300 min ⁇ 1 .
  • cutting force is generated at intervals of 0.009 seconds, and there is time for the tip 121 to idle, so intermittent cutting force is applied.
  • the result of discrete Fourier transform of FIG. 5B is shown in FIG. 5C.
  • the fundamental wave is a frequency of 110 Hz (3300 min ⁇ 1 / 60 ⁇ 2 blades) corresponding to the tool rotation speed of 3300 min ⁇ 1 , and harmonics that are an integral multiple of the fundamental wave are generated. Harmonics are generated because the cutting force is intermittent and there are discontinuous portions.
  • FIG. 6 shows an example in which the diameter cutting amount h is large, and the diameter cutting amount h is equal to the diameter of the machining tool 104. Since there is no free running time of the tip 122, the cutting force is continuous. It can be seen that only the signal of the fundamental wave 110 Hz is generated in the frequency conversion result.
  • FIG. 5B is a waveform obtained by removing the waveform of the idle running period of the chip 121 from the graph of FIG. 6B. Therefore, the waveform shown in FIG. 5B can be obtained by applying a window function that makes the signal effective only during the time when the chip 121 is cut into the work material 105.
  • a method for deriving a relational expression between the cutting depth h and the cutting waveform and Fourier transform will be described with reference to FIGS. 7A, 7B, and 7C.
  • Fig. 7A shows the window function.
  • the window function is a rectangular wave having a magnitude of 1, with a period of fc and a rectangular wave width of s ⁇ fc.
  • the rectangular ratio s is a value related to the idle time of the chip 121, and takes a value of 0 ⁇ s ⁇ 1.
  • FIG. 7B shows a cutting force waveform with the same diameter cutting amount as FIG. 6B.
  • the maximum cutting force is F
  • the period is fc as with the window function.
  • FIG. 7C is a waveform obtained by multiplying the window function (FIG. 7A) and the cutting force waveform (FIG. 7B), and corresponds to the waveform of FIG. 5B.
  • Equation 1 The window function M (t) in FIG. 7A is expressed by Equation 1.
  • 2 ⁇ fc.
  • Equation 2 is a mathematical expression of the cutting force waveform when two chips 121 are attached to the rotary shaft 122 at equal intervals, and is determined by the number of chips, the interval between the chips, and the size of the rotary shaft.
  • Equation 4 When the radius of the machining tool 104 is r and the number of the chips 121 is N, the relationship between the rectangular ratio s and the diameter cutting amount h is expressed by Equation 4.
  • the magnitude of the harmonic component is a function of the diameter cutting amount h, and the diameter cutting amount h can be calculated from the harmonic ratio.
  • the fundamental frequency corresponding to the number of tool rotations is F0
  • the first harmonic is F1
  • the nth harmonic is Fn.
  • F1 / F0, F2 / F0,..., Fn / F0 are functions of the diameter cutting amount h, and other parameters (for example, the shaft cutting amount, machining tool 104, work material, etc. It can be seen that it is not influenced by the rigidity of 105.
  • the rectangular ratio s is obtained from the actual measurement values P1 / P0 using Equation 7, and the cutting amount h can be calculated from Equation 4.
  • general techniques such as Runge-Kutta method, Euler method, simulation, and the like can be used.
  • the harmonic ratios calculated from Equation 3 are P1s / P0s, P2s / P0s,..., Pns / P0s, and the harmonic ratios calculated from the actually measured values are P1m / P0m, P2m / P0m,.
  • Expression 8 is defined as an error function and the cutting amount h is used as a parameter, the cutting amount h having the smallest error function may be obtained.
  • the rectangular ratio s that minimizes the error function in Equation 8 may be obtained.
  • n only needs to be calculated up to a sufficiently high-order term, and does not need to be calculated up to infinity.
  • a general technique such as Runge-Kutta method, Euler method, simulation, or the like can be used.
  • harmonic ratios P1 / P0, P2 / P0,..., Pn ⁇ P0
  • the harmonic ratio P1m / P0m, P2m / P0m,..., Pnm / P0m
  • the rectangular ratio s that minimizes the error function (Equation 9) is selected. The calculation accuracy can be improved as the number of divisions of the rectangular ratio s is increased.
  • the abnormality detection threshold value calculation (S5) performed by the abnormality detection threshold value calculation unit in the controller 107 will be described.
  • the magnitude F of the cutting force in Equation 3 depends on the rigidity, diameter cutting amount, and shaft cutting amount of the machining tool 104 and the work material 105.
  • the parameters that can be changed during machining are the diameter cutting amount and the shaft cutting amount, so that a threshold value using these two parameters as parameters is provided in a table as shown in FIG. 9A so that it can be referred to.
  • the threshold value for each condition is derived in advance by a simulation or experiment and is set in a table according to the magnitude of the cutting force. Since the diameter cutting amount and the harmonic ratio have the relationship of Equation 3, a table in which the diameter cutting amount is replaced with the harmonic ratio as shown in FIG. 9B may be used.
  • the cutting force F can be obtained from Expression 12, and a value obtained by adding a margin D to F can be used as an abnormality detection threshold value.
  • the cutting force calculation unit in the controller 107 obtains the magnitude of the cutting force by performing inverse Fourier transform on the frequency component of the cutting force separated in the cutting force component extraction (S3).
  • the abnormality determination (S7) the abnormality determination unit in the controller 107 detects the cutting abnormality by comparing the cutting force obtained in S6 with the abnormality detection threshold value obtained in S5.
  • the present embodiment it is possible to provide a method for dynamically setting an abnormality detection threshold in a machining path in which the diameter cutting changes from moment to moment, so that it is possible to avoid generation of defective products due to machining failure. Can contribute to the reduction of manufacturing costs.
  • FIG. 10 is a configuration diagram for explaining an embodiment of a portion related to processing abnormality detection in the controller 107 of the processing apparatus.
  • the MPU of the controller 107 includes a cutting state quantity measuring unit 11 and a frequency converting unit 12, a cutting force component extracting unit 13, a cutting force calculating unit 14, an abnormality determining unit 15, a cutting amount calculating unit 16, and an abnormality detection threshold value calculating unit.
  • the memory includes a machining condition storage unit 18, a cutting amount conversion coefficient storage unit 19, a threshold conversion coefficient storage unit 20, a machining condition input unit 21, a threshold conversion coefficient calculation unit 23, a threshold value.
  • a condition input unit 25 is included.
  • the cutting state quantity measuring unit 11 includes sensors such as a force sensor, a driving current value of a spindle motor, an acceleration sensor, an acoustic sensor, and acoustic emission, and is a means for measuring a change in a signal accompanying a cutting force or mechanical vibration.
  • the force sensor can be installed by being built in the table 106 or the spindle stage 102 or by being disposed so as to be sandwiched between the work material 105 and the table 106. Since the driving current value of the spindle motor is a value proportional to the force applied to the machining tool 104, the machining load can be measured.
  • the acceleration sensor and the acoustic emission are mainly attached to the casing 101, the spindle stage 102, and the table 106, and measure the vibration of the apparatus.
  • the acoustic signal collects sound generated by the vibration of the apparatus with a microphone or the like.
  • the frequency converter 12 is a means for converting the frequency of the sensor signal output from the cutting state quantity measuring unit 11.
  • the frequency conversion method general techniques such as discrete Fourier transform and fast Fourier transform can be used.
  • the cutting force component extraction unit 13 is a means for separating the cutting force component using the natural frequency of the processing tool 104 or the frequency of the cutting force.
  • the cutting amount calculation unit 16 is a means for calculating the diameter cutting amount from the harmonic ratio of the cutting force component separated by the cutting force component extraction unit 13.
  • the cutting amount calculation unit 16 acquires the coefficient of the equation for calculating the diameter cutting amount from the harmonic ratio or the conversion table from the cutting amount conversion coefficient storage unit 19, and calculates the diameter cutting amount. Since the expression for calculating the cutting amount is determined by the number of chips, the interval between the chips, and the rotation axis size, these pieces of information are acquired from the cutting amount conversion coefficient storage unit 19.
  • the abnormality detection threshold value calculation unit 17 uses the information in the machining condition storage unit 18 and the threshold conversion coefficient storage unit 20 and uses the calculation formula or the conversion table from the cut amount calculated by the cut amount calculation unit 16. Means for determining an abnormality detection threshold.
  • the threshold conversion coefficient storage unit 20 the machining conditions set by the machining condition setting unit 23, the cutting amount, and the threshold value are stored in association with each other.
  • the cutting force calculation unit 14 is a means for calculating a cutting force by performing inverse frequency conversion on the cutting force component separated by the cutting force component extraction unit 13. General techniques such as inverse discrete Fourier transform and inverse fast Fourier transform can be used.
  • the abnormality determination unit 15 determines abnormality by comparing the cutting force output from the cutting force calculation unit 14 with the threshold value output from the abnormality detection threshold value calculation unit 17.
  • FIG. 11 is a schematic diagram showing an example of an input screen 1001 for inputting a machining condition setting method.
  • FIG. 12 is a diagram illustrating an embodiment of a file format of the library information described in FIG.
  • the library information includes, for example, a library number 1005 and a library item 1006 such as a spindle rotation speed input method.
  • Display items 1002 are displayed on the input screen 1001 in FIG. 11 based on the library information in FIG. 12, and a condition to be used for each item is selected by pressing a radio button 1003. After selecting all the items, pressing the enter button 1004 terminates the input and stores the selected items in the processing condition storage unit 18.
  • the cutting force component extraction unit 13 extracts the cutting force component using the spindle rotation speed acquired by the controller 107 from the machining apparatus 100.
  • the spindle speed of the program stored in the machining apparatus 100 or the controller 107 is acquired.
  • the machining program is composed of several steps, and it is desirable to acquire the spindle rotation speed for each step.
  • FIG. 13 shows an example of the file information when “Acquire from file” is selected as the axis cutting amount input method.
  • file information for example, a library number 1007, a library first item 1008, and a library second item 1009 are included. Enter the path number or program step number as the first library item, and enter the axis cut amount as the second library item to set the axis cut amount corresponding to each pass or each program step number. Can do.
  • FIG. 14 is a schematic diagram showing an example of an input screen 1040 for inputting an abnormality detection threshold value input method.
  • An input method can be selected with a radio button 1003.
  • FIG. 15 shows an example of the outline of the input screen 1041 that transitions when “Acquire from table” is pressed.
  • the vertical axis of the threshold setting table 1045 is the axis cutting amount
  • the horizontal axis is the harmonic ratio or the diameter cutting amount
  • the horizontal axis is linked to the radio button 1003 selected in FIG.
  • FIG. 15 is an example of a screen when “Acquire from table (harmonic ratio conversion)” is selected in FIG.
  • the number of parameters and the range of the threshold setting table 1045 are determined by the numerical values input to the parameter setting table 1044. For each item, a lower limit value, an upper limit value, and a step amount are input, and when the setting button 1043 is pressed, the number of parameters and numerical values in the threshold setting table 1045 are determined and displayed according to the input values. After inputting a numerical value in the threshold value input field 1046, the input is terminated by pressing the enter button 1004.
  • the threshold value and the parameter may be input from a file. In this case, data can be input by designating a file to be read into the threshold setting table 1045 using the file name input unit 1047 and pressing a read button 1048.
  • FIG. 17 shows an example of the input screen 1011 that transitions when “Acquire from cutting force coefficient” is selected as the abnormality detection threshold value input method.
  • setting items 1012 based on the library information of FIG.
  • FIG. 19 shows an example of an input screen that is transitioned to when “Acquire from machining specifications” is selected.
  • setting items 1022 based on the library information of FIG. 20 are displayed, and information is input.
  • the threshold conversion coefficient calculation unit 23 sets the threshold value of the file format item shown in FIG. Then, threshold setting table information in which the input fixed value is set is created and stored in the threshold conversion coefficient storage unit 20.
  • the threshold setting table input in FIG. 15 is stored in the threshold conversion coefficient storage unit.
  • An abnormality detection threshold value is determined by multiplying the calculated cutting force by a threshold setting magnification.
  • a threshold value is calculated while changing the vertical axis and horizontal axis cutting amounts and harmonic ratios shown in the embodiment of FIG. 16, and data including the file information shown in FIG. 16 is created to convert the threshold value.
  • Store in the coefficient storage unit 20 Values stored in advance are used for the lower and upper limits of the values on the vertical axis and the horizontal axis, and for the step. Alternatively, an input screen may be provided.
  • the present embodiment it is possible to provide a means for dynamically setting an abnormality detection threshold in a machining path in which the diameter cutting changes from moment to moment, so that generation of defective products due to machining failure can be avoided. Can contribute to the reduction of manufacturing costs.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)
  • Machine Tool Sensing Apparatuses (AREA)

Abstract

 切削工具を回転させる加工に伴う切削状態量を測定し、測定した信号から基本波及び高調波を含む切削力成分を抽出し、切削力成分の基本波と高調波との比率である高調波比率に基いて異常判定のしきい値を算出し、抽出した切削力成分から切削力を算出し、算出した切削力と前記算出したしきい値とに基いて異常を判定する。

Description

加工異常検知方法および加工装置
 本発明は、機械加工中に加工状態をモニタリングし、異常を検知する方法および加工装置に関するものである。
 機械加工は様々な金属加工に用いられる一般的な加工方法であり、回転工具に取り付けられた切刃を被削材に切り込み、材料を除去することで様々な形状に加工する。複雑な形状を持つ部品を加工する場合、除去量が多くなるため、切込量や送り速度、工具回転速度を大きくする等して、高能率化を図っている。
 切込量や工具回転数を上げると、切刃にかかる力が大きくなるため、工具の振動や切刃の摩耗、折損等の加工トラブルが発生しやすい。加工トラブルが発生すると、加工部分の表面粗さが悪化したり、傷ついたりするため、材料を廃棄しなければならず、廃棄コストがかかる。そこで、加工状態をモニタリングし、異常が発生する直前に加工条件を変更したり、加工を停止することができるシステムを構築する技術が不可欠となっている。
 従来、工具摩耗の検知方法として、主軸回転に用いるモータの駆動電流地量を測定することでモータ負荷を推定し、あらかじめ設定したしきい値と比較することによって異常を検知する方法が知られている。このとき、しきい値の設定方法として、あらかじめ実験やシミュレーションによりモータ駆動電流値の変化パターンを把握しておき、この変化パターンから加工パス毎にしきい値を設定する発明が特許文献1(特開平5-337790)に開示されている。
特開平5-337790号公報
 しかしながら、加工パス毎にあらかじめしきい値を設定する方法は、一つの加工パスにおける切込量が一定のときのみ適用可能であり、切込量が変化して加工負荷が変化する場合には適用できない。また、複雑な3次元形状の加工では、短い加工パスを多数分割する必要があり、それぞれの加工パスに対して、しきい値を設定することは困難である。
 本発明の目的は、切込量が時々刻々と変化する加工パスにおいても、動的に切削力異常検知しきい値を決定可能とする方法を提供することにある。
 上記課題を解決するために、例えば特許請求の範囲に記載の構成を採用する。本願は上記課題を解決する手段を複数含んでいるが、その一例を挙げるならば、回転切削によって発生した信号を測定し、測定した信号から基本波と高周波を含む切削力成分を抽出し、この切削力成分の基本波と高周波との比率に基いて異常検知のしきい値を算出し、切削力成分に基いて切削力を算出し、切削力としきい値とを比較して加工異常を判定する。
 本発明によれば、切込量の変化に応じて動的に切削力異常検知しきい値を決定することができるため、異常検知しきい値の設定精度を向上させるとともに、加工精度を向上させることができる。
本発明の実施形態1の加工異常検知方法を説明するフローチャートである。 本発明の実施例にかかる加工装置の構成を説明するための図である。 径切込量の変化が小さい加工パスにおいて異常判定する方向の決定方法を説明するための図である。 径切込量の変化が大きい加工パスにおいて異常判定する方向の決定方法を説明するための図である。 径切込量が小さい場合の加工状態を示す図である。 切削力を示す図である。 切削力の周波数変換の一例を説明するための図である。 径切込量が大きい場合の加工状態を示す図である。 切削力を示す図である。 切削力の周波数変換の一例を説明するための図である。 径切込量の変化を定式化する方法を説明するための図である。 径切込量の変化を定式化する方法を説明するための図である。 径切込量の変化を定式化する方法を説明するための図である。 切削力の周波数変換後の高調波成分を説明するための図である。 本発明の実施形態1の異常検知しきい値を決定する方法を説明するための図である。 本発明の実施形態1の異常検知しきい値を決定する方法を説明するための図である。 本発明の実施形態1の加工装置の構成を示した図である。 加工条件設定方法を入力する入力画面の一例を示す概略図である。 図11に記載のライブラリ情報のファイルフォーマットの一実施形態を表す図である。 ファイル情報の一実施例を示す図である。 異常検知しきい値入力方法を入力する入力画面の一例を示す概略図である。 遷移する入力画面の概要の一実施例を示す図である。 ファイルフォーマット情報の一実施例を示す図である。 遷移する入力画面の一実施例を示す図である。 ライブラリ情報に基づいた設定項目の一例を示す図である。 遷移する入力画面の一実施例を示す図である。 ライブラリ情報に基づいた設定項目の表示の一例を示す図である。 しきい値変換係数算出部の詳細について説明する図である
 以下、本発明が適用された実施形態の例について、図面を用いて説明する。ただし、以下の説明において、同一構成要素には同一符号を付し繰り返しの説明は省略する。
 図1~図9Cを用いて第一の実施例を説明する。図2に本実施例で用いる一般的な機械加工装置の装置構成を示す。本実施例では3軸制御の機械加工装置を例に説明するが、制御軸数や装置構成はこれに限られるものではない。機械加工装置100は、筐体101と加工工具104,加工工具104を保持して回転させる主軸103、主軸103を軸方向に移動させる主軸ステージ102、被削材105,被削材を保持して移動させるテーブル106、機械加工装置100を制御するコントローラ107で構成される。コントローラ内のMPU(図示せず)は、ソフトウェアを実行することにより、後述する周波数変換部,切削力成分抽出部、切削力算出部,異常判定部、切込量算出部、異常検知しきい値算出部として機能し、メモリ(図示せず)は加工条件記憶部、切込量変換係数記憶部、しきい値変換係数記憶部を有する。機械加工装置100は、加工工具104を回転させて被削材105に切り込み、被削材105を除去することによって、被削材105の形状を加工するものである。加工工具104は被削材105から受ける力により、加工工具104や筐体101等が振動し、加工面の表面粗さが低下したり、加工工具104が折損する等の不具合が発生する。
 図1に加工異常判定方法の処理フローを示す。まず、切削状態量測定(S1)を実施し、測定した信号の周波数変換(S2)を実施する。次に切削力成分抽出(S3)を実施し、さらに抽出した信号から切込量算出(S4)を実施する。次に算出した高調波比率から異常検知しきい値算出(S5)を実施後、切削力成分抽出(S3)で抽出した切削力成分を逆周波数変換して切削力を算出する切削力算出(S6)を実施する。最後に、切削力算出(S6)で算出した切削力と異常検知しきい値算出(S5)で算出したしきい値を比較することによって異常状態を判定する異常判定(S7)を実施する。
 切削状態量測定(S1)では、センサ(図示せず)を用いて、切削状態量を測定する。切削状態量として、一般的に力センサ信号や主軸モータの駆動電流値、加速度センサ信号、音響信号、アコースティックエミッション等のセンサ出力を使用することができる。力センサは、テーブル106や主軸ステージ102に内蔵したり、被削材105とテーブル106の間に挟み込むように配置する等して、設置することができる。主軸モータの駆動電流値は、加工工具104を回転させる力に比例した値となるため、加工負荷を測定することが可能である。加速度センサとアコースティックエミッションは主に筐体101や主軸ステージ102、テーブル106に取り付けられ、装置の振動を測定する。音響信号は装置の振動に伴って発生する音をマイクロフォン等で収集するものである。
 図3、図4を用いて、信号を解析するときの軸方向について説明する。加工工具104は、回転軸122に切刃を形成したチップ121を取り付けた構造である。加工工具104は回転中心Cを中心として回転し、チップ121を被削材105に切り込んで加工する。図3、図4ではチップ121が2枚取り付けられている例を示しているが、チップ枚数は工具によって異なってもよい。
 信号解析に使用する軸方向として、軸切込方向(紙面垂直方向)と加工工具104の送り方向、それらに垂直な径切り込み方向を考える。図3に示した例のように、工具送り方向Xがほぼ一定の方向であり、回転軸122の回転中心が描く軌跡31の移動平均線32がほぼ直線となる場合には、工具送り方向Xを固定して考えることができる。また、図4のように、工具送り方向Xが大きく変化し、回転軸122の回転中心が描く軌跡31の移動平均線32が曲線となる場合は、現在の回転中心位置における移動平均線32の接線方向をFx、垂線方向をFyとなるよう、測定信号を座標変換すればよい。
 異常検知においては、必ずしも3方向の異常判定をする必要はなく、代表的な方向、例えば径切り込み方向の信号成分Fyを用いて判定すれば十分である。あるいは、切削状態量の変動が顕著に表れる方向の信号成分で判定してもよい。切削状態量の変動が顕著に表れる方向は、チップ121の取り付け角度や工具移動方向等によって決まる。
 周波数変換(S2)では、コントローラ107内の周波数変換部が、測定した切削状態量測定値を周波数変換する。周波数変換方法としては、離散フーリエ変換やファーストフーリエ変換等の一般的な技術を用いることができる。切削力成分抽出(S3)では、コントローラ107内の切削力成分抽出部が、切削力に関する周波数成分を抽出する。力センサ出力を例に取ると、測定した信号には、被削材を除去するときに発生する切削力と工具振動等により発生する振動力が混在している。この信号を周波数変換することによって、工具回転数と切刃数から決まる切削力周波数(例えば、切刃2枚の加工工具104を工具回転数3300min―1で回転させると切削力周波数は110Hz(=3300min―1/60×2枚刃)となる)と、加工工具104の固有振動数で決まる振動力周波数に分離することができる。すなわち、切削力成分抽出(S3)では、主軸モータの回転数に基いて加工工具の回転数を算出し、その回転数に刃数を乗じた値に対応する周波数を基本波とする。そして、基本波及びその整数倍付近の周波数を、測定した信号から、切削力成分として抽出する。
 切込量算出(S4)では、コントローラ107内の切込量算出部が、径方向の切込量を算出する。これについて、図5A、図5B、図5C、図6A,図6B、図6Cを用いて説明する。図5A~図5Cは径切込量hが小さい場合であり、径切込量hが加工工具104の半径と同程度な場合の例である。
 図5Bに工具回転数3300min-1で回転させたときの切削力信号の例を示す。工具回転数に対応し、0.009秒間隔で切削力が発生し、チップ121が空走する時間があるため、断続的な切削力がかかっている。図5(b)を離散フーリエ変換した結果を図5Cに示す。工具回転数3300min-1に対応した周波数110Hz(3300min-1/60×2枚刃)を基本波とし、基本波の整数倍の高調波が発生している。高調波は、切削力が断続的であり、不連続部分があるために発生するものである。図6は径切込量hが大きい場合であり、径切込量hが加工工具104の直径と同等な場合の例である。チップ122の空走時間がないため、連続的な切削力となっている。周波数変換結果においても、基本波110Hzの信号のみ発生していることがわかる。
 図6Bの切削力信号は余弦波で近似され、図5Bは図6Bのグラフからチップ121の空走期間の波形を取り除いた波形となっている。したがって、チップ121が被削材105に切り込んでいる時間だけ信号を有効とするような窓関数を図6B波形にかけることによって、図5Bの波形を得ることができる。図7A、図7B、図7Cを用いて切込量hと切削波形の関係式導出方法、およびフーリエ変換について説明する。
 図7Aに窓関数を示す。窓関数は、大きさ1の矩形波であり、周期をfc、矩形波幅をs・fcとおく。矩形比率sはチップ121の空走時間と関係する値であり、0≦s≦1の値をとる。図7Bは図6Bと同じ径切込量での切削力波形を示している。最大切削力をF、周期は窓関数と同じくfcとする。図7Cは窓関数(図7A)と切削力波形(図7B)を掛け合わせた波形であり、図5B波形に相当する。
 図7Aの窓関数M(t)は、式1で表される。以下、簡単のため、ω=2πfcの関係式より、角周波数ωを用いて説明する。
Figure JPOXMLDOC01-appb-M000001
 また、図7Bの切削力波形G(t)は式2で表される。式2は、2枚のチップ121が回転軸122に等間隔に取り付けられている場合の切削力波形を数式化したものであり、チップ枚数や各チップ間の間隔、回転軸寸法によって決まる。
Figure JPOXMLDOC01-appb-M000002
 径切り込みが小さい場合である図7Cの切削力波形H(t)は、式3で表される。
Figure JPOXMLDOC01-appb-M000003
 加工工具104の半径をr、チップ121の枚数をNとおくと、矩形比率sと径切込量hの関係は、式4で表される。
Figure JPOXMLDOC01-appb-M000004
 式3,式4より、高調波成分の大きさは、径切込量hの関数となっており、高調波の比率から、径切込量hを算出することができる。
 高調波比率から径切込量を算出する方法の一例を説明する。図8に示すように、工具回転数に対応した基本周波数をF0、一次高調波をF1,n次高調波をFnとする。式3、式4から、F1/F0,F2/F0、…、Fn/F0は、径切込量hの関数であり、他のパラメータ(例えば、軸切込量や加工工具104、被削材105の剛性)に左右されないことが分かる。基本波F0(t),第一高調波F1(t)を式3から求めると、式5、6となる。
Figure JPOXMLDOC01-appb-M000005
Figure JPOXMLDOC01-appb-M000006
 フーリエ変換のパワースペクトルをそれぞれP0,P1とすると、
P0=|F0(t)|、P1=|F1(t)|であることから、式5,式6より
P1/P0は式7で表される。
Figure JPOXMLDOC01-appb-M000007
 実測値P1/P0から式7を使用して、矩形比率sを求め、式4から切込量hを算出することができる。式7から矩形比率sを求める方法としては、ルンゲクッタ法やオイラー法、シミュレーション等の一般的な技術を用いることができる。
 高調波比率から径切込量を算出する方法の別の一例について説明する。式3から算出される高調波比率をP1s/P0s、P2s/P0s、…、Pns/P0s、とし、実測値から算出される高調波比率をP1m/P0m、P2m/P0m、…、Pnm/P0m、とする。誤差関数として式8を定義し、切込量hをパラメータとしたとき、最も誤差関数が小さくなる切込量hを求めればよい。式4の関係から、式8において最も誤差関数が小さくなる矩形比率sを求めてもよい。また、nは十分高次の項まで計算すればよく、無限大まで計算する必要はない。式8から矩形比率sを求める方法としては、ルンゲクッタ法やオイラー法、シミュレーション等の一般的な技術を用いることができる。
Figure JPOXMLDOC01-appb-M000008
 高調波比率から径切込量を算出する方法の別の一例について説明する。複数個の矩形比率sについて、高調波比率(P1/P0,P2/P0,、、Pn・P0)をシミュレーションや実験によりあらかじめ求めて記憶しておく。次に、実測データの高調波比率(P1m/P0m、P2m/P0m、…、Pnm/P0m、)を求める。最後に、誤差関数(式9)が最も小さくなる矩形比率sを選択する。矩形比率sの分割数を多くするほど、算出精度を向上することができる。
Figure JPOXMLDOC01-appb-M000009
 軸切込量を算出する方法の一例について説明する。切削力の大きさFは、加工工具104と被削材105の剛性で決まる定数Cと軸切込量wを用いて、F=C・wと表される。式3より直流成分はF・s/2であり、C・w・s/2と表される。切削力の直流成分の実測値をLとおくと、式10が得られる。定数Cをあらかじめシミュレーションや実験により求めておけば、直流線分の実測値Lと、式7または式8または式9で求めた矩形比率sを用いて、軸切込量wを式11から算出することができる。
Figure JPOXMLDOC01-appb-M000010
Figure JPOXMLDOC01-appb-M000011
 コントローラ107内の異常検知しきい値算出部が行う異常検知しきい値算出(S5)について説明する。式3における切削力の大きさFは、加工工具104と被削材105の剛性、および径切込量、軸切込量に依存する。この中で加工中に変更可能なパラメータは、径切込量と軸切込量であるため、この二つをパラメータとしたしきい値を図9Aに示すようなテーブルに持たせ、参照できるようにする。各条件におけるしきい値は、あらかじめシミュレーションや実験によって切削力を導出し、その切削力の大きさに応じてテーブルに設定しておく。径切込量と高調波比率は式3の関係があることから、図9Bのように径切込量を高調波比率に置き換えたテーブルを使用してもよい。
 また、式12から切削力Fを求め、このFにマージンDを加えた値を異常検知しきい値とすることができる。
Figure JPOXMLDOC01-appb-M000012
 コントローラ107内の切削力算出部が、切削力算出(S6)では、切削力成分抽出(S3)で分離した切削力の周波数成分を、逆フーリエ変換することで切削力の大きさを求める。コントローラ107内の異常判定部が、異常判定(S7)ではS6で求めた切削力とS5で求めた異常検知しきい値を比較することによって、切削異常を検知する。
 本実施例によれば、径切り込みが時々刻々と変化する加工パスにおいて、動的に異常検知しきい値を設定する方法を提供することが可能なため、加工失敗による不良品発生を回避することができるとともに、製造コストの削減に寄与する。
 図10は、加工装置のコントローラ107うち、加工異常検知にかかる部分の一実施例を説明する構成図である。コントローラ107のMPUは、切削状態量測定部11と周波数変換部12,切削力成分抽出部13、切削力算出部14,異常判定部15、切込量算出部16、異常検知しきい値算出部17として機能し、メモリは、加工条件記憶部18、切込量変換係数記憶部19、しきい値変換係数記憶部20、加工条件入力部21、しきい値変換係数算出部23、しきい値条件入力部25を有する。
 切削状態量測定部11は、力センサや主軸モータの駆動電流値、加速度センサ、音響センサ、アコースティックエミッション等のセンサを備え、切削力や機械振動に伴う信号の変化を測定する手段である。力センサは、テーブル106や主軸ステージ102に内蔵されたり、被削材105とテーブル106の間に挟み込むように配置する等して、設置することができる。主軸モータの駆動電流値は、加工工具104にかかる力に比例した値となるため、加工負荷を測定することが可能である。加速度センサとアコースティックエミッションは主に筐体101や主軸ステージ102、テーブル106に取り付けられ、装置の振動を測定する。音響信号は装置の振動に伴って発生する音をマイクロフォン等で集音するものである。
 周波数変換部12は、切削状態量測定部11から出力されるセンサ信号を周波数変換する手段である。周波数変換方式としては、離散フーリエ変換やファーストフーリエ変換、等の一般的な技術を用いることができる。切削力成分抽出部13は、加工工具104の固有振動数や、切削力の振動数を用いて、切削力成分を分離する手段である。切込量算出部16は、切削力成分抽出部13で分離された切削力成分の高調波比率から径切込量を算出する手段である。切込量算出部16は切込量変換係数記憶部19より、高調波比率から径切込量を算出する式の係数、あるいは、変換テーブルを取得し、径切込量を算出する。切込量を算出する式は、チップ枚数や各チップ間の間隔、回転軸寸法によって決まるため、これらの情報を切込量変換係数記憶部19から取得する。
 異常検知しきい値算出部17は、加工条件記憶部18としきい値変換係数記憶部20の情報を用い、切込量算出部16で算出した切込量から、計算式あるいは変換テーブルを用いて異常検知しきい値を決定する手段である。しきい値変換係数記憶部20には、加工条件設定部23で設定した加工条件と切込量としきい値が関連づけられて記憶されている。
 切削力算出部14は、切削力成分抽出部13で分離された切削力成分を逆周波数変換することにより、切削力を算出する手段である。逆離散フーリエ変換や逆ファーストフーリエ変換等の一般的な技術を用いることができる。異常判定部15は、切削力算出部14から出力される切削力と異常検知しきい値算出部17から出力されるしきい値を比較することによって異常を判定する。
 加工条件入力部21の詳細について、図11~図13を用いて説明する。図11は加工条件設定方法を入力する入力画面1001の一例を示す概略図である。図12は、図11に記載のライブラリ情報のファイルフォーマットの一実施形態を表す図である。ライブラリ情報には、例えばライブラリ番号1005と、主軸回転速度入力方法等のライブラリ項目1006が含まれる。図12のライブラリ情報に基づいて図11の入力画面1001に表示項目1002を表示し、項目毎に使用する条件を、ラジオボタン1003を押下することによって選択する。全ての項目を選択した後、決定ボタン1004を押下することによって、入力を終了し、加工条件記憶部18に選択した項目を記憶する。主軸回転速度入力方法で「装置から取得」を選択した場合には、機械加工装置100からコントローラ107が取得した主軸回転速度を用いて、切削力成分抽出部13で切削力成分を抽出する。また、「プログラムから取得」を選択した場合には、機械加工装置100またはコントローラ107に保存されているプログラムの主軸回転速度を取得する。一般的に加工プログラムは数ステップで構成されており、各ステップ毎に主軸回転速度を取得することが望ましい。軸切込み量入力方法で「ファイルから取得」を選択した場合のファイル情報の一実施例を図13に示す。ファイル情報として、例えばライブラリ番号1007とライブラリ第一項目1008、ライブラリ第二項目1009が含まれる。ライブラリ第一項目として、パス番号、またはプログラムのステップ番号を入力し、ライブラリ第二項目として、軸切込み量を入力することで、各パスまたは各プログラムステップ番号に対応した軸切込み量を設定することができる。
 しきい値条件入力部25の詳細について、図14~図20用いて説明する。図14は異常検知しきい値入力方法を入力する入力画面1040の一例を示す概略図である。ラジオボタン1003により、入力方法を選択できるようにする。「テーブルから取得」を押下したときに遷移する入力画面1041の概要の一実施例を図15に示す。しきい値設定テーブル1045の縦軸は軸切込み量、横軸は高調波比率または径切込み量であり、横軸は図11で選択したラジオボタン1003に連動して、高調波比率または径切込み量の表示が切り替わる。図15は、図11において、「テーブルから取得(高調波比率変換)」を選択したときの画面の一例である。しきい値設定テーブル1045のパラメータ数および範囲は、パラメータ設定テーブル1044に入力された数値で決定される。各項目について、下限値と上限値およびステップ量を入力し、設定ボタン1043を押下すると、入力された値に応じてしきい値設定テーブル1045のパラメータ数と数値が決定され、表示される。しきい値入力欄1046に数値を入力後、決定ボタン1004を押下することによって、入力を終了する。しきい値およびパラメータの入力は、ファイルから読み込む方法でもよい。この場合、ファイル名入力部1047でしきい値設定テーブル1045に読み込むファイルを指定し、読込みボタン1048を押下することで、データを入力することができる。しきい値設定テーブル1045に読み込むファイルのファイルフォーマット情報の一実施例を図16に示す。ファイル情報として、縦軸の項目名と横軸の項目名、縦軸の下限値と上限値、およびステップ、横軸の下限値と上限値、およびステップ、しきい値で構成される。しきい値は、縦軸のステップ数mと横軸のステップ数nの積であるmn個のデータを持つ。異常検知しきい値入力方法として、「切削力係数から取得」を選択した場合に遷移する入力画面1011の一実施例を図17に示す。入力画面1011には、図18のライブラリ情報に基づいた設定項目1012を表示し、必要な情報を入力する。また、「加工諸元から取得」を選択した場合に遷移する入力画面の一実施例を図19に示す。入力画面1021には、図20のライブラリ情報に基づいた設定項目1022を表示し、情報を入力する。
しきい値変換係数算出部23の詳細について、図21を用いて説明する。図14に示す異常検知しきい値入力方法で、「固定値入力」のラジオボタン1003を選択した場合、しきい値変換係数算出部23は、図16に示したファイルフォーマットの項目のしきい値に、入力した固定値を設定したしきい値設定テーブル情報を作成し、しきい値変換係数記憶部20に記憶する。「テーブルから取得」のラジオボタン1003を選択した場合、図15で入力したしきい値設定テーブルをしきい値変換係数記憶部に記憶する。「切削力係数から算出」または「加工諸元から算出」のラジオボタン1003を選択した場合、図17または図19で入力した値をもとにシミュレーションを実施し、工具磨耗量0μmにおける切削力を算出する。算出した切削力にしきい値設定倍率を乗算して、異常検知しきい値を決定する。図16の実施例に示した縦軸、横軸としての軸切込み量と高調波比率を変えながらしきい値を算出し、図16に示すファイル情報を含むデータを作成して、しきい値変換係数記憶部20に記憶する。縦軸、横軸の値の下限値および上限値、ステップはあらかじめ記憶しておいた値を使用する。あるいは、入力画面を設けてもよい。
 本実施例によれば、径切り込みが時々刻々と変化する加工パスにおいて、動的に異常検知しきい値を設定する手段を提供することが可能なため、加工失敗による不良品発生を回避することができるとともに、製造コストの削減に寄与する。
 以上、本発明者によってなされた発明を、前記発明の実施形態に基づき具体的に説明したが、本発明は、前記発明の実施形態に限定されるものではなく、その要旨を逸脱しない範囲において変更可能であることは勿論である。
 101・・・筐体、102・・・主軸テーブル、103・・・主軸、104・・・加工工具、105・・・被削材、106・・・テーブル、107・・・コントローラ、121・・・チップ、122・・・回転軸。

Claims (23)

  1.  切削工具を回転させる加工に伴う切削状態量を測定する測定ステップと、
     前記測定した信号から、基本波及び高調波を含む切削力成分を抽出する抽出ステップと、
     前記切削力成分の基本波と高調波との比率である高調波比率に基いて、異常判定のしきい値を算出するしきい値算出ステップと、
     前記抽出した切削力成分から、切削力を算出する切削力算出ステップと、
     前記算出した切削力と、前記算出したしきい値とに基いて、異常を判定する異常判定ステップと、
     を含む加工異常検知方法。
  2.  請求項1において、
     前記抽出ステップでは、前記測定した信号を周波数変換して切削力成分を抽出し、
     前記切削力算出ステップでは、前記周波数変換されて抽出された切削力成分を、逆周波数変換することにより切削力を算出することを特徴とする加工異常検知方法。
  3.  請求項1または2において、
     前記しきい値算出ステップでは、前記高調波比率に基いて、径切込量を算出し、当該径切込量に基いて、しきい値を算出することを特徴とする加工異常検知方法。
  4.  請求項1乃至3のいずれかにおいて、
     軸切込量を算出する軸切込量算出ステップを有し、
     前記しきい値算出ステップでは、前記高調波比率または径切込み量と、前記軸切込量とに基いて、前記しきい値を設定することを特徴とする加工異常検知方法。
  5.  請求項1乃至4のいずれかにおいて、
     前記測定ステップでは、前記切削状態量として、被加工材の振動、加工装置の振動、前記加工工具を回転させるモータの電流、振動に伴って発生する音のいずれかを測定することを特徴とする加工異常検知方法。
  6.  請求項1乃至5のいずれかにおいて、
     測定した信号を、前記切削工具の回転中心が描く軌跡の移動平均線の接線方向成分と、
    垂直方向成分に座標変換し、前記抽出ステップでは当該垂直方向成分を用いることを特徴とする加工異常検知方法。
  7.  請求項3において、
     前記しきい値算出ステップでは、前記測定信号の、一次高調波の振幅F1と基本波の振幅F0との比である高調波比率と切込量を関連付けて記憶した変換テーブルまたは数式を用いて、前記径切込量を算出することを特徴とする加工異常検知方法。
  8.  請求項7において、
     前記しきい値算出ステップは、前記測定信号の、一次高調波からn次高調波の振幅F1からFnと基本波の振幅F0との比を算出するステップと、シミュレーションあるいは数式から算出される信号の一次高調波からn次高調波の振幅F1からFnと基本波の振幅F0との比を算出するステップと、各高調波比の差を最小化する切込量を算出するステップと
     を有することを特徴とする加工異常検知方法。
  9.  切削工具と、前記切削工具を回転させるモータと、制御を行う制御手段と、を備えた加工装置において、
     切削工具を回転させる加工に伴う切削状態量を測定する測定手段を有し、
     前記制御手段は、
     前記測定した信号から、基本波及び高調波を含む切削力成分を抽出する抽出部と、
     前記切削力成分の基本波と高調波との比率である高調波比率に基いて、異常判定のしきい値を算出するしきい値算出部と、
     前記抽出した切削力成分から、切削力を算出する切削力算出部と、
     前記算出した切削力成分と、前記算出したしきい値とに基いて、異常を判定する異常判定分と、
     とを有することを特徴とする加工装置。
  10.  請求項9において、
     前記抽出部は、前記測定した信号を周波数変換して切削力成分を抽出し、
     前記切削力算出部は、前記周波数変換されて抽出された切削力成分を、逆周波数変換することを特徴とする加工装置。
  11.  請求項9または10において、
     前記しきい値算出部では、前記高調波比率に基いて、径切込量を算出し、当該径切込量に基いて、しきい値を算出することを特徴とする加工装置。
  12.  請求項9乃至11のいずれかにおいて、
     軸切込量を算出する軸切込量算出部を有し、
     前記しきい値算出部は、前記高調波比率または径切込み量と、前記軸切込量とに基いて、前記しきい値を設定することを特徴とする加工装置。
  13.  請求項9乃至12のいずれかにおいて、
     前記測定手段は、前記切削状態量として、被加工材の振動、加工装置の振動、前記加工工具を回転させるモータの電流、振動に伴って発生する音のいずれかを測定することを特徴とする加工装置。
  14.  請求項9乃至13のいずれかにおいて、
     前記しきい値算出部は、前記高調波と前記基本波の比率と切込量を関連付けるテーブルあるいは数式を用いて前記しきい値を算出することを特徴とする加工装置。
  15.  請求項9乃至14のいずれかにおいて、
     前記しきい値算出部は、切込量と加工条件情報と異常検知しきい値を関連付けるテーブルあるいは数式に基いてを用いて前記しきい値を算出する
     をことを特徴とする加工装置。
  16.  請求項9乃至15のいずれかにおいて、
     前記切削工具の回転軸の回転中心が描く軌跡の移動平均線の接線方向と垂線方向に測定値を座標変換する手段を有することを特徴とする加工装置。
  17.  請求項9乃至16のいずれかにおいて、
     加工条件記憶部から加工条件を取得し、シミュレーションまたは数式により、切込量変換係数を算出する手段を有することを特徴とする加工装置。
  18.  請求項15または17において、
     前記加工条件情報として、チップ枚数およびチップ取り付け位置を有することを特徴とする加工装置。
  19.  請求項15または17において、
     前記加工条件情報として、チップ枚数およびチップ取り付け位置を有することを特徴とする加工装置。
  20.  切削工具を回転させる加工に伴う切削状態量を測定し加工異常検知する加工装置におけるデータ入力を支援するデータ入力支援装置であって、
     異常検知しきい値を算出するための加工条件のライブラリ項目と前記加工条件のライブラリ項目をユーザに対して提示しユーザからの加工条件のライブラリ項目の指定を受ける加工条件入力部と、
     異常検知しきい値を算出するためのしきい値のライブラリ項目と前記しきい値のライブラリ項目をユーザに対して提示しユーザからのしきい値のライブラリ項目の指定を受けるしきい値条件入力部と、
     前記ユーザから指定を受けたしきい値から異常検知しきい値を算出するしきい値変換係数算出部と、
     前記しきい値変換係数算出部で算出したしきい値変換係数を記憶するしきい値変換係数記憶部と、
     を備えたことを特徴とするデータ入力支援装置。
  21.  請求項20において、
    前記しきい値変換係数算出部はユーザが入力したしきい値条件からシミュレーションによってしきい値変換係数を算出することを特徴とするデータ入力支援装置。
  22. 請求項20において、
    前記しきい値条件入力部で選択した入力項目に応じて、しきい値変換係数の算出方法を変更することを特徴とするデータ入力支援装置。
  23.  請求項20において、
     前記しきい値変換係数算出部は、高調波比率と軸切込み量、および異常検知しきい値を関連づけたデータを生成することを特徴とするデータ入力支援装置。
PCT/JP2012/066102 2011-09-02 2012-06-25 加工異常検知方法および加工装置 WO2013031353A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US14/126,198 US20140288882A1 (en) 2011-09-02 2012-06-25 Processing Abnormality Detection Method and Processing Device
JP2013531140A JP5740475B2 (ja) 2011-09-02 2012-06-25 加工異常検知方法および加工装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2011-191257 2011-09-02
JP2011191257 2011-09-02

Publications (1)

Publication Number Publication Date
WO2013031353A1 true WO2013031353A1 (ja) 2013-03-07

Family

ID=47755865

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/066102 WO2013031353A1 (ja) 2011-09-02 2012-06-25 加工異常検知方法および加工装置

Country Status (3)

Country Link
US (1) US20140288882A1 (ja)
JP (1) JP5740475B2 (ja)
WO (1) WO2013031353A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013150905A1 (ja) * 2012-04-04 2013-10-10 株式会社日立製作所 切削加工システム及び方法
JP2019072806A (ja) * 2017-10-17 2019-05-16 オムロン株式会社 切削加工装置
DE102019117684A1 (de) 2018-07-18 2020-01-23 Hitachi, Ltd. System zur Bestimmung eines Werkzeugmaschinenzustands und Verfahren zur Bestimmung des Werkzeugmaschinenzustands
JP2020160830A (ja) * 2019-03-27 2020-10-01 ブラザー工業株式会社 数値制御装置、工作機械、制御プログラム、及び記憶媒体

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6578195B2 (ja) * 2015-11-26 2019-09-18 Dmg森精機株式会社 切削工具の固有振動数導出方法及び安定限界曲線作成方法、並びに切削工具の固有振動数導出装置
KR101867136B1 (ko) * 2015-12-07 2018-06-12 현대위아 주식회사 절삭 작업을 위한 공구의 마모 및 파손 감지 방법
CN106863009B (zh) * 2017-01-20 2018-01-09 西北工业大学 基于刀杆两点变形的切削力测量方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62193749A (ja) * 1986-02-19 1987-08-25 Ichiro Inazaki 多刃工具損傷検出装置
JPH05337790A (ja) * 1992-06-01 1993-12-21 Ntn Corp 工具異常検出装置
JPH08323585A (ja) * 1995-05-31 1996-12-10 Fanuc Ltd 異常負荷検出方式
JP2000263377A (ja) * 1999-03-17 2000-09-26 Ntn Corp 金型加工装置
JP2007276031A (ja) * 2006-04-05 2007-10-25 Kitagawa Iron Works Co Ltd 加工対象物把握手段の監視方法および加工対象物把握手段の監視装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3025421B2 (ja) * 1995-06-14 2000-03-27 三菱電機株式会社 制御システムの異常検知装置
KR100579083B1 (ko) * 2002-12-30 2006-05-12 두산인프라코어 주식회사 공작기계의 공구 이상 검출장치 및 검출방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62193749A (ja) * 1986-02-19 1987-08-25 Ichiro Inazaki 多刃工具損傷検出装置
JPH05337790A (ja) * 1992-06-01 1993-12-21 Ntn Corp 工具異常検出装置
JPH08323585A (ja) * 1995-05-31 1996-12-10 Fanuc Ltd 異常負荷検出方式
JP2000263377A (ja) * 1999-03-17 2000-09-26 Ntn Corp 金型加工装置
JP2007276031A (ja) * 2006-04-05 2007-10-25 Kitagawa Iron Works Co Ltd 加工対象物把握手段の監視方法および加工対象物把握手段の監視装置

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013150905A1 (ja) * 2012-04-04 2013-10-10 株式会社日立製作所 切削加工システム及び方法
JP2019072806A (ja) * 2017-10-17 2019-05-16 オムロン株式会社 切削加工装置
US12017316B2 (en) 2017-10-17 2024-06-25 Omron Corporation Cutting processing apparatus
DE102019117684A1 (de) 2018-07-18 2020-01-23 Hitachi, Ltd. System zur Bestimmung eines Werkzeugmaschinenzustands und Verfahren zur Bestimmung des Werkzeugmaschinenzustands
US10990085B2 (en) 2018-07-18 2021-04-27 Hitachi, Ltd. Machine-tool-state determination system and machine-tool-state determination method
DE102019117684B4 (de) * 2018-07-18 2021-06-24 Hitachi, Ltd. System zur Bestimmung eines Werkzeugmaschinenzustands und Verfahren zur Bestimmung des Werkzeugmaschinenzustands
JP2020160830A (ja) * 2019-03-27 2020-10-01 ブラザー工業株式会社 数値制御装置、工作機械、制御プログラム、及び記憶媒体
JP7131454B2 (ja) 2019-03-27 2022-09-06 ブラザー工業株式会社 数値制御装置、工作機械、制御プログラム、及び記憶媒体

Also Published As

Publication number Publication date
JPWO2013031353A1 (ja) 2015-03-23
JP5740475B2 (ja) 2015-06-24
US20140288882A1 (en) 2014-09-25

Similar Documents

Publication Publication Date Title
JP5740475B2 (ja) 加工異常検知方法および加工装置
JP5732325B2 (ja) 振動判別方法、及び振動判別装置
JP5686760B2 (ja) 振動判別方法、及び振動判別装置
JP5793200B2 (ja) 工作機械の切削力検出装置、切削力検出方法、加工異常検出方法、および加工条件制御システム
TWI472402B (zh) Tool flutter monitoring method
JP5609739B2 (ja) 加工びびり振動検出装置、及び工作機械
US20020146296A1 (en) Method and device for avoiding chatter during machine tool operation
JP5622626B2 (ja) 回転速度表示装置
WO2015015987A1 (ja) 軸受装置の振動解析方法、軸受装置の振動解析装置、および転がり軸受の状態監視装置
CN105204443A (zh) 无须进行繁杂的准备工作便可算出能够抑制再生颤振的主轴稳定转数的装置
JP5301380B2 (ja) 回転刃具の寿命予測方法
WO2013088849A1 (ja) 切削加工装置及びそれを用いた加工方法
JP6354349B2 (ja) 振動検出装置と工作機械
JP6511573B1 (ja) 転がり軸受の異常診断方法及び異常診断装置、異常診断プログラム
CN114536103A (zh) 用于检测机床中的刀具状态的装置
WO2022090846A1 (en) Diagnostic apparatus, machining system, diagnostic method, and recording medium
JP5637840B2 (ja) 振動検出方法
JP5660850B2 (ja) 振動表示装置
JP5631792B2 (ja) 工作機械のモニタ装置
KR101865081B1 (ko) 가공 정밀도 향상을 위한 공작 기계의 채터 검출방법
EP4538658A1 (en) Data extraction device and abnormality monitoring device
CN119891857A (zh) 无传感器的颤动检测
JP7538930B2 (ja) 加工状態検出方法、加工状態検出プログラム、および加工状態検出装置
WO2022080505A1 (ja) 作業機械の工具損傷判定方法及びシステム
JP2022145537A (ja) 判定装置、判定方法、プログラムおよび加工システム

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12828534

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2013531140

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 14126198

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 12828534

Country of ref document: EP

Kind code of ref document: A1