DEBUT();

sig400=DEFI_FONCTION(
                     NOM_PARA='EPSI',
                     NOM_RESU='SIGM',
                     VALE=(0.01,40,
                           0.02,51,
                           0.03,58,
                           0.04,65,
                           0.05,67,
                           0.06,70,
                           0.096,80,
                           ),
                     PROL_DROITE='LINEAIRE',
                     PROL_GAUCHE='LINEAIRE',);


MA=DEFI_MATERIAU(ELAS=_F(E=4000,
                         NU=0.3,),
                 TRACTION=_F(SIGM=sig400,),);

MAIL=LIRE_MAILLAGE(FORMAT='MED',);


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


MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
                   AFFE=_F(TOUT='OUI',
                           MATER=MA,),);



CHAR=AFFE_CHAR_MECA(MODELE=MODE,
                    DDL_IMPO=(
                        _F(GROUP_MA=('kotei1',),
                           DX=0.0,
                           DY=0.0,
                           DZ=0.0,),
                        _F(GROUP_MA='kotei2',
                           DX=0.0,
                           DY=0.0,
                           DZ=0.0,),
                        _F(GROUP_MA='kotei3',
                           DX=0.0,
                           DY=0.0,
                           DZ=0.0,),
                        _F(GROUP_MA='kotei4',
                           DX=0.0,
                           DY=0.0,
                           DZ=0.0,),
                        _F(GROUP_MA='kotei5',
                           DX=0.0,
                           DY=0.0,
                           DZ=0.0,),
                        _F(GROUP_MA='kotei6',
                           DX=0.0,
                           DY=0.0,
                           DZ=0.0,),
                       _F(GROUP_MA='x',
                           DX=0.0,),
                       _F(GROUP_MA='y',
                           DY=0.0,),
                       _F(GROUP_MA='z',
                           DZ=-50,),
                        ),
                    
                  LIAISON_UNIF=(_F(GROUP_MA=('lia','saika1','saika2','saika3','saika4','saika5','saika6',),
                                 DDL=('DZ',),),),
                    );



inst=DEFI_LIST_REEL(DEBUT=0,
                    INTERVALLE=_F(JUSQU_A=1,
                                  NOMBRE=50,),);


fonc1=DEFI_FONCTION(
                    NOM_PARA='INST',
                    VALE=(0,0,
                          1,1,
                          ),);

RESU=STAT_NON_LINE(COMPORTEMENT=(_F(RELATION='ELAS',
                                    DEFORMATION='GROT_GDEP',
                                    ),
                                 ),
                   MODELE=MODE,
                   CHAM_MATER=MATE,
                   EXCIT=_F(CHARGE=CHAR,
                            FONC_MULT=fonc1,),
                   INCREMENT=_F(LIST_INST=inst,),
                   NEWTON=_F(PREDICTION='TANGENTE',
                             MATRICE='TANGENTE',
                             REAC_ITER=1,),
                   CONVERGENCE=_F(RESI_GLOB_RELA=0.01,
                                  ITER_GLOB_MAXI=100,),
                   ARCHIVAGE=_F(LIST_INST=inst,
                                CHAM_EXCLU='VARI_ELGA',),);



RESU=CALC_CHAMP(CRITERES=('SIEQ_ELNO',),
                CONTRAINTE=('SIEF_ELNO',),
                reuse=RESU,
               MODELE=MODE,
               CHAM_MATER=MATE,
               RESULTAT=RESU,
               EXCIT=_F(CHARGE=CHAR,),);


RESU=CALC_CHAMP(CRITERES=('SIEQ_NOEU',),
                  CONTRAINTE=('SIGM_NOEU','SIEF_NOEU',),
                  reuse=RESU,
             RESULTAT=RESU,);

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=_F(MAILLAGE=MAIL,
                  RESULTAT=RESU,
                  NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),),);

FIN();
#CHECKSUM:accf509e4411ddb0702f1d47bd602908  -:FIN CHECKSUM

