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

jo · 2008-08-07 17:51:00

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)

jo · 2008-08-19 17:53:30

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: 426

todd_alan_martin · 2008-08-20 01:18:23

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,),),);

Thomas DE SOZA · 2008-08-20 10:14:52

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

Fil · 2008-08-23 11:57:21

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,

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

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

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

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

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

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

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

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

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

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

Board footer