Code_Aster forums

jeanpierreaubry · 2012-08-02 16:31:06

hello

in 11.2, when i do the following

Code:

stat=CALC_CHAMP(reuse =stat,RESULTAT=stat,
  CONTRAINTE=(
    'SIEF_ELNO','SIPO_ELNO','SIPM_ELNO',
    'EFGE_ELNO',
  ),
  CRITERES=('SIEQ_ELNO',),
  FORCE=('REAC_NODA'),
);

statsup=POST_CHAMP(
  GROUP_MA=('panel',), RESULTAT=stat,
  EXTR_COQUE=_F(
    NUME_COUCHE=1,
    NIVE_COUCHE='SUP',
    NOM_CHAM='SIEQ_ELNO',
  ),
);

'panel' being a group of DKT
i get the following

Code:

   !------------------------------------------------------------------------!
   ! <EXCEPTION> <CALCULEL2_15>                                             !
   !                                                                        !
   ! Erreur utilisateur dans la commande POST_CHAMP :                       !
   !  On demande l'extraction pour des champs n'ayant pas de "sous-points". !
   !------------------------------------------------------------------------!

what i want is the VMIS criteria at the top face of the elements

looks like i am missing something!
but i cannot figure out what!

jean pierre aubry

Last edited by jeanpierreaubry (2012-08-02 17:42:34)

jeanpierreaubry · 2012-08-02 17:42:03

i am answering myself,
it is perfectly explained here,
http://www.code-aster.org/forum2/viewtopic.php?id=17451
thanks to the guys invaolved

once again, one never read the documentation, and the forum, enough!

pilafa · 2012-08-06 14:03:54

Hi Jean-Pierre,

Since you're working with COQUE_3D, maybe you can help me with an unresolved post I did recently.

I am not able to project the results on the original 8 nodes mesh (with PROJ_CHAMP).

It is working for the displacements, but not for the stress (I have only element results for now)

I obtain this message:

Code:

   !------------------------------------------------------------------------------------------------------------!
   ! <A> <CALCULEL_54>                                                                                          !
   !                                                                                                            !
   ! Le paramètre:  PCONTRR  de l'option:  SIGM_ELNO  n'est pas connu des TYPE_ELEM du LIGREL:  modellin.MODELE !
   !  Champ : salomres.018.000000                                                                               !
   !                                                                                                            !
   !                                                                                                            !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette                                             !
   ! alarme, vous pouvez obtenir des résultats inattendus !                                                     !
   !------------------------------------------------------------------------------------------------------------!
   
   
   !----------------------------------------------------------------------!
   ! <A> <CALCULEL4_63>                                                   !
   !                                                                      !
   ! On n'a pas pu projeter le champ SIGM_ELNO de la sd_resultat Solution !
   !  vers la sd_resultat salomres pour le numéro d'ordre 1               !
   !                                                                      !
   !                                                                      !
   ! Ceci est une alarme. Si vous ne comprenez pas le sens de cette       !
   ! alarme, vous pouvez obtenir des résultats inattendus !               !
   !----------------------------------------------------------------------!

here's the comm:

Code:

#
#Pour VERSION  10
#

DEBUT();

#
# PREPROC
#

#quad8 mesh
mesh_1=LIRE_MAILLAGE(FORMAT='MED',VERI_MAIL=_F(VERIF='OUI'));



#Dummy #QUAD8 model to postproc
        # Will produce warnings: 
        #<CALCULEL2_64> 
        #<CALCULEL2_63> 
modellin=AFFE_MODELE(MAILLAGE=mesh_1,
                 AFFE=(_F(GROUP_MA='plate',
                          PHENOMENE='MECANIQUE',
                          MODELISATION='3D',),),);   


mesh=CREA_MAILLAGE(MAILLAGE=mesh_1,
                    MODI_MAILLE=_F(GROUP_MA='plate',OPTION='QUAD8_9',),);


mesh=MODI_MAILLAGE(reuse =mesh,
                       MAILLAGE=mesh,
                       ORIE_PEAU_2D=_F(GROUP_MA=('plate',),),                     
                       ); 


#Create model for 3D shell structure
model=AFFE_MODELE(MAILLAGE=mesh,
                 INFO=2,
                 AFFE=(_F(GROUP_MA='plate',
                          PHENOMENE='MECANIQUE',
                         MODELISATION='COQUE_3D',
                          ),),);


steel=DEFI_MATERIAU(ELAS=_F(E=210.0e9,
                            NU=0.3,
                            RHO=7800.0),);


material=AFFE_MATERIAU(MAILLAGE=mesh,
                      AFFE=_F(TOUT='OUI',MATER=steel,),);


realcst=AFFE_CARA_ELEM(
                   MODELE=model,
                   COQUE=(
                         _F(GROUP_MA='plate',EPAIS=0.02,VECTEUR=(1.0,0.0,0.0)),
                        ),
                        );


bchold=AFFE_CHAR_MECA(MODELE=model,
                       DDL_IMPO=(
                                 _F(GROUP_NO=("h1"),DX=0.0,DY=0.0,DZ=0.0,DRX=0.0,DRY=0.0,DRZ=0.0,),  #FIXED
                                 _F(GROUP_NO=("h2"),DX=0.0,DY=0.0,DZ=0.0,DRX=0.0,DRY=0.0,DRZ=0.0,),
                                 _F(GROUP_NO=("s1"),DRX=0.0,),  #symetry
                                 _F(GROUP_NO=("s2"),DRY=0.0,),       
                           #      _F(GROUP_NO=("nmid"),DRZ=0.0),  #rigid motion                                                                                                                        
                                 ),
                        PRES_REP=(
                                  _F(GROUP_MA=("plate"),PRES=-0.5E4,), 
                                  ),
                                  );

#
# SOLUTION
#


Solution=MECA_STATIQUE(MODELE=model,
                   CHAM_MATER=material,
                   CARA_ELEM=realcst,
                   EXCIT=_F(CHARGE=bchold,),);    

#
# POSTPTOC
#

Solution=CALC_ELEM(reuse = Solution,
                 RESULTAT=Solution,
                 INFO=2,
                 REPE_COQUE=_F(NIVE_COUCHE='INF',
                               VECTEUR=(1.0,0.0,0.0,),),                 
                 TYPE_OPTION='TOUTES',
                 OPTION=('SIGM_ELNO','SIEQ_ELNO',),);

     

#
# Project and print
#

# determine / project the result fields (champ) from modified model modmod to the initial model modinit
salomres=PROJ_CHAMP(RESULTAT=Solution,
                   MODELE_1=model,        #project modes/results of this model
                   MODELE_2=modellin,);   #...to this one 



IMPR_RESU(FORMAT='MED',
           RESU=(_F(
                    NOM_CHAM=('SIGM_ELNO','SIEQ_ELNO','DEPL',),
                 #   RESULTAT=salomres,          
                    RESULTAT=Solution,                                                  
                    MAILLAGE=mesh,),   
                    ),
                    );  

FIN();

Thanks
Pierre

jeanpierreaubry · 2012-08-06 14:52:24

hello

it looks like you are trying too fast to calculate the stress in the bottom face of the plates
you have to calculate in the middle plane as a first step
and then calculate another concept for the top plate, second step
and the eventually a third concept for Von Mises and the like
this is explained in
http://www.code-aster.org/forum2/viewtopic.php?id=17451

another point it seems better to use CALC_CHAMP instead of CALC_ELEM in 11.* version

jean pierre aubry

pilafa · 2012-08-06 15:19:44

Hi Jean-Pierre.

Some comments :

- Let say that I don't want the Von Mises stress for now;

- I thought that the complicated procedure explained in the post you are referring to (I actually initiated this post) was "temporary" because version 11.2 is still in development. I switched to version 10 for this reason.

- I can calculate the correct value of stress with SIGM_ELNO (for all three position (inf,mid,sup)).
It is just a problem of projecting it back to the 8 nodes elements.
I need this because paraVIS cannot postptoc 9 nodes correctly.

Thanks
Pierre

jeanpierreaubry · 2012-08-06 15:45:27

i wonder whether
the ORIE_PEAU command should not be before the quad mesh creation

and, more important,
in IMPR_RESU in MED format you should have
MAILLAGE=mesh_1, which is the original mesh
instead of MAILLAGE=mesh, which is the modified quadratic mesh
to project the results on the original mesh

pilafa · 2012-08-08 14:15:15

Thanks for the suggestions.
However, putting ORI_PEAU before doesn't change anything.

You are certainly right about mesh_1, but the error message is emitted just after the PROJ_CHAMP command, so correcting this doesn't solve the problem either.

Strange...
Thanks anyway!
PL

Thomas DE SOZA · 2012-08-10 00:06:58

pilafa wrote:
- I thought that the complicated procedure explained in the post you are referring to (I actually initiated this post) was "temporary" because version 11.2 is still in development. I switched to version 10 for this reason.

- I can calculate the correct value of stress with SIGM_ELNO (for all three position (inf,mid,sup)).
It is just a problem of projecting it back to the 8 nodes elements.
I need this because paraVIS cannot postptoc 9 nodes correctly.

Note that the problem you have drawn our attention to in this other post (http://www.code-aster.org/forum2/viewtopic.php?id=17451) has been fixed in the unstable branch (11.2.x). So if you can afford the upgrade you should be all set.

TdS

pilafa · 2012-08-10 11:13:38

Thanks for the info Thomas. I will try it.
For now I am playing around with version 10.
Pierre

Code_Aster forums

#1 2012-08-02 16:31:06

[SOLVED]calculating VMIS in DKT in 11.2

Code:

Code:

#2 2012-08-02 17:42:03

Re: [SOLVED]calculating VMIS in DKT in 11.2

#3 2012-08-06 14:03:54

Re: [SOLVED]calculating VMIS in DKT in 11.2

Code:

Code:

#4 2012-08-06 14:52:24

Re: [SOLVED]calculating VMIS in DKT in 11.2

#5 2012-08-06 15:19:44

Re: [SOLVED]calculating VMIS in DKT in 11.2

#6 2012-08-06 15:45:27

Re: [SOLVED]calculating VMIS in DKT in 11.2

#7 2012-08-08 14:15:15

Re: [SOLVED]calculating VMIS in DKT in 11.2

#8 2012-08-10 00:06:58

Re: [SOLVED]calculating VMIS in DKT in 11.2

pilafa wrote:

#9 2012-08-10 11:13:38

Re: [SOLVED]calculating VMIS in DKT in 11.2

Board footer