
DEBUT();

sig1=DEFI_FONCTION(
                   NOM_PARA='EPSI',
                   NOM_RESU='SIGM',
                   VALE=(0.001 ,210.0 ,
                         0.4 ,400.0 ,),
                   PROL_DROITE='LINEAIRE',
                   PROL_GAUCHE='LINEAIRE',);

sig2=DEFI_FONCTION(
                   NOM_PARA='EPSI',
                   NOM_RESU='SIGM',
                   VALE=(0.001 ,1000000.0 ,
                         0.4 ,400000.0 ,),
                   PROL_DROITE='LINEAIRE',
                   PROL_GAUCHE='LINEAIRE',);


MA=DEFI_MATERIAU(ELAS=_F(E=210000.0,
                         NU=0.3,
                         RHO=7.87e-9,),
                 TRACTION=_F(SIGM=sig1,),
                 PRAGER=_F(C=318.2,),);


MA1=DEFI_MATERIAU(ELAS=_F(E=1000000000.0,
                         NU=0.0,
                         RHO=0.0,),
                 TRACTION=_F(SIGM=sig2,),
                 PRAGER=_F(C=318.2,),);


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(GROUP_MA='ita',
                           MATER=MA1,),
                         _F(GROUP_MA='kou',
                           MATER=MA,),),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
                    DDL_IMPO=(_F(GROUP_MA='x',
                                 DX=0.0,),
                              _F(GROUP_MA='y',
                                 DY=0.0,),
                              _F(GROUP_MA='x1',
                                 DX=0.0,),
                              _F(GROUP_MA='y1',
                                 DY=0.0,),
                              _F(GROUP_MA='z1',
                                 DZ=0.0,),),
                    FORCE_FACE=_F(GROUP_MA='saika',
                                  FZ=220.0,),
                    LIAISON_UNIF=_F(GROUP_MA=('lia','lia1',),
                                    DDL=('DZ',),),
);


inst=DEFI_LIST_REEL(DEBUT=0.0,
                    INTERVALLE=(_F(JUSQU_A=0.08181818,
                                   NOMBRE=210,),
                                _F(JUSQU_A=1.8,
                                   NOMBRE=70,),
                                _F(JUSQU_A=1.96363636,
                                   NOMBRE=220,),
                                _F(JUSQU_A=5.4,
                                   NOMBRE=140,),
                                _F(JUSQU_A=5.56363636,
                                   NOMBRE=420,),
                                _F(JUSQU_A=9.0,
                                   NOMBRE=140,),
                                _F(JUSQU_A=9.16363636,
                                   NOMBRE=420,),
                                _F(JUSQU_A=12.6,
                                   NOMBRE=140,),
                                _F(JUSQU_A=12.76363636,
                                   NOMBRE=420,),
                                _F(JUSQU_A=16.2,
                                   NOMBRE=140,),
                                _F(JUSQU_A=16.36363636,
                                   NOMBRE=420,),
                                _F(JUSQU_A=19.8,
                                   NOMBRE=140,),
                                _F(JUSQU_A=19.96363636,
                                   NOMBRE=420,),
                                _F(JUSQU_A=23.4,
                                   NOMBRE=140,),
                                _F(JUSQU_A=23.56363636,
                                   NOMBRE=420,),
                                _F(JUSQU_A=27.0,
                                   NOMBRE=140,),
                                _F(JUSQU_A=27.16363636,
                                   NOMBRE=420,),
                                _F(JUSQU_A=30.6,
                                   NOMBRE=140,),
                                _F(JUSQU_A=30.76363636,
                                   NOMBRE=420,),
                                _F(JUSQU_A=34.2,
                                   NOMBRE=140,),
                                _F(JUSQU_A=34.36363636,
                                   NOMBRE=420,),
                                _F(JUSQU_A=37.8,
                                   NOMBRE=140,),),);

dinst=DEFI_LIST_INST(DEFI_LIST=_F(METHODE='AUTO',
                                  LIST_INST=inst,
                                  PAS_MINI=1e-05,
                                  PAS_MAXI=0.5,
                                  NB_PAS_MAXI=1000000,),);

fonc1=DEFI_FONCTION(
                    NOM_PARA='INST',
                    VALE=(0.0 ,0.0 ,
                          0.08181818 ,0.95454545 ,
                          1.8 ,0.97851664 ,
                          1.96363636 ,-0.89046 ,
                          5.4 ,-1.02171664 ,
                          5.56363636 ,0.927914618 ,
                          9.0 ,0.97851664 ,
                          9.16363636 ,-0.89046 ,
                          12.6 ,-1.02171664 ,
                          12.76363636 ,0.927914618 ,
                          16.2 ,0.97851664 ,
                          16.36363636 ,-0.89046 ,
                          19.8 ,-1.02171664 ,
                          19.96363636 ,0.927914618 ,
                          23.4 ,0.97851664 ,
                          23.56363636 ,-0.89046 ,
                          27.0 ,-1.02171664 ,
                          27.16363636 ,0.927914618 ,
                          30.6 ,0.97851664 ,
                          30.76363636 ,-0.89046 ,
                          34.2 ,-1.02171664 ,
                          34.36363636 ,0.927914618 ,
                          37.8 ,0.97851664 ,),);

RESU=DYNA_NON_LINE(MODELE=MODE,
                   CHAM_MATER=MATE,
                   EXCIT=_F(CHARGE=CHAR,
                            FONC_MULT=fonc1,),
                   AMOR_RAYL_RIGI='ELASTIQUE',
                   COMPORTEMENT=_F(RELATION='VMIS_ECMI_TRAC',
                                   DEFORMATION='PETIT_REAC',),
                   INCREMENT=_F(LIST_INST=dinst,),
                   SCHEMA_TEMPS=_F(SCHEMA='NEWMARK',
                                   BETA=0.25,
                                   GAMMA=0.5,
                                   FORMULATION='DEPLACEMENT',),
                   NEWTON=_F(PREDICTION='TANGENTE',
                             MATRICE='TANGENTE',
                             REAC_ITER=1,),
                   CONVERGENCE=_F(RESI_GLOB_RELA=0.001,
                                  ITER_GLOB_MAXI=10,),
                   ARCHIVAGE=_F(LIST_INST=inst,
                                CHAM_EXCLU='VARI_ELGA',),);

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

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

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

FIN();
#CHECKSUM:25886a0c5aa5a31820b763ef3e308c38  -:FIN CHECKSUM
