Herezh++

oui,
effectivement l'intégrale s'effectue sur le volume de l'élément support, c-a-d ligne pour les éléments 1D, surface pour 2D et volume pour 3D.
Maintenant si on veut le volume interne, par exemple entre le plan xy et la surface il faut utiliser une fonction nD qui a pour argument la coordonnée Z du point et qui donne comme résultat Z

C'est ce que tu as appelé dans ton fichier test "le moment geo / y" ?

J'ai créé dans l'ordre logique :

la fonction X1

fct_ddl_X1 FONCTION_EXPRESSION_LITTERALE_nD
deb_list_var_ X1 fin_list_var_
fct= X1
fin_parametres_fonction_expression_litterale_

l'intégrale de volume

integrale_sur_volume_   #--------# 
E_tout  une_fonction_nD_ fct_ddl_X1

et la fonction nD

fct_volume_2D_plan_zy FONCTION_EXPRESSION_LITTERALE_nD
deb_list_var_ int_vol_E_to_fct_nD_fct_ddl_X1 fin_list_var_
fct= int_vol_E_to_fct_nD_fct_ddl_X1
fin_parametres_fonction_expression_litterale_

J'utilise cette fonction dans un chargement pression bidon pour récupérer la valeur

F_to PRESSION 1 Fonction_nD_CHARGE: fct_volume_2D_plan_zy

et j'obtiens à chaque nœuds 1 (voir le fichier Gmsh modele_pression_ext_ddl_etendu_noe_Gmsh), ce qui n'est pas la valeur attendue.

OK pour la fonction nD:
> la plaque fait 10x1 et les x valent 10, donc l'espace entre la plaque et le plan yz est 10 x 10 = 100
> on lit dans le fichier .BI que l'intégrale vaut 99.99999... donc 100, -> donc a priori on obtient la bonne valeur de l'intégrale

Concernant l'utilisation de l'intégrale, il y a une erreur de mise en donnée:
-> la variable globale générée par l'intégrale se nomme : "int_vol_E_tout_ddl_etendu_TEMP" et non "int_vol_E_to_ddl_etendu_TEMP" !! Actuellement, Herezh ne signale pas le pb (je vais essayer de modifier cela dans le futur).
Deux solutions pour avoir le nom des variables généré: soit regarder dans le .BI (pour les utilisateurs curieux et averti), ou simplement lorsqu'herezh est en interactif, on peut avoir le nom des variables globale (cas 9).
Cependant, avec un delta_t de 1, le fonctionnement ne va pas être celui que tu attends car la fonction va utiliser la valeur à t, or l'intégrale à t = 0 !!
Donc par exemple une solution est d'utiliser plusieurs incrément de temps, et ici dès le second incrément l'intégrale utilisée sera correcte et la pression aussi.

Bonjour Frank,

1) Le calcul du volume s'effectue de la manière suivante:
- la surface de l'élément est décomposée en surfaces élémentaires triangulaires, cette décomposition dépend du type d'élément et d'interpolation,
- ensuite pour chaque triangle le code suivant est exécuté, sachant que A correspond aux coordonnées du point A (premier point du triangle), idem B pour le deuxième point, et C pour le troisième point:

" Coordonnee AB; Coordonnee AC;
Coordonnee G=1./3.*(A+B+C); // coordonnées du centre de gravité
Coordonnee so=0.5*Util::ProdVec_coor(AB,AC);
so(1) *=G(1);so(2) *=G(2);so(3) *=G(3);
volumePlan += so;
"
donc : AB = un coté du triangle, AC = l'autre coté, d'où so = la surface du triangle
G est la position du centre de gravité
d'où so(3) = Z(G) * (la surface projetée sur z => so(3)) c-a-d le volume entre le triangle et le plan xy
et de même pour les deux autres volumes.

Ensuite tous les volumes sont sommés dans une grandeur globale. Le résultat est en relatif, ce qui permet, lorsqu'il s'agit d'une surface fermée d'avoir le volume interne de la surface indépendamment du plan de projection.

2) Variation: dans mon cas j'obtiens un volume de 100. et ensuite si je fais une rotation autour de z, la valeur change ce qui est normal.
Pour 45° j'obtiens -4.75 par rapport à xz et yz
Pour une rotation de 90°, cela donne 0 pour yz et 100 pour xz,

Donc, j'ai l'impression qu'herezh donne les résultats attendus (aux approximations de la facettisation près).

3) Dans le cas du calcul d'une intégrale, elle s'effectue non pas à partir de la facettisation, mais via les points d'intégration utilisés pour le calcul d'équilibre. On peut donc obtenir un résultat un peu différent du précédent, suivant le nombre de points d'intégration utilisé. A priori (à vérifier), l'intégrale via les points d'intégration devrait être (?) un peu plus précise (mais à confirmer).

Bonjour Gérard,

Je viens aux nouvelles concernant l'ajout dans Herezh de la possibilité de récupérer les coordonnées de la normal aux faces d'éléments pour une utilisation dans des fonctions nD.

Merci
Frank