DEBUT(LANG='FR')


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

mesh1 = MODI_MAILLAGE(identifier='1:1',
                      reuse=mesh1,
                      MAILLAGE=mesh1,
                      ORIE_NORM_COQUE=_F(GROUP_MA=('web', 'flange1', 'flange2')))

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

elemprop = AFFE_CARA_ELEM(identifier='3:1',
                          COQUE=(_F(COQUE_NCOU=1,
                                    EPAIS=12.0,
                                    GROUP_MA=('web', ),
                                    INER_ROTA='OUI',
                                    VECTEUR=(0.0, 0.0, 1.0)),
                                 _F(COQUE_NCOU=1,
                                    EPAIS=21.0,
                                    GROUP_MA=('flange1', ),
                                    INER_ROTA='OUI'),
                                 _F(COQUE_NCOU=1,
                                    EPAIS=21.0,
                                    GROUP_MA=('flange2', ),
                                    INER_ROTA='OUI')),
                          MODELE=model)

mater = DEFI_MATERIAU(identifier='4:1',
                      ELAS=_F(E=206000.0,
                              NU=0.3))

fieldmat = AFFE_MATERIAU(identifier='5:1',
                         AFFE=_F(MATER=(mater, ),
                                 TOUT='OUI'),
                         MAILLAGE=mesh1,
                         MODELE=model)

load = AFFE_CHAR_MECA(identifier='6:1',
                      DDL_IMPO=(_F(DX=0.0,
                                   DY=0.0,
                                   GROUP_MA=('kotei1', 'kotei2')),
                                _F(DX=0.0,
                                   DY=0.0,
                                   DZ=0.0,
                                   GROUP_MA=('kasetukotei1', 'kasetukotei2'))),
                      FORCE_COQUE=_F(FZ=-1.0,
                                     GROUP_MA=('flange1', )),
                      LIAISON_MAIL=(_F(GROUP_MA_ESCL=('web1', ),
                                       GROUP_MA_MAIT=('flange1', ),
                                       TYPE_RACCORD='COQUE'),
                                    _F(GROUP_MA_ESCL=('web2', ),
                                       GROUP_MA_MAIT=('flange2', ),
                                       TYPE_RACCORD='COQUE')),
                      LIAISON_UNIF=(_F(DDL=('DRX', 'DRY', 'DRZ'),
                                       GROUP_MA=('web1', )),
                                    _F(DDL=('DRX', 'DRY', 'DRZ'),
                                       GROUP_MA=('web2', ))),
                      MODELE=model)

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

unnamed1 = CREA_CHAMP(identifier='8: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='9:1',
           CARA_ELEM=elemprop,
           CHAM_MATER=fieldmat,
           CHARGE=(load, ),
           MATR_ASSE=(_F(MATRICE=CO('SAmech'),
                         OPTION='RIGI_MECA'),
                      _F(MATRICE=CO('SAgoeom'),
                         OPTION='RIGI_GEOM',
                         SIEF_ELGA=unnamed1)),
           MODELE=model,
           NUME_DDL=CO('nddl'))

IMPR_RESU(identifier='10:1',
          FORMAT='MED',
          RESU=_F(CARA_ELEM=elemprop,
                  NOM_CHAM=('DEPL', ),
                  RESULTAT=reslin),
          UNITE=80)

table = INFO_MODE(identifier='13:1',
                  CHAR_CRIT=(-5.0, 5.0),
                  MATR_RIGI=SAmech,
                  MATR_RIGI_GEOM=SAgoeom,
                  TYPE_MODE='MODE_FLAMB')

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

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

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

unnamed4 = CALC_CHAMP(identifier='17:1',
                      CARA_ELEM=elemprop,
                      CHAM_MATER=fieldmat,
                      CONTRAINTE=('SIGM_NOEU', ),
                      MODELE=model,
                      RESULTAT=unnamed3)

IMPR_RESU(identifier='18:1',
          FORMAT='MED',
          RESU=(_F(CARA_ELEM=elemprop,
                   NOM_CHAM=('DEPL', ),
                   RESULTAT=modes),
                _F(CARA_ELEM=elemprop,
                   NOM_CHAM=('SIGM_NOEU', ),
                   RESULTAT=unnamed4)),
          UNITE=2)

FIN()