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

You are not logged in.

#1 2008-08-07 17:51:00

jo
Member
From: Ingenieurbüro für Mechanik, DE
Registered: 2008-01-29
Posts: 146
Website

[RESOLU]How can I display top-and bot-stress of 2D-elements in Salome?

Bonjour,

I make a lin static analysis with linear 4-noded DKT-elements.

My impression is (but am not shure) that

1)  for fields of  type *_ELGA_*, only the middle-plane-stresses  can  be grafically displayed in Salome. Is that true?

2)  fields of the type *_ELNO_* cannot be grafically displayed in Salome. True?

Thus, I use to calculate the nodal stress field *_NOE_* and this I can display in Salome. In order to get the bottom- and the top-2D-stress, I make twice CALC_ELEM and CALC_NO (first with NIVE_COUCHE='INF', then with NIVE_COUCHE='SUP'). I try to write the results in 2 steps into the same .resu.med-file. First I write the fields: NOM_CHAM=('DEPL','SIGM_NOEU_DEPL','EQUI_NOEU_SIGM',) where the stresses are bottom-stresses. Then I write the fields: NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM',) where the stresses are top-stresses. But Aster ends with a fatal error, it seems that I should not write a certain field twice.

Question 3) What is the easiest way to postprocess top- AND bottom-stresses in Salome? Must I write in 2 different .resu.med-files or is there an easier way?

Thank You for any help. I list below part of the .mess-file
Regards, Johannes


    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0001            CONCEPT DE TYPE :
    #  -------------                  -----------------
   DEBUT(RESERVE_CPU=_F(BORNE=900),
        PAR_LOT='OUI',
        IMPR_MACRO='NON',
        );

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0002            CONCEPT DE TYPE :
    #  -------------                  -----------------
    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0003            CONCEPT DE TYPE : maillage_sdaster
    #  -------------                  -----------------
   MeshLin=LIRE_MAILLAGE(INFO=2,
                        VERI_MAIL=_F(APLAT=1.E-3,
                                     VERIF='OUI'),
                        UNITE=20,
                        FORMAT='ASTER',
                        );

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0004            CONCEPT DE TYPE : maillage_sdaster
    #  -------------                  -----------------
   MeshLin=DEFI_GROUP(reuse = MeshLin,
                     MAILLAGE=MeshLin,
                     CREA_GROUP_NO=_F(GROUP_MA='PLATE_T1',
                                      CRIT_NOEUD='TOUS'),
                     ALARME='OUI',
                     );



    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0005            CONCEPT DE TYPE : modele_sdaster
    #  -------------                  -----------------
   FEM=AFFE_MODELE(MAILLAGE=MeshLin,
                  INFO=1,
                  AFFE=_F(PHENOMENE='MECANIQUE',
                          TOUT='OUI',
                          MODELISATION='DKT'),
                  VERIF='MAILLE',
                  );


SUR LES            3 MAILLES DU MAILLAGE MeshLin
    ON A DEMANDE L'AFFECTATION DE            3
    ON A PU EN AFFECTER                      3

    MODELISATION      ELEMENT FINI      TYPE MAILLE          NOMBRE
                      MEBODKT           SEG2                      2
    DKT               MEDKQU4           QUAD4                     1
 
    #  FIN COMMANDE NO : 0005   DUREE TOTALE:        0.01s (SYST:        0.01s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0006            CONCEPT DE TYPE : mater_sdaster
    #  -------------                  -----------------
   elastic=DEFI_MATERIAU(INFO=1,
                        ELAS=_F(ALPHA=9.9999999999999995E-07,
                                E=2.E5,
                                NU=0.0,
                                RHO=5.0000000000000002E-05),
                        );


MATERIAU : elastic
RELATION DE COMPORTEMENT: ELAS           


 
    #  FIN COMMANDE NO : 0006   DUREE TOTALE:        0.01s (SYST:        0.00s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0007            CONCEPT DE TYPE : cham_mater
    #  -------------                  -----------------
   Mat=AFFE_MATERIAU(VARC_M_ZIRC=_F(),
                    INFO=1,
                    VARC_TEMP=_F(),
                    MAILLAGE=MeshLin,
                    VARC_HYDR=_F(),
                    VARC_CORR=_F(),
                    VARC_EPSA=_F(),
                    VARC_IRRA=_F(),
                    VARC_NEUT1=_F(),
                    VARC_NEUT2=_F(),
                    AFFE=_F(MATER=elastic,
                            TOUT='OUI'),
                    VARC_M_ACIER=_F(),
                    VARC_SECH=_F(),
                    );

 
    #  FIN COMMANDE NO : 0007   DUREE TOTALE:        0.00s (SYST:        0.00s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0008            CONCEPT DE TYPE : cara_elem
    #  -------------                  -----------------
   Cara=AFFE_CARA_ELEM(INFO=1,
                      COQUE=_F(MODI_METRIQUE='NON',
                               EPAIS=1.0,
                               GROUP_MA='PLATE_T1',
                               COEF_RIGI_DRZ=1.0000000000000001E-05,
                               COQUE_NCOU=1),
                      MODELE=FEM,
                      );


     LE MODELE FEM      CONTIENT :
                                        3 ELEMENT(S) COQUE
 
    #  FIN COMMANDE NO : 0008   DUREE TOTALE:        0.01s (SYST:        0.00s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0009            CONCEPT DE TYPE : char_meca
    #  -------------                  -----------------
   fixed=AFFE_CHAR_MECA(INFO=1,
                       VERI_NORM='OUI',
                       DDL_IMPO=_F(LIAISON='ENCASTRE',
                                   GROUP_MA='FIXED_EL'),
                       LIAISON_XFEM='NON',
                       MODELE=FEM,
                       );

 
    #  FIN COMMANDE NO : 0009   DUREE TOTALE:        0.03s (SYST:        0.01s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0010            CONCEPT DE TYPE : char_meca
    #  -------------                  -----------------
   Tensio=AFFE_CHAR_MECA(FORCE_ARETE=_F(GROUP_MA='LOADED_E',
                                       FZ=50.0),
                        INFO=1,
                        VERI_NORM='OUI',
                        LIAISON_XFEM='NON',
                        MODELE=FEM,
                        );

 
    #  FIN COMMANDE NO : 0010   DUREE TOTALE:        0.01s (SYST:        0.00s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0011            CONCEPT DE TYPE : char_meca
    #  -------------                  -----------------
   Traver=AFFE_CHAR_MECA(INFO=2,
                        FORCE_ARETE=_F(FY=-50.0,
                                       GROUP_MA='LOADED_E'),
                        VERI_NORM='OUI',
                        LIAISON_XFEM='NON',
                        MODELE=FEM,
                        );

 
    #  FIN COMMANDE NO : 0011   DUREE TOTALE:        0.00s (SYST:        0.00s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0012            CONCEPT DE TYPE : evol_elas
    #  -------------                  -----------------
   Res=MECA_STATIQUE(EXCIT=(_F(TYPE_CHARGE='FIXE',
                              CHARGE=fixed),
                           _F(TYPE_CHARGE='FIXE',
                              CHARGE=Traver)),
                    INFO=1,
                    OPTION='SIEF_ELGA_DEPL',
                    SOLVEUR=_F(RENUM='METIS',
                               STOP_SINGULIER='OUI',
                               METHODE='MULT_FRONT',
                               NPREC=8),
                    INST=0.0,
                    CARA_ELEM=Cara,
                    MODELE=FEM,
                    CHAM_MATER=Mat,
                    );

--- NOMBRE TOTAL DE NOEUDS :  28 DONT :  24 NOEUDS "LAGRANGE"
--- NOMBRE TOTAL D'EQUATIONS :  48
--- NOMBRE DE COEFFICIENTS NON NULS DANS LA MATRICE :  480
   CHAMP STOCKE :             DEPL INSTANT :  0.00000E+00  NUMERO D'ORDRE :     1
 
    #  FIN COMMANDE NO : 0012   DUREE TOTALE:        0.22s (SYST:        0.17s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0013            CONCEPT DE TYPE : evol_elas
    #  -------------                  -----------------
   Res=CALC_ELEM(reuse = Res,
                REPE_COQUE=_F(NUME_COUCHE=1,
                              NIVE_COUCHE='INF',
                              ANGLE=0,
                              PLAN='MAIL'),
                INFO=1,
                OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
                TYPE_OPTION='TOUTES',
                PRECISION=1.E-3,
                RESULTAT=Res,
                CRITERE='RELATIF',
                CARA_ELEM=Cara,
                );

 
    #  FIN COMMANDE NO : 0013   DUREE TOTALE:        0.05s (SYST:        0.01s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0014            CONCEPT DE TYPE : evol_elas
    #  -------------                  -----------------
   Res=CALC_NO(reuse = Res,
              CRITERE='RELATIF',
              RESULTAT=Res,
              OPTION=('EQUI_NOEU_SIGM','SIGM_NOEU_DEPL',),
              PRECISION=1.E-3,
              );

 
    #  FIN COMMANDE NO : 0014   DUREE TOTALE:        0.01s (SYST:        0.01s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0015            CONCEPT DE TYPE :
    #  -------------                  -----------------
   IMPR_RESU(INFO=1,
            RESU=_F(INFO_MAILLAGE='NON',
                    RESULTAT=Res,
                    SOUS_TITRE='Depl and bottom-stress',
                    NOM_CHAM=('DEPL','SIGM_NOEU_DEPL','EQUI_NOEU_SIGM',)),
            UNITE=80,
            MODELE=FEM,
            FORMAT='MED',
            );

RESULTAT           : Res     
CHAMP              : DEPL           
   ==> NOM MED DU CHAMP : Res_____DEPL____________________

  . ECRITURE DES  6 *        4 VALEURS POUR LE NUMERO D'ORDRE :      1
     ( PAS DE TEMPS NUMERO :    1, T =    0.0000     )
 
---------------------------------------------------------------------------------

RESULTAT           : Res     
CHAMP              : SIGM_NOEU_DEPL 
   ==> NOM MED DU CHAMP : Res_____SIGM_NOEU_DEPL__________

  . ECRITURE DES  6 *        4 VALEURS POUR LE NUMERO D'ORDRE :      1
     ( PAS DE TEMPS NUMERO :    1, T =    0.0000     )
 
---------------------------------------------------------------------------------

RESULTAT           : Res     
CHAMP              : EQUI_NOEU_SIGM 
   ==> NOM MED DU CHAMP : Res_____EQUI_NOEU_SIGM__________

  . ECRITURE DES  6 *        4 VALEURS POUR LE NUMERO D'ORDRE :      1
     ( PAS DE TEMPS NUMERO :    1, T =    0.0000     )
 
---------------------------------------------------------------------------------

 
    #  FIN COMMANDE NO : 0015   DUREE TOTALE:        0.05s (SYST:        0.03s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0016            CONCEPT DE TYPE : evol_elas
    #  -------------                  -----------------
   Res=CALC_ELEM(reuse = Res,
                REPE_COQUE=_F(NUME_COUCHE=1,
                              NIVE_COUCHE='SUP',
                              ANGLE=0,
                              PLAN='MAIL'),
                INFO=1,
                OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
                TYPE_OPTION='TOUTES',
                PRECISION=1.E-3,
                RESULTAT=Res,
                CRITERE='RELATIF',
                CARA_ELEM=Cara,
                );

   
   !----------------------------------------------------------------!
   ! <A> <UTILITAI8_31>                                             !
   !                                                                !
   !                                                                !
   !   Option deja calcule:  option  SIGM_ELNO_DEPL  NUME_ORDRE  1 !
   !   On la recalcule car les donnes peuvent etre diffrentes     !
   !                                                                !
   !                                                                !
   !                                                                !
   !                                                                !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette !
   ! alarme, vous pouvez obtenir des rsultats inattendus !         !
   !----------------------------------------------------------------!
   
   
   !----------------------------------------------------------------!
   ! <A> <UTILITAI8_31>                                             !
   !                                                                !
   !                                                                !
   !   Option deja calcule:  option  EQUI_ELNO_SIGM  NUME_ORDRE  1 !
   !   On la recalcule car les donnes peuvent etre diffrentes     !
   !                                                                !
   !                                                                !
   !                                                                !
   !                                                                !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette !
   ! alarme, vous pouvez obtenir des rsultats inattendus !         !
   !----------------------------------------------------------------!
   
 
    #  FIN COMMANDE NO : 0016   DUREE TOTALE:        0.06s (SYST:        0.02s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0017            CONCEPT DE TYPE : evol_elas
    #  -------------                  -----------------
   Res=CALC_NO(reuse = Res,
              CRITERE='RELATIF',
              RESULTAT=Res,
              OPTION=('EQUI_NOEU_SIGM','SIGM_NOEU_DEPL',),
              PRECISION=1.E-3,
              );

   
   !----------------------------------------------------------------------!
   ! <A> <PREPOST5_1>                                                     !
   !                                                                      !
   !                                                                      !
   !  L'option EQUI_NOEU_SIGM est deja calculee pour le numero d'ordre 1. !
   !  On la recalcule car les donnees peuvent etre differentes.           !
   !                                                                      !
   !                                                                      !
   !                                                                      !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette       !
   ! alarme, vous pouvez obtenir des rsultats inattendus !               !
   !----------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------!
   ! <A> <PREPOST5_1>                                                     !
   !                                                                      !
   !                                                                      !
   !  L'option SIGM_NOEU_DEPL est deja calculee pour le numero d'ordre 1. !
   !  On la recalcule car les donnees peuvent etre differentes.           !
   !                                                                      !
   !                                                                      !
   !                                                                      !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette       !
   ! alarme, vous pouvez obtenir des rsultats inattendus !               !
   !----------------------------------------------------------------------!
   
 
    #  FIN COMMANDE NO : 0017   DUREE TOTALE:        0.01s (SYST:        0.00s)
    #  ---------------------------------------------------------------------------

    #  ---------------------------------------------------------------------------
    #  COMMANDE NO :  0018            CONCEPT DE TYPE :
    #  -------------                  -----------------
   IMPR_RESU(INFO=1,
            RESU=_F(INFO_MAILLAGE='NON',
                    RESULTAT=Res,
                    SOUS_TITRE='2D-top-stress',
                    NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM',)),
            UNITE=80,
            MODELE=FEM,
            FORMAT='MED',
            );

RESULTAT           : Res     
CHAMP              : SIGM_NOEU_DEPL 
   ==> NOM MED DU CHAMP : Res_____SIGM_NOEU_DEPL__________

   
   !-------------------------------------------------------------------------------------------!
   ! <A> <PREPOST5_24>                                                                         !
   !                                                                                           !
   !                                                                                           !
   !  fichier med :  fort.80 champ :  Res_____SIGM_NOEU_DEPL__ retour de mdexch : existc =  3  !
   !                                                                                           !
   !                                                                                           !
   !                                                                                           !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                            !
   ! alarme, vous pouvez obtenir des rsultats inattendus !                                    !
   !-------------------------------------------------------------------------------------------!
   
   
   !--------------------------------------------!
   ! <F> <PREPOST_73>                           !
   !                                            !
   !                                            !
   !  le champ est crit dans le fichier        !
   !                                            !
   !                                            !
   !                                            !
   ! Cette erreur est fatale. Le code s'arrete. !
   !--------------------------------------------!
   



DEPL and

Last edited by jo (2008-09-09 13:22:27)


[url=http://www.code-aster.de]www.code-aster.de[/url]                       Training & Support NASTRAN & CODE_ASTER

Offline

#2 2008-08-19 17:53:30

jo
Member
From: Ingenieurbüro für Mechanik, DE
Registered: 2008-01-29
Posts: 146
Website

Re: [RESOLU]How can I display top-and bot-stress of 2D-elements in Salome?

Bonjour,

I can't believe that I am the only person interested in postprocessing stresses of 2D-elements in Salome. So, could someone competent please respond to my questions 1), 2) and 3)?

I tryed an answer to my question 3), but I am still not shure if this is the easiest to handle way to display TOP and BOTTOM stress in Salome.

I create now an additional resultat using CREA_CHAMP and CREA_RESU. The new resultat contains the 2D-Top-stresses. But It seem that the new resultat must contain not only the nodal stresses but must beginn containing displacements, then Gauss-stresses, element-nodal stresses, then (at the happy end) the nodal stresses.

Question 4)
Is there an easier way to put only the nodal stresses into the new result?


I attach all files and hope to get into discussion.
Regards, Johannes


Attachments:
2D_stress_Top-and-Bottom.zip, Size: 28.95 KiB, Downloads: 430

[url=http://www.code-aster.de]www.code-aster.de[/url]                       Training & Support NASTRAN & CODE_ASTER

Offline

#3 2008-08-20 01:18:23

todd_alan_martin
Member
Registered: 2008-03-06
Posts: 131

Re: [RESOLU]How can I display top-and bot-stress of 2D-elements in Salome?

Hi

Perhaps you need to ask questions about using Salome in the Salome forum.

Regarding printing top and bottom stresses in a plate, for some reason code-aster doesn't like you re-using the same result object in CALC_ELEM. Instead create two new results and output them.

  Bottom=CALC_ELEM(
                REPE_COQUE=_F(NUME_COUCHE=1,
                              NIVE_COUCHE='INF',
                              PLAN='MAIL'),
                OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
                TYPE_OPTION='TOUTES',
                RESULTAT=Res,
                );

  Top=CALC_ELEM(
                REPE_COQUE=_F(NUME_COUCHE=1,
                              NIVE_COUCHE='SUP',
                              PLAN='MAIL'),
                OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
                TYPE_OPTION='TOUTES',
                RESULTAT=Res,
                );


IMPR_RESU(RESU=(_F(RESULTAT=Top,),
                _F(RESULTAT=Bottom,),
                _F(RESULTAT=Res,),),);

Offline

#4 2008-08-20 10:14:52

Thomas DE SOZA
Guru
From: EDF
Registered: 2007-11-23
Posts: 3,066

Re: [RESOLU]How can I display top-and bot-stress of 2D-elements in Salome?

Hi,

For question 2), fields of type ELNO (in Code_Aster terminology) that is fields that aren't continuous can only be displayed in Salomé 4.1.3 and above.

If you have an older version, you can't still view such fields in Gmsh for example (which as of version 2.2 can read MED format).

TdS

Offline

#5 2008-08-23 11:57:21

Fil
Member
From: Bayonne (64)
Registered: 2008-06-02
Posts: 5

Re: [RESOLU]How can I display top-and bot-stress of 2D-elements in Salome?

Hello,

Is there any way to obtain this latest version of Salome, or at least, what kind of improvment has been added compared to the 3.2.9 version.

Thanks,

Offline