Courbe de Prédominance d'un Acide
<sxh python; title : courbe_predominance_acide.py> #!/usr/bin/env python # -*- coding: utf-8 -*- # travail de KH, ba2 chimie 2012-2013
# Courbe de Prédominance d'un Acide # from math import * import matplotlib.pyplot as plt from Tkinter import *
def equations (n,C,a,b,c,d,e,f,g,pKa1,pKa2,pKa3): #initialiser (valeur 0)
pH = 0. #pH de départ
if n==3:
Ka1=10**(-pKa1)
Ka2=10**(-pKa2)
Ka3=10**(-pKa3)
while pH < 14: #pH va jusque 14
H = 10**(-pH) #équation pour avoir les concentrations en H30+ aux différents pH demandé
O = (10**(-14)) / H #équation pour avoir les concentrations en OH- aux différents pH demandé
K = 1 + H/Ka3 + (H**2)/(Ka2*Ka3) + (H**3)/(Ka1*Ka2*Ka3) #équation qui détermine la valeur de alpha H
a0= 1/K # équation pour avoir alpha0
a1=a0*(H/Ka3) # équation pour avoir alpha1
a2=a1*(H/Ka2) # équation pour avoir alpha2
a3=a2*(H/Ka1) # équation pour avoir alpha3
H3Y = a3 * C # équation pour avoir la [H3Y]
H2Y = a2 * C # équation pour avoir la [H2Y-]
HY2 = a1 * C # équation pour avoir la [HY2-]
Y3 = a0 * C # équation pour avoir la [Y3-]
# listes auxquels nous ajouterons les différentes valeurs calculées au dessus pour chaque pH demandé :
a.append(pH)
b.append (H)
c.append (O)
d.append (H3Y)
e.append (H2Y)
f.append (H2Y)
g.append (Y3)
pH = pH + 0.1
elif n==2:
Ka1=10**(-pKa1)
Ka2=10**(-pKa2)
while pH < 14:
H = 10**(-pH) #équation pour avoir les concentration en H3O+ aux différents pH demandé
O = (10**(-14)) / H # équation pour avoir les concentration en OH - aux différents pH demandé
K = 1 + H/Ka2 + (H**2)/(Ka2*Ka1) #équation pour avoir la valeur de alpha H
a0 = 1/K #équation pour avoir alpha 0
a1= a0*(H/Ka1) #équation pour avoir alpha 1
a2 = a1*(H/Ka1) #équation pour avoir alpha 2
H2Y = a2*C #équation pour avoir la [H2Y]
HY=a1*C #équation pour avoir la [HY-]
Y2=a0*C # équation pour avoir la [Y2-]
# listes auxquels nous ajouterons les différentes valeurs calculées au dessus pour chaque pH demandé :
a.append(pH)
b.append(H)
c.append(O)
d.append(H2Y)
e.append(HY)
f.append(Y2)
pH = pH + 0.1
elif n==1:
Ka1=10**(-pKa1)
while pH < 14 :
H = 10**(-pH) #équation pour avoir les concentration en H3O+ aux différents pH demandé
O = (10**(-14)) / H # équation pour avoir les concentration en OH - aux différents pH demandé
K = 1 + H/Ka1 #équation pour avoir la valeur de alpha H
a0 = 1/K #équation pour avoir alpha 0
a1= a0*(H/Ka1) #équation pour avoir alpha 1
HY= a1*C #équation pour avoir la [HY]
Y=a0*C #équation pour avoir la [Y-]
# listes auxquels nous ajouterons les différentes valeurs calculées au dessus pour chaque pH demandé :
a.append(pH)
b.append (H)
c.append (O)
d.append (HY)
e.append (Y)
pH = pH + 0.1
#2eme fenêtre dans laquelle on rentrera les paramètres propre à la molécule étudiée. def fenetre2():
root2 = Tk() n = float (champn.get ()) # transformation en nombre flottant de la valeur entrée pour n
txtn = Label(root2, text ='pKa1 : ', fg = 'red')
txtn.grid(row=2)
champ1 = Entry(root2) #demander le pKa1
champ1.grid(row=3) # ligne à laquelle apparaitra le champs ou l'on donnera le pKa1
if n > 1 : #paramètre demander si la molécule possède 2 H+
txt2 = Label(root2, text ='pKa2 : ', fg = 'red')
txt2.grid(row=4)
champ2 = Entry (root2) #demande le pKa2
champ2.grid(row=5) # ligne à laquelle apparaitra le champs ou l'on donnera le pKa2
if n == 3 : #paramètre demander si la molécule possède 3 H+
txt3 = Label(root2, text ='pKa3 : ', fg = 'red')
txt3.grid(row=6)
champ3 = Entry (root2) #demande le pKa3
champ3.grid(row=7) #ligne à laquelle apparaitra le champs dans lequel on donnera le pKa3
def graphe():
a=[]
b=[]
c=[]
d=[]
e=[]
f=[]
g=[]
n = float(champn.get())
C = float(champc.get())
pKa1 = float(champ1.get())
pKa2 = 0
pKa3 = 0
if n > 1:
pKa2 = float (champ2.get())
if n == 3:
pKa3 = float (champ3.get())
equations(n,C,a,b,c,d,e,f,g,pKa1 = pKa1, pKa2 = pKa2, pKa3 = pKa3)
plt.xlabel("pH") #nom donné à l'axe des abscisses
plt.ylabel("Concentrations") #nom donné à l'axe des ordonnées
# les différentes droites qui seront présentes sur le graphique semilog :
c1 = plt.semilogy(a,b,'r') #courbe donnant le concentration en H3O+ en fonction du pH
c2 = plt.semilogy(a,c, 'b') #courbe donnant le concentration en OH- en fonction du pH
c3 = plt.semilogy(a,e,'g') #courbe donnant le concentration en l'acide qui à tout ces H+ en fonction du pH
c4 = plt.semilogy(a,d,'m')
if len(f)>1:
c5 = plt.semilogy(a,f,'c') #courbe donnant le concentration en l'acide qui a perdu 2H+ en fonction du pH
if len(g)>1:
c6 = plt.semilogy(a,g,'y') #courbe donnant le concentration en l'acide qui a perdu 3H+ en fonction du pH
plt.title("Courbe representant la predominance")
plt.show()
b2=Button(root2,text="Tracer le Graphique",command=graphe, bg = 'black', fg = 'red') #bouton qui tracera la courbe de prédominance de l'acide choisit.
b2.grid(row=8) #ligne à laquelle se situe le bouton2
bfin=Button(root2,text="Quitter",command=root2.quit, bg ='green', fg = 'blue') bfin.grid(row=9) #ligne à laquelle se trouve le bouton de fin.
root2.mainloop()
# éliminer la fenêtre : root2.destroy()
#fenêtre n°1 dans laquelle on donnera le nombre de proton de la molécule afin de pouvoir donner les valeurs nécessaires à ceux-ci.
root1= Tk()
txtn = Label(root1, text = “Combien de proton peut perdre la molecule ?”, fg = 'red') txtn.grid(row=0) champn=Entry(root1) #demande le nombre de H+ de la molécule (1, 2 ou 3) champn.grid(row=1) #ligne à laquelle apparaitra le champ à complèter
txtc = Label(root1, text = “Concentration (mol/L) :”, fg ='orange') txtc.grid(row=2) champc = Entry(root1) #demande la concentration de la molécule. champc.grid(row=3) #ligne à laquelle apparaîtra le champ où l'on donnera la concentration
b1 = Button(root1, text = 'Valeurs des pKa ?',command = fenetre2,fg = 'red', bg = 'black') #bouton qui ouvrira la fenêtre n°2 afin de rentrer les autres paramètres.
b1.grid(row = 4)
b2 = Button(root1, text = 'Quitter', command = root1.quit,fg = 'red', bg = 'black') b2.grid (row = 5) root1.mainloop()
# éliminer la fenêtre : root1.destroy()
#Référence:
# http://w3.umons.ac.be/perso/Villers.Didier/wiki/doku.php?id=floss:python # http://www.siteduzero.com/informatique/python/tutoriels # https://moodle.umons.ac.be/course/view.php?id=107 # http://www.inforef.be/swi/download/apprendre_python3.pdf
</sxh>