Méthode FORM-SORM

La méthode FORM-SORM est une méthode utilisée en fiabilité pour déterminer la probabilité de défaillance d'un composant. FORM est l'acronyme de First Order Reliability Method, méthode de fiabilité du premier ordre ; SORM est l'acronyme de Second Order Reliability Method, méthode de fiabilité du second ordre. Il s'agit de trouver le point de conception, c'est-à-dire le cas de défaillance du maximum de vraisemblance. La solution est obtenue par une procédure d'optimisation.

Présentation du problème

Une conception « au plus juste » — compromis entre le coût de fabrication, le coût de service (maintenance, garantie) et la satisfaction de l'utilisateur (sécurité, disponibilité du système) — nécessite d'avoir recours aux statistiques ; c'est le domaine de la fiabilité. On veut connaître la probabilité que la durée de vie T d'un système soit supérieure à une valeur t, R(t) = P(T ≥ t), ce qui revient à connaître la probabilité de défaillance avant t, F(t) = P(T ≤ t) ; F = 1 - R.

En particulier, on veut garantir que cette probabilité de défaillance est inférieure à une limite α, choisie en fonction des incertitudes et des conséquences d'une situation de défaillance : F(t) ≤ α, c'est-à-dire R(t) ≥ 1 - α. On se place en général à la limite F(t) = α. La notion de durée de vie implique celle de processus, c'est-à-dire de variable aléatoire indexée sur le temps, ce que la méthode FORM présentée ne prend pas en compte. Cette durée de vie dépend de plusieurs facteurs : des facteurs de charge, ou contraintes S (stress) : c'est ce que subit la machine, sollicitations mécaniques, humidité, température, intensité du courant électrique… des facteurs de fabrication, ou résistance R : dimensions des composants, propriétés de la matière, … Tous ces facteurs, contraintes ou résistances, peuvent s'exprimer par des grandeurs chiffrées xi. Ces valeurs sont des réalisations de variables aléatoires Xi : tous les systèmes, même produits à la chaîne, sont différents les uns des autres, et ils ne subissent pas non plus les mêmes sollicitations.

Dans un certain nombre de cas, on peut définir une fonction de défaillance g(x1, x2, …, xn) pour exprimer la condition sur la durée de vie : F(t) ≤ α ⇔ g(x1, x2, …, xn) ≥ 0. L'espace à n dimensions (x1, x2, …, xn) est donc coupé en deux demi-espaces : le domaine de défaillance F (failure), vérifiant l'équation g(x1, x2, …, xn) < 0 ; le domaine de bon fonctionnement, vérifiant l'équation g(x1, x2, …, xn) > 0 ; la frontière entre les deux est une hypersurface d'équation g(x1, x2, …, xn) = 0 qui définit l'état limite.

Chaque variable aléatoire Xi suit une loi de probabilité dont la fonction de densité de probabilité est notée ƒi. La densité de probabilité d'avoir un n-uplet donné (x1, x2, …, xn) vaut donc y = ƒ1(x1)׃2(x2)ׅ׃n(xn) La fonction y(x1, x2, …, xn) forme une hypersurface de l'espace (x1, …, xn, ƒ). Le maximum de y, le « sommet » de l'hypersurface, est le point le plus probable.

Pour déterminer la probabilité de défaillance Pf du système, il faut intégrer y sur la zone de défaillance F

La valeur de Pf est le « volume » — l'hypervolume — délimité par les hypersurfaces y et g = 0.

La solution analytique de cette intégrale est en général complexe, voire impossible. Une manière d'estimer Pf consiste à utiliser une méthode de Monte-Carlo : on génère un grand nombre de valeurs aléatoires (x1, x2, …, xn) suivant les lois statistiques connues, et l'on compte le nombre de cas pour lesquels g est négatif.

Toutefois, la majeure partie de l'information est contenue dans une petite zone de l'espace, dite « zone critique », qui est l'ensemble des n-uplet (x1, x2, …, xn) pour lesquels la densité de probabilité de défaillance est importante ; le point où la densité de probabilité de défaillance est maximale est appelée « point de conception » et noté x* (c'est un n-uplet). Or, si la conception est robuste, cette zone critique est loin du sommet de l'hypersurface y, puisque l'on veut que la plupart des situations soient des situations de bon fonctionnement. Donc, la méthode de Monte-Carlo génère beaucoup de n-uplet proche du maximum de y, c'est-à-dire loin de la zone d'intérêt.

Il faut donc un nombre considérable de simulations pour avoir une bonne estimation de Pf.

La méthode FORM-SORM est une méthode d'approximation consistant à trouver le point de conception, et de s'en servir pour déterminer la probabilité de défaillance Pf.

Méthode générale

La méthode FORM-SORM consiste à : Remplacer les variables aléatoires par des lois normales centrées réduites. Trouver le n-uplet (x1, x2, …, xn) de l'hypersurface d'état limite le plus proche de l'origine ; c'est alors le point de conception x*, le n-uplet de la zone de défaillance ayant la densité de probabilité maximale. Faire une approximation linéaire (FORM) ou quadratique (SORM) de g pour déterminer la probabilité de défaillance Pf.

Hypothèses

On suppose que les variables aléatoires ont une distribution normale.

Si les variables sont normales d'espérance μi et d'écart type σi, on les transforme également en variables de loi centrées réduites Ui = (Xi - μi)/σi Xi = σiUi + μi

Si les variables aléatoires Xi ne sont pas normales, on les transforme pour avoir des variables Ui suivant lois normales centrées réduites ; cette opération est la « transformation de Nataf et Rosenblatt » : Ui = Φ-1(Fi(Xi)) et réciproquement Xi = F-1(Φi(Ui)) où Φ (lettre grecque Phi) est la fonction de répartition de la loi normale centrée réduite, et Fi la fonction de répartition de Xi.

Méthode générale

Plus on est proche du sommet de l'hypersurface y, plus la densité de probabilité est élevée. Les fonctions de densité ƒi étant des courbes en cloche symétriques, y forme une « colline » : plus on s'éloigne du sommet, plus la densité de probabilité y diminue. Donc, le point le plus critique, le point de conception x*, est le point de la zone de défaillance F le plus proche du sommet.

La méthode FORM-SORM consiste à trouver le point de la zone frontière g(x1, …, xn) = 0 le plus proche du sommet. Dans l'espace centré réduit, le sommet de y se trouve à l'origine, le point de coordonnées (0, 0, …, 0). La fonction de défaillance devient alors g* : F(t) = α ⇔ g*(u1, …, un) ≥ 0. La distance à l'origine est donc

et ainsi, la méthode consiste à minimiser d*(u1, …, un) avec pour contrainte g*(u1, …, un) = 0. C'est une méthode d'optimisation quadratique : en effet, minimiser d* revient à minimiser d*2 (puisque la fonction carré est strictement croissante sur [0 ; +∞[), on doit donc minimiser une fonction quadratique ∑xi2 avec une contrainte linéaire.

La distance minimale trouvée est appelée indice de fiabilité et est noté β (lettre grecque bêta) : β = d*(x*).

Méthode FORM On fait l'hypothèse que la fonction g est linéaire. Si ce n'est pas le cas, la méthode revient finalement à en prendre le développement limité du premier ordre.

La probabilité de défaillance vaut alors Pf = Φ(-β).

Méthode SORM

La méthode SORM est identique à la méthode FORM, mais si g* n'est pas linéaire, alors on utilise une approximation du second ordre pour déterminer Pf.

On définit les courbures principales κi (lettre grecque kappa) :

Elles sont au nombre de n - 1, puisque l'équation g* = 0 définit une hypersurface de dimension n - 1.

La probabilité de défaillance s'écrit alors : Pf = S1 + S2 + S3 avec :

où Re est la partie réelle, i est l'imaginaire, et φ est la densité de probabilité de la loi normale centrée réduite.

La fonction S1 est l'approximation asymptotique de Pf lorsque β tend vers +∞ — c'est donc le terme prépondérant si la conception est robuste (le point de conception est loin de l'origine) —, et S2 et S3 sont des termes correctifs.

Application à un problème à deux variables normales

Présentation du problème

vignette|300px|Distribution de la contrainte S et de la résistance R : toutes les situations (R, S) ont une densité de probabilité (surface en niveaux de gris). La zone où la marge m = R - S est positive est la zone où le système survit (R > S).

vignette|300px|Même situation vue de dessus.

vignette|300px|Point de conception : maximum de la courbe m = 0.

Dans la méthode contrainte-résistance, la condition sur la durée de vie F(t) ≤ α s'exprime par une résistance R supérieure à une contrainte S (stress) : R ≥ S ou plutôt en définissant la marge m = R - S : m ≥ 0. Cette marge m est donc une fonction de défaillance : g(R, S) = m.

Si l'on dispose des distributions (densités de probabilité) ƒR et ƒS, on peut alors déterminer cette probabilité : F(t) = P(R ≤ S) = P(m ≤ 0) La probabilité d'avoir une couple de valeurs dans un intervalle ([S, S + dS], [R, R+dR]) est donné par P([S, S + dS]∩[R, R+dR]) = P([S, S + dS])×P([R, R+dR]) = ƒS׃R×dS×dR. On peut tracer la surface y = ƒS׃R, ainsi que la ligne de séparation m = 0 : cette ligne marque l'état limite entre la situation de survie et la situation de défaillance. Une méthode simple d'estimer la probabilité P(m ≤ 0) consiste à utiliser une méthode de Monte-Carlo : on effectue N tirages aléatoires de couples (R, S), et l'on dénombre les n situations pour lesquelles m ≤ 0 : P(m ≤ 0) ≃ n/N. Si l'on veut un risque α faible, cela signifie que le sommet de la surface doit être loin de l'état limite m = 0.

Le point que l'on veut le mieux connaître est le maximum de la courbe m = 0 ; c'est notre « point de conception ». Si le risque α est faible, alors le point de conception est loin du sommet de la surface. Cela signifie que lors du tirage aléatoire, peu de points seront proches du point de conception ; il faudra donc de très nombreux tirages — N devra être très grand — si l'on veut caractériser ce point.

Exemple

Nous avons représenté ci-contre une situation pour laquelle : la contrainte S suit une loi normale d'espérance μS = et un écart type σS = ; la résistance R suit une loi normale d'espérance μR = et un écart type σR = . Graphiquement, on constate que le maximum de l'état limite se trouve à . Ce point est en général proche de l'intersection des courbes ƒS et ƒR.

Une méthode de Monte-Carlo avec un million de tirages et un intervalle de confiance de 5 %l'erreur pour un risque α avec N tirages vaut soit donne 18,1 % ± 0,2 % de défaillance.

Le logiciel LibreOffice Calc est un tableur libre et gratuit.

La fonction ALEA (ou RAND en anglais)Fonctions mathématiques (aide lLibreOffice) donne un nombre aléatoire compris entre 0 et 1, mais avec une loi uniforme. Pour avoir une loi normale, il faut appliquer la fonction LOI.NORMALE.INVERSE (NORMINV)Fonctions statistiques - Quatrième partie (aide LibreOffice 3.3). On compare simplement deux valeurs avec l'opérateur >, ce qui donne un booléen ; Calc calcule la moyenne (commande MOYENNE, MEAN) de booléens en prenant 1 pour « vrai » et 0 pour « faux ». Il faut ensuite étirer une boîte de… un million de lignes, ce qui n'est pas raisonnable, et touche les limites du logicielQuel est le nombre maximum de cellules dans une feuille de calcul LibreOffice ?, FAQ LibreOffice Calc.

On se limite ici à mille lignes, ce qui donne un résultat à ± 5 %.

A B C D 1 R S 2 mu 250 230 3 sigma 17 14 4 5 n° Valeur R Valeur S S>R 6 1 =LOI.NORMALE.INVERSE(ALEA() ; B$2 ; B$3) =LOI.NORMALE.INVERSE(ALEA() ; C$2 ; C$3) =C6>B6 … 1006 =MOYENNE(D6:D1005) 1007 1008 epsilon =RACINE(LN(2/0,005)/(2*1001))

Le logiciel R est un logiciel libre et gratuit, sous la forme d'un langage de programmation interprété, qui permet d'effectuer des calculs statistiques. La méthode de Monte-Carlo se fait de la manière suivante : On crée deux matrices R et S de un million de nombres tirés aléatoirement avec une loi normale (commandes rnorm()). On crée une matrice booléenne Mbool, contenant des « vrai » si la contrainte S est supérieure à la résistance R (opérateur >). On calcule la moyenne de cette matrice, les « vrai » étant comptés comme des 1 et les « faux » comme des 0 (commande mean()).

N = 1e6 # nombre d'événements (tirs) alpha = 0.005 # risque