Herezh++

Merci Julien pour ta réponse (et Gérard de le me l'avoir transmise). J'ai également faut des tests sur des formes simples et j'en suis arrivé à la même conclusion pour des profils courbes (genre 1/4 de cercle -> demi sphère). Par contre sur les profils droits comme un cylindre ou un tronc de cône cela ne fonctionne pas avec la formule initiale. Dans le cas d'un cône à 45° le résultat est même nulle mais c'est normal puisque le produit scalaire est nulle !

J'ai donc en effet pu réutiliser les fonctions nD du modèle de couplage volume-pression en axi et cela fonctionne bien. Un regret concernant ce modèle axi : la loi d'orthotropie utilisée pour intégrer le comportement uniaxial des rubans rend le calcul très très lent, ce qui limite fortement l'intérêt de passer en axi. Par contre plus de problème de distorsion de maillage même dans le cas sol !

Frank

Frank, pourrais-tu expliquer ou transmettre un fichier test sur les cas qui ne fonctionnent pas, car je ne comprends pas bien ce qui se passe et je suis intéressé pour regarder.
Pour la rapidité de la loi d'orthotropie en axi, là aussi je suis preneur d'un test de base...

L'exemple de Julien sur le cylindre m'a fait comprendre le problème : les profils axi que j'ai testés ne sont pas fermés ! En fermant mes figures les calculs sont correctes (cylindre et cône). Pour ma demi-sphère le profil était naturellement fermé (stamm). Il faut donc que je fasse de même pour mes ballons. Ils sont quasi fermés mais il doit quand même y avoir une petite erreur...
Concernant le sens des normales cela change uniquement le signe du résultat. Il faut juste s'assurer qu'elles sont toutes orientées dans le même sens.

Grand merci pour votre aide à tous les 2.
Gérard, Le ticket peut être fermé.

Frank

merci pour cet épilogue. Effectivement, je me suis trompé pour l'histoire du produit scalaire. X.N n'est pas du tout nul dans le cas du cone à 45°, oups...

Je réponds à Gérard concernant ma remarque sur la lenteur de la loi PROJECTION_ANISOTROPE_3D. Je constate que le temps par itération pour un même maillage dépend très fortement des difficultés de convergence. Or ma loi est complexe car elle contient :

TYPE_DE_CRITERE_ PLISSEMENT_MEMBRANE
...
LOI_CONTRAINTES_PLANES
...
LOI_ADDITIVE_EN_SIGMA

ISOHYPER3DFAVIER3    
      +
        PROJECTION_ANISOTROPE_3D sur loi ISOELAS non linéaire par fct nD 
...

Je pense que le temps passé se trouve réparti dans l'ensemble de ces biques.

Frank