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

You are not logged in.

#1 2021-10-21 17:07:11

CASIO
Member
Registered: 2021-10-05
Posts: 15

I get an error in the POU_D_E model, please tell me what is wrong

I created a model with reference to framabook, but I get the following error.

!------------------------------------------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <FACTOR_10>                                                                                                               
   !                                                                                                                                                                                 
   ! 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 43 qui correspond au degré de liberté donné ci-dessus.                     
   !   Le nombre de décimales "perdues" pour cette ligne est : 9.                                                                           
   !   Ce nombre de décimales perdues est lié au degré de singularité de la matrice.                                             
   !   Plus il est grand, plus le système est singulier.                                                                         
   !   Quand on a perdu plus de 8 décimales, on estime que l'on a perdu 50% de la précision                                       
   !   des nombres réels (qui ont 15 décimales environ).                                                                         
   !                                                                                                                             
   ! 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.                                                                             
   !------------------------------------------------------------------------------------------------------------------------------!

Using AFFE_CARA_MECA's ORIENTATION did not eliminate the error
Please tell me where I made a mistake in Geometry, Mesh or Aster Study.

I pasted the .comm, .med and .hdf file, so please answer


Attachments:
POU_D_E.zip, Size: 29.77 KiB, Downloads: 31

Offline

#2 2021-10-24 10:50:41

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

Re: I get an error in the POU_D_E model, please tell me what is wrong

hello

i am the author of what you call the Framabook
i had like to be called by my author name

1- the group of discrete element hinge is not defined in the the model
2- each hinge should be only one element in its length

i suggest you try to run the example in my book strictly as it is written to understand what is going on
then, and then only you can try to invent things

jean pierre aubry


consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/

Offline

#3 2021-10-25 15:06:58

CASIO
Member
Registered: 2021-10-05
Posts: 15

Re: I get an error in the POU_D_E model, please tell me what is wrong

thank you for your reply.

I reread the framabook and tried to make a model that included hinges, but I got a similar error again.

It will bother you again, but please tell me what is wrong


Attachments:
hinge2.zip, Size: 41.06 KiB, Downloads: 24

Offline

#4 2021-10-25 15:24:40

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

Re: I get an error in the POU_D_E model, please tell me what is wrong

i told you in the previous post

2- each hinge should be only one element in its length
think of a folding ruler!

try to follow my instruction

and the the framabook is "Beginning with code_aster" by jean-Pierre Aubry


consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/

Offline

#5 2021-11-02 08:29:03

CASIO
Member
Registered: 2021-10-05
Posts: 15

Re: I get an error in the POU_D_E model, please tell me what is wrong

Thank you for your reply.

I tried it with a different version of salome-meca and it succeeded without any errors.

But the answer doesn't match the actual value.
I think there is a problem with Geometry, but is it meshing this structure as one model or as two structures, mast and topbeam?

If it is a 3D model, the former shares nodes at the interface and the latter does not share nodes at the interface.

I'm sorry I asked you many times.
I would appreciate it if you could answer.

Offline

#6 2021-11-02 08:57:56

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

Re: I get an error in the POU_D_E model, please tell me what is wrong

hello

how is made the geometry ?
how  is it meshed then ?

jean pierre aubry


consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/

Offline

#7 2021-11-03 02:41:08

CASIO
Member
Registered: 2021-10-05
Posts: 15

Re: I get an error in the POU_D_E model, please tell me what is wrong

Thank you for your reply.

I'm sorry.
I couldn't accurately reproduce the contents of your book because there was something I couldn't do with the version of salome-meca I'm using now.

The file was posted in the previous question.
Write briefly what I did

1 Decide a node and connect it with a line.
2 Make the connected lines into one model using Partition.
3 Use Create Group to define mast, topbeam, load, etc.
4 mesh the model (the model that was unified using Partition)

Offline

#8 2021-11-03 09:53:29

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

Re: I get an error in the POU_D_E model, please tell me what is wrong

there can be very large differences in between what one believes one is doing and what one is actually doing
particularly if the one is a beginner
so i cannot say more than what i said in post #4

But the answer doesn't match the actual value.

what do you mean?
please try to be clear!!


consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/

Offline

#9 2021-11-04 04:22:43

CASIO
Member
Registered: 2021-10-05
Posts: 15

Re: I get an error in the POU_D_E model, please tell me what is wrong

thank you for your reply.

Carefully, I start over from the beginning.

Finally, I have a question.
In my salome-meca version, typing CARA ='K_TR_D_L', VALE = (1e6,1e6,1e6,0,1e9,1e9,)) gives an error, probably because this is DRX = 0. So, when I tried to change it into DRX = 1, the error disappeared.

Also, when you answered another questioner about hinge, you tiped CARA ='K_TR_D_L', VALE = (1e3,1e3,1e3,1e0,1e0,1e0,1e0).
I think it's a mistake to have 7 values instead of 6, but why is it different from what is written in the book?

What kind of value should I enter to get the correct value?
I want you to answer if you would like.

Thank you for answering the questions so far.
Also, if I have any questions, I would appreciate it if you could answer them.

Offline

#10 2021-11-04 09:29:30

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

Re: I get an error in the POU_D_E model, please tell me what is wrong

with such a simple (isostatic) system DRX=0 make the system unstable so DRX=1 is a good value

the seven values is a typing error it should be 6, of course


consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/

Offline

#11 2021-11-05 09:07:54

CASIO
Member
Registered: 2021-10-05
Posts: 15

Re: I get an error in the POU_D_E model, please tell me what is wrong

Thank Koyu for reply.

I hope we can meet and you are able to assist me once more.

Offline