FR 3029624 A1 20160610 FR 1462000 A 20141205 1 Procédé de suivi de la fabrication de pièces basé sur l'analyse d'indicateurs statistiques en situation d'allègement de contrôle DOMAINE DE L'INVENTION L'invention concerne le domaine des indicateurs statistiques en milieu industriel, par exemple dans l'industrie aéronautique, en particulier pour faciliter le suivi et le contrôle de la fabrication d'un flux de pièces. ETAT DE LA TECHNIQUE La fabrication de pièces, notamment mécaniques, en milieu industriel se heurte à deux contraintes opposées, à savoir l'augmentation des cadences et volumes de fabrication d'une part, et l'augmentation des niveaux de qualité requis d'autre part, ce qui est en particulier vrai dans le domaine aéronautique. Il est aujourd'hui impossible d'effectuer des contrôles qualité pour la totalité des pièces fabriquées sauf à nuire considérablement à la cadence de fabrication. On utilise donc en général des indicateurs statistiques de la fabrication, permettant de déduire une information globale sur la qualité de l'ensemble des pièces fabriquées à partir d'information spécifiques de la qualité d'un nombre fini de pièces prélevées en tant qu'échantillon. Outre les contrôles en fin de production qui peuvent être faits sur des échantillons ayant un nombre limité de pièces, il est généralement également effectué des contrôles en cours de production de manière à pouvoir éventuellement piloter le flux de production, c'est-à-dire ajuster les conditions de fabrication pour s'assurer que les pièces fabriquées continuent de répondre aux critères de qualités requis. Dans certains cas, ces contrôles statistiques en cours de production peuvent amener à stopper complètement la production, notamment si les pièces produites présentent des défauts de qualité trop importants et que le flux de fabrication doit être réinitialisé entièrement. Les contrôles qualité sont faits en relation à une ou plusieurs dimensions caractéristiques des pièces qui sont fabriquées. Cette dimension caractéristique peut être par exemple une cote particulière de la pièce, sa masse, ou toute autre caractéristique mesurable desdites pièces. Pour effectuer les contrôles statistiques, on prélève successivement plusieurs échantillons, chaque échantillon comprenant plusieurs pièces du flux de fabrication, et on mesure ensuite la dimension caractéristique de chaque pièce de l'échantillon prélevé. A partir des différentes mesures de la dimension caractéristique des pièces de l'échantillon 3029624 2 prélevé, on calcule la valeur d'un indicateur statistique préalablement choisi pour suivre la qualité du flux de fabrication. Il existe divers indicateurs statistiques pouvant être utilisés afin de suivre l'évolution d'un flux de fabrication de pièces, chaque indicateur statistique donnant des 5 informations différentes permettant d'ajuster les conditions de fabrication d'une manière ou d'une autre. La plupart des indicateurs statistiques utilisés pour le suivi d'un processus industriel de fabrication sont calculés à partir d'une moyenne p et d'un écart-type a de la dimension caractéristique mesurée sur plusieurs pièces. Plus précisément, p correspond 10 à la moyenne du décentrage mesuré pour la dimension caractéristique par rapport à la valeur de référence pour cette dimension caractéristique. Le processus de fabrication peut par exemple être piloté en étudiant des indices de capabilité qui caractérisent la performance réelle du procédé par rapport à la performance souhaitée. De tels indices permettent en effet de mesurer la capacité du 15 processus de fabrication à réaliser des pièces dont les dimensions caractéristiques sont comprises dans un intervalle de tolérance IT souhaité. L'intervalle de tolérance IT est l'écart entre les valeurs extrêmes admissibles de la dimension caractéristique, étant donc calculé comme la différence entre la tolérance supérieure TS et la tolérance inférieure TI de la dimension caractéristique mesurée, soit IT = TS - TI. 20 On peut par exemple se référer à l'indice de capabilité Cpk qui représente la dispersion mais également le centrage de la production par rapport aux limites de la tolérance. Dans ce cas, lorsque l'indice de capabilité Cpk est élevé, cela signifie que la production est répétable et qu'elle est en outre bien centrée dans l'intervalle de tolérance IT, c'est-à-dire qu'il y aura moins de risque que des pièces soient fabriquées en dehors 25 des tolérances. L'indice de capabilité Cpk est défini généralement par la formule : Min(TS - ; µ - TI) Cpk = 3a Le principal intérêt du Cpk est que lorsque la répartition de la production suit une loi normale, la valeur du Cpk de la population produite est directement liée à la proportion de non-conformités produites. Plus précisément, lorsque l'indice de capabilité mesuré pour l'ensemble d'une production donnée est égal à une valeur Cpk, la proportion 30 (exprimée en parties par million, notés ppm) de pièces produites avec une valeur supérieure à TS ou inférieure à TI pour la dimension caractéristique, c'est-à-dire la proportion de pièces en non-conformité, vaut au maximum : Pppm = 106[1 - LOI. NORMALE. STANDARD(3 * Cpk) + LOI. NORMALE. STANDARD(-3 * Cpk)] 3029624 3 Par exemple, un indice de capabilité Cpk égal à 1 correspond à environ 2700 ppm de pièces en non-conformité, tandis qu'un indice de capabilité Cpk égal à 1.33 correspond à environ 66 ppm de pièces en non-conformité. Comme indiqué plus haut, comme les contrôles sont coûteux et sans valeur 5 ajoutée sur la pièce, on ne contrôle en général qu'une partie de la production globale. Pour ce faire, on prélève un échantillon de n pièces, on mesure la dimension caractéristique sur ces n pièces, et calcule sur ces n valeurs mesurées un indice de capabilité Cpk qu'on espère être représentatif de celui de la population produite toute entière. Bien entendu, plus le nombre de pièces prélevées n est grand, et plus le degré de 10 confiance dans cette représentativité est élevé. Pour quantifier cette confiance, il existe des formules approchées donnant en fonction du nombre de pièces n mesurées et de l'indice de capabilité Cpk mesuré, que l'on note Cpk, les bornes de l'intervalle de confiance sur cette mesure d'indice de capabilité à un seuil de confiance a donné. 15 La plus couramment utilisée est la formule de Kushler, qui permet de calculer les bornes inférieures de l'intervalle de confiance sur l'indice de capabilité Cpk à un risque unilatéral a près (a E [0; 1]). Selon cette formule, donné un indice de capabilité mesuré sur n pièces, noté Cpk, on court un risque a que l'indice de capabilité Cpk de la population réelle soit inférieur 20 à Cpk [1 + za I9nCpk2 + 2(n où za représente le dème percentile de la loi normale centrée réduite. Par exemple, on a 5% de chance que l'indice de capabilité Cpk réel de la population soit inférieur à Cpk [1 - 1.645 i9nCpk2 2(n-1) Dans la pratique, on compare généralement la valeur calculée de la borne 25 inférieure de l'intervalle de confiance unilatéral à 5% sur l'indice de capabilité mesuré à la valeur de l'indice de capabilité minimum acceptable, noté Cpkmini, spécifiée par le concepteur ou les procédures de qualité de l'entreprise. Dans ces conditions, étant donnée la mesure d'un indice de capabilité noté Cpk sur un échantillon de n pièces, le niveau de qualité de la production est jugé satisfaisant, 30 au regard des incertitudes pesant sur la représentativité de l'échantillon prélevé à la condition que : Cpk [1 - 1.645 i9n-Cpk2 2(n-1)] > LpKmini. Ce mode de contrôle permet de fournir un cadre clair et simple sur la manière de procéder pour prendre en compte les incertitudes liées à l'échantillonnage dans la sanction finale (acceptable ou non acceptable) qui sera rendue à l'issue de l'analyse des 35 mesures effectuées, et permettra ou non la poursuite de l'allègement de contrôle, tout en participant au pilotage de l'outil de production. 3029624 4 Au moyen d'une sévérisation (le remplacement de la valeur Cpk par la borne inférieure de l'intervalle de confiance unilatéral à 5% sur cette mesure Cpk [1. - 1.645 \i9n-Clpk2 + 2(n1-1) on se prémunit à chaque analyse contre un risque de 5% que l'échantillon prélevé paraisse acceptable alors que la production réelle pourrait être non 5 acceptable. En d'autres termes, on se prémunit contre un risque de non-représentativité d'un échantillon qui fournirait une image de la production plus optimiste que la réalité effective des pièces livrées sans avoir été mesurées. Le problème de ce mode de fonctionnement est qu'il conduit à imposer des sévérisations très importantes. Par exemple lorsque la valeur minimale de l'indice de 10 capabilité Cpkmini vaut 1 (ce qui signifie que l'on tolère 2700 ppm de pièces en non- conformité) et que l'on souhaite avoir un allègement à 50 pièces par période (ce qui représente tout de même un nombre important de contrôles), la valeur de l'indice de capabilité mesuré Cpk suivant ce mode de fonctionnement sera jugée acceptable sous la condition que Cpk [1 - 1.645 i9n-Cpk2 + 2(n-1)1 > 1, soit Cpk > 1.22, ce qui correspond à 15 252 ppm de pièces en non-conformité. Ainsi, alors que la spécification tolère 2700 ppm de pièces en non-conformité, on devra produire continuellement à un taux de non-conformité dix fois inférieur pour pouvoir rester durablement allégeable au taux pourtant relativement élevé de 50 pièces par période. Un but de la présente invention est donc de fournir un procédé de suivi d'un flux 20 de fabrication de pièces basé sur l'analyse d'au moins un indicateur statistique, qui permet de remédier aux inconvénients précités. Plus spécifiquement, un but de la présente invention est de fournir un procédé de suivi d'un flux de fabrication de pièces basé sur l'analyse d'au moins un indicateur statistique qui permet un allègement de contrôle fiable et plus précis. 25 Encore un but de la présente invention est de fournir un procédé de suivi d'un flux de fabrication de pièces basé sur l'analyse d'au moins un indicateur statistique avec une sévérisation réduite, tout en garantissant un allègement de contrôle au moins aussi fiable. EXPOSE DE L'INVENTION 30 A cette fin, on propose un procédé de suivi d'un flux de fabrication d'une population de pièces basé sur l'analyse d'au moins un indicateur statistique représentatif d'une dimension caractéristique des pièces, selon lequel : a) On prélève dans le flux de fabrication de pièces un échantillon comprenant un nombre n de pièces ; 35 b) On mesure la dimension caractéristique de chaque pièce de l'échantillon, et on calcule pour l'échantillon une valeur mesurée de l'indicateur statistique ; 3029624 5 c) On calcule une espérance mathématique de la proportion de pièces en non-conformité par rapport à une spécification sur la dimension caractéristique, ledit calcul étant effectué à partir de la valeur mesurée de l'indicateur statistique pour l'échantillon prélevé et du nombre n de pièces de l'échantillon ; 5 d) On compare l'espérance mathématique de la proportion de pièces en non- conformité calculée à une valeur seuil de proportion de pièces en non-conformité ; e) On pilote le flux de production des pièces en fonction des résultats de comparaison de l'étape d). 10 Des aspects préférés mais non limitatifs de ce procédé, pris seuls ou en combinaison, sont les suivants : - à l'étape c), on calcule en outre un intervalle de confiance sur l'indicateur statistique à un risque unilatéral a près, où a est compris entre 0 et 1, ledit calcul étant effectué à partir de la valeur mesurée de l'indicateur statistique pour 15 l'échantillon prélevé et du nombre n de pièces de l'échantillon ; et à l'étape d), on compare en outre la borne inférieure de l'intervalle de confiance calculé à une valeur seuil de l'indicateur statistique correspondant à la valeur seuil de proportion de pièces en non-conformité. - à l'étape e) : 20 o Si la borne inférieure de l'intervalle de confiance calculé est supérieure à la valeur seuil de l'indicateur statistique, alors on n'effectue aucune modification au flux de fabrication de pièces ; o Si la borne inférieure de l'intervalle de confiance calculé est inférieure à la valeur seuil de l'indicateur statistique et que l'espérance mathématique de 25 la proportion de pièces en non-conformité calculée pour l'échantillon est supérieure à la valeur seuil de proportion de pièces en non-conformité, alors on effectue des actions correctives sur le flux de fabrication de pièces et/ou on augmente le nombre de pièces prélevées lors d'un prélèvement ultérieur d'échantillon ; 30 o Si l'espérance mathématique de la proportion de pièces en non-conformité calculée pour l'échantillon est inférieure à la valeur seuil de proportion de pièces en non-conformité, alors on stoppe le flux de fabrication de pièces. - le au moins un indicateur statistique est un indice de capabilité Cpk défini par la formule : Cpk = 3a Min(TS - ; µ - TI) 35 où : 3029624 6 o p est la moyenne des dimensions caractéristiques mesurées sur les pièces de l'échantillon prélevé ; o a est l'écart-type des dimensions caractéristiques mesurées sur les pièces de l'échantillon prélevé ; 5 o TS est une tolérance supérieure de la dimension caractéristique mesurée ; o Tl est une tolérance inférieure de la dimension caractéristique mesurée. - à l'étape c), on prend pour hypothèse que la population des pièces suit une loi normale, et on effectue les calculs suivants : cl) Calcul de la fonction de répartition inverse Gppm de la proportion de pièces en 10 non-conformité que l'on risque d'accepter en ne connaissant que la valeur mesurée Cpk de l'indice de capabilité sur la dimension caractéristique pour l'échantillon de n pièces, selon la formule : 1 1 = 106 [Gppm 1 - LOI. NORMALE. STANDARD 3 * Cpk 1 + z 2(n 1)1) - 9nCpk2 [+ LOI. NORMALE. STANDARD(-3 * Cpk 1 + z 1 1 9nCpk2+ 2(n - 1) ) où a E [0; 1] et za représente le aème percentile de la loi normale standard ; c2) Calcul de l'espérance mathématique de la proportion de pièces en non- 15 conformité en prenant l'intégrale entre 0 et 1 de la fonction de répartition inverse Gppm selon a. - à l'étape c), on prend pour hypothèse que la population des pièces suit une loi normale, et on détermine l'espérance mathématique de la proportion de pièces en non-conformité à partir d'un calcul de la densité de probabilité de la proportion de 20 pièces en non-conformité, où : o On considère que le risque d'avoir X0 ppm de pièces en non-conformité dans la population globale correspond à un indice de capabilité donné par la formule : Cpko = 3 LOI. NORMALE. STANDARD. INVERSE(1 - 2xi°06)) 25 o On détermine la densité de probabilité de la proportion de pièces en non- conformité en se basant sur un percentile p de l'intervalle de confiance calculé à partir de la valeur mesurée Cpk de l'indice de capabilité sur la dimension caractéristique pour l'échantillon de n pièces selon la formule : Cpko - Cpk 1 1 + 9nCpk2 2(n - 1) p = LOI. NORMALE.STANDARD( 3029624 7 - à l'étape c), on prend pour hypothèse que la population des pièces suit une loi normale, et on effectue les calculs suivants : cl) Simulation d'une pluralité de tirages de population réelle à partir d'une moyenne ri et d'un écart-type 6 de la dimension caractéristique mesurée sur 5 l'échantillon de n pièces en prenant pour hypothèse que la variable aléatoire r) suit une loi de Student à n - 1 degrés de liberté, et que la variable (32 aléatoire (n -1)p est distribuée suivant une loi du X2 à n - 1 degrés de liberté, où pt. et a sont respectivement la moyenne réelle et l'écart-type réel de la dimension caractéristique pour la population ; puis 10 c2) Calcul de la proportion de pièces en non-conformité pour chaque tirage simulé, et calcul de la moyenne de la proportion de pièces en non-conformité sur l'ensemble des tirages simulés pour en déduire l'espérance mathématique de la proportion de pièces en non-conformité. - le flux de fabrication de pièces correspond à un lot ayant un nombre de pièces No 15 déterminé, où : o A l'étape a), le nombre de pièces n prélevées pour former l'échantillon est inférieur au nombre No de pièces du lot ; o A l'étape c), on calcule en outre une moyenne de la proportion de pièces en non-conformité dans le lot en pondérant le nombre de pièces de 20 l'échantillon mesurées en conformité par la valeur 0, en pondérant le nombre de pièces de l'échantillon mesurées en non-conformité par la valeur de l'espérance mathématique calculée, et en pondérant le nombre de pièces du lot non-prélevées dans l'échantillon par la valeur de l'espérance mathématique calculée ; 25 o A l'étape d), on compare la moyenne de la proportion de pièces en non- conformité dans le lot à la valeur seuil de proportion de pièces en non-conformité ; o A l'étape e), on accepte ou on refuse le lot de pièces produites en fonction des résultats de comparaison de l'étape d). 30 DESCRIPTION DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative et doit être lue en regard des dessins annexés, sur lesquels : 35 - la figure 1 est un tableau présentant les valeurs seuils en termes d'indice de capabilité Cpk en fonction de la taille de l'échantillon, pour une spécification d'indice de capabilité minimum égal à 1. 3029624 8 - la figure 2 est un tableau présentant les valeurs seuils en termes d'indice de capabilité Cpk en fonction de la taille de l'échantillon, pour une spécification d'indice de capabilité minimum égal à 1.33. - La figure 3 est un graphique illustrant l'évolution des valeurs seuils en termes 5 d'indice de capabilité Cpk en fonction de la taille de l'échantillon, pour une spécification d'indice de capabilité minimum égal à 1. - La figure 4 est un graphique illustrant l'évolution des valeurs seuils en termes d'indice de capabilité Cpk en fonction de la taille de l'échantillon, pour une spécification d'indice de capabilité minimum égal à 1.33. 10 - La figure 5 est un graphique illustrant la fonction de répartition inverse de la proportion de pièces en non-conformité pour un indice de capabilité mesuré de 1.51 sur un échantillon de 50 pièces ; - La figure 6 est un graphique illustrant la densité de probabilité de la proportion de pièces en non-conformité pour un indice de capabilité mesuré de 1.51 sur un 15 échantillon de 50 pièces ; - La figure 7 est un graphique illustrant la densité de probabilité de la proportion de pièces en non-conformité pour un indice de capabilité mesuré de 1.1 sur un échantillon de 50 pièces ; - La figure 8 est un graphique représentant les simulations de plusieurs tirages 20 représentant dans le plan (p ; a) les positions potentielles de la population réelle. DESCRIPTION DETAILLEE DE L'INVENTION Le principe de l'invention se base sur le fait que ce n'est pas parce qu'un échantillon ne respecte pas strictement le critère statistique sévérisé requis pour une 25 dimension caractéristique, que l'ensemble de la population ne pourra pas en moyenne vérifier la spécification imposant de produire une proportion de pièces en non-conformité inférieure à un certain seuil imposé. Reprenons l'exemple plus haut de l'art antérieur où, avec des échantillons de 50 pièces, on impose d'avoir un indice de capabilité mesuré supérieur à 1.22 pour être 30 certain que l'ensemble de la population de pièces satisfasse l'exigence d'un indice de capabilité minimum égal à 1. Des simulations, de type tirage de Monte-Carlo par exemple, des populations réelles que l'on est susceptible d'avoir produites lorsque la mesure d'un échantillon de 50 pièces donne un indice de capabilité mesuré d'une valeur de 1.06 (c'est-à-dire bien inférieur à la valeur minimale de 1.22 tolérée par l'état de l'art) montrent que 35 l'acceptation répétée de populations ayant un indice de capabilité mesuré de 1.06 sur 50 pièces permet de vérifier en moyenne, c'est-à-dire sur de longues périodes, la spécification imposant de produire moins de 2700 ppm de pièces en non-conformité. 3029624 9 Parfois la proportion de pièces en non-conformité dans la population livrée se trouvera supérieure à 2700, mais ces situations sont compensées par les périodes au cours desquelles la proportion de non-conformités est inférieure. Selon cet exemple, il n'est pas dangereux vis-à-vis de la spécification tolérant 2700 ppm de non-conformités d'accepter 5 de manière répétée des lots de pièces sur la base de la mesure d'un indice de capabilité supérieur ou égal à 1.06 sur des échantillons de 50 pièces, puisque les non-conformités éventuelles seront compensées en moyenne. Ainsi, la valeur pour l'indice de capabilité mesuré de 1.22 qui faisait référence dans l'état de l'art antérieur peut par conséquent être abaissée à 1.06 (lorsque le critère minimal de l'indice de capabilité sans allègement de 10 contrôle est égale à 1), ce qui représente un saut quantitatif important. On propose donc de suivre la production des pièces en se fondant sur un calcul d'espérance de la proportion de pièces en non-conformité. On compare ensuite cette espérance de la proportion de pièces en non-conformité à un seuil de pièces en non-conformité acceptables pour piloter la production en conséquence. 15 Ainsi, le suivi d'un flux de fabrication d'une population de pièces basé sur l'analyse d'au moins un indicateur statistique représentatif d'une dimension caractéristique des pièces, peut se dérouler selon les étapes successives suivantes : a) On prélève dans le flux de fabrication de pièces un échantillon comprenant un nombre n de pièces ; 20 b) On mesure la dimension caractéristique de chaque pièce de l'échantillon, et on calcule pour l'échantillon une valeur mesurée de l'indicateur statistique ; c) On calcule une espérance mathématique de la proportion de pièces en non-conformité par rapport à une spécification sur la dimension caractéristique, ledit calcul étant effectué à partir de la valeur mesurée de l'indicateur statistique pour 25 l'échantillon prélevé et du nombre n de pièces de l'échantillon ; d) On compare l'espérance mathématique de la proportion de pièces en non-conformité calculée à une valeur seuil de proportion de pièces en non-conformité ; e) On pilote le flux de fabrication des pièces en fonction des résultats de comparaison de l'étape d). 30 L'espérance mathématique d'une variable aléatoire réelle correspond à la valeur que l'on s'attend à trouver, en moyenne, si l'on répète un grand nombre de fois la même expérience aléatoire. Elle correspond à une moyenne pondérée des valeurs que peut prendre cette variable. On peut également prévoir de suivre la production en ayant deux valeurs 35 comparatives pour affiner le pilotage du flux de production. La première valeur peut ainsi correspondre, comme dans l'art antérieur, à la borne inférieure de l'intervalle de confiance sur l'indicateur statistique à un risque unilatéral a près, qui si elle est respectée garantit 3029624 10 que la production en cours est conforme. La deuxième valeur correspond comme proposé à l'espérance mathématique de la proportion de pièces en non-conformité calculée, qui permet d'agir plus finement sur le pilotage de la production. Ainsi, dans le processus de suivi de la production, on peut à l'étape c), calculer en 5 outre un intervalle de confiance sur l'indicateur statistique à un risque unilatéral a près, où a est compris entre 0 et 1, ledit calcul étant effectué à partir de la valeur mesurée de l'indicateur statistique pour l'échantillon prélevé et du nombre n de pièces de l'échantillon. Dans ce cas, à l'étape d), on compare en outre la borne inférieure de l'intervalle de confiance calculé à une valeur seuil de l'indicateur statistique correspondant à la valeur 10 seuil de proportion de pièces en non-conformité. Le pilotage de la production selon ce mode de réalisation pourrait être le suivant : - Si la borne inférieure de l'intervalle de confiance calculé est supérieure à la valeur seuil de l'indicateur statistique, alors on n'effectue aucune modification au flux de fabrication de pièces (mode de pilotage 1) ; 15 - Si la borne inférieure de l'intervalle de confiance calculé est inférieure à la valeur seuil de l'indicateur statistique et que l'espérance mathématique de la proportion de pièces en non-conformité calculée pour l'échantillon est supérieure à la valeur seuil de proportion de pièces en non-conformité, alors on effectue des actions correctives sur le flux de fabrication de pièces et/ou on augmente le nombre de 20 pièces prélevées lors d'un prélèvement ultérieur d'échantillon (mode de pilotage 2) , - Si l'espérance mathématique de la proportion de pièces en non-conformité calculée pour l'échantillon est inférieure à la valeur seuil de proportion de pièces en non-conformité, alors on stoppe le flux de fabrication de pièces (mode de 25 pilotage 3). Tout comme pour les intervalles de confiance, la valeur de l'espérance mathématique de la proportion de pièces en non-conformité dépend du nombre de pièces contrôlées : plus ce nombre est élevé et plus l'incertitude liée à la représentativité de l'échantillonnage est faible, et donc moins il sera nécessaire de sévériser la valeur 30 minimale de l'indicateur statistique à prouver sur l'échantillon. Le reste de la description est fait en référence à l'indice de capabilité Cpk comme indicateur statistique utilisé dans le suivi de la production du flux de pièces, mais l'enseignement de l'invention n'est pas limité à cet indicateur statistique et pourra être 35 adapté à d'autres indicateurs statistiques classiquement utilisés dans le suivi de productions industrielles. 3029624 11 Le tableau de la figure 1 présente les valeurs seuils en termes d'indice de capabilité Cpk (ce tableau pourrait également être présenté en fonction des proportions seuil de pièces en non-conformité), la valeur CPkintervalle étant la valeur seuil d'indice de capabilité selon l'intervalle de confiance à 5% et la valeur CPkespérance étant la valeur seuil 5 d'indice de capabilité selon l'espérance mathématique de la proportion de pièces en non- conformité. La valeur minimale d'indice de capabilité Cpki',n, pour l'ensemble de la population est fixée ici à 1. On constate sur le tableau de la figure 1 que les valeurs seuil évoluent en fonction du nombre de pièces n dans l'échantillon. 10 Le tableau de la figure 2 est semblable au tableau de la figure 1 mais présente les valeurs seuil correspondant à une valeur minimale d'indice de capabilité Cpki',n, pour l'ensemble de la population est fixée ici à 1.33. Les tableaux des figures 1 et 2 illustrent également différentes zones facilitant le pilotage de du flux de production comme décrites ci-dessus. Ainsi, les zones Z1, Z2, et Z3 15 correspondent par exemple respectivement au mode de pilotage 1, mode de pilotage 2, et mode de pilotage 3. On peut également représenter ces sévérisations sous forme de graphique illustrant, en fonction du nombre de pièces mesurées, la valeur d'indice de capabilité minimale à prouver pour être en zone 1 ou zone 2. 20 Le graphique de la figure 3 correspond à une telle représentation pour une spécification d'indice de capabilité globale Cpkmini = 1. Le graphique de la figure 4 correspond à une telle représentation pour une spécification d'indice de capabilité globale Cpkmini = 1.33. Selon le procédé proposé, on effectue donc un calcul pour déterminer l'espérance 25 mathématique de la proportion de pièces en non-conformité potentiellement engendrée par une situation dans laquelle on a mesuré un indice de capabilité donné sur un échantillon contenant un nombre de pièces donné. Pour un prélèvement donné caractérisé par un nombre de pièces n prélevées dans l'échantillon et un indice de capabilité calculé, on cherche donc à quantifier, en 30 nombre de pièces en non-conformité engendrées (en ppm), l'espérance du risque pris lorsqu'on accepte la population toute entière. L'enjeu est d'établir une sanction intermédiaire (correspondant au mode de pilotage 2 par exemple) lorsque l'espérance du risque engendré sera inférieure à la proportion de pièces en non-conformité tolérées par la spécification. 35 Il existe plusieurs possibilités pour calculer cette espérance mathématique. Nous allons ci-dessous décrire trois méthodes différentes, qui peuvent être mises en oeuvre 3029624 12 indépendamment les unes des autres, ou en parallèle si l'on souhaite avoir une estimation moyennée de l'espérance mathématique. On entend par densité de probabilité la fonction permettant de quantifier la probabilité d'effectuer un tirage dans une plage de valeurs donnée. Plus cette fonction 5 possède une valeur élevée dans une plage de données considérée, et plus les valeurs tirées seront densément réparties dans la plage de données en question après un grand nombre de tirages. Dans la suite de cet exposé nous la noterons f. On entend par fonction de répartition la fonction associant à chaque valeur possible de résultat la probabilité d'effectuer un tirage inférieur à cette valeur. Nous la 10 noterons dans la suite F. Elle est égale à la primitive de f qui s'annule en -09, et se o calcule par conséquent de la façon suivante : F(x0) = L xoo f (t)dt. Enfin, on entend par inverse de la fonction de répartition (ou fonction de répartition inverse) la fonction associant à chaque probabilité p la plus petite valeur de tirage X0 telle que la probabilité de tirer une valeur inférieure à X0 soit égale à p. Nous la noterons G. 15 Elle est définie par la relation : G (t) = inf {x E R : F (x) t} . Pour une variable aléatoire U de loi uniforme sur [0; 1] et F une fonction de répartition donnée, on sait que la variable aléatoire G (U) a pour fonction de répartition F. En corollaire, si X suit une loi de probabilité caractérisée par une densité de probabilité fx, une fonction de répartition Fx et une fonction de répartition inverse G x, alors 20 on a l'égalité : Espérance(X) = fô Gx(t)dt. La première méthode de calcul proposée est basée sur le calcul fonction de répartition inverse Gppm de la proportion de pièces en non-conformité que l'on risque d'accepter en ne connaissant que la valeur mesurée Cpk de l'indice de capabilité sur la 25 dimension caractéristique pour l'échantillon de n pièces. Nous allons décrire cette méthode sur un exemple particulier, qui sera aisément généralisable à une taille d'échantillon quelconque et un indice de capabilité quelconque. Selon cet exemple, on a mesuré sur 50 pièces une valeur Cpk de l'indice de capabilité égale à 1.51. La population globale des pièces du flux de production est supposée 30 normale. La proportion maximale de non-conformités autorisées est fixée à 66 ppm (Cpk = 1.33). Sous l'hypothèse de normalité de la population, on sait convertir un indice de capabilité Cpk en une proportion de pièces en non-conformité engendrées dans l'hypothèse pénalisante d'une loi centrée : 35 Pppm 106[1 - LOI. NORMALE. STANDARD (3 * Cpk) + LOI. NORMALE. STANDARD (-3 * Cpk)]) et réciproquement : 3029624 13 131919 Cpk = 1 -3 LOI. NORMALE. STANDARD. INVERSE(1 2.10 1n 6) La formule de Kushler permet de calculer les bornes inférieures de l'intervalle de confiance sur l'indice de capabilité Cpk au risque unilatéral a près (a E [0;1]). Etant donné un indice de capabilité Cpk mesuré sur n pièces, on court un risque a que l'indice de capabilité Cpk de la population réelle soit inférieur à Cpk [1 + zai 2 + ], où 9nCpk 2(n-1) 5 za représente le aème percentile de la loi normale standard. Par exemple, on a 5% de chance que l'indice de capabilité Cpk réel de la population soit inférieur à Cpk [1 - 1.645 \i9nCpk2 2(n-1) On en déduit que la fonction de répartition inverse de la proportion de pièces en non-conformité que l'on risque d'accepter lorsque l'on tolère la livraison d'un lot duquel on 10 ne connaît qu'un indice de capabilité Cpk mesuré sur n pièces est la fonction qui à a E [0; 1] associe : Gppm = 106 [1 - LOI. NORMALE. STANDARD (3 * Cpk[1+ Z a ,\/9nCpk2 2(n-1) + KI LOI. NORMALE. STANDARD(-3 * Cpk[1+ Z a ,\/9nCpk2 2(n-1)]) Dans le cas de notre exemple où Cpk = 1.51 et n = 50, on peut représenter 15 graphiquement cette fonction de répartition inverse comme illustré à la figure 5. Conformément au corollaire présenté ci-dessus, l'espérance recherchée se calcule comme étant égale à la valeur de l'intégrale entre 0 et 1 de la fonction de répartition inverse de notre variable aléatoire. Pour notre exemple particulier, l'espérance mathématique calculée selon cette 20 première méthode vaut environ 34 ppm. La deuxième méthode de calcul proposée est basée sur un calcul de la densité de probabilité de la proportion de pièces en non-conformité. A partir de la formule de Kushler, il est également possible de représenter l'allure 25 de la densité de probabilité de la variable aléatoire étudiée, en effectuant un raisonnement inverse à celui développé pour la première méthode. Un risque d'avoir X0 ppm de pièces en non-conformité dans la population globale correspond à Cpko = LOI. NORMALE. STANDARD. INVERSE(1 - Le percentile de la loi normale correspondant dans la formule de Kushler vaut : Cpko- Cpk 30 Zo = 1 1 gncpk2+2(n-1) On en déduit qu'on se situe dans le percentile de l'intervalle de confiance p = LOI. NORMALE. STANDARD(zo) ; et il est ensuite possible à partir de nombreux calculs de 3029624 14 valeurs différentes pour p de différencier p par rapport au nombre de ppm engendrés pour obtenir la valeur de la densité de probabilité recherchée. Dans le cas de notre exemple où Cpk = 1.51 et n = 50, on peut représenter graphiquement cette densité de probabilité comme illustré à la figure 6. 5 Pour notre exemple particulier, l'espérance mathématique calculée selon cette deuxième méthode vaut environ 31 ppm. La différence de 3 ppm observée entre les deux méthodes est à mettre sur le compte des imprécisions d'intégration et sur le fait que seule une plage finie de ppm engendrés a été prise en compte dans la seconde méthode, regroupant seulement un 10 peu plus de 99% de l'ensemble des densités de probabilité (ce qui signifie que l'aire sous la courbe de la figure 6, dans le domaine sur lequel elle est tracée, est seulement égale à 0.9904). A titre d'illustration, et pour faciliter la compréhension du phénomène, nous pouvons également représenter la densité de probabilité de notre variable aléatoire dans 15 le cas Cpk = 1.1 et n = 50. Cette représentation est illustrée graphiquement à la figure 7. Dans ce cas précis, si la spécification tolère 2700 ppm de pièces en non-conformité, cette population aurait été rejetée dans l'état antérieur de la technique, car le 95ème percentile se situe à 6843 ppm, bien au-delà du maximum toléré. Or on voit bien que l'espérance reste malgré cela bien inférieure à 2700 ppm (d'ailleurs la médiane et le 20 mode sont encore plus faibles) La troisième méthode de calcul proposée ne repose plus sur la formule de Kushler, mais sur des tirages aléatoires. On considère la mesure d'une moyenne et d'un écart-type a donnés sur un 25 échantillon de 50 pièces issues d'une population supposée normale. En notant respectivement Fi et G la moyenne réelle et l'écart-type réel de la population globale, la variable aléatoire Vil suit une loi de Student à n - 1 degrés de a a2 liberté, et la variable aléatoire (n - 1) a2 est distribuée suivant une loi du x2 à n - 1 degrés de liberté. 30 Ainsi, on peut de simuler plusieurs tirages de population réelle à partir de la seule connaissance que l'on a de l'échantillon. Le graphe de la figure 8 représente ces simulations effectuées sur 10 000 tirages représentant dans le plan 1.1.; G les positions potentielles de la population réelle. Ici la population contrôlée sur 50 pièces a pour moyenne 0 et pour écart-type 0.1167, avec des 35 tolérances à +/-0.53, ce qui représente bien un Cpk de 1.51. Sur chacune des dix mille populations réelles que l'on a tirées, on peut alors calculer la proportion de non-conformités, puis en calculer la moyenne. 3029624 15 Chaque tirage de 10 000 populations donnera un résultat différent, mais ce nombre moyen de la proportion de pièces en non-conformité engendrées varie entre 37ppm et 42ppm. La différence avec les deux précédentes méthodes tient au fait que la formule de 5 Kushler sur lesquelles elles reposent ne fournit qu'une approximation, réputée excellente, des bornes de intervalles de confiance sur l'indice de capabilité. En particulier, on voit que les bornes données par la formule de Kushler sont symétriques, ce qui néglige notamment l'asymétrie de la répartition des écarts-types, qui sont plus dispersés vers les fortes valeurs que vers les petites valeurs. 10 Les résultats des calculs d'espérance mathématique de la proportion de pièces en non-conformité selon ces trois méthodes peuvent se synthétiser selon le tableau suivant : Mesure Borne inférieure de l'intervalle de confiance à Espérance du Espérance Espérance du brute 5% sur la nombre de du nombre nombre de mesure : ppm de ppm engendrés calculée par la deuxième ppm engendrés méthode engendrés calculée par la calculée par la première troisième méthode méthode Cpk 1.51 1.25 N/A N/A N/A ppm 5.9 182 31 34 39 engendrés Sachant que selon l'exemple, la cible en terme de pièces en non-conformité 15 tolérées est fixée à 66 ppm (Cpk = 1.33), l'étude réalisée sur le prélèvement de 50 pièces qui a donné pour mesure un indice de capabilité Cpk nous permet d'affirmer que : - L'espérance du nombre de pièces en non-conformité que l'on risque d'accepter en tolérant la livraison de la population est bien supérieure à ce que le seul indice de capabilité Cpk donné par la mesure brute sans prise en compte des incertitudes 20 laissait présager (environ 35 ppm contre 5.9 ppm pour la mesure brute) ; - Si l'on prend pour référence la borne inférieure de l'intervalle de confiance à 5% sur la mesure de l'indice de capabilité Cpk, on trouve un nombre de pièces en non-conformité engendrés égal à 182 ppm, ce qui est bien supérieur à la spécification de 66ppm. Cela laisse penser que la population possède un risque 25 supérieur à 5% d'être non-conforme à la spécification. - L'espérance du nombre de pièces en non-conformité que l'on risque d'accepter en tolérant la livraison de la population est malgré tout bien inférieure à ce que l'on 3029624 16 envisage si l'on prend pour référence la borne inférieure de l'intervalle de confiance à 5% sur la mesure de l'indice de capabilité Cpk (environ 35 ppm contre 182 ppm). L'espérance du nombre de pièces en non-conformité que l'on risque d'accepter en 5 tolérant la livraison de la population est inférieure à la spécification (environ 35 ppm contre 66 ppm), ce qui permet d'affirmer qu'en tolérant de façon répétée la livraison de telles populations, la spécification serait respectée sur le volume global de livraison. Pour toutes ces raisons, il nous semble raisonnable d'accepter la livraison de ce 10 lot, mais selon un mode de pilotage intermédiaire. Dans l'optique d'une gestion dynamique des prélèvements, il pourra par exemple être demandé pour la période suivante le contrôle d'un nombre plus important de pièces, en vue de mieux caractériser la population livrée et de faire diminuer les risques encourus (avec un niveau accru de prélèvement et à qualité du prélèvement constante, les valeurs 15 renseignées dans les quatre dernières colonnes du tableau ci-dessus diminueront). En conclusion, au lieu de sanctionner une population caractérisée par prélèvement uniquement sur la base du calcul du 95ème percentile le plus défavorable de l'intervalle de confiance sur l'indice de capabilité Cpk, c'est désormais l'ensemble de sa fonction de répartition qui est prise en compte et rentre en jeu dans le calcul d'une espérance 20 permettant de statuer sur son acceptabilité ou non, et ainsi piloter le processus de production. Une application supplémentaire du calcul d'espérance proposé peut également être de permettre la libération d'un lot. 25 Reprenons l'exemple illustré par le graphique de la figure 7, où on a contrôlé 50 pièces et on a trouvé un Cpk égal à 1.1. Le calcul d'espérance exposé plus haut permet d'établir que l'espérance de la proportion de pièces en non-conformité parmi les pièces non contrôlées est de 1906ppm. Si le lot contient 75 pièces (50 contrôlées et 25 non contrôlées) et que l'on n'a 30 détecté aucune non-conformité parmi les pièces contrôlées, alors l'espérance de la proportion totale de pièces non conformes dans le lot est de (50*0+25*1906)/(50+25) = 635ppm. Si la spécification exige d'avoir moins de 700 ppm de non-conformités parmi les pièces livrées (ce qui correspond à un indice de capabilité Cpk de 1.14), on en déduit que 35 le lot pourra être libéré car le risque en jeu, exprimé en termes de non-conformités produites, est inférieur à celui toléré par la spécification. 3029624 17 Ainsi, cela permet d'alléger les contrôles, alors même que le Cpk mesuré (1.1) est inférieur à celui spécifié (1.14). Dans l'état antérieur de l'art, ce lot aurait bien entendu été refusé sans contrôle supplémentaire. Ainsi, dans un flux de fabrication de pièces correspond à un lot ayant un nombre 5 de pièces No déterminé, on peut adapter le procédé de suivi présenté de la manière suivante: - A l'étape a), le nombre de pièces n prélevées pour former l'échantillon est inférieur au nombre No de pièces du lot ; - A l'étape c), on calcule en outre une moyenne de la proportion de pièces en non- 10 conformité dans le lot en pondérant le nombre de pièces de l'échantillon mesurées en conformité par la valeur 0, en pondérant le nombre de pièces de l'échantillon mesurées en non-conformité par la valeur de l'espérance mathématique calculée, et en pondérant le nombre de pièces du lot non-prélevées dans l'échantillon par la valeur de l'espérance mathématique calculée ; 15 - A l'étape d), on compare la moyenne de la proportion de pièces en non-conformité dans le lot à la valeur seuil de proportion de pièces en non-conformité ; - A l'étape e), on accepte ou on refuse le lot de pièces produites en fonction des résultats de comparaison de l'étape d). Ce dernier exemple illustre bien la puissance, l'utilité et les gains avérés pour le 20 producteur d'une caractérisation adéquate et précise des risques associés à une situation d'allègement de contrôle, dans laquelle il convient de prendre en compte les incertitudes liées au fait que seule une partie de la production aura été mesurée. REVENDICATIONS1. Procédé de suivi d'un flux de fabrication d'une population de pièces basé sur l'analyse d'au moins un indicateur statistique représentatif d'une dimension caractéristique des pièces, selon lequel : a) On prélève dans le flux de fabrication de pièces un échantillon comprenant un nombre n de pièces ; b) On mesure la dimension caractéristique de chaque pièce de l'échantillon, et on calcule pour l'échantillon une valeur mesurée de l'indicateur statistique ; c) On calcule une espérance mathématique de la proportion de pièces en non- conformité par rapport à une spécification sur la dimension caractéristique, ledit calcul étant effectué à partir de la valeur mesurée de l'indicateur statistique pour l'échantillon prélevé et du nombre n de pièces de l'échantillon ; d) On compare l'espérance mathématique de la proportion de pièces en non-conformité calculée à une valeur seuil de proportion de pièces en non-conformité ; e) On pilote le flux de fabrication des pièces en fonction des résultats de comparaison de l'étape d). 2. Procédé selon la revendication 1, dans lequel : - à l'étape c), on calcule en outre un intervalle de confiance sur l'indicateur statistique à un risque unilatéral a près, où a est compris entre 0 et 1, ledit calcul étant effectué à partir de la valeur mesurée de l'indicateur statistique pour l'échantillon prélevé et du nombre n de pièces de l'échantillon ; - à l'étape d), on compare en outre la borne inférieure de l'intervalle de confiance calculé à une valeur seuil de l'indicateur statistique correspondant à la valeur seuil de proportion de pièces en non-conformité. 3. Procédé selon la revendication 2, dans lequel à l'étape e) : - Si la borne inférieure de l'intervalle de confiance calculé est supérieure à la valeur seuil de l'indicateur statistique, alors on n'effectue aucune modification au flux de fabrication de pièces ; - Si la borne inférieure de l'intervalle de confiance calculé est inférieure à la valeur seuil de l'indicateur statistique et que l'espérance mathématique de la proportion de pièces en non-conformité calculée pour l'échantillon est supérieure à la valeur seuil de proportion de pièces en non-conformité, alors on effectue des actions correctives sur le flux de fabrication de pièces et/ou on augmente le nombre de pièces prélevées lors d'un prélèvement ultérieur d'échantillon ; Si l'espérance mathématique de la proportion de pièces en non-conformité calculée pour l'échantillon est inférieure à la valeur seuil de proportion de pièces en non-conformité, alors on stoppe le flux de fabrication de pièces. 3029624 19 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel le au moins un indicateur statistique est un indice de capabilité Cpk défini par la formule : Min(TS - t; - TI) Cpk = 3o- où : - p est la moyenne des dimensions caractéristiques mesurées sur les pièces de 5 l'échantillon prélevé ; - a est l'écart-type des dimensions caractéristiques mesurées sur les pièces de l'échantillon prélevé ; - TS est une tolérance supérieure de la dimension caractéristique mesurée ; - Tl est une tolérance inférieure de la dimension caractéristique mesurée. 10 5. Procédé selon la revendication 4, dans lequel à l'étape c), on prend pour hypothèse que la population des pièces suit une loi normale, et on effectue les calculs suivants : cl) Calcul de la fonction de répartition inverse Gppm de la proportion de pièces en non-conformité que l'on risque d'accepter en ne connaissant que la valeur mesurée Cpk de l'indice de capabilité sur la dimension caractéristique pour l'échantillon de 15 n pièces, selon la formule : = 106 [Gppm 1 - LOI. NORMALE. STANDARD 3 * Cpk 1+ za 1 1 LOI. NORMALE. STANDARD(-3 * Cpk 1+ za [+ 9nCpk2+ 2(n - 1)1) 1 1 9nCpk2+ 2(n - 1)1)1 où a E [0; 1] et za représente le aème percentile de la loi normale standard ; c2) Calcul de l'espérance mathématique de la proportion de pièces en non-conformité en prenant l'intégrale entre 0 et 1 de la fonction de répartition inverse Gppm selon a. 20 6. Procédé selon la revendication 4, dans lequel à l'étape c), on prend pour hypothèse que la population des pièces suit une loi normale, et on détermine l'espérance mathématique de la proportion de pièces en non-conformité à partir d'un calcul de la densité de probabilité de la proportion de pièces en non-conformité, où : On considère que le risque d'avoir X0 ppm de pièces en non-conformité dans la 25 population globale correspond à un indique de capabilité donné par la formule : Cpko = -31 LOI.NORMALE.STANDARD.INVERSE (1 - 2xi°06)) - On détermine la densité de probabilité de la proportion de pièces en non- conformité en se basant sur un percentile p de l'intervalle de confiance calculé à partir de la valeur mesurée Cpk de l'indice de capabilité sur la dimension 30 caractéristique pour l'échantillon de n pièces selon la formule : 3029624 20 p = LOI.NORMALE.STANDARD( 1_L 1 9nCik2 2(n - 1) 7. Procédé selon la revendication 4, dans lequel à l'étape c), on prend pour hypothèse que la population des pièces suit une loi normale, et on effectue les calculs suivants : cl) Simulation d'une pluralité de tirages de population réelle à partir d'une moyenne ri et d'un écart-type 6 de la dimension caractéristique mesurée sur l'échantillon de n 5 pièces en prenant pour hypothèse que la variable aléatoire Vil, r) suit une loi (32 de Student à n - 1 degrés de liberté, et que la variable aléatoire (n - 1)p est distribuée suivant une loi du X2 à n - 1 degrés de liberté, où pt. et a sont respectivement la moyenne réelle et l'écart-type réel de la dimension caractéristique pour la population ; puis 10 c2) Calcul de la proportion de pièces en non-conformité pour chaque tirage simulé, et calcul de la moyenne de la proportion de pièces en non-conformité sur l'ensemble des tirages simulés pour en déduire l'espérance mathématique de la proportion de pièces en non-conformité. 8. Procédé selon l'une quelconque des revendications des revendications 1 à 7, dans 15 lequel le flux de fabrication de pièces correspond à un lot ayant un nombre de pièces No déterminé, où : - A l'étape a), le nombre de pièces n prélevées pour former l'échantillon est inférieur au nombre No de pièces du lot ; - A l'étape c), on calcule en outre une moyenne de la proportion de pièces en non- 20 conformité dans le lot en pondérant le nombre de pièces de l'échantillon mesurées en conformité par la valeur 0, en pondérant le nombre de pièces de l'échantillon mesurées en non-conformité par la valeur de l'espérance mathématique calculée, et en pondérant le nombre de pièces du lot non-prélevées dans l'échantillon par la valeur de l'espérance mathématique calculée ; 25 - A l'étape d), on compare la moyenne de la proportion de pièces en non-conformité dans le lot à la valeur seuil de proportion de pièces en non-conformité ; - A l'étape e), on accepte ou on refuse le lot de pièces produites en fonction des résultats de comparaison de l'étape d). Cpko - Cpk