DEBUT();

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

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

urakome=DEFI_MATERIAU(ELAS=_F(E=16800000.0,  #tekitou(Pa)
                               NU=0.33,
                               RHO=1.50e3,  #kg/m^3:14.7kN/m^3
			                   AMOR_ALPHA=0.00004,),); #tekitou

isi=DEFI_MATERIAU(ELAS=_F(E=36750000000.0,  #heikin(Pa):(1.5~6)*10^5 kg/cm^2
                          NU=0.3,
                          RHO=2.50e3,  #kg/m^3
                          AMOR_ALPHA=0.00002,),); #tekitou

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

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

MACRO_MATR_ASSE(MODELE=MODELE,
                CHAM_MATER=CHMAT,
                CHARGE=BLOCAGE,
                NUME_DDL=CO('NUMEDDL'),
                MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
                              OPTION='RIGI_MECA',),
                           _F(MATRICE=CO('MASSE'),
                              OPTION='MASS_MECA',),
						   _F(MATRICE=CO('amore'),
						      OPTION='AMOR_MECA',),
                           ),
                );

loadV=DEFI_FONCTION(NOM_PARA='INST',VALE=(0.0,1.0,
                                          0.079,1.0,
                                          0.080,100.0, #tekitou:jijuu no 100 bai
                                          0.081,1.0,
					                      0.200,1.0,
					                      ),);

stepTR=DEFI_LIST_REEL(DEBUT=0.0,
                      INTERVALLE=_F(JUSQU_A=0.20,
		                    PAS=0.001,),);

MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
                       MATR_B=MASSE,
                       CALC_FREQ=_F(OPTION='PLUS_PETITE',
                                    NMAX_FREQ=5,),
                       );

dynaTR=DYNA_LINE_TRAN(MODELE=MODELE,
                      CHAM_MATER=CHMAT,
					  MATR_MASS=MASSE,
					  MATR_RIGI=RIGIDITE,
					  MATR_AMOR=amore,
					  NEWMARK=_F(DELTA=0.5,),
					  EXCIT=_F(CHARGE=BLOCAGE,
					           FONC_MULT=loadV,),
					  INCREMENT=_F(LIST_INST=stepTR,),);

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

IMPR_RESU(MODELE=MODELE,
          FORMAT='MED',
		  RESU=_F(RESULTAT=dynaTR,
		          NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),),);

FIN();
