How to choose an interpolation method for seismic ground motion signals I am working with seismic ground motion records from the PEER database. They are acceleration measurements that are band limited to 50 Hz and sampled at 200 Hz.

I use these records to drive an explicit finite element analysis of a structure to determine its response under the earthquake. This analysis runs with a 1e-5 time step (100,000 Hz). If I were to input the sampled records directly the FEA program would linearly interpolate the time steps between the samples, which would distort the signal. 

Therefore I need to reconstruct the original signal and resample it at 100,000 Hz (or high enough that the interpolation doesn't affect accuracy much).

One property that I would like to preserve is the double integral of the signal, which gives displacement. This keeps the structure from drifting away over time. I do not care about the calculation time, as it only needs to be done once and is trivial compared to the FEA calculation time.

A previous answer discussed several methods for interpolation but not how to choose from among them.

I have experimented with sinc interpolation using a Kaiser window with an artificially generated signal that I sampled and then reconstructed and I see decent results, but there are visible differences in the time domain.

Is there some way to determine what interpolation method (and parameters in the case of a windowed method) would give the most accurate results?
If there are differences then some likely sources might be  that the window on the Sinc wasn't wide enough, that the data wasn't really sufficiently bandlimited, or too much quantization noise in the samples.

The double integration of anything with system offsets or quantization noise is unlikely to be zero. (See random walk) One solution would be to compute the double integral of the first pass interpolated result and then bias the interpolator (noise filtered if needed) until the resulting double integral goes to zero.