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

You are not logged in.

#1 2019-10-16 02:49:24

nirmaljoshi
Member
From: Japan
Registered: 2018-10-12
Posts: 161

[solved] reading rmed file for post processing

Hi,
Which command should i use to read result.rmed file (created after very long processing time) so that i can post process it (CALC_CHAMP) without re-running the analysis. In the previous analysis i forgot to calculate SIGM_ELNO and the hdf file has broken. I only have the mesh file and "rmed", the result file.

The CREA_RESU has no option to read external files.

Please help.

Last edited by nirmaljoshi (2019-10-17 04:24:50)

Offline

#2 2019-10-16 19:05:41

sb1966
Member
Registered: 2015-03-16
Posts: 128

Re: [solved] reading rmed file for post processing

Hi,
I suppose you may add another stage (in AsterStudy) and use POURSUITE instead of Debut. There are examples discussed in the 2nd edition of Aubry's book, with reference to astk.

Hope this helps.

Regards.

Offline

#3 2019-10-17 04:24:30

nirmaljoshi
Member
From: Japan
Registered: 2018-10-12
Posts: 161

Re: [solved] reading rmed file for post processing

Thanks. I could do the calculation.
For completeness, i post here the steps.
#1. Read the old result.rmed file by LIRE_MAILLAGE
#2. Use LIRE_RESU to result result variable (you should know the result variable)
#3. Define same material properties and assign to model and also while using CALC_CHAMP
------

mesh = LIRE_MAILLAGE(FORMAT='MED',
                     UNITE=20)

model = AFFE_MODELE(AFFE=_F(MODELISATION=('3D', ),
                            PHENOMENE='MECANIQUE',
                            TOUT='OUI'),
                    MAILLAGE=mesh)

steel = DEFI_MATER_GC(ACIER=_F(D_SIGM_EPSI=10000000000.0,
                               E=2e+11,
                               NU=0.3,
                               SY=400000000.0),
                      RHO=7850.0)

conMazMa = DEFI_MATERIAU(ELAS=_F(E=37303280300.0,
                                 NU=0.2,
                                 RHO=2400.0),
                         MAZARS=_F(AC=1.26865552,
                                   AT=0.9,
                                   BC=1681.54688,
                                   BT=12688.1906,
                                   EPSD0=7.88134441e-05,
                                   K=0.7,
                                   SIGM_LIM=39000000.0))

fieldmat = AFFE_MATERIAU(AFFE=(_F(GROUP_MA=('group_concrete', ),
                                  MATER=(conMazMa, )),
                               _F(GROUP_MA=('group_mainrebar', 'group_shearrebars'),
                                  MATER=(steel, ))),
                         MAILLAGE=mesh,
                         MODELE=model)

oldres = LIRE_RESU(FORMAT='MED',
                   FORMAT_MED=_F(NOM_CHAM='DEPL',
                                 NOM_RESU='resnonl0'),
                   MODELE=model,
                   TOUT_ORDRE='OUI',
                   TYPE_RESU='MODE_MECA',
                   UNITE=20)

stress = CALC_CHAMP(CHAM_MATER=fieldmat,
                    CONTRAINTE=('SIGM_ELNO', ),
                    DEFORMATION=('EPSI_ELNO', ),
                    MODELE=model,
                    RESULTAT=oldres)

IMPR_RESU(FORMAT='MED',
          RESU=_F(RESULTAT=stress),
          UNITE=80)

Offline