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

You are not logged in.

#1 2020-06-15 07:07:29

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

[Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Hi all,

I am new to code_aster, and 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.

This is code of the case setting:

DEBUT(LANG='EN')

mesh = LIRE_MAILLAGE(identifier='0:1',
                     UNITE=20)

model1 = AFFE_MODELE(identifier='1:1',
                     AFFE=_F(MODELISATION=('3D', ),
                             PHENOMENE='MECANIQUE',
                             TOUT='OUI'),
                     MAILLAGE=mesh)

Cu_L_ISO = DEFI_MATERIAU(identifier='2:1',
                         ECRO_COOK=_F(A=0.0631779,
                                      B=0.283815,
                                      C=0.025,
                                      EPSP0=1e-06,
                                      M_PUIS=1.09,
                                      N_PUIS=0.31,
                                      TMELT=4.708333,
                                      TROOM=1.0),
                         ELAS=_F(ALPHA=5.903e-08,
                                 COEF_AMOR=0.01,
                                 E=106.92,
                                 NU=0.31,
                                 RHO=1.0),
                         THER=_F(LAMBDA=0.006422,
                                 RHO_CP=0.959))

fieldma0 = AFFE_MATERIAU(identifier='3:1',
                         AFFE=_F(MATER=(Cu_L_ISO, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh,
                         MODELE=model1)

listr = DEFI_LIST_REEL(identifier='4:1',
                       DEBUT=0.0,
                       INTERVALLE=_F(JUSQU_A=10.0,
                                     NOMBRE=2500),
                       TITRE='Time_step')

timeline = DEFI_LIST_INST(identifier='5:1',
                          DEFI_LIST=_F(LIST_INST=listr),
                          METHODE='MANUEL')

BC_Fix = AFFE_CHAR_MECA(identifier='6:1',
                        DDL_IMPO=_F(DX=0.0,
                                    DY=0.0,
                                    DZ=0.0,
                                    GROUP_NO=('Substrate_Bottom', )),
                        MODELE=model1)

contact = DEFI_CONTACT(identifier='7:1',
                       MODELE=model1,
                       ZONE=_F(APPARIEMENT='MAIT_ESCL',
                               GROUP_MA_ESCL=('Particle_Shell', ),
                               GROUP_MA_MAIT=('Substrate_Shell', ),
                               NORMALE='MAIT_ESCL',
                               VECT_ESCL='AUTO',
                               VECT_MAIT='AUTO'))

init_vel = CREA_CHAMP(identifier='8:1',
                      AFFE=_F(GROUP_MA=('Particle', ),
                              GROUP_NO=('Particle_Surface', ),
                              NOM_CMP=('DZ', ),
                              VALE=(-1.765, )),
                      MODELE=model1,
                      OPERATION='AFFE',
                      TYPE_CHAM='NOEU_DEPL_R')

solver = DYNA_NON_LINE(identifier='9:1',
                       CHAM_MATER=fieldma0,
                       CONTACT=contact,
                       ETAT_INIT=_F(VITE=init_vel),
                       INCREMENT=_F(LIST_INST=timeline,
                                    NUME_INST_FIN=10,
                                    NUME_INST_INIT=0),
                       MODELE=model1,
                       SCHEMA_TEMPS=_F(FORMULATION='DEPLACEMENT',
                                       MODI_EQUI='OUI',
                                       SCHEMA='HHT'))
FIN()

However, during execution of the study, the following error happened:

.. __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,
                                      NUME_INST_INIT=0,
                                      NUME_INST_FIN=10,
                                      PRECISION=1.E-06,),
                         SCHEMA_TEMPS=_F(SCHEMA='HHT',
                                         MODI_EQUI='OUI',
                                         ALPHA=-0.3,
                                         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 365355 éléments
     Relation                             : ELAS
     Déformation                          : PETIT
     Nombre total de variables internes   : 1
            V1 : VIDE
  Le système linéaire à résoudre a 1658619 degrés de liberté:
   - 1558995 sont des degrés de liberté physiques
     (ils sont
portés par 519665 noeuds du maillage)
   - 99624 sont les couples de paramètres de Lagrange associés
     aux 49812 relations
linéaires dualisées.
  La matrice est de taille 1658619 équations.
  Elle contient 67789758 termes non nuls si elle est symétrique et 133920897 termes
non nuls si elle n'est pas symétrique.
  Soit un taux de remplissage de   0.005 %.
    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.

   
   !---------------------------------------------------------------!
   ! <EXCEPTION> <FACTOR_78>                                       !
   !                                                               !
   ! Solver MUMPS:                                                 !
   !  Overflow of the terms of the matrix. One detected at least:  !
   !  - 593685 too small terms,                                    !
   !  - 5 too large terms.                                         !
   !---------------------------------------------------------------!
   
   !-----------------------------------------------------------------!
   !                                                                 !
   !    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                      !
   !                                                                 !
   !-----------------------------------------------------------------!


Criterion of destruction of the file ( 1.00%) associated with the base VOLATILE exceeded 7.89%
Nombre of records used:  4965

Occupied disc volume:  3879 Mo.
Maximum number of records:  62914


  Ouverture en écriture du fichier ./vola.1

Validation of the concept  solver.
<S> ERREUR UTILISATEUR RECUPEREE PAR LE SUPERVISEUR


Could you please explain the problem, and show me a solution to it?
Thanks!

Last edited by hkboondoggle (2020-06-29 04:55:28)


Attachments:
stderr_command_salome.log, Size: 66.67 KiB, Downloads: 70

Offline

#2 2020-06-15 17:42:49

Johannes_ACKVA
Member
From: Ingenieurbüro für Mechanik, DE
Registered: 2009-11-04
Posts: 701
Website

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

hello,
the field which you create and which should contain the initial velocity is uncomplete because you assigned only one component (DZ), but it is a vector field. So you must complete the components DX and DY. Example with any numbers (1111., 2222.):

init_vel = CREA_CHAMP(identifier='8:1',
                      AFFE=_F(GROUP_MA=('Particle', ),
                              GROUP_NO=('Particle_Surface', ),
                              NOM_CMP=('DX','DY','DZ', ),
                              VALE=(1111., 2222., -1.765, )),
                      MODELE=model1,
                      OPERATION='AFFE',
                      TYPE_CHAM='NOEU_DEPL_R')

Regards
Johannes_ACKVA

______________________________________________________________________
CODE-ASTER-courses at Ingenieurbüro für Mechanik, Germany

ONLINE-COURSES STARTING AT AUGUST 2020


Ingenieurbüro für Mechanik
D 91717 Wassertrüdingen / Germany

www.code-aster.de                                                Training & Support for NASTRAN and CODE-ASTER

Offline

#3 2020-06-16 04:18:06

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Hi Johannes, Thanks for your feedback.

I've corrected the the initial velocity field as instructed.

init_vel = CREA_CHAMP(identifier='8:1',
                      AFFE=_F(GROUP_MA=('Particle_Shell', 'Particle'),
                              GROUP_NO=('Particle_Surface', ),
                              NOM_CMP=('DX', 'DY', 'DZ'),
                              VALE=(0.0, 0.0, -1.765)),
                      MODELE=model1,
                      OPERATION='AFFE',
                      TYPE_CHAM='NOEU_DEPL_R')

However, the same problem of <EXCEPTION> <FACTOR_78> still persists. Could you please advice?

Thanks!

Offline

#4 2020-06-16 10:19:50

Johannes_ACKVA
Member
From: Ingenieurbüro für Mechanik, DE
Registered: 2009-11-04
Posts: 701
Website

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

hard to say what is the reason for the exception

You defined a nonlin material, but do a linear elastic analysis (no comportement defined, only DEFI_MATERIAU(ELAS=_F(E=106.92,NU=0.31,RHO=1.0)) matters), so it is already simple.

You can just try some modifications/simplifications to detect if there is any physical or numerical reason: add displacement constraints, eliminate the contact, do a static analysis with the model,..

Regards
Johannes_ACKVA

______________________________________________________________________
CODE-ASTER-courses at Ingenieurbüro für Mechanik, Germany

ONLINE-COURSES STARTING AT AUGUST 2020


Ingenieurbüro für Mechanik
D 91717 Wassertrüdingen / Germany

www.code-aster.de                                                Training & Support for NASTRAN and CODE-ASTER

Offline

#5 2020-06-17 06:58:36

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Johannes, would the exception be mesh related as well?
I'm wondering this because the log of the study has multiple warning about the mesh which I don't quite understand:

.. __stg1_cmd0:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0002 Concept of the type:  maillage_sdaster
  # ------------------------------------------------------------------------------------------
  mesh = LIRE_MAILLAGE(UNITE=20,
                       INFO_MED=1,
                       FORMAT='MED',
                       VERI_MAIL=_F(VERIF='OUI',
                                    APLAT=1.E-3,),
                       INFO=1,)

Vérification du maillage.
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M22 et M1 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M23 et M2 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M24 et M3 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M25 et M4 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !-----------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                              !
   !                                                                 !
   ! Les mailles M26 et M5 sont doubles (même noeuds en support).    !
   !     Type de maille : SEG3                                       !
   !                                                                 !
   !                                                                 !
   ! This is a warning. If you do not understand the meaning of this !
   !  warning, you can obtain unexpected results!                    !
   !-----------------------------------------------------------------!
   
   
   !--------------------------------------------------------------------------------------!
   ! <A> <CATAMESS_41>                                                                    !
   !                                                                                      !
   ! The warning message  MODELISA8_13 was emitted 5 time, it will not be displayed more. !
   !                                                                                      !
   !                                                                                      !
   ! This is a warning. If you do not understand the meaning of this                      !
   !  warning, you can obtain unexpected results!                                         !
   !--------------------------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------------------!
   ! <A> <MODELISA4_8>                                                             !
   !                                                                               !
   !  - > Phase of checking of the mesh: presence of double meshes (or triple,')   !
   !  The multiple meshes are meshes of different names                            !
   ! which have same connectivity                                                  !
   !  (they lean on the same nodes).                                               !
   !                                                                               !
   !  - > Risks & advices:                                                         !
   !  The risk is to modelize 2 times (or more)                                    !
   ! space. One can for example have                                               !
   !  a model 2 times too heavy or 2 times too rigid.                              !
   !  Note: the meshes concerned are printed in the                                !
   ! file" message".                                                               !
   !  On this mesh, it is imprudent to affect quantities with key word TOUT='OUI'. !
   !                                                                               !
   !                                                                               !
   ! This is a warning. If you do not understand the meaning of this               !
   !  warning, you can obtain unexpected results!                                  !
   !-------------------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M19392 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000429, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M20423 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000138, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M20942 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000007, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M21441 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000005, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_15>                                                !
   !                                                                   !
   ! La maille M22848 possède des noeuds confondus géométriquement.    !
   !     Rapport entre le plus petit côté et le plus grand = 0.000695, !
   ! type de maille : TETRA10                                          !
   !                                                                   !
   !                                                                   !
   ! This is a warning. If you do not understand the meaning of this   !
   !  warning, you can obtain unexpected results!                      !
   !-------------------------------------------------------------------!
   
   
   !--------------------------------------------------------------------------------------!
   ! <A> <CATAMESS_41>                                                                    !
   !                                                                                      !
   ! The warning message  MODELISA8_15 was emitted 5 time, it will not be displayed more. !
   !                                                                                      !
   !                                                                                      !
   ! This is a warning. If you do not understand the meaning of this                      !
   !  warning, you can obtain unexpected results!                                         !
   !--------------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------!
   ! <A> <MODELISA4_9>                                                    !
   !                                                                      !
   !  - > Phase of checking of the mesh: presence of degenerated meshes.  !
   !  The relationship between the smallest edge on largest is            !
   ! lower than 0,001                                                     !
   !                                                                      !
   !  - > Risks & advices:                                                !
   !  Check your mesh. The presence of such meshes can lead to            !
   !  problems of convergence and harm                                    !
   ! the quality of the results.                                          !
   !                                                                      !
   !                                                                      !
   ! This is a warning. If you do not understand the meaning of this      !
   !  warning, you can obtain unexpected results!                         !
   !----------------------------------------------------------------------!
   

------------ MAILLAGE mesh     - IMPRESSIONS NIVEAU  1 ------------

MA-16-JUIN-2020 09:31:03                                                       

NOMBRE DE NOEUDS                      519665

NOMBRE DE MAILLES                     365355
                              SEG3                  898
                              TRIA6                  28
                              QUAD8                8281
                              TETRA10            347867
                              PYRAM13              8281

NOMBRE DE GROUPES DE NOEUDS                6
                              GR_1_Particle_Surface             2341
                              Particle_Surface                  2341
                              GR_1_Substrate_Surface            4961
                              Substrate_Surface                 4961
                              GR_1_Substrate_Bottom             4961
                              Substrate_Bottom                  4961

NOMBRE DE GROUPES DE MAILLES               4
                              Substrate                       336925
                              Particle                         19223
                              Substrate_Shell                   7520
                              Particle_Shell                     789

--------------------------------------------------------------------------------

# Memory (Mo):    216.99 /    162.55 /    163.38 /    107.86 (VmPeak / VmSize / Optimum / Minimum)
  # Fin orders No:  0002 user+syst:          2.77S (syst:          0.16S, Elaps:          0.00S)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd1:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0004 Concept of the type:  modele_sdaster
  # ------------------------------------------------------------------------------------------
  model1 = AFFE_MODELE(MAILLAGE=mesh,
                       AFFE=_F(TOUT='OUI',
                               PHENOMENE='MECANIQUE',
                               MODELISATION=('3D', ),),
                       INFO=1,
                       DISTRIBUTION=_F(PARTITIONNEUR='METIS',
                                       METHODE='SOUS_DOMAINE',),
                       VERI_JACOBIEN='OUI',)


Sur les 365355 mailles du maillage mesh, on a demandé l'affectation de 365355, on a pu en affecter 365355
Modélisation     Type maille  Élément fini     Nombre
                  SEG3         MECA_ARETE3      898
                  TRIA6        MECA_FACE6       28
                  QUAD8        MECA_FACE8       8281
3D               PYRAM13      MECA_PYRAM13     8281
3D               TETRA10      MECA_TETRA10     347867
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9245 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9292 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9588 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9591 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !------------------------------------------------------------------!
   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9661 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.                                               !
   !                                                                  !
   !  Risks & advices:                                                !
   !  Was the mesh produced by a mesh generator?                      !
   !  Does                                                            !
   ! connectivity respect Aster convention strictly?                  !
   !                                                                  !
   !                                                                  !
   ! This is a warning. If you do not understand the meaning of this  !
   !  warning, you can obtain unexpected results!                     !
   !------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------------------------!
   ! <A> <CATAMESS_41>                                                                   !
   !                                                                                     !
   ! The warning message  MAILLAGE1_1 was emitted 5 time, it will not be displayed more. !
   !                                                                                     !
   !                                                                                     !
   ! This is a warning. If you do not understand the meaning of this                     !
   !  warning, you can obtain unexpected results!                                        !
   !-------------------------------------------------------------------------------------!

.. __stg1_cmd7:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0010 Concept of the type:  char_contact
  # ------------------------------------------------------------------------------------------
  contact = DEFI_CONTACT(MODELE=model1,
                         STOP_INTERP='NON',
                         LISSAGE='NON',
                         VERI_NORM='OUI',
                         ITER_GEOM_MAXI=10,
                         RESI_GEOM=1.E-2,
                         ALGO_RESO_GEOM='POINT_FIXE',
                         REAC_GEOM='AUTOMATIQUE',
                         ITER_CONT_MULT=4,
                         STOP_SINGULIER='OUI',
                         NB_RESOL=10,
                         ITER_GCP_MAXI=0,
                         RECH_LINEAIRE='ADMISSIBLE',
                         PRE_COND='SANS',
                         ZONE=_F(APPARIEMENT='MAIT_ESCL',
                                 GROUP_MA_MAIT='Substrate_Shell',
                                 GROUP_MA_ESCL='Particle_Shell',
                                 NORMALE='MAIT_ESCL',
                                 VECT_MAIT='AUTO',
                                 VECT_ESCL='AUTO',
                                 GLISSIERE='NON',
                                 TYPE_PROJECTION='ORTHOGONALE',
                                 DIST_POUTRE='NON',
                                 DIST_COQUE='NON',
                                 DIST_APPA=-1.0,
                                 TOLE_PROJ_EXT=0.5,
                                 RESOLUTION='OUI',
                                 ALGO_CONT='CONTRAINTE',),
                         FORMULATION='DISCRETE',
                         FROTTEMENT='SANS',)

T
T
   
   !------------------------------------------------------------------------------------------------------------------------------------!
   ! <A> <CONTACT_8>                                                                                                                    !
   !                                                                                                                                    !
   ! Formulation discrète du contact.                                                                                                   !
   !  -> Il y a des éléments de type QUAD8 sur les surfaces de contact. Ces éléments ne permettent pas                                  !
   ! de respecter exactement la condition de contact.                                                                                   !
   !     Afin d'empêcher une pénétration intempestive des surfaces, on a procédé à des                                                  !
   ! liaisons cinématiques (LIAISON_DDL) entre les noeuds milieux et les noeuds sommets, sur les deux surfaces (maître et esclave).     !
   !                                                                                                                                    !
   !                                                                                                                                    !
   ! Risques et conseils :                                                                                                              !
   !    - Ces liaisons supplémentaires peuvent provoquer des incompatibilités avec les conditions limites en                            !
   ! particulier dans le cas de symétries, ce qui se traduira par une matrice singulière.                                               !
   !      Dans ce cas, il est possible de                                                                                               !
   ! n'appliquer les conditions aux limites concernées que sur les noeuds sommets (on pourra utiliser la commande DEFI_GROUP pour créer !
   ! le groupe de noeuds adéquat).                                                                                                      !
   !    - Il est toutefois préférable d'utiliser sur les surfaces de contact des éléments de type QUAD9.                                !
   ! Pour cela, transformez les éléments volumiques HEXA20 en HEXA27 ou PENTA15 en PENTA18 (par exemple avec la commande                !
   ! CREA_MAILLAGE).                                                                                                                    !
   !                                                                                                                                    !
   !                                                                                                                                    !
   ! This is a warning. If you do not understand the meaning of this                                                                    !
   !  warning, you can obtain unexpected results!                                                                                       !
   !------------------------------------------------------------------------------------------------------------------------------------!

Offline

#6 2020-06-18 10:22:00

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Hello, could some one please help me out?

The links of the mesh file are attached FYI in case it helps the 'dissection'.
I am using mega.nz to host them since they exceed the size limit of attachment.

Mesh file: https: // mega.nz / file / HcoXHIxD#tp3AeMwVUqwFjUy2hEwiFzFT2v4aAWEU92Wb7nuR334  (38.7MB)

I am running code aster version 14.4.0 on Windows-10-10.0.17134-SP0.

Many thanks.

Offline

#7 2020-06-18 10:40:35

mecour
Member
From: Ostrava (Czech)
Registered: 2011-04-04
Posts: 158

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

There are doubled nodes and elements in your mesh. Try merge nodes and merge elements.

mecour

Offline

#8 2020-06-19 09:23:23

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Hi mecour, appreciate for your feedback!

I've merged the nodes (tolerance = 5e-4) as well as the mesh and gave it a shot.

Yet, the same problem of <EXCEPTION> <FACTOR_78> still persists. Any thoughts?

The Message Log of the latest run is attached FYI.

Thanks!


Attachments:
stderr_command_salome.log, Size: 54.29 KiB, Downloads: 65

Offline

#9 2020-06-19 17:59:44

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Hello,

Quality of your  mesh !

   ! <A> <MAILLAGE1_1>                                                !
   !                                                                  !
   !  Error of mesh:                                                  !
   !  The mesh M9212 of the type MECA_TETRA10 is distorted too much.  !
   !  The jacobian of the geometrical transformation                  !
   ! does not have the same sign on all                               !
   !  the Gauss points.     

For these kind of element: probably degenerated "flat" elements


Code_Asterの開発者

Offline

#10 2020-06-20 04:57:02

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Hello,
there are a few mistakes in your model and the main issue is that you're meshing solids AND faces. You must simply define groups in the geometry and not duplicate faces. Moreover be careful, you mesh has quadrangles only on the surface. I corrected a few points in your comm files, it's converging nicely.

Ball

Konyaro

Last edited by konyaro (2020-06-20 06:31:56)


Attachments:
ball.gif, Size: 712.03 KiB, Downloads: 325

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

Offline

#11 2020-06-20 04:58:18

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Attached comm and med.


Attachments:
ball.zip, Size: 731.2 KiB, Downloads: 82

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

Offline

#12 2020-06-20 11:15:26

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Excellent Konrayo.

Your Remark " Moreover be careful, you mesh has quadrangles only on the surface. I corrected a few points in your comm files  "

Quadrangles on surface for All type of Contact Analysis Studies.

Is above true.

Last edited by sameer21101970 (2020-06-20 11:18:32)

Offline

#13 2020-06-21 06:14:36

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Moreover be careful, you mesh has quadrangles only on the surface

I just wanted to point out that one may think that a Netgen mesh with the option "Quad-dominated" produces a nice hexahedral mesh. This is not the case, the elements are pyramids connected to tetrahedrons.
Fake Hexas

Last edited by konyaro (2020-06-21 06:15:04)


Attachments:
Mesh_fake_hexa.png, Size: 56.16 KiB, Downloads: 235

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

Offline

#14 2020-06-21 16:14:55

Anirudh
Member
Registered: 2017-09-11
Posts: 257

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

@hkboondogle
Please introduce yourself in the introduction section. Your username is not enough to understand what you do.

Regards
Anirudh

Offline

#15 2020-06-22 04:52:20

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

Wow, I wasn't expecting to receive such attention to my help request throughout the weekend.
Thanks for all of your enthusiastic feedback every one! I really appreciate it!

First of all, I have to admit that I'm a complete doolie for the Salome-Meca / Code_aster platform.
So during the geometry creation & meshing stage, I basically followed this video tutorial:
youtube .com / watch?v=eOvjusbQ_RE

Now, with the disclaimer done, may I inquire the following:

1.

konyaro wrote:

there are a few mistakes in your model and the main issue is that you're meshing solids AND faces. You must simply define groups in the geometry and not duplicate faces.

Konyaro, could you please clarify? Because aren't face mesh is required for "DEFI_CONTACT" ?

Or do you mean having one should assign the contact rule to the 3D elements of the ball & substrate instead?
I was approaching this way but encounter problems (see: viewtopic.php?id=24974). You did mention that "I need surfaces or edges in 3D"  for "DEFI_CONTACT" in that post. So it bears the question: how do one obtain surfaces or edges in 3D?


2.

AsterO'dactyle wrote:

For these kind of element: probably degenerated "flat" elements.

konyaro wrote:

you mesh has quadrangles only on the surface. I corrected a few points in your comm files.

I just wanted to point out that one may think that a Netgen mesh with the option "Quad-dominated" produces a nice hexahedral mesh. This is not the case, the elements are pyramids connected to tetrahedrons.

Given that, should I use Netgen ONLY for meshing the ball? Or should I avoid using Netgen at all?


@Anirudh: I'm attempting to simulate a spray particle colliding & attaching itself to a substrate. High strain rate deformation on both entities are expected upon the collision. I'll introduce myself there right away.


Many thanks guys!

Last edited by hkboondoggle (2020-06-22 04:52:44)

Offline

#16 2020-06-27 06:47:20

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

Re: [Solved]Advice needed on <EXCEPTION> <FACTOR_78> (DYNA_NON_LINE)

So it bears the question: how do one obtain surfaces or edges in 3D?

You must create groups, not new elements. The easiest way is to use "Create groups from geometry" in the MESH module. If you select a face it will create a group of 2D elements, if you select an edge it will create a group of 1D elements.

Given that, should I use Netgen ONLY for meshing the ball? Or should I avoid using Netgen at all?

You can use whatever you want. If you don't know what to choose I recommend quadratic tetrahedrons.

Konyaro


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

Offline