×

Notice

The forum is in read only mode.

Problem Getting The Code Aster Tutorial 1 Working

  • Gordon Ellison
  • Topic Author
  • Offline
  • New Member
  • New Member
More
13 years 1 month ago #6200 by Gordon Ellison
I am a very new user of CAELinux. I am trying it on a LiveUSB drive. My first attempt is using the Perforated Cube, known as Tutorial 1. I can make through everything, including solving the problem with Code Aster, but I don't get results. I am hoping someone will see something in the following .mess file. There is an exit code at the end of the file.

Thanks in advance.

Version parallèle de Code_aster compilée avec MPI
Exécution sur le processeur de nom caelinux de rang 0
Nombre de processeurs utilisés 1
==========================================
==========================================
#Linear Statics with 3D linear solid elements
#
#template by J.Cugnoni, CAELinux.com, 2005

DEBUT();

#Read MED MESH File
#
#First command LIRE_MAILLAGE is used to read the MED mesh file generated by Salome.
#
#To Do:
#Enter the name of your Salome Mesh in b_format_med->NOM_MED

MeshLin=LIRE_MAILLAGE(UNITE=20,
FORMAT='MED',
NOM_MED='Mesh_1',
INFO_MED=1,);

#Assigns a physical model to geometric entities.
#
#Here we assume that all the geometric entities (TOUT=OUI) are used for mechanical simulation (PHENOMENE=MECANIQUE) with 3D solid elements.
#
#To Do (Optional):
#you can assign other physics or element types (like shells for example) to some of the elements by replacing TOUT=OUI in AFFE with GROUP_MA = TheElementGroupYouWantToModel
#

FEMLin=AFFE_MODELE(MAILLAGE=MeshLin,
AFFE=_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),);

#Material properties
#
#To Do:
#Enter your material properties in this section
#if necessary Copy/Paste the DEFI_MATERIAU command to add a second material

Steel=DEFI_MATERIAU(ELAS=_F(E=2.1e11,
NU=0.27,
RHO=7800.0,),);

#Assign Material properties to Elements
#
#TO DO:
#If you need more than one material, you need to enter pairs of element group <-> materials by duplicating the AFFE option.

Mat=AFFE_MATERIAU(MAILLAGE=MeshLin,
MODELE=FEMLin,
AFFE=_F(TOUT='OUI',
MATER=Steel,),);

#Boundary conditions
#
#This section defines the boundary conditions of the FEA, use DDL_IMPO on selected groups to impose displacements and FORCE_* / PRESSION to apply forces / pressures to selected groups
#
#To Do:
#for each boundary conditions, you need to choose the appropriate option, for example DDL_IMPO for imposed displacements, and assign this option to a selected region of the mesh by using the GROUP_NO option for a node group or the GROUP_MA for face / volume groups.

BCnd=AFFE_CHAR_MECA(MODELE=FEMLin,
DDL_IMPO=_F(GROUP_MA='depl',
DX=0.0,
DY=0.0,
DZ=0.0,),
PRES_REP=_F(GROUP_MA='press',
PRES=1,),);

#Finite Element Solution

Solution=MECA_STATIQUE(MODELE=FEMLin,
CHAM_MATER=Mat,
EXCIT=_F(CHARGE=BCnd,),
SOLVEUR=_F(METHODE='MUMPS',
RENUM='AUTO',),);

#Compute VonMises Stress / Strain

Solution=CALC_ELEM(reuse =Solution,
RESULTAT=Solution,
OPTION=('SIEF_ELNO','SIEQ_ELNO',),);

Solution=CALC_NO(reuse =Solution,
RESULTAT=Solution,
OPTION=('REAC_NODA','SIEQ_NOEU','SIEF_NOEU',),);

#Write Results to MED file

IMPR_RESU(MODELE=FEMLin,
FORMAT='MED',
UNITE=80,
RESU=_F(RESULTAT=Solution,
TOUT_CHAM='OUI',
TOUT_CMP='OUI',),);

FIN(FORMAT_HDF='OUI',);

==========================================
==========================================
EXECUTION_CODE_ASTER_EXIT_5666-caelinux=139

<S>_ABNORMAL_ABORT exit code = 1
More
13 years 1 month ago #6201 by kwou
Hoi GordonEll

is that the complete mess file?
If not, pls post it completely.

kind regards - kees

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

--
kind regards - kees
Moderators: catux
Time to create page: 0.113 seconds
Powered by Kunena Forum