DEBUT();

empty=DEFI_MATERIAU(ELAS=_F(E=1.0E-20,              #Empty
                          NU=0.0,),);

mat=DEFI_MATERIAU(ELAS=_F(E=29e2,              #Carbon
                          NU=0.4,),);

fib=DEFI_MATERIAU(ELAS=_F(E=28e4,              #Carbon
                          NU=0.2,),);

#mat=DEFI_MATERIAU(ELAS=_F(E=400,                #Carbon
#                          NU=0.05,),);

#fib=DEFI_MATERIAU(ELAS=_F(E=2e5,                #Steel
#                          NU=0.25,),);


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

MAIL=MODI_MAILLAGE(reuse =MAIL,
                   MAILLAGE=MAIL,
                   ORIE_PEAU_3D=_F(GROUP_MA=('ZY_Side','XY_Side',),),);

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

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
                   AFFE=(_F(GROUP_MA='Matrix',
                            MATER=mat,),
                         _F(GROUP_MA='Fibers',
                            MATER=fib,),
                         _F(GROUP_MA='Fiber',
                            MATER=empty,),
            #             _F(GROUP_MA='fs1',
            #                MATER=empty,),
            #             _F(GROUP_MA='fs2',
            #                MATER=mat,),
                        ),
                   );

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
                    DDL_IMPO=(
                              _F(GROUP_MA='YX_Side',
                                 DZ=0.0,),
                              _F(GROUP_MA='YZ_Side',
                                 DX=0.0,),
                              _F(GROUP_MA='ZX_Bott',
                                 DY=0,),
                              _F(GROUP_MA='XZ_Top',
                                 DY=0,),
                 #             ),
                 #   PRES_REP=(
                              _F(GROUP_MA='XY_Side',
                 #                PRES=1,),
                                  DZ=-1.0),
                              _F(GROUP_MA='ZY_Side',
                 #                PRES=1,),
                                  DX=1.0),
                              ),
                    );

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','EPSI_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','SIGM_ELNO',),),);

IMPR_RESU(FORMAT='RESULTAT',
          UNITE=22,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM='SIGM_NOEU',
                  FORM_TABL='OUI',
                  GROUP_MA='Matrix',
                  IMPR_COOR='OUI',),);

IMPR_RESU(FORMAT='RESULTAT',
          UNITE=23,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM='SIEQ_NOEU',
                  FORM_TABL='OUI',
                  GROUP_MA='Matrix',
                  IMPR_COOR='OUI',),);

IMPR_RESU(FORMAT='RESULTAT',
          UNITE=24,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM='SIGM_NOEU',
                  FORM_TABL='OUI',
                  GROUP_MA='Fibers',
                  IMPR_COOR='OUI',),);

IMPR_RESU(FORMAT='RESULTAT',
          UNITE=25,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM='SIEQ_NOEU',
                  FORM_TABL='OUI',
                  GROUP_MA='Fibers',
                  IMPR_COOR='OUI',),);

IMPR_RESU(FORMAT='RESULTAT',
          UNITE=26,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM='SIGM_NOEU',
                  FORM_TABL='OUI',
                  IMPR_COOR='OUI',),);

IMPR_RESU(FORMAT='RESULTAT',
          UNITE=27,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM='SIEQ_NOEU',
                  FORM_TABL='OUI',
                  IMPR_COOR='OUI',),);

IMPR_RESU(FORMAT='RESULTAT',
          UNITE=28,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM='SIEQ_NOEU',
                  FORM_TABL='OUI',
                  GROUP_MA='Fiber',
                  IMPR_COOR='OUI',),);
       
FIN();