METHOD FOR TESTING AND MANUFACTURING SPIRAL SPRINGS FOR A TIMEPIECE

20240069496 · 2024-02-29

    Inventors

    Cpc classification

    International classification

    Abstract

    A method for inspecting a spiral or a spiral blank arranged to form a spiral, the spiral having to exhibit at least one predetermined resonant frequency, the inspecting method including the following steps: a. applying to the spiral or to the spiral blank a time-varying vibrational excitation to cover a predetermined frequency range, b. identifying at least one characteristic of a resonant frequency, such as a resonance peak, of the spiral or of the spiral blank during the vibrational excitation over the predetermined frequency range, c. submitting to a predicting machine the resonant frequency characteristic identified in the step b. to determine a stiffness of the spiral or of the spiral blank and/or determine whether or not a dimensional correction of the spiral or of the spiral blank is necessary to obtain the predetermined resonant frequency.

    Claims

    1. A method for inspecting a spiral or a spiral blank arranged to form a spiral, the spiral having to exhibit at least one predetermined resonant frequency, the inspecting method including the following steps: a. applying to the spiral or to the spiral blank a time-varying vibrational excitation to cover a predetermined frequency range, b. identifying at least one characteristic of a resonant frequency, such as a resonance peak, of the spiral or of the spiral blank during the vibrational excitation over the predetermined frequency range, c. submitting to a predicting machine the resonant frequency characteristic identified in the step b. to determine a stiffness of the spiral or of the spiral blank and/or determine whether or not a dimensional correction of the spiral or of the spiral blank is necessary to obtain the predetermined resonant frequency.

    2. The inspecting method as claimed in claim 1, wherein the frequency range is predetermined to encompass at least one frequency range: centered on the predetermined resonant frequency, and of an extent of at least 30% of the predetermined resonant frequency.

    3. The inspecting method as claimed in claim 1, the spiral having at least two predetermined resonant frequencies, wherein the frequency range is predetermined to cover at least the two predetermined resonant frequencies.

    4. The inspecting method as claimed in claim 1, wherein the step b is based on a measurement over time of an amplitude or of a speed or of an acceleration of displacement of at least one point of the spiral or of the spiral blank, preferably carried out at least partially during the step a.

    5. The inspecting method as claimed in claim 1, the spiral or the spiral blank being contained in a base plane, wherein the step b comprises: a step b of measuring an amplitude or a speed or an acceleration of displacement of at least one point of the spiral or of the spiral blank along a direction normal to the base plane, and/or a step b of measuring an amplitude or a speed or an acceleration of displacement of at least one point of the spiral or of the spiral blank along a direction contained in the base plane.

    6. The inspecting method as claimed in claim 4, wherein the step b comprises: a step of identifying a resonance peak of the spiral or of the spiral blank as a function of an amplitude or of a speed of displacement of at least one point of the spiral or of the spiral blank.

    7. The inspecting method as claimed in claim 6, wherein the resonant frequency is identified on the basis of the width of the resonance peak, at mid-height of the maximum value of the resonance peak.

    8. The inspecting method as claimed in claim 1, wherein, if a dimensional correction is necessary, then the method comprises a step of: d. computing, with the predicting machine, the dimensional modification to be applied based on the resonant frequency characteristic identified in the step b.

    9. The inspecting method as claimed in claim 1, wherein the predicting machine implements a polynomial formula to predict whether or not a dimensional correction is necessary.

    10. The inspecting method as claimed in claim 1, wherein the predicting machine implements a classification carried out for example by a neural network to predict whether or not a dimensional correction is necessary.

    11. The inspecting method as claimed in claim 1, the spiral blank being formed on a wafer comprising a plurality of spiral blanks distributed over several sectors of the wafer, wherein the step b comprises a step consisting in identifying at least one characteristic of a resonant frequency of at least one spiral blank for each sector, and wherein the step c comprises a step consisting in determining for the spiral blanks of each sector a stiffness and whether or not a dimensional correction is necessary.

    12. The inspecting method as claimed in claim 1, comprising a preliminary step consisting in taking into account the material of the spiral or of the spiral blank, and in adjusting a maximum amplitude of the vibrational excitation and/or a range of frequency of the predetermined frequency range as a function of the material of the spiral or of the spiral blank.

    13. The inspecting method as claimed in claim 1, wherein the frequency range extends over a range of frequencies ranging from 0 Hz to 100 kHz, preferably of 0 Hz to 50 kHz, more preferably of 0 Hz to 40 kHz, and most preferably of 10 kHz to 35 kHz.

    14. A method for fabricating a spiral having at least one predetermined resonant frequency comprising the steps consisting in: forming at least one spiral or one spiral blank having dimensions contained within predetermined tolerances necessary to obtain the predetermined resonant frequency, inspecting the spiral or the spiral blank according to the inspecting method of one of the preceding claims.

    15. The fabricating method as claimed in claim 14, comprising a step consisting in: correcting at least one dimension of the spiral blank formed during the step a., according to the computation of step d. of claim 8, in order to obtain a spiral having the predetermined resonant frequency.

    16. The fabricating method as claimed in claim 14, wherein the spiral blank is formed on a wafer, with a plurality of other spiral blanks.

    17. A method for training a predicting machine for implementing the step c of the inspecting method of claim 13, comprising the steps consisting in: iforming spirals or spiral blanks, iiapplying to each of the spirals or to each of the spiral blanks a time-varying vibrational excitation to cover a predetermined frequency range, iiiidentifying at least one characteristic of a resonant frequency of each spiral or of each spiral blank during the application of the predetermined frequency range, ivsetting up a plurality of spirals or of spiral blanks in an oscillating mechanism having a predetermined inertia such as to measure for each spiral or each spiral blank a free oscillation frequency or a stiffness, and/or ivmodeling in a simulation tool a plurality of spirals or of spiral blanks in an oscillating mechanism having a predetermined inertia such as to measure for each spiral or each spiral blank a free oscillation frequency or a stiffness vsupplying to the predicting machine, and for each spiral or each blank: the characteristic of the resonant frequency identified in the step iii; the free oscillation frequency or the stiffness(es) measured in the step ivand/or computed in the step iv.

    18. The training method as claimed in claim 17, wherein the step iiicomprises a preliminary phase of identifying reference measuring points with: the measurement of a displacement of a plurality of predetermined points of the spiral or of the spiral blank, the identification of nodes from among the plurality of predetermined points, which have at least one frequency or resonance peak a displacement amplitude which is zero or less than a first threshold peak value, selecting reference points to be measured during the inspection from among the plurality of predetermined points, which are different from the identified nodes, and which preferably each have a displacement amplitude peak greater than a second threshold peak value.

    Description

    BRIEF DESCRIPTION OF THE DRAWINGS

    [0119] Other details of the invention will appear more clearly on reading the following description, given with reference to the appended drawings wherein:

    [0120] FIG. 1 shows the uncorrected stiffness dispersion curves for the spirals on three different wafers,

    [0121] FIG. 2 shows the centering of the mean stiffness over one wafer around a nominal value,

    [0122] FIGS. 3A-3F are a simplified representation of a method for fabricating a mechanical resonator, here a spiral, on a wafer,

    [0123] FIG. 4 shows a device for evaluating the torque of a spiral,

    [0124] FIG. 5 schematically represents the implementation of the evaluation of the stiffness of a spiral by vibrational analysis,

    [0125] FIG. 6 shows an example of frequencies applied to a silicon wafer bearing spiral blanks, to impose a vibrational excitation,

    [0126] FIG. 7 shows an example of measurement of the amplitudes of displacement of a point of a spiral blank, in response to the imposed frequency range in FIG. 6,

    [0127] FIG. 8 shows in retail a resonance peak identified at a particular frequency in FIG. 7,

    [0128] FIG. 9 shows the measured and superimposed resonance peaks for the particular frequency of FIG. 8,

    [0129] FIG. 10 represents an example of a prediction model constructed based on data extracted from FIG. 9.

    EMBODIMENT OF THE INVENTION

    [0130] FIGS. 3A-3F are a simplified representation of a method for fabricating a mechanical resonator 100 on a wafer 10. The resonator is in particular intended to equip a regulating member of a part for a timepiece and, according to this example, is in the form of a silicon spiral spring 100 which is intended to equip a balance of a mechanical movement for a timepiece.

    [0131] The wafer 10 is illustrated in FIG. 3A as a wafer SOI (silicon on insulator) wafer and comprises a substrate or handler 20 carrying from a sacrificial silicon oxide (SiO.sub.2) layer 30 and a monocrystalline silicone layer 40. By way of example, the substrate 20 can have a thickness of 500 m, the sacrificial layer 30 can have a thickness of 2 m and the silicon layer 40 can have a thickness of 120 m. The monocrystalline silicon layer 40 can have any crystalline orientation.

    [0132] A step of lithography is shown in FIGS. 3B and 3C. The term lithography should be understood to mean all the operations making it possible to transfer an image or pattern onto or above the wafer 10 toward the latter. With reference to FIG. 3B, in this exemplary embodiment, the layer 40 is covered with a protective layer 50, for example made of curable resin. This layer 50 is structured, typically by a step of photolithography using an ultraviolet light source as well as, for example, a photomask (or another type of exposure mask) or a stepper and reticle system. This structuring by lithography forms the patterns for the plurality of resonators in the layer 50, as illustrated in FIG. 3C.

    [0133] Subsequently, in the step of FIG. 3D, the patterns are machined, particularly etched, to form the plurality of resonators 100 in the layer 40. The etching can be carried out by a DRIE (Deep Reactive Ion Etching) technique. After the etching, the remaining part of the protective layer 50 is subsequently eliminated.

    [0134] In FIG. 3E, the resonators are released from the substrate 20 by locally removing the sacrificial layer 30 or by etching all or part of the silicon of the substrate or handler 20. A smoothing (not illustrated) of the etched surfaces can also take place before the release step, for example by a step of thermal oxidization followed by a step of deoxidization, constituted for example of hydrofluoric acid (HF) wet etching.

    [0135] In the last step of the fabricating method in FIG. 3F, the windings 110 of the silicon resonator 100 are covered with a silicon oxide (SiO2) layer 120, typically by a step of thermal oxidization to produce a thermocompensated resonator. The formation of this layer 120, which generally has a thickness of 2-5 m, generally affects the final stiffness of the resonator and must therefore be taken into account during the preceding steps to obtain vibrational characteristics of the spiral leading to the obtainment of a particular natural frequency of the spiral-balance pair in a given watch mechanism. As indicated above, in the stage preceding the production of the thermocompensating layer, the different resonators formed in the wafer generally present a considerable geometric dispersion between them and therefore a considerable dispersion between therefore a considerable dispersion between their stiffnesses, even though the steps of forming the patterns and machining/etching through these patterns are the same for all the resonators.

    [0136] Moreover, this stiffness dispersion is even greater between the spirals of two wafers etched at different times even if the same method specifications are used.

    [0137] The description above relates to silicon resonators 100, but it can be envisioned to make the resonators out of glass, ceramic, carbon nanotubes, or metal.

    [0138] To center the mean of the stiffnesses of the resonators on different wafers in relation to a nominal stiffness value as illustrated in FIG. 2, the resonators obtained in the step 3E on the wafer 10 in question can be deliberately formed with dimensions d that are different to the required dimensions (for example greater) for the obtainment of a nominal or target stiffness. Thus, it is possible to put in place an inspection method intended to estimate the vibrational characteristics of the resonators (natural frequency and/or resonant frequencies) to deduce therefrom the stiffness and/or the actual dimensions of the resonators 100 to correct the dimensions thereof, which will lead to the obtainment of the natural frequency of the desired resonatorbalance pair.

    [0139] This invention makes provision for determining on the basis of at least one characteristic of a resonant frequency of a sample of resonators 100 on the wafer in the step 3E whether or not a geometrical correction of the resonators is necessary. If so, this invention makes provision for accurately computing the thickness of material to be modified (to be removed or added), around each spiral, to obtain the dimensions leading to the obtainment of the vibrational characteristics of the resonators (natural frequency and/or resonant frequencies, and/or stiffness) corresponding to target values, according to a more effective method than the methods of the prior art.

    [0140] Thus, the invention makes provision for determining at least one characteristic of a resonant frequency of a sample of resonators by vibrational measurement and applying a predictive method (for example a computer model or a classification or categorization method) to relate the result of said vibrational measurement to the necessary geometrical correction.

    [0141] One thus makes use of the modal properties of the spiral attached to the wafer. During a training phase, and by an analytic and numerical approach, it is possible to introduce a predicting machine by establishing a predictive model relating the dimensions (particularly thickness) and/or stiffness at certain frequencies (natural frequency or resonant frequencies associated with a resonance peak or with a mid-height width) specifically chosen.

    [0142] Once the training phase is finished (once the modes to be made use of as well as the excitation frequencies have been determined), it is possible to pass on to a predicting phase and to use the predicting machine by making use of the predictive model to inspect the resonators of a produced wafer, in order to predict whether or not a correction of the dimensions is necessary, and where applicable, compute or predict the exact correction to be made to the dimensions of the resonators (by removal if the blank is made with dimensions greater than the final dimensions required, or by addition of material if the blank is made with dimensions less than the final dimensions required, for example).

    [0143] Thus, it is possible to incorporate the inspecting method into a fabricating method to correct, if necessary, the vibrational characteristics of the resonators (natural frequency and/or resonant frequencies, and/or stiffness) to obtain a particular and predetermined oscillation natural frequency, once the resonators are each paired with a balance of a given watch mechanism.

    Vibrational Excitation

    [0144] The measurement of the vibrational excitation of the resonators makes it possible to deduce at least one characteristic of a resonant frequency, such as for example a value of a resonant frequency. To describe in detail: a vibrational excitation must first be imposed on the wafer. Several options are available: [0145] a. Measurements in the frequency domain: [0146] 1Using a piezo-electric source (or any other source making it possible to induce or impose an acoustic excitation) on the edge of the wafer, on, or under the spiral blank 200 to be specifically excited (preferred blank) which excites at a particular frequency f.sub.0 (continuous single-frequency excitation). In this variant, the excitation is sustained. [0147] 2In a variant, the piezo-electric source (or any other source making it possible to induce or impose an acoustic excitation) can also be used on the edge of the wafer, on, or under the spiral blank 200 to be specifically excited (preferred spiral blank) which excites at a time-varying frequency to cover a predetermined frequency range, ranging for example from 0 to 100 kHz, preferably from 0 to 75 kHz, preferably from 0 to 50 kHz, preferably from 5 kHz to 50 kHz, and preferably from 10 to 35 kHz. The entirety of the frequency range can be scanned or covered in a time interval which can range from a fraction of a second to a few seconds. For example, provision can be made for scanning or covering the range of frequencies of the frequency range in less than 0.5 s, less than 1 s, or less than 1.5 s. In this variant, the excitation frequency changes continuously. [0148] b. Time-domain measurements: use an excitation hammer (or any other source making it possible to induce a pulsed acoustic excitation) on the edge of the wafer, on, or under the spiral to be specifically excited (preferred spiral) which gives an acoustic pulse which is as short as possible (pulsed multi-frequency excitation). In this variant, the excitation is a point excitation and not sustained.

    [0149] Moreover, the measurements can be taken following a particular sampling, for example in a sampling range of 4, 2 to 1 Hz. Specifically, the resolution for processing the acquisition data, for example according to a Fourier transform, depends directly on the duration of this acquisition.

    [0150] Moreover, a signal sampling frequency of at least 100 kHz may be chosen if the frequency range extends to 50 kHz for example.

    [0151] In general, provision can finally be made for changing the direction of excitation, i.e. the direction of the movements imposed by the source (vibrations can be imposed along one or more axial directions, and this direction or directions can be varied over time). In the case where a wafer comprising a plurality of resonators is excited, provision can be made for setting the direction of the vibrations such as to point to one or the other of the resonators, as a function of the displacement amplitude measurements described below.

    [0152] Finally, provision can be made for coupling the acoustic source with a divergent cone directed toward the resonators to be excited, and for setting the acoustic source to emit an excitation signal with sufficient amplitude to be detected and measured in an accurate manner by the chosen measurement instruments.

    Measurement of Amplitude or Speed, or Acceleration of Displacement

    [0153] During the excitation, a recording is made, via a suitable measuring means, of the amplitude and the phase (with respect to the excitation source) of oscillation in the 3 directions X, Y (in the plane) and Z (out of the plane) of the specifically excited spiral. Without limitation, the following possible measuring means can be cited: [0154] Optical methods by interferometry: [0155] a. By 3D Doppler effect (laser Doppler vibrometer), [0156] b. Holographic, [0157] Stroboscopic optical methods, [0158] High time-domain resolution chromatic confocal profilometry, [0159] Optical reflectometry: [0160] a. Vibrational analysis by beam deflection on multi-dial detector or camera, [0161] b. Analysis by time-domain analysis of TCSPC type, [0162] Acoustic methods by Doppler ultrasound.

    [0163] FIG. 5 schematically represents a silicon wafer 25 on which are formed a plurality of spiral blanks 200. A vibrational excitation source 400 is coupled with the wafer 25, so as to be able to impose a vibrational excitation. Consequently, each spiral blank 200 will be set to vibrate, and a laser vibrometer 300, here focused on a point of the right spiral blank 200 will be able to measure the amplitudes of vibration of the measurement point over time. Provision can be made for measuring the displacements along a direction normal to the plane of the wafer 25, but it is equally possible to measure the displacements along one or more directions contained in the plane of the wafer 25.

    [0164] Once a particular point has been studied, the laser vibrometer 300 can be displaced to another measuring point of the spiral blank 200, or one can pass on to another spiral blank 200 of the wafer 25. Of course, the spiral blank 200 can alternatively be displaced in relation to the laser vibrometer.

    [0165] FIG. 6 shows an example of vibrational excitation over time. In the given example, the excitation frequency varies over time, between 0 Hz and 50 kHz, and a succession of rising edges can be imposed, each spaced apart by an idle period without excitation. For each measurement point of the spiral blank 200, a plurality of rising edges can be imposed (between 2 rising edges and 60 rising edges), each lasting between 0.5 s and 2 s for example.

    Selection of Reference Points to be Measured

    [0166] As regards the displacement amplitude measurement, during the training phase, provision can be made for a step consisting in identifying points of the resonator for which the vibrational response is significant. Specifically, in the case of a spiral on which a vibration is imposed, especially if the frequency is time-varying, the vibration response will cause the appearance on the spiral of nodes, i.e. particular points of the spiral, the displacement amplitude of which is low or zero. If a measurement of the displacement is taken on a point of the spiral which proves to be a node with one of more particular frequencies, the identification of resonant frequency characteristics will be negatively affected.

    [0167] Thus, it is advantageous to make provision for a preliminary step of measuring displacements at a plurality of predetermined points of the spiral, for example at least ten predetermined points, preferably at least twenty or so predetermined points, and most preferably at least thirty predetermined points. Provision can be made for selecting the predetermined points arranged on an orthonormal frame of reference X-Y in the plane of the spiral.

    [0168] At the end of this preliminary step of measuring amplitudes at the predetermined points, provision can be made for identifying resonant frequencies for each measurement point, and next a step of selecting reference points for which the measurement of displacement amplitude during excitation shows that they are not nodes at these resonant frequencies. In other words, the identified nodes exhibit, at least one resonant frequency, a displacement amplitude which is zero or less than a first threshold peak value, and these points forming nodes are excluded from the reference points to be considered for subsequent measurements. It can also be noted that the reference points differ as a function of the position of the spiral blank 200 on the wafer 25.

    [0169] Typically, it can be considered that at least two reference points will be selected, and preferably at least four reference points will be selected. In the case where the resonator has a radius Ra and is secured or fixed on the wafer by its outer pinning end, four reference points can preferably be chosen and located: [0170] in a first area at less than 0.20Ra (for example on the central ferrule), or [0171] in a second area between 0.05Ra and 0.30Ra (for example on the second winding starting from the ferrule), or [0172] in a third area between 0.35Ra and 0.65Ra (for example on a winding located in the middle of the spiral), or [0173] in a fourth area between 0.65Ra and 0.85Ra (for example on a winding located at three-quarters of the length of the spiral).

    [0174] Thus, the reference points are distant from the part secured to the wafer and naturally have a considerable oscillatory displacement capability, which ensures better accuracy of the displacement measurement.

    [0175] Moreover, one can also measure the displacements of one point of the body of the wafer, and/or one point of the excitation source, to identify or measure, for example, a phase difference or a vibrational attenuation, or else a resonance from a vibrational coupling or else from the wafer. These additional measurements make it possible to ensure that the identified peaks are indeed those of the spiral alone. One can also synchronize the displacement amplitude measurement and the vibrational excitation.

    Determination of the Vibrational Characteristics

    [0176] Several scenarios then arise according to the field previously chosen for the excitation: [0177] a. Measurements in the frequency domain [0178] 1variant with sustained excitation: [0179] i. Integrating over time the amplitude and phase of oscillation over a long enough time to have a good spectral resolution at the excitation frequency f.sub.0, [0180] ii. Offsetting the oscillation frequency by delta f to excite at the frequency f.sub.0+f and repeat the integrating step i, [0181] iii. Reconstructing amplitude and oscillation phase spectra as a function of the excitation frequency (possibly with several peaks at several frequencies). [0182] 2variant with excitation of time-varying frequency: [0183] i. Recording in the time domain the amplitude and the oscillation phase during the frequency scanning of the frequency range, [0184] ii. Repeating the step iat least once, preferably at least three times, [0185] iii. Reconstructing the amplitude and oscillation phase spectra as a function of the excitation frequency (possibly with several frequency peaks). [0186] b. Measurements in the time domain: [0187] i. Recording the displacement over time of the winding along X, Y and Z over a long enough time, in such a way as to obtain a sufficiently representative signal, such as for example a few seconds. [0188] ii. One may choose to record the signal to make it into a reference signal to be compared with other signals measured at other parts. One may also choose to subject the signal to processing of Fourier transform type to identify resonant frequencies in the recorded signal.

    [0189] Consequently, at least one resonance peak can be identified for each excited resonator, and provision is made for determining the resonant frequency not on the basis of the resonance peak apex, i.e. on the maximum amplitude, but rather over an area of the curve located between 25% and 75% of the maximum amplitude value of the resonance peak, for example based on its mid-height width. This is because this processing method, which focuses on a part of the curve between 25% and 75% of the maximum amplitude value of the resonance peak, makes it possible to limit errors due to the singularity of the maximum amplitude point and due to the approximation calculations to reconstruct the apicial part of the resonance peak. The area of the curve located between 25% and 75% of the maximum amplitude value of the resonance peak has a better accuracy than the part above 75% (typically the peak), which offers better accuracy on the exact determined resonant frequency. One may for example take the middle of the segment connecting the two points at mid-height of the resonance peak to determine the resonant frequency associated with the peak in question.

    [0190] FIG. 7 shows an example of a vibrational spectrum for a point of a spiral blank 200 of FIG. 5, reconstructed based on displacement amplitude measurements of the measurement point under consideration in response to the vibrational excitation of FIG. 6, between 10 kHz and 15 kHz. Note the presence of three amplitude peaks, at approximately 11 kHz, 12.3 kHz, and 13.7 kHz. Although this is not shown, between 10 and 30 amplitude peaks can typically be identified if the vibrational excitation scans a frequency range between 0 Hz and 50 kHz. Each amplitude peak has a resonant frequency, and the maximum amplitudes vary greatly.

    [0191] FIG. 8 shows in detail the processing that can be done on an amplitude peak, that at 11 kHz for example. The aim is to find the resonant frequency and give it as accurate a value as possible. Instead of basing this processing on the maximum value of the peak, the applicant has noticed that a better accuracy could be reached by determining the length of the segment connecting the rising part and the falling part of the curve, at mid-height of the peak. The resonant frequency is typically the value in the middle of this segment. However, one may carry out an interpolation on points in the vicinity of the resonance peak to improve the accuracy, and offset the chosen point on the segment, which will not be the middle, in particular if the actual position of the resonance peak is offset, for example due to the chosen sampling frequency.

    [0192] FIG. 9 shows, for the example of an amplitude peak at approximately 10 kHz, the amplitude peaks constructed for ten tested spiral blanks 200. Note that from one spiral blank to the other, the frequency position of the amplitude peak varies (from 9.8 kHz to 10.02 kHz approximately), and the maximum displacement amplitude varies in a ratio of approximately 1 to 5. Since the amplitude peak apices are not truly synthetic, it seems advisable to determine the resonant frequency on the basis of the width of the peak at mid-height.

    [0193] For these tests of FIG. 9, the following resonant frequencies were identified:

    TABLE-US-00001 Spiral no. Resonant frequency (Hz) 2 9824 9 9824 3 9840 8 9840 7 9848 4 9863 10 10020 5 10121 1 10129 6 10148
    Determination of the Stiffness and/or Actual Dimensions of the Bar of Tested Resonators

    [0194] To establish a prediction model that can receive as input the vibrational characteristics (typically a resonant frequency) and give as output a stiffness and/or a dimensional correction, it is necessary, during the training phase, to supply the data relating to the actual stiffness and/or dimensions of the bar of tested resonators. For this purpose, provision can be made for practically measuring a natural frequency of a spiralbalance system in an environment similar to that of a particular watch mechanism.

    [0195] Two alternatives can be implemented. As a first alternative it is possible to couple a predetermined balance directly onto the resonator still attached to the wafer, and measure a natural oscillation frequency of the resonatorbalance pair to compare this natural frequency with an expected natural frequency and above all compute the actual stiffness or the actual dimensions based on the equations 1 to 3 above. As a second alternative, it is possible to finish fabricating the tested resonators, in order to set them up or couple them with a balance individually, here again to measure a natural oscillation frequency of the resonatorbalance pair.

    [0196] In the two alternatives above, it is possible to proceed via an intermediate step of determining the stiffness of each resonator, and then to determine the actual dimensions of the bar of the tested resonators. In other words, it is possible to determine the natural frequency or a resonant frequency and then the stiffness or the dimensions of the bar of the resonator by analyzing the free oscillations of a spiral coupled with a reference balance. In this approach, a laser pointed at the arms of the balance or at the spiral carrier records the times at which the arms of the balance or a safeguard device pass. One then deduces therefrom an estimate of the period, then of the frequency and finally the stiffness. The data collected are essentially scatter plots of the times of passage.

    [0197] Specifically, to assess the stiffness of a spiral on the wafer, several solutions are available, as described by M. Vermot et al, in the Trait de construction horlogre (2011) on pages 178-179. For example, a dynamic evaluation can be made, by coupling the spiral with a reference balance, the inertia of which is known. The measurement of the frequency of the assembly makes it possible to deduce the stiffness of the spiral, accurately. This evaluation can be carried out on the wafer or by detaching the spiral from the wafer. The references and priorities given above provide details of this method.

    [0198] Similarly, the stiffness can also be deduced from a measurement of the reaction torque at the ferrule by means of a rheometer. The acquired signal represents the variation of the torque as a function of amplitude. The analysis of the gradient of this curve for low amplitudes (linear part) makes it possible to deduce the stiffness, and then the dimensions of the bar of the resonator. The dimensions of the bar of the spiral can then be determined.

    [0199] Moreover, provision can be made for estimating by simulation a natural frequency and/or a resonant frequency and/or the stiffness for each tested resonator on the wafer. For this purpose, dimensional measurements of each tested resonator can be taken to reconstruct the resonator by numerical modeling in order to simulate its vibrational response to the imposed spectrum by numerical calculation, and to also find the stiffness of the resonator.

    [0200] An approach by 3D tomography by high-resolution X-ray would make it possible to extract scatter plots giving the 3D material density of the spirals, and, subject to a reconstruction of the adapted images, a mapping of the section of the spiral. These different types of data can be used to deduce the dimensions of the bar and to estimate the stiffness of the spiral by a geometrical approach.

    [0201] Another approach consists in analyzing the driven oscillations of a spiral on a reference balance with an escapement. A laser measurement of the time of passage of the arms of the balance (scatter plots), as described above, can be used to measure the frequency and to deduce the stiffness therefrom. An alternative can be envisioned based on an acoustic acquisition (Witschi-type microphone) which records the shocks of the different operating phases of the escapement/securing system. The measured data are either scatter plots of the times of passage of the arms of the balance, or the variation over time of the acoustic pressure level. These types of experimental data can be used to deduce the period, then the frequency, then the stiffness and finally the dimensions of the bar of the resonator.

    [0202] Returning to the tests discussed above in FIG. 9, a measurement of the stiffness has been taken by coupling each spiral blank 200 with a reference balance, and the stiffnesses below were deduced:

    TABLE-US-00002 Spiral no. Measured stiffness (10.sup.7 N .Math. mm) 2 3.89 9 3.88 3 3.92 8 3.90 7 3.91 4 3.95 10 4.111 5 4.135 1 4.119 6 4.196

    Establishment of the Prediction Model

    [0203] During the training phase, the oscillation amplitude measurements are carried out on physical resonators, and resonant frequencies are identified. In order to be able to subsequently relate the resonant frequencies measured on resonators to stiffnesses and/or corrections of dimension (thickness) to be made, provision must be made for a correlating phase during which a predictive model is constructed.

    [0204] The operations described above (vibrational measurements, identification of the resonance peaks, mid-height bandwidth and its middle or corrected value, determination of the stiffness and/or dimensions of the bar) are used to supply a database which can relate the position of the spiral on the wafer, of the spectra or oscillation periods or mid-height bandwidth and its middle or corrected value with the stiffnesses and/or effective dimensions of the bar of the spiral. As seen above, this database can be constructed from numerical simulations on a finite element model of a spiral. These simulations make it possible to generate reference spectra or oscillation periods associated with the stiffnesses. This database can also be completed by experimental measurements by measuring vibration spectra, oscillation periods and the positions of spirals on the wafer along with their associated stiffnesses. One of the advantages of this approach lies in the fact that the training database is enriched gradually as the tests proceed. This can make it possible to have an adaptive model according to the wafers and the spirals and contributes to the reduction of the standard deviation on wafer stiffness.

    [0205] This database can be used to build a prediction model, and several solutions are available.

    [0206] A numerical model, for example polynomial, can be constructed to compute, as a function of a resonant frequency value, an actual thickness, a dimensional correction or an actual stiffness.

    [0207] A categorization can also be made by making a partitioning into k-means of input data (the results of vibrational measurements, typically the frequency of the resonance peaks) and output data (stiffness, and/or dimensions of the bar of the resonator) and relate them to one another to establish a correspondence.

    [0208] Provision can also be made for processing the images of the resonance peaks by a neural network, for example a perceptron, to make a classification by stiffness or bar dimensions, the classes being able to be defined by increments of values.

    [0209] In summary, the training phase comprises a test phase (excitation of resonators with measurement of the vibrational characteristics to reconstruct a vibration spectrum and identify resonant frequencies). A phase of measuring stiffnesses and/or dimensions of the bar of the resonators is also carried out. Once the input data (resonant frequencies) and output data (stiffnesses and/or dimensions of the bar) for a significant sample are available, the phase of construction of the prediction model can be carried out.

    [0210] To return to the example addressed and described in relation to FIG. 9, the collected data are as follows:

    TABLE-US-00003 Resonant frequency Measured stiffness Spiral no. (Hz) (10.sup.7 N .Math. mm) 2 9824 3.89 9 9824 3.88 3 9840 3.92 8 9840 3.90 7 9848 3.91 4 9863 3.95 10 10020 4.111 5 10121 4.135 1 10129 4.119 6 10148 4.196

    [0211] A linear regression model was made for the data above for the first six rows, and the relationship below was able to be established: [0212] R=0.0015 F10.894, [0213] With [0214] R for the stiffness in 10.sup.7 N.Math.mm [0215] F for the resonant frequency in Hz.

    [0216] The stiffness can therefore be predicted and compared with the actual stiffness measured as shown in the table below, with for the first six rows the data used to build or train the linear regression, and for the last four rows, a prediction only:

    TABLE-US-00004 R R F (10.sup.7 N .Math. mm) (10.sup.7 N .Math. mm) No (Hz) measured predicted deviation 2 9824 3.89 3.84 1.20% 9 9824 3.88 3.84 1.10% 3 9840 3.92 3.87 1.40% 8 9840 3.9 3.87 0.90% 7 9848 3.91 3.88 0.70% 4 9863 3.95 3.9 1.20% Test 10 10020 4.111 4.136 0.60% 5 10121 4.135 4.288 3.70% 1 10129 4.119 4.3 4.40% 6 10148 4.196 4.328 3.20%

    [0217] A maximum error of 4.40% was able to be measured, and the FIG. 10 shows the linear regression curve for the values of the first six rows.

    [0218] Note that it is advantageous to check that the established prediction model has a good sensitivity, i.e. for two different input values, the model gives two distinct output values. The applicant noticed that the sensitivity of the prediction model was not the same for all the resonance peaks. In particular, if one refers to the prediction formula established and shown in FIG. 10, the slope coefficient is 0.0015 10.sup.7 N.Math.mm/Hz. On the one hand, the applicant noticed that the slope coefficient could be greater for high resonant frequencies, which procures a better prediction sensitivity, to predict distinct stiffness values or dimensional correction values, even based on similar resonant frequency values. It is advantageous to make provision, during the training phase, for a step of comparing the sensitivity of the prediction to check/confirm that it is preferable to consider and choose certain resonance peaks at high frequencies (for example above 5 kHz) to then predict as accurately as possible a stiffness and/or a dimensional correction as a function of the measured vibrational response.

    [0219] On the other hand, the applicant also noticed that even for similar resonant frequencies, the resonance modes (particularly modes of deformation and/or displacement of the resonators) could significantly differ, which can also affect the stiffness and/or dimensional correction prediction sensitivity. It is advantageous to make provision, during the training phase, for a step of comparing the sensitivity of the prediction to choose considering subsequently one particular resonant frequency to then predict as accurately as possible a stiffness and/or a dimensional correction as a function of the vibrational response.

    [0220] From the remarks above relating to the study of the prediction sensitivity, provision can be made, during the training phase, for classing the different identified resonance peaks according to the stiffness and/or dimensional correction prediction sensitivity. Provision can then be made for defining the excitation frequency range (which will be applied during a pure prediction phase) to include at least one or more resonant frequency peaks which give the best sensitivity. Thus, imposing a vibrational excitation which varies over the frequency range thus predetermined will guarantee the ability to make an accurate prediction for the identified resonance peak or predictions for each of the identified resonance peaks, which overlap or confirm one another.

    [0221] In general, the training phase makes it possible to choose either resonance peaks at high frequencies and/or resonance peaks which correspond to particular modes of resonance used to predict accurate and reliable values, and the frequency range will be predetermined to include at least one resonance peak and preferably several, to be able to make either a single prediction which is as accurate as possible, or several predictions (one per resonance peak deemed of interest) to then make overlaps, means or recalibrations of the predicted values.

    [0222] Provision can for example be made for predicting several values of stiffness or of dimensional corrections based on several peaks or resonant frequencies, and then computing a final value, by taking, based on the predicted values, a weighted mean by assigning weights to each predicted value, each weight being determined as a function of the sensitivity identified for each peak or corresponding resonant frequency.

    [0223] Alternatively and preferably, provision can be made for having only a single model which takes all the peaks of resonant frequencies as input and which returns the stiffness or the dimensional correction, the training phase of the model service serving precisely to compute the weightings on the input peaks or resonant frequencies.

    Prediction Phase

    [0224] Once the training phase is finished, it is possible to pass on to a prediction phase, for example during a method for inspecting resonators. The inspecting method can typically be carried out on spiral blanks made on a wafer and still attached to this wafer, such as to estimate the stiffness and/or dimensions of the bar of the spirals of the sample, in order to determine whether or not a dimensional correction is to be made.

    [0225] Once the model has been trained, the inspecting procedure to be deployed is as follows: [0226] 1) Marking the position of the spiral on the wafer, vibrational measurement of the spectra or oscillation period (as described above), [0227] 2) Predicting the stiffness and/or dimensions of the bar of the spiral by applying the predictive model, [0228] 3) Determining whether or not a dimensional correction is necessary to achieve the target natural frequency or stiffness.

    [0229] During the inspecting method, it is also possible to quantify the exact correction to be made, such that the fabricating method can include, in addition to the inspection above: [0230] 1) Knowing the effective stiffness of the spiral(s) estimated according to the model and the target stiffness and/or the target dimensions of the bar: applying the necessary corrective measure.

    [0231] Repeating step 1) and step 2) of the inspecting method to control the stiffness/dimensions of the spiral and confirm that the target values are achieved, to the nearest tolerance threshold, or repeating these steps and the dimensional correction until the stiffness/dimension predicted by the model reaches the target values.

    Sampling

    [0232] It is known to fabricate several hundreds of spirals on one wafer and that the dimensions of the bar of the spirals made can vary according to the regions of the wafer. If the evaluation of the stiffness can be done on a single spiral, in practice, it will be carried out on a sample of spirals, distributed across the wafer.

    [0233] Based on the evaluations made, the corrections can be made for any wafer uniformly, or else differentiated by region, if the obtained results vary from one spiral to another. It is thus possible to reduce the standard deviation on the stiffness dispersion. Moreover, if the stiffnesses of all the spirals by application of the model are known, it is possible to determine the optimal correction for reducing the overall dispersion.

    [0234] It can even be envisioned to evaluate all the spirals of the wafer, particularly with a vibrational evaluation, since this is very fast to carry out.

    [0235] Although the examples above have been predetermined mainly on the basis of a fabrication of spirals having initial bar dimensions which are greater than the target bar dimensions, provision can also be made for making spirals having initial bar dimensions which are smaller than the target bar dimensions. The correcting step then consists in adding material, as for example described in the abovementioned document EP3181939.

    [0236] The method, consisting in identifying resonant frequencies by imposing a vibrational excitation on the spiral blanks alone makes it possible to quickly obtain measuring data, without for example having to perform operations of assembly of a balance, while limiting measurement errors since only the spiral blank is tested (there is no error that can be related to the balance, such as its mass, its assembly position etc.)