Perte de charge de type Borda au franchissement d’un pont; ne peut êter posée que sur un nœud de rivière.
La section d’ouverture du pont est donnée par une courbe tabulée b(Z).
Le champ discharged for headloss computation permet de considérer un contournement potentiel de l'ouvrage par le lit majeur sans obstacle :
si l'option full section est activée la perte de charge est calculée par rapport au débit total au noeud amont.
si l'option river only est activé la perte de charge est calculée par rapport au débit transitant dans le lit mineur seul.
L'option full section est activée par défaut. L'option river only ne doit être utilisé que dans le cas d'un lit majeur ne présentant aucun obstacle.
La loi de perte de charge dépend des conditions d'écoulement au pas de temps considéré, fonction de la position de la cote critique z c z_cz c par rapport à la cote de voute de l'ouvrage z v z_vz v et de la cote de la ligne d'eau en aval de l'ouvrage z 2 z_2z 2 .
Le régime d'écoulement est déterminé comme suit:
si z c < z v z_c < z_vz c < z v et z 2 < z c z_2 < z_cz 2 < z c : régime déversoir dénoyé (DD),
si z c > z v z_c > z_vz c > z v et z 2 < z v z_2 < z_vz 2 < z v : régime vanne dénoyée (VD),
dans les autres cas : régime noyé.
On détermine la hauteur critique sous le pont par :
v b g h c = 1 \frac{v_b} {\sqrt{gh_c}} = 1
g h c v b = 1
avec
h c h_ch c : hauteur critique
v b v_bv b : vitesse sous le pont pour la hauteur critique, prise égale à v b = Q B m o y h c v_b = \frac{Q}{B_{moy}h_c}v b = B m o y h c Q
Q QQ : débit passant sous le pont,
B m o y = S b 0 z v − z r B_{moy} = \frac{S_{b0}}{z_v-z_r}B m o y = z v − z r S b 0 , largeur moyenne sous le pont,
z r z_rz r : cote de radier du pont,
z v z_vz v : cote de la voûte du pont,
S b 0 S_{b0}S b 0 : section sous le pont.
Il s'en suit que la cote critique z c z_cz c sous le pont s'exprime par :
z c = z r + h c = z r + ( Q B m o y g ) 2 / 3 z_c = z_r + h_c = z_r + \left(\frac {Q} {B_{moy}\sqrt{g}}\right)^{2/3}
z c = z r + h c = z r + ( B m o y g Q ) 2 / 3
Q QQ étant le débit sous le pont.
Soit z 2 z_2z 2 : niveau d'eau à l'aval du pont,
E 1 − E 2 = 1 2 g ( u b − u 2 ) 2 ( 2.1 ) E_1 - E_2 = \frac{1}{2g}(u_b-u_2)^2 \ \ \ \ \ \ \ \ \ \ \ (2.1)
E 1 − E 2 = 2 g 1 ( u b − u 2 ) 2 ( 2 . 1 )
où:
E 1 = z 1 + u 1 2 g E_1 = z_1 + \frac{u_1}{2g}E 1 = z 1 + 2 g u 1 : énergie spécifique à l'amont de la singularité
E 2 = z 2 + u 2 2 g E_2 = z_2 + \frac{u_2}{2g}E 2 = z 2 + 2 g u 2 : énergie spécifique à l'amont de la singularité
z 1 z_1z 1 : niveau d'eau à l'amont du pont,
z 2 z_2z 2 : niveau d'eau à l'aval du pont,
u b = Q m i n c c s b u_b = \frac{Q_{min}}{c_c s_b}u b = c c s b Q m i n : vitesse d'écoulement dans la section rétrécie au droit du bridge
s b s_bs b : section mouillée du bridge
c c c_cc c : coefficient de contraction de l'écoulement au passage dans le bridge
u 2 = Q m i n s 2 u_2 = \frac{Q_{min}}{s_2}u 2 = s 2 Q m i n : vitesse d'écoulement dans la section aval du lit mineur
L'équation (2.1) est réécrite comme suit:
z 1 − z 2 = B Q m i n 2 2 g s b 2 ( 2.2 ) z_1 - z_2 = \frac{B Q_{min}^2} {2gs_b^2} \ \ \ \ \ \ \ \ \ \ \ (2.2)
z 1 − z 2 = 2 g s b 2 B Q m i n 2 ( 2 . 2 )
où
B = ( 1 c c − s b s 2 ) 2 + β 2 ( s b s 2 ) 2 − β 1 ( s b s 1 ) 2 ( 2.3 ) B = {(\frac{1}{c_c} - \frac{s_b}{s_2})}^2 + \beta_2 {(\frac{s_b}{s_2})}^2 - \beta_1{(\frac{s_b}{s_1})}^2 \ \ \ \ \ \ \ \ \ \ \ (2.3)
B = ( c c 1 − s 2 s b ) 2 + β 2 ( s 2 s b ) 2 − β 1 ( s 1 s b ) 2 ( 2 . 3 )
β 1 \beta_1β 1 et β 2 \beta_2β 2 sont des coefficients correcteurs en cas de vitesse élevée, précisés plus bas.
Dans le cas d'une modélisation 1D le débit nodal est le débit total. Il faut alors corriger l'équation (2.2) si la singularité n'intercepte que le débit du lit mineur:
z 1 − z 2 = α q 2 B Q 2 2 g s b 2 ( 2.4 ) z_1 - z_2 = \frac{\alpha_q^2 B Q^2} {2gs_b^2} \ \ \ \ \ \ \ \ \ \ \ (2.4)
z 1 − z 2 = 2 g s b 2 α q 2 B Q 2 ( 2 . 4 )
où α q = Q m i n / Q \alpha_q = Q_{min}/Qα q = Q m i n / Q
Q m i n Q_{min}Q m i n : débit dans le lit mineur à l'amont de la singularité,
Q QQ : débit total.
La loi d'écoulement en régime noyé est résolue avec l'équation (2.4) réécrite comme suit:
Q = β ( z 1 − z 2 ) 0.5 . 2 g . s b = c n ( z 1 − z 2 ) 0.5 ( 2.5 ) Q \ =\ \beta(z_1-z_2)^{0.5}.\sqrt{2g}.s_b \ =\ c_n(z_1-z_2)^{0.5}\ \ \ \ \ \ \ \ \ \ \ (2.5)
Q = β ( z 1 − z 2 ) 0 . 5 . 2 g . s b = c n ( z 1 − z 2 ) 0 . 5 ( 2 . 5 )
avec β = 1 α q B \beta=\frac{1}{\alpha_q \sqrt B}β = α q B 1 , c n = β 2 g . s b c_n = \beta \sqrt{2g}.s_bc n = β 2 g . s b
Soient:
q qq le débit dans le lit mineur non influencé (aval),
s 1 s_1s 1 la section mouillée à l'amont du pont,
b 1 b_1b 1 la largeur au miroir à l'amont du pont,
u 1 u_1u 1 la vitesse à l'amont du pont: u 1 = q / s 1 u_1 = q/s_1u 1 = q / s 1
On définit le nombre de Froude par F r 1 = u 1 . b 1 g s 1 F_{r1} = u_1.\sqrt{ \frac {b_1} {gs_1} }F r 1 = u 1 . g s 1 b 1 qu'on borne à 0.7
On exprime le terme d'énergie cinétique e c = u 1 2 2 g e_c = \frac {u_1^2 }{2g}e c = 2 g u 1 2 :
e c = F r 1 2 s 1 2 b 1 e_c = F_{r1}^2 \frac{s_1}{2b_1}
e c = F r 1 2 2 b 1 s 1
Le coefficient β 1 \beta_1β 1 est exprimé par:
β 1 = e c . 2 g ( s 1 q ) 2 \beta_1 = {e_c}. 2g\left( {\frac{s_1}{q}}\right)^2
β 1 = e c . 2 g ( q s 1 ) 2
Si F r 1 < 0.7 F_{r1} < 0.7F r 1 < 0 . 7 alors β 1 = 1 \beta_1 = 1β 1 = 1
Le coefficient β 2 \beta_2β 2 est exprimé de la même manière à partir des paramètres en aval du pont.
Il est donné en fonction du rapport s b / s 1 s_b/s_1s b / s 1 selon le tableau ci-dessous, où:
s 1 s_1s 1 est la section mouillée correspondant à la cote amont z 1 z_1z 1 ,
s b s_bs b est la section mouillée sous le pont sous la même cote z 1 z_1z 1 .
Tableau c c c_cc c
s b / s 1 s_b/s_1s b / s 1
0
0.6
0.8
0.9
0.95
1
c c c_cc c
0.6
0.6
0.75
0.82
0.87
0.999
L'équation en régime déversoir dénoyé est écrite comme suit:
Q p o n t = 2 3 μ d 2 g B m o y ( E 1 − Z r ) 3 / 2 (2.7) Q_{pont} = \frac{2}{3} \mu_d \sqrt{2g} B_{moy}(E_1-Z_r)^{3/2} \tag{2.7}
Q p o n t = 3 2 μ d 2 g B m o y ( E 1 − Z r ) 3 / 2 ( 2 . 7 )
avec
Q p o n t Q_{pont}Q p o n t : débit sous le pont,
E 1 E_1E 1 : énergie amont (voir plus bas),
Z r Z_rZ r : cote de radier du pont,
B m o y = S b 0 z v − z r B_{moy} = \frac{S_{b0}}{z_v-z_r}B m o y = z v − z r S b 0 , largeur moyenne sous le pont,
S b 0 S_{b0}S b 0 : section pleine sous le pont,
z v z_vz v : cote de la voûte du pont
μ d = 0.56 \mu_d = 0.56μ d = 0 . 5 6
Dans le cas d'une modélisation 1D le débit nodal est le débit total. Il faut alors corriger l'équation (2.7) si la singularité n'intercepte que le débit du lit mineur:
Q = 1 α q Q p o n t Q = \frac{1}{\alpha_q} Q_{pont}
Q = α q 1 Q p o n t
où α q = Q m i n / Q \alpha_q = Q_{min}/Qα q = Q m i n / Q
Q m i n Q_{min}Q m i n : débit dans le lit mineur à l'amont de la singularité,
Q QQ : débit total.
L'équation en régime vanne dénoyée est écrite comme suit:
Q p o n t = μ d 2 g B m o y ( Z v − Z r ) ( E 1 − Z r ) 1 / 2 ( 2.8 ) Q_{pont} = \mu_d \sqrt{2g} B_{moy}(Z_v-Z_r)(E_1-Z_r)^{1/2} \ \ \ \ \ \ \ \ \ \ \ (2.8)
Q p o n t = μ d 2 g B m o y ( Z v − Z r ) ( E 1 − Z r ) 1 / 2 ( 2 . 8 )
avec
Q p o n t Q_{pont}Q p o n t : débit sous le pont,
E 1 E_1E 1 : énergie amont (voir plus bas),
Z r Z_rZ r : cote de radier du pont,
B m o y = S b 0 z v − z r B_{moy} = \frac{S_{b0}}{z_v-z_r}B m o y = z v − z r S b 0 , largeur moyenne sous le pont,
S b 0 S_{b0}S b 0 : section pleine sous le pont,
z v z_vz v : cote de la voûte du pont
μ d = 0.56 \mu_d = 0.56μ d = 0 . 5 6
De même que pour le régime dénoyé il faut corriger l'équation (2.8) si la singularité n'intercepte que le débit du lit mineur:
Q = 1 α q Q p o n t Q = \frac{1}{\alpha_q} Q_{pont}
Q = α q 1 Q p o n t
E 1 = Z 1 + u 1 2 2 g E_1 = Z_1 + \frac{u_1^2}{2g}
E 1 = Z 1 + 2 g u 1 2
Elle est réexprimée en fonction du nombre de Froude.
Pour des raisons de stabilité numérique le nombre de Froude F r 1 F_{r1}F r 1 est borné par F r m a x = 0.7 F_{rmax} = 0.7F r m a x = 0 . 7
E 1 = Z 1 + s 1 2 b m 1 F r 1 2 E_1 = Z_1 + \frac{s_1}{2 b_{m1}} F_{r1}^2
E 1 = Z 1 + 2 b m 1 s 1 F r 1 2
s 1 s_1s 1 : section mouillée amont,
b m 1 b_{m1}b m 1 : largeur au miroir amont,
F r 1 F_{r1}F r 1 : nombre de Froude à l'amont borné par F r m a x = 0.7 F_{rmax} = 0.7F r m a x = 0 . 7
Si la cote amont z 1 z_1z 1 dépasse la cote de la chaussée z r o a d z_{road}z r o a d on ajoute un débit de surverse sous la forme d'une loi de déversoir:
en régime dénoyé si ( z 2 − z r o a d ) < 2 3 ( z 1 − z r o a d ) (z_2-z_{road}) < \frac {2}3{}(z_1-z_{road})( z 2 − z r o a d ) < 3 2 ( z 1 − z r o a d ) :
Q r o a d = 2 3 μ d 2 g . b r o a d ( z 1 − z r o a d ) 3 / 2 a v e c μ d = 0.58 ( 2.12 ) Q_{road} = \frac {2}{3} \mu_d \sqrt{2g}.b_{road} (z_1 - z_{road})^{3/2}\ \ \ \ \ \ \ \ \ avec\ \ \ \ \ \ \ \ \mu_d = 0.58\ \ \ \ \ \ \ \ (2.12)
Q r o a d = 3 2 μ d 2 g . b r o a d ( z 1 − z r o a d ) 3 / 2 a v e c μ d = 0 . 5 8 ( 2 . 1 2 )
en régime noyé dans le cas contraire :
Q r o a d = 2.6 2 3 μ d 2 g . b r o a d ( z 1 − z 2 ) 1 / 2 ( z 2 − z r o a d ) ( 2.13 ) Q_{road} = 2.6\ \frac {2}{3} \mu_d \sqrt{2g}.b_{road} (z_1 - z_2)^{1/2}(z_2 - z_{road})\ \ \ \ \ \ \ \ \ \ (2.13)
Q r o a d = 2 . 6 3 2 μ d 2 g . b r o a d ( z 1 − z 2 ) 1 / 2 ( z 2 − z r o a d ) ( 2 . 1 3 )
Il faut corriger ces équations si la singularité n'intercepte que le débit du lit mineur:
Q = 1 α q Q p o n t Q = \frac{1}{\alpha_q} Q_{pont}
Q = α q 1 Q p o n t
Cet objet est de type singularité . Il ne peut être posé que sur un container de type noeud de rivière ( River node ).