Welcome to the forums. Please post in English or French.

You are not logged in.

#1 2022-03-04 17:42:17

Groguiguy
Member
From: Centrale Lille
Registered: 2017-11-07
Posts: 94

Force unique de réaction sur une face

Bonjour à tous,

Je souhaiterai obtenir l’effort de réaction orthogonal à une face contrainte en DDL, l'effort globale mesuré dans des mords d'un essai de traction simulé par exemple.
J'ai déjà lu ce sujet (code-aster.org/forum2/viewtopic.php?id=19987) et le livre de Jean-Pierre Aubry qui mentionne l'utilisation de :

#U4.81.21
sum_reac=POST_RELEVE_T (
ACTION=_F (
INTITULE=’sum reactions’ ,
GROUP_NO = ( ’groundS’ , ’groundN’ , ) ,
RESULTAT=stat ,
NOM_CHAM=’REAC_NODA’ ,
TOUT_ORDRE=’OUI’ ,
RESULTANTE = ( ’DX’ , ’DY’ , ’DZ’ , ) ,
MOMENT = ( ’DRX’ , ’DRY’ , ’DRZ’ , ) ,
POINT = ( 0 , 0 , 0 , ) ,
OPERATION=’EXTRACTION’ ,
               ) ,
) ;
#U4.91.03
IMPR_TABLE ( TABLE=sum_reac , FORMAT_R=’1PE12.3’ , )

Cependant cette méthode sort une valeur de force de réaction pour chaque noeud de la face sélectionnée.
Hors, la résultante n'est pas la somme de ces efforts, notamment parce que des efforts vont être négatifs et positifs dans un même élément (différence entre DLOAD et DSLOAD de mémoire).
Si j'en crois la doc u4/u4.81.04 il faut appliquer TOUT=OUI à la définition du groupe observé, mais dans ce cas je ne comprends pas la différence entre Reac_Noda et Force_Noda puisque le calcul est effectué partout.

J'ai forcément mal compris quelque chose!
En vous remerciant!

Offline

#2 2022-03-22 17:28:25

Groguiguy
Member
From: Centrale Lille
Registered: 2017-11-07
Posts: 94

Re: Force unique de réaction sur une face

Bonjour à tous,

Je n'ai toujours pas trouvé de solution pour calculer l'effort de réaction sur une face en déplacement imposée.
Dois je prévoir un programme Python prenant l'effort à chaque nœud et la pondération selon les lois de la MEF au sujet des efforts distribués ?
Il n'y a pas d'option de Code_Aster pour une réponse sur toute la face ?

En vous remerciant!

Offline

#3 2022-03-23 09:30:47

jeanpierreaubry
Guru
From: nantes (france)
Registered: 2009-03-12
Posts: 3,917

Re: Force unique de réaction sur une face

bonjour

il doit y avoir une erreur quelque part, tellement grosse qu'elle passe inaperçue (ça arrive souvent !)
car je procède régulièrement de cette façon et j'obtiens le résultat sur la face
par exemple, pour une face comportant 1681 noeuds :

NOMBRE  DE GROUPES DE NOEUDS CREES :      1
=======================================

               ------------------------------------------------------
               !         NOM DU GROUPE         ! NBRE DE NOEUDS DU  !
               !            NOEUDS             !      GROUPE_NO     !
               ------------------------------------------------------
               !  mesh_face_xy_0               !      1681          !
               ------------------------------------------------------

##ASTER 15.02.00 CONCEPT 0000001c CALCULE LE 23/03/2022 A 08:20:06 DE TYPE        
#TABLE_SDASTER                                                                   
INTITULE         RESU     NOM_CHAM          INST          RESULT_X      RESULT_Y      RESULT_Z      MOMENT_X      MOMENT_Y      MOMENT_Z    
sum reactions    0000001b REAC_NODA            0.000E+00    -3.464E-13     9.780E-14     5.252E+03    -6.238E-11     7.887E-11    -5.080E-13

consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/

Offline

#4 2022-04-01 15:38:32

Groguiguy
Member
From: Centrale Lille
Registered: 2017-11-07
Posts: 94

Re: Force unique de réaction sur une face

C'est bon, j'ai fini par y arriver.
Effectivement il y avait bien une erreur dans mon .com pour ne pas sortir convenablement les valeurs, mais lorsque j'eus corrigé l'erreur j'avais alors le classique :

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee.      !
   ! Les bases sont fermees.                                    !
   ! Type de l'exception : error                                !
   !                                                            !
   !  le calcul de l'option :  FORC_NODA                        !
   !  n'est possible pour aucun des types d'éléments du LIGREL. !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Même en définissant convenablement CREA_GROUP_NO=_F( GROUP_MA='nom_du_groupe_de_maille_cible',),

Le problème est ailleurs :
Il faut dans CALC_CHAMP définir un groupe d'éléments volumiques, et non 2D. C'est trompeur sans avoir créé de groupe 3D dans l'outil de maillage Aster Study nous propose naturellement les groupes 2D qui servent à la définition des conditions limites.

Offline

#5 2022-04-01 16:05:12

jeanpierreaubry
Guru
From: nantes (france)
Registered: 2009-03-12
Posts: 3,917

Re: Force unique de réaction sur une face

serait ce un de méfaits de la définition avec TOUT='OUI' dans les défintions de groupes, de matériau, etc..
je mets en garde contre cela au chapitre 20.12 de mon livre
et je ne l'utilise plus jamais en tout cas


consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/

Offline