====== Graphiques des pressions partielles de systèmes non-idéaux ======

<sxh  python; title : pressions_partielles_systemes_non_ideaux.py>
#!/usr/bin/env python
# -*- coding: utf-8 -*-
"""
Graphiques des pressions partielles de systèmes non-idéaux
Basé sur le travail de ML et VM, ba2 chimie 2013-2014
"""

import numpy as np
import matplotlib.pyplot as plt
from Tkinter import*

#données
R = 0.082                  #cste des gaz parfaits
T = 293                    # temperature cste (k)
U01 = np.array([343,330,328,313,308,295,290,284,275,274,263,262,254,252,250,246,242,241,233,232,227,225,217,207,180,124,121,110,103,86,73,62,0]) #[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]  #potientiel chimique standard, donnee = exemple
U02 = np.array([0,110,122,191,206,259,272,284,323,328,357,361,381,382,390,394,403,405,419,420,427,428,438,447,465,490,489,491,491,495,502,501,513])
a1 = np.array([0,0.061,0.071,0.121,0.133,0.185,0.198,0.209,0.275,0.288,0.35,0.356,0.405,0.414,0.447,0.452,0.492,0.498,0.569,0.572,0.612,0.615,0.67,0.72,0.828,0.918,0.923,0.934,0.939,0.954,0.961,0.969,1])#f represente le rapport des constituants, o etant le constituant 1 et p le consituant 2 par ex
a2 = np.ones(1)-a1

x1 = np.array([0,0.1891,0.3851,0.5945,0.8243,1])
x2 = np.ones(1)-x1
U03 = np.array([80])
U04 = np.array([45])

def graph1():
    U1= (np.log(a1)*R*T)+U01
    U2= (np.log(a2)*R*T)+U02
    plt.figure(1)
    plt.plot(a1,U1,'r''x''-')
    plt.plot(a1,U2,'g''o''-')
    plt.plot(a1, U1+U2,'b')
    plt.xlabel ('fraction molaire des constituants')
    plt.ylabel ('pression en mmHg')
    plt.title ('diagramme de l acetone et du CS2')
    plt.legend(("acetone","CS2","somme des 2"),'best')
    plt.show()

def graph2():
    U3= (np.log(x1)*R*T)+U03
    U4= (np.log(x2)*R*T)+U04
    plt.figure(2)
    plt.plot(x2,U3,'b''-')
    plt.plot(x2,U4,'g''x''-')
    plt.plot(x2, U3+U4,'k''o''-')
    plt.xlabel ('fraction molaire des constituants')
    plt.ylabel ('pression en mmHg')
    plt.title ('diagramme du benzene et du toluene')
    plt.legend(("benzene","toluene","somme des 2"),'best')
    plt.show()

def graph3 ():
    U1= (np.log(a1)*R*T)+U01
    U2= (np.log(a2)*R*T)+U02
    U3= (np.log(x1)*R*T)+U03
    U4= (np.log(x2)*R*T)+U04
    
    plt.figure(3)

    plt.subplot(121)
    plt.plot(a1,U1,'r''x''-')
    plt.plot(a1,U2,'g''o''-')
    plt.plot(a1, U1+U2,'b')
    plt.xlabel ('fraction molaire des constituants')
    plt.ylabel ('pression en mmHg')
    plt.title ('diagramme de l acetone et du CS2')
    plt.legend(("acetone","CS2","somme des 2"),'best')

    plt.subplot(122)
    plt.plot(x2,U3,'b''-')
    plt.plot(x2,U4,'y''x''-')
    plt.plot(x2, U3+U4,'k''o''-')
    plt.xlabel ('fraction molaire des constituants')
    plt.ylabel ('pression en mmHg')
    plt.title ('diagramme du benzene et du toluene')
    plt.legend(("benzene","toluene","somme des 2"),'best')

    plt.show()

def abientot ():
    print "en esperant vous revoir bientôt !"
    fen1.destroy()

fen1 = Tk()
fen1.title("diagramme des cas non ideaux et ideaux")

#création widget

Label(fen1, text="choisissez votre diagramme parmis les propositions ci dessous :", bg="orange", fg="black").grid(row=2, column=1,columnspan=2, padx=3, pady=3)
Button(fen1, text="acetone CS2", command=graph1, bg="yellow", fg="black").grid(row=3, column=0, columnspan=2)
Button(fen1, text="benzene toluene", command=graph2, bg="blue", fg="white").grid( row=3, column=2, columnspan=2)
Button(fen1, text="comparatif des 2 diagrammes", command=graph3, bg="green",fg="black").grid( row=3,column=4,columnspan=2)
Button(fen1, text="quitter", command=abientot, bg="red", fg="black").grid(row=5, column=2, columnspan=2)

fen1.mainloop()

</sxh>