Assistance #335
Lois de comportement ORTHOELA2D_C et fonctions nD
Description
Bonjour Gérard,
Dans mon modèle axi de BSO je cherche à faire dépendre les coefficients d'une loi ORTHOELA2D_C d'une fonction nD pour obtenir un comportement non linéaire de type ISO_ELAS_ESPO1D (modèle non linéaire des rubans proposé par Julien).
J'y arrive dans le cas d'un cylindre (en axi) et avec une fonction nD de la forme :
- Fonction raideur ruban
f_E_ruban FONCTION_COURBE1D
courbe= f_eps_loi1D
un_argument= EPS22
fin_fonction_courbe1D_
où f_eps_loi1D et la loi 1D non linéaire identifiée par Julien.
La fonction nD intervient dans la loi de comportement :
loi_film ORTHOELA2D_C
E1= E1_fonction_nD: f_E_ruban
.....
J'ai besoin de complexifier cette dépendance pour 2 raisons :
1. Le repère d'orthotropie (repère local de l'élément) ne coïncide pas avec le repère global du tenseur EPS, il faudrait remplacer la variable EPS22 par Almansi_local11. Cela n'est pas accepté par Herezh, pourquoi ?
2. La proportion de ruban par rapport au film dans la largeur du fuseau est variable.
J'ai essayé de construire une unique fonction nD f_E_ruban de type FONCTION_EXPRESSION_LITTERALE_nD. La dépendance avec les composantes du vecteurs normal N_surf_1,... est possible mais les composantes globales et/ou locales du tenseur des déformations (Almansi) ne sont pas acceptées. Est-ce une limitation actuelle de Herezh ou une impossibilité de principe ?
Est-ce que la grandeur Almansi_local11 est bien la composante locale du tenseur d'Almansi ? Dans ce cas elle devrait dans mon cas coïncider avec EPS22, ce qui n'est pas le cas (voir le fichier .Maple).
Ci-joint mon cas test élémentaire à 1 élément axi.
Merci d'éclairer ma lanterne !
Frank
Fichiers