FR 2686714 A1 19930730 FR 9201137 A 19920124 PROCEDE DE SIMULATION D'UN PROCESSUS INDUSTRIEL ET UTILISATION POUR TESTER LE FONCTIONNEMENT D'UN AUTOMATISME La présente invention a pour objet principal un procédé de simulation d'un processus industriel piloté par un automatisme de commande La complexité croissante des processus mis en oeuvre dans les installations industrielles actuelles, implique une complexité croissante des automatismes de commande qui les pilotent. Aussi est-il devenu impératif de tester le bon fonctionnement de tout automatisme de commande qui est soit à logique cablée, soit à logique programmée du type contrôleur programmable industriel, avant la mise en service de l'installation industrielle qu'il contrôle, afin d'éviter des coûts importants de mise au point sur site. Un procédé connu pour tester la validité d'un automatisme de commande, consiste à relier audit automatisme un simulateur qui reproduit fidèlement le processus de l'installation industrielle, de telle sorte que les entrées de l'un communiquent avec les sorties de l'autre et vice et versa. Le procédé de simulation de processus qui est utilisé par de tels simulateurs, qui sont la plupart du temps des systèmes basés sur des moyens informatiques, consiste à gérer des temporisations qui sont déclenchées à la réception par le simulateur d'une commande émise par l'automatisme Chaque fois qu'une temporisation touche à sa fin, ce qui correspond dans la réalité au changement d'état d'un capteur de l'installation industrielle, le simulateur rend compte de l'état du processus audit automatisme Ainsi, la simulation du processus industriel se fait en temps réel. L'inconvénient majeur d'un tel procédé de simulation réside dans la perte de temps consacrée à attendre la fin d'une temporisation; sachant que plus le processus simulé est lent, plus cette perte de temps est importante. Le but principal de la présente invention est donc de proposer un procédé de simulation d'un processus, qui peut être mis en oeuvre automatiquement, et qui élimine ces temps d'attente, tout en tenant compte du temps réel de l'évolution du processus. Cet objectif est atteint par un procédé de simulation d'un processus industriel qui fait varier la valeur de chaque grandeur physique caractéristique de son état, en fonction de l'état des signaux de commande émis par un automatisme de commande et qui rend compte de son état émis par des capteurs en fonction de la valeur de la grandeur physique qu'ils contrôlent, caractérisé en ce qu'il comprend les étapes suivantes: a) définir un modèle fonctionnel du contrôle du processus industriel, b) définir un modèle fonctionnel de la commande du processus, c) initialiser le processus industriel, d) calculer pour chaque grandeur physique Pi, une vitesse de variation Vi qui est fonction de l'état des signaux de commande, e) déterminer, pour chaque grandeur physique Pl dont la vitesse de variation n'est pas nulle, l'intervalle de temps Ati pendant lequel la valeur de la grandeur physique peut varier sans que l'un quelconque des capteurs contrôlant cette grandeur Pl change d'état, f) calculer l'intervalle de temps Lt, défini comme étant le plus petit des intervalles de temps At i déterminés à l'étape précédente, g) mettre à jour la valeur de chaque grandeur physique Pl en l'incrémentant de la valeur du produit de la vitesse de variation Vi par l'intervalle de temps ht, h) déterminer en fonction de la valeur des grandeurs physiques mises à jour à l'étape précédente, les capteurs qui changent d'état, et mettre à jour les signaux d'état correspondants, i) attendre un laps de temps, dit temps d'attente, dont la valeur est paramétrable, avant de recommencer le procédé à l'étape d. Un premier mode particulier de réalisation du procédé de simulation consiste à décomposer l'étape a précitée en la succession d'étapes suivantes: j) définir les capteurs qui contrôlent chaque grandeur physique jj) associer à chaque capteur un signal d'état, jjj) décomposer la plage de variation de chaque grandeur physique Pl en une succession de secteurs contigus; chaque secteur S i correspondant à une plage de variation de la grandeur physique Pl pendant laquelle aucun capteur contrôlant la grandeur Pl ne change d'état, jjjj) associer à chaque secteur Si 1 i une borne inférieure, une borne supérieure, et les états des capteurs contrôlant la grandeur Pi. Un deuxième mode particulier de réalisation du procédé de simulation consiste à décomposer l'étape b précitée, en la succession d'étapes suivantes: j) associer à chaque grandeur physique Pl les signaux de commande qui font varier la valeur de la grandeur Pi, jj) associer à chaque signal de commande défini en (j), une vitesse de variation de la grandeur Pl en fonction de l'état binaire dudit signal de commande. D'autre part de façon plus particulière, à partir des deux modes particuliers de réalisation précédents, les étapes c à i du procédé de simulation peuvent être effectuées automatiquement à l'aide de la succession des étapes suivantes 1) choisir une valeur pour le temps d'attente 2) choisir pour chaque grandeur physique Pl un secteur courant Si, parmi les secteurs qui composent la plage de variation de la grandeur physique Pl et initialiser une valeur pi qui est comprise entre les bornes du secteur courant Si, 3) calculer pour chaque grandeur physique Pl une vitesse de variation Vi à partir des vitesses de variation associées conformément à l'étape b (jj), à chaque signal de commande en fonction de son état, 4) déterminer pour chaque grandeur physique Pi, l'intervalle de temps A ti en fonction des bornes du secteur courant Si, de la vitesse de variation Vi et de la valeur pi de la grandeur physique. ) calculerat conformément à l'étape f, 6) mettre à jour la valeur de chaque grandeur physique P. conformément à l'étape g, 7) calculer pour chaque grandeur physique dont on a déterminé un intervalle de temps t ti égal à à t, un nouveau secteur courant en fonction du signe de la vitesse de variation Vi, 8) déterminer parmi les capteurs contrôlant une grandeur physique Pi, pour laquelle on a calculé un nouveau secteur courant à l'étape précédente, les capteurs dont l'état associé au nouveau secteur courant est différent de l'état associé à l'ancien secteur courant, 9) mettre à jour les signaux d'état correspondant aux capteurs déterminés à l'étape précédente, ) déclencher une temporisation dont la valeur initiale correspond à la valeur du temps d'attente choisie à l'étape 1, et attendre la fin de cette temporisation avant de reprendre automatiquement le procédé à l'étape 3 Enfin la présente invention propose une application du procédé de simulation permettant de tester le fonctionnement d'un automatisme de commande qui pilote un processus industriel, selon laquelle on met en oeuvre les étapes suivantes: A) simuler une première fois le processus industriel dans sa totalité avec un temps d'attente plus grand que le temps d'exécution estimé de l'automatisme de commande ( 1) afin de valider le séquencement logique des signaux de commande émis par ledit automatisme, B) réitérer la simulation avec des temps d'attente de plus en plus petits, jusqu'à obtenir un temps d'attente minimum pour lequel l'automatisme de commande séquence correctement les signaux de commande. Le procéde de simulation, peut-être appliqué à un automatisme de commande constitué par un automate programmable industriel Dans ce cas, on choisit la valeur du temps d'attente à l'étape A, de manière à ce qu'il soit supérieur au temps de cycle automate. D'autres caractéristiques et avantages de l'invention, ressortiront de la description suivante de modes particuliers de réalisation, donnée en référence aux dessins annexés sur lesquels: la figure 1 schématise un processus industriel qui reçoit des signaux de commande émis par un automatisme de commande et qui émet en retour des signaux d'état, la figure 2 représente un processus industriel particulier mis en oeuvre dans une installation industrielle comportant une cuve et un axe mécanique, les figures 3 a et 3 b sont des logigrammes d'état des signaux d'état en fonction de la position de l'axe mécanique et du niveau de liquide dans la cuve, et les figures 4 a, 4 b, 4 c et 4 d schématisent le processus industriel de la figure 2 à différents instants caractéristiques du procédé de simulation. Si l'on se reporte à la figure 1, on voit que le procédé conformément à l'invention vise à simuler un processus industriel ( 2) qui réagit à des signaux de commande ( 3) émis par un automatisme de commande ( 1), en faisant varier des grandeurs physiques caractéristiques de son état, et qui rend compte de son état audit automatisme par l'intermédiaire de signaux d'état ( 4). Un dispositif de mise en oeuvre du procédé selon l'invention est un système informatique dont les entrées sont reliées aux sorties de l'automatisme de commande, et vice et versa Un tel dispositif est déjà connu et ne fera pas l'objet d'une description plus détaillée L'automatisme de commande ( 1) peut être soit à logique cablée, soit à logique programmée du type contrôleur programmable industriel, encore appelé automate programmable industriel Cette automatisme de commande peut aussi être un système informatique simulant la logique de fonctionnement d'un automate programmable industriel ou d'une armoire à relais. Les grandeurs physiques caractéristiques d'un processus industriel sont de natures très diverses On peut citer comme grandeurs physiques sans que cela soit limitatif, la position d'un axe mécanique, le niveau de liquide dans une cuve, le degré de température dans un four, le niveau de la pression dans une chambre hermétique, la grandeur temporelle mise en oeuvre par exemple dans des temporisations, la vitesse de rotation d'un moteur ou enfin le poids d'un objet. La figure 2 représente un processus industriel particulier mis en oeuvre dans une installation industrielle simple, qui comprend un axe mécanique ( 21), piloté par un moteur ( 210) et une cuve ( 22) dont le remplissage dépend de l'ouverture d'une vanne ( 220) Les grandeurs physiques caractéristiques de l'état de ce processus, sont la position de l'axe mécanique ( 21) qui est repérée par une came ( 211) de largeur 10-2 m et le niveau de liquide dans la cuve ( 22) qui est repéré par la position de la surface ( 221) du liquide. Dans la suite de la description, on appelera Pl la position de l'axe et P 2, le niveau de liquide Conformément à l'invention, dans une première étape on définit un modèle fonctionnel du contrôle du processus (étape a) Pour cela, selon un mode particulier de réalisation, on commence par définir les capteurs qui contrôlent chaque grandeur physique (étape a(j)). Ainsi pour le processus particulier de la figure 2 le contrôle de la position de l'axe mécanique est réalisé à l'aide des capteurs ( 212) et ( 213) et le contrôle du niveau de liquide dans le cuve, à l'aide des capteurs ( 222) et ( 223). Ensuite on associe les capteurs ( 212), ( 213), ( 222) et ( 223) aux signaux d'état qu'ils délivrent (étape a(jj)) et qui sont respectivement les signaux ( 41), ( 42), ( 43) et ( 44). Les capteurs qui contrôlent le processus, ont une position fixe dans l'installation industrielle De cette position et de la largeur de la came ( 211), on peut déduire facilement les logigrammes des figures 3 a et 3 b qui représentent l'état respectivement des signaux ( 41), ( 42), et des signaux ( 43) et ( 44) en fonction de la valeur en mètre des grandeurs physiques Pl et P 2 A partir des logigrammes de la figure 3 a, on peut facilement repérer les plages de variation de la grandeur physique P 1, pendant lesquelles les capteurs ( 212) et ( 213) ne changent pas d'état De même les logigrammes de la figure 3 b font apparaître lesdites plages de variation de la grandeur physique P 2 pendant lesquelles les capteurs ( 222) et ( 223) ne changent pas d'état. Ainsi on décompose la plage de variation de la grandeur physique Pl en la succesion de secteurs contigus S 1, 51,2 et 51,3, et la plage de variation de la grandeur physique P 2, en la succession de secteurs contigus, 52,1 X 522 et 52,3 (étape a(jjj)). Ensuite, on associe au secteur S une borne inférieure qui est 0, une borne supérieure qui est 0 01, l'état du capteur ( 212) qui vaut 1, et l'état du capteur ( 213) qui vaut O On procède de même pour chacun des secteurs de la grandeur physique Pl Ensuite on associe au secteur 52,1 une borne inférieure qui est O, une borne supérieure qui est 2, l'état du capteur ( 221) qui vaut O et l'état du capteur ( 223) qui vaut O On procède de même pour chacun des secteurs de la grandeur physique P 2 (étape a(jjj)). Il résulte de cette succession d'étapes, un modèle fonctionnel du contrôle du processus Ce modèle peut être décrit dans un fichier informatique dit fichier de contrôle: FICHIER DE CONTROLE Première grandeur physique: Pl premier capteur: ( 212) signal d'état: ( 41) deuxième capteur: ( 213) signal d'état: ( 42) premier secteur: 51,1 borne inférieure: O borne supérieure: O état: Deuxième grandeur physique premier capteur: 1 deuxième capteur:0 deuxième secteur: 51,2 borne inférieure: 0 01 borne supérieure: 1 état: premier capteur: O deuxième capteur:0 troisième secteur: 51,3 borne inférieure: 1 borne supérieure: 1 01 état: premier capteur: O deuxième capteur:1 P 2 premier capteur: ( 222) signal d'état: ( 43) deuxième capteur: ( 223) signal d'état: ( 44) premier secteur: 52,1 borne inférieure O borne supérieure 2 état: premier capteur: O deuxième capteur:0 deuxième secteur: S 2 borne inférieure: 2 borne supérieure: 4 état: premier capteur: 1 deuxième capteur: O troisième secteur: 52,3 borne inférieure: 4 borne supérieure: 5 état: premier capteur: 1 deuxième capteur: 1 Ce fichier contient, selon le mode de réalisation particulier toutes les données relatives au fonctionnement du contrôle du processus industriel que l'on doit simuler Ces données seront utilisées ultérieurement pour simuler l'évolution du processus industriel Dans un mode de réalisation équivalent, on pourrait se contenter de spécifier, pour chaque secteur, dans le modèle fonctionnel du contrôle du processus, exclusivement les capteurs actifs (respectivement inactifs), un capteur non spécifié, pour un secteur donné, étant par déduction inactif (respectivement actif). Dans une deuxième étape, on définit un modèle fonctionnel de la commande du processus (étape b). Si on se référe à la figure 2, on voit que la commande du processus industriel ( 2) comprend un moteur ( 210) pour l'entraînement de l'axe mécanique et d'une vanne ( 220) pour la commande du remplissage de la cuve Le moteur ( 210) qui permet de piloter le mouvement de l'axe mécanique ( 21) et donc de faire varier la grandeur physique Pl est commandé par l'automatisme de commande ( 1) par l'intermédiaire des signaux de commande ( 31) et ( 32) De même, l'ouverture de la vanne ( 220) qui permet de remplir la cuve ( 22), et donc de faire varier la grandeur physique P 2, est commandée par ledit automatisme, à l'aide du signal de commande ( 33). Selon un deuxième mode particulier de réalisation, on associe à la grandeur physique P 1, les signaux de commande ( 31) et ( 32), et à la grandeur physique P 2, le signal de commande ( 33) (étape b(j)). Selon que le signal de commande ( 31) ou que le signal de commande ( 32) est actif, c'est à dire lorsque leur état vaut 1, le moteur ( 210) fait se déplacer l'axe mécanique ( 21) vers la droite à une vitesse de 0 05 ms, ou vers la gauche avec la même vitesse. Lorsque les signaux ( 31) et ( 32) sont inactifs, le moteur ( 210) est à l'arrêt, l'axe mécanique ( 21) est donc immobile De même, si le signal ( 33) est actif, la vanne ( 220) est ouverte et la cuve ( 22) se remplit avec une vitesse de lm/s Si le signal ( 33) est inactif, la vanne est fermée et le remplissage de la cuve s'arrête. Ainsi on peut associer au signal de commande ( 31), la vitesse de variation + 0 05 m/s ou O m/s selon que son état vaut 1 ou O, au signal de commande ( 32), la vitesse de variation 0 05 mis ou O m/s selon que son état 1 ou O, et au signal de commande ( 33), la vitesse de variation 1 m/s ou O mis selon que son état vaut 1 ou O (étape b(jj)). La différence de signe entre les vitesses associées aux états actifs des signaux de commande ( 31) et ( 32), s'explique par le fait que si le signal de commande ( 31) est actif la valeur de la grandeur physique Pl augmente, alors que si le signal de commande ( 32) est actif, la valeur de la grandeur physique Pl diminue Il résulte de cette succession d'étapes un modèle fonctionnel de la commande du processus Celui-ci peut être décrit dans un fichier informatique dit fichier de commande FICHIER DE COMMANDE première grandeur physique: Pl premier signal de commande: ( 31) vitesse état O: O vitesse état 1: + 0 05 m/s deuxième signal de commande: ( 32) vitesse état O: O vitesse état 1: 0 05 m/s deuxième grandeur physique: P 2 premier signal de commande: ( 33) vitesse état 0: O vitesse état: 1 m/s Une fois les modèles fonctionnels du contrôle et de la commande du processus définis, on peut démarrer la simulation du processus conformément à l'invention. Pour cela, on doit initialiser le processus dont l'état initial est représenté à la figure 4 a; l'axe mécanique ( 21) étant en position 0, et la cuve ( 22) contenant un liquide dont le niveau est à 2 mètres, l'état des capteurs ( 212) et ( 222) vaut 1 alors que l'état des capteurs ( 213) et ( 223) vaut O et les valeurs des grandeurs physiques Pl et P 2 sont respectivement O et 2 mètres. L'état des capteurs et les valeurs de Pl et P 2, permettent de dire que la valeur de la grandeur physique Pl est située dans la plage de variation correspondant au secteur S 1 et que la valeur de la grandeur physique P 2 est située dans la plage de variation correspondant au secteur S 2 ' 51 1 et S 2 étant définis précédemment Ainsi on définit pour Pl un secteur courant 51 qui est le secteur 511 et pour P 2 un secteur courant 52 qui est le secteur S 22 et on initialise la valeur pl de la grandeur physique Pl à O mètre et la valeur P 2 de la grandeur physique P 2 à 2 mètres (étape 2) A ce stade de la simulation on a donc: pour la grandeur physique Pl Pl = O m 51 51,1 valeur du signal d'état ( 41) = 1 valeur du signal d'état ( 42) = O pour la grandeur physique P 2 P 2 = 2 m il 52 = S 2,2 valeur du signal d'état ( 43) = 1 valeur du signal d'état ( 44) = O L'étape suivante conformément à l'invention consiste à calculer pour chacune des grandeurs physiques Pl et P 2 des vitesses de variation, appelées respectivement V 1 et V 2 (étape 3) Cette étape est la première étape d'une boucle itérative. Si l'on se référe à la figure 4 a, on remarque que les signaux de commande ( 31) et ( 33) schématisés par des fléchas en trait plein sont actifs alors que le signal de commande ( 32) schématisé par une fléche en pointillé est inactif Conformément à ce qui a été décrit, dans la partie de la description concernant la commande du processus et compte tenu de l'état des signaux de commande ( 31), ( 32) et ( 33) on déduit facilement que V 1 vaut + 0 05 m/s et V 2 vaut lm/s En effet, pour obtenir Vi, il suffit d'additionner la vitesse correspondant à l'état actif du signal de commande ( 31), et la vitesse correspondant à l'état inactif du signal de commande ( 32), ces vitesses ayant été spécifiées dans le fichier de commande conformément au deuxième mode particulier de réalisation De même V 2 est la vitesse correspondant à l'état actif du signal de commande ( 33). Ensuite on détermine pour la grandeur physique Pl (respectivement pour P 2) l'intervalle de temps Ct 1 (respectivement A t 2) qui correspond à l'intervalle de temps pendant lequel la grandeur physique Pl (respectivement P 2) peut varier sans que l'un des capteurs ( 212) et ( 213) (respectivement ( 222) et ( 223)) change d'état. Conformément au premier mode particulier de réalisation, A t 1 correspond au temps nécessaire au processus pour faire varier la valeur p 1 de la grandeur physique Pl de telle sorte que p 1 devienne égale soit à la borne supérieure du secteur courant S si la vitesse de variation 51 est positive, soit à la borne inférieure du secteur courant 51, si la vitesse de variation est négative. (borne inférieure de 51) Pl Ainsi si V 1 est négative, t 1 = V 1 (borne supérieure de 51) Pl et si V 1 est positive, t 1- V 1 Etant donné que V 1 est positive, que p 1 vaut O, et que le secteur courant de la grandeur physique Pl est le secteur 51 qui a pour borne supérieure 0 01, on trouve que Lt 1 vaut 0 2 s. De même étant donné que V 2 est positive et vaut lm/s, que P 2 vaut 2 m, et que le secteur courant de P 2 est le secteur S 2 qui a pour borne supérieure 4 m on trouve que At 2 vaut 2 s. On calcule alors l'intervalle de temps Lst qui est le plus petit des intervalles de temps àt 1 et At 2 (étape 5) Ainsi on trouve que At vaut At 1, soit 0 2 s. On met alors à jour la valeur de la grandeur physique Pl (respectivement P 2) en l'incrémentant de la valeur du produit de t et de V 1 (respectivement V 2), (étape 6). Ainsi pl qui valait O m, après un intervalle de temps /St de 0 2 s, vaut 0 01 m et P 2 qui valait 2 m, vaut 2 2 m. Le nouvel état du processus industriel est représenté à la figure 4 b. Conformément à l'invention, on calcule pour la grandeur physique P 1, un nouveau secteur courant Ce secteur est choisi automatiquement parmi les secteurs S 1, 51,2 et 51,3 qui sont associés à la grandeur physique P 1, comme étant le secteur dont la borne soit inférieure, soit supérieure vaut p 1, selon que la vitesse de variation V 1 est positive ou négative Ainsi dans le cas présent, le nouveau secteur courant de P,est le secteur qui a pour borne inférieure, la nouvelle valeur de la grandeur physique Pl V 1 étant positive Par conséquent le nouveau secteur courant de Pl est le secteur 51,2 (étape 7). Ainsi pour la grandeur physique Pl on a Pl = 0 01 m 51 = 51,2 et pour la grandeur physique P 2 on a: P 2 2 2 m 52 = 52,2 On détermine ensuite, parmi les capteurs ( 212) et ( 213), les capteurs dont l'état associé au nouveau secteur courant 51 est différent de l'état associé à l'ancien secteur courant (étape 8). Dans le fichier de contrôle, on voit que l'état du capteur ( 213) associé au secteur 51,1 est identique à l'état associé au 1,1 secteur 51,2 et vaut O alors que l'état du capteur ( 212) associé au secteur 51,1 est différent de celui associé au secteur 51,2 On en déduit que seul le capteur ( 212) change d'état en devenant inactif, et on modifie en conséquence le signal d'état ( 41) qui lui est associé dans le fichier de contrôle et qui vaut maintenant O (étape 9). Quant aux deux capteurs ( 222) et ( 223), ils ne changent pas d'état car l'intervalle de temps At est strictement inférieur à l'intervalle de temps At 2. A ce stade du procédé de simulation on a: pour la grandeur physique Pl: Pl = 0 01 m 51 = 51,2 V 1 =+ O 05 m/s valeur du signal d'état ( 41) = O valeur du signal d'état ( 42) = O pour la grandeur physique P 2 P 2 = 2 2 m 52 = 52,2 V 2 = + l'm/s valeur du signal d'état ( 43) = 1 valeur du signal d'état ( 44) = O L'automatisme de commande ( 1) est alors informé du changement d'état du signal d'état ( 41). On marque alors une pause dans le procédé de simulation, en déclenchant une temporisation (étape 10) afin de laisser le temps à l'automatisme de commande, de calculer s'il'y a lieu de nouveaux signaux de commande ( 31), ( 32) ou ( 33) Le choix de la valeur de la temporisation (étape 1) joue un rôle important dans l'utilisation du procédé de simulation pour tester le fonctionnement d'un automatisme Celle-ci sera décrite ultérieurement. Quand la temporisation, qui est généralement de l'ordre de quelques dizaines de millisecondes touche à sa fin, ce qui marque la fin de la première itération, on reprend le procédé tel que décrit à partir du calcul des vitesses de variation V 1 et V 2; ce qui correspond à la première étape de la boucle itérative précitée. Les après une itération, ( 31), ( 32) figures 4 c, et 4 d, montrent l'évolution du processus deuxième et une troisième itération A la deuxième on voit sur la figure 4 c que les signaux de commande et ( 33) n'ont pas changé d'état On a donc toujours: V = + O 05 m/s En appliquant précédemment, on et V 2 = + 1 m/s. le procédé de simulation tel que décrit trouve alors: Atl = 19 8 S; t 2 = 1 8 S; et = tt 2 = 1 8 s. on en déduit les nouvelles valeurs de Pl et P suivantes Pl = 0 1 m; P 2 = 4 m. Seul le capteur ( 223) change d'état en devenant actif alors la valeur des signaux d'état ( 41) ( 42) et ( 43) inchangée Quant au signal d'état ( 44), il est à présent Si on se reporte maintenant à la figure 4 d, on l'automatisme de commande, suite au changement d'état c ( 44), a changé la valeur du signal de commande ( 31 on a qui est actif. voit que Je signal 3) en le désactivant, ce qui a pour effet d'arrêter le remplissage de la cuve ( 22) en fermant la vanne ( 220). La troisième itération du procédé de simulation conduit alors à faire évoluer le processus industriel vers un nouvel état qui est représenté à la figure 4 d; l'intervalle de temps a t correspondant à cette nouvelle itération valant 18 s. Ainsi par itérations successives, on peut simuler le processus industriel dans sa totalité en fonction des signaux de commande ( 3) émis par l'automatisme de commande ( 1). L'avantage d'un tel procédé de simulation est de simuler le processus en accéléré, et non pas en temps réel, tout en gardant une trace de ce temps réel par l'intermédiaire de la valeur de l'intervalle de temps At à chaque itération. Les deux modes particuliers de réalisation des modèles fonctionnels de la commande et du contrôle du processus, permettent de réaliser automatiquement le procédé Tout autre procédé de modélisation de la commande et du contrôle du processus visant à obtenir le même résultat, et permettant notamment un calcul automatique de l'intervalle de temps ât, ne sortirait pas du cadre de l'invention. D'autre part le procédé de simulation selon ses modes particuliers de réalisation, peut facilement être mis en oeuvre par des moyens informatiques tels qu'un microcalculateur personnel relié à un automatisme de commande En outre, ces moyens informatiques peuvent avantageusement être aussi aptes à simuler la logique de fonctionnement de l'automatisme de commande du processus. Une autre caractéristique importante du procédé selon l'invention est le choix d'une valeur de temporisation que l'on a appelée temps d'attente La description qui va suivre a pour but de montrer une application du procédé pour valider le fonctionnement d'un automatisme de commande,qui peut être soit à logique cablée soit à logique programmée, en utilisant la caractéristique précitée. En effet, supposons que l'automatisme de commande ( 1), soit un automate programmable industriel, dont le programme est conçu de telle sorte qu'en fonction des signaux d'état ( 41) ( 42) ( 43) ( 44), il calcule les signaux de commande ( 31) ( 32) ( 33) Il est intéressant de détailler ce qui se passe si l'on simule le processus industriel vis à vis del'automate, avec un temps d'attente inférieur à un cycle d'automate En effet, si on applique le procédé de simulation précédemment décrit, avec un temps d'attente nul, les deux premières itérations ne sont pas altérées car les signaux de commande ( 31) ( 32) et ( 33) ne changent pas En revanche, à la troisième itération, du fait d'un temps d'attente nul, l'automate n'a pas le temps de désactiver le signal de commande ( 33), afin d'arrêter le remplissage de la cuve On en conclut, que P 2 ne vaut alors plus 4 m à la fin de la troisème itération, mais vaut 22 m; la cuve continuant à se remplir à la vitesse de 1 m/s Cette caractéristique du procédé de simulation conduit à spécifier un temps d'attente, dont la valeur, lorsque l'automatisme de commande est un automate programmable industriel, doit être au moins de un cycle d'automate En effet, pendant ce cycle, l'automate lit toutes ses entrées, ce qui correspond à une prise en compte des signaux d'état ( 4), puis effectue le calcul de ses sorties à partir de ses entrées et de bits internes, et enfin rafraîchit l'état de ses sorties; ce qui correspond à l'émission des signaux de commande ( 3) Il faut donc au moins un cycle pour que l'automate prenne en compte le dernier changement du signal d'état ( 44), afin de mettre à jour le signal de commande ( 33) qui arrête le remplissage de la cuve. Ainsi pour tester le fonctionnement d'un automate programmable industriel, on simule dans un premier temps (étape A) le processus industriel ( 2) qu'il pilote, avec un temps d'attente, plus important que son temps de cycle A titre indicatif, on choisit un temps d'attente qui équivaut à cinq fois le temps de cycle automate. Ainsi dans le cas présent, on laisse le temps à l'automate à la troisième itération, de calculer un nouveau signal de commande ( 33) En analysant à chaque itération successive les signaux de commande émis par l'automate, on valide le séquencement logique desdits signaux. Une fois cette validation effectuée, dans un deuxième temps (étape B), on réitère plusieurs fois la simulation globale du processus, mais avec des temps d'attente de plus en plus petits, jusqu'à l'observation d'une différence dans le séquencement des signaux de commande Ainsi dans le cas du remplissage de la cuve, cette divergence apparaît lorsque le signal de commande ( 33) n'est pas désactivé à temps On peut alors déterminer un temps d'attente minimum, qui correspond au temps de calcul maximum de l'automate. Ce temps d'attente minimum dépend évidement de l'automate et de son temps de cycle, mais aussi de la structure du programme, notamment dans la gestion des bits internes intervenant dans les calculs intermédiaires Il ressort de cette description que le temps d'attente est dans le meilleur des cas un cycle automate. Dans certaines applications industrielles, ce renseignement sur le temps de calcul réel de l'automate peut avoir une importance capitale, soit au regard de la sécurité, soit au regard de la précision du processus industriel mis en oeuvre. L'application particulière de cette utilisation du procédé de simulation pour tester le fonctionnement d'un automate programmable est donnée à titre indicatif On peut de la même manière utiliser le procédé pour valider le fonctionnement de tout automatisme de commande, en choisissant pour valeur initiale du temps d'attente à l'étape A, un temps supérieur au temps estimé d'exécution de l'automatisme de commande. Il apparaît donc que le procédé de simulation et son utilisation selon l'invention trouvent leur application dans l'étape importante de l'automatisation d'installations industrielles, qui consiste à certifier la conformité de l'automatisme avec le cahier de charges de l'installation Ils permettent en effet, dans un premier temps de valider le fonctionnement logique de l'automatisme avant la mise en oeuvre du processus industriel et dans un deuxième temps, de quantifier le temps d'exécution de l'automatisme; ce qui permet de certifier le respect des contraintes temporelles liées à une installation industrielle donnée. REVENDICATIONS 1 Procédé de simulation d'un processus industriel ( 2) qui fait varier la valeur de chaque grandeur physique caractéristique de son état, en fonction de l'état des signaux de commande ( 3) émis par un automatisme de commande ( 1) et qui rend compte de son état audit automatisme par l'intermédiaire de signaux d'état ( 4) émis par des capteurs en fonction de la valeur de la grandeur physique qu'ils contrôlent, caractérisé en ce qu'il comprend les étapes suivantes: a) définir un modèle fonctionnel du contrôle du processus industriel, b) définir un modèle fonctionnel de la commande du processus, c) initialiser le processus industriel, d) calculer, pour chaque grandeur physique Pi, une vitesse de variation Vi qui est fonction de l'état des signaux de commande ( 3), e) déterminer, pour chaque grandeur physique Pl dont la vitesse de variation n'est pas nulle, l'intervalle de temps t ti pendant lequel la valeur de la grandeur physique peut varier sans que l'un quelconque des capteurs contrôlant cette grandeur Pl change d'état, f) calculer l'intervalle de temps L t, défini comme étant le plus petit des intervalles de temps At i déterminés à l'étape précédente, g) mettre à jour la valeur de chaque grandeur physique Pl en l'incrémentant de la valeur du produit de la vitesse de variation Vi par l'intervalle de temps L Vt, h) déterminer en fonction de la valeur des grandeurs physiques mises à jour à l'étape précédente, les capteurs qui changent d'état, et mettre à jour les signaux d'état correspondants, i) attendre un laps de temps, dit temps d'attente, dont la valeur est paramétrable, avant de recommencer le procédé à l'étape d). 2 Procédé de simulation selon la revendication 1 caractérisé en ce qu'il consiste à décomposer l'étape a) en la succession d'étapes suivantes: j) définir les capteurs qui contrôlent chaque grandeur physique, il) associer à chaque capteur un signal d'état, jjj) décomposer la plage de variation de chaque grandeur physique Pl en une succession de secteurs contigus, chaque secteur S J correspondant à une plage de variation de la grandeur physique Pl pendant laquelle aucun capteur contrôlant la grandeur Pl ne change d'état, jjjj) associer à chaque secteur Si i une borne inférieure, une borne supérieure, et les états des capteurs contrôlant la grandeur Pi. 3 Procédé de simulation selon la revendication 1 ou 1; revendication 2 caractérisé en ce qu'il consiste à décompose: l'étape b) en la succession d'étapes suivantes: j) associer à chaque grandeur physique Pl les signaux de commande qui font varier la valeur de la grandeur Pi, jj) associer à chaque signal de commande défini en (j), uni vitesse de variation de la grandeur Pl en fonction de l'étal binaire dudit signal de commande. a r e t 4 Procédé de simulation selon les revendications 2 et 3 caractérisé en ce que les étapes c à i, sont effectuées automatiquement à l'aide de la succession des étapes suivantes 1) choisir une valeur pour le temps d'attente, 2) choisir pour chaque grandeur physique Pl un secteur courant Si, parmi les secteurs qui composent la plage de variation de la grandeur P et initialiser une valeur pi qui est comprise entre les bornes du secteur courant Si, 3) calculer pour chaque grandeur physique Pl une vitesse de variation Vi à partir des vitesses de variation associées conformément à l'étape b (jj), à chaque signal de commande en fonction de son état, 4) déterminer pour chaque grandeur physique Pi, l'intervalle de temps bti en fonction des bornes du secteur courant Si, de la vitesse de variation Vi et de la valeur pi de la grandeur physique, 5) calculer At conformément à l'étape f, 6) mettre à jour la valeur de chaque grandeur physique P. conformément à l'étape g), i 7) calculer pour chaque grandeur physique dont on a déterminé un intervalle de temps A ti égal à Et, un nouveau secteur courant en fonction du signe de la vitesse de variation V, 8) déterminer parmi les capteurs contrôlant une grandeur physique Pl pour laquelle on a calculé un nouveau secteur courant à l'étape précédente, les capteurs dont l'état associé au nouveau secteur courant est différent de l'état associé à l'ancien secteur courant, 9) mettre à jour les signaux d'état correspondant aux capteurs déterminés à l'étape précédente, 10) déclencher une temporisation dont la valeur initiale correspond à la valeur du temps d'attente choisie à l'étape 1, et attendre la fin de cette temporisation avant de reprendre automatiquement le procédé à l'étape 3, 5 Procédé de simulation selon l'une quelconque des revendications 1 à 4, pour tester le fonctionnement d'un automatisme de commande ( 1) pilotant un processus industriel ( 2), caractérisé en ce que l'on met en oeuvre les étapes suivantes: A) simuler une première fois le processus industriel ( 2) dans sa totalité avec un temps d'attente plus grand que le temps d'exécution estimé de l'automatisme de commande ( 1) afin de valider le séquencement logique des signaux de commande ( 3) émis par ledit automatisme, B) réitérer la simulation avec des temps d'attente de plus en plus petits, jusqu'à obtenir un temps d'attente minimum pour lequel l'automatisme de commande séquence correctement les signaux de commande ( 3). 6 Procédé de simulation selon la revendication 5, appliqué à un automatisme de commande constitué par un automate programmable industriel, caractérisé en ce que le temps d'attente à l'étape A, est choisi de manière à ce qu'il soit supérieur au temps de cycle automate.