diff --git a/abaqueCuve.py b/abaqueCuve.py index 3cb92d8..2508aff 100644 --- a/abaqueCuve.py +++ b/abaqueCuve.py @@ -2,7 +2,7 @@ """ Created on Mon Feb 10 16:29:16 2025 -@author: Maël Dbq +@authors: Maël Dbq & Achille Toupin """ import numpy as np @@ -11,6 +11,15 @@ import matplotlib.dates as mdates from scipy.optimize import root_scalar import datetime +# Données + +rayon_cuve = 0.6 # Rayon de la cuve en m +longueur_cuve = 1.77 # Longueur de la cuve en m +nb_grad_fleche = 25 # Nombre de graduations sur l'axe des hauteurs +nb_grad_volume = 21 # Nombre de graduations sur l'axe des volumes +date_debut = datetime.date.today() +duree_graphe = 24 # Durée du graphe en mois + # Fonction pour calculer l'aire du segment de disque pour une flèche donnée def segment_area(f, R): if f == 0: @@ -29,46 +38,63 @@ def find_f(A, R): except ValueError: return None # Retourne None si aucun zéro trouvé -# Demande du rayon et de la longueur avec conversion de la virgule -rayon_cuve = float(input("Entrez le rayon de la cuve (m) : ").replace(",", ".")) -longueur_cuve = float(input("Entrez la longueur de la cuve (m) : ").replace(",", ".")) +# Création d'une échelle temporelle +dates = [date_debut.replace(day=1) + datetime.timedelta(days=30*i) for i in range(duree_graphe)] # Capacité maximale de la cuve (cylindre entier) volume_max = np.pi * rayon_cuve**2 * longueur_cuve -# Liste des volumes espacés de 0.1 m³ -volumes = np.arange(0, volume_max + 0.1, 0.1) +# Liste des volumes pour chaque date +volumes = np.linspace(0, volume_max, len(dates)) -# Calcul des hauteurs de flèche associées -flèches = [find_f(v / longueur_cuve, rayon_cuve) for v in volumes] +# Création de la liste de flèches +flèches = np.linspace(0, 2 * rayon_cuve, len(dates)) + +# Création des graduations de l'axe de droite +right_ticks = np.linspace(0, 2 * rayon_cuve, nb_grad_fleche) +right_tick_positions = [find_f(v / longueur_cuve, rayon_cuve) for v in right_ticks * volume_max / (2 * rayon_cuve)] # Supprimer les valeurs None (au cas où) volumes = np.array([v for i, v in enumerate(volumes) if flèches[i] is not None]) flèches = np.array([f for f in flèches if f is not None]) -# Création d'une échelle temporelle : une valeur par mois sur un an -dates = [datetime.date.today() + datetime.timedelta(days=30*i) for i in range(len(volumes))] - # Espacement en flèche pour avoir des graduations équidistantes en volume -nb_grad = 10 # Nombre de graduations voulues -volumes_grad = np.linspace(0, volume_max, nb_grad) +volumes_grad = np.linspace(0, volume_max, nb_grad_fleche) flèche_grad = [find_f(v / longueur_cuve, rayon_cuve) for v in volumes_grad] +plt.rcParams.update({'font.size': 22}) # Taille du texte + # Tracer le graphe -fig, ax1 = plt.subplots(figsize=(10, 6)) +fig, ax1 = plt.subplots(figsize=(29.7-1, 21-1)) # Format A4 avec une marge de 1 cm # Axe principal : Volume en fonction des mois -ax1.plot(dates, volumes, label="Volume (m³)", color='b', linewidth=2) -ax1.set_xlabel("Temps (mois)") -ax1.set_ylabel("Volume (m³)", color='b') + +volumes_L = volumes * 1000 # Conversion en litres +ax1.plot(dates, volumes_L, label="Volume (L)", color='b', linewidth=2, alpha=0) +ax1.set_ylabel("Volume (L)", color='b') ax1.tick_params(axis='y', labelcolor='b') -ax1.xaxis.set_major_formatter(mdates.DateFormatter("%b %Y")) # Format mois année + +ax1.yaxis.set_major_locator(plt.MaxNLocator(nb_grad_volume)) # Configuration du nombre de graduations +ax1.set_ylim(bottom=0, top=volume_max * 1000) # Limite inférieure à 0 et supérieure à la capacité maximale + +# Dates en abscisse +ax1.xaxis.set_major_formatter(mdates.DateFormatter("%m/%Y")) # Format mois (numérique) année ax1.xaxis.set_major_locator(mdates.MonthLocator()) # Une graduation par mois plt.xticks(rotation=45) # Deuxième axe pour les hauteurs de flèche ax2 = ax1.twinx() -ax2.plot(dates, flèches, label="Hauteur de remplissage (m)", color='r', linestyle='dashed', linewidth=2) + +# Application des ticks et labels personnalisés +ax2.set_yticks(right_tick_positions) # Positions réelles des ticks + +right_ticks_cm = [f * 100 for f in right_ticks] # Valeurs en cm +right_ticks_text = [f"{f:.0f}" for f in right_ticks_cm] # Texte des labels +ax2.set_yticklabels(right_ticks_text) # Labels affichés + +ax2.set_ylim(bottom=0, top=2 * rayon_cuve) # Limite inférieure à 0 et supérieure à la hauteur maximale + +ax2.plot(dates, flèches, label="Hauteur de remplissage (cm)", color='r', linestyle='dashed', linewidth=2) ax2.set_ylabel("Hauteur de remplissage (m)", color='r') ax2.tick_params(axis='y', labelcolor='r') @@ -76,7 +102,13 @@ ax2.tick_params(axis='y', labelcolor='r') ax2.set_yticks(flèche_grad) # Aligner les graduations de l'axe de droite sur les hauteurs ax2.grid(visible=True, linestyle='--', alpha=0.6, which='both') # Quadrillage en pointillé -# Légende et titre -fig.suptitle("Évolution du remplissage de la cuve au fil du temps") +# Ajouter un quadrillage qui colle aux volumes +ax1.set_yticks(volumes_grad) # Aligner les graduations de l'axe de gauche sur les volumes +ax1.grid(visible=True, linestyle='--', alpha=0.6, which='both') # Quadrillage en pointillé + +# Show the plot fig.tight_layout() -plt.show() \ No newline at end of file +plt.show() + +# Enregistrement de la figure +fig.savefig("graphes/abaque.png") \ No newline at end of file