SHLL101 harmonic analysis

Hello,

i'm trying to perform an harmonic analysis using code_aster.
I've been trying to do the SHLL101 validation case but i don't manage to obtain the same results.

Please find below what i did (.comm file).
I used DYNA_LIM_HARM.

Thks in advance.

Thajup

DEBUT();

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

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

MAT=DEFI_MATERIAU(ELAS=_F(E=1.658,
NU=0.3,
RHO=1.3404106104,
AMOR_ALPHA=0.001,
AMOR_BETA=0,),);

CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
MATER=MAT,),);

elem=AFFE_CARA_ELEM(MODELE=MODELE,
POUTRE=_F(GROUP_MA='mailla',
SECTION='GENERALE',
CARA=('A','IY','IZ','JX',),
VALE=(3.439e-3,1.377e-5,1.377e-5,2.754e-5,),),);

BLOCAGE=AFFE_CHAR_MECA(MODELE=MODELE,
DDL_IMPO=_F(GROUP_NO='A',
LIAISON='ENCASTRE',),);

fono=AFFE_CHAR_MECA(MODELE=MODELE,
FORCE_NODALE=_F(GROUP_NO='B',
FX=0,
FY=3000,
FZ=0,),);

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

jkh=CALC_VECT_ELEM(OPTION='CHAR_MECA',
CHARGE=fono,
CHAM_MATER=CHMAT,);

jjd=ASSE_VECTEUR(VECT_ELEM=jkh,
NUME_DDL=NUMEDDL,);

MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
MATR_B=MASSE,
METHODE='TRI_DIAG',
TYPE_RESU='DYNAMIQUE',
CALC_FREQ=_F(OPTION='PLUS_PETITE',
NMAX_FREQ=15,),);

MACRO_PROJ_BASE(BASE=MODES,
NB_VECT=999,
PROFIL='DIAG',
MATR_ASSE_GENE=(_F(MATRICE=CO('M_gen'),
MATR_ASSE=MASSE,),
_F(MATRICE=CO('K_gen'),
MATR_ASSE=RIGIDITE,),),
VECT_ASSE_GENE=_F(VECTEUR=CO('V_gen'),
VECT_ASSE=jjd,),);

freqs=DEFI_LIST_REEL(DEBUT=0.01,
INTERVALLE=_F(JUSQU_A=500,
PAS=10,),);

hamor=DYNA_LINE_HARM(CHAM_MATER=CHMAT,
CARA_ELEM=elem,
MATR_MASS=M_gen,
MATR_RIGI=K_gen,
LIST_FREQ=freqs,
TOUT_CHAM='OUI',
EXCIT=_F(VECT_ASSE=V_gen,
COEF_MULT=1,),);

teqt=REST_GENE_PHYS(RESU_GENE=hamor,
NOM_CHAM='DEPL',);

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

IMPR_RESU(FORMAT='RESULTAT',
RESU=_F(RESULTAT=teqt,
NOM_CHAM='DEPL',
LIST_FREQ=freqs,
NOM_CMP='DY',),);

resu1=RECU_FONCTION(RESULTAT=teqt,
LIST_FREQ=freqs,
NOM_CHAM='DEPL',
NOM_CMP='DY',
GROUP_NO='B',
NOM_PARA='FREQ',);

FIN();