FR 2950178 A1 20110318 FR 0956286 A 20090914 La présente invention relève du domaine des structures. Il concerne en particulier les structures de type panneaux raidis, et encore plus particulièrement de tels panneaux, renforcés par des raidisseurs. L'invention s'intéresse au calcul de résistance de telles structures soumises à des charges combinées. Art antérieur et problème posé Les structures fines raidies représentent la plupart des structures primaires des aéronefs commerciaux. Les panneaux sont généralement renforcés avec des raidisseurs perpendiculaires les uns aux autres et définissant des zones rectangulaires sur la peau du panneau, limitées par les raidisseurs, et appelées poches. Une structure d'avion est ainsi conçue avec un squelette de raidisseurs et une peau : -• des raidisseurs longitudinaux (généralement appelés longerons): ils permettent de soutenir la structure dans la direction principale des charges -• des raidisseurs transversaux (généralement nommés "frame" ou "rib") : leur fonctionnalité principale est de fournir un support pour les longerons -• un panneau (généralement nommé peau): en règle générale, il 20 reprend le chargement dans le plan (membrane) Les longerons et les lisses sont orientés à 90° les uns des autres et définissent sur la peau des poches rectangulaires. Cependant, dans les années 1950 et 60, pour les structures spatiales, 25 la NASA a développé un nouveau concept de structure raidie appelée "Isogrid" (voir figure 1). Une telle structure raidie est ainsi composée d'une peau renforcée avec un réseau de raidisseurs orientés à 0° (0 =60°, dans les structures envisagées par la NASA) entre eux. Les raidisseurs sont en forme de lame et 30 sont intégrés au panneau. En raison de sa géométrie, cette configuration possède un comportement orthotrope (isotrope lorsque O = 60 °) et les poches formées sur la peau sont triangulaires. On utilise dans la suite de la description les termes de structure raidie par des poches triangulaires ou panneau raidi par des poches triangulaires pour définir des structures ou panneaux renforcés par des raidisseurs croisés formant des poches triangulaires. Des données limitées sont disponibles dans la littérature pour calculer la résistance et la stabilité d'une telle structure raidie par des poches triangulaires. Etat de l'art des méthodes de calcul de panneau raidis par des poches triangulaires Un procédé de calcul analytique de panneau raidi par des poches triangulaires équilatérales est décrit dans le document Rapport de Contrat NASA "Isogrid" design handbook" (NASA-CR-124075, 02/1973) Cette méthode est bien documentée, mais présente de sérieuses limitations : utilisation de triangles équilatéraux uniquement : angle = 60°, calcul de contraintes appliquées mais pas de calcul de contraintes admissibles, coefficient de Poisson du matériau valant uniquement 1/3. La méthode antérieure présente de nombreuses limitations et ne prend pas en compte tous les problèmes qui se présentent sur une structure d'aéronef, notamment en ce qui concerne les conditions aux limites et la plasticité. Elle n'est donc pas utilisable de façon fiable pour le calcul analytique de structure de panneaux raidis par des poches triangulaires Objectifs de l'invention Afin d'effectuer une analyse structurelle des panneaux raidis par des poches triangulaires, une méthode d'analyse structurale a été développée, basée sur une théorie de plaque composite et tenant compte de ses modes d'effondrement spécifiques. Cette méthode, s'applique aux panneaux plats réalisés dans un matériau à propriétés isotropes. Le procédé décrit ici envisage une modification de l'angle de base entre les raidisseurs (qui est de 60 °dans les structures "Isogrid"). Cela signifie que l'isotropie du panneau est n'est plus garantie. Exposé de l'invention L'invention vise à cet effet un procédé de dimensionnement par méthode analytique d'un panneau sensiblement plan en matériau homogène et isotrope, le panneau étant composé d'une peau renforcée par un ensemble (dit "grille") de trois faisceaux parallèles de raidisseurs intégrés au panneau, les poches déterminées sur la peau par ces familles de raidisseurs étant triangulaires, les raidisseurs étant en forme de lame, le panneau raidi devant satisfaire un cahier des charges de résistance mécanique à des charges externes prédéterminées, les angles entre faisceaux de raidisseurs étant tels que les poches triangulaires sont de forme isocèle quelconque. Selon une mise en oeuvre avantageuse, le procédé comporte des étapes : Etape 2 - de calcul des contraintes appliquées dans la peau et les raidisseurs, ainsi que de flux dans la peau et charges dans les raidisseurs, à partir de la géométrie du panneau raidi, et des charges externes, supposées situées dans le plan du panneau et appliquées au centre de gravité de la section (du panneau), le panneau raidi étant représenté par un assemblage de deux plaques orthotropes, la grille de raidisseurs étant représentée par un panneau équivalent. Etape 3 - de calcul des charges internes dans le panneau raidi, Etape 4 - d'analyse de résistance comportant un calcul des facteurs de réserve du matériau à charge limite et à charge ultime, Etape 5 û de calcul de contraintes locales admissibles, Préférentiellement, le procédé prend en compte la redistribution des contraintes appliquées entre le panneau et la grille de raidisseurs due : au post-flambage des raidisseurs, par la définition d'une section droite efficace pour chaque type de raidisseur (0°, +0 or -0): Aoô t , A±est et A_est, au post-flambage de la poche à travers le calcul d'une épaisseur efficace du panneau : ts_eff, à la plasticité des charges externes appliquées, grâce à un procédé itératif sur les différentes propriétés du matériau, notamment modules d'Young et coefficients de Poisson : Eoô t, E+est, E-est pour les raidisseurs et EXS, Eys et vepst pour la peau, utilisant la loi de Ramberg û Osgood. Selon un mode préféré de mise en oeuvre du procédé selon l'invention, celui-ci comporte une étape de correction des charges appliquées pour tenir compte de la plasticité, utilisant un procédé de calcul itératif de contraintes plastiques, réalisé jusqu'à ce que les cinq paramètres du matériau (Eo.st E+est E est, Epeau, vep) entrés au début du processus soient sensiblement égaux aux mêmes paramètres obtenus après calcul de contrainte plastique. Selon une mise en oeuvre avantageuse, le procédé comporte une étape 4, d'analyse de résistance comportant un calcul des facteurs de réserve du matériau à charge limite et à charge ultime, réalisée par comparaison des charges appliquées calculées dans les composants du panneau raidi, et des contraintes maximums admissibles du matériau, les charges appliquées étant corrigées pour tenir compte de la plasticité du panneau raidi. Selon une mise en oeuvre avantageuse, le procédé comporte une étape 5 de calcul de contraintes locales admissibles, laquelle comporte une sous-étape 5A de calcul de flux de flambage admissible, et de facteur de réserve pour des poches triangulaires isocèles, les contraintes appliquées à prendre en compte pour le calcul de facteur de réserve étant les contraintes agissant dans la peau uniquement, les flux externes utilisés étant des flux de peau ne correspondant pas au chargement complet du panneau raidi. Dans ce cas, l'étape 5A de calcul de flux de flambage admissible, et de facteur de réserve pour des poches triangulaires isocèles comporte favorablement deux sous-étapes : d'abord de calcul de valeurs admissibles pour des plaques soumises à des cas de chargement pur (compression selon deux directions dans le plan, cisaillement) en utilisant une méthode par éléments finis, puis de calcul des courbes d'interaction entre ces cas de chargement pur. Encore plus précisément, le calcul de valeurs admissibles comporte des sous-étapes suivantes : • Création d'un modèle paramétrique FEM d'une plaque triangulaire • Exécution de nombreuses combinaisons différentes pour obtenir des résultats de flambage, • Obtention de paramètres compatible avec une formulation analytique polynomiale. Dans un mode particulier de mise en oeuvre, dans le cas d'un chargement pur, les courbes d'interaction sont définies par des sous-étapes suivantes : - de création de modèles éléments finis d'une pluralité de plaques triangulaires présentant des angles isocèles différents, l'angle isocèle (8) étant défini comme l'angle de base du triangle isocèle, - pour chaque angle isocèle : 1/ de calcul par Modèle Eléments Finis pour déterminer le flux admissible de plissement (sans correction plastique) pour diverses épaisseurs de plaque. 2/ de tracé d'une courbe de flux de flambage admissible en fonction du rapport D (D raideur de plaque, h hauteur du triangle), cette courbe étant déterminée pour des petites valeurs du rapport D , par une équation au second degré fonction de ce rapport, dont les coefficients K1 et K2 dépendent de l'angle et du cas de charge considéré, 3/ de tracé de l'évolution des coefficients de l'équation polynomiale K1 et K2 en fonction de l'angle de base du triangle isocèle, ces coefficients étant tracés en fonction de l'angle des plaques triangulaires considérées, et d'interpolation pour déterminer une équation polynomiale qui permet de calculer ces constantes quel que soit l'angle isocèle. Toujours dans le cas de calcul de flux de flambage admissible, et de facteur de réserve pour des poches triangulaires isocèles, selon une mise en oeuvre avantageuse, dans le cas d'un chargement combiné, on utilise l'hypothèse suivante : si certains composants de la charge combinée sont en tension, alors ces composants ne sont pas pris en compte pour le calcul, et les courbes d'interaction sont définies par des sous-étapes suivantes : - de création de modèles éléments finis d'une pluralité de plaques triangulaires présentant des angles isocèles différents, l'angle isocèle (8) étant défini comme l'angle de base du triangle isocèle, - pour chaque angle, 1/ de calcul par Modèle Eléments Finis (FEM) pour déterminer la valeur propre de flambage correspondant à différentes répartitions de charges externes. 2/ de tracé des courbes d'interaction, pour chaque angle et chaque combinaison de charges, puis d'approximation de ces courbes avec une unique équation couvrant toutes les combinaisons : aPP RexA +ReyB +Rsc =1 (où Ri = N, i = cX, cY ou s), A, B, C étant des coefficients empiriques. Avantageusement, le procédé comporte en outre une sous-étape de calcul de facteurs de réserve, par résolution de l'équation suivante : iR 'A iR 'B iR \c cY + cX + s =1 R~ R, R, N aPP N aPP 1 cX = s = N cric cric cric RF cYcomb cXcomb scomb Selon une mise en oeuvre avantageuse, le procédé utilise, pour le calcul des contraintes admissibles corrigées de plasticité, un facteur de correction de plasticité r, défini par : 25 - pour tous les cas de chargement (pur et combiné) à l'exception du cisaillement, - pour les cas de chargement en cisaillement pur,: 20 aPP avec R = NcY (1+ve) Esec (1+v) E, la correction de plasticité étant calculée en utilisant la contrainte élastique équivalente de Von Mises. Selon une mise en oeuvre avantageuse, dans le cas de plaques triangulaires isocèles en appui simple ou encastrées, en cas de chargement 3 combiné, on utilise une courbe d'interaction : RX +RY +Rsz =1, pour tous les cas de chargement. Selon une mise en oeuvre avantageuse, le procédé comporte une étape 5 de calcul de contraintes locales admissibles, laquelle comporte une sous-étape 5B de calcul de contrainte de flambage admissible et de facteur de réserve de l'âme du raidisseur, considéré comme un panneau rectangulaire, les contraintes appliquées pour les calculs de facteur de réserve étant les contraintes dans les âmes des raidisseurs seulement. Selon une mise en oeuvre avantageuse du procédé, celui-ci comporte une étape 6, de calcul d'instabilité générale, fournissant des données de flux admissible de flambage, et des facteurs de réserve, pour un panneau plat raidi, dans des conditions de chargement pur ou combiné, les flux appliqués, à prendre en compte pour le calcul de facteur de réserve, étant les flux externes du panneau raidi. Dans ce cas, plus précisément, le procédé comporte avantageusement des sous-étapes : - d'utilisation d'une loi de comportement générale (équation 6-8), définissant les relations entre flux et moments, d'une part, et déformations, d'autre part, un état de contraintes plan étant considéré, - d'utilisation d'équations d'équilibre général (équations 6-9 et 6-10) d'un élément du panneau raidi, reliant les flux, moments et la densité de forces de surface, - de résolution d'une équation différentielle générale (équation 6-17) entre les flux de contraintes, la densité de forces de surface, les déplacements et les raideurs de flexion. '76 = Selon une mise en oeuvre favorable, le procédé comporte une étape d'itération, permettant de modifier les valeurs de charges appliquées, ou les valeurs dimensionnelles des panneaux, selon les résultats de l'une au moins des étapes 3 à 6. L'invention vise sous un autre aspect un produit-programme d'ordinateur comportant une série d'instructions adaptées à mettre en oeuvre un procédé tel qu'exposé, lorsque cette suite d'instructions est exécutée sur un ordinateur. Brève description des figures La description qui va suivre, donnée uniquement à titre d'exemple d'un mode de réalisation de l'invention, est faite en se référant aux figures annexées qui représentent : Figure 1 Exemple d'un panneau plat raidi par des poches triangulaires, Figure 2 ù Définition du chargement et du système de coordonnées, Figure 3 ù Définition géométrique d'un panneau Figure 4 ù Noeud dans une structure raidie par des poches triangulaires, Figure 5 ù Exemple d'instabilité générale d'un panneau raidi par des poches triangulaires, Figure 6 ù Théorie de la largeur efficace, Figure 7 ù Organigramme général du procédé selon l'invention, Figure 8 ù Décomposition de la grille en triangles élémentaires, Figure 9 ù Triangle isocèle élémentaire utilisé dans le calcul de masse du panneau, Figure 10 ù Triangle rectangle élémentaire utilisé dans le calcul de masse du panneau, Figure 11 ù Forme élémentaire de grille de raidisseurs dans un panneau raidi par des poches triangulaires, Figure 12 ù Cas de chargements purs de la plaque raidie, Figure 13 ù Schématisation de charges sur un raidisseur, Figure 14 ù Expression des coefficients Kc en fonction des cas de conditions aux limites, Figure 15 ù Panneau de raidisseurs considéré comme assemblage de deux plaques orthotropes, Figure 16 ù Charges sur une forme élémentaire de grille de raidisseurs pour panneau raidi par des poches triangulaires, Figure 17 ù Procédé de calcul de charges appliquées corrigées de la plasticité, Figure 18 ù Conventions de notation du triangle isocèle élémentaire, Figure 19 ù Interpolation linéaire ou quadratique du coefficient K, Figure 20 ù Cas de chargement combiné, Figure 21 ù Conventions de flux et moments, Figure 22 - Valeur de h(â) en fonction des différentes conditions aux limites, pour un cas de compression, Figure 23 ù Coefficient de flambage en cisaillement pour une configuration à quatre bords en appui simple, Figure 24 ù Tableau de valeurs du coefficient de flambage en cisaillement, Figure 25 ù Coefficient de flambage en cisaillement pour une configuration à quatre bords encastrés, Figure 26 ù Evolution de la constante K1 selon l'angle isocèle pour une plaque triangulaire en appui simple, Figure 27 ù Evolution de la constante K2 selon l'angle isocèle pour une plaque triangulaire en appui simple, Figure 28 ù Evolution de la constante K1 selon l'angle isocèle pour une plaque triangulaire encastrée, Figure 29 - Evolution de la constante K2 selon l'angle isocèle pour une plaque triangulaire encastrée. Description détaillée d'un mode de réalisation de l'invention Le procédé d'analyse de résistance de panneau métallique raidi par des poches triangulaires, principalement plan, décrit est destiné à être mis en oeuvre sous forme de programme sur un ordinateur de type connu en soi. Le procédé est destiné à être utilisé pour une structure principalement plane (raidisseurs et peau). Le procédé décrit ici s'applique exclusivement au calcul de paramètres structurels typiques ("typical structural settings") avec les limitations suivantes : Les extrémités de la zone étudiée ne bordent pas d'ouverture. Aucun raidisseur ne déborde de la zone étudiée. Chaque coupe doit être limitée par des raidisseurs. Toutes les poches triangulaires sur la peau sont supposées avoir la même épaisseur. Tous les raidisseurs sont supposés de mêmes dimensions. Ce procédé est utilisé pour calculer des panneaux construits à partir de matériau homogène et isotrope (par exemple mais non limitativement du métal) pour lequel les courbes monotones croissantes décrivant (a, c) peuvent être idéalisées par le moyen de formules telles que R&O (voir plus loin). L'organigramme simplifié du procédé selon l'invention est illustré figure 7. Deux types d'effondrements (dont la survenue est évaluée dans les étapes 4 et 6 du procédé) peuvent se produire sur une structure raidie par des poches triangulaires: une défaillance du matériau (qui fait l'objet de l'étape 4): les contraintes appliquées atteignent les contraintes maximums admissibles du matériau (Ftu or Fm), effondrement global : un flambage général (comprenant la grille de raidisseurs) se produit sur la totalité du panneau (cette vérification fait l'objet de l'étape 6). Par ailleurs, deux types d'instabilité (objet de l'étape 5) affaiblissent la raideur globale de la structure raidie par des poches triangulaires mais ne causent pas l'effondrement global de la structure complète : Instabilité du panneau : flambage des poches triangulaires Instabilité des raidisseurs : flambage des âmes des raidisseurs Les parties flambées peuvent supporter seulement une partie de la charge qu'elles auraient supportée avant leur flambage. De ce fait, les charges appliquées sont redistribuées dans la structure. On note que dans la présente invention, les calculs post-flambage ne sont pas traités. De ce fait, les deux types de flambage cités ci-dessus sont considérés comme des modes d'effondrement. Notation et Unités Les conventions de notations et de systèmes d'axes sont explicitées dans la figure 2 Un système local de coordonnées est défini pour chaque raidisseur. Un axe X est défini dans le plan de la section droite du raidisseur, c'est l'axe sortant, dans la direction de la dimension principale du raidisseur. Un axe Z est défini comme l'axe normal au plan de la peau, dans la direction du raidisseur. Enfin, un axe Y est le troisième axe dans un système de coordonnées droit. Pour les forces et charges, un signe négatif sur une force selon l'axe X signifie une compression du raidisseur, un signe positif signifie une tension. Un moment de flexion positif induit une compression dans la peau et une tension dans les raidisseurs. On utilise les notations générales définies dans le tableau suivant. Symbole Unité Description 1 mm2 mm4 Surface Inertie Constante de torsion Raideur normale (tension/compression) d'une plaque Rigidité de flexion ("bending" en anglais) d'une plaque N/mm2 Contrainte ("stress" en anglais) Déformation ("strain" en anglais) Déformation hors du plan Facteur de correction plastique z mm Coordonnée sur l'axe z k Coefficient de flambage Caractéristiques géométriques Les caractéristiques géométriques d'un panneau considéré ici à titre d'exemple nullement limitatif sont donnés figure 3. Pour la suite de la description, on utilise plusieurs hypothèses. On suppose que l'axe Z est un plan de symétrie pour la section droite d'un raidisseur. De même, les dimensions a et h sont définies à propos de la fibre neutre d'un raidisseur. Par ailleurs, le panneau raidi par des poches triangulaires envisagé n'a pas de raidisseurs aux deux bords définis par : X=0 et X=Lx Suffixes: ........ ............ Symbole Unite ............................ g a E K mm longueur Ly mm largeur mm Longueur du côté d'un triangle Angle du triangle mm Hauteur d'un triangle (= 2 . tane ) mm Epaisseur de la peau mm Hauteur de la toile de raidisseurs mm Épaisseur de la toile de raidisseurs mm Rayon de filet de noeud mm Rayon de poche mm2 Section droite du raidisseur selon l'axe x Décalage ("offset") du panneau entre son centre de vp mm gravité et le point d'origine du système local de coordonnées Décalage du réseau de raidisseurs entre son centre de mm gravité et le point d'origine du système local de coordonnées Matériaux Symbole Unité Designation MPa Limite élastique du matériau en compression : MPa Résistance de tension ultime du matériau : MPa Résistance de cisaillement ultime du matériau : MPa Contrainte de référence : - Déformation ("strain") plastique ultime (=e%) : - Coefficient de Poisson élastique : - Coefficient de Poisson plastique (= 0.5) Ftu su cuit ve vp Coefficient de Poisson elasto-plastique MPa Module d'Young élastique en compression MPa Module d'Young élastique en tension MPa Module sécant : MPa Module Tangent : Coefficient de Ramberg Osgood (R&O) en compression MPa Module de cisaillement : MPa Module de cisaillement sécant kg/mm Densite du materiau Contraintes Symbol Unité Désignation Charge normale avec i égal à : • 0° pour la charge appliquée sur un raidisseur à 0° • x° pour la charge appliquée sur un raidisseur à x ....................... N/mm Flux M N.mm Moment de flexion : MPa Contrainte de cisaillement : 6 MPa Contrainte normale . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . : . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6` . MPa Contrainte de flambage d'un panneau i en compression z` . MPa Contrainte de flambage d'un panneau i en cisaillement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . : . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6` MPa Contraintes appliquées sur l'élément i app : : RF ù Facteur de réserve V sec Etan nec Gsec Taux de charge en compression Taux de charge en cisaillement Taux de charge dû à la pression Charge limite Charge ultime Définitions On définit pour la suite de la description les termes suivants. Dans une structure raidie par des poches triangulaires, on appelle grille le réseau complet de raidisseurs seuls On définit de même un noeud est une intersection de plusieurs raidisseurs dans une structure raidie par des poches triangulaires (voir figure 4). En pratique, il s'agit d'un élément au design complexe comportant des rayons de courbure dans les deux directions. Quand une structure (soumise uniquement à des charges dans son plan) subit des déplacements visibles significatifs transversaux au plan des charges, elle est dite en flambage ("buckle" en langue anglo-saxonne). La figure 5a illustre un tel cas d'instabilité locale de panneau raidi par des poches triangulaires. Le phénomène de flambage peut être mis en évidence en pressant les cotés opposés d'un feuille plane de carton l'un vers l'autre. Pour de petites charges, le flambage est élastique (réversible) puisque le flambage disparaît lorsque la charge est supprimée. Un flambage local (ou instabilité locale) de plaques ou coques ("shells") est indiqué par l'apparition de bosses, ondulations ou vagues et est couramment rencontré pour les plaques composants des structures minces. Lorsqu'on considère des panneaux raidis, un flambage local, par opposition au flambage général, décrit une instabilité dans laquelle le panneau entre les LL UL longerons (raidisseurs) flambe, mais les raidisseurs continuent à soutenir les panneaux et ne montrent pas de déformation significative hors du plan. La structure peut alors présenter deux états d'équilibre : • Stable : dans ce cas les déplacements augmentent de façon contrôlée lorsque les charges augmentent, c'est-à-dire que la capacité de la structure à supporter des charges supplémentaires est maintenue, ou • Instable : dans ce cas, les déformations augmentent instantanément, la capacité à supporter des charges diminue rapidement Un équilibre neutre est également possible en théorie durant le flambage, cet état est caractérisé par une augmentation de la déformation sans modification de la charge Si les déformations de flambage deviennent trop grandes, la structure s'effondre ("fails" en langue anglo-saxonne). Si un composant ou une partie d'un composant risque de subir un flambage, alors sa conception doit satisfaire à la fois les contraintes de résistance et de flambage. On nomme instabilité générale le phénomène qui apparaît lorsque les raidisseurs ne sont plus capables de contrecarrer les déplacements du panneau hors du plan lors du flambage. La figure 5b montre un exemple de flambage global en compression d'une structure raidie par des poches triangulaires, lorsque le panneau atteint son premier mode de flambage général. De ce fait, il est nécessaire de vérifier si les raidisseurs agissent comme de simples supports du panneau (en compression, en cisaillement, et en combinaison). Si cette condition n'est pas remplie, on doit supposer que l'assemblage du panneau et des raidisseurs flambe de façon globale dans un mode d'instabilité, ce qui doit être évité dans une structure destinée à un usage aéronautique. Un effondrement général (ou global) se produit lorsque la structure n'est plus capable de supporter de charge supplémentaire. On dit alors que la structure a atteint la charge d'effondrement ("failure loading" en langue anglo- saxonne) ou la charge maximale admissible ("loading capacity" en langue anglo-saxonne) L'effondrement général couvre tous les types d'effondrement : • Effondrement dû à une instabilité (instabilité générale, post- flambage...) • Effondrement dû au dépassement de la charge maximum supportable par le matériau (par exemple après un flambage local) La largeur efficace (ou largeur travaillante) de la peau d'un panneau est définie comme la portion de la peau supportée par un longeron dans une structure de panneau raidi qui ne flambe pas lorsqu'elle est soumise à une charge de compression axiale. Le flambage de la peau seule ne constitue pas un effondrement du panneau; en fait, le panneau va supporter des charges additionnelles jusqu'à la contrainte pour laquelle la colonne formée par le raidisseur et le panneau efficace ("effective panel") commence à être détruite ("fait"). Lorsque la contrainte dans le raidisseur augmente au delà de la contrainte de flambage de la peau, la peau adjacente au raidisseur supporte une contrainte additionnelle du fait du support fourni par les raidisseurs. Cependant, la contrainte au centre du panneau ne dépassera pas la contrainte initiale de flambage, quel que soit la contrainte atteinte au niveau du raidisseur. La peau est plus efficace autour de la position des raidisseurs parce qu'il existe localement un support contre le flambage. A un niveau de contrainte donné, inférieur au flambage local de la peau, la largeur efficace est égale à la largeur du panneau. La théorie de la largeur efficace est illustrée par la figure 6 Idéalisation du matériau On rappelle ici que jusqu'à la contrainte de limite élastique ("yield stress" en langue anglaise) (Fcy), la courbe contrainte-déformation ("stressstrain curve" en langue anglaise) du matériau est idéalisée par la loi connue de Ramberg et Osgood (dite "formule R&O" dans la suite de la description) : e = 6 + 0.002 • E Fcy Equation 0-1 On en déduit les expressions suivantes: Module sécant: R&Olaw >E = 1 Equation 0-2 sec / V (ncù1) Module tangent: 1 + 0.002 E F e ev 6 F cY / 1 a(e) ù ne + 1ù ne R&Olaw a(6) Esec Ec Etan ( (ncù1) 1 0.002 6 + n E F F y cy 1 Equation 0-3 1 Esec Ec Coefficient de Poisson: v=Esec .v + e Ecv, Avec vp =0.5 Equation 0-4 Il est à noter que, avec le rapport de R&O (paramètre n ou n corrigé), connu de l'homme de l'art, ces équations sont correctes uniquement dans la zone [0 ; Fcy]. Pour la suite de cette étude, cette zone doit être étendue de Fcy à Ftu. Au delà de Fcy, différentes courbes peuvent être utilisées jusqu'à la contrainte ultime, notamment : formule de R&O utilisant un coefficient n modifié, ou méthode elliptique. Dans la suite, la formule de R&O utilise un coefficient modifié. La continuité entre les deux courbes est maintenue. Le coefficient n modifié de la 20 formule R&O est calculé par : Equation 0-5 Avec : £2p = eult £1p = 0.002 62 =Ftu a = F Il est à noter que pour utiliser cette formule, les critères suivants doivent être respectés : Ftu > Fcy et su,t>0.002 Dans la méthode elliptique, au dessus de Fcy, une autre courbe est utilisée jusqu'à la contrainte ("stress") ultime: la courbe d'extension elliptique. Naturellement, la continuité entre la courbe de R&O et la courbe d'extension elliptique est assurée. Les relations contrainte-déformation de l'extension elliptique sont : 19 n= b2 Equation 0-6 avec: û ellipse paremeters \'ob D) ±D) r x.1 2D E ret-... R nE ,et Plasticité Toujours à titre de rappel, il est connu que les facteurs de correction de plasticité dépendent du type de chargement et des conditions aux limites. 20 Les facteurs de correction de plasticité pour des panneaux rectangulaires plats sont présentés dans la table 1 ci-dessous. Chargement Conditions aux Equation limites Bride ("flange") avec 2 un bord charnière E ("hinged") non = sec 1ùv, chargé (bords appuyésùlibres) h1 E 1ùvz Bride avec un bord .......................................................................... .......................................................................... .......................................................................... .......................................................................... .... Compression et fixe non chargé = 0.33+0.335. 1+3. Etan sec flexion (encastré - bords libres) • 172 Plaque avec bords .......................................................................... .......................................................................... .......................................................................... .......................................................................... .... charnières ("hinged") Etan non chargés (bords r73 =r7i • 0.5+0.25 1 + 3 appuyés) Esec .......................................................................... .......................................................................... .......................................................................... .......................................................................... .... Plaque avec bords E fixes non chargés 174 =17i • 0.352+0.324 l+3 tan (bords encastrés) Esec . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Compression Colonne = tan 175 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . toutes conditions Cisaillement Gsec ~6= 6eq =T . 1l G . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Cisaillement bords retendus i =r7i E. 0.83+0.17 tan 6e4 = . 2 Esec / Dans le cas particulier du cisaillement, la courbe contrainte-5 déformation de compression du matériau est également utilisée comme suit : • Calcul de la contrainte normale équivalente : 6eq = z \fi • Calcul des valeurs correspondantes ES et v à partir de cette contrainte : 6eq Etape 1 ù Module d'entrée de données : ciéométrie, matériau, chargement Le procédé comporte une première phase, d'entrée des données relative au panneau raidi par des poches triangulaires considéré et au chargement appliqué à ce panneau. Ces données sont entrées par des moyens connus en soi et mémorisées dans des bases de données de type également connu. Les paramètres d'entrée du procédé de calcul analytique de panneau raidi par des poches triangulaires, comportent notamment : Dimensions générales : panneau rectangulaire (dimensions: LX, Ly) Section droite des raidisseurs : dimensions de l'âme : b, d Epaisseur constante du panneau (t) Charges aux limites du panneau : NX, Ny, NXy Calcul de masse Cette partie est destinée au calcul de la masse du panneau raidi par des poches triangulaires complet, incluant la prise en compte des rayons de filet et de noeud. Cette étape de calcul de masse est indépendante du reste du procédé décrit ici. La masse est calculée de façon connue en soi à partir de la définition géométrique du panneau. Les données d'entrée pour ce procédé sont la géométrie du panneau incluant les rayons des poches et des noeuds (Rn et Rf). La donnée de sortie est la masse du panneau. La masse est calculée en additionnant la masse de la peau et des 30 longerons. Les rayons des filets entre deux longerons, et entre la peau et les longerons sont également pris en compte. Le calcul de masse est basé sur y~ G, /1+ve Esee /1ùv • 16 =ù _ G 1+v Ee ~1ùVe/ ~~ deux triangles élémentaires: triangle isocèle et triangle rectangle (voir figures 8, 9 et 10). Etape 2 û Calcul des charqes appliquées Cette étape permet de calculer les contraintes appliquées dans la peau et les raidisseurs à partir de la géométrie du panneau raidi par des poches triangulaires, et des charges externes. La méthode prend en compte une correction de plasticité des charges appliquées, réalisée par un processus itératif. Elle permet de prendre en compte le post-flambage des raidisseurs et des poches. Il s'agit d'un progrès substantiel vis à vis du NASA "Isogrid" design handbook" (NASA-CR-124075, 02/1973) en ce qu'il prend notamment en compte les points suivants : grille de raidisseurs avec 8#60, panneau raidi par des poches triangulaires considéré comme assemblage de deux plaques orthotropes. Les données d'entrée de cette étape sont : • Données géométriques : • 0 : angle de base du triangle, o a : base du triangle, o Aist: section droite du raidisseur, i = 0°, 0 or - O. o ts : épaisseur de la peau, o tg : épaisseur du panneau équivalent à la grille Données sur le matériau : o EXS , Eys : module d'Young de la peau, o GXys : module de cisaillement de la peau, o vxy5, vyxs : coefficient de Poisson de la peau, o Est : module d'Young des raidisseurs, o vst : coefficient de Poisson des raidisseurs o Données de matériau (n : coefficient de Ramberg & Osgood, Fcy, Ftu, vp,ast=0.5) • Charges appliquées à la structure (NX°, Nÿ , NXy°) Les données obtenues en sortie de cette étape sont : 30 • NXS, Nys, Nxys : flux dans la peau, • o , ayS, ti)ys : contraintes dans la peau, • aoo, 6e, 6_e : contraintes dans les raidisseurs, • Fo., Fo, F_e : charges dans les raidisseurs. Dans la suite de la description, la peau est supposée être un matériau isotrope. Le procédé fournit des entrées pour : o L'analyse de résistance (étape 4) : contraintes dans la peau et dans les raidisseurs o l'analyse de flambage de poche (étape 5.1) : contraintes dans la peau o l'analyse du flambage du raidisseur (étape 5.2) : contraintes dans les raidisseurs o l'analyse d'instabilité générale (étape 6) : contraintes dans la peau et dans les raidisseurs pour calculer la raideur en flexion du panneau raidi par des poches triangulaires. Le procédé de calcul nécessite des données d'entrée de post-flambage des raidisseurs: Aost A+est et A est et de post-flambage des poches: ts eff Le procédé prend en compte la redistribution des contraintes appliquées ("applied stresses") entre le panneau et la grille de raidisseurs due en premier lieu, au post-flambage des raidisseurs, par la définition d'une section droite efficace pour chaque type de raidisseur (0°, +0 or -0): Ao st , A+est et A est en second lieu, au post-flambage de la poche à travers une épaisseur efficace du panneau : ts eff, enfin, à la plasticité des charges externes appliquées, grâce à un procédé itératif sur les différentes propriétés du matériau : Eost E+est, E_est pour les raidisseurs et ERS, Eys et vepst pour la peau. La charge externe est supposée être dans le plan du panneau et appliquée au centre de gravité de la section : A B B C30 Nx donc : s#0 et x=0 4 {N}= [A] • {e} avec {N}= Ny Nxy En conséquence, les contraintes dans la peau ne sont pas fonction de l'épaisseur de ladite peau et de la position dans le plan. De même, les contraintes dans les raidisseurs ne sont pas fonction de la position sur la section du raidisseur, mais seulement de l'angle du raidisseur. La définition géométrique de la grille de raidisseurs utilisée pour réaliser les calculs est définie figure 11 : Pour obtenir un panneau raidi par des poches triangulaires, cette forme élémentaire est associée à la peau et est répétée autant de fois que nécessaire. De ce fait, ce procédé ne prend pas en compte le concept de géométrie aux bords. Pour chaque raidisseur, la section réelle (A1st avec i: 0°, +e ou - e) est donnée par la relation : Aist = %AZ x Ast (dans le présent exemple nullement limitatif, seul le cas %AZ =1 est envisagé). La section droite des raidisseurs comprend la section du rayon de poche (2 . R f2 1ù 4 ùj). Quelle que soit leur position sur la grille, les contraintes et déformations sont identiques pour chaque type de raidisseur (0°, +e, -e) . Pour prendre en compte la plasticité qui peut survenir dans chaque raidisseur, le module d'Young est spécifique à chaque type de raidisseur (0°, +e, -e): E0ost East E a st La matrice "matériau" E est définie par : iEoost 0 0 Eg = 0 East 0 0 0 East) Etape 3 - Calcul des charries internes 3.1 Plaque équivalente aux raidisseurs 3.1.1 Relation entre les déformations globales et les déformations des raidisseurs Les notations et conventions géométriques sont illustrées par la figure 16. On recherche la relation entre (EX, sy, Exy) et (so., se, s_e) . Les déformations générales sont définies par les formules suivantes : en, = n.e.n Donc on a pour les déformations : eo. = ioo.e.ioo = (1 0 0). eB = iB.e.ie = (cos 9 sin 9 0 e_B = = (cos B ù sin 9 0 Et finalement : ~eo. 1 0 £B = cosy 8 sine 9 2cos2 G 2sin2 G eyz 0 e zz 0/ exz 'cos 0' EyZ sin e ~zz/\ 0 ) xy exz cos9 eyy eyZ ùsin 9 e e 0 zy zz/\ 0 2sinGcose -2sinGcose eyx eYY e zY \£xYJ Equation 3-1 La matrice ci-dessus est notée Z: 1 0 0 Z= cos2e sine8 2sin0cos8 2 cos 2 0 2 sin 2 0 ù 2 sin 0 cos 0) 3.1.2 Relation entre contraintes et déformations Comme on l'a dit, les notations et conventions géométriques sont illustrées par la figure 16. Les charges dans les raidisseurs sont données par l'expression suivante : Pz = EZ.Es`.As`Z.il , (i = 0°, e, -e) Equation 3-2 Ainsi l'élément de base ci-dessus est soumis à : Poo,PB,PB (Poo compté deux fois puisque la dimension de l'élément de base selon l'axe Y est 2h, donc le raidisseur correspondant à 0 ° doit aussi être pris en compte). Selon l'axe x: Pxg = P.x = (2P0o + Pa + P a).x = 2Eoos`Aoos`Eoo + EastAastEa cos 8 + E_es`A-astE_a cos 8 Selon l'axe y: Pyg = EastAa"ea sin 9 + E_as`A_as`E_a sin Cisaillement dans le plan (x, y) : Aa stEa sine ù E_as` A_a st5_a sine Equation 3-3 Pour obtenir les contraintes, la charge est divisée par la surface d'un élément de base. La section de l'élément de base sur une surface normale selon l'axe X est 2htg = atanO.tg La section de l'élément de base sur une surface normale selon l'axe Y est de même atg. En termes de contraintes on a: g 1 2Eo0s`Ao0s` stA-st cos2 8 E stA-st cos2 8 x atg tan 8 e e sine e e sin 8 Pour ay et tixy, on obtient avec la même méthode : eo. g p g 1 (2E0ost Aoo st Ea st ast cos 8 E_a st 6x = _ A A_ast cos O) Ea 2htg 2htg eo. /e0. 6yg = 1 (0 EastAast sin 9 E_as`A_as` sin e) £e ~£-e J /soc (0 East Aast sin 9 ù E_0stA_0st sin e) £e atg Les mêmes résultats peuvent être présentés sous forme matricielle : atg 1 y g Y J / 2EoostAo0st EstAst cos2 9 tan 9 a a sin 9 0 EastAast sin 9 0 EastAast sin 9E ~stA pst cos2 Gv sin 9 E astAgst sin 9 ù E astAast sin 9 J Equation 3-4 La matrice ci-dessus est notée T: / st st 2Eoo A00 st st cos2 9 st st cos2 9 tan 8 EgA9 sin 8 EgAg _ sin 9 T = 1 0 EastAast sin 9 E_astAast sin 9 atg 0 EastAast sin 9 ù E_0 st Aa st sin 9 En utilisant alors l'Equation 3-1 et la notation de la matrice Z : / g 6x £x 6yg = T.Z. £y g ~z J £xy J Equation 3-5 La matrice qui précède est notée W: W =T.Z Cette relation (Equation 3-5) signifie que le comportement du panneau équivalent aux raidisseurs est similaire à un matériau anisotrope (la matrice W peut être complète : toutes ses cellules ont des valeurs non nulles). 3.2 Panneaux raidis On utilise l'hypothèse de Kirchhoff : des sections planes restent planes après déformation. Le réseau de raidisseurs est modélisé par un panneau équivalent avec une matrice de comportement W (voir l'Equation 3-5). Cette 28 disposition de modélisation du panneau raidi par des poches triangulaires, par deux plaques orthotropes est illustrée par la figure 15. Pour le calcul de vyx, on a: xys = yx Ex Eys 3.2.1 Contraintes et charges des raidisseurs Flux dans le panneau équivalent aux raidisseurs L'expression générale des flux est : Naa Equation 3-6 Le flux selon l'axe X s'exprime comme suit : s v s h --+t f 2 h 2 1ùvxysv en utilisant l'Equation 3ù5: h sdz+f h xxgdz +ts 2 Exsts +t .W1,1 1ùvxYsV s g vsEst xy y s 1ù v. s v, s + tg .W 1,2 y + tg .W 1,3 £xY Nx = .£x + Et, en utilisant la même méthode pour les flux Ny et Nxy: .£x + sE~st Vyxvs == s +tg.W2,1 xy yxEyS S + tg .W 2,2 1ùv v xy yx .£y + tg .W2,3 .exY = tg .W3,2 .£x + tg .W3,3 .£y + (2Gts + tg .W3,1 £xy Equation 3-7 Ces expressions montrent clairement la répartition du flux entre la peau et le panneau équivalent aux raidisseurs. Dans la peau, la relation entre les flux 20 et déformations est : i s Ex ts v , SE Y sts /Ns~ iN Ng x x x NY s g = NY NY N g ~NxY / N } 1ùvxysvyxs v SEst x s 1_v S S xY Est y s 0 Equation 3-8 On note X la matrice précédente : f Exsts 1 ù vx,,svys v sEst X sVsEst XY Y s 1û vxYsvys Est y s 0 0 X= 10 1 ù vxYsVYxs 1 ù 0 0 2Gx,Sts Alors : ~N /Ng` /Ng x x x Ng + 1 .X.W . Y N g tg J Equation 3-9 En inversant cette relation, les flux dans la grille sont exprimés en fonction des flux appliqués globaux : N g' f-N x x NYg =V. NY g Nxy / Nxy/ Equation 3-10 NY = NxY Ng Y J 0 0 £Y 2GxYsts J Avec: V = 1 =_ù1 Id +ù.X.W tg (Id est la matrice identité) Contraintes et charges dans les raidisseurs 15 Les flux dans le panneau équivalent aux raidisseurs peuvent être exprimés par : Ngv £0 x 0 Nyg = tg.T. ea N g xy % £ùBJ Equation 3-11 /Ng x Ng Y ==1 = tg .T.Eg . Y J Equation 3-12 En utilisant la notation suivante : U = tg.T.Eg On a : 600 / Nxg \ 6B = U-1 Ng y x Equation 3-13 Finalement, les charges et contraintes dans les raidisseurs s'expriment en fonction des flux de charges externes : /Nx v /Fov / e st 6 v/6v 00 Ao. 0o 6B = U 1.V. Ny et Fg ABst 6e st ~6 e/ ~Nxy~ ~F e/ A-e } 6-a) Equation 3-14 3.2.2 Flux et contraintes dans la peau Comme montré par l'Equation 3-8, les flux dans la peau s'expriment 15 comme ci-dessous : /NS\ x Nys N s xy //N /Ng\ x x Ny NY g y g N N XY i ~xy i Equation 3-15 ainsi, les contraintes dans la peau sont exprimées par : s 6x Nx s 1 s 6Y ù tY ~xyJ J Equation 3-16 3.3 Le procédé de calcul de charges appliquées corrigées de la plasticité est présenté ici avec les références aux matrices introduites dans la description (figure 17). On note que la théorie utilisée pour calculer la plasticité suppose l'isotropie du matériau formant la peau. La solution pour les charges appliquées corrigées pour la plasticité est fournie par une méthode itérative. Un procédé de convergence doit être réalisé jusqu'à ce que les cinq paramètres du matériau (Eoost ,E+est, E-est, Epeau, vep) entrés au début du processus itératif soient égaux aux mêmes paramètres calculés en sortie (après calcul de contrainte plastique). Dans la figure 17, déjà citée, les paramètres de convergence sont indiqués par un fond grisé. Plus précisément, dans ce processus itératif, les entrées initiales sont les charges appliquées dans la grille des raidisseurs, et dans la peau. Pour la grille des raidisseurs, les données de la ième itération de Module d'Young des raidisseurs Eoô t, E+est, E_est en 3 directions : 0°, +0, - 0, permettent, en association avec la valeur de l'angle 0 et de la géométrie, de calculer la matrice [T] (équation 3-4). Les valeurs d'angle 0 et de géométrie fournissent la matrice [Z] (équation 3-1). Les matrices [T] et [Z] donnent la matrice [W] (équation 3-5). Pour la peau, les données de matériau de peau isotrope Epeau, Vep) permettent de calculer la matrice [X] (équation 3-8). Les matrices [W] et [X] permettent de calculer les matrices [U] (équation 3-13) et [V] (équation 3-10). Les résultats tirés de ces matrices comprennent : les flux, les efforts élastiques, les corrections de plasticité sur les efforts, les valeurs de la i+1 ème itération des Modules d'Young des raidisseurs et de la peau corrigés, et du coefficient de Poisson de la peau corrigé, et les charges dans les raidisseurs. On comprend que le calcul est itéré jusqu'à ce que les modules d'Young et coefficients de Poisson varient au cours d'une itération d'une valeur inférieure à un seuil prédéterminé. Impact de la correction de plasticité sur le calcul de flambage général : Naturellement, la correction de plasticité altère la matrice de loi de 5 comportement calculée dans le module d'instabilité générale à travers les 5 paramètres de matériau (voir la partie sur l'instabilité générale). / Nxx "cc' £ B i£o 411 `412 A13 B11 B12 13 .)CC 0 A B £yy0 Al2 A22 A23 B12 B22 B23 eYY 0 0 Nxy £xy A13 A23 `411 B13 B23 B33 -M Kxx = B11 B12 B13 C11 C12 C13 K.x Myy B C Kyy B12 B11 B23 C12 C22 C23 KYY ù Mxy % J \Kxy % \B13 B23 B33 C13 C23 C33 } Kxy Equation 3-17 De ce fait, la correction de plasticité altère également les coefficients ni 10 (i=1..3) utilisés pour le calcul du flambage général. La correction de plasticité dans l'analyse de flambage général est fournie par ces coefficients modifiés 52;. 3.4 Exemple: distribution de charges sous bi-compression et cisaillement 15 Dans l'étape de calcul de contraintes appliquées, le rayon de filet de la poche est pris en compte pour le calcul de la section du raidisseur. Par ailleurs, il n'y a pas de post-flambage, on peut donc écrire : %Ao0 f = %A+est = %A_esf = 1000/0 ts_eff = ts 20 Dans le présent exemple décrit ici à titre nullement limitatif, la géométrie du panneau raidi par des poches triangulaires est définie par : LX = 1400.45 mm a = 198 mm rayon de noeud : Rn = 9 mm Ly = 685.8 mm t = 3.64 mm rayon de poche : Rf= 4 mm 8=58° b=2.5mm d=37.36mm 25 On considère un matériau isotrope. La loi élasto-plastique utilisée est celle de Ramberg & Osgood. Les charges externes appliquées (à charge ultime) sont : Nx = -524.65N/mm Ny = -253.87N/mm Nxy = 327.44N/mm Le procédé de calcul de charges internes et de charges appliquées 5 incluant la correction de plasticité s'écrit sous forme matricielle : w = T7 1518 504 0 epsx 504 1290 0 epsy 0 0 504 gamxy X 312000 93600 0 epsx 93600 312000 0 epsy 0 0 109200 gsmxy 0.152 0.005 0.000 Nx 0.011 0.130 0.000 Ny 0.000 0.000 0.147 Nxy u-1 On obtient alors comme résultats les flux, contraintes et charges dans 10 les raidisseurs : Nxg = -81.1N/mm 600 = -101.14MPa Foo = -10441N Nyg = -39N/mm 6+0 = 44.25MPa F+B = 4437N Nxyg = 48.1N/mm 6 e = -135.07MPa F_0 = -13537N Ainsi que les flux et contraintes dans la peau : .......................... ........................ sigxg sigyg tiuxyg Nx-Nxg Ny-Nyg Nxy-Nxyg Nxq Nxyg .......... ............ Nxs = -443.5N/mm Nys = -214.87N/mm Nxys = 279.3N/mm = -121.85MPa = -59.03MPa = 76.74MPa o-xs o-ys 15 Dans cet exemple, les contraintes appliquées demeurent dans le domaine élastique. La distribution de flux entre la peau et la grille de raidisseurs est résumée par le tableau suivant : Flux Flux Distribution de charges - Distribution de charges - (N / mm) externes Flux Pourcentages Grille Peau Grille Peau Nx -524.65 -81.13 -443.52 15.46% 84.54% Ny -253.87 -39.00 -214.87 15.36% 84.64% Nxy +327.44 +48.12 +279.32 14.7% 85.3% Etape 4 ù Module d'analyse de résistance : Cette phase a pour but le calcul de facteurs de réserve (RF) par comparaison des charges appliquées calculées dans les composants du panneau raidi par des poches triangulaires, et les contraintes maximums 10 admissibles du matériau. Les contraintes admissibles dans le matériau à charge ultime sont déterminés par : Ftu (Résistance de tension ultime du matériau) comparé aux contraintes appliquées dans les âmes des raidisseurs, Ftu comparé aux contraintes principales appliquées dans la peau, Fsu (Résistance de 15 cisaillement ultime du matériau) comparé à la contrainte de cisaillement maximal dans la peau. L'analyse de résistance consiste à calculer les facteurs de réserve du matériau à la charge limite et à la charge ultime. Les charges appliquées proviennent des charges dans le plan (compression, cisaillement) ou hors du 20 plan (pression). Les données d'entrée pour ce calcul sont : • Valeurs admissible pour le matériau : Fty, Fcy, Fsy, Ftu, Fsu • Contraintes appliquées à la structure : o Contraintes sur la peau (axs, ays et tixys) 25 o Contraintes normales dans les raidisseurs Note: les contraintes appliquées sont corrigées pour la plasticité dans le procédé de calcul de contrainte appliquée, comme on l'a dit plus haut. Les données de sortie sont les facteurs de réserve. L'analyse de résistance dans le plan du panneau se base sur les hypothèses suivantes. Les contraintes dans la peau ne sont pas une fonction de l'épaisseur de la peau et de la position dans le plan. Les contraintes dans un raidisseur ne sont pas fonction de la position sur la section du raidisseur, mais seulement de l'angle du raidisseur. Ces hypothèses ne sont plus valables lorsque le post-flambage et le comportement en dehors du plan sont également pris en compte. Dans ces cas, des fonctions max/min, de type connu en soi, doivent être mise en oeuvre pour tenir compte de ces phénomènes. Calcul des contraintes principales : Pour calculer le facteur de réserve de la peau, les contraintes principales (6max, Emin et tmax) sont utilisées : z 6xs + 6ys + 6xs ûa 2 2 6max s 2 + 2xys Emin s 2 6xs + 6ys _ 6xs ù am Z 2 2 + 2xys 6max s ù Emin s 2max 2 Equation 4-1 La valeur 6max utilisée dans le calcul du facteur de réserve est définie comme le maximum absolu entre 6max_s et 6min_s calculés dans l'Equation 4-1. Facteur de réserve à charge limite (LL) : Facteur de réserve sur les âmes des raidisseurs : F Y blade 0 in plane LL ablade 0 LL blade _+B _ in plane LL F Y ablade +0 LL blade _ùB _ in plane LL F Y ablade -0 LL 10 36 Facteur de réserve sur la peau : Cisaillement maximal : RFinkipla eeazL û 2max LL Dans cette formule, si Fsy n'est pas connu, Fm / -q3 peut être utilisé RF = m111{RFblade_0 RFblade+B RFbladeùB RF skin shear • RF skin { material_LL in plane LL' in plane LL' in plane LL' in_plane _LL' in_plane _LL Equation 4-2 Facteur de réserve à charge ultime UL : • Facteur de réserve sur les âmes des raidisseurs: Contrainte principale : skin _ y F RF in plane LL a max LL • Un facteur de réserve enveloppe est calculé à charge limite RF _ blade0 _ l~ m plane UL Fa ablade 0 UL 15 blade+B = Fu in _ plane _ UL ablade +B UL bladeùB _ in _ plane _ UL ablade ùB UL Note: si Fcu n'est pas connu, Fcy ou Ftu peut être utilisé • Facteur de réserve sur la peau: Cisaillement maximal: RF skin _shear = Fsu in _ plane _ UL 2max UL 20 Contrainte principale : RF.skin = in_plane _UL Fa 6max UL • Un facteur de réserve enveloppe est calculé à charge ultime : RF D blade0 DD blade+B DD blade_ùB DD skin shear DD skin materialUL = min RinplaneUL, in_plane_UL' in_plane_UL; plane UL n_plane_UL L'analyse de résistance hors du plan du panneau sort du cadre de la présente invention. Dans un exemple de mise en oeuvre de cette partie du procédé, d'analyse de résistance dans le plan, la même géométrie et le même cas de charges que dans les parties précédentes est étudié. Dans un premier temps, le ratio de facteur de charges est supposé égal à 1. Facteur de réserve des âmes des raidisseurs : RFblade0 in ù 4.71 plane UL Dblade_+B _ Rl in_planeUL -11'14 bladeùB = 3.63 inplane UL Facteur de réserve de la peau : Calcul des contraintes principales: !17 6max s ù~7.52MPa am. s p -1 / 3.35M17a Z = 82.91MPa max RF skin shear ù 3.41 in plane UL RFinaterial UL 2.82 20 Dans le cas de calcul plastique, un calcul itératif est réalisé jusqu'à ce que les charges appliquées atteignent la charge de résistance d'effondrement (de manière à effectuer sur elles une correction de plasticité). A chaque boucle d'itération, les mêmes calculs que décrits précédemment sont réalisés. 25 Etape 5 ù Calcul de contraintes locales admissibles Deux types d'instabilité affaiblissent la raideur globale de la structure raidie par des poches triangulaires, mais ne causent pas l'effondrement global de la structure complète Instabilité du panneau : flambage des poches triangulaires 15 amax skip UL = 1 / 3.35M1 a Cisaillement maximal : Contrainte : RF iSkin =2.82 principale nplane UL Facteur de réserve enveloppe RF à charge ultime (UL) Instabilité des raidisseurs : flambage des âmes des raidisseurs Les parties flambées peuvent supporter seulement une partie de la charge qu'elles auraient supportée avant leur flambage. De ce fait, les charges appliquées sont redistribuées dans la structure. On note que dans la présente invention, les calculs post-flambage ne sont pas traités. De ce fait, les deux types de flambage cités ci-dessus sont considérés comme des modes d'effondrement. 5A ù Calcul de flambage local du panneau : Dans les panneaux raidis par des poches triangulaires, les poches sont des plaques triangulaires soumises à des chargements combinés dans le plan. De manière à calculer des flux de flambage admissible et des facteurs de réserve sous un chargement pur du panneau raidi complet, une méthode basée sur un modèle par éléments finis (FEM) est utilisée. Cette partie fournit un calcul de flux de flambage de poche admissible pour des plaques triangulaires isocèles : la base du triangle peut varier entre toutes valeurs, et, dans le présent exemple non limitatif, entre 45° et 70°. Le calcul de flux est réalisé avec deux types de conditions aux limites : en appui simple ("simply supported,") et encastrées ("clamped") Les contraintes appliquées à prendre en compte pour le calcul de facteur de réserve sont les contraintes agissant dans la peau uniquement, et déterminés dans la partie, décrite plus haut, de calcul de contraintes appliquées. Les données d'entrée de cette partie sont : • données géométriques (base du triangle, angle isocèle, épaisseur de la peau) • données du matériau (linéaire (E, v) et non-linéaire (Fcy, tu,e%,n~ o matériau isotrope uniquement o les valeurs de flux de flambage plastique admissible sont pertinentes jusqu'à Fcy seulement • Conditions aux limites: en appui simple ou encastrées • Charges appliquées à la peau Il est à noter que tous les flux externes utilisés dans cette partie sont des flux de peau et ne correspondent pas au chargement complet du panneau raidi. Par ailleurs, la hauteur du triangle (h) utilisée comme longueur de référence pour le calcul de flambage est réduite de la demi-épaisseur des âmes des raidisseurs. Dans cette partie, la formule utilisée pour la hauteur du triangle est : h=h red h red = a _Red • tan e 2 a red =aû b• 1 1 tane sine Les données de sortie sont les suivantes : • Flambage de poche admissible. • Facteur de réserve Cette partie a pour but de calculer des flux admissibles pour des plaques triangulaires isocèles. Elle comporte deux parties : 1/ Calcul de valeurs admissibles pour des plaques triangulaires en appui simple (partie 5A.5), 2/ Calcul de valeurs admissibles pour des plaques triangulaires encastrées (partie 5A.6). Les deux parties suivent la même démarche : d'abord un calcul de valeurs admissibles pour des plaques triangulaires soumises à des cas de chargement pur (compression selon X, compression selon Y et cisaillement), puis un calcul des courbes d'interaction entre les trois cas de chargement pur 5A.1 Principe de calcul Les cas de chargement pur envisagés sont présentés dans la figure 12. Il existe dans la littérature plusieurs méthodes de formulation analytique du flambage local de triangles. La comparaison entre ces méthodes montre de grandes différences entre les contraintes de flambage prédites. Plus encore, certains des paramètres utilisés dans ces méthodes sont dérivés de calculs par éléments finis, les tests sont souvent empiriques et certaines méthodes ne fournissent aucune donnée pour des angles différents de 60° Le développement d'une théorie complète de ce problème étant long, le procédé tel que décrit ici à titre nullement limitatif met en oeuvre une méthode basée entièrement sur un Modèle Eléments Finis (FEM): • Création d'un modèle paramétrique FEM d'une plaque triangulaire (paramètres: angle de base, épaisseur, hauteur du triangle, conditions aux limites), • Exécution de nombreuses combinaisons différentes pour 10 obtenir des résultats de flambage linéaires, • Obtention de paramètres qui puissent être utilisés dans une formulation analytique (coefficients K). Les effets de plasticité induits doivent également être pris en compte dans le calcul des valeurs admissibles. Les charges appliquées sont soit des 15 charges simples, soit des combinaisons de ces charges simples. 5A.2 Cas de chargement pur Des courbes d'interaction sont définies comme suit. Six modèles éléments finis de plaques triangulaires ont été créés avec des angles compris 20 entre 45° et 70°, dans le présent exemple nullement limitatif en ce qui concerne les angles du triangle. Dans cette partie, l'angle isocèle (8) est défini comme l'angle de base du triangle isocèle (voir figure 18). Pour chaque angle isocèle et pour chaque cas de chargement pur, l'étude est organisée en trois points : 1/ Calcul par Modèle Eléments Finis (FEM) 25 Des calculs linéaires de flambage local des triangles par modèle éléments finis (FEM) de type connu, ont été réalisés pour déterminer le flux admissible de plissement (sans correction plastique) pour diverses épaisseurs de plaque et donc diverses rigidités. On note que le premier mode observé présente toujours un seul flambage (une seule bosse). 30 2/ tracé de la courbe de flux de flambage admissible en fonction de D 41 De façon générale, dans la littérature, le flux de flambage admissible est exprimé comme suit : Ncrit = K.- Kest une constante, 121-v2)' h est la hauteur du triangle : h= 2.tane. Mais cette étude démontre que, dans le cas de plaques triangulaires, et pour de petites valeurs du rapport D , il n'est pas pertinent d'exprimer le flux de flambage admissible avec une équation au premier degré en D . Des 10 résultats meilleurs sont obtenus avec une équation au second degré de la forme suivante : Nerit = K~ ~ +K2. D (valeur élastique) Equation 5-1 Les constantes K1 et K2 dépendent de l'angle et du cas de charge 15 considéré. Donc, pour chaque cas et chaque angle, une valeur des constantes K1 et K2 est obtenue. 3/ tracé de l'évolution de K, et K2 en fonction de l'angle de base du triangle isocèle K1 et K2 sont tracés en fonction de l'angle et une interpolation est 20 réalisée pour déterminer une équation polynomiale qui permet de calculer ces constantes quel que soit l'angle compris entre 45° et 70°. La figure 19 illustre une interpolation linéaire ou quadratique pour les coefficients K. Il est clair que cette fonction permet également d'extrapoler à des valeurs en dehors mais proches du domaine allant de 45 à 70°. Ainsi, en connaissant l'angle isocèle et 25 les conditions aux limites, il est possible de calculer directement le flux de flambage admissible de la plaque triangulaire étudiée. D la raideur de plaque : D = E't3 5A.3 Cas de chargement combiné On utilise dans ce cas l'hypothèse suivante : si certains composants de la charge combinée sont en tension, alors ces composants sont ramenés à zéro (ne sont pas pris en compte pour le calcul). Il est en effet conservatif de considérer que des composants en tension n'ont pas d'effet vis à vis du flux de flambage étudié, et n'améliorent pas la contrainte de flambage de la plaque. Par exemple, si Nx =+200N/mm (qui met en évidence une tension) et Nsapp=300N/mm, la charge combinée admissible est réduite à la charge admissible de cisaillement pur. La présentation des cas de chargements envisagés est illustrée figure 21. Dans cette partie, trois modèles d'éléments finis ont été utilisés : trois plaques triangulaires isocèles avec des angles égaux à 45°, 60° et 70°. Pour chaque angle, l'étude est organisée en deux points : 1/ Calcul par Modèle Eléments Finis (FEM) Pour toutes les combinaisons présentées ci-dessus, des calculs linéaires par méthode d'éléments finis ont été réalisés pour déterminer la valeur propre ("eigenvalue") de flambage correspondant à différentes répartitions de charges externes. On constate que le premier mode observé présente toujours un seul cloquage. Il apparaît que la courbe d'interaction est peu dépendante de la valeur de D. 2/Tracé des courbes d'interaction Pour chaque angle et chaque combinaison de charges, la courbe d'interaction est tracée. Ensuite ces différentes courbes ont été approximées avec des courbes classiques dont l'équation a la forme suivante : RIA +R2B =1 N ' Avec Ri = `rNä , i = cX, cY ou s. Les résultats et les choix effectués ont montré que les équations des courbes d'interaction ne dépendent pas de l'angle de base du triangle isocèle et sont donc compatibles et peuvent être unifiées avec une unique équation couvrant toutes les combinaisons avec la forme suivante : RDA+RcyB +Rsc = 1 A partir de cette équation, pour déterminer les facteurs de réserve, on résout l'équation suivante : R A ,R B R cY + + s =1 R R R N aPP N aPP N aPP 1 avec R = cY = cX = s = N cric N cric N cric RF cYcomb cXcomb scomb 5A.4 Facteur de correction de plasticité L'obtention d'une correction de plasticité pour les cas de chargement pur, fonction de l'angle isocèle et des conditions aux limites est très complexe. En fait, pour des plaques triangulaires, les fonctions de déflexion sont complexes et de nombreux problèmes d'intégration numériques se posent. En conséquence, il est décidé d'utiliser un facteur rt conservatif, sur la base du document NACA Report 898 ("A Unified Theory of Plastic Buckling of Columns and Plates", juillet 1947). Ce facteur est défini pour tous les cas de chargement (pur et combiné) à l'exception du cisaillement, par : Et, pour les cas de chargement en cisaillement pur, par : (1 + V ) Esec (1+v) Ec La correction est calculée en utilisant la contrainte élastique équivalente de Von Mises : _ crit crit crit crit 2 crit 6VM 6xcomb + aycomb ù 6xcomb 6y_comb + 3 xycomb Ainsi, les contraintes admissibles corrigées peuvent être calculées : Etan plast crit x comb ù '! 6x comb plast y~ crat 6y _ comb = '! 6y comb plast ù crit Zxy comb ù '! Zxy comb Pour les cas de chargement pur (compression selon X, compression selon Y ou cisaillement), la correction de plasticité est également appliquée à la contrainte de Von Mises, donc pour les cas de charge de cisaillement pur, la contrainte corrigée est :: zxV . 5A.5 Plaques triangulaires isocèles en appui simple ("simply supported") Cas de chargement pur Dans les figures 26 et 27 sont présentées les évolutions des constantes K1 et K2 en fonction de l'angle isocèle. Les équations de ces courbes sont (avec e en degrés): Klex = -0.0000002417 •94 + 0.0000504863 •93 ù 0.0039782194 .92 + 0.1393226958 •9 -1.8379213492 K1 = -0.0000007200 •94 + 0.0001511407 •93 -0.0119247778 .92 + 0.4177844180 .0ù5.4796530159 Kas = -0.0000018083.94 +0.0003804181. 0' -0.0300743972 .02 +1.0554840265.0ù13.8695053175 K2,X = 0.0029565 .93 -0.4291321. 02 + 21.1697836 •9 ù 291.6730902 . = 0.0068664 . 93 -1.0113413 . 92 +51.3462358 • 9 -852.1945224 K2, = 0.013637 •93 ù 2.017207 .02 +102 .120039 .0ù1674 .287384 Que ce soit pour K1 ou K2, leurs valeurs en compression pure X et pure Y sont égales pour un angle isocèle de 60°. Le point d'intersection à 60° témoigne du comportement isotrope de la structure raidie par des poches triangulaires à 60°, en termes de flambage local de la peau. Cas de chargement combiné On réalise dans les cas de chargement combiné une analyse par modèle éléments finis de calcul linéaire de flambage pour des plaques triangulaires en appui simple. On choisit alors une courbe d'interaction conservative, proche des courbes d'interaction calculées, mais de formulation simple, qui devient la courbe utilisée dans le procédé décrit ici. Son équation N apP est: RIA + R2B =1. Avec Ri = r.t , i = cX, cY ou s. Interaction compression X + compression Y (cas 1) Dans ce cas de chargement, pour les angles compris entre 45° et 70°, on choisit une courbe d'interaction conservative. C'est-à-dire une courbe déclarant une interaction à une valeur inférieure de la somme des compressions Rcx sur X et Rcy sur Y, vis-à-vis de toutes les courbes d'interaction calculées pour des valeurs angles compris entre 45° et 70°. Cette courbe est définie par l'équation suivante : Re), =1 Interaction compression X + cisaillement (cas 2) Dans ce cas de chargement, pour les angles compris entre 45° et 70°, on choisit une courbe d'interaction conservative, vis à vis des différentes courbes d'interaction selon les angles compris entre 45° et 70°, définie par l'équation suivante : Re, + R, Î =1 Interaction compression Y + cisaillement (cas 3) Pour les angles compris entre 45° et 70°, pour déterminer le facteur de réserve dans le cas de chargement combiné de compression selon Y et de cisaillement, on choisit une équation conservative de formule : Re), +RS2 =1 . De manière à aboutir à une seule équation couvrant tous les cas de chargement, on choisit d'utiliser une autre courbe encore plus conservative Re), + Rs z =1 Interaction compression X + compression Y + cisaillement (cas 4) L'équation choisie pour ce cas de charge est : Rx+Re), +Rsz =1. Cette unique équation est utilisée pour tous les cas de chargement combiné 5A.6 Plaques triangulaires isocèles encastrées Cas de chargement pur Dans les figures 28 et 29 sont présentées les évolutions des constantes K1 et K2 en fonction de l'angle isocèle. Les équations de ces courbes sont (avec e en degrés): Klex = -0.0000018547.94 +0.0003940252 .0' -0.0314632778 .02 +1.1143937831.0ù14.8040153968 K1 = -0.0000027267 •94 + 0.0005734489 •93 ù 0.0453489667 •92 +1.5921323016 •9 -20.9299676191 K1 = -0.0000069990 .94 +0.0014822211 .93 -0.1179080417 •92 + 4.1617623127 .9 -54.9899559524 K2,X = 0.0110488.93 -1.6258419.92 +81.8278420.0ù1254.6580819 K2,P = 0.0158563 •93 ù 2.3439723 •92 +119.5038876 .0ù1970.9532998 K2, = 0.0252562 .93 -3.7563673 .02 +191.0642156 .0ù3113.4527806 Cas de chargement combiné On réalise dans les cas de chargement combiné une analyse par modèle éléments finis de calcul linéaire de flambage pour des plaques triangulaires encastrées. On choisit alors une courbe d'interaction conservative, proche des courbes d'interaction calculées, mais de formulation simple, qui devient la courbe utilisée dans le procédé décrit ici. Son équation est : RDA +R2B =1, avec Ri = , i = cX, cY ou s. Interaction compression X + compression Y (cas 1) Dans ce cas de chargement, pour les angles compris entre 45° et 70°, on choisit une courbe d'interaction conservative, vis à vis des différentes courbes d'interaction selon les angles compris entre 45° et 70°, définie par l'équation suivante : R~X + Re), =1 Interaction compression X + cisaillement (cas 2) Dans ce cas de chargement, pour les angles compris entre 45° et 70°, on choisit une courbe d'interaction conservative, vis à vis des différentes courbes d'interaction selon les angles compris entre 45° et 70°, définie par l'équation suivante : RX + RS z =1 Interaction compression Y + cisaillement (cas 3) Pour les angles compris entre 45° et 70°, pour déterminer le facteur de réserve dans le cas de chargement combiné de compression selon Y et de cisaillement, on choisit une équation conservative de formule : Re), +RS2 =1 . De manière à aboutir à une seule équation couvrant tous les cas de chargement, on choisit d'utiliser une autre courbe encore plus conservative 3 Re), +Rsz =1 . Interaction : compression X + compression Y + cisaillement (cas 4) L'équation, utilisée pour ce cas de charge, est Rex+Rey+Rs =1 Cette unique équation est utilisée pour tous les cas de chargement combiné. 5B Calcul de flambage local de raidisseur : Ce module calcule la contrainte de flambage et le facteur de réserve de l'âme ("blade") du raidisseur, considéré comme un panneau rectangulaire avec diverses conditions aux limites à définir par l'utilisateur. Les contraintes appliquées à prendre en compte pour les calculs de facteur de réserve sont les contraintes dans les âmes des raidisseurs seulement, provenant du module de calcul des contraintes appliquées. Sur la grille de raidisseurs, un ou plusieurs types d'âmes des raidisseurs sont chargés en compression. De ce fait, la contrainte de compression admissible doit être calculée. Les données d'entrée de ce module sont : • Données géométriques : dimensions des âmes des raidisseurs (longueur, hauteur, épaisseur), • Données de matériau (linéaire (E, v) et non-linéaire (Fby, Ftu, suit, ne)) . On considère dans le présent exemple un matériau isotrope seulement, • conditions aux limites (quatre sont disponibles), • charges appliquées sur les âmes des raidisseurs. Les données de sortie sont le flambage admissible de l'âme du raidisseur et un facteur de réserve. La contrainte de flambage admissible de l'âme du raidisseur est (voir figure 13 pour les conventions de notation) : . z z ablade k E~ ~ b blade c 1 2.( 1û) Equation 5-2 Avec: b: épaisseur de l'âme des raidisseurs d: hauteur de l'âme des raidisseurs Lb: longueur de l'âme des raidisseurs Eb: Module d'Young en compression ve: Coefficient de Poisson dans le domaine élastique kb: Facteur de flambage local (dépendant des conditions aux limites et de la géométrie) n: Facteur de correction de plasticité 25 30 Note: la longueur d'une âme de raidisseur est donnée par (Lb) = a (pour les âmes des raidisseurs dans la direction X) (Lb) = 2 . J1+tant 9 (pour les âmes des raidisseurs dans les directions transverses) De nombreuses conditions aux limites peuvent être appliquées sur l'âme de raidisseur en fonction de la structure environnante. (voir figure 14). On note que si Lb/d est supérieur à la valeur Lim, alors kc vaut kc infini. Le facteur de flambage conservatif recommandé pour les calculs basés sur de nombreuses analyses par éléments finis est le cas 2 (2 bords encastrés û 1 bord en appui simple û 1 bord libre). Selon les conditions aux limites citées plus haut et selon la table 1 de facteur de correction plastique pour des plaques rectangulaires, le facteur de correction de plasticité utilisé dans ce cas est : 2 ~_= ECC 1ûv2 La formule pour le calcul de facteur de réserve pour le flambage de l'âme de raidisseur est valide pour tous les types d'âmes de raidisseurs utilisés dans le présent panneau raidi par des poches triangulaires (0°, +0, -0): crit RF blade = 6blade buck 6 app blade L'exemple suivant est basé sur la même géométrie que dans les parties précédentes. La géométrie des raidisseurs est : b = 2.5 mm, et d = 37.36 mm. La longueur des âmes des raidisseurs est (Lb) : Lb = a = 198mm pour les âmes des raidisseurs dans la direction X et Lb = 2..I1+tan29 = 25 186.82 mm pour les âmes des raidisseurs transverses. Les conditions aux limites utilisées sont : 2 bords encastrés û 1 bord en appui simple û 1 bord libre. 20 5 d s\ 2 Ainsi = 4.143 • ù +0.384 = 0.5315 pour les âmes des Lb ids\ 2 raidisseurs dans la direction X , et k, e° = 4.143. ù + 0.384 = 0.5497 pour les Lb âmes des raidisseurs transverses Et la charge de flambage pour chaque raidisseur est : 2 6biadeo° ke 0° E = e b =167.78MPa 12.1ùve) vd Au flambage de l'âme de raidisseur, le facteur de correction de plasticité est : 2 2 7J.E'e .Tt- b J 6biadee° = kea° 12 1û v2 d = 173.52MPa 2 Esec 1 ùVe 17 _= E 1ûv2 c=1 (élastique) 10 Les charges appliquées dans les raidisseurs sont : 60° = -101.14MPa 6+0 = 44.25MPa 6 e = -135.07MPa Les résultats du calcul de facteur de réserve sont les suivants : crit RF blade _ 0° _ ablade _ 0° =1.61 buck 6app blade 0° crit buake +e° - 6blade e° pas _ de _ gauchissement 6app blade +B° crit R.blade _-a° _ abladeB° buck =1.28 6app blade -B° Etape 6 ù Calcul d'instabilité ciénérale : 15 Cette étape fournit des données de flux admissible de flambage pour un panneau plat raidi par des poches triangulaires, dans des conditions de chargement pur ou combiné. Les formules sont basées sur le flambage de plaques orthotropes. Deux ou quatre conditions aux limites sont possibles selon le cas de chargement (4 bords en appui simple, 4 bords encastrés, 2 bords chargés en appui simple et 2 bords latéraux encastrés, 2 bords chargés encastrés et 2 bords latéraux en appui simple). Les flux appliqués, à prendre en compte pour le calcul de facteur de réserve sont les flux externes du panneau raidi par des poches triangulaires, qui sont des données d'entrée. Les données d'entrée sont les suivantes : • Données géométriques: o Lx : longueur du panneau équivalent à la grille, o Ly : largeur du panneau équivalent à la grille, o ts : épaisseur de la peau, o tg : épaisseur du panneau équivalent à la grille, • Données sur le matériau : o Exs , Eys : modules d' Young de la peau, o Gxys : module de cisaillement de la peau, o vxys, vyxs : coefficient de Poisson de la peau, o Exg , Eyg : modules d'Young de la grille, o Gxyg : module de cisaillement de la grille, o vxys, vyxs : coefficient de Poisson de la grille, 20 • charges appliquées à la structure (Nx°, Nÿ , N,y , pz) • conditions aux limites (2 ou 4 sont possibles selon le type de chargement) Les données de sortie sont : • Nxc, Nyc, Nxyc : flux de flambage admissibles, 25 • Nxccomb, NyCcomb, NxyCcomb : flux de flambage combinés admissibles, • Facteurs de réserve On utilise l'hypothèse de Kirchhoff : des sections planes demeurent principalement planes après déformation. La grille (de raidisseurs) est ici 30 modélisée par un panneau équivalent. La peau et le panneau équivalent à la grille sont considérés comme des plaques de comportement orthotrope. 15 vxy v yx Les paramètres du matériau vérifient la relation suivante : E = E1 y Les conventions de flux et moments sont illustrées par la figure 21. 6.1.1 Déplacements Le vecteur U représente le déplacement d'un point M(x,y) de la surface médiane : t if> [u, v, w] = u(x, y)z + v(x, y) ÿ + w(x, y)z Les différentes variables ne dépendent pas de z parce qu'un état de contraintes plan a été envisagé (a- = 0). 6.1.2 Déformations L'expression générale des déformations dans une section de la plaque située à une distance z de l'axe médian est : £xx = £° + z.Kxx £ =£°+z.iç yy y Equation 6-1 Avec: £o =v + w -+ 1-(wy)2 1 =~(uy+vx)+-w Equation 6-2 Les termes e° , eyy et e° représentent la contribution en déformation 20 dans le plan de la plaque. Les termes 1 (wx)2 1 (wy)2 et 1 wx.wy représentent la contribution non-linéaire en déformation dans le plan de la plaque. Le terme 15 53 R représente le rayon de la coque ("shell"), mais on considère ici une plaque plane, donc R = 0 . Equation 6-3 Les termes z.Kxx , z.Kyy et z.Kxy représentent la contribution en déformation due au changement de courbe de la plaque (z est la distance depuis l'axe médian de la plaque). 6.1.3 Lois de comportement La peau et le panneau équivalent aux raidisseurs sont considérés comme des plaques orthotropes. De ce fait les relations entre les contraintes et les déformations sont : ai yy v Ex xy y 1ùv~,vyx 1ùvxyvyx v` E` E` yx y 1-Vxyvyx 1- vxyvyx 0 0 0 0 2Gxy J Equation 6-4 avec i=(s, g) (indice s pour les valeurs relatives à la peau et indice g pour les valeurs relative à la grille de raidisseurs). 6.1.4 Flux et moments Les expressions des flux et moments par unité de longueur sont : Nad=fo-a/3dz h Maki = f aafl.zdz h Equation 6-5 avec (a,/3)=(x,y). • Flux: h h +t, Naft = J o dz = J h o-asfi dz + J `h 6c-edz h 2 2+t, En utilisant l'Equation 6-5 et la relation h = ts +tg on trouve : 0 £xx + Est + Egx tg 1ùVx~,V; 1ùVgVg/ s g g VxYEst ys + Vx7,EY tg £0 1ùVgVg 1ùVgVg yY \ xY Yx xY Yx / 1- tt • Eg Ex \1ùVgVg 1ù Vx,vs xY yx yx / gEg s Es V xY Y vs y - - vg V y' 1ùVgVg / + 2 s g Kxx + 2-t s tg YY N + / vs Est Vg Egtg v1ùvgvg 1ùVgVg/ • vg Eg v; Es - tt 2 s g 1ù VgVg 1ù VgVg g + Ey tg £o vs vs 1_ vg / vg YY û vg 1ùVx~,VY 0 + Kxx + Ests -t t 2 s g i Eg Y s y + + Nxy = 2(GsYts + vgtg ), + [tstg (Gg ù Gg )]Kxy Equation 6-6 • Moments par unité de longueur_ h h +t, Ma/3 = J6aflzdz = J h (Tala zdz+ h o zdz h 2 2+t, En utilisant l'Equation 6-5 et la relation h =ts +tg on trouve : • Eg 2 s Es ttg\1ûVgvg 1ùv vyx 1 Ï Exts(ts2+3tg2) Eg x tg(tg2+3ts2Y 12 1_ vs v;, 1ùVgVg / Ï vs Ests (ts2 + 3tg2) VgEgtg (tg2 + 3ts2 1ùVgVg 1ùVgVg 1 12 gEg s Es VxYY VxYY _tst _ 2 g1ùvgv 1ûVgvg 0 £xx + 0 £YY 1 vyy x vs Es tst 2 g\1ùvg y., yy vg 1ùvgvg/ 1 Eg ùtst y 2 glùVgvg M = yy 0 £xx + + + /v; Exts(ts2+3tg2) vgEgtg(tg2+3ts2Y ù VgVg ù vgvg / /Ests(ts2+3tg2) Egtg(tg2+3t52)A 1ûvx s yx ù M _ [tstg (Gg GXy )° o +6 (Gxts (t s2 ±32)± Ggtg (tg2 +3t52)) Equation 6-7 Dès lors, la loi de comportement générale (entre flux et moments d'une 5 part, et déformations, d'autre part), est obtenue : /All 42 0 Al2 A22 O O O A11 B11 B12 O B12 B11 O Kxy ) O O B33 / N xx C B11 B12 0 CI 1 Cu. O O B A B B12 O C~2 10 Equation 6-8 Les matrices A, B et C sont symétriques. 6.2 Equations d'équilibre Les équations d'équilibre général d'un élément du panneau (ou coque) sont données par les expressions suivantes, reliant les flux, moments et la densité de forces de surface : y+Px =0 x + p y =0 M ,+ yy,y2 ù +pz =0 R 15 Equation 6-9 où f = pxx + pY ÿ + pzz est la densité de forces de surface agissant sur l'élément de coque. La densité de force de surface agit seulement le long de la direction radiale Z, et non dans les autres directions. Donc, px = pY = 0. Par ailleurs, on considère ici le cas d'une plaque plate, donc R = 0. De ce fait, on 5 obtient des équations d'équilibre simplifiées : Equation 6-10 6.3 Pour la résolution générale de ces équations, on définit les 10 vecteurs suivants : YY,Y + -. yx,x Nwxz +NB,w, Yz +2Nx.,,wxY + pz = =0 =0 YY,Y2 + 2MXy,xY [w [M Equation 6-11 Les charges usuelles appliquées sur la plaque sont : • Flux de compression uniforme le long de l'axe x : ûN.,° 15 • Flux de compression uniforme le long de l'axe y : û N° • Flux de cisaillement uniforme dans le plan x-y : û • Pression uniforme le long de l'axe z : pzz De ce fait, comme les charges appliquées définies ci-dessus sont uniformes, il peut être déduit que les deux premières équations de l'Equation 20 6-10 sont vérifiées: x = ~ Y = NxY Y = N x =o. Expression des Moments: En utilisant les Equation 6-8 et Equation 6-11, les relations suivantes sont trouvées : [N] = A.[E]+B.[K] [M] = B.[E]+C.[K] Equation 6-12 Comme les flux appliqués sont uniformes, les relations suivantes sont obtenues : [N]x2 =0 A.[E]2 +B.[K]x2 =0 [E] x2 -1 = A .B.[W] 2 [N] Y2 = 0 A.[E],Y2 + B.[K]4Y2 = 0 [E] Y2 = A .B.[W], 2 [N]x, =0 Ale] xY+B.[K]XY =0 [d, = A .B.[W],,y Equation 6-13 Avec le fait que : [w]= On a donc pour les moments : /__ _ _ [ =B.[E]x2 +C.[K]x2 = B.A .B -C [w = -D.[W [M] = B.[E],Y2 + C.[K],Y2 = B.A .B -C .[W],Y2 = -D.[W],Y2 [M ],xY =B.[E]xy+C.[K]xy = B.A .B - C .[WL =-D.[W],xY Equation 6-14 avec D = C ù B.A .B Equation 6-15 D est la matrice de raideur globale, et est symétrique. D12 0 D22 0 0 D33/ D= Les dérivées des moments de l'Equation 6-10 sont alors obtenues : M 2=D11w x4+D12w x2y2 xx,x -M 2=D12w22+D22w4 YY,Y x Y Y MxY,Xy D33wx2Y2 Equation 6-16 L'expression de l'Equation 6-10 en termes de déplacements donne alors l'équation différentielle générale : ùN°wxz ùN°w z -2N°ywxy + pz = S21wx, +122w , +S23wxz Y2 Equation 6-17 avec : S21 =D11 n2 = D22 S23 =2.( 12+D33) Equation 6-18 Dans la partie suivante, le panneau raidi par des poches triangulaires est modélisé avec ses trois raideurs de flexion (521, S22 et S23) de façon à calculer les flux de flambage d'une plaque orthotrope. Ici encore, on utilise dans ce cas l'hypothèse suivante : si certains composants de la charge combinée sont en tension, alors ces composants ne sont pas pris en compte pour le calcul. Il est en effet conservatif de considérer que des composants en tension n'ont pas d'effet vis à vis du flux de flambage étudié, et n'améliorent pas la contrainte de flambage de la plaque. 6.4 Flux de flambage admissibles 6.4.1 Flux de compression longitudinale (compression selon X) La plaque est soumise à un flux de compression longitudinal (selon l'axe X) uniforme . De ce fait : N° = Nom, = pz = 0. L'équation différentielle générale (Equation 6-17) s'exprime alors : ûNo wxz =S21wx, +S22w +S23wx2 2 ,y , Equation 6-19 On considère dans un premier temps une plaque en appui simple (les conditions aux limites sont généralisées plus loin): Pour x=0 et x=Lx : w=0 et =0 Pour y = 0 et y = Ly : w=OetMyy=O L'expression suivante pour le déplacement w satisfait toutes les conditions aux limites détaillées plus haut : / m x ly w(x, y) = C. sin / sin n2(m, n) E N2 L L x Y Equation 6-20 L'expression précédente pour w doit satisfaire l'équation différentielle générale (Equation 6-19), donc on obtient : / 1 ml 4 +SZ2/ net +3 L x y L met x L 2 =S 4 Lx / 2/ 1 2 net Ly (m,n)E N2 N° x La valeur minimale de N° correspond à la valeur du flux admissible de flambage général Nx N. On démontre que cette valeur est : 2/ / Vn1~2 + ç3 Ly 2 Equation 6-21 Cette formule peut être généralisée pour différentes conditions aux limites (bord chargés et bords latéraux en appui simple ou encastrés) : Nx =2 Nx= avec Equation 6-22 kc =h(â)+q./3 a_ Lx 4 Ç2 Ly SZl //~~_ 3 N - 2.nIS22 h(â) q=2 La figure 22 montre la valeur de h(â) en fonction des différentes conditions aux limites (le cas des quatre bords simplement supportés est la courbe du bas). 6.4.2 Flux de compression transversale (compression selon Y) La plaque est soumise à un flux de compression uniforme transversal (selon l'axe Y):ûN. De ce fait : N° = N =p, = 0. La résolution est analogue à celle qui a été décrite dans la partie précédente. Le flux de flambage admissible Nÿ s'exprime par : 2 Ny = ./122121 Equation 6-23 avec kc = h(â) + q.,(3 ù~ +a, sitar1 sinon h(â) q=2 6.4.3 Flux de cisaillement La plaque est soumise à un flux de cisaillement uniforme :ûN . De ce fait : N° = N° = pz = O. On note que dans ce paragraphe, les formules suivantes sont valides seulement pour Ly 61 bords encastrés) (source : S. G. Lekhnitskii: Anisotropic Plates. Gordon and Breach). Les données d'entrée de la table et du graphe sont définies avec: 6.4.4 Flux de compression bi-axiale La plaque est soumise à un chargement combiné : un flux uniforme de compression longitudinale (selon l'axe x) et un flux uniforme de compression transversal (selon l'axe y) : û N°omb et û Ny omb . De ce fait : Nom, = pz = 0. On 10 définit a par : Ny omb = 2,N° comb . L'équation différentielle générale s'exprime -N°combw,x2 -NYcombw,Y2 _ ~iwx4 +Ç w 4 +n3W Y2 Equation 6-25 • Conditions aux limites : quatre bords en appui simple : 15 Pour x=0 et x=Lx : w=0 et ùM =0 Pour y=0ety=LY: w=OetMB,=O Expression du déplacement: L'expression suivante du déplacement w satisfait toutes les conditions aux limites détaillées plus haut : 20 w(x, y) = C. sin m sin n (m, n)E N2 L L x y Equation 6-26 Flux de flambage admissibles (NC omb , Nycomb )_ L'expression précédente pour w doit satisfaire l'équation différentielle générale (Equation 6-25), et on obtient donc : / 2 m7L o L + Nycomb x / 2 nTC ù mTc L L 1\ Y x / 4 / 2/ 2 net + mit net \ Ly s Lx Ly 4 N° comb (m,n)E N2 L'expression de N° omb fonction de a doit satisfaire : 2 , 2 SZI/Ly m4 +522 Lx n4 +SZ3m2n2 Lx \LY 0 Nxcomb ù Ly2m2 + x /LL2n2 7c2 N2 ainsi, les flux de flambage admissibles sont obtenus : 7c2 L2m2 +2L 2n2 y 2 , 2 /L Y m4 +522 L Lx y L c Nx comb =Min SZ1 n4+SZ Ne Y comb Equation 6-27 Conditions aux limites : quatre bords encastrés : Pour x=0 et x=Lx : w=0 et aw=0 y Pour y=0 et y=Ly : w=0 et az =0 Expression du déplacement: L'expression suivante du déplacement w satisfait toutes les conditions 15 aux limites détaillées plus haut : x comb W(x,y)=Cmn 1 ù cos /n2TCy Ly N2 / m•2.z.x ù cos Lx Equation 6-28 Flux de flambage admissible (Nxcomb , Nxcomb )= L'expression de N° comb fonction de a doit satisfaire : 4,t2 N0 xcomb ù L 2m2 + x2n2 2 , 2 SZ~/~y m4+522 Lx n43S23m2n2 x y 20 (m,n)E N2 ainsi, les flux de flambage admissibles sont obtenus : 4,t-2 2 + /ÎLx 2n 2 z /L _ y 4 +~2 Lx n4 + S23m2n2 Lx ~Ly 3 Nxc comb =Min SZ1 (m,n)E N2 Ne = Y comb Equation 6-29 6.4.5 Flux de compression longitudinale et de cisaillement La plaque est soumise à des chargement combiné : flux de compression longitudinale uniforme (selon l'axe X) et de cisaillement : ù N°COmb et ùN 0 xY comb De ce fait : N° = pz = 0 . Equation d'interaction : L'équation d'interaction pour des flux combinés de compression longitudinale et de cisaillement est : R x + R 1.75 =1 xv Equation 6-30 Avec N Rx = N°mb rapport de flux de compression longitudinal x rapport de flux de cisaillement où Nxc et Nxye sont les flux de flambage admissibles calculés ci-dessus pour un chargement uni-axial. 20 6.4.6 Flux de compression transversale et de cisaillement La plaque est soumise à un flux de compression transversale uniforme (selon l'axe y) et à un flux de cisaillement : ù Nyet ùN° comb xYcomb De ce fait: N° = pz = 0 . L'équation d'interaction pour les flux combinés de compression 25 transversale et de cisaillement est : C x comb15 R +R 1.75 =1 Y xY Equation 6-31 Avec : NC Y comb Nÿ transverse rapport de flux de compression rapport de flux de cisaillement où Nÿ et Ne sont les flux de flambage admissibles calculés plus haut pour un chargement uni-axial. 6.4.7 Flux de compression Bi-axiale et de cisaillement La plaque est soumise à des chargement combiné : un flux de compression longitudinal uniforme (selon l'axe x) et un flux de compression transversal uniforme (selon l'axe y) en plus d'un flux de cisaillement : ù N°comb ùN° et ù Y comb De ce fait: pz = 0 . L'équation d'interaction est obtenue en deux étapes. D'abord, on détermine un facteur de réserve RFb; correspondant au flux de compression biaxial : RFbi = No ù No x comb Y comb Equation 6-32 Puis cette valeur est utilisée dans une équation d'interaction pour des flux combinés de compression bi-axiale et de cisaillement : Rb + RXY1.75 =1 Equation 6-33 25 Avec : Rbi = b 1 rapport de flux de compression bi-axial o comb Ne c Nx comb NY comb20 N° xY comb c rapport de flux de cisaillement où N est le flux de flambage admissible calculé en cisaillement pur. Le procédé selon l'invention comporte également une boucle d'itération (voir figure 7). Cette boucle permet de modifier les valeurs de charges appliquées, ou les valeurs dimensionnelles des panneaux raidis par des poches triangulaires considérés, selon les résultats de l'une au moins des étapes 3 à 6. La procédé, tel qu'il a été décrit, peut être mis en oeuvre au moins partiellement sous forme d'une macro sur un logiciel de type tableur. Celle-ci utilise alors, par exemple, comme entrées des données de matériau et de géométrie stockées sur une zone dédiée, ainsi que divers cas de charges considérés, et de conditions aux limites, et fournit en sortie des valeurs de masse du panneau, de facteur de réserve à charge ultime concernant notamment les poches triangulaires, les raidisseurs, l'effondrement général. Ces données de sortie mettent ainsi en évidence les cas de charges ou de dimensionnement incompatibles avec les facteurs de réserve souhaités. Avantages de l'invention On comprend que le procédé NASA antérieurement connu, a été substantiellement étendu dans le cadre de la présente invention pour prendre en compte les spécificités du domaine aéronautique : - Valeurs locales admissibles pour les raidisseurs (destruction, instabilité latérale etc.) pour la compression selon la direction X ou Y et le cisaillement, - Valeurs locales admissibles de la peau triangulaire pour la compression selon la direction X ou Y et le cisaillement, - Correction de plasticité, - Calcul de masse préliminaire, - Calcul de flambage général pour une compression selon la direction X ou Y et un cisaillement. Les améliorations principales sont le calcul des contraintes admissibles pour les différents types de flambage et les calculs des facteurs de réserve adaptés. Extensions: Cas de chargement - double-compression (pour le flambage local et global) - Chargement combiné : compression et cisaillement Extensions: améliorations des paramètres de la méthode - Pas de limitation sur le coefficient de Poisson du matériau - Variation de l'angle de la grille (différent de 60°) - Plasticité lorsque la structure raidie par des poches triangulaires est considérée comme un panneau équivalent raidi - Conditions aux limites (attache û"clamping"- ou conditions aux limites intermédiaires) sur le flambage local ou global 15 Un des avantages les plus significatifs du procédé de dimensionnement selon l'invention est la possibilité d'installer des panneaux raidis par des poches triangulaires, en lieu et place de panneaux précédemment réalisés avec deux familles perpendiculaires de raidisseurs 20 sous ( longerons et traverses), résultant, à résistance mécanique égale, en un gain de masse atteignant 30% sur certaines pièces. Variantes de l'invention La portée de la présente invention ne se limite pas aux détails des 25 formes de réalisation ci-dessus considérées à titre d'exemple, mais s'étend au contraire aux modifications à la portée de l'homme de l'art. On a cité dans la présente description des angles de base de triangles isocèles compris entre 45° et 70°, qui correspondent au besoin courant dans les structures aéronautiques. Il est clair cependant qu'un procédé analogue 30 peut être mis en oeuvre pour toute valeur d'angle isocèle dans des panneaux raidis par des poches triangulaires. 10 REVENDICATIONS1. Procédé de dimensionnement par méthode analytique d'un panneau sensiblement plan en matériau homogène et isotrope, le panneau étant composé d'une peau renforcée par un ensemble (dit "grille") de trois faisceaux parallèles de raidisseurs intégrés au panneau, les poches déterminées sur la peau par ces familles de raidisseurs étant triangulaires, les raidisseurs étant en forme de lame, le panneau raidi devant satisfaire un cahier des charges de résistance mécanique à des charges externes prédéterminées, caractérisé en ce que les angles entre faisceaux de raidisseurs sont tels que les poches triangulaires sont de forme isocèle quelconque. 10 2. Procédé selon la revendication 1, caractérisé en ce que le procédé prend en compte la redistribution des contraintes appliquées entre le panneau et la grille de raidisseurs due : au post-flambage des raidisseurs, par la définition d'une section droite 15 efficace pour chaque type de raidisseur (0°, +0 or -0): Aoô t , A+est et A_esf, au post-flambage de la poche à travers le calcul d'une épaisseur efficace du panneau : ts eff, à la plasticité des charges externes appliquées, grâce à un procédé itératif sur les différentes propriétés du matériau, notamment modules d'Young 20 et coefficients de Poisson : Eoô t, E+est, E_est pour les raidisseurs et EXS, Eys et vepst pour la peau, et utilisant la loi de Ramberg û Osgood. 3. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte une étape de correction des charges 25 appliquées pour tenir compte de la plasticité utilisant un procédé de calcul itératif de contraintes plastiques, réalisé jusqu'à ce que les cinq paramètres du matériau (Eo.st ,E+est, E-est, Epeau, vep) entrés au début du processus soient sensiblement égaux aux mêmes paramètres obtenus après calcul de contrainte plastique. 30 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte une étape 5A de calcul de flux de flambage admissible, et de facteur de réserve pour des poches triangulaires isocèles, les contraintes appliquées à prendre en compte pour le calcul de facteur de réserve étant les contraintes agissant dans la peau uniquement, les flux externes utilisés étant des flux de peau ne correspondant pas au chargement complet du panneau raidi. 5. Procédé selon la revendication 4, caractérisé en ce qu'il comporte deux sous-étapes : d'abord de calcul de valeurs admissibles pour des plaques soumises à des cas de chargement pur (compression selon deux directions dans le plan, cisaillement) en utilisant une méthode par éléments finis, puis de calcul des courbes d'interaction entre ces cas de chargement pur. 6. Procédé selon la revendication 5, caractérisé en ce que le calcul de valeurs admissibles comporte des sous-étapes suivantes : • Création d'un modèle paramétrique FEM d'une plaque triangulaire • Exécution de nombreuses combinaisons différentes pour obtenir des résultats de flambage, • Obtention de paramètres compatible avec une formulation analytique polynomiale. 7. Procédé selon la revendication 6, caractérisé en ce que, dans le cas d'un chargement pur, les courbes d'interaction sont définies par des sous- étapes suivantes : - de création de modèles éléments finis d'une pluralité de plaques triangulaires présentant des angles isocèles différents, l'angle isocèle (0) étant défini comme l'angle de base du triangle isocèle, - pour chaque angle isocèle :1/ de calcul par Modèle Eléments Finis pour déterminer le flux admissible de plissement (sans correction plastique) pour diverses épaisseurs de plaque. 2/ de tracé d'une courbe de flux de flambage admissible en fonction du rapport D (D raideur de plaque, h hauteur du triangle), cette courbe étant déterminée pour des petites valeurs du rapport D , par une équation au second degré fonction de ce rapport, dont les coefficients K1 et K2 dépendent de l'angle et du cas de charge considéré, 3/ de tracé de l'évolution des coefficients de l'équation polynomiale K1 et K2 en fonction de l'angle de base du triangle isocèle, ces coefficients étant tracés en fonction de l'angle des plaques triangulaires considérées, et d'interpolation pour déterminer une équation polynomiale qui permet de calculer ces constantes quel que soit l'angle isocèle. 8. Procédé selon l'une quelconque des revendications 6 à 7, caractérisé en ce que, dans le cas d'un chargement combiné, on utilise l'hypothèse suivante : si certains composants de la charge combinée sont en tension, alors ces composants ne sont pas pris en compte pour le calcul, et en ce que les courbes d'interaction sont définies par des sous-étapes suivantes : - de création de modèles éléments finis d'une pluralité de plaques triangulaires présentant des angles isocèles différents, l'angle isocèle (0) étant défini comme l'angle de base du triangle isocèle, - pour chaque angle, 1/ de calcul par Modèle Eléments Finis (FEM) pour déterminer la valeur propre de flambage correspondant à différentes répartitions de charges externes. 2/ de tracé des courbes d'interaction, pour chaque angle et chaque combinaison de charges, puis d'approximation de ces courbes avec une unique équation couvrant toutes les combinaisons :R~XA B c où R~ = N aPP +Rey +Rs = 1 (où , i = cX, cY ou s), A, B, C étant des coefficients empiriques. 9. Procédé selon l'une quelconque des revendications 6 à 8, caractérisé en ce que, dans le cas de plaques triangulaires isocèles en appui 5 simple ou encastrées, en cas de chargement combiné, on utilise une courbe 3 d'interaction : R~X +R~Y +Rsz =1, pour tous les cas de chargement. 10. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte une étape 6, de calcul d'instabilité générale, 10 fournissant des données de flux admissible de flambage, et des facteurs de réserve, pour un panneau plat raidi, dans des conditions de chargement pur ou combiné, les flux appliqués, à prendre en compte pour le calcul de facteur de réserve, étant les flux externes du panneau raidi. 15 11. Produit-programme d'ordinateur comportant une série d'instructions adaptées à mettre en oeuvre un procédé selon l'une quelconque des revendications précédentes, lorsque cette suite d'instructions est exécutée sur un ordinateur.