Discussion :

[matt394]

Bonjour,

J'ai un modèle d'éléments finis déjà maillé (je connais les coordonnées des noeuds et les noeuds qui composent les éléments), mais je voudrai pouvoir rajouter des elements dans cette geometrie.

Pour ça, il faut que je coupe en deux certains éléments dejà existant (ceux qui sont sur la trajectoire de la face que je veux rajouter). Quelqu'un a-t-il une idée pour savoir quels sont les éléments que je dois couper? (je ne connais que les coordonnées des sommets).

Merci d'avance

PS : pour info si c'est utile, je code sous MatLab

[Jerome Briot]

(ceux qui sont sur la trajectoire de la face que je veux rajouter).

Tu peux nous donner un peu plus d'informations pour qu'on visualise bien ton problème.

[matt394]

Tu peux nous donner un peu plus d'informations pour qu'on visualise bien ton problème.

Le souci c'est que je ne peux pas poster d'images parce que c'est sur un ordi sans réseau... Alors je vais essayer de vous faire comprendre autrement.

Il faut imaginer 2 plans verticaux et un plan horizontal (dans mon modèle, ce ne sont pas tous des plans verticaux, il y a aussi des surfaces courbes), tous les 3 maillés avec des éléments irréguliers (un mélange de carrés et de triangles). Je sais, c'est pas pratique, mais c'est comme ça... Ce que je voudrai faire c'est couper tout ce beau monde par un plan vertical. Donc je voudrai pouvoir récupérer dans un tableau les numéros des éléments coupés par le plan, pour pouvoir ensuite les couper en deux. En fait ce que je cherche, c'est un méthode permettant de savoir si un élément dont je ne connait que les coordonnées des sommets est coupé par le plan ou pas (sachant que ce plan est lui-même défini par 4 points, j'ai bien dit 4 et pas 3...). Après je m'attaquerai au maillage du nouveau plan, mais c'est la prochaine étape.

[pseudocode]

bah il suffit de projeter tes mesh sur un plan normal a ton plan de coupe.

Ensuite choisir un repere sur ce plan normal, par exemple le repere dans lequel l'axe "X" est l'intersection plan normal / plan de coupe.

Dans ce cas, les mesh "a decouper" sont ceux pour lesquels les coordonnées en "Y" des sommets sont a la fois positives ET negatives.

[matt394]

Oula oui... c'est vraiment pas compliqué en fait... On va mettre ça sur le compte de la fatigue...

Merci beaucoup en tout cas.

[pseudocode]