Hi Jingjing,
In order to calculate a moving load, say of a train, on a bridge, here is what I do :
==================================================================================================
V=0.1; # Speed of the train
L_TRAIN=1.; # Length of the train
P=1000.; # Weight of the train
FONC_T = FORMULE(NOM_PARA=('INST','V','P','L_TRAIN','X'),VALE="""P*(((V*INST-X)>0)-((V*INST-L_TRAIN-X)>0))*(V*INST-X)*(V*INST-L_TRAIN-X)""");
FORCE_T=AFFE_CHAR_MECA_F(MODELE=MODELE,
FORCE_NODALE=_F(GROUP_NO='BRIDGE',
FY=FONC_T)
);
==================================================================================================
This results in a parabolic load in the Y direction moving at a speed V on the bridge in the X direction. All you have to do is to apply this load on the bridge and run a dynamic calculation.
I think it should help you.
AA