Ajouté par Frank Petitjean il y a plus de 6 ans. Mis à jour il y a environ 6 ans.
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
| cube.her (4,34 ko) cube.her | Frank Petitjean, 17/06/2018 12:10 | ||
| cylindre.her (7,6 ko) cylindre.her | Frank Petitjean, 17/06/2018 12:10 | ||
| modele.CVisu (7,35 ko) modele.CVisu | Frank Petitjean, 17/06/2018 12:10 | ||
| modele.info (2,81 ko) modele.info | Frank Petitjean, 17/06/2018 12:10 | ||
| modele.info (3,02 ko) modele.info | Gérard Rio, 19/06/2018 18:16 | ||
| modele.info (3,02 ko) modele.info | Frank Petitjean, 21/06/2018 09:46 | ||
| modele.log (160 ko) modele.log | Frank Petitjean, 21/06/2018 09:46 | ||
| modele_deplace_Gmsh.pos (15,9 ko) modele_deplace_Gmsh.pos | Frank Petitjean, 21/06/2018 09:46 |
avec pièces jointes !
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.
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
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 !
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