Contrainte de cisaillement (HU)

Traduction anglaise : Shear stress

Mot en chantier

Dernière mise à jour : 07/02/2024

Contrainte mécanique appliquée de manière tangentielle à une face d'un matériau ; comme toute contrainte, elle s'exprime sous la forme du rapport d'une force à une surface (voir figure 1).

$ τ = \frac{F}{A} \qquad (1) $

Avec :

• $ τ $ : contrainte de cisaillement (Pa) ;
• $ F $ : Force tangentielle à la face (N) ;
• $ A $ : Surface de la face (m²) ;

Figure 1 : La contrainte de cisaillement appliquée sur la surface $ A $, parallèle à la force $ F $ est égale à $ τ $.

Contrainte de cisaillement à l'intérieur des fluides

Une contrainte de cisaillement se développe sur chaque particule fluide chaque fois que l'on met un fluide en mouvement sur une surface fixe. Celle-ci est due au gradient de vitesse, également appelé taux de cisaillement, qui s'établit à l'intérieur du fluide à partir du fond (également des parois ou des rives) (voir figure 2). Elle est particulièrement importante à proximité de la fine couche laminaires formée au contact des aspérités du fond (radier de canal, collecteur ou canalisation ou substratum rocheux d’un cours d’eau) et/ou de la surface d’un dépôt de particules solides. Ces contraintes de cisaillement provoquent une déformation angulaire du fluide, dont l'importance dépend de sa viscosité (Oms, 2004).

Figure 2 : Le gradient de vitesse qui se développe à partir du fond (ou éventuellement des parois), supposé immobile, est à l'origine de contraintes de cisaillement qui s'appliquent sur toutes les particules fluides ainsi que sur les particules solides présentes dans l'écoulement ou sur le fond ; c'est ce dernier aspect qui nous intéresse le plus dans cet article.

Von Karman a montré que dans un écoulement turbulent sur fond fixe rugueux de grande largeur (pour ne pas avoir à tenir compte des parois), la vitesse variait de façon logarithmique en fonction de la hauteur d'eau (relation (2)) :

$ \frac{u}{u_*} = \frac{1}{K}ln{\frac{h}{h_0}}+B_s\qquad(2) $

avec :

• $ K $ : constante de Von Karman (sans dimension), dont la valeur est comprise entre 0,36 et 0,42 (normalement 0,41) ;
• $ u $ : vitesse moyenne du fluide à la hauteur $ h $ (m/s) ;
• $ h_0 $ : hauteur de rugosité (hauteur assimilée à celle pour laquelle $ u $ devient nulle) (m) ;
• $ h $ : hauteur dans la veine liquide (m) ;
• $ u_* $ : vitesse de frottement (m/s) ;
• $ B_s $ : constante d'intégration fonction du nombre de Reynolds et dont la valeur est comprise entre 6 et 10 selon le type d'écoulement (sans dimension).

Contrainte de cisaillement sur les parois et lien avec les pertes de charge

Considérons un canal de forme quelconque et un écoulement uniforme à surface libre (figure 3).

Figure 3 : Les forces qui s'appliquent sur le volume de contrôle (en bleu) sont le poids, $ W $, les forces de frottement sur le fond et les parois, $ F $, et les forces de pression (qui s'annulent dans le cas d'un régime uniforme).

Dans ce cas, un volume de contrôle est simplement soumis à son poids et aux forces de frottement (il n'y a aucune accélération et les forces de pression s'annulent du fait de la répartition hydrostatique de cette dernière). La composante du poids parallèle à l'écoulement équilibre donc exactement les forces de frottement. On peut donc écrire la relation (3), en projetant ces forces sur un axe parallèle à l'écoulement :

$ F + W_x = 0 \qquad (3) $

La composante $ W_x $ du poids s’exprime par la relation (4) :

$ W_x = ρ.g.L.S.sin(α) \qquad (4) $

Avec :

• $ ρ $ : masse volumique du fluide (kg/m³),
• $ g $ : accélération gravitationnelle (m/s²),
• $ S $ : section mouillée (m²),
• $ L $ : longueur du volume de contrôle (m),
• $ α $ : angle entre le fond du canal et l'horizontale.

La surface de frottement, $ A $, est égale au produit du périmètre mouillé, $ P $ (longueur du contact entre le fluide et les parois dans une section mouillée) par la longueur de frottement, $ L $, soit :

$ A = P.L \qquad (5) $

Comme l'angle $ α $ est petit, on peut assimiler $ sin(α) $ à la pente $ I $ du canal et obtenir la relation (6) qui permet le calcul de la contrainte de cisaillement sur les parois ou sur le fond, $ τ_w $ (en assainissement on parle de contrainte de cisaillement au radier) :

$ τ_w = \frac{F}{A} = \frac{ρ.g.L.S.sin(α)}{ P.L} = ρ.g.R_h.I\qquad (6) $

avec :

• $ R_h $ : rayon hydraulique ($ R_h = \dfrac{S}{P} $) (m) ;
• $ P $ : périmètre mouillé (m).

Cette relation peut être rapprochée de la formule permettant le calcul des pertes de charge, $ J $ (voir Coefficient de rugosité (HU)), ce qui montre le lien direct entre les deux notions :

$ J = \frac{λ.V^2}{8.g.R_h}\qquad (7) $

Avec :

• $ λ $ : coefficient de pertes de charge (sans dimension) ;
• $ R_h $ : rayon hydraulique de la conduite (m) ;
• $ g $ : accélération de la pesanteur (m/s²) ;
• $ J $ : pertes de charge par unité de longueur (m/m) ;
• $ V $ : vitesse moyenne de l'écoulement (m/s).

Cas du régime uniforme

En régime uniforme la ligne d'eau est parallèle au fond et les pertes de charges par unité de longueur sont égales à la pente du radier : $ I = J $. Il existe donc une relation directe entre la vitesse moyenne de l'écoulement et la contrainte de cisaillement au voisinage des parois (relations (8) ou (9)) :

$ V = \sqrt{\frac{8.τ_w}{λ.ρ}}\qquad (8) $

$ τ_w = \frac{λ.ρ}{8}. V^2\qquad (9) $

Généralisation au cas des écoulements graduellement variés

La relation précédente peuvent être généralisées au cas des écoulements graduellement variés.

La seconde relation des équations de équations de Barré de Saint-Venant, issue du théorème d'Euler, traduit la conservation de la quantité de mouvement. En supposant une pente faible et un fluide incompressible, on obtient la formulation suivante (10) :

$ \frac{\partial Q}{\partial t}+\frac{\partial (Q^2/S)}{\partial x}+g.S.\frac{\partial h}{\partial x}=g.S.(I-J)+k.q.\frac{Q}{S} \quad(10) $

Avec :

• q : perte ou apport latéral de débit (m³/s) ;
• $ I $ : pente du canal (m/m) ;
• $ J $ : pente de la ligne d'énergie ou pertes de charge par unité de longueur (m/m) ;
• $ k $ : coefficient compris entre -1 et +1 ; avec

• $ k $ positif si le débit latéral est entrant (apport de quantité de mouvement) et négatif s'il est sortant (export de quantité de mouvement) ;
• $ k $ proche de 0 si le débit d'apport (ou de fuite) est approximativement orthogonal à $ x $, de façon à ne pas apporter (ou emporter) en même temps une quantité de mouvement, tandis $ k=1 $ suppose que le débit d'apport (ou de fuite) apporte (emporte) avec lui sa quantité de mouvement avec une vitesse égale à celle de l'écoulement principal.

Si on néglige la quantité de mouvement apportée ou emportée par les apports ou les fuites latérales de débit, l'expression se simplifie, et en divisant par S, on obtient la relation classique (11) :

$ \frac{\partial V}{\partial t}+V.\frac{\partial V}{\partial x}+g.\frac{\partial h}{\partial x}=g.(I-J) \quad(11) $

$ g.I $ et $ - g.J $ correspondent respectivement au poids et à la résultante des forces de frottement.

On peut aussi exprimer ces forces de frottement, parallèlement à l’axe x d’écoulement, par le produit de la contrainte de frottement au fond et sur les parois :

$ g.J = τw.P.dx $

soit

$ g.J = τw.\frac{S}{Rh.dx $

Ce qui donne en divisant, comme plus haut, par ρ.S.dx, comme pour l’établissement de l’équation (11) :

$ - τw.P.dx / ρ.S.dx = - τw / ρ. Rh $

On peut donc écrire

- τw / ρ. Rh = - g. J

soit :

$ τ_w = ρ.g.R_h.J\qquad (6) $

Utilisation en transport solide

Une contrainte de cisaillement s'applique également pour toute particule ou objet présent dans l'écoulement et conditionne donc le transport solide. Si l'objet est posé sur le fond il sera mis en mouvement dès que cette contrainte, généralement notée $ τ_w $, deviendra supérieure à une contrainte critique, $ τ_* $, appelée contrainte critique d'entrainement ou seuil de cisaillement.

La contrainte de cisaillement près du fond, $ τ_w $, dépend de la vitesse de l'écoulement et plus précisément de la vitesse près des parois. Elle peut ainsi être rapprochée d'une grandeur particulière, homogène à une vitesse, appelée vitesse de frottement (relation (10)) :

$ τ_w = ρ.u_*^2\qquad(10) $

avec :

• $ τ_w $ : contrainte de cisaillement près du fond (N/m²) ;
• $ ρ $ : masse volumique du fluide (kg/m³).

La valeur de la contrainte critique d'entrainement, $ τ_* $, peut être déterminée en écrivant le bilan des forces qui s'appliquent sur une particule solide posée au fond (voir figure 4).

Figure 4 : Équilibre des forces s'exerçant sur une particule solide posée sur le fond d'un canal soumis à un écoulement ; Source Dufrene et Isenmann (2018).

Cette approche théorique est cependant très difficile à appliquer en pratique, en particulier du fait de la diversité des formes des particules, lesquelles conditionnent la valeur des forces en présence (voir Contrainte critique d’entrainement (HU)) et de la diversité des modes de transport. On se contente donc le plus souvent de modèles empiriques (qui reposent d'ailleurs souvent sur des particules idéalisées, par exemple sphériques), mettant en œuvre l'analyse dimensionnelle, le plus connu étant le diagramme de Shields.

Bibliographie :

• Oms, C. (2004) : Localisation, nature et dynamique de l'interface eau-sédiment en réseau d'assainissement unitaire ; Thèse ENPC ; Partie II. Étude du taux de cisaillement en réseau d’assainissement unitaire ; 40p. ; disponible sur https://pastel.archives-ouvertes.fr/tel-00005725/file/partie_II_Cisaillement.pdf

Pour en savoir plus :

• Dufrene, M., Isenmann, G. (2018) : Transport solide ; cours ENGEES ; disponible sur : https://drive.google.com/file/d/1gO97J6d9Xzg8tx17HPaReEpGVL7sZmMY/view
• https://pastel.hal.science/tel-00005725/file/partie_II_Cisaillement.pdf