Human Vibration

The OROS Body Vibration tool allow you to evaluate the affect of vibration on the human body according to standards ISO 2651 and ISO 5349. These standard define measurement practice and vibration signal analysis to evaluate the affect on health and and comfort of environmental and equipment vibrations on the Human body.

The ISO 2651 describe the affect on health and comfort of vibration on the whole-body in transportation system, and the ISO 5349 the affect on health of vibration on hands and arms when manipulating machine-tools or vibrating objects. In the following, we will see how to use OROS to evaluate these effect.

This is a post analysis toolkit, operating alongside NVGate after the recording of the signal and will help you to calculate time-weighted signal of acceleration and specific indicators both defined in the standards.

File:OROS BodyVib2.png

Download & install

Download the program here (version from 03/12/2020) : https://orossas.sharepoint.com/:u:/g/support/EXCVJE_bLGZMgHg1mr4p1hwB5ICqdqKBXTK3kH9stxozKg?e=UJg72W

Once downloaded, you can unzip the folder and lunch the installer program Setup_OROS_BodyVibration_Tool_vXX.exe. Follow the step of the program to properly install the software.

The program will be installed in "C:\OROS\Programs\ExternalTools". Shortcuts are created on the main desktop, in the Windows Start program and in the "Link" repertory of NVGate. By that way you can directly run the program from NVGate.

How to use

1. Launch NVGate
   
   
2. Launch OROS_BodyVibration_Tool.exe
   
   
3. Open a project in NVGate
   


4. In the tool, Click on list measurement to list all the measurement present in the project (only the measurement where the signal files are actually on the disk will be displayed. Ensure you have downloaded signal from the OR3X disk)
   
   

5. Add the measurement you want to analyse in the "selected list" area using the ">>" button. Use the "<<" button to remove file if needed.
   
   

6. Click on "Select Directions" to set the direction and sensor for each channel (This step is not necessary if you only want the time-weighted signals).
   
   

7. Select the weighting you want to apply on the time signal
   
   

8. Select the weighting you want to apply for each direction X, Y and Z used for the indicators
   
   

9. Select the metrics you want to calculate
   
   

10. Click on "Start Processing"
    
    

11. The processing is completed when a the following pop-up window appear. The results and weighted signal have added to the current NVGate project. See the next sections for details.
    
    

Toolkit results

Time-weighted signal

Studies of the affect of vibrations on the human body suggest that part of the body don't have the same response to vibrations, and will be sensitive to specific frequencies and transient component of vibration signal. In order to represent the sensitivity of the human boy, ISO standards defined time-weighting that must applied to the signal according to the environmental conditions. The following part present you how to apply time-weighting filters to your signal using the toolkit.

Time-weighting filters

As defined in the standards, specific time weighting must be applied to the acceleration signal in order to represent the effect of vibrations on health and comfort. Here is the list of the different time-weighting filters implemented in the toolkit :

Time-weighting filters

$W_{k}$	Time weighting for the Z axis for whole-body measurement (ISO 2651-1)
$W_{d}$	Time weighting for the X and Y axis for whole-body measurement (ISO 2651-1)
$W_{h}$	Time weighting for the hand-arms measurement in any direction (ISO 5349-1)
$W_{f}$	Time weighting for motion sickness measurement in the vertical direction (ISO 2651-1)
$W_{c}$	Time weighting for the X axis for whole-body measurement (ISO 2651-1)
$W_{e}$	Time weighting for all rotational directions for whole-body measurement (ISO 2651-1)
$W_{j}$	Time weighting for the Z axis for head comfort measurement (ISO 2651-1)
$W_{m}$	Time weighting for all directions for building's vibration comfort measurement (ISO 2651-2)

Filters have been validated from 0.5 to 10kHz at 25.6kSample/second.

Usage in the toolkit

The toolkit allow you to calculate the filtered raw signal with any of the above weighting filters. This will allow you to display spectrums and calculate RMS via post-analysis in NVGate.

To calculate the raw weighted signal, simply select the filters you want in the "Time signal filters" section :

Once the processing is completed, a new measurement containing the filtered signal for each selected filter is created in NVGate Project Manager :

You can then load the new signals in the player as any NVGate signal. Each track of the signal is filtered with the corresponding filter.

Health and comfort indicator

In addition to the time-weighting filters, the ISO standards present a series of specific metrics to evaluate the affect of vibrations on the Human body. The following section will present you how to use the toolkit to access these metrics.

Indicators

Body vibration indicators

$A_{w}$	RMS of the weighted acceleration magnitude.
$a_{v}$	Total RMS value calculated as a quadratic average of the three direction for each tri-axial sensor.
$A_{h}$	RMS of the $W_{h}$ weighted acceleration magnitude.
$Peak$	Amplitude Peak of the vibration; maximal amplitude of the signal from the 0 (Calculated from Post processing in NVGate)
$Crestfactor$	Ratio between the Peak level and the RMS of the weighted signal(Calculated from Post processing in NVGate)
$MTVV$	Maximum Transient Vibration Value, represent the maximal RMS value of the signal
$VDV$	Vibration Dose Value, taking into account the temporal shocks in the signal
$MSDV$	Motion Sickness Dose Value, representing the comfort in transportation measurement
$D_{k}$	The Acceleration Dose, representing the affect of the vibration on the spine
$Se$	The equivalent Static compression stress on the spine, caused by repeated shocks

Aw

The $A_{w}$ is the RMS of the weighted acceleration signal. The weighting is applied in accordance with the direction set for each channel and the weighting filter defined for the direction.

For each channel, the $A_{w}$ and $A_{w}(T)$ will be calculated. The $A_{w}(T)$ is the Daily exposure value, defined as :

$$A_{w}(T)=A_{w}*k_{i}*{\sqrt {\frac {T}{T_{m}}}}$$

With $T_{m}$ the duration of the measurement and $T$ the total exposure duration represented by the "Reference time" parameter. the $k_{i}$ factor is defined in the ISO 2631 as $k_{X}=k_{Y}=1.4$ and $k_{Z}=1$.

When comparing result from one measurement point in the three direction, the maximum of the Daily exposure value must be used as the total daily exposure value at that point :

$$A_{w}(T)=max(A_{wX}(T),A_{wY}(T),A_{wZ}(T))$$

av

The $a_{v}$ is the total RMS vibration value of the weighted acceleration signal. This value will only be calculated if a sensor have been defined in the definition of direction window.

For each tri-axial sensor, the $a_{v}$ will be calculated as a quadratic average of the weighted RMS value in the three directions :

$$a_{v}={\sqrt {{a_{vX}}^{2}*{k_{X}}^{2}+{a_{vY}}^{2}*{k_{Y}}^{2}+{a_{vZ}}^{2}*{k_{Z}}^{2}}}$$

With $k_{X,Y,Z}$ a factor defined in the ISO 2631 as $k_{X}=k_{Y}=1.4$ and $k_{Z}=1$

Ah

The $A_{h}$ is the RMS of the $W_{h}$ weighted acceleration signal for hand-arms vibration measurement. For each channel, the $A_{h}$ and $A_{h}(T)$ will be calculated. The $A_{h}(T)$ is the Daily exposure value, defined as :

$$A_{h}(T)=A_{h}*{\sqrt {\frac {T}{T_{m}}}}$$

With $T_{m}$ the duration of the measurement and $T$ the total exposure duration represented by the "Reference time" parameter.

If any sensor is defined, the total $A_{h}$ and $A_{h}(T)$ will be calculated for each sensor as a quadratic average of the three direction of the sensor :

$$A_{hTotal}={\sqrt {{a_{hX}}^{2}*{k_{X}}^{2}+{a_{hY}}^{2}*{k_{Y}}^{2}+{a_{hZ}}^{2}*{k_{Z}}^{2}}}$$

and :

$$A_{hTotal}(T)={\sqrt {{a_{hX}}^{2}*{k_{X}}^{2}+{a_{hY}}^{2}*{k_{Y}}^{2}+{a_{hZ}}^{2}*{k_{Z}}^{2}}}*{\sqrt {\frac {T}{T_{m}}}}$$

With $T_{m}$ the duration of the measurement and $T$ the total exposure duration represented by the "Reference time" parameter, and $k_{X}=k_{Y}=k_{Z}=1$

Peak

This indicator is the peak amplitude of the vibration signal. It represent the maximal amplitude of the weighted signal. This indicator can be obtained in NVGate by post-analysing the weighted signal with the TDA plugin of NVGate.

MTVV

The MTVV is the Maximal Transient Vibration Value, is the maximum of the running RMS of the weighted acceleration. This value is calculated for each channel and is expressed in m/s².

The MTVV Tau is the integration time, expressed in seconds. This value can be edited before the calculation, but the ISO 2631 recommand to use 1 second. Therefore this is the default value of this parameter.

VDV

The VDV is the forth power Vibration Dose Value. It is calculated as a cumulative sum of the weighted acceleration at the forth power. The forth power is used to give more weight to the transient event of the signal over the periodic event. Therefore, this indicator must be used to asses the effect of shock on human health and comfort. The VDV is calculated for each channel and expressed in m/s^(1.75).

If the "VDVTexp" box is checked, an evaluation of the VDV value over a longer time period will be calculated of each channel as :

$$VDV_{X,Y,Z}(T)=VDV_{X,Y,Z}*k_{X,Y,Z}*{\sqrt[{4}]{\frac {T}{T_{m}}}}$$

With $T_{m}$ the duration of the measurement and $T$ the total exposure duration represented by the "Reference time" parameter, and $k_{X}=k_{Y}=1.4$, and $k_{Z}=1$

MSDV

This metric is the Motion Sickness Dose Value, used to evaluate the ride comfort in in-vehicle measurement. It is expressed in m.s^(1.5).

Although this value will be calculated for each direction, only the Z-direction will have a prominent relevance to assess ride comfort.

Dk

$D_{k}$ is the acceleration dose. It is use to evaluate the effect of vibrations on the human spine, and calculated as the 6th root of the sum of the peaks in the acceleration signal raised at the 6th power. it is expressed in m/s². Peaks are defined as the absolute maximum of the acceleration signal between two zeros crossing. In the Z direction, only the positive peak will be considered, were all positive and negative peaks are used for the X and Y directions.

Specific filters are used to weight the acceleration signal before calculating the $D_{k}$. These filter are different for the vertical and horizontal direction, and are defined for a sampling rate of 160 Sample/second. Therefore, we advise to record the signal with a low sampling rate in NVGate.

Assessing the affect of the vibration on human spine most is important when dealing with transportation vehicular, especially from driver seat. for a better evaluation, the ISO 2631-5 advise to calculate the daily acceleration dose, defined as :

$$D_{X,Y,Z}(T)=D_{X,Y,Z}*{\sqrt[{6}]{\frac {T}{T_{m}}}}$$

This value will be calculated if the box "Dk Texp" is checked.

Se

The Se is the compressive stress on the human spine, expressed in MPa. It is calculated by the combination the Dk in the three direction. Therefore, it will only be calculated if sensors are defined in the "Select Direction" window. One value of Se will be calculated for each sensor by the formula :

$$Se={\sqrt[{6}]{\sum _{k=X,Y,Z}(D_{k}*m_{k})^{6}}}$$

Where $m_{X}=0,015$ MPa/(m/s²), $m_{Y}=0,035$ MPa/(m/s²) and $m_{Z}=0,032$ MPa/(m/s²).

To evaluate the lumbar stress over a daily exposition time, the Se(T) will also be calculated as :

$$Se(T)={\sqrt[{6}]{\sum _{k=X,Y,Z}(D_{k}(T)*m_{k})^{6}}}$$

T being the total exposure duration represented by the "Reference time" parameter.

In the toolkit

Before calculating the indicators in the toolkit, you must define a filter for each directions. Filters will be selected in accordance with the wanted result (Hand-arms; Whole-Body, ..) and the direction of the measurement as specified in the table describing the filters. Filters will be applied to the signal before calculating the indicators, but the time-filtered signal will not be saved.

You will then need to enter the reference time you want to use to estimate the parameters on a greater duration. This duration must be entered in second. The default value is 288500 second, corresponding to a 8 hours measurement:

You will finally need to enter the MTVV Tau value. The MTVV Tau is the integration time in seconds used for the calculation of the MTVV. The ISO 2651 recommand to use 1sec. This is the default value of the setting.

Once the filters selected and the parameters competed, you can select the metrics you want to calculate by checking the different boxes.



You can then start the processing to calculate the indicators.

Here, we are calculating all the indicator for the "Seat_Pad_Measurement" signal. This signal contain 3 Tracks, one for each channel X, Y and Z of a three axis accelerometer. As we want to calculate the multi-axis metrics, we also define a sensor on the three tracks, called "Seat_Pas_3axis".


Once completed, the result can be find the currently opened project in the NVGate Project manager. A new measurement called "NAME_OF_THE_Signal_out" containing the result is created.

Each indicator is calculated for each direction and for each track of the signal. The name of the result contain the name of the track, the considered direction and the filter used for that direction (except for the $D_{k}$ which use specific filters). The following picture show the indicators for the first track (SeatPadX) of the signal :

File:Result SeatPad.png

Here this track represent the X direction. We can remove the results in irrelevant directions (Y and Z) for this track by a right click, and repeat this operation for each channel :

/Insert image of the result with only relevant direction when bug solve/

You can then display the result in NVGate by double clicking.

Calculating the other indicators with NVGate

$A_{w}$

$A_{w}$ is the daily vibration dose. It is calculated from an average of the RMS values form the weighted signal and represent the "quantity" of vibration the human body will receive when using of the tested equipment.

This value is defined as : $A_{w}={\sqrt {(k_{x}*a_{w_{x}})^{2}+(k_{y}*a_{w_{y}})^{2}+(k_{z}*a_{w_{z}})^{2}}}$, with $k_{i}$ defined in the standard and $a_{w_{x}}$ the RMS value of the weighted acceleration signal in accordance to the direction, as defined in the table of the precedent section.

Calculation in NVGate

To calculate the $A_{w}$ in NVGate, calculmate the RMS level of the weighted signal in Post analysis. Here, we will calculate the $A_{w}$ for the SeatPad measurement as defined for a whole-body measurement in the ISO 2651. According to the [Human_Vibration#Time-weighting_filters|table]], we need to calculate the RMS level of the $W_{k}$ weighted signal for the Z direction and the $W_{d}$ weighted signal for the X and Y directions. You then must display the three RMS level in the same window :

Practical information about measurements

Durée de mesure Type de mesure Placement des capteurs Exemple de capteur : [Pad accelerometer] Capteur avec adaptateurs hand arms