clc
clear

a=80*10^-3 //alésage m
c=84*10^-3 //course m

cylunit=((%pi*a^2)/4)*c //cylindrée unitaire
cyl=cylunit*4

//cherchons V1 et V2 connaissant les deux formules de Pic et de la cyl
V1=1782713*10^-9 //m^3
V2=93827*10^-9 //m^3

Tc=V1/V2 //taux compression

//car Diesel
V4=V1

//on a P1=Patm, T1 = temp ambiante (en K)

Pa=101325
Y=1.4
Cp=1000
Cv=714
P1=101325
T1=15+273.15 //en K

//après calcul

P2=(P1*V1^Y)/(V2^Y) //Pascals, 62.5 bars

T2=(T1*V1^(Y-1))/V2^(Y-1)

P3=P2


//trouvons T3
r=286.8
PCi=43*10^6
rapport_stoechio=14.7
m_air=(P1*V1)/(r*T1)
m_fuel=m_air/rapport_stoechio
Qin=m_fuel*PCi
T3=Qin/((m_air+m_fuel)*Cp)+T2

V3=V2*T3/T2

P4=(P3*V3^Y)/(V4^Y)

T4=P4*T1/P1

Td=V4/V3 //taux détente


rendement=1-(1/Y)*( 1/Td^Y-1/Tc^Y)/(1/Td-1/Tc)



disp ("Qin="+string(Qin)+" J")

//Verifions...

Q2_3=m_air*Cp*(T3-T2)
Q4_1=m_air*Cv*(T1-T4)
Qtot=Q2_3+Q4_1


W2_3=-P2*(V3-V2)
W1_2=(P2*V2-P1*V1)/(Y-1)
W3_4=(P4*V4-P3*V3)/(Y-1)

Wtot=W2_3+W1_2+W3_4

Qtot

rpm=2300

//puissancce moyenne lors de la détendre =MEP
//MEP=Mild effective pressure


MEP=-Wtot/(V4-V3)


P=MEP*cyl*rpm/120

xset('window',1)


// de 1 à 2

cst1=P1*(V1^Y)
volume12=[V2:0.00000001:V1]
pression12=cst1./(volume12^Y)
plot(volume12,pression12,'r')


// de 2 à 3

volume23=[V2,V3]
pression23=[P2,P3]
plot(volume23,pression23,'r')

// de 3 à 4
cst2=P3*(V3^Y)
volume34=[V3:0.00000001:V4]
pression34=cst2./(volume34^Y)
plot(volume34,pression34,'r')

// de 4 à 1
pression41=[P1,P4]
volume41=[V1,V4]
plot(volume41,pression41,'r')



disp ("Qtot="+string(Qtot)+" J")
disp ("Wtot="+string(Wtot)+" J")



//graphique de la puissance en fonction de la température//////////////////////////////////////////////////


xset('window',2)
i=210
while (i<320)

    i=i+0.25
    T1=i
    T2=(T1*V1^(Y-1))/V2^(Y-1)
     m_air=(P1*V1)/(r*T1)
    m_fuel=m_air/rapport_stoechio

    Qin=m_fuel*PCi
    T3=(Qin/(m_air+m_fuel)*Cp)+T2
    V3=V2*T3/T2
    P3=P2
    P4=(P3*V3^Y)/(V4^Y)
    T4=P4*T1/P1


    W3_4=(P4*V4-P3*V3)/(Y-1)
    W1_2=(P2*V2-P1*V1)/(Y-1)
    W2_3=-P2*(V3-V2)
    Wtot=W2_3+W1_2+W3_4
    MEP=-Wtot/(V4-V3)
    P=MEP*cyl*rpm/120
    plot (i,P,".r")
end


//remise a niveau////////////////////////////////////////////////////////////////////////////////////////
clc
clear

a=80*10^-3 //alésage m
c=84*10^-3 //course m

cylunit=((%pi*a^2)/4)*c //cylindrée unitaire
cyl=cylunit*4

//cherchons V1 et V2 connaissant les deux formules de Pic et de la cyl
V1=1782713*10^-9 //m^3
V2=93827*10^-9 //m^3

Tc=V1/V2 //taux compression

//car Diesel
V4=V1

//on a P1=Patm, T1 = temp ambiante (en K)

Pa=101325
Y=1.4
Cp=1000
Cv=714
P1=101325
T1=15+273.15 //en K

//après calcul

P2=(P1*V1^Y)/(V2^Y) //Pascals, 62.5 bars

T2=(T1*V1^(Y-1))/V2^(Y-1)

P3=P2


//trouvons T3
r=286.8
PCi=43*10^6
rapport_stoechio=14.7
m_air=(P1*V1)/(r*T1)
m_fuel=m_air/rapport_stoechio
Qin=m_fuel*PCi
T3=Qin/((m_air+m_fuel)*Cp)+T2

V3=V2*T3/T2

P4=(P3*V3^Y)/(V4^Y)

T4=P4*T1/P1

Td=V4/V3 //taux détente


rendement=1-(1/Y)*( 1/Td^Y-1/Tc^Y)/(1/Td-1/Tc)



disp ("Qin="+string(Qin)+" J")

//Verifions...

Q2_3=m_air*Cp*(T3-T2)
Q4_1=m_air*Cv*(T1-T4)
Qtot=Q2_3+Q4_1


W2_3=-P2*(V3-V2)
W1_2=(P2*V2-P1*V1)/(Y-1)
W3_4=(P4*V4-P3*V3)/(Y-1)

Wtot=W2_3+W1_2+W3_4

Qtot

rpm=2300

//puissancce moyenne lors de la détendre =MEP
//MEP=Mild effective pressure


MEP=-Wtot/(V4-V3)


P=MEP*cyl*rpm/120








//graphique de la puissance en fonction de la pression ///////////////////////////////////////////////////////////////////////////////////////


xset('window',3)

i=50000

while i<(152000)


    i=i+400
    P1=i
    P2=(P1*V1^Y)/(V2^Y)

    T2=(T1*V1^(Y-1))/V2^(Y-1)
     m_air=(P1*V1)/(r*T1)
    m_fuel=m_air/rapport_stoechio

     Qin=m_fuel*PCi
    T3=(Qin/(m_air+m_fuel)*Cp)+T2
    V3=V2*T3/T2
    P3=P2
    P4=(P3*V3^Y)/(V4^Y)
    T4=P4*T1/P1

    W3_4=(P4*V4-P3*V3)/(Y-1)
    W1_2=(P2*V2-P1*V1)/(Y-1)
    W2_3=-P2*(V3-V2)
    Wtot=W2_3+W1_2+W3_4
    MEP=-Wtot/(V4-V3)
    P=MEP*cyl*rpm/120
    plot (i,P,".r")
end

//remise a niveau ///////////////////////////////////////////////////////////////////////////////////////
clc
clear

a=80*10^-3 //alésage m
c=84*10^-3 //course m

cylunit=((%pi*a^2)/4)*c //cylindrée unitaire
cyl=cylunit*4

//cherchons V1 et V2 connaissant les deux formules de Pic et de la cyl
V1=1782713*10^-9 //m^3
V2=93827*10^-9 //m^3

Tc=V1/V2 //taux compression

//car Diesel
V4=V1

//on a P1=Patm, T1 = temp ambiante (en K)

Pa=101325
Y=1.4
Cp=1000
Cv=714
P1=101325
T1=15+273.15 //en K

//après calcul

P2=(P1*V1^Y)/(V2^Y) //Pascals, 62.5 bars

T2=(T1*V1^(Y-1))/V2^(Y-1)

P3=P2


//trouvons T3
r=286.8
PCi=43*10^6
rapport_stoechio=14.7
m_air=(P1*V1)/(r*T1)
m_fuel=m_air/rapport_stoechio
Qin=m_fuel*PCi
T3=Qin/((m_air+m_fuel)*Cp)+T2

V3=V2*T3/T2

P4=(P3*V3^Y)/(V4^Y)

T4=P4*T1/P1

Td=V4/V3 //taux détente


rendement=1-(1/Y)*( 1/Td^Y-1/Tc^Y)/(1/Td-1/Tc)



disp ("Qin="+string(Qin)+" J")

//Verifions...

Q2_3=m_air*Cp*(T3-T2)
Q4_1=m_air*Cv*(T1-T4)
Qtot=Q2_3+Q4_1


W2_3=-P2*(V3-V2)
W1_2=(P2*V2-P1*V1)/(Y-1)
W3_4=(P4*V4-P3*V3)/(Y-1)

Wtot=W2_3+W1_2+W3_4

Qtot

rpm=2300

//puissancce moyenne lors de la détendre =MEP
//MEP=Mild effective pressure


MEP=-Wtot/(V4-V3)


P=MEP*cyl*rpm/120