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

You are not logged in.

#1 2019-10-06 13:44:51

marcelo
Member
Registered: 2017-06-20
Posts: 44

REAC_NODA error?!

Hi,

I have performed a three element study (solid, plate and beam) and I am not confident of the results presented by the REAC_NODA analysis (link: encurtador.com.br/tGTVX). I don't know if the iteration I made between connecting the elements was correct. Below is the .comm file:

DEBUT(LANG='EN')

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

mesh = DEFI_GROUP(reuse=mesh,
                  CREA_GROUP_MA=_F(NOM='TOUT',
                                   TOUT='OUI'),
                  CREA_GROUP_NO=_F(TOUT_GROUP_MA='OUI'),
                  MAILLAGE=mesh)

model = AFFE_MODELE(AFFE=(_F(MODELISATION=('3D', ),
                             PHENOMENE='MECANIQUE',
                             TOUT='OUI'),
                          _F(GROUP_MA=('Plate', ),
                             MODELISATION=('DKT', ),
                             PHENOMENE='MECANIQUE'),
                          _F(GROUP_MA=('Beam', ),
                             MODELISATION=('POU_D_T', ),
                             PHENOMENE='MECANIQUE')),
                    MAILLAGE=mesh)

elempro = AFFE_CARA_ELEM(COQUE=_F(EPAIS=10.0,
                                  GROUP_MA=('Plate', )),
                         MODELE=model,
                         POUTRE=_F(CARA=('HY', 'HZ'),
                                   GROUP_MA=('Beam', ),
                                   SECTION='RECTANGLE',
                                   VALE=(100.0, 10.0)))

mater = DEFI_MATERIAU(ELAS=_F(E=210000.0,
                              NU=0.3,
                              RHO=7.8e-09))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh)

Fix = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0,
                                 DY=0.0,
                                 DZ=0.0,
                                 GROUP_MA=('Fix', )),
                     MODELE=model)

ForceZ = AFFE_CHAR_MECA(FORCE_NODALE=_F(FZ=2000.0,
                                        GROUP_NO=('Load', )),
                        MODELE=model)

Solid1 = AFFE_CHAR_MECA(LIAISON_SOLIDE=_F(GROUP_NO=('Solid1', )),
                        MODELE=model)

Solid2 = AFFE_CHAR_MECA(LIAISON_SOLIDE=_F(GROUP_NO=('Solid2', )),
                        MODELE=model)

reslin = MECA_STATIQUE(CARA_ELEM=elempro,
                       CHAM_MATER=fieldmat,
                       EXCIT=(_F(CHARGE=Fix),
                              _F(CHARGE=ForceZ),
                              _F(CHARGE=Solid1),
                              _F(CHARGE=Solid2)),
                       MODELE=model)

reslin = CALC_CHAMP(reuse=reslin,
                    CONTRAINTE=('SIEF_ELNO', 'SIEF_ELGA', 'SIPO_ELNO', 'SIPO_NOEU', 'SIGM_ELNO', 'SIGM_NOEU', 'EFGE_NOEU'),
                    CRITERES=('SIEQ_ELNO', 'SIEQ_NOEU'),
                    FORCE=('REAC_NODA', ),
                    RESULTAT=reslin)

IMPR_RESU(FORMAT='MED',
          RESU=_F(INFO_MAILLAGE='OUI',
                  RESULTAT=reslin,
                  TOUT_CMP='OUI'),
          UNITE=80)

FIN()

I would like to know if I made the right choices and if there is any way to improve.


Attachments:
message, Size: 44.32 KiB, Downloads: 15

Offline

#2 2019-10-09 13:53:30

marcelo
Member
Registered: 2017-06-20
Posts: 44

Re: REAC_NODA error?!

Can anybody help me?

Offline