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

You are not logged in.

#1 2019-01-16 16:38:38

pironman
Member
Registered: 2018-05-15
Posts: 50

Simulate a mesh composed by more solids

Hi to everyone.
I have a component made by 3 separate parts, connected by bolts. I modelled it as a compound, then in the Mesh module of Salome-Meca I create a mesh composed by 3 sub-meshes of the three component.
Then, I load that mesh on the Aster module, apply constraints, loads and so on.
When I simulate it I have every time an error. In the first, the error was a lot of orphane nodes. After it's deletion, an error about "mailles doubles", and so on. Now I'm in a situation where I don't know how to correct the error.
Then, I try to simulate this system with only 1 body, instead of three, and all works correctly.

For this reason I'm thinking that my error is about the mesh, but I don't know how to correct it.

Any help?

Thank you very much,

Mattia

Offline

#2 2019-01-16 18:32:41

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

Re: Simulate a mesh composed by more solids

Hello,
This error is very common. You need to explode the compound into its constituent solid. Now make a mesh on compound without any hypothesis. Now create 3 sub meshes on the main mesh, each time cycling over the exploded solids. Algorithm is your wish, Netgen 1D-2D-3D works fine. Go ahead and compute the main mesh. Now create volume groups by selecting the top level mesh and using "create groups from geometry" icon. Use this mesh in Aster study.

Regards
Anirudh

Offline

#3 2019-01-17 11:26:21

pironman
Member
Registered: 2018-05-15
Posts: 50

Re: Simulate a mesh composed by more solids

Thankyou for your answer.
I repeat the simulation, but I have the same problem.

This is the .comm file:

DEBUT(LANG='EN')


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

model = AFFE_MODELE(identifier=u'1:1',
                    AFFE=(_F(GROUP_MA=('SUPP_2', 'SUPP_1', 'PIASTRA'),
                             MODELISATION=('3D', ),
                             PHENOMENE='MECANIQUE'),
                          _F(GROUP_MA=('TA1', 'TA2', 'TA3', 'IS3', 'IS2', 'IS1', 'TV1', 'TV3', 'TV2'),
                             MODELISATION=('DIS_TR', ),
                             PHENOMENE='MECANIQUE')),
                    MAILLAGE=mesh)

elemprop = AFFE_CARA_ELEM(identifier=u'2:1',
                          DISCRET=(_F(CARA='M_TR_D_N',
                                      GROUP_MA=('TA1', 'TA2', 'TA3'),
                                      VALE=(0.029, 191.5, 292.8, 209.4, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)),
                                   _F(CARA='M_TR_D_N',
                                      GROUP_MA=('TV1', 'TV3', 'TV2'),
                                      VALE=(0.02, 208.2, 150.9, 150.6, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)),
                                   _F(CARA='M_TR_D_N',
                                      GROUP_MA=('IS3', 'IS2', 'IS1'),
                                      VALE=(0.0011, 3.3, 3.3, 0.75, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0))),
                          MODELE=model)

steel = DEFI_MATERIAU(identifier=u'3:1',
                      ELAS=_F(AMOR_HYST=0.02,
                              E=190000.0,
                              NU=0.3,
                              RHO=7.85e-09))

fieldmat = AFFE_MATERIAU(identifier=u'4:1',
                         AFFE=_F(MATER=(steel, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh,
                         MODELE=model)

fix = AFFE_CHAR_MECA(identifier=u'5:1',
                     DDL_IMPO=_F(GROUP_NO=('FIX2', 'FIX1'),
                                 LIAISON='ENCASTRE'),
                     MODELE=model)

link = AFFE_CHAR_MECA(identifier=u'6:1',
                      LIAISON_SOLIDE=(_F(GROUP_NO=('TV1_BASE', 'TV1')),
                                      _F(GROUP_NO=('TV2_BASE', 'TV2')),
                                      _F(GROUP_NO=('TV3', 'TV3_BASE')),
                                      _F(GROUP_NO=('TA1', 'TA1_BASE')),
                                      _F(GROUP_NO=('TA2', 'TA2_BASE')),
                                      _F(GROUP_NO=('TA3_BASE', 'TA3')),
                                      _F(GROUP_NO=('IS1', 'IS1_BASE')),
                                      _F(GROUP_NO=('IS2', 'IS2_BASE')),
                                      _F(GROUP_NO=('IS3', 'IS3_BASE')),
                                      _F(GROUP_NO=('FORO_SUPP_1', 'FORO_VITE_1')),
                                      _F(GROUP_NO=('FORO_VITE_2', 'FORO_SUPP_2')),
                                      _F(GROUP_NO=('FORO_VITE_3', 'FORO_SUPP_3')),
                                      _F(GROUP_NO=('FORO_SUPP_4', 'FORO_VITE_4')),
                                      _F(GROUP_NO=('FORO_SUPP_5', 'FORO_VITE_5')),
                                      _F(GROUP_NO=('FORO_PIASTRA_6', 'FORO_VITE_6')),
                                      _F(GROUP_NO=('FORO_VITE_7', 'FORO_PIASTRA_7')),
                                      _F(GROUP_NO=('FORO_PIASTRA_8', 'FORO_VITE_8')),
                                      _F(GROUP_NO=('FORO_VITE_9', 'FORO_PIASTRA_9')),
                                      _F(GROUP_NO=('FORO_PIASTRA_10', 'FORO_VITE_'))),
                      MODELE=model)

ASSEMBLAGE(identifier=u'7:1',
           CARA_ELEM=elemprop,
           CHAM_MATER=fieldmat,
           CHARGE=(fix, link),
           MATR_ASSE=(_F(MATRICE=CO('M'),
                         OPTION='MASS_MECA'),
                      _F(MATRICE=CO('K'),
                         OPTION='RIGI_MECA')),
           MODELE=model,
           NUME_DDL=CO('nddl'))

modes = CALC_MODES(identifier=u'10:1',
                   CALC_FREQ=_F(NMAX_FREQ=10),
                   MATR_MASS=M,
                   MATR_RIGI=K,
                   OPTION='PLUS_PETITE',
                   TYPE_RESU='DYNAMIQUE',
                   VERI_MODE=_F(STOP_ERREUR='NON'))

IMPR_RESU(identifier=u'11:1',
          FORMAT='MED',
          RESU=_F(RESULTAT=modes),
          UNITE=80)

FIN()

And this is the message file:

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


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

model = AFFE_MODELE(identifier=u'1:1',
                    AFFE=(_F(GROUP_MA=('SUPP_2', 'SUPP_1', 'PIASTRA'),
                             MODELISATION=('3D', ),
                             PHENOMENE='MECANIQUE'),
                          _F(GROUP_MA=('TA1', 'TA2', 'TA3', 'IS3', 'IS2', 'IS1', 'TV1', 'TV3', 'TV2'),
                             MODELISATION=('DIS_TR', ),
                             PHENOMENE='MECANIQUE')),
                    MAILLAGE=mesh)

elemprop = AFFE_CARA_ELEM(identifier=u'2:1',
                          DISCRET=(_F(CARA='M_TR_D_N',
                                      GROUP_MA=('TA1', 'TA2', 'TA3'),
                                      VALE=(0.029, 191.5, 292.8, 209.4, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)),
                                   _F(CARA='M_TR_D_N',
                                      GROUP_MA=('TV1', 'TV3', 'TV2'),
                                      VALE=(0.02, 208.2, 150.9, 150.6, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0)),
                                   _F(CARA='M_TR_D_N',
                                      GROUP_MA=('IS3', 'IS2', 'IS1'),
                                      VALE=(0.0011, 3.3, 3.3, 0.75, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0))),
                          MODELE=model)

steel = DEFI_MATERIAU(identifier=u'3:1',
                      ELAS=_F(AMOR_HYST=0.02,
                              E=190000.0,
                              NU=0.3,
                              RHO=7.85e-09))

fieldmat = AFFE_MATERIAU(identifier=u'4:1',
                         AFFE=_F(MATER=(steel, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh,
                         MODELE=model)

fix = AFFE_CHAR_MECA(identifier=u'5:1',
                     DDL_IMPO=_F(GROUP_NO=('FIX2', 'FIX1'),
                                 LIAISON='ENCASTRE'),
                     MODELE=model)

link = AFFE_CHAR_MECA(identifier=u'6:1',
                      LIAISON_SOLIDE=(_F(GROUP_NO=('TV1_BASE', 'TV1')),
                                      _F(GROUP_NO=('TV2_BASE', 'TV2')),
                                      _F(GROUP_NO=('TV3', 'TV3_BASE')),
                                      _F(GROUP_NO=('TA1', 'TA1_BASE')),
                                      _F(GROUP_NO=('TA2', 'TA2_BASE')),
                                      _F(GROUP_NO=('TA3_BASE', 'TA3')),
                                      _F(GROUP_NO=('IS1', 'IS1_BASE')),
                                      _F(GROUP_NO=('IS2', 'IS2_BASE')),
                                      _F(GROUP_NO=('IS3', 'IS3_BASE')),
                                      _F(GROUP_NO=('FORO_SUPP_1', 'FORO_VITE_1')),
                                      _F(GROUP_NO=('FORO_VITE_2', 'FORO_SUPP_2')),
                                      _F(GROUP_NO=('FORO_VITE_3', 'FORO_SUPP_3')),
                                      _F(GROUP_NO=('FORO_SUPP_4', 'FORO_VITE_4')),
                                      _F(GROUP_NO=('FORO_SUPP_5', 'FORO_VITE_5')),
                                      _F(GROUP_NO=('FORO_PIASTRA_6', 'FORO_VITE_6')),
                                      _F(GROUP_NO=('FORO_VITE_7', 'FORO_PIASTRA_7')),
                                      _F(GROUP_NO=('FORO_PIASTRA_8', 'FORO_VITE_8')),
                                      _F(GROUP_NO=('FORO_VITE_9', 'FORO_PIASTRA_9')),
                                      _F(GROUP_NO=('FORO_PIASTRA_10', 'FORO_VITE_'))),
                      MODELE=model)

ASSEMBLAGE(identifier=u'7:1',
           CARA_ELEM=elemprop,
           CHAM_MATER=fieldmat,
           CHARGE=(fix, link),
           MATR_ASSE=(_F(MATRICE=CO('M'),
                         OPTION='MASS_MECA'),
                      _F(MATRICE=CO('K'),
                         OPTION='RIGI_MECA')),
           MODELE=model,
           NUME_DDL=CO('nddl'))

modes = CALC_MODES(identifier=u'10:1',
                   CALC_FREQ=_F(NMAX_FREQ=10),
                   MATR_MASS=M,
                   MATR_RIGI=K,
                   OPTION='PLUS_PETITE',
                   TYPE_RESU='DYNAMIQUE',
                   VERI_MODE=_F(STOP_ERREUR='NON'))

IMPR_RESU(identifier=u'11:1',
          FORMAT='MED',
          RESU=_F(RESULTAT=modes),
          UNITE=80)

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

  Ouverture en écriture du fichier ./vola.1

 <INFO> Démarrage de l'exécution.
                                                                                
                                                                                
              -- CODE_ASTER -- VERSION:  EXPLOITATION (stable) --               
                                                                                
                       Modified 13.6.0 version 07/08/2018                       
                    the revision 3c51e6c237ea - branch  v13                     
                        Copyright EDF R & D 1991 - 2019                         
                                                                                
                    Exécution of:  Thu Jan 17 11:15:27 2019                     
                          Name of the machine:  debian                          
                               Structure:  64bit                                
                               CPU type:  x86_64                                
               Operating system:  Linux debian 9.6  4.9.0-8-amd64               
                     Language of the messages:  en (UTF-8)                      
                                                                                
                             Python version: 2.7.10                             
                             NumPy version: 1.12.1                              
                             Parallelism MPI: idle                              
                           OpenMP parallelism: Nombre                           
                          credit of processes used:  2                          
                      Version of the HDF5 library : 1.8.14                      
                       Version of the MED library : 3.3.1                       
                     Version de la librairie MFront : 3.0.0                     
                          MUMPS library version: 5.1.1                          
                          PETSc library: not available                          
                     Version of the SCOTCH library : 6.0.4                      

                  Mémoire limite pour l'exécution : 4096.00 Mo                  
                   consommée par l'initialisation : 546.44 Mo                   
                  par les objets du jeu de commandes : 9.83 Mo                  
                 reste pour l'allocation dynamique : 3539.60 Mo                 
                   Cut limiting files of exchange:  12.00 Go                    
                                                                                

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


 Initial value of maximum time CPU = 900 second
 Valeur of the maximum time CPU placed to the orders = 810 second
 Réserve CPU envisaged = 90 seconds 

  Ouverture en écriture du fichier ./glob.1

  Ouverture en écriture du fichier ./vola.1

  Ouverture en lecture du fichier /home/mattia/salome_meca/V2018.0.1_public/tools/Code_aster_stable-v136_smeca/lib/aster/elem.1

     Nom de la base                          :  ELEMBASE
     Créée avec la version                   :  13.06.00
     Nombre d'enregistrements utilisés       :  44
     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/mattia/salome_meca/V2018.0.1_public/tools/Code_aster_stable-v136_smeca/lib/aster/elem.1


     Nom de la base                          :  ELEMBASE
     Nombre d'enregistrements utilisés       :  44
     Nombre d'enregistrements maximum        :  512
     Longueur d'enregistrement (octets)      :  819200
     Nombre total d'accès en lecture         :  62
     Volume des accès en lecture             :         48.44 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 %

 Second reading of the catalogues of the elements made. 

 End of reading (lasted 0.038937 S.)   
  # Memory (Mo):    601.64 /    593.80 /     44.67 /     23.16 (VmPeak / VmSize / Optimum / Minimum)
  # Fin orders No:  0001 user+syst:          0.05S (syst:          0.04S, Elaps:          0.07S)
  # ------------------------------------------------------------------------------------------

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

Vérification du maillage.
   
   !-------------------------------------------------------------------!
   ! <A> <MODELISA8_13>                                                !
   !                                                                   !
   ! Les mailles M6687 et M6618 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 M6886 et M6663 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 M6984 et M6633 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 M6998 et M6976 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 M7024 et M6902 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!                                    !
   !---------------------------------------------------------------------------------!
   

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

JE-17-JANV-2019 11:15:27                                                        

NOMBRE DE NOEUDS                      234875

NOMBRE DE MAILLES                     294685
                              POI1                    9
                              SEG3                17393
                              TRIA6              114210
                              TETRA4             163073

NOMBRE DE GROUPES DE NOEUDS               40
                              TV1                                  1
                              TV3                                  1
                              TA1                                  1
                              TA2                                  1
                              TA3                                  1
                              IS3                                  1
                              IS2                                  1
                              IS1                                  1
                              TV2                                  1
                              FORO_SUPP_1                         66
                              FORO_SUPP_2                         66
                              FORO_SUPP_3                         66
                              FORO_SUPP_4                         66
                              FORO_SUPP_5                         66
                              FIX1                               328
                              FORO_VITE_1                         48
                              FORO_VITE_2                         48
                              FORO_VITE_3                         48
                              FORO_VITE_4                         48
                              FORO_VITE_5                         48
                              FORO_VITE_6                         48
                              FORO_VITE_7                         48
                              FORO_VITE_8                         48
                              FORO_VITE_9                         48
                              FORO_VITE_                          48
                              TV3_BASE                           192
                              TV2_BASE                           192
                              TV1_BASE                           192
                              TA1_BASE                           192
                              TA2_BASE                           192
                              TA3_BASE                           192
                              IS1_BASE                            72
                              IS2_BASE                            72
                              IS3_BASE                            72
                              FIX2                               328
                              FORO_PIASTRA_6                      66
                              FORO_PIASTRA_7                      66
                              FORO_PIASTRA_8                      66
                              FORO_PIASTRA_9                      66
                              FORO_PIASTRA_10                     66

NOMBRE DE GROUPES DE MAILLES              43
                              PIASTRA                          87828
                              SUPP_1                           37024
                              SUPP_2                           38221
                              TV3                                  1
                              TV2                                  1
                              TV1                                  1
                              IS3                                  1
                              IS2                                  1
                              IS1                                  1
                              TA1                                  1
                              TA2                                  1
                              TA3                                  1
                              FORO_SUPP_1                         16
                              FORO_SUPP_2                         16
                              FORO_SUPP_3                         16
                              FORO_SUPP_4                         16
                              FORO_SUPP_5                         16
                              FIX1                                80
                              FORO_VITE_1                         16
                              FORO_VITE_2                         16
                              FORO_VITE_3                         16
                              FORO_VITE_4                         16
                              FORO_VITE_5                         16
                              FORO_VITE_6                         16
                              FORO_VITE_7                         16
                              FORO_VITE_8                         16
                              FORO_VITE_9                         16
                              FORO_VITE_                          16
                              TV3_BASE                            64
                              TV2_BASE                            64
                              TV1_BASE                            64
                              TA1_BASE                            64
                              TA2_BASE                            64
                              TA3_BASE                            64
                              IS1_BASE                            24
                              IS2_BASE                            24
                              IS3_BASE                            24
                              FIX2                                80
                              FORO_PIASTRA_6                      16
                              FORO_PIASTRA_7                      16
                              FORO_PIASTRA_8                      16
                              FORO_PIASTRA_9                      16
                              FORO_PIASTRA_10                     16

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

  # Memory (Mo):    655.41 /    631.80 /     97.11 /     58.74 (VmPeak / VmSize / Optimum / Minimum)
  # Fin orders No:  0002 user+syst:          1.28S (syst:          0.07S, Elaps:          1.28S)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd1:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0003 Concept of the type:  modele_sdaster
  # ------------------------------------------------------------------------------------------
  model = AFFE_MODELE(MAILLAGE=mesh,
                      AFFE=(_F(GROUP_MA=('SUPP_2', 'SUPP_1', 'PIASTRA'),
                               PHENOMENE='MECANIQUE',
                               MODELISATION=('3D', ),),
                            _F(GROUP_MA=('TA1', 'TA2', 'TA3', 'IS3', 'IS2', 'IS1', 'TV1', 'TV3', 'TV2'),
                               PHENOMENE='MECANIQUE',
                               MODELISATION=('DIS_TR', ),),
                            ),
                      INFO=1,
                      VERI_JACOBIEN='OUI',
                      DISTRIBUTION=_F(PARTITIONNEUR='METIS',
                                      METHODE='SOUS_DOMAINE',),
                      )


Sur les 294685 mailles du maillage mesh, on a demandé l'affectation de 163082, on a pu en affecter 163082
 Modélisation     Type maille  Élément fini     Nombre
 DIS_TR           POI1         MECA_DIS_TR_N    9
 3D               TETRA4       MECA_TETRA4      163073
#2      Calculs elementaires et assemblages               CPU (USER+SYST/SYST/ELAPS):      0.10      0.01      0.09
  # Memory (Mo):    688.44 /    682.11 /     97.11 /     68.85 (VmPeak / VmSize / Optimum / Minimum)
  # Fin orders No:  0003 user+syst:          0.27S (syst:          0.08S, Elaps:          0.19S)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd2:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0004 Concept of the type:  cara_elem
  # ------------------------------------------------------------------------------------------
  elemprop = AFFE_CARA_ELEM(MODELE=model,
                            DISCRET=(_F(CARA='M_TR_D_N',
                                        GROUP_MA=('TA1', 'TA2', 'TA3'),
                                        VALE=(0.029, 191.5, 292.8, 209.4, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0),
                                        REPERE='GLOBAL',
                                        SYME='OUI',),
                                     _F(CARA='M_TR_D_N',
                                        GROUP_MA=('TV1', 'TV3', 'TV2'),
                                        VALE=(0.02, 208.2, 150.9, 150.6, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0),
                                        REPERE='GLOBAL',
                                        SYME='OUI',),
                                     _F(CARA='M_TR_D_N',
                                        GROUP_MA=('IS3', 'IS2', 'IS1'),
                                        VALE=(1.1E-3, 3.3, 3.3, 0.75, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0),
                                        REPERE='GLOBAL',
                                        SYME='OUI',),
                                     ),
                            INFO=1,)


     LE MODELE model    CONTIENT : 
                                        9 ELEMENT(S) DISCRET         
                                   163073 ELEMENT(S) MASSIF          
#2      Calculs elementaires et assemblages               CPU (USER+SYST/SYST/ELAPS):      0.00      0.00      0.00
  # Memory (Mo):    784.41 /    750.81 /    192.22 /    118.72 (VmPeak / VmSize / Optimum / Minimum)
  # Fin orders No:  0004 user+syst:          2.58S (syst:          0.18S, Elaps:          2.59S)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd3:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0006 Concept of the type:  mater_sdaster
  # ------------------------------------------------------------------------------------------
  steel = DEFI_MATERIAU(ELAS=_F(E=1.90000E5,
                                AMOR_HYST=0.02,
                                RHO=7.85E-09,
                                NU=0.3,
                                COEF_AMOR=1.0,),
                        INFO=1,)


 MATERIAU : steel   
 RELATION DE COMPORTEMENT: ELAS                            


  # Memory (Mo):    784.41 /    750.81 /    192.22 /    118.72 (VmPeak / VmSize / Optimum / Minimum)
  # Fin orders No:  0006 user+syst:          0.00S (syst:          0.00S, Elaps:          0.01S)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd4:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0007 Concept of the type:  cham_mater
  # ------------------------------------------------------------------------------------------
  fieldmat = AFFE_MATERIAU(MAILLAGE=mesh,
                           MODELE=model,
                           AFFE=_F(MATER=(steel, ),
                                   TOUT='OUI',),
                           INFO=1,
                           )

  # Memory (Mo):    784.41 /    750.81 /    192.22 /    118.72 (VmPeak / VmSize / Optimum / Minimum)
  # Fin orders No:  0007 user+syst:          0.01S (syst:          0.00S, Elaps:          0.00S)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd5:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0008 Concept of the type:  char_meca
  # ------------------------------------------------------------------------------------------
  fix = AFFE_CHAR_MECA(MODELE=model,
                       DDL_IMPO=_F(LIAISON='ENCASTRE',
                                   GROUP_NO=('FIX2', 'FIX1'),),
                       INFO=1,
                       VERI_NORM='OUI',)

  # Memory (Mo):    784.41 /    759.93 /    192.22 /    118.72 (VmPeak / VmSize / Optimum / Minimum)
  # Fin orders No:  0008 user+syst:          0.39S (syst:          0.02S, Elaps:          0.40S)
  # ------------------------------------------------------------------------------------------

.. __stg1_cmd6:1
  # ------------------------------------------------------------------------------------------
  # Commands No:  0009 Concept of the type:  char_meca
  # ------------------------------------------------------------------------------------------
  link = AFFE_CHAR_MECA(MODELE=model,
                        LIAISON_SOLIDE=(_F(GROUP_NO=('TV1_BASE', 'TV1'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('TV2_BASE', 'TV2'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('TV3', 'TV3_BASE'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('TA1', 'TA1_BASE'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('TA2', 'TA2_BASE'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('TA3_BASE', 'TA3'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('IS1', 'IS1_BASE'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('IS2', 'IS2_BASE'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('IS3', 'IS3_BASE'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_SUPP_1', 'FORO_VITE_1'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_VITE_2', 'FORO_SUPP_2'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_VITE_3', 'FORO_SUPP_3'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_SUPP_4', 'FORO_VITE_4'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_SUPP_5', 'FORO_VITE_5'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_PIASTRA_6', 'FORO_VITE_6'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_VITE_7', 'FORO_PIASTRA_7'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_PIASTRA_8', 'FORO_VITE_8'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_VITE_9', 'FORO_PIASTRA_9'),
                                           NUME_LAGR='NORMAL',),
                                        _F(GROUP_NO=('FORO_PIASTRA_10', 'FORO_VITE_'),
                                           NUME_LAGR='NORMAL',),
                                        ),
                        INFO=1,
                        VERI_NORM='OUI',)

   
   !-----------------------------------------------------------------------------------!
   ! <EXCEPTION> <CHARGES2_31>                                                         !
   !                                                                                   !
   !  Erreur utilisateur:                                                              !
   !     On cherche à imposer une condition aux limites sur le ddl DX du noeud N57001. !
   !     Mais ce noeud ne porte pas ce ddl.                                            !
   !                                                                                   !
   !     Conseils :                                                                    !
   !      - vérifiez le modèle et les conditions aux limites :                         !
   !         - le noeud incriminé fait-il partie du modèle ?                           !
   !         - le noeud porte-t-il le ddl que l'on cherche à contraindre ?             !
   !-----------------------------------------------------------------------------------!
   

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

.. __stg1_txt91
  # ------------------------------------------------------------------------------------------
  # Commands No:  0010 Concept of the type:  -
  # ------------------------------------------------------------------------------------------
  FIN(RETASSAGE='NON',
      INFO_RESU='OUI',
      FORMAT_HDF='NON',
      PROC0='OUI',
      )

   
   !-----------------------------------------------------------------!
   ! <I> <CATAMESS_89>                                               !
   !                                                                 !
   !  List of warnings emitted during the execution of computation.  !
   !                                                                 !
   !  Warnings which you chose to ignore of are preceded by (*).     !
   !  Number of occurrences for each warning:                        !
   !              MODELISA4_8          emitted    1 time             !
   !              MODELISA8_13         emitted 5331 time             !
   !-----------------------------------------------------------------!
   
 ---------------------------------------------------------------------------
 Concepts de la base: G
    Nom       Type                 Taille (Mo)         Nombre      Nombre de
                                                     d'objets       segments
    TOTAL                               158.41            150            370
  
    mesh       MAILLAGE_SDASTER          36.67             36            140
    model      MODELE_SDASTER            16.95              8             13
    elemprop   CARA_ELEM                 68.02             27             52
    steel      MATER_SDASTER              0.00              5              5
    fieldmat   CHAM_MATER                 0.00              6             11
    fix        CHAR_MECA                  9.14             23             43
    &FOZERO                               0.00              2              2
    &&_NUM_C                              0.00              1              1
    &CATA.AC                              0.00              2              4
    &CATA.CL                              0.86              1              3
    &CATA.GD                              0.14              4             11
    &CATA.ME                              0.20              2              4
    &CATA.OP                              0.32              4             19
    &CATA.PH                              0.00              1              1
    &CATA.PR                              0.00              2              4
    &CATA.TE                             26.07             17             42
    &CATA.TH                              0.01              2              4
    &CATA.TM                              0.01              7             11
 ---------------------------------------------------------------------------


     Nom de la base                          :  GLOBALE
     Nombre d'enregistrements utilisés       :  221
     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        :  235
     Volume des accès en écriture            :        183.59 Mo.
     Nombre d'identificateurs utilisés       :  346
     Taille maximum du répertoire            :  2000
     Pourcentage d'utilisation du répertoire :  17 %


     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        :  4
     Volume des accès en écriture            :          3.12 Mo.
     Nombre d'identificateurs utilisés       :  236
     Taille maximum du répertoire            :  2000
     Pourcentage d'utilisation du répertoire :  11 %

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

 The memory requested from launching is over-estimated, it is of 4096.00 Mo. 
 The peak report used is of 784.41 Mo. 
 


     Nom de la base                          :  ELEMBASE
     Nombre d'enregistrements utilisés       :  44
     Nombre d'enregistrements maximum        :  512
     Longueur d'enregistrement (octets)      :  819200
     Nombre total d'accès en lecture         :  62
     Volume des accès en lecture             :         48.44 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 :                     118.72 Mo
 <I> <FIN> MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION :                     192.22 Mo
 <I> <FIN> MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS LORS DE L'EXECUTION :     784.41 Mo
  
 <I>       FERMETURE DES BASES EFFECTUEE
  
   STATISTIQUES CONCERNANT L'ALLOCATION DYNAMIQUE :
     TAILLE CUMULEE MAXIMUM            :                  192  Mo.
     TAILLE CUMULEE LIBEREE            :                  188  Mo.
     NOMBRE TOTAL D'ALLOCATIONS        :              1236411
     NOMBRE TOTAL DE LIBERATIONS       :              1236391
     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 :     118.72 Mo
     - IMPOSE DE NOMBREUX ACCES DISQUE
     - RALENTIT LA VITESSE D'EXECUTION
   MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION :     192.22 Mo
     - LIMITE LES ACCES DISQUE
     - AMELIORE LA VITESSE D'EXECUTION
   MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS     :     784.41 Mo
     - COMPREND LA MEMOIRE CONSOMMEE PAR  JEVEUX, 
       LE SUPERVISEUR PYTHON, LES LIBRAIRIES EXTERNES
  
 <I>       FIN D'EXECUTION LE : JE-17-JANV-2019 11:15:32

 ********************************************************************************
 * COMMAND                  :       USER :     SYSTEM :   USER+SYS :    ELAPSED *
 ********************************************************************************
 * init (jdc)               :       0.54 :       0.06 :       0.60 :       0.61 *
 *  . compile               :       0.00 :       0.00 :       0.00 :       0.01 *
 *  . exec_compile          :       0.11 :       0.01 :       0.12 :       0.11 *
 *  . report                :       0.01 :       0.00 :       0.01 :       0.01 *
 *  . build                 :       0.00 :       0.00 :       0.00 :       0.00 *
 * DEBUT                    :       0.01 :       0.04 :       0.05 :       0.07 *
 * LIRE_MAILLAGE            :       1.21 :       0.07 :       1.28 :       1.28 *
 * AFFE_MODELE              :       0.19 :       0.08 :       0.27 :       0.19 *
 * AFFE_CARA_ELEM           :       2.40 :       0.18 :       2.58 :       2.59 *
 * DEFI_MATERIAU            :       0.00 :       0.00 :       0.00 :       0.01 *
 * AFFE_MATERIAU            :       0.01 :       0.00 :       0.01 :       0.00 *
 * AFFE_CHAR_MECA           :       0.37 :       0.02 :       0.39 :       0.40 *
 * AFFE_CHAR_MECA           :       0.07 :       0.15 :       0.22 :       0.39 *
 * FIN                      :       0.03 :       0.11 :       0.14 :       0.17 *
 *  . part Superviseur      :       0.58 :       0.10 :       0.68 :       0.73 *
 *  . part Fortran          :       4.23 :       0.50 :       4.73 :       4.83 *
 ********************************************************************************
 * TOTAL_JOB                :       4.81 :       0.60 :       5.41 :       5.56 *
 ********************************************************************************



  <I> Informations sur les temps d'exécution
      Temps cpu user total                    4.81 s
      Temps cpu systeme total                 0.60 s
      Temps elapsed total                     5.56 s
      Temps restant                         804.44 s


 The signup of the saved base is (with the address 0): 
  28025dd0c4e3bacb86c00b8e7d34065cde707a6b 

 ********************************************************************************
 * COMMAND                  :       USER :     SYSTEM :   USER+SYS :    ELAPSED *
 ********************************************************************************
 * init (jdc)               :       0.54 :       0.06 :       0.60 :       0.61 *
 *  . compile               :       0.00 :       0.00 :       0.00 :       0.01 *
 *  . exec_compile          :       0.11 :       0.01 :       0.12 :       0.11 *
 *  . report                :       0.01 :       0.00 :       0.01 :       0.01 *
 *  . build                 :       0.00 :       0.00 :       0.00 :       0.00 *
 * DEBUT                    :       0.01 :       0.04 :       0.05 :       0.07 *
 * LIRE_MAILLAGE            :       1.21 :       0.07 :       1.28 :       1.28 *
 * AFFE_MODELE              :       0.19 :       0.08 :       0.27 :       0.19 *
 * AFFE_CARA_ELEM           :       2.40 :       0.18 :       2.58 :       2.59 *
 * DEFI_MATERIAU            :       0.00 :       0.00 :       0.00 :       0.01 *
 * AFFE_MATERIAU            :       0.01 :       0.00 :       0.01 :       0.00 *
 * AFFE_CHAR_MECA           :       0.37 :       0.02 :       0.39 :       0.40 *
 * AFFE_CHAR_MECA           :       0.07 :       0.15 :       0.22 :       0.39 *
 * FIN                      :       0.03 :       0.11 :       0.14 :       0.17 *
 *  . part Superviseur      :       0.58 :       0.10 :       0.68 :       0.73 *
 *  . part Fortran          :       4.23 :       0.50 :       4.73 :       4.83 *
 ********************************************************************************
 * TOTAL_JOB                :       4.81 :       0.60 :       5.41 :       5.56 *
 ********************************************************************************

JDC.py : ERREUR WITH THE EXECUTION - INTERRUPTION
>> JDC.py: DEBUT CR of execution of JDC in MIXTE
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception user raised but not interceptee.                            !
   ! The bases are fermees.                                                    !
   ! Type of the exception: error                                              !
   !                                                                           !
   !  Erreur utilisateur:                                                      !
   !     On cherche à imposer une condition aux limites sur le ddl DX du noeud !
   ! N57001.                                                                   !
   !     Mais ce noeud ne porte pas ce ddl.                                    !
   !                                                                           !
   !     Conseils :                                                            !
   !      - vérifiez le modèle et les conditions aux limites :                 !
   !         - le noeud incriminé fait-il partie du modèle ?                   !
   !         - le noeud porte-t-il le ddl que l'on cherche à contraindre ?     !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fine CR of execution of JDC in MIXTE

 REPORT
 >> JDC.py: FIN REPORT
EXECUTION_CODE_ASTER_EXIT_1834-0002-debian=1

 <I>_EXIT_CODE = 1

I'm stopped here. The last time, in the mesh module, I use the command "modification -> transformation -> merge elements" in order to remove the warning "Les mailles M**** et M**** sont doubles", but the error "Erreur utilisateur: On cherche à imposer une condition aux limites sur le ddl DX du noeud N*****. Mais ce noeud ne porte pas ce ddl" remain, then I don't understand what is wrong.


P.S. I think that the problem is in the "link" BC and load, I try to switch from "GROUP_NO" to "GROUP_MA", but nothing change

Last edited by pironman (2019-01-17 12:18:07)

Offline

#4 2019-01-17 14:57:43

jeanpierreaubry
Guru
From: nantes (france)
Registered: 2009-03-12
Posts: 3,569

Re: Simulate a mesh composed by more solids

hello

the .comm file is just not enough
without the mesh it is impossible to identify what is wrong

jean pierre aubry

Offline

#5 2019-01-17 15:08:56

pironman
Member
Registered: 2018-05-15
Posts: 50

Re: Simulate a mesh composed by more solids

Hi, here there is an archive with the .hdf file.

https://drive dot google dot com/file/d/1dpCnry2a76gR4aeaI8FzkGr8xNqDUzGM/view?usp=sharing

Thank you

P.S. Another one, in which I have modified the group of nodes, but the results are the same (doesn't work)

https://drive dot google dot com/file/d/1G2LxzGznt7_RbhGSYgz96Rtvy3W-k1ve/view?usp=sharing

Last edited by pironman (2019-01-17 15:23:38)

Offline

#6 2019-01-22 10:49:21

pironman
Member
Registered: 2018-05-15
Posts: 50

Re: Simulate a mesh composed by more solids

Right now, I solved by adding some 3D-parts that connect the two sub-meshes, and works. In this case is ok, but in some case I need to connect two submeshes, and I still don't know how to do that.

Any help?

Thank you

Offline

#7 2019-01-22 16:58:32

nirmaljoshi
Member
From: Japan
Registered: 2018-10-12
Posts: 146

Re: Simulate a mesh composed by more solids

Use Operation>Partition method after forming the compound
(referece- docs.salome-platform.org/latest/gui/GEOM/partition_explanation.html)

Offline

#8 2019-01-25 10:06:24

pironman
Member
Registered: 2018-05-15
Posts: 50

Re: Simulate a mesh composed by more solids

Thank you, but I want to connect two solid not in contact between them, then I can't use partition.

Offline

#9 2019-01-28 10:05:28

pironman
Member
Registered: 2018-05-15
Posts: 50

Re: Simulate a mesh composed by more solids

Now I'm trying with LIAISON_MAIL, and I obtain this error message:

   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee.                       !
   ! Les bases sont fermees.                                                     !
   ! Type de l'exception : error                                                 !
   !                                                                             !
   !  les 12 mailles imprimées ci-dessus n'appartiennent pas au modèle           !
   !  et pourtant elles ont été affectées dans le mot-clé facteur : LIAISON_MAIL !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Offline