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

You are not logged in.

#1 2022-01-13 14:34:36

Teddy
Member
From: Deutschland
Registered: 2022-01-13
Posts: 5

Error in STAT_NON_LINE

Hello dear community,

On the first day of my accession, I already have a question, I have a problem with a shell element that I have previously calculated statically.
If I switch to STAT_NON_LINE, then I get an error that I can't quite follow.

Maybe someone can help me with this?

P.S. I'm a newbie, for tips on literature etc. please contact me smile



Mesh3 = LIRE_MAILLAGE(identifier='0:1',

                      FORMAT='MED',

                      UNITE=2)



model = AFFE_MODELE(identifier='1:1',

                    AFFE=_F(GROUP_MA=('Gesamt', ),

                            MODELISATION=('DST', ),

                            PHENOMENE='MECANIQUE'),

                    MAILLAGE=Mesh3)



#comment: elemprop = AFFE_CARA_ELEM(identifier='2:1',

#comment:                           COQUE=_F(EPAIS=2.0,

#comment:                                    GROUP_MA=('Partition_1', ),

#comment:                                    VECTEUR=(1.0, 1.0, 1.0)),

#comment:                           MODELE=model)



mater = DEFI_MATERIAU(identifier='3:1',

                      ELAS=_F(E=210000.0,

                              NU=0.3,

                              RHO=7.85e-09))



fieldmat = AFFE_MATERIAU(identifier='4:1',

                         AFFE=_F(MATER=(mater, ),

                                 TOUT='OUI'),

                         MODELE=model)



func = DEFI_FONCTION(identifier='5:1',

                     NOM_PARA='INST',

                     PROL_DROITE='CONSTANT',

                     PROL_GAUCHE='CONSTANT',

                     VALE=(0.0, 0.0, 1.0, 1.0))



listr = DEFI_LIST_REEL(identifier='6:1',

                       DEBUT=0.0,

                       INTERVALLE=_F(JUSQU_A=1.0,

                                     NOMBRE=6))



load = AFFE_CHAR_MECA(identifier='7:1',

                      DDL_IMPO=_F(GROUP_MA=('Schraubenauflage', ),

                                  LIAISON='ENCASTRE'),

                      MODELE=model)



load0 = AFFE_CHAR_MECA(identifier='8:1',

                       MODELE=model,

                       PRES_REP=_F(GROUP_MA=('Druckbereich', ),

                                   PRES=-0.084))



#comment: load1 = AFFE_CHAR_MECA(identifier='9:1',

#comment:                        DDL_IMPO=_F(DY=0.0,

#comment:                                    GROUP_MA=('Fix', )),

#comment:                        MODELE=model)



#comment: reslin = MECA_STATIQUE(identifier='10:1',

#comment:                        CARA_ELEM=elemprop,

#comment:                        CHAM_MATER=fieldmat,

#comment:                        EXCIT=(_F(CHARGE=load),

#comment:                               _F(CHARGE=load0)),

#comment:                        MODELE=model,

#comment:                        SOLVEUR=_F(METHODE='MUMPS',

#comment:                                   NPREC=8,

#comment:                                   RESI_RELA=0.001))



resnonl = STAT_NON_LINE(identifier='11:1',

                        CHAM_MATER=fieldmat,

                        COMPORTEMENT=_F(),

                        EXCIT=(_F(CHARGE=load),

                               _F(CHARGE=load0)),

                        INCREMENT=_F(LIST_INST=listr),

                        MODELE=model)



resnonl = CALC_CHAMP(identifier='12:1',

                     reuse=resnonl,

                     CRITERES=('SIEQ_ELNO', ),

                     DEFORMATION=('EPSI_ELNO', ),

                     RESULTAT=resnonl)



#comment: unnamed0 = POST_CHAMP(identifier='13:1',

#comment:                       EXTR_COQUE=_F(NIVE_COUCHE='INF',

#comment:                                     NOM_CHAM=('SIGM_ELNO', 'SIGM_ELGA'),

#comment:                                     NUME_COUCHE=1),

#comment:                       RESULTAT=reslin)



#comment: unnamed1 = CALC_CHAMP(identifier='14:1',

#comment:                       CRITERES=('SIEQ_NOEU', 'SIEQ_ELNO'),

#comment:                       RESULTAT=unnamed0)



IMPR_RESU(identifier='15:1',

          FORMAT='MED',

          RESU=_F(MAILLAGE=Mesh3,

                  NOM_CHAM=('DEPL', 'SIEQ_NOEU'),

                  RESULTAT=resnonl),

          UNITE=80)



Erreur utilisateur :                                                           !
   !   Vous essayez de faire un calcul non-linaire mcanique ou un post-traitement !
   ! sur un modle dont les                                                         !
   ! lments                                                                       !
   !   ne sont pas programms pour cela.                                            !
   !   On arrte le calcul.                                                         !
   ! Risques & conseils :                                                           !
   !   Vous devriez changer de modlisation.


Be your own Hero

Offline

#2 2022-01-13 16:35:10

mihe
Member
From: ENGCALengineeringcalculations
Registered: 2008-04-08
Posts: 259

Re: Error in STAT_NON_LINE

Hi,
CARA_ELEM=elemprop, is missing for your STAT_NON_LINE case.
BR/Micke


Salome Meca 2019 on Ubuntu 18

Offline

#3 2022-01-13 20:12:05

Teddy
Member
From: Deutschland
Registered: 2022-01-13
Posts: 5

Re: Error in STAT_NON_LINE

Hi Mike,

Unfortunately that wasn't the problem, I've changed the analysis so many times now that I forgot the elemprops in the last analysis.

Do you perhaps know what the error message could mean?

Regards, Theo


Mesh3 = LIRE_MAILLAGE(identifier='0:1',

                      FORMAT='MED',

                      UNITE=2)



model = AFFE_MODELE(identifier='1:1',

                    AFFE=_F(GROUP_MA=('Gesamt', ),

                            MODELISATION=('DST', ),

                            PHENOMENE='MECANIQUE'),

                    MAILLAGE=Mesh3)



elemprop = AFFE_CARA_ELEM(identifier='2:1',

                          COQUE=_F(EPAIS=2.0,

                                   GROUP_MA=('Partition_1', ),

                                   VECTEUR=(1.0, 1.0, 1.0)),

                          MODELE=model)



func = DEFI_FONCTION(identifier='3:1',

                     NOM_PARA='INST',

                     PROL_DROITE='CONSTANT',

                     PROL_GAUCHE='CONSTANT',

                     VALE=(0.0, 0.0, 1.0, 1.0))



listr = DEFI_LIST_REEL(identifier='4:1',

                       DEBUT=0.0,

                       INTERVALLE=_F(JUSQU_A=1.0,

                                     NOMBRE=6))



FSI = DEFI_FONCTION(identifier='5:1',

                    NOM_PARA='EPSI',

                    PROL_DROITE='LINEAIRE',

                    PROL_GAUCHE='LINEAIRE',

                    VALE=(0.002, 420.0, 0.004, 500.0))



mater = DEFI_MATERIAU(identifier='6:1',

                      ELAS=_F(E=210000.0,

                              NU=0.3,

                              RHO=7.85e-09),

                      TRACTION=_F(SIGM=FSI))



fieldmat = AFFE_MATERIAU(identifier='7:1',

                         AFFE=_F(MATER=(mater, ),

                                 TOUT='OUI'),

                         MAILLAGE=Mesh3,

                         MODELE=model)



load = AFFE_CHAR_MECA(identifier='8:1',

                      DDL_IMPO=_F(GROUP_MA=('Schraubenauflage', ),

                                  LIAISON='ENCASTRE'),

                      MODELE=model)



load0 = AFFE_CHAR_MECA(identifier='9:1',

                       MODELE=model,

                       PRES_REP=_F(GROUP_MA=('Druckbereich', ),

                                   PRES=-0.084))



#comment: load1 = AFFE_CHAR_MECA(identifier='10:1',

#comment:                        DDL_IMPO=_F(DY=0.0,

#comment:                                    GROUP_MA=('Fix', )),

#comment:                        MODELE=model)



#comment: reslin = MECA_STATIQUE(identifier='11:1',

#comment:                        CARA_ELEM=elemprop,

#comment:                        CHAM_MATER=fieldmat,

#comment:                        EXCIT=(_F(CHARGE=load),

#comment:                               _F(CHARGE=load0)),

#comment:                        MODELE=model,

#comment:                        SOLVEUR=_F(METHODE='MUMPS',

#comment:                                   NPREC=8,

#comment:                                   RESI_RELA=0.001))



resnonl = STAT_NON_LINE(identifier='12:1',

                        CARA_ELEM=elemprop,

                        CHAM_MATER=fieldmat,

                        COMPORTEMENT=_F(DEFORMATION='PETIT',

                                        RELATION='VMIS_ISOT_TRAC',

                                        TOUT='OUI'),

                        EXCIT=(_F(CHARGE=load),

                               _F(CHARGE=load0)),

                        INCREMENT=_F(LIST_INST=listr),

                        METHODE='NEWTON',

                        MODELE=model)



resnonl = CALC_CHAMP(identifier='13:1',

                     reuse=resnonl,

                     CRITERES=('SIEQ_ELNO', ),

                     DEFORMATION=('EPSI_ELNO', ),

                     RESULTAT=resnonl)



#comment: unnamed0 = POST_CHAMP(identifier='14:1',

#comment:                       EXTR_COQUE=_F(NIVE_COUCHE='INF',

#comment:                                     NOM_CHAM=('SIGM_ELNO', 'SIGM_ELGA'),

#comment:                                     NUME_COUCHE=1),

#comment:                       RESULTAT=reslin)



#comment: unnamed1 = CALC_CHAMP(identifier='15:1',

#comment:                       CRITERES=('SIEQ_NOEU', 'SIEQ_ELNO'),

#comment:                       RESULTAT=unnamed0)



IMPR_RESU(identifier='16:1',

          FORMAT='MED',

          RESU=_F(MAILLAGE=Mesh3,

                  NOM_CHAM=('DEPL', 'SIEQ_NOEU'),

                  RESULTAT=resnonl),

          UNITE=80)


<S> Exception user raised but not interceptee.                                 !
   ! The bases are fermees.                                                         !
   ! Type of the exception: error                                                   !
   !                                                                                !
   ! Erreur utilisateur :                                                           !
   !   Vous essayez de faire un calcul non-linaire mcanique ou un post-traitement !
   ! sur un modle dont les                                                         !
   ! lments                                                                       !
   !   ne sont pas programms pour cela.                                            !
   !   On arrte le calcul.                                                         !
   ! Risques & conseils :                                                           !
   !   Vous devriez changer de modlisation.                                        !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fine CR of execution of JDC in MIXTE


Be your own Hero

Offline

#4 2022-01-13 21:43:48

mihe
Member
From: ENGCALengineeringcalculations
Registered: 2008-04-08
Posts: 259

Re: Error in STAT_NON_LINE

No i do not. Maybe change to linear elastic material first and then troubleshoot from there.
BR/Micke


Salome Meca 2019 on Ubuntu 18

Offline

#5 2022-01-13 22:26:09

AsterO'dactyle
Administrator
Registered: 2007-11-29
Posts: 413

Re: Error in STAT_NON_LINE

Hello,

DST elements are not allowed in non-linear analysis.


Code_Asterの開発者

Offline