An Analytical Modelling of Demand for Driving Torque of a Wheelchair with Electromechanical Drive

: This study aimed at analysing the inﬂuence of the position of the centre of gravity variability and the movement velocity on the demand for a torque and work time coverage of a wheelchair with an electromechanical drive. The variable parameter in the study was the conﬁguration of the wheelchair, namely changes in the position of the batteries which changed the weight distribution. An analytical model describing the demand for torque was used in the analysis. The set of equations was introduced into the numerical calculation software. Simulations were carried out which allowed it to analyse selected parameters of the wheelchair dynamics. An increase in the torque demand was observed due to the increase in the mass of the system from 427.7 N to 533.1 N, ranging from 6.1% to 31.6% at the simulated velocity v7 = 4.2 m/s. The increase in the demand for torque due to the increase in velocity of the wheelchair from v2 = 1.05 m/s to v7 = 4.2 m/s ranged from 25.9% to 31.6% compared to the reference velocity v1 = 0.525 m/s. The centre of gravity of the wheelchair structure localization has a non-linear impact on the analysed values. At the same time, it was not possible to deﬁne its nature—this issue remains open and requires further research. of rolling resistance for wheel e 1 in the opposite direction to the of the The frictional F t opposite sign F and the of the arm of rolling resistance for larger wheel e 2 . the same the ground reaction forces F 1 and F 2 which load the individual wheel axles. The sum written in this way allows to determine the value of the necessary driving force F n , from which the demand for driving torque is then calculated. individual forces are described in more detail into physical quantities that deﬁne their value. wheelchair wheels’ axes. For increasing values of the a parameter, the load on the larger wheel axle increases. Moving the CG towards the axis of the larger wheel causes an increase in the demand for torque, as it signiﬁcantly increases the second part of Equation (2), which is F 2 f 2 · r d 2 . Placing the battery at the front part of the wheelchair results in a more even distribution of weight by adding load to the front axle, thus reducing the need for driving torque.


Introduction
Devices supporting the mobility of the disabled and infirm increase their independence [1] and psychological comfort [2], bringing great benefits to the whole society [3,4]. One of the oldest and most commonly used means of supporting the mobility of disabled people are wheelchairs. Wheelchair's structures are dynamically developed in many directions, e.g., in order to increase the range of overcoming terrain obstacles (thresholds [5], stairs [6][7][8], elevations [9], ditches characterized by large differences in height [10]), both weight and dimensions reduction [11], frame reinforcements [12], technologies supporting propulsion [13], systems enabling the transmission of driving force from hands to wheels of a wheelchair with various degrees of dysfunction [14], structural improvements increasing the safety of movement [15,16], systems facilitating manual driving of the wheelchair, e.g., gears [8,17]. The group that has been most developed in recent years consists of electric wheelchairs with the possibility of supporting the manual propulsion [18][19][20]. This technique is beneficial for the health of people with disabilities, as it provides them with physical activity in the range of their abilities resulting from the driving force demand [21]. The disadvantages of structures equipped with an electric drive that facilitates drive are increased resistance to motion [22], greater weight [23], problems with learning to control such a wheelchair [24][25][26][27], and difficulties arising in transporting such wheelchair by another means of transportation [28]. However, the benefits resulting from the possibility of independent movement, along with the possibility of physical exertion by a disabled person during such activities, are much more important than the above-mentioned disadvantages. The authors believe that the technology that uses electric drives to assist the movement of disabled people will continue to develop dynamically [29]. The developed structures will be faced with ever higher requirements concerning the limitation of the energy consumption of the driving process realized by electric drives [30,31], and the reduction of both manufacturing and operation of such wheelchairs [32][33][34]. Moreover, the mechanisms and electrical systems used in them should also be characterized by the lowest possible energy consumption [35]. An electric assisted wheelchair is, in principle, a parallel hybrid (manual-electric) drive [36]. For this reason, when designing such devices, the mechanical parameters [37] should be taken into account, such as the location of the centre of gravity [38], the weight of the structure [39], the efficiency of mechanisms [40], as well as electrical parameters such as, for example, battery capacity [39], power and efficiency of drive units [41], etc. This paper deals with the analysis of the impact of changing the position of the centre of gravity resulting from the distribution of the wheelchair's elements (primarily the battery packs) on the torque demand and energy consumption during movement. In a wheelchair, different types of wheels are used on the front and rear axles, which have different resistance to motion characteristics also at the same time depending on the load distribution [40]. The location of the batteries in wheelchairs is also important for accessibility during the charging process [42].
The aim of the presented research work was to determine the influence of mass and spatial arrangement of the elements of the electric wheelchair drive system on the demand for drive torque and on the range of its movement. The range of the wheelchairs movement was determined in the context of the time that the electric system can work for the adopted weight variant and the arrangement of the batteries during the movement at a given velocity. The set of such data can be used to more efficiently select electromechanical drives and the elements included in such systems.
The results of the work presented in this paper are a continuation of author's research on wheelchairs. Previously, the influence of changes in the position of the human body while driving a wheelchair on the position of the centre of gravity of the entire system was analysed [38,43]. It can therefore be stated that this work is rather focused on analysing changes occurring in the system as a result of its configuration, apart from the human factor.

Research Work Conception
The entirety of work carried out can be presented as a diagram which shows a sequence in which its successive stages were performed ( Figure 1). First, a system of equations describing the forces acting on a wheelchair was produced. The basis for the development of this system is the sum of the forces acting on the wheelchair, projected on the direction of movement. A wheelchair with a hybrid electric-manual drive was used for the tests and analysis. It is a certain development of the concept of an electric wheelchair, but in principle it has a mode of operation that allows it to be classified into this group. The above-mentioned steps became the basis for carrying out measurements. First, the geometry of the tested system was measured, as it was crucial to obtain a number of parameters necessary for the simulation. Then, variants of variable parameters in the simulation were defined and adopted. They were associated with different masses and the arrangement of batteries supplying the drives and control system of the wheelchair. As a result, it was possible to perform a number of measurements of the variable values-that is, the masses and the location of the centre of gravity of the entire system. A series of collected data was introduced into the numerical computing environment. As a result, simulation tests were carried out for all adopted variants. The determined values were then related to the characteristics of the drive system and power supply, which allowed for the development of a set of results. This provided the basis for compiling and comparing the effects of the calculations performed for all the adopted variants. ing the effects of the calculations performed for all the adopted variants.
A detailed description of the methodology for the implementation of the indicated stages of research work has been described in the following subchapters of the "Methods" section. For the sake of clarity, the order indicated in this subsection was followed. A complete list of all notations and physical quantities used in this paper can be found in Appendix A in Table A1. Diagram illustrating the course of research works; I-analytical development of the force system acting on the wheelchair, II-development of a wheelchair with a hybrid drive, III-adaptation of simulation variants for the purpose of measuring variable parameters at the test stand, IV-measurements of constant parameters on the test stand, V-measurements of variable parameters on the test stand, VI-introducing the input parameters to the simulation, VII-performing the simulation in the environment for numerical calculations for all the adopted variants, VIIIcalculation of the work time coverage based on the simulation results and the characteristics of the engine and the supply system, IX-summary of the results for all variants.
The novelty of the presented research process results from the use of the mathematical model used classically to describe the movement of a motor vehicle to describe the movement of a wheelchair. It required introducing changes resulting from differences between these systems. In the case of a wheelchair, the location of the centre of gravity is of particular importance, which results from both the human factor (analysed in the work [38] as well as the configuration of the wheelchair itself. In the presented study, this matter was analysed on the example of an electric wheelchair with a different arrangement of the batteries which altered both the weight of the system as well as its concentration in relation to the wheel's axles. The high significance of the position in the case of a wheelchair results from a much lower weight ratio of a wheelchair and a human being compared to a motor vehicle.
Limitations of the research result mainly from the elements used in the adopted method. The results of the measurements refer to specific variants of the wheelchair resulting from its dimensions, the choice of the type of battery and electric drives and its placement. Although it should be expected that the nature of the observed changes will be similar, on the basis of the presented data it is impossible to directly relate them to, e.g., a wheelchair equipped with a different type of electric drive. In turn, the advantage of this approach is that the developed mathematical model is so general that its implementation in different cases remains possible after minor changes or modifications are made. In its present form, the model describes the movement of a wheelchair in a straight line.

Analytical Basics of Modelling
The distribution of forces acting on the wheelchair along with the user, projected into the direction of movement, is shown in Figure 2. The presented mathematical model is based on the description of the sum of forces acting on the wheelchair during its movement on a hard, flat surface at a given angle . From a mechanical point of view, the system is considered to be flat, therefore the forces are applied at the points of operation on a two-dimensional plane. The inertia force of the system is applied in the centre of Figure 1. Diagram illustrating the course of research works; I-analytical development of the force system acting on the wheelchair, II-development of a wheelchair with a hybrid drive, III-adaptation of simulation variants for the purpose of measuring variable parameters at the test stand, IV-measurements of constant parameters on the test stand, V-measurements of variable parameters on the test stand, VI-introducing the input parameters to the simulation, VII-performing the simulation in the environment for numerical calculations for all the adopted variants, VIII-calculation of the work time coverage based on the simulation results and the characteristics of the engine and the supply system, IX-summary of the results for all variants.
A detailed description of the methodology for the implementation of the indicated stages of research work has been described in the following subchapters of the "Methods" section. For the sake of clarity, the order indicated in this subsection was followed. A complete list of all notations and physical quantities used in this paper can be found in Appendix A in Table A1.
The novelty of the presented research process results from the use of the mathematical model used classically to describe the movement of a motor vehicle to describe the movement of a wheelchair. It required introducing changes resulting from differences between these systems. In the case of a wheelchair, the location of the centre of gravity is of particular importance, which results from both the human factor (analysed in the work [38] as well as the configuration of the wheelchair itself. In the presented study, this matter was analysed on the example of an electric wheelchair with a different arrangement of the batteries which altered both the weight of the system as well as its concentration in relation to the wheel's axles. The high significance of the COG position in the case of a wheelchair results from a much lower weight ratio of a wheelchair and a human being compared to a motor vehicle.
Limitations of the research result mainly from the elements used in the adopted method. The results of the measurements refer to specific variants of the wheelchair resulting from its dimensions, the choice of the type of battery and electric drives and its placement. Although it should be expected that the nature of the observed changes will be similar, on the basis of the presented data it is impossible to directly relate them to, e.g., a wheelchair equipped with a different type of electric drive. In turn, the advantage of this approach is that the developed mathematical model is so general that its implementation in different cases remains possible after minor changes or modifications are made. In its present form, the model describes the movement of a wheelchair in a straight line.

Analytical Basics of Modelling
The distribution of forces acting on the wheelchair along with the user, projected into the direction of movement, is shown in Figure 2. The presented mathematical model is based on the description of the sum of forces acting on the wheelchair during its movement on a hard, flat surface at a given angle α. From a mechanical point of view, the system is considered to be flat, therefore the forces are applied at the points of operation on a two-dimensional plane. The inertia force of the system F b is applied in the centre of gravity of the system at point W, as is the gravity force (weight of the system) G. The aerodynamic drag force F a was applied in the centre of the system height (half of the height) at point V. The frictional force of the smaller wheel F t1 was applied at the moment arm of rolling resistance for smaller wheel e 1 in the opposite direction to the movement of the wheelchair. The frictional force of the larger wheel F t2 was summed (with the opposite sign to counteract the movement) with the propulsion force F n and applied to the moment of the arm of rolling resistance for larger wheel e 2 . In the same places, the ground reaction forces F 1 and F 2 were applied, which load the individual wheel axles. The sum written in this way allows to determine the value of the necessary driving force F n , from which the demand for driving torque is then calculated. Ultimately, individual forces are described in more detail into physical quantities that define their value. right pair of wheels. Symmetry can be disturbed by many facto sure differences in the larger wheels of a wheelchair. With that b metry assumption is acceptable because the worst-case criterion culating and selecting a drive. Knowing the demand for driving metrical) case, it is easy to determine the demand for driving torq loading of either side of the wheelchair. The advantage arising o however, the possibility of considering only half of the wheelcha down to a coplanar system (i.e., one that is contained on a two-dim significantly simplifies the calculations. As a result, the values o are determined for one of the two larger wheels of the wheelcha where is to be considered total friction force while -longit ing physical quantities are: and are the moments of inert of rotation for respectively smaller and larger wheels. The expres eration. The other values depend on the geometry of the system are the dynamic radiuses of the wheels for respectively sma and are consistently the rolling radiuses of the smaller and la The case analysis was based on the assumption of symmetry of the system. This is, of course, a certain approximation of reality, because neither the weight of the wheelchair nor the weight of the user has to be symmetrically distributed with respect to the left and right pair of wheels. Symmetry can be disturbed by many factors, including small pressure differences in the larger wheels of a wheelchair. With that being said, the ideal symmetry assumption is acceptable because the worst-case criterion is usually used when calculating and selecting a drive. Knowing the demand for driving torque in an ideal (symmetrical) case, it is easy to determine the demand for driving torque in the event of uneven loading of either side of the wheelchair. The advantage arising of such an assumption is, however, the possibility of considering only half of the wheelchair, which basically comes down to a coplanar system (i.e., one that is contained on a two-dimensional plane)-which significantly simplifies the calculations. As a result, the values of the demand for torque are determined for one of the two larger wheels of the wheelchair.
The sum of the projections of forces on the wheelchair travel axis can be written as: .. ..
where F t is to be considered total friction force while F w -longitudinal force. The following physical quantities are: J k1 and J k2 are the moments of inertia with respect to the axis of x(t) is linear acceleration. The other values depend on the geometry of the system are defined as r d1 and r d2 are the dynamic radiuses of the wheels for respectively smaller and larger ones; r t1 and r t2 are consistently the rolling radiuses of the smaller and larger wheels.
This expression, after describing the component forces in detail and being multiplied by the dynamic radius of the larger wheel, allows to obtain the dependence on the torque T n necessary to drive the wheelchair under given conditions: where f 1 is rolling resistance coefficient of the smaller wheel, f 2 is rolling resistance coefficient of the larger wheel, G is weight of the system, α is inclination angle, A p is the surface area of the front face of the system, C x is a dimensionless factor of aerodynamic drag, . ϕ(t) is angular velocity, and .. ϕ(t) is angular acceleration, g is acceleration due to gravity. In Equation (2) the value F 1 represents a load on the smaller wheel axle and F 2 -on the larger wheel axle which can be described as: . .
where: d is vertical position of system's centre of height, h is vertical position system's centre of gravity, a is the distance from the axle of the smaller wheel to the centre of gravity, and l is the distance between the axles of the wheelchair's wheels. Detailed derivation of the model described by Equations (1)-(3) can be found in the previous publication of the authors [44].
The demand for the mechanical power P M that must be met by the wheelchair drive is described by Equation (4): where T n is the driving torque determined according to Equation (2) based on previous assumptions, n is rotational speed in (rpm) of the wheelchair's larger wheel. If into Equation (4) the value of T n will be substituted in (Nm) and the value of n in (rpm), then the resulting unit for P M will be (kW), which is favourable for further calculations. Assuming that the drive will be utilizing a DC motor, the electric power will be defined by the formula (5): where U is the voltage and I stands for electric current. Both power values P M and P E are related to each other by the total efficiency of the system. The total efficiency η of a system is the product of the efficiency of its elements and mechanisms. A number of such elements can be indicated in the system, but the specification of all individual efficiency values is extremely difficult, therefore the most important ones were adopted in further considerations, and the others were simplified according to Equation (6): where η em is the efficiency of electric motor, η b is the efficiency of bearings and η r is the efficiency of the rest of the system (losses resulting from, among others, electrical resistance, heat losses and power consumption of the control system). Formulas (4) and (5) must be multiplied by time in order to receive units as kilowatt-hour, yielding: where t M is work time coverage (work time of the mechanical part), and in turn t E is work time of the electric system determined on the basis of the knowledge of the ampere-hours parameter of the batteries. Substituting (4) and (5) to (7) and transforming those equations allows for obtaining the relationship (8) which relates the quantity that defines the battery capacity, the operating parameters of electric motors, efficiency, but also the demand for driving torque: The value of T n in turn is dependent, as Equation (2) shows, on a number of parameters resulting from the geometry and weight of the wheelchair. In this study, the variables (i.e., the analysed parameters) were the mass of the system and its distribution, and thus also the resulting displacement of the centre of gravity of the wheelchair itself.
Substituting successive values of mass m and the location of the center of gravity (defined by parameter a), and for various cases, it is possible to determine the operating time of the system (work time coverage) t M , which results not only from the capacity of the batteries used, but also from their location in relation to the plane of the wheelchair movement.

Hybrid Electric-Manual Wheelchair
A wheelchair with a hybrid electric-manual drive was used in the research. It was created by equipping the classic Breezy 90 wheelchair (Sunrise Medial S. L.) with the MagicPie 5 electric drive units by Goldenmotor. The invention has a control system that allows it to work in various modes-it can use electric drives to completely propel the system, but also only to support the manual drive carried out by the user. The system is equipped with control levers, manual brakes and an LCD screen so as to simplify its operation as much as possible from the user's point of view. The whole is powered by a Li-ION battery, which has been permanently attached to the frame at the rear of the wheelchair. The developed device is shown in Figure 3a,b. It was created on the basis of the PAT 2102 patent; a detailed description of this wheelchair can be found in [45]. The developed drive system of the presented device constitutes an electro-mechanical drive in its entirety. equations allows for obtaining the relationship (8) which relates the quantity that the battery capacity, the operating parameters of electric motors, efficiency, but a demand for driving torque: The value of in turn is dependent, as Equation (2) shows, on a number of eters resulting from the geometry and weight of the wheelchair. In this study, the bles (i.e., the analysed parameters) were the mass of the system and its distributio thus also the resulting displacement of the centre of gravity of the wheelchair itsel Substituting successive values of mass and the location of the center of (defined by parameter ), and for various cases, it is possible to determine the op time of the system (work time coverage) , which results not only from the cap the batteries used, but also from their location in relation to the plane of the whe movement.

Hybrid Electric-Manual Wheelchair
A wheelchair with a hybrid electric-manual drive was used in the research created by equipping the classic Breezy 90 wheelchair (Sunrise Medial S. L.) w MagicPie 5 electric drive units by Goldenmotor. The invention has a control syste allows it to work in various modes-it can use electric drives to completely pro system, but also only to support the manual drive carried out by the user. The sy equipped with control levers, manual brakes and an LCD screen so as to simplify eration as much as possible from the user's point of view. The whole is powered b ION battery, which has been permanently attached to the frame at the rear of the chair. The developed device is shown in Figure 3a,b. It was created on the basis of t 2102 patent; a detailed description of this wheelchair can be found in [45]. The dev drive system of the presented device constitutes an electro-mechanical drive in its e

Test and Simulation Variants
In order to carry out calculations and simulations, nine variants of the whe configuration were adopted -. They were mainly created by different mass d tions, in the form of batteries within an electric wheelchair. The diagram of the ge of the adopted variants is shown in Figure 4. The drawings are made to scale, the sions of symmetrically placed elements are not repeated, and the defined dimensi not given again, unless they have changed. The dimensions were drawn with resp fixed point, which was always in the same place and was the reference for the coo system (the origin of the coordinate system), which was marked in blue. Referring

Test and Simulation Variants
In order to carry out calculations and simulations, nine variants of the wheelchair configuration were adopted A 1 -A 9 . They were mainly created by different mass distributions, in the form of batteries within an electric wheelchair. The diagram of the geometry of the adopted variants is shown in Figure 4. The drawings are made to scale, the dimensions of symmetrically placed elements are not repeated, and the defined dimensions are not given again, unless they have changed. The dimensions were drawn with respect to a fixed point, which was always in the same place and was the reference for the coordinate system (the origin of the coordinate system), which was marked in blue. Referring to the coordinate system introduced in Figure 2, the axis in the direction of movement is the x axis and the second axis in the plane of movement is the z axis, since the y axis remains upward (see Figure 2a). Introduced axis designations of the Cartesian coordinate system are consistently preserved throughout the text. The first variant A 1 is a classic Breezy 90 wheelchair, while the next variant A 2 is a modified wheelchair, equipped with MagicPie 5 electric drives. The variant A 3 is a modified wheelchair, with a battery in the rear part, on the left-hand side from the user's perspective. This is the case of the actual system described in Section 2.3. The variant A 4 it has a battery located centrally under the seat. A 5 and A 6 have two symmetrically placed batteries, the first of which is located at the rear of the wheelchair, and the second under the seat. The variants form A 7 to A 9 are modifications of A 5 created by adding weights simulating batteries successively by adding one in the front part of the wheelchair, adding two (one after the other) behind the seat and finally adding two symmetrically placed on the sides of the stroller under the seat. The variants A 1 and A 2 were adopted for the purpose of comparison and experimental determination of the mass of electric drives and batteries. The small dimensional deviations in relation to the commercially available Breezy 90 wheelchair result from modifications made to its structure in order to attach the cooperating elements constituting the invention.
Energies 2021, 14, x FOR PEER REVIEW 7 of 20 coordinate system introduced in Figure 2, the axis in the direction of movement is the axis and the second axis in the plane of movement is the axis, since the axis remains upward (see Figure 2a). Introduced axis designations of the Cartesian coordinate system are consistently preserved throughout the text. The first variant is a classic Breezy 90 wheelchair, while the next variant is a modified wheelchair, equipped with MagicPie 5 electric drives. The variant is a modified wheelchair, with a battery in the rear part, on the left-hand side from the user's perspective. This is the case of the actual system described in Section 2.3. The variant it has a battery located centrally under the seat. and have two symmetrically placed batteries, the first of which is located at the rear of the wheelchair, and the second under the seat. The variants form to are modifications of created by adding weights simulating batteries successively by adding one in the front part of the wheelchair, adding two (one after the other) behind the seat and finally adding two symmetrically placed on the sides of the stroller under the seat. The variants i were adopted for the purpose of comparison and experimental determination of the mass of electric drives and batteries. The small dimensional deviations in relation to the commercially available Breezy 90 wheelchair result from modifications made to its structure in order to attach the cooperating elements constituting the invention. --designations of individual variants, -wheelchair outline; red points represent the geometric centre of the battery, blue points represent the origin of the coordinate system, light blue points represent the geometric centre of the rear wheels.  Before the test, the larger (pneumatic) wheels of the wheelchair were inflated to the nominal pressure value of 3 ± 0.1 bar suggested by the manufacturer. The configurations from A 2 to A 9 used the same wheelchair; the position and amount of additional mass were changed without changing its position and other parameters, including tire pressure. Therefore, the pressure in the pneumatic wheels did not affect the position of the centre of gravity in the test conducted on dynamometer (detailed dynamometer description is presented in the next section).

Constant Parameters of the Simulation
The values of the geometrical and physical quantities of the system were either adopted on the basis of the literature or measured directly, as it is shown in Table 1. The geometry was determined during measurements made on a wheelchair immobilized on a test stand (dynamometer).

Variable Parameters of the Simulation
Conducting the simulation required the determination of the mass and the location of the centre of gravity of a wheelchair that was not modified and equipped with electric drives and a control system with batteries. For this purpose, a wheelchair dynamometer was used. It is equipped with a measurement system that uses the HBM PMX data acquisition system. Mass measurements are carried out thanks to the Z6FC6 100 kg transducers, also by HBM. The measurement consists of registering the weight of the weighing pan with the wheelchair in four support places R 1 − R 4 , and its general idea is presented in Figure 5. The knowledge of systems dimensions allows for determining the location of centre of gravity (CG) at the intersection of lines drawn form points P 1 − P 4 . Characteristic dimensions k 12 , k 14 , k 23 , k 43 that define the points P 1 − P 4 are determined using statics equations for individual planes of the system. A detailed description of the principle of determining the centre of gravity according to this method can be found in [43]. The values of dimensions k 01 and k 02 define the location of strain gauges in relation to the adopted coordinate system (as shown on Figure 4).
As a result of the measurements, a cloud of 2500 points was obtained defining the location of the CG of the wheelchair for a given variant with respect to the origin of the x coordinate system as shown in Figure 5. The masses of the wheelchair in individual configurations were also recorded A 1 − A 9 . For the estimator of the searched values (the location of CG on x axis and the mass of the system for a given variant) the arithmetic mean was adopted, and the estimator error was the standard deviation of the arithmetic mean.
Determining the distance from the axle of the smaller whee for individual cases equals to the difference in dimensions g variant , that being the value of 810 mm and the result of the p ment on axis.

Numerical Simulation
The developed system of Equations (2), (3) and (8) was in Simulink numerical computing environment (MAthWorkks, v.R were conducted with a time step of Δ = 0.005 s, assuming a kn the sake of simplicity, it has been assumed that the wheelchair velocity with a constant acceleration for the duration of 10 s, (so dynamic loads). After reaching a set value, it will be maintained ulation, as shown in Figure 6. Determining the distance from the axle of the smaller wheel to the centre of gravity a for individual cases equals to the difference in dimensions given in Figure 4 for the variant A 1 , that being the value of 810 mm and the result of the position of CG measurement on x axis.

Numerical Simulation
The developed system of Equations (2), (3) and (8) was introduced to the MatLab Simulink numerical computing environment (MAthWorkks, v.R2015a). The simulations were conducted with a time step of ∆t = 0.005 s, assuming a known velocity profile. For the sake of simplicity, it has been assumed that the wheelchair will accelerate to the set velocity with a constant acceleration for the duration of 10 s, (so as not to cause too high dynamic loads). After reaching a set value, it will be maintained until the end of the simulation, as shown in Figure 6. Determining the distance from the axle of the smaller wheel to the centre of gravity for individual cases equals to the difference in dimensions given in Figure 4 for the variant , that being the value of 810 mm and the result of the position of measurement on axis.

Numerical Simulation
The developed system of Equations (2), (3) and (8) was introduced to the MatLab Simulink numerical computing environment (MAthWorkks, v.R2015a). The simulations were conducted with a time step of Δ = 0.005 s, assuming a known velocity profile. For the sake of simplicity, it has been assumed that the wheelchair will accelerate to the set velocity with a constant acceleration for the duration of 10 s, (so as not to cause too high dynamic loads). After reaching a set value, it will be maintained until the end of the simulation, as shown in Figure 6.

Characteristics of the Electric Drive
The parameters related to the motors driving the wheelchair and electric energy batteries were read directly from the data provided by the manufacturer of the applied solutions. The graph in Figure 7 shows the performance of the Golden Motor MagicPie 5 system.
Energies 2021, 14, x FOR PEER REVIEW

Characteristics of the Electric Drive
The parameters related to the motors driving the wheelchair and elect teries were read directly from the data provided by the manufacturer of th tions. The graph in Figure 7 shows the performance of the Golden Motor M tem. in Equation (8). Weights simulating auxiliary b similar mass (determined in detail in the further part of the study). The pr pendence of capacity to battery mass was assumed, therefore the use of aux assumes a proportional increase in value of ⋅ . It is worth noting here th ship between the mass of the battery and its capacity in general is not nec and depends on a number of factors, including its type. However, due to th ences in masses shown later, the assumption of a linear relationship does n significant error in the calculations. The introduction of the results obtain tions (2) and (3) and the values read from the characteristics of the used allows for determination the working time of the mechanical system u (8). It essentially is the work time coverage, which directly translates int travel range (in the context of the maximum possible distance achievabl conditions for a given velocity).
In order to perform numerical calculations, the efficiency curve of the was approximated by (9) according to the dependency:  For the calculated values of the torque demand T n of the wheelchair, depending on the simulation variant, appropriate values of the supply current I, rotational speed n, and efficiency η em were assumed for further calculations. The battery [54] used in the solution had the nominal voltage of U = 36 V. Its capacity equals to 13 Ah and it corresponds to the term: I·t E in Equation (8). Weights simulating auxiliary batteries had a similar mass (determined in detail in the further part of the study). The proportional dependence of capacity to battery mass was assumed, therefore the use of auxiliary batteries assumes a proportional increase in value of I·t E . It is worth noting here that the relationship between the mass of the battery and its capacity in general is not necessarily linear and depends on a number of factors, including its type. However, due to the slight differences in masses shown later, the assumption of a linear relationship does not introduce a significant error in the calculations. The introduction of the results obtained from Equations (2) and (3) and the values read from the characteristics of the used electric drive allows for determination the working time of the mechanical system t M using Equation (8). It essentially is the work time coverage, which directly translates into the vehicle's travel range (in the context of the maximum possible distance achievable under given conditions for a given velocity).
In order to perform numerical calculations, the efficiency curve of the electric motor was approximated by (9) according to the dependency: with a 1 − a 4 being the parameters of this equation that were calculated as values: a 1 = −3.84·10 −1 , a 2 = 4.95·10, a 3 = 4.49·10 −1 and a 4 = 8.48·10 −3 . For the adopted values, the correlation coefficient was 0.99.

Centre of Gravity and Mass Measurements Results
The measurement results of the location of CG for individual variants are presented in Figure 8, the specific values are collected in Table 2. It also presents the determined values of the weight G and distance from the axle of the smaller wheel to the centre of gravity a for individual cases. They were adopted for further calculations. The difference in the masses of selected variants makes it possible to easily determine the masses of the system components. Thus, the weight of the wheelchair itself is m 0 = 17.55 kg (for traditional and standard wheelchair, such mass should be considered typical, it is consistent with the results of other researchers [55,56]), the mass of the drives with accessories and mounting is m 1 = 26.0 kg, the mass of a single battery is m 2 = 2.66 kg (this is typical mass for batteries in such a range of solutions, other researchers report it in the range between 0.7 kg and 7 kg [57,58]), and the averaged mass of a single weight is m 3 = 2.72 kg. Therefore, considering the earlier assumption of a linear mass to capacity ratio of the battery, auxiliary batteries capacity equals to I·t E = 13.3 Ah (compared to nominal batteries capacity I·t E = 13 Ah).
Energies 2021, 14, x FOR PEER REVIEW 11 of 2 system components. Thus, the weight of the wheelchair itself is = 17.55 kg (for tradi tional and standard wheelchair, such mass should be considered typical, it is consisten with the results of other researchers [55,56]), the mass of the drives with accessories and mounting is = 26.0 kg, the mass of a single battery is = 2.66 kg (this is typical mas for batteries in such a range of solutions, other researchers report it in the range between 0.7 kg and 7 kg [57,58]), and the averaged mass of a single weight is = 2.72 kg. There fore, considering the earlier assumption of a linear mass to capacity ratio of the battery auxiliary batteries capacity equals to ⋅ = 13.3 Ah (compared to nominal batteries ca pacity ⋅ = 13 Ah). The analysis of the results shows that the of the wheelchair is not ideally located in its geometric axis of symmetry. This may be due to a number of factors. The manufac turing accuracy of the individual elements of the system is finite, so obtaining perfectl equal masses and symmetry of their mounting is not achievable [59]. In any case, thi parameter does not affect the simulation results because the position of the on the axis is not analyzed in this paper. At this point however, it is worth noting that, in genera the position of the on the axis affects the load on the individual axles of wheelchai [60] and may be an important parameter in some cases-in particular during curvilinea motion [61,62].   The analysis of the results shows that the CG of the wheelchair is not ideally located in its geometric axis of symmetry. This may be due to a number of factors. The manufacturing accuracy of the individual elements of the system is finite, so obtaining perfectly equal masses and symmetry of their mounting is not achievable [59]. In any case, this parameter does not affect the simulation results because the position of the CG on the z axis is not analyzed in this paper. At this point however, it is worth noting that, in general, the position of the CG on the z axis affects the load on the individual axles of wheelchair [60] and may be an important parameter in some cases-in particular during curvilinear motion [61,62].

Work Time Coverage and Torque Demand
The values obtained as a result of the simulation are presented in Table 3. The calculation results related to the determination of the work time coverage value as a function of the wheelchair's velocity and the total system weight are shown in Figure 9a. As it results from the analysis of the presented graph, with the increase in velocity, the possible working time of the system decreases because the electric drive requires greater quantities of electricity to achieve and maintain it. The changes caused by the weight alteration are of the opposite nature, and also occur with much larger increments. This seemingly illogical effect is the result of the fact that the mass added to the system by additional batteries results in an increase in their capacity. Apparently, the increase in capacity progresses more rapidly than the increase in torque demand caused by the extra weight. This parameter however has an upper limit resulting mainly from the available space that the batteries can occupy. It is also indirectly related to the strength (in terms of load capacity) of the wheelchair frame structure. Table 3. Values obtained in course of conduction the simulation.

Velocity
Variant The effect of weight on the tested parameters is approximately linear, while the effect of velocity is a quadratic function, as can be seen more clearly in Figure 9b. On the other hand, the position of the centre of gravity on the x axis, represented by the parameter a, remains in a much more complicated relation to the demand for driving torque, as shown by the plane shown in Figure 10. The value of the parameter a seems to have a non-linear effect on the demand for driving torque [63]. The location of the centre of gravity depends on the mass, but also on the geometry of the entire system [64], and this relationship in the general case remains a complicated function resulting from the location of the centre of mass in three dimensions. The non-linear relationship between the tested parameters and the value of the parameter a is shown in Figure 11a,b. Analytical determination of the equations determining the exact relationship between the position of the centre of gravity and the distribution of mass remains theoretically possible. However, performing such analysis for a system with so many elements would be inefficient as it would be extremely time-consuming process. In practice, the position of the CG of complex systems (such as vehicles) is determined experimentally, for example, by means of strain gauges [65,66], similar to the methodology described in Section 2.6, or the measurement of the moment of inertia [67]. An unambiguous and direct determination of impact nature of the parameter a on the studied quantities is problematic and needs more study. analysis for a system with so many elements would be inefficient as it would be extremely time-consuming process. In practice, the position of the of complex systems (such as vehicles) is determined experimentally, for example, by means of strain gauges [65,66], similar to the methodology described in Section 2.6, or the measurement of the moment of inertia [67]. An unambiguous and direct determination of impact nature of the parameter on the studied quantities is problematic and needs more study.   of inertia [67]. An unambiguous and direct determination o eter on the studied quantities is problematic and needs m      Next, a comparison of the position of the system's additional mass in relation to the geometrical centre of the wheelchair was analysed. Thus, two groups of configurations were determined. The "Front" configuration including variants: A 4 , A 6 , A 8 that have the greater part of the battery weight shifted towards the front of the wheelchair (towards increasing values on the x axis). On the other hand, configuration variants "Rear", being: A 3 , A 5 , and A 9 were determined as the opposite. A similar arrangement was analysed for motor vehicles in [68]. The additional mass was concentrated in the direction of decreasing quantities on the axis according to the adopted coordinate system, therefore they had a greater value of the parameter a because this quantity is measured from the axis of rotation of the smaller wheel. These configurations properly arranged in pairs (as shown in Figure 12) have similar weights; the slight difference is due to measurement inaccuracy. The analysis of this diagram shows that the shift of the CG towards the rear wheels increases the torque demand. Most likely, this is due to the uneven load on the wheelchair wheels' axes. For increasing values of the a parameter, the load on the larger wheel axle increases. Moving the CG towards the axis of the larger wheel causes an increase in the demand for torque, as it significantly increases the second part of Equation (2), which is F 2 f 2 ·r d2 . Placing the battery at the front part of the wheelchair results in a more even distribution of weight by adding load to the front axle, thus reducing the need for driving torque.
gies 2021, 14, x FOR PEER REVIEW Next, a comparison of the position of the system's addition geometrical centre of the wheelchair was analysed. Thus, two were determined. The "Front" configuration including variants: greater part of the battery weight shifted towards the front of increasing values on the axis). On the other hand, configuratio , , and were determined as the opposite. A similar ar for motor vehicles in [68]. The additional mass was concentrat creasing quantities on the axis according to the adopted coordin had a greater value of the parameter because this quantity is of rotation of the smaller wheel. These configurations proper shown in Figure 12) have similar weights; the slight differenc inaccuracy. The analysis of this diagram shows that the shift of wheels increases the torque demand. Most likely, this is due to wheelchair wheels' axes. For increasing values of the parame wheel axle increases. Moving the towards the axis of the la crease in the demand for torque, as it significantly increases the (2), which is ⋅ . Placing the battery at the front part of th more even distribution of weight by adding load to the front axl for driving torque.   In order to show the scale of the changes in the value of torque demand and work time coverage, the proportions were arranged according to Equation (10): where x 1 is the quantity after the change and x 0 is the quantity before the change (reference quantity). For the assumed relationship, increases will be positive, and decreases negative. For the analysis of the torque demand, the reference value was the value for a wheelchair with electric drives but without a battery (A 2 variant for velocity of v1), while for the analysis of work time coverage, the value for a wheelchair with one battery in the "Rear" position was adopted (A 2 variant for velocity of v1). The set of obtained results is presented in Figures 13 and 14.
Energies 2021, 14, x FOR PEER REVIEW sign, because in the relation to the variant, the addition of auxiliary batter ficiently low simulated wheelchair velocity results in a large increase in their ity. For those cases it results in an increase in work time coverage.

•
As the weight of the system increases, the demand for drive torque inc work time coverage is increasing, with an approximately linear trend; • As the weight of the system increases, so does the work time coverage increase in driving torque demand; • The increase in the simulated velocity causes changes in the demand torque and work time coverage, the change being approximately a squar • The set of recorded data does not allow for an unambiguous and direct det

•
As the weight of the system increases, the demand for drive torque in work time coverage is increasing, with an approximately linear trend; • As the weight of the system increases, so does the work time coverage increase in driving torque demand; • The increase in the simulated velocity causes changes in the demand torque and work time coverage, the change being approximately a squar • The set of recorded data does not allow for an unambiguous and direct de of impact nature of the parameter value (distance from the axle of The analysis of the presented charts confirms the earlier observations. The increase in velocity and at the same time, the rise of the system weight along with its shift towards the larger wheels of the wheelchair increases the demand for drive torque. Simultaneously, adding auxiliary batteries to the system (although it increases the overall mass and affects the CG), for selected ranges of variability of said parameters, results in an increase in work time coverage, effectively increasing the battery capacity. Conversely, an increase in the wheelchair velocity causes a decrease in the maximum capacity of the battery as it results in an increased demand for power. Part of the graph's bars (on Figure 14) change sign, because in the relation to the A 2 variant, the addition of auxiliary batteries at a sufficiently low simulated wheelchair velocity results in a large increase in their total capacity. For those cases it results in an increase in work time coverage.

•
As the weight of the system increases, the demand for drive torque increases and work time coverage is increasing, with an approximately linear trend; • As the weight of the system increases, so does the work time coverage despite the increase in driving torque demand; • The increase in the simulated velocity causes changes in the demand for driving torque and work time coverage, the change being approximately a square function; • The set of recorded data does not allow for an unambiguous and direct determination of impact nature of the parameter value a (distance from the axle of the smaller wheel to the center of gravity) on the torque demand and work time coverage. This issue remains open and requires further research; • The most advantageous case is one where the greatest possible mass of the battery is placed as close as possible to the front wheels of the wheelchair.

Informed Consent Statement: Not applicable.
Data Availability Statement: Not applicable.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
A complete list of all notations and physical quantities used in this paper can be found in Appendix A in Table A1.   x(t) -Linear acceleration