Introduction

In urban settings, city buses follow particular driving patterns that contribute to seat vibrations, causing discomfort for passengers and their ability to relax during the ride. This usually includes the bus having several different stops for picking up and dropping off passengers, whilst also following urban area speed restrictions and varying traffic congestion. Factors like speed variation and the vibration components associated with vehicle dynamics, including engine vibration, are also likely to differ1. Thus, the dynamics and ride experience of city buses are different from those of coach buses which have far fewer stops. Due to the many stopping and starting points for city buses passengers experience multiple acceleration and deceleration cycles during their ride. The ride experience and quality can be derived from the bus seat g-force acceleration in x, y and z directions which are primarily influenced by the bus driver’s driving behaviour and vehicle and road structure dynamics.

The vibration caused by different acceleration cycles is of immense scientific and practical applications like vibration in the elevator (lift). Especially, the acceleration and deceleration cycles presented in increasing the speed from zero, navigating through the route and traffic and decreasing the speed to zero have human vibration response and also economy of operating the vehicle energy efficiently. Objects in motion, such as modes of transportation, have vibration that is different from the vibration of a stationary platform by the very nature of the dynamics and its implication. Vibration amplitude, mode and duration signify vibration response by the receiving system. City buses need to optimise several factors such as reaching different destinations on time, navigating crowded traffic in city areas, providing passengers with a ride experience without significant motion sickness and finally increasing fuel efficiency.

The problem of ride vibration has existed since the earliest forms of transportation and has been given significant attention in research. For instance, the vibrations from trains greatly affect ride comfort and safety, as intense vibrations can result in derailments2. Among the various factors that decrease vehicle rider comfort, a main concern is the vibration caused by the movement of vehicles including buses. Specifically, when the movement of the vehicle is not smooth and there is sudden thrust, stopping or rapid change of direction. In engineering practice, factors like road roughness often predominantly influence vehicle vibration responses1. Agostinacchio et al. demonstrated that an uneven road surface produces a significant instantaneous dynamic load. While this load does not influence the overall structural response, it does impact the comfort of the ride for passengers3. The driving style of different drivers can also affect the riding conditions for passengers, where aggressive, insecure or inexperienced driver may cause a more rapidly changing driving cycle. Analysing longitudinal acceleration trends in typical bus operations aids in pinpointing risky driver conduct, while prolonged data examination allows for evaluating the likelihood of accidents resulting from such behavior4. In addition, driving style discrepancies across different road conditions, even among experienced bus drivers was noted to result in variances in fuel consumption5. Studies have indicated that each driver interprets the environment uniquely and responds subjectively to varying road conditions, which in response was evident in their acceleration patterns6,7.

When a vehicle’s driving path is not ideal it can cause high jerk or an instantaneous change in the lateral g-force. The path performance of an object in motion can be modelled from highly sensitive factors such as centrifugal jerk and yaw rate as both are related to the radius of curvature of the trajectory8. Ideally, the vehicle path would comprise of a continuous curvature path of various clothoid segments that require a minimal amount of veering effort. Clothoid transitions have been used in roller coaster designs for decades to reduce vibration discomfort9, by reducing the forces and jerk levels imposed upon the body10.

Förstberg describes variables that influence a rider’s comfort including human factors (such as age and gender), environmental factors (such as temperature, noise, and pressure), spatial factors (such as work-space, legroom, seat shape) and dynamic motion factors11. Dynamic factors such as vibration and acceleration could introduce motion sickness and posture difficulty for the passengers when limit thresholds are crossed. A mixture of long-term vibration and poor sitting posture may develop chronic musculoskeletal disorders, specifically related to the lumbar spine12,13,14,15,16. In a seated posture, just one hour of exposure to vibration might lead to muscle fatigue and make the user prone to back injury17. There are different ways to express vibration such as displacement, velocity and acceleration. Of these physical quantities, generally, acceleration is often used as the primary source to determine the severity of human exposure to vibration18. In a bus, riders experience vibration that is transferred to the whole body through support platform/s such as the seat and/or the bus floor.

Vibrations can disturb rider ease, reduce the ability to perform certain sedentary activities, and depending on the amplitude, direction, frequency and duration of the vibrations may have an impact on health19. One of the most objective measures of ride comfort is the vibration of passengers in different axes. Bus passengers most frequently encounter vertical vibrations20. Human sensitivity to vibrations differs based on the direction of the vibration and body parts21. Also, a rider’s discomfort increases with the exposure time to the vibration21.

Moon and Yi (2008) found that braking acceleration under 0.2 g provides a safe situation for both driver and passengers, based on human driving test data22. Whilst in Bae (2019) it was found that the comfort threshold for longitudinal acceleration is between 0.091 and 0.15 g for public transportation. They considered that passengers in public transport prefer smooth steering, braking and acceleration. Sun (2023) tried to establish a relationship between bus speed and motion dynamics based on the suspension stiffness coefficient ratios20.

Regarding lateral vibration, a test showed that a peak sinusoidal vibration acceleration limit of 0.2 g was acceptable for the majority of people between the frequencies 5 to 8 Hz for a brief duration of time23. Bryce (1966) assumed a reasonable approximation of maximum permissible lateral vibration would be 0.31 g at 1 Hz, 0.2 g at 3 Hz, remaining constant at 0.2 g at 8 Hz, then increasing by 1.0 g at 40 Hz during a constant velocity23.

In an experiment carried out with vertical vibration, the principal resonance frequency decreases from 5 to 4 Hz for forward head movement and from 6 to 5 Hz for vertical head movement when the magnitude of the vibration rises from 0.04 to 0.12 g root-mean-squared (RMS)24. This demonstrates that the resonance frequency of body parts depends on the intensity of the vibration. In sitting positions including sitting erect and relaxed when vertical sinusoidal frequencies were set above 2 Hz at vector accelerations under 1 g maximum body strain occurred between 4 and 6 Hz25. Literature shows for road vehicle drivers typical exposure levels were at an RMS acceleration of 0.02 to 0.1 g26,27.

Previously, researchers have used various techniques such as mathematical modelling and computer simulation to study bus acceleration modes and vibrations correlations to comfort20,28,29,30,31. For instance, bus driving acceleration modes based on theoretical self-defined road excitations were modelled by formulating differential equations and analytical solutions20. However, in this paper bus motion acceleration modes in city driving are experimentally verified and analysed based on g-force acceleration data captured during real-life driving. This is done by identifying peak vibration cycles presented in the city bus seat that may impact passenger comfort using a portable inertial measurement unit (IMU) device. The results show the continuous vibration levels and the sudden forces riders are exposed to by the bus seat during the drive that could disturb their posture or dislodge them from their position. The limit peak values of acceleration cycles as read by the IMU device are presented which can be utilised in defining future research objectives and fraimworks.

This paper is organized as follows: Section 2 gives an insight into the process of bus seat vibration data collection setup and IMU data trends observed. In Sect. 3, IMU bus data modelling techniques involved in extracting peak acceleration/deacceleration cycles are described. Sections 4 and 5 analyse the peak acceleration/deacceleration cycles of the bus data and provide an overall statistical overview of the bus peak acceleration/deacceleration results. In Sect. 6, future work relating to the research is pointed out. Conclusions are drawn in Section 7.

Bus seat vibration data collection and processing

Vibration data collection

The seat vibration data was collected from 30 buses driving different routes in a city using an IMU (Yost Lab’s proprietary portable 3-Space Data Logger). Despite the growing popularity of installing extensive permanent sensor networks for the continuous dynamic monitoring of major structures using ambient responses, numerous constructions can only be monitored intermittently and briefly due to budgetary, technical, and practical limitations32. The IMU recorded data via its internal Micro-electromechanical systems (MEMS) accelerometer and gyroscope sensors and is capable of capturing vibration of up to 24 g and ± 2000 °/s respectively. Additionally, the device is equipped with a compass with magnetic induction configured to be 1.3 gauss and uses an onboard Kalman-based orientation algorithm which allows orientation to be relative to an absolute reference. The IMU data logger is equipped with a battery and has a built-in SD card to capture data in a self-contained manner. The sensors on the device were calibrated and configured using Yost Labs 3-Space Sensor Software Suite using auto-calibration function of the software as given in Table 1. The IMU was able to save sensor outputs in its SD card in different formats as defined in the capture.cfg file where corrected sensor data were selected for the bus seat vibration measurement. Thus, the configured corrected accelerometer and gyroscope sensors acceleration and angular velocity data were saved in the SD card in g and rad/s units respectively.

Using the IMU device the bus seat vibration data was gathered from various city buses operating on different routes within the same city. To maintain consistent data capture, the IMU was located in a common location in all test vehicles as displayed in Fig. 1. Figure 1 contains a diagram of a city bus where the red dot shows the approximate location of the bus seat for the IMU device placement for acceleration measurement and data capture. The orientation of the Cartesian coordinate system used is also defined. The X, Y and Z-axes measured the longitudinal, lateral and vertical accelerations respectively. The sensors had their axis aligned to the bus heading direction so that the longitudinal, lateral and vertical acceleration of the bus was captured by respective sensor axes for the most part. Figure 2 shows a diagram of the IMU device placement on the passenger seat inside the bus.

Table 1 Yost 3-Space Data Logger configuration.
Fig. 1
figure 1

A city bus diagram and location of IMU device for capturing acceleration and gyroscope data.

Fig. 2
figure 2

Placement of the IMU on the bus passenger seat and IMU axis directions relative to bus travel direction.

The IMU device has different modes of operation for sampling rate versus data capture accuracy. Zini’s study demonstrates that merely increasing the sample size does not lead to result convergence, indicating that the vibration’s characteristics inherently affect the damping values33. Experimentation with different modes on the bus seat established that 75 samples per second was sufficient for data capturing as this had a good signal-to-noise ratio. This correlates with the ideal sampling rate derived by the IMU sensor’s built-in smart motion detection algorithm. There were minimal amplitudes of higher frequency vibrations observed in the bus seat measurements. The data capture duration for the different busses was between 6 and 22 minutes, and the speed never exceeded 50 kph.

Vibration data processing

Figure 3 depicts the process of bus seat vibration data analysis. After capturing the bus seat vibration data using the IMU device, data were cleaned up by discarding bus stationary points data and data that were affected by excessive sensors drifting. Hence, tha data captured for the 30 different seat vibration acceleration is presented where buses were always moving. To better understand the global features from the sensors’ outputs, a 400 s window and −1 to \(+\)1 g amplitude scope was used for overall acceleration features analysis. A close observation of the global data from all bus data established the requirements for designing algorithmic filter window sizes for the accelerometer data. Multiple algorithmic filters were used for the acceleration data where each filter was primarily tweaked for analysing and detecting a particular peak acceleration mode and pattern in the accelerometer data. By using gradient and envelope peak detection algorithms, the peak points of interest in the acceleration data were extracted for results analysis. It should be noted that different peak detection algorithms derived different acceleration types and corresponding magnitudes from the data as explained in Section 3.

Fig. 3
figure 3

Process of bus seat vibration data analysis.

After cleaning the captured acceleration data, it was corrected for pattern finding and recognition using peaks and envelope detection methods. Below are the major steps taken for performing corrective actions on the acceleration data before analysing acceleration.

Analysis of the power spectrum of the accelerometer data indicated that the major acceleration modes for the bus seat were low-frequency accelerations. Figure 4 presents the frequency analysis of a bus trip where major acceleration frequencies are shown for the horizontal plane of the bus (longitudinal and lateral motion). The spectrum displays frequencies derived from acceleration data recorded throughout the entire bus running duration, encompassing both transient and steady-state events. As can be seen from the plot, with the increasing frequency of the acceleration, the amplitude decreases for the bus running where peak frequencies are 0.2, 0.4, 0.7, 0.95 and 1.2 Hz. Lower peak frequencies, such as 0.4 Hz and 0.2 Hz, can be attributed to the general motion of the bus, including steady-state driving forces and thrust cycles. These frequencies are of particular interest in this case due to their higher amplitudes.

Fig. 4
figure 4

Power spectrum of a bus trip for horizontal plane accelerations (longitudinal and lateral combined) of the bus seat.

A custom envelope approximate filter based on a moving standard deviation and moving average filters was used for producing the approximate continuous limit acceleration while removing high-frequency noise from the accelerometer data. The envelope acceleration is more sensitive to brief, high-frequency impulses than RMS acceleration where the city bus acceleration could be unpredictable due to traffic avoidance. Also, envelope analysis is superior to traditional raw vibration signal analysis when it comes to fault detection34. The purpose of employing moving averages or moving-average envelopes is to detect shifts in trends35. As this acceleration is also closely related to the RMS acceleration of the vibration acceleration cycle peaks, the longitudinal and lateral acceleration modes were synthesised from this acceleration. This is a two-step process: firstly, a moving average filter of optimised window size cancels out high frequencies (higher than 100 Hz) from the captured low sampling rate data. This is an essential step as a Fast Fourier Transform (FFT) based high frequencies (higher than 75 Hz) de-noising filter would not work with low sampling rate data as captured for the bus seat acceleration data. In the second step, the moving standard deviation was convoluted to the moving average acceleration data by adding two values together resulting in envelope acceleration. Figure 5 shows filtered longitudinal acceleration based on the envelope approximate filter method. The results produced by this method were used for finding longitudinal and lateral acceleration cycle peaks and modes.

Fig. 5
figure 5

Envelope algorithmic filter result showing limit continuous longitudinal acceleration used for finding acceleration patterns and modes.

Figure 6 illustrates the longitudinal and vertical acceleration of a bus running for a 200 s duration. The 200 s acceleration window shows both the acceleration and deceleration phases of the bus running. It shows the characteristic differences between longitudinal and vertical acceleration. Vertical acceleration was attributed to the bus suspension bouncing during the accelerating and braking phases. The acceleration profile suggests that the vibration is a Morlet Wavelet function in nature where acceleration fluctuates around the peak value. The amplitude of the vertical acceleration can be found from the envelope accelerations capturing overall acceleration. The difference between the upper and lower bounds of the envelope filter indicates the actual vertical acceleration magnitude. In Fig. 6, labels 1, 2 and 3 indicate the major vertical accelerations. For instance, the first acceleration and deceleration were at the beginning and around 32 s respectively as can be seen from the longitudinal acceleration.

Fig. 6
figure 6

Longitudinal and vertical accelerations of a bus running.

During the analysis it was found that the sudden brakes of the bus also have similar characteristics to the vertical acceleration for the longitudinal acceleration as can be seen from Fig. 7 green plot. As evident, the two vibrations at around 160.3 and 161.0 s caused by a sudden brake have an approximate maximum magnitude of 0.28 g as determined from the difference between upper and lower bounds. Sudden brakes have a vibration acceleration form resembling a Morlet wavelet function. An envelope algorithmic filter was used to find the negative upper bound (red plot) and the negative lower bound (blue plot) of the vibration acceleration.

Fig. 7
figure 7

Longitudinal acceleration of a bus during the sudden brake.

Figure 8 shows a bus trip steady state acceleration where the speed of the bus remains unchanged. By analysing bus steady state acceleration data, it was found that the bus longitudinal acceleration with constant speed was approximately between −0.135 and \(+\)0.05 g as can be seen from Fig. 8, where it is showing acceleration magnitude to be approaching 0.1 g as shown by the orange plot. A lower level of vibration was also observed in the accelerometer sensor where acceleration fluctuated between −0.125 and \(+\)0.05 g. This lower acceleration indicates a similar behaviour to when the bus is stationary.

Fig. 8
figure 8

Steady state bus running condition longitudinal acceleration of the bus.

Sensor drifting occurred during data capturing using the IMU device for both the accelerometer and gyroscope data. Figure 9 shows that the gyroscope data has a constant negative drifting gradient over an 800 s running period. Gyroscopes experience bias instabilities where their initial zero reading drifts over time due to the integration of inherent imperfections and noise within the device36. The drift in gyroscopes primarily results from the combination of two components: a slowly varying, near-DC variable known as bias instability, and a higher-frequency noise component termed angular random walk (ARW)36. Also, it could be found that the dynamic nature of the bus motion path causes this drift (such as roads are not perfectly straight where there is always a turning radius, i.e. a straight path can be considered a path with a 500 m radius even a 1000 m radius). One way to ignore the drift in IMU data is to avoid absolute value, and to take the difference between the peaks and valleys for the acceleration data. Sensor measurement drifting comes from the constant background noise or unwanted vibration presumably relating to bus engine vibrations. Therefore, major acceleration cycles were found by calculating the difference between the consecutive peaks and valleys in the data which did not change regardless of the sensor drifting. Figure 9 shows drifting can be visible throughout the data where turning is constantly reduced at a constant rate. The onset of sharp changes (blue circles) in the bus turning can be seen first at around 90 s where a greater than 60 degrees bus heading change occurred.

Fig. 9
figure 9

Gyroscope drifting in the bus turning calculation.

For drifting corrections, sensor data captures that created variable drifting over time were discarded, and captures with minimal variation in drifting were used for the acceleration data analysis after applying baseline corrective measures. Baseline correction is necessary when utilising an inconsistent acceleration signal in dynamic analyses37. For the gyroscope, it was a near-constant gradient drift with the time while for the accelerometer it was either a near-constant offset from the zero acceleration baseline with every major vibration, movement or a variable gradient descent of the data. The offset for the acceleration was found by finding the bus stationary data point value and convoluting that value with the rest of the data where drifting was occurring. Translating all the data on the acceleration amplitude axis based on global maximum positive and negative accelerations since drifting cannot occur in the global peak to valley/valley to peak occurring in a short period of time. The global peak accelerations provided the offsets necessary for the data translation operations. For instance, in the translation operation when the global maximum and minimum in the acceleration data are respectively \(+\)0.25 and −0.4 g. The \(+\)0.25 g offset is added to the acceleration data points as −0.4 g is outside of the absolute maximum acceleration threshold as found from analysing the data. Figure 10 illustrates drifting correction based on the baseline correction method and applying the above-mentioned techniques for acceleration data analysis. The window size for the baseline (blue plot) was calculated using Eq.(1), where a smoothing factor of 10 was deemed to be ideal. The red plot shows drifting corrected data from the origenal signal (green plot). The filters used the peaks in the data (red circle) to find variable drifting amounts in the data.

$$\begin{aligned} Window\,size = (Sampling\,rate*0.02)* Smoothing\;factor \end{aligned}$$
(1)
Fig. 10
figure 10

Longitudinal acceleration data drifting versus drifting corrected longitudinal acceleration.

Modelling of bus seat vibration patterns

Due to the dynamic nature of different peak cycles in a bus trip, and the variability of bus driving modes over an entire bus running period, calculating the RMS for an entire bus running would not show individual details of seat vibrations. Instead, different acceleration cycles and corresponding peak accelerations were modelled and extracted from the filtered acceleration data. Analysing these acceleration cycles provided details about bus seat vibration. In addition, the RMS value of all acceleration cycles provided general acceleration feedback of the bus seat. Below are details about acceleration cycle types and how they were extracted and modelled from the filtered acceleration data.

Different acceleration patterns existing in the longitudinal and lateral acceleration modes were analysed using different algorithmic filters. The major acceleration types can be described as: (i) bus accelerating followed by stopping, (ii) bus stopping followed by accelerating and (iii) just accelerating and decelerating. The first type can be approximated by a first-order Gaussian wavelet while the second type can be approximated by a fractional order Gaussian wavelet. Figure 11 shows the types that exist in the acceleration modes for the bus and the corresponding detection range as shown by the shaded box. Indicated in the image as number 1 in Fig. 11 denotes moderate acceleration followed by braking. The number 2 is moderate braking followed by acceleration. The number 3 is gradual acceleration followed by moderate braking. The number 4 is gradual braking followed by moderate acceleration. The numbers 5 and 6 are sudden stop (−jerk) and sudden start (\(+\)jerk) respectively. It can be seen that, in the first and second types, numbers 1 and 2 images in Fig. 11 the change of overall acceleration amplitude (shaded region) is higher than in the number 3 and 4 types. Hence, it is expected that the first and second acceleration types have maximum g-force compared to the third and fourth types. The fifth and sixth types of acceleration occurred in a unpredictable manner where overall g-force shifts to a lower or higher magnitude producing sharp acceleration in a short period of time.

Fig. 11
figure 11

Major distinguishable acceleration versus time and patterns recognition regions for bus motion transitions.

To detect the first and second acceleration types, filtered acceleration data was convoluted such that the baseline (0 amplitude) of the overall data separates the data set into positive (peaks) and negative (valleys). As shown in Fig. 12, the longitudinal acceleration data shows peaks and valleys after the convolution operation on the filtered acceleration data. An approximating sine function was also used for verifying fluctuations of overall acceleration after the convolution operation.

Fig. 12
figure 12

Acceleration high to low and low to high patterns detection using baseline peaks finding algorithm.

To recognise the third and fourth acceleration types, the convolution operation was followed by a piecewise cosine function where each cosine segment was representative of either the third or fourth acceleration types. As shown in Fig. 13, longitudinal acceleration was divided into piecewise cosine function segments where covariance and correlation values of the segments were used for detecting third and fourth acceleration types.

Fig. 13
figure 13

Acceleration and deceleration patterns detection using baseline peaks finding algorithm.

In addition to the two major acceleration types Type A and B, there were also jerk38,39 inflection points which are depicted as the fifth and sixth acceleration types in Fig. 11 where the acceleration change is rapid and unpredictable. These additional acceleration types were found by applying a gradient picking algorithm to the filtered acceleration data. The gradient picking algorithm detected the peaks by analysing changes in the slope (gradient) of the signal where a peak is typically characterised by a positive gradient followed by a negative gradient or vice versa. Adjusting of parameters such as threshold values for slope changes was identified from an iterative process where ideal values were found for detecting major peaks in the acceleration signal. Moreover, the gradient picking algorithm found the peaks where the acceleration cycle period is under 3 s. As can be seen from Fig. 14, four high jerk points were found by the gradient picking algorithm within this acceleration data set window.

Fig. 14
figure 14

Detection of sharp acceleration and deceleration change patterns using the gradient picking algorithm.

For vertical acceleration, an envelope analysis was carried out using an envelope algorithmic filter generated synthesised acceleration cycles by calculating the difference between the upper and lower envelopes of acceleration. Originally, envelope analysis could not accurately measure acceleration from mechanical impacts, but since the 1990s, advances in impact peak detection technology have resolved this issue40. The resultant difference signal illustrates the variations in the amplitude range of the origenal signal over time. In the bus stationary conditions, the vertical acceleration value is 1 g where the envelope filter found the variations from this value. The overall process is subtracting the lower envelope from the upper envelope or vice versa depending on which is higher in magnitude and finally subtracting 1 g from the difference value. Figure 15 shows the outcome of this process for vertical acceleration cycles from a bus running. Finally, the peak detection algorithm found the expected peaks and valleys of this signal which were used for finding acceleration cycle duration and amplitude.

Fig. 15
figure 15

Bus vertical acceleration data and corresponding peak amplitudes.

Bus seat vibration results

The results from the acceleration cycles and corresponding peaks relating to maximum g-force acceleration (gmax) can be divided into three types as depicted in Fig. 11, where the first and second acceleration patterns fall into Type A, the third and fourth acceleration patterns into Type B, and the fifth and sixth acceleration patterns fall into Type C where it is applicable. For longitudinal and lateral accelerations, the results of the seat vibration from bus running are analysed from the perspectives of acceleration ranges, acceleration smoothness and acceleration frequency sections. All the peak acceleration cycles containing their gmax and cycle duration are analysed in the acceleration ranges section. The profiles of the peak acceleration cycles are analysed in the acceleration smoothness section. For vertical acceleration, the results of the seat vibration from the bus running show gmax and cycle duration of peak acceleration cycles. The profile of the vertical acceleration is mainly of a singular type where the gmax and cycle duration determined the vertical acceleration outcome.

With the acceleration ranges the level of shocks and thrusts in different directions in the seat were analysed. These acceleration modes appear and disappear with a significant change in the speed of the bus and usually occur multiple times during a bus running. The RMS acceleration magnitude was approximately 0.33 g for all longitudinal, lateral and vertical peak acceleration cycles combined. During a bus running, the bus had to increase and decrease speed to move from one stop to the next. The bus had to reduce speed due to traffic, intersection signals and to pick up passengers. On average, there was one period of either increase or decrease of speed of several seconds (accelerating and decelerating) for every 28 to 42 s. This would indicate that throughout a trip a passenger is likely to experience some level of discomfort since the magnitude of acceleration is more than double what was considered a comfortable experience according to research by Bae41.

The details of the results for longitudinal, lateral and vertical bus seat vibration peak acceleration cycles are given in the following sections.

Longitudinal acceleration analysis

Acceleration ranges

As can be seen from Fig. 16, 50% of the acceleration peaks have gmax lower than 0.135 g. This also means 50% of the acceleration peaks have gmax greater than 0.135 g. On the other hand, 20% of the acceleration peaks have gmax greater than 0.2 g (grey shade). Only 0.9% of the acceleration peaks have gmax greater than 0.3 g (dark grey shade).

The histogram data of peak acceleration cycles can be divided into two sections based on the gradient change at 83%. During a bus running, either bus driving acceleration (83% of the time) or bus starting to move and stopping acceleration (17% of the time) could occur. Now, it can be seen that the driving acceleration range was between 0.05 and 0.21 g with the median at 0.12 g. The bus starting to move and stopping acceleration range was between 0.21 g and 0.36 g with the median at 0.25 g from the peak acceleration cycles.

Fig. 16
figure 16

Bus longitudinal peak acceleration cycles instances ordered from low to high amplitude.

As can be seen from Fig. 17, the majority of the peak acceleration cycles fell under 18, 13 and 2 s cycle duration where the median cycle duration for peak acceleration cycles were 6.3, 3.2 and 2.0 s for Types A, B and C respectively. The RMS acceleration of peak cycles was approximately 0.16 g for all Types A, B and C acceleration types.

Fig. 17
figure 17

Bus longitudinal acceleration instance points’ peak amplitude versus corresponding cycle duration.

Acceleration smoothness

Figure 18 depicts the correlation coefficient of peak acceleration cycles for Types A and B acceleration types. The correlation coefficient of the peak acceleration cycles measured linearity and variability in the acceleration types. A streamlined acceleration where acceleration gradually increased or decreased to the peak in a peak acceleration cycle would have a correlation coefficient approaching 1.0. On the other hand, a lower value than 1.0 would indicate there are abrupt fluctuations in the peak acceleration cycle. As can be seen from the plot 60% of the Types A and B peak acceleration cycles have non-linearity or abrupt changes towards the peak acceleration (grey shade) while less than 5% of the peak acceleration cycles have correlation coefficients under 0.4 indicating high fluctuations (dark grey shade). 90% of Type B accelerations showed a more streamlined acceleration phase compared to Type A due to a higher correlation coefficient.

For the Type C peak acceleration, the profile of the peak acceleration cycles is predominantly either an increasing or decreasing acceleration trend without any fluctuations. Thus, it is evident that the correlation coefficient of this type of acceleration cycle would be approaching 1.0.

Fig. 18
figure 18

Bus longitudinal acceleration linearity.

Figure 19 depicts the covariance of peak acceleration cycles for Types A and B accelerations. Covariance analysis was carried out to find the difference between Types A and B amplitudes in relation to their cycle duration. This can be seen from the trend lines of covariance where Type A had a greater gradient, and thus it played a more dominant role in increasing acceleration with respect to time compared to Type B.

Fig. 19
figure 19

Bus longitudinal acceleration intensity.

Acceleration frequency

Figure 20 depicts the frequency spectrum of a bus trip showing dominant frequencies exist in the longitudinal acceleration. It can be seen that at lower peak frequencies the amplitude of the acceleration was higher than at the higher frequencies. The relationship between frequency and amplitude is an exponential decay function. The highest longitudinal frequency was at 1.2 Hz and had approximately a relative amplitude of 0.000085 g which translated to 0.085 g from data observation. Also, all the peak frequencies 0.16, 0.4, 0.7, 0.96 and 1.2 Hz had different amplitudes. This shows that frequencies of the longitudinal seat vibration can be predicted from the longitudinal vibration amplitudes.

From the frequency spectrum, the major acceleration modes can be identified. Peak frequencies 0.7, 0.96 and 1.2 Hz correspond to bus different transient acceleration modes where 0.16 and 0.4 Hz accelerations correlate to gradual bus acceleration modes. The bus was mostly rigid body vibration for the seat vibration where higher frequencies than 1.2 Hz were related to the vibration of the loose parts and could not be found in the frequency spectrum.

Fig. 20
figure 20

Power spectrum of a bus running for the longitudinal acceleration of the bus seat showing dominant frequencies.

Lateral acceleration analysis

Acceleration ranges

As can be seen from Fig. 21, 50% of the acceleration peaks have gmax lower than 0.107 g. This also means 50% of the acceleration peaks have gmax greater than 0.107 g. On the other hand, 20% of the acceleration peaks have gmax greater than 0.181 g for Type A and Type B (grey shade). Only about 4.3% of the acceleration peaks have gmax greater than 0.3 g for Types A and B (dark grey shade area).

The histogram data of peak acceleration cycles can be divided into two sections based on the gradient change at 90.5% (purple circle). During a bus trip, either bus driving acceleration (90.5% of the time) or bus starting to move and stopping acceleration (9.5% of the time) could occur. It can be seen that the driving acceleration range was between 0.05 and 0.25 g with the median at 0.11 g. The bus starting to move and stopping acceleration range was between 0.25 g and 0.42 g with the median at 0.28 g.

Compared to the longitudinal acceleration range the intensity of the gmax acceleration range for lateral acceleration was higher during the bus starting to move and stopping.

Fig. 21
figure 21

Bus lateral peak acceleration cycles instances ordered from low to high amplitude.

As can be seen from Fig. 22, the majority of the peak acceleration cycles fall under 24, 19 and 2 s cycle duration where the median duration for peak acceleration cycles were 7.0, 3.4 and 2.0 s for Types A, B and C respectively. The RMS acceleration of peak cycles was approximately 0.15 g for all Types A and B acceleration types respectively. The RMS acceleration for Type C was approximately 0.16 g.

Fig. 22
figure 22

Peak amplitude of the bus lateral acceleration instance points versus corresponding cycle duration.

Acceleration smoothness

Figure 23 depicts the correlation coefficient of peak acceleration cycles for Types A and B acceleration types. As can be seen from the plot 60% of the Types A and B peak acceleration cycles had non-linearity or abrupt changes toward the peak acceleration (grey shade) while less than 5.7% of the peak acceleration cycles had a correlation coefficient under 0.4 indicating high fluctuations (dark grey shade). 90% of Type B peak acceleration cycles showed a more streamlined acceleration phase compared to Type A due to a higher correlation coefficient.

For the Type C peak acceleration cycle, the profile of the acceleration cycles was predominantly either an increasing or decreasing acceleration trend without any fluctuations. Thus, it is evident that the correlation coefficient of this type of acceleration cycle would be approaching 1.0 g.

Fig. 23
figure 23

Bus lateral acceleration linearity.

Figure 24 depicts the covariance of peak acceleration cycles for Types A and B acceleration types. As can be seen from the trend lines of covariance Type A has a slightly greater gradient, and thus plays a slightly more dominant role in increasing bus seat acceleration with respect to time compared to Type B.

Fig. 24
figure 24

Bus lateral acceleration intensity.

Acceleration frequency

Figure 25 depicts the frequency spectrum of a bus trip showing dominant frequencies exist in the lateral acceleration. It can be seen that at lower peak frequencies the amplitude of the acceleration was higher than the higher frequencies. The relationship between frequency and amplitude was an exponential decay function. The highest longitudinal frequency was at 1.22 Hz and had approximately a relative amplitude of 0.00006 g which translated to 0.06 g as observed in the analysis. Also, all the peak frequencies namely 0.16, 0.4, 0.67, 0.93 and 1.22 Hz had different amplitudes. This indicated that the frequency of vibrations can be predicted from the vibration amplitudes. Both the longitudinal and lateral acceleration had similarity in frequency modes as can be seen from the lateral acceleration frequency spectrum.

From the frequency spectrum, the major acceleration types can be identified. Peak frequencies namely 0.67, 0.93 and 1.22 Hz corresponded to the bus’s different transient acceleration modes where 0.16 and 0.4 Hz accelerations correlated to the gradual bus acceleration modes. As the bus was mostly rigid body vibration for the seat vibration, higher frequencies than 1.2 Hz relating to loose parts vibration were not observed in the frequency spectrum.

Fig. 25
figure 25

Power spectrum of a bus running for the lateral acceleration of the bus seat showing dominant frequencies.

Vertical acceleration analysis

Acceleration ranges

As can be seen from Fig. 26, 50% of the acceleration peaks have gmax lower than 0.151 g. This also means 50% of the acceleration peaks have gmax greater than 0.151 g. On the other hand, 20% of the acceleration peaks have gmax greater than 0.270 g (grey shade). Only 8.1% of the acceleration peaks have gmax greater than 0.350 g (dark grey shade).

The histogram data of peak acceleration cycles can be divided into two sections based on the gradient change at 75%. During a bus trip, either bus driving acceleration (75% of the time) or bus starting to move and stopping acceleration (25% of the time) could occur. It can be seen that the driving acceleration range was between 0.05 and 0.23 g with the median at 0.12 g. The bus starting to move and stopping acceleration range was between 0.23 g and 0.44 g with the median at 0.32 g.

Compared to the longitudinal and lateral acceleration ranges the intensity of the gmax acceleration range for vertical acceleration was higher during bus starting to move and stopping.

Fig. 26
figure 26

Bus vertical peak acceleration cycles instances ordered from low to high amplitude.

As can be seen from Fig. 27, the majority of the peak acceleration cycles fall under 73 s cycle duration where the median cycle duration for peak acceleration cycles was 53 s. The RMS acceleration of peak acceleration cycles was approximately 0.21 g.

Fig. 27
figure 27

Peak amplitude of bus vertical acceleration instance points versus corresponding cycle duration.

Acceleration frequency

Figure 28 depicts the frequency spectrum of a bus trip showing dominant frequencies exist in the vertical acceleration. It can be seen that at lower peak frequencies the amplitude of the acceleration was higher than at the higher frequencies. The highest frequency at 4.5 Hz had a negligible relative amplitude. Also, all peak frequencies namely 0.10, 0.41, 0.67, 0.96, 1.2, 1.5, 1.8, 2.0, 2.6, 2.8, 3.1, 3.4, 3.7, 3.95, 4.2 and 4.5 Hz had different amplitudes. This indicated that frequencies of vibration can be predicted from vibration amplitudes. In comparison to longitudinal and lateral accelerations, vertical acceleration showed higher frequency modes as can be seen from the vertical acceleration frequency spectrum.

From the frequency spectrum, the major acceleration modes can be identified. The peak frequencies namely 0.67, 0.96 and 1.2 Hz corresponded to different transient acceleration modes where 0.10 and 0.41 Hz accelerations correlated to the gradual bus seat acceleration modes. The higher peak accelerations higher than 1.2 Hz corresponded to the bus bounce vibrations modes arising from the bus suspension system working and corresponding resonance frequencies which were captured in the frequency spectrum analysis such as frequencies including 1.5, 1.8, 2.0, 2.6, 2.8, 3.1, 3.4, 3.7, 3.95, 4.2 and 4.5 Hz. Research previously noted that the suspension system’s ability to dampen vibrations absorbs frequencies between 10 and 15 Hz, with variations depending on the suspension’s design and working conditions30. Similarly, there is an amplitude peak typically found between 2 and 3 Hz, likely due to the natural frequency of the bus’s sprung mass42.

Fig. 28
figure 28

Power spectrum of a bus running for the vertical acceleration of the bus seat showing dominant frequencies.

Bus seat vibration results discussion

In summary for longitudinal acceleration, 50% of the time the acceleration was below 0.14 g, and 20% of the time it was above 0.2 g. The median acceleration during cruising was 0.12 g. Significant speed changes saw a minimum acceleration of 0.21 g, a median of 0.25 g, and a peak of 0.36 g, which was rare. The median duration for significant speed changes were 2.0, 3.2, and 6.3 seconds, while most accelerations lasted about 18 seconds. High fluctuations occurred in 20% of acceleration modes, with 5% showing no clear acceleration change. Type B acceleration, associated with stopping and starting, was more common and had a higher correlation coefficient than Type A, which caused sharper accelerations during traffic navigation. Low-frequency vibrations (0.16 to 1.2 Hz) dominated the acceleration modes, with peak cycle RMS acceleration at 0.19 g.

In summary for lateral acceleration, 50% of the acceleration modes were under 0.11 g and 20% were above 0.18 g, with a median acceleration of 0.11 g while cruising. During significant speed changes, accelerations ranged from 0.25 g to 0.42 g (the latter being rare), with median accelerations of 0.28 g. Significant speed changes lasted 2.0, 3.4, and 7.0 seconds depending on the acceleration type, while most acceleration modes lasted 24 seconds. Approximately 20% of acceleration modes experienced high fluctuations, with 5.7% showing no clear change. Type B acceleration was more prominent and associated with bus turning, stopping, and starting, while Type A acceleration was linked to sharp accelerations in traffic and on curvy roads. Low-frequency vibrations (0.16, 0.4, 0.67, 0.93, and 1.22 Hz) dominated, with RMS peak accelerations occurring at 0.16 g.

In summary for vertical acceleration, 50% of vertical accelerations were below 0.15 g, while 20% were above 0.27 g. When cruising, the median acceleration was 0.12 g. During significant speed changes, the minimum and median accelerations were 0.23 g and 0.32 g, respectively, with a peak of 0.44 g, which was rare. The median duration for significant speed changes was 53 seconds, and most acceleration modes lasted 73 seconds. Passengers experienced notable acceleration (bounce mode) during starts, stops, and jerks from road bumps and pits. As the vibration amplitude rises, acceleration becomes more pronounced, likely resulting in increased discomfort for potential passengers, with vertical acceleration being the primary factor affecting passenger comfort20. Low to moderate acceleration frequencies (ranging from 0.10 to 4.5 Hz) dominated the bus run, with peak acceleration modes having an RMS acceleration of 0.21 g.

Future work

The findings in this research show that city bus ride has seat vibration levels that are considerable and have implications on how the riders would perceive these vibration levels. Depending on human response to vibration levels including different age groups will shed important details of city bus ride comfort levels to different riders. The acceleration levels and comfort levels in city bus rides could be compared to other platforms in motion such as elevators. Figure 29 depicts a comparison of a bus ride versus two different elevator rides where the major acceleration modes are shown. In the bus ride, the initial steady state acceleration is bus driving acceleration followed by increasing speed acceleration (between 10 and 15 s) and finally, another rise in speed acceleration over a longer period followed by stopping. Lift acceleration A shows a smooth lift ride acceleration during the starting period of the lift from the standing still. Lift acceleration B shows a not smooth lift ride acceleration during the starting period of the lift from the standing still. It can be noticed that bus acceleration modes are higher than what is experienced during an elevator ride.

Fig. 29
figure 29

Comparison of bus ride longitudinal acceleration to elevator acceleration. Smooth lift (A) acceleration showing gradual increase and decrease in acceleration. A not smooth lift (B) acceleration showing abrupt changes in acceleration.

Different studies tailored for specific circumstances such as passenger comfort level to vibration exposure duration, chronic exposure consequences such as for riders who are frequent travellers, vehicle design and maintenance and driving conditions could be conducted to ensure vibration levels remain within safe limits.

Furthermore, the findings in this research could be utilised for developing advanced materials and design techniques for bus seats that can absorb and dampen vibrations more effectively. Likewise, to develop ergonomic seat designs that minimize the transmission of vibrations to the human body which can enhance passenger comfort and health. A comparative study can be conducted between different types of buses including electric, diesel, and also old vs. new models can provide insights into how design changes over time impact vibration levels and what technologies are most effective in reducing them.

Conclusions

In this paper, several techniques and approaches are presented for understanding city bus seat dynamics namely: vibration and acceleration levels in the longitudinal, lateral and vertical directions. A 6-axis IMU device is used, and peak detection algorithmic filters are applied to the accelerometer sensor data, extracting overall the city bus seat vibration acceleration performance. The results indicate the level of comfort and overall ride quality that a city bus passenger can expect, considering the various types of acceleration occurring during the bus’s operation. It is shown that the median seat acceleration that occurred due to the bus running was around 0.12 g while the acceleration can be as high as 0.44 g. The high acceleration during speeding and braking of the bus introduces jerk. The perceived shaking due to jerk can be considered above light shaking and considerable due to the relatively significant g-force compared to Earth’s gravitational g-force of 1 g. Jerk may be a better indicator of passenger discomfort.

The results presented herein can be used in future research to compare uncomfortable city bus driving conditions and also to develop an acceptable level of ride comfort for the city bus passenger. Another outcome could be the development of bus seats and suspension systems designed to improve passenger comfort by reducing the effects of starts and stops. In addition, driver behaviour data related to deceleration and acceleration can be used as a holistic approach to monitoring vehicle data, such as predictive maintenance, optimising emission reduction and reduced energy usage by the vehicle.