ANNEXE 4: MODELE DU STICK-SLIP Annexe 4: Modèle du stick-slip La vitesse axiale du baladeur varie beaucoup durant le processus de changement de vitesses. En effet, le baladeur démarre de la position engagée dans une des vitesses, s’accélère, puis se ralentit puis s’arrête pour la synchronisation. Après cela, il accélère de nouveau, déplace la bague de synchronisateur, engage la vitesse suivante avec un choc, puis s’arrête définitivement. Durant tout ce processus, il glisse sur des cannelures et subit une force tangentielle variable. Ainsi, on peut supposer, que du phénomène de stick-slip intervient lors du déplacement. Pour décrire cela, on utilise le modèle suivant, d’après Thomsen [41]. Le baladeur est assimilé à une masse posée sur un tapis roulant, et relié au mur fixe par un ressort et un amortisseur (Fig. A4-1). F représente la force tangentielle, la vitesse du tapis celle axiale du baladeur. Fig. A4-1 Modèle pour étudier le stick-slip [41] Les équations de mouvement du système, en forme adimensionnée: Pour le glissement: .. . .⎛ ⎞x+2β x+x+µ x−vb =0 ⎜ ⎟⎝ ⎠Pour le collage: .. .x=0 ; x+2βvb<µs ; x=vb Les paramètres adimensionnés sont: ~- t=ω0t le temps, K- ω0= la fréquence propre du système, MX- x= le déplacement de la masse, L 254F- L= le facteur de déplacement, KVb- vb= la vitesse d’excitation, ω0Lc- 2β = l’amortissement KM Fig. A4-2 Fonction du coefficient de frottement en fonction de la vitesse relative ...
La vitesse axiale du baladeur varie beaucoup durant le processus de changement de vitesses. En effet, le baladeur démarre de la position engagée dans une des vitesses, s’accélère, puis se ralentit puis s’arrête pour la synchronisation. Après cela, il accélère de nouveau, déplace la bague de synchronisateur, engage la vitesse suivante avec un choc, puis s’arrête définitivement. Durant tout ce processus, il glisse sur des cannelures et subit une force tangentielle variable. Ainsi, on peut supposer, que du phénomène de stickslip intervient lors du déplacement. Pour décrire cela, on utilise le modèle suivant, d’après Thomsen [41]. Le baladeur est assimilé à une masse posée sur un tapis roulant, et relié au mur fixe par un ressort et un amortisseur (Fig. A41).Freprésente la force tangentielle, la vitesse du tapis celle axiale du baladeur.
Fig. A41 Modèle pour étudier le stickslip [41] Les équations de mouvement du système, en forme adimensionnée: Pour le glissement:
Pour le collage:
.. . . ⎛ ⎞ = x+2βx+x+µ⎜x−vb⎟0⎝ ⎠
.. . x=0;x+2βvb<µs;x=vb
Les paramètres adimensionnés sont: t=0tle temps,
K ω0=la fréquence propre du système, M
X x=le déplacement de la masse, L
254
F L=le facteur de déplacement, K Vb vb=la vitesse d’excitation, ω0L c 2β=l’amortissement KM
Fig. A42 Fonction du coefficient de frottement en fonction de la vitesse relative [41] e Le frottement est décrit par une fonction de 3 degré (Fig. A42): 3 ( ) (v)v v µvr=signr⋅µs−κ1r+κ3r
où
. vr=x−vb la vitesse de glissement relative.
µs−µm 3 κ1= ⋅2vm µs−µm 1 κ3= ⋅3 2vm Les paramètres de la fonction sont: µsle coefficient de frottement statique,
µmle coefficient de frottement minimal,
vmla vitesse relative appartenant au coefficient de frottement minimal. Pour faciliter l’étude du mouvement, on placera l’origine du système de coordonnées au point où la force de frottement et la force du ressort sont en équilibre. Les équations décrivant ce point=:
L’équation de l’équilibre des forces:
. .. =:x=0 ;x=0
3 x=−µ(−v)=µ−κv+κvbb s1b3
A ce point, les vitesses de déplacement et d’excitation sont égales: vr=vb
255
En connaissant le point d’équilibre, on introduit la nouvelle variable pour le déplacement: u(t)=x(t)−x
En remplaçant cela dans les équations du mouvement, on obtient la nouvelle forme de l’équation de mouvement: Pour le glissement:
Pour le collage:
.. . u+u+εh(u)=0
.. . 3 uε κbκ εβb;vb u=0 ;+(−1v+3vb)+2v≤0u=
où la fonctionhest la suivante:
2 3 . . . . . . . . 2 h(u)=2βu+µ(u−v)−µ(−v)=2βu+µ(1+sign(u−v))+(−κ+κvb)u−3κvbu+κub b s b1 3 3 3 Le facteurεest théoriquement petit:ε<<1, en pratique on supposeε≈1.
. En fonction de la vitesse relativeu, on distingue trois états de fonctionnement:
•
•
•
. u=glissement sans oscillation, avec une vitesse d’excitation0 : vb.
. u<vb: glissement avec oscillation.
. u≤vb: glissement avec oscillation et collage (stickslip).
A41. Description du glissement avec oscillation
. La condition de base:u<vb. L’équation du mouvement:
La fonctionhest la suivante:
Les coefficients:
.. . u+u+εh(u)=0
2 3 . . . . h(u)=h1u+h2u+h u3
2 h1=2β−κ1+3κ3vb
h2=−3v 3κb
h3=3
256
La solution de l’équation de mouvement a la forme suivante: u=Asinψ
où
ψ=t+θ(t) la phase,
θ(t)– le déphasage initial. Calcul deA(t)et deθ(t)avec la méthode de la moyenne standard:
La solution non triviale:
. 1 3 2 A=−εA(h1+h3A)2 4
. Aθ=0
−4h1 A=3h3 =constEn remplaçant cela dans l’équation de la solution, on obtient une solution périodique. La solution n’est pas stable, siu=0 pourh1<0. Cela donne la condition suivante pour la vitesse d’excitation:
4βvm b<vb1=vm− v13(µs−µm)
Cela veut dire qu’en absence d’amortissement visqueux, l’équilibre est instable pour des vitesses d’excitation inférieures àvm. Donc, l’amortissement visqueux a un effet stabilisateur sur le système. Siβdépasse un seuil, le système est toujours stable: 3 β>(µs−µm) 4 Les conditions de l’existence et de la stabilité de l’oscillation périodique sont: h1<0; h3>0 Preuve de la vérité de la deuxième condition: Etant donnéµs>µm, la deuxième condition est vraie: µs−µm 1 h3=κ3= ⋅ >0 3 2m v L’équation modifiée de l’amplitude de déplacement:
2 vb4βv2 2 m A=2v1−( )− =2vb1−vbm vm3(µs−µm)
Ainsi, on peut définir des valeurs limites pour l’oscillation pure. L’oscillation pure est présente si la vitesse d’excitationvbest entre les valeurs limites suivantes: vb0<vb<vb1
257
Quand la vitesse d’excitationvb diminue à partir d’une valeur grande, le glissement sans oscillation se transforme en oscillation pure à la valeurvb1. L’oscillation pure dure jusqu’à ce que la vitesse d’excitation atteintvb0, vitesse de l’apparition du stickslip. Similairement, l’amplitude de l’oscillation commence à augmenter aprèsvb1. Il atteint sa valeur maximale à vb0. A ce point:
4 vb0=vb1 5 De cela, l’amplitude de l’oscillation au début de la phase de stickslip:
A A v =max=b0= vb=vb0
4βvm 4 v1−m 5 3(µs−µm)
A42. Description de l’oscillation stickslip
La condition de base:vb<vb0. Dans la discussion, on s’occupe séparément de la phase glissement et de la phase collage. A421. Phase glissement Pour décrire la phase de glissement, on choisit pour moment initial la fin de la phase de
. collage, où la vitesse de la masse est égale à celle d’excitation:u=vb. A ce moment, la force de frottement et celle du ressort sont égales. L’équation du mouvement au moment initial:
. 3 t=0:u(0)=−2εβvb+ε(κ1v−κ3vb) ;u(0)=vbb L’équation de mouvement sera la même que pour le mouvement d’oscillation pure:
.. . u+u+εh(u)=0
Avec ces conditions initiales, l’équation est nonlinéaire. Pour la résoudre, on fait appel à des méthodes approchées. Appliquons la méthode des perturbations. Décomposons la variable: u(t)u0(t)εu1(t)
quandt∈[0,t], oùtsest le début de la nouvelle phase de collage, etε<<1un paramètre. s
. Remplaçons l’équationu(t)=v dans celle précédente, et arrangeons les membres pour s b
obtenirts:
. ⎛3π⎞ vb−u⎜ ⎟2 ⎛ ⎞ 3π π ⎝2⎠ ⎛3⎞ ⎜ ⎟ t O t s'+ ⎜= + s− ⎟⎜ ⎟ .. ⎛3⎞ ⎝2⎠ 2π ⎝ ⎠ u⎜ ⎟ ⎝2⎠
. .. Remplaçons les équations deu(t) etu(t) , et négligeons les puissances dets: s s
259
1 13 v+c−c b1 3 3π 2 8 t'= +s 3 2 v−c−2c b1 3 4 Selon Thomsen [41], le pourcentage d’erreur de cette approximation est inférieur à 0,3%. Toutefois, l’erreur peut augmenter en cas de(µs- µm) plus petit. Pour compenser cela, on effectue le même développement en série, mais autour du pointts=2π:
2πc 1 t"=2π−s 3π −v−c−2c b1 3 4 Pour une phase de collage court,ts’est plus précise. Pour une phase de collage plus longue, ts’’ est mieux adaptée. Pour être précis en chaque cas, considérons que la solution la plus précise est la combinaison des deux approximations:
v−v b b0 κ t=t'+e(t"−t')s s s s
4βv m oùκ=. 3(µ−µ) s m L’erreur de cette approximation pourtsest inférieur à 1% pour les conditions suivantes: µ−µ s m vb<vb0;∈[0,11;1];∈[0,07;0,61]; v m β=0,05;µs=0,25;µm=0,25;vm=0,5.
A423. Calcul de la période et des amplitudes Les conditions de départ sont:
. .. t=t;u(t)=v;u(t)=0 s s b s
L’équation du mouvement pourt∈]t,T[: s
où
u(t)=u(t)+v(t−t) s b s
u(0)−u(t) s T=t+la période du mouvement stickslip. s v b
260
Fig. A43 Types du mouvement stickslip [41] On ne peut pas définir une amplitude de la même façon que pour un mouvement harmonique. Continuons la discussion en traitant les deux phases de façon séparée (Fig. A43).
. Pour le collage:u=v=constb
. . .. Pour le glissement:u≠const, etu=maxsiu(t)=0 m
Une solution approchée pourtmsoit:
t m
t+εtm0m1
oùε<<1est un paramètre. Remplaçons cela dans l’équation du mouvement développée en série de Taylor:
On obtient la solution pourtm1de façon similaire. Remplaçonstm0ettm1dans l’équation detm:
261
π
73 c+c 1 3 32 t=π−m v b Remplaçons le temps dans l’équation de la vitesse:
. . .. . ⎛ ⎞2⎛π⎞2 u(t)=u0(t)+εt u0(t)+u1(t)+O(ε)= −v+ε⎜c⎟ +O(ε)m m0⎜m1m0m0⎟b1 ⎝ ⎠ ⎝2⎠ Comme les extrema de vitesse ne sont pas symétriques par rapport àu=0, l’amplitude est donnée par la définition suivante: . 1⎛ ⎞ A= ⎜v−u(t)⎟v b m 2⎝ ⎠
. En remplaçantu(ton obtient pour) , vb<vb0: m
2 ⎛ ⎞ ⎛ ⎞ ⎛π⎞3πv5v ⎜ ⎟ b b ⎜ ⎟ A= ⎜1−β⎟v+(µ−µ)1−v b s m ⎜ ⎟ 2 8v4v ⎝ ⎠m⎝m⎠ ⎝ ⎠ L’amplitude du déplacement est calculée de façon semblable. Les approximations pour les
valeurs de temps (Fig. A43):
3π1 52 = − t−c1+c3m 2 2 32 π1 52 t= +c−cm+1 3 2 2 32 Les approximations pour les valeurs de déplacement (Fig. A43): 3π = u(t−)−v+2c3+c1m b 4
( )+2−c u tm+=vbc3 14 L’amplitude de déplacementvb<vb0:
2 ⎛ ⎞ π πv⎛v⎞ 1 3b⎜5b⎟ πβ µ µ A=(u(tm)−u(tm))=vb−c=(1−)vb+(s−m)1−⎜ ⎟+ −1 ⎜ ⎟ 2 2 4v4v m⎝m⎠ ⎝ ⎠ Le maximum de la vitesse :
⎧ µ µ 4 (s−m) 3 v=v, si≤(1−(π−3)β) ⎪b0b1 5v8 ⎪ m ⎪ ⎛ ⎞ * v=b⎨ ⎜ ⎟ v−πβ ⎪4m(1 ) ⎜ ⎟ v1+ailleurs m ⎪π 15⎜3⎟ (µ−µ) ⎪ s m ⎩ ⎝4⎠
262
Siµ−µ est petit, le maximum de l’oscillation apparaît pourvb=vb0. Dans le cas de s m
µ−µle maximum apparaît quelque part pour grande, v∈]0,v[A44). Pour (Fig. s m b b0
µ−µpetit, on a: s m
A=A=v= maxb0 vb=vb0
4 v m 5
4βv m 1− =Av,max 3(µ−µ) s m
donc l’amplitude du déplacement et de la vitesse sont les mêmes à la vitesse maximale d’oscillation .
La fréquence de base du mouvement stickslip est la suivante (Fig. A45): 2π2πvb ω= = <1=ωss0 T v t+u(0)−u(t) b s s
Fig. A44 La variation de l’amplitude de l’oscillation en fonction de la vitesse d’excitation
Fig. A45 La variation de la fréquence de l’oscillation en fonction de la vitesse d’excitation