Method for controlling the quality of industrial processes, in particular laser-welding processes

Information

  • Patent Grant
  • 8571823
  • Patent Number
    8,571,823
  • Date Filed
    Tuesday, March 16, 2010
    14 years ago
  • Date Issued
    Tuesday, October 29, 2013
    11 years ago
Abstract
A system and a method for controlling the quality of an industrial process, of the type that comprises the steps of: providing one or more reference signals for the industrial process; acquiring one or more real signals that are indicative of the quality of said industrial process; and comparing said one or more reference signals with said one or more real signals in order to identify defects in said industrial process. According to the invention, the method moreover comprises the operations of: obtaining a transformed signal from said reference signal; obtaining a transformed signal from said real signal; and calculating energies of said transformed reference signal and said real signal, respectively, said comparison operation comprising; comparing with one another said energies of said transformed reference signal and said transformed real signal, respectively, in order to extract corresponding time-frequency distribution for selected frequency values; calculating energies of said time-frequency distributions; and comparing the energies of said time-frequency distributions with threshold values in order to identify energy values associated to defects.
Description
BACKGROUND OF THE INVENTION

The present invention relates to methods for controlling the quality of an industrial process, comprising the steps of:

    • providing one or more reference signals for an industrial process;
    • acquiring one or more real signals that are indicative of the quality of said industrial process; and
    • comparing said one or more reference signals with said one or more real signals in order to identify defects in said industrial process.


Monitoring of the defects in industrial processes assumes increasing economic importance on account of its impact on the analysis of quality of industrial products. The possibility of obtaining an on-line and automatic assessment of the quality of an industrial process presents many advantages both from the economic point of view and from the standpoint of the speed of the process. Desirable characteristics of the system are hence:

    • on-line and real-time processing; and
    • capacity for accurate recognition of the main defects in production.


Currently the problem of recognition of the quality of an industrial process, and consequently the identification of defects, is approached by an inspection carried out off line by skilled staff, or else using automatic methods, which, by means of sensors, identify only some of the defects listed above, in a way that is far from satisfactory and is moreover sensitive to the different settings of the machine.


There are known methods and systems for controlling the quality in industrial processes, for example applied to on-line monitoring of the laser-welding process, in particular in the case of welding of sheet metal. The control system is able to assess the presence of porosities in the welding area or else, in the case of butt-welded thin sheet metal, the presence of defects due to overlapping or poor jointing of the sheet metal.


The above used systems base quality control on a comparison between the signals detected during the process and one or more predetermined reference signals indicating a good-quality weld. Said reference signals, which usually range in number between two and ten, are arranged starting from a number of samples of good-quality welds. Obviously, said mode of procedure implies the presence of a skilled operator who is able to certify the goodness of the weld at the moment of creation of the reference signals, and involves expenditure in terms of time and sometimes also in terms of waste of material (used for making the samples necessary for obtaining the reference signals). In some cases there are also pre-arranged reference signals indicating a defective weld, this, however, involving additional problems and difficulties.


From the European patent application No. EP-A-1275464 filed in the name of the present applicant, it is known to split into blocks the signal acquired via photodiode, which collects the radiation emitted by a welding spot, calculating the mean of the signal in each block sampled and considering the blocks having a value smaller than or equal to the offset of the photodiode as indicating the presence of a defect. Said method eliminates the need for the reference signal; however, it enables only a very approximate detection of the defects.


SUMMARY OF THE INVENTION

The aim of the present invention is to overcome all the aforesaid drawbacks.


In order to achieve such aim, the object of the present invention is a method for controlling the quality of industrial processes which has the characteristics indicated at the beginning and is further characterized in that said method further comprises the operations of:

    • obtaining a transformed signal from said reference signal;
    • obtaining a transformed signal from said real signal; and
    • calculating energies of said transformed reference signal and said real signal,


said comparison operation comprising:

    • comparing with one another said energies of said transformed reference signal and said transformed real signal for extracting corresponding time-frequency distributions for selected frequency values;
    • calculating energies of said time-frequency distributions; and
    • comparing the energies of said time-frequency distributions with threshold values in order to identify the energy values associated to defects.


In the preferred embodiment, said steps of obtaining a transformed signal from said reference signal and of obtaining a transformed signal from said real signal comprise a filtering operation by means of the application of a discrete wavelet transform (DWT), whilst said operation of comparing said energies of said transformed reference signal and said transformed real signal for obtaining corresponding time-frequency distributions comprises performing a calculation of the conjugate of the Fourier transform of the envelope of the real signal and of the envelope of the normalized signal, to obtain a real conjugated transformed signal and a reference conjugated transformed signal, respectively, as well as comparing the energies of the reference signal and of the real signal, extracting the frequency values for which the energy of the real signal is greater than that of the reference signal.


Of course, a further object of the invention is the system for controlling the quality of industrial processes that implements the method described above, as well as the corresponding computer product, directly loadable into the memory of a computer, such as a processor, and comprises software code portions for performing the method according to the invention when the product is run on a computer.





BRIEF DESCRIPTION OF THE DRAWINGS

Further characteristics and advantages of the invention will emerge from the ensuing description, with reference to the annexed drawings, which are provided purely by way of non-limiting example and in which:



FIG. 1 is a block diagram representing a system implementing the method according to the invention;



FIG. 2 shows a detail of the system of FIG. 1;



FIGS. 3, 4 and 5 are flowcharts representing operations of the method according to the invention; and



FIG. 6 is a diagram of quantities processed by the method according to the invention.





DETAILED DESCRIPTION OF THE INVENTION

The method according to the invention will now be exemplified with reference to a laser-welding method. Said laser-welding method constitutes merely, however, just one non-limiting example of an industrial process to which the method for controlling the quality of industrial processes according to the invention can be applied.


With reference to FIG. 1, the reference number 1 indicates, as a whole, a system for controlling the quality of a laser-welding process. The example relates to the case of two pieces of sheet metal 2, 3, which are welded by means of a laser beam. The number 4 designates as a whole the focussing head, including a lens 5, at which there arrives the laser beam originated by a laser generator (not illustrated) and reflected by a half-reflecting mirror 6, after passage through a lens L. The radiation E emitted by the welding area passes through the half-reflecting mirror 6 and is detected by a sensor 7 constituted by a photodiode that is able to send its outgoing signal to an electronic control and processing unit 8, associated to a personal computer 9.


In a concrete embodiment, the half-reflecting mirror 6 used is a mirror made of ZnSe, with a diameter of 2 ins and a thickness of 5 mm. The sensor 7 is made up of a photodiode with a spectral response of between 190 nm and 1100 nm and an active area of 1.1×1.1 mm and a quartz window.



FIG. 2 illustrates in greater detail the electronic control and processing unit 8 associated to the personal computer 9. Said processing unit 8 comprises an anti-aliasing filter 11, which operates on the signal sent by the sensor 7. There is then envisaged an acquisition card 12 equipped with an analog-to-digital converter, which samples the filtered signal and converts it numerically. Said acquisition card 12 is preferably directly associated to the personal computer 9.


Once again in the case of a concrete embodiment, the acquisition card 12 is a data-acquisition card of the type PC card NI 6110E, with a maximum frequency of acquisition of 5 Msamples/s.


The anti-aliasing filter 11 carries out a filtering of the signal by means of a low-pass filter (for example a Butterworth IIR filter).


In the personal computer 9, according to the invention there is implemented a method for quality control, which is based upon a comparison of a real signal xreal, acquired via the photodiode 7 and a reference signal xref, representing a defective weld, stored in said personal computer 9.


The reference signal, designated by xref(t) is acquired at a frequency of acquisition fs, and hence, according to Nyquist's theorem, has associated to it a frequency band of the signal having the value of fs/2, whilst the number of samples acquired for the reference signal xref(t) is N.



FIG. 3 illustrates a flowchart representing the operations performed on the reference signal xref(t).


In a first step 100, an operation of filtering of the reference signal xref(t) is performed by means of the application of a discrete wavelet transform (DWT). At output from the step 100 there is thus obtained a signal xrefDWT having N/2 samples in the band 0:fs/4.


Subsequently, a Hilbert-transform operation is applied to the signal xrefDWT in a step 101, to obtain a complex analytical signal xrefHIL, having N/2 samples and null negative frequencies.


A normalization operation is applied to said analytical signal xrefHIL, in a step 102, which produces at output a normalized signal xrefnorm.


On said normalized signal xrefnorm there is then performed, in a step 103, an operation of calculation of an envelope of the normalized signal, designated by xrefinvnorm, whilst in a step 104, a fast-Fourier-transform (FFT) operation is applied to said envelope of the normalized signal xrefinvnorm, to obtain a transformed envelope Xrefinvnorm.


Finally, in a step 105, an operation of calculation of the energy of the reference signal, designated by Eref, is performed by applying the relation:

∫|xrefinvnorm(t)|2dt=∫|Xrefinvnorm(f)|2df  (1)


As regards the real signal xreal(t), also this is acquired at a frequency of acquisition fs, and thus, according to Nyquist's theorem, has associated to it a frequency band of the signal having a value of fs/2, whilst the number of samples acquired for the real signal xreal(t) is N.



FIG. 4 illustrates a flowchart representing the operations performed on the real signal xreal(t).


In particular, represented in FIG. 4 is a first step 200, in which an operation of filtering of the real signal xreal(t) is performed by the application of a DWT. At output from step 200, there is thus obtained a signal xrealDWT having N/2 samples in the band 0:fs/4.


A fast-Fourier-transform operation is performed on said signal xrealDWT, in a step 211, to obtain a transformed signal FFTreal, which, subsequently, in a step 212, is normalized, to obtain a transformed normalized signal FFTrealnorm.


In a step 250, an operation of calculation of a mean frequency f0 is performed on the transformed normalized signal FFTrealnorm, according to the relation:

f0=∫f*FFTrealnorm(f)*FFTrealnorm(f)df  (2)


In a step 251, an operation of calculation of a standard deviation B is performed, according to the relation:

B=(∫f2*FFTrealnorm*FFTrealnormdf−f02)1/2  (3)


In a step 252, there are then calculated a lower band F_Sn==(f0−B/2) and an upper band F_Dx=(f0+B/2).


In parallel, in a step 201, a Hilbert-transform operation is applied to the signal xrealDWT, to obtain a complex analytical signal xrealHIL, which has N/2 samples and having null negative frequencies.


In a step 202, a normalization operation is applied to said analytical signal xrealHIL, which produces at output a normalized signal xrealnorm.


On said normalized signal xrealnorm there is then performed, in a step 203, an operation of calculation of the envelope, designated by xrealinvnorm, whilst, in a step 204, a fast-Fourier-transform operation (FFT) is applied to said envelope of the normalized signal xrealinvnorm, to obtain a transformed envelope Xrealinvnorm.


Finally, in a step 205, an operation of calculation of an energy of the real signal Ereal is performed by applying the following relation:

∫|xrealinvnorm(t)|2dt=∫|Xrealinvnorm(f)|2df  (4)


The operations of calculation of the energies Ereal and Eref are performed in a band delimited between the lower band F_Sn and the upper band F_Dx calculated in step 252. In greater detail, the calculation is performed on the band thus delimited, considering frequency steps, for example of one hertz, i.e.,

Fs|Xrealinvnorm(f)|2df∫stepFDX|Xrealinvnorm(f)|2df
Fs|Xrefinvnorm(f)|2df∫stepFDX|Xrefinvnorm(f)|2df


In this way, the operation of calculation of the energies Eref and Ereal produces two respective vectors, namely a vector of energies of the reference signal Energy_Ref_step (1, . . . k), and a vector of energies of the real signal Energy_Real_step (1, . . . k), both comprising k frequency values.


Subsequently, a procedure of calculation of the quadratic time-frequency distributions is performed, illustrated in the flowchart of FIG. 5, which comprises the following operations:

    • in a step designated by 300, calculating the conjugates of the fast-Fourier transforms (FFTs) of the envelope of the real signal Xrealinvnorm(f) and of the envelope of the reference signal Xrefinvnorm(f), to obtain conjugated transformed signals, namely, the real conjugated transformed signal X*realinvnorm(f) and the reference conjugated transformed signal X*refinvnorm(f), respectively;
    • in a step 301, considering the energies of the reference signal Eref and of the real signal Ereal, represented, respectively, by the vector of energies of the reference signal Energy_Ref_step (1, . . . k) and the vector of energies of the real signal Energy_Real_step (1, . . . k), and, for each element k of said two vectors, assessing whether the following criterion is satisfied:

      Energy_Real_step (1, . . . k)>Energy_Ref_step (1 . . . k)  (5)


This operation can be appreciated also with reference to the graph of FIG. 6, which shows the amplitudes of the energy of the reference signal Eref and of the energy of the real signal Ereal (thicker line) as a function of the frequency;

    • if the criterion (5) is satisfied, in a step 302 carrying out an operation of extraction of the frequency value for which said criterion (5) is verified, said value being indicated as f_e; according to how many times the condition is satisfied, obtaining up to a maximum of k frequency values f_e; FIG. 6 shows the regions corresponding to the frequency values f_e for which the criterion (5) is satisfied;
    • in a step 303, constructing a matrix M, the rows of which are represented by the extracted frequency values f_e, whilst the columns are represented by N/2 temporal values t1 . . . tN/2 of the signal at output from the DWT operation 200;
    • in a step 304, for each row of the matrix M, calculating a quadratic time-frequency distribution both for the reference signal, designated by Tfdref, and for the real signal, designated by Tfdreal, using the Margenau-Hill relation, namely

      Tfdreal=Real(xrealDWT(tXrealinvnorm*(f)·and−j2πf)  (6)
      Tfdref=Real(xrefDWT(tXrefinvnorm*(f)·and−j2πf)  (7)
    • then, in a step 305, calculating, for both the reference signal and the real signal, energies associated to the distributions for each instant of time, designated as Etref and Etreal, respectively; and
    • then, in a step 306, calculating a maximum value of the energy max_Tfdref for the time-frequency distribution of the reference Tfdref.


Finally, in order to obtain an estimate of the defects, in a step 307 each temporal value of the energy Etreal of the quadratic time-frequency distribution of the real signal Tfdreal is compared with the maximum value of the energy max_Tfdref.


If said value of energy of the quadratic time-frequency distribution of the real signal Tfdreal exceeds the maximum value of the energy max_Tfdref, this means that there is a defect on that time co-ordinate.


In this way, it is possible to locate the defects in time.


Of course, without prejudice to the principle of the invention, the details of construction and the embodiments may vary widely with respect to what is described and illustrated herein purely by way of example, without thereby departing from the scope of the present invention.

Claims
  • 1. A system for controlling quality of a manufacturing process for recognizing defects in result of a production, comprising: sensor which detects one or more process parameters; andan electronic control and processing unit configured to process the signals issued by said sensor and configured to control the quality of the manufacturing process,wherein said electronic control and processing unit is configured to acquire, by detecting one or more quantity of the production of said manufacturing process, one or more real signals produced by a machine involved in the manufacturing process, where the one or more real signals are indicative of the quality of said manufacturing process;wherein said electronic control and processing unit comprises storing means for storing one or more reference signals corresponding to said one or more real signals and indicative of a good-quality process for the manufacturing process and is configured to provide from the storing means said one or more reference signals corresponding to said one or more real signals and indicative of a good-quality process for the manufacturing process;wherein said electronic control and processing unit is further configured by a computer-executable software code loaded in a memory of said electronic control and processing unit to obtain a reference transformed signal from said reference signal storied in the storing means and to acquire a transformed real signal from said real signal, and is configured to convert each of the transformed real signal and the transformed reference signal into a respective vector of energies, where the vector of energies comprises an integer number k of frequency values; andwherein said electronic control and processing unit is further configured by a computer-executable software code loaded in the memory of said electronic control and processing unit to compare said one or more reference signals with said one or more real signals to identify the defects in the result of the production, wherein said electronic and control processing unit is configured to compare by: comparing said vector of energies of said transformed reference signal with said vector of energies of said transformed real signal and selecting frequency values among said k frequency values on the basis of said comparison, respectively, for extracting corresponding time-frequency distributions for said selected frequency values;calculating energies of said time-frequency distributions; andcomparing the energies of said time-frequency distributions with threshold values in order identify energy values associated with said defects, and the corresponding time coordinates of said defects.
  • 2. The system according to claim 1, wherein said industrial process is a laser-welding process.
  • 3. A system for controlling quality of a manufacturing process for recognizing defects in result of a production, comprising: means for detecting one or more process parameters; andmeans for processing the signals issued by said detecting means,wherein said processing means are means for controlling the quality of the manufacturing process such that: said processing means are means for acquiring, by detecting one or more quantity of the production of said manufacturing process, one or more real signals produced by a machine involved in the manufacturing process, where the one or more real signals are indicative of the quality of said manufacturing process;said processing means are means for providing one or more reference signals corresponding to said one or more real signals and indicative of a good-quality process for the manufacturing process;said processing means are means for obtaining a transformed signal from said reference signal; for obtaining a transformed signal from said real signal; for converting each of the transformed real signal and the transformed reference signal into a respective vector of energies, where the vector of energies comprises an integer number k of frequency values; andsaid processing means are means for comparing said one or more reference signals with said one or more real signals in order to identify the defects in the result of the production,wherein said comparison means are means for comparing said vector of energies of said transformed reference signal with said vector of energies of said transformed real signal and selecting frequency values among said k frequency values on the basis of said comparison, respectively, for extracting corresponding time-frequency distributions for said selected frequency values, means for calculating energies of said time-frequency distributions, and means for comparing the energies of said time-frequency distributions with threshold values in order identify energy values associated with said defects, and the corresponding time coordinates of said defects.
  • 4. The system according to claim 3, wherein said industrial process is a laser-welding process.
  • 5. An apparatus which controls quality of a manufacturing process and which recognizes defects in result of a production, the apparatus comprising: a memory which stores process parameters detected by at least one sensor and software modules that process the detected parameters and that control the quality of the manufacturing process; anda processor which executes the software modules,wherein the software modules comprise: a module configured for acquiring, by detecting one or more quantity of the production of said manufacturing process, one or more real signals produced by a machine involved in the manufacturing process, where the one or more real signals are indicative of the quality of said manufacturing process;a module configured for one or more reference signals corresponding to said one or more real signals and indicative of a good-quality process for the manufacturing process;a module configured for obtaining a transformed signal from said reference signal;a module configured for obtaining a transformed signal from said real signal;a module configured for converting each of the transformed real signal and the transformed reference signal into a respective vector of energies, where the vector of energies comprises an integer number k of frequency values; anda module configured for comparing said one or more reference signals with said one or more real signals in order to identify the defects in the result of the production, wherein the comparing module is further configured for:comparing said vector of energies of said transformed reference signal with said vector of energies of said transformed real signal and for selecting frequency values among said k frequency values on the basis of said comparison, respectively, for extracting corresponding time-frequency distributions for said selected frequency values;calculating energies of said time-frequency distributions; andcomparing the energies of said time-frequency distributions with threshold values in order identify energy values associated with said defects, and the corresponding time coordinates of said defects.
  • 6. A system for controlling quality of a manufacturing process for recognizing defects in result of a production, comprising: sensor configured to detect one or more process parameters; andan electronic control and processing unit configured to process the signals issued by said sensor and configured to control the quality of the manufacturing process,wherein the processing unit is configured to acquire, by detecting one or more quantity of the production of said manufacturing process, one or more real signals produced by a machine involved in the manufacturing process, where the one or more real signals are indicative of the quality of said manufacturing process;wherein the processing unit is configured to provide one or more reference signals corresponding to said one or more real signals and indicative of a good-quality process for the manufacturing process;wherein the processing unit is configured to obtain a transformed signal from said reference signal; to obtain a transformed signal from said real signal; to convert each of the transformed real signal and the transformed reference signal into a respective vector of energies, where the vector of energies comprises an integer number k of frequency values; andwherein the processing unit is configured to compare said one or more reference signals with said one or more real signals in order to identify the defects in the result of the production, wherein the processing unit is configured to compare by: comparing said vector of energies of said transformed reference signal with said vector of energies of said transformed real signal and selecting frequency values among said k frequency values on the basis of said comparison, respectively, for extracting corresponding time-frequency distributions for said selected frequency values;calculating energies of said time-frequency distributions; andcomparing the energies of said time-frequency distributions with threshold values in order identify energy values associated with said defects, and the corresponding time coordinates of said defects.
Priority Claims (1)
Number Date Country Kind
TO2004A0013 Jan 2004 IT national
CROSS REFERENCE PARAGRAPH

This is a Divisional of application Ser. No. 11/033,498, filed on Jan. 12, 2005 now U.S. Pat. No. 7,728,254, which claims the benefit of priority from Italian Patent Application No. TO2004-A000013, filed on Jan. 13, 2004. The entire disclosures of the prior applications are considered part of the disclosure of the accompanying divisional application and are hereby incorporated by reference in their entirety.

US Referenced Citations (27)
Number Name Date Kind
3659974 Neugroschl May 1972 A
4813288 Gamberini Mar 1989 A
4877940 Bangs et al. Oct 1989 A
5121339 Jenuwine et al. Jun 1992 A
5410406 Webster Apr 1995 A
5486677 Maischner et al. Jan 1996 A
5517420 Kinsman et al. May 1996 A
5553614 Chance Sep 1996 A
5651903 Shirk Jul 1997 A
5659479 Duley et al. Aug 1997 A
5698120 Kurosawa et al. Dec 1997 A
5874708 Kinsman et al. Feb 1999 A
5986252 Kawamura Nov 1999 A
6181888 Scheuer et al. Jan 2001 B1
6335504 Ling et al. Jan 2002 B1
6400741 Matsunaga et al. Jun 2002 B1
6549022 Cole, Jr. et al. Apr 2003 B1
6661820 Camilleri et al. Dec 2003 B1
6674044 Bolognese et al. Jan 2004 B2
6914215 Davis et al. Jul 2005 B2
7107118 Orozco et al. Sep 2006 B2
7124063 D'Angelo et al. Oct 2006 B2
7440870 D'Angelo et al. Oct 2008 B2
7640125 D'Angelo et al. Dec 2009 B2
7702476 D'Angelo et al. Apr 2010 B2
20020144984 Mori et al. Oct 2002 A1
20090001262 Visser et al. Jan 2009 A1
Foreign Referenced Citations (6)
Number Date Country
19962967 Jul 2001 DE
1238744 Sep 2002 EP
1275464 Jan 2003 EP
1371443 Dec 2003 EP
40-9057476 Mar 1997 JP
2003-053569 Feb 2003 JP
Related Publications (1)
Number Date Country
20100228377 A1 Sep 2010 US
Divisions (1)
Number Date Country
Parent 11033498 Jan 2005 US
Child 12724724 US