DEBUT();

urakome=DEFI_MATERIAU(ELAS=_F(E=16800000.0,  #tekitou(Pa)
                         NU=0.33,),);

isi=DEFI_MATERIAU(ELAS=_F(E=36750000000.0,  #heikin(Pa):(1.5~6)*10^5 kg/cm^2
                         NU=0.3,),);

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

#MAIL=MODI_MAILLAGE(reuse=MAIL,
#                   MAILLAGE=MAIL,
#                   ORIE_PEAU_3D=_F(GROUP_MA=('sita',),),
#                   );

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

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
                   AFFE=(_F(GROUP_MA='urakome',
                            MATER=urakome,),
        	         _F(GROUP_MA='isi',
	                    MATER=isi,),),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
                    DDL_IMPO=(
                        _F(GROUP_MA='hidari',
                           DY=0.0,),
                        _F(GROUP_MA='migi',
                           DY=0.0,),
                        _F(GROUP_MA='usiro',
						   DX=0.0,),
                        _F(GROUP_MA='sita',
                           DZ=0.0,),
                        ),
                    FORCE_INTERNE=(             #tannitaiseki
                        _F(GROUP_MA='urakome',
                           FZ=-14700.0,),       #jijuu(N/m^3):gamma=14.7kN/m^3
						_F(GROUP_MA='isi',
						   FZ=-24500.0,),       #jijuu(N/m^3):rho=2.50g/cm^3
                        ),
                    );

RESU=MECA_STATIQUE(MODELE=MODE,
                   CHAM_MATER=MATE,
                   EXCIT=_F(CHARGE=CHAR,),);

RESU=CALC_ELEM(reuse=RESU,
               MODELE=MODE,
               CHAM_MATER=MATE,
               RESULTAT=RESU,
               OPTION=('SIGM_ELNO','SIEQ_ELNO',),
               EXCIT=_F(
               CHARGE=CHAR,),);

RESU=CALC_NO(reuse=RESU,
             RESULTAT=RESU,
             OPTION=('SIGM_NOEU', 'SIEQ_NOEU', ),);

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

FIN();
