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

You are not logged in.

#1 2021-07-22 10:09:19

rmoortgat
Member
From: Belgium
Registered: 2016-08-30
Posts: 37
Website

SOLVED Simple Contact

I'm trying to simulate a simple contact between two solids.
I'm getting the following:
Tryed to figure it out but I'm stuck.
Any help would be appreciated.

Roger

Summary
Output messages for stage Stage_1


Output messages for stage Stage_1 - [top / bottom]

# ------------------------------------------------------------------------------------------
# Impression du contenu du fichier de commandes à exécuter :
# ------------------------------------------------------------------------------------------
DEBUT(LANG='FR')


Mup = LIRE_MAILLAGE(identifier='0:1',
                    FORMAT='MED',
                    UNITE=20)

Mlow = LIRE_MAILLAGE(identifier='1:1',
                     FORMAT='MED',
                     UNITE=2)

mesh = ASSE_MAILLAGE(identifier='2:1',
                     MAILLAGE_1=Mup,
                     MAILLAGE_2=Mlow,
                     OPERATION='SUPERPOSE')

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

steel = DEFI_MATERIAU(identifier='4:1',
                      ELAS=_F(E=210000.0,
                              NU=0.3,
                              RHO=7.8e-09))

fieldmat = AFFE_MATERIAU(identifier='5:1',
                         AFFE=_F(MATER=(steel, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh)

list = DEFI_LIST_REEL(identifier='6:1',
                      DEBUT=0.0,
                      INTERVALLE=_F(JUSQU_A=1.0,
                                    PAS=0.1))

func = DEFI_FONCTION(identifier='7:1',
                     NOM_PARA='INST',
                     PROL_DROITE='CONSTANT',
                     PROL_GAUCHE='CONSTANT',
                     VALE=(0.0, 0.0, 1.0, 1.0))

FixLeft = AFFE_CHAR_MECA(identifier='8:1',
                         DDL_IMPO=_F(DX=0.0,
                                     DY=0.0,
                                     DZ=0.0,
                                     GROUP_MA=('FixLeft', )),
                         MODELE=model)

FixRight = AFFE_CHAR_MECA(identifier='9:1',
                          DDL_IMPO=_F(DX=0.0,
                                      DY=0.0,
                                      DZ=0.0,
                                      GROUP_MA=('FixRight', )),
                          MODELE=model)

contact = DEFI_CONTACT(identifier='10:1',
                       FORMULATION='DISCRETE',
                       MODELE=model,
                       ZONE=_F(GROUP_MA_ESCL=('ContactLow', ),
                               GROUP_MA_MAIT=('ContactUp', )))

pressure = AFFE_CHAR_MECA(identifier='11:1',
                          FORCE_FACE=_F(FZ=-100.0,
                                        GROUP_MA=('pressure', )),
                          MODELE=model)

result = STAT_NON_LINE(identifier='12:1',
                       CHAM_MATER=fieldmat,
                       COMPORTEMENT=_F(DEFORMATION='PETIT',
                                       RELATION='ELAS',
                                       TOUT='OUI'),
                       CONTACT=contact,
                       CONVERGENCE=_F(ITER_GLOB_MAXI=10,
                                      RESI_GLOB_MAXI=0.0001,
                                      RESI_GLOB_RELA=0.0001),
                       EXCIT=(_F(CHARGE=FixLeft),
                              _F(CHARGE=FixRight),
                              _F(CHARGE=pressure,
                                 FONC_MULT=func)),
                       INCREMENT=_F(LIST_INST=list),
                       METHODE='NEWTON',
                       MODELE=model,
                       NEWTON=_F(MATRICE='TANGENTE',
                                 REAC_INCR=1,
                                 REAC_ITER=1),
                       SOLVEUR=_F(METHODE='MULT_FRONT'))

result = CALC_CHAMP(identifier='13:1',
                    reuse=result,
                    CONTRAINTE=('SIGM_ELNO', 'SIGM_NOEU'),
                    CRITERES=('SIEQ_ELNO', 'SIEQ_NOEU'),
                    FORCE=('REAC_NODA', ),
                    RESULTAT=result)

IMPR_RESU(identifier='14:1',
          FORMAT='MED',
          RESU=_F(MAILLAGE=mesh,
                  NOM_CHAM=('DEPL', 'SIGM_ELNO', 'SIEQ_NOEU'),
                  RESULTAT=result),
          UNITE=80)

IMPR_RESU(identifier='15:1',
          FORMAT='RESULTAT',
          RESU=_F(MAILLAGE=mesh,
                  NOM_CHAM=('DEPL', 'SIGM_ELNO', 'SIEQ_NOEU'),
                  RESULTAT=result),
          UNITE=8)

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

  Ouverture en écriture du fichier ./vola.1

<INFO> Démarrage de l'exécution.
                                                                               
                                                                               
              -- CODE_ASTER -- VERSION : EXPLOITATION (stable) --               
                                                                               
                     Version 14.6.0 modifiée le 11/06/2020                     
                     révision b3490fa3b76c - branche 'v14'                     
                         Copyright EDF R&D 1991 - 2021                         
                                                                               
                    Exécution du : Thu Jun 17 09:42:39 2021                     
                       Nom de la machine : roger-W65-67SZ                       
                              Architecture : 64bit                             
                          Type de processeur : x86_64                           
Système d'exploitation : Linux-5.4.0-74-generic-x86_64-with-debian-bullseye-sid
                        Langue des messages : fr (UTF-8)                       
                                                                               
                           Version de Python : 3.6.5                           
                           Version de NumPy : 1.15.1                           
                           Parallélisme MPI : inactif                           
                          Parallélisme OpenMP : actif                           
                        Nombre de processus utilisés : 2                       
                     Version de la librairie HDF5 : 1.10.3                     
                      Version de la librairie MED : 4.0.0                       
                     Version de la librairie MFront : 3.2.1                     
                     Version de la librairie MUMPS : 5.1.2                     
                        Librairie PETSc : non disponible                       
                     Version de la librairie SCOTCH : 6.0.4                     

                  Mémoire limite pour l'exécution : 2048.00 Mo                 
                   consommée par l'initialisation : 451.18 Mo                   
                  par les objets du jeu de commandes : 1.30 Mo                 
                 reste pour l'allocation dynamique : 1595.39 Mo                 
                Taille limite des fichiers d'échange : 48.00 Go                 
                                                                               

.. __stg1_txt1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0001            Concept de type : -
  # ------------------------------------------------------------------------------------------
  DEBUT(LANG='FR',
        PAR_LOT='OUI',
        IMPR_MACRO='NON',
        DEBUG=_F(JXVERI='NON',
                 SDVERI='NON',
                 JEVEUX='NON',
                 HIST_ETAPE='NON',),
        MESURE_TEMPS=_F(NIVE_DETAIL=1,
                        MOYENNE='NON',),
        MEMOIRE=_F(TAILLE_BLOC=800.0,
                   TAILLE_GROUP_ELEM=1000,),
        RESERVE_CPU=_F(BORNE=900,),
        INFO=1,)


  Valeur initiale du temps CPU maximum =   900 secondes
  Valeur du temps CPU maximum passé aux commandes =   810 secondes
  Réserve CPU prévue = 90 secondes

  Ouverture en écriture du fichier ./glob.1

  Ouverture en écriture du fichier ./vola.1

  Ouverture en lecture du fichier
/home/roger/salome_meca_2020/V2020.0.1_universal_universal/tools/Code_aster_stable-1460/lib/aster/elem.1

     Nom de la base                          :  ELEMBASE
     Créée avec la version                   :  14.06.00
     Nombre d'enregistrements utilisés       :  43
     Nombre d'enregistrements maximum        :  512
     Longueur d'enregistrement (octets)      :  819200
     Nombre d'identificateurs utilisés       :  123
     Taille maximum du répertoire            :  300
     Pourcentage d'utilisation du répertoire :  41 %

  Ouverture en lecture du fichier
/home/roger/salome_meca_2020/V2020.0.1_universal_universal/tools/Code_aster_stable-1460/lib/aster/elem.1


     Nom de la base                          :  ELEMBASE
     Nombre d'enregistrements utilisés       :  43
     Nombre d'enregistrements maximum        :  512
     Longueur d'enregistrement (octets)      :  819200
     Nombre total d'accès en lecture         :  60
     Volume des accès en lecture             :         46.88 Mo.
     Nombre total d'accès en écriture        :  0
     Volume des accès en écriture            :          0.00 Mo.
     Nombre d'identificateurs utilisés       :  123
     Taille maximum du répertoire            :  300
     Pourcentage d'utilisation du répertoire :  41 %

  Relecture des catalogues des éléments faite.

  Fin de lecture (durée  0.544247  s.)
  # Mémoire (Mo) :   497.35 /   489.65 /    44.29 /    22.92 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0001   user+syst:        0.12s (syst:        0.08s, elaps:        0.58s)
  # ------------------------------------------------------------------------------------------

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

Vérification du maillage.

------------ MAILLAGE Mup      - IMPRESSIONS NIVEAU  1 ------------

JE-17-JUIN-2021 09:42:39                                                       

NOMBRE DE NOEUDS                         712

NOMBRE DE MAILLES                        870
                              SEG2                  150
                              TRIA3                  20
                              QUAD4                 700

NOMBRE DE GROUPES DE NOEUDS                3
                              pressure                           121
                              ContactUp                          121
                              FixLeft                            121

NOMBRE DE GROUPES DE MAILLES               3
                              pressure                           100
                              FixLeft                            100
                              ContactUp                          100

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

  # Mémoire (Mo) :   497.35 /   490.70 /    44.29 /    22.92 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0002   user+syst:        0.03s (syst:        0.01s, elaps:        0.20s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd1:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0004            Concept de type : maillage_sdaster
  # ------------------------------------------------------------------------------------------
  Mlow = LIRE_MAILLAGE(FORMAT='MED',
                       UNITE=2,
                       INFO_MED=1,
                       VERI_MAIL=_F(VERIF='OUI',
                                    APLAT=1.E-3,),
                       INFO=1,)

Vérification du maillage.

------------ MAILLAGE Mlow     - IMPRESSIONS NIVEAU  1 ------------

JE-17-JUIN-2021 09:42:40                                                       

NOMBRE DE NOEUDS                         712

NOMBRE DE MAILLES                        870
                              SEG2                  150
                              TRIA3                  20
                              QUAD4                 700

NOMBRE DE GROUPES DE NOEUDS                2
                              FixRight                           121
                              ContactLow                         121

NOMBRE DE GROUPES DE MAILLES               2
                              FixRight                           100
                              ContactLow                         100

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

  # Mémoire (Mo) :   497.35 /   490.83 /    44.29 /    22.92 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0004   user+syst:        0.02s (syst:        0.00s, elaps:        0.01s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd2:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0006            Concept de type : maillage_sdaster
  # ------------------------------------------------------------------------------------------
  mesh = ASSE_MAILLAGE(MAILLAGE_1=Mup,
                       MAILLAGE_2=Mlow,
                       OPERATION='SUPERPOSE',)

  # Mémoire (Mo) :   497.35 /   490.99 /    44.29 /    22.92 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0006   user+syst:        0.01s (syst:        0.00s, elaps:        0.05s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd3:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0007            Concept de type : modele_sdaster
  # ------------------------------------------------------------------------------------------
  model = 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 1740 mailles du maillage mesh, on a demandé l'affectation de 1740, on a pu en affecter 1740
   
   !-------------------------------------------------------------------------------------!
   ! <A> <MODELE1_63>                                                                    !
   !                                                                                     !
   !   -> La maille M1 porte un élément fini de bord, mais elle ne borde                 !
   !      aucun élément ayant une "rigidité".                                            !
   !                                                                                     !
   !   -> Risque & Conseil :                                                             !
   !      Cela peut entraîner des problèmes de "pivot nul" lors de la résolution.        !
   !      Si la résolution des systèmes linéaires ne pose pas de problèmes, vous         !
   !      pouvez ignorer ce message.                                                     !
   !      Sinon, vérifier la définition du modèle (AFFE_MODELE) en évitant l'utilisation !
   !      de l'opérande TOUT='OUI'.                                                      !
   !                                                                                     !
   !                                                                                     !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                      !
   ! alarme, vous pouvez obtenir des résultats inattendus !                              !
   !-------------------------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------------------------!
   ! <A> <MODELE1_63>                                                                    !
   !                                                                                     !
   !   -> La maille M2 porte un élément fini de bord, mais elle ne borde                 !
   !      aucun élément ayant une "rigidité".                                            !
   !                                                                                     !
   !   -> Risque & Conseil :                                                             !
   !      Cela peut entraîner des problèmes de "pivot nul" lors de la résolution.        !
   !      Si la résolution des systèmes linéaires ne pose pas de problèmes, vous         !
   !      pouvez ignorer ce message.                                                     !
   !      Sinon, vérifier la définition du modèle (AFFE_MODELE) en évitant l'utilisation !
   !      de l'opérande TOUT='OUI'.                                                      !
   !                                                                                     !
   !                                                                                     !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                      !
   ! alarme, vous pouvez obtenir des résultats inattendus !                              !
   !-------------------------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------------------------!
   ! <A> <MODELE1_63>                                                                    !
   !                                                                                     !
   !   -> La maille M3 porte un élément fini de bord, mais elle ne borde                 !
   !      aucun élément ayant une "rigidité".                                            !
   !                                                                                     !
   !   -> Risque & Conseil :                                                             !
   !      Cela peut entraîner des problèmes de "pivot nul" lors de la résolution.        !
   !      Si la résolution des systèmes linéaires ne pose pas de problèmes, vous         !
   !      pouvez ignorer ce message.                                                     !
   !      Sinon, vérifier la définition du modèle (AFFE_MODELE) en évitant l'utilisation !
   !      de l'opérande TOUT='OUI'.                                                      !
   !                                                                                     !
   !                                                                                     !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                      !
   ! alarme, vous pouvez obtenir des résultats inattendus !                              !
   !-------------------------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------------------------!
   ! <A> <MODELE1_63>                                                                    !
   !                                                                                     !
   !   -> La maille M4 porte un élément fini de bord, mais elle ne borde                 !
   !      aucun élément ayant une "rigidité".                                            !
   !                                                                                     !
   !   -> Risque & Conseil :                                                             !
   !      Cela peut entraîner des problèmes de "pivot nul" lors de la résolution.        !
   !      Si la résolution des systèmes linéaires ne pose pas de problèmes, vous         !
   !      pouvez ignorer ce message.                                                     !
   !      Sinon, vérifier la définition du modèle (AFFE_MODELE) en évitant l'utilisation !
   !      de l'opérande TOUT='OUI'.                                                      !
   !                                                                                     !
   !                                                                                     !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                      !
   ! alarme, vous pouvez obtenir des résultats inattendus !                              !
   !-------------------------------------------------------------------------------------!
   
   
   !-------------------------------------------------------------------------------------!
   ! <A> <MODELE1_63>                                                                    !
   !                                                                                     !
   !   -> La maille M5 porte un élément fini de bord, mais elle ne borde                 !
   !      aucun élément ayant une "rigidité".                                            !
   !                                                                                     !
   !   -> Risque & Conseil :                                                             !
   !      Cela peut entraîner des problèmes de "pivot nul" lors de la résolution.        !
   !      Si la résolution des systèmes linéaires ne pose pas de problèmes, vous         !
   !      pouvez ignorer ce message.                                                     !
   !      Sinon, vérifier la définition du modèle (AFFE_MODELE) en évitant l'utilisation !
   !      de l'opérande TOUT='OUI'.                                                      !
   !                                                                                     !
   !                                                                                     !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                      !
   ! alarme, vous pouvez obtenir des résultats inattendus !                              !
   !-------------------------------------------------------------------------------------!
   
   
   !------------------------------------------------------------------------------!
   ! <A> <CATAMESS_41>                                                            !
   !                                                                              !
   ! Le message d'alarme 'MODELE1_63' a été émis 5 fois, il ne sera plus affiché. !
   !                                                                              !
   !                                                                              !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette               !
   ! alarme, vous pouvez obtenir des résultats inattendus !                       !
   !------------------------------------------------------------------------------!
   
   
   !--------------------------------------------------------------------------------------------------!
   ! <A> <MODELE1_64>                                                                                 !
   !                                                                                                  !
   !   -> Le modèle model n'a pas d'éléments sachant calculer la rigidité.                            !
   !                                                                                                  !
   !   -> Risque & Conseil :                                                                          !
   !      Ce modèle ne pourra donc pas (en général) être utilisé pour faire des calculs.              !
   !      Vérifier la définition du modèle (AFFE_MODELE) et assurez-vous que les                      !
   !      types de mailles du maillage (SEG2, TRIA3, QUAD4, ...) sont compatibles avec votre          !
   !      modélisation.                                                                               !
   !      Exemples d'erreur :                                                                         !
   !        * affecter une modélisation "3D" sur un maillage formé de facettes.                       !
   !        * affecter une modélisation qui ne sait pas traiter tous les types de mailles du maillage !
   !          (par exemple 'PLAN_DIAG' en thermique, 'AXIS_SI' en mécanique)                          !
   !                                                                                                  !
   !                                                                                                  !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                                   !
   ! alarme, vous pouvez obtenir des résultats inattendus !                                           !
   !--------------------------------------------------------------------------------------------------!
   
Modélisation     Type maille  Élément fini     Nombre
                  SEG2         MECA_ARETE2      300
                  TRIA3        MECA_FACE3       40
                  QUAD4        MECA_FACE4       1400
#2      Calculs elementaires et assemblages               CPU (USER+SYST/SYST/ELAPS):      0.00      0.00      0.00
  # Mémoire (Mo) :   524.41 /   524.41 /    44.29 /    33.96 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0007   user+syst:        0.11s (syst:        0.06s, elaps:        0.17s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd4:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0008            Concept de type : mater_sdaster
  # ------------------------------------------------------------------------------------------
  steel = DEFI_MATERIAU(ELAS=_F(E=2.10000E5,
                                NU=0.3,
                                RHO=7.8E-09,
                                COEF_AMOR=1.0,),
                        INFO=1,)

Propriétés définies dans le matériau: ELAS.
  # Mémoire (Mo) :   524.41 /   524.41 /    44.29 /    33.96 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0008   user+syst:        0.02s (syst:        0.00s, elaps:        0.04s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd5:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0009            Concept de type : cham_mater
  # ------------------------------------------------------------------------------------------
  fieldmat = AFFE_MATERIAU(MAILLAGE=mesh,
                           AFFE=_F(TOUT='OUI',
                                   MATER=(steel, ),),
                           INFO=1,)

  # Mémoire (Mo) :   524.41 /   524.41 /    44.29 /    33.96 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0009   user+syst:        0.01s (syst:        0.00s, elaps:        0.03s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd6:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0010            Concept de type : listr8_sdaster
  # ------------------------------------------------------------------------------------------
  list = DEFI_LIST_REEL(DEBUT=0.0,
                        INTERVALLE=_F(JUSQU_A=1.0,
                                      PAS=0.1,),
                        INFO=1,)

  # Mémoire (Mo) :   524.41 /   524.41 /    44.29 /    33.96 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0010   user+syst:        0.00s (syst:        0.00s, elaps:        0.01s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd7:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0011            Concept de type : fonction_sdaster
  # ------------------------------------------------------------------------------------------
  func = DEFI_FONCTION(NOM_PARA='INST',
                       VALE=(0.0, 0.0, 1.0, 1.0),
                       PROL_DROITE='CONSTANT',
                       PROL_GAUCHE='CONSTANT',
                       NOM_RESU='TOUTRESU',
                       INTERPOL=('LIN', ),
                       VERIF='CROISSANT',
                       INFO=1,)

  # Mémoire (Mo) :   524.41 /   524.41 /    44.29 /    33.96 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0011   user+syst:        0.00s (syst:        0.00s, elaps:        0.00s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd8:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0012            Concept de type : char_meca
  # ------------------------------------------------------------------------------------------
  FixLeft = AFFE_CHAR_MECA(MODELE=model,
                           DDL_IMPO=_F(GROUP_MA=('FixLeft', ),
                                       DX=0.0,
                                       DY=0.0,
                                       DZ=0.0,),
                           VERI_NORM='OUI',
                           DOUBLE_LAGRANGE='OUI',
                           INFO=1,)

  # Mémoire (Mo) :   524.88 /   524.63 /    44.29 /    33.96 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0012   user+syst:        0.02s (syst:        0.01s, elaps:        0.05s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd9:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0013            Concept de type : char_meca
  # ------------------------------------------------------------------------------------------
  FixRight = AFFE_CHAR_MECA(MODELE=model,
                            DDL_IMPO=_F(GROUP_MA=('FixRight', ),
                                        DX=0.0,
                                        DY=0.0,
                                        DZ=0.0,),
                            VERI_NORM='OUI',
                            DOUBLE_LAGRANGE='OUI',
                            INFO=1,)

  # Mémoire (Mo) :   525.09 /   524.85 /    44.29 /    33.96 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0013   user+syst:        0.01s (syst:        0.00s, elaps:        0.01s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd10:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0014            Concept de type : char_contact
  # ------------------------------------------------------------------------------------------
  contact = DEFI_CONTACT(MODELE=model,
                         FORMULATION='DISCRETE',
                         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(GROUP_MA_MAIT='ContactUp',
                                 GROUP_MA_ESCL='ContactLow',
                                 GLISSIERE='NON',
                                 APPARIEMENT='MAIT_ESCL',
                                 NORMALE='MAIT',
                                 VECT_MAIT='AUTO',
                                 VECT_ESCL='AUTO',
                                 TYPE_PROJECTION='ORTHOGONALE',
                                 DIST_POUTRE='NON',
                                 DIST_COQUE='NON',
                                 DIST_APPA=-1.0,
                                 TOLE_PROJ_EXT=0.5,
                                 RESOLUTION='OUI',
                                 ALGO_CONT='CONTRAINTE',),
                         FROTTEMENT='SANS',)

  # Mémoire (Mo) :   525.09 /   524.93 /    44.29 /    33.96 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0014   user+syst:        0.02s (syst:        0.00s, elaps:        0.04s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd11:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0015            Concept de type : char_meca
  # ------------------------------------------------------------------------------------------
  pressure = AFFE_CHAR_MECA(MODELE=model,
                            FORCE_FACE=_F(GROUP_MA=('pressure', ),
                                          FZ=-100.0,),
                            VERI_NORM='OUI',
                            DOUBLE_LAGRANGE='OUI',
                            INFO=1,)

  # Mémoire (Mo) :   525.14 /   524.95 /    44.29 /    33.96 (VmPeak / VmSize / Optimum / Minimum)
  # Fin commande No : 0015   user+syst:        0.01s (syst:        0.00s, elaps:        0.04s)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd12:1
  # ------------------------------------------------------------------------------------------
  # Commande No :  0016            Concept de type : evol_noli
  # ------------------------------------------------------------------------------------------
  result = STAT_NON_LINE(MODELE=model,
                         CHAM_MATER=fieldmat,
                         CONTACT=contact,
                         METHODE='NEWTON',
                         EXCIT=(_F(CHARGE=FixLeft,
                                   TYPE_CHARGE='FIXE_CSTE',),
                                _F(CHARGE=FixRight,
                                   TYPE_CHARGE='FIXE_CSTE',),
                                _F(CHARGE=pressure,
                                   FONC_MULT=func,
                                   TYPE_CHARGE='FIXE_CSTE',),
                                ),
                         COMPORTEMENT=_F(TOUT='OUI',
                                         RELATION='ELAS',
                                         DEFORMATION='PETIT',
                                         RESI_CPLAN_RELA=1.E-06,
                                         RESI_INTE_RELA=1.E-06,
                                         ITER_INTE_MAXI=20,
                                         ITER_INTE_PAS=0,
                                         ITER_CPLAN_MAXI=1,
                                         SYME_MATR_TANG='OUI',
                                         PARM_THETA=1.0,),
                         INCREMENT=_F(LIST_INST=list,
                                      PRECISION=1.E-06,),
                         CONVERGENCE=_F(RESI_GLOB_MAXI=1.E-4,
                                        RESI_GLOB_RELA=1.E-4,
                                        ITER_GLOB_MAXI=10,
                                        ITER_GLOB_ELAS=25,
                                        ARRET='OUI',),
                         SOLVEUR=_F(METHODE='MULT_FRONT',
                                    RENUM='METIS',
                                    NPREC=8,
                                    ELIM_LAGR='NON',
                                    STOP_SINGULIER='OUI',),
                         NEWTON=_F(REAC_INCR=1,
                                   MATRICE='TANGENTE',
                                   REAC_ITER=1,
                                   REAC_ITER_ELAS=0,
                                   MATR_RIGI_SYME='NON',),
                         MESURE=_F(TABLE='NON',),
                         ARCHIVAGE=_F(PRECISION=1.E-06,
                                      CRITERE='RELATIF',),
                         )

   
   !-------------------------------------------------------------------------------------------------------------!
   ! <A> <CALCULEL_54>                                                                                           !
   !                                                                                                             !
   !  Le paramètre:  PCOMPOR  de l'option:  FULL_MECA  n'est pas connu des TYPE_ELEM du LIGREL:  model   .MODELE !
   !  Champ : &&CRCMEL.CEL1                                                                                      !
   !                                                                                                             !
   !                                                                                                             !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                                              !
   ! alarme, vous pouvez obtenir des résultats inattendus !                                                      !
   !-------------------------------------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <MECANONLINE_3>                                                                                    !
   !                                                                                                                !
   ! Erreur utilisateur :                                                                                           !
   !   Vous essayez de faire un calcul non-linéaire mécanique ou un post-traitement sur un modèle dont les éléments !
   !   ne sont pas programmés pour cela.                                                                            !
   !   On arrête le calcul.                                                                                         !
   ! Risques & conseils :                                                                                           !
   !   Vous devriez changer de modélisation.                                                                        !
   !----------------------------------------------------------------------------------------------------------------!
   

  Validation du concept 'result'.
<S> ERREUR UTILISATEUR RECUPEREE PAR LE SUPERVISEUR

.. __stg1_txt111
  # ------------------------------------------------------------------------------------------
  # Commande No :  0017            Concept de type : -
  # ------------------------------------------------------------------------------------------
  FIN(FORMAT_HDF='NON',
      RETASSAGE='NON',
      INFO_RESU='OUI',
      PROC0='OUI',
      )

   
   !-------------------------------------------------------------------------!
   ! <I> <CATAMESS_89>                                                       !
   !                                                                         !
   !     Liste des alarmes émises lors de l'exécution du calcul.             !
   !                                                                         !
   !     Les alarmes que vous avez choisies d'ignorer sont précédées de (*). !
   !     Nombre d'occurrences pour chacune des alarmes :                     !
   !            CALCULEL_54          émise    1 fois                         !
   !            MODELE1_63           émise 1740 fois                         !
   !            MODELE1_64           émise    1 fois                         !
   !-------------------------------------------------------------------------!
   
---------------------------------------------------------------------------
Concepts de la base: G
    Nom       Type                 Taille (Mo)         Nombre      Nombre de
                                                     d'objets       segments
    TOTAL                                28.30            266            465
 
    Mup        MAILLAGE_SDASTER           0.10             36             63
    Mlow       MAILLAGE_SDASTER           0.10             36             61
    mesh       MAILLAGE_SDASTER           0.20             35             66
    model      MODELE_SDASTER             0.11              8             13
    steel      MATER_SDASTER              0.00              5              5
    fieldmat   CHAM_MATER                 0.00              6             11
    list       LISTR8_SDASTER             0.00              6              6
    func       FONCTION_SDASTER           0.00              4              4
    FixLeft    CHAR_MECA                  0.19             23             43
    FixRight   CHAR_MECA                  0.19             23             43
    contact    CHAR_CONTACT               0.09             29             29
    pressure   CHAR_MECA                  0.01              8             13
    result     EVOL_NOLI                  0.00              2              2
    &FOZERO                               0.00              2              2
    &&_NUM_C                              0.00              1              1
    &CATA.AC                              0.00              2              4
    &CATA.CL                              0.60              1              3
    &CATA.GD                              0.14              4             11
    &CATA.ME                              0.20              2              4
    &CATA.OP                              0.30              4             19
    &CATA.PH                              0.00              1              1
    &CATA.PR                              0.00              2              4
    &CATA.TE                             26.03             17             42
    &CATA.TH                              0.01              2              4
    &CATA.TM                              0.01              7             11
---------------------------------------------------------------------------


     Nom de la base                          :  GLOBALE
     Nombre d'enregistrements utilisés       :  47
     Nombre d'enregistrements maximum        :  62914
     Longueur d'enregistrement (octets)      :  819200
     Nombre total d'accès en lecture         :  7
     Volume des accès en lecture             :          5.47 Mo.
     Nombre total d'accès en écriture        :  63
     Volume des accès en écriture            :         49.22 Mo.
     Nombre d'identificateurs utilisés       :  459
     Taille maximum du répertoire            :  2000
     Pourcentage d'utilisation du répertoire :  22 %


     Nom de la base                          :  VOLATILE
     Nombre d'enregistrements utilisés       :  3
     Nombre d'enregistrements maximum        :  62914
     Longueur d'enregistrement (octets)      :  819200
     Nombre total d'accès en lecture         :  0
     Volume des accès en lecture             :          0.00 Mo.
     Nombre total d'accès en écriture        :  6
     Volume des accès en écriture            :          4.69 Mo.
     Nombre d'identificateurs utilisés       :  305
     Taille maximum du répertoire            :  2000
     Pourcentage d'utilisation du répertoire :  15 %

<FIN> Arrêt normal dans "FIN".
<I> <FIN> ARRET NORMAL DANS "FIN" PAR APPEL A "JEFINI".

La mémoire demandée au lancement est surestimée, elle est de 2048.00 Mo.
Le pic mémoire utilisée est de 526.26 Mo.



     Nom de la base                          :  ELEMBASE
     Nombre d'enregistrements utilisés       :  43
     Nombre d'enregistrements maximum        :  512
     Longueur d'enregistrement (octets)      :  819200
     Nombre total d'accès en lecture         :  60
     Volume des accès en lecture             :         46.88 Mo.
     Nombre total d'accès en écriture        :  0
     Volume des accès en écriture            :          0.00 Mo.
     Nombre d'identificateurs utilisés       :  123
     Taille maximum du répertoire            :  300
     Pourcentage d'utilisation du répertoire :  41 %
 
<I> <FIN> MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION :                      33.96 Mo
<I> <FIN> MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION :                      44.29 Mo
<I> <FIN> MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS LORS DE L'EXECUTION :     526.26 Mo
 
<I>       FERMETURE DES BASES EFFECTUEE
 
   STATISTIQUES CONCERNANT L'ALLOCATION DYNAMIQUE :
     TAILLE CUMULEE MAXIMUM            :                   44  Mo.
     TAILLE CUMULEE LIBEREE            :                    7  Mo.
     NOMBRE TOTAL D'ALLOCATIONS        :                11391
     NOMBRE TOTAL DE LIBERATIONS       :                11371
     APPELS AU MECANISME DE LIBERATION :                    0
     TAILLE MEMOIRE CUMULEE RECUPEREE  :                    0  Mo.
     VOLUME DES LECTURES               :                    0  Mo.
     VOLUME DES ECRITURES              :                    0  Mo.
 
   MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION :      33.96 Mo
     - IMPOSE DE NOMBREUX ACCES DISQUE
     - RALENTIT LA VITESSE D'EXECUTION
   MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION :      44.29 Mo
     - LIMITE LES ACCES DISQUE
     - AMELIORE LA VITESSE D'EXECUTION
   MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS     :     526.26 Mo
     - COMPREND LA MEMOIRE CONSOMMEE PAR  JEVEUX,
       LE SUPERVISEUR PYTHON, LES LIBRAIRIES EXTERNES
 
<I>       FIN D'EXECUTION LE : JE-17-JUIN-2021 09:42:40


  <I> Informations sur les temps d'exécution
      Temps cpu user total                    0.76 s
      Temps cpu systeme total                 0.22 s
      Temps elapsed total                     2.56 s
      Temps restant                         807.44 s


  La signature de la base sauvegardée est (à l'adresse 0) :
    abbcdcc2fe9c9a3f56f95063f8a7933d06adf522

********************************************************************************
* COMMAND                  :       USER :     SYSTEM :   USER+SYS :    ELAPSED *
********************************************************************************
* init (jdc)               :       0.49 :       0.04 :       0.53 :       1.00 *
*  . compile               :       0.00 :       0.00 :       0.00 :       0.01 *
*  . exec_compile          :       0.14 :       0.01 :       0.15 :       0.45 *
*  . report                :       0.01 :       0.00 :       0.01 :       0.00 *
*  . build                 :       0.00 :       0.00 :       0.00 :       0.00 *
* DEBUT                    :       0.04 :       0.08 :       0.12 :       0.58 *
* LIRE_MAILLAGE            :       0.02 :       0.01 :       0.03 :       0.20 *
* LIRE_MAILLAGE            :       0.02 :       0.00 :       0.02 :       0.01 *
* ASSE_MAILLAGE            :       0.01 :       0.00 :       0.01 :       0.05 *
* AFFE_MODELE              :       0.05 :       0.06 :       0.11 :       0.17 *
* DEFI_MATERIAU            :       0.02 :       0.00 :       0.02 :       0.04 *
* AFFE_MATERIAU            :       0.01 :       0.00 :       0.01 :       0.03 *
* DEFI_LIST_REEL           :       0.00 :       0.00 :       0.00 :       0.01 *
* DEFI_FONCTION            :       0.00 :       0.00 :       0.00 :       0.00 *
* AFFE_CHAR_MECA           :       0.01 :       0.01 :       0.02 :       0.05 *
* AFFE_CHAR_MECA           :       0.01 :       0.00 :       0.01 :       0.01 *
* DEFI_CONTACT             :       0.02 :       0.00 :       0.02 :       0.04 *
* AFFE_CHAR_MECA           :       0.01 :       0.00 :       0.01 :       0.04 *
* STAT_NON_LINE            :       0.05 :       0.02 :       0.07 :       0.15 *
* FIN                      :       0.03 :       0.02 :       0.05 :       0.05 *
*  . part Superviseur      :       0.57 :       0.12 :       0.69 :       1.80 *
*  . part Fortran          :       0.20 :       0.10 :       0.30 :       0.77 *
********************************************************************************
* TOTAL_JOB                :       0.76 :       0.22 :       0.98 :       2.56 *
********************************************************************************

JDC.py : ERREUR A L'EXECUTION - INTERRUPTION
>> JDC.py : DEBUT RAPPORT
CR d'execution de JDC en MIXTE
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee.                          !
   ! Les bases sont fermees.                                                        !
   ! Type de l'exception : error                                                    !
   !                                                                                !
   ! Erreur utilisateur :                                                           !
   !   Vous essayez de faire un calcul non-linéaire mécanique ou un post-traitement !
   ! sur un modèle dont les éléments                                                !
   !   ne sont pas programmés pour cela.                                            !
   !   On arrête le calcul.                                                         !
   ! Risques & conseils :                                                           !
   !   Vous devriez changer de modélisation.                                        !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR d'execution de JDC en MIXTE

>> JDC.py : FIN RAPPORT
EXECUTION_CODE_ASTER_EXIT_3584-0001-roger-W65-67SZ=1

<I>_EXIT_CODE = 1

Last edited by rmoortgat (2021-07-26 08:07:40)


Attachments:
SimpleContact.tar.gz, Size: 12.83 KiB, Downloads: 6

Offline

#2 2021-07-22 10:45:36

mf
Member
Registered: 2019-06-18
Posts: 227

Re: SOLVED Simple Contact

Hello Roger,

your two bodies are hollow (only 2D surface meshes) and you assigned 3D elements to them, that doesn't work.

If this should be a contact between two 3D bodies, they'll need a 3D mesh (3D meshing of solids, not only 2D mesh of their faces).

Mario.

Last edited by mf (2021-07-22 10:46:02)


Attachments:
Bildschirmfoto vom 2021-07-22 11-42-37.png, Size: 120.68 KiB, Downloads: 8

Offline

#3 2021-07-22 10:54:24

mf
Member
Registered: 2019-06-18
Posts: 227

Re: SOLVED Simple Contact

One more thing: instead of using ASSE_MAILLAGE in Asterstudy (I also did this for a long time...) it is easier to Build a compound mesh in the mesh module. The result is more or less the same, but additionally you'll be able to click on your groups in Asterstudy instead of typing them in. With contact use 'Build Compound' like in the attached image (without 'Merge coincident nodes..', otherwise possibility of merging nodes in contact-->errors)
With build compound it is also easier to find errors in the mesh, e.g. if Asterstudy tells you N324543 is singular, you are easily able to find this node in the compound mesh in 'Mesh Information'.


Mario.

Last edited by mf (2021-07-22 11:00:26)


Attachments:
Bildschirmfoto vom 2021-07-22 11-50-30.png, Size: 251.03 KiB, Downloads: 17

Offline

#4 2021-07-23 13:45:45

rmoortgat
Member
From: Belgium
Registered: 2016-08-30
Posts: 37
Website

Re: SOLVED Simple Contact

Thank you Mario!

Roger

Offline