Discussion :

[phy4me]

Bonjour,
je cherche un mailleur qui génère un maillage triangulaire non structuré; j'ai trouvé pas mal de mailleur écrits en C++, mais puisque je suis expeeert en C, j'ai pas pu les tourner . Svp si quelqu'un en connaît un, je lui serais reconnaissant de me diriger.

Merci

[FR119492]

Salut.

Si je t'ai bien compris, tu cherches à faire un préprocesseur pour un programme d'éléments finis. Dans ce cas, on peut envisager deux approches distinctes:

La première est utilisée dans un grand nombre de programmes du commerce. Elle se base sur la méthode de discrétisation de Delaunay: le contour du domaine est approximé par une ligne polygonale fermée. On choisit le moins mal possible deux sommets éloignés que l'on relie par une autre ligne polygonale (ouverte), qui peut parfaitement être constituée de segments alignés, définissant ainsi deux sous-domaines. On recommence avec chacun de ces sous-domaines, et ainsi de suite jusqu'à ce que tous les sous-domaines soient des triangles. Ensuite, en général, on renumérote les noeuds, afin que la matrice du système ait la forme d'une bande aussi étroite que possible. En principe, ça marche toujours et c'est utilisable par quelqu'un qui n'a que peu d'idées sur la méthode. Mais c'est très loin d'être optimal du point de vue du temps d'exécution.

C'est pourquoi, personnellement, j'ai opté pour une autre approche: dans tous les problèmes industriels que j'ai eu à résoudre, la plus grande partie du domaine avait une géométrie simple et pouvait être maillée de manière systématique; il ne restait donc que peu à mailler à la main.

J'ai écrit entre 65 et 69 un programme basé sur cette méthode et qui permettait jusqu'à 3000 noeuds sur une machine ne disposant que de 128k de mémoire vive. Ce programme n'existe évidemment plus. Par la suite, il y a eu diverses versions entre 74 et 87, puis encore d'autres par la suite au gré de mes pérégrinations professionnelles. Actuellement, il n'y a aucune version opérationnelle, mais je ne désespère pas de réactiver ce projet quand j'en trouverai le temps, dans l'idée de le mettre ensuite à disposition sur le site. Je ne peux donc rien t'offrir pour le moment, mais, à long terme, pourquoi pas?

Jean-Marc Blanc

[phy4me]

Bonjour,
je te remercie de m'avoir repondu, et pour les éxplications données sur la méthode Delauny

[feynman]

Salut,
Si tu veux utiliser la methode des elements finis, avec un maillage ainsi que la solution, n utilise pas le fortran mais je te conseil d utiliser Freefem++

[phy4me]

Merci beaucoup pour l'info,
un mailleur n'utilise ni les Elémenet Finis ni Volume Finis ... c'est une étape préliminaire dont l'objectif est de préparer le domaine (le mailler ) sur lequel on veut résoudre nos équations que ce soit par EF ou VF (que j'utilise). donc ça n'a aucune relation avec les EDP

Merci pour FreeFem++

[feynman]

Mais , avec Freefem , tu mailles ton domaine avant de resoudre ton EDP!
et tu voit ton domaine plain avec les triangles avant la resolution,
Va voire les exemples sur http://www.freefem.org/ff++/ftp/manual.pdf