
DEBUT(PAR_LOT='NON',);

sig400=DEFI_FONCTION(
                     NOM_PARA='EPSI',
                     NOM_RESU='SIGM',
                     VALE=(0.01 ,40.0 ,
                           0.02 ,80.0 ,
                           0.04 ,160.0 ,),
                     PROL_DROITE='LINEAIRE',
                     PROL_GAUCHE='LINEAIRE',);

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

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

MAILf=CREA_MAILLAGE(MAILLAGE=MAIL,
                    MODI_MAILLE=_F(TOUT='OUI',
                                   OPTION='TRIA6_7',),);

MODE=AFFE_MODELE(MAILLAGE=MAILf,
                 AFFE=_F(TOUT='OUI',
                         PHENOMENE='MECANIQUE',
                         MODELISATION='COQUE_3D',),);

cara=AFFE_CARA_ELEM(MODELE=MODE,
                    COQUE=_F(GROUP_MA='Group_3',
                             EPAIS=0.1,
                             COQUE_NCOU=1,),);

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

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
                    DDL_IMPO=(_F(GROUP_MA='Group_3',
                                 DZ=0.0,),
                              _F(GROUP_NO='Group_2',
                                 DX=-60.0,
                                 DZ=0.0,),
                              _F(GROUP_NO='Group_1',
                                 DX=0.0,
                                 DY=0.0,
                                 DZ=0.0,),),);

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

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

RESU=STAT_NON_LINE(MODELE=MODE,
                   CHAM_MATER=MATE,
                   CARA_ELEM=cara,
                   EXCIT=_F(CHARGE=CHAR,
                            FONC_MULT=fonc1,),
                   COMPORTEMENT=_F(RELATION='ELAS',
                                   DEFORMATION='GROT_GDEP',),
                   INCREMENT=_F(LIST_INST=inst,),
                   NEWTON=_F(PREDICTION='TANGENTE',
                             MATRICE='TANGENTE',
                             REAC_ITER=1,),
                   CONVERGENCE=_F(RESI_GLOB_RELA=1.0,
                                  ITER_GLOB_MAXI=50,),
                   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='SIGM_ELNO',
                CRITERES='SIEQ_ELNO',);

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

PROJ=PROJ_CHAMP(RESULTAT=RESU,
                MODELE_1=MODE,
                MODELE_2=MODE2,
                NOM_CHAM=('DEPL','SIEF_NOEU','SIEQ_NOEU',),);

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=_F(RESULTAT=PROJ,),);

FIN();
#CHECKSUM:2cc8c994362772d1dff54a18b9b060f8  -:FIN CHECKSUM
