Assistance #179
Contact glue_contact
Description
Gérard,
J'ai positionné les noeuds des bielles très proches des faces des quadrangles et j'ai activé le contact glue_contact. Dès la 2e itération les forces de contact sont énormes et repoussent les noeuds des bielles. Le calcul diverge de suite.
Même chose bien sûr avec glue_contact_init_gap_zero puisque les noeuds sont déjà sur les faces.
Avec les paramètres de contact par défaut même chose.
Désolé mais je ne serais pas allé bien loin tout seul ce matin...
Je vais demander à Hugo des exemples.
Frank
Fichiers
Mis à jour par Gérard Rio il y a environ 6 ans
- Fichier modele_5.info modele_5.info ajouté
- Statut changé de Nouveau à Résolu
- % réalisé changé de 0 à 80
Des éléments de réponse:
1) par défaut avec la définition de zone de collage, il y a prise en compte d'une pénalisation de frottement: on remarque que dans l'affichage des forces tangentielles cela donne des grandeurs non nulles et énormes. Pour neutraliser le frottement tangentiel, il faut indiquer une pénalisation nulle par exemple pour le cas traité j'ai mis:
PENALISATION_TANGENTIELLE 0.
TYPE_PENALISATION_TANGENTIELLE 8
TANGENTIELLE_CONTACT_MAXI 0.1
TANGENTIELLE_BORNE_REGULARISATION 0.1
FORCE_TANGENTIELLE_NOEUD_MAXI 100000
2) En fait le type de pénalisation qu'il faut a priori employer c'est le type 8
TYPE_PENALISATION_PENETRATION 8 #2 par def.
3) j'ai également diminué fortement le coefficient de pénalisation sinon les forces transitoires sont énormes et créent des oscillations et finalement une divergence:
PENALISATION_PENETRATION 0.001# 1 #0.1 par def.
4) j'ai passé la précision globale à 1.e-3 (au lieu de 5e-3) histoire de vérifier la bonne convergence
avec ces modifs, cela fonctionne: je mets le fichier qui m'a servi.
Mis à jour par Frank Petitjean il y a environ 6 ans
Merci Gérard pour ce nouveau dépannage. Mes différents cas tests tournent bien maintenant.
Il faudra mettre à jour la doc avec TYPE_PENALISATION_PENETRATION 8 !!
Il y a une situation qui n'est pas correctement gérée : c'est lorsque un noeud esclave est dans le cône normal à deux facettes adjacentes. Dans ce cas le noeud n'est pas projeté sur une des faces.
Ce cas c'est produit plusieurs fois car j'ai des maillages simples l'un généré par stamm pour le 1/4 cylindre et l'autre par Python pour le tendon. Je suis obligé de prévoir un offset.
Dans l'exemple que tu as traité les bielles sont tournées de 5° ( rotation_= 5. 0. 0 en_degre_ )
Je poursuis en essayant de faire glisser les bielles !!!
Frank
Mis à jour par Frank Petitjean il y a environ 6 ans
Complément :
Si les noeud esclaves sont sur les noeuds des éléments maitres cela fonctionne très bien puisque pas de projection...
Mis à jour par Frank Petitjean il y a environ 6 ans
- Fichier modele.CVisu modele.CVisu ajouté
- Fichier modele.info modele.info ajouté
- Fichier plaque01.her plaque01.her ajouté
- Fichier tendon_droit01.her tendon_droit01.her ajouté
Voici le test dans lequel le fuseau démarre à plat.
Je croise les doigt...
Frank