Procédure
Construction d’un modèle TELEMAC 2D
1. Création d’un maillage sur BlueKenue
Ouvrir BlueKenue 64
1.1. Création d’un maillage linéaire (type Channel Mesher)
• Importer les couches dessinées préalablement sur QGIS dans BlueKenue :
Les couches qui doivent être importée sur BlueKenue sont :
o L’emprise du modèle (polygone ou ligne) ;
o Les lignes de contraintes (lignes) ;
o Les digues et merlons (lignes) ;
o Le lit mineur (lignes) ;
o Le Bâti (polygones ou lignes)
o Les autres singularités : seuils … (lignes)
Maillage des éléments linéaires singuliers :
o Lit mineur
o Bras de décharge
o Routes, voies ferrées, digues
On maille ces éléments à partir de 2 lignes « parallèles ».
• Importer les couches SHP dans Bluekenue :
File > Import > ArcView Shape File
• Visualiser les entités importées dans BlueKenue :
Faire un glisser-déposer de l’entité à visualiser de « Data Items » vers « 2D View »
• Création des « Channel Mesher » :
Avant de générer un maillage entre deux lignes parallèles (lit mineur, digues, bras secondaires …), il faut enregistrer
les entités .shp en fichier .i2s (extension de Bluekenue).
- Création de deux « rives » (left and right banks) :
On réalise une manipulation rapide afin de créer la structure des fichiers i2s dans lesquels on viendra ajouter les
lignes correspondant à nos 2 « rives ».
Pour créer un fichier i2s polyligne : Cliquer sur l’icône entourée en rouge, puis dessiner une polyligne temporaire,
puis échappe, et la nommer RD (ex: rive droite).
Supprimer la polyligne temporaire. BluKenue ne
permet pas la suppression d’une entité. Pour
contourner cela : clic droit > cut.
Ajouter la ligne de rive réelle dans le fichier i2s.
Dans le SHP chargé, clic droit sur la polyligne de
rive réelle > copy
Sur la couche RD i2s créée, clic droit > paste.
Enregistrer la couche de rive droite (RD) ainsi créée en fichier i2s : sélection de la couche > file > save
Répéter la même opération pour la rive gauche (RG).
- Génération du maillage « Channel Mesher » :
A ce stade on peut générer un maillage par un « Channel
Mesher »
File > New > T3 Channel Mesher …
Valider par OK
Un new T3 Channel Mesh vient d’être crée et apparait dans
l’arborescence à gauche.
Il faut ensuite faire glisser dans ce répertoire les entités
nécessaires :
- RD précédemment créée
- RG précédemment créee
Accéder au paramètres du T3 Channel Mesh par un double-clic sur l’élément. Les différents paramètres à fixer sont :
- Le nombre de lignes de nœuds en travers à générer dans le channel : Cross Channel Node Count
- L’intervalle entre deux nœuds d’une ligne de nœuds : Along Channel Interval
Cross Channel Node Count
Along Channel Interval
Valider le paramétrage, puis Run, puis Ok.
Un Maillage (Mesh) vient d’être créé. Le glisser-déposer dans la 2D View pour le voir. File > Save pour l’enregistrer.
Répéter cette procédure de Channel Mesher pour toutes les autres singularités que l’on veut sous-mailler (voies
ferrées, routes, digues …)
1.2. Création d’un maillage global (type Mesh Generator)
Pour finaliser le maillage global du modèle, on ajoute toutes les singularités à faire figurer dans le modèle (lignes de
contraintes, bâti …) ainsi que les sous maillages précédemment générés et l’emprise globale du modèle.
Les fichiers importés en SH (bâti, lignes de contraintes, emprise du modèle …) doivent être enregistrés en i2s, puis
ré-échantillonnés afin de s’adapter à la densité du maillage souhaité (Resample).
• Pour ré-échantillonner un objet :
Clic droit dur la couche visée > Resample
Dans la fenêtre qui apparait (ci-contre), choisir
la méthode et le pas d’échantillonnage.
Valider par OK.
Tous les éléments sont désormais disponibles
pour générer le maillage global.
• Génération du maillage global :
Créer un Mesh Generator : File > New > T3 Mesh Generator
L’arborescence suivante apparait dans le volet de navigation à
gauche de la page :
Faire glisser dans chaque sous-groupe les éléments associés :
- Outline : Emprise du modèle
- HardLines : Glisser les lignes « majeures », indispensables à représenter fidèlement dans le maillage. Lignes
pour lesquelles on souhaite que BK maille obligatoirement sur chacun des sommets échantillonnés le long
(bâtis, ponts, …)
- SoftLines : Glisser les lignes plus « facultatives ». Lignes pour lesquelles BK respectera la position sans pour
autant mailler sur chacun des noeuds échantillonnés le long.
- SubMeshes : Glisser les sous-maillages précédemment générés (lits mineurs, digues, routes …)
- Density : Glisser la couche de zone de densités de maillage
- HardPoints : Par analogie avec le HardLines, couche de points sur lesquels on souhaite que le maillage
s’appuie obligatoirement.
Pour lancer la génération du maillage : double-clic sur new T3Mesh dans le volet à gauche. Paramétrer :
- La taille des mailles (default Edge Length)
- La densité maximale souhaitée (Density).
- Resample outline = OK
- Edge Growth Ration = 1,2 (conserver la valeur ou faire
varier très légèrement)
Lancer la génération par RUN.
Glisser le maillage généré dans le viewer afin de contrôler qu’il
correspond aux attentes.
1.3. Correction du maillage dans des endroits spécifiques
Afin de corriger / modifier localement le maillage, on procède à un remaillage :
- Suppression du maillage dans la zone à remailler :
o Dessiner un polygone autour de la zone à remailler en utilisant le bouton de création
« New Closed Line »
o Cliquer une fois sur le maillage dans le volet de navigation
à gauche puis Edit > T3Mesh > Delete element in polygon
…
o Choisir dans la liste le polygone précédemment créé.
Le maillage est ainsi supprimé à l’intérieur des nœuds entourés
par le polygone.
- Création du nouveau maillage pour la zone à remailler :
o Extraire le contour de la zone précédemment découpée.
Cela permet ensuite de créer un sous-maillage dont la
limite sera exactement identique à celle de la zone
découpée.
Sélectionner le maillage dans le volet à gauche puis Tools
> T3mesh > Extract Edges (Shorelines)
o Sélectionner le polygone de contour correspondant à la
zone découpée et le Copier / Coller dans un nouveau
fichier i2s (création d’une ligne fermée temporaire > cut,
puis copy / paste de la ligne de contour vers ce nouveau
fichier)
o Générer le maillage comme décrit précédemment sans activer le ré-échantillonnage de l’Outline.
L’enregistrer en t3s.
- Fusion du sous-maillage avec le maillage découpé :
Sélectionner l’un des deux maillages à fusionner dans le volet à gauche.
Puis Edit > T3 Mesh > Merge Mesh
Choisir le second maillage à fusionner dans la liste.
Les deux maillages sont ainsi fusionnés. La fusion fonctionne parce que le
contour du sous-maillage est exactement identique à celui de la zone
découpée (même sommets et même lignes).
2. Edition du maillage avec TriGrid
L’objectif est d’ouvrir le maillage créé par bluekenue dans TriGrid afin de réaliser des modifications dessus (ajout de
lignes, déplacement de nœuds, division de mailles, swapping de lignes, détection de voisins trop nombreux …)
2.1. Prérequis, conversion du format de maillage
• Création d’un fichier Selafin avec BlueKenue :
Dans BK : File > New > SELAFIN Object. Faire glisser le maillage dans « newSelafin » et l’enregistrer (extension .slf).
• Convertir le fichier en .ngh :
Il faut utiliser la petite routine « [Link] » (copier la routine dans le dossier où est rangé le .slf).
Ouvrir la routine, renseigner le nom du fichier de maillage d’origine (.slf), puis le nom du fichier de
maillage de sortie (.ngh). Attention à bien écrire les extensions
2.2. Edition du maillage avec TriGrid
Pour éviter le plantage du TriGrid clic droit sur .exe, Puis
Propriétés, Compatibilité, Exécuter ce programme en
compatibilité pour Windows XP
Lancer le programme TriGrid.
Pour accéder à l’éditeur : Application > Editor
Pour charger le maillage souhaité : File > Open neighbour file (Text) puis sélectionner le fichier .ngh
• Les différents outils de navigation > Display :
o Zoomer / Dézoomer : > In / > Out (> Out (multiple) permet de dézoomer sur l’ensemble du maillage).
o Afficher les mailles : > Draw all edges
o N’afficher que les contours : > Draw Outlines
• Les outils de recherche de nœuds / éléments > Info :
o Chercher un nœud par index : > Node, puis cliquer sur INDEX (bas à droite de l’écran), puis entrer le numéro
du nœud dans la fenêtre qui apparait. Valider par OK.
• Les outils d’analyse du maillage > Test :
o > Vertices : Permet de localiser les points selon leur nombre de voisin. Dans NBGT/NBLT renseigner un
nombre max/min de nœuds voisins. Les nœuds ayant plus/moins de voisins que la valeur sont encadrés en
rouge.
Il faut tester : NBGT : 10 et NBLT : 3 et supprimer les nœuds concernés
• Les outils d’édition du maillage > Edit :
o > Add node : à éviter
o > Add line : permet d’ajouter une ligne entre deux nœuds. Clic sur le premier nœud puis clic sur le second.
Valider par OK.
o > Del node ou > Del line : permet de supprimer un nœud ou une ligne
o > Move Node : Déplacer un nœud. Clic sur le nœud à déplacer puis clic sur l’endroit où le repositionner.
o Attention : il faut forcément réaliser les 2 clics sinon TriGrid peut planter.
o Attention : s’il y a un nœud proche du point du 2eme clic, le nœud déplacer viendra s’y fusionner.
o > Insert node (on line) : permet d’ajouter un nœud sur une ligne. Clic sur la ligne à l’endroit où l’on souhaite
ajouter le nœud. Ce nœud est automatiquement connecté aux nœuds opposés à la ligne dans les mailles.
Nœud à ajouter sur la Ajout de 2 lignes de
ligne centrale jonction avec les nœuds
opposés à la ligne centrale
o > Exchange line (in quadrilateral) : permet d’échanger l’orientation d’une ligne entre deux mailles. Clic sur la
ligne à échanger
• Autres fonctionnalités de TriGrid :
o Rafraichir la vue : > Refresh !
o La dernière commande utilisée reste en mémoire à l’extrémité droite de la barre d’outil.
Une fois les modifications finies, enregistrer le fichier par File > Save As (Text), au format .ngh
3. Corrections du maillage – DETECTRIPLAT
La Moulinette [Link] permet d’analyser le maillage pour détecter les mailles à angles trop faibles ou les
segments trop courts, pouvant générer des erreurs de calcul par la suite.
Pour lancer la moulinette, placer l’exécutable dans le répertoire du maillage. Lancer l’exe et se laisser guider.
Renseigner les seuils de détection souhaités :
- Angles : usuellement 1° [pour des maillages lâches] à 0,1° [pour des maillages plus resserrés (urbain …)].
- Longueurs : usuellement de 1 m à 0,1 m.
Le programme génère 2 fichiers .txt finaux :
- Un fichier récapitulatif des mailles avec angles sous le seuil fixé
- Un fichier récapitulatif des segments à longueur inférieure au seuil fixé
Selon les « erreurs » détectées, corriger dans TriGrid ou dans BlueKenue les nœuds/mailles/segments identifiés.
4. Opérations sous Fudaa
Ouvrir le logiciel Fudaa ([Link]). La page qui s’ouvre constitue le Navigateur de fichiers de Fudaa.
Se placer dans le répertoire, là où est enregistré le fichier de maillage.
4.0. Vérification de la conformité du maillage avec Fudaa
Fudaa peut détecter les erreurs sur le maillage et indiquer les Id des éléments problématiques.
Dans le navigateur Fudaa, clic droit sur le fichier de maillage > Ouvrir avec > Post.
A l’ouverture, le message d’erreur que Fudaa affiche indique le nœud/l’élément concerné. Il faut ensuite repasser
par TriGrid pour corriger les erreurs sur le .ngh.
Si le maillage ne présente pas d’erreur, fudaa indique :
Il est conseillé de réaliser cette vérification en amont des étapes suivantes afin d’éviter de nouvelles bascules de
formats de fichiers (geo vers ngh puis ngh vers geo).
4.1. Renumérotation des nœuds du maillage
Suite aux modifications réalisées sur le maillage, notamment avec TriGrid, la renumérotation du maillage est
indispensable.
Dans le Navigateur Fudaa, naviguer vers le dossier de travail contenant le fichier .ngh corrigé final.
o Clic droit sur le fichier de maillage > Traitement du maillage
o Dans la fenêtre qui apparait renseigner : le fichier de maillage d’entrée (.ngh) et son format (Trigrid (ngh))
o Cocher Renumérotation
o Fichier de maillage de sortie : T2GEO_*.geo (préférentiel)
o Et le Format de fichier de sortie de maillage : Serafin
o Exécuter
Remarque :
Si le fichier ne s’ouvre pas ou que la renumérotation ne se lance pas, il y a probablement une erreur dans le maillage
d’entrée. Dans ce cas, vérifier que les vérifications préalables du maillage ont bien été réalisées.
4.2. Création d’un projet TELEMAC
Le fichier du maillage (.geo) est désormais crée. On crée le projet telemac à partir de ce maillage :
o Clic droit sur le fichier maillage > Créer projet Telemac
o Dans la fenêtre qui apparait :
- Fichier sérafin : déjà renseigné
- Fichier des paramètres : pré rempli mais peut être adapté (T2DCAS_*.cas)
- Charger un fichier dico : Aller chercher le dico pour la version T2D utilisée (pour la v7p2 utiliser le dico
de la v7p1)
- Valider
o Une nouvelle page de Fudaa apparait avec un message : « Les conditions initiales vont être initialisées »,
valider ce message, le maillage apparait.
• Pour accéder aux paramètres généraux du projet telemac
> Bouton Maison bleue à gauche
o Renseigner le nom du fichier de résultat : T2DRES_*.res
o Renommer si besoin le fichier des conditions aux limites :
T2DCLIM_*.conlim
• Pour visualiser le maillage
> Bouton Mesh rouge à gauche
Parcourir visuellement le maillage pour vérifier que les nœuds de bord sont bien placés comme prévu.
Attention notamment à ce qu’il n’y ait pas de « trous » dans le maillage (notamment des trous triangles assimilables
à des mailles). Ce type de bords peut générer des erreurs dans les calculs du logiciel par la suite.
• Gérer les variables des Noeuds
Ajouter des variables dans les nœuds par clic droit sur Nœuds > Gérer les
propriétés nodales
- Topographie du modèle (Bathymétrie)
- Surface libre du modèle
- Coefficient de frottement de fond
- Vitesses
- RUISSELLEMENT
ATTENTION : C’est à ce moment que l’on peut intégrer la variable RUISSELLEMENT pour les modèles de ruissellement.
Attention, il faut le faire avant d’activer les conditions initiales, sinon la variable ne s’inscrit que dans le PRE et pas
dans le GEO et il ne pleut pas dans le modèle (même si ça ne plante pas au démarrage). Par contre, si on veut
appliquer un coeff de ruissellement différent de 1, il ne faut pas activer les conditions intiales.
• Gérer les conditions initiales
Les conditions initiales sont stockées dans un fichier [Link]. Les différentes données stockées sont :
- La topographie pour le modèle (Bathymétrie)
- Le coefficient de frottement de fond
- La surface libre : c’est la hauteur d’eau + la cote au sol
Activer les conditions initiales : Ce n’est pas obligatoire
Clic droit sur CI: Nœuds > Activer les conditions initiales
ATTENTION, il ne faut pas activer les conditions initiales dans le cas d’un
modèle de ruissellement où on veut modifier le coeff. de ruissellement
• Appliquer la topographie sur le modèle
Dans un premier temps on applique la topographie générale sur tout le modèle. Puis localement écrase / remplace
par de la donnée topo particulière sur des emprises définies préalablement sur QGIS : lit mineur, digues, merlons, …
o Préparation de la donnée topo
2 méthodes possibles selon que l’on dispose de points topo ou d’un lidar
Méthode 1 : Utilisation d’une donnée Raster
Si la donnée topo est une donnée raster type MNT (RGE Alti, Lidar, …).
- Sur BK, extraire les nœuds du maillage sous la forme d’un fichier de points (.xyz), ouvrir ce fichier sous
QGis et l’enregistrer en SHP. Préparer 3 champs : X / Y / Z.
- Sur Qgis, utiliser l’outil « Prélever les valeurs rasters » pour associer à chaque point la valeur du pixel
dans le Raster altimétrique.
- Copier la table attributaire dans une feuille excel et l’organiser en 3 colonnes sans en-têtes :
Colonne 1 Colonne 2 Colonne 3
Coordonnées X Coordonnées Y Altitudes Z
- Si on a appliqué un décalage sur le modèle, penser à appliquer le X+dX et Y+dY aux coordonnées
- Enregistrer le fichier en CSV (séparateur : point-virgule) (*.csv)
Méthode 2 : Utilisation de données sous forme de points topo
Si les points sont en .dwg ou .dxf dans un format classique de bloc points topo (ALT, MAT, …) :
• Dans Autocad : Extraire les points à exporter et les enregistrer en .dxf
• Dans Mapinfo :
- importer le fichier .dxf
- cocher « Préserver les données des attributs »
- dans les options cocher « Stocker Altitude »
- Valider, un fichier .tab est créé (ouvrable ensuite sous QGis)
• Dans QGis :
- Ouvrir le .tab
- Ajouter 2 attributs X et Y pour récupérer les coordonnées de chaque point :
- X = $x et Y = $y (type = réel / longueur = 10 / précision = 3 à 5)
- Copier la table d’attributs
- Voir Préparation du fichier topo CSV
• Sur Excel, créer un CSV comme expliqué dans la préparation de la Topo Raster.
o Application de la donnée topographique dans Fudaa :
- Maillage -> Nœud.
- Faire la sélection des nœuds souhaités (ils se mettent en
surbrillance bleue)
- Clic droit -> « Initialiser depuis des données géographiques » ->
cocher « Modifier à partir de fichiers » puis « Suivant »
- « Ajouter » -> sélectionner le fichier CSV préparé -> Ouvrir
- « Charger » (les composantes autres que X/Y s’affichent) -> Suivant
-
- Choisir Z pour la bathymétrie et la
surface libre
- « Suivant » puis « Terminer »
Remarque : lorsqu’on souhaite appliquer de la données topographique sur certains nœuds du maillage, il est possible
de faire une sélection des nœuds à partir de couches SIG préalablement chargées dans Fudaa. Dans ce cas, ne pas
oublier de cocher l’option « Modifier les objets sélectionnés uniquement » au risque d’écraser les données topo déjà
appliquées sur le maillage.
o Extraction des éléments du maillage avec fudaa
Parfois, il peut être utile d’afficher les éléments de son maillage sous QGis
(contour, éléments, nœuds).
BlueKenue permet d’enregistrer le maillage sous forme de points pouvant
ensuite via plusieurs manips être transformés en SHP sous QGIS, mais il ne
permet pas d’enregistrer le contour ou les éléments du maillage en shp
directement sans passer par d’autres outils (MapInfo par exemple).
Avec Fudaa :
- Dans le visualisateur de maillage, « Fichier »
-> « Exporter » -> « Maillage »
- Choisir le nom du fichier et le type SHP
- Cocher la/les donnée(s) à extraire
• Appliquer les coefficients de frottement de fond sur le maillage
o Possibilité 1 : On n’a pas beaucoup de secteurs de Ks
différents dans le modèle.
Dans ce cas on peut utiliser une couche de zones qui permet de
sélectionner successivement les nœuds où l’on veut appliquer tel
ou tel coefficient de frottement.
Ensuite :
- Sélection de nœuds souhaités
- Clic droit sur « Nœuds » et « Editer »
- Cocher « valeurs agrégées » pour changer toutes
les valeurs à la fois (et non nœud par nœud)
On peut utiliser d’autres variables contenues dans le maillage pour
calculer des valeurs à appliquer via « f(x).
• Les Conditions aux limites
Une fois les conditions aux limites appliquées, on peut passer aux conditions initiales :
- Le débit sur les limites « amont »
- La hauteur d’eau sur les limites « aval »
o Sélectionner les limites amont/aval :
Sélectionner les nœuds sur lesquels on veut placer une frontière liquide.
Clic droit -> insérer des nœuds liquides
Choisir Hauteur imposée ou Débit imposé
Pour les débits, dans la valeur de débit renseigner un numéro en incrémentant de 1 à chaque nouvelle frontière.
Cela va permettre de connaitre la numérotation des frontières par le modèle (plus tard, selon ordre affiché dans le
fichier CAS)
Renseigner une courbe de tarage :
Remarque : les deux courbes de tarage et du débit doivent être jointe au dossier où se trouve le fichier CAS, puis il
faut enlever les extensions (.txt) et les faire appeler sur le fichier cas.