# PATRON 
# ANALYSE STATIQUE
# 3D

DEBUT();
# MISE EN DONNEES

MAIL=LIRE_MAILLAGE();

MAIL=DEFI_GROUP(reuse =MAIL,
                MAILLAGE=MAIL,
                CREA_GROUP_MA=_F(NOM='TOUT',
                                 TOUT='OUI',),);

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

ACIER=DEFI_MATERIAU(ELAS=_F(E=2.E11,
                            NU=0.3,),);

CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
                    AFFE=_F(TOUT='OUI',
                            MATER=ACIER,),);
# CONDITIONS AUX LIMITES ET CHARGEMENT
#POINT A : ENCASTREMENT

BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODMECA,
                        DDL_IMPO=_F(GROUP_MA='ENCAST',
                                    LIAISON='ENCASTRE',),);
#POINT B : FORCE IMPOSEE

CHARGT=AFFE_CHAR_MECA(MODELE=MODMECA,
                      FORCE_NODALE=_F(GROUP_NO='PB',
                                      FY=1000.0,),);
# RESOLUTION DU PROBLEME ELASTIQUE

RESU1=MECA_STATIQUE(MODELE=MODMECA,
                    CHAM_MATER=CHMAT,
                    EXCIT=(_F(CHARGE=BLOQUAGE,),
                           _F(CHARGE=CHARGT,),),);
# CALCUL DES EFFORTS

RESU1=CALC_ELEM(reuse =RESU1,
                RESULTAT=RESU1,
                OPTION='EFGE_ELNO_DEPL',);
# RESULTATS FORMAT GMSH

IMPR_RESU(MODELE=MODMECA,
          FORMAT='GMSH',
          RESU=_F(MAILLAGE=MAIL,
                  RESULTAT=RESU1,
                  NOM_CHAM=('DEPL','EFGE_ELNO_DEPL',),),);

FIN();