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

You are not logged in.

#1 2018-12-22 10:20:23

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

[solved] matrix factorization error and duplicate nodes

I am getting this error. I tried the recommended solutions but none of them works. Any hint for the solution? Mesh file is attached.   Comm file is below.


   !------------------------------------------------------------------------------------------------------------------------------!
   ! <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 1277 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.                                                                             !
   !------------------------------------------------------------------------------------------------------------------------------!

------COMM  FILE.......................


DEBUT(LANG='EN')

mesh = LIRE_MAILLAGE(FORMAT='MED',
                     UNITE=20)

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

concrete = DEFI_MATERIAU(ELAS=_F(E=2.1e+11,
                                 NU=0.2,
                                 RHO=7850.0))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(concrete, ),
                                 TOUT='OUI'),
                         MODELE=model)

fload = DEFI_FONCTION(INTERPOL=('LIN', ),
                      NOM_PARA='INST',
                      PROL_DROITE='CONSTANT',
                      PROL_GAUCHE='CONSTANT',
                      VALE=(0.0, 0.0, 10.0, 1.0))

listr = DEFI_LIST_REEL(DEBUT=0.0,
                       INTERVALLE=_F(JUSQU_A=10.0,
                                     PAS=1.0))

support = AFFE_CHAR_MECA(DDL_IMPO=(_F(DX=0.0,
                                      DY=0.0,
                                      DZ=0.0,
                                      GROUP_NO=('supportnode', )),
                                   _F(DX=0.0,
                                      GROUP_NO=('symmetryFace', ))),
                         MODELE=model)

load0 = AFFE_CHAR_MECA(FORCE_NODALE=_F(FZ=-50000.0,
                                       GROUP_NO=('loadingNode', )),
                       MODELE=model)

reslin = MECA_STATIQUE(CHAM_MATER=fieldmat,
                       EXCIT=(_F(CHARGE=load0,
                                 FONC_MULT=fload),
                              _F(CHARGE=support)),
                       MODELE=model,
                       SOLVEUR=_F(METHODE='MUMPS'))

mass = POST_ELEM(MASS_INER=_F(TOUT='OUI'),
                 TITRE='mass')

IMPR_TABLE(FORMAT='TABLEAU',
           TABLE=mass,
           UNITE=8)

FIN()

Last edited by nirmaljoshi (2019-01-10 10:23:07)


Attachments:
half_beam.hdf, Size: 629.6 KiB, Downloads: 21

Offline

#2 2018-12-22 11:25:05

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

Re: [solved] matrix factorization error and duplicate nodes

I tried by appling load and support to faces instead of node, but still same error.


Attachments:
halfbeam.zip, Size: 49.5 KiB, Downloads: 18

Offline

#3 2018-12-23 22:07:03

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

Re: [solved] matrix factorization error and duplicate nodes

Hello Nirmaljoshi,
There is one DOF which is not constrained: the X rotation around the node 'supportnode'.

Even after adding a second node to 'supportnode' in order to constrain that DOF, the matrix remains singular.

I tried to re-mesh your geometry and it works. I don't know why your mesh doesn't work, there are some elements with poor aspect ratio (>6) but I don't think that should make the matrix singular.

I decreased your nodal force from 50000 to 50 because the deformation was huge.

I attach the old mesh, the modified mesh and the modified comm file.

Any ideas why the first mesh doesn't work are welcome.

Good luck,

Konyaro


Attachments:
2018_12_23_beam.zip, Size: 75.33 KiB, Downloads: 14

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

Offline

#4 2018-12-24 08:29:51

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

Re: [solved] matrix factorization error and duplicate nodes

hello both

it looks like the mesh contains a few hundreds of double node
the number depends on the precision chosen

how was  it made?

jean pierre aubry

Offline

#5 2019-01-08 06:11:59

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

Re: [solved] matrix factorization error and duplicate nodes

Yes, I had also suspected that the error is due to duplicate nodes.

a) How can we avoid creation of such duplicate nodes in the Mesh model? Or (b) merge them if already present in mesh?

Offline

#6 2019-01-08 09:43:28

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

Re: [solved] matrix factorization error and duplicate nodes

a) How can we avoid creation of such duplicate nodes in the Mesh model? Or (b) merge them if already present in mesh?

answer to b)
in the salome_meca Mesh module
menu Modification > Transformation > Merge nodes
the in the dialog box Manual
and looks what is going on

and answer to a)
but i am afraid this is more or less hopeless as the geometry is probably very badly defined to start with

Offline

#7 2019-01-09 07:01:42

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

Re: [solved] matrix factorization error and duplicate nodes

Thanks Jean!

I found that I was using Compound method to join the objects. It was a mistake and perhaps this created multiple nodes. I used Operation>Partition method after forming the coumpund and the problem of duplicate nodes has been solved. (referece- docs.salome-platform.org/latest/gui/GEOM/partition_explanation.html)

But the matrix factorization error still occurs.  Konyaro's mesh works fine. But when I create new mesh, it dosent work.

What is the problem?

Offline

#8 2019-01-09 07:03:50

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

Re: [solved] matrix factorization error and duplicate nodes

the same file with new mesh-saying factoriztion error


Attachments:
half_beam.hdf, Size: 1.08 MiB, Downloads: 6

Offline

#9 2019-01-09 08:39:20

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

Re: [solved] matrix factorization error and duplicate nodes

there is no .comm file !

Offline

#10 2019-01-09 11:07:40

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

Re: [solved] matrix factorization error and duplicate nodes

here is the file.


Attachments:
meshAndComm.zip, Size: 30.8 KiB, Downloads: 8

Offline

#11 2019-01-09 11:45:53

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

Re: [solved] matrix factorization error and duplicate nodes

group_supportNode consist of one single node that belongs to a 3D (solid) mesh
it carries no rotation DOF so the model is able to rotate freely around the 3 axes
this is basic modeling and is not specific to code_aster

Offline

#12 2019-01-10 10:22:38

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

Re: [solved] matrix factorization error and duplicate nodes

Yes Jean, thanks. That was my silly mistake. I completely forgot about the rotational constrainsts.

Offline

Board footer