DEBUT();

MA=DEFI_MATERIAU(ELAS_ORTH=_F(E_L=3.0776e2,
                               E_T=3.0776e2,
                               E_N=7.694e3,
                               G_LT=5.12933e2,
                               G_TN=5.12933e2,
                               G_LN=5.12933e2,
                               NU_LT=0.016,
                               NU_TN=0.016,
                               NU_LN=0.016,),);

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


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='kotei',
                           DX=0.0,
                           DY=0.0,
                           DZ=0.0,),
                    FORCE_NODALE=
                        _F(GROUP_NO='saika_n',
                           FY=-8.26446e-7,),
                    );

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();
