DEBUT();

#Read MED mesh file


Mc1=LIRE_MAILLAGE(UNITE=20,
                  FORMAT='MED',
                  NOM_MED='Mcyl',
                  INFO_MED=3,
                  INFO=1,);

Mb1=LIRE_MAILLAGE(UNITE=21,
                  FORMAT='MED',
                  NOM_MED='Mblock',
                  INFO_MED=3,
                  INFO=1,);

Massy=ASSE_MAILLAGE(MAILLAGE_1=Mb1,
                    MAILLAGE_2=Mc1,
                    OPERATION='SUPERPOSE',);

#Assign to the Model
Blmod=AFFE_MODELE(MAILLAGE=Massy,
                  AFFE=_F(TOUT='OUI',
                          PHENOMENE='MECANIQUE',
                          MODELISATION='3D',),);

#Define material properties

steel=DEFI_MATERIAU(ELAS=_F(E=2.1e5,
                            NU=.28,
                            RHO=0.00785,),);

#Assign material to mesh

Blmat=AFFE_MATERIAU(MAILLAGE=Massy,
                    AFFE=_F(TOUT='OUI',
                            MATER=steel,),);

# define boundary conditions and loads
# boundary conditions - load

load=AFFE_CHAR_MECA(MODELE=Blmod,
                    DDL_IMPO=(_F(GROUP_MA='BPbot',
                                 DX=0.0,
                                 DY=0.0,
                                 DZ=0.0,),
                              _F(GROUP_MA='CPtop',
                                 DX=0.1,
                                 DY=0.0,
                                 DZ=-0.01,),),
                    LIAISON_MAIL=_F(GROUP_MA_MAIT='Block',
                                    GROUP_MA_ESCL='CPbot',
                                    TYPE_RACCORD='MASSIF',
                                    #DDL_MAIT='DNOR',
                                    #DDL_ESCL='DNOR',
                                    CENTRE=0.0,),);

#Execute analysis

result=MECA_STATIQUE(MODELE=Blmod,
                     CHAM_MATER=Blmat,
                     EXCIT=_F(CHARGE=load,),);

result=CALC_ELEM(reuse=result,
                  MODELE=Blmod,
                  RESULTAT=result,
                  TOUT='OUI',
                  OPTION=('SIGM_ELNO','SIEQ_ELNO'),);
                  #OPTION=('SIGM_ELNO_DEPL','EQUI_ELNO_SIGM','SIEF_ELNO_ELGA',),);

result=CALC_NO(reuse =result,
               RESULTAT=result,
               OPTION=('SIEQ_NOEU','SIGM_NOEU'),);
#               OPTION=('SIGM_NOEU_DEPL','EQUI_NOEU_SIGM',),);
#               OPTION=('SIGM_NOEU_DEPL','SIGM_NOEU_DEPL','EQUI_NOEU_SIGM',),);

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=_F(MAILLAGE=Massy,
                  RESULTAT=result,
                  NOM_CHAM=('SIGM_ELNO','SIGM_NOEU','SIEQ_ELNO','DEPL'),),);
                  #NOM_CHAM=('SIEQ_NOEU','DEPL',),),);

FIN();