Effects of the performance parameters of a wheelchair on the changes in the position of the centre of gravity of the human body in dynamic condition

Purpose The aim of this research is to establish whether, and to what extent, the tilt angle, gear ratio of the propulsion system and propulsion frequency of a wheelchair influence the position of the centre of gravity. Furthermore, it verifies the usefulness of such research using an original test stand. Materials and methods The article presents the effects of three operational parameters of a wheelchair on the position of the centre of gravity of the human body. The study included 27 wheelchair propulsion tests of a wheelchair with pushrim propulsion using the following variable parameters: gear ratio of the propulsion system, propulsion frequency and wheelchair tilt angle. The position of the centre of gravity of the human body was measured in dynamic conditions at 100 Hz. The results were represented with ellipses defining the region of variability of the position of the centre of gravity of the human body. The coordinates of the centre of gravity were measured in relation to the reference system, with the start point at the centre of the axis of rotation of the rear wheelchair wheels. The measurements were taken in a horizontal plane in relation to the base on which the test stand was positioned. Results The research carried out shows that the inclination angle of the wheelchair has the greatest influence on position of the ellipse describing the position of the centre of gravity of the human body. By controlling the change in the inclination angle value in the range from 0° to 5.4°, the standard deviation of the length of the horizontal half-axis of the ellipse (SD a) equal to 31.2 mm was obtained. For comparison, by changing the frequency of pushes (40 to 50 pushes per minute) of the wheelchair at a constant inclination angle, the standard deviation of the horizontal half-axis length (SD a) equal to 8 mm was recorded. The results of the study show a change in the position of the centre of gravity of the human body in dynamic conditions. They are relative to the contact points of the wheelchair wheels with the ground. Using the dimensions of the plotted ellipses, one can determine the values of pressure that affect the wheelchair’s individual wheels. Conclusions–Increasing the value of each aforementioned parameter resulted in the increase of strength required by the operator to propel the wheelchair. It directly influenced the position of the centre of gravity during the test.


Materials and methods
The article presents the effects of three operational parameters of a wheelchair on the position of the centre of gravity of the human body. The study included 27 wheelchair propulsion tests of a wheelchair with pushrim propulsion using the following variable parameters: gear ratio of the propulsion system, propulsion frequency and wheelchair tilt angle. The position of the centre of gravity of the human body was measured in dynamic conditions at 100 Hz. The results were represented with ellipses defining the region of variability of the position of the centre of gravity of the human body. The coordinates of the centre of gravity were measured in relation to the reference system, with the start point at the centre of the axis of rotation of the rear wheelchair wheels. The measurements were taken in a horizontal plane in relation to the base on which the test stand was positioned.

Results
The research carried out shows that the inclination angle of the wheelchair has the greatest influence on position of the ellipse describing the position of the centre of gravity of the human body. By controlling the change in the inclination angle value in the range from 0˚to 5.4˚, the standard deviation of the length of the horizontal half-axis of the ellipse (SD a) equal to 31.2 mm was obtained. For comparison, by changing the frequency of pushes (40 to 50 pushes per minute) of the wheelchair at a constant inclination angle, the standard deviation of the horizontal half-axis length (SD a) equal to 8  analyse how the wheelchair's operational parameters affect the centre of gravity, but also how it changes its position in dynamic conditions related to manual propelling. Determination of changes in the centre of gravity is important, because it allows to determine the place of application of the load resulting from human body weight. In addition, the position of the centre of gravity of the human is important when determining the stability of the wheelchair, both static and dynamic.
The study of the impact of operational parameters required the use of a specialised test stand that allows the measurement of dynamic [18] and biomechanical parameters. It also allows simultaneous simulation of various operational parameters of the wheelchair [19]. The test stand was designed and manufactured as part of the Lider VII project "Study of the biomechanics of manually propelled wheelchair for innovative manual and hybrid drives" (LIDER/ 7/0025/L-7/15/2016) financed by the National Centre for Research and Development. The aim of this research is to establish whether, and to what extent, the tilt angle, gear ratio of the propulsion system and propulsion frequency of a wheelchair influence the position of the centre of gravity. Furthermore, it verifies the usefulness of such research using an original test stand.

Materials and methods
The method of determining the position of the centre of gravity can be divided into three stages: preparation for the test, the measurement test and the analytical processing of the obtained data. As part of the preparations, the wheelchair was weighed, then along with the examined person, it was placed on the test bench. After securing the wheelchair, its position was measured to determine the exact position relative to the reference system adopted in the study. The examined person performed 30 propulsion cycles with a constant frequency (equal time intervals between pushes). Each measurement test was carried out with its unique and constant operational parameters, which include: the frequency of pushes, the inclination of the wheelchair and the gear ratio.
The wheelchair was tested using a wheelchair dynamometer [17] with a weighing scale with four strain gauges W i (Fig 1). The weighing scale was used to determine the response in the supports of the weighing scale platform on which the wheelchair and its operator were positioned. The wheelchair on the weighing scale platform can be propelled by the operator and a wheel torque is compensated by the traction rollers in the weighing scale platform.
Based on the measured response of the supports of the weighing scale platform, the position of the centre of gravity f ij on four planes π1-π4, in relation to the system with a start point at the weighing scale platform W1 (Fig 2), was determined. Based on the position, the coordinates of the position of the centre of gravity on one plane, parallel to the floor, were determined using Eqs (1) and (2).
The centre of gravity was determined every 0.01 seconds in each test, giving a set of points. The set of points representing the position of the centre of gravity of the human body in time was replaced with an ellipse (Fig 3), with the central point in the plot being an average of the measured coordinates of the position of the centre of gravity, and the inclination angle α of the  semi-major axis of the ellipse corresponded to the slope of a trend line for the measured points.
The dimensions of the ellipse that replaced the set of points of the measured centre of gravity of the human body were determined in such a way that at least 75% of the measured points were in the area drawn by the ellipse. Knowing that the measured coordinates of the centre of gravity need not be characterized by a normal distribution, the Chebyshev inequality was used in which for k = 2 a maximum of 25% of the measured points was rejected. On this basis, the dimensions of the semi-axes of the ellipse a and b were determined by calculating the standard deviation σ x and σ y multiplied by the adopted factor k according to Eqs (3) and (4).
A test stand for simulating operating conditions and measuring dynamic parameters of a wheelchair [17] was used in the tests (Fig 4). The basic component of the test stand is a weighing scale platform ( Fig 4A) with four strain gauges ( Fig 4B). The rear drive wheels are supported on the traction rollers ( Fig 4D) coupled with two electric motors (Fig 4C) simulating the resistance to motion. The front section of the wheelchair is connected to an arm tilting the wheelchair (Fig 4F). Kinematic parameters were measured using two incremental encoders ( Fig 4E) coupled with the rear drive wheels. A manually propelled wheelchair with multi-gear transmission (3 gears) was tested [18] ( Fig 5). The wheelchair was propelled by the operator using three constant propulsion frequencies at three different wheelchair tilt angles and three gear ratios (Table 1). Propulsion frequency was controlled by a metronome. The operator's task was to execute individual propulsion phases according to the rhythm set by this device. The drive wheels were loaded with the same constant anti-torque for all tests. A total of 27 tests were carried out. During the test, the position of the centre of gravity of the wheelchair was measured in constant time intervals.
The results were compared with only one of the operational parameters of the wheelchair variable. For example, changes in the propulsion frequency were measured at a constant wheelchair tilt angle and constant gear ratio. The result of the test procedure was 27 tests, for which the region of variability of the position of the centre of gravity of the human body was determined in relation to the gear ratio, propulsion frequency and wheelchair tilt angle. To recognise the defined goal as accomplished, it took one person, because the research was qualitative not quantitative. The key task was to establish whether the change in the constructed devices (research stand and drive transmission with epicyclic gear) results in a measurable change of the position of the centre of gravity. Currently, the relations between the quantitative factors of separate parameters are being examined, such as muscular effort and the position of the centre of gravity. It is, however, beyond the scope of this research, which has a finite character. To enable statistical analysis of the collected results, it was decided to use a large number of propulsion cycles.
One female subject, aged 25 years, with a body weight of 61 kg and the BMI value of 20,4 kg/m 2 was examined. The examined person is a student at the University of Physical  Wieczorek B. The authors obtained written consent of the examined person for the publication of research results with her participation. The data was presented in such a way as to ensure her complete anonymity. The set of registered results describes only one patient. The analysis does not concern a larger number of people, and therefore they are not categorised.
On this basis, it should be stated that there are no ethical concerns.

Results
During the tests, the effects of the selected operational parameters of the wheelchair on the region of variability of the position of the centre of gravity of the human body in the sagittal plane was determined. Variability of the position of the centre of gravity was defined with ellipses covering 80-90% of all measured points of the position of the centre of gravity of the human body for a single test, in which the wheelchair operator carried out 30 propulsion cycles. The test aimed to analyse the effects of the following parameters: • gear ratio at constant propulsion frequency and constant wheelchair tilt angle; • propulsion frequency at constant gear ratio and constant wheelchair tilt angle; • wheelchair tilt angle at constant propulsion frequency and constant gear ratio.
The central point on the axis of rotation of the rear wheels was selected as a starting point of the system of coordinates representing the changes in the position of the centre of gravity of the human body. The ellipses were drawn in the sagittal plane and defined using three parameters: inclination angle α of the semi-major axis a in relation to the axis x, corresponding to the travel direction, semi-major axis a defining the length of the ellipse in the direction corresponding to the travel direction, and semi-major axis b defining the width of the ellipse, perpendicular to the semi-major axis a. Charts showing the differences in registered positions of the centre of gravity for different conditions (Figs 6-8). Tables 2-4 show the parameters of the ellipses describing the region of variability of the position of the centre of gravity for selected wheelchair operational parameters.

Discussion
Changes in the gear ratio of the propulsion system affected the load on the muscular system of the upper limbs and the shoulder girdle of the operator. The higher the gear ratio of the multigear transmission of the tested wheelchair, the higher the load on the muscular system. Analysis of the variability of the region of the position of the centre of gravity depending on the gear ratio for different operating conditions (Fig 6, Table 2) shows an slight increase in ellipse size with an increase in gear ratio. Semi-major axis a of the ellipse describes the variability of the position of the centre of gravity changes to the highest degree. Dimension changes for the semi-axis a are 7.6 mm while for the semi-axis b are 1.3 mm.
The change in load affects the kinematics of the upper limb propelling the pushrim [4]. The load acting on the propulsion system affects the motion of the upper limbs and torso and changes the dimensions of the region of variability of the position of the centre of gravity of the human body. In the tests, an increase in the propulsion force [19] resulted in more chaotic movements of the upper limb and elongation of the hand motion path to the front of the wheelchair [20]. It resulted in an increase in ellipse dimensions, with an increase in the gear ratio of the wheelchair's propulsion system.
Analyses of the results of the dimensions of the regions of variability of the centre of gravity at different gear ratios show that the dimensions of the ellipses differ slightly despite a growing tendency ( Table 2). An increase in semi-major axis a, in relation to semi-major axis b with an increase in gear ratio, can be observed for all operating conditions. Another variable parameter was the propulsion frequency expressed as the number of push motions per minute. The wheelchair velocity was controlled by changing the propulsion frequency. Based on the measured position of the centre of gravity of the human body, the ellipses describing the variability of the centre of gravity at three propulsion frequencies at constant wheelchair tilt angle were determined (Fig 7, Table 3). Analysis of the results shows that with an increase in the propulsion frequency, semi-major axis a of the ellipse lengthens, which translates into an increase in the region of variability of the position of the centre of gravity of the human body. The results also show that the propulsion frequency increases the scatter of the points of the position of the centre of gravity along the travel direction (axis x), measured dimensional difference for extreme propulsion frequency values was 12.7 mm. It results from  can be explained by the biomechanics of wheelchair propulsion, in which an increase in propulsion frequency translates to changes in the shift of the human body segments in the direction corresponding to the travel direction only (axis x).
In reference to the actual operating conditions, the wheelchair tilts as a result of uneven terrain. To reflect those conditions in the tests, the wheelchair's frame was tilted by lifting the front castor wheels to determine the effect of wheelchair tilt angle β on the distribution of the position of the centre of gravity of the human body at constant propulsion frequency and constant gear ratio (Fig 8, Table 4).  Table 2

. Parameters of the ellipses describing the region of variability of the position of the centre of gravity for selected operational parameters of the wheelchair and three gear ratios, where i-Gearbox ratio, test-Independent measuring test, α-Inclination of the axis of the ellipse, a-Horizontal axis of the ellipse, b-Vertical axis of the ellipse, M-Average value, SD-Standard deviation, x CG -Average position of the centre of gravity on x-axis, y CG -Average position of the centre of gravity on y-axis, ρ-The number of measured points inside the ellipse.
p1-i = 1.96 p2-i = 1.00 p3-i = 0.51   ellipses is due to an increase in torso flexion and an increase in the wheelchair tilt angle. High increase in semi-axis a dimension (of 30.7 mm), resulted from high torso flexion angles during wheelchair propulsion at high seat [22] tilt angle. Comparison of the average dimensions of ellipses shows that the biggest difference between the dimensions of the region of variability of the position of the centre of gravity of the human body can be observed for changes in the wheelchair tilt angle. For the semi-axis a the increase was 30.7 mm and for the semi-axis b 6.8 mm. Also, the propulsion frequency affects the dimensions of the ellipses. In this case, the semi-axis dimension a increased by 12.7 mm as the propulsion frequency increased. The analysis of the effect of the gear ratio on the region of variability of the position of the centre of gravity of the human body shows a minor effect of this parameter. In the case of semi-axis a the difference in dimensions was 7.6 mm and for the semi-axis b the dimensional difference was 1.3 mm.
Analysis of the changes in the dimensions of the semi-axles shows that any changes in operating parameters affect major on the dimensions of semi-axis a and the ellipse inclination angle in relation to axis x. The tested configurations show that the operating parameters affect the dimension of semi-axis b of the ellipse defining the distribution of the centre of gravity in relation to the left and right wheelchair wheel to a negligible degree.

Summary
The tests verified the effects of the operational parameters of a wheelchair on the position of the centre of gravity of the body of an operator propelling a manual wheelchair. During the test, the gear ratio of the propulsion drive, propulsion frequency and wheelchair tilt angle were variable. For all tests, an increase in aforementioned parameters resulted in any of the parameters resulted in an increase in effort required by the operator. The effort was compensated by the wheelchair operator changing their body kinematics, which translated into changes in the dimensions of the regions of variability of the position of the centre of gravity defined with ellipses. It can be assumed that an increase in wheelchair resistance to motion [23] translates β -0˚β-1.5˚β-5.4α  The tests represent the effects of operational parameters on the distribution of the position of the centre of gravity in dynamic conditions. Since the position of the centre of gravity of the wheelchair is constant, the position of the centre of gravity of the entire system is affected by the variable position of the centre of gravity of the human body. The position of the centre of gravity of the human-wheelchair system is the result of the kinematics of the human body propelling the manual wheelchair. The design of innovative and advanced wheelchairs requires thorough knowledge of the biomechanical parameters affecting the wheelchair operation. The position of the centre of gravity is a major biomechanical parameter affecting wheelchair stability [24,25]. Stability is a key safety factor. The position of the centre of gravity also affects changes in resistance to motion and thus the muscular activity of the upper limbs and shoulder girdle. Determining the position of the centre of gravity of the human body in the horizontal plane is important. It allows to define the place of load application to the wheelchair structure (the mass force that the support frame has to withstand). The innovation of the described method is to perform measurements in dynamic conditions, which allows to significantly refine the load values assumed in calculations regarding the strength of the structure and its stability. The novelty of the results obtained in the presented method is their determination based on various parameters, in dynamic conditions. The realised experiments represent a selected part of the ongoing research. In order to consider them as complete, it is necessary to measure the position of the centre of gravity also on the vertical axis. This is necessary in order to determine the dynamic stability. This measurement is carried out within the described research project, but with the use of another method not included in this article.
The presented test results allow the assessment of the change of the position of the centre of gravity depending on the movement conditions of the manually operated wheelchair. Therefore, it is possible to determine the pressure values (reaction forces) to each wheelchair wheel, among others. It is a crucial parameter which determines, for instance, the torque required to move a wheelchair. This research is unique also due to the use of an original research stand and innovative transmission of manual drive to the wheelchair wheels, which was submitted as a patent application [12,17]. Work related to the development of innovative designs of manual wheelchairs is important because the physical activity associated with the propelling of the wheelchair enriches the rehabilitation process. Wheelchairs equipped with electric drive do not have this advantage. The developed constructional solutions and the results of the conducted research remain an added value in the field of both theoretical research and practical design guidelines for wheelchairs and their drive systems. Therefore, it broadens the scope of the design of adjusting the wheelchair to individual needs of disabled people.

Conclusions
According to the research and analysis, it can be established that the tilt angle, gear ratio of the propulsion system and propulsion frequency of a wheelchair influence the position of the centre of gravity. The biggest observed change is the result of the largest wheelchair tilt angle. Most likely, it is a result of the increased effort that the operator had to put in when propelling the wheelchair. It directly influenced the position of the centre of gravity during the test. It can be concluded that the research stand used, enables to examine in a measurable way the influence of individual parameters on the position of the centre of gravity. Establishing the quantitative influence on this parameter is the aim for further research.