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

You are not logged in.

#1 2016-09-12 19:42:59

GuillermoBarraza
Member
Registered: 2016-05-24
Posts: 49

[SOLVED]Thermal contact between two concentric cylinders ECHANGE_PAROI

Hello, everyone.

I'm having troubles running an analysis of two concentric cylinders sharing a common face. I don't know if this has been posted before but here it goes: I have an imposed temperature on the bottom part of the cylinders, a convection condition in the external face of the outer cylinder and the heat transfer between internal faces. When running the case the following error comes up:

JDC.py : ERREUR A L'EXECUTION - INTERRUPTION
>> JDC.py : DEBUT RAPPORT
CR d'execution de JDC en MIXTE
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee.                            !
   ! Les bases sont fermees.                                                          !
   ! Type de l'exception : error                                                      !
   !                                                                                  !
   ! 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 31 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.                                                              !
   !                                                                                  !
   !    * 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.                                 !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR d'execution de JDC en MIXTE

I've found the previous alarm several times, some of them being a lack of essential boundary conditions. I have the idea that the usage of EXANGE_PAROI demands two congruent meshes with a clearance between both, but I don't know how this could be done in this particular case. I attach the .hdf, .comm and .mess files and a couple of imaegs hopping that someone can see what's going on.

best,

Guillermo

Last edited by GuillermoBarraza (2016-11-17 23:34:22)


Attachments:
Cylinders.tar.gz, Size: 1.69 MiB, Downloads: 287

Offline

#2 2016-09-12 21:12:08

AsterO'dactyle
Administrator
Registered: 2007-11-29
Posts: 416

Re: [SOLVED]Thermal contact between two concentric cylinders ECHANGE_PAROI

Hello,

I think your problem is not from boundary conditions but from your mesh.

You have this (very important) alarm:
MODELE1_63

That means you have "skin" elements without volumic ones with same nodes at surface. You "skin" elements are "alone"

Check your mesh in Salome


Code_Asterの開発者

Offline

#3 2016-09-14 19:42:43

GuillermoBarraza
Member
Registered: 2016-05-24
Posts: 49

Re: [SOLVED]Thermal contact between two concentric cylinders ECHANGE_PAROI

Thanks for the reply.

I checked my mesh and didn't find anything unusual (at least at first sight). Although, I've found a similar post here:

http://code-aster.org/forum2/viewtopic.php?id=20134

With the same error as mine, but with mechanic loads, Jean Pierre Aubry recomends to add the surface mesh elements (with GROUP_MA) to AFFE_MODELE and avoid the implementation of TOUT='OUI'. Did it but still have the same result.

I don't know what could be the problem.

Guillermo

Offline

#4 2016-11-17 23:33:29

GuillermoBarraza
Member
Registered: 2016-05-24
Posts: 49

Re: [SOLVED]Thermal contact between two concentric cylinders ECHANGE_PAROI

It is solved. I just had to add a line to each of the contact surfaces in order to partition both in halves cylinders. It turns out that there's a bug in curved surfaces regardless the congruent mesh in the interface.

Offline