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: 17

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: 19

Offline

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

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

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: 16

Offline

#3 2020-07-08 08:49:38

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

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: 15

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

Offline

#4 2020-07-17 06:52:22

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

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

Thanks Konyaro for your inspection.

I have shorten the calculation time from 10 to 3.24 (<DVP_1> would happen at 3.26), and the calculation does converge well.
Although the result is physically incorrect again with the ball re-bounces backward after the collision, whereas the reference I'm following [1] has the ball undergoes plastic deformation and stays attached to the substrate.

The material properties and initial velocity of the ball in my code-aster simulation follow the presented data in [1] and are the normalized with respect to the density of copper, the speed of sound, the diameter of the ball, and the room temperature.

After inspection, I found that the RELATION of behavior of DYNA_NON_LINE by default is 'ELASTIC'.
Thinking that this might be the reason of the discrepancy of the calculated result, I changed it to 'VMIS_JOHN_COOK'.
Other changes included halving the time increment to 5e-3, and enabling Deformation in CLAC_CHAMP.

Now it appears to be having trouble in convergence, as time increment automatically down-sizing itself to the order of -16 once the ball impacts the substrate &  'VMIS_JOHN_COOK' kicks into action, which ultimately leads to <MECANONLINE9_3>.
The comm, log & med files of this latest trial are attached FYI.

I also tried switching DEFORMATION of behavior of DYNA_NON_LINE from 'GROT_GDEP' to 'SIMO_MIEHE', as high strain rate deformation is expected upon the collision.
The calculation refused to launch and instantly gave out <F>_ABNORMAL_ABORT with EXECUTION_CODE_ASTER_EXIT_8627=0.

Any suggestions / ideas in handling the latest issue?

[1] T. Schmidt, F. Gärtner, H. Assadi, et al., "Development of a generalized parameter window for cold spray deposition", Acta Materialia, 2006, 54, pp 729-742.

Last edited by hkboondoggle (2020-07-17 06:55:10)


Attachments:
Case_Files.zip, Size: 635.26 KiB, Downloads: 16

Offline

#5 2020-07-29 08:53:52

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

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

Hello, could some one please help me out?

Offline

#6 2020-07-29 09:20:39

sameer21101970
Member
Registered: 2019-09-06
Posts: 177

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

You may try below, hope it will work.

Apply Liasion Mail Command to ball & surface.
Apply Velocity or acceleration to ball.

This will ensure your ball remain stick to the surface and you can study deformation.

Offline

#7 2020-07-31 09:25:00

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

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

Hello sameer, thank you for your reply.
May I ask a few questions about the 'Liaison_Mail' approach:

1. I was bumped with the following problem when trying out 'Liaison_Mail'. Can you please explain what the matter is?

.. __stg1_cmd6:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0009 Concept of the type:  char_meca
  # ------------------------------------------------------------------------------------------
  BC_Fix = AFFE_CHAR_MECA(MODELE=model1,
                          DDL_IMPO=_F(GROUP_NO=('Substrate_Bottom', ),
                                      DX=0.0,
                                      DY=0.0,
                                      DZ=0.0,),
                          LIAISON_MAIL=_F(GROUP_MA_MAIT=('Substrate_Top', ),
                                          GROUP_MA_ESCL=('Particle_Face', ),
                                          TYPE_RACCORD='MASSIF',
                                          ELIM_MULT='NON',),
                          VERI_NORM='OUI',
                          DOUBLE_LAGRANGE='OUI',
                          INFO=1,)

   
   !-----------------------------------------------------------------------------!
   ! <EXCEPTION> <CALCULEL4_55>                                                  !
   !                                                                                                          !
   !  Pas de mailles à projeter ou en correspondance.                            !
  !                                                                                                          !
   !  Dans le cas de l'utilisation de AFFE_CHAR_MECA / LIAISON_MAIL, les mailles !
   ! maîtres                                                                                        !
   !  doivent avoir la même dimension que l'espace de modélisation :             !
   !  - mailles volumiques pour un modèle 3D                                     !
   !  - mailles                                                                                           !
   ! surfaciques pour un modèle 2D                                                     !
   !-----------------------------------------------------------------------------!
   
   !-----------------------------------------------------------------!
  !                                                                                                          !
   !    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                                         !
  !                                                                                                          !
   !-----------------------------------------------------------------!

Destruction of the concept  BC_Fix.
<S> ERREUR UTILISATEUR RECUPEREE PAR LE SUPERVISEUR


2. Would you please elaborate on your suggestion "Apply Velocity or acceleration to ball" ? 'Cause I have already set the initial velocity of the ball in the following command:

init_vel = CREA_CHAMP(TYPE_CHAM='NOEU_DEPL_R',
                        OPERATION='AFFE',
                        MODELE=model1,
                        AFFE=_F(GROUP_NO=('Particle_Solid', ),
                                NOM_CMP=('DX', 'DY', 'DZ'),
                                VALE=(0.0, 0.0, -1.765),),
                        INFO=1,)

solver = DYNA_NON_LINE(MODELE=model1,
                         CHAM_MATER=fieldma0,
                         CONTACT=contact,
                         COMPORTEMENT=_F( ... ),
                         ETAT_INIT=_F(VITE=init_vel,
                                      PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         INCREMENT=_F( ... ),
                         SCHEMA_TEMPS=_F( ... ),
                         SOLVEUR=_F( ... ),
                         AMOR_RAYL_RIGI='TANGENTE',
                         METHODE='NEWTON',
                         NEWTON=_F( ... ),
                         CONVERGENCE=_F( ... ),
                         MESURE=_F(TABLE='NON',),
                         ARCHIVAGE=_F( ... ),
                         INFO=1,)


3. How does setting up 'Liaison_Mail' help tackling the convergence issue of 'DYNA_NON_LINE' with 'VMIS_JOHN_COOK' enabled?

Thanks again!

Last edited by hkboondoggle (2020-07-31 09:25:30)

Offline

#8 2020-07-31 11:17:12

sameer21101970
Member
Registered: 2019-09-06
Posts: 177

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

odélisation :             !
   !  - mailles volumiques pour un modèle 3D                                     !
   !  - mailles                                                                                           !
   ! surfaciques pour un modèle 2D         


odelization:!
   ! - voluminal meshes for a 3D model!
   ! - stitches!
   ! surface for a 2D model

GROUP_MA_MAIT WILL BE 3D VOLUME
AND  SELECT GROUP_NO_ESCL = 2D SURFACE

Offline

#9 2020-07-31 16:38:39

sameer21101970
Member
Registered: 2019-09-06
Posts: 177

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

i hv tested, it is working under your given material law.

Offline

#10 2020-08-03 08:03:44

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

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

sameer, would you mind sharing the details of your test run?

'Cause I have managed to get 'Liaison_Mail' implemented in my latest trial, but it is still giving out  <MECANONLINE9_3> under 'VMIS_JOHN_COOK'.

Also, switching DEFORMATION of behavior of DYNA_NON_LINE from 'GROT_GDEP' to 'SIMO_MIEHE' still instantly leads to <F>_ABNORMAL_ABORT with EXECUTION_CODE_ASTER_EXIT_8627=0.


Attachments:
stderr_command_salome.log, Size: 446.85 KiB, Downloads: 7

Offline

#11 2020-08-03 12:36:34

sameer21101970
Member
Registered: 2019-09-06
Posts: 177

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

Whether you hv made changes to your geometry. No gap between ball and striking face.

Offline