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

You are not logged in.

#1 2020-06-29 06:13:27

hkboondoggle
Member
Registered: 2020-05-29
Posts: 13

Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE)

Hi all,

I am trying to simulate an impact scenario with a ball travelling on negative z-axis direction and being about to collide on a static plate.

With the aid of konyaro, Johannes, mecour, and others (I really appreciate it guys!), I have set up the case which details can be found in the attached zip file.
It includes the .py script for geometry & mesh generation (.med too large), the .comm file and .log of the problematic Run case.

The following error happened as the Run Case was being launched:

.. __stg1_cmd9:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0012 Concept of the type:  evol_noli
  # ------------------------------------------------------------------------------------------
  solver = DYNA_NON_LINE(MODELE=model1,
                         CHAM_MATER=fieldma0,
                         CONTACT=contact,
                         ETAT_INIT=_F(VITE=init_vel,
                                      PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INCREMENT=_F(LIST_INST=timeline,
                                      PRECISION=1.E-06,),
                         SCHEMA_TEMPS=_F(SCHEMA='HHT',
                                         ALPHA=-0.3,
                                         MODI_EQUI='NON',
                                         FORMULATION='DEPLACEMENT',
                                         COEF_MASS_SHIFT=0.0,),
                         AMOR_RAYL_RIGI='TANGENTE',
                         METHODE='NEWTON',
                         NEWTON=_F(REAC_INCR=1,
                                   MATRICE='TANGENTE',
                                   REAC_ITER=1,
                                   REAC_ITER_ELAS=0,
                                   MATR_RIGI_SYME='NON',),
                         CONVERGENCE=_F(ITER_GLOB_MAXI=10,
                                        ITER_GLOB_ELAS=25,
                                        ARRET='OUI',),
                         SOLVEUR=_F(RENUM='AUTO',
                                    NPREC=8,
                                    ELIM_LAGR='LAGR2',
                                    STOP_SINGULIER='OUI',
                                    TYPE_RESOL='AUTO',
                                    ACCELERATION='AUTO',
                                    LOW_RANK_SEUIL=0.0,
                                    PRETRAITEMENTS='AUTO',
                                    POSTTRAITEMENTS='AUTO',
                                    PCENT_PIVOT=20,
                                    RESI_RELA=-1.0,
                                    GESTION_MEMOIRE='AUTO',
                                    FILTRAGE_MATRICE=-1.0,
                                    MIXER_PRECISION='NON',
                                    MATR_DISTRIBUEE='NON',
                                    METHODE='MUMPS',),
                         MESURE=_F(TABLE='NON',),
                         ARCHIVAGE=_F(PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INFO=1,)


Comme vous n'avez pas défini explicitement le comportement, tout le modèle est supposé élastique en petites perturbations.
   Liste des comportements
   Affecté sur 186954 éléments
     Relation                             : ELAS
     Déformation                          : PETIT
     Nombre total de variables internes   : 1
            V1 : VIDE
  Le système linéaire à résoudre a 734008 degrés de liberté:
   - 734008 sont des degrés de liberté physiques
     (ils sont portés
par 244274 noeuds du maillage)
   - 0 sont les couples de paramètres de Lagrange associés
     aux 0 relations linéaires dualisées.
  La matrice est de taille 734008 équations.
  Elle contient 31072249 termes non nuls si elle est symétrique et 61410490 termes non
nuls si elle n'est pas symétrique.
  Soit un taux de remplissage de   0.011 %.
    Lecture de l'état initial
      Il n'y a pas d'état initial défini. On prend un état initial nul.
      Le champ <DEPL> est initialisé a zéro
      Le champ <SIEF_ELGA> est initialisé a zéro
      Le champ <VARI_ELGA> est initialisé a zéro
   
   !-----------------------------------------------------------------!
   ! <A> <ALGELINE7_20>                                              !
   !                                                                 !
   !  On ne peut pas remplir la composante DX du noeud numéro N1     !
   !  pour le champ "&&NMETL2.CHAMP.CONV"                            !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------------!
   ! <A> <MECANONLINE_2>                                                                    !
   !                                                                                        !
   !  Lors de la recopie du champ &&NMETL2.CHAMP.CONVER donné dans ETAT_INIT de la commande !
   !  STAT_NON_LINE vers le champ                                                           !
   ! &&NMCH1P.VITMOI, certaines composantes de &&NMCH1P.VITMOI                              !
   !  ont du être mises à zéro.                                                             !
   !                                                                                        !
   !  Ce problème survient lorsque le champ donné                                           !
   ! dans ETAT_INIT ne comporte                                                             !
   !  pas assez de composantes, on complète donc par des zéros.                             !
   !                                                                                        !
   !                                                                                        !
   ! This is a warning. If you do not understand the meaning of this                        !
   !  warning, you can obtain unexpected results!                                           !
   !----------------------------------------------------------------------------------------!
   
      Le champ <VITE> est lu dans ETAT_INIT, par un champ donné explicitement
      Le champ <ACCE> est initialisé a zéro

One considers an acceleration initial.

The initial state does not have acceleration given.
It is computed.


Filing of the initial state

Filing of the fields
    Field stored  DEPL at time  0.000000000000e+00 for the sequence number  0
    Field stored  SIEF_ELGA at time  0.000000000000e+00 for the sequence number  0
    Field stored  VARI_ELGA at time  0.000000000000e+00 for the sequence number  0
    Field stored  COMPORTEMENT at time  0.000000000000e+00 for the sequence number  0
    Field stored  VITE at time  0.000000000000e+00 for the sequence number  0
    Field stored  ACCE at time  0.000000000000e+00 for the sequence number  0
    Field stored  CONT_NOEU at time  0.000000000000e+00 for the sequence number  0
   Il y a 0 points initialement en contact et 0 points exclus.
-----------------------------------------------------------------------------------------------------------------------------------------

Time of computation:   4.000000000000e-03
-----------------------------------------------------------------------------------------------------------------------------------------
|     CONTACT    |     NEWTON     |     RESIDU     |     RESIDU     |     OPTION     |     CONTACT    |     CONTACT    |     CONTACT    |
|    BCL. GEOM.  |    ITERATION   |     RELATIF    |     ABSOLU     |   ASSEMBLAGE   |   NEWTON GENE  |    PRESSURE    |     CRITERE    |
|    ITERATION   |                | RESI_GLOB_RELA | RESI_GLOB_MAXI |                |   VARI. CONT.  |    ERROR       |    VALEUR      |
-----------------------------------------------------------------------------------------------------------------------------------------
Traceback returned by pdbhelp:
[0] print_trace_
[1] utmess_core_
[2] utmess_
[3] assert_
[4] mmcaln_
[5] mmvppe_
[6] te0365_
[7] te0000_
[8] calcul_
[9] nmelcv_
[10] nmfocc_
[11] nmprta_
[12] nmpred_
[13] nmnewt_
[14] op0070_
[15] execop_
[16] expass_
[17] initaster_fonctions
[18] PyMethodDef_RawFastCallKeywords
[19] PyMethodDef_RawFastCallKeywords
[20] PyEval_EvalFrameDefault
[21] PyMethodDef_RawFastCallKeywords
[22] PyEval_EvalFrameDefault
[23] PyMethodDef_RawFastCallKeywords
[24] PyEval_EvalFrameDefault
   
   !-----------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <DVP_1>                                                                           !
   !                                                                                               !
   ! Program error.                                                                                !
   !                                                                                               !
   ! Condition not met:                                                                            !
   !     ASTER_FALSE                                                                               !
   ! File                                                                                          !
   ! /home/siavelis/PROJETS/aster/windows/build/codeaster-src/bibfor/cont_elem/mmcaln.F90, line 63 !
   !                                                                                               !
   ! --------------------------------------------                                                  !
   ! Contexte du message :                                                                         !
   !    Option         : CHAR_MECA_CONT                                                            !
   !    Type d'élément : COT6T6                                                                    !
   !                                                                                               !
   ! Maillage       : mesh                                                                         !
   !    Maille         : XXX                                                                       !
   !    Type de maille : TRIA66                                                                    !
   !    Cette maille appartient aux groupes de mailles suivants                                    !
   ! :                                                                                             !
   !                                                                                               !
   !    Position du centre de gravité de la maille :                                               !
   !       x=nan y=0.000000 z=0.000000                                                             !
   !                                                                                               !
   !                                                                                               !
   !                                                                                               !
   ! Il y a probablement une erreur dans la programmation.                                         !
   ! Veuillez contacter votre assistance technique.                                                !
   !-----------------------------------------------------------------------------------------------!
   
   !-----------------------------------------------------------------!
   !                                                                 !
   !    FATAL ERROR detected in Code_Aster for Windows version       !
   !                                                                 !
   !    If an ABNORMAL_ABORT occured, that is not reproductible      !
   !    on official Linux version, or if you need some help to       !
   !    understand the message above, please report issues at:       !
   !                                                                 !
   !                      support@simulease.com                      !
   !                                                                 !
   !-----------------------------------------------------------------!

Could you please advice in dealing this problem?
Thanks!

I have attached the mesh info below FYI:

50056298003_4ee4aa89f8_o.jpg

50057117222_c947c2115b_o.jpg


Attachments:
Case_Files.zip, Size: 13.4 KiB, Downloads: 6

Offline

#2 2020-07-03 09:40:59

hkboondoggle
Member
Registered: 2020-05-29
Posts: 13

Re: Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE)

Thinking that the problem might be caused by the stretching of the mesh, I have replace the whole mesh set up with an uniform setting shown below:

50071164116_1540215690_o.jpg

However, the same exception still occurs.
The comm, log & med files of this latest trial are attached FYI.

Last edited by hkboondoggle (2020-07-03 09:41:48)


Attachments:
Trial_2_Study_Files.zip, Size: 1.27 MiB, Downloads: 4

Offline

#3 Yesterday 08:49:38

konyaro
Member
From: Genève
Registered: 2016-04-02
Posts: 246

Re: Advice needed on <EXCEPTION> <DVP_1> (DYNA_NON_LINE)

Hello,
your case runs fine. It requires a lot of RAM as there are many nodes and you're using a direct solver. I guess that's the source of your problem. The Windows version doesn't depict the memory error messages. You should try to increase the amount of memory or reduce the size of the mesh.

In the attached files I switched to larger linear elements and changed a few things in your comm file: GROT_GDEP, contact, reuse for CALC_CHAMP, less steps as I don't want to wait for hours. It does converge well.

Konyaro


Attachments:
ball2.zip, Size: 427.24 KiB, Downloads: 1

失敗は成功のもと (L'échec est la base de la réussite)

Offline