Brevet FR2654849A1 — publié le 19910524

Les procédés industriels mettent en jeu un nombre important de variables qui sont, notamment, des paramètres de composition des matières premières utilisées et/ou des facteurs de réglage caractéristiques de la mise en oeuvre. Généralement, on cherche à privilégier une propriété propre au procédé ou au produit et on tente de trouver la meil- leure combinaison de variables permettant de parvenir au but recherché La situation se complique bien évidemment à mesure que croît le nombre de variables mises en jeu pour peu que croisse également le nombre de propriétés concernées Si on souhaite quadriller l'espace des variables pour évaluer inci- dence sur plusieurs propriétés, on arrive très vite, même avec l'aide de l'informatique, à des temps de calcul rendant cette étude impossible A titre d'information, pour obtenir une pré- cision de 10-3 sur le domaine de variation, et avec l'aide d'un microordinateur dont la vitesse de calcul est de 5/100 de seconde par point calculé, on constate que pour un calcul com- plet mettant en jeu 4 variables, il faut effectuer 1012 calculs élémentaires ( 1000 niveaux à chacune des variables puis combi- naison de tous les niveaux possibles de toutes les variables) soit 1600 ans d'utilisation du microordinateur. L'invention propose une technique de détermination de la relation optimale variables/propriétés, cette technique étant telle qu'on ne sorte pas par construction du domaine de calcul préalablement choisi. L'invention propose également une technique permettant de fixer à l'avance le nombre de calculs élémentaires et, en conséquence, le temps nécessaire pour parvenir à la solution désirée. D'autres avantages de la technique faisant l'objet de l'invention apparaîtront à la lecture de ce qui suit. L'invention a donc pour objet une technique de déter- mination de la relation optimale variable/propriété dans un procédé et/ou une composition mettant en jeu k variables X in- fluant sur p propriétés Y, technique comprenant les opérations suivantes: 2 2654849 à partir d'un plan d'essai et/ou par l'analyse de données, établissement des équations de réponse xi = f (ui) choix d'une valeur Yj calcul de la valeur minimum de j 62 = z (xi Vi*)2 formules dans lesquelles Yj désigne la valeur numérique de la propriété Y de numéro j Yj: la valeur pondérée de Yj (valeur visée) Xj: la valeur pondérée de Xj donnée par le modèle de régression f(ui) vj*: la valeur pondérée codée de Yj* vj: la valeur pondérée codée de Yj ui: la valeur pondérée codée de la variable de valeur numérique X de numéro i (Xi) les pondérations étant obtenues par les équations vj = 2 Yj -Yj Y max Y min et xi xi ui = 2 X max X min Yj désignant la valeur moyenne de la propriété et Xi la valeur moyenne de la variable, Y max et Y min désignant les valeurs réelles extrêmes des propriétés dans le domaine expérimental, X max et X min désignant les valeurs réelles extrêmes des variables dans le domaine expérimental, Yj désignant la valeur réelle de la propriété et Xi la valeur réelle de la variable, cette technique de détermination étant caractérisée en ce qu'elle consiste à créer un hypercube des variables codées constituant le domaine expérimental, à effectuer les calculs de vj pour les valeurs de ui correspondant aux 2 k sommets de l'hypercube dont les cotes sont une fraction des cotes des variables codées au 1 stade puis de 62, le sommet donnant la plus faible valeur de 62 devenant le centre d'un nouvel hyper- cube dont les cotes sont à nouveau une fraction des cotes de l'hypercube précédent, et ainsi de suite jusqu'à obtenir la valeur minimum de 62. Pour des raisons évidentes de simplification, et pour faciliter la compréhension, l'hypercube du domaine expérimental peut avoir pour cotes -1 et + 1, le premier hypercube de cal- cul pouvant avoir pour cotes -1, + 1, le 2 hypercube de cal- 2 2 cul pouvant avoir pour cotes -1, + 1 et ainsi de suite, en pre- 4 4 nant chaque fois pour cotes la moitié des cotes de l'hypercube précédant. Dans la mise en place de la technique de détermination conforme à l'invention, les valeurs successives de ui peuvent être calculées au moyen de l'algorithme (A 1): D * 1 r E n-1 1 E E-1 i =ci+* _ i + *cos E i 1 + (p-2) E i*) 2 Rp Npi pi p_ o: ci désigne les coordonnées du centre de l'hypercube R o se déroule le calcul; D désigne l'amplitude du calcul partiel dans l'ité- 2 R ration de rang R, D désignant la distance de calcul à partir du centre de l'hypercube expérimental, R correspondant à l'ordre des itérations successives et prenant par conséquent les valeurs 1, 2, 3 Rmax; 1 r n 1 -LE + 1 désigne le dimensionnement Nl des hypercubes gigognes dans l'espace R, N désignant le nombre d'hypercubes gigognes, E indiquant qu'on ne retient que la par- tie entière de la fraction qui suit, k désignant le nombre de variables prises en considération, N représentant le pas de calcul dans chaque hypercube R, et variant de 1 à Npk, p dési- 4 2654849 gnant le nombre de points calculés sur chaque arête de l'hypercube; E(E nl n-l os (E i 1 + (p-2) E -i) * p 1 désignant les coordonnées de chaque point exploré dans les hypercubes gigognes, E ayant la signification donnée précédem- ment et i désignant le numéro de la variable, variant de 1 à k. Si on effectue les calculs à partir du centre du domaine expérimental, on retiendra Ci = O pour la première ité- ration (R = 1). Si on effectue les calculs à partir d'un autre point, on retiendra: BS + BI Ci = Xmax Xmin avec Xmax et Xmin pour désigner les valeurs de variables constituant les limites du domaine expérimental et BS et BI les bornes du domaine de calcul. Si le calcul est strictement limité au domaine expéri- mental, D = 1 On pourra de la même manière explorer des solu- tions plus proches du centre (D 1. Il ressort des explications précédentes que, sous une forme simplifiée (cas o on ne crée pas d'hypercubes gigognes et o on ne calcule que les 2 k sommets des différents hypercu- bes) Al peut s'écrire (A 2): 1 r n-i ui = ci + cos LE 2 i rj PR 2 i_ 1 De même si on forme des hypercubes gigognes en ne retenant que les sommets desdits hypercubles ( 2 points par arête), l'algorithme pourra s'écrire (A 3): 1 1 E n-l m n-l ui =ci N L 2 k + 1 * cos * 2 R N 2 k i 2 i-1 Z De même si on veut étudier l'influence d'un très grand nombre de variables en n'explorant que les bornes du domaine expérimental, en procédant à 2 k calculs, on pourra utiliser l'algorithme (A 4): r n-l m Ui = Cos LE i_ 1 r Les symboles utilisés dans A 2, A 3 et A 4 ont les signi- fications données précédemment. Le lecteur aura compris que l'optimisation de la rela- tion variable/propriété peut consister soit dans la mise en évidence d'une ou de la solution exacte (cas o le nombre de propriétés est égal ou inférieur au nombre de variables, soit dans la mise en évidence du meilleur compromis (cas o le nombre de propriétés est supérieur au nombre de variables). La technique objet de l'invention sera illustrée à l'aide des figures. Le problème posé est la détermination du meilleur couple de variables V 1 et V 2 permettant d'obtenir deux proprié- tés ayant respectivement pour valeur P 1 l = 50 et P 2 = 15 (valeurs visées). Les figures 1 à 3 sont données à titre comparatif, en s'appuyant sur des calculs traditionnels. Sur la figure 1, on indique les courbes de niveau de Pl en fonction de V 1 et V 2: on constate qu'il existe un grand nombre de couples V 1 et V 2 donnant Pl = 50. La figure 2 indique les courbes de niveau de P 2 et montre également une pluralité de couples donnant P 2 15 La superposition des figures 1 et 2 (figure 3) permet d'établir les valeurs V 1 R et V 2 R conduisant au résultat cherché. On constate d'une part que l'établissement des diverses courbes de niveau a nécessité un grand nombre de cal- culs élémentaires et, d'autre part, que cette technique est d'autant plus difficile à mettre en oeuvre que le nombre de variables est élevé, la superposition d'un grand nombre de gra- phiques devenant impossible. Les figures 4 à 8 illustrent la technique de l'inven- tion. La figure 4 permet de visualiser les valeurs choisies (P* = 50; P 2 = 15): (point C). La figure 5 illustre le 1 stade du calcul: le carré représentant le domaine expérimental, on calcule, pour quatre points situés à mi-distance entre le centre du domaine et les quatre sommets du carré, les valeurs des propriétés On les compare aux propriétés visées (C) et retient la solution la plus proche de C. Au 2 stade (figure 6), on part de la meilleure solu- tion trouvée à l'étape précédente (sommet du carré la plus proche de C) et renouvelle le calcul avec le pas de calcul pré- cédent divisé par 2 On forme un nouveau carré dont l'un des sommets est plus proche de C qu'à la deuxième étape. Au 30 stade (figure 7), on effectue à nouveau le cal- cul à partir du nouveau sommet retenu et avec un pas de calcul à nouveau divisé par 2. La figure 8 illustre globalement l'algorithme en représentant 4 itérations, sachant que 9 itérations permettent de parvenir à une précision de 10-3. Les figures 9 et 10 sont des représentations de cubes gigognes (figure 9), avec balayage des centres des arêtes et des centres des faces (figure 10). La méthode de détermination de la relation optimale variables/propriétés telle que définie ci-avant, et qui fait appel à la succession d'hypercubes et plus précisément qui met en oeuvre l'algorithme (A 1, A 2, A 3 ou A 4) permet, ainsi qu'il a été précisé, de parvenir à cette relation optimale en fixant, par construction, le domaine du calcul et le nombre de calculs, c'est-à-dire le temps nécessaire pour parvenir au résultat sou- haité A titre purement indicatif, avec le microordinateur men- tionné plus haut, la mise en oeuvre de 4 variables nécessite un temps de calcul de 8 secondes (à comparer aux 1600 ans), tou- jours pour une précision de 10-3 Ce nombre total de calculs élémentaires est donné par la formule générale: Q = R max * N * pk dans laquelle R max représente le nombre d'itérations (par exemple R max = 9 correspond à une précision de 10-3 sur la variation du domaine), N, p et k ayant été définis précédem- ment. Cette méthode est utilisable dans un grand nombre de domaines et d'applications industriels, qu'il s'agisse de la mise en oeuvre de procédés, qu'il s'agisse de la préparation de compositions ou de formulations, qu'il s'agisse d'optimiser les conditions de fabrication ou de traitement de matériaux divers, d'étudier l'influence d'impuretés, par exemple en analyse spec- troscopique ou encore dans la transformation de produits divers, par exemple de matière plastique en particulier en optimisant les paramètres de réglage d'une extrudeuse compte tenu des propriétés souhaitées pour l'extrudat. Les exemples suivants sont donnés à titre purement indicatif. EXEMPLE 1 L'objectif est de réaliser un gris dont les coordon- nées trichromatiques ont pour valeur: L = 63: L = brillance variant de O (noir) à 100 (blanc) a = 2 05: a = axe vert-rouge (; +) b = 3 83: b = axe bleu-jaune (; +) Une telle teinte peut être obtenue en mélangeant 2 pigments: un noir et un bleu. Pour atteindre rapidement l'objectif, on met en oeuvre un plan 22 avec les niveaux suivants: pigment noir: niveau mini (-1) = 50 ppm niveau maxi (+ 1) = 500 ppm pigment bleu: niveau mini (-1) = 200 ppm niveau maxi (+ 1) = 2000 ppm On mesure pour chaque condition d'essai les valeurs de 3 propriétés concernées: Essai Noir Bleu L a b (N) (B) 1 66 01 1 27 0 87 2 + 62 88 3 25 7 01 3 + 53 27 0 55 0 41 4 + + 51 24 -1 44 -2 96 Le calcul classique à partir du plan d'essai conduit aux 3 modèles cidessous dont l'écriture en minuscule rappelle qu'il s'agit des variables explicatives codées: L = 58 35 6 095 N 1 290 b + 0 275 nb a = 1 63 + 0 633 N 0 718 b + 0 273 nb b = 2 81 + 1 128 N 2 173 b + 0 898 nb On voit que chacune des composantes trichromatiques est influencée par les 2 pigments utilisés Etant donné qu'on met en jeu 3 propriétés pour 2 variables, on sait qu'il n'existe pas de solution exacte (cela ressort de la superposition des 3 isoréponses visées, tracées dans le plan des variables figure 11: sur cette figure, le point fléché désigne l'objectif). Le calcul d'optimisation donne en 2 secondes la meilleure solution avec seulement 36 calculs (R max = 9, N = 1, p = 2 et k = 2): Valeur visée Valeur de compromis L 63 62 895 a 2 05 2 110 b 3 83 3 684 Illa 9 2654849 On constate, au niveau du calcul, que la moins mauvaise solution est acceptable puisqu'elle conduit à un écart de teinte A E = O 19, alors que l'oeil n'a pas cette sensibilité de perception. Les niveaux des variables calculées pour la technique de l'invention, respectivement 69 ppm pour le noir et 1257 ppm pour le bleu, ont été retenus pour réaliser la teinte, jugée acceptable par le coloriste. EXEMPLE 2 L'objectif à atteindre est défini par les coordonnées trichromatiques: L = 50 44 a = 27 24 b = 29 21 Trois pigments sont utilisés pour réaliser la teinte souhaitée (marron soutenu), avec les taux mini et maxi suivants: Taux mini Taux maxi Noir 10 ppm 30 ppm Rouge 200 ppm 400 ppm Jaune 1 500 ppm 1000 ppm Etant donné le coût des essais qui impliquent un mélange, une granulation puis un malaxage sur cylindres, on a décidé de les réduire le plus possible, et on a donc mis en oeuvre un plan fractionnaire 2 ( 3-1), dont les résultats sont: Il Essai N R J L a b 1 + + + 47 27 24 31 25 62 2 + 52 49 21 68 28 30 3 + 56 90 31 21 40 10 4 + 52 45 34 96 29 24 La modélisation conduit aux équations codées: L = 52 28 2 40 N 2 42 r 0 19 j a = 20 04 5 05 N + 1 59 r 0 28 j b = 30 82 3 85 N 3 39 r + 2 04 j Pour le calcul d'optimisation, on a retenu p= 3 et N= 3 en vue d'obtenir une meilleure précision de calcul, et au terme de 729 calculs élémentaires on aboutit aux résultats suivants: Valeurs visées Valeurs calculées L 50 44 50 435 a 27 24 27 238 b 29 21 29 214 Les quantités de pigments à utiliser pour atteindre l'objectif, sont d'après le calcul d'optimisation, en valeurs brutes: Noir: 23 ppm Rouge: 345 ppm Jaune: 871 ppm La réalisation effective a donné un résultat très proche de ce qui était visé, soit: L = 50 45 a = 27 18 b = 28 81 -1 On poursuit le calcul en ne faisant évoluer que l'axe rouge seul, les deux autres coordonnées trichromatiques restant constantes. Valeurs visées Valeurs trouvées N R J L a b L a b 50.44 25 5 29 21 50 539 25 515 29 029 25 316 848 26.0 " 50 440 26 000 29 210 25 326 879 26.5 50 440 26 499 29 210 24 334 875 27.0 50 440 27 001 29 211 23 341 871 27.5 50 440 27 499 29 211 22 349 861 28.0 50 440 27 998 29 209 21 356 863 28.5 50 514 28 503 29 070 21 361 859 29.0 50 460 28 998 29 170 20 370 847 L'examen du parvenir à l'objectif tableau qui est qui précède montre que pour ici d'augmenter la seule valeur de a de 26 à 28 5, il faut simultanément: diminuer le noir de 25 à 21 ppm augmenter le rouge de 326 à 361 ppm diminuer le jaune de 879 à 859 ppm (si l'on désire demeurer à l'intérieur du domaine expérimen- tal) Au-delà des limites définies par 26 dans le domaine expérimental. Il faut noter qu'une telle teinte est difficile à réaliser, puisque des variations minimes de pigments font bouger cette teinte au délà de ce que décèle l'oeil: celui-ci perçoit un écart de 0 5 sur A a, ce qui correspond aux varia- tions suivantes: Noir: 1 ppm environ ( 4 %) Rouge: 7 ppm environ ( 2 %) Jaune: 4 ppm environ ( 0 5 %) il Le tracé des courbes isoréponses (figure 12) pour l'objectif initial, dans le plan Noir Rouge pour jaune = 870 ppm, illustre bien à la fois: la sensibilité de la teinte aux faibles variations de concentrations la nécessité de calculer avec une extrême précision puisque l'angle entre les plans b et L est faible. EXEMPLE 3 Dans le but d'optimiser 4 propriétés d'un tissu enduit, on a pratiqué un plan d'expérience mettant en jeu 3 variables qui sont respectivement: A = pourcentage d'additif de traitement T = température de traitement H = humidité du tissu à traiter Les propriétés à optimiser sont les suivantes ADH = adhérence du revêtement DCH = résistance du tissu, sens chaîne DTR = résistance du tissu, sens trame RF = résistance au frottement (note sur 20) Pour simplifier, les valeurs réelles des variables sont remplacés par leurs niveaux codés (-1, + 1). Plan d'essai et résultats obtenus VARIABLES PROPRIETES A T H ADH DCH DTR RF + 1 -1 -1 5 3 1 00 1 08 14 -1 -1 -1 5 6 1 15 1 40 4 + 1 + 1 -1 7 3 0 90 1 18 18 -1 + 1 -1 6 0 1 18 1 40 12 + 1 -1 + 1 8 2 0 95 1 30 12 -1 -1 + 1 5 7 1 20 1 55 6 + 1 + 1 + 1 8 5 0 85 1 30 16 + 1 + 1 6 1 1 35 1 50 10 On notera, à la lecture du tableau ci-dessus, que la valeur maximale de chaque propriété correspond chaque fois à une combinaison différente des variables. On va rechercher progressivement des solutions de plus en plus performantes, pour l'ensemble des 4 propriétés mesurées. Les modèles retenus comportent à la fois les actions principales des variables et leurs interactions deux à deux. ADH = 6 59 + 0 74 a + 0 39 t + 0 54 h + 0 19 at + 0 49 ah 0 21 th DCH = 1 07 0 15 a 0 003 t + 0 016 h 0 047 at 0 041 ah + 0 016 th DTR = 1 34 0 12 a + 0 006 t + 0 074 h + 0 018 at + 0 012 ah 0 018 th RF = 11 5 + 3 5 a + 2 5 t 0 5 h 0 5 at 0 5 ah 0 5 th Le tableau suivant donne les résultats des différents calculs pour les objectifs successifs. Afin de vérifier la cohérence des données, on a également calculé l'optimisation pour les valeurs moyennes, identiques aux constantes des équations centrées qui précèdent. On retrouve bien sur dans ce cas que le niveau de toutes les variables est égal à 0. VALEURS VISEES VALEURS TROUVEES VARIABLES ADH DCH DTR RF ADH DCH DTR RF A T H 0 9 1 0 6 5 5 1 04 1 14 11 8 0 63 -1 -1 5.5 0 95 1 05 8 5 5 1 04 1 14 11 8 0 63 -1 -1 6 1 0 1 1 10 5 8 1 02 1 15 12 1 0 72 -1 -0 75 6.5 1 05 1 15 12 6 3 1 00 1 18 12 5 0 84 -1 -0 33 7 1 1 1 2 14 6 6 1 01 1 26 14 3 0 24 0 53 -0 75 7.5 1 15 1 25 15 7 1 1 04 1 33 14 4 0 13 + 1 O 8 1 2 1 3 16 7 4 1 06 1 37 13 6 0 14 + 1 0 72 8.5 1 25 1 35 18 7 5 1 06 1 38 13 4 0 18 + 1 1 9 1 3 1 4 20 7 6 1 05 1 38 13 5 0 21 + 1 1 VALEURS MOYENNES 6.581 1 071 1 34111 516 581 07 1 1 34111 51 O I I O I 16.581 lo 7 1 l 41 1 810 Il 4 l jo j O O On peut tirer du tableau qui précède les conclusions pratiques suivantes lorsque l'on cherche à améliorer simulta- nément l'ensemble des 4 propriétés: il convient à la fois pour les variables: de mettre peu d'additif d'augmenter la température de traitement de travailler avec le taux maximal d'humidité pour les propriétés, on notera: que la résistance au déchirement dans le sens de la chaîne demeure pratiquement constante 2654849 que s'il est possible d'augmenter simultanément la résistance au déchirement dans le sens de la trame et l'adhérence, propriétés pour lesquelles on peut s'approcher des valeurs maximales expérimentales, il n'en va pas de même pour la résistance au froissement: on est limité à une valeur de 13.5 du fait du compromis souhaité. 16 2654849 R E V E N D I C A T I O N 8 1./ Technique de détermination de la relation optimale variables/propriétés dans un procédé et/ou une composition mettant en jeu k variables X influant sur p propriétés Y, technique comprenant les opérations suivantes: à partir d'un plan d'essai et/ou par l'analyse de données, établissement des équations de réponse: Yi = f (ui) choix d'une valeur Yj calcul de la valeur minimum de j 62 = z (xi vi*)2 formules dans lesquelles Yj désigne la valeur numérique de la propriété Y de numéro j Yj*: la valeur pondérée de Yj (valeur visée) Xj: la valeur pondérée de Xj donnée par le modèle de régression f(ui) vj: la valeur pondérée codée de Yj* vj: la valeur pondérée codée de Xj ui: la valeur pondérée codée de la variable de valeur numérique X de numéro i (Xi) les pondérations étant obtenues par les équations vj = 2 Yj i Y max Y min et Ui2 Xi -Xîi ni = 2 X max X min Yj désignant la valeur moyenne de la propriété et Xi la valeur moyenne de la variable, Y max et Y min désignant les valeurs réelles extrêmes des propriétés dans le domaine expérimental, X max et X min désignant les valeurs réelles extrêmes des variables dans le domaine expérimental, Yj désignant la valeur 17 2654849 réelle de la propriété et Xi la valeur réelle de la variable, cette, technique de détermination étant caractérisée en ce qu'elle consiste à créer un hypercube des variables codées constituant le domaine expérimental, à effectuer les calculs de vj pour les valeurs de ui correspondant aux 2 k sommets de l'hypercube dont les cotes sont une fraction des cotes des variables codées au 1 stade puis de 62, le sommet donnant la plus faible valeur de 62 devenant le centre d'un nouvel hyper- cube dont les cotes sont à nouveau une fraction des cotes de l'hypercube précédent, et ainsi de suite jusqu'à obtenir la valeur minimum de 62. 2./ Technique selon la revendication 1, caractérisée en ce que l'hypercube du domaine expérimental est coté -1, + 1, le premier hypercube de calcul est noté -1, + 1, le deuxième 2 2 -1, + 1 et ainsi de suite. 4 4 3./ Technique selon l'une quelconque des revendica- tions 1 et 2, caractérisée en ce que les valeurs successives de ui sont déterminées au moyen de l'algorithme A 1: i = Ci+DT n * I cos SE 1 + (p-2) E * 2 RNPk pi 1 pi Pl o: ci désigne les coordonnées du centre de l'hypercube R o se déroule le calcul; D désigne l'amplitude du calcul partiel dans l'ité- 2 R ration de rang R, D désignant la distance de calcul à partir du centre de l'hypercube expérimental, R correspondant à l'ordre des itérations successives et prenant par conséquent les valeurs 1, 2, 3 Rmax; 1 n-1 E + 1 désigne le dimensionnement lE 7 pk 18 2654849 des hypercubes gigognes dans l'espace R, N désignant le nombre d'hypercubes gigognes, E indiquant qu'on ne retient que la par- tie entière de la fraction qui suit, k désignant le nombre de variables prises en considération, N représentant le pas de calcul dans chaque hypercube R, et variant de 1 à Npk, p dési- gnant le nombre de points calculés sur chaque arête de l'hypercube; n-1 nl X cos (E + (p-2) E) * - pi-l pi p-l désignant les coordonnées de chaque point exploré dans les hypercubes gigognes, E ayant la signification donnée précédem- ment et i désignant le numéro de la variable, variant de 1 à k. 4./ Technique selon l'une quelconque des revendica- tions 1 à 3, caractérisée en ce que les valeurs successives de ui sont déterminées au moyen de l'algorithme A 2 1 n-l - ui = ci + -R cos E 2 i-1 R_ 2 il 1 * 5./ Technique selon l'une quelconque des revendications 1 à 3, caractérisée en ce que les valeurs successives de ui sont déterminées au moyen de l'algorithme A 3 1 lE n-1 E n-1 l u = Ci À + 1 *+ * cos E 2 r * 2 R N 2 k 2 i-1 6./ Technique selon l'une quelconque des revendica- tions 1 à 3, caractérisée en ce que les valeurs successives de ui sont déterminées au moyen de l'algorithme A 4: n-l ui = cos LE i-i * 2 i_ 7./ Utilisation de la technique selon l'une quelconque des revendications 1 à 6 dans la détermination de la relation 19 2654849 optimale variables/propriétés dans des procédés, formulations, compositions, traitements, analyses, transformation de produits industriels. 8./ L'algorithme utilisé dans les techniques décrites dans l'une quelconque des revendications 3 à 6.