Hi, i'm having some issues with LIAISON_MAIL for my 3D thermal analysis.
For my project i'm using YACS so everything is automated but i'm trying to solve the issue directly on code_aster. The geometry is quite complex and is composed of 29 solids. I have a python script that detect all the contacts between the parts and creates a surface group on the GEOM part for each contact detected. Then the parts are meshed, then coumpound mesh with all the surface and volume groups for the contact. The meshes are non-conformal that's why i use LIAISON_MAIL.
Unfortunately i'm having this error message :
<EXCEPTION> <FACTOR_11> ║
║ ║
║ 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 29151 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. ║
║
The error message indicates wich nodes causes the error so i went to see the location of the error and nothing in particular except the nodes that cause issues are usually located on a border and really close to a "master" node.
I don't really understand why it happens, i tried isolating random nodes and it works but some nodes cause the singularity error.
Do you have more info about the error ? And do you have any idea how to correct it ?
Feel free to ask questions if you want any more info.
Thanks in advance !
