I have made relative displacement calculations for a base mass excitation.
Here is an approach:
1. Perform a "frequency response fuction" analysis over the frequency domain from 0 to 2000 Hz, or even higher as needed. Use a constant frequency increment for the response.
2. Output the displacement FRFs at the base and at a few carefully select response nodes.
3. Manipulate the FRFs to form a transmissibility function relating the relative displacement PSD to the base input PSD. Typically, the base input is an acceleration PSD.
4. Multiply the base input PSD by the appropriate transmissibility function to obtain the relative displacement PSD at the selected node.
5. Integrate the relative displacement PSD and then take the square root to get the relative displacement RMS.
I have performed this analysis using NE/Nastran.
This method requires some postprocessing steps. There are several advantages, however. One is that the transmissibility functions yield additional insight into the behavior of the system. Another is that the transmissibilty functions can be used very efficiently if there are several different PSD inputs.
Tom Irvine