,

Une nouvelle méthode de calcul des reseaux maillés

A new mesh network calcl/lation method
André Divenot
Compagnie Générale des Eaux, Paris

Introduction f(x) = 0 (1 )

Au cours des quinze dernières années, la méthodo- la solution est recherchée en appliquant à l'inconnue
logie du calcul des réseaux maillés de distribution d'eau des corrections successives, de valeurs absolues décrois-
en état permanent a beaucoup évolué, en liaison avec santes, à partir d'une valeur approchée x o . A l'itération i,
les possibilités nouvelles offertes par les ordinateurs. une approximation meilleure est obtenue ainsi:
Vers le milieu des années soixante dix, la Compagnie
Générale des Eaux a souhaité disposer de nouveaux (2a)
programmes, car ceux qui étaient en exploitation pré- avec
sentaient quelques restrictions' qui n'étaient plus
admises. (2b)
11 est rapidement apparu que, pour les raisons qui
seront exposées ci-après, aucune des méthodes connues
ce qui peut aussi s'écrire:
ne pouvait donner pleinement satisfaction. Les recher-
ches entreprises en 1975 ont débouché sur une méthode
originale de génération et de traitement du système
d'équations des mailles. Le premier programme à l'uti-
liser est le produit-programme PROGRES (pROG pour A l'itération i + l, xi est pris comme nouvelle valeur
programme et RES pour réseaux) qui tient plus que ses approchée de x. Le processus est répété jusqu'à l'itéra-
promesses depuis sa mise en exploitation, au cours de tion p où un critère de précision est satisfait. Le critère
l'été] 976. peut être appliqué à la valeur absolue de DX p ou à la
Afin- de bien situer le problème, il est nécessaire de valeur absolue de f(x P _ I ) ou encore aux deux. Alors,
rappeler tout d'abord les principes de la technique de il est admis que x = x p'
Newton-Raphson et des méthodes de résolution Dans le cas d'un système de 11 équations non linéaires
publiées. à 11 inconnues, et en notation matricielle, soit le vec-
teur X des 11 variables inconnues ou variables d'ité-
ration:
La méthode de Newton-Raphson (3)
Calculer un réseau en état permanent revient à ré- 11 s'agit de résoudre le système:
soudre un système de n équations algébriques non
linéaires à f1 inconnues, 11 pouvant atteindre des valeurs
très élevées, de l'ordre du millier. Parmi les techniques
itératives utilisables, c'est celle de Newton-Raphson La solution est recherchée en ajoutant au vecteur
qui a la faveur des spécialistes du calcul des réseaux des variables inconnues des vecteurs de correction
maillés de distribution d'eau. successifs, de normes décroissantes, à partir d'une
Dans le cas simple d'une seule équation à une seule valeur approchée X o. A l'itération i, une approxima-
inconnue,x: tion meilleure est obtenue ainsi:

LA HOUILLE BLANCHE/N° 6-1980

Article published by SHF and available at http://www.shf-lhb.org or http://dx.doi.org/10.1051/lhb/1980034

366 LA HOUILLE BLANCHE/N° 6-1980

(5a) où Qij est le débit du tronçon ij et Ci la consommation

avec prélevée au nœud i. Qij est nul si le tronçon ij n'existe
pas. Tout débit s'éloignant du nœud est positif.
(5b) 2) La somme des pertes de charge des tronçons
appartenant à toute maille m est nulle. C'est la loi de
(5c) Kirchoff aux mailles, ou loi de conservation de
l'énergie:
T
ce qui peut aussi s'écrire:
lm L h ij = 0 (m = l, ... ,M) (9)
(5d) ij€m
ij= 1
Ji __ 1 est la matrice jacobienne correspondant au
vecteur Xi _ l' Afin d'éviter les indices multiples, et où h ij est la perte de charge du tronçon ij, positive si
maintcnant que le processus est bien établi, la réfé- le sens du débit du tronçon coïncide avec le sens de la
rence à l'itération i est abandonnée. L'équation 5d maille, lequel est arbitraire mais doit être parfaitement
s'écrit alors: défini. La relation fonctionnelle qui lie la perte de charge
et le débit est, sous une forme générale:
JI:!.X =- F(X) (5e)
h ij = Hi -- H j = R ij Qij IQijIO:-1 (lOa)
et J peut être défini par l'élément de la ligne i et de la
colonne j : ou, réciproquement
1 1--0.:
J=aF(X)=[a / i (X 1 , ... ,X n )] (5f)
Qij = R ij (Hi H j ) IHi HIJ 0: (lOb)
ax aXj 0:

où Hi et H j sont les niveaux piézométriques des nœuds i

La similitude entre les équations 2a, 2b, 2c d'une
et j. R ij est la résistance hydraulique du tronçon ij.
part et les équations 5a, 5c, 5d d'autre part n'est pas
L'exposant 0: est habituellement compris entre 1,8 et 2.
fortuite. Comme dans le cas d'une seule équation, le
Le système d'équations non linéaires qu'il faut
calcul est poursuivi jusqu'à ce qu'un ou plusieurs cri-
résoudre pour équilibrer le réseau peut avoir pour dimen-
tères de précision soient satisfaits. La différence est
sions le nombre de nœuds N, ou le nombre de mailles
que le test d'un vecteur peut porter sur la valeur ab-
M, ou encore le nombre de tronçons T. L'évaluation
solue de chacun de ses éléments ou, globalement, sur
initiale porte, dans le premier cas, sur les niveaux piézo-
sa norme. métriques inconnus et, dans le deuxième cas, sur les
débits des tronçons. Dans le troisième cas, aucune éva-
luation initiale n'est nécessaire. Toutes les méthodes
Les bases du calcul des réseaux maillés en état qui vont maintenant être présentées s'appuient sur l'un
permanent ou l'autre des trois systèmes d'équations.

Afin d'alléger l'exposé, les pseudo-mailles ou mailles

ouvertes sont ignorées dans cette première partie, ce
qui revient à dire que, par convention, un réseau ne Système d'équations des noeuds
comprend qu'un seul nœud de niveau piézométrique
déterminé, appelé racine. La question sera reprise en La loi de Kirchoff aux mailles est constamment
détail plus loin. satisfaite et l'objectif est de satisfaire la continuité des
Pour un réseau de T tronçons, N nœuds et M mailles, débits aux nœuds, en corrigeant les niveaux piézomé-
on peut alors écrire: triques inconnus. Le rapprochement entre les équa-
tions 8 et lOb donne:
T=M+N--- 1 (6) yN _1.. 1-0:

Comme, dans les réseaux réels, T est couramment fi = '-' R ij 0: (Hi - H j ) IHi - Hjl-O:- + Ci =0
supérieur à N de 40 à 50 %, il vient, en première appro- joFi
ximation : j=l (lI)

N::::::2M (7) Le système est alors obtenu en appliquant les équa-

T ::::::3M tions 5d et 5f :

Un réseau maillé en état permanent est équilibré

lorsque deux ensembles de conditions sont satisfaits :
ail ail
- -- 8H 1 - Il
1) En tout nœud i autre que la racine r, la somme aH l aHN
des débits est nulle. C'est la loi de Kirchoff aux nœuds (12)
ou loi de continuité du débit: '
N alN N al
fi = L Qi! + Ci 0 Ci =1= r ;i = l, ... ,N) (8)
-
aH l
" . --
aH
8HN ·-IN
=1 N
joFi
j=l avec
 A. DIVENOT 367

l-a
a appliquer immédiatement la correction de débit de
(i =1= j) chaque maille, avant de passer aux mailles suivantes,
(13) faire en sorte que la partie commune à deux mailles
et soit de résistance hydraulique minimale.
N
Y
..... -a/i (14)
On le verra, les deux propositions sont toujours per-
tinentes. La deuxième a été formulée à nouveau par
j*i aHj
Voyles et Wilke [16] en 1962.
j=l
Epp et Fowler [6] ont noté que la matrice est symé-
La matrice jacobienne est définie positive symétrique trique et ont décrit un algorithme d'arrangement en
et elle est de rang (N - 1) puisque le nœud de niveau matrice-bande. La contribution de Kesavan et Chan-
piézométrique déterminé n'y figure pas. drashekar [9] porte sur la génération des systèmes
En 1936, Cross [3] établissait les relations de base d'équations (des nœuds ou des mailles) et celle de
pour le traitement de chaque nœud considéré indivi- Jeppson et Davis [7] sur l'inclusion des réducteurs de
duellement, ce qui conduisait à des convergences parti- pression (également avec les systèmes des nœuds et
culièrement difficiles. En 1963, Martin et Peters [14] des tronçons).
étaient les premiers à utiliser la technique de Newton-
Raphson. Plus tard, Shamir et Howard [15] ont montré
que, tout en gardant le même nombre d'équations et Système d'équations des tronçons
d'inconnues, il est possible de substituer à une partie
de ces dernières des consommations aux nœuds ou des Les T équations à T inconnues sont obtenues en
résistances de tronçons. Liu [13] essaie de ne conserver écrivant (N - 1) équations de continuité des débits
que la diagonale principale de la matrice jacobienne. aux nœuds et M équations de mailles. Conformément
Donachie [4] insiste sur la possibilité d'ajuster les résis- à l'équation 6, le nombre total d'équations est alors
tances pour mieux caler le modèle. bien égal à T. Les variables d'itération sont les débits
Chandrashekar et Stewart [2] tirent parti du fait que des tronçons.
la matrice jacobienne est particulièrement clairsemée Les (N - 1) équations aux nœuds sont directement
et ne calculent pas, explicitement, son inverse. Afin déduites de l'équation 8. Ce sont des équations linéaires.
d'éviter également les inversions répétées, Lam et Wolla Les M équations de mailles sont linéarisées en réécrivant
[10, Il] utilisent une inverse approchée, retouchée à l'équation 1Da sous la forme:
chaque itération. Lemieux [12] propose une autre
modification de l'algorithme de Newton-Raphson, h ij = RijQijlQijla-1 = KijQij (16)
basée sur les propriétés des matrices-bandes symétriques, et les K ij sont recalculés à chaque itération avec les
définies positives. L'arrangement en matrice-bande, nouveaux débits.
dans le but de réduire l'encombrement, est également Aucune évaluation initiale n'est nécessaire. Plus
pratiqué par d'autres auteurs [7, 14]. exactement, Wood et Charles [17] qui ont décrit cette
technique obtiennent de bons résultats en calculant
les Kij initiaux de l'équation 16 pour un débit égal à
Système d'équations des mailles l'unite. Jeppson et Tavallaee [8] ont ensuite explicité
les modifications nécessaires pour inclure les pompes
Les lignes qui suivent exposent seulement le principe et réservoirs.
de la méthode. D'autres détails seront donnés plus loin.
La loi de Kirchoff aux nœuds est satisfaite dès la répar-
tition initiale des débits et le reste constamment. L'ob- Limites des méthodes
jectif est de satisfaire la loi de conservation de l'énergie
en corrigeant les débits des mailles. Chaque correction Ce qui précède n'est qu'un aperçu des contributions,
ôQm est appliquée à tous les tronçons de la maille m. la plupart très intéressantes, publiées au cours des der-
Le système résulte de l'application des équations 5d nières années. D'une façon générale, il est fait appel, et

l
et 5f aux équations 9 et IDa: avec succès, à tout l'arsenal mathématique dans le but de
diminuer l'encombrement de la matrice et la durée des
ail
aQI
ail
--
aQM
r ÔQI
-t,
calculs.
Toutefois, des problèmes subsistent. Plusieurs auteurs
(15) signalent que la méthode de Newton-Raphson ne con-
verge pas lorsque les valeurs initiales sont mal choisies
al alM [4, 10, Il, 12, 15, 17]. C'est probablement pour cette
-M
- -- ÔQM -lM J raison que le système des nœuds a la faveur. Il est en
aQI aQM
effet plus facile de programmer la génération de niveaux
ta 1936, Cross [3] présentait l'algorithme de calcul piézométriques plausibles que de programmer la géné-
de chaque maille prise individuellement, toutes les ration de débits de tronçons plausibles.
corrections étant appliquées seulement à la fin de l'ité- Malheureusement, la méthode des nœuds présente
ration. En 1947, Dubin [5] remarquait que la conver- deux difficultés. D'une part, les tronçons de faible perte
gence est souvent médiocre et énonçait deux principes de charge entraînent des dérivées partielles très grandes.
pour l'améliorer: D'autre part, les oscillations près de la solution sont
368 LA HOUILLE BLANCHE/N° 6-1980

fréquentes. Cependant, il convient de noter que bien des métrique déterminé n'est pas obligatoirement un noeud
auteurs annoncent le plein succès, grâce aux diverses de niveau piézométrique constant quel que soit le débit.
améliorations apportées à la méthode de Newton - II suffit que la relation entre le débit et le niveau piézo-
Raphson. métrique soit connue. Ce peut être, par exemple, l'en-
Par ailleurs, les méthodes matricielles sont théo- semble constitué par un réservoir et une pompe de
riquement - moins bien adaptées que les méthodes refoulement.
maille par maille ou nœud par nœud au traitement des Le nombre total de mailles d'un réseau, M, s'écrit
clapets. Ces organes de coupure asymétriques ont long- finalement:
temps été ignorés ou pris en considération seulement
par le biais du réducteur de pression [7] ou de la station M = (T - N + 1) + (Z - 1) = T +Z - N (17)
de pompage [Il, 15]. Chandrashekar [1] a récemment
où M est le nombre de mailles indépendantes strictement
fait le point en la matière, dans le cadre du système
nécessaire et suffisant pour équilibrer le réseau. La rela-
d'équations des noeuds.
tion est exacte pour les réseaux palmés, maillés ou
A chaque itération, toute méthode matricielle con-
mixtes. Ainsi, si T = 263, N = 264 et Z = 1, le réseau
sidère globalement un système d'équations qui corres-
est palmé. Si T = 180, N = 145 et Z = 9, le réseau
pond à un état bien précis des clapets: chacun d'entre
comprend 44 mailles dont 8 ouvertes et il est moyen-
eux est soit ouvert soit fermé. Application faite des cor-
nement maillé car le ratio (T + Z)/N est égal à 1,3. Pour
rections, il est facile de tester la validité des hypothèses
les réseaux de distribution d'eau fortement maillés, ce
mais s'il en résulte le moindre changement d'état, il n'est
ratio peut dépasser 1,5. Il est évidemment fonction de la
généralement pas possible d'en déduire à coup sûr le sys-
proportion de noeuds de niveaux piézométriques déter-
tème d'équations approprié pour l'itération suivante,
minés, mais aussi de la proportion de noeuds poin ts de
faute d'algorithme de prédiction des interactions en
jonction d'au moins trois tronçons, ainsi que de l'impor-
chaîne : l'ouverture ou la fermeture d'un clapet peut
tance relative des antennes.
entraîner une cascade imprévisible d'autres modifi-
La répartition entre mailles ouvertes et mailles
cations. La bonne fm du calcul ne peut donc pas être
fermées peut être modifiée par un type d'appareil qui,
garantie [1, Il, 15] bien que, dans la pratique, il semble
tel Janus, a deux visages. C'est le réducteur ou stabili-
qu'elle soit atteinte sans difficultés majeures.
sateur de pression qui est un noeud de consommation
Toutes ces raisons ont motivé la recherche d'une
variable à l'amont et un nœud de niveau piézométrique
nouvelle méthode. Entre les divers systèmes d'équations,
déterminé à l'aval, ce qui détruit la symétrie de la ma-
la préférence allait aux mailles: en appliquant à un gros
trice jacobienne [1,7]. Si par exemple une zone basse est
réseau les approximations grossières données précédem-
alimentée par l'intermédiaire de n réducteurs de pres-
ment (équation 7), il est logique de choisir 400 équa-
sion, alors (n - 1) mailles ouvertes supplémentaires se
tions de mailles plutôt que 800 équations de noeuds ou
substituent à (n - 1) mailles fermées.
1 200 équations de tronçons. Mais, de proche en proche,
II est maintenant possible de généraliser l'équation 9
une propriété beaucoup plus intéressante a été mise en
sous la forme:
évidence. Tandis que les systèmes des tronçons ou des
nœuds sont lourdement figés, un système d'équations
de mailles peut être sélectionné et, surtout, dynamisé. fm = Lh t - L. hp - (HA - H B ) =0
C'est ce qui va maintenant être exposé. tEm pEm

(m = 1, ... , M) (18)
Compléments sur le système d'équations des où les deux sommes s'appliquent aux éléments appar-
mailles tenant à la maille m. La perte de charge ht dans le
tronçon t peut s'écrire:
L'équation 6 donne le nombre M de mailles fermées (19)
dans un réseau connexe de N noeuds et T tronçons ne
comprenant qu'un seul noeud de niveau piézométrique sans référence à quelque formule de perte de charge que
déterminé. Chacune des (T - N + 1) mailles est un ce soit. La résistance hydraulique R t est fonction des
cheminement fermé de tronçons tous différents ; cha- caractéristiques physiques de la canalisation (diamètre,
que tronçon appartient à au moins une maille; chaque rugosité, longueur) et l'exposant Œ est habituellement
maille a au moins un tronçon qui lui est propre. Les compris entre 1,8 et 2. Mais, en généralisant, R t peut
antennes, qui sont des cheminements ouverts issus du aussi être fonction de la vitesse de circulation du fluide
réseau et aboutissant à des noeuds de niveaux piézo- ou résulter de la combinaison de plusieurs éléments, en
métriques indéterminés, sont naturellement ignorées au série ou en parallèle, et l'exposant Œ peut varier d'un
cours de la recherche. tronçon à l'autre. Il est parfaitement possible, par
Si, parmi les N noeuds, il en est Z, et non plus un seul, exemple, de combiner des pertes de charge singulières
dont les niveaux piézométriques sont déterminés, le - fonction du carré du débit - avec d'autres pertes de
nombre de mailles fermées est inchangé mais il faut alors charge calculées avec la formule de Colebrook - White.
ajouter (Z - 1) pseudo-mailles ou mailles ouvertes. Dans l'équation 19, le débit Qt est positif lorsqu'il est
Chaque maille ouverte joint un noeud de niveau piézo- orienté dans le sens de la maille.
métrique déterminé à un autre noeud de niveau piézomé- La deuxième somme de l'équation 18 est étendue aux
trique déterminé ; chacun des Z noeuds appartient à au stations de refoulement ou de surpression de la maille.
moins une maille ouverte. Un noeud de niveau piézo- La hauteur de pompage hp s'exprime habituellement
 A. DIVENOT 369

- ce n'est pas une obligation - par un polynôme :

B C
hp = bo + bIQp + b2Q~ + ... + bnQ; (20)
Le débit de pompage Qp et la hauteur de pompage h p
sont normalement positifs lorsque le sens du pompage
coïncide avec le sens de la maille.
EnfIn, toujours dans l'équation 18, les niveaux piézo-
métriques HA et Hs sont ceux des nœuds d'origine (A)
et d'extrémité (B) de la maille. Les deux nœuds ne sont A
distincts que si la maille est ouverte.
figure 1

Figure 1 - Mini-réseau exemple

Calcul de deux mailles liées

Le mini-réseau de la fIgure 1 sera pris en exemple. Les

Donnée Tronçon 1 Tronçon 2 Tronçon 3
nœuds A, B et C sont des réservoirs de niveaux piézomé-
triques déterminés. Le point de jonction des trois tron-
çons est un nœud ordinaire sans consommation. Appli- 3,75.10- 5 3,75.10- 6 1,25.10- 6
rési stance
cation de l'équation 17, le nombre strictement nécessaire
et suffIsant demaillesindépendantes.M.estégalà2.et débit initial 200 100 100
les deux mailles sont ouvertes. Il est supposé que le débit exa~t 400 200 200
choix s'est porté sur les mailles formées par les tronçons h exact 6 0,15 0,05
1 et 2 d'une part, et par les tronçons 1 et 3 d'autre part,
ce qui sera noté 1 + 2 et 1 + 3. Il est important de
remarquer que l'équilibrage du réseau nécessite la résolu-
tion d'un système de deux équations indépendantes mais Tableau 1 - Données pour le calcul du réseau de la figure 1.
que le choix est possible entre le système retenu, mailles Unités: I/s et m
1 + 2 et 1 + 3 (le tronçon 1 est la partie commune aux
deux mailles) ou les mailles 1 + 2 et 2 + 3 (le tronçon 2 avec
est la partie commune) ou enfin les mailles 1 + 3 et
2 + 3 (le tronçon 3 est la partie commune).
Les flèches de la fIgure 1 indiquent le sens des mailles, afI2 = IR Q +2R 2Q2
aQI2 - 1 1
qui est également celui des débits positifs. Les débits
seront d'ailleurs constamment positifs, ce qui évitera = 7,5 x 10- 5 (QI + 0,1 Q2)
toute référence aux valeurs absolues dans l'équation 19
et dans l'expression des dérivées. Les unités sont le Ils (23)
pour les débits et le m pour les pertes de charge. Pour af12 = afl3 = 7,5 x 10- 5 QI
simplifier, l'exposant ex de l'équation 19 est pris égal à 2. aQl3 aQ12
Les résistances, les débits initiaux, les débits exacts et
les pertes de charge à la solution sont fournis dans le
tableau 1. Pour donner un simple ordre de grandeur,
af13 = 7,5 x 10- 5 (QI +~ Q3)
aQI3 3
les trois tronçons pourraient avoir la même longueur
de 1 km et des diamètres respectifs de 500, 800 et Un calcul est terminé lorsque la valeur absolue de
1000 mm. chaque élément du vecteur F est inférieure à 0,5 mm et
Ci-après, dans l'application des équations 15, 18 et que la valeur absolue de chaque élément du vecteur ~Q
19, le double indice se rapporte aux deux tronçons est inférieure à 0,05 Ils. Les quantités des tableaux 2 et
concernés. Ainsi fi} est ia fermeture en cote de la maille 3 ont été arrondies individuellement.
ij, c'est-à-dire de la maille i + j. Le système choisi est La partie supérieure du tableau 2 montre le dérou-
le suivant: lement du calcul selon la méthode de Newton - Raphson.
Les tests sont satisfaits à l'itération 5. La convergence est
f 12 = 3,75 x 10- 5 Qî \ rapide et la précision fmale excellente. Le contraste est
+ 3,75 x 10- 6 Q~ - 6,15 =0 (
frappant avec la partie inférieure du même tableau 2,
où le réseau est calculé maille par maille avec application
(21)
f 13 = 3,75 ~ 10- 5
Qî \ immédiate de la correction. Seuls les deux éléments dia-
gonaux de la matrice jacobienne de l'équation 22 sont
+ 1,25 x 10- 6 Q~ - 6,05 =0 J pris en considération, l'un après l'autre. La maille 1 + 2
génère la correction oQ 12 qui est appliquée aux débits
af12 af12 QI et Q2 avant de traiter la maille 1 + 3. Les tests de
-- --
aQI2 aQl3 fIn de calcul ne sont satisfaits qu'à l'itération 101 et les
(22) débits Q2 et Q 3 ne sont équilibrés qu'à 0,2 Ils près.
afl3 af13 Avec un test moins sévère, 1 Ils par exemple, le calcul
---- aurait cessé à l'itération 38, mais avec un déséquilibre
aQI2 aQl3
de l'ordre de 141/s.
370 LA HOUILLE BLANCHE/N° 6-1980

li
itération f f 00 12
12 13 0°13 °1 °2 °3

a 200,00 100,00 100,00

1 -4,6125 -4,5375 150,00 150 ,00 500,00 250,00 250 ,00
2 3,4594 3,4031 -45,00 -45,00 410,00 205,00 205,00
3 0,3113 0,3063 -4,94 -4,94 400,12 200,06 200,06
4 0,0038 0,0037 -0,06 -0,06 400,00 200,00 200,00
5 0,0000 0,0000 -0,00 -0,00 400,00 200,00 200,00

0 200,00 100,00 100,00

1 -4,6125 3,0716 292,86 -82,54 410 ,32 392,86 17,46
2 0,7424 -0,3953 -22,02 13,55 401,86 370,84 31,02
3 0,4216 -0,3727 -12,81 12,74 401,79 358,03 43,76
4 0,3846 -0,3417 -11,72 11 ,64 401,71 346,32 55,39
5 0,3512 -0,3137 -10,73 10 ,65 401,63 335,59 66,04
--------- - --
10 0,2287 -0,2086 -7,07 7,00 401,25 293,65 107,60
30 0,0518 -0,0488 -1,63 1,61 400,37 223,48 176,88
38 0,0299 -0,0283 -0,95 0,93 400,22 213,77 186,45
50 0,0134 -0,0127 -0,42 0,42 400,10 206,23 193,87
101 0,0005 -0,0004 -0,01 0,01 400,00 200,22 199,78

Tableau 2 - Calcul du réseau de la figure 1. Le tronçon 1 est commun aux deux mailles. En
haut, Newton-Raphson. En bas, maille par maille. Unités: Ils et m.

En supposant que le calcul démarre avec les débits vant être la partie commune des deux mailles ouvertes.
de l'itération 10, ou avec ceux d'une itération de rang Un examen rapide des matrices jacobiennes correspon-
supérieur, il apparaît que la convergence resterait lente dantes montre que le choix fait à l'origine est le pire et
avec des débits initiaux pourtant proches de la solution. qu'à l'opposé le meilleur système est celui obtenu en
L'observateur attentif acquiert la certitude que les deux prenant le tronçon 3 comme partie commune, ce qui
mailles s'équilibreront difficilement dès l'itération 3, et revient à substituer /23 à /13 dans l'équation 21, d'où:
au plus tard à l'itération 4 ou 5. Cette certitude est basée
sur le fait symptomatique que les deux fermetures en /23 =- 3,75 X 10-6 Q~
cote d'une part et les deux corrections de débit d'autre
part sont de signes opposés et proches en valeurs
+ 1,25 x 10- 6 Q~ + 0,10 = ° (26)
absolues. /13 3,75 x lO- s Qi
La conclusion (provisoire) est évidente : il faut pros-
crire à jamais le calcul maille par maille car il n'est pas
possible d'ignorer les éléments non diagonaux de la
+ 1,25 x 10- 6 Q~ - 6,05 = °
matrice jacobienne. Et en effet, pour l'itération 6 par Comme les modifications sont évidentes, les équa-
exemple, on obtient, avec les débits de la fin de l'ité- tions 22 et 23 ne seront pas réécrites, mais il est inté-
ration 5 : ressant de donner la nouvelle matrice jacobienne de
l'itération 6 :

l27
326
Js = 10- 4
[ 301
301]
303
(24)
l' = 10-
4 r 2] (27)
s 2 303
Tous les éléments sont du même ordre de grandeur
La différence avec la matrice relative au maillage
mais, ainsi que le montre bien l'équation 23, parce que
d'origine, équation 24, est éloquente. Il est maintenant
chacun d'entre eux contient:
tout à fait normal d'ignorer les éléments non diagonaux,
puisqu'ils sont effectivement négligeables par rapport
aux éléments diagonaux. La partie supérieure du tableau
qui est la dérivée partielle prépondérante. Or, il en est 3 le prouve : la solution est atteinte en 4 itérations
ainsi parce que le tronçon 1 est commun aux deux nouvelles.
mailles. Un autre calcul, qui ne figure pas sur le tableau 3
Cette dernière observation permet de sortir de l'im- pour des raisons d'encombrement, montre - ce qui est
passe. Le maillage d'origine ayant fait la preuve de son tout à la fois intuitif et facile à vérifier - que la situation
inadaptation, rien ne s'oppose à un changement de mail- est aussi énergiquement redressée si la maille 2 + 3 est, à
lage à partir de l'itération 6 (par hypothèse, car le rang partir de l'itération 6, ajoutée au système d'origine, et
de l'itération est secondaire). Il a déjà été remarqué que non pas substituée à l'une des équations du système
le choix était possible entre trois systèmes d'équations, d'origine. Un procédé aussi peu orthodoxe (système de
chacun des trois tronçons du réseau de la figure 1 pou- 3 équations non indépendantes à 2 iIlconnues) ne se
 A. DIVENOT 371

itération f f 13 0°23 0°13

23 °1 °2 °3

6 -0,3169 0,0413 118,14 -1,35 400 ,28 217,44 182,84

7 -0 ,0355 0,0082 17,01 -0,27 400,01 200 ,43 199,58
8 -0 ,0009 0,0002 0,43 -0,01 400,00 200 ,00 200,00
9 -0 ,0000 0,0000 0,00 -0,00 400 ,00 200 ,00 200 ,00

a 200 ,00 100,00 100 ,00

1 0,0750 -4,5492 -75,00 302,02 502,02 175,00 327,02
2 0,1188 3,4930 -55,79 -91,13 410 ,89 230,79 180,10
3 -0 ,0592 0,3349 27,14 -la ,69 400 ,20 203,65 196,55
4 -0,0072 0,0062 3,58 -0,20 400,00 200 ,07 199,93
5 -0,0001 0,0000 0,07 -0,00 400 ,00 200 ,00 200,00
6 -0 ,0000 0,0000 0,00 -0,00 400 ,00 200,00 200,00

Tableau 3 - Calcul maille par maille du réseau de la figure 1. Le tronçon 3 est commun aux
deux mailles. En haut, reprise du calcul du bas du tableau 2. En bas, calcul normal. Unités:
I/s et m.

justifie pas dans le cas simple de l'exemple mais il est étrangers à Cross. Même sans aller jusqu'à la dynami-
bien utile pour traiter les réseaux réels. sation, le fait que Cross soit muet sur le choix des mailles
La partie inférieure du tableau 3 est relative au dérou- rend déjà ambiguë l'expression, pourtant souvent uti-
lement du calcul lorsque le maillage est bien choisi dès lisée : "le calcul de tel réseau par la méthode de Cross
l'origine. Le rapprochement avec la partie supérieure du demande tant d'itérations". Au fait, quelles sont les
tableau 2 conduit à la conclusion que, dans ce cas, la mailles Cross conformes du réseau de la figure 1 ?
méthode de Newton - Raphson est à peine plus efficace
(en nombre d'itérations, par en durée de calcul) que le
calcul maille par maille. Pour un petit réseau comme Principe de la méthode
celui-ci, il n'est d'ailleurs par rare de constater le même
nombre d'itérations. Le raisonnement ne s'applique plus au mini-réseau
D'autres calculs, qu'il n'est pas possible de reproduire de la figure 1 mais, très généralement, à un réseau
ici mais qui sont faciles à réaliser, démontrent en com- important. Toute maille du réseau a au moins un tron-
plément: çon qui lui est propre (les exceptions qui peuvent se
le cas présenté n'a rien d'exceptionnel. En première produire avec les mailles ouvertes ne modifient pas le
approximation (les débits interviennent également), raisonnement). Soient Mi et M j les ensembles d'élé-
plus les résistances sont différenciées, plus le choix du ments appartenant respectivement aux mailles i et j.
maillage est important; Ces ensembles sont la réunion des parties propres des
peu importe la formule de perte de charge utilisée, deux mailles, Pi et Pj' et de leur partie commune Pc :
même si l'exposant Œ est différent de 2, même si la
formule est celle de Colebrook - White; Mi = Pc + Pi (28)
il est toujours préférable d'appliquer immédiatement M j = Pc + Pj
la correction de chaque maille. L'avantage est d'au-
tant plus grand que la convergence est moins rapide. Quatre éléments de la matrice jacobienne sont re-
latifs à ces deux mailles: les éléments diagonaux aii et
En résumé, le calcul maille par maille peut, à la limite, ajj , et les éléments non diagonaux aij et aji (équations
ne pas demander plus d'itérations que l'application de 15, 22, 23). Soit, en désignant par <P", <Pi et <Pi les déri-
la méthode de Newton - Raphson. La répartition initiale vées partielles de premier ordre des parties Pc, Pi
des débits est secondaire. En revanche, le choix des etPj :
mailles, qui est fondamental, peut être optimisé et
dynamisé. aii = <Pc + <Pi
Dans ce chapitre, aucune référence n'a été faite à la a jj = <Pc+ <Pj (29)
méthode de Cross [3]. Pourtant, l'auteur n'oublie pas aij = aji = <Pc
que Cross a ouvert la voie à tous ceux qui ont apporté
leur contribution dans ce domaine. Ainsi, appliquer la Il a été démontré qu'il faut éviter de former des
méthode de Newton - Raphson, c'est traiter simulta- paires de mailles conduisant à un <Pc grand par rap-
nément des équations de Cross. Ces mêmes équations port à <Pi et <p .• Fait favorable, la matrice jacobienne
sont à la base de la méthode proposée ici. Mais l'algo- contient beauc6up d'éléments nuls car Pc est fréquem-
rithme de sélection des mailles qui sera explicité plus ment vide, ce qui entraîne <Pc = O. Les deux mailles
loin et plus encore le concept de dynamisation du sont alors totalement indépendantes et l'équilibrage de
maillage qui est l'aspect essentiel - sont totalement l'une n'interfère pas du tout avec l'équilibrage de l'autre.
372 LA HOUILLE BLANCHE/N° 6-1980

S'il pouvait en être amsl pour toute paire de mailles, L'expérience acquise avec le mini-réseau de la figure 1
tous les éléments non diagonaux seraient nuls et le est généralisable. 11 n'est pas du tout nécessaire de con-
calcul du réseau ne poserait aucun problème. Mais il naître la matrice jacobienne, qui reste parfaitement im-
est vrai que le réseau ne serait pas maillé, au sens habi- plicite. L'attitude de deux mailles dont le 'Pc n'est plus
tuel du terme. négligeable est un message suffisamment clair. Les deux
Faute de pouvoir annuler tous les 'Pc' l'objectif de la mailles s'attardent, convergent mal, leurs fermetures en
méthode est de faire en sorte que, pour toute paire de cote et corrections de débit sont de signes opposés et
mailles, 'Pc soit négligeable devant 'Pi et 'PI' Alors, tous proches en valeurs absolues. La formation de la nouvelle
les éléments non diagonaux sont négligeables devant la maille, k, nécessite l'identification de la partie com-
diagonale principale, et le calcul maille par maille est mune par intersection (équation 28) :
justifié. Dans la pratique, l'objectif n'est pas pleinement
atteint mais il l'est suffisamment pour commencer le Pc =Mj nMI (32)
calcul dans de bonnes conditions. Plus tard, le maillage puis l'union des deux parties propres:
est amélioré, si nécessaire et juste de ce qui est
nécessaire. Mk = (Mi - Pc) U (MI - Pc) (33)
Le discriminant utilisé pour évaluer les poids respec-
tifs des éléments de la matrice jacobienne est la résis- L'ensemble M k peut être substitué à Mi ou MI' mais
tance cumulée des tronçons constitutifs. Chacune des il est plus simple de ne pas bouleverser le vecteur du
quantités 'Pc' 'Pj et 'PI est une somme de dérivées de maillage (voir ci-après) et d'ajouter M k au maillage
pertes de charge (équations 18 et 19) : existant. Selon l'option exercée, l'amélioration est dite
être réalisée par modification du maillage ou par mail-
dh lage supplémentaire. Dans le premier cas, le nombre
_t = aR IQ 1",-1 (30)
dQt t t strictement nécessaire et suffisant de mailles indépen-
dantes est préservé. Dans le second cas, les mailles sup-
ou de hauteurs de pompage (équations 18 et 20) : plémentaires sont redondantes.
L'opération peut être répétée plusieurs fois au cours
d'un calcul. Le maillage est alors véritablement dyna-
misé, ce qui permet le traitement approprié des interac-
tions entre les mailles, dès qu'elles se manifestent. Le
La valeur de a dans l'équation 30 est constante ou concept est extrêmement puissant. L'équilibrage des
évolue dans une fourchette étroite, de l'ordre de 1,8 à gros réseaux demande à peine plus d'itérations qu'avec
2. Le débit Qt peut varier très largement d'un tronçon la méthode de Newton-Raphson, mais moins de temps
à l'autre mais généralement beaucoup moins que la car l'exécution de chaque itération est beaucoup plus
résistance R t . Celle-ci est extrêmement différenciée, rapide. Bien qu'aucune tentative n'ait encore été réa-
ce. qui en fait le bon critère : à longueur et rugosité lisée, il n'est pas impossible que la méthode puisse être
identiques, elle est en effet inversement proportion- appliquée dans d'autres domaines.
nelle à la puissance 5,25 du diamètre, environ (l'expo-
sant dépend de la formule de perte de charge utilisée). Arborescence
Ainsi par exemple, si les diamètres sont dans le rapport
1 à 5, les résistances sont dans le rapport 5 000 à l, Conformément à ce qui vient d'être développé, le
environ. but est d'extraire du graphe un maillage proche de
S'il est fait usage de la formule de Colebrook-White, l'optimum. L'arborescence (le terme est utilisé dans un
la résistance R t est fonction du débit Qt' mais une sens large) est particulièrement bien adaptée au pro-
simple approximation suffit puisqu'il s'agit seulement blème car elle permet, selon le sens de l'exploration,
d'évaluer l'ordre de grandeur relatif des dérivées, et que le calcul des débits initiaux ou celui des niveaux piézo-
l'erreur est admise. A ce stade, il est de même possible métriques équilibrés.
d'assimiler, même grossièrement, chaque courbe de Les tronçons sont regroupés en arêtes, échelons
pompage à la courbe de perte de charge d'un tronçon intermédiaires entre les tronçons et les parties propres
fictif. ou communes des mailles. Une arête est une séquence
Comme toutes les dérivées sont fonction des débits, de tronçons, non orientée tant qu'elle n'est pas intégrée
et que les débits initiaux sont très éloignés de la solution, dans l'arborescence. Ses nœuds courants sont de degré 2,
les éléments de la matrice jacobienne implicite évoluent ses nœuds d'extrémité de degré différent de 2. La fi-
au cours du calcul. Certains éléments non diagonaux gure 2 donne le détail de l'arête 1-20 de la figure 3. Les
peuvent en conséquence devenir gênants, soit parce que nombres supérieurs sont les numéros de nœuds, les
l'algorithme de sélection initiale des mailles - qui sera nombres inférieurs les numéros de tronçons. L'arête
explicité plus loin - n'est pas parfait, soit parce que les compte 3 tronçons. Les nœuds 14 et 5 sont courants.
débits corrigés renforcent singulièrement le poids relatif
de quelques éléments non diagonaux. Il est alors néces-

6--3--'.4--a-----:>-----2K
saire et parfois indispensable d'améliorer le maillage car
il suffit que deux mailles deviennent farouchement anta-
gonistes pour que le réseau soit rendu incalculable.
Qu'elle soit nécessaire ou indispensable, l'amélioration ligure 2
est toujours largement bénéfique. Figure 2 - Détail de l'arête 1-20 du réseau de la figure 3.
 A. DIVENOT 373

Les degrés des nœuds d'extrémité, 1 et 20, sont respec- Maillage

tivement de 1 et 4. La résistance d'une arête est le cumul
des résistances hydrauliques des tronçons constitutifs. L'examen des caractéristiques de chaque chemin
La racine de départ est l'un des nœuds de niveau permet de savoir s'il génère une maille et, dans l'affir-
piézométrique déterminé. La méthode prévoyait à l'ori- mative, si elle est ouverte ou fermée. Comme les nœuds
gine le choix de la racine par le programme (meilleure 1, 3 et 6 ont leurs niveaux piézométriques déterminés,
conductance au point d'entrée dans le réseau) avec la le premier chemin de la figure 4 génère la maille ouverte
1-20-9-3. Simultanément, le programme laisse une trace
3 6
aux nœuds 20 et 9.
3
Le deuxième chemin (9-16-6) est l'amorce d'une
maille ouverte qui, outre le chemin lui-même, peut
9 comprendre 9-3 ou 9-20-1. C'est grâce à la trace laissée
par la première maille que le bouclage s'effectue sur le
nœud de niveau piézométrique déterminé d'écart mini-
mal, tel que le cumul des résistances à partir du nœud
figu re 2 confluent (le nœud 9) soit minimal. D'où la nature de
la trace: lors de l'examen de la première maille ouverte,
le programme a noté que, par exemple, la distance 20-1
est moindre que la distance 20-9-3, et la distance 9-3
figure 3 figure 4 moindre que la distance 9-20-1. Les distances sont
Figure 3 - Exemple de réseau. Les niveaux piézométriques des mesurées en résistances cumulées. Lors de l'examen
nœuds 1,3 et 6 sont déterminés.
du deuxième chemin, la trace laissée en 9 indique que
Figure 4 - Arborescence du réseau de la figure 3. la maille ouverte doit être complétée par 9-3, c'est-à-
dire en prenant la direction de 3, et non en prenant
possibilité pour l'utilisateur de forcer le choix dans l'en- celle de 1. Il est évident que le problème posé par le
semble des nœuds de niveaux piézométriques déter- réseau de la figure 1 est ainsi radicalement réglé, mais
minés. L'expérience a montré que c'était une précau- le procédé va bien au-delà de ce cas simple.
tion inutile qui d'ailleurs faisait, au moins partielle- La formation des mailles fermées est immédiate :
ment, double emploi avec le nœud d'écart minimal, au troisième chemin (20-16), il suffit d'ajouter 20-9-
autre notion qui sera bientôt explicitée. 16, en remontant vers la racine.
Chaque arête n'est parcourue qu'une seule fois dans Enfin, le dernier chemin (20-7) est une antenne
un sens précis. La progression s'effectue en marquant qui ne génère pas de maille. Lors de la formation de
les arêtes intégrées et les nœuds de passage. A chaque l'arborescence, les arêtes les moins résistantes ont
extrémité d'arête, l'arête choisie pour poursuivre le constamment été sélectionnées de façon préférentielle.
chemin est celle de résistance minimale. Le chemin Les arêtes de forte résistance n'ont été intégrées que
est totalement constitué lorsque le nœud d'extrémité lorsqu'il ne pouvait plus en être autrement. Il est main-
de la dernière arête incorporée est de degré 1 ou déjà tenant clair qu'il s'agissait, autant que possible pour
marqué. Chaque fois que ce deuxième cas se produit, un algorithme d'exécution rapide, de parcourir préfé-
l'arborescence est "fermée", ce qui est contraire à la rentiellement les arêtes susceptibles d'appartenir à des
définition habituelle. La première arête du chemin parties communes, et au contraire de rejeter en parties
suivant est l'arête pendante (l'une de ses extrémités, propres les arêtes de résistances relativement élevées.
au moins, est marquée) de résistance minimale. La re· Ainsi, il est facile de vérifier que la maille fermée 20-9-
cherche est poursuivie jusqu'à ce que toutes les arêtes 16-20 a pour partie propre l'arête 20-16 dont la résis-
soient intégrées dans l'arborescence. Si le graphe n'est tance est supérieure à celles des arêtes 20-9 et 9-16.
pas connexe, ses composantes sont traitées succes- Il est également devenu évident qu'il est indispensable
sivement. de former les arêtes car la considération des seuls tron-
Sur la figure 4, qui est à rapprocher de la figure 3, çons pourrait conduire à de graves mécomptes.
on distingue le chemin 1-20-9-3, le chemin 9-16-6, Comme pour l'arborescence, l'organisation du vec-
le chemin 20-16, le chemin 20-7. Au fur et à mesure teur du maillage est fonction du programme. Beaucoup
de l'avancement, un tableau à une dimension, le vec- de temps sera gagné si toutes les informations relatives
teur de l'arborescence, est constitué. Son contenu à chaque maille sont directement accessibles au cours
dépend de l'organisation du programme, par exemple du processus itératif : identification des tronçons, des
(figure 2) pour le début du premier chemin, 1-3-14-8- pompes de refoulement ou de surpression, des clapets,
5-7-20 si toutes les informations, alternativement des stabilisateurs de pression, etc.
nœuds et tronçons, figurent dans le vecteur. Le maillage issu de l'arborescence est toujours bon.
Les débits initiaux sont déterminés en remontant Toutefois, il a été exposé qu'il peut ne pas être optimal
le vecteur de l'arborescence, dans le sens contraire à et que la surveillance du déroulement des itérations
son orientation. Les consommations sont cumulées est partie intégrante de la méthode, afin de repérer les
de proche en proche dans les tronçons, jusqu'à la racine. mailles qui ont tendance à s'attarder. 11 ne peut s'agir
Le calcul des niveaux piézométriques finals sera pra- que d'un nombre restreint de mailles, d'où un rema-
tiqué dans l'autre sens, lorsque le réseau sera équilibré. niement limité du maillage. Supposons, figure 4, que
Mais, auparavant, le maillage est déduit de ce même soient en cause les mailles 3-9-16-6 et 20-9-16-20.
vecteur. Conformément aux équations 32 et 33, la maille nou-
374 LA HOUILLE BLANCHE/N° 6-1980

velle, 3-9-20-16-6, est formée par union après élimi- caractère différentiel, en performances hydrauliques
nation de la partie commune. Elle pourrait être subs- et en coûts. Les différences de performances hydrau-
tituée à l'une ou l'autre des mailles génératrices mais liques risquent d'être gravement altérées si les résultats
il est beaucoup plus simple et tout aussi efficace de de chaque calcul sont approximatifs.
l'ajouter en surnombre aux mailles existantes, ce qui se La dernière colonne du tableau 4 est produite à
traduit par une simple extension du vecteur du maillage. simple titre d'information. Si le principe de la dynami-
sation repose sur une base solide, son application est
encore empirique et peut vraisemblabement être amé-
Application liorée. Le thème est soumis aux chercheurs et aux pra-
ticiens. Il convient de remarquer que la minimisation
PROGRES, qui est un programme basé sur la mé- de la durée du calcul doit normalement avoir la prio-
thode présentée, a déjà permis la réalisation de plusieurs rité sur la minimisation du nombre d'itérations et que
milliers de calculs dans d'excellentes conditions, aussi le point d'équilibre dépend, au moins partiellement,
bien sur de gros ordinateurs que sur des minis. Quel de l'ordinateur utilisé. Comme l'auteur n'a pas de con-
que soit le réseau et quelle que soit la version du pro- tribution fondamentale à présenter en la matière, on
gramme, la méthode n'a jamais été mise en défaut. voudra bien l'excuser de conserver durant quelque temps
Afin de produire quelques résultats tout en évitant l'avantage de l'expérience.
de longs développements, des exemples ont été recher- Pour les réseaux d'essai de 10 à 20 mailles du ta-
chés dans les références bibliographiques. Parmi les bleau 4, le nombre d'itérations en maillage dynamique
huit réseaux que l'on y trouve, quatre n'ont pas pu et haute précision est de 10 à 20. Dans les mêmes con-
être exploités soit parce que les courbes de pompage ditions, des réseaux réels de 50 à 150 mailles demandent
ne sont pas fournies [l, 2,4,8] soit parce que les débits 15 à 40 itérations, selon les données particulières du
aux nœuds ne paraissent pas cohérents [exemple 2 calcul. Un réseau de 1 000 tronçons et 250 mailles en-
de 17]. Les résultats obtenus avec les quatre autres viron (selon le calcul) est couramment équilibré en 25
réseaux sont reportés dans le tableau 4. à 45 itérations. Le même gros réseau a été traité en
Les nombres d'itérations des calculs effectués en maillage statique, afin de pouvoir citer les résultats
maillage statique sont aisément vérifiables puisqu'il ici : il a fallu 9 itérations pour la précision courante
suffit d'appliquer ce qui a été dit sur l'arborescence (1 lis), 36 itérations pour 0,05 lis, 73 itérations pour
et le maillage. La précision de 1 lis est une précision la grande précision, 0,05 lis et 0,5 mm, en réalité
courante. Le calcul cesse lorsque toutes les corrections 8.10- 4 I/s et 0,5 mm.
de débit des mailles sont, en valeurs absolues, infé-
rieures à Il/s. La grande précision, qui est la précision
souhaitable, n'est atteinte que lorsque toutes les valeurs Conclusion
absolues des corrections de débit sont inférieures à
0,05 lis et toutes les valeurs absolues des fermetures Le maillage optimal serait celui qui rendrait les élé-
en cote inférieures à 0,5 mm. La deuxième contrainte ments non diagonaux de la matrice jacobienne cons-
est habituellement la plus sévère et il n'est pas excep- tamment négligeables. L'objectif ne peut pas être
tionnel que le millième de lis soit nécessaire pour la atteint à coup sûr car le poids relatif des éléments
satisfaire. évolue au cours du calcul. D'où la double nécessité
Aucune comparaison n'est possible avec les nombres de soigneusement sélectionner les mailles initiales
d'itérations cités dans les articles référencés car les pré- puis de dynamiser le maillage.
cisions y sont omises ou ne sont pas transposables. Il La méthode présentée est un tout dont les mots-
est naturel de choisir la grande précision quand le clés sont arborescence, maillage, retouches strictement
programme permet de l'atteindre aisément pour une indispensables. Le mode de formation de l'arbores-
légère augmentation de la durée du calcul. Mais une cence et le concept du nœud de niveau piézométrique
bonne précision est toujours souhaitahle parce que déterminé de moindre écart conduisent à un maillage
les enseignements tirés des calculs sont souvent de pour le moins satisfaisant. Les imperfections éven-

Nombre de ~1a i 11 age sta ti que ~'a i 11 age dynami que

Référence
mai 11 es 1 l/s 0,05 l/s - 0,5 mm 0,05 l/s - 0,5 mm

[9]-[ 13 1 16 15 29 19
[11 ] - [14] 20 3 36 16
[12 ] 14 10 27 17
[17 ] 8 5 28 11

Tableau 4 - Nombres d'itéractions obtenus avec PROGRES.

 A. DIVENOT 375

tuelles sont éliminées au cours du processus itératif - stations de pompage, clapets, stabilisateurs de pres-
en prenant en considération un nombre restreint de sion - que les singularités de calcul. L'aspect essen-
mailles supplémentaires. tiel de )a méthode est probablement d'être à l'écoute
L'ensemble est efficace, performant, et permet de d'une matrice jacobienne qui reste parfaitement
traiter sans difficulté aussi bien tous les ouvrages implicite.

Liste des principaux symboles

niveau piézométrique du nœud i. résistance hydraulique du tronçon t.
perte de charge dans le tronçon t. nombre de tronçons du réseau.
matrice jacobienne. nombre de nœuds de niveau piézométrique
nombre de mailles fermées du réseau. déterminé du réseau.
nombre total de mailles du réseau. exposant du débit dans la loi de perte de
charge.
ensemble des éléments appartenant à la
maille i. correction de débit appliquée aux éléments de
la maille m.
N nombre de nœuds du réseau.
dérivée partielle de premier ordre de Pc,
partie commune de deux mailles.
dérivée partielle de premier ordre de Pi'
partie propre de la maille i.
débit du tronçon t.

Liste bibliographique
Il) CHANDRASHEKAR M. - Extended set of components equations, Journal of the Hydraulics Division, ASCE,
in pipe networks, Journal of the Hydraulics Division, février 1972, p. 335-344.
ASCE, janvier 1980, p. 133-149. (11) LAM C.F. and WOLLA M.L. - Computer analysis of
(2) ÇHANDRASHEKAR M. and STEWART K.H. - Sparsity water distributiom systems : Part II - Numerical solu-
oriented analysis of large pipe networks, Journal of the tion, Journal of the Hydraulics Division, ASCE, mars
Hydraulics Division, ASCE, avril 1975, p. 341-355. 1972, p. 447-460.
(3)'CROSS H. - Analysis of f10w in networks of conduits or (12] LEMIEUX P.F. - Efficient algorithm for distribution
conductors, University of Illinois Engineering Experiment networks, Journal of the Hydraulics Division, A SCE,
Station, Bulletin n° 286, 1936. novembre 1972, p. 1911-1920.
(4) DONACHIE R.P. - Digital program for water network [13] LIU K.T.H. - The numerical analysis of water supply
analysis, Journal of the Hydraulics Division, ASCE, mars networks by digital computer, Proceedings, 13th Congress
1974, p. 393-403. of the International Association for Hydraulic Research,
(5) DUBIN C. - Le calcul des réseaux maillés, La Houille Kyoto, Japan, Vol. l, Subject A, septembre 1969,
Blanche, n° 3,1947, p. 213-227. p.35-42.
(6) EPP R. and FOWLER A.G. - Efficient code for steady- (14] MARTIN D.W. and PETERS G. - The application of
state f10ws in networks, Journal of the Hydraulics Divi- Newton 's method to network analysis by digital com-
sion, ASCE, Janvier 1970, p. 43-56. puter, Journal of the Institution of Water Engineers,
(7) JEPPSON R.W. and DAVIS A.L. - Pressure reducing Vol. 17, 1963,p.115-129.
valves in pipe network analyses, Journal of the Hydraulics [15) SHAMIR U. and HOWARD C.D.D. - Water distribution
Division, ASCE, juillet 1976, p. 987-1001. systems analysis, Journal of the Hydraulics Division,
[8) JEPPSON R.W. and TAVALLAEE A. - Pumps and reser- ASCE, janvier 1968,p. 219-234.
voirs in networks by linear theory, Journal of the Hy- (16) VOYLES C.F. and W1LKE H.R. - Selection of circuit
draulics Division, ASCE, mars 1975, p. 576-580. arrangements for distribution network analysis by the
[9) KESAVAN H.K. and CHANDRASHEKAR M. - Graph- Hardy Cross method, Journal of the American Water
theoretic models for pipe network analysis, Journal of the Works Association, mars 1962, p. 285-290.
Hydraulics Division, ASCE, février 1972, p. 345-364. [171 WOOD D.1. and CHARLES C.O.A. - Hydraulic network
(10) LAM C.F. and WOLLA M.L. - Computer analysis of analysis using linear theory, Journal of the Hydraulics
water distribution systems : Part 1 - Formulation of Division, ASCE, juillet 1972, p. 1157-1170.