[SOLVED] hey need help with a 2D beam problm .comm
- Mihir
- Topic Author
- Offline
- Junior Member
Less
More
- Posts: 29
- Thank you received: 0
11 years 7 months ago - 11 years 6 months ago #6849
by Mihir
[SOLVED] hey need help with a 2D beam problm .comm was created 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();
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.
- Claus
- Offline
- Moderator
Less
More
- Posts: 670
- Thank you received: 34
11 years 7 months ago #6850
by Claus
Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
Replied by Claus on topic Re: hey need help with a 2D beam problm .comm file
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
/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
Less
More
- Posts: 29
- Thank you received: 0
11 years 7 months ago #6851
by Mihir
Replied by Mihir on topic Re: hey need help with a 2D beam problm .comm file
I ran the Aster in the Salome-Meca 2011.2 module, cant figure out what I'm missing tho , been doing a tutorial
- Claus
- Offline
- Moderator
Less
More
- Posts: 670
- Thank you received: 34
11 years 7 months ago #6852
by Claus
Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
Replied by Claus on topic Re: hey need help with a 2D beam problm .comm file
Please post your .comm file - don't copy/paste it, but upload the file
/C
/C
Code_Aster release : STA11.4 on OpenSUSE 12.3 64 bits - EDF/Intel version
- Mihir
- Topic Author
- Offline
- Junior Member
Less
More
- Posts: 29
- Thank you received: 0
11 years 7 months ago - 11 years 7 months ago #6853
by Mihir
Replied by Mihir on topic Re: hey need help with a 2D beam problm .comm file
This is the command file...
Last edit: 11 years 7 months ago by Mihir.
- kwou
- Offline
- Moderator
11 years 7 months ago - 11 years 7 months ago #6854
by kwou
Interest: structural mechanics, solar energy (picture at 'my location' shows too little pv panels)
--
kind regards - kees
Replied by kwou on topic Re: hey need help with a 2D beam problm .comm file
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'
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