FR 3068136 A1 20181228 FR 1755709 A 20170622 OBTENTION D'UN VITRAGE RESISTANT A UN IMPACT La présente invention a trait à un procédé non destructif de validation qu'un vitrage installé dans une structure résiste à un impact dynamique, tel qu'un impact à l'oiseau pour un vitrage d'aéronef. L'invention a également trait à un procédé de fabrication d'un vitrage pour qu'il résiste à un impact dynamique. Au sens de l'invention, on entend par "vitrage" une structure comprenant au moins une feuille de verre. De plus, on entend par "vitrage feuilleté" une structure de vitrage comprenant un empilement d'au moins une feuille de verre et un intercalaire polymérique, y compris un empilement d'une seule feuille de verre et d'un seul intercalaire polymérique assemblés entre eux. Dans le cadre de l'invention, une feuille de verre est un substrat transparent rigide, qui peut être en verre minéral ou en verre organique. Un intercalaire polymérique est un feuillet intercalaire à base de matériau polymère, notamment thermoformable ou sensible à la pression, c'est-à-dire le type de feuillet qu'on utilise comme intercalaire dans les vitrages feuilletés. Il peut s'agir d'un intercalaire monolithique ou d'un intercalaire composite constitué par l'assemblage de plusieurs éléments polymériques sous forme de couches, résines ou films. Il est connu de tester la résistance à l'impact des vitrages de véhicules par des tests destructifs normalisés avec des impacteurs représentatifs de situations réelles, tels qu'un test à l'oiseau dans le domaine aéronautique pour des vitrages d'avions ou d'hélicoptères, un test au pavé pour des vitrages de trains, un test balistique pour des vitrages de véhicules blindés. De tels tests destructifs sont coûteux, notamment en ce qu'ils nécessitent de la main d'ceuvre et la mise au rebut des vitrages testés en cas de casse. De plus, ils ne permettent pas d'optimiser la composition et l'intégration des vitrages dans une structure, notamment en fonction de la raideur des matériaux du vitrage ou des systèmes d'attache. Il en résulte souvent un surdimensionnement des vitrages, c'est-à-dire des épaisseurs trop importantes des feuilles de verre, et éventuellement des intercalaires polymériques dans le cas de vitrages feuilletés, par rapport aux niveaux de contraintes mécaniques susceptibles d'être appliqués aux vitrages. Dès lors, le coût et la masse des vitrages ne sont pas optimisés. C'est à ces inconvénients qu'entend plus particulièrement remédier l'invention en proposant un procédé non destructif de validation qu'un vitrage installé dans une structure résiste à un impact dynamique, et un procédé de fabrication d'un vitrage garantissant l'obtention d'un vitrage optimisé à la fois en termes de masse, de coût et de résistance à un impact. A cet effet, l'invention a pour objet un procédé non destructif de validation 5 qu'un vitrage installé dans une structure résiste à un impact dynamique, tel qu'un impact à l'oiseau pour un vitrage installé dans un aéronef, le vitrage comprenant au moins une feuille de verre, caractérisé en ce qu'il comprend des étapes dans lesquelles : - à l'aide d'un modèle numérique par éléments finis du vitrage installé dans 10 la structure et soumis à l'impact, en utilisant une loi de comportement du matériau constitutif de chaque feuille de verre, on calcule l'enveloppe de contrainte maximale sur au moins une feuille de verre critique du vitrage (de préférence on calcule l'enveloppe de contrainte maximale sur chaque feuille de verre du vitrage) ; 15 - pour au moins la feuille de verre critique du vitrage (de préférence pour chaque feuille de verre du vitrage), on compare l'enveloppe de contrainte maximale à la contrainte à rupture de la feuille de verre obtenue selon une méthode expérimentale de détermination de la contrainte à rupture d'une feuille de verre sélectionnée en fonction du type d'impact, en particulier on vérifie si 20 l'enveloppe de contrainte maximale est strictement inférieure à la contrainte à rupture. Selon un aspect de l'invention, le procédé est un procédé non destructif de validation qu'un vitrage feuilleté installé dans une structure résiste à un impact dynamique, le vitrage feuilleté comprenant un empilement d'au moins une feuille 25 de verre et un intercalaire polymérique, le procédé comprenant des étapes dans lesquelles : - à l'aide d'un modèle numérique par éléments finis du vitrage feuilleté installé dans la structure et soumis à l'impact, en utilisant une loi de comportement du matériau constitutif de chaque feuille de verre, une loi de comportement du 30 matériau constitutif de chaque intercalaire polymérique, et une loi de comportement de chaque interface entre une feuille de verre et un intercalaire polymérique, on calcule l'enveloppe de contrainte maximale sur au moins une feuille de verre critique du vitrage feuilleté (de préférence on calcule l'enveloppe de contrainte maximale sur chaque feuille de verre du vitrage feuilleté) ; - pour au moins la feuille de verre critique du vitrage feuilleté (de préférence pour chaque feuille de verre du vitrage feuilleté), on compare l'enveloppe de contrainte maximale à la contrainte à rupture de la feuille de verre obtenue selon une méthode expérimentale de détermination de la contrainte à rupture d'une 5 feuille de verre sélectionnée en fonction du type d'impact, en particulier on vérifie si l'enveloppe de contrainte maximale est strictement inférieure à la contrainte à rupture. L'invention a également pour objet un procédé de fabrication d'un vitrage pour qu'il résiste à un impact dynamique lorsqu'il est installé dans une structure, 10 tel qu'un impact à l'oiseau pour un vitrage installé dans un aéronef, le vitrage comprenant au moins une feuille de verre, caractérisé en ce qu'il comprend des étapes dans lesquelles : - à l'aide d'un modèle numérique par éléments finis du vitrage installé dans la structure et soumis à l'impact, en utilisant une loi de comportement du matériau 15 constitutif de chaque feuille de verre, on calcule l'enveloppe de contrainte maximale sur au moins une feuille de verre critique du vitrage (de préférence on calcule l'enveloppe de contrainte maximale sur chaque feuille de verre du vitrage), en fonction des dimensions du vitrage ; - on ajuste les caractéristiques du vitrage parmi ses dimensions, le matériau 20 constitutif de chaque feuille de verre, de telle sorte que l'enveloppe de contrainte maximale calculée pour au moins la feuille de verre critique du vitrage (de préférence pour chaque feuille de verre du vitrage) est strictement inférieure à la contrainte à rupture de la feuille de verre obtenue selon une méthode expérimentale de détermination de la contrainte à rupture d'une feuille de verre 25 sélectionnée en fonction du type d'impact, tout en ayant un dimensionnement optimisé du vitrage ; - on prépare et on assemble chaque feuille de verre du vitrage avec les caractéristiques ajustées. Selon un aspect de l'invention, le procédé est un procédé de fabrication 30 d'un vitrage feuilleté pour qu'il résiste à un impact dynamique lorsqu'il est installé dans une structure, le vitrage feuilleté comprenant un empilement d'au moins une feuille de verre et un intercalaire polymérique, le procédé comprenant des étapes dans lesquelles : - à l'aide d'un modèle numérique par éléments finis du vitrage feuilleté installé dans la structure et soumis à l'impact, en utilisant une loi de comportement du matériau constitutif de chaque feuille de verre, une loi de comportement du matériau constitutif de chaque intercalaire polymérique, et une loi de 5 comportement de chaque interface entre une feuille de verre et un intercalaire polymérique, on calcule l'enveloppe de contrainte maximale sur au moins une feuille de verre critique du vitrage feuilleté (de préférence on calcule l'enveloppe de contrainte maximale sur chaque feuille de verre du vitrage feuilleté), en fonction des dimensions du vitrage feuilleté ; 10 - on ajuste les caractéristiques du vitrage feuilleté parmi ses dimensions, le matériau constitutif de chaque feuille de verre, le matériau constitutif de chaque intercalaire polymérique, la nature de chaque interface entre une feuille de verre et un intercalaire polymérique, de telle sorte que l'enveloppe de contrainte maximale calculée pour au moins la feuille de verre critique du vitrage feuilleté (de 15 préférence pour chaque feuille de verre du vitrage feuilleté) est strictement inférieure à la contrainte à rupture de la feuille de verre obtenue selon une méthode expérimentale de détermination de la contrainte à rupture d'une feuille de verre sélectionnée en fonction du type d'impact, tout en ayant un dimensionnement optimisé du vitrage feuilleté ; 20 - on prépare et on assemble chaque feuille de verre et chaque intercalaire polymérique du vitrage feuilleté avec les caractéristiques ajustées. Dans le cadre de l'invention, on entend par "feuille de verre critique" une feuille de verre du vitrage dont on sait qu'elle est la plus susceptible de casser lors de l'impact dynamique, par exemple en raison de son épaisseur, de son matériau 25 constitutif, de son positionnement vis-à-vis de l'impacteur, etc. L'invention peut alors être mise en oeuvre en calculant l'enveloppe de contrainte maximale sur cette feuille de verre critique uniquement. Toutefois, dans un mode de réalisation préféré de l'invention, l'enveloppe de contrainte maximale est calculée sur chaque feuille de verre du vitrage. 30 Au sens de l'invention, on entend par "dimensions du vitrage" non seulement ses dimensions périphériques, par exemple dans le cas d'un vitrage rectangulaire, sa largeur et sa longueur, mais également les épaisseurs de sa ou ses feuilles de verre et, dans le cas d'un vitrage feuilleté, de son ou ses intercalaires polymériques constitutifs. De plus, on entend par "dimensionnement optimisé du vitrage" le fait d'avoir une épaisseur de verre et éventuellement une épaisseur d'intercalaire polymérique dans le vitrage ajustées pour minimiser la masse et/ou le coût du vitrage. Dans le cadre de l'invention, la ou chaque feuille de verre du vitrage peut 5 être une feuille en verre minéral, notamment un verre d'oxyde tel qu'un silicate, borate, sulfate, phosphate, ou autre. Chaque feuille de verre du vitrage qui est en verre minéral est avantageusement renforcée, notamment par un procédé de trempe thermique ou par un procédé d'échange ionique également appelé "trempe chimique". 10 De manière connue, les procédés de trempe thermique et de trempe chimique permettent d'augmenter la résistance aux chocs des feuilles de verre minéral, en créant une zone superficielle en compression et une zone centrale en tension. Dans le cas de la trempe chimique, la substitution superficielle d'un ion de la feuille de verre (généralement un ion alcalin tel que le sodium ou le lithium) par 15 un ion de rayon ionique plus grand (généralement un ion alcalin, tel que le potassium ou le sodium) permet de créer en surface de la feuille de verre des contraintes résiduelles de compression, jusqu'à une certaine profondeur. Dans l'ensemble du texte, une profondeur correspond, selon une coupe transversale, à une distance entre un point considéré et une surface de la feuille de verre, 20 mesurée selon une normale à ladite surface. Les contraintes sont parallèles à la surface de la feuille de verre, et sont des contraintes d'épaisseur, au sens où, exception faite des zones de bord, la moyenne des contraintes sur toute l'épaisseur de la feuille de verre est nulle. Les contraintes superficielles en compression sont en effet équilibrées par la présence d'une zone centrale en 25 tension. Il existe donc une certaine profondeur à laquelle se produit la transition entre compression et tension, cette profondeur étant appelée "profondeur de compression P" dans la suite du texte. Dans le cadre de l'invention, la ou chaque feuille de verre du vitrage peut également être une feuille en verre organique, contenant une ou plusieurs 30 substances polymères organiques de haut poids moléculaire, par exemple en polycarbonate (PC) ou en polyméthacrylate de méthyle (PMMA). Par ailleurs, dans le cas d'un vitrage feuilleté, le ou chaque intercalaire polymérique du vitrage feuilleté peut être un feuillet thermoformable ou sensible à la pression, notamment à base de polybutyral de vinyle (PVB), de polyuréthane (PU), d'éthylène-acétate de vinyle (EVA), de polyéthylène téréphtalate (PET), de polychlorure de vinyle (PVC). Selon l'invention, on considère un impact dynamique entre un vitrage et un impacteur, où la vitesse relative entre le vitrage et l'impacteur est comprise entre 5 15 m/s et 1500 m/s. L'impacteur peut être de nature variée, en particulier l'impacteur peut être un élément dur tel qu'une bille en acier, un projectile de type pavé, un projectile balistique, ou bien l'impacteur peut être un élément mou tel qu'un oiseau. Il peut également s'agir d'un impacteur de type fluide, par exemple un gaz sous pression 10 dans le cas d'un vitrage soumis à un impact explosif, ou encore un volume d'eau projeté dans le cas d'un vitrage soumis à un impact avec un paquet d'eau, notamment pour des applications marines. A titre d'exemples non limitatifs d'impacts envisageables dans le cadre de l'invention, on peut citer les impacts suivants, correspondant à des essais 15 normalisés : - un impact à l'oiseau, utilisé pour tester des vitrages d'aéronefs (avions, hélicoptères), où l'impacteur est un poulet de 0,5 kg à 2 kg et la vitesse relative entre le vitrage et l'impacteur est comprise entre 50 m/s et 200 m/s; - un impact au projectile UIC (Union Internationale des Chemins de fer) 20 utilisé pour tester des vitrages de train selon la norme européenne des chemins de fer, où l'impacteur est le projectile UIC et la vitesse relative entre le vitrage et l'impacteur est comprise entre 20 m/s et 220 m/s; - un impact à la bouteille en verre, utilisé pour tester des vitrages de train, où l'impacteur est une bouteille en verre et la vitesse relative entre le vitrage et 25 l'impacteur est comprise entre 25 m/s et 180 mis; - un impact de type gravillonnage, utilisé pour tester des vitrages de train, où l'impacteur est un élément de 20 g en aluminium ayant une tête pointue et la vitesse relative entre le vitrage et l'impacteur est comprise entre 25 m/s et 150 m/s, 30 - un impact aux grêlons, utilisé pour tester des vitrages d'aéronefs (avions, hélicoptères), où l'impacteur est un ensemble de coton et d'eau gelée de diamètre prédéfini (10 mm ; 12,7 mm ; 25,4 mm ; 50,8 mm) et la vitesse relative entre le vitrage et l'impacteur est comprise entre 40 m/s et 260 m/s ; - un impact balistique, utilisé pour tester des vitrages de véhicules blindés, où l'impacteur est un projectile balistique pouvant être de différentes formes et de différents calibres et la vitesse relative entre le vitrage et l'impacteur est comprise entre 350 m/s et 1000 m/s. 5 Selon une caractéristique, le modèle numérique par éléments finis est obtenu en réalisant un maillage de modèles géométriques de l'impacteur, d'une part, et du vitrage avec les éléments environnants qui le maintiennent en position dans la structure, d'autre part. Ces modèles géométriques peuvent notamment être réalisés à l'aide d'un 10 logiciel de dessin assisté par ordinateur (DAO) ou de conception assistée par ordinateur (CAO), tel que les logiciels AUTOCAD, CATIA, PROENGINEER/CREO, SOLIDWORKS. De manière avantageuse, le maillage des modèles géométriques de l'impacteur et du vitrage avec ses éléments environnants, ainsi que le calcul de 15 l'enveloppe de contrainte maximale sur chaque feuille de verre du vitrage, sont réalisés à l'aide d'un logiciel d'analyse par éléments finis, tel que les logiciels ABAQUS, ANSYS, RADIOSS. Selon un aspect de l'invention, on fournit en entrée, pour le calcul par éléments finis, les propriétés des matériaux de l'impacteur, du vitrage et de ses 20 éléments environnants, sur au moins les gammes de fréquences et de températures caractéristiques de l'impact. Selon un autre aspect de l'invention, on fournit en entrée, pour le calcul par éléments finis, les caractéristiques de l'impact dynamique, notamment le lieu et l'angle d'impact de l'impacteur sur le vitrage, la vitesse relative entre le vitrage et 25 l'impacteur, la masse de l'impacteur, la température de chaque composant. Quel que soit le mode de réalisation, un procédé selon l'invention comprend : - la mise en place d'un modèle numérique par éléments finis de l'impacteur et du vitrage avec les éléments environnants qui le maintiennent en position dans 30 la structure, - l'injection dans le modèle numérique par éléments finis des propriétés des matériaux de l'impacteur, du vitrage et de ses éléments environnants, sur au moins les gammes de fréquences et de températures caractéristiques de l'impact, - l'injection dans le modèle numérique par éléments finis des caractéristiques de l'impact dynamique, notamment le lieu et l'angle d'impact de l'impacteur sur le vitrage, la vitesse relative entre le vitrage et l'impacteur, la masse de l'impacteur, la température de chaque composant, 5 - la sélection d'une méthode expérimentale de détermination de la contrainte à rupture d'une feuille de verre correspondant au type d'impact, et l'obtention de la contrainte à rupture de chaque feuille de verre du vitrage selon la méthode sélectionnée, afin de la comparer avec l'enveloppe de contrainte maximale sur cette feuille de verre calculée à l'aide du modèle numérique par 10 éléments finis. Une étape importante de l'invention est la sélection d'une méthode expérimentale de détermination de la contrainte à rupture d'une feuille de verre correspondant au type d'impact, c'est-à-dire dans laquelle l'impacteur sollicite la feuille de verre de manière similaire à ce qui se passe lors de l'impact réel. En 15 particulier, la méthode sélectionnée doit être représentative du type de sollicitations des défauts critiques présents dans la feuille de verre, qui peut dépendre notamment de la composition du verre, du type de traitement appliqué au verre (trempe thermique, trempe chimique, etc.), du type d'impacteur, de la vitesse d'impact. Ainsi, par exemple, la méthode sélectionnée ne sera pas la 20 même pour un vitrage d'avion soumis à un impact à l'oiseau, pour un vitrage de train soumis à un impact au pavé, ou encore pour un vitrage automobile soumis à un impact balistique. Des exemples de méthodes expérimentales de détermination de la contrainte à rupture d'une feuille de verre comprennent, notamment : un test 25 d'impact à la tour de chute ; un test de flexion par anneau sur tripode sans indentation ; un test de flexion par anneau sur tripode après indentation. Dans le test d'impact à la tour de chute, un impact d'un impacteur rigide, qui est une bille d'acier, est réalisé sur une éprouvette de la feuille de verre préalablement instrumentée de jauges de déformation. La bille d'acier est 30 positionnée à différentes hauteurs, jusqu'à obtenir la rupture de l'éprouvette. Parallèlement, un modèle numérique par éléments finis de l'essai est réalisé, en modélisant les jauges de déformation dans le modèle numérique. Sur l'essai réel, les déformations en dynamique sont relevées grâce aux jauges de déformation afin de valider le modèle numérique, ce qui permet de déterminer la contrainte au départ de casse, qui correspond à la contrainte à rupture de la feuille de verre. En pratique, deux jauges de déformation sont utilisées pour chaque essai, ce qui permet de voir le centrage de la bille vis-à-vis du centre de l'éprouvette. En s'appuyant sur les calculs du modèle numérique par éléments finis, on trace un 5 graphe donnant la contrainte au départ de casse en fonction de la hauteur de chute pour différents diamètres de billes. On peut ainsi déduire, en fonction des résultats des essais, la probabilité de rupture pour une contrainte donnée. De manière avantageuse, dans le test d'impact à la tour de chute, la feuille de verre est sollicitée de manière dynamique, ce qui est analogue à ce qui se passe lors 10 d'un impact à l'oiseau, un impact au projectile UIC, un impact balistique par exemple. Selon des modes de réalisation particuliers de l'invention : - le vitrage est un vitrage destiné à être installé dans un aéronef, l'impact est un impact à l'oiseau et la méthode sélectionnée pour déterminer la contrainte à 15 rupture de la feuille de verre est un test d'impact à la tour de chute, - le vitrage est un vitrage destiné à être installé dans un train, l'impact est un impact au projectile UIC et la méthode sélectionnée pour déterminer la contrainte à rupture de la feuille de verre est un test d'impact à la tour de chute, - le vitrage est un vitrage destiné à être installé dans un véhicule 20 automobile, l'impact est un impact balistique et la méthode sélectionnée pour déterminer la contrainte à rupture de la feuille de verre est un test d'impact à la tour de chute. Dans le test de flexion par anneau sur tripode sans indentation, une force croissante est appliquée sur une éprouvette de la feuille de verre placée entre 25 trois billes et un anneau, jusqu'à la rupture de l'éprouvette. Ce test permet de déterminer la contrainte à la rupture de défauts à la surface d'éprouvettes en s'affranchissant des bords du fait de la localisation de la contrainte maximale située sous l'anneau. De plus, la contrainte est constante et isotrope (égale dans toutes les directions) sous l'anneau de chargement. En pratique, une face 30 quelconque de chaque éprouvette est revêtue d'un film adhésif sur une face qui sera mise par la suite en compression. Le rôle de ce film est de permettre la localisation de l'origine de rupture. Le test de flexion par anneau sur tripode est réalisé par exemple à l'aide d'une machine Instron 5567, réglée avec une vitesse de descente de la traverse de 2 mm/min, instrumentée avec un capteur de force 10 de 10 kN, d'un anneau de 10 mm de diamètre avec un tore de 1 mm de rayon, fixée en bout de la machine lnstron, et d'un socle sur lequel sont collées trois billes de rayon 5 mm, disposées à 1200 sur un cercle de 20 mm de rayon et dont le centre est confondu avec le centre de l'anneau. L'éprouvette est placée entre 5 ces trois billes et l'anneau. On applique alors sur l'anneau une force croissante jusqu'à la rupture de l'éprouvette. Seules les éprouvettes dont l'origine de rupture est sous l'anneau sont comptabilisées. La contrainte à rupture en fonction de la force à rupture et de l'épaisseur de l'éprouvette est donnée par la formule suivante, le résultat étant pris comme la moyenne de cinq essais : 0,847 x Force(N) a(mPa) = épaisseur(m) 10 Dans le test de flexion par anneau sur tripode après indentation, le test de flexion est mis en oeuvre comme précédemment, à la différence que les éprouvettes ont été préalablement soumises à une indentation, réalisée sur la face opposée au film adhésif à l'aide de poids posés par-dessus une pointe Vickers. Pour l'indentation, chaque éprouvette est positionnée sous la pointe de telle sorte 15 que l'indentation soit réalisée au milieu de l'éprouvette, à 1 mm près. La descente de la pointe sur l'éprouvette est effectuée par exemple à l'aide d'une machine lnstron munie d'un capteur de force de 5 kN. En position initiale, la pointe est placée entre 2 et 5 mm au-dessus de l'éprouvette. Puis la pointe est approchée du verre à une vitesse de 10 mm/min. Après contact entre la pointe et le verre, la 20 force appliquée par la machine devient nulle et seuls les poids posés sur la pointe provoquent l'indentation du verre. L'indentation dure 20 secondes, puis la pointe est remontée par la machine. Le verre est alors entreposé pendant au moins 12 heures afin de stabiliser la propagation des fissures. En cas de rupture après indentation mais avant le test de flexion, la contrainte à rupture en flexion est 25 déclarée nulle. Pour le test de flexion, l'éprouvette est ensuite placée entre les trois billes et l'anneau de telle manière que la marque d'indentation soit alignée avec le centre de l'anneau, à 1 mm près. Pour une feuille de verre minéral, l'évolution de la probabilité de rupture en fonction de la contrainte obtenue selon le test d'impact à la tour de chute ou le test 30 de flexion par anneau sur tripode après ou sans indentation peut dépendre du volume de l'éprouvette testée. En effet, la résistance mécanique de la feuille de verre minéral est fixée par le défaut le plus important dans la zone sollicitée. En 11 changeant les conditions de sollicitation ou en prenant une éprouvette plus grande, on a une probabilité plus grande de rencontrer un défaut plus important. Pour prendre en compte ce phénomène, on utilise alors une méthode statistique reposant sur le modèle de Weibull. 5 Dans ce modèle, on exprime une probabilité cumulée de rupture PRi en fonction de la contrainte appliquée o-R. En effectuant N fois le même type d'essai, et si l'on ordonne les contraintes de rupture par ordre croissant : o-R1 1), en considérant que le test de flexion par anneau sur tripode après indentation est trop conservateur, étant donné qu'on ajoute un défaut. La valeur de ce facteur X est choisie par rapport à l'expérience et aux observations faites sur des tests destructifs. Dans un mode de réalisation, le vitrage est un vitrage feuilleté d'avion 20 constitué d'un empilement de trois feuilles de verre et de deux intercalaires polymériques intercalés entre les feuilles de verre. Une telle structure à trois plis de verre est une structure classique de vitrage frontal, latéral avant ou latéral arrière d'avion. En particulier, un exemple de structure classique pour un vitrage feuilleté 25 d'avion est l'empilement suivant : Verre minéral (3 mm) / PU (5,3 mm) / Verre minéral (8 mm) / PVB (2 mm) / Verre minéral (8 mm). Dans un autre mode de réalisation, le vitrage est un vitrage feuilleté d'hélicoptère constitué d'un empilement comprenant au moins une feuille de verre et un intercalaire polymérique. 30 En particulier, un exemple de structure classique pour un vitrage feuilleté d'hélicoptère est l'empilement suivant : Verre minéral (0,7 mm) / PU (2,5 mm) / PMMA (7 mm), ou encore un empilement d'une seule feuille de verre et d'un seul intercalaire polymérique composite : Verre minéral (3 mm) / PU (3,56 mm) + PET (0,18 mm). 13 Selon un aspect de l'invention, dans le cas d'un vitrage feuilleté, la loi de comportement du matériau constitutif de chaque intercalaire polymérique du vitrage feuilleté est un modèle viscoélastique déterminé à partir de mesures de DMA (analyse mécanique dynamique ou "Dynamic Mechanical Analysis"). La 5 DMA est utilisée pour caractériser la réponse d'un matériau en température et en fréquence, en appliquant de petites déformations cycliques. En pratique, à partir de résultats de DMA sur un échantillon de l'intercalaire polymérique, on étudie les propriétés en cisaillement du matériau de l'intercalaire en établissant : - la courbe d'évolution du module de conservation G' du matériau en 10 fonction de la fréquence pour différentes températures, notamment la fréquence est comprise entre 5 Hz et 285 Hz et la température est comprise entre -60°C et +60°C, - la courbe d'évolution du module de perte G" du matériau en fonction de la fréquence pour différentes températures, notamment la fréquence est comprise 15 entre 5 Hz et 285 Hz et la température est comprise entre -60°C et +60°C. A partir de ces données G'(f) et G"(f), on construit une courbe maîtresse pour les modules de conservation G' et de perte G" sur au moins les gammes de fréquences et de températures caractéristiques de l'impact, en utilisant par exemple la loi d'équivalence fréquence/température établie par la méthode WLF 20 (Williams-Landel-Ferry). On peut ensuite utiliser le modèle de Maxwell généralisé, qui permet de décrire les paramètres des matériaux selon une distribution en temps de relaxation, sous forme de série de Prony. La courbe maîtresse établie précédemment permet d'identifier les paramètres du modèle viscoélastique du 25 matériau constitutif de l'intercalaire polymérique, par identification (ou "fitting") d'une série de Prony sur la courbe maîtresse, sous la forme : t G(t) = Go 1 - gk (1 - e Tic) k=1 avec Go le module instantané (haute fréquence ou basse température), gk les modules relatifs, Tk les temps de relaxation. Un autre objet de l'invention est un vitrage obtenu par le procédé de 30 fabrication tel que décrit ci-dessus de manière à résister à un impact dynamique donné lorsqu'il est installé dans une structure donnée. 14 Dans un mode de réalisation, au moins une partie des étapes du procédé non destructif de validation qu'un vitrage installé dans une structure résiste à un impact dynamique tel que décrit précédemment, ou du procédé de fabrication d'un vitrage pour qu'il résiste à un impact dynamique lorsqu'il est installé dans une 5 structure tel que décrit précédemment, sont déterminées par des instructions de programmes d'ordinateurs. En conséquence, l'invention a également pour objet un programme d'ordinateur sur un support d'enregistrement, ce programme étant susceptible d'être mis en oeuvre dans un terminal, ou plus généralement dans un ordinateur, 10 ce programme comportant des instructions adaptées à la mise en oeuvre de tout ou partie des étapes d'un procédé tel que décrit précédemment. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code source, code objet, ou de code intermédiaire entre code source et code objet, tel que dans une forme partiellement compilée. 15 L'invention a aussi pour objet un support d'enregistrement lisible par un ordinateur, et comportant des instructions d'un programme d'ordinateur tel que mentionné ci-dessus. Le support d'enregistrement peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un 20 moyen de stockage, tel qu'une mémoire morte, une mémoire non volatile réinscriptible, par exemple une clé USB, une carte SD, une EEPROM, ou encore un moyen d'enregistrement magnétique, par exemple un disque dur. Le support d'enregistrement peut aussi être un circuit intégré dans lequel le programme est incorporé, le circuit étant adapté pour exécuter ou pour être utilisé 25 dans l'exécution du procédé. Le support d'enregistrement peut être un support transmissible tel qu'un signal électrique ou optique, qui peut être acheminé via un câble électrique ou optique, par radio ou par d'autres moyens. Le programme selon l'invention peut être en particulier téléchargé sur un réseau de type Internet. 30 L'invention a également pour objet un terminal comportant un module de traitement configuré pour : - calculer par analyse par éléments finis l'enveloppe de contrainte maximale sur chaque feuille de verre d'un vitrage installé dans une structure et soumis à un impact dynamique, où le vitrage comprend au moins une feuille de verre, à l'aide 15 d'un modèle numérique par éléments finis du vitrage installé dans la structure et soumis à l'impact, en utilisant une loi de comportement du matériau constitutif de chaque feuille de verre, ainsi que, dans le cas d'un vitrage feuilleté comprenant un empilement d'au moins une feuille de verre et un intercalaire polymérique, une loi 5 de comportement du matériau constitutif de chaque intercalaire polymérique et une loi de comportement de chaque interface entre une feuille de verre et un intercalaire polymérique, et - comparer l'enveloppe de contrainte maximale calculée pour chaque feuille de verre du vitrage à une valeur de contrainte à rupture de la feuille de verre 10 obtenue selon une méthode expérimentale de détermination de la contrainte à rupture du verre sélectionnée en fonction du type d'impact. Selon un mode de réalisation, le module de traitement du terminal est également configuré pour : - calculer par analyse par éléments finis l'enveloppe de contrainte maximale 15 sur chaque feuille de verre du vitrage en fonction des dimensions du vitrage, et - ajuster les dimensions du vitrage de telle sorte que l'enveloppe de contrainte maximale calculée pour chaque feuille de verre du vitrage est strictement inférieure à une valeur de contrainte à rupture de la feuille de verre, obtenue selon une méthode expérimentale de détermination de la contrainte à 20 rupture du verre sélectionnée en fonction du type d'impact, tout en ayant un dimensionnement optimisé du vitrage. Selon un aspect, le module de traitement du terminal selon l'invention comporte un programme d'ordinateur tel que mentionné ci-dessus, ce programme étant enregistré sur un support d'enregistrement conforme à l'invention et 25 constitué par une mémoire non volatile réinscriptible du terminal, les instructions du programme étant interprétables par un processeur du terminal. Le terminal, le programme d'ordinateur et le support d'enregistrement présentent, selon l'invention, les mêmes caractéristiques que le procédé selon l'invention. L'invention peut être mise en oeuvre avec tout type de terminal, par 30 exemple un ordinateur portable ou non. Enfin, l'invention a pour objet un système de validation par analyse par éléments finis qu'un vitrage installé dans une structure résiste à un impact dynamique, où le vitrage comprend au moins une feuille de verre, le système comprenant : 16 - une interface graphique, configurée pour afficher un modèle de l'impacteur et un modèle du vitrage avec ses éléments environnants, fournir des données en entrée pour l'analyse par éléments finis et afficher des résultats de l'analyse par éléments finis ; 5 - un module de modélisation des matériaux de l'impacteur et de chaque feuille de verre, et éventuellement de chaque intercalaire polymérique dans le cas d'un vitrage feuilleté, et des éléments environnants, pour définir les propriétés de ces matériaux sur au moins les gammes de fréquences et de températures caractéristiques de l'impact ; 10 - un module de modélisation de l'impact, pour définir notamment le lieu et l'angle d'impact de l'impacteur sur le vitrage, la vitesse relative entre le vitrage et l'impacteur, la masse de l'impacteur, la température de chaque composant ; - un module de traitement, pour préparer le modèle numérique par éléments finis du vitrage installé dans la structure et soumis à l'impact, réaliser 15 l'analyse par éléments finis, et calculer l'enveloppe de contrainte maximale sur chaque feuille de verre du vitrage. Dans un tel système, le module de traitement utilise les données définies dans le module de modélisation des matériaux et dans le module de modélisation de l'impact. 20 Les caractéristiques et avantages de l'invention apparaîtront dans la description qui va suivre d'un exemple de mise en oeuvre d'un procédé selon l'invention pour l'obtention d'un vitrage feuilleté d'avion qui résiste à un impact à l'oiseau, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés dans lesquels : 25 - la figure 1 est une vue schématique de face d'un cockpit d'avion comportant plusieurs vitrages feuilletés ou "glaces", respectivement en position frontale l(F), latérale avant l(FL) et latérale arrière 1(BL) ; - la figure 2 est une coupe schématique partielle d'un vitrage feuilleté latéral arrière 1(BL) d'avion et de ses éléments environnants lorsqu'il est installé dans la 30 structure d'un cockpit d'avion, le vitrage feuilleté étant constitué d'un empilement de trois feuilles de verre et deux intercalaires polymériques intercalés entre les feuilles de verre ; - la figure 3 est un diagramme schématique montrant les étapes successives d'un procédé selon l'invention, mis en oeuvre pour évaluer si le 17 vitrage feuilleté de la figure 2 installé dans la structure de l'avion résiste à un impact à l'oiseau ; - la figure 4 est un ensemble maillé destiné à être utilisé dans un modèle numérique par éléments finis réalisé à partir de modèles DAO d'un oiseau et du 5 vitrage feuilleté de la figure 2 avec ses éléments environnants, lors d'un impact de l'oiseau au centre du vitrage feuilleté ; - la figure 5 est un graphe montrant l'identification (ou "fitting") d'une série de Prony sur la courbe maîtresse du module de conservation G'(f), obtenue selon l'invention pour l'intercalaire polymérique en PU du vitrage feuilleté de la figure 2; 10 - la figure 6 est un graphe montrant l'enveloppe de contrainte maximale en fonction du temps pour une feuille de verre d'épaisseur 3 mm du vitrage feuilleté de la figure 2, telle que calculée à l'aide du modèle numérique par élément finis, dans une première configuration d'impact à l'oiseau sur le vitrage feuilleté ; - la figure 7 est un graphe montrant l'enveloppe de contrainte maximale en 15 fonction du temps pour une feuille de verre d'épaisseur 3 mm du vitrage feuilleté de la figure 2, telle que calculée à l'aide du modèle numérique par élément finis, dans une deuxième configuration d'impact à l'oiseau sur le vitrage feuilleté ; - la figure 8 est un graphe représentatif de la probabilité de rupture d'une feuille de verre de même composition verrière et même renforcement par trempe 20 chimique que les feuilles de verre du vitrage feuilleté de la figure 2 en fonction de la contrainte, obtenu selon le test d'impact à la tour de chute, où la probabilité de rupture dépend du volume de la feuille de verre et est donnée sur le graphe pour une feuille de verre de mêmes dimensions que la feuille de verre d'épaisseur 3 mm du vitrage feuilleté ; 25 - la figure 9 est un graphe représentatif de la probabilité de rupture d'une feuille de verre de même composition verrière et même renforcement par trempe chimique que les feuilles de verre du vitrage feuilleté de la figure 2 en fonction de la contrainte, obtenu selon le test de flexion par anneau sur tripode après indentation, où la probabilité de rupture est indépendante du volume de la feuille 30 de verre ; - la figure 10 est un diagramme schématique montrant les étapes successives d'un procédé de fabrication selon l'invention, mis en oeuvre pour obtenir le vitrage feuilleté de la figure 2 avec un dimensionnement optimisé lui 18 permettant de résister à un impact à l'oiseau lorsqu'il est installé dans la structure, tout en ayant une masse et/ou un coût minimisé ; et - la figure 11 est un diagramme schématique d'un système permettant de mettre en oeuvre un procédé selon l'invention. 5 Le procédé selon l'invention est mis en oeuvre pour vérifier qu'un vitrage feuilleté 1(BL), destiné à être intégré dans un cockpit d'avion en tant que glace latérale arrière, résiste à un impact à l'oiseau dans deux configurations différentes (exemples 1 et 2). Comme bien visible sur la figure 2, le vitrage feuilleté 1 est constitué d'un 10 empilement de trois feuilles de verre 11, 13, 15 et deux intercalaires polymériques 12, 14 intercalés entre les feuilles de verre. Une telle structure à trois plis de verre est une structure classique de vitrage feuilleté de cockpit d'avion. Chaque feuille de verre 11, 13, 15 est une feuille en verre alumino silicate qui a été renforcée par un procédé de trempe chimique. Pour chaque feuille de 15 verre 11, 13, 15, la profondeur de compression P résultant de la trempe chimique est de l'ordre de 200 à 250 pm. L'intercalaire polymérique 12 est un feuillet intercalaire à base de polyuréthane (PU). L'intercalaire polymérique 14 est un feuillet intercalaire à base de 20 polybutyral de vinyle (PVB). Les épaisseurs h_i des feuilles de verre i = 11, 13, 15 et h] des intercalaires polymériques j = 12, 14 sont les suivantes : Verre (h11 = 3 mm) / PU (h12 = 5,3 mm) / Verre (h13 = 8 mm) / PVB (h14 = 2 mm) / Verre (h15 = 8 mm). La figure 2 montre le vitrage feuilleté 1 avec ses éléments environnants 3, 25 5, 7 qui maintiennent le vitrage feuilleté 1 en position dans la structure de l'avion. Le vitrage feuilleté 1 est lié à la structure 7 (ou fuselage) de l'avion par l'intermédiaire d'un joint périphérique 3 en silicone, qui se comporte mécaniquement comme une rotule entre la structure 7 et le vitrage feuilleté 1. Une cale 5 en composite verre époxy est également fournie sur la périphérie intérieure 30 du vitrage feuilleté 1, afin de respecter l'entrefer défini par la structure 7. Le vitrage feuilleté 1 et ses éléments environnants 3, 5, 7 sont les mêmes pour les deux exemples 1 et 2, qui ne diffèrent l'un de l'autre que par les caractéristiques de l'impact à l'oiseau. La mise en oeuvre du procédé selon l'invention, pour vérifier que le vitrage feuilleté 1 intégré dans le cockpit d'avion 19 résiste à un impact à l'oiseau, est la même pour les deux exemples 1 et 2. Le procédé comprend les étapes montrées dans le diagramme de la figure 3 et décrites ci-après. Il est à noter que l'ordre des étapes de la figure 3 n'est pas impératif et peut faire l'objet de toute modification techniquement possible. 5 Aux étapes 110 et 120, on fournit respectivement un modèle géométrique CAD GLZ du vitrage feuilleté 1 avec ses éléments environnants 3, 5, 7, et un modèle géométrique CAD_BRD de l'oiseau 9. Pour les exemples 1 et 2, les modèles géométriques CAD_BRD et CAD_GLZ ont été réalisés à l'aide du logiciel CATIA. 10 A l'étape 130, on réalise un maillage des modèles géométriques CAD_BRD et CAD GLZ et on obtient un modèle numérique par élément finis FE IMP du vitrage feuilleté 1 installé dans la structure 7 de l'avion et soumis à l'impact avec l'oiseau 9. Pour les exemples 1 et 2, le maillage des modèles géométriques CAD_BRD et CAD GLZ a été réalisé à l'aide de l'outil de maillage HYPERMESH 15 et le codage a été réalisé à l'aide du code de calcul par éléments finis ABAQUS EXPLICIT. La figure 4 montre un exemple de modèle obtenu à l'étape 130, comprenant les représentations maillées de l'oiseau 9 et du vitrage feuilleté 1 avec ses éléments environnants 3, 7. A l'étape 140, on fournit en entrée du modèle numérique par élément finis 20 FE IMP les propriétés des matériaux des composants maillés : - une loi de comportement MAT_BRD du matériau constitutif de l'oiseau, à savoir dans cet exemple une loi de comportement hydrodynamique sous la forme d'une équation d'état, qui se trouve dans la littérature scientifique ; - une loi de comportement MAT_i du matériau constitutif de chaque feuille 25 de verre i = 11, 13, 15, qui sont les mêmes pour les trois feuilles de verre, à savoir dans cet exemple une masse volumique p = 2450 kg/m3, un module d'Young E = 72 GPa, un coefficient de Poisson y = 0,23; - une loi de comportement VISCMOD] du matériau constitutif de chaque intercalaire polymérique j = 12, 14, à savoir dans cet exemple un modèle 30 viscoélastique qui a été déterminé pour chaque intercalaire polymérique ; - une loi de comportement (INT_ij) de chaque interface entre une feuille de verre i = 11, 13, 15 et un intercalaire polymérique j = 12, 14, en particulier dans cet 20 exemple on considère à chaque fois une adhésion et un contact parfaits entre les surfaces, de sorte que le maillage est coïncident ; - une loi de comportement LAW_SEAL du joint 3 et de la cale 5, à savoir dans cet exemple, pour le joint 3, une loi de comportement hyper-élastique de 5 type NEO HOOKE ou VAN DER WALLS et, pour la cale 7, une loi de comportement élastique. Dans les exemples 1 et 2, pour chaque intercalaire polymérique j = 12, 14, le modèle viscoélastique VISCMOD] du matériau constitutif de l'intercalaire a été déterminé en procédant aux étapes suivantes : 10 - on a établi, à partir de résultats de DMA sur un échantillon de l'intercalaire j, la courbe d'évolution du module de conservation G'(f) et la courbe d'évolution du module de perte G"(f) du matériau de l'intercalaire j, pour une fréquence comprise entre 5 Hz et 285 Hz et différents isothermes compris entre -60°C et +60°C ; - on a construit, à partir des données G'(f) et G"(f), une courbe maîtresse 15 pour les modules de conservation G' et de perte G", sur des gammes de fréquences et de températures allant du plateau vitreux au plateau caoutchoutique du matériau, en utilisant la loi d'équivalence fréquence/température établie par la méthode WLF (Williams-Landel-Ferry) ; - on a identifié les paramètres du modèle viscoélastique du matériau 20 constitutif de l'intercalaire j, par identification d'une série de Prony sur la courbe maîtresse, sous la forme : / N t G(t) = Go 1 gk (1 - e Tic) - k=1 avec Go le module instantané, gk les modules relatifs, Tk les temps de relaxation. La figure 5 montre un exemple de "fitting" d'une série de Prony sur la 25 courbe maîtresse du module de conservation G'(f) du matériau de l'intercalaire polymérique 12 en PU du vitrage feuilleté 1. En pratique, les données matériaux ont été définies sous le format du logiciel ABAQUS, par exemple pour l'intercalaire polymérique 12 en PU (les valeurs sont données en unités SI) : 30 *Material, name=PU Visco *Density 21 1070e-09, *Elastic, moduli=LONG TERM 7e6, 0.49 *Viscoelastic, time=PRONY 5 0.085, 0., le-10 0.075, 0., le-09 0.06, 0., le-08 0.063, 0., le-07 0.06, 0., le-06 10 0.054, 0., le-05 0.052, 0., 0.0001 0.037, 0., 0.001 0.024, 0., 0.01 0.01, 0., 0.1 15 0.0066, 0., 1. 0.0019, 0., 10. 0.0010, 0., 100. 0.00097, 0., 1000. 20 A l'étape 150, on fournit en entrée du modèle numérique par élément finis FE_IMP les caractéristiques de l'impact, notamment : - lieu et angle d'impact de l'oiseau : au centre du vitrage feuilleté, avec l'oiseau qui se déplace parallèlement à la trajectoire de l'avion ; - vitesse relative entre le vitrage feuilleté et l'oiseau : 152,8 m/s pour 25 l'exemple 1 et 187,6 m/s pour l'exemple 2; - masse de l'oiseau : 1,807 kg pour l'exemple 1 et 1,812 kg pour l'exemple 2; - température de chaque composant : > pour l'exemple 1, T°C ambiante : 25°C, T°C face interne du vitrage 1 : 20,4°C, 30 T°C face externe du vitrage 1 : 23,9°C, > pour l'exemple 2, T°C ambiante : 23°C, T°C face interne du vitrage 1 : 23°C, T°C face externe du vitrage 1 : 24,3°C. A l'étape 160, on calcule par analyse par éléments finis à l'aide du modèle numérique FE_IMP l'enveloppe de contrainte maximale o-m_i sur chaque feuille de verre i = 11, 13, 15 du vitrage feuilleté 1. En pratique, pour les exemples 1 et 2, on a utilisé le solveur du logiciel ABAQUS pour calculer les champs de contraintes et de déformations induits par l'impact à l'oiseau sur chaque feuille de verre i = 11, 13, 15 et chaque intercalaire polymérique j = 12, 14 du vitrage feuilleté 1. Les figures 6 et 7 montrent, respectivement pour l'exemple 1 et pour l'exemple 2, l'enveloppe de contrainte maximale o-m_11 en fonction du temps calculée à l'aide du modèle numérique FE_IMP pour la feuille de verre 11 d'épaisseur 3 mm du vitrage feuilleté 1. 22 A l'étape 170, on compare l'enveloppe de contrainte maximale o-m_i sur chaque feuille de verre i = 11, 13, 15 du vitrage feuilleté 1 avec la contrainte à rupture o-r_i de la feuille de verre obtenue selon une méthode sélectionnée pour correspondre au type de sollicitations propres à un impact à l'oiseau, et on déduit 5 si le vitrage feuilleté 1 résiste à l'impact dans la configuration d'impact à l'oiseau considérée. Comme expliqué précédemment, les inventeurs ont mis en évidence que le test d'impact à la tour de chute sollicite la feuille de verre de manière similaire à ce qui se passe lors d'un impact à l'oiseau. Dès lors, pour les exemples 1 et 2, pour 10 chaque feuille de verre i du vitrage feuilleté 1, l'enveloppe de contrainte maximale o-m_i peut être comparée avec les résultats de probabilité de rupture d'une feuille de verre minéral de même composition verrière, même renforcement par trempe chimique et même volume que la feuille de verre i en fonction de la contrainte, obtenu selon le test d'impact à la tour de chute. 15 Ainsi, pour la feuille de verre 11 d'épaisseur 3 mm du vitrage feuilleté 1, on peut utiliser le graphe de la figure 8 montrant la probabilité de rupture d'une feuille de verre minéral de mêmes dimensions, même composition verrière et même renforcement que la feuille 11 en fonction de la contrainte, obtenu selon le test d'impact à la tour de chute. On choisit comme valeur de la contrainte à rupture 20 o-r_11, à laquelle on va comparer l'enveloppe de contrainte maximale o-m_11, la valeur de contrainte à 10% de probabilité de rupture sur le graphe de la figure 8, soit o-r_11 = 590 MPa. Pour l'exemple 1, par comparaison de la figure 6 avec la figure 8, on constate que l'enveloppe de contrainte maximale o-m_11 calculée à l'aide du 25 modèle numérique FE_IMP pour la feuille de verre 11 reste au cours du temps toujours strictement inférieure à une valeur de 450 MPa, ce qui est strictement inférieur à la contrainte à rupture o-r_11 d'une feuille de verre d'épaisseur 3 mm obtenue selon le test d'impact à la tour de chute. Cela permet de valider que la feuille de verre 11 du vitrage feuilleté 1 ne présente pas de risque de casse lors de 30 l'impact à l'oiseau selon l'exemple 1. En outre, une analyse analogue (non illustrée sur les figures) réalisée pour les feuilles de verre 13 et 15 d'épaisseur 8 mm du vitrage feuilleté 1 permet de valider que les feuilles de verre 13 et 15 du vitrage feuilleté 1 ne présentent pas non plus de risque de casse lors de l'impact à 23 l'oiseau selon l'exemple 1. Le vitrage feuilleté 1 est donc considéré comme résistant à l'impact à l'oiseau selon l'exemple 1 (résultat : OK). A contrario, pour l'exemple 2, par comparaison de la figure 7 avec la figure 8, on constate que l'enveloppe de contrainte maximale o-m_11 calculée à l'aide du 5 modèle numérique FE_IMP pour la feuille de verre 11 atteint au cours du temps une valeur de 600 MPa, ce qui est supérieur à la contrainte à rupture o-r_11 d'une feuille de verre d'épaisseur 3 mm obtenue selon le test d'impact à la tour de chute. Il en résulte que la feuille de verre 11 du vitrage feuilleté 1 présente un risque de casse lors de l'impact à l'oiseau selon l'exemple 2 et le vitrage feuilleté 1 est 10 considéré comme ne résistant pas à l'impact à l'oiseau selon l'exemple 2 (résultat : NOK). En variante, pour les exemples 1 et 2, pour chaque feuille de verre i du vitrage feuilleté 1, il est également possible de comparer l'enveloppe de contrainte maximale o-m_i avec les résultats de probabilité de rupture d'une feuille de verre 15 minéral de même composition verrière et même renforcement par trempe chimique que la feuille de verre i et de volume quelconque en fonction de la contrainte, obtenus selon le test de flexion par anneau sur tripode après indentation, en choisissant une profondeur d'indentation supérieure à la taille de défaut maximale du verre et inférieure à la profondeur de compression P résultant 20 de la trempe chimique. En particulier, pour chacune des feuilles de verre 11, 13, 15 du vitrage feuilleté 1, on peut utiliser le graphe de la figure 9 montrant la probabilité de rupture d'une feuille de verre minéral de même composition verrière et même renforcement en fonction de la contrainte, obtenu selon le test de flexion 25 par anneau sur tripode après indentation avec une profondeur d'indentation comprise entre 5 et 25 pm. De manière très avantageuse, le graphe de la figure 9 peut être utilisé directement pour chacune des trois feuilles de verre 11, 13, 15 du vitrage feuilleté 1, car dans ce cas la probabilité de rupture est indépendante du volume de la feuille de verre. On choisit comme valeur de la contrainte à rupture 30 o-r_i, à laquelle on va comparer l'enveloppe de contrainte maximale o-m_i, la valeur de contrainte à 10% de probabilité de rupture sur le graphe de la figure 9, soit o-r_i = 450 MPa. Ici encore, pour l'exemple 1, par comparaison de la figure 6 avec la figure 9, on constate que l'enveloppe de contrainte maximale o-m_11 calculée à l'aide du 24 modèle numérique FE_IMP pour la feuille de verre 11 reste au cours du temps strictement inférieure à 450 MPa, qui est la valeur de la contrainte à rupture o-r_11 obtenue selon le test de flexion par anneau sur tripode après indentation. Cela permet de valider que la feuille de verre 11 du vitrage feuilleté 1 ne présente pas 5 de risque de casse lors de l'impact à l'oiseau selon l'exemple 1. Une analyse analogue peut être réalisée pour les feuilles de verre 13 et 15. Le vitrage feuilleté 1 est donc considéré comme résistant à l'impact à l'oiseau selon l'exemple 1 (résultat : OK). A contrario, pour l'exemple 2, par comparaison de la figure 7 avec la figure 10 9, on constate que l'enveloppe de contrainte maximale o-m_11 calculée à l'aide du modèle numérique FE_IMP pour la feuille de verre 11 atteint au cours du temps une valeur de 600 MPa, ce qui est supérieur à la contrainte à rupture o-r_11 obtenue selon le test de flexion par anneau sur tripode après indentation. Il en résulte que la feuille de verre 11 du vitrage feuilleté 1 présente un risque de casse 15 lors de l'impact à l'oiseau selon l'exemple 2 et le vitrage feuilleté 1 est considéré comme ne résistant pas à l'impact à l'oiseau selon l'exemple 2 (résultat : NOK). Ces résultats sont bien conformes aux essais réels de tests destructifs à l'oiseau effectués sur le vitrage feuilleté 1. Il ressort des exemples précédents que la valeur de contrainte à rupture 20 déterminée selon le test de flexion par anneau sur tripode après indentation fournit un critère plus sévère que celle déterminée selon le test d'impact à la tour de chute. En effet, avec le test de flexion par anneau sur tripode après indentation, on ajoute un défaut au verre, ce qui conduit à surévaluer la probabilité de rupture du verre. Il est alors possible d'ajouter un facteur X (X>1) à la valeur de contrainte à 25 10% de probabilité de rupture, ce facteur X étant choisi de manière empirique, en comparant aux observations faites sur des tests destructifs. La figure 10 montre les étapes d'un procédé de fabrication du vitrage feuilleté 1 pour qu'il résiste à un impact à l'oiseau lorsqu'il est installé dans un cockpit d'avion. Ici encore, l'ordre des étapes de la figure 10 n'est pas impératif et 30 peut faire l'objet de toute modification techniquement possible. Les étapes 220, 230, 240, 250 du procédé de la figure 10 sont identiques respectivement aux étapes 120, 130, 140, 150 du procédé de la figure 3. Le procédé de la figure 10 diffère de celui de la figure 3 en ce que : 25 - à l'étape 210, on fournit un modèle géométrique CAD_GLZ_p du vitrage feuilleté 1 avec ses éléments environnants 3, 5, 7 qui est un modèle paramétré, en définissant en tant que paramètres du modèle l'épaisseur h_i des feuilles de verre i = 11, 13, 15 et l'épaisseur h] des intercalaires polymériques j = 12, 14 du vitrage 5 feuilleté 1 ; - à l'étape 260, on calcule par analyse par éléments finis à l'aide du modèle numérique FE_IMP l'enveloppe de contrainte maximale o-m_i_p sur chaque feuille de verre i = 11, 13, 15 du vitrage feuilleté 1 en fonction de l'épaisseur h_i des feuilles de verre et de l'épaisseur h] des intercalaires polymériques ; 10 - à l'étape 270, on ajuste l'épaisseur h_i des feuilles de verre et l'épaisseur h] des intercalaires polymériques de telle sorte que l'enveloppe de contrainte maximale o-m_i_p sur chaque feuille de verre du vitrage feuilleté est strictement inférieure à la contrainte à rupture o-r_i de la feuille de verre obtenue selon une méthode sélectionnée pour correspondre au type de sollicitations d'un impact à 15 l'oiseau ; - à l'étape 280, une fois que les valeurs d'épaisseurs h_i et h] ajustées ont été calculées, on prépare les feuilles de verre i = 11, 13, 15 et les intercalaires polymériques j = 12, 14 à ces épaisseurs h_i, h] ajustées, et on les assemble de manière à former le vitrage feuilleté 1. 20 Un tel procédé de fabrication de vitrage feuilleté garantit l'obtention d'un vitrage feuilleté 1 optimisé à la fois en termes de masse, de coût et de résistance à l'impact dans la configuration d'impact à l'oiseau définie à l'étape 250. La figure 11 montre un système 30 selon l'invention, susceptible d'être utilisé pour mettre en oeuvre le procédé décrit précédemment en lien avec la figure 25 3 pour vérifier que le vitrage feuilleté 1 intégré dans un cockpit d'avion en tant que glace latérale arrière résiste à un impact à l'oiseau, et/ou le procédé décrit ci-dessus en lien avec la figure 10 pour fabriquer le vitrage feuilleté 1 pour qu'il résiste à un impact un impact à l'oiseau lorsqu'il est installé dans un cockpit d'avion. 30 Le système 30 comprend une interface graphique utilisateur 31, un module de modélisation des matériaux 32, un module de modélisation de l'impact 33 et un module de traitement 34. En référence aux figures 3 et 10, l'interface graphique 31 est configurée pour afficher les modèles de l'impacteur et du vitrage feuilleté avec ses éléments 26 environnants obtenus aux étapes 110/210, 120/220, 130/230; fournir les données d'entrée aux étapes 140/240, 150/250; afficher les résultats de l'analyse par éléments finis à l'étape 170/270. Le module de modélisation des matériaux 32 est configuré pour stocker et gérer les données matériaux fournies à l'étape 140/240. 5 Le module de modélisation de l'impact 33 est configuré pour stocker et gérer les données caractéristiques de l'impact fournies à l'étape 150/250. Le module de traitement 34 est configuré pour préparer le modèle numérique par éléments finis FE IMP à l'étape 130/230; réaliser l'analyse par éléments finis et calculer l'enveloppe de contrainte maximale o-m_i, o-m_i_p à l'étape 160/260 et, dans le 10 cas du procédé de la figure 10, les valeurs d'épaisseurs h_i et h] ajustées à l'étape 270. Pour l'analyse par éléments finis, le module de traitement 34 utilise les données définies dans le module de modélisation des matériaux 32 et le module de modélisation de l'impact 33. L'invention n'est pas limitée aux exemples décrits et représentés. En 15 particulier, l'invention a été illustrée avec des exemples d'impacts à l'oiseau sur un vitrage feuilleté d'avion, mais il est bien entendu qu'elle est applicable pour tout type d'impact dynamique et tout type de vitrage comprenant au moins une feuille de verre, qu'il s'agisse d'un vitrage feuilleté ou non. Une condition importante pour la bonne mise en oeuvre de l'invention est la sélection, pour déterminer la 20 contrainte à rupture d'une feuille de verre, d'une méthode représentative des sollicitations associées au type d'impact considéré et d'une REVENDICATIONS1. Procédé non destructif de validation qu'un vitrage (1) installé dans une structure (7) résiste à un impact dynamique, le vitrage (1) comprenant au moins une feuille de verre (11, 13, 15), caractérisé en ce qu'il comprend des étapes dans lesquelles : - à l'aide d'un modèle numérique par éléments finis (FE_IMP) du vitrage (1) installé dans la structure (7) et soumis à l'impact, en utilisant une loi de comportement (MAT_i) du matériau constitutif de chaque feuille de verre (11, 13, 15), on calcule l'enveloppe de contrainte maximale (o-m_i) sur au moins une feuille de verre critique (11, 13, 15) du vitrage ; - pour au moins la feuille de verre critique (11, 13, 15) du vitrage, on compare l'enveloppe de contrainte maximale (o-m_i) à la contrainte à rupture (o-r_i) de la feuille de verre obtenue selon une méthode expérimentale de détermination de la contrainte à rupture d'une feuille de verre sélectionnée en fonction du type d'impact. 2. Procédé selon la revendication 1, caractérisé en ce que le vitrage (1) est un vitrage feuilleté comprenant un empilement d'au moins une feuille de verre (11, 13, 15) et un intercalaire polymérique (12, 14), le procédé comprenant des étapes dans lesquelles : - à l'aide d'un modèle numérique par éléments finis (FE_IMP) du vitrage feuilleté (1) installé dans la structure (7) et soumis à l'impact, en utilisant une loi de comportement (MAT_i) du matériau constitutif de chaque feuille de verre (11, 13, 15), une loi de comportement (VISCMOD]) du matériau constitutif de chaque intercalaire polymérique (12, 14), et une loi de comportement (INT_ij) de chaque interface entre une feuille de verre et un intercalaire polymérique, on calcule l'enveloppe de contrainte maximale (o-m_i) sur au moins une feuille de verre critique (11, 13, 15) du vitrage feuilleté ; - pour au moins la feuille de verre critique (11, 13, 15) du vitrage feuilleté, on compare l'enveloppe de contrainte maximale (o-m_i) à la contrainte à rupture (o-r_i) de la feuille de verre obtenue selon une méthode expérimentale de détermination de la contrainte à rupture d'une feuille de verre sélectionnée en fonction du type d'impact. 28 3. Procédé de fabrication d'un vitrage (1) pour qu'il résiste à un impact dynamique lorsqu'il est installé dans une structure (7), le vitrage (1) comprenant au moins une feuille de verre (11, 13, 15), caractérisé en ce qu'il comprend des étapes dans lesquelles : 5 - à l'aide d'un modèle numérique par éléments finis (FE_IMP) du vitrage (1) installé dans la structure (7) et soumis à l'impact, en utilisant une loi de comportement (MAT_i) du matériau constitutif de chaque feuille de verre (11, 13, 15), on calcule l'enveloppe de contrainte maximale (o-m_i_p) sur au moins une feuille de verre critique (11, 13, 15) du vitrage en fonction des dimensions (h_i, 10 h_j) du vitrage ; - on ajuste les caractéristiques du vitrage (1) parmi ses dimensions (h_i, h_j), le matériau constitutif de chaque feuille de verre (11, 13, 15)' de telle sorte que l'enveloppe de contrainte maximale (o-m_i_p) sur au moins la feuille de verre critique (11, 13, 15) du vitrage est strictement inférieure à la contrainte à rupture 15 (o-r_i) de la feuille de verre obtenue selon une méthode expérimentale de détermination de la contrainte à rupture d'une feuille de verre sélectionnée en fonction du type d'impact, tout en ayant un dimensionnement optimisé du vitrage ; - on prépare et on assemble chaque feuille de verre (11, 13, 15) du vitrage avec les caractéristiques ajustées (h_i, h]). 20 4. Procédé selon la revendication 3, caractérisé en ce que le vitrage (1) est un vitrage feuilleté comprenant un empilement d'au moins une feuille de verre (11, 13, 15) et un intercalaire polymérique (12, 14), le procédé comprenant des étapes dans lesquelles : - à l'aide d'un modèle numérique par éléments finis (FE_IMP) du vitrage 25 feuilleté (1) installé dans la structure (7) et soumis à l'impact, en utilisant une loi de comportement (MAT_i) du matériau constitutif de chaque feuille de verre (11, 13, 15), une loi de comportement (VISCMOD]) du matériau constitutif de chaque intercalaire polymérique (12, 14), et une loi de comportement (INT_ij) de chaque interface entre une feuille de verre et un intercalaire polymérique, on calcule 30 l'enveloppe de contrainte maximale (o-m_i_p) sur au moins une feuille de verre critique (11, 13, 15) du vitrage feuilleté en fonction des dimensions (h_i, h_j) du vitrage feuilleté ; - on ajuste les caractéristiques du vitrage feuilleté (1) parmi ses dimensions (h_i, h_j), le matériau constitutif de chaque feuille de verre (11, 13, 15), le matériau 29 constitutif de chaque intercalaire polymérique (12, 14), la nature de chaque interface entre une feuille de verre et un intercalaire polymérique, de telle sorte que l'enveloppe de contrainte maximale (o-m_i_p) sur au moins la feuille de verre critique (11, 13, 15) du vitrage feuilleté est strictement inférieure à la contrainte à 5 rupture (o-r_i) de la feuille de verre obtenue selon une méthode expérimentale de détermination de la contrainte à rupture d'une feuille de verre sélectionnée en fonction du type d'impact, tout en ayant un dimensionnement optimisé du vitrage feuilleté ; - on prépare et on assemble chaque feuille de verre (11, 13, 15) et chaque 10 intercalaire polymérique (12, 14) du vitrage feuilleté avec les caractéristiques ajustées (h_i, h]). 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on obtient le modèle numérique par éléments finis (FE_IMP) en réalisant un maillage de modèles géométriques de l'impacteur (CAD_BRD) et 15 du vitrage avec ses éléments environnants (CAD_GLZ). 6. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le maillage de modèles géométriques de l'impacteur (CAD BRD) et du vitrage (CAD_GLZ) et le calcul de l'enveloppe de contrainte maximale (o-m_i) sur chaque feuille de verre (11, 13, 15) du vitrage sont réalisés à 20 l'aide d'un logiciel d'analyse par éléments finis. 7. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on fournit en entrée, pour le calcul par éléments finis, les propriétés des matériaux de l'impacteur, du vitrage et des éléments environnants (MAT_BRD, MAT_i, VISCMOD], LAW_SEAL) sur au moins les gammes de 25 fréquences et de températures caractéristiques de l'impact. 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'on fournit en entrée, pour le calcul par éléments finis, les caractéristiques de l'impact, notamment le lieu et l'angle d'impact de l'impacteur sur le vitrage, la vitesse relative entre le vitrage et l'impacteur, la masse de 30 l'impacteur, la température de chaque composant. 9. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la méthode sélectionnée pour déterminer la contrainte à rupture de la feuille de verre est un test d'impact à la tour de chute, un test de 30 flexion par anneau sur tripode sans indentation ou un test de flexion par anneau sur tripode après indentation. 10. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le vitrage (1) comprend au moins une feuille de verre 5 minéral renforcée par trempe chimique (11, 13, 15) et la méthode sélectionnée pour déterminer la contrainte à rupture de la feuille de verre est un test de flexion par anneau sur tripode après indentation. 11. Procédé selon la revendication 10, caractérisé en ce que la profondeur d'indentation est choisie supérieure à la plus grande taille de défaut du verre et 10 inférieure à la profondeur de compression résultant de la trempe chimique. 12. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce que le vitrage (1) est un vitrage feuilleté d'avion constitué d'un empilement de trois feuilles de verre (11, 13, 15) et de deux intercalaires polymériques (12, 14) intercalés entre les feuilles de verre. 15 13. Procédé selon l'une quelconque des revendications 1 à 11, caractérisé en ce que le vitrage est un vitrage feuilleté d'hélicoptère constitué d'un empilement comprenant au moins une feuille de verre et un intercalaire polymérique. 14. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le vitrage (1) est un vitrage feuilleté comprenant un 20 empilement d'au moins une feuille de verre (11, 13, 15) et un intercalaire polymérique (12, 14) et la loi de comportement du matériau constitutif de chaque intercalaire polymérique (12, 14) du vitrage (1) est un modèle viscoélastique (VISCMOD j) déterminé en procédant aux étapes suivantes : - on établit, à partir de résultats de DMA sur un échantillon de l'intercalaire 25 polymérique (12, 14), la courbe d'évolution du module de conservation G'(f) du matériau en fonction de la fréquence pour différentes températures et la courbe d'évolution du module de perte G"(f) du matériau en fonction de la fréquence pour différentes températures ; - à partir des données G'(f) et G"(f), on construit une courbe maîtresse pour 30 les modules de conservation G' et de perte G" sur au moins les gammes de fréquences et de températures caractéristiques de l'impact, en utilisant par exemple la loi d'équivalence fréquence/température établie par la méthode WLF (Williams-Landel-Ferry) ; 31 - on identifie les paramètres du modèle viscoélastique du matériau constitutif de l'intercalaire polymérique (12, 14), par identification d'une série de Prony sur la courbe maîtresse, sous la forme : t G(t) = Go 1 - gk (1 - e Tic) k=1 avec Go le module instantané, gk les modules relatifs, Tk les temps de 5 relaxation. 15. Vitrage (1) destiné à résister à un impact dynamique donné lorsqu'il est installé dans une structure donnée, caractérisé en ce qu'il est obtenu par le procédé de l'une quelconque des revendications précédentes. 16. Programme d'ordinateur comportant des instructions pour l'exécution de 10 tout ou partie des étapes d'un procédé selon l'une quelconque des revendications 1 à 14 lorsque ledit programme est exécuté par un ordinateur. 17. Support d'enregistrement lisible par un ordinateur sur lequel est enregistré un programme d'ordinateur comprenant des instructions pour l'exécution de tout ou partie des étapes d'un procédé selon l'une quelconque des 15 revendications 1 à 14. 18. Terminal comportant un module de traitement (33) configuré pour : - calculer par analyse par éléments finis l'enveloppe de contrainte maximale (o-m_i, o-m_i_p) sur chaque feuille de verre d'un vitrage installé dans une structure et soumis à un impact dynamique, où le vitrage comprend au moins une feuille de 20 verre, à l'aide d'un modèle numérique par éléments finis (FE_IMP) du vitrage installé dans la structure et soumis à l'impact, en utilisant une loi de comportement (MAT_i) du matériau constitutif de chaque feuille de verre (11, 13, 15), et - comparer l'enveloppe de contrainte maximale (o-m_i, o-m_i_p) calculée pour chaque feuille de verre du vitrage à une valeur de contrainte à rupture (o-r_i) 25 de la feuille de verre obtenue selon une méthode expérimentale de détermination de la contrainte à rupture du verre sélectionnée en fonction du type d'impact. 19. Terminal selon la revendication 18, dans lequel le module de traitement (33) est également configuré pour : - calculer par analyse par éléments finis l'enveloppe de contrainte maximale 30 (o-m_i_p) sur chaque feuille de verre du vitrage en fonction des dimensions (h_i, h_j) du vitrage, et 32 - ajuster les dimensions (h_i, h_j) du vitrage de telle sorte que l'enveloppe de contrainte maximale (o-m_i_p) calculée pour chaque feuille de verre du vitrage est strictement inférieure à la contrainte à rupture (o-r_i) de la feuille de verre, obtenue selon une méthode expérimentale de détermination de la contrainte à 5 rupture du verre sélectionnée en fonction du type d'impact, tout en ayant un dimensionnement optimisé du vitrage. 20. Système (30) de validation par analyse par éléments finis qu'un vitrage (1) installé dans une structure (7) résiste à un impact dynamique, où le vitrage (1) comprend au moins une feuille de verre (11, 13, 15), le système comprenant : 10 - une interface graphique (31), configurée pour afficher des modèles de l'impacteur (9) et du vitrage (1) avec ses éléments environnants (3, 5, 7), fournir des données en entrée pour l'analyse par éléments finis et afficher des résultats de l'analyse par éléments finis ; - un module de modélisation des matériaux (32) de l'impacteur (9), de 15 chaque feuille de verre (11, 13, 15) du vitrage (1), et des éléments environnants (3, 5, 7), pour définir les propriétés de ces matériaux sur au moins les gammes de fréquences et de températures caractéristiques de l'impact ; - un module de modélisation de l'impact (33), pour définir notamment le lieu et l'angle d'impact de l'impacteur (9) sur le vitrage (1), la vitesse relative entre le 20 vitrage et l'impacteur, la température de chaque composant ; - un module de traitement (34), pour préparer le modèle numérique par éléments finis (FE_IMP) du vitrage (1) installé dans la structure (7) et soumis à l'impact, réaliser l'analyse par éléments finis, calculer l'enveloppe de contrainte maximale (o-m_i, o-m_i_p) sur chaque feuille de verre (11, 13, 15) du vitrage. 25 21. Système selon la revendication 20, dans lequel le module de traitement (34) utilise les données définies dans le module de modélisation des matériaux (32) et le module de modélisation de l'impact (33).