Etude cinématique – statique d’un
moteur thermique a deux temps
{
⃗
z 1=−cos ( α ) ⃑z +sin (α ) ⃗y
⃗y 1=sin(α )⃗z −cos( α )
⃗
AB=l ⃗
z 1=lsin ( α ) ⃗y−l cos ( α ) ⃗z
⃗
BC =b ¿
⃗ 2
BC=b−2b cos ( α ) ⃗z + sin ( 2 α ) ⃗y
→
⃗
AC =⃗ BC =( b−2 b cos 2 ( α )−l cos ( α ) ) ⃗z + ( l sin ( α ) +sin ( 2 α ) ) ⃗y
AB+ ⃗
Etude cinématique
{ |}{ } { |}
α̇ ⃗x 0 α̇ ⃗x 0
{ } 0
ϑ ( ) A= 0 0
1
0 0 A
1
ϑ ( ) B= 0 0
2
0 0 B
{ |}{ } { | }
ω23 ⃗x 0 0 0
{ } 2
ϑ ( ) C= 0 0
3
0 0 C
3
ϑ ( ) C = ω30 ⃗y V 30
0
0 0 C
Liaisons L(1/2) dans A :
⃗ ⃗ ( B ∈1/2 ) +ω 12 ^ ⃗
V ( A ∈1/2 )=V BA
¿ 0⃗ + α̇ ⃗x ^ (−lsin ( α ) ⃗y+l cos ( α ) ⃗z )
¿−l α̇ sin ( α ) ⃗z −l α̇ cos ( α ) ⃗y
{| }
α̇ ⃗x 0
{ } 1
ϑ ( ) A = 0 −l α̇ cos ( α ) ⃗y
2
0 −l α̇ sin ( α ) ⃗z A
Liaison L 2/3 dans A :
⃗
V ( A ∈2/3 ) =⃗
V ( c ∈ 2/3 )+ ω23 ^ ⃗
CA
¿ 0⃗ + ω23 ( b−2b cos 2 ( α ) −lcos ( α ) ) ⃗y −ω23 ( lsin ( α ) +sin ( 2 α ) ) ⃗y
{ | }
ω 23 ⃗x 0
{ } 2
ϑ ( ) A=
3 0 ω23 ( b−2 b cos ( α )−l cos ( α ) ) ⃗y
0
2
ω 23 ( l sin ( α )+ sin ( 2 α ) ) ⃗z A
Liaison L 3/0 dans A :
⃗
V ( A ∈3 /0 )=⃗
V ( c ∈ 3/0 ) + ω30 ^ ⃗
CA
V ( c ∈ 3/0 ) + ω30 ( b−2 b cos 2 ( α ) −lcos ( α ) ) ⃗x
¿⃗
{ | }
ω 23 ⃗x ω30 ( b−2 b cos2 ( α )−l cos ( α ) ) ⃗x
{ } 2
ϑ ( ) A= 0
3
0
⃗
V ( c ∈ 3/0 )
0 A
{ϑ ( 01 )} +{ϑ ( 12 )} +{ϑ ( 23 )} + {ϑ ( 32 )} ={0⃗ }
A A A A
On déduire le système suivant :
{ }
2 α̇ +ω 23=0
ω30=0
0=0
ω30 ( b−2 b cos ( α ) −lcos ( α ) ) =0
2
⇒
⃗ ( c ∈3 /0 ) +ω23 ( b−2 b cos2 ( α )−l cos ( α ) )−l α̇ cos ( α )=0
V
ω23 ( l sin ( α ) +sin ( 2 α ) ) −l α̇ sin ( α )=0
En résume qu’on a un système de Nc=4 équations avec Ic=5 inconnues la
mobilité du système est donc :
m=I c −N c =5−4=1
{ }
2 α̇ +ω 23=0
ω30=0
⇒ ⃗ ( c ∈3 /0 ) +ω23 ( b−2 b cos2 ( α )−l cos ( α ) )−l α̇ cos ( α )=0
V
ω23 ( l sin ( α ) +sin ( 2 α ) ) −l α̇ sin ( α )=0
Le degré d’hyper statisme de système est donc de H=2
H = Mc +6 – Ic= 2+6-8= 0
Loi entre sortie :
AC =( b−2 b cos 2 ( α )−l cos ( α ) ) ⃗z + ( lsin ( α ) +sin ( 2 α ) ) ⃗y
⃗
AC‖=√ ( b−2 b cos ( α )−l cos ( α ) ) + ( lsin ( α ) +sin ( 2 α ) )
Ou y=‖⃗ 2 2 2
Loi cinématique entrée sortie :
dα
ω=
dt
La loi d’entrée sortie entre la vitesse de rotation ω et la vitesse de position V
se déduit en dérivant la loi géométrique :
3
dy sin (α )(8 b ² cos(α )−16 b ² cos (α )+2 bl−12 bl cos ² (α )+ 4 l)
V ( t )= =
dt
√ 2
2 ( b−2 b cos ( α )−l cos ( α ) ) + ( l sin ( α )+ sin (2 α ) )
2 2
Etude statique :
{ | } { ( )} { | }
X 01 0 X 12 0
{ ( )}
τ
0
1
A = Y 01 M 01
Z 01 N 01 A
τ
1
2
B = Y 12 M 12
Z 12 N 12 B
{ | }{ } { | }
X 23 0 X 30 L30
{ } 2
τ ( ) C = Y 23 M 23
3
Z 23 N 23 C
3
τ ( ) C= 0 0
0
Z 30 N 30 C
L’équation d’équilibre pour chaque pièce :
{ 1} { 2}
(
τ
0
) + τ ( ) + { τ ( effort exterieur ) } = ⃗0
A
1
A A
{ 2} { 3}
(
τ
1
) + τ ( ) + { τ ( effort exterieur ) } =⃗0
B
2
B B
{ ( )} { ( )}
τ
2
3 C + τ
3
0 C + { τ ( effort exterieur ) } C = ⃗0
Equilibre de vilebrequin 1 au point A
{ 1} { 2}
(
τ
0
) + τ ( ) + { τ ( effort exterieur ) } = ⃗0
A
1
A A
{|}
X 12 0
{ ( )}
τ
1
2
B = Y 12 M 12
Z 12 N 12 B
Avec : ⃗
AB=l ⃗
z 1=lsin ( α ) ⃗y−l cos ( α ) ⃗z
⃗
M ( A ∈ 1/2 )=⃗
M ( B ∈ 1/2 )+ R 12 ^ ⃗
BA
{| }
X 12 L cos(α )Y 12+l sin ( α ) Z12
{ ( )}
τ
1
2 A= Y 12
Z12
M 12−lcos ( α ) X 12
N 12−l sin(α ) X 12 A
{| }{| } {| }
X 01 0 X 12 L cos( α )Y 12+l sin ( α ) Z12 0 −C
Y 01 M 01 + Y 12 M 12−lcos ( α ) X 12 + 0 0 =0⃗
Z 01 N 01 A Z12 N 12−l sin(α ) X 12 A
0 0 A
{ }
X 01+ X 12=0
Y 01+Y 12=0
Z 01+ Z 12=0
⇒ L cos (α)Y 12+ lsin ( α ) Z 12−C=0
M 01+ M 12−l cos ( α ) X 12=0
N 01 + N 12−lsin ( α ) X 12=0
Equilibre de la bielle 2 au point B :
{ | }{ } { | }
X 12 0 X 23 0
{ ( )}
τ
1
2
B = Y 12 M 12
Z 12 N 12 B
2
τ ( ) C = Y 23 M 23
3
Z 23 N 23 C
Avec : ⃗ 2
BC =b−2 b cos ( α ) ⃗z +sin ( 2 α ) ⃗y
⃗
M ( B ∈2 /3 )=⃗
M (C ∈2/ 3 ) + R 23 ^ ⃗
BC
{| }
X 23 b−2 b cos 2 ( α ) Y 23−sin ( 2 α ) Z23
{ ( )}
τ
2
3 B = Y 23
Z 23
2
M 23−b−2 b cos ( α ) X 23
N 12−sin ( 2 α ) X 12 B
{| }{ | }
X 12 0 X 23 b−2 b cos 2 ( α ) Y 23−sin ( 2 α ) Z 23
Y 12 M 12 + Y 23 2
M 23−b−2b cos ( α ) X 23 = ⃗0
Z 12 N 12 B Z 23 N 12−sin ( 2 α ) X 12 B
{ }
X 12 + X 23=0
Y 12 +Y 23=0
Z 12 +Z 23=0
2
⇒ b−2 b cos ( α ) Y 23−sin ( 2 α ) Z 23=0
2
M 12+ M 23−b−2 b cos ( α ) X 23=0
N 12 + N 23−sin ( 2 α ) X 12=0
Equilibre de la piston 3 au point 3 :
{ 3} { 0}
(
τ
2
) + τ ( ) + { τ ( effort exterieur ) } = ⃗0
C
3
C C
{ | } { | } {| }
X 23 0 X 30 L30 0 0
M
Y 23 23 + 0 0 + F =0⃗
0 −⃗
Z 23 N 23 C Z 30 N 30 C 0 0 C
{ }
X 23+ X 30=0
Y 23+ Y 30=0
Z 23+ Z 30=0
⇒ L30=0
M 23−⃗ F =0
N 30 + N 23=0
Étude thermodynamique
Dans cette étude on va utiliser un octane de formule brute C 8 H 18 de masse
volumique ρ=720 Kg ∙ m−1
C 8 H 18 +12.5 O2 → 8 C O2 +9 H 2 O
Enthalpie de la liaison :
On réalise un cycle thermodynamique faisant intervenir la relation
d’atomisation de l’octane
On donne :
Enthalpie de liaison (kJ/mol) : C-C : 415 ; C-H : 345 ;
Enthalpie standards de formation (kJ/mol): carbone C(g) : 716,7 ; H(g) : 218 ;
∆ atome H °=7 ∆liaison H ° ( C−C )+18 ∆liaison H ° ( C−H )
¿ 9885 Kj /mol
∆ f H ° ( C 8 H 18) =8 ∆ f H ° ( C ) +18 ∆ f H ° ( H )=9657 Kj /mol
−1
∆ f H ° ( octane )=9657.6−9885=−227.4 Kj ∙ mol
Enthalpie standard de combustion de l’octane :
C 8 H 18 +12.5 O2 → 8 C O2 +9 H 2 O
3
∆ comb H °=−∆ f H ° ( octane ) +12.5 ∆ f H ° ( O2 ) +8 ∆f H ° ( C O2 ) +9 ∆ f H ° ( H 2 O )=−5.1 ∙10 Kj /mol
Chaleur libérée par la combustion d’un Cm3 d’octane :
La mass d’un Cm3 d’octane m=0.72 g
Et sa mass molaire est M =114 g/mol
0.72 −3
Quantité de matière dans 1Cm3 d’octane 114 =6.31∙ 10 mol
−3 3 −3
Q=6.31∙ 10 ∙5.1 ∙10 =32.2 Kj ∙ cm
Calcule thermo dynamique du moteur
deux temps
Le mélange carburé est ( air – essence [C 8 H 18] – huile )
En suppose qu’on brule que d’essence on donne :
Cylindrée de moteur est :V A =101.021 cm3
3 5
T A=300 K V B =6.507 cm P a=10
Le course du piston est de 105mm
Vitesse moyen de piston sur un cycle
Un tour de vilebrequin correspend a un cycle thermodynamique
A 6000 tr/min
60 −3
soit 1 tour effectué en : 6000 =9.99 ∙10 s=9.99 ms c’est la durée d’un cycle
Un cycle correspond à un aller-retour soit alors :
−2
2 ∙105=2.1∙ 10 m
−2
2.1∙ 10 −1
Vmoy = −3
=2.1 m∙ s
9.99 ∙ 10
Vitesse moyenne quadratique des molécules :
√ V −2= √ 3∙ K b ∙ T
m
avec : m=
M
=
0.029
N A 6.02 ∙ 1026
=4.81 ∙10−26
( )
−23 1
3 ∙1 , 38 ∙10 ∙ 300 2
−26
⇒ 4.81 ∙ 10
Valeur très supérieure à la vitesse du piston. La transformation est quasi-
statique et on peut la considérer comme réversible.
Transformation isotopique :
S'il n'y a pas d'échange de chaleur lors des transformations CD et AB, celles-ci
peuvent être modélisée par des adiabatiques réversibles. PV δ =Cte Pour un gaz
parfait et pour une transformation adiabatique réversible. Dans une
transformation adiabatique, il n'y a pas d’échange de chaleur donc l'entropie
échangée est nulle.si la transformation est en plus réversible, l'entropie crée
est nulle.la variation d'entropie est donc nulle : transformation isentropique.
Ce résultat est indépendant de la nature du système
Taux de compression
( )
1
δ VB
δ PB δ
PB ∙ V =P A ∙ V et a= =
C B =0.064
VA PA
Travail du moteur au cour d’un cycle :
La détente CD est adiabatique réversible donc :
W CD =∆ U=n∙ C Vn ∙ ∆ T =n ∙C Vn ∙(T D−T C )
La compression AB est adiabatique réversible donc :
W AB =∆ U =n ∙C Vn ∙ ∆ T =n ∙ CVn ∙(T B−T A )
W Cycle =n ∙C Vn ∙(T D−T C +T B−T A )
n∙R
¿ :n ∙ CVn= et P A ∙V D =n ∙ R ∙ T A
( δ−1 )
P A ∙V D
W Cycle =
( ( δ−1 ) ∙T A )(T D −T C +T B −T A )
Chaleur libérée per la combustion
Pour la transformation BC(W=0) et ∆ U =Qchaud =n∙ C Vn ∙(T c −T b)
PA ∙ V D
Qchaud =
( ( δ−1 ) ∙ T A )(T C −T B)
Calcule de rendement :
travail fourni a l ' exterieur
η=
chaleur liberee par la combustion
−W Cycle (T D −T C +T B −T A ) 1−(T D −T A )
η= = =
Qchaud T c −T b T C −T B
δ P∙V
Ou pour une adiabatique reversible PV =Cte et T =Cte
' ( δ−1 )
d ou T ∙V =Cte
( )VC
VD
δ−1
=
T D T A T D −T A
= =
T C T B T C −T B
1
η=1− (δ−1)
a
Calcule de chaleur libérée :
η∙ Qchaud
W traval =η ∙Q chaud ∨P=
durée du cycle
Calcule de puissance et de couple moteur :
δ
δ δ P A ∙V A
P A ∙ V =P B ∙V
A B P B= δ
⇒ VB
PB =1344782.24 Pa
δ
δ δ PA ∙ V A
PB ∙ V B =PC ∙ V C PC = δ
⇒ VB
Avec ça en peux déduire la puissance moteur :
PC
F C= =6916.66 N
Sp
Puissance en (KW et Cv.din)
F ∙ V moy
Pmoteur = =14.5 Kw
1000
¿ 19.02 Cv . din
Couple moteur :
Pmoteur 6000 tr
C= avec : N =
2∙ π ∙ N min
C=23.10 N ∙m