Abstract

Motive forces by muscles are applied to different parts of the human body in a periodic fashion when walking at a uniform rate. In this study, the whole human body is modeled as a multidegree of freedom (MDOF) system with seven degrees of freedom. In view of the changing contact conditions with the ground due to alternating feet movements, the system under study is considered piecewise time invariant for each half-period when one foot is in contact with the ground. Forces transmitted from the body to the ground while walking at a normal pace are experimentally measured and numerically simulated. Fourth-order Runge-Kutta method is employed to numerically simulate the forces acting on different masses of the body. An optimization problem is formulated with the squared difference between the measured and simulated forces transmitted to the ground as the objective function, and the motive forces on the body masses as the design variables to solve.

1. Introduction

A human body may essentially be considered a multidegree of freedom (MDOF) vibrating system. Whole body vibration studies and the related adverse health effects have been a topic of great interest and many researchers have attempted to model the parameters of MDOF models of the whole body and body parts [1–3]. Such models have been successfully employed in ride vibrations, ride comfort studies, and driver-seat dynamics [4–7]. Some researchers have studied comparative evaluation of different models of the same body part as well [8, 9]. In these models, the mass information is rather simple; and the stiffness and damping values of the various elements in the MDOF system have been reported based on experimental studies [1, 10]. However, there has virtually been no attempt to identify the muscular motive forces acting on different parts of the human body when the body itself performs simple physical activities instead of being excited by external sources.

Parameter identification can be carried out by subjecting the human body to sinusoidal excitation on a vibrating platform and curve fitting the analytical model results onto the experimental results [11, 12]. However, such vibration excitation can cause undue reactions and discomfort to the body. Even though such a technique can be safely used on a healthy and normal individual, it cannot be used on subjects with poor health or with a medical condition for obvious reasons. For example, R. B. Bhat and P. P. Bhat [13] studied the dynamics of a foetus, uterus, and body model of a pregnant woman in order to study the dynamic force transmitted to the cervix from the foetus when the mother’s body is subjected to ride vibrations. The parameters had to be chosen approximately from available information, as there is no way to obtain these parameters through vibration tests on the pregnant mother.

In this study, the body of a walking human is modeled as a 7 degrees of freedom system, taking into consideration the changing contact conditions due to the single foot contact with the ground during each half-period. Previous studies by Bhat [14] assumed that the sinusoidal excitation was exerted by the ground on the person while walking. McGibbon [15], Lewis and Ferris [16], and Garcia et al. [17] have chosen rather more simplistic models, with only two rigid links connected at the hip and some minor modifications of such model in order to simulate different methods of actuation. A major limitation of these simple bipedal models is that they cannot incorporate hip torques and forces. Furthermore, the contributions of other parts of the human body are completely disregarded in these models. In this study, the 7 DOF model employed does not incorporate any rotation around body joints and motion out of the saggittal plane. However, it does include different parts of the body modeled as simple spring-mass-damper systems. The excitation is considered as the sum total of self-excitation due to the change in the foot contact condition with the ground, the excitation by the muscles, and the joint reaction forces represented by the lumped force acting on each mass.

The force transmitted from the walking human to the ground is simulated in terms of motive forces applied by the muscles on different parts of the human body while walking, which is then curve fitted to the experimentally measured results using MATLAB [18] optimization module. Different motive forces are applied on different parts of the body since the masses, spring constants, and damping coefficients are all different owing to dissimilar morphological characteristics. The motive force defined in the paper does not only result from a specific muscle or muscles in the body, but is also a resultant force coming from active forces from muscles, and passive forces from ligaments, tendons and interactions between bones at joints called joint reaction forces. There are other models that relate motive forces to specific muscles in the body, but the model used in this paper does not consider the details of the anatomical structure of the muscles and other soft tissues. This paper presents the optimized approximations for these musculoskeletal motive forces using the method of least squares by matching the experimentally measured forces transmitted to the ground with simulated results.

2. MDOF Model

The 7 degrees of freedom model of a human body employed to study the walking dynamics is shown in Figure 1(a). Masses π‘š1, π‘š2 and π‘š3 represent the discretized head, trunk which includes the upper limbs, and the hip (or pelvis) parts, while π‘š4 and π‘š5 represent the thighs and π‘š6 and π‘š7 represent the legs including feet on right and left sides of the body, respectively. The π‘˜ξ…žπ‘–π‘  and π‘ξ…žπ‘–π‘  indicate the effective springs and dampers representing the viscoelastic characteristics of the muscles and joints. The springs and dampers π‘˜7, π‘˜8, 𝑐7, and 𝑐8 are varying as constant functions of time, depending on whether both right and left legs are in contact with the ground (double leg support) or one leg is in swing phase (single leg support). These stiffness and damping coefficients in different time periods of a gait cycle are considered as follows:π‘˜7β‰ 0,𝑐7β‰ 0,π‘˜8β‰ 0,𝑐8π‘˜β‰ 0,indoublelegsupport,7β‰ 0,𝑐7β‰ 0,π‘˜8=0,𝑐8π‘˜=0,inrightsinglesupport,7=0,𝑐7=0,π‘˜8β‰ 0,𝑐8β‰ 0,inleftsinglesupport,(1) Figure 1(b) shows the free body diagram of a representative mass π‘š7. The equation of motion for the MDOF whole body system is written as[𝑀][𝐢][𝐾]𝐹{̈π‘₯}+(𝑑){Μ‡π‘₯}+(𝑑){π‘₯}=𝑔+ξ€½πΉπ‘šξ€Ύ(𝑑),(2) where {𝐹𝑔} and {πΉπ‘š(𝑑)} are the gravity and muscular forces, respectively. Gravity forces are easily computed as the weight of each of the masses. Most of the stiffness and damping coefficients are taken from published literature [14], and reasonable values are assumed for the rest. These are presented in Table 1. The undamped natural frequencies (i.e., square roots of eigen values) and corresponding normal modes (eigenvectors) for this 7 DOF model are presented in Table 2 from a MATLAB run.

3. Experimental Measurements

Ground reaction forces (GRF’s) are measured when a human subject walks on force platforms fitted with transducers. During the experiment, the subject walks at a steady pace towards the force platforms for a considerable distance in order to achieve steady state walking conditions. These measured GRF’s can thus be reasonably assumed to be the forces transmitted from the human body to the ground under steady state walking conditions. As the time period in double leg support during walking is very short, it is assumed that only one foot is in contact with the ground at any time and hence the full body weight will be carried by the right leg or the left leg for a half period in a gait cycle. In this experiment, the subject stepped on a set of four force platforms in steady state walking condition, which were arranged in a square pattern as shown in Figure 2. The forces transmitted by the right and the left feet contacting in succession were measured on two force platforms. At the instant of peak transmissibility, the whole body would be standing on only one foot, with the other being entirely in air; however, the measured force transmitted to the force platform would be slightly higher than the observed body weight of the subject owing to the dynamic effects of the body.

Figure 3 shows the results of measurements in one typical run of the experiment. The first peak in the measured signal indicates the force exerted on the platform when the heel hits the platform. The second peak, which is slightly lower than the first, indicates the corresponding force exerted when the toe hits the platform and the heel begins to rise. The trough in between the two peaks is indicative of the straddling motion of the sole between the two contacts. It should also be noted that the signals for the right and the left legs overlap slightly over time due to the fact that the heel of the second foot would have already hit the platform in the front while the toes of the first foot leave contact with the platform behind. The forces exerted by the two feet on the force platforms were measured for a period of approximately 5 s, which is the time required to complete one cycle of walk on the two platforms in line.

4. Formulation and Optimization

The forces transmitted to the ground are expressed as𝑇1=π‘˜7(𝑑)π‘₯6(𝑑)+𝑐7(𝑑)Μ‡π‘₯6𝑇(𝑑),2=π‘˜8(𝑑)π‘₯7(𝑑)+𝑐8(𝑑)Μ‡π‘₯7(𝑑),(3) where 𝑇1 and 𝑇2 represent the forces transmitted through the right and the left feet, respectively.

Casting the equations of motion in the state-space form, fourth-order Runge-Kutta integration is carried out in MATLAB in order to obtain the responses and forces transmitted to the ground. The set of forces {πΉπ‘š(𝑑)} are identified as design variables, and the squared error between the measured and computed forces transmitted to the ground is the objective function to be minimized. During walking the human upper body is usually maintained in a neutral position in which the muscle activities in the upper body for maintaining the posture are trivial compared to the muscle activity in the lower extremities [19]. Therefore, it is assumed in this study that the muscular forces acting on the head mass (π‘š1) and the trunk mass (π‘š2) are zero during walking. The forces at all the other degrees of freedom are periodic. It is also assumed that the forces on the degrees of freedom 4 and 6 are identical to those on 5 and 7, respectively, with a phase shift of 180∘.This is because an ideal human body is completely symmetric along the vertical midline, and one half-side can be taken as an exact mirror image of the other.

From Figure 3 it can be seen that one complete period for the right and left foot contact in succession is 1.2167 s and the transmissibility consists of four peaks corresponding to the successive right and left heel and the toe contact to the ground. Assuming that the four peaks are caused by four force pulses, the forces {πΉπ‘š(𝑑)} are taken asξ€½πΉπ‘šξ€Ύ=𝐴(𝑑)1𝑓1𝐴(𝑑)+2𝑓2𝐴(𝑑)+3𝑓3𝐴(𝑑)+4𝑓4(𝑑),(4) where 𝑓𝑖(𝑑), 𝑖 = 1, 2, 3, and 4 are predefined pulse forms and their amplitudes are given by the vectors {𝐴𝑖}. In general 𝑓𝑖(𝑑) can be of different durations within the period, but in this study they are assumed to be of equal durations and are rectangular, as shown in Figure 4. They are further expanded in Fourier series in order to assign smooth, varying sinusoidal waveforms for the force pulses 𝑓𝑖(𝑑) as below: 𝑓1=π‘Ž1,0+𝑁𝑛=1ξ‚€π‘Ž1,𝑛cos2π‘›πœ‹π‘‘π‘‡ξ‚+𝑏1,𝑛sin2π‘›πœ‹π‘‘π‘‡,𝑓2=π‘Ž2,0+𝑁𝑛=1ξ‚€π‘Ž2,𝑛cos2π‘›πœ‹π‘‘π‘‡ξ‚+𝑏2,𝑛sin2π‘›πœ‹π‘‘π‘‡,𝑓3=π‘Ž3,0+𝑁𝑛=1ξ‚€π‘Ž3,𝑛cos2π‘›πœ‹π‘‘π‘‡ξ‚+𝑏3,𝑛sin2π‘›πœ‹π‘‘π‘‡,𝑓4=π‘Ž4,0+𝑁𝑛=1ξ‚€π‘Ž4,𝑛cos2π‘›πœ‹π‘‘π‘‡ξ‚+𝑏4,𝑛sin2π‘›πœ‹π‘‘π‘‡.(5)

Using (4), the muscular motive forces can be obtained as follows:ξ€½πΉπ‘šξ€Ύ=𝛼(𝑑)0ξ€Ύ+𝑁𝑛=1𝛽𝑛cos2π‘›πœ‹π‘‘π‘‡ξ‚+𝛾𝑛sin2π‘›πœ‹π‘‘π‘‡ξƒ°,𝛼where0ξ€Ύ=𝐴1ξ€Ύπ‘Ž1,0+𝐴2ξ€Ύπ‘Ž2,0+𝐴3ξ€Ύπ‘Ž3,0+𝐴4ξ€Ύπ‘Ž4,0,𝛽𝑛=𝐴1ξ€Ύπ‘Ž1,𝑛+𝐴2ξ€Ύπ‘Ž2,𝑛+𝐴3ξ€Ύπ‘Ž3,𝑛+𝐴4ξ€Ύπ‘Ž4,𝑛,𝛾𝑛=𝐴1𝑏1,𝑛+𝐴2𝑏2,𝑛+𝐴3𝑏3,𝑛+𝐴4𝑏4,𝑛.(6)

The resulting four shapes, 𝑓𝑖(𝑑), 𝑖 = 1, 2, 3, and 4 are shown in Figure 5.

Nonlinear iterative least square curve fitting from MATLAB [18] optimization toolbox is used, and the error between the measured and computed transmitted forces to the ground is calculated using the following equation:1𝐸=2𝑁𝑖=1𝑇1,𝑖𝐴1𝑇,𝐴2𝑇,𝐴3𝑇,𝐴4ξ€Ύπ‘‡ξ‚„ξ‚βˆ’π‘‡exp,𝑖2,(7) where N is the size of the vector {𝑇1} and {𝑇exp}, and the notations {𝑇1,𝑖} and {𝑇exp,𝑖} stand for the 𝑖th element of the vectors {𝑇1} and {𝑇exp}, respectively. Note that {𝑇1} is a function of {𝐴1}, {𝐴2}, {𝐴3}, and {𝐴4} and is not multiplied by them.

Figure 6 shows a workflow schematic of the numerical integration and optimization process. Only 𝑇1 is considered in numerical computation assuming the symmetrical nature of the problem as indicated before.

5. Results and Discussion

Least squares curve fitting resulted in a set of coefficients that indicated the transmitted forces as shown in Figure 7. The measured forces transmitted to the ground compare closely with those obtained by simulation.

Figure 8 shows the force profile through one period of the force {𝐹3(𝑑)} acting on the hip. Figure 9 shows the plots of {𝐹4(𝑑)} and {𝐹5(𝑑)}, and Figure 10 shows those of {𝐹6(𝑑)} and {𝐹7(𝑑)}. It can be seen from both figures that these forces are similar with a shift of one half-period, which is originally set as a constraint in the simulation. Moreover, Figure 10 shows that during the foot contact, precisely represented from time 0.6 s to 1.2 s of {𝐹6(𝑑)} corresponding to the step involving the spring and damper π‘˜7 and 𝑐7, the lumped muscular force acting on the mass π‘š6 is positive when the heel is landing and the load is all directed to the ground. In the next half-period, the muscle force applied to the foot is lifting the whole body to move upward. This force is almost equal to the weight of the body.

Comparing the force profiles in Figures 8, 9, and 10, {𝐹3(𝑑)} shows the highest overall value. This result is in agreement with previous observations reported in other published literature [16, 17, 20, 21], that the force acting on the hips is quite significant in comparison to that acting on the thighs and legs during walking or running activity. This significant force is assumedly due to the greatly increased muscle contraction of the hip abductor muscles for maintaining the upper body in a vertical position during walking with a single leg support.

The implications of this conclusion are quite significant to those suffering from musculoskeletal disorders, osteoarthritis, or diabetes in advanced stages, especially in the hips and the lower limbs. People suffering from joint pains in the lower limbs may use a less tiring mechanism with smaller step length during walking. Conversely, people with hip pain can walk with an exaggerated lower limb movement. Such recommendations have also been made by Mueller et al. [21] and Lewis and Sahrmann [22]. Based on the reason proposed above, a walking pattern in which a small lateral bending of upper body to the ipsilateral side of the supporting leg during walking might be able to reduce the hip abductor muscle contraction so that the load applied to the hip joint could be significantly reduced. This assumption can be investigated in further studies with this 7 DOF model. The effect of arm swing on ground reaction forces in gait also contributes substantially to the vertical ground reaction forces, and an 8 or 9 DOF model with the arms separated from the trunk could be a better representation for including the arm swing effect [23, 24].

The 7 DOF model incorporates contributions of different parts of the body on force distributions in the 5 relevant degrees of freedom that are involved in the physical process of walking. While such a detailed model has rarely been studied hitherto in gait studies, a major limitation of this 7 DOF model is that it considers the motion of hips, thighs, and legs only in the vertical direction in the saggittal plane. This simplified 7 DOF model also tends to ignore the rotational movements of different body parts around their joints while walking. The upper body parts and lower limbs may also have small movements laterally in the frontal plane, which will affect the muscular forces in the vertical direction. Furthermore, the relative movement between the legs and feet is not incorporated as separate masses with additional DOF’s in this model.

In the normal walking process, the arms usually swing back and forth which would contribute to the inertia forces of the body. However, in this model the upper limbs are incorporated into the trunk mass (π‘š2) and their movements with respect to the trunk during walking are completely ignored. The effect of this simplification on inertia force of arms depends on the degree of arm swinging. In this study, the simplification was justified by the fact that the degree of arm swinging was found small and the effect of inertia force induced by the arm movement in vertical direction could be largely mitigated.

It is noted that all the degrees of freedom in this model are taken as lumped masses, as it is difficult to come up with an exact formulation of distributed mass system of different body parts. Such data is largely non-existent in biomedical literature, and moreover, a universal distributed mass system is extremely difficult to formulate owing to different musculoskeletal characteristics of individual people. This model is at best simply a rough approximation of relevant parts of the human body with a focus on predicting and comparing force distributions in chosen discretized masses while walking at a normal pace and with a normal gait.

6. Conclusions

Experiments were carried out to measure forces transmitted to force platforms when a human subject is walking. The human body is modeled as a 7 degrees of freedom system. The forces transmitted to the ground are measured and the muscle forces acting on the masses are estimated by fitting least squares curve to the measured forces transmitted to the ground. The results show that the force acting on the hip is quite significant in comparison to forces acting on the thighs and legs. Relevant conclusions are drawn about the suitability of this 7 DOF model and the implications of these observations for people with musculoskeletal disabilities or ailments.