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

You are not logged in.

#1 2021-01-14 15:10:40

VonPire
Member
Registered: 2021-01-14
Posts: 3

Is it possible get SIEQ_ELNO and SIEG_ELGA from Coque_3D

Hello,

I have a problem with Coque_3d post processing with Paravis.

I've read that there is two possible ways to print stresses in Paravis from COQUE_3D model:

1. Duplicate Quadratic (QUAD8 or TRIA6 elements) mesh with CREA_MAILLE before changing the mesh to QUAD9 or TRIA7. Then at IMPR_RESU requesting results for the copied quadratic mesh GROUP_MA.

2. Project results from QUAD9/TRIA7 model to QUAD8/TRIA6 model.

I've managed to get DEPL, SIGM_NOEU and SIEQ_NOEU visible in Paravis, but not SIEQ_ELNO, SIEQ_ELGA or SIGM_ELNO, SIGM_ELGA with the procedure number 2 (projecting):
------------------------------------------------------------------------------------------------
resnonl = STAT_NON_LINE(identifier='14:1',
                        CARA_ELEM=elempro0,
                        CHAM_MATER=fieldmat,
                        COMPORTEMENT=_F(DEFORMATION='GROT_GDEP',
                                        RELATION='VMIS_ISOT_TRAC'),
                        CONVERGENCE=_F(ITER_GLOB_MAXI=100),
                        EXCIT=(_F(CHARGE=FIXL),
                               _F(CHARGE=load1,
                                  FONC_MULT=func)),
                        INCREMENT=_F(LIST_INST=times),
                        METHODE='NEWTON',
                        MODELE=model,
                        NEWTON=_F(MATRICE='TANGENTE',
                                  PREDICTION='TANGENTE',
                                  REAC_ITER=1))

proj0 = PROJ_CHAMP(identifier='15:1',
                   MODELE_1=model,
                   MODELE_2=model1,
                   RESULTAT=resnonl,
                   TOUT_CHAM='OUI',
                   TOUT_ORDRE='OUI')

resnonl = CALC_CHAMP(identifier='16:1',
                     reuse=resnonl,
                     CONTRAINTE=('SIGM_ELNO', ),
                     RESULTAT=resnonl,
                     TOUT_ORDRE='OUI')

ResINF = POST_CHAMP(identifier='17:1',
                    EXTR_COQUE=_F(NIVE_COUCHE='INF',
                                  NOM_CHAM=('SIGM_ELNO', 'SIEF_ELGA', 'SIEQ_ELNO'),
                                  NUME_COUCHE=1),
                    RESULTAT=resnonl,
                    TOUT_ORDRE='OUI')

ResINF = CALC_CHAMP(identifier='18:1',
                    reuse=ResINF,
                    CONTRAINTE=('SIGM_NOEU', 'SIGM_ELNO'),
                    CRITERES=('SIEQ_ELNO', 'SIEQ_NOEU'),
                    RESULTAT=ResINF,
                    TOUT_ORDRE='OUI')

projINF = PROJ_CHAMP(identifier='19:1',
                     MODELE_1=model,
                     MODELE_2=model1,
                     RESULTAT=ResINF,
                     TOUT_CHAM='OUI',
                     TOUT_ORDRE='OUI')

ResSup = POST_CHAMP(identifier='20:1',
                    EXTR_COQUE=_F(NIVE_COUCHE='SUP',
                                  NOM_CHAM=('SIGM_ELNO', 'SIEF_ELGA', 'SIEQ_ELNO'),
                                  NUME_COUCHE=1),
                    RESULTAT=resnonl)

ResSup = CALC_CHAMP(identifier='21:1',
                    reuse=ResSup,
                    CONTRAINTE=('SIGM_NOEU', 'SIGM_ELNO'),
                    CRITERES=('SIEQ_ELNO', 'SIEQ_NOEU'),
                    RESULTAT=ResSup)

projSUP = PROJ_CHAMP(identifier='22:1',
                     MODELE_1=model,
                     MODELE_2=model1,
                     RESULTAT=ResSup,
                     TOUT_CHAM='OUI',
                     TOUT_ORDRE='OUI')

IMPR_RESU(identifier='23:1',
          FORMAT='MED',
          RESU=(_F(RESULTAT=projSUP,
                   TOUT_CHAM='OUI',
                   TOUT_ORDRE='OUI'),
                _F(RESULTAT=projINF,
                   TOUT_CHAM='OUI',
                   TOUT_ORDRE='OUI'),
                _F(RESULTAT=proj0,
                   TOUT_CHAM='OUI',
                   TOUT_ORDRE='OUI')),
          UNITE=80)

FIN()
------------------------------------------------------------------------------------------------------

Are there anyone who had the same problem and managed to get it work ?

It is also a little bit unclear what TOUT_ORDRE does. Is this only for modal analysis ?

I have a simple test model for this where 2mm sheet metal is fixed from the bottom edge and load is applied to the upper edge. The .comm file is as attachment.


Attachments:
RunCase_3_Stage_1.comm, Size: 6.41 KiB, Downloads: 3

Offline

#2 2021-01-14 17:03:16

VonPire
Member
Registered: 2021-01-14
Posts: 3

Re: Is it possible get SIEQ_ELNO and SIEG_ELGA from Coque_3D

Here is also the .med file as attachment


Attachments:
_ExportedFromSalomeObject_0_1_2_3.med, Size: 68.24 KiB, Downloads: 1

Offline