WO2013037880A1 - Procédé de mesure des propriétés dynamiques d'une structure mécanique - Google Patents

Procédé de mesure des propriétés dynamiques d'une structure mécanique Download PDF

Info

Publication number
WO2013037880A1
WO2013037880A1 PCT/EP2012/067930 EP2012067930W WO2013037880A1 WO 2013037880 A1 WO2013037880 A1 WO 2013037880A1 EP 2012067930 W EP2012067930 W EP 2012067930W WO 2013037880 A1 WO2013037880 A1 WO 2013037880A1
Authority
WO
WIPO (PCT)
Prior art keywords
measurement signal
frequency
transient phase
excitation
profile
Prior art date
Application number
PCT/EP2012/067930
Other languages
English (en)
Inventor
Thibault GOUACHE
Original Assignee
Gouache Thibault
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 Gouache Thibault filed Critical Gouache Thibault
Publication of WO2013037880A1 publication Critical patent/WO2013037880A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0066—Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by exciting or detecting vibration or acceleration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00—Vibration-testing of structures; Shock-testing of structures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H11/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by detecting changes in electric or magnetic properties
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H13/00—Measuring resonant frequency

Definitions

  • the present invention belongs to the field of calculating structures.
  • it relates to a method for measuring the dynamic properties of a mechanical structure, possibly associated with a method for detecting a modification of these properties greater than a given threshold.
  • the invention relates to a method for measuring the dynamic properties of a mechanical structure, the method using an acceleration sensor transmitting acceleration (or speed or displacement) data at regular intervals to a computing device, in the form of a measurement signal transmitted through a network.
  • the method for measuring dynamic properties preferably comprises:
  • a step 320 of identifying the real transient phase profile with a reference transient phase profile extracted from a database the identification being carried out according to a determined distance criterion.
  • a time-frequency map of the measurement signal is produced, then zones of interest are identified in the time-frequency map. which are the zones of maximum amplitude of the measurement signal and a frequency of interest is defined as the frequency of an area of interest.
  • the filtering of the measurement signal around a frequency of interest is carried out so as to retain in the filtered measurement signal only the mode or the resonance around the frequency of interest.
  • this filtering of the signal in the frequency domain is carried out by transforming the fourrier, then by multiplying the Fourier transform of the measurement signal by the Fourier transform of the filter to be applied, and then by the inverse Fourier transform of the product.
  • the method comprises, prior to the step of calculating the real transient phase profile, steps in which:
  • the dominant frequency of the filtered measurement signal is identified, advantageously by the search for a maximum of the Fourier transform of said filtered measurement signal,
  • the local extremums of the filtered measurement signal (output of 130) are identified, this identification comprising a local extremum by half-oscillation of the filtered measurement signal,
  • the phase of the filtered measurement signal is calculated:
  • the method comprises a step in which:
  • a basic model is chosen to carry out the simulations, which is, in a particular implementation, a harmonic oscillator model with a degree of freedom that represents either an oscillator with a degree of freedom or an isolated mode. of a more complex structure.
  • the solution of the equation is simulated numerically at a degree of freedom (or isolated mode) at the excitation as defined by the equations defining the chosen base model.
  • the simulation is done by numerically solving the differential equation of the oscillator with a damped degree of freedom subjected to a transient excitation.
  • the identification of the real transient phase profile to a reference transient phase profile is by a least squares method by exploring the entire database.
  • the identification of the real transient phase profile to a reference transient phase profile is done by a gradient method by generating the reference transient phase profiles on the fly.
  • the method comprises a step in which:
  • the basic model is a harmonic oscillator model with one degree of freedom (that is to say one seeks to identify a vibration mode or a system with a degree of freedom)
  • these parameters are :
  • the method comprises steps in which:
  • FIG. 1 is a flowchart of the steps of an exemplary implementation of the method according to the invention. Detailed description of an embodiment of the invention
  • the method is implemented using an acceleration sensor. Such a sensor is known per se, and is therefore not detailed further here. Alternatively, the method uses a speed or displacement sensor.
  • This acceleration sensor transmits data at regular intervals to a computing device in the form of a measurement signal transmitted through a network known per se.
  • the signal is acquired according to a sampling frequency at least twice that of the frequency of interest (Shannon rule).
  • the measurement signal contains acceleration values. It can also contain values in time. If it does not contain time values, the sampling frequency must be known.
  • This computing device is, in the present example in no way limiting, PC-type microcomputer. It has conventional user interfaces, as well as data storage means.
  • the method for measuring dynamic properties comprises, firstly, in the present example, which is in no way limitative, after a preliminary step 50 of reception by the computing device of a measurement signal emitted by the sensor acceleration, a phase 100 of conditioning the measurement signal.
  • This phase 100 of conditioning the measurement signal comprises several steps:
  • Step 1 In the case where the excitation frequency, called the frequency of interest, around which we want to measure modal parameters of the structure is not known a priori, we realize a time-frequency map of the signal of measured. This operation can be performed in particular by calculating a sliding window Fourier transform, or by wavelet analysis. These two techniques are known per se.
  • modal parameters that one wishes to identify for the structure are for example:
  • Step 120 In a next step, areas of interest are identified in the time-frequency map. These are, in this example, the areas of maximum amplitude of the measurement signal. A frequency of interest is then defined as the frequency of an area of interest. Step 130 After determining a frequency of interest (previously known, or calculated in steps 1 10-120), the measurement signal is filtered around this frequency of interest.
  • this filtering is performed to retain in the measurement signal, only the mode or the resonance around the frequency of interest.
  • This filtering is done numerically. The simplest is to realize it in the frequency domain, that is to say by transform of fourrier, then by multiplication of the Fourier transform of the measurement signal by the Fourier transform of the filter to be applied, then by transform of Fourier inverse of the product. This example is a method among other digital filtering methods.
  • the method of measuring dynamic properties here comprises a phase 200 for calculating the transient phase.
  • This phase 200 for calculating the transient phase is subdivided into four stages.
  • Step 210 the dominant frequency, denoted f, of the filtered measurement signal (output of step 130) is identified, generally by searching for a maximum of the Fourier transform of said filtered measurement signal.
  • Step 220 In parallel, the local extremums of the filtered measurement signal (output of 130) are identified. This identification comprises a local extremum by half-oscillation of the filtered measurement signal. Step 230 Then, based on local extremums identified in step
  • the phase of the filtered measurement signal which is pseudo-harmonic through the filtering step 130, is calculated by the following formula:
  • Step 240 the difference between the phase of the filtered pseudo-harmonic measurement signal and the phase of a purely harmonic dummy signal which has the same frequency as the dominant frequency identified in step 210 is calculated. That is, subtracted from the phase calculated in step 230, 2 * ⁇ * f * time. This difference is called real transient phase profile.
  • a third phase 300 of the method consists of an identification of the real transient phase profile.
  • identification is meant the association of this real transient phase profile with a reference transient phase profile extracted from a previously generated database.
  • Step 310 In a first step of this phase 300, responses of the structure whose modal parameters are to be determined are simulated to form a database of responses.
  • a harmonic oscillator model with a degree of freedom that represents either an oscillator with one degree of freedom is an isolated mode of a more complex structure.
  • This model makes it possible to generate reference curves (for example acceleration / time curves, speed / time, displacement / time or transient reference phase profile) used in the subsequent identification.
  • the template can be written as follows:
  • m is the mass (in the case of a system with one degree of freedom) or the modal mass (in the case of a more complex structure), the modal mass being defined as the partial mass of the structure associated with own mode considered, it is the damping (case with a degree of freedom) or the modal damping, k is the stiffness (case with a degree of freedom) or the modal stiffness,
  • x is the position, x is the speed, x is the acceleration.
  • x is the acceleration.
  • a more complex model may be chosen, in particular with parameters of the model that vary over time.
  • This step is performed independently of the measurement acquisition phase 50 and the measurement signal conditioning phase 100.
  • the solution of the equation is simulated numerically at one degree of freedom (or isolated mode) at the excitation as defined by the three equations above.
  • This simulation is done by a numerical resolution of the differential equation of the oscillator with a damped degree of freedom subjected to a transient excitation. This is for example carried out using a commercial digital calculation software. The number of responses to be simulated depends on the accuracy of the desired identification.
  • Step 320 In the transient phase profile database thus created, a transient profile closest to that identified in step 240 ("real" transient phase profile) is identified according to a determined distance criterion.
  • This identification can be done in different ways: least squares method by exploring the whole database, gradient method by generating transient phase profiles on the fly, genetic search algorithm, learning, etc. These different methods are known to those skilled in the art and are therefore not detailed further here.
  • Step 330 This transient phase profile, identified as being closest to the transient phase profile of the measurement signal, is completely defined by a set of parameters, which are combinations of modal parameters to be identified and parameters of the excitation.
  • the excitation frequency is actually the frequency of interest. It should be noted that the words pulsation and frequency are used interchangeably in the present description, even if there is a difference of 2 ⁇ between the 2.
  • a fourth phase 400 of the method consists of an extension of the measurement. This phase, not always necessary according to the applications, has two stages.
  • Step 420 With the data obtained in step 410, it is possible to simulate a response of the measured structure or mode with arbitrary amplitude excitation.
  • modal mass by adjusting the amplitude of the arbitrary excitation so that the response of the simulated system is of the same amplitude as the signal studied.
  • the adjustment is simply proportional between the output amplitude and the amplitude of the excitation.
  • the aim is to measure the modal parameters obtained in the steps 320, 330, 410 and / or 420. These parameters are then compared to a specification, and according to the conformity or not In that specification, design changes are made to the structure or manufacturing process of the design.
  • the goal is to follow the evolution of the parameters measured in steps 320, 330 and / or 410 (very rarely in step 420). We repeat the process several times and we follow the evolution of the measured parameters.
  • the evolution of the parameters over time is then analyzed according to one or more criteria, so as to detect, for example, a passage of a parameter under a predetermined threshold, or a progressive drift of a slope greater than a predetermined slope. Such events then lead to further verification of the structure or maintenance action.
  • This process allows, from the recordings of displacements, velocities or accelerations of a structure, to obtain a measurement of the properties of the structure (eigenfrequencies, modal damping, modal forms) as well as a measurement of the excitation that it has undergone (frequency of the excitation, amplitude and transitory character).
  • the monitoring of structures (civil engineering, aeronautics, energy, railway) over time by a regular measurement of the properties of the structure thanks to the invention.
  • the method then makes it possible to monitor the health of the structure, its compliance with standards of mechanical strength, and where appropriate to identify a degradation of the structure (or a modification of the properties of the structure).

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

L'invention vise un procédé de mesure des propriétés dynamiques d'une structure mécanique, le procédé utilisant un capteur d'accélération transmettant des données d'accélération à intervalles régulier à un dispositif de calcul, sous forme d'un signal de mesure émis à travers un réseau. Le procédé de mesure de propriétés dynamiques comprend préférentiellement : - une étape 50 de réception par le dispositif de calcul d'un signal de mesure émis par le capteur d'accélération, - une étape 130 de conditionnement du signal de mesure par filtrage du signal de mesure autour d'une fréquence d'intérêt, - une étape 240 de calcul de la différence, nommée profil de phase transitoire réel, entre la phase du signal de mesure filtré, pseudo-harmonique, et la phase d'un signal fictif purement harmonique qui a la même fréquence que la fréquence dominante du signal de mesure filtré, - une étape 320 d'identification du profil de phase transitoire réel à un profil de phase transitoire de référence extrait d'une base de données, l'identification étant réalisée selon un critère de distance déterminé.

Description

Procédé de mesure des propriétés dynamiques d'une structure mécanique
La présente invention appartient au domaine du calcul des structures. Elle vise en particulier un procédé de mesures des propriétés dynamiques d'une structure mécanique, éventuellement associé à un procédé de détection d'une modification de ces propriétés supérieure à un seuil donné.
Préambule et art antérieur
La connaissance des propriétés précises d'une structure mécanique est rendue nécessaire pour des raisons évidentes de sécurité.
Dans le domaine du génie civil, notamment, les structures de type ouvrage d'art (ponts etc.) subissent au cours du temps un vieillissement qui altère leurs propriétés mécaniques. Ces modifications de santé de ces structures peuvent entraîner des dégradations potentiellement dangereuses.
Il est donc souhaitable de mesurer ponctuellement ou de suivre l'évolution de quelques paramètres dynamiques de ces structures.
Exposé de l'invention
L'invention vise un procédé de mesure des propriétés dynamiques d'une structure mécanique, le procédé utilisant un capteur d'accélération transmettant des données d'accélération (ou de vitesse ou de déplacement) à intervalles régulier à un dispositif de calcul, sous forme d'un signal de mesure émis à travers un réseau.
Le procédé de mesure de propriétés dynamiques comprend préférentiellement :
- une étape 50 de réception par le dispositif de calcul d'un signal de mesure émis par le capteur d'accélération,
- une étape 130 de conditionnement du signal de mesure par filtrage du signal de mesure autour d'une fréquence d'intérêt,
- une étape 240 de calcul de la différence, nommée profil de phase transitoire réel, entre la phase du signal de mesure filtré, pseudo-harmonique, et la phase d'un signal fictif purement harmonique qui a la même fréquence que la fréquence dominante du signal de mesure filtré,
- une étape 320 d'identification du profil de phase transitoire réel à un profil de phase transitoire de référence extrait d'une base de données, l'identification étant réalisée selon un critère de distance déterminé. Dans un mode de mise en œuvre particulier du procédé, préalablement à l'étape 130 de conditionnement du signal de mesure, on réalise une carte temps- fréquence du signal de mesure, puis on identifie dans la carte temps-fréquence des zones d'intérêt qui sont, les zones d'amplitude maximale du signal de mesure et on définit une fréquence d'intérêt comme la fréquence d'une zone d'intérêt.
Dans une mise en œuvre préférée, le filtrage du signal de mesure autour d'une fréquence d'intérêt est réalisé pour ne retenir dans le signal de mesure filtré, que le mode ou la résonnance autour de la fréquence d'intérêt.
Avantageusement, on réalise ce filtrage du signal dans le domaine fréquentiel, par transformée de fourrier, puis par multiplication de la transformée de fourrier du signal de mesure par la transformée de Fourier du filtre à appliquer, puis par transformée de Fourier inverse du produit.
Dans une mise en œuvre préférée, le procédé comprend, préalablement à l'étape de calcul du profil de phase transitoire réel, des étapes dans lesquelles :
- 210 on identifie la fréquence dominante du signal de mesure filtré, avantageusement par la recherche d'un maximum de la transformée de Fourrier dudit signal de mesure filtré,
- 220 on identifie les extremums locaux du signal de mesure filtré (sortie de 130), cette identification comprenant un extremum local par demi-oscillation du signal de mesure filtré,
- 230 sur la base des extremums locaux identifiés dans l'étape 220, on calcule la phase du signal de mesure filtré :
arcsin (valeur du signal de mesure filtré / extremum local), si la dérivée locale (pente du signal) du signal de mesure filtré est positive, et
π - arcsin (valeur du signal de mesure filtré / extremum local), si la dérivée locale du signal de mesure filtré est négative.
Dans une mise en œuvre préférée, le procédé comprend une étape dans laquelle :
- 310 on simule des réponses de la structure dont on cherche à déterminer des paramètres modaux, pour constituer une base de données de réponses, puis on calcule pour chaque réponse simulée un profil de phase transitoire dit " de référence". Préférentiellement, on choisit, pour réaliser les simulations, un modèle de base, lequel est, dans une mise en œuvre particulière, un modèle d'oscillateur harmonique à un degré de liberté qui représente soit un oscillateur à un degré de liberté soit un mode isolé d'une structure plus complexe.
Afin de réaliser les profils de phase transitoire de référence, on simule numériquement la solution de l'équation à un degré de liberté (ou du mode isolé) à l'excitation comme défini par les équations définissant le modèle de base choisi.
La simulation se fait en résolvant numériquement l'équation différentielle de l'oscillateur à un degré de liberté amorti soumis à une excitation transitoire.
Dans une mise en œuvre préférée, l'identification du profil de phase transitoire réel à un profil de phase transitoire de référence se fait par une méthode de moindres carrés en explorant toute la base de données.
Dans une mise en œuvre alternative, l'identification du profil de phase transitoire réel à un profil de phase transitoire de référence se fait par une méthode par gradient en générant les profils de phase transitoire de référence à la volée. Dans une mise en œuvre préférée, le procédé comprend une étape dans laquelle :
- 330 on extrait du profil de phase transitoire de référence le plus proche du profil de phase transitoire réel, des paramètres définissant complètement ledit profil de phase transitoire de référence, ces paramètres de définition du profil, combinaisons de paramètres modaux à identifier et de paramètres de l'excitation initiale, étant ensuite considérés comme ceux du profil de phase transitoire réel.
Dans le cas où le modèle de base est un modèle d'oscillateur harmonique à un degré de liberté (c'est-à-dire qu'on cherche à identifier un mode de vibration ou un système à un degré de liberté), ces paramètres sont :
ω
- rapport — de la fréquence d'excitation ω de la structure sur la fréquence ωθ
propre ωθ de la structure,
- rapport de l'amortissement modal c sur le produit de la masse modale m.coO
m par la fréquence propre ωθ, et - rapport de l'amortissement modal c sur le produit de la masse modale m m.a
par la constante de temps transitoire a.
Dans une mise en œuvre avantageuse, le procédé comprend des étapes dans lesquelles :
- 410 sur la base des paramètres de définition du profil (identifiés dans l'étape 330), et la valeur de la fréquence d'excitation de la structure (fréquence dominante du signal de mesure filtré, identifiée dans l'étape 210), on calcule la fréquence propre ωθ, le rapport (amortissement modal / masse modale) = — et la constante de temps m
transitoire a,
- 420 on simule une réponse de la structure ou du mode mesuré avec une excitation d'amplitude arbitraire, puis on calcule un rapport amplitude de l'excitation réelle / masse modale en ajustant l'amplitude de l'excitation arbitraire pour que la réponse du système simulé soit de même amplitude que le signal étudié, l'ajustement étant proportionnel entre l'amplitude de sortie et l'amplitude de l'excitation.
Présentation des figures
Les caractéristiques et avantages de l'invention seront mieux appréciés grâce à la description qui suit, description qui expose les caractéristiques de l'invention au travers d'un exemple non limitatif d'application.
La description s'appuie sur les figures annexées dans lesquelles :
La figure 1 est un organigramme des étapes d'un exemple de mise en œuvre du procédé selon l'invention. Description détaillée d'un mode de réalisation de l'invention
La suite de la description est réalisée dans le cas particulier de la surveillance d'une structure de type pont.
Il est cependant clair que le procédé peut s'appliquer à tous types de structures.
Le procédé est mis en œuvre en utilisant un capteur d'accélération. Un tel capteur est connu en soi, et n'est donc pas détaillé plus avant ici. En variante, le procédé utilise un capteur de vitesse ou de déplacement.
Ce capteur d'accélération transmet des données à intervalles régulier, à un dispositif de calcul, sous forme d'un signal de mesure émis à travers un réseau connu en soi. Conformément à une règle classique, le signal est acquis selon une fréquence d'échantillonnage égale au moins au double de celle de la fréquence d'intérêt (règle de Shannon).
Le signal de mesure contient des valeurs d'accélérations. Il peut aussi contenir des valeurs en temps. S'il ne contient pas de valeurs en temps, il faut que la fréquence d'échantillonnage soit connue.
Ce dispositif de calcul est, dans le présent exemple nullement limitatif, de type micro-ordinateur PC. Il dispose d'interfaces utilisateur classiques, ainsi que de moyens de stockage de données. Comme on le voit sur la figure 1 , le procédé de mesure de propriétés dynamiques comprend en premier lieu, dans le présent exemple nullement limitatif, après une étape préliminaire 50 de réception par le dispositif de calcul d'un signal de mesure émis par le capteur d'accélération, une phase 100 de conditionnement du signal de mesure.
Cette phase 100 de conditionnement du signal de mesure comporte plusieurs étapes :
Etape 1 10 Dans le cas où la fréquence d'excitation, dite fréquence d'intérêt, autour de la quelle on veut mesurer des paramètres modaux de la structure n'est pas connue a priori, on réalise une carte temps-fréquence du signal de mesure. Cette opération peut être réalisée notamment par calcul d'une transformée de Fourrier à fenêtre glissante, ou par une analyse par ondelettes. Ces deux techniques sont connues en soi.
Les paramètres modaux que l'on souhaite identifier pour la structure sont par exemple :
- la fréquence de résonnance.
- l'amortissement modal.
- la masse modale.
- la raideur modale.
- la déformée modale.
Etape 120 Dans une étape suivante, on identifie dans la carte temps- fréquence des zones d'intérêt. Celles-ci sont, dans le présent exemple, les zones d'amplitude maximale du signal de mesure. On définit alors une fréquence d'intérêt comme la fréquence d'une zone d'intérêt. Etape 130 Après la détermination d'une fréquence d'intérêt (préalablement connue, ou calculée dans les étapes 1 10-120), on filtre le signal de mesure autour de cette fréquence d'intérêt.
Physiquement, ce filtrage est réalisé pour ne retenir dans le signal de mesure, que le mode ou la résonnance autour de la fréquence d'intérêt. Ce filtrage est réalisé numériquement. Le plus simple est de le réaliser dans le domaine fréquentiel, c'est-à- dire par transformée de fourrier, puis par multiplication de la transformée de fourrier du signal de mesure par la transformée de Fourier du filtre à appliquer, puis par transformée de Fourier inverse du produit. Cet exemple est une méthode parmi d'autres méthodes numériques de filtrage.
Après cette phase 100 de conditionnement du signal de mesure, le procédé de mesure de propriétés dynamiques comprend ici une phase 200 de calcul de la phase transitoire. Cette phase 200 de calcul de la phase transitoire se subdivise en quatre étapes.
Etape 210 Dans un premier temps, on identifie la fréquence dominante, notée f, du signal de mesure filtré (sortie de l'étape 130), en général par la recherche d'un maximum de la transformée de Fourrier dudit signal de mesure filtré.
Etape 220 En parallèle, on identifie les extremums locaux du signal de mesure filtré (sortie de 130). Cette identification comprend un extremum local par demi-oscillation du signal de mesure filtré. Etape 230 Puis, sur la base des extremums locaux identifiés dans l'étape
220, on calcule la phase du signal de mesure filtré, qui est pseudo-harmonique grâce à l'étape de filtrage 130, par la formule suivante :
arcsin (valeur du signal de mesure filtré / extremum local), si la dérivée locale (pente du signal) du signal de mesure filtré est positive, et
77 - arcsin (valeur du signal de mesure filtré / extremum local), si la dérivée locale du signal de mesure filtré est négative.
Il est à noter que des formules similaires en arccosinus sont également utilisables dans le procédé de la présente invention.
Etape 240 Enfin, on calcule la différence entre la phase du signal de mesure filtré, pseudo-harmonique, et la phase d'un signal fictif purement harmonique qui a la même fréquence que la fréquence dominante identifiée dans l'étape 210. C'est-à-dire on soustrait, à la phase calculée dans l'étape 230, 2*π *f*temps. Cette différence est nommée profil de phase transitoire réel.
Une troisième phase 300 du procédé consiste en une identification du profil de phase transitoire réel. On entend par identification l'association de ce profil de phase transitoire réel à un profil de phase transitoire de référence extrait d'une base de données préalablement générée.
Etape 310 Dans une première étape de cette phase 300, on simule des réponses de la structure dont on cherche à déterminer des paramètres modaux, pour constituer une base de données de réponses.
On choisit, pour réaliser les simulations ou pour définir les paramètres à identifier, un modèle de base, qui est dans l'exemple de mise en œuvre décrit ici mais non limitatif, un modèle d'oscillateur harmonique à un degré de liberté qui représente soit un oscillateur à un degré de liberté soit un mode isolé d'une structure plus complexe. Ce modèle permet de générer des courbes de références (par exemple courbes accélération / temps, vitesse / temps, déplacement / temps ou profil de phase transitoire de référence) utilisées dans la suite de l'identification.
Le modèle peut-être écrit comme suit :
mx' + ex + kx— input(t)
inp'ut(t)— F · a(t) · cos( t)
Figure imgf000009_0001
dans lequel m est la masse (dans le cas d'un système à un degré de liberté) ou la masse modale (dans le cas d'une structure plus complexe), la masse modale étant définie comme la masse partielle de la structure associée au mode propre considéré, c est l'amortissement (cas à un degré de liberté) ou l'amortissement modal, k est la raideur (cas à un degré de liberté) ou la raideur modale,
f représente le temps,
x est la position, x est la vitesse, x est l'accélération. En ce qui concerne l'excitation :
F est l'amplitude de l'excitation,
le choix d'un cosinus ou d'un sinus avec ou sans décalage cos/sin(wt + décalage) ne change rien,
a constante de temps transitoire. La forme de a(t) est également sans influence sur le résultat. L'essentiel est que a module la pente de l'amplitude de l'excitation. On pourrait également choisir ici une rampe dont la pente est égale à a puis un plateau, ωο est égale à la racine carré de k/m,
la relation entre la fréquence et ω s'exprime par ω = 2πΐ.
On peut donc différencier les paramètres de la structure (un degré de liberté ou mode) et les paramètres de l'excitation. En fonction des applications les paramètres de l'excitation sont ou ne sont pas d'intérêt.
Il est à noter qu'un des points forts du présent procédé est d'être capable de remonter à l'excitation.
Pour des applications particulières, un modèle plus complexe peut être choisi, avec notamment des paramètres du modèle variables dans le temps.
On applique à ces réponses simulées, formatées à cet effet comme des signaux de mesure, les étapes 210 à 240. On génère ainsi une base de données de profils de phase transitoire de référence.
Cette étape est réalisée indépendamment de la phase 50 d'acquisition des mesures et à la phase 100 de conditionnement du signal de mesure. Afin de réaliser les profils de phase transitoire de référence, on simule numériquement la solution de l'équation à un degré de liberté (ou du mode isolé) à l'excitation comme défini par les trois équations ci-dessus.
Cette simulation se fait par une résolution numérique de l'équation différentielle de l'oscillateur à un degré de liberté amorti soumis à une excitation transitoire. Celle-ci est par exemple réalisée en utilisant un logiciel de calcul numérique du commerce. Le nombre de réponse à simuler dépend de la précision de l'identification voulue.
Par exemple si on veut identifier ω/ωο à une valeur A près, il faut simuler un cas tous les ω/ωο=Α. La base de données de réponses n'a évidement pas besoin d'être générée à chaque fois mais peut-être conservée en mémoire. Il est possible de complexifier le modèle tout en conservant exactement le même processus. L'exemple le plus simple est de faire varier les paramètres de l'excitation ou de l'oscillateur à un degré de liberté / mode dans le temps.
Il est aussi possible d'entreprendre la même démarche en prenant un modèle à deux degrés de liberté couplés (ou N, avec N entier relatif) et d'identifier les réponses des deux degrés de liberté en même temps (à partir de deux ou N signaux).
Etape 320 On identifie dans la base de données de profils de phase transitoire ainsi créée, un profil transitoire le plus proche de celui identifié dans l'étape 240 (profil de phase transitoire "réel") selon un critère de distance déterminé.
Cette identification peut se faire de différentes manières : méthode de moindres carrés en explorant toute la base de données, méthode par gradient en générant les profils de phase transitoire à la volée, algorithme de recherche génétique, apprentissage, etc. Ces différentes méthodes sont connues de l'homme du métier et ne sont donc pas détaillées plus avant ici.
Etape 330 Ce profil de phase transitoire, identifié comme étant le plus proche du profil de phase transitoire du signal de mesure, est complètement défini par un ensemble de paramètres, qui sont des combinaisons de paramètres modaux à identifier et des paramètres de l'excitation.
Ces paramètres sont, dans le cas où on cherche à identifier un mode de vibration ou un système à un degré de liberté :
ω
- rapport — de la fréquence d'excitation ω de la structure sur la fréquence ωθ
propre ωθ de la structure. La fréquence d'excitation est en réalité la fréquence d'intérêt. Il faut noter que les mots pulsations et fréquence sont utilisés indifféremment dans la présente description, même s'il y a une différence de 2π entre les 2.
- rapport de l'amortissement modal c sur le produit de la masse modale m.coO
m par la fréquence propre ωθ, et
- rapport de l'amortissement modal c sur le produit de la masse modale m m.a
par la constante de temps transitoire a.
Ces paramètres sont alors considérés comme la mesure des mêmes paramètres pour le profil de phase transitoire du signal de mesure. Une quatrième phase 400 du procédé consiste en une extension de la mesure. Cette phase, non toujours nécessaire selon les applications, comporte deux étapes.
Dans le cas du « suivi dans le temps des paramètres / application génie civile», on peut faire le suivi dans le temps sur les trois rapports identifiés dans l'étape 330 ou sur les paramètres identifiés dans l'étape 400. Dans le cas de vérification de « conformité à un cahier des charges / application spatiale», les paramètres à identifier ou d'intérêt dépendent de l'application et du client (en fonction des paramètres dans le cahier des charges ou d'intérêt pour le client). Etape 410 Avec les valeurs des trois paramètres identifiés dans l'étape 330 et la valeur de la fréquence d'excitation (de façon équivalente la fréquence dominante f : ω ou f, différents de 2ττ) de la structure, identifiée dans l'étape 210, on calcule la c fréquence propre ωθ, le rapport (amortissement modal / masse modale) = — et la m constante de temps transitoire a.
Etape 420 Avec les données obtenues dans l'étape 410, il est possible de simuler une réponse de la structure ou du mode mesuré avec une excitation d'amplitude arbitraire.
F
On obtient alors un rapport — c'est à dire (amplitude de l'excitation réelle / m
masse modale) en ajustant l'amplitude de l'excitation arbitraire pour que la réponse du système simulé soit de même amplitude que le signal étudié. L'ajustement est simplement proportionnel entre l'amplitude de sortie et l'amplitude de l'excitation.
Dans le cas d'une application au domaine spatial, le but est de mesurer les paramètres modaux obtenus dans les étapes 320, 330, 410 et/ou 420. Ces paramètres sont alors comparés à un cahier des charges, et selon la conformité ou non audit cahier des charges, des modifications de conception sont apportées à la structure ou au procédé de fabrication de celle-ci. Dans le cas d'une application à la surveillance d'infrastructure, le but est de suivre l'évolution des paramètres mesurés dans les étapes 320, 330 et/ou 410 (très rarement dans l'étape 420). On répète donc le processus plusieurs fois et on suit l'évolution des paramètres mesurés. L'évolution des paramètres dans le temps est alors analysée selon un ou plusieurs critères, de manière à détecter, par exemple, un passage d'un paramètre sous un seuil prédéterminé, ou une dérive progressive d'une pente supérieure à une pente prédéterminée. De tels événements entraînent alors une action de vérification plus approfondie de la structure ou une action de maintenance.
Avantages de l'invention
Ce procédé, permet, à partir des enregistrements des déplacements, des vitesses ou des accélérations d'une structure, d'obtenir une mesure des propriétés de la structure (fréquences propres, amortissements modaux, formes modales) ainsi qu'une mesure de l'excitation qu'elle a subit (fréquence de l'excitation, amplitude et caractère transitoire).
Deux applications principales sont envisagées:
- la surveillance de structures (génie civile, aéronautique, énergétique, ferroviaire) dans le temps par une mesure régulière des propriétés de la structure grâce à l'invention. Le procédé permet alors de surveiller la santé de la structure, sa conformité à des normes de tenue mécanique, et le cas échéant d'identifier une dégradation de la structure (ou une modification des propriétés de la structure).
- la caractérisation des propriétés d'une structure et de l'excitation qu'elle subit à un instant donné (si cette excitation est transitoire). Cette identification permet de vérifier la conformité à un cahier des charges. Elle permet également de mettre à jour les propriétés modales de la structure dans les modèles utilisés comme référence ou pour simuler la réponse de la structure à une charge exceptionnelle ou à une charge de fonctionnement.

Claims

REVENDICATIONS
1 . Procédé de mesure des propriétés dynamiques d'une structure mécanique, le procédé utilisant un capteur d'accélération transmettant des données d'accélération à intervalles régulier à un dispositif de calcul, sous forme d'un signal de mesure émis à travers un réseau,
caractérisé en ce que le procédé de mesure de propriétés dynamiques comprend:
- une étape 50 de réception par le dispositif de calcul d'un signal de mesure émis par le capteur d'accélération,
- une étape 130 de conditionnement du signal de mesure par filtrage du signal de mesure autour d'une fréquence d'intérêt,
- une étape 240 de calcul de la différence, nommée profil de phase transitoire réel, entre la phase du signal de mesure filtré, pseudo-harmonique, et la phase d'un signal fictif purement harmonique qui a la même fréquence que la fréquence dominante du signal de mesure filtré,
- une étape 320 d'identification du profil de phase transitoire réel à un profil de phase transitoire de référence extrait d'une base de données, l'identification étant réalisée selon un critère de distance déterminé.
2. Procédé selon la revendication 1 , caractérisé en ce que, préalablement à l'étape 130 de conditionnement du signal de mesure, on réalise une carte temps- fréquence du signal de mesure, puis on identifie dans la carte temps-fréquence des zones d'intérêt qui sont, les zones d'amplitude maximale du signal de mesure et on définit une fréquence d'intérêt comme la fréquence d'une zone d'intérêt.
3. Procédé selon la revendication 2, caractérisé en ce que le filtrage du signal de mesure autour d'une fréquence d'intérêt est réalisé pour ne retenir dans le signal de mesure filtré, que le mode ou la résonnance autour de la fréquence d'intérêt.
4. Procédé selon la revendication 3, caractérisé en ce qu'on réalise ce filtrage du signal dans le domaine fréquentiel, par transformée de fourrier, puis par multiplication de la transformée de fourrier du signal de mesure par la transformée de Fourier du filtre à appliquer, puis par transformée de Fourier inverse du produit.
5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce qu'il comprend, préalablement à l'étape de calcul du profil de phase transitoire réel, des étapes dans lesquelles :
- 210 on identifie la fréquence dominante du signal de mesure filtré, avantageusement par la recherche d'un maximum de la transformée de Fourrier dudit signal de mesure filtré,
- 220 on identifie les extremums locaux du signal de mesure filtré (sortie de 130), cette identification comprenant un extremum local par demi-oscillation du signal de mesure filtré,
- 230 sur la base des extremums locaux identifiés dans l'étape 220, on calcule la phase du signal de mesure filtré :
arcsin (valeur du signal de mesure filtré / extremum local), si la dérivée locale (pente du signal) du signal de mesure filtré est positive, et
π - arcsin (valeur du signal de mesure filtré / extremum local), si la dérivée locale du signal de mesure filtré est négative.
6. Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il comprend, une étape dans laquelle :
- 310 on simule des réponses de la structure dont on cherche à déterminer des paramètres modaux, pour constituer une base de données de réponses, puis on calcule pour chaque réponse simulée un profil de phase transitoire dit " de référence".
7. Procédé selon la revendication 6, caractérisé en ce qu'on choisit, pour réaliser les simulations, un modèle de base, lequel est un modèle d'oscillateur harmonique à un degré de liberté qui représente soit un oscillateur à un degré de liberté soit un mode isolé d'une structure plus complexe.
8. Procédé selon la revendication 7, caractérisé en ce que, afin de réaliser les profils de phase transitoire de référence, on simule numériquement la solution de l'équation à un degré de liberté (ou du mode isolé) à l'excitation comme défini par les équations définissant le modèle de base choisi, la simulation se faisant en résolvant numériquement l'équation différentielle de l'oscillateur à un degré de liberté amorti soumis à une excitation transitoire.
9. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que l'identification du profil de phase transitoire réel à un profil de phase transitoire de référence se fait par une méthode de moindres carrés en explorant toute la base de données.
10. Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que l'identification du profil de phase transitoire réel à un profil de phase transitoire de référence se fait par une méthode par gradient en générant les profils de phase transitoire de référence à la volée.
1 1. Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce qu'il comprend une étape dans laquelle :
- 330 on extrait du profil de phase transitoire de référence le plus proche du profil de phase transitoire réel, des paramètres définissant complètement ledit profil de phase transitoire de référence, ces paramètres de définition du profil, combinaisons de paramètres modaux à identifier et de paramètres de l'excitation initiale, étant ensuite considérés comme ceux du profil de phase transitoire réel.
12. Procédé selon la revendication 1 1 , caractérisé en ce que les paramètres de définition du profil sont, dans le cas où on cherche à identifier un mode de vibration ou un système à un degré de liberté :
ω
- rapport — de la fréquence d'excitation ω de la structure sur la fréquence ωθ
propre ωθ de la structure,
- rapport de l'amortissement modal c sur le produit de la masse modale m.coO
m par la fréquence propre ωθ, et
- rapport de l'amortissement modal c sur le produit de la masse modale m m.a
par la constante de temps transitoire a.
13. Procédé selon l'une quelconque des revendications 1 1 à 12, caractérisé en ce qu'il comprend des étapes dans lesquelles :
- 410 sur la base des paramètres de définition du profil (identifiés dans l'étape 330), et la valeur de la fréquence d'excitation de la structure (fréquence dominante du signal de mesure filtré, identifiée dans l'étape 210), on calcule la fréquence propre ωθ, le rapport (amortissement modal / masse modale) = — et la constante de temps m
transitoire a,
- 420 on simule une réponse de la structure ou du mode mesuré avec une excitation d'amplitude arbitraire, puis on calcule un rapport amplitude de l'excitation réelle / masse modale en ajustant l'amplitude de l'excitation arbitraire pour que la réponse du système simulé soit de même amplitude que le signal étudié, l'ajustement étant proportionnel entre l'amplitude de sortie et l'amplitude de l'excitation.
PCT/EP2012/067930 2011-09-13 2012-09-13 Procédé de mesure des propriétés dynamiques d'une structure mécanique WO2013037880A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1158151 2011-09-13
FR1158151A FR2980010B1 (fr) 2011-09-13 2011-09-13 Procede de mesure des proprietes dynamiques d'une structure mecanique

Publications (1)

Publication Number Publication Date
WO2013037880A1 true WO2013037880A1 (fr) 2013-03-21

Family

ID=46881048

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2012/067930 WO2013037880A1 (fr) 2011-09-13 2012-09-13 Procédé de mesure des propriétés dynamiques d'une structure mécanique

Country Status (2)

Country Link
FR (1) FR2980010B1 (fr)
WO (1) WO2013037880A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11660049B2 (en) 2015-12-22 2023-05-30 The University Of Sheffield Apparatus and methods for determining force applied to the tip of a probe

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4164149A (en) * 1978-05-09 1979-08-14 Shigeo Okubo Method and system for monitoring the angular deformation of structural elements
WO1991019173A1 (fr) * 1990-06-01 1991-12-12 Technology Integration And Development Group Incorporated Procede pour evaluer l'integrite structurale de structures composites
US20080011091A1 (en) * 2006-06-27 2008-01-17 Abnaki Systems, Inc. Method for measuring loading and temperature in structures and materials by measuring changes in natural frequencies
US20110036180A1 (en) * 2009-08-17 2011-02-17 Fdh Engineering, Inc. Method for Determining Tension in a Rod

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4164149A (en) * 1978-05-09 1979-08-14 Shigeo Okubo Method and system for monitoring the angular deformation of structural elements
WO1991019173A1 (fr) * 1990-06-01 1991-12-12 Technology Integration And Development Group Incorporated Procede pour evaluer l'integrite structurale de structures composites
US20080011091A1 (en) * 2006-06-27 2008-01-17 Abnaki Systems, Inc. Method for measuring loading and temperature in structures and materials by measuring changes in natural frequencies
US20110036180A1 (en) * 2009-08-17 2011-02-17 Fdh Engineering, Inc. Method for Determining Tension in a Rod

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11660049B2 (en) 2015-12-22 2023-05-30 The University Of Sheffield Apparatus and methods for determining force applied to the tip of a probe

Also Published As

Publication number Publication date
FR2980010A1 (fr) 2013-03-15
FR2980010B1 (fr) 2013-11-22

Similar Documents

Publication Publication Date Title
Zhang et al. Allan variance analysis on error characters of MEMS inertial sensors for an FPGA-based GPS/INS system
Brincker Some elements of operational modal analysis
CN109298457B (zh) 适用于原子干涉重力仪的振动噪声修正补偿方法
KR100994247B1 (ko) 간섭 신호의 고조파 성분분석을 이용한 가속도계의 위상 감도 평가 방법 및 장치
Filipot et al. A unified deep‐to‐shallow water wave‐breaking probability parameterization
FR2976353A1 (fr) Procede d'estimation simplifie de l'orientation d'un objet et centrale d'attitude mettant en oeuvre un tel procede
CN107764257A (zh) 一种惯性器件数值模拟方法
CN108225374A (zh) 一种融合遗传算法的Allan方差分析法
FR2832801A1 (fr) Procede pour evaluer la frequence instantanee d'une excitation mecanique exercee sur une roue d'un vehicule automobile, et applications
Hadianfard et al. Analysis of modal frequencies estimated from frequency domain decomposition method
FR2502341A1 (fr) Appareil pour mesurer les deplacements de differents points de structures excitees par des forces artificielles ou naturelles
FR3069316B1 (fr) Procede d'estimation du mouvement d'un objet evoluant dans un champ magnetique
WO2013037880A1 (fr) Procédé de mesure des propriétés dynamiques d'une structure mécanique
Damadipour et al. Seismic modal identification using a new approach based on weighted transmissibility
Kholikov et al. Meridional-flow measurements from global oscillation network group data
EP2210059B1 (fr) Procede d'integration d'un signal vibratoire a partir de differentielles du signal et procede correspondant de correction d'un signal acquis par un capteur d'images
CN103398799A (zh) 法布里珀罗干涉圆环图像处理方法
CA3182712A1 (fr) Procede de determination des composantes de la vitesse du vent au moyen d'un capteur de teledetection par laser et au moyen d'une coherence temporelle
WO2020038985A1 (fr) Procédé de détection d'un fluide et système associé
Kendall et al. Comparison of point and array-computed rotations for the TAIGER explosions of 4 March 2008
Kijewski et al. Identification of dynamic properties of a tall building from full-scale response measurements
Pham et al. Characterizing Stochastic Errors of MEMS–Based Inertial Sensors
CN103543294B (zh) 一种基于附加质量的微米光栅加速度计测试方法
CA3028856A1 (fr) Procede de caracterisation du sous-sol d'une region utilisant des signaux sismiques passifs, et systeme correspondant
EP2488838B1 (fr) Procédé et dispositif d'estimation de la pression acoustique dans un fluide en mouvement et programme d'ordinateur

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: 12761720

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

WPC Withdrawal of priority claims after completion of the technical preparations for international publication

Ref document number: 1158151

Country of ref document: FR

Date of ref document: 20140311

Free format text: WITHDRAWN AFTER TECHNICAL PREPARATION FINISHED

122 Ep: pct application non-entry in european phase

Ref document number: 12761720

Country of ref document: EP

Kind code of ref document: A1