lecture du maillage dans Code-Aster

  • bevilacqua
  • Topic Author
  • Offline
  • New Member
  • New Member
More
11 years 7 months ago - 11 years 7 months ago #6446 by bevilacqua
lecture du maillage dans Code-Aster was created by bevilacqua
Bonjour,

Je débute dans l'utilisation de Salomé et de code_aster. C'est en cherchant un code sur la toile que j'ai découvert CAELinux que je trouve très convivial avec le duo Salomé/code_aster. Les exemples proposés sont très intéressant et permettent d'entrer dans l'utilisation du code. Un grand merci à ceux qui ont fait ce site et qui le maintienne.
J'ai commencé par étudier les exemples sur le site CAELinux qui m'aident beaucoup dans l'apprentissage de ce code. Actuellement j'étudie l'exemple de la contribution de Kees Wouters : prestressmodal, que je trouve très intéressant car il aborde l'analyse modale sans et avec une force appliquée et le flambage.

J'ai commencé par la création du cylindre :
Création du point null (centrale bas (0,0,0) du cylindre)
Création d'un cercle puis une extrusion pour obtenir le cylindre finale.
Création des groupe géométriques Cbot et Ctop.
J'ai continué par la création du maillage :
Quadrangle (Mapping)
Création des groupe de mailles - Edge MCbot et MCtop - Face MSurfin
Puis création du fichier commande pour Aster et lancement du calcul.
Mais le processus s'arrête juste après la lecture du maillage sur une erreur :
EXECUTION_CODE_ASTER_EXIT_0010-2845-albert-desktop=139 et je n'arrive pas à comprendre le problème et y remédier pour continuer.

Pouvez-vous m'aider sur ce sujet afin que j'aille jusqu'au bout de l'exercice.

Merci

PS : Je vous joins les fichier hdf, med, comm et mess.

English version

Hello,

I started to use Salomé and Code_Aster. It is in looking for a mechanical calculation code on the web that I discovered CAELinux which is very friendly with the duo Salomé / Code_Aster. The examples are very interesting and permit to enter in the code. A big thank you to all who have made this site and continue to maintain it.
I started by studying the examples on the site CAELinux that help me a lot in learning the code. Currently I am studying the example of the contribution of Kees Wouters: prestressmodal, which I find very interesting because it talks about modal analysis with and without an applied force and buckling.

I started by creating the cylinder:
null point is created (bottom center (0,0,0) of the cylinder)
Creating a circle and extruding of it to get the final cylinder.
Creating of geometric group Cbot and Ctop.
I continued with the creation of the mesh:
Quadrangle (Mapping)
Creating mesh group - Edge MCbot and MCtop - Face MSurfin
And control file creation and launch of Aster calculation.
But the process stops after reading the mesh on an error:
EXECUTION_CODE_ASTER_EXIT_0010-2845-albert-desktop = 139 and I can not seem to understand the problem and fix it to continue.

Can you help me on this subject so I go to the end of the year.

thank you

PS: I attache the file hdf, med and comm mess.
Attachments:
Last edit: 11 years 7 months ago by bevilacqua.

Please Log in or Create an account to join the conversation.

More
11 years 6 months ago #6447 by kwou
Replied by kwou on topic Re: lecture du maillage dans Code-Aster
Bonjour Albert - vingt quatre

first of all: I will need to update the examples. If I have time and when the site allows it once again.

A few remarks:
COQUE_3D needs quadratic elements. You have QUAD4 elements:

NOMBRE DE NOEUDS 1640

NOMBRE DE MAILLES 1720
SEG2 120
QUAD4 1600

NOMBRE DE GROUPES DE MAILLES 3
MCbot 40
MCtop 40
MSurfin 1600

Coque_3d can not create an element from this.
You get this error message:
Alarme utilisateur dans CREA_MAILLAGE/MODI_MAILLE :
Occurrence du mot clé facteur MODI_MAILLE : 1.
Vous avez demandé la transformation de type QUAD8_9.
Mais il n'y a aucune maille à transformer.

Then it returns by FIN(), at least in my calculation.

Second remark:
meshmod=CREA_MAILLAGE(MAILLAGE=meshinit,
MODI_MAILLE=(_F(TOUT='OUI',
OPTION='QUAD8_9',),
_F(TOUT='OUI',
OPTION='TRIA6_7',),),);

In the upgrade process of C-A you can only change elements that are available, in your case QUAD8 (after making quad4 qaudratic). Then leave out the _F(TOUT='OUI',OPTION='TRIA6_7',),),) part because you donot have any triangles.

meshmod=CREA_MAILLAGE(MAILLAGE=meshinit,
MODI_MAILLE=(_F(TOUT='OUI',
OPTION='QUAD8_9',),),);

Third: change GROUP_MA='Surfin' to GROUP_MA='MSurfin':
meshmod=MODI_MAILLAGE(reuse =meshmod,
MAILLAGE=meshmod,
ORIE_NORM_COQUE=_F(GROUP_MA='MSurfin',),);

Fourth, roughly the same as third:
clamp=AFFE_CHAR_MECA(MODELE=modelc,
DDL_IMPO=(_F(GROUP_MA='MCbot',
DX=0.0,
DY=0.0,
DZ=0.0,
DRX=0.0,
DRY=0.0,
DRZ=0.0,),
_F(GROUP_MA='MCtop',
DX=0.0,
DY=0.0,
DRX=0.0,
DRY=0.0,
DRZ=0.0,),),
FORCE_ARETE=_F(GROUP_MA='MCtop',FZ=(Fax / Lcyl),),);



Lastly I modified the postprocessing a bit, then it runs.

Note that the retrieval of stresses, strains etc has been completely overhauled in the latter versions of C-A.
So if you have a working version, pls let me know. I will try to update the files.

So far, so good.
If you have further problems, pls let us know.

Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)

--
kind regards - kees

Please Log in or Create an account to join the conversation.

More
11 years 6 months ago #6448 by kwou
Replied by kwou on topic Re: lecture du maillage dans Code-Aster

keeswouters wrote: Bonjour Albert - vingt quatre

first of all: I will need to update the examples. If I have time and when the site allows it once again.

A few remarks:
COQUE_3D needs quadratic elements. You have QUAD4 elements:

NOMBRE DE NOEUDS 1640

NOMBRE DE MAILLES 1720
SEG2 120
QUAD4 1600

NOMBRE DE GROUPES DE MAILLES 3
MCbot 40
MCtop 40
MSurfin 1600

Coque_3d can not create an element from this.
You get this error message:
Alarme utilisateur dans CREA_MAILLAGE/MODI_MAILLE :
Occurrence du mot clé facteur MODI_MAILLE : 1.
Vous avez demandé la transformation de type QUAD8_9.
Mais il n'y a aucune maille à transformer.

Then it returns by FIN(), at least in my calculation.

Second remark:

meshmod=CREA_MAILLAGE(MAILLAGE=meshinit,
MODI_MAILLE=(_F(TOUT='OUI',
OPTION='QUAD8_9',),
_F(TOUT='OUI',
OPTION='TRIA6_7',),),);


In the upgrade process of C-A you can only change elements that are available, in your case QUAD8 (after making quad4 qaudratic). Then leave out the _F(TOUT='OUI',OPTION='TRIA6_7',),),) part because you donot have any triangles.

meshmod=CREA_MAILLAGE(MAILLAGE=meshinit,
MODI_MAILLE=(_F(TOUT='OUI',
OPTION='QUAD8_9',),),);

Third: change GROUP_MA='Surfin' to GROUP_MA='MSurfin':
meshmod=MODI_MAILLAGE(reuse =meshmod,
MAILLAGE=meshmod,
ORIE_NORM_COQUE=_F(GROUP_MA='MSurfin',),);

Fourth, roughly the same as third:
clamp=AFFE_CHAR_MECA(MODELE=modelc,
DDL_IMPO=(_F(GROUP_MA='MCbot',
DX=0.0,
DY=0.0,
DZ=0.0,
DRX=0.0,
DRY=0.0,
DRZ=0.0,),
_F(GROUP_MA='MCtop',
DX=0.0,
DY=0.0,
DRX=0.0,
DRY=0.0,
DRZ=0.0,),),
FORCE_ARETE=_F(GROUP_MA='MCtop',FZ=(Fax / Lcyl),),);



Lastly I modified the postprocessing a bit, then it runs.

Note that the retrieval of stresses, strains etc has been completely overhauled in the latter versions of C-A.
So if you have a working version, pls let me know. I will try to update the files.

So far, so good.
If you have further problems, pls let us know.


Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)

--
kind regards - kees

Please Log in or Create an account to join the conversation.

  • bevilacqua
  • Topic Author
  • Offline
  • New Member
  • New Member
More
11 years 6 months ago #6449 by bevilacqua
Replied by bevilacqua on topic Re: lecture du maillage dans Code-Aster
Hi Kees,

Thank you for your answer.

I have taken into account your corrections and I have understood them. Nevertheless, the problem of the mesh reading (LIRE_MAILLAGE) remains.
I have CAELinux Salome-Meca 2011.2 I had installed in a computer with a quadricorps cpu and 8 Go of memory and I have a hard disk with 500 Go of capcity.
For me, it is not a problem of capacity of my computer.
Up to now, I have done the exercices on both the volumique cube and the holed plate and I have had no problem at all.
I still do not understand what is the problem with the meshing cylinder.


PS : Here attached my corrected files hdf, med, comm and mess.
Attachments:

Please Log in or Create an account to join the conversation.

More
11 years 6 months ago - 11 years 6 months ago #6451 by kwou
Replied by kwou on topic Re: lecture du maillage dans Code-Aster
Hoi Albert

you still need to update the name of your element/surface group in read_maillage, see error message below:

CR d'execution de JDC en MIXTE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! <S> Exception utilisateur levee mais pas interceptee. !
! Les bases sont fermees. !
! Type de l'exception : error !
! !
! le GROUP_MA Surfin ne fait pas partie du maillage : meshmod !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fin CR d'execution de JDC en MIXTE


Info:
MAILLAGE meshinit - IMPRESSIONS NIVEAU 2

VE-28-SEPT-2012 10:39:13
NOMBRE DE NOEUDS 3080
NOMBRE DE MAILLES 1105
SEG3 105
QUAD8 1000
NOMBRE DE GROUPES DE MAILLES 3
MCbot 40
MCtop 40
MSurfin 1000



shellch=AFFE_CARA_ELEM(MODELE=modelc,
COQUE=_F(GROUP_MA='Surfin', <<<<<< change to MSurfin
EPAIS=0.10,
ANGL_REP=(0.10,0.3,),
COQUE_NCOU=1,
EXCENTREMENT=0.000,
INER_ROTA='OUI',),);

Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)

--
kind regards - kees
Last edit: 11 years 6 months ago by kwou.

Please Log in or Create an account to join the conversation.

  • bevilacqua
  • Topic Author
  • Offline
  • New Member
  • New Member
More
11 years 6 months ago #6452 by bevilacqua
Replied by bevilacqua on topic Re: lecture du maillage dans Code-Aster
Hi Kees,

Excuse me for this forgetting.

But, the problem in my version on CAELinux Salome remains and I can't find the solution.
I have changed my model with a construction like you have done.
I performed the example and the problem still remains.


PS : Here attached my corrected files hdf, med, comm and mess.
Attachments:

Please Log in or Create an account to join the conversation.

Moderators: catux
Time to create page: 0.137 seconds
Powered by Kunena Forum