Un modèle numérique, permettant d’obtenir de manière automatique des informations sur la propagation des ondes de coup de bélier dans les réseaux maillés de conduites, est présenté. Le modèle tient compte des bifurcations et des dérivations dans les réseaux de conduites ainsi que des pertes de charges. Il est constitué d’un système de deux équations aux dérivées partielles non-linéaires de type hyperbolique résolu par la méthode des caractéristiques. L’algorithme numérique ainsi construit fournit une estimation des pressions maximales dans le fluide et des contraintes maximales dans les parois dues à la fermeture rapide de vannes. Dans certains cas, la contrainte maximale peut devenir supérieure à la contrainte admissible et provoque la rupture des conduites. La dangerosité d’un défaut de type cratère de corrosion a été analysée en déterminant la distribution des contraintes en tête de ce défaut. Les résultats obtenus permettent de calculer le facteur d’intensité de contraintes d’entaille appliqué. Cette grandeur est insérée dans un diagramme Intégrité-Rupture de type SINTAP et les nœuds en situation critique sont déterminés. Il y a risque de rupture si le facteur de sécurité est inférieur à 2. Les résultats de calcul montrent que presque tous les nœuds du réseau analysé sont situés en dehors du domaine d’intégrité.

Dans cette étude, on examine, par simulation numérique, la possibilité de réduire les pressions provoquées par le phénomène du coup de bélier dans un réseau de conduites quasi-rigides en remplaçant l'une des conduites du réseau par une conduite viscoélastique en polymère. Le modèle numérique développé est constitué d'un système de deux équations aux dérivées partielles non-linéaires de type hyperbolique résolu par la méthode des caractéristiques avec mise en mémoire et interpolation. Le calcul au niveau des nœuds est effectué par la résolution d'un système linéaire avec prise en considération des conditions aux limites imposées. La loi de comportement de la conduite viscoélastique est décrite par le modèle de Kelvin-Voigt. Afin de simplifier les calculs, on se limite à l'élément élastique de ce modèle. Pour les applications, un code de calcul écrit en langage FORTRAN a été élaboré. Les résultats obtenus montrent bien l'effet de l'élasticité des parois sur l'évolution et l'amortissement des ondes de pression dans les réseaux de conduites.