FR 3078161 A1 20190823 FR 1851515 A 20180222 Méthode de simulation de la puissance optique d'un verre feuilleté La présente invention concerne une méthode de simulation de la puissance optique d'un verre feuilleté. Elle est particulièrement adaptée à la détermination d'une valeur minimale d'indice de qualité optique à laquelle une combinaison de feuilles de verre doit satisfaire pour former un verre feuilleté de forme donnée avec une contrainte donnée de qualité optique. L'invention a également pour objet un système de traitement de données comprenant des moyens de mettre en oeuvre ladite méthode. Les niveaux de qualité optique exigés pour les verres feuilletés utilisés comme vitrage dépendent généralement de l'application. Les déviations angulaires que subit un faisceau lumineux en traversant un défaut optique présent dans un vitrage n'ont pas la même criticité selon que ce vitrage est destiné à procurer d'abord de la visibilité ou de la luminosité. Par exemple, un pare-brise de véhicule aérien ou terrestre requiert un niveau de qualité optique plus élevé que celui d'un vitrage destiné à procurer davantage de lumière dans un bâtiment. Dans le premier cas, aucun défaut optique ne doit avoir une taille et/ou une position susceptible de réduire le confort visuel du conducteur ou de gêner la conduite du véhicule car il en va de la sécurité de ses passagers. En revanche, dans le second cas, des défauts optiques de plus grandes tailles ou situés au milieu du vitrage peuvent être tolérés tant qu'ils ne dégradent pas l'esthétique et les performances thermiques globales du vitrage. De même, un vitrage utilisé dans un dispositif d'affichage numérique tel qu'un écran à plasma ou un écran à cristaux liquides ne doit pas présenter de défaut optique susceptible de perturber l'affichage des informations et leur perception par un utilisateur. En ce sens, de nombreuses méthodes de contrôle de la qualité optique des vitrages ont été développées de manière à pouvoir isoler et élimer ceux dont les caractéristiques optiques sont impropres à l'usage auquel ils étaient destinés. Ces méthodes sont souvent implémentées pendant ou après fabrication. Elles sont décrites en détail dans l'état de la technique. Par exemple, les brevets EP 0463940 B1 et EP0342127 B1 décrivent des procédés automatisables de contrôle d'un vitrage dans lequel les niveaux de déformation optique du vitrage sont déterminés à partir d'une image ombroscopique puis comparés à des valeurs seuils préalablement définies. Les demandes de brevets W09817993 et GB2152210 ainsi que le brevet EP 1061357 divulguent des procédés de détection des anomalies optiques d'une feuille transparente par analyse de l'image d'un motif géométrique réfléchi ou transmis par la feuille. Toutefois, ces méthodes de contrôle, qu'elles soient mises en oeuvre pendant ou après la production, aussi rapides qu'elles puissent être, ne peuvent pas éviter les pertes de 35 production. 2 Certaines demandes de brevet proposent des méthodes de simulation numérique pour anticiper d'éventuelles anomalies optiques liées à la mise en forme du vitrage. Par exemple, les méthodes de simulation fondées sur le lancer de rayons décrite dans les demandes de brevet 5 CN 105138759 et JP2008070946 peuvent être adaptées pour évaluer l'effet, sur la lumière transmise, des gradients d'indice de réfraction d'un vitrage après mise en forme. La demande de brevet WO 2021028630 divulgue un procédé de simulation du formage d'un vitrage exploitant une base de données comprenant des informations sur les altérations de forme d'une feuille de verre selon les distributions spatio-temporelles de température qu'elle peut o subir au cours du formage. Chaque forme peut ensuite faire l'objet d'une évaluation de sa qualité optique. Toutefois, ces méthodes ne permettent pas d'anticiper les difficultés techniques inhérentes à la fabrication de certaines formes de verre feuilleté de qualité optique donnée. Une démarche heuristique est alors aujourd'hui souvent utilisée afin de déterminer la qualité 15 optique minimale requise pour les feuilles de verre utilisées lors du formage. En effet, en- dehors des aléas ou des imperfections liées au procédé de fabrication, la qualité optique d'un verre feuilleté est conditionnée à la fois par sa forme et la qualité optique des feuilles de verre utilisées pour le fabriquer. Il arrive que les feuilles de verre d'une certaine qualité optique utilisées pour produire un vitrage d'une forme et d'une qualité optique données ne puissent 20 pas être directement utilisées pour produire un vitrage d'une autre forme et d'une autre qualité optique. La démarche heuristique consiste alors généralement, pour une forme de vitrage de qualité optique donnée, à expérimenter cette forme sur des feuilles de verre de qualité optique différentes jusqu'à obtenir le résultat recherché. La perte en temps et en matériaux n'est généralement pas négligeable. 25 D'autre part, il serait techniquement et économiquement peu pertinent de fabriquer un vitrage de forme et de qualité optique données à l'aide de feuilles de verre dont la qualité optique serait supérieure à l'usage auquel le vitrage est destiné. La solution qui consisterait à utiliser systématiquement des feuilles de verre de qualité optique supérieure pour produire tout type de verre feuilleté, y compris ceux pour lesquels certains types de défauts optiques sont 30 tolérables, reviendrait à ajouter des contraintes techniques inutiles sur les procédés de fabrication. En revanche, il serait techniquement plus avantageux de pouvoir déterminer a priori, avant toute fabrication, le niveau de qualité optique minimal que les feuilles de verre doivent posséder pour la fabrication d'un verre feuilleté de forme et de qualité optique donnée. Non seulement le gain en temps et en matériau serait plus important, mais aussi les difficultés techniques liées à la mise en forme du vitrage pourraient être anticipées. Or, les - 3 - méthodes de contrôle ou de simulation actuelles ne permettent pas d'obtenir une telle information. La présente invention résout ces problèmes. Elle concerne une méthode de simulation, mise en oeuvre par ordinateur, de la puissance optique d'un verre feuilleté de forme donnée 5 susceptible d'être obtenu en feuilletant au moins deux feuilles de verre et au moins un intercalaire de feuilletage, ledit verre feuilleté étant susceptible d'avoir un bord, une première face principale, et une deuxième face principale, et ladite méthode comprenant les étapes suivantes : al. la modélisation géométrique d'une première surface susceptible de 10 correspondre à la première face principale du verre feuilleté à partir d'un modèle numérique de la forme dudit verre feuilleté pour former une première surface modélisée ; 131. la modélisation géométrique d'une deuxième surface obtenue par translation de la surface modélisée dans l'étape (al) selon la direction de plus petite dimension du verre feuilleté d'une distance E correspondant à l'épaisseur du verre feuilleté, et selon la direction 15 de l'une des deux plus grandes dimensions du verre feuilleté d'une distance 8 définie, le rapport de la distance 8 sur la plus grande dimension du verre feuilleté étant supérieur à 0 et inférieur ou égal à 0,01, et ladite modélisation permettant d'obtenir une deuxième surface modélisée ; cl. le calcul du décalage géométrique t, en chacun des points des surfaces 20 modélisées, par soustraction de la première surface modélisée et de la deuxième surface modélisée ; dl. le calcul de l'épaisseur géométrique E', en chacun des points des surfaces modélisées, par addition de la distance E et du décalage géométrique T; el. le calcul de l'épaisseur topologique if, en chacun des points des surfaces 25 modélisées, par soustraction des profils topographiques des surfaces des feuilles de verre susceptibles de correspondre à la première face principale et deuxième face principale du verre feuilleté ; fi. le calcul de la puissance optique PO, en chacun des points des surfaces modélisées, à partir de la valeur du Laplacien de la somme des épaisseurs géométrique EI et 30 topologique Tc. Un premier avantage de la méthode de l'invention est qu'elle permet de contrôler, avant fabrication, toute combinaison des effets de forme du verre feuilleté et de la qualité optique des surfaces des feuilles de verre utilisées en première et deuxième faces 4 principales. Les combinaisons potentiellement défectueuses peuvent être anticipées et les pertes de production être ainsi avantageusement réduites. Un deuxième avantage de la méthode est sa rapidité d'exécution. Elle est, en termes de vitesse d'exécution, environ 50 fois plus rapide à mettre en oeuvre que les 5 méthodes de simulation actuelles fondées sur le lancer de rayons. Par exemple, à capacités comparables de traitement des données, une méthode de simulation fondée sur le lancer de rayons requiert 30 minutes pour obtenir la puissance optique d'un verre feuilleté. La méthode de simulation de l'invention permet d'obtenir la puissance optique de 50 verres feuilletés dans le même temps. Cet avantage rend possible l'utilisation en temps réel de la 10 méthode sur des lignes de production de verre feuilleté où des changements dans le planning de production sont susceptibles d'avoir lieu. En particulier, si une forme particulière de verre feuilleté doit être produite en urgence, la méthode de l'invention permet de contrôler la puissance optique susceptible d'être obtenue en quelques dizaines de secondes. 15 La modélisation géométrique des étapes (b1) et (cl) peut être réalisée par tout moyen adapté. Cette modélisation géométrique peut avantageusement être réalisée à l'aide des outils de la conception assistée par ordinateur (CAO) qui est une technique très répandue dans l'industrie. Ces outils permettent de créer des modèles numériques de pièces mécaniques, et sont particulièrement adaptés à la modélisation de la forme d'un 20 verre feuilleté. A l'étape (cl) de la méthode de simulation de l'invention, la deuxième surface modélisée est obtenue par translation de la première surface modélisée. Les première et deuxième surfaces modélisées sont généralement parallèles entre elles. Le rapport de la distance ô sur la plus grande dimension du vitrage est supérieur à 0 et inférieur ou égal à 25 0,01. Le rapport est, de préférence, compris entre 0,001 et 0,01. A l'étape (el) de la méthode de simulation de l'invention, les profils topographiques des surfaces de feuilles de verre susceptibles de correspondre à la première face principale et à la deuxième face principale du verre feuilleté sont utilisés pour le calcul de l'épaisseur topologique IT. Il a été constaté expérimentalement que la contribution des 30 profils topographiques des surfaces des feuilles de verre en contact de l'intercalaire de feuilletage, et la contribution de l'intercalaire lui-même étaient généralement négligeables. En effet, la comparaison des valeurs de puissance optique obtenues à l'aide de la méthode de simulation de l'invention avec celles obtenues à l'aide de méthode tenant compte des profils topographiques des surfaces des feuilles de verre en contact de l'intercalaire étaient relativement proches. Dans un verre feuilleté l'indice de réfraction de l'intercalaire de feuilletage est identique, sinon proche, de celui des feuilles de verre. Dans certains modes de réalisation de l'invention, il est néanmoins possible de tenir compte de la contribution des profils topographiques des surfaces des feuilles de verre en contact de l'intercalaire de 5 feuilletage, et de la contribution de l'intercalaire dans le calcul de l'épaisseur topologique IT. En raison de leur forme, notamment en raison de leur courbure générale, les verres feuilletés utilisés comme vitrage peuvent être affectés de distorsions optiques ayant des conséquences sur la puissance optique. En ce sens, dans un mode particulier de réalisation du procédé de l'invention, les épaisseurs géométrique El et topologique IT sont chacune 10 pondérées par un facteur d'amplification optique dépendant de l'angle, 61, formé entre la direction normale à la première surface et une direction de référence d'observation du verre feuilleté. A titre d'exemple, cette direction de référence d'observation peut être la direction d'observation d'un conducteur et/ou des passagers dans un véhicule terrestre ou aérien dans lequel le verre feuilleté est utilisé comme pare-brise. 15 La puissance optique en chacun des points peut par exemple être calculée à l'aide de l'expression suivante : PO = -(n - 1) * A[A,(61)n- + AE, (6)E1 où n est l'indice de réfraction du verre et A, et AE, (61) sont des facteurs d'amplification calculés à l'aide de formules suivantes : (Vn2 - sin2 0 - cost9) A,(61) = n - 1 (Vn2 - sin261 - cost9) 4(6) = (n - 1)cos2 0 20 Les profils topographiques des surfaces de feuilles de verre susceptibles de correspondre à la première face principale et à la deuxième face principale du verre feuilleté peuvent généralement être obtenus de deux façons différentes. Selon une première façon, les profils sont des profils mesurés par méthodes 25 profilométriques avec ou sans contact. Une méthode profilomètrique avec contact consiste, par exemple, à balayer la surface à l'aide d'un stylet en contact physique avec ladite surface et en en mesurant les déplacements verticaux. Dans le domaine de l'industrie du verre, les méthodes profilométriques sans contact sont généralement plus usuelles. Il s'agit souvent de méthodes profilomètriques optiques matriciels ou à balayage comme l'ombroscopie, l'interférométrie optique, ou encore l'holographie numérique. En particulier, les méthodes ombroscopiques sont très utilisées. Elles sont décrites avec en détail dans l'état de la technique, par exemple dans les demandes de brevet EP0463940 ou EP 0342127. Dans le cadre de la fabrication d'un verre feuilleté, les profils topographiques sont généralement 5 fournis par le fabriquant des feuilles de verre. Selon une deuxième façon, les profils topographiques sont des profils simulés numériquement. Il peut être avantageux de générer les profils topographiques à l'aide de fonctions mathématiques afin de contrôler, par anticipation, l'effet d'un ou plusieurs profils topographiques sur les valeurs de la puissance optiques d'un verre feuilleté avec une forme 10 donnée à l'aide de la méthode de l'invention. Une telle démarche permet, par exemple, d'établir des spécifications techniques sans avoir à recourir à des données expérimentales sur les profils topographiques. Ces spécifications techniques peuvent ensuite être transmises à un fournisseur de feuilles de verre avant tout approvisionnement afin qu'il puisse vérifier que les feuilles de verre qu'il est susceptible de fournir répondent aux besoins. La génération des 15 profils topographiques par des fonctions mathématiques peut également être réalisée à partir d'une bibliothèque de profils topographiques expérimentaux utilisés comme modèles. Dans certaines applications, plutôt que d'obtenir des valeurs de puissance optique en tout point du verre feuilleté, il peut être davantage pratique d'utiliser un paramètre synthétique représentatif de la puissance optique globale du verre feuilleté. En ce sens, dans un mode de 20 réalisation de l'invention, la méthode de simulation comprend en outre, après l'étape (f1), une étape de détermination d'un indice de fluctuation optique cp à partir des valeurs des puissances optiques PO obtenues en chacun des points des surfaces modélisées. Cet indice de fluctuation optique peut avoir la dimension d'une puissance optique. En particulier, cet indice de fluctuation optique cp peut être l'écart-type des valeurs des puissances optiques PO ou la 25 valeur maximale des puissances optiques PO. Lorsque le verre feuilleté est utilisé comme parebrise dans le domaine des véhicules terrestres ou aériens, en particulier dans le domaine automobile, il est commun d'utiliser un taux de variation de la puissance optique (« Rate of change » en anglais). Dans un mode particulier de réalisation du procédé de l'invention, l'indice de fluctuation optique cp peut 30 correspondre à ce taux de variation. Il peut alors être défini comme étant la valeur maximale d'un ensemble de valeurs correspondant aux différences entre les valeurs maximales et minimales des puissances optiques PO dans une fenêtre d'échantillonnage glissante. Dans la plupart des applications, la fenêtre d'échantillonnage peut avantageusement être une fenêtre 7 d'échantillonnage carrée ou rectangulaire dont les dimensions latérales sont comprises entre 2mm et 100mm, de préférence entre 5 et 50mm. Dans le cas particulier des pare-brise, seules la ou les régions du verre feuilleté correspondant à la ou les zones de vision du conducteur et/ou des passagers doivent 5 généralement avoir une qualité optique telle que la visibilité n'est pas affectée lors de la conduite du véhicule. Ainsi, dans un mode de réalisation de la méthode de simulation de l'invention, la détermination de l'indice de fluctuation optique cp est réalisée dans une région limitée des faces principales du verre feuilleté. Cette région limitée peut être définie conformément à la directive 92/22/EEC de l'Union Européenne ou conformément au o règlement 43 de la commission économique pour l'Europe des Nations Unies, en particulier au point 9.2.5.2.3 dudit règlement. Au sens de l'invention, cette région limitée peut correspondre à la région délimitée par l'intersection des faces principales dudit verre feuilleté avec une pyramide à base carrée dont l'apex est situé à une distance définie des faces principales et les deux angles au sommet entre les faces latérales en vis-à-vis de ladite pyramide 15 sont respectivement compris entre 100 et 20°, et entre 5 et 15°. La méthode de simulation de l'invention décrite ci-dessus est particulièrement adaptée à la détermination d'une valeur minimale d'indice de qualité optique à laquelle une combinaison de feuilles de verre doit satisfaire pour former un verre feuilleté avec une contrainte donnée de qualité optique. La méthode de simulation peut alors être intégrée comme étapes 20 intermédiaires dans une méthode permettant cette détermination. Ainsi, l'invention a également pour objet une méthode de détermination, mise en oeuvre par ordinateur, d'un intervalle d'indices de qualité optique requis pour une combinaison de feuilles de verre susceptibles de former un verre feuilleté de forme donnée susceptible d'être obtenu en feuilletant au moins deux feuilles de verre et au moins un intercalaire de feuilletage, ledit verre 25 feuilleté ayant une contrainte donnée de qualité optique, 0, et ladite méthode comprenant les étapes suivantes : a2. la sélection d'une combinaison de feuilles de verre parmi un ensemble de feuilles de verre auxquelles sont associés au moins un identifiant, un profil topographique et un indice de qualité optique ; 30 b2. le calcul de l'indice de qualité optique de la combinaison de feuilles de verre à partir des indices de qualité optique des deux feuilles de verre de la combinaison, l'une des surfaces de chacune desdites deux feuilles de verre est susceptible de correspondre respectivement à la première face principale ou à la deuxième face principale du verre feuilleté ; c2. la détermination d'un indice de fluctuation optique cp d'un verre feuilleté susceptible de comprendre la combinaison des feuilles de verre à l'aide d'une méthode de simulation selon l'un des quelconques modes de réalisation décrits ci-dessus ; d2. la répétition des étapes (a2) à (c2) avec une nouvelle combinaison de feuilles 5 de verre de manière à obtenir un ensemble de valeurs d'indices de fluctuation optique cp pour les différentes combinaisons possibles ; e2. la comparaison des valeurs des indices de fluctuation optique, cp, à la valeur de la contrainte de qualité optique O; f2. la sélection de l'intervalle des indices de qualité optique des combinaisons des 10 feuilles de verre satisfaisant le résultat de la comparaison de l'étape (e2). Cette méthode permet de déterminer, a priori, avant toute production, le niveau de qualité optique minimal qu'une combinaison de feuilles de verre doit posséder pour fabriquer un verre feuilleté de forme et de qualité optique données. L'avantage est un gain de temps et en matériaux non négligeable par rapport à une démarche heuristique expérimentale et une 15 anticipation plus efficace des difficultés liées à la mise en forme du verre feuilleté. La contrainte de qualité optique 0 peut avoir la dimension d'une puissance optique. Afin d'optimiser le temps nécessaire à l'exécution des étapes de ladite méthode, il n'est pas nécessaire de réaliser les étapes (al) à (dl) de la méthode de simulation de l'étape (c2) pour chaque combinaison de feuilles de verre sélectionnée à l'étape (a2). Pour chaque forme de 20 verre feuilleté pour laquelle est recherché l'intervalle d'indice de qualité optique des combinaisons de feuilles de verre permettant de satisfaire à la contrainte de puissance optique 0, les étapes (al) à (dl) de la méthode de simulation de l'étape (c2) peuvent être exécutées une seule fois. Le résultat peut être ensuite conservé pour l'exécution des autres étapes avec chaque combinaison de feuille de verre. 25 L'identifiant, le profil topographique et l'indice de qualité optique associés à chacune des feuilles de verre peuvent être sélectionnés à partir d'une base de données numérique comprenant des entités représentant les feuilles de verre. Chaque entité représentant les feuilles de verre possède alors, parmi ses attributs, un identifiant, un profil topographique et un indice de qualité optique. 30 Cette base de données peut être alimentée par les données transmises par les fabricants ou fournisseurs de feuilles de verre. Elle peut également comprendre des entités représentant des feuilles de verre dont les profils topographiques ont été obtenus par simulation numérique. En pratique, l'indice de qualité optique des combinaisons de feuilles de verre est un paramètre synthétique permettant de comparer les combinaisons entre elles. Cet indice peut notamment être une valeur numérique permettant de les trier. Dans un mode particulier de réalisation de la méthode de détermination d'un intervalle d'indices de qualité optique, l'indice 5 de qualité optique des combinaisons des feuilles de verre est la moyenne quadratique des indices moyens de qualité optique des deux feuilles de verre de la combinaison, l'une des surfaces de chacune desdites feuilles de verre étant susceptible de correspondre respectivement à la première face principale et à la deuxième face principale du verre feuilleté, lesdits indices moyens de qualité optiques étant obtenus par des méthodes ombroscopiques. 10 Par exemple, un indice de qualité optique, noté NO, peut être calculé à l'aide de la formule suivante : (NOFv1)2 + (NOFv2)2 NO 2 où NOFvi et NOFv2 sont respectivement les indices moyens de qualité optique la première feuille de verre et la deuxième feuille de verre obtenues à l'aide d'une méthode ombroscopique. Des exemples de méthodes ombroscopiques sont décrits en détail dans l'état 15 de la technique, par exemple dans les demandes de brevet EP0463940 ou EP 0342127. L'unité de mesure des grandeurs NO, NOnti et NOFv2 sont généralement la dioptrie. La méthode de simulation de la puissance optique et la méthode de détermination de la qualité optique de la combinaison de feuilles de verre sont mises en oeuvre par ordinateur. L'invention a également pour objet un programme informatique comprenant des instructions 20 exécutables par un ordinateur permettant l'exécution des étapes de ces deux méthodes de l'invention dans tous leurs modes de réalisation possibles. Tout type de langage de programmation compilé vers une forme binaire ou directement interprété peut être utilisé pour implémenter les étapes des méthodes par une suite d'instructions arithmétiques ou logiques exécutables par un ordinateur ou tout système de traitement de l'information 25 programmable. Le programme informatique peut faire partie d'un logiciel, c'est-à-dire d'un ensemble d'instructions exécutables et/ou d'un ou plusieurs jeux de données ou de bases de données. Le programme informatique peut être stocké sur un support de stockage déchiffrable par ordinateur. Ce support de stockage est de préférence une mémoire informatique non 30 volatile ou rémanente, par exemple une mémoire de masse magnétique ou à semi-conducteur (solid state drive, flash memory). Elle peut être amovible ou intégrée à l'ordinateur qui en déchiffre le contenu et en exécute les instructions. - 10 - Le support de stockage peut être intégré à un ordinateur distant, appelé « serveur », différent de celui qui exécute les instructions, appelé le « client ». Pour exécuter les instructions contenues dans le support de stockage, l'ordinateur « client » accède à l'espace mémoire de l'ordinateur « serveur » dans laquelle est enregistré le programme informatique à l'aide d'un 5 moyen de télécommunication physique et/ou aérien approprié. L'ordinateur « serveur » peut aussi déchiffrer le support de stockage sur lequel est stocké le programme d'ordinateur et communiquer les instructions sous forme binaire à l'ordinateur « client » par tout moyen de télécommunication. Il peut être avantageux que le support de stockage soit un support amovible ou soit 1 o accessible à distance par un moyen de télécommunication de manière à faciliter la diffusion de l'invention dans les lieux où elle est susceptible d'être utilisée. La méthode de détermination d'un intervalle d'indices de qualité optique requis pour une combinaison de feuilles de verre susceptibles de former un verre feuilleté avec une contrainte donnée de qualité optique, 0, peut être intégré à un procédé de fabrication d'un 15 verre feuilleté avec une contrainte de qualité optique. La présente invention concerne également un procédé de fabrication d'un verre feuilleté comprenant au moins deux feuilles de verre et au moins intercalaire, ledit verre feuilleté de forme donnée ayant un bord, une première face principale et une deuxième face principale, ledit verre feuilleté ayant une contrainte donnée de qualité optique, 0, ledit procédé comprend les étapes suivantes : 20 a. la sélection, mise en oeuvre par ordinateur, d'une combinaison de feuilles de verre parmi un ensemble de feuilles de verre auxquelles sont associés au moins un identifiant, un profil topographique et un indice de qualité optique, ladite combinaison de feuilles de verre étant susceptible de former ledit verre feuilleté ; 25 b. le calcul, mis en oeuvre par ordinateur, de l'indice de qualité optique de la combinaison de feuilles de verre à partir des indices de qualité optique des deux feuilles de verre de la combinaison, l'une des surfaces de chacune desdites deux feuilles de verre est susceptible de correspondre respectivement à la première face principale ou à la deuxième face principale du verre feuilleté ; 30 c. la modélisation géométrique, mise en oeuvre par ordinateur, d'une première surface correspondant à la première face principale du verre feuilleté à partir d'un modèle numérique dudit verre feuilleté pour former une première surface modélisée ; d. la modélisation géométrique, mise en oeuvre par ordinateur, d'une deuxième surface obtenue par translation de la surface modélisée dans l'étape a selon la direction de plus petite dimension du verre feuilleté d'une distance £ correspondant à l'épaisseur du verre feuilleté, et selon la direction de l'une des 5 deux plus grandes dimensions du verre feuilleté d'une distance ô définie, le rapport de la distance ô sur la plus grande dimension du verre feuilleté étant supérieur à 0 et inférieur ou égale à 0,01, et ladite modélisation permettant d'obtenir une deuxième surface modélisée ; e. le calcul, mis en oeuvre par ordinateur, du décalage géométrique -C obtenu par 10 soustraction de la première surface modélisée et de la deuxième surface modélisée en chacun des points des surfaces ; f. le calcul, mis en oeuvre par ordinateur, de l'épaisseur géométrique El par addition de la distance £ et du décalage géométrique -c en chacun des points des surfaces ; 15 g. le calcul, mis en oeuvre par ordinateur, de l'épaisseur topologique IT obtenu par soustraction des profils topographiques des surfaces des feuilles de verre susceptibles de correspondre à la première face principale et deuxième face principale du verre feuilleté en chacun des points des faces ; h. le calcul, mis en oeuvre par ordinateur, de la puissance optique PO à partir de la 20 valeur du Laplacien de la somme des épaisseurs géométrique El et topologique IT en chacun des points des surfaces modélisées ; i. le calcul, mis en oeuvre par ordinateur, d'un indice de fluctuation optique cp correspondant à la valeur maximale d'un ensemble de valeurs correspondant aux valeurs maximales des différences entre les valeurs maximales et minimales 25 des puissances optiques PO dans une fenêtre d'échantillonnage glissante ; j. la répétition des étapes (a) à (i) avec une nouvelle combinaison de feuilles de verre de manière à obtenir un ensemble de valeurs d'indices de fluctuation optique cp pour chaque combinaison possible ; k. la comparaison, mise en oeuvre par ordinateur, des valeurs des indices de 30 fluctuation optique, cp, à la valeur de la contrainte de qualité optique O; 1. la sélection, mise en oeuvre par ordinateur, de l'intervalle des indices de qualité optique des combinaisons des feuilles de verre satisfaisant le résultat de la comparaison de de l'étape (j) ; - 12 - m. la fabrication d'un verre feuilleté comprenant une combinaison des feuilles de verre dont l'indice de qualité optique est compris dans l'intervalle déterminé à l'étape (k). 5 La méthode de simulation de la puissance optique et la méthode de détermination de la qualité optique d'une combinaison de feuilles de verre peut être mise en oeuvre à l'aide d'un système de traitement de données comprenant des moyens permettant cette mise en oeuvre. Ce système de traitement de données peut comprendre un ordinateur intégrant une ou plusieurs unités centrales de traitement (Central Processing Unit) adaptées à la mise en oeuvre 1 o de la méthode de simulation de l'invention. L'ordinateur peut comprendre d'autres composants électroniques nécessaires à son fonctionnement, tels que des interfaces d'entrée-sortie, des systèmes de stockages volatiles et/ou rémanents et des BUS, nécessaires au transfert des données entre les unités centrales de traitement et à la communication avec des systèmes extérieurs. L'une des interfaces d'entrée-sortie peut être avantageusement une 15 interface utilisateur permettant des interactions humain-machine. En particulier, afin de faciliter la représentation et la communication des données obtenues à l'aide du système de l'invention, cette interface utilisateur peut être graphique. A titre d'exemple, un système de traitement de données permettant de simuler la puissance optique d'un verre feuillé de forme donnée susceptible d'être obtenu en feuilletant 20 au moins deux feuilles de verre et au moins un intercalaire de feuilletage, ledit verre feuilleté étant susceptible d'avoir un bord, une première face principale, et une deuxième face principale, peut comprendre les moyens suivants : a. un moyen de modélisation géométrique d'une première surface correspondant à la première face principale du verre feuilleté à partir d'un modèle numérique 25 de la forme dudit verre feuilleté pour former une première surface modélisée ; b. un moyen de modélisation géométrique d'une deuxième surface obtenue par translation de la surface modélisée dans l'étape a selon la direction de plus petite dimension du verre feuilleté d'une distance £ correspondant à l'épaisseur du vitrage, et selon la direction de l'une des deux plus grandes dimensions du 30 verre feuilleté d'une distance ô définie, le rapport de la distance ô sur la plus grande dimension du verre feuilleté étant supérieur à 0 et inférieur ou égal à 0,01, et ladite modélisation permettant d'obtenir une deuxième surface modélisée ; - 13 - c. un moyen de calcul du décalage géométrique -c, en chacun des points des surfaces modélisées, par soustraction de la première surface modélisée et de la deuxième surface modélisée ; d. un moyen de calcul de l'épaisseur géométrique E', en chacun des points des 5 surfaces modélisées, par addition de la distance £ et du décalage géométrique T; e. un moyen de calcul de l'épaisseur topologique 71", en chacun des points des surfaces modélisées, par soustraction des profils topographiques des surfaces des feuilles de verre susceptibles de correspondre à la première face principale 10 et deuxième face principale du verre feuilleté ; f. un moyen de calcul de la puissance optique PO, en chacun des points des surfaces modélisées, à partir de la valeur du Laplacien de la somme des épaisseurs géométrique El et topologique IT. 15 Dans un mode particulier de réalisation, le système de traitement des données peut comprendre en outre un moyen de détermination des indices de fluctuation optique (p. A titre d'exemple, un système de détermination d'un intervalle d'indices de qualité optique requis pour une combinaison de feuilles de verre susceptibles de former un verre feuilleté de forme donnée susceptible d'être obtenu en feuilletant au moins deux feuilles de 20 verre et au moins un intercalaire de feuilletage, ledit verre feuilleté ayant une contrainte donnée de qualité optique, 0, peut comprendre les moyens suivants : a. un moyen de sélection d'une combinaison de feuilles de verre parmi un ensemble de feuilles de verre auxquelles sont associés au moins un identifiant, un profil topographique et un indice de qualité optique ; 25 b. un moyen de calcul de l'indice de qualité optique de la combinaison de feuilles de verre à partir des indices de qualité optique des deux feuilles de verre de la combinaison, l'une des surfaces desdites deux feuilles de verre est susceptible de correspondre respectivement à la première face principale ou à la deuxième face principale du verre feuilleté ; 30 c. un système de traitement de données permettant de déterminer un indice de fluctuation optique cp d'un verre feuilleté susceptible de comprendre la combinaison des feuilles de verre, d. un moyen de comparaison des valeurs des indices de fluctuation optique cp à la valeur de la contrainte de qualité optique O; - 14 - e. un moyen de sélection de l'intervalle des indices de qualité optique des combinaisons des feuilles de verre satisfaisant le résultat de la comparaison. Dans un mode de réalisation, ledit système peut comprendre en outre un moyen de télécommunication directe ou indirecte avec un support de stockage déchiffrable par 5 ordinateur comprenant une base de données contenant pour chaque feuille de verre de l'ensemble des feuilles de verre, un identifiant, un profil topographique et un indice de qualité optique. Dans les deux exemples de système de traitement, les moyens de calcul, de modélisation géométrique, de comparaison et de sélection, peuvent être une ou plusieurs 10 unités de calculs. Ils peuvent aussi être des moyens virtuels instanciés sous la forme d'objets à partir de classes dans la mémoire vive d'un ordinateur, éventuellement assistée par une mémoire virtuelle, par un programme informatique ou un logiciel informatique. Un programme informatique de l'invention et/ou un système de traitement de données permettant de mettre oeuvre les étapes de la méthode de simulation et/ou de la 15 méthode de détermination d'un intervalle d'indices de qualité optique peuvent être intégrés dans une infrastructure informatique de type « cloud computing » ou infonuage. Un ordinateur « client » peut servir de plateforme d'accès à l'infonuage. Un tel mode de réalisation a, par exemple, pour avantage de permettre à plusieurs acteurs mettant en oeuvre la ou les méthodes de l'invention de mutualiser des moyens techniques et à centraliser l'information. Par exemple, 20 plusieurs sites de fabrication de verres feuilletés avec des formes et des contraintes de qualité optique différentes peuvent transmettre les informations pertinentes quant à ces formes et contraintes de qualité optique au travers de ladite infrastructure. Une fois les informations traitées par ledit programme ou le ou les dits systèmes, les résultats relatifs à la puissance optique ou à l'intervalle de qualité optique peuvent être transmis aux sites de fabrication au 25 travers de la même infrastructure. Un système de détermination d'un intervalle d'indices de qualité optique requis pour chacune des feuilles de verre d'une combinaison de feuilles de verre susceptibles de former un verre feuilleté comprenant au moins deux feuilles de verre et au moins un intercalaire, ledit verre feuilleté ayant une contrainte donnée de qualité optique, 0, peut en outre comprendre 30 une base données contenant pour chaque feuille de verre d'un ensemble de verre feuilleté, un identifiant, un profil topographique et un indice de qualité optique. Cette base de données peut également être comprise dans une infrastructure informatique de type « cloud computing » ou infonuage. - 15 - Ce mode de réalisation est particulièrement avantageux dans le cas où un fournisseur ou un fabriquant de feuilles de verre met à disposition dans une infrastructure informatique une base de données dans laquelle sont contenues des informations techniques sur les feuilles de verre qu'il est susceptible de fournir. Si les informations comprennent, pour chaque feuille 5 de verre, un indice de qualité optique et un profil topographique, un acteur ou un client, qui met en oeuvre l'invention, peut accéder à cette base de données et vérifier avant toute demande d'approvisionnement si ledit fournisseur ou ledit fabriquant est en capacité immédiate de fournir des feuilles de verre ayant une qualité optique suffisante pour la fabrication d'un verre feuilleté de forme donnée et une contrainte donnée de qualité optique. 1 o Les caractéristiques et les avantages de l'invention sont illustrés par les figures ci- dessous. La figure 1 est une représentation schématique d'un exemple de verre feuilleté avec une forme donnée. La figure 2 est une représentation schématique de la coupe transverse selon le plan A- 15 A du détail II du verre feuilleté de la figure 1. La figure 3 est une représentation sous forme de diagramme logique de la méthode de simulation de l'invention. La figure 4 est une représentation sous forme de diagramme logique d'un mode de réalisation de la méthode de détermination d'un intervalle d'indices de qualité optique requis 20 pour une combinaison de feuilles de verre susceptibles de former un verre feuilleté. La figure 5 est une représentation schématique d'un mode de réalisation d'un système de traitement de données comprenant des moyens de mettre en oeuvre les étapes d'une méthode de simulation de l'invention. La figure 6 est une modélisation géométrique de la forme du verre feuilleté de la figure 25 1. La figure 7 est une représentation graphique tridimensionnelle des exemples de profils topographiques des surfaces de feuilles de verre susceptible de correspondre à la première face principale et à la deuxième face principale d'un verre feuilleté. La figure 8 est une représentation schématique d'une délimitation d'une région du 30 verre feuilleté dans laquelle un indice de fluctuation cp est susceptible d'être calculé. La figure 9 est une représentation simplifiée de la première et deuxième surfaces modélisées du verre feuilleté avec la translation d'une distance ô et l'épaisseur géométrique E. La figure 10 est une représentation graphique de la puissance optique du verre feuilleté de la figure 1 obtenue à l'aide de la méthode du lancer de rayons. - 16 - La figure 11 est une représentation graphique de la puissance optique du verre feuilleté de la figure 1 obtenue à l'aide de la méthode de simulation de l'invention avec une valeur ô = 0,04. La figure 12 une représentation graphique de la puissance optique du verre feuilleté de 5 la figure 1 obtenue à l'aide de la méthode de simulation de l'invention avec une valeur ô = 0,004. La figure 13 est une représentation graphique de la variation des indices de fluctuation optique cp d'un verre feuilleté comprenant deux feuilles de verre et un intercalaire en fonction de l'indice de qualité optique NO des combinaisons de deux feuilles de verre choisies parmi un 1 o ensemble de 10 feuilles de verre, pour deux formes données 1 et 2 de verre feuilleté et une contrainte donnée de puissance optique, O. La figure 1 représente schématiquement, dans un repère orthonormé XYZ, un exemple de verre feuilleté 1000 avec une première face principale 1001, une deuxième face 15 principale 1002 et un bord 1003. La figure 2 est un schéma de la coupe transverse selon le plan A-A du détail II du verre feuilleté de la figure 1. A titre d'exemple, le verre feuilleté est formé d'une première feuille de verre 2001, une deuxième feuille de verre 2002 et un intercalaire de feuilletage 2003. La première face principale 1001, une deuxième face principale 1002 et le bord 1003 sont 20 également indiqué sur le schéma. La méthode de simulation de l'invention est représentée sous forme de diagramme logique sur la figure 3. Elle comprend les étapes suivantes : a. la modélisation géométrique E3001 d'une première surface susceptible de correspondre à la première face principale 1001 du verre feuilleté à partir d'un 25 modèle numérique de la forme dudit verre feuilleté pour former une première surface modélisée S1 ; b. la modélisation géométrique E3002 d'une deuxième surface obtenue par translation TE,8(Si) de la surface modélisée dans l'étape (a) selon la direction de plus petite dimension du verre feuilleté d'une distance £ correspondant à 30 l'épaisseur du vitrage, et selon la direction de l'une des deux plus grandes dimensions du verre feuilleté d'une distance ô définie, le rapport de la distance ô sur la plus grande dimension du verre feuilleté étant supérieur à 0 et inférieur ou égal à 0,01, et ladite modélisation permettant d'obtenir une deuxième surface modélisée S2 ; - 17 - c. le calcul E3003 du décalage géométrique -c, en chacun des points des surfaces modélisées, par soustraction de la première surface modélisée Si et de la deuxième surface modélisée S2 ; d. le calcul E3004 de l'épaisseur géométrique E', en chacun des points des 5 surfaces modélisées, par addition de la distance £ et du décalage géométrique T; e. le calcul E3005 de l'épaisseur topologique 7 1 - , en chacun des points des surfaces modélisées, obtenu par soustraction des profils topographiques PT1 et PT2 des surfaces de feuilles de verre susceptibles de correspondre à la première face 10 principale 1001 et à la deuxième face principale 1002 du verre feuilleté ; f. le calcul E3006 de la puissance optique PO, en chacun des points des surfaces modélisées, à partir de la valeur du Laplacien de la somme des épaisseurs géométrique El et topologique I T . 15 La figure 4 est une représentation sous forme de diagramme logique d'un mode de réalisation de la méthode de détermination, mis en oeuvre par ordinateur, d'un intervalle d'indices de qualité optique requis pour une combinaison de feuilles de verre susceptibles de former un verre feuilleté de forme donnée susceptible d'être obtenu en feuilletant au moins deux feuilles de verre et au moins un intercalaire de feuilletage, ledit verre feuilleté ayant une 20 contrainte donnée de qualité optique, O. Le verre feuilleté 1000 comprend une première face principale 1001, une deuxième face principale 1002 et un bord 1003. Pour une forme donnée de verre feuilleté, la méthode, mise en oeuvre par ordinateur, comprend les étapes suivantes : a. la modélisation géométrique E4000 d'une première surface susceptible de correspondre à la première face principale 1001 du verre feuilleté à partir d'un 25 modèle numérique de la forme dudit verre feuilleté pour former une première surface modélisée Si ; b. la modélisation géométrique E4001 d'une deuxième surface obtenue par translation TES (S1) de la surface modélisée dans l'étape (a) selon la direction de plus petite dimension du verre feuilleté d'une distance £ correspondant à 30 l'épaisseur du vitrage, et selon la direction de l'une des deux plus grandes dimensions du verre feuilleté d'une distance ô définie, le rapport de la distance ô sur la plus grande dimension du verre feuilleté étant supérieur à 0 et inférieur - 18 - ou égal à 0,01, ladite modélisation permettant d'obtenir une deuxième surface modélisée S2; c. le calcul E4002 du décalage géométrique -c, en chacun des points des surfaces modélisées, par soustraction de la première surface modélisée S1 et de la 5 deuxième surface modélisée S2 ; d. le calcul E4003 de l'épaisseur géométrique E', en chacun des points des surfaces modélisées, par addition de la distance £ et du décalage géométrique T; e. la génération E4004 d'un ensemble 3"k(E) de k feuilles de verre choisies 10 parmi un ensemble E de feuilles de verres auxquelles sont associés au moins un identifiant, un profil topographique et une indice de qualité optique ; f. l'initialisation 14005 d'une variable itérative 1; g. la sélection E4006 d'une combinaison C1 dans l'ensemble 3"k (E) ; h. le calcul E4007 l'épaisseur topologique i 1 , en chacun des points des surfaces 15 modélisées, obtenu par soustraction des profils topographiques PT1 et PT2 des surfaces de feuilles de verre susceptibles de correspondre à la première face principale 1001 et à la deuxième face principale 1002 du verre feuilleté susceptible de comprendre la combinaison C ; i. le calcul E4008 de la puissance optique PO, en chacun des points des surfaces 20 modélisées, à partir de la valeur du Laplacien de la somme des épaisseurs géométrique El et topologique I T ; j. le calcul E4008 de l'indice de fluctuation optique cp correspondant à la valeur maximale de l'ensemble des différences entre les valeurs maximales et minimales des puissances optiques PO dans une fenêtre d'échantillonnage 25 glissante ; k. le calcul E4009 de l'indice de qualité optique NO(ci) de la combinaison C1 de feuilles de verre à partir des indices de qualité optique des deux feuilles de verre de la combinaison, l'une des surfaces de chacune desdites deux feuilles de verre est susceptible de correspondre respectivement à la première face 30 principale ou à la deuxième face principale du verre feuilleté ; 1. l'incrémentation 140010 de la variable itérative 1 d'une unité ; la répétition des étapes (g) à (j) avec, à chaque itération, une nouvelle combinaison C1 de feuilles de verre ; la variable itérative 1 étant incrémentée d'une unité jusqu'à ce qu'elle - 19 - devienne égale à la valeur du coefficient binomial Ccard(E)k avec Card(E) la cardinalité de l'ensemble E de feuilles de verres ; m. la construction E4011 d'un ensemble cri de valeurs d'indices de fluctuation optique cp pour chaque combinaison c1 dans l'ensemble 3"k (E) ; 5 n. la sélection E4012 de l'intervalle 1 des indices de qualité optique 1O(c1) des combinaisons c1 de feuilles de verre satisfaisant le résultat de la comparaison de de l'étape (1). La figure 5 est une représentation schématique d'un mode de réalisation d'un système 10 de traitement de données comprenant des moyens de mettre en oeuvre les étapes d'un mode de réalisation de la méthode de détermination d'un intervalle d'indices de qualité optique requis pour une combinaison de feuilles de verre susceptibles de former un verre feuilleté de forme donnée susceptible d'être obtenu en feuilletant au moins deux feuilles de verre et au moins un intercalaire de feuilletage, ledit verre feuilleté ayant une contrainte donnée de qualité 15 optique, O. Le système comprend un ordinateur 5001 en télécommunication avec une infrastructure informatique 5003 de type « cloud computing » ou infonuage. Les moyens sont intégrés dans l'infrastructure informatique 5003. La dite infrastructure informatique 5003 comprend les moyens suivants : a. un moyen de sélection 5003a d'une combinaison de feuilles de verre parmi un 20 ensemble de feuilles de verre auxquelles sont associés au moins un identifiant, un profil topographique et un indice de qualité optique ; b. un moyen de calcul E5003b de l'indice de qualité optique de la combinaison de feuilles de verre à partir des indices de qualité optique des deux feuilles de verre de la combinaison, l'une des surfaces de chacune desdites deux feuilles 25 de verre est susceptible de correspondre respectivement à la première face principale ou à la deuxième face principale du verre feuilleté ; c. un système de traitement de données 5003c détermination d'un indice de fluctuation optique cp d'un verre feuilleté susceptible de comprendre la combinaison des feuilles de verre ; 30 d. un moyen de comparaison 5003d des valeurs des indices de fluctuation optique cp à la valeur de la contrainte de qualité optique O; e. un moyen de sélection 5003e de l'intervalle des indices de qualité optique des combinaisons des feuilles de verre satisfaisant le résultat de la comparaison. - 20 - L'infrastructure informatique comprend en outre une base de données 5003f contenant pour chaque feuille de verre d'un ensemble de verre feuilleté, un identifiant, un profil topographique et un indice de qualité optique. L'ordinateur comprend une interface graphique d'entrée-sortie 5002 permettant des interactions humain-machine. 5 Un exemple de modélisation 6000 de la forme du verre feuilleté de la figure 1 est représenté sur la figure 6. La modélisation a été obtenue par maillage de la surface. Une première surface modélisée 6001 est obtenue par la modélisation géométrique de la première face principale 1001 du verre feuilleté. La figure 7 est une représentation graphique tridimensionnelle d'exemples de profils 10 topographiques PT1 et PT2 des surfaces des feuilles de verre susceptibles de correspondre respectivement à la première face principale 1001 et à la deuxième face principale 1002 d'un verre feuilleté. Plus particulièrement, elle représente la variation en um du relief des surfaces des feuilles de verre. Sur la figure, l'axe Z est l'échelle de la variation du relief et les axes X et Y sont dans le plan des feuilles de verre. Pour des raisons de concision, seuls les reliefs d'une 15 région des surfaces sont représentés. La figure 8 est une représentation schématique d'une délimitation d'une région 8001 du verre feuilleté 1000 dans laquelle un indice de fluctuation cp est susceptible d'être calculé. En particulier, cette région est délimitée par l'intersection des faces principales dudit verre feuilleté 1000 avec une pyramide à base carrée (non représenté) dont l'apex est situé dans la 20 direction de l'axe Z à une distance définie des faces principales 1001 et 1002, les deux angles au sommet entre les faces latérales en vis-à-vis de ladite pyramide étant respectivement compris entre 100 et 20°, et entre 5 et 15°. La forme de la région est généralement rectangulaire ou carrée. La figure 9 est une représentation simplifiée de la deuxième surface modélisée 9001 25 du verre feuilleté obtenue par translation de la première surface modélisée 6001 selon la direction de plus petite dimension du verre feuilleté d'une distance £ correspondant à l'épaisseur du vitrage, et selon la direction de l'une des deux plus grandes dimensions du verre feuilleté d'une distance ô définie, le rapport de la distance ô sur la plus grande dimension du verre feuilleté étant supérieur à 0 et inférieur ou égal à 0,01. Dans le repère orthonormé XYZ 30 de la figure, la translation d'une distance ô est réalisée selon la direction des axes X ou Y ou selon une direction combinaison linéaire des axes X et Y. La translation de la distance £ est réalisée selon la direction de l'axe Z. La figure 10 est une représentation graphique de la puissance optique d'une région du verre feuilleté de la figure 1 obtenue à l'aide la méthode du lancer de rayons. L'échelle de - 21 - niveaux de gris représente la puissance optique exprimée en milidioptries (mD). La région du verre feuilleté dans laquelle la puissance optique a été calculée est une région du type de celle de la figure 8. Les axes X et Y de la figure correspondent aux axes de la figure 1. L'unité des axes et la position du centre du repère sont arbitraires. La méthode du lancer de rayons est une 5 méthode connue de l'état de la technique, fréquemment utilisée dans le domaine technique. La figure 11 est une représentation graphique de la puissance optique d'une région du verre feuilleté de la figure 1 obtenue à l'aide de la méthode de simulation de l'invention avec une valeur ô = 0,04. Pour la simulation, les profils topographiques des feuilles de verre du type de ceux de la figure 7 ont été utilisés. L'échelle de niveaux de gris représente la puissance 1 o optique exprimée en milidioptrie (mD). La région du verre feuilleté dans laquelle la puissance optique a été calculée est une région du type de celle de la figure 8. Les axes X et Y de la figure correspondent aux axes de la figure 1. L'unité des axes et la position du centre du repère sont arbitraires. La figure 12 est une représentation graphique de la puissance optique d'une région du 15 verre feuilleté de la figure 1 obtenue à l'aide de la méthode de simulation de l'invention avec une valeur ô = 0,004. Pour la simulation, les profils topographiques des feuilles de verre du type de ceux de la figure 7 ont été utilisés. L'échelle de niveaux de gris représente la puissance optique exprimée en milidioptries (mD). La région du verre feuilleté dans laquelle la puissance optique a été calculée est une région du type de celle de la figure 8. Les axes X et Y de la figure 20 correspondent aux axes de la figure 1. L'unité des axes et la position du centre du repère sont arbitraires. La figure 13 est une représentation graphique de la variation des indices de fluctuation optique cp d'un verre feuilleté susceptible de comprendre deux feuilles de verre et un intercalaire en fonction de l'indice de qualité optique NO des combinaisons de deux feuilles de 25 verre choisies parmi un ensemble de 10 feuilles de verre. Les indices de fluctuation optique cp ont été calculés pour deux formes données 1 et 2 du verre feuilleté et une contrainte donnée de qualité optique O. Le nombre de combinaisons pour chacune des deux formes est de 45. Cette représentation graphique a été obtenue à l'aide de la méthode de l'invention permettant de déterminer l'intervalle d'indices de qualité optique requis pour une combinaison de feuilles 30 de verre susceptibles de former un verre feuilleté comprenant au moins deux feuilles de verre et au moins un intercalaire, ledit verre feuilleté ayant une contrainte donnée de puissance optique, O. L'indice de qualité optique NO des combinaisons des deux feuilles de verre est la moyenne quadratique des indices moyens de qualité optique des deux feuilles de verre de la combinaison. Lesdits indices moyens de qualité optiques ont été obtenus par des méthodes - 22 - ombroscopiques. Sur la figure la contrainte de qualité optique, 0, est représentée par la ligne en pointillés. Les figurés en forme de ronds pleins représentent les résultats obtenus pour la forme 1 du verre feuilleté ; ceux en forme de carrés vides représentent les résultats obtenus pour la forme 2 du verre feuilleté. 5 Exemple 1 Le premier exemple compare les valeurs de puissance optique d'un verre feuilleté à l'aide de la méthode de simulation de l'invention et à l'aide de la méthode du lancer de rayons, connue de l'état de la technique et fréquemment utilisée dans le domaine technique. 1 o Le verre feuilleté comprend deux feuilles de verre sodo-calcique et un intercalaire de feuilletage. Sa forme et sa structure sont représentées schématiquement sur les figures 1 et 2 respectivement. Le verre feuilleté 1000 comprend une première feuille de verre 2001, une deuxième feuille de verre 2002 et un intercalaire de feuilletage 2003 de manière à former une première face principale 1001, une deuxième face principale 1002 et un bord 1003. 15 L'épaisseur des feuilles de verre est de 2,1mm. L'indice de réfraction du verre est 1,51. Des extraits des profils topographiques PT1 et PT2 des surfaces des feuilles de verre correspondant à la première face principale 1001 et à la deuxième face principale 1002 du verre feuilleté respectivement sont représentés sur le figure 7. Les profils topographiques ont été acquis par ombroscopie. 20 Les deux méthodes ont été mises en oeuvre à l'aide d'un système de traitement de données comprenant un ordinateur de performances standards. La figure 10 est une représentation graphique de la puissance optique d'une région du verre feuilleté obtenue à l'aide la méthode du lancer de rayons. La figure 11 est la représentation graphique de la puissance optique d'une région du verre feuilleté obtenue à 25 l'aide de la méthode de simulation de l'invention avec une valeur ô = 0,04. La figure 12 représente la puissance optique obtenue avec une valeur ô = 0,004. L'échelle en niveaux de gris représente la puissance optique exprimée en milidioptries (mD). Sur les trois figures, l'unité des axes X et Y, et la position du centre du repère sont arbitraires. La comparaison de figures 10 et 11 montre que les résultats obtenus à l'aide de la 30 méthode de simulation de l'invention avec une valeur ô supérieure à 0,01 sont différents de ceux obtenus avec la méthode du lancer de rayon. Il a, en outre, été constaté expérimentalement que lorsque la valeur de ô est supérieure à 0,01 les résultats sont incorrects. - 23 - La comparaison de figures 10 et 12 montre que les résultats obtenus à l'aide de la méthode de simulation de l'invention avec une valeur ô inférieure à 0,01 sont identiques, sinon très proches, de ceux obtenus avec la méthode du lancer de rayon. En revanche, la méthode de simulation de l'invention est, en terme vitesse d'exécution, environ 50 fois plus 5 rapide à mettre en oeuvre que la méthode du lancer de rayons. Cette rapidité de mise en oeuvre rend possible l'utilisation en temps réel de la méthode sur des lignes de production de verre feuilleté où des changements dans le planning de production sont susceptibles d'avoir lieu. Exemple 2 1 o Le second exemple illustre les avantages de la méthode de simulation de l'invention lorsqu'elle est utilisée dans une méthode d'un intervalle d'indices de qualité optique requis pour une combinaison de feuilles de verre susceptibles de former un verre feuilleté de forme donnée susceptible d'être obtenu en feuilletant au moins deux feuilles de verre et au moins un intercalaire de feuilletage, ledit verre feuilleté ayant une contrainte donnée de qualité optique, 15 0. L'intervalle d'indices de qualité optique requis pour une combinaison de deux feuilles de verre susceptibles de former un verre feuilleté comprenant deux feuilles de verre (Verre 1 et Verre 2) et un intercalaire de feuilletage a été déterminé pour deux formes 1 et 2 de verre feuilleté. La combinaison de deux feuilles de verre est choisie parmi un ensemble de 10 verres. 20 L'indice de qualité optique NO des combinaisons de deux feuilles de verre est la moyenne quadratique des indices moyens de qualité optique des deux feuilles de verre de la combinaison. Lesdits indices moyens de qualité optiques des feuilles de verre ont été obtenus par des méthodes ombroscopiques. Le tableau ci-dessous regroupe les 45 combinaisons évaluées. Dans la première ligne 25 (Verre 1) et la première colonne (Verre 2) du tableau sont indiquées en gras les valeurs des indices de qualité optique de chacun des dix verres. Chaque intersection des lignes et colonnes représente une combinaison de deux verres Verre 1 et Verre 2. La valeur à chaque intersection est la valeur de l'indice de qualité NO de la combinaison, moyenne quadratique des indices de qualité optique du Verre 1 et du Verre 2. 30 - 24 - Verre 1 9,54 11,32 12,97 14,06 15,29 16,31 18,23 19 20,14 21,86 Verre 2 9,54 11,32 10,47 - 12,97 11,38 12,17 - 14,06 12,01 12,76 13,53 15,29 12,74 13,45 14,18 14,69 16,31 13,36 14,04 14,73 15,23 15,81 18,23 14,55 15,17 15,82 16,28 16,82 17,30 19 15,03 15,64 16,27 16,71 17,25 17,71 18,62 20,14 15,76 16,34 16,94 17,37 17,88 18,33 19,21 19,58 - 21,86 16,87 17,41 17,97 18,38 18,86 19,29 20,13 20,48 21,02 - Pour chaque combinaison et chacune des deux forme du verre feuilleté, un indice de fluctuation optique cp a été calculé. Dans l'exemple, l'indice de fluctuation optique cp est défini comme la valeur maximale d'un ensemble de valeurs correspondant aux différences entre les 5 valeurs maximales et minimales des puissances optiques dans une fenêtre d'échantillonnage glissante carrée. La dimension latérale de la fenêtre d'échantillonnage est 80mm. Le calcul des indices de fluctuation optique a été réalisé dans une région du verre feuilleté. Cette région est la région délimitée par l'intersection des faces principales dudit verre feuilleté avec une pyramide à base carrée dont l'apex est situé à une distance définie des faces principales et les 1 o deux angles au sommet entre les faces latérales en vis-à-vis de ladite pyramide sont respectivement compris entre 100 et 20°, et entre 5 et 15°. Les dimensions latérales de cette région sur le verre feuilleté sont 200mm x 1000mm. L'unité de cp est le milidioptrie (mD). Les résultats sont représentés sur la figure 13. La contrainte de puissance optique, 0, est représentée par la ligne en pointillés. Elle est définie à 63. Les figurés en forme de ronds 15 pleins représentent les résultats obtenus pour la forme 1 du verre feuilleté ; ceux en forme de carrés vides représentent les résultats obtenus pour la forme 2 du verre feuilleté. Pour la forme 1, la contrainte de qualité optique 0 est satisfaite lorsque la valeur de l'indice de qualité optique NO de la combinaison de deux feuilles de verre est inférieure à 13 ; pour la forme 2 lorsqu'elle est inférieure à 17. 20 La figure 13 montre clairement que ladite méthode permet de déterminer un intervalle d'indices REVENDICATIONS1. Méthode de simulation, mise en oeuvre par ordinateur, de la puissance optique d'un verre feuilleté de forme donnée susceptible d'être obtenu en feuilletant au moins deux feuilles de verre et au moins un intercalaire de feuilletage, ledit verre feuilleté étant susceptible d'avoir un bord, une première face principale, et une deuxième face principale, et ladite méthode comprenant les étapes suivantes : al. la modélisation géométrique d'une première surface susceptible de correspondre à la première face principale du verre feuilleté à partir d'un modèle numérique de la forme dudit verre feuilleté pour former une première surface modélisée ; bl. la modélisation géométrique d'une deuxième surface obtenue par translation de la surface modélisée dans l'étape (al) selon la direction de plus petite dimension du verre feuilleté d'une distance e correspondant à l'épaisseur du verre feuilleté, et selon la direction de l'une des deux plus grandes dimensions du verre feuilleté d'une distance 8 définie, le rapport de la distance 8 sur la plus grande dimension du verre feuilleté étant supérieur à 0 et inférieur ou égal à 0,01, et ladite modélisation permettant d'obtenir une deuxième surface modélisée ; cl. le calcul du décalage géométrique r, en chacun des points des surfaces modélisées, par soustraction de la première surface modélisée et de la deuxième surface modélisée ; dl. le calcul de l'épaisseur géométrique Et , en chacun des points des surfaces modélisées, par addition de la distance E et du décalage géométrique r; el. le calcul de l'épaisseur topologique ir, en chacun des points des surfaces modélisées, par soustraction des profils topographiques des surfaces des feuilles de verre susceptibles de correspondre à la première face principale et deuxième face principale du verre feuilleté ; fi. le calcul de la puissance optique PO, en chacun des points des surfaces modélisées, à partir de la valeur du Laplacien de la somme des épaisseurs géométrique Et et topologique m. 2. Méthode de simulation selon la revendication 1, telle que les épaisseurs géométrique EI et topologique Tt sont chacune pondérées par un facteur d'amplification optique dépendant de l'angle, 0, formé entre la direction normale à la première surface et une direction de référence d'observation du verre feuilleté. - 26 - 3. Méthode de simulation selon l'une quelconque des revendications 1 à 2, telle que les profils topographiques des surfaces de feuilles de verre susceptibles de correspondre à la première face principale et à la deuxième face principale du verre feuilleté sont des profils mesurés par méthodes profilométriques avec ou sans contact. 5 4. Méthode de simulation selon l'une quelconque des revendications 1 à 2, telle que les profils topographiques des surfaces de feuilles de verre susceptibles de correspondre à la première face principale et à la deuxième face principale du verre feuilleté sont des profils simulés numériquement. 5. Méthode de simulation selon l'une quelconque des revendications 1 à 4, telle qu'elle 10 comprend en outre, après l'étape (f1) une étape de détermination d'un indice de fluctuation optique cp à partir des valeurs des puissances optiques PO obtenues en chacun des points des surfaces modélisées. 6. Méthode de simulation selon la revendication 5, telle que l'indice de fluctuation optique cp est l'écart-type des valeurs des puissances optiques PO ou à la valeur 15 maximales des puissances optiques. 7. Méthode de simulation selon la revendication 5, telle que l'indice de fluctuation optique cp est la valeur maximale d'un ensemble de valeurs correspondant aux différences entre les valeurs maximales et minimales des puissances optiques PO dans une fenêtre d'échantillonnage glissante. 20 8. Méthode de simulation selon la revendication 7, telle que la fenêtre d'échantillonnage glissante est une fenêtre d'échantillonnage carrée ou rectangulaire dont les dimensions latérales sont comprises entre 2mm et 100mm, de préférence entre 5 et 50mm. 9. Méthode de simulation selon l'une quelconque des revendications 5 à 8, telle que la détermination de l'indice de fluctuation optique (t:i est réalisée dans une région limitée 25 des surfaces modélisées. 10. Méthode de simulation selon la revendication 9, telle que la région limitée du verre feuilleté est la région délimitée par l'intersection des faces principales dudit verre feuilleté avec une pyramide à base carrée dont l'apex est situé à une distance définie des faces principales et les deux angles au sommet entre les faces latérales en vis-à-vis 30 de ladite pyramide sont respectivement compris entre 100 et 20°, et entre 5 et 150 . 11. Méthode de détermination, mise en oeuvre par ordinateur, d'un intervalle d'indices de qmlité optique requis pour une combinaison de feuilles de verre susceptibles de former un verre feuilleté de forme donnée susceptible d'être obtenu en feuilletant au moins deux feuilles de verre et au moins un intercalaire de feuilletage, ledit verre - 27 - feuilleté ayant une contrainte donnée de qualité optique, 0, et ladite méthode comprenant les étapes suivantes : a2. la sélection d'une combinaison de feuilles de verre parmi un ensemble de feuilles de verre auxquelles sont associés au moins un identifiant, un profil topographique 5 et un indice de qualité optique ; b2. le calcul de l'indice de qualité optique de la combinaison de feuilles de verre à partir des indices de qualité optique des deux feuilles de verre de la combinaison, l'une des surfaces de chacune desdites deux feuilles de verre est susceptible de correspondre respectivement à la première face principale ou à la deuxième face principale du verre 10 feuilleté ; c2. la détermination d'un indice de fluctuation optique (I) d'un verre feuilleté susceptible de comprendre la combinaison des feuilles de verre à l'aide d'une méthode de simulation selon l'une quelconque des revendications 5 à 10; d2. la répétition des étapes (a2) à (c2) avec une nouvelle combinaison de feuilles 15 de verre de manière à obtenir un ensemble de valeurs d'indices de fluctuation optique cp pour les différentes combinaisons possibles ; e2. la comparaison des valeurs des indices de fluctuation optique, cp, à la valeur de la contrainte de qualité optique O; f2. la sélection de l'intervalle des indices de qualité optique des combinaisons des 20 feuilles de verre satisfaisant le résultat de la comparaison de l'étape (e2). 12. Méthode selon la revendication 11, telle que l'indice de qualité optique des combinaisons des feuilles de verre est la moyenne quadratique des indices moyens de qualité optique des deux feuilles de verre de la combinaison dont l'une des surfaces est susceptible de correspondre respectivement à la première face principale et à la 25 deuxième face principale du verre feuilleté, lesdits indices moyens de qualité optiques étant obtenus par des méthodes ombroscopiques. 13. Programme informatique comprenant des instructions exécutables par un ordinateur permettant l'exécution des étapes d'une méthode selon l'une quelconque des revendications 1 à 12. 30 14. Support de stockage déchiffrable par un ordinateur comprenant un programme informatique selon la revendication 13. 15. Système de traitement de données comprenant des moyens de mettre en oeuvre les étapes d'une méthode selon l'une quelconque des revendications 1 à 12. - 28 - 16. Procédé de fabrication d'un verre feuilleté comprenant au moins deux feuilles de verre et au moins intercalaire, ledit verre feuilleté de forme donnée ayant un bord, une première face principale et une deuxième face principale, ledit verre feuilleté ayant une contrainte donnée de qualité optique, 0, ledit procédé comprend les étapes suivantes : 5 a. la sélection, mise en oeuvre par ordinateur, d'une combinaison de feuilles de verre parmi un ensemble de feuilles de verre auxquelles sont associés au moins un identifiant, un profil topographique et un indice de qualité optique, ladite combinaison de feuilles de verre étant susceptible de former ledit verre feuilleté ; 10 b. le calcul, mis en oeuvre par ordinateur, de l'indice de qualité optique de la combinaison de feuilles de verre à partir des indices de qualité optique des deux feuilles de verre de la combinaison, l'une des surfaces de chacune desdites deux feuilles de verre est susceptible de correspondre respectivement à la première face principale ou à la deuxième face principale du verre feuilleté ; 15 c. la modélisation géométrique, mise en oeuvre par ordinateur, d'une première surface correspondant à la première face principale du verre feuilleté à partir d'un modèle numérique dudit verre feuilleté pour former une première surface modélisée ; d. la modélisation géométrique, mise en oeuvre par ordinateur, d'une deuxième 20 surface obtenue par translation de la surface modélisée dans l'étape a selon la direction de plus petite dimension du verre feuilleté d'une distance £ correspondant à l'épaisseur du verre feuilleté, et selon la direction de l'une des deux plus grandes dimensions du verre feuilleté d'une distance 8 définie, le rapport de la distance 5 sur la plus grande dimension du verre feuilleté étant 25 supérieur à 0 et inférieur ou égal à 0,01, et ladite modélisation permettant d'obtenir une deuxième surface modélisée ; e. le calcul, mis en oeuvre par ordinateur, du décalage géométrique T obtenu par soustraction de la première surface modélisée et de la deuxième surface modélisée en chacun des points des surfaces ; 30 f. le calcul, mis en oeuvre par ordinateur, de l'épaisseur géométrique el par addition de la distance e et du décalage géométrique T en chacun des points des surfaces ; g. le calcul, mis en oeuvre par ordinateur, de l'épaisseur topologique Tt obtenu par soustraction des profils topographiques des surfaces des feuilles de verre - 29 - susceptibles de correspondre à la première face principale et deuxième face principale du verre feuilleté en chacun des points des faces ; h. le calcul, mis en oeuvre par ordinateur, de la puissance optique PO à partir de la valeur du Laplacien de la somme des épaisseurs géométrique E' et topologique 5 if en chacun des points des surfaces modélisées.. i. le calcul, mis en oeuvre par ordinateur, d'un indice de fluctuation optique cp correspondant à la valeur maximale d'un ensemble de valeurs correspondant aux valeurs maximales des différences entre les valeurs maximales et minimales des puissances optiques PO dans une fenêtre d'échantillonnage glissante ; 10 j. la répétition des étapes (a) à (i) avec une nouvelle combinaison de feuilles de verre de manière à obtenir un ensemble de valeurs d'indices de fluctuation optique cp pour chaque combinaison possible ; k. la comparaison, mise en oeuvre par ordinateur, des valeurs des indices de fluctuation optique, cp, à la valeur de la contrainte de qualité optique ; 15 1. la sélection, mise en oeuvre par ordinateur, de l'intervalle des indices de qualité optique des combinaisons des feuilles de verre satisfaisant le résultat de la comparaison de de l'étape (j) ; m. la fabrication d'un verre feuilleté comprenant une combinaison des feuilles de verre dont l'indice de qualité optique est compris dans l'intervalle déterminé à 20 l'étape (k).