DEBUT(LANG='FR')


mesh = LIRE_MAILLAGE(identifier='0:1',
                     FORMAT='MED',
                     UNITE=20)

model = AFFE_MODELE(identifier='1:1',
                    AFFE=_F(MODELISATION=('DKT', ),
                            PHENOMENE='MECANIQUE',
                            TOUT='OUI'),
                    MAILLAGE=mesh)

elemprop = AFFE_CARA_ELEM(identifier='2:1',
                          COQUE=_F(COQUE_NCOU=1,
                                   EPAIS=12.0,
                                   GROUP_MA=('zentai', ),
                                   INER_ROTA='OUI'),
                          MODELE=model)

kou = DEFI_MATERIAU(identifier='3:1',
                    ELAS=_F(E=206000.0,
                            NU=0.3))

fieldmat = AFFE_MATERIAU(identifier='4:1',
                         AFFE=_F(MATER=(kou, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh,
                         MODELE=model)

load = AFFE_CHAR_MECA(identifier='5:1',
                      DDL_IMPO=(_F(DX=0.0,
                                   DY=0.0,
                                   DZ=0.0,
                                   GROUP_MA=('kotei2', )),
                                _F(DX=0.0,
                                   DZ=0.0,
                                   GROUP_MA=('kotei1', ))),
                      FORCE_ARETE=_F(FY=100.0,
                                     GROUP_MA=('saika', )),
                      MODELE=model)

reslin = MECA_STATIQUE(identifier='6:1',
                       CARA_ELEM=elemprop,
                       CHAM_MATER=fieldmat,
                       EXCIT=_F(CHARGE=load),
                       MODELE=model,
                       OPTION='SIEF_ELGA',
                       SOLVEUR=_F(METHODE='MUMPS',
                                  RESI_RELA=0.006))

unnamed = CREA_CHAMP(identifier='7:1',
                     NOM_CHAM='SIEF_ELGA',
                     NUME_ORDRE=(1, ),
                     OPERATION='EXTR',
                     RESULTAT=reslin,
                     TYPE_CHAM='ELGA_SIEF_R',
                     TYPE_MAXI='MINI',
                     TYPE_RESU='VALE')

ASSEMBLAGE(identifier='8:1',
           CARA_ELEM=elemprop,
           CHAM_MATER=fieldmat,
           CHARGE=(load, ),
           MATR_ASSE=(_F(MATRICE=CO('SAmech'),
                         OPTION='RIGI_MECA'),
                      _F(MATRICE=CO('SAgeom'),
                         OPTION='RIGI_GEOM',
                         SIEF_ELGA=unnamed)),
           MODELE=model,
           NUME_DDL=CO('nddl'))

IMPR_RESU(identifier='9:1',
          FORMAT='MED',
          RESU=_F(NOM_CHAM=('DEPL', ),
                  RESULTAT=reslin),
          UNITE=4)

table = INFO_MODE(identifier='12:1',
                  CHAR_CRIT=(-5.0, 80.0),
                  MATR_RIGI=SAmech,
                  MATR_RIGI_GEOM=SAgeom,
                  TYPE_MODE='MODE_FLAMB')

modes = CALC_MODES(identifier='13:1',
                   CALC_CHAR_CRIT=_F(CHAR_CRIT=(-5.0, 80.0),
                                     TABLE_CHAR_CRIT=table),
                   MATR_RIGI=SAmech,
                   MATR_RIGI_GEOM=SAgeom,
                   OPTION='BANDE',
                   SOLVEUR_MODAL=_F(DIM_SOUS_ESPACE=80,
                                    METHODE='SORENSEN',
                                    NMAX_ITER_SOREN=80),
                   TYPE_RESU='MODE_FLAMB',
                   VERI_MODE=_F(SEUIL=1e-06,
                                STOP_ERREUR='NON'))

unnamed0 = CALC_CHAMP(identifier='14:1',
                      CARA_ELEM=elemprop,
                      CHAM_MATER=fieldmat,
                      CONTRAINTE=('SIGM_ELNO', ),
                      MODELE=model,
                      RESULTAT=modes)

unnamed2 = POST_CHAMP(identifier='15:1',
                      EXTR_COQUE=_F(NIVE_COUCHE='SUP',
                                    NOM_CHAM=('SIGM_ELNO', ),
                                    NUME_COUCHE=1),
                      RESULTAT=unnamed0)

unnamed1 = CALC_CHAMP(identifier='16:1',
                      CARA_ELEM=elemprop,
                      CHAM_MATER=fieldmat,
                      CONTRAINTE=('SIGM_NOEU', ),
                      MODELE=model,
                      RESULTAT=unnamed2)

IMPR_RESU(identifier='17:1',
          FORMAT='MED',
          RESU=(_F(CARA_ELEM=elemprop,
                   NOM_CHAM=('DEPL', ),
                   RESULTAT=modes),
                _F(CARA_ELEM=elemprop,
                   NOM_CHAM=('SIGM_NOEU', ),
                   RESULTAT=unnamed1)),
          UNITE=80)

FIN()