Partage
  • Partager sur Facebook
  • Partager sur Twitter

Listes dans une boucle for

    24 mai 2024 à 14:46:37

    Bonjour, 

    Je viens vers vous pour récolter de l'aide.

    Mon objectif est de mettre des couleurs différentes pour chaque courbe qu'il y a sur mon graphique.

    Sur mon code je trace mon graphique dans ma boucle for, cependant je n'arrive pas à mettre une couleur différente pour chaque courbe ce qui est normal. Voir même pouvoir créer une légende pour identifier les courbes.

    Ma première idée était de créer des listes pour récolter les données qui m'intéressent pour ainsi tracer les courbes avec les couleurs que je veux. Or il faudrait que je créer une nouvelle liste à chaque fois que T1 prend une nouvelle valeur pour sauvegarder les données de ntag.

    Mais je ne sais comment faire.

    Cordialement, 

    import CoolProp.CoolProp as CP
    import matplotlib.pyplot as plt
    
    P2 = [2*1e5, 3*1e5, 4*1e5, 5*1e5, 6*1e5, 7*1e5, 8*1e5, 9*1e5, 10*1e5, 11*1e5, 12*1e5]
    
    for T1 in [357.7, 306.6, 268.3, 238.5] :
        for P2 in [2*1e5, 3*1e5, 4*1e5, 5*1e5, 6*1e5, 7*1e5, 8*1e5, 9*1e5, 10*1e5, 11*1e5, 12*1e5] :
            # Cycle Turbine à gaz
    
            # Au point 1 : Entrée Compresseur
            P1 = 1*1e5 # (Pascal)
            h1 = CP.PropsSI('H','P',P1,'T',T1,'air')
            s1 = CP.PropsSI('S','P',P1,'T',T1,'air')
    
            # Au point 2 isentropique
            s2is = s1 # (kJ/kg.K)
            T2is = CP.PropsSI('T','P',P2,'S',s2is,'air')
            h2is = CP.PropsSI('H','P',P2,'S',s2is,'air')
    
            # Au point 2 : Sortie Compresseur / Entrée Réchauffeur
            n2 = 0.85 # Rendement isentropique compresseur
            h2 = (h2is-h1)/n2 + h1
            s2 = CP.PropsSI('S','P',P2,'H',h2,'air')
            T2 = CP.PropsSI('T','P',P2,'H',h2,'air')
    
            # Au point 3 : Sortie Réchauffeur / Entrée Turbine à gaz 
            T3 = 800 + 273.15
            P3 = P2
            h3 = CP.PropsSI('H','P',P3,'T',T3,'air')
            s3 = CP.PropsSI('S','P',P3,'T',T3,'air')
    
            # Au point 4 isentropique
            P4 = P1
            s4is = s3
            T4is = CP.PropsSI('T','P',P4,'S',s4is,'air')
            h4is = CP.PropsSI('H','P',P4,'S',s4is,'air')
    
            # Au point 4 : Sortie Turbine à gaz
            n4 = 0.85 # Rendement isentropique turbine
            h4 = n4*(h4is-h3)+h3
            s4 = CP.PropsSI('S','P',P4,'H',h4,'air')
            T4 = CP.PropsSI('T','P',P4,'H',h4,'air')
            
            Pre = 10*1e6       
            mair = Pre/(h3 - h2)
            Ptag = mair*(h4 - h3 + h2 - h1)
            ntag = (abs(Ptag))/Pre 
            plt.plot(P2, ntag, 'b.-' )
    
     

    -
    Edité par FlorianFerri 24 mai 2024 à 14:49:56

    • Partager sur Facebook
    • Partager sur Twitter
      25 mai 2024 à 10:47:38

      Pour tracer des courbes avec des couleurs différentes et ajouter une légende sur votre graphique, il est possible d'utiliser un tableau de couleurs prédéfini et une liste pour stocker les valeurs de ntag pour chaque valeur de T1.

      • Partager sur Facebook
      • Partager sur Twitter

      Celui qui trouve sans chercher est celui qui a longtemps cherché sans trouver.(Bachelard)
      La connaissance s'acquiert par l'expérience, tout le reste n'est que de l'information.(Einstein)

      Listes dans une boucle for

      × Après avoir cliqué sur "Répondre" vous serez invité à vous connecter pour que votre message soit publié.
      • Editeur
      • Markdown