Aspect physique » History » Version 10

Pierre Guinault, 03/27/2013 10:27 AM

1 1 Pierre Guinault
h1. Aspect physique
2 1 Pierre Guinault
3 2 Pierre Guinault
h2. Définitions utiles dans le domaine du quadricoptère :
4 2 Pierre Guinault
5 2 Pierre Guinault
* Angle de calage : angle entre la corde et le plan de rotation de l'hélice
6 2 Pierre Guinault
* Traînée : résistance de l'air sur l'hélice sur une hélice en rotation.
7 2 Pierre Guinault
* Vitesse induite : vitesse communiquée à la molécule d'air lors de son passage au travers du disque rotor.
8 2 Pierre Guinault
9 2 Pierre Guinault
h2. Effet à prendre en compte pour le vol d'un drone.
10 2 Pierre Guinault
11 2 Pierre Guinault
_Cette partie relate des contraintes de l'environnement extérieur au drone et qui agisse dessus._
12 2 Pierre Guinault
13 2 Pierre Guinault
* L'effet de couple : Du fait que le l'hélice tourne dans un sens, le moteur tend à tourner dans le sens opposé en raison de la traînée induite par les pales.	Il en résulte une légère inclinaison du moteur concerné. L'effet de couple dépend de la masse du l'hélice, de son pas, de son diamètre et de son régime, mais aussi et surtout à la puissance disponible : plus la puissance est grande, plus l'effet de couple est grand.
14 2 Pierre Guinault
	L'inclinaison de l'axe du moteur en dû à l'effet de couple dépend de la résistance qui lui est opposé. Pour réduire cet effet, on peut appliquer deux types de « contre-mesures » :
15 2 Pierre Guinault
-	les mesures passives qui ont lieu lors de la conception.
16 2 Pierre Guinault
-	les mesures actives qui sont celle que l'on peut induire durant le vol (comme le report de poids d'un côté ou de l'autre, un réglage asymétrique du Trim, une réduction de poussée...). 
17 2 Pierre Guinault
18 2 Pierre Guinault
* La précession gyroscopique : Une hélice qui tourne agit comme un gyroscope, qui initialement tend à résister à toute force qui pourrait lui faire changer d'axe de rotation.
19 2 Pierre Guinault
20 2 Pierre Guinault
!precession_gyroscopique.png!
21 3 Pierre Guinault
_F la force instantanée appliquée au système en rotation, et R la résultante à 90 degrés dans le sens de rotation que le système va réellement "ressentir"._
22 3 Pierre Guinault
23 3 Pierre Guinault
* La poussée asymétrique des pales : celle-ci a trois facteurs qui sont :
24 3 Pierre Guinault
-	Elle se manifeste lorsque l'hélice n'est pas dans le plan horizontal, chacune des pales lorsqu'elle descend est entraînée vers l'avant, et lorsqu'elle monte vers l'arrière (phénomène dû bien entendu à la force de pesanteur).
25 3 Pierre Guinault
D'où en descente : vitesse air = vitesse sol + vitesse de rotation + vitesse d'entraînement.
26 3 Pierre Guinault
Et en montée : vitesse air = vitesse sol + vitesse de rotation – vitesse d'entraînement. 
27 1 Pierre Guinault
!poussee_asymetrique_pales.png!
28 4 Pierre Guinault
_Le poids de la pale P vient dans un cas s'additionner à la trainé f, et dans l'autre cas s'en soustraire. Il y a donc asymétrie dès lors que le drone n'est pas à l'horizontal._
29 4 Pierre Guinault
-	Du fait de l'inclinaison du moteur, l'angle d'attaque de chaque pale est plus élevé dans la phase de descente que dans la phase de montée. Cela provoque plus de poussée d'un côté de l'hélice que de l'autre.
30 4 Pierre Guinault
-	Les deux facteurs précédents induisent le fait que le côté où la pale descend est exposée à un flux d'air dit « propre » alors que l'autre subit un masquage aggravant. Cela provoque une poussée accrue du coté de descente, et réduit de l'autre, ce qui s'ajoute encore aux facteurs précédents, et on a donc une boucle.
31 1 Pierre Guinault
32 4 Pierre Guinault
A noter que la principale différence entre la précession gyroscopique et la poussée asymétrique des pales est que la première est une force momentanée qui disparaît dès le prochain changement d'assiette ; alors que la seconde est à peu de chose près permanente, de sorte que la direction réelle de l'appareil est toujours décalée par rapport à la direction souhaitée.
33 4 Pierre Guinault
34 4 Pierre Guinault
Il est également intéressant de noter que la poussée asymétrique des pales et la précession gyroscopique sont rapidement contrariées par une réduction de puissance et par l'évitement de changement d'assiette trop brut.
35 4 Pierre Guinault
36 4 Pierre Guinault
h2. Analyse mécanique
37 4 Pierre Guinault
38 4 Pierre Guinault
_Afin de déterminer quelles sont les contraintes internes au quadricoptère, il nous a fallu réaliser une  analyse mécanique, et déterminer, si ces contraintes engendrent des problèmes,  si une approche mécanique (classique) pouvait les paliers._
39 2 Pierre Guinault
40 5 Pierre Guinault
h3. Théorie de l'élément de la pale (Froude)
41 5 Pierre Guinault
	
42 5 Pierre Guinault
La théorie de l'élément de la pale a été initialement développée par Froude qui assimile dans cette théorie l'élément de la pale à une plaque mince, soumise à un vent relatif. Cette théorie a été améliorée depuis par Drzewiecki qui calcule les caractéristiques aérodynamiques des hélices.
43 5 Pierre Guinault
44 5 Pierre Guinault
La pale est ici découpée en un nombre fini de tranches (éléments) à l'aide de surfaces cylindrique.
45 5 Pierre Guinault
46 5 Pierre Guinault
On suppose que l'écoulement dans un anneau, limité par deux surfaces voisines, est indépendant de l'écoulement dans les autres anneaux. Il est ainsi possible d'obtenir de façon indépendante les forces de traînée dD, et les forces de portance dL appliquée sur chaque tranche, sans tenir compte de l'écoulement dans les éléments voisins. 
47 5 Pierre Guinault
48 5 Pierre Guinault
Enfin, en intégrant les efforts élémentaires, on peut calculer es caractéristique aérodynamique intégrales du rotor. 
49 6 Pierre Guinault
!froude.png!
50 6 Pierre Guinault
51 6 Pierre Guinault
Pour obtenir l'effort élémentaire, chaque élément de la pale est donc représenté comme une aile cylindrique (aile à allongement infini), avec une surface dS = c*dr.
52 6 Pierre Guinault
L'aile est soumise au vent résultant W, qui est une composition de la vitesse du vent à l'infini Vinf et celle de rotation U = Ωr.
53 6 Pierre Guinault
<pre>
54 6 Pierre Guinault
W = √(Vinfini² + Ω²r²)    , c’est la vitesse relative en m/s (pas le travail)
55 6 Pierre Guinault
</pre>
56 6 Pierre Guinault
57 10 Pierre Guinault
_Vinfini est la vitesse à l’amont en m/s.
58 6 Pierre Guinault
Ω = vitesse angulaire en rad /r.
59 6 Pierre Guinault
r = rayon local, correspond au r sur le schéma précédent
60 10 Pierre Guinault
_
61 6 Pierre Guinault
L'angle d'incidence α est calculé à partir de l'angle d'écoulement φ = arctan(Vinf / U), et l'angle de calage de la pale β : α =(φ-β).
62 6 Pierre Guinault
La portance et la traînée peuvent maintenant être calculées à partir des formules :
63 6 Pierre Guinault
64 6 Pierre Guinault
<pre>
65 6 Pierre Guinault
dL=½ .ρ .W^2."CL" .α .c .dr
66 6 Pierre Guinault
dS= ½ .ρ .W² ."CD" .α .c .dr
67 6 Pierre Guinault
</pre>
68 6 Pierre Guinault
69 6 Pierre Guinault
Où  CL  et  CD  sont les coefficients, respectivement, de traînée et de portance du profil d'aile.
70 6 Pierre Guinault
71 6 Pierre Guinault
Les valeurs de ces coefficients dépendent de l'angle d'incidence entre la corde du profil et la vitesse résultante. De façon générale, ces coefficients sont obtenus à partir d'essais de profils ou des simulations numériques prenant en compte la viscosité du milieu ; dans notre cas la méthode de recherche des coefficients par essais nous a semblé être plus facilement réalisable que la seconde méthode, qui nécessiterait des moyens d'envergures à son application. 
72 6 Pierre Guinault
73 7 Pierre Guinault
!pale.png!
74 6 Pierre Guinault
Afin d'obtenir la force tangentielle et la force axiale, appliquée à chaque élément de la pale, on  projette les forces de traînée et de portance sur le plan de rotation et sur l'axe de rotation :
75 7 Pierre Guinault
<pre>
76 1 Pierre Guinault
dFa = ½ .ρ .W² .[ "CL" .α .cos β + "CD" .α .sin β ] .c .r .dr 
77 7 Pierre Guinault
dFt= ½ .ρ .W² .[ "CL" .α .sin β - "CD" .α .cos β ] .c .r .dr 
78 7 Pierre Guinault
</pre> 
79 7 Pierre Guinault
80 8 Pierre Guinault
Finalement l'intégration de ces forces le long de la pale, du pied à l'extrémité permet d'obtenir la force axiale et la puissance du rotor.
81 8 Pierre Guinault
En prenant en compte le nombre de pale N et la vitesse angulaire du rotor Ω, on obtient la force axiale Fa et la poussée P (en Watts) :
82 5 Pierre Guinault
83 8 Pierre Guinault
<pre>
84 8 Pierre Guinault
Fa = ½ .ρ . N .∫ [rp,R] ( W² .[ "CL" .α .cos β + "CD" .α .sin β ] .c .r .dr )
85 8 Pierre Guinault
P = ½ .ρ .N .∫ [rp,R] ( W² .[ "CL" .α .sin β - "CD" .α .cos β ] .c .r .dr )
86 8 Pierre Guinault
</pre>
87 8 Pierre Guinault
88 8 Pierre Guinault
La théorie de l'élément de la pale donne des résultats satisfaisant à condition que les vitesses induites soient faibles.
89 9 Pierre Guinault
h3. 
90 9 Pierre Guinault
Théorie de l'hélice
91 9 Pierre Guinault
92 9 Pierre Guinault
Comme la théorie de l'élément de la pale ne nous convient pas totalement nous nous sommes alors penchés sur la théorie de l'hélice, plus générale, et qui ne prend pas en compte le profil des pales.
93 9 Pierre Guinault
Soient (i) l'angle d'incidence de la pale, (θ) l'angle de pas et (β) l'angle de conicité. 
94 9 Pierre Guinault
!helice1.png!
95 9 Pierre Guinault
!helice2.png!
96 1 Pierre Guinault
!helice3.png!
97 1 Pierre Guinault
98 10 Pierre Guinault
Plus l'angle d'incidence augmente, plus ρ us la portance augmente, jusqu'à une valeur limite où le profil décroche. Au décrochage, la portance s'écroule et la traînée croît brutalement.
99 10 Pierre Guinault
La vitesse en chaque point de la pale augment proportionnellement à l'éloignement de l'axe de rotation, et la portance croît en fonction du carré de la vitesse.
100 10 Pierre Guinault
La vitesse de rotation de la pale vaut : ωR, V est la vitesse de la pale et Vt la vitesse de translation.
101 10 Pierre Guinault
102 10 Pierre Guinault
<pre>
103 10 Pierre Guinault
La résultante aérodynamique en hauteur Rz est : Rz =  ½ . ρ . S . V² . CZ . ẑ
104 10 Pierre Guinault
</pre>
105 10 Pierre Guinault
106 10 Pierre Guinault
_ρ = masse volumique du fluide du milieu en Kg/m3
107 10 Pierre Guinault
S = maître couple en m²
108 10 Pierre Guinault
V = vitesse en m/s
109 10 Pierre Guinault
CZ = coefficient portance_
110 10 Pierre Guinault
111 10 Pierre Guinault
h3. Conclusion de l'analyse
112 10 Pierre Guinault
 	
113 10 Pierre Guinault
Dans le cadre de notre projet, nous avons préféré substitué la résolution mécanique de l'instabilité du drone par une résolution mathématique utilisant un asservissement PID.
114 10 Pierre Guinault
115 10 Pierre Guinault
Ce choix est dû à la quantité de traitement qu'il n'est pas possible de réaliser avec le matériel présent sur le drone ou avec le budget qu'il nous est fourni. De plus la complexité des équations demanderait trop de temps pour pouvoir avoir un vol stabilisé dans les délais impartis.
116 10 Pierre Guinault
117 10 Pierre Guinault
Dans la prochaine partie nous aborderons donc la solution de stabilisation par PID que nous avons implémenté de façon concrète, et nous étudierons également la recherche des coefficients nécessaires à sa mise en place.
118 5 Pierre Guinault
119 5 Pierre Guinault
120 1 Pierre Guinault
_Dernière mise à jour : 27/03/2013 10h13
121 1 Pierre Guinault
Auteur : Benjamin RIVRAY_