Assistance #177
Contact bielle sur quadrangle
Description
Gérard,
J'ai écrit un petit test élémentaire pour tester le contact entre des éléments 1D (tendon) et des éléments 2D (film).
L'animation Gmsh montre que le contact a lieu entre les incréments 8 et 9 puis les bielles passent au travers des faces des quadrangles.
Le fichier log montre que F_N max 0. J'ai fait une demande de sortie CVisu des grandeurs de contact et je n'ai rien. Je ne comprends pas pourquoi le contact n'est pas activé. Mon jeu de données est-il incomplet ?
Merci pour ton support.
Frank
Fichiers
Mis à jour par Frank Petitjean il y a environ 6 ans
- Fichier m_cylindre01.her m_cylindre01.her ajouté
- Fichier m_tendon01.her m_tendon01.her ajouté
- Fichier modele.CVisu modele.CVisu ajouté
- Fichier modele.info modele.info ajouté
Avec les fichiers...
Mis à jour par Gérard Rio il y a environ 6 ans
- Statut changé de Nouveau à Résolu
- % réalisé changé de 0 à 100
Lorsqu'il y a plusieurs maillages dans le calcul, les maillages sont organisés de la manière suivante:
- maillage(s) esclaves
- maillage(s) mixtes (peuvent-être esclaves et maîtres à la fois)
- maillage(s) maîtres
Le contact s'étudie entre les noeuds esclaves et les frontières du maître.
Dans l'exemple, nous avons 2 maillages et on indique qu'il y a deux maillages esclaves. Donc finalement il n'y a pas de maître, d'où pas de contact a étudier.
exécution: Herezh affiche qu'il y a 2 mailles esclaves et 0 maillage maître
Modifs effectués:
------------------
De manière à avoir au final:
- le tendon = le maillage esclave
- la toile = le maillage maître
1) on change l'ordre d'apparition des maillages
2) on indique un seul maillage esclave
3) pour la visualisation: une sortie tous les 100 itérations
4) le contact introduit des raideurs supplémentaires qui réduisent la stabilité, on augmente le lambda: 4 -> 8
-> on observe du contact et une déformée qui commence à ressembler à quelque chose
Mis à jour par Frank Petitjean il y a environ 6 ans
Tout est dans la doc Herezh que j'ai pourtant lu mais bon...
Merci Gérard pour ces explications. Le contact a en effet lieu maintenant. Il me reste à ajuster les paramètres pour obtenir quelque chose de correcte.
Frank
Mis à jour par Frank Petitjean il y a environ 6 ans
- Fichier m_cylindre01.her m_cylindre01.her ajouté
- Fichier m_tendon01.her m_tendon01.her ajouté
- Fichier modele.CVisu modele.CVisu ajouté
- Fichier modele.info modele.info ajouté
Gérard,
J'ai affiné le test et essayé de "jouer" avec les paramètres sans réel succès : le contact n'est pas bon dans cet exemple élémentaire mais représentatif de mon problème.
1. Le fichier log indique bien que le contact a lieu à partir de l'incrément 4 (F_N max > 0) et se poursuit ensuite. Pourtant les sorties Gmsh (FORCE_CONTACT) et Maple affiche des valeurs non nulles à l'incr 4 puis 0. Le contact est relâché et la bielle passe au travers !
2. Le paramètre AJOUT_EXTRA_BOITE_PRELOCALISATION n'est pas dans la doc or il joue un rôle important. Avec la valeur 0.1 que tu as proposée le calcul ne converge pas.
Il me faut absolument maitriser ce test élémentaire mais représentatif avant de passer à un BPL avec ses tendons glissants !
Je suis en poste tous les jours.
Frank
Mis à jour par Gérard Rio il y a environ 6 ans
des éléments de réponses:
- les facettes en contact sont orientées. Dans le test ici, a priori il faut inverser l'orientation des facettes de manière à avoir un contact sur les noeuds esclaves du centre et non des bords,
- pour un contact tendon-surface, il est préférable de définir une zone de contact de type "glue_contact" (ou "glue_contact_init_gap_zero") ce qui aura comme conséquence d'imposer systématiquement le contact. De plus il est sans doute préférable d'utiliser un contact avec "TYPE_PENALISATION_PENETRATION 7 "
- pour étudier l'influence du contact lorsque l'on ne comprend pas ce qui se passe, il est intéressant d'utiliser momentanément les options suivantes:
1) "mode_debug_= 1" # dans le fichier .info -> sortie des infos de post-traitement à chaque itération
2) "FORCE_GENE_EXT FORCE_GENE_INT" en post-traitement -> sortie des forces externes à chaque maillage, => les forces de contact pour les zones en contact (visualiser les isovaleurs sous forme de vecteurs aux noeuds).
3) "NIVEAU_COMMENTAIRE_CONTACT " dans le .info -> donne des infos sur les éléments de contact en cours et en recherche.
Mis à jour par Gérard Rio il y a environ 6 ans
j'ai fait une erreur dans mes remarques: c'est le type de contact 8 au lieu du type 7 (indiqué) qu'il faut utiliser !!