pre-pos_processing_gmsh-herezh

Bonjour Gérard,

Je suis toujours dans la définition de références pour une pièce importée en STEP.

Je te rapporte 2 choses :

• Je n'arrive pas à récupérer les références associées aux Physical Groups que j'ai introduits dans Gmsh.
Via msh2her puis l'optimisation du maillage.
Les références n'apparaissent pas dans mon .lis
Je pense que c'est à cause de la nature importée de la géométrie ?
Le .geo contient en effet juste un Merge.

• Du coup, j'utilise ton utilitaire "informations avec plus creation_reference"
Et je rencontre un autre souci avec la condition de définition d'un cylindre.
L'utilitaire ne considère pas les valeurs négatives pour les coordonnées d'un point du cylindre.
Je n'arrive donc à référencer qu'un trou sur les 2 (celui qui est dans les X positifs).
Voir l'image dans l'archive (c'est complet) et la capture du shell relativement au problème (en PJ).
Là je suis bloqué, sauf à opérer une translation vers les X positifs de toute la pièce...

NB :
• Pour "attraper" une arête linéaire (ça marche) de noeuds entre 2 points, j'ai utilisé 2 conditions :
- l'inclusion dans un cylindre de faible diamètre (d'axe défini par les 2 points),
+ l'inclusion entre 2 plans // contenant ces 2 points.
-> je me demandais s'il n'y avait pas plus simple ?

Merci pour ton retour,
Christian

Bonjour Christian,

je répond tout d'abord à la première méthode, que je viens de simuler et qui fonctionne.
1) je récupère le .step dans gmsh > ok
2) je crée une physical sur une ligne que je nomme titi. Gmsh me demande si je veux créer un .geo, réponse "Create new '.geo' file "
3) je crée une physical sur tout le volume ->> obligatoire car sinon dans le fichier de sortie il n'y aura "que la ligne titi !!". La raison est que dès que l'on définit une physical, le maillage produit ne contient que les physical !!. Je valide ("e") puis "q" pour sortir des physical
4) je maille (ici par défaut)
5) je sauvegarde au format ascii 2
6) j'utilise msh2her.pl pour récupérer le format her --> DATA.her
7) je visualise le résultat avec: hz_visuMail.pl DATA, on voit qu'il y a bien une référence de ligne: "A_titi"
--> on peut donc ainsi définir des ref de ligne, de surface, de volume

Ensuite pour la mise en données d'Herezh, on aura donc:
- des éléments volumiques tétraédriques qui permettront de calculer les efforts internes, c-a-d ceux correspondants aux actions des contraintes internes. Ils dépendent de la loi de comportement adopté
- des éléments linéiques, pour lesquels on utilisera une loi "LOI_RIEN1D" car ces éléments ne doivent pas intervenir dans la raideur interne. Par contre on pourra leur appliquer des efforts externes ce qui permettra d'imposer les efforts linéiques que l'on veut. On peut ainsi suivant le même modèle, imposer des efforts ponctuels, surfaciques et volumiques. De même on peut imposer des conditions cinématiques sur des ref de noeuds.

Est-ce que cela répond à tes questions ?

a+

re...
suite avec la seconde méthode. --> il y avait un mauvais fct dans Herezh. Je ne dis pas bug car c'est un truc bizarre que j'ai un peu de mal à comprendre.
Au moment de la lecture via un getline, je suis "en général" obligé de supprimer le premier caractère que je lis ce se trouve être un caractère non exploitable, mais dans certain cas il faut que je réintègre ce caractère en particulier si c'est un chiffre. J'ai cherché partout l'explication ... et à part trouvé des gens qui avait la même interrogation que moi, je n'ai pas trouvé d'explication.
Bref, ici dans ton cas (et c'était vrai également pour d'autres questions/réponses dans Herezh) le premier '-' était supprimé car ce n'était pas un chiffre , et cela je ne l'avais pas pris en compte. Du coup, j'ai modifié et maintenant cela fonctionne avec version 6.954 (et suivantes).
Je met à jour les versions mac.
Peux-tu tester pour me dire si c'est ok pour toi ?
Merci d'avance

Gérard,
Le calcul tourne depuis 2h sur mon Mac (le maillage est sans doute trop fin) mais seulement si je déclare une loi élastique pour le matériau.
J'avais lancé initialement le calcul avec une loi élasto-plastique et Herezh ne semblait pas converger (+ message d'erreur).
Je pensais à cet égard utiliser la loi de PRANDTL & REUSS (?).
Je mets le .info en PJ. Merci si tu peux déjà y jeter un oeil !
A+
Christian