# Mathematical Modelling of Biomechanical Interactions between Backpack and Bearer during Load Carriage

Authors: Lei Ren; David Howard; Richard K. Jones
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

This paper proposes a three-dimensional mathematical model of the biomechanical interactions between backpack and bearer during load carriage. The model considers both the coupled pack motions, which follow the torso, and also the longitudinal compliance and damping in the backpack suspension. The pack interaction forces and moments, acting on the bearer, are determined from kinematic relationships, equations of motion, and a dynamic pack suspension model. The parameters of the pack suspension model were identified from test data obtained using a load carriage test rig. Output from the load carriage mathematical model has been compared with measurement data during human gait and conclusions drawn with regard to the validity of the proposed approach.


---

## Body

## 1. Introduction

Backpacks are common devices for increasing human load carriage performance but when heavily loaded may lead to excessive joint loadings, muscle fatigue, or even injury [1–5]. Most studies of personal load carriage systems have focused on physiological and biomechanical aspects [6–10]. Biomechanical studies have concentrated on experimental gait analysis, including the effects of load carriage on electromyographic activity [3, 8, 11], gait and posture [6, 7, 10], and ground reactions [6, 8, 9]. Almost all of these studies consider the backpack to be part of the trunk segment, with no relative motion between the two. Little is known about the pack interaction forces and moments acting on the bearer’s trunk.During load carriage, the pack interaction forces exerted on the torso relate directly to perceived discomfort, fatigue, and the risk of injury, for example, rucksack palsy and back problems [12, 13]. A better understanding of these interaction forces would help to improve the design of future load carriage systems. Unfortunately, in contrast to the contact pressure distribution [14–16], the interaction forces and moments between pack and torso cannot be measured directly; however mathematical modelling and simulation offer an alternative approach.Very little has been published on the modelling of load carriage biomechanics. While examining the effects of pack load on back muscle EMG, Bobet and Norman [3] considered the pack to be part of the trunk segment for inverse dynamics purposes. Recently, Pelot et al. [17] developed a two-dimensional backpack model to investigate the forces at the shoulder and hip belts. However, their model assumes no relative movement between pack and torso and can only be used to evaluate static forces. To properly model backpack dynamics, it is necessary to include the relative motions between pack and torso, which is capable of representing the dynamic impacts on the bearer rather than static interactions only.This paper proposes a mathematical load carriage model, which represents the three-dimensional biomechanical interactions between pack and torso. All pack motions are directly coupled to those of the bearer’s trunk, with the exception of longitudinal motion up and down the bearer’s back, which depends on a pack suspension model, the parameters of which were identified from test data obtained using a load carriage test rig [18]. The validity of the proposed approach is investigated and conclusions drawn with regard to future developments in the modelling of load carriage biomechanics.


## 2. Methods

### 2.1. Pack-Trunk Kinematics

Throughout this study, the typical backpack system with two shoulder straps and a waist belt is addressed. If deformations of the backpack are neglected, the pack can be considered to be a rigid body (Figure1). Due to the actions of the shoulder straps and hip belt, holding the pack against the bearer’s back, pack rotations relative to the trunk in the sagittal and transverse planes, and translation perpendicular to the back are all constrained. For simplicity, relative pack rotation in the frontal plane and lateral pack translation are also considered to be negligible. In other words, pack motion relative to the trunk only occurs in a direction which coincides with the trunk’s longitudinal axis. So, the backpack is modelled as a rigid body that slides up and down the bearer’s back while all other degrees of freedom follow the motions of the trunk.Figure 1
The 3D pack and trunk motions during load carriage and the three coordinate systems.OXYZ is the global (inertial) coordinate system, otxtytzt is the trunk coordinate system, and opxpypzp is the backpack coordinate system.A global coordinate systemOXYZ, a trunk system otxtytzt, and a backpack system opxpypzp are defined (Figure 1). The global system’s x-axis lies in the sagittal plane and points in the direction of forward progress, the y-axis also lies in the sagittal plane and points upwards, and the z-axis lies in the frontal plane and points to the right. The trunk coordinate system otxtytzt is attached to the bearer’s trunk with its origin located at the trunk centre of mass (CoM). The yt-axis is parallel to the trunk longitudinal axis (midpoint between C7 and T8 spinous process to the midpoint between xiphoid process and jugular notch) and points upwards. The xt-axis is perpendicular to the bearer’s back and points in the anterior direction. The zt-axis is perpendicular to both xt-axis and yt-axis and points to the right. The local pack coordinate system opxpypzp is attached to the backpack with its origin at the pack CoM, and its axes are parallel to those of the trunk coordinate system. In the analysis below, we have followed the common practice of not expressing vector equations in any particular coordinate frame. Where vector equations have been expanded, to derive the individual components, the resulting scalar equations are expressed in a particular coordinate frame, which is stated.Based on the rigid body pack model proposed above, the pack’s kinematic relationship with the trunk can be described by(1a)a⃑p=a⃑t+a⃑ptn+a⃑ptτ+a⃑r+a⃑rc,(1b)ω⃑p=ω⃑t,(1c)α⃑p=α⃑t,wherea⃑p is the absolute acceleration of the pack CoM; a⃑t is the absolute acceleration of the trunk CoM; a⃑ptn  and  a⃑ptτ are the centripetal and tangential acceleration vectors, relative to the trunk CoM, of a point fixed in the trunk coordinate frame and instantaneously coincident with the pack CoM; a⃑r is the relative pack acceleration vector along the back; and a⃑rc is the Coriolis acceleration vector.The various acceleration terms in (1a) are given by(2a)a⃑ptτ=α⃑p×r⃑pt,(2b)a⃑ptn=ω⃑t×(ω⃑t×r⃑pt),(2c)a⃑r=u⃑¨,(2d)a⃑rc=2ω⃑p×u⃑˙,wherer⃑pt=(d⃑+u⃑) is the position vector of the pack CoM relative to the trunk CoM, u⃑ is the pack translation relative to the trunk, and d⃑ is the position of the pack CoM, relative to the trunk CoM, in the unloaded condition.Substituting  (2a), (2b), (2c), and (2d) into (1a) and considering each component in the backpack coordinate system, the absolute translational accelerations of the pack CoM in all three directions can be derived as(3a)apx=atx+αtydz-αtz(dy+u)+ωtxωty(dy+u)-(ωty2+ωtz2)dx+ωtxωtzdz-2ωtzu˙,(3b)apy=aty+αtzdx-αtxdz+ωtyωtzdz-(ωtx2+ωtz2)(dy+u)+ωtxωtydx+u¨,(3c)apz=atz+αtx(dy+u)-αtydx+ωtxωtzdx-(ωtx2+ωty2)dz+ωtyωtz(dy+u)+2ωtxu˙.From (3a), (3b), and (3c), it can be seen that the backpack’s absolute accelerations are not only influenced by trunk translation but also strongly coupled with trunk rotation and relative pack motion.


### 2.2. Backpack’s Equation of Motion

If the backpack is regarded as a separate object and isolated from the torso (Figure2), the pack interaction forces and moments can be assessed. In this case, the interaction between the pack and the body is described by the resultant force and moment (F⃑p and M⃑p) acting on the pack and at its CoM. This equivalent force system represents the net effects of the various forces at the pack-human interface as they act on the pack.Figure 2
The human-pack interaction forceF⃑p (Fpx, Fpy, Fpz) and moment M⃑p acting at the pack mass centre, and the associated reactions F⃑pack and M⃑pack acting at the trunk mass centre.Applying the Newton-Euler equations leads to(4a)mpa⃑p=F⃑p+mpg⃑,(4b)Jpα⃑p+ω⃑p×(Jpω⃑p)=M⃑p.Equation (4a) can be expressed in the backpack coordinate system as(5a)Fpx=mpapx+mpg(sinψcosϕ+cosψcosθsinϕ),(5b)Fpy=mpapy+mpg(sinψsinϕ+cosψcosθcosϕ),(5c)Fpz=mpapz-mpgcosψsinθ,where the orientation of the backpack with respect to the global coordinate systemOXYZ is defined by Euler angles (ψ,θ,ϕ) in ZXZ sequence.


### 2.3. Backpack Suspension Model

In this study, a general non-linear pack suspension model is employed, which relates the pack interaction force along the trunk longitudinal axis to the relative pack motion. The relationship can be written as(6)Fpy=Fpe+Fpd+Fpi,
where Fpe is the elastic component of the interaction force, Fpd is the damping component, and Fpi is an inertial coupling term. Cubic polynomials are used to describe these non-linear properties [19] as follows:
(7)Fpe=a3u3+sign(u)a2u2+a1u,Fpd=b3u˙3+sign(u˙)b2u˙2+b1u˙,Fpi=c3u¨3+sign(u¨)c2u¨2+c1u¨,
where a1, a2, a3, b1, b2, b3, c1, c2, c3 are model parameters, which are constant for a given type of pack and specified working conditions (i.e., pack load, load distribution, strap and belt tensions, etc.).For a particular backpack, the parameters of the suspension model can be identified from dynamic test data, obtained using the hydraulically driven load carriage test rig shown in Figure3 [19]. The mannequin was covered in a neoprene-like material to mimic the soft tissues. The hydraulic ram drives the mannequin up and down with different frequency and amplitude inputs, which allows non-linear frequency response testing over a range of frequencies and amplitudes. The harmonic analysis method was used to identify the model parameters [18].Figure 3
The dynamic load carriage test rig. The hydraulic ram drives the mannequin up and down at different frequencies and amplitudes. Accelerometers measure the motion of the backpack and the mannequin, and a load cell measures the dynamic force propelling the mannequin.


### 2.4. Pack-Bearer Interaction Forces during Walking

Using (1a), (1b), and (1c) to (7) and measured trunk motion data, the pack interaction forces and moments can be derived as follows. Substituting the expression for acceleration apy (3b) into the pack’s equation of motion (5b), and using the backpack suspension model (6) and (7) to substitute for Fpy, leads to the following non-linear second-order differential equation, which describes the dynamic pack behaviour along the bearer’s back:
(8)c3u¨3+c2sign(u¨)u¨2+(c1+mp)u¨+b3u˙3+b2sign(u˙)u˙2+b1u˙+a3u3+a2sign(u)u2+(a1-mp(ωx2+ωz2))u+mp(aty+αzdx-αxdz+ωyωzdzaty+αzdx-αxdz+ωyωzdz-(ωx2+ωz2)dy+ωxωydx+g(sinψsinϕ+cosψcosθcosϕ)aty+αzdx-αxdz+ωyωzdz)=0.As this differential equation is of a highly non-linear form and some of the terms are time varying, a numerical solution is required (4th-order Runge-Kutta). Because the numerical integration time step is normally smaller than the gait measurement interval, cubic interpolation is used to provide trunk motion data at the necessary frequency. The initial values of the relative pack displacement and velocity were set to zero, and the numerical integration algorithm executed until a steady-state cyclic pack motion was achieved. Steady state was defined as the initial and final values of the state variables for one gait cycle being equal within an acceptable tolerance (1.0e-6 was used in this study).Thus, given the trunk motion data (aty, ψ, θ, ϕ, ωx, ωy, ωz, αx, αz), the relative pack motion can be derived (u, u˙, and u¨). Then the pack force along the bearer’s back Fpy can be obtained by substituting the calculated relative motion (u, u˙, and u¨) into (6) and (7).The expressions for the normal pack forceFpx and lateral pack force Fpz are obtained by substituting the expressions for accelerations apx and apz ((3a) and (3c)) into the pack’s equations of motion ((5a) and (5c)), which leads to
(9)Fpx=mp(atx+αydz-αz(dy+u)+ωxωy(dy+u)-(ωy2+ωz2)dx+ωxωzdz-2ωzu˙+g(sinψcosϕ+cosψcosθsinϕ)αz(dy+u)),Fpz=mp(atz+αx(dy+u)-αydx+ωxωzdx-(ωx2+ωy2)dz+ωyωz(dy+u)+2ωxu˙-gcosψsinθ(dy+u)).Thus, using (6), (7), and (9), the three components of the pack force, acting on the pack CoM, can be determined from the trunk motion data and the calculated relative pack motion. The resultant pack moment, acting about the pack CoM, can be calculated directly from the trunk rotations, given that the pack rotates with the trunk (see (1b), (1c), and (4b)).The pack interaction force and moment can also be represented by an equivalent force system (F⃑pack, M⃑pack) acting on the trunk CoM (Figure 2). According to Newton’s third law, the interaction force and moment acting at the trunk CoM can be expressed as
(10)M⃑pack=F⃑p×r⃑tp-M⃑p,F⃑pack=-F⃑p.The pack force and moment can be used to assess the mechanical loads imposed on the human body during load carriage. Moreover, when input into a three-dimensional human gait model, joint loads and mechanical energy expenditure can also be evaluated.The proposed dynamic load carriage model has been implemented in the MATLAB programming environment (see Figure4 for a schematic diagram). Measured 3D trunk rotations and translations over a complete gait cycle, with a loaded backpack, were used as input data.Figure 4
The schematic diagram for the pack-bearer interaction force calculation based on pack-trunk kinematics, backpack’s equation of motion, and backpack suspension model.


### 2.5. Experimental Model Assessment

To assess the validity of the proposed modelling approach, three-dimensional gait measurement was conducted to capture trunk and pack motions whilst carrying a backpack. Two healthy male subjects (age:30±2 yrs, weight: 75±3.2 kg) were selected from a population of postgraduate students. Prior to participation, the subjects provided informed consent in accordance with the policies of Salford University’s Ethical Advisory Committee. The subjects walked in trainers along a walkway, inside a gait laboratory, while motion data was collected at 100 Hz using a 6-camera Qualisys motion analysis system (Qualisys AB, Savedalen, Sweden). The subjects walked at two different speeds, normal and fast, with two backpack loads, 11.5 kg and 23.0 kg (see Figure 5). Each experimental condition was measured six times to ensure that a repeatable data set for a complete walking cycle was obtained.Figure 5
The 3D gait measurement setup for capturing the trunk and pack motions during load carriage. A specially designed plastic plate carrying four markers, firmly attached to the thorax, was used to capture the position and orientation of the trunk. Eight markers were attached to the corners of the pack to capture its position and orientation.To capture trunk motions, a specially designed plastic plate, carrying four reflective markers, was firmly attached to the thorax [20]. Four anatomical landmarks (C7 and T8 spinous processes, xiphoid process, and jugular notch) were used to define the trunk coordinate system. Before the walking trials, a static calibration procedure was used to locate these anatomical landmarks based on the CAST methodology (calibrated anatomical system technique) [21]. The position of the trunk CoM was estimated from the anatomical landmarks using anthropometric data from the literature [22].To minimize the deformation of the backpack, the test pack contained a specially designed plastic foam block with cylindrical cavities to allow metal bars to be firmly inserted. By putting different metal bars in different cavities, different pack loads and load distributions were possible. During walking trials, eight reflective markers were attached to the corners of the backpack, and the geometric centre of the pack was assumed to be equidistant from these markers. The pack’s CoM position and inertia matrix were then calculated based on this geometric centre and the load distribution within the pack.The raw 3D marker data were processed using a custom-written software package, SMAS (Salford Motion Analysis Software), which has been developed for 3D kinematic and kinetic analysis of general biomechanical multibody systems [20, 23]. The data were filtered using a low pass 4th-order Butterworth digital filter with a cut-off frequency of 15 Hz. The orientations of the backpack and trunk segments were derived using an optimal estimation algorithm [24]. The pack and trunk angular velocities, angular accelerations, and CoM linear accelerations were calculated using a finite difference method.Based on these measurement data, the interaction force and moment acting at the pack CoM were estimated as follows:(11)F⃑pm=mp(a⃑pm-g⃑),M⃑pm=Jpα⃑pm+ω⃑pm×(Jpω⃑pm).These data for the whole gait cycle were then used to validate the 3D load carriage model.


## 2.1. Pack-Trunk Kinematics

Throughout this study, the typical backpack system with two shoulder straps and a waist belt is addressed. If deformations of the backpack are neglected, the pack can be considered to be a rigid body (Figure1). Due to the actions of the shoulder straps and hip belt, holding the pack against the bearer’s back, pack rotations relative to the trunk in the sagittal and transverse planes, and translation perpendicular to the back are all constrained. For simplicity, relative pack rotation in the frontal plane and lateral pack translation are also considered to be negligible. In other words, pack motion relative to the trunk only occurs in a direction which coincides with the trunk’s longitudinal axis. So, the backpack is modelled as a rigid body that slides up and down the bearer’s back while all other degrees of freedom follow the motions of the trunk.Figure 1
The 3D pack and trunk motions during load carriage and the three coordinate systems.OXYZ is the global (inertial) coordinate system, otxtytzt is the trunk coordinate system, and opxpypzp is the backpack coordinate system.A global coordinate systemOXYZ, a trunk system otxtytzt, and a backpack system opxpypzp are defined (Figure 1). The global system’s x-axis lies in the sagittal plane and points in the direction of forward progress, the y-axis also lies in the sagittal plane and points upwards, and the z-axis lies in the frontal plane and points to the right. The trunk coordinate system otxtytzt is attached to the bearer’s trunk with its origin located at the trunk centre of mass (CoM). The yt-axis is parallel to the trunk longitudinal axis (midpoint between C7 and T8 spinous process to the midpoint between xiphoid process and jugular notch) and points upwards. The xt-axis is perpendicular to the bearer’s back and points in the anterior direction. The zt-axis is perpendicular to both xt-axis and yt-axis and points to the right. The local pack coordinate system opxpypzp is attached to the backpack with its origin at the pack CoM, and its axes are parallel to those of the trunk coordinate system. In the analysis below, we have followed the common practice of not expressing vector equations in any particular coordinate frame. Where vector equations have been expanded, to derive the individual components, the resulting scalar equations are expressed in a particular coordinate frame, which is stated.Based on the rigid body pack model proposed above, the pack’s kinematic relationship with the trunk can be described by(1a)a⃑p=a⃑t+a⃑ptn+a⃑ptτ+a⃑r+a⃑rc,(1b)ω⃑p=ω⃑t,(1c)α⃑p=α⃑t,wherea⃑p is the absolute acceleration of the pack CoM; a⃑t is the absolute acceleration of the trunk CoM; a⃑ptn  and  a⃑ptτ are the centripetal and tangential acceleration vectors, relative to the trunk CoM, of a point fixed in the trunk coordinate frame and instantaneously coincident with the pack CoM; a⃑r is the relative pack acceleration vector along the back; and a⃑rc is the Coriolis acceleration vector.The various acceleration terms in (1a) are given by(2a)a⃑ptτ=α⃑p×r⃑pt,(2b)a⃑ptn=ω⃑t×(ω⃑t×r⃑pt),(2c)a⃑r=u⃑¨,(2d)a⃑rc=2ω⃑p×u⃑˙,wherer⃑pt=(d⃑+u⃑) is the position vector of the pack CoM relative to the trunk CoM, u⃑ is the pack translation relative to the trunk, and d⃑ is the position of the pack CoM, relative to the trunk CoM, in the unloaded condition.Substituting  (2a), (2b), (2c), and (2d) into (1a) and considering each component in the backpack coordinate system, the absolute translational accelerations of the pack CoM in all three directions can be derived as(3a)apx=atx+αtydz-αtz(dy+u)+ωtxωty(dy+u)-(ωty2+ωtz2)dx+ωtxωtzdz-2ωtzu˙,(3b)apy=aty+αtzdx-αtxdz+ωtyωtzdz-(ωtx2+ωtz2)(dy+u)+ωtxωtydx+u¨,(3c)apz=atz+αtx(dy+u)-αtydx+ωtxωtzdx-(ωtx2+ωty2)dz+ωtyωtz(dy+u)+2ωtxu˙.From (3a), (3b), and (3c), it can be seen that the backpack’s absolute accelerations are not only influenced by trunk translation but also strongly coupled with trunk rotation and relative pack motion.


## 2.2. Backpack’s Equation of Motion

If the backpack is regarded as a separate object and isolated from the torso (Figure2), the pack interaction forces and moments can be assessed. In this case, the interaction between the pack and the body is described by the resultant force and moment (F⃑p and M⃑p) acting on the pack and at its CoM. This equivalent force system represents the net effects of the various forces at the pack-human interface as they act on the pack.Figure 2
The human-pack interaction forceF⃑p (Fpx, Fpy, Fpz) and moment M⃑p acting at the pack mass centre, and the associated reactions F⃑pack and M⃑pack acting at the trunk mass centre.Applying the Newton-Euler equations leads to(4a)mpa⃑p=F⃑p+mpg⃑,(4b)Jpα⃑p+ω⃑p×(Jpω⃑p)=M⃑p.Equation (4a) can be expressed in the backpack coordinate system as(5a)Fpx=mpapx+mpg(sinψcosϕ+cosψcosθsinϕ),(5b)Fpy=mpapy+mpg(sinψsinϕ+cosψcosθcosϕ),(5c)Fpz=mpapz-mpgcosψsinθ,where the orientation of the backpack with respect to the global coordinate systemOXYZ is defined by Euler angles (ψ,θ,ϕ) in ZXZ sequence.


## 2.3. Backpack Suspension Model

In this study, a general non-linear pack suspension model is employed, which relates the pack interaction force along the trunk longitudinal axis to the relative pack motion. The relationship can be written as(6)Fpy=Fpe+Fpd+Fpi,
where Fpe is the elastic component of the interaction force, Fpd is the damping component, and Fpi is an inertial coupling term. Cubic polynomials are used to describe these non-linear properties [19] as follows:
(7)Fpe=a3u3+sign(u)a2u2+a1u,Fpd=b3u˙3+sign(u˙)b2u˙2+b1u˙,Fpi=c3u¨3+sign(u¨)c2u¨2+c1u¨,
where a1, a2, a3, b1, b2, b3, c1, c2, c3 are model parameters, which are constant for a given type of pack and specified working conditions (i.e., pack load, load distribution, strap and belt tensions, etc.).For a particular backpack, the parameters of the suspension model can be identified from dynamic test data, obtained using the hydraulically driven load carriage test rig shown in Figure3 [19]. The mannequin was covered in a neoprene-like material to mimic the soft tissues. The hydraulic ram drives the mannequin up and down with different frequency and amplitude inputs, which allows non-linear frequency response testing over a range of frequencies and amplitudes. The harmonic analysis method was used to identify the model parameters [18].Figure 3
The dynamic load carriage test rig. The hydraulic ram drives the mannequin up and down at different frequencies and amplitudes. Accelerometers measure the motion of the backpack and the mannequin, and a load cell measures the dynamic force propelling the mannequin.


## 2.4. Pack-Bearer Interaction Forces during Walking

Using (1a), (1b), and (1c) to (7) and measured trunk motion data, the pack interaction forces and moments can be derived as follows. Substituting the expression for acceleration apy (3b) into the pack’s equation of motion (5b), and using the backpack suspension model (6) and (7) to substitute for Fpy, leads to the following non-linear second-order differential equation, which describes the dynamic pack behaviour along the bearer’s back:
(8)c3u¨3+c2sign(u¨)u¨2+(c1+mp)u¨+b3u˙3+b2sign(u˙)u˙2+b1u˙+a3u3+a2sign(u)u2+(a1-mp(ωx2+ωz2))u+mp(aty+αzdx-αxdz+ωyωzdzaty+αzdx-αxdz+ωyωzdz-(ωx2+ωz2)dy+ωxωydx+g(sinψsinϕ+cosψcosθcosϕ)aty+αzdx-αxdz+ωyωzdz)=0.As this differential equation is of a highly non-linear form and some of the terms are time varying, a numerical solution is required (4th-order Runge-Kutta). Because the numerical integration time step is normally smaller than the gait measurement interval, cubic interpolation is used to provide trunk motion data at the necessary frequency. The initial values of the relative pack displacement and velocity were set to zero, and the numerical integration algorithm executed until a steady-state cyclic pack motion was achieved. Steady state was defined as the initial and final values of the state variables for one gait cycle being equal within an acceptable tolerance (1.0e-6 was used in this study).Thus, given the trunk motion data (aty, ψ, θ, ϕ, ωx, ωy, ωz, αx, αz), the relative pack motion can be derived (u, u˙, and u¨). Then the pack force along the bearer’s back Fpy can be obtained by substituting the calculated relative motion (u, u˙, and u¨) into (6) and (7).The expressions for the normal pack forceFpx and lateral pack force Fpz are obtained by substituting the expressions for accelerations apx and apz ((3a) and (3c)) into the pack’s equations of motion ((5a) and (5c)), which leads to
(9)Fpx=mp(atx+αydz-αz(dy+u)+ωxωy(dy+u)-(ωy2+ωz2)dx+ωxωzdz-2ωzu˙+g(sinψcosϕ+cosψcosθsinϕ)αz(dy+u)),Fpz=mp(atz+αx(dy+u)-αydx+ωxωzdx-(ωx2+ωy2)dz+ωyωz(dy+u)+2ωxu˙-gcosψsinθ(dy+u)).Thus, using (6), (7), and (9), the three components of the pack force, acting on the pack CoM, can be determined from the trunk motion data and the calculated relative pack motion. The resultant pack moment, acting about the pack CoM, can be calculated directly from the trunk rotations, given that the pack rotates with the trunk (see (1b), (1c), and (4b)).The pack interaction force and moment can also be represented by an equivalent force system (F⃑pack, M⃑pack) acting on the trunk CoM (Figure 2). According to Newton’s third law, the interaction force and moment acting at the trunk CoM can be expressed as
(10)M⃑pack=F⃑p×r⃑tp-M⃑p,F⃑pack=-F⃑p.The pack force and moment can be used to assess the mechanical loads imposed on the human body during load carriage. Moreover, when input into a three-dimensional human gait model, joint loads and mechanical energy expenditure can also be evaluated.The proposed dynamic load carriage model has been implemented in the MATLAB programming environment (see Figure4 for a schematic diagram). Measured 3D trunk rotations and translations over a complete gait cycle, with a loaded backpack, were used as input data.Figure 4
The schematic diagram for the pack-bearer interaction force calculation based on pack-trunk kinematics, backpack’s equation of motion, and backpack suspension model.


## 2.5. Experimental Model Assessment

To assess the validity of the proposed modelling approach, three-dimensional gait measurement was conducted to capture trunk and pack motions whilst carrying a backpack. Two healthy male subjects (age:30±2 yrs, weight: 75±3.2 kg) were selected from a population of postgraduate students. Prior to participation, the subjects provided informed consent in accordance with the policies of Salford University’s Ethical Advisory Committee. The subjects walked in trainers along a walkway, inside a gait laboratory, while motion data was collected at 100 Hz using a 6-camera Qualisys motion analysis system (Qualisys AB, Savedalen, Sweden). The subjects walked at two different speeds, normal and fast, with two backpack loads, 11.5 kg and 23.0 kg (see Figure 5). Each experimental condition was measured six times to ensure that a repeatable data set for a complete walking cycle was obtained.Figure 5
The 3D gait measurement setup for capturing the trunk and pack motions during load carriage. A specially designed plastic plate carrying four markers, firmly attached to the thorax, was used to capture the position and orientation of the trunk. Eight markers were attached to the corners of the pack to capture its position and orientation.To capture trunk motions, a specially designed plastic plate, carrying four reflective markers, was firmly attached to the thorax [20]. Four anatomical landmarks (C7 and T8 spinous processes, xiphoid process, and jugular notch) were used to define the trunk coordinate system. Before the walking trials, a static calibration procedure was used to locate these anatomical landmarks based on the CAST methodology (calibrated anatomical system technique) [21]. The position of the trunk CoM was estimated from the anatomical landmarks using anthropometric data from the literature [22].To minimize the deformation of the backpack, the test pack contained a specially designed plastic foam block with cylindrical cavities to allow metal bars to be firmly inserted. By putting different metal bars in different cavities, different pack loads and load distributions were possible. During walking trials, eight reflective markers were attached to the corners of the backpack, and the geometric centre of the pack was assumed to be equidistant from these markers. The pack’s CoM position and inertia matrix were then calculated based on this geometric centre and the load distribution within the pack.The raw 3D marker data were processed using a custom-written software package, SMAS (Salford Motion Analysis Software), which has been developed for 3D kinematic and kinetic analysis of general biomechanical multibody systems [20, 23]. The data were filtered using a low pass 4th-order Butterworth digital filter with a cut-off frequency of 15 Hz. The orientations of the backpack and trunk segments were derived using an optimal estimation algorithm [24]. The pack and trunk angular velocities, angular accelerations, and CoM linear accelerations were calculated using a finite difference method.Based on these measurement data, the interaction force and moment acting at the pack CoM were estimated as follows:(11)F⃑pm=mp(a⃑pm-g⃑),M⃑pm=Jpα⃑pm+ω⃑pm×(Jpω⃑pm).These data for the whole gait cycle were then used to validate the 3D load carriage model.


## 3. Results

Experimental and simulated data for one subject walking at two different speeds, normal (1.6 m/s) and fast (1.9 m/s), with two backpack loads, 11.5 kg and 23.0 kg, are presented in Figures6 to 11. Results are shown for repeated experimental trials and the corresponding simulations.Figure 6
The simulated relative pack displacement along the back (mean (solid red line)± one standard deviation (dashed red lines), (a)), compared with measured data (mean (solid blue line) ± one standard deviation (dashed blue lines), (a)), for six repeated trials for one subject (age: 28 yrs, weight: 71.8 kg, height: 178 cm), and measured vertical trunk CoM displacement (mean (solid blue line) ± one standard deviation (dashed blue lines), (b)). The average walking speed was 1.63±0.03 m/s, and the pack load was 11.5 kg. The stance phase for the left limb is from 0 to 60%, and the swing phase is from 60% to 100%. The double support phase is from 0 to 10% and from 50 to 60%.Figures6 and 7 show the vertical trunk CoM motion, the simulated relative pack displacement, and the pack interaction forces, over a complete gait cycle, for normal walking with an 11.5 kg load. Figures 8 and 9 show the same data for fast walking and an 11.5 kg load. Figures 10 and 11 show the same data for normal walking with a 23 kg load.Figure 7
The simulated pack interaction forces in the pack coordinate system mean ((solid red line))± one standard deviation (dashed red lines): (a) longitudinal pack force, (b) normal pack force (c) and lateral pack force, compared with measured data (mean (solid blue line) ± one standard deviation (dashed blue lines)), for six repeated trials for one subject (age: 28 yrs, weight: 71.8 kg, height: 178 cm). The average walking speed was 1.63±0.03 m/s, and the pack load was 11.5 kg. The stance phase for the left limb is from 0 to 60%, and the swing phase is from 60% to 100%. The double support phase is from 0 to 10% and from 50 to 60%.Figure 8
The simulated relative pack displacement along the back (mean (solid red line)± one standard deviation (dashed red lines), (a)), compared with measured data (mean (solid blue line) ± one standard deviation (dashed blue lines), (a)), for six repeated trials for one subject (age: 28 yrs, weight: 71.8 kg, height: 178 cm), and measured vertical trunk CoM displacement (mean (solid blue line) ± one standard deviation (dashed blue lines), (b)). The average walking speed was 1.93±0.03 m/s, and the pack load was 11.5 kg. The stance phase for the left limb is from 0 to 59%, and the swing phase is from 59% to 100%. The double support phase is from 0 to 9% and from 50 to 59%.Figure 9
The simulated pack interaction forces in the pack coordinate system (mean (solid red line)± one standard deviation (dashed red lines)): (a) longitudinal pack force, (b) normal pack force (c), and lateral pack force, compared with measured data (mean (solid blue line) ± one standard deviation (dashed blue lines)), for six repeated trials for one subject (age: 28 yrs, weight: 71.8 kg, height: 178 cm). The average walking speed was 1.93±0.03 m/s, and the pack load was 11.5 kg. The stance phase for the left limb is from 0 to 59%, and the swing phase is from 59% to 100%. The double support phase is from 0 to 9% and from 50 to 59%.Figure 10
The simulated relative pack displacement along the back (mean (solid red line)± one standard deviation (dashed red lines), (a)), compared with measured data (mean (solid blue line) ± one standard deviation (dashed blue lines), (a)), for six repeated trials for one subject (age: 28 yrs, weight: 71.8 kg, height: 178 cm), and measured vertical trunk CoM displacement (mean (solid blue line) ± one standard deviation (dashed blue lines), (b)). The average walking speed was 1.62±0.04 m/s, and the pack load was 23.0 kg. The stance phase for the left limb is from 0 to 60.3%, and the swing phase is from 60.3% to 100%. The double support phase is from 0 to 10.3% and from 50 to 60.3%.Figure 11
The simulated pack interaction forces in the pack coordinate system (mean (solid red line)± one standard deviation (dashed red lines)): (a) longitudinal pack force, (b) normal pack force (c), and lateral pack force, compared with measured data (mean (solid blue line) ± one standard deviation (dashed blue lines)), for six repeated trials for one subject (age: 28 yrs, weight: 71.8 kg, height: 178 cm). The average walking speed was 1.62±0.04 m/s, and the pack load was 23.0 kg. The stance phase for the left limb is from 0 to 60.3%, and the swing phase is from 60.3% to 100%. The double support phase is from 0 to 10.3% and from 50 to 60.3%.


## 4. Discussion

A 3D load carriage model has been proposed, which simulates the dynamic interactions between pack and bearer. This considers pack motions that are directly coupled to the bearer’s trunk movement and also the longitudinal compliance and damping in the backpack suspension. The 3D pack interaction forces and moments are determined from kinematic relationships, equations of motion, and a dynamic backpack suspension model. The parameters of the pack suspension model were identified from test data collected using a load carriage test rig [18, 19]. Output from the load carriage simulation model was compared with measurement data obtained from a motion analysis system at two different pack loads, and during both normal and fast walking.Referring to Figures6(a), 8(a), and 10(a), the simulated relative pack displacements (longitudinally along the bearer’s back) show reasonably good agreement with the measured data at both pack loads, and at both walking speeds. They reproduce a sinusoidal fluctuation about 20 mm with a very similar pattern to vertical displacement of trunk CoM (Figures 6(b), 8(b), and 10(b)), which is consistent with previous load carriage measurements [10]. However, the results appear to contain an additional higher-frequency harmonic that is not present in the measurement data, especially noticeable at the normal walking speed.Referring to Figures7(b), 9(b), and 11(b), the simulated longitudinal pack interaction forces also show reasonably good agreement with the measured data in terms of amplitudes and general trends. The largest pack forces occur during double stance because of the high accelerations at heel strike and push-off. As one would expect, the pack forces reach their minima in midstance when the vertical accelerations are negative. Both peak and mean values of the longitudinal pack forces increase significantly with increased pack load, which is consistent with the results of previous experimental load carriage studies [10]. However, the simulated longitudinal pack force also exhibits an additional harmonic which is not present in the measured data.It is likely that the additional harmonic seen in the longitudinal results is due to the limitations of the load carriage test-rig [18, 19]. Despite the fact that the mannequin was covered in a neoprene-like material to mimic the soft tissues, it is unlikely that there was sufficient compliance; resulting in a suspension model with a high natural frequency. Furthermore, neither the mannequin nor the simulation model can alter torso posture in response to changes in pack interaction forces. In practise, the bearer may dynamically adjust their posture to damp the pack motion and thereby decrease high-frequency fluctuations in pack force. Furthermore, the test rig only generates vertical displacement of the mannequin, rather than the true 3D trunk motion seen during walking. Hence, the pack suspension model is limited to one degree of freedom and the suspension model parameters do not properly account for coupling with other motions.Although the pack model neglects all relative pack motions except for longitudinal translation up and down the back, the simulated normal pack forces (Figures7(a), 9(a), and 11(a)) were in reasonable agreement with the experimental data with the exception of unexplained discrepancies around 40% to 60% of the gait cycle. In percentage terms, there are major discrepancies between the model predictions and the measurement data in the lateral pack forces (Figures 7(c), 9(c), and 11(c)). However, these forces are of small amplitude and the experimental data exhibit significant variability.In conclusion, whilst promising results have been achieved, the proposed methodology still needs some improvements especially the technique to identify pack suspension parameters [18, 19]. Although the load carriage test-rig has the advantage that full non-linear frequency response testing can be conducted, it does not adequately represent the effective compliance of the torso. Therefore, we suggest that future work should focus on techniques for identifying backpack suspension models by combining both load carriage test-rig data and in vivo data captured in the gait laboratory. This would have the advantage that the small but significant rotations (in the plane of the back) and lateral translations of the pack, relative to the torso, could be taken into account. Including these relative motions, in addition to the longitudinal translation, may further improve the accuracy of the predicted pack forces.



---

*Source: 349638-2013-04-29.xml*