# ------------------------------------------------------------------------------ Execution of code_aster # ------------------------------------------------------------------------------ Prepare environment in /tmp/run_aster_8_qv4h3g/proc.0 # ------------------------------------------------------------------------------ Command file #1 / 1 Content of the file to execute: # coding=utf-8 # original filename __file__ = r"/out/M02_oring.comm" # EXPORT CONTROLLED — ITAR USML Category XX (Ocean Guard M02) # Code_Aster STAT_NON_LINE hyperelastic O-ring/groove squeeze — AS568A-333 EPDM 90A. # Auto-generated by oring_hyperelastic_comm.py. Do not hand-edit. # temporarly added for compatibility with code_aster legacy from math import * import code_aster from code_aster.Commands import * from code_aster import CA DEBUT(LANG='EN') MESH = LIRE_MAILLAGE(FORMAT='MED', UNITE=20) # D_PLAN (plane strain), not AXIS — see module docstring 2026-07-09 note # (ELASHYPER_97: this Aster build forbids ELAS_HYPER+AXIS outright). MODEL = AFFE_MODELE( MAILLAGE=MESH, AFFE=_F(TOUT='OUI', PHENOMENE='MECANIQUE', MODELISATION='D_PLAN'), ) # EPDM 90A hyperelastic Mooney-Rivlin (dual-oring-sealing skill FEM baseline). MAT = DEFI_MATERIAU( ELAS_HYPER=_F(C10=5.500000e+05, C01=1.400000e+05, C20=0.0, K=2.000000e+09), ) FIELDMAT = AFFE_MATERIAU(MAILLAGE=MESH, AFFE=_F(TOUT='OUI', MATER=MAT)) # Groove floor fixed. Gland face ramps to the full squeeze displacement # (0.9800 mm = 21.9% squeeze on CS=4.47mm). # GROOVE_FLOOR alone (full DX=DY=0 over a non-degenerate 2D edge) already # removes all 3 planar rigid-body modes of this D_PLAN model — see # generate_oring_comm_file() docstring, 2026-07-10 CONTACT-BC fix note. LOAD_FIXED = AFFE_CHAR_MECA( MODELE=MODEL, DDL_IMPO=_F(GROUP_MA='GROOVE_FLOOR', DX=0.0, DY=0.0), ) RAMP = DEFI_FONCTION( NOM_PARA='INST', VALE=(0.0, 0.0, 1.0, 1.0), PROL_DROITE='CONSTANT', PROL_GAUCHE='CONSTANT', ) LOAD_SQUEEZE = AFFE_CHAR_MECA( MODELE=MODEL, DDL_IMPO=_F(GROUP_MA='GLAND_FACE', DY=9.800000e-04), ) LINST = DEFI_LIST_REEL(DEBUT=0.0, INTERVALLE=_F(JUSQU_A=1.0, NOMBRE=40)) DLIST = DEFI_LIST_INST( METHODE='AUTO', DEFI_LIST=_F(LIST_INST=LINST), ECHEC=_F(ACTION='DECOUPE', SUBD_METHODE='MANUEL', SUBD_PAS=4, SUBD_NIVEAU=6), ) # STAT_NON_LINE — large displacement + hyperelastic (ELAS_HYPER). # NOTE (2026-07-08 fix, COMPOR1_44): Code_Aster 17.4 rejects # DEFORMATION='GROT_GDEP' combined with RELATION='ELAS_HYPER' # (" Le type de deformation choisi est # incompatible avec le comportement ") — box-confirmed real # run, first real STAT_NON_LINE attempt on this .comm. Same class of fix # already applied in gate_codeaster_nonlinear.py's generate_comm_file() # (that file's own 2026-07-03 COMPOR5_22 note): Aster's own exception text # recommends GREEN_LAGRANGE as the large-strain-compatible alternative for # hyperelastic materials on this mesh/element combination. # NOTE (2026-07-10 fix, MECANONLINE9_4): CONTACT= removed — see module # docstring "CONTACT-BC FIX" section. GROOVE_FLOOR + GLAND_FACE's own # independent Dirichlet BCs fully determine this model without it. RESU = STAT_NON_LINE( MODELE=MODEL, CHAM_MATER=FIELDMAT, EXCIT=( _F(CHARGE=LOAD_FIXED), _F(CHARGE=LOAD_SQUEEZE, FONC_MULT=RAMP), ), COMPORTEMENT=_F(RELATION='ELAS_HYPER', DEFORMATION='GREEN_LAGRANGE', TOUT='OUI'), INCREMENT=_F(LIST_INST=DLIST), NEWTON=_F(MATRICE='TANGENTE', REAC_ITER=1), CONVERGENCE=_F(RESI_GLOB_RELA=1.0e-5, ITER_GLOB_MAXI=30), ) RESU = CALC_CHAMP(reuse=RESU, RESULTAT=RESU, CONTRAINTE=('SIEF_NOEU',), FORCE=('REAC_NODA',)) IMPR_RESU(FORMAT='RESULTAT', RESU=_F(RESULTAT=RESU, NOM_CHAM='DEPL')) # Reaction-force contact-pressure proxy (2026-07-10 fix — replaces the # redundant DEFI_CONTACT/LAGR_C extraction; see module docstring # "CONTACT-PRESSURE EXTRACTION — REVISED"). GLAND_FACE already carries a # complete, independent DDL_IMPO(DY=squeeze) Dirichlet BC — STAT_NON_LINE # computes REAC_NODA (the reaction force needed to hold that prescribed # displacement) natively at every constrained node, no CONTACT block # required. In a compressed-slab idealization (this D_PLAN cross-section: # fixed GROOVE_FLOOR / prescribed-displaced GLAND_FACE), the reaction force # at the prescribed face divided by its contact width IS the interface # contact pressure. GLAND_FACE's real in-plane width equals the O-ring # cross-section CS_MM (4.47mm — bbox-confirmed in # KD-CODEASTER-ORING-CONTACT-SINGULAR-2026-07-09.md: GLAND_FACE spans # x=[147.765,152.235]mm = 4.47mm = CS_MM exactly). D_PLAN reports REAC_NODA # as a force-per-unit-out-of-plane-depth (Aster's plane-strain convention # treats the model as 1m thick); dividing the summed DY reaction by the # real in-plane width converts N/m -> Pa, i.e. contact pressure. TAB_REAC = POST_RELEVE_T( ACTION=_F(INTITULE='GLAND_REACTION', GROUP_NO='GLAND_FACE', RESULTAT=RESU, NOM_CHAM='REAC_NODA', NOM_CMP='DY', OPERATION='SOMME'), ) # Deterministic Python-level marker (unchanged format from Round-1 — see # module docstring) so downstream parsing needs no changes. _gland_width_m = 4.470000e-03 _max_contact_pa = 0.0 try: _arr = TAB_REAC.EXTR_TABLE().Array('INTITULE', 'DY') _reac_dy_n = abs(float(_arr[0][1])) if len(_arr) else 0.0 _max_contact_pa = _reac_dy_n / _gland_width_m if _gland_width_m > 0.0 else 0.0 except Exception: import re as _re _raw = str(TAB_REAC) _nums = _re.findall(r"[-+]?\d+\.\d+E[-+]\d+", _raw) _vals = [abs(float(n)) for n in _nums] _reac_dy_n = max(_vals) if _vals else 0.0 _max_contact_pa = _reac_dy_n / _gland_width_m if _gland_width_m > 0.0 else 0.0 print("OG_CONTACT_PRESSURE_PA=%.6E" % _max_contact_pa) print("OG_CONTACT_GATE_MPA=%.4f" % 30.0000) FIN() # ------------------------------------------------------------------------------ Command line #1: ulimit -c unlimited ; ulimit -t 108000 ; ( /opt/spack/opt/spack/linux-zen2/python-3.11.14-crxgllx4ijrk74sukxvebk3ulzuavesl/bin/python3 -m mpi4py "M02_oring.comm.changed.py" --last --link="F::mmed::/work/hardware/mechanical/fem_meshes/ORING_HYPERELASTIC_STANDARD.med::D::20" --tpmax 86400 ; echo $? > _exit_code_ ) 2>&1 | tee -a fort.6 setting '--memory' value to 3686.40 MB (keyword RESERVE_MEMOIRE) copying '/work/hardware/mechanical/fem_meshes/ORING_HYPERELASTIC_STANDARD.med' to 'fort.20'... checking MPI initialization... using COMM_WORLD. MPI is initialized. Ouverture en écriture du fichier ./vola.1 Démarrage de l'exécution. -- CODE_ASTER -- VERSION : EXPLOITATION (stable) -- Version 17.4.0 modifiée le 02/01/2026 révision spack-local - branche 'v17.4.0' Copyright EDF R&D 1991 - 2026 Exécution du : Sun Aug 2 23:46:30 2026 Nom de la machine : 2138298f84ef Architecture : 64bit Type de processeur : x86_64 Système d'exploitation : Linux-6.8.0-124-generic-x86_64-with-glibc2.36 Langue des messages : en (UTF-8) Version de Python : 3.11.14 Version de NumPy : 1.26.4 Parallélisme MPI : actif Rang du processeur courant : 0 Nombre de processeurs MPI utilisés : 1 Parallélisme OpenMP : actif Nombre de processus OpenMP utilisés : 1 Version de la librairie HDF5 : 1.10.7 Version de la librairie MED : 4.2.0 Version de la librairie MFront : 5.0.0 Version de la librairie MUMPS : 5.6.2 Version de la librairie PETSc : 3.24.0p0 Version de la librairie SCOTCH : 7.0.10 starting the execution... Valeur initiale du temps CPU maximum = 86400 secondes Valeur du temps CPU maximum passé aux commandes = 77760 secondes Réserve CPU prévue = 8640 secondes Ouverture en écriture du fichier ./glob.1 Ouverture en écriture du fichier ./vola.1 Ouverture en lecture du fichier /opt/spack/opt/spack/linux-zen2/code-aster-17.4.0-ecm2bfgr5obydotnlte6xilvggppqnap/lib64/aster/elem 1 Nom de la base : ELEMBASE Créée avec la version : 17.04.00 Nombre d'enregistrements utilisés : 46 Nombre d'enregistrements maximum : 512 Nombre d'enregistrements par fichier : 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 /opt/spack/opt/spack/linux-zen2/code-aster-17.4.0-ecm2bfgr5obydotnlte6xilvggppqnap/lib64/aster/elem 1 Nom de la base : ELEMBASE Nombre d'enregistrements utilisés : 46 Nombre d'enregistrements maximum : 512 Nombre d'enregistrements par fichier : 512 Longueur d'enregistrement (octets) : 819200 Nombre total d'accès en lecture : 64 Volume des accès en lecture : 50.00 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.219288 s.) Mémoire limite pour l'allocation dynamique : 4254.17 Mo ajouté à l'initialisation du processus : 647.56 Mo Limite cible du processus : 4901.73 Mo Taille limite des fichiers d'échange : 2048.00 Go # Mémoire (Mo) : 647.56 / 639.44 / 207.61 / 183.67 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0001 user+syst: 0.08s (syst: 0.55s, elaps: 0.76s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt18 # ---------------------------------------------------------------------------------------------- # Commande #0002 de M02_oring.comm.changed.py, ligne 18 MESH = LIRE_MAILLAGE(FORMAT='MED', INFO=1, INFO_MED=0, PARTITIONNEUR='SANS', UNITE=20, VERI_MAIL=_F(APLAT=0.001, VERIF='OUI')) Nom MED du maillage : oring_axis ------------ MAILLAGE 00000001 - IMPRESSIONS NIVEAU 1 ------------ NOMBRE DE NOEUDS 214 NOMBRE DE MAILLES 400 SEG2 48 TRIA3 352 NOMBRE DE GROUPES DE MAILLES 5 -------------------------------------------------------------------------------- Vérification du maillage. # Résultat commande #0002 (LIRE_MAILLAGE): MESH ('<00000001>') de type # Mémoire (Mo) : 647.56 / 640.15 / 207.61 / 183.67 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0002 user+syst: 0.02s (syst: 0.01s, elaps: 0.09s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt22 # ---------------------------------------------------------------------------------------------- # Commande #0003 de M02_oring.comm.changed.py, ligne 22 MODEL = AFFE_MODELE(AFFE=_F(MODELISATION='D_PLAN', PHENOMENE='MECANIQUE', TOUT='OUI'), DISTRIBUTION=_F(METHODE='SOUS_DOMAINE', PARTITIONNEUR='METIS'), INFO=1, MAILLAGE=MESH, VERI_FORMULATION_IFS='OUI', VERI_JACOBIEN='OUI', VERI_NORM_IFS='OUI', VERI_PLAN='OUI') Sur les 400 mailles du maillage 00000001, on a demandé l'affectation de 400, on a pu en affecter 400. Modélisation Formulation Type maille Élément fini Nombre D_PLAN _ TRIA3 MEDPTR3 352 _ _ SEG2 MEPLSE2 48 #2 Calculs elementaires et assemblages CPU (USER+SYST/SYST/ELAPS): 0.00 0.00 0.01 # Résultat commande #0003 (AFFE_MODELE): MODEL ('<00000002>') de type # Mémoire (Mo) : 768.88 / 768.88 / 207.61 / 197.08 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0003 user+syst: 0.02s (syst: 0.01s, elaps: 0.08s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt28 # ---------------------------------------------------------------------------------------------- # Commande #0004 de M02_oring.comm.changed.py, ligne 28 MAT = DEFI_MATERIAU(ELAS_HYPER=_F(C01=140000.0, C10=550000.0, C20=0.0, K=2000000000.0, RHO=0.0), INFO=1) # Résultat commande #0004 (DEFI_MATERIAU): MAT ('<00000003>') de type # Mémoire (Mo) : 768.88 / 768.88 / 207.61 / 197.08 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0004 user+syst: 0.14s (syst: 0.00s, elaps: 0.13s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt32 # ---------------------------------------------------------------------------------------------- # Commande #0005 de M02_oring.comm.changed.py, ligne 32 FIELDMAT = AFFE_MATERIAU(AFFE=_F(MATER=MAT, TOUT='OUI'), INFO=1, MAILLAGE=MESH) # Résultat commande #0005 (AFFE_MATERIAU): FIELDMAT ('<00000004>') de type # Mémoire (Mo) : 768.88 / 768.88 / 207.61 / 197.08 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0005 user+syst: 0.14s (syst: 0.01s, elaps: 0.14s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt39 # ---------------------------------------------------------------------------------------------- # Commande #0006 de M02_oring.comm.changed.py, ligne 39 LOAD_FIXED = AFFE_CHAR_MECA(DDL_IMPO=_F(DX=0.0, DY=0.0, GROUP_MA='GROOVE_FLOOR'), DOUBLE_LAGRANGE='OUI', INFO=1, MODELE=MODEL, SYNTAXE='NON', VERI_AFFE='OUI', VERI_NORM='OUI') # Résultat commande #0006 (AFFE_CHAR_MECA): LOAD_FIXED ('<00000005>') de type # Mémoire (Mo) : 768.88 / 768.88 / 207.61 / 197.08 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0006 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt44 # ---------------------------------------------------------------------------------------------- # Commande #0007 de M02_oring.comm.changed.py, ligne 44 RAMP = DEFI_FONCTION(INFO=1, INTERPOL='LIN', NOM_PARA='INST', NOM_RESU='TOUTRESU', PROL_DROITE='CONSTANT', PROL_GAUCHE='CONSTANT', VALE=(0.0, 0.0, 1.0, 1.0), VERIF='CROISSANT') # Résultat commande #0007 (DEFI_FONCTION): RAMP ('<00000006>') de type # Mémoire (Mo) : 768.88 / 768.88 / 207.61 / 197.08 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0007 user+syst: 0.00s (syst: 0.00s, elaps: 0.01s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt50 # ---------------------------------------------------------------------------------------------- # Commande #0008 de M02_oring.comm.changed.py, ligne 50 LOAD_SQUEEZE = AFFE_CHAR_MECA(DDL_IMPO=_F(DY=0.00098, GROUP_MA='GLAND_FACE'), DOUBLE_LAGRANGE='OUI', INFO=1, MODELE=MODEL, SYNTAXE='NON', VERI_AFFE='OUI', VERI_NORM='OUI') # Résultat commande #0008 (AFFE_CHAR_MECA): LOAD_SQUEEZE ('<00000007>') de type # Mémoire (Mo) : 768.88 / 768.88 / 207.61 / 197.08 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0008 user+syst: 0.01s (syst: 0.00s, elaps: 0.01s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt55 # ---------------------------------------------------------------------------------------------- # Commande #0009 de M02_oring.comm.changed.py, ligne 55 LINST = DEFI_LIST_REEL(DEBUT=0.0, INFO=1, INTERVALLE=_F(JUSQU_A=1.0, NOMBRE=40)) # Résultat commande #0009 (DEFI_LIST_REEL): LINST ('<00000008>') de type # Mémoire (Mo) : 768.88 / 768.88 / 207.61 / 197.08 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0009 user+syst: 0.00s (syst: 0.00s, elaps: 0.00s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt57 # ---------------------------------------------------------------------------------------------- # Commande #0010 de M02_oring.comm.changed.py, ligne 57 DLIST = DEFI_LIST_INST(ADAPTATION=_F(CRIT_COMP='LE', EVENEMENT='SEUIL', MODE_CALCUL_TPLUS='FIXE', NB_INCR_SEUIL=2, NOM_PARA='NB_ITER_NEWTON', PCENT_AUGM=100.0), DEFI_LIST=_F(LIST_INST=LINST, NB_PAS_MAXI=1000000), ECHEC=_F(ACTION='DECOUPE', EVENEMENT='ERREUR', SUBD_METHODE='MANUEL', SUBD_NIVEAU=6, SUBD_PAS=4, SUBD_PAS_MINI=0.0), INFO=1, METHODE='AUTO') Gestion automatique de la liste d'instants. Paramètres de la gestion automatique de la liste d'instants. Pas minimum: 1E-12 Pas maximum: 1 Nombre de pas maximum: 1000000 La liste d'instants contient 41 pas de temps et le pas de temps minimum vaut 0.025. Il y a 1 événements. L'événement 1 est pour capturer les erreurs. Si cet événement se déclenche, on découpe le pas de temps. La découpe du pas de temps est manuelle. On découpe le pas de temps en 4 incréments jusqu'à ce que le niveau atteigne 6. Il y a 1 schémas d'adaptations du pas de temps. Le schéma d'adaptation 1 du pas de temps se déclenche pour un seuil donné. Le seuil est franchi quand, 2 fois de suite, on fait exactement ou moins de 0 itérations de Newton. Le mode de calcul de l'instant suivant est fixe. Le pas de temps suivant sera modifié de 100.00 %. # Résultat commande #0010 (DEFI_LIST_INST): DLIST ('<00000009>') de type # Dépend de : # - LINST ('<00000008>') de type # Mémoire (Mo) : 768.88 / 768.88 / 207.61 / 197.08 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0010 user+syst: 0.02s (syst: 0.00s, elaps: 0.02s) # ---------------------------------------------------------------------------------------------- .. _stg1_txt76 # ---------------------------------------------------------------------------------------------- # Commande #0011 de M02_oring.comm.changed.py, ligne 76 RESU = STAT_NON_LINE(AFFICHAGE=_F(INFO_RESIDU='NON', INFO_TEMPS='NON'), ARCHIVAGE=_F(CHAM_EXCLU=('RESI_NOEU', 'RESI_RELA_NOEU'), PAS_ARCH=1), CHAM_MATER=FIELDMAT, COMPORTEMENT=_F(DEFORMATION='GREEN_LAGRANGE', ITER_CPLAN_MAXI=10, ITER_INTE_MAXI=20, PARM_THETA=1.0, REGU_VISC='NON', RELATION='ELAS_HYPER', RESI_CPLAN_RELA=1e-06, RESI_INTE=1e-06, TOUT='OUI'), CONVERGENCE=_F(ARRET='OUI', ITER_GLOB_ELAS=25, ITER_GLOB_MAXI=30, RESI_GLOB_RELA=1e-05, VERIF='TOUT'), EXCIT=(_F(CHARGE=LOAD_FIXED, TYPE_CHARGE='FIXE_CSTE'), _F(CHARGE=LOAD_SQUEEZE, FONC_MULT=RAMP, TYPE_CHARGE='FIXE_CSTE')), INCREMENT=_F(LIST_INST=DLIST, NUME_INST_INIT=0), MESURE=_F(TABLE='NON'), METHODE='NEWTON', MODELE=MODEL, NEWTON=_F(MATRICE='TANGENTE', MATR_RIGI_SYME='NON', REAC_INCR=1, REAC_ITER=1, REAC_ITER_ELAS=0), SOLVEUR=_F(ACCELERATION='AUTO', ELIM_LAGR='LAGR2', FILTRAGE_MATRICE=-1.0, GESTION_MEMOIRE='AUTO', LOW_RANK_SEUIL=0.0, MATR_DISTRIBUEE='NON', METHODE='MUMPS', MIXER_PRECISION='NON', NB_RHS=1, NPREC=8, PCENT_PIVOT=35, POSTTRAITEMENTS='AUTO', PRETRAITEMENTS='AUTO', REDUCTION_MPI=0, RENUM='AUTO', RESI_RELA=-1.0, STOP_SINGULIER='OUI', TYPE_RESOL='AUTO')) Liste des comportements Affecté sur 400 éléments Relation : ELAS_HYPER Déformation : GREEN_LAGRANGE Pas de régularisation visqueuse Pas de post-traitement à la fin du pas de temps Nombre total de variables internes : 1 V1 : VIDE Le système linéaire à résoudre a 506 degrés de liberté: - 428 sont des degrés de liberté physiques (ils sont portés par 214 noeuds du maillage) - 78 sont les couples de paramètres de Lagrange associés aux 39 relations linéaires dualisées. La matrice est de taille 506 équations. Elle contient 3171 termes non nuls si elle est symétrique et 5836 termes non nuls si elle n'est pas symétrique. Soit un taux de remplissage de 2.279 %. Il n'y a pas d'état initial défini. On prend un état initial nul. Le champ est initialisé a zéro Le champ est initialisé a zéro Le champ est initialisé a zéro Archivage de l'état initial Archivage des champs Champ stocké à l'instant 0.000000000000e+00 pour le numéro d'ordre 0 Champ stocké à l'instant 0.000000000000e+00 pour le numéro d'ordre 0 Champ stocké à l'instant 0.000000000000e+00 pour le numéro d'ordre 0 Champ stocké à l'instant 0.000000000000e+00 pour le numéro d'ordre 0 --------------------------------------------------------------------- Instant de calcul: 2.500000000000e-02 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 6.250000000000e-03) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.073 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.021 s (1 intégrations) * Temps total factorisation matrice : 0.034 s (1 factorisations) * Temps construction second membre : 0.003 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.015 s Mémoire (Mo) : 771.65 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 6.250000000000e-03 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 1.562500000000e-03) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.020 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.010 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.002 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 1.562500000000e-03 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 3.906250000000e-04) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.018 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.009 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.002 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 3.906250000000e-04 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 9.765625000000e-05) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.016 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.007 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 9.765625000000e-05 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 2.441406250000e-05) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.018 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.009 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 2.441406250000e-05 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 6.103515625000e-06) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.017 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 6.103515625000e-06 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 1.525878906250e-06) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.017 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.007 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.002 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 1.525878906250e-06 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 3.814697265625e-07) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.014 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.007 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 3.814697265625e-07 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 9.536743164063e-08) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.016 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.002 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 9.536743164063e-08 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 2.384185791016e-08) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.018 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.002 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 2.384185791016e-08 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 5.960464477539e-09) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.016 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.007 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.002 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 5.960464477539e-09 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 1.490116119385e-09) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.018 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 1.490116119385e-09 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 3.725290298462e-10) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.019 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.002 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.006 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 3.725290298462e-10 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 9.313225746155e-11) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.017 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 9.313225746155e-11 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 2.328306436539e-11) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.019 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.009 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.005 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 2.328306436539e-11 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 5.820766091347e-12) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.018 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.009 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.002 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.005 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 5.820766091347e-12 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 1.455191522837e-12) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.016 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.007 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 1.455191522837e-12 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 3.637978807092e-13) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.017 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 3.637978807092e-13 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 9.094947017729e-14) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.017 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.004 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 9.094947017729e-14 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 2.273736754432e-14) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.017 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 2.273736754432e-14 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 5.684341886081e-15) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.016 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 5.684341886081e-15 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 1.421085471520e-15) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.016 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.008 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 1.421085471520e-15 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 3.552713678801e-16) On découpe le pas de temps. Temps CPU consommé dans ce pas de temps : 0.019 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.007 s (1 intégrations) * Temps total factorisation matrice : 0.005 s (1 factorisations) * Temps construction second membre : 0.002 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) --------------------------------------------------------------------- Instant de calcul: 3.552713678801e-16 --------------------------------------------------------------------- | NEWTON | RESIDU | RESIDU | OPTION | | ITERATION | RELATIF | ABSOLU | ASSEMBLAGE | | | RESI_GLOB_RELA | RESI_GLOB_MAXI |(+ TYPE MATRICE)| --------------------------------------------------------------------- Degré de liberté physique associé au noeud 6 et à la composante DX. Problème : la matrice est singulière ou presque singulière : Lors de la factorisation de la matrice, on a rencontré un problème (pivot nul ou presque nul) à la ligne 19 qui correspond au degré de liberté donné ci-dessus. Risques et conseils : * Si la ligne correspond a un degré de liberté physique, il s'agit probablement d'un mouvement de corps rigide mal bloqué. Vérifiez les conditions aux limites. Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le contact. Vérifiez également les caractéristiques matériaux (module d'Young, ...). * Si la ligne correspond a un degré de liberté de Lagrange, il s'agit sans doute d'une condition limite redondante. En particulier, il se peut que la relation linéaire surabondante provienne des conditions de contact. Peut-être devriez vous exclure certains noeuds des conditions de contact (mots clés SANS_NOEUD et SANS_GROUP_NO). * Si le solveur utilisé est LDLT ou MULT_FRONT, vous pouvez utiliser le solveur MUMPS car celui-ci est le seul à pouvoir factoriser les matrices qui ne sont pas définies positives. * Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop pessimiste ! Il reste néanmoins souvent possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 la valeur du mot-clé NPREC) ou en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de processeurs. * Il se peut aussi que ce phénomène soit tout à fait normal avec X-FEM si la fissure passe très près d'un noeud. Si le nombre de décimales perdues n'est pas trop grand (max 10 décimales), vous pouvez relancer le calcul en augmentant le nombre de décimales perdues autorisé : mot-clé NPREC du mot clé facteur SOLVEUR. Sinon, contactez l'équipe de développement. | 0 E | | |TANGENTE (SYM)| --------------------------------------------------------------------- La matrice du système est singulière On essaie de découper le pas de temps. On utilise la découpe manuelle. Découpe uniforme à partir de l'instant 0.000000000000e+00 en 4 pas de temps. (soit un incrément constant de 8.881784197001e-17) Le pas de temps minimum 0.000000000000e+00 (SUBD_PAS_MINI) est atteint. <Échec> Échec dans la tentative de découper le pas de temps. Temps CPU consommé dans ce pas de temps : 0.016 s * Nombre d'itérations de Newton : 1 * Temps total intégration comportement : 0.007 s (1 intégrations) * Temps total factorisation matrice : 0.003 s (1 factorisations) * Temps construction second membre : 0.001 s * Temps total résolution K.U=F : 0.000 s (0 résolutions) * Temps assemblage matrice : 0.000 s * Temps autres opérations : 0.004 s Mémoire (Mo) : 771.81 / 771.65 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) ╔════════════════════════════════════════════════════════════════════════════════════════════════╗ ║ ║ ║ ║ ║ Arrêt pour cause de matrice non inversible. ║ ║ La base globale est sauvegardée. Elle contient les pas archivés avant l'arrêt. ║ ║ ║ ║ Conseils : ║ ║ - Vérifiez vos conditions aux limites. ║ ║ - Vérifiez votre modèle, la cohérence des unités. ║ ║ - Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le ║ ║ contact. ║ ║ ║ ║ - Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop ║ ║ pessimiste ! Il reste néanmoins souvent ║ ║ possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 ║ ║ la valeur du mot-clé NPREC) ou ║ ║ en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de ║ ║ processeurs. ║ ╚════════════════════════════════════════════════════════════════════════════════════════════════╝ Validation du concept '0000000a'. # Résultat commande #0011 (STAT_NON_LINE): RESU ('<0000000a>') de type # Mémoire (Mo) : 771.81 / 770.32 / 207.61 / 197.63 (VmPeak / VmSize / Optimum / Minimum) # Fin commande #0011 user+syst: 0.54s (syst: 0.06s, elaps: 0.71s) # ---------------------------------------------------------------------------------------------- Traceback (most recent call last): File "", line 198, in _run_module_as_main File "", line 88, in _run_code File "/opt/spack/var/spack/environments/simvia_env/.spack-env/view/lib/python3.11/site-packages/mpi4py/__main__.py", line 7, in main() File "/opt/spack/var/spack/environments/simvia_env/.spack-env/view/lib/python3.11/site-packages/mpi4py/run.py", line 230, in main run_command_line(args) File "/opt/spack/var/spack/environments/simvia_env/.spack-env/view/lib/python3.11/site-packages/mpi4py/run.py", line 47, in run_command_line run_path(sys.argv[0], run_name='__main__') File "", line 291, in run_path File "", line 98, in _run_module_code File "", line 88, in _run_code File "M02_oring.comm.changed.py", line 76, in RESU = STAT_NON_LINE( ^^^^^^^^^^^^^^ File "/opt/spack/opt/spack/linux-zen2/code-aster-17.4.0-ecm2bfgr5obydotnlte6xilvggppqnap/lib64/aster/code_aster/Supervis/ExecuteCommand.py", line 183, in run return cmd.run_(**kwargs) ^^^^^^^^^^^^^^^^^^ File "/opt/spack/opt/spack/linux-zen2/code-aster-17.4.0-ecm2bfgr5obydotnlte6xilvggppqnap/lib64/aster/code_aster/Supervis/ExecuteCommand.py", line 226, in run_ self.exec_(keywords) File "/opt/spack/opt/spack/linux-zen2/code-aster-17.4.0-ecm2bfgr5obydotnlte6xilvggppqnap/lib64/aster/code_aster/Supervis/ExecuteCommand.py", line 516, in exec_ self._call_oper(syntax) File "/opt/spack/opt/spack/linux-zen2/code-aster-17.4.0-ecm2bfgr5obydotnlte6xilvggppqnap/lib64/aster/code_aster/Supervis/ExecuteCommand.py", line 282, in _call_oper return libaster.call_oper(syntax, 0) ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ libaster.SolverError: ╔════════════════════════════════════════════════════════════════════════════════════════════════╗ ║ ║ ║ ║ ║ Arrêt pour cause de matrice non inversible. ║ ║ La base globale est sauvegardée. Elle contient les pas archivés avant l'arrêt. ║ ║ ║ ║ Conseils : ║ ║ - Vérifiez vos conditions aux limites. ║ ║ - Vérifiez votre modèle, la cohérence des unités. ║ ║ - Si vous faites du contact, il ne faut pas que la structure ne "tienne" que par le ║ ║ contact. ║ ║ ║ ║ - Parfois, en parallèle, le critère de détection de singularité de MUMPS est trop ║ ║ pessimiste ! Il reste néanmoins souvent ║ ║ possible de faire passer le calcul complet en relaxant ce critère (augmenter de 1 ou 2 ║ ║ la valeur du mot-clé NPREC) ou ║ ║ en le débranchant (valeur du mot-clé NPREC=-1) ou en relançant le calcul sur moins de ║ ║ processeurs. ║ ╚════════════════════════════════════════════════════════════════════════════════════════════════╝ Exception: Trying to close the database after an uncaught exception... Publishing the result of the current command STAT_NON_LINE... No database in results, objects not saved on processor #0 ╔════════════════════════════════════════════════════════════════════════════════════════════════╗ ║ ║ ║ ║ ║ 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 : ║ ║ aucune alarme ║ ╚════════════════════════════════════════════════════════════════════════════════════════════════╝ ARRET NORMAL DANS "FIN" PAR APPEL A "JEFINI". MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION : 197.63 Mo MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION : 207.61 Mo MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS LORS DE L'EXECUTION : 779.02 Mo FERMETURE DES BASES EFFECTUEE STATISTIQUES CONCERNANT L'ALLOCATION DYNAMIQUE : TAILLE CUMULEE MAXIMUM : 208 Mo. TAILLE CUMULEE LIBEREE : 10 Mo. NOMBRE TOTAL D'ALLOCATIONS : 53452 NOMBRE TOTAL DE LIBERATIONS : 52919 APPELS AU MECANISME DE LIBERATION : 0 TAILLE MEMOIRE CUMULEE RECUPEREE : 0 Mo. VOLUME DES LECTURES : 0 Mo. VOLUME DES ECRITURES : 1 Mo. MEMOIRE JEVEUX MINIMALE REQUISE POUR L'EXECUTION : 197.63 Mo - IMPOSE DE NOMBREUX ACCES DISQUE - RALENTIT LA VITESSE D'EXECUTION MEMOIRE JEVEUX OPTIMALE REQUISE POUR L'EXECUTION : 207.61 Mo - LIMITE LES ACCES DISQUE - AMELIORE LA VITESSE D'EXECUTION MAXIMUM DE MEMOIRE UTILISEE PAR LE PROCESSUS : 779.02 Mo - COMPREND LA MEMOIRE CONSOMMEE PAR JEVEUX, LE SUPERVISEUR PYTHON, LES LIBRAIRIES EXTERNES FIN D'EXECUTION LE : DI-02-AOUT-2026 23:46:33 INFO './glob.1' deleted Deleting './glob.2': No such file or directory INFO './vola.1' deleted Deleting './vola.2': No such file or directory -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD Proc: [[9583,0],0] Errorcode: 1 NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes. You may or may not see output from other processes, depending on exactly when Open MPI kills them. -------------------------------------------------------------------------- EXECUTION_CODE_ASTER_EXIT_48=1 execution ended (command file #1): _NO_CONVERGENCE # ------------------------------------------------------------------------------ Content of /tmp/run_aster_8_qv4h3g/proc.0 after execution: .: total 176 -rw-r--r-- 1 root root 259 Aug 2 23:46 48.export -rw-r--r-- 1 root root 5511 Aug 2 23:46 M02_oring.comm.changed.py drwxr-xr-x 2 root root 4096 Aug 2 23:46 REPE_IN drwxr-xr-x 2 root root 4096 Aug 2 23:46 REPE_OUT -rw-r--r-- 1 root root 38065 Jul 3 02:01 fort.20 -rw-r--r-- 1 root root 115399 Aug 2 23:46 fort.6 -rw-r--r-- 1 root root 0 Aug 2 23:46 fort.8 -rw-r--r-- 1 root root 0 Aug 2 23:46 fort.9 REPE_OUT: total 0 # ------------------------------------------------------------------------------ Copying results copying 'fort.6' to '/out/M02_oring.mess'... WARNING: file not found: fort.80 # ------------------------------------------------------------------------------ Execution summary cpu system cpu+sys elapsed -------------------------------------------------------------------------------- Preparation of environment 0.00 0.00 0.00 0.01 Execution of code_aster 6.42 1.92 8.34 14.60 Copying results 0.00 0.00 0.00 0.01 -------------------------------------------------------------------------------- Total 6.42 1.92 8.34 14.62 -------------------------------------------------------------------------------- ------------------------------------------------------------ --- DIAGNOSTIC JOB : _NO_CONVERGENCE ------------------------------------------------------------