Regulation method for a chemical medium

Abstract: Plant permitting a changing chemical medium to be controlled and regulated to arrive at an optimum value of a predetermined function of physical or chemical parameters representing the medium, the function Y being related to the values of the parameters by a relationship of the form: <IMAGE> the said plant comprising: - sensors (5, 6) intended to supply the instantaneous value of each of the said parameters, - a scanning unit (26) allowing the value of the group of parameters to be recorded at uniform intervals, - a calculating unit (10) making it possible to determine the modification to be made to the said parameters to give an optimum value to the function Y, and - means for adjusting the said parameters controlled by the said calculating unit (10) to bring the value of the said parameters to the optimum. <IMAGE>

Description: Procédé et installation permettant dune Dart. de connaitre d avance à un instant donné le résultat auauel conduit nécessairement un milieu chimique évolutif donné et d'autre Dart. d'assurer la régulation de ce milieu oour Darvenir à un résultat déterminé-fixé à l'avance. L'invention a pour objet un procédé permettant, d'une part, de connaitre d'avance à un instant donné le résultat auquel conduit nécessairement un milieu chimique évolutif donné et, d'autre part, d'assurer la régulation de ce milieu pour parvenir à un résultat déterminé fixé à 1 avance. On rencontre dans l'industrie de nombreux milieux chimiques qui sont évolutifs de par la variation des paramètres qui les caractérisent, à savoir concentration en les divers constituants, pH, température et autres. A titre d'exemples de tels milieux. on peut citer - les bains acides de conversion chimique, - les bains alcalins ou acides de dégraissage des métaux, - l'eau contenue dans des chaudières et installations de chauffage, et autres. L'utilisation de ces milieux conduit à des résultats souhaités dans certains cas et non souhaités dans d'autres et qui peuvent se traduire dans le cas des exemples susdits respectivement - par des propriétés plus ou moins favorables des revetements obtenues par conversion chimique, - par un état de propreté plus ou moins poussé des surfaces dégraissées, par - par des états d'entartrage variables des chau- dières ou installations de chauffage. Or, il importe à l'industriel de savoir à tout instant si le résultat obtenu à l'instant en question par l'utilisation du milieu correspondant est bien celui qui est attendu. Mais la réponse est souvent longue et difficultueuse à obtenir Ainsi, par exemple dans le cas d un bain de conversion chimique. - il faut disposer de 5 ou 6 semaines pour savoir si la résistance au brouillard salin du revêtement obtenu est supérieure à 1000 heures ou non, - il faut disposer d un matériel coûteux et encombrant (microscope électronique à balayage) pour connaitre la structure cristalline du revêtement obtenu. Et il est bien certain qu'une fois ces tests effectués. c'est-à-dire par exemple 6 semaines plus tard, le bain a pu évoluer de manière importante, sa composition pouvant alors etre très différente de celle qu'il avait au moment où avait été obtenu le revêtement testé. Il sten- suit que la modification qu'on aurait éventuellement du apporter au bain à 1 instant en question pour faire en sorte que le résultat obtenu soit celui souhaité, se trouve en fait dépassée au moment où le résultat est connu à l'issue des tests en question Il existe déjà des procédés de régulation de milieux chimiques évolutifs. Mais, dans ces procédes, il est d usage de mesurer de façon répétitive au moins l'un des paramètres d'un tel milieu, de comparer chaque fois la valeur mesurée à une valeur de consigne et d agir sur le paramètre en question de manière à rapprocher la valeur mesurée suivante de la valeur de consigne. par exemple en ramenant, si le paramètre mesuré est un constituant du milieu évolutif, la concentration en ce constituant à la valeur de consigne par apport de produits de régénération appropries Et il est bien certain que. dans un milieu chimique évolutif du genre en question, ces opérations ne peuvent garantir la constance des résultats que 1 on souhaite obtenir. Le problème que la Société Demanderesse s'est donc proposé de résoudre est celui de mettre au point un procédé propre à pérmettre à l'utilisateur - d'une part, de connaître, à partir de la mesure de divers paramètres évolutifs d'un milieu chimique évolutif, les particularités du résultat auquel conduit le milieu en question à l'instant choisi. - d'autre part, au cas où la valeur trouvée par cette mesure montre que le résultat en question ne correspond pas à celui souhaité, de déterminer tout d'abord un ensemble de modifications des variables ou paramètres du milieu. propre à faire en sorte que le résultat auquel conduira le milieu ainsi modifié corresponde à celui souhaité, puis apporter cet ensemble de modifications au milieu chimique évolutif en question, par action sur les paramètres concernés. Et elle a eu le mérite d'élaborer un procédé permettant d'atteindre ce but, ce procédé étant caractérisé par le fait - que, dans le cas d'un milieu chimique évolutif donné et pour au moins l'un des résultats auxquels conduit ce milieu. on établit une équation dite prévisionnelle du type Y = b0 + b1x1 + b2x2 ... + bkxk + ... + b12x1x2 + ... I 2 k n n 12 i 2 + bk-1kxk-1xk + ... + bkkxkê +...+ bnnxnê + nn n dans laquelle - Y représente le résultat auquel conduit le milieu. - X1, X2 ... Xk, ... Xn représentent les variables par paramètres du milieu chimique évolutif dont dépend le résultat Y et - b0, b1 ...bk,..., bnn sont des coefficients par rapport auxquels s'exprime la linéarité de l'équation. la valeur de ces coefficients étant déterminée par des calculs reposant sur des méthodes de régression en ellesmêmes connues, - que l'on mesure de manière répétitive chacun desdits variables ou paramètres par autant de sondes ou capteurs et - que. d'une part, à un instant t donné. c'est-àdire à l'issue d'une des mesures répétitives, on calcule la valeur de chacune des équations Y en utilisant les valeurs mesurées pour les différents variables à l'instant t, - que. d'autre part, éventuellement et au cas où la valeur trouvée pour l'une ou l'autre des équations Y ne correspond pas au résultat souhaité, on détermine par le calcul un ensemble de variations ou modifications à impartir aux variables ou paramètres, cet ensemble de variations donnant pour Y la valeur souhaitée et constituant, parmi toutes les solutions possibles de l'équation prévisionnelle Y, celle qui est la plus économique et la plus facile à réaliser, et - que, d'autre part encore, on agit individuellement sur autant d'organes de réglage que de paramètres, fournissant à chacun de ces organes de réglage un signal propre à amener la variable ou paramètre qui lui est associé à la valeur correspondant à la susdite solution de l'équation prévisionnelle Y la plus économique et la plus facile à obtenir Dans un mode de réalisation avantageux, on effectue les calculs susdits à l'aide d'une calculatrice numérique que l'on munit de bases de données appropriées. L'invention vise encore d'autres dispositions qui s'utilisent de préférence en meme temps et dont il sera plus explicitement question ci-après. Et elle pourra. de toute façon, étre bien comprise à l'aide du complément de description, du dessin y relatif et des exemples. lesdits complément de description et exemples concernant des modes de réalisation avantageux. La figure I de ces dessins est un diagramme relatif à une équation Y et la figure 2 est une représentation schématique d'une installation propre à la mise en oeuvre du procédé selon l'invention. Les figures 3a et 3b prises dans leur ensemble, d'une part, et 4a. 4b et 4c prises dans leur ensemble, d'autre part, représentent chacune un logigramme relatif à la mise en oeuvre du procédé conforme à l'invention. La figure 5 est un schéma synoptique d'un exemple d'installation permettant de mettre en oeuvre l'invention. Se proposant, par conséquent, de se donner les moyens permettant, d'une part, de connaitre à l'avance, à un instant donné, le résultat auquel conduit nécessairement un milieu chimique évolutif et, d'autre part, d'assurer la régulation de ce milieu pour parvenir à un résultat déterminé fixe à l'avance, on s'y prend comme suit ou de façon équivalente. Pour un milieu chimique évolutif donné, il peut être nécessaire de connaitre à tout instant le résultat ou propriété Y1 auquel il conduit, ce résultat ou propriété étant la fonction d'un certain nombre de variables X1, X2 ..... Xi, du milieu et dont dépend le résultat, Y1 et les variables X1, X2, ... Xi étant liés- par l'équation prévisionnelle Y1 = b0 + b1x1 +...+ bixi +...+ bijxixj + bi,ixiê . dans laquelle b0 b1 ... bu ... sont des coefficients, au départ inconnus et déterminés de la manière qui sera indiquée plus loin, la valeur Y1 effective traduisant à un moment donné le résultat auquel conduit le milieu pris dans sa composition correspondant à ce meme moment. En admettant que le résultat en question doit être par exemple égal ou supérieur à la valeur N, il faut que. pour que le milieu à l'instant envisagé puisse être consi déré comme convenable, qu'à cet instant la mesure des variables X1, X2 ... X. fournisse pour Y1 une valeur Y1 # N la mesure des variables X1, X2 ... X. pouvant etre effectuée par autant de capteurs ou sondes plongeant dans le milieu et affectés respectivement à chacune des variables x x .. .x.. 1' 2 i Si Y1 ainsi mesurée est 5 N, le milieu conduit au résultat souhaite ; sinon il faut apporter des modifications au milieu. Mais avant toute chose, il faut bien entendu établir l'équation Y et, par conséquent, - selectionner les variables X. que l'on considère comme etant caracteristiques du milieu (il peut s'agir de la concentration en certains constituants, de la température, du pH et autres) - déterminer les valeurs des coefficients b0 ... b Le choix des variables est déterminant pour l'ob- tention d'une équation prévisionnelle de bonne qualité. l'oubli d'un paramètre influent conduisant à l'établisse- ment d'une équation erronée. Pour chaque problème posé, le choix des paramètres résulte des connaissances antérieures de l'homme de l'art. La détermination des coefficients b nécessite le recours à des calculs reposant sur des méthodes de régression classiques (sur lesquelles il ne sera pas insisté ici mais dont le principe peut; être trouvé dans I ouvrage de L LEBART. A. MORINEAU, J.P. FENELON, "Traitement des données statistiques, méthodes et programmes", DUNOD) en utilisant les données fournies par un certain nombre d'expériences. Du point de vue pratique, on procède pour ce faire de la manière suivante - on choisit la ou les propriétés à mesurer, - on choisit les paramètres expérimentaux ou variables qui vont décrire l'équation, - on choisit pour chaque paramètre un domaine de variations qui est fonction du problème posé et de la connaissance antérieure de l'expérimentation, - on choisit un modèle d'équation simple. le premier modele choisi pouvant être du type linéaire Y1 = b0 + b1X1 + ....+ bi Xi +....+ bnXn ou linéaire avec interactions Y1 = b0 + b1x1 +...+ bixi + bjxj +..+ bnxn +..+ bijxixj+.. - on réalise le nombre d'expériences défini par le choix du modèle et on relève pour chacune d'elles la valeur de la réponse Y obtenue, - on calcule les coefficients ..bi... du modèle. Le modèle est alors testé. Ce modèle étant choisi a priori", il est en effet indispensable de le tester pour l'accepter et l'utiliser, ou bien au contraire le refuser. Pour cela, on réalise quelques expériences supplémentaires pour lesquelles les valeurs de chacun des paramètres sont fixées simultanément à égale distance des extrémités de leur domaine de variations et la réponse obtenue est comparée à la réponse calculée par le modèle. Si les deux résultats sont voisins, le modèle est accepté. Dans le cas contraire, on complique le modèle par adjonction de termes supplémentaires, par exemple des termes quadratiques. L'équation prévisionnelle représentative du nouveau modèle s'écrit alors Y1 b b0 + b1 X1 +. + biXi + bjXj + ij i j +biixiê + bjjxjê + On réalise le supplément d'expériences nécessaires au calcul des coefficients des termes quadratiques et le modèle est de nouveau testé. Ainsi, de proche en proche, par une démarche séquentielle, on affine le modèle jusqu'à ce que la prévision soit de bonne qualité. L'équation ainsi établie est représentée par une hypersurface dans un espace de n+1 dimensions pour n facteurs. Pour matérialiser l'incidence des divers facteurs, on projette l'hypersurface successivement sur des sousespaces à trois dimensions dont deux sont constituées par deux variables que l'on désire étudier et la troisième par la réponse Y, les variables restantes étant fixées à des valeurs constantes ; ceci se traduit par des courbes de niveau dans le plan des deux variables, la réponse étant matérialisée par les valeurs des courbes de niveau ; ainsi on peut envisager de projeter l'hypersurface représentative de l'équation Y1 sur le plan des paramètres X3 et X4; à titre d'exemple, on peut imaginer qu'Y1 représente la résistance au brouillard salin d'une surface métallique qui a été soumise au préalable à un traitement de conversion chimique et à un revêtement par peinture, X3 représentant par exemple la teneur du milieu évolutif, c'est-àdire du bain de conversion, en ions Zn et X4 la teneur du même bain en ions PO43-, les variables restantes X1 et X2 étant fixées à des valeurs constantes. Dans ce cas, l'hypersurface se traduit dans le plan defini par X3, X4 par un certain nombre de courbes dites "courbes isoréponses" ou courbes de niveau et dési gnees par C1, C2 C , chacune de ces courbes correspon n dant à une valeur donnée de la résistance au brouillard salin et chaque point de chacune des courbes représentant un couple de valeurs de X3 et X4 qui conduit à une surface métallique ayant la resistance au brouillard salin correspondant à la courbe. A la figure f, an a représenté un diagramme (projection de l'hypersurfaces sur le plan X3, X4) montrant les courbes isorèponses suivantes courbe C1 : resistance au brouillard salin de 80 h courbe C2 . résistance au brouillard salin de 100 h courbe C3 : resistance au brouillard salin de 400 h courbe C4 : resistance au brouillard salin de 500 h courbe C5 : resistance au brouillard salin de 800 h courbe C6 : resistance au brouillard salin de 900 h courbe C7 : résistance au brouillard salin de 950 h courbe C8 : résistance au brouillard salin de 1000 h. Si on admet que, pour être satisfaisante. la résistance au brouillard salin des surfaces traitées doit être de 950 à 1000 h, tous les couples de valeurs de X3 et Xq correspondant à des points situés sur l'une des courbes C7 et C8 ou dans la surface comprise entre ces courbes conduisent au résultat recherché. En d autres termes si, à un moment donné, la mesure de X3 et X4 conduit à un point P1 répondant à cette condition, la résistance au brouillard salin de la pièce traitée sera satisfaisante ; si, par contre, elle conduit à un point P2 ne répondant pas à cette condition, c'est-àdire situé en dehors du domaine défini par les courbes C7 et C8, le résultat ne sera pas satisfaisant et il s'agira d'apporter au milieu évolutif, c'est-à-dire au bain de conversion, une modification des concentrations X3 et X4 qui ramène à l'intérieur dudit domaine le point situé en dehOS de celui-ci. Et la modification retenue parmi toutes les modifications possibles conduisant à ce résultat devra etre celle qui soit la plus facile et la plus économique à réaliser. Cette détermination est réalisée à chaque instant en effectuant les calculs nécessaires pour déterminer pour chacune des différentes variables concernées et non pas seulement X3 et X4, de chacune des équations prévisionnel les Y1, Y2 .. Y ... la modification la plus simple et la n plus économique parmi toutes les modifications possibles conduisant à un résultat satisfaisant. Par exemple, on peu-t effectuer ces calculs à l'aide d'une calculatrice munie des bases de données né nécessaires, notamment les données relatives aux équations prévisionnelles Y1, Y2 an... . n Pour pouvoir effectuer l'ensemble des calculs, il convient de disposer à l'instant donné pour chaque équation prévisionnelle de l'ensemble des valeurs des paramètres correspondants, ce qui peut être réalisé dans le cadre des mesures répétitives dont il a été question plus haut par des capteurs ou sondes plongeant dans le milieu chimique évolutif considéré. La modification est alors réalisée par l'envoi à un ensemble de dispositifs aptes à agir sur les divers paramètres, de signaux propres à déclencher un ensemble de réponses de la part de ces dispositifs, réponses qui se traduisent au niveau de chaque paramètre par un changement précis, l'ensemble de ces changements conduisant à la susdite modification, sélectionnée et déterminée par les cal culs Pour la commodité de l'exposé, on admet dans ce qui suit et en particulier dans les exemples que le milieu chimique évolutif pris en considération est une séquence de traitements de phosphatationde bain de conversion chimique au zinc, étant bien entendu que les principes qui vont être développés en rapport avec cet exemple particulier de milieu chimique évolutif se transposent aisement à tout autre milieu chimique évolutif en tenant compte des considérations plus générales qui précédent. Il est connu que les bains de conversion en question comprennent des ions nickel. zinc, phosphate et sont accélérés par des nitrites et/ou des chlorates. Par ailleurs, la séquence de traitements peut comprendre de façon classique, par exemple, un traitement de conditionnement préalable de la surface et un traitement subséquent consistant en un rinçage passivant ; il va de soi que la séquence peut comprendre d'autres traitements, mais ceux-ci ne sont pas pris en considération ici pour la commodité de la description. La figure 2 montre schématiquement une installation propre à mettre en oeuvre le procédé conforme à l'in Invention Cette installation comprend - trois enceintes représentées schématiquement en 1, 2 et 3 et renfermant respectivement les bains de conditionnement, de conversion proprement dite et de rinçage passivant. - six sondes ou capteurs représentés schematique- ment en s à 9, constituées par des pH-mètres, des électro des sélectives, des thermomètres et autres et plongeant respectivement, pour ce qui est de la sonde 4, dans l'enceinte 1, pour ce qui est des sondes 5, 6 et 7, dans l'enceinte 2 et, pour ce qui est des sondes 8 et 9, dans l'enceinte 3, - une calculatrice 10 dans laquelle ont été d'une part introduites sous forme de programme les algorithmes de résolution des diverses équations Y envisagées, d'au n tre part mémorisées les données permanentes et coefficients, cette calculatrice étant reliée aux sondes 4 à 9 dont elle peut par conséquent recevoir les valeurs mesurées par l'intermédiaire d'interfaces comportant des convertisseurs A/N, - une chaine de commande 11 reliée à la calculatrice 10 et propre à commander, sous l'influence de signaux émis par la calculatrice et par l'intermédiaire de dispositifs schématiquement représentés par des flèches 12, 13 et 14. l'alimentation en produits notamment chimiques appropriés des diverses enceintes 1, 2 et 3 et - eventuellement une imprimante 15 fournissant en clair les variations des diverses équations Y sur les n quelles on entend jouer et dont les données ont été fournies à la calculatrice. En rapport avec la susdite calculatrice, on a in diqué sur les figures 3a-3c et 4avec des logigrammes montrant l'une des approches possibles pour, tout d'abord, charger en mémoire de la calculatrice des courbes d'isoréponses sous forme discrète, d'autre part, rechercher l'optimum de la réponse Y par la méthode du gradient. Tracé et mémorisation des courbes isoréoonses Le cas illustré par le logigramme des figures 3a3c correspond à celui d'équations polynomiales dont certains termes sont des fonctions d'un des paramètres et d'autres des fonctions de deux paramètres. On commence par indiquer, en mémoire de travail de la calculatrice, le nombre N de facteurs expérimentaux X(i) qui interviennent dans la détermination de Y On introduit ensuite I et J, qui sont les deux valeurs de i constituant chacun des couples pour lequel on va tracer des courbes d'isoréponses. Au cours du déroulement du programme, I et J prendront évidemment tous les couples de valeurs possibles entre I et N, à l'exception des couples pour lesquels I=J, et pour chaque couple (X(i). X(J)), on affectera, aux N-2 variables restantes, tous les jeux de valeurs constantes Une fois ce chargement effectué, on introduit (cadre 20) l'équation de la fonction réponse que, pour simplifier, on a supposé être de la forme linéaire avec interactions EMI12.1 Avant de tracer les courbes, on va déterminer la valeur de Y pour tous les points d'une grille rectangulaire correspondant à des valeurs de X(i > et X(j) régulièrement réparties Pour que les grilles correspondant à tous les couples de valeurs aient le même nombre de points. il faut effectuer une normalisation et un centrage tenant compte de la plage de variations possibles pour chacun des paramètres X. Cela revient à faire varier le nombre représentant chaque paramètre X entre -1 et +1, ce qui implique d'introduire un facteur d'échelle Pour cela, on introduit en mémoire les valeurs limites des paramètres (cadre 22) puis on détermine. pour chaque couple de paramètres correspondant à une courbe d'isoréponse, le facteur d'échel- le Le cadre 24 correspond à une grille à P points pour X(I) et Q points pour X(J), correspondant à un premier jeu de valeurs constantes de tous les X autres que X(I) et X(J). Le tracé des courbes peut alors être effectué à partir des valeurs minimales prévues pour X(I) et X(J) (par exemple température ambiante, d'une part, concentra tion nulle, d'autre part). Après initialisation à ces valeurs en 26, on effectue le calcul de tous les points en séquence, en mettant en oeuvre deux boucles imbriquées. Le cadre montré en 28 correspond à la détermination de tous les points de la grille correspondant à une valeur donnée d'un des paramètres avant de passer à une nouvelle valeur de ce paramètre pour recommencer le calcul. On détermine ainsi les valeurs prises par Y pour tous les points de la grille. On trace ensuite les courbes d'isoréponses par une interpolation dont les paramètres sont introduits en 30. Le pas d'écartement des isoréponses (valeurs de Y pour lesquelles on veut disposer des courbes) sera choisi en fonction de la finesse d'ajustement que l'on recherche. L'opération sera effectuée à partir des points de la grille, en utilisant un programme d'interpolation qui peut être classique. A l'issue de l'opération, schématisée en 32, toutes les valeurs d'isoréponses seront affichées et mémorisées. La même opération sera répétée pour tous les couples X(I), X(J) avec tous les jeux de valeurs constantes affectées aux paramètres fixes. Une fois les données ainsi stockées dans la calculatrice, cette dernière peut être utilisée pour optimiser la fonction Y par la méthode du gradient. Ootimisation de Y Dar la méthode du aradient Le logigramme montré en figures 4a à 4c reprend une fraction du logigramme précédent. Mais il ne met pas en oeuvre cette partie commune pour la totalité du domaine de variation des paramètres X mais seulement autour du point expérimental, au lieu de repartir des valeurs minimalles. Une installation fonctionnant suivant le logigramme des figures 4a à 4c constitue donc un ensemble autonome. Ce logigramme prévoit des contraintes, visant à diminuer la durée des calculs et à éviter des réglages trop fréquents des paramètres, Le parcours des boucles de calcul est arrêté lorsque le gain sur Y obtenu par un parcours supplémentaire est inférieur à une valeur prédéterminée. De plus, les valeurs des réponses doivent rester à l'intérieur d'un domaine de validité défini par un rayon R. La figure 4a montre les opérations d'initialisation du programme On introduit M et l'équation de la fonctionrreponse Y. c'est-à-dire des données identiques à celles utilisées au début du programme de tracé de courbes. Pour appliquer la méthode du gradient, on introduit ensuite les dérivées partielles D[X(i)] = #Y;/#Xi de Y; par rapport à un premier paramètre Xi (cadre 34) puis les valeurs limites et l'algorithme du gradient (cadre 36 > avec une valeur prédéterminée, la même pour tous les X. On détermine ensuite en 38 (figure 4b) le sens de l'influence sur Y d'une première variation de 1 ensemble des paramètres avant d'effectuer la recherche de l'optimum, par application de l'algorithme du gradient. il est important de remarquer que le logigramme réduit le coefficient multiplicateur appliqué à la dérivée lorsque la dernière itération précédente amène à dépasser le rayon de validité R et recommence le même calcul, mais avec la valeur réduite Dans l'exemple représenté, cette réduction est faite dans un rapport Z. Si le nouveau calcul conduit encore à dépas- ser le rayon R, une nouvelle réduction dans le même rapport est appliquée : la précision est ainsi très accrue. Le logigramme limite par ailleurs le nombre de parcours de la boucle en arrêtant les itérations dès que l'on arrive à proximité de l'optimum, ce qui se traduit par le fait que la différence entre deux résultats successifs est inférieure à un seuil DY ou que. à la fois, on a dépassé le rayon R et que le coefficient L appliqué à la dérivée dans l'algorithme du gradient est tombé au-dessous d'une valeur déterminée (test 44. figure 4c) La figure 5 montre, à titre de simple exemple, une installation susceptible de mettre en oeuvre le procédé qui vient d'être défini. Sur cette figure. les organes correspondant à ceux de la figure 2 sont désignés par le même numéro de référence. Dans la réalisation montrée en figure 5, un circuit 46 de scrutation périodique échantillonne, par exemple une fois par heure, les paramètres X, qui dans le cas illustré sont au nombre de deux et fournis respectivement par les sondes 5 et 6. La valeur échantillonnée est appliquée à un convertisseur analogique-numérique 48 et envoyée à un organe de calcul 10, qui peut être un microprocesseur ou un automate programmable. Les nouvelles valeurs de consigne des paramètres X élaborées par 1 organe de calcul 10 sont envoyées à des cartes 50 et 52 de réglage des paramètres. Dans les systèmes à cinétique lente, se prêtant mal à une régulation, la carte 50 pourra être prévue pour élaborer, à partir de données connues sur le système, de la valeur de consigne qui vient d'être fournie par 1 organe de calcul 10 et de la valeur réelle fournie par la sonde 5, la quantité de réactif à ajouter au système pour obtenir la valeur de consigne. L'opération de réglage peut alors être effectuée par commande d'ouverture d'une vanne 52 pendant une durée déterminée. Dans le cas où la cinétique est plus rapide, on peut utiliser une carte 52 qui reçoit également, en vue d'une régulation, un signal de sortie provenant du convertisseur 48. le circuit de scrutation 46 fournissant alors en permanence la valeur du paramètre correspondant. afin de permettre d'en suivre l'évolution. La carte 52 arrête alors l'ouverture de la vanne correspondante 54 dès que la valeur de consigne est atteinte. C'est en mettant en oeuvre les procédé et installation conformes à l'invention qu'on a réalisé les exemples suivants. EXEMPLE 1 - Détermination des équations prévisionnelles représentatives de deux propriétés d'une sur face métallique traitée dans un bain de con version chimique. Dans cet exemple, il a été choisi de modéliser et donc de prévoir deux propriétés d'un revêtement de conversion couramment mesurées en contrôle dans les traitements de phosphatation, à savoir : - le poids de couche exprimé en g/mê, - la densité de structure cristalline exprimée en nombre de cristaux par unité de surface et mesurée à l'aide d'une photographie prise au microscope électronique à balayage avec un grossissement de 1500 fois. Cinq variables X1 à X5 du bain de conversion ont été prises en considération X1 : concentration en ions nickel X2 : concentration en ions chlorate X3 : concentration en ions zinc concentration en ions phosphate X5 : valeur de l'acidité libre ou AcL étant entendu que, par AcL, on désigne la valeur obtenue exprimée en ml de NaOH N/10 sur un échantillon de 10 ml du bain dosé par NaOH N/10 au virage de l'hélianthine. Des éprouvettes en tôle d'acier laminé à froid de qualité courante sont traitées suivant une gamme au trempé à base de produits commercialisés par la Société Demanderesse et couramment employés dans l'industrie Automobile pour la préparation des carrosseries avant peinture. Les étapes et produits mis en oeuvre comprennent - un dégraissage alcalin à chaud avec utilisation des agents dégraissants à base de sels alcalins et de tensio-actifs commercialisés par la Société Demanderesse respectivement sous les marques de fabrique "RIDOLINE 1550 CF/2" et "RIDOSOL 550 CF", - un rinçage froid à l'eau courante, - un conditionnement de surface avec utilisation d'un agent d'affinage à base de sel de titane commercialisé par la Société Demanderesse sous la marque de fabrique "FIXODINE 5" - un traitement de conversion chimique avec utilisation d une solution acide à base de phosphate de zinc et autres ions classiques commercialisée par la Société Demanderesse sous la marque de fabrique "GRANODINE 908", - un rinçage à froid, - un rinçage passivant avec utilisation d'un agent à base d'ions chrome commercialisé par la Société Demanderesse sous la marque de fabrique "DEOXYLYTE 41", - un étuvage. Les cinq variables du bain de conversion chimique qui sont prises en considération sont mesurées par tout moyen classique habituellement utilisé : dosage chimique, pH-métrie, potentiométrie, conductivité ionique et autres. On effectue 16 expériences pour déterminer les constantes de l'équation prévisionnelle correspondante qui s'écrit Y1 = b0 + b1X1 + b2X2 + b3X3 = b4X4 + b5X5 + b12X1X2 + b13X1X3 + b14X1X4 + b15X1X5 + b23X2X3 + b24X2X4 + b25X2X5 + b34X3X4 + b35X3X5 + b45X4X5 et on mesure pour chaque expérience - le poids de couche de chaque éprouvette traitée, par dissolution du revêtement dans une solution d'acide chromique à chaud (à 10 % 60C, 15 minutes), - le nombre de cristaux présents à l'intérieur d'une fenetre de 7,5 cm de diamètre de chaque éprouvette traitée, ces cristaux étant comptés à partir d'une photographie prise au microscope électronique à balayage avec un grossissement de 1500 fois. Pour déterminer les caractéristiques de ces expériences de façon à pouvoir déterminer les susdites constantes, on procède comme suit Pour chacun des cinq facteurs pris en considéra tion, on définit un domaine de variations comme suit ++ 0,8 g/1 # Zn # 1,4 g/1 15 g/1 # PO43- # 23 g/1 0,3 g/l # ClO - 0,9 g/l 0,5 g/l ( Ni # 0,8 g/l 0,9 ml NaOH N/10 ( AcL # 1,3 ml NaOH N/10. Compte tenu de l'équation précédente, il faut déterminer p (p+1) coefficients liés aux termes du 1er degré 2 et aux termes rectangles (p représentant le nombre de pa- ramètres avec, dans le cas présent p = 5), auxquels vient s'ajouter le terme constant Cela conduit à réaliser les 16 expériences figurant dans le tableau suivant pour lesquelles chacun des paramètres ne prend que les valeurs extrêmes de son domaine de variation. Exp. Ni++ ClO3- Zn++ PO43- ACL &gamma;1 p.d.c. &gamma;1 n.c. nb1 (g/m2) 1 0,5g/1 0,3g/1 0,8g/1 15g/1 1,3ml 1,74 147 2 0,5 g/1 0,9 g/1 0,8 g/1 15g/1 0,8ml 1,14 460 3 0,8g/l 0,9g/1 0,8g/1 15g/l 1,3ml 1,19 168 4 0,5g/1 0,3g/1 1,4g/1 15g/1 0,8ml 1,61 117 5 0,5g/1 0,9g/1 1,4g/1 15g/1 1,3ml 1,78 91 6 0,8g/1 0,3g/1 0,8g/1 23g/1 1,3ml 1,39 128 7 0,5g/1 0,3g/1 1,4g/1 23g/1 1,3ml 1,86 101 8 0,8g/1 0,9g/1 1,4g/1 15g/1 0,8ml 1,29 81 9 0,8g/1 0,9g/1 0,8g/1 23g/1 0,8ml 1,19 205 10 0,8g/1 0,9g/1 1,4g/1 23g/1 1,3ml 1,41 66 11 0,5g/1 0,9g/1 1,4g/1 23g/1 0,8ml 1,61 88 12 0,5g/1 0,9g/1 0,8g/1 23g/1 1,3ml 1,36 157 13 0,5g/1 0,3g/1 0,8g/1 23g/1 0,8ml 1,34 126 14 0,8g/1 0,3g/1 0,8g/1 15g/1 0,8ml 1,39 181 t5 0,8g/l 0,3g/1 1,4g/1 23g/l 0,8ml 1,81 97 16 0,8g/1 0,3g/l 0,8g/1 15g/l 0,8ml 1,91 90 Pour s'affranchir des unités de facteurs, pour simplifier les calculs et les interprétations et pour comparer les effets des facteurs entre eux, les valeurs de chaque facteur sont centrées et réduites. Le calcul des coefficients b. se fait alors par la méthode des moindres carrés décrite en détail dans l'ouvrage cité précédemment. Seuls sont pris en compte dans l'équation, les coefficients les plus influents. Après traitement des données, les équations des deux propriétés du revêtement, à savoir respectivement le poids de couche pdc et le nombre de cristaux nc par unité de surface, s'écrivent Y pdc = 1,50 + 0,16 X3 - 0,13 X2 + 0,079 X5 - 0,070 X4X5 (I) Y nc = 119 - 39 X3 -17 X2 X3 14 X1 3 2 4 (Il) EXEMPLE 2 - Détermination de l'équation prévisionnelle représentative d'une propriété d'une surface métallique traitée dans un bain de conversion chimique. La propriété prise en considération est toujours le poids de couche. Pour la détermination de l'équation correspondante, on peut prendre en compte non seulement des constituants chimiques mais aussi des facteurs physiques tels que la température, la vitesse d'agitation et autres, ou bien des paramètres discrets tels que la nature ou l'état de surface d'un substrat. Dans le cas présent, on a déterminé l'équation Ypdc à partir de sept variables. Six de ces variables, respectivement X1, X2 X3, Xs, X5 et X7, concernent le bain de conversion chimique, à savoir - X1 : concentration en ions zinc X X2 : concentration en ions chlorate - X3 : concentration en ions phosphate - X4 : taux d'accélérateur (nitrite de sodium) - X5 : température du bain - X7 : valeur de l'acidité libre la septième, X6, représentant la concentration du bain de conditionnement de surface relatif à l'étape qui précède le traitement de conversion chimique proprement dit. En opérant comme à l'exemple 1 quant au mode opératoire, l'équation Y pdc obtenue s'écrit Y pdc = 2,15 - 0,19 X6 + 0,16 X1 - 0,13 X2 + 0,13 X5 + O,11 X3 - 0.10 X4 + 0,07 X7 (III). EXEMPLE 3 - Exemples pratiques d application. Un revêtement de conversion chimique tel que celui obtenu avec I agent de conversion commercialisé par la Société Demanderesse sous la marque "GRANODINE 908" doit réaliser l'adéquation suivante - son épaisseur ou poids de couche doit être suf 2 fisamment faible, de l'ordre de 2 à 2,5 g/m , pour permet- tre une bonne adhérence de la peinture, - sa structure cristalline doit etre la plus dense et la plus régulière possible pour minimiser la porosite et pour améliorer la protection qu'il confère contre la corrosion. C'est cette adéquation "adhérence peinture-densite de structure cristalline qu'il est proposé de résoudre à l'aide de l'équation indiquée à l'exemple 2. De façon générale, les paramètres d'un tel bain varient suivant les fourchettes ci-dessous 0,8 g/l < Zn < 1,4 g/l 15 g/l < POs < 23 g/l 0,3 g/l < Cl03 < 0,9 g/l 0,5 g/l < Ni < 0,8 g/l 0,9 ml NaOH N/10 < AcL < 1,3 ml NaOH N/10 50 C < température du bain de conversion < 60 C 1,0 ml < taux d'accélérateur < 2,0 ml 1,0 g/l < concentration du bain de conditionnement de surface "FIXODINE 5" < 3,0 g/l. On mesure tous les paramètres X1 -X7 à des instants successifs, par exemple à intervalles d'une demi-heure. Par exemple a) Mesure à l'instant t Concentration du bain de conditionnement de surface: 2 g/l Paramètres du bain de conversion chimique ++ Zn : 0,8 g/l ClO3- : 0,6 g/1 PO43- : 19 g/1 taux d'accélérateur : 1,5 ml température : 55' C acidité libre : 0,9 ml. On introduit cette série de valeurs dans l'équation (III) Y pdc, après les avoir centrées et réduites Y pdc = 2,15 - 0,19 X6 + 0,16 X1 - 0,13 X2 + 0,13 X5 + 0,11 X3 - 0,10 X4 + 0,07 XT (IV) On obtient alors, pour l'instant t : 2 Y pdc = 2,06 g/m Cette valeur de poids de couche qu'on obtiendra se 2 situe bien dans la fourchette comprise entre 2 et 2,5 g/m qui a été préconisée précédemment pour l'obtention d'une bonne adhérence de la peinture. Il n'y a donc pas lieu d'intervenir sur les paramètres. b) Mesure à l'instant t : Concentration du bain de conditionnement de surface: 1 g/l Paramétres du bain de conversion chimique ++ Zn : 1,4 g/l (X1) ClO3- : 0,3 g/1 (X2) PO43- : 19 g/1 (X3) taux d'accélérateur : 1,0 ml température : 55' C acidité libre : 1,0 ml. En introduisant cette série de valeurs dans l'équation (IV), on trove Y pdc = 2,77 g/mê. Cette valeur calculée est hors de la fourchette acceptable. Elle traduit une dérive d'un ou de plusieurs paramètres qu'il convient de corriger pour l'ajustement le plus économique. Pour cela, on utilisera les projections de l'hypersurface à 8 dimensions représentant l'équation (III) sur les- plans correspondant chacun à deux paramètres. Par exemple, la figure 6 montre, pour des valeurs constantes de X2, X3, X;, X5 et X7 correspondant à la mesure à l'instant t2-, les courbes d'isoréponse, ou courbes de niveau, dans le plan X1. X6, où X1 est la teneur en ions Zn++ du bain de conversion. X6 est la teneur en FIXODINE 5 du bain de condition nement de surface. Chaque courbe de niveau correspond à une même valeur du poids de couche Y (indiqué en g/mê). Les valeurs des autres paramètres sont X2 : 0,3 g/1 de ClO3 X3 : 19 g/1 de PO43 X4 : 1 ml d'accélérateur X5 : 55-C X7 . 1,0 ml d'acidité libre. La calculatrice va alors intervenir pour determi ner, à partir des valeurs de X, les variations optimales à provoquer dans- le bain de conversion chimique et/ou le bain -de conditionnement pour ramener Y dans la plage ac 2 ceptable, par exemple à une valeur proche de 2 g/m Pour cela, à chaque scrutation, toutes les valeurs X sont introduites en mémoire de travail, soit automatiquement, soit manuellement à la fin d'une analyse d'échan- tillon du bain. La calculatrice cherche alors, par la méthode du gradient, les variations optimales de chacun des paramètres conduisant à l'obtention de la valeur de poids de couche voulue. On rappelle que la méthode du gradient permet d'effectuer des déplacements sur l'hypersurface en chaque point perpendiculairement aux courbes isoréponses, le gra dient étant la direction qui permet de se déplacer le plus rapidement vers un optimum. On appelle gradient d'une fonction u = f (x,y,z), le vecteur dont les projections sur les axes de coordonnées sont les dérivées partielles correspondantes de la fonction donnée. du i + du i + du k grad u = dx dy dz i, j, k représentant les vecteurs unitaires du système de référence choisi. Le gradient d'une fonction de trois variables est en chaque point dirigé suivant la normale à la surface passant par ce point. La direction du gradient d'une fonction en un point donné est la direction suivant laquelle la fonction a la plus grande vitesse de croissance en ce point, c'est-à-dire pour 1 = grad u la dérivée du atteint dl sa plus grande valeur qui est égale à (du)ê+ (du)ê + (du)ê dx dy dz De meme on définit le gradient d une fonction à n variables u = f (x,y,z,...,n) grad u = du i + du j + du k + ... + du n grad u dx dy dz dn Celui-ci est en chaque point dirigé suivant la normale à l'hypersurface passant par ce point. Après traitement, la calculatrice fournit les valeurs optimales X1 : 1,09 g/l de Zn X2 . 0,56 g/l de ClO3 X3 : 16 g/l de P043 X4 : 1,33 ml X5 : 51 C X6 : 2,25 g/l X7 : 0,91 ml correspondant à Y = 1,94 g/mê(choisie inférieure à l'op- timum étant donné la poursuite de l'évolution). On supposera maintenant que l'ajustement doit se faire par - apport de FIXODINE 5 dans le bain de conditionnement. sous forme d'un prémélange aqueux contenant 3 g/1 de constituants actifs, - apport en produit phosphatant réalisé à partir de deux solutions aqueuses concentrées : I une A permettant la régénération du bain de con- version en fonctionnement normal, contenant des ions Zn++, PO43-, Ni++, par example dans la formu lation suivante (Z en poids) ++ 3,0 ( Zn < 6,5 30 < H3PO4 < 45 0,1 < Ni++ < 1,9 dont l'acidité libre est comprise entre 3 et 11 ml de NaOH N (prise d'essai (p.e.) 1 ml) et l'acidité totale est comprise entre 9,8 et 13.4 ml de NaOH N (prise d'essai : 1 ml), l'autre B permettant de corriger les augmentations de l'acidité libre et de la concentration en zinc du bain et contenant des ions PO43-, Ni++, par exemple, dans la formulation suivante 30 < NaH2PO4 < 40 0,1 < Ni++ < 1 i,2 < pH < 3,8 3 < AcT < 5 ml NaOH N (p.e. 1 ml) - apport en accélérateurs dans le bain de conversion chi mique réalisé séparément à partir de solutions aqueuses concentrées l'une C contenant des ions nitrites dans une formu lation comprise, par exemple, entre 16 et 20 Z, l'autre D contenant des ions chlorates dans une formulation comprise. par exemple, entre 20 et 30Z L'ajustage du bain de conditionnement de surface peut s'effectuer de façon très simple. Le calculateur 10 envoie à la carte de commande correspondante une valeur représentative de l'augmentation de teneurs à obtenir, c'est-à-dire 2,25-1= 1,25 g/l. La carte comporte un microprocesseur ou simplement une mémoire morte de codage qui détermine, à partir de cette valeur et du volume du bain, le volume de prémélange à ajouter et ouvre une électrovanne à débit constant d'amenée de prémélange pendant la durée nécessaire. L'ajustage du bain de conversion est un peu plus complexe. La calculatrice détermine les apports de solutions A, B, C et D requis pour amener le bain à la composition requise et envoie les informations correspondantes aux cartes de commande 56. Chacune de ces cartes comporte un circuit de détermination de la durée d'apport nécessaire à partir d'un réservoir correspondant 58 et un circuit de puissance qui ouvre l'électrovanne appropriée 60 pendant le temps requis. Une autre carte 62, adressable comme les cartes 56, reçoit l'information de température et commande un circuit de commutation 64 réglant la puissance électrique appliquée par une source 66 à une résistance chauffante 68, dans le sens de la diminution dans le cas considéré. * * * Il est également possible de réguler un bain du genre en question en effectuant un compromis entre plusieurs propriétés. Ainsi, on peut utiliser simultanément les équations prévisionnelles Y pdc (I) et Y nc (II) de l'exemple 1. Le système permet alors de faire évoluer les paramètres simultanément pour obtenir le meilleur compromis entre le poids de couche et la densité de structure. Comme il va de soi et comme il résulte d'ailleurs déjà de ce qui précède, l'invention ne se limite nullement à ceux de ses modes d'application et de réalisation qui ont été plus spécialement envisagés ; elle en embrasse, au contraire, toutes les variantes.

Claims: Method adapted to enable: on the one hand, to know, from the measurement of various developing parameters of a developing chemical medium, the particularities of the result to which the medium concerned leads at the moment selected, on the other hand, in the case where the value found by this measurement shows that the result concerned does not correspond to that desired, which must be comprised within a range of values limited by at least one value N fixed in advance, to determine first of all a set of modifications of the variables or parameters of the medium, the said set being proper to arrange so that the result to which the thus modified medium will lead, corresponds to that desired, then to apply this set of modifications to the developing chemical medium concerned, by action on the parameters concerned, this method, which comprises the steps: of, on the one hand, in a first time, establishing for the result to which this medium leads, a so-called predictional equation Y, of, on the other hand, the equation Y being established, of the use of said equation for the regulation of the developing chemical medium in proceeding to the repetitive measuring of the value of each of the said variables or parameters by as many probes or detectors and to the calculation after each of said repetitive measurings of the value of Y by introduction of the values measured for the various variables or parameters in the predictional equation, being characterized by the fact the predictional equation is of the type Y = b₀ + b₁X₁ + b₂X₂ ... + b k X k + ... + b n X n + b₁₂X₁X₂ + ... + b k-1,k X k-1 X k + ... + b kk X k ² + ... + b nn X n ² + ... in which Y represents the result to which the medium leads, X₁, X₂ .. X k , ... X n represent the adjustable variables or parameters of the developing chemical medium from which the result Y depends and b₀, b₁ ... b k , ..., b nn are coefficients with respect to which the linearity of the equation is expressed, the value of these coefficients being determined by calculations based on regression methods known in themselves, and by the fact that it is proceeded: to the determination by calculation, if the value found for Y does not correspond to the desired result, of all sets of variations or modifications to be imparted to the variables or parameters of the predictional equation adapted to provide the desired value for Y and to select, among all these sets of variations or modifications, that is to say among all of the possible corresponding solutions of the predictional equation Y, the one which is the most economical and the easiest to carry out having recourse to the gradient method, that is to say in carrying out successively: the initialization operation of the program by introduction of the equation Y and of the value N in respect to which Y is defined, the introduction of the partial derivatives D[X( i )= δY/δX i of Y with respect to a first parameter X i , then the limit values and the algorithm of the gradient with a predetermined value, the same for all the Xs, the determination of the direction of the influence on Y of a first variation of the set of parameters before carrying out the search for the optimum, by the application of the algorithm of the gradient, being understood that the logigram, on the one hand, reduces the multiplier coefficient applied to the derivative when the preceding latter iteration leads to exceeding the validity radius and recommences the same calculation, but with the reduced value and, on the other hand, limits the number of passages through the loop by stopping the iterations as soon as Y is brought back into the range of values limited by N, and to the action of as many regulating members as parameters, providing each of these regulating members with a signal adapted to bring the variable or parameter associated therewith to the value corresponding to the abovesaid solution of the predictional equation Y, which is the most economical and the most easy to obtain. Method according to claim 1, characterized by the fact that the set of the calculations is carried out by means of a digital computer which is supplied with suitable data bases. Method according to one of claims 1 and 2, characterized by the fact that the modifications to be imparted to the variables or parameters are calculated by the process represented in Figures 4a-4c.

Inventors: HOESSLER CHRISTIAN, KEN VICTOR, SCHAPIRA JOSEPH, HOESSLER, CHRISTIAN, KEN, VICTOR, SCHAPIRA, JOSEPH
Applicants: PRODUITS IND CIE FSE, COMPAGNIE FRANCAISE DE PRODUITS INDUSTRIELS
Application year: 1986
Main classification: B01J19
IPC codes: C23C22, G05B13, B01J19