pressure over COQUE_3D

The option:

OPTION='TRIA6_7',),),);

Does not work with Aster 9.04, so if you comment that out, you should be fine. I think Kees mentions it on the page from where you got the .comm template.

/C

Hoi Claus, Borja

I hope that by now I updated most of the shell (coque_3d) examples for use with Code Aster 10+. The command CREA_MAILLAGE should be adjusted according to the mesh: only TRIA6, only QUAD8 or a mixture of both. See extract below:


=='''Note on Code Aster 10'''==
Code Aster 10+ is stricter then CA9- with respect to the MODI_MAILLAGE keyword in CREA_MAILLAGE. In case no TRIA6 or QUAD8 elements are available in the mesh, change near line xyz in the command file <command>.comm:

[code:1]
change
COQmesh=CREA_MAILLAGE(MAILLAGE=meshinit,
MODI_MAILLE=(_F(TOUT='OUI',OPTION='QUAD8_9',),
_F(TOUT='OUI',OPTION='TRIA6_7',),),);

to, for QUAD8 mesh only:

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

or, for TRIA6 mesh only:

COQmesh=CREA_MAILLAGE(MAILLAGE=meshinit,
MODI_MAILLE=(_F(TOUT='OUI',OPTION='TRIA6_7',),),);
[/code:1]

But if I missed some examples, please mention them here.

Though I am not sure whether this your problem now Borja, since the calculation runs until the end.

Added thursday 4 march.
Borja, replace the command MECA_STATIQUE() to FIN() by this part and you obtain the requested stresses.
In the CALC_ELEM() you need to add the layer where you want to extract the stresses eg NIVE_COUCHE='SUP'. ('SUP' can be either 'SUP','MOY' and 'INF' [outer, centre and inner layer]). You can also add resuINF, resuMOY if you want.
This is where I had a glance:
www.caelinux.org/wiki/index.php/Contrib:..._stresses_-_commands

[code:1]
resu=MECA_STATIQUE(MODELE=modmod,
CHAM_MATER=material,
CARA_ELEM=shellch,
EXCIT=_F(CHARGE=load_bc,),
OPTION='SIEF_ELGA_DEPL',);


resu=CALC_ELEM(reuse=resu,
MODELE=modmod,
CHAM_MATER=material,
RESULTAT=resu,
REPE_COQUE=_F(NIVE_COUCHE='SUP',),#this determines layer of EQUIvalent stresses 'EQUI_ELNO_SIGM'
OPTION=('SIGM_ELNO_DEPL',),); #includes data for 'DEPL' and 'SIGM_NOEU_DEPL'
#EXCIT=_F(CHARGE=bc_force,),);

resuSUP=CALC_ELEM(MODELE=modmod,
CHAM_MATER=material,
RESULTAT=resu,
REPE_COQUE=_F(NIVE_COUCHE='SUP',),
OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),);
#EXCIT=_F(CHARGE=bc_force,),);

resuSUP=CALC_NO(reuse=resuSUP,
RESULTAT=resuSUP,
OPTION=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM'),);

IMPR_RESU(UNITE=80,
FORMAT='MED',
MODELE=modmod,
RESU=(_F(MAILLAGE=meshmod,
RESULTAT=resu,
NOM_CHAM=('DEPL',),),
_F(MAILLAGE=meshmod,
RESULTAT=resuSUP,
NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM',),),),);

FIN();
[/code:1]

In the med file the displacements DEPL, the equivalents stresses EQUI_NOEU_SIGM at layer SUP and the local stresses SIGM_NOEU_DEPL at layer SUP are present in the corresponding fields.

kind regards - kees<br /><br />Post edited by: Kees Wouters, at: 2010/03/04 15:35