Spannungen für Top und Buttom bei 2D Elementen

Hallo,

wie kann ich mir die Top und Buttom Spannungen bei 2D Systemen ausgeben lassen? Mein jetziges Modell gibt die Buttom Spannungen aus. Ich habe aber auch Interesse an den Top Spannungen.


Gruß und Dank Markus

Hallo,

ist das der prinzipielle Weg oder gibt es eine andere Möglichkeit?

[code:1]
The stresses in the plate - shell elements

Depending on which layer you would like to have the stresses, you need to define in the command file:

* for layer NIVE_COUCHE='SUP' (superieur, top)
* for layer NIVE_COUCHE='MOY' (moyenne, center)
* for layer NIVE_COUCHE='INF' (inferieur, bottom)

e.g to define the stresses in the top layer you can use the following sequence of commands:

# define top layer by 'SUP'
top=CALC_ELEM(REPE_COQUE=_F(NUME_COUCHE=1,NIVE_COUCHE='SUP',),
OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
RESULTAT=result,);
....
top=CALC_NO(reuse=top,
RESULTAT=top,
OPTION='SIGM_NOEU_DEPL',);
...
# because Salome uses a different mesh (quad8 iso coque_3d/quad9) we have to project the stress fields to the initial Salome mesh *)
salomT=PROJ_CHAMP(RESULTAT=top,
MODELE_1=modelc, # project fields of this model to
MODELE_2=modinit,); # field of modele_2
...
# and finally print the results.
# stresses in , and layer,
# displacements in salomR
IMPR_RESU(FORMAT='MED',
UNITE=21,
RESU=(_F(MAILLAGE=meshinit,RESULTAT=salomT,
NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM','EQUI_ELNO_SIGM','DEPL','SIGM_ELNO_COQU',),),
_F(MAILLAGE=meshinit,RESULTAT=salomB,
NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM','EQUI_ELNO_SIGM','DEPL','SIGM_ELNO_COQU',),),
_F(MAILLAGE=meshinit,RESULTAT=salomC,
NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM','EQUI_ELNO_SIGM','DEPL','SIGM_ELNO_COQU',),),
_F(MAILLAGE=meshinit,RESULTAT=salomR,
NOM_CHAM=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM','EQUI_ELNO_SIGM','DEPL','SIGM_ELNO_COQU',),),),);

*) Added july 2010: In later versions op C-Aster (vs10.x) this in not necessary any more. You can directly write the results on a quadratic mesh (quad8 or tria6 elements). I have not tried linear elements.
[/code:1]


oder kann man wie folgt arbeiten?

[code:1]
Bottom=CALC_ELEM(
REPE_COQUE=_F(NUME_COUCHE=1,
NIVE_COUCHE='INF',
PLAN='MAIL'),
OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
TYPE_OPTION='TOUTES',
RESULTAT=Res,
);

Top=CALC_ELEM(
REPE_COQUE=_F(NUME_COUCHE=1,
NIVE_COUCHE='SUP',
PLAN='MAIL'),
OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM',),
TYPE_OPTION='TOUTES',
RESULTAT=Res,
);


IMPR_RESU(RESU=(_F(RESULTAT=Top,),
_F(RESULTAT=Bottom,),
_F(RESULTAT=Res,),),);
[/code:1]

Gruß und Dank Markus<br /><br />Post edited by: MGolbs, at: 2011/02/15 19:23