×

Notice

The forum is in read only mode.

[SOLVED] hey need help with a 2D beam problm .comm

  • Mihir
  • Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
11 years 7 months ago - 11 years 6 months ago #6849 by Mihir
Hey! as I'm new here, i was trying to get started with a simple 2D beam fixed at both ends and having a uniformly distributed load. So I created the geometry and mesh using Salome-Meca , and then generated the .comm file using Eficas . I noticed in the .comm file that the AFFE_CARA_ELEM icon is in yellow, while all the rest are in green (i get a feeling that something's not right here-so I've attached the source file of the .comm file at the end).

When i try to run the Aster file later, it does not do the computation; it comes up with a message the command file is not in use.

Kindly tell me where I'm going wrong as I've been going through tutorials and web content for the past 2 days, and I really cant proceed any further without help.

Thanking in advance.




DEBUT();

MA=DEFI_MATERIAU(ELAS=_F(E=210000000000.0,
NU=0.3,
RHO=8e-9,),);

MAIL=LIRE_MAILLAGE(FORMAT='MED',);

MAIL=MODI_MAILLAGE(reuse =MAIL,
MAILLAGE=MAIL,
ORIE_PEAU_2D=_F(GROUP_MA='Group_1',),);

MODE=AFFE_MODELE(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='C_PLAN',),);

SD_4=AFFE_CARA_ELEM(MODELE=MODE,
POUTRE=_F(GROUP_MA='Group_1',
SECTION='RECTANGLE',
VARI_SECT='CONSTANT',
CARA=('HY','HZ',),
VALE=(0.04,0.02,),),);

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
MATER=MA,),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
PESANTEUR=_F(GROUP_MA='Group_1',
GRAVITE=10,
DIRECTION=(0,0,-1,),),
DDL_IMPO=_F(GROUP_MA='Group_1',
DX=0.0,
DY=0.0,),
PRES_REP=_F(GROUP_MA='Group_1',
PRES=0.0,),
FORCE_POUTRE=_F(GROUP_MA='Group_1',
FZ=-100,),);

RESU=MECA_STATIQUE(MODELE=MODE,
CHAM_MATER=MATE,
EXCIT=_F(CHARGE=CHAR,),);

RESU=CALC_ELEM(reuse =RESU,
MODELE=MODE,
CHAM_MATER=MATE,
RESULTAT=RESU,
OPTION=('SIGM_ELNO','SIEQ_ELNO',),
EXCIT=_F(CHARGE=CHAR,),);

RESU=CALC_NO(reuse =RESU,
RESULTAT=RESU,
OPTION=('SIGM_NOEU','SIEQ_NOEU',),);

IMPR_RESU(FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=MAIL,
RESULTAT=RESU,
NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),),);

FIN();
Last edit: 11 years 6 months ago by Mihir.
More
11 years 7 months ago #6850 by Claus
What version of Code_Aster are you using? A yellow icon indicates that you're missing a field in the CALC_ELEM group. I must admit that I can't remember the requirements in my head right now, but compare to a working .comm file.

/C

Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
The following user(s) said Thank You: Mihir
  • Mihir
  • Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
11 years 7 months ago #6851 by Mihir
I ran the Aster in the Salome-Meca 2011.2 module, cant figure out what I'm missing tho :(, been doing a tutorial
More
11 years 7 months ago #6852 by Claus
Please post your .comm file - don't copy/paste it, but upload the file :)

/C

Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
  • Mihir
  • Topic Author
  • Offline
  • Junior Member
  • Junior Member
More
11 years 7 months ago - 11 years 7 months ago #6853 by Mihir
This is the command file...
Attachments:
Last edit: 11 years 7 months ago by Mihir.
More
11 years 7 months ago - 11 years 7 months ago #6854 by kwou
Hoi D_Riel

I havenot used beams in 2D, but you may need to add a wall thickness.
In 3D I use:

Beam1=AFFE_CARA_ELEM(MODELE=BeamMod,
POUTRE=(_F(GROUP_MA='beam1',
SECTION='RECTANGLE',
CARA=('HY','HZ','EP',),
VALE=(0.0200,0.0100,0.0025,),),),
ORIENTATION=(_F(GROUP_MA='beam1',
CARA='VECT_Y',
VALE=(0.0,0.0,1.0,),),),)

u4.42.01.pdf says:

|'H'|'EP',
or
|‘HY’|'HZ'|'EPY'|'EPZ'

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

--
kind regards - kees
Last edit: 11 years 7 months ago by kwou.
The following user(s) said Thank You: Mihir
Moderators: catux
Time to create page: 0.133 seconds
Powered by Kunena Forum