DEBUT();

DSDE__=200.

SIGY__ = 1.

YOUN__=8.0e+04


# DSDE__=DEFI_CONSTANTE(VALE=DSDE__,);
# 
# SIGY_=DEFI_CONSTANTE(VALE=SIGY__,);
# 
# YOUN__=DEFI_CONSTANTE(VALE=YOUN__,);


NU=DEFI_CONSTANTE(VALE=0.38,);

MA=LIRE_MAILLAGE();

MA=DEFI_GROUP(reuse =MA,
              MAILLAGE=MA,
              CREA_GROUP_NO=_F(TOUT_GROUP_MA='OUI',),);


ACIER=DEFI_MATERIAU(ECRO_LINE=_F(D_SIGM_EPSI=DSDE__,
                                 SY=SIGY__,),
                    ELAS=_F(NU=0.3,
                            E=YOUN__,),);

# ACIER=DEFI_MATERIAU(ELAS_FO=_F(E=YOUN__,
#                                NU=NU,),
#                     ECRO_LINE_FO=_F(D_SIGM_EPSI=DSDE__,
#                                     SY=SIGY_,),);

CHMAT=AFFE_MATERIAU(AFFE=_F(MATER=ACIER,
                            TOUT='OUI',),
                    MAILLAGE=MA,);

MO=AFFE_MODELE(AFFE=_F(TOUT='OUI',
                       PHENOMENE='MECANIQUE',
                       MODELISATION='AXIS',),
               MAILLAGE=MA,);

TRACTION=AFFE_CHAR_MECA(DDL_IMPO=(_F(DY=5.E-3,
                                     GROUP_NO='CD',),
                                  _F(DY=0.0,
                                     GROUP_NO='AB',),),
                        MODELE=MO,);


RAMPE=DEFI_FONCTION(NOM_PARA='INST',
                    VALE=(0.0,0.0,1.0,1.0),);

INSTANTS=DEFI_LIST_REEL(INTERVALLE=_F(JUSQU_A=1.0,
                                      NOMBRE=10,),
                        DEBUT=0.0,);

EVOL=STAT_NON_LINE(CHAM_MATER=CHMAT,
                   MODELE=MO,
                   ARCHIVAGE=_F(LIST_INST=INSTANTS,
                                ARCH_ETAT_INIT='OUI',),
                   CONVERGENCE=_F(ITER_GLOB_MAXI=10,
                                  RESI_GLOB_RELA=1.E-05,),
                   COMP_INCR=_F(RELATION='VMIS_ISOT_LINE',),
                   INCREMENT=_F(LIST_INST=INSTANTS,
                                ),
                   NEWTON=_F(REAC_ITER=1,
                             REAC_INCR=1,),
                   EXCIT=_F(CHARGE=TRACTION,
                            FONC_MULT=RAMPE,),
                   );

EVOL=CALC_ELEM(reuse =EVOL,
               RESULTAT=EVOL,
               CHAM_MATER=CHMAT,
               MODELE=MO,
               OPTION=('SIEF_ELNO_ELGA','VARI_ELNO_ELGA'),);

EVOL=CALC_ELEM(reuse =EVOL,
               RESULTAT=EVOL,
               CHAM_MATER=CHMAT,
               MODELE=MO,
               OPTION=('SIEF_ELNO_ELGA','VARI_ELNO_ELGA'),
               );

REPONSE1=POST_RELEVE_T(ACTION=_F(OPERATION='EXTRACTION',
                       INTITULE = 'SIGYY',
                       RESULTAT = EVOL,
                       NOM_CHAM = 'SIEF_ELNO_ELGA',
                       NOM_CMP  = 'SIYY',
                       GROUP_NO = 'A',),);

REPONSE2=POST_RELEVE_T(ACTION=_F(OPERATION='EXTRACTION',
                       INTITULE = 'V1',
                       RESULTAT = EVOL,
                       NOM_CHAM = 'VARI_ELNO_ELGA',
                       NOM_CMP  = 'V1',
                       GROUP_NO = 'A',),);

FIN();
