Pyzo

Pour ne pas se faire piéger par différentes unités utilisées dans le code, on prendra l'habitude de les indiquer dans le nom de la variable (exemple au lieu de prendre comme nom de variable t, alpha, Wmot... on indiquera t_s, Alpha_deg, Wmot_radpars

Exemple :

import numpy as np

import matplotlib.pyplot as plt

## 1. Récupération des valeurs EXPérimentales

xEXP_mm,alphaEXP_deg,bêtaEXP_deg=np.loadtxt("P:\\Documents\\TP SII\\valeursEXP.csv", skiprows=7, unpack=True, usecols=(1,3,8))

##

## 2. Récupération des valeurs issues des simulations NUMériques

xNUM_mm,alphaNUM_deg,bêtaNUM_deg=np.loadtxt("P:\\Documents\\TP SII\\valeursNUM.csv", skiprows=2, unpack=True, usecols=(0,1,2))

##

## 3. Détermination des valeurs issues des simulations MANuelles

# Paramètres caractéristiques

L1_mm=50

L2_mm=35

L3_mm=90

#NB : toutes les expressions des variables doivent être précisées littéralement

xMAN_mm=np.arange(20,55,1)

alphaMAN_rad=-np.cos((L3_mm**2-xMAN_mm**2)/(L1_mm*L2_mm))

bêtaMAN_rad=-alphaMAN_rad+np.sin((L2_mm**2-xMAN_mm**2)/(L3_mm*xMAN_mm))

alphaMAN_deg=360/2/np.pi*alphaMAN_rad

bêtaMAN_deg=360/2/np.pi*bêtaMAN_rad

##

## 4. Tracés

plt.subplot(211)

plt.plot(xEXP_mm,alphaEXP_deg,label="expérimental")

plt.plot(xNUM_mm,alphaNUM_deg,label="simul. numérique")

plt.plot(xMAN_mm,alphaMAN_deg,label="simul. manuelle")

plt.legend(fontsize=15)

plt.xlabel("x en mm")

plt.ylabel("alpha en degré")

plt.title("Alpha en fonction de x")

plt.subplot(212)

plt.plot(xEXP_mm,bêtaEXP_deg,label="expérimental")

plt.plot(xNUM_mm,bêtaNUM_deg,label="simul. numérique")

plt.plot(xMAN_mm,bêtaMAN_deg,label="simul. manuelle")

plt.legend(fontsize=15)

plt.xlabel("x en mm")

plt.ylabel("bêta en degré")

plt.title("Bêta en fonction de x")

plt.show()

##