Anomalie #173
Problème contact 3D, test élémentaire
Description
Gérard,
J'ai constitué un cas test beaucoup plus simple et pour lequel le contact n'a pas lieu.
Il s'agit d'une portion de cylindre mince sous pression (20 quadrangles ) que je viens écraser avec la face d'un hexaèdre.
Chargement 1 : mise en pression
Chargement 2 : descente du cube en 4 incréments de charge.
Le test dure moins de 1 seconde !
Pourquoi est-ce que le contact n'est pas actif ?
Frank
Fichiers
Mis à jour par Frank Petitjean il y a plus de 6 ans
- Fichier cube.her cube.her ajouté
- Fichier cylindre.her cylindre.her ajouté
- Fichier modele.CVisu modele.CVisu ajouté
- Fichier modele.info modele.info ajouté
avec pièces jointes !
Mis à jour par Gérard Rio il y a plus de 6 ans
- Fichier modele.info modele.info ajouté
- Statut changé de Nouveau à En cours
- % réalisé changé de 0 à 80
Le cas de contact proposé ne peut pas fonctionner en l'état.
1) il manque des paramètres de contrôles qui permettent de préciser la recherche du contact. Je joins un fichier avec des paramètres qui fonctionnent. Il existe surement d'autres jeux de paramètres possibles et le jeu proposé n'est sans doute pas optimisé.
2) par défaut (cf. doc) le paramètre TYPE_PENALISATION_TANGENTIELLE est 2, ce qui signifie que la pénalisation est fonction de la raideur des éléments maîtres. Dans l'exemple fourni la cinématique du maître est entièrement fixée d'où l'utilisation d'une loi rien3D ce qui est correct, mais qui évidemment ne fournit pas de raideur donc le calcul de la pénalisation n'est pas possible avec le type 2. Par contre en utilisant le type 1 c'est ok.
Mis à jour par Frank Petitjean il y a plus de 6 ans
- Fichier modele.info modele.info ajouté
- Fichier modele.log modele.log ajouté
- Fichier modele_deplace_Gmsh.pos modele_deplace_Gmsh.pos ajouté
Bonjour Gérard,
Quelle version de Herezh utilises-tu ? Avec la dernière publiée, la V6.839, le test que tu m'as envoyé ne fonctionne pas : lorsque je visualise la déformée le prisme traverse le cylindre sans entrer en contact !!!
Je te joins les fichier log et Gmsh
Mis à jour par Gérard Rio il y a plus de 6 ans
OK, effectivement j'utilise la version 6.847 ...
bon, normalement je devrais la mettre en ligne prochainement. Je suis en train de finir les validations des tests existants.
affaire à suivre !
Mis à jour par Frank Petitjean il y a plus de 6 ans
Bonjour Gérard,
Merci d'avoir publié la version 6.847. J'ai pu vérifier que mon cas test était maintenant correctement calculé. Par contre je n'arrive toujours pas à traiter le cas du BSO sous largueur (le premier modèle fourni), même avec des hexaèdre comme faces de contact. Avec les paramètres proposée pour le cas test les nœuds du maillage esclave (fuseau) sont bien entrainés lors du contact mais ils rentrent très largement dans le corps de contact.
J'ai augmenté la PENALISATION_PENETRATION mais le calcul explose si ce paramètre devient trop élevé. J'ai essayé de modifier d'autres paramètres mais sans succès.
Pour l'étude que je dois rendre... après-demain je vais conserver le contact type 1 qui fonctionne avec la version 6.812 (mais je ne bénéficie pas du calcul d'intégrale pour le calcul du volume).
Il va falloir que je travaille sérieusement ces problèmes de contact car c'est au cœur du lot 5 de la R&T. C'est pas gagné !
Bon we
Frank