Herezh++

Le problème n'est pas finalement du côté d'herezh. En fait, il est dû à une mauvaise constitution du maillage.
Globalement vous avez pour la plupart des éléments, deux éléments superposés à la place d'un.
Avant d'étudier le contact, le programme détermine les noeuds frontière des maillages esclaves, et c'est seulement ces noeuds esclaves qui vont ensuite être testés pour le contact. Pour déterminer si un noeud est frontière ou non, la technique retenue dans Herezh est d'étudier si le noeud appartient à une face qui est une frontière commune à deux éléments, si oui ce n'est pas un noeud frontière.
Dans le cas de deux éléments facettes superposés, la face de la facette est commune à deux éléments donc aucun des noeuds de la facette n'est considéré comme noeud frontière !
Donc voici ce que j'ai fait pour mener votre calcul à terme, correctement:
1) suppression des éléments superposés à l'aide de l'utilitaire d'Herezh,
2) j'ai utilisé un facteur de pénalisation de 0.01 et un lambda de 1.
3) j'ai utilisé un pas de temps de 1., donc le calcul est effectué en 1 incrément,
--> on obtient une déformée correcte à l'autocontact près qui n'est pas géré ici.
@+

Bonjour,

J'ai bien essayé ces solutions mais cela me conduit à un calcul sans fin, de plus j'ai eu un problème avec le fichier CVisu je ne peux donc pas voir ce qui se passe.

Je ne comprends pas bien d'où peut venir le problème...

Adrien

Il faut penser à modifier le lambda, la pénalisation etc.
bon... voici les fichiers qui fonctionnent avec moi ce qui vous permettra de démarrer, maintenant à vous de jouer pour la suite !