I am putting together some simple ride comfort models that use road surface PSDs from ISO 8606 (road surface profiles) as inputs to a quarter car suspension model to produce body acceleration spectra. These are then weighted according to ISO 2631 (human exposure to whole body vibration).
The road surface data is vertical displacement PSDs referencing spatial frequency (wave number; inverse of wavelength). The suspension model is simply a gain of body displacement per unit of ground displacement over a range of temporal frequencies. The input to the weighting process needs to be an RMS acceleration with respect to temporal frequency. So, the problem is to convert the road surface profile equation from a PSD against spatial frequency to an RMS against temporal frequency.
The following is the process I am using:
- Convert road surface displacement PSD to acceleration PSD (given in ISO 8608)
- Convert road surface acceleration PSD to acceleration RMS (still with respect to spatial frequency)
- Convert road surface acceleration RMS to displacement RMS by dividing by square of spatial frequency
- Convert road surface displacement RMS to acceleration RMS by multiplying by square of temporal frequency
- Multiply by suspension gain to get sprung mass acceleration
- Weight sprung mass acceleration using filter described in ISO 2631.
However, the sprung mass acceleration results (pre or post weighting) are rather high, so I just wanted to check whether the above process was sound. The units seem to check out so that gives me some confidence.
Can anyone confirm for me that steps 2-4 in the above process are correct?
Thanks.