Hello, I am trying to load a component with a harmonic vibration from the outside and would like to see how the component behaves with regard to vibrations.
This type of simulation in Salome is completely new to me. I have seen a video on dynamic simulations and tried to transfer what was shown to a simple beam. The beam is coupled with a node. I would like to load this node with a force as a vibration and later read out the resulting vibrations in the component. unfortunately the calculation does not work as I would like it to, I would be grateful for help and tips on the project!
Am I perhaps going down the wrong path or are there simpler options? the whole thing has to be done with a solid!

Many thanks in advance!

DEBUT(LANG='EN')

PI10 = 10 * pi

PI3 = 3 * pi

Beam = LIRE_MAILLAGE(FORMAT='MED',
UNITE=20)

mesh0 = CREA_MAILLAGE(CREA_POI1=_F(GROUP_NO=('load', ),
NOM_GROUP_MA=('load', )),
MAILLAGE=Beam)

Modell = AFFE_MODELE(AFFE=(_F(GROUP_MA=('Beam', ),
MODELISATION='3D',
PHENOMENE='MECANIQUE'),
_F(GROUP_MA=('load', ),
MODELISATION='DIS_TR',
PHENOMENE='MECANIQUE')),
MAILLAGE=mesh0)

elempro0 = AFFE_CARA_ELEM(DISCRET=(_F(CARA='K_T_D_N',
GROUP_MA=('load', ),
VALE=(0.1, 0.1, 0.1)),
_F(CARA='M_T_D_N',
GROUP_MA=('load', ),
VALE=1.0)),
MODELE=Modell)

Stahl = DEFI_MATERIAU(ELAS=_F(E=210000.0,
NU=0.3,
RHO=5e-09))

fieldmat = AFFE_MATERIAU(AFFE=_F(MATER=(Stahl, ),
TOUT='OUI'),
MODELE=Modell)

sinf = FORMULE(NOM_PARA=('INST', ),
VALE='sin(2pi5.0*INST)')

listr = DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=PI3,
PAS=0.1))

sinvals = CALC_FONC_INTERP(FONCTION=sinf,
LIST_PARA=listr)

bc = AFFE_CHAR_MECA(DDL_IMPO=F(BLOCAGE=('DEPLACEMENT', ),
GROUP_MA=('FIX', )),
LIAISON_SOLIDE=F(GROUP_NO=('Rigid_con', )),
MODELE=Modell)

load = AFFE_CHAR_MECA(FORCE_NODALE=_F(FZ=1000.0,
GROUP_NO=('load', )),
MODELE=Modell)

ASSEMBLAGE(CARA_ELEM=elempro0,
CHARGE=(bc, ),
MATR_ASSE=(_F(GROUP_MA=('load', ),
MATRICE=CO('mass'),
OPTION='MASS_MECA'),
F(GROUP_MA=('load', ),
MATRICE=CO('rigi'),
OPTION='RIGI_MECA')),
MODELE=Modell,
NUME_DDL=CO('nddl'),
VECT_ASSE=F(CHARGE=(load, ),
OPTION='CHAR_MECA',
VECTEUR=CO('fz')))

resharm = DYNA_VIBRA(BASE_CALCUL='PHYS',
CARA_ELEM=elempro0,
EXCIT=F(FONC_MULT=sinf,
VECT_ASSE=fz),
INCREMENT=F(LIST_INST=listr),
MATR_MASS=mass,
MATR_RIGI=rigi,
MODELE=Modell,
TYPE_CALCUL='TRAN')

resharm = CALC_CHAMP(reuse=resharm,
CONTRAINTE=('SIEF_NOEU', ),
CRITERES=('SIEQ_NOEU', ),
GROUP_MA=('Beam', ),
RESULTAT=resharm)

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