import numpy as np

import matplotlib.pyplot as plt

# Constantes

rho = 1000 # kg/m^3, densité de l'eau

U = 2 # m/s, vitesse du courant

S = 42 # m², surface de l'aile (2b x c)

# Angles d’attaque en degrés

alpha_deg = np.linspace(-15, 15, 300)

alpha_rad = np.radians(alpha_deg) # Conversion en radians

# Coefficient de portance

Cl = 2 * np.pi * np.sin(alpha_rad)

# Force de portance (en Newtons)

L = 0.5 * rho * S * Cl * U**2

# Tracé

plt.figure(figsize=(8,5))

plt.plot(alpha_deg, L, label='Force de portance L(α)', color='blue')

plt.xlabel("Angle d’attaque α (°)")

plt.ylabel("Force de portance L (N)")

plt.title("Portance d'une aile NACA0015 en fonction de l’angle d’attaque")

plt.grid(True)

plt.axhline(0, color='black', lw=0.5)

plt.legend()

plt.tight_layout()

plt.show()