
DEBUT();

#MA=DEFI_MATERIAU(ELAS=_F(E=1.0,
#                         NU=0.4,),);

MA=DEFI_MATERIAU(ELAS_ORTH=_F(E_L=0.04,        #Ex
                              E_T=0.04,        #Ey
                              E_N=1.0,         #Ez
                              NU_LT=0.016,     #nu_xy 
                              NU_LN=0.016,     #nu_xz
                              NU_TN=0.016,     #nu_yz
                              G_LT=0.0667,     #Gxy
                              G_LN=0.0667,     #Gxz
                              G_TN=0.0667,),); #Gyz

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

#MAIL=MODI_MAILLAGE(reuse=MAIL,
#                   MAILLAGE=MAIL,
#                   ORIE_PEAU_3D=_F(GROUP_MA=('Face_6',),),
#                   );

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

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

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
                    DDL_IMPO=(_F(GROUP_MA='DX',
                                 DX=0.0,),
                              _F(GROUP_MA='DY',
                                 DY=0.0,),
                              _F(GROUP_MA='DZ',
                                 DZ=0.0,),
                              ),
                    FORCE_FACE=(
                                _F(GROUP_MA='FZ',
                                FZ=-0.1,),
                                #_F(GROUP_MA='FX',
                                #FX=-0.004,),
                                #_F(GROUP_MA='FY',
                                #FY=-0.004,),
                                ),
                    );

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

RESU=CALC_CHAMP(reuse =RESU,
                RESULTAT=RESU,
                CONTRAINTE=('SIGM_ELNO','SIGM_NOEU',),
                DEFORMATION='EPSI_NOEU',
                CRITERES=('SIEQ_ELNO','SIEQ_NOEU',),);

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

FIN();
