Brevet FR2938648A1 — publié le 20100521

Procédé d'évaluation du vieillissement d'un ensemble électronique L'invention concerne un procédé pour évaluer l'état de vieillissement appelé dommage par fatigue en vibration dans un système électronique quelconque. Il peut s'agir de cartes électroniques individuelles ou d'un ensemble de cartes électroniques. L'évaluation repose sur un dispositif de mesure et de traitement numérique de grandeurs physiques relatives à la sollicitation qui est à l'origine du dommage. Des dispositifs similaires existent actuellement pour des pièces et structures mécaniques, mais ne disposent pas de certaines caractéristiques nécessaires pour être appliqués aux cartes électroniques. Ces caractéristiques sont la capacité à traiter en temps réel, et de façon autonome avec toute la capacité de calcul embarquée, l'estimation du dommage sans avoir recours ni à de la télémesure, ni à du post-traitement. La télémesure et le post-traitement de données sont couramment utilisés pour surveiller des aérostructures et des ouvrages d'art, mais les dispositifs correspondants ne répondent pas simultanément aux critères d'embarquabilité et d'autonomie de toute la puissance de calcul. L'invention repose sur l'association d'un ou plusieurs capteurs cinématiques de la vibration placés en des points particuliers, d'un observateur d'état pour évaluer dynamiquement des contraintes mécaniques localisées en des points critiques et d'un calcul de dommage en ces points critiques, intégrables tous les trois dans un dispositif embarqué et autonome quant à sa fonction. Les points particuliers de mesure et les points critiques de calcul de dommage sont distincts ou confondus selon les applications. La notion d'observateur d'état trouve sa définition dans les disciplines de l'automatique et de la théorie de l'information. Elle s'applique à des systèmes dynamiques représentés par des grandeurs d'état. Lorsqu'un état d'un système n'est pas accessible à la mesure, on construit un modèle dynamique qui permet de reconstruire l'état non mesurable à partir de mesures d'autres grandeurs accessibles. Ce modèle dynamique s'appelle observateur d'état de la grandeur. Cependant la construction d'un observateur d'état n'est pas toujours possible, et la possibilité d'en construire un s'appelle la condition d'observabilité de cette grandeur. Dans la présente invention, des grandeurs cinématiques (une accélération ou une vitesse ou une position ou une déformation) de points particuliers de l'ensemble électronique à surveiller sont accessibles à la mesure à l'aide de capteurs, par exemple d'accéléromètres. À l'opposé, des contraintes mécaniques, et plus encore des états de dommage, à des points critiques de l'assemblage électronique sont presque toujours inaccessibles à la mesure, sauf en laboratoire. Un observateur d'état des contraintes mécaniques aux points critiques utilisant la mesure d'une ou plusieurs grandeurs cinématiques en des points particuliers peut être établi, sous réserve d'observabilité. Le caractère observable des contraintes mécaniques dépend directement du ou des points particuliers où sont placés le ou les capteurs cinématiques. En général, il suffit de placer le ou les capteurs cinématiques hors des noeuds de vibration, soit de le(s) placer par exemple sur la fixation de l'ensemble électronique. Le besoin de savoir évaluer correctement le dommage en vibration sur des cartes électroniques embarquées est de plus en plus prononcé. Ce besoin se rapporte à un certain nombre de points critiques de ces cartes électroniques, auxquels la soumission à des vibrations est susceptible d'entrainer la défaillance de l'ensemble électronique. Par exemple, les points critiques à surveiller sont des joints de brasure de certains composants électroniques. Dans le cas de telles cartes électroniques, des systèmes de surveillance embarqués existent. Ils peuvent être munis d'accéléromètres ou de jauges de déformation pour mesurer des accélérations ou des déformations subies par cette carte. Toutefois, ces capteurs ne peuvent presque jamais être placés aux points critiques de la carte à cause des très petites dimensions des composants, et plus encore, des joints brasés. Les systèmes existants ne peuvent donc pas accéder à la donnée des contraintes mécaniques en ces points critiques. On sait mesurer un état de vieillissement général de la carte, mais on ne sait pas mesurer d'état de fatigue local d'un composant, faute de pouvoir y mettre un capteur qui renseigne sur les contraintes subies par ce composant. Des méthodes ont été proposées par le laboratoire CALCE, Center for Advanced Life Cycle Engineering, pour faire le lien entre mesure et contraintes aux points critiques. Mais la grande lacune de ces méthodes est de ne pas prendre en compte les effets dynamiques des vibrations. En effet, la méthode proposée par le CALCE néglige la plupart des modes propres de vibration et repose uniquement sur des calculs statiques de courbure. À l'opposé, la présente invention repose sur des observateurs d'état qui rendent compte des effets dynamiques liés aux forces d'inertie et aux modes propres de vibration de l'ensemble électronique. Le système proposé par l'invention exploite beaucoup plus largement les mesures cinématiques pour évaluer le plus fidèlement possible les contraintes mécaniques aux points critiques. Le caractère multiaxial des contraintes mécaniques et la dynamique de l'ensemble électronique sont respectés. L'avantage d'évaluer précisément ces contraintes mécaniques est de pouvoir calculer beaucoup plus précisément le dommage subi par l'ensemble électronique, à l'endroit où sont placés le ou les composants critiques. L'utilisation d'un tel dispositif permet d'améliorer les prévisions de maintenance, car l'indication de dommage qu'il délivre sert à anticiper les défaillances potentielles. L'approche proposée permet d'appliquer les principes de la surveillance des structures mécaniques et le traitement numérique des mesures pour connaitre l'état de dommage en plusieurs endroits, grâce à un dispositif embarqué, avec un ou plusieurs capteurs. Le dispositif intégrant le procédé complet est un petit système électronique embarqué de surveillance (PHM pour Prognostics and Health Management) pourvu d'un ou plusieurs capteurs cinématiques (accélération ou vitesse ou position ou déformation), et d'une unité de calcul. Le système de surveillance est destiné à estimer le dommage lié à une défaillance potentielle des cartes électroniques constituant l'ensemble électronique à surveiller. Le ou les capteurs mesurent une grandeur cinématique en un ou respectivement plusieurs points de l'ensemble électronique à surveiller. Une des fonctions du système de surveillance est l'évaluation, à partir de la ou des mesures du ou des capteurs, de la contrainte mécanique en vibration qui est à l'origine d'une défaillance potentielle. Cette fonction est réalisée par un observateur d'état qui permet de reconstruire la dynamique de la contrainte mécanique à partir d'un modèle dynamique et des mesures cinématiques. La fonction s'opère en temps réel. Le ou les points de mesure et les points d'évaluation des contraintes ne sont pas forcément les mêmes. L'utilisation d'un observateur d'état ne requiert qu'une puissance de calcul de l'unité assez limitée, et le dispositif complet peut ainsi être architecturé autour d'un simple microcontrôleur par exemple. C'est donc le procédé même d'évaluation du vieillissement par fatigue qui permet l'intégration dans un petit dispositif embarqué disposant d'une capacité de calcul autonome. En pratique, l'observateur d'état peut être réalisé par exemple par une opération de convolution de l'accélération mesurée avec une réponse impulsionnelle en contrainte mécanique. De cette façon, l'ensemble de la bande de fréquence sollicitée est pris en compte, restituant ainsi tous les effets dynamiques. De plus, le caractère multiaxial de la contrainte mécanique est respecté dès lors qu'on utilise une réponse impulsionnelle pour chacune des composantes du tenseur des contraintes. L'observateur d'état pourrait être réalisé autrement que par une opération de convolution : par exemple avec un schéma de résolution d'équation différentielle ou encore en utilisant un modèle à dérivée non-entière (thèse de doctorat de Xavier Moreau, No 95BOR10512, 1995). L'estimation du dommage repose sur la contrainte mécanique donnée par l'observateur d'état. Par exemple, l'estimation du dommage peut être réalisée par comptage rainflow (décrit dans la norme NF-A03-406, � Fatigue sous sollicitations d'amplitude variable, méthode rainflow de comptage des cycles �, novembre 1993, ISSN 0335-3931), et par application de la règle de Palmgren-Miner. La puissance de calcul requise pour ces opérations est très faible, ce qui permet aussi l'intégration dans un dispositif embarqué disposant d'une capacité de calcul autonome. L'invention a donc pour objet un procédé d'évaluation du vieillissement d'un ensemble électronique, typiquement une carte électronique, soumis à un mouvement vibratoire, procédé dans lequel : ù on mesure une ou plusieurs grandeurs cinématiques, par exemple l'accélération ou la déformation, en un ou plusieurs points particuliers subissant le mouvement vibratoire, par exemple à l'endroit ou aux endroits de sa fixation à un ou des supports, ù on établit un modèle dynamique reliant les grandeurs cinématiques mesurées et des contraintes mécaniques subies en des points critiques au regard de la fatigue en vibration, ce modèle étant appelé observateur d'état, ù on déduit des contraintes mécaniques calculées selon l'observateur d'état un état de dommage de chacun de ces points critiques, états de dommage caractérisant le vieillissement recherché, ù on intègre la chaîne mesures cinématiquesûobservateur d'état calcul des états de dommage dans un dispositif embarqué disposant d'une capacité de calcul autonome. L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. Celles ci ne sont présentées qu'à titre indicatif et nullement limitatif de l'invention. Les figures montrent : ù figure 1 : un schéma synoptique du procédé de l'invention, ù figure 2 : un exemple pratique d'utilisation de l'invention, ù figure 3 : un exemple d'observateur d'état d'une contrainte mécanique d'un joint brasé d'un composant électronique critique, tiré d'un cas réel. La figure 1 montre qu'un ensemble électronique, par exemple une carte électronique 7, est soumis par un environnement 6 à des vibrations mécaniques. Selon l'invention, le procédé de surveillance, et donc d'évaluation du vieillissement, est mis en oeuvre à l'aide d'un système 3 comportant un capteur cinématique 1, une fonction dite observateur d'état 4 qui traite les signaux y produits par le capteur 1 pour restituer une estimation & de la contrainte mécanique vraie a en un point critique de l'ensemble 7, et une fonction 5 de calcul de dommage en fatigue accumulé en ce point critique de l'ensemble 7 prenant comme variable les résultats ô- produits par l'observateur d'état 4. L'observateur d'état 4 est capable de simuler les conséquences de l'excitation de l'environnement 6 à un endroit choisi alors que le capteur 1 est situé à un autre endroit a priori. Par exemple, figure 2, l'ensemble 7 est une carte électronique munie de composants utiles 8 à 10, en particulier de composants critiques au sens où leur fonctionnement doit être garanti, et donc surveillé, et sensibles au sens où des vibrations sont susceptibles d'entrainer une défaillance. La carte 7 est fixée dans son support d'utilisation par des fixations telles que 11 à 12. Des capteurs 1 à 2, des accéléromètres ici, mesurent les effets d'excitations el à e2 appliquées aux fixations 11 à 12. Selon l'invention, avec l'observateur d'état 4, on est capable d'évaluer les contraintes mécaniques aux endroits sur la carte 7 où sont placés les composants 8 à 10. Pour simplifier l'explication, on retiendra que les excitations el à e2 sont orientées dans une même direction, perpendiculairement à la carte 7, mais ne sont pas nécessairement synchrones. Toutefois, il serait possible de construire un observateur d'état même si cette obligation n'était pas respectée, mais les vibrations selon des directions dans le plan de la carte 7 causent des déformations de beaucoup plus faibles amplitudes qui causent moins de dommage aux composants électroniques. En principe c'est le support qui vibre et qui transmet ses vibrations à la carte par les fixations 11 à 12. Dans un exemple la vibration perpendiculairement à la carte 7 entraine majoritairement des contraintes de traction/compression, ou contraintes normales, sur des joints de brasure 0 d'un composant électronique critique 8 (type BGA, Bali Grid Array) qui conduisent potentiellement à la rupture en traction/compression de ce joint et à la défaillance de la carte électronique 7. Dans l'invention, en particulier si l'accéléromètre est un accéléromètre à trois dimensions, on peut calculer les dommages résultant des trois directions simultanément. Pour simplifier l'explication, on retiendra aussi que l'ensemble électronique 7 réagit de façon linéaire aux excitations vibratoires. Cela signifie que toute grandeur physique relative à la vibration se comporte de façon additive lorsqu'on superpose plusieurs excitations vibratoires appliquées en des endroits différents ou confondus. Cette restriction n'est pas limitative, car presque tous les ensembles électroniques se comportent ainsi dans la pratique. Un observateur d'état de la contrainte de traction/compression a au point critique 0 est réalisé par exemple par une opération de convolution entre les mesures d'accélération y, et 72 des points particuliers 1 et 2, et des réponses impulsionnelles s, et s2 dont la détermination est précisée plus loin. L'opération de convolution réalisant l'observateur d'état s'écrit : n-1 n-1 5-(t) 51(kT)'yl(t ù kT) +~ s2(kT)ry2(t ù kT) k=O k=O n étant un nombre de points nécessaires et suffisants pour représenter les réponses impulsionnelles s, et s2, k étant un nombre entier prenant successivement toutes les valeurs de 0 à n-1, i étant un pas de temps d'échantillonnage des mesures cinématiques, et en notant â (t) la contrainte normale évaluée au point critique 0 à l'instant présent t, y4(tûkt) la k-ième dernière valeur mesurée de y, suivant l'échantillonnage au pas de temps t, y2(tûkt) la k-ième dernière valeur mesurée de y2 suivant l'échantillonnage au pas de temps i, S1(kt) la k-ième valeur tabulée de la réponse impulsionnelle s, selon l'échantillonnage au pas de temps i, s2(kt) la k-ième valeur tabulée de la réponse impulsionnelle s2 selon l'échantillonnage au pas de temps i, et en utilisant la notation bien connue E désignant une sommation discrète selon une variable discrète et des bornes indiquées en indice et exposant de cette lettre. La variable indiquant le temps présent t est elle-même échantillonnée au pas de temps t et prend donc des valeurs discrètes multiples de T. L'observateur présenté ici délivre donc à chaque instant t multiple de t une estimation ô de la contrainte mécanique vraie a au point critique 0 dont la mesure directe est presque toujours impossible. La figure 3 montre un exemple d'observateur d'état sous forme de convolution tiré de l'étude d'une carte électronique réelle. Dans cet exemple, il n'y a qu'une seule excitation vibratoire el qui s'applique de façon synchrone sur tous les points de fixation de la carte dans la direction perpendiculaire au plan de la carte. Cette excitation est donc caractérisée par un nombre de degrés de liberté égale à 1. Un seul capteur cinématique est nécessaire et suffisant pour construire un observateur d'état pour chaque composante du tenseur des contraintes en un point critique, sous réserve d'observabilité liée aux positions respectives du capteur et du point critique et aux caractéristiques du capteur. Dans l'exemple de la figure 3, le capteur cinématique est un accéléromètre qui mesure l'accélération y perpendiculairement au plan de la carte, et la contrainte mécanique évaluée & est la contrainte normale dans un joint brasé d'un composant électronique de type BGA. Dans cet exemple, le nombre n de points représentant la réponse impulsionnelle est 512, le pas de temps est i = 0,5 ms et donc la convolution porte sur une plage de temps de 0,256 s. La dynamique de vibration de la carte électronique est prise en compte sur une plage de fréquence de 0 Hz à 1000 Hz, compte tenu du pas de temps d'échantillonnage t = 0,5 ms. Dans cet exemple tiré d'un cas réel, le nombre d'opérations arithmétiques par seconde à réaliser par l'unité de calcul du dispositif de surveillance (PHM) pour la fonction d'observateur d'état du point critique 0 est 1 024 000 multiplications par seconde et 1 024 000 additions par seconde. La capacité de calcul requise pour la fonction d'observateur d'état est donc assez modeste, selon les architectures électroniques actuelles. Le dispositif se passe donc complètement de tout post-traitement ou télémesure, car la technique de calcul utilisée fournit le résultat au prix d'un nombre limité et assez faible de calculs. On peut remarquer que le nombre de capteurs nécessaires doit généralement être supérieur ou égal au nombre de degrés de liberté d'excitation. Si le nombre de capteur est inférieur, la condition d'observabilité n'est généralement pas vérifiée. Ceci éclaire l'exemple donné en figure 2 avec deux capteurs cinématiques 1 et 2 pour deux excitations vibratoires indépendantes el et e2. Comme annoncé précédemment, un exemple de détermination des réponses impulsionnelles s1 et s2 est donné ci-dessous pour l'exemple de la figure 2. Dans la suite, chaque grandeur portant un tilde (-) désigne la transformée de Fourier d'une grandeur temporelle. La pulsation, qui est la variable d'une transformée de Fourier, est notée co et est homogène à l'inverse d'un temps. L'unité de co est usuellement le radian par seconde. Ainsi, la fonction 'l(w) résulte de la transformation de Fourier de la fonction y1(t), accélération délivrée par le capteur 1. La fonction 'y2(w) résulte de la transformation de Fourier de la fonction y2(t), accélération délivrée par le capteur 2. La fonction éi(w) résulte de la transformation de Fourier de la fonction el (t), excitation s'appliquant à la fixation 11. La fonction e2(w) résulte de la transformation de Fourier de la fonction e2(t), excitation s'appliquant à la fixation 12. La fonction 5-o(w) résulte de la transformation de Fourier de la fonction ao(t), contrainte mécanique au point critique 0. On note que les transformées de Fourier sont des grandeurs scalaires complexes. Quatre fonctions H11(co), H12(o), H21(o), H22(o) de la pulsation co sont définies chacune comme la réponse d'un des points de mesure en terme d'accélération, à une excitation unitaire de pulsation co en un des points d'excitation, toutes les autres excitations étant maintenues identiquement nulles. Par la linéarité précisée ci-dessus, ces fonctions s'écrivent de façon générale comme le rapport d'une accélération sur une excitation non nulle qui en est la cause, les autres excitations étant maintenues identiquement nulles. Ce rapport est indépendant l'excitation appliquée à cause de la linéarité précisée plus haut. De façon explicite : Hii(w) = 7i(w) ei (w) e2=0 e2=0 H21(w) ù 72(w) ei (w) e1=0 H22(W) ù ï2(W) e2(W) La notation en indice e1=0 signifie que le rapport est calculé pour l'excitation e2 appliquée seule. Ces quatre fonctions sont indépendantes des excitations appliquées. Elles peuvent être déterminées expérimentalement, par exemple en procédant successivement aux excitations unitaires dont il est question ci- dessus, en balayant successivement des valeurs de pulsation co. La méthode des éléments finis, technique de calcul très usuelle et reconnue dans l'industrie, le génie civil et les sciences, offre également des moyens multiples pour déterminer ces quatre fonctions, par exemple en calculant numériquement la procédure expérimentale qui permet de déterminer ces quatre fonctions. De façon analogue, deux fonctions G1(co) et G2(co) de la pulsation co sont définies chacune comme la réponse du point critique 0 en terme de contrainte mécanique, à une excitation unitaire de pulsation co en un des points d'excitation, toutes les autres excitations étant maintenues identiquement nulles. Par linéarité, comme précédemment : Gl(W) = 5-0(w) et (W ) e1=0 La notation en indice e1=0 ou e2=0 signifie que le rapport est calculé pour une excitation appliquée seule. Les deux fonctions G1 et G2 peuvent être déterminées expérimentalement ou par la méthode des éléments finis, comme pour les quatre fonctions précédentes. Les transformées de Fourier 'ïl(W) et `}'2(W) des deux accélérations y1(t) et y2(t) mesurées par les capteurs 1 et 2 s'écrivent en superposant linéairement les réponses aux deux excitations el et e2 : J ;fi. (W) = H11(W)ê1(W) + H12(w)ê2(w) Y2(W) = H21(W)él(W) + H22(W)é2(W) 9 La notation en indice e2=0 signifie que le rapport est calculé pour l'excitation el appliquée seule. H12(w) = ïl(Lu) e2(W) e1=0 e2=0 G2(W) = 5-0(W) e2(W) Ce système scalaire s'écrit sous forme matricielle : }1(W) H11(W) H12(w) é1(W) 2( H21(W) H22(W) éz(W) Le système matriciel inverse s'écrit : 10 où l'exposant -1 indique l'inverse de la matrice. Soit, sous forme de système scalaire : el(W) _ ( H11(W) H12(W) 1 ê2 (w) / ù \ Hz 1(W) H22 (W) ) \ Î'z (Wé1(W) = H22(W)'i'1(W) ù H1z(W);ï2(W) H11(W)H22(W) ù H12(W)H21(W) ùH21(W)`}'1(W) + H11(W);ï2(W) H11(W)H22(W) ù H12(W)H21(W) 15 La transformée de Fourier 5-o(W) de la contrainte mécanique 6o au point critique 0 s'écrit en superposant linéairement les réponses aux deux excitation e, et e2 : Qo((,J) = Gl((,J)él((,J) + G2(W)é2(W) Cette expression s'écrit alors selon les transformées de Fourier des accélérations : Gl(W)Hzz(W) ù G2(W)H21(W) %Y1(W) H11(W)H22(W) ù H12(W)H21(W) ùG1(W)H12(W) + G2((,J)Hll(W) +H11((,J)Hzz(W) ù H12(W)H21(w)Î'2(W) On définit alors les deux facteurs devant '1(W), respectivement 72(W) dans l'expression qui précède comme deux fonctions s1(W) et respectivement sz(w) : 20 25 'é-o(W) = S1(w) _ G1(w)H224 ) ù G2(w)H21(w) H11(i)H22(ù) ù H12(w)H21(w) S2(w) G1(w)H12(w) + G2((,J)Hll((,J) H11(i)H22(ù) ù H12(w)H21(i) Les réponses impulsionnelles si(t) et s2(t) sont simplement les transformées de Fourier inverses de s1(w) et s2(w) respectivement. On note que si(t) et s2(t) sont des grandeurs scalaires réelles, par leur nature physique. Il y a une raison algébrique à cela qui tient aux propriétés de la transformation de Fourier. si(t) et s2(t) sont les deux fonctions qui, convoluées respectivement avec les mesures y1(t) et y2(t), restituent la contrainte mécanique ao(t) au point critique 0 de l'ensemble électronique 7 : 0-o(t) = f t si(u)'Yi(t ù u)du +1 t s2(u)'Y2(t ù u)du o o La variable u sous le signe intégrale est la variable muette d'intégration. Une version discrétisée de cette équation intégrale est l'expression donnée plus haut de l'observateur d'état -(t). Pour obtenir une telle version discrétisée, il suffit par exemple de substituer les intégrales par des sommations discrètes et de n'utiliser que les valeurs de s1 et s2 prises aux points multiples du pas de temps i: nù1 nù1 Q(t) sl(kT)'yl(t ù kT) + s2kT)-ï2(t ù kT) k=0 k=0 La discrétisation n'est valable qu'avec un pas de temps suffisamment petit pour représenter correctement les intégrales et les réponses s1 et s2, comme l'illustre la figure 3. Bien entendu, dans tout ce qui précède, le dénominateur H11(w)H22(w) û H12(w)H21(w) ne doit s'annuler pour aucune des valeurs de pulsation co dans la plage de fréquence des vibrations. Si c'était le cas, la condition d'observabilité ne serait pas vérifiée : il ne serait pas possible de déduire les excitations à partir des mesures cinématiques, et il faudrait positionner les capteurs cinématiques différemment. Le calcul précédent n'est qu'une illustration d'une façon de mettre au point un observateur d'état pour une contrainte mécanique d'un point critique. Elle est immédiatement généralisable à un seul capteur ou à un nombre de capteurs supérieur à deux, toujours sous réserve d'observabilité, donc sous réserve que le nombre de voies de mesure soit au moins égal au nombre de degrés de liberté des excitations. Le calcul par éléments finis repose sur le maillage de la structure géométrique de l'ensemble électronique en éléments polyédriques et sur l'écriture des équations de la mécanique en vibration selon la méthode de Galerkin, appelée aussi formulation faible. La méthode des éléments finis est abondamment documentée dans les revues de mathématiques appliquées et de mécanique. Le calcul peut être réalisé avec l'un des nombreux logiciels du marché. Une présentation de cette technique est disponible dans un article intitulé � Durability Modelling of a BGA Component under Random Vibration � et présenté au congrès EuroSimE à Freiburg-im-Breisgau les 21- 23 avril 2008 et accessible sur le site Internet http://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=4525047&isnumber=4 525005. Le calcul 5 du dommage est réalisable par exemple selon la méthode rainflow, bien connue (norme NF-A03-406), lorsque la contrainte mécanique évaluée par l'observateur d'état est scalaire. Dans cette méthode, tous les extréma locaux de la contrainte mécanique sont identifiés et un algorithme permet d'identifier des cycles de contrainte mécanique, chaque cycle se définissant par l'un des minima locaux et par l'un des maxima locaux de la contrainte mécanique. Chaque cycle contribue au dommage du point critique 0. Un exemple de calcul du dommage accumulé consiste à utiliser l'hypothèse de Palmgren-Miner, selon laquelle le dommage accumulé est une grandeur scalaire s'incrémentant à chaque cycle d'une valeur qui ne dépend que des caractéristiques du cycle considéré, indépendamment des cycles qui l'ont précédé. Lorsque le dommage attend une valeur seuil, on peut considérer qu'une défaillance du point critique 0 a une probabilité notable de se produire et le dispositif de surveillance est alors en mesure de générer une alarme. Alors, des mesures préventives peuvent être prises avant que la défaillance considérée ne se produise. Lorsque la contrainte mécanique est plus généralement tensorielle, il existe plusieurs méthodes. L'une d'elle est présentée dans l'article � Durability Modelling of a BGA Component under Random Vibration � précité. Autrement, il existe aussi des méthodes rainflow généralisées aux grandeurs tensorielles. Enfin, on peut citer des critères multiaxiaux de Dang Van et de Sines : voir � Mécanique des matériaux solides �, J. Lemaître, J.-L. Chaboche, ISBN 2-10-001397-1. REVENDICATIONS1 - Procédé d'évaluation du vieillissement par fatigue en vibration d'un ensemble électronique, typiquement une carte électronique, dans lequel - on mesure une ou plusieurs grandeurs cinématiques, par exemple l'accélération ou la déformation, en un ou plusieurs points particuliers subissant le mouvement vibratoire, par exemple à l'endroit ou aux endroits de sa fixation à un ou des supports, - on établit un modèle dynamique reliant les grandeurs cinématiques mesurées et des contraintes mécaniques subies en des points critiques au regard de la fatigue en vibration, ce modèle étant appelé observateur d'état, - on déduit des contraintes mécaniques calculées selon l'observateur d'état un état de dommage de chacun de ces points critiques, états de dommage caractérisant le vieillissement recherché, - on intègre la chaîne mesures cinématiquesûobservateur d'étatûcalcul des états de dommage dans un dispositif embarqué disposant d'une capacité de calcul autonome. 2 - Procédé selon la revendication 1, caractérisé en ce que - un capteur cinématique est placé en un endroit de l'ensemble électronique qui n'est ni une fixation de l'ensemble électronique au support, ni une point critique à surveiller. 3 - Procédé selon la revendication 1, caractérisé en ce que l'observateur d'état est déterminé par calcul selon la méthode des éléments finis. 4 - Procédé selon l'une des revendications 1 à 3, caractérisé en ce qu'on place le ou les capteurs en un ou des endroits qui permettent l'observabilité des contraintes mécaniques aux points critiques en se fondant sur l'analyse vibratoire de l'ensemble électronique au moyen de la méthode des éléments finis. 5 - Procédé selon l'une des revendications 1 à 4, caractérisé en ce que - le nombre de points de mesure est supérieur ou égal au nombre de degrés de liberté des excitations.