Biomechanical Analysis of Sit To Stand Movement in Normal and Obese Subjects
Biomechanical Analysis of Sit To Stand Movement in Normal and Obese Subjects
www.elsevier.com/locate/clinbiomech
a
Department of Bioengineering, Politecnico di Milano, p.zza Leonardo da Vinci 32, 20133 Milano, Italy
b
Ospedale S. Giuseppe, Piancavallo, Verbania, Italy
Received 16 January 2003; accepted 17 June 2003
Abstract
Objective. Main purpose of this study was to develop a biomechanical model for the analysis of sit-to-stand movement in normal
and obese subjects.
Design. A biomechanical model describing sit-to-stand was developed using kinetic and kinematic experimental data. Trunk
flexion, feet movement, knee and hip joint torques were assumed as sensible indexes to discriminate between normal and obese
subjects.
Background. Sit-to-stand is a functional task that may become difficult for certain patients. The analysis of its execution provides
useful biomechanical information on the motor ability of selected subjects.
Methods. Sit-to-stand was recorded using an optoelectronic system and a force platform in 40 obese patients and 10 normal
subjects. A biomechanical model was developed using inverse dynamics equations.
Results. Kinematic and kinetic indexes evidenced differences in motion strategy between normal and obese subjects. Obese
subjects rise from the chair limiting trunk flexion (mean value: 73.1°) and moving their feet backwards from initial position (mean
deviation: 50 mm). Normal subjects, instead, show a higher trunk flexion (mean value: 49.2°, a lower angular value between trunk
and the horizontal means increased flexion) and fixed feet position (mean deviation: 5 mm). As for kinetics, obese patients show knee
joint torque higher than hip torque (maximum knee torque: 0.75 Nm/kg; maximum hip torque: 0.59 Nm/kg), while normal subjects
show opposite behaviour (maximum knee torque: 0.38 Nm/kg; maximum hip torque: 0.98 Nm/kg).
Relevance
We found differences in motion strategy between normal and obese subjects performing sit-to-stand movement, which may be
used to plan and evaluate rehabilitative treatments.
Ó 2003 Elsevier Ltd. All rights reserved.
2. Methods
modelled through a hemisphere positioned on the ab- Fourth segment––FOOT (Fig. 2):
domen. The hemisphere was supposed to be homoge-
neous and fat density value was set at 0.99 g/cm3 . An Mankle þ ðD G4 ÞKðP4 Þ þ ðD P ÞKðFy Þ ¼ 0 ð8Þ
ideal BMI value was fixed at 23 kg/m2 (this value coin-
cides with the mean value of CG); from this value an where Mankle ¼ unknown moment; P4 ¼ weight vec-
‘‘ideal weight’’ for each obese subject was calculated as tor ¼ mðfootÞ g, ðD G4 Þ ¼ G4 centre of mass position
follows: 23 ¼ ideal weight [kg]/(effective height2 ) [m2 ] and with respect to D, ðD P Þ ¼ P ground reaction position
the difference between the effective weight of each obese with respect to D, Fy ¼ vertical component of ground
patient and her ideal weight represented the fat mass that reaction force, inertial vector not significant.
had to be modelled through the hemisphere. Since all our Third segment––LEG (Fig. 2):
obese patients are women, the fat mass was distributed Mknee Mankle þ ðD G3 ÞKðP3 Þ
40% on the abdomen and 60% on the thighs (‘‘pear
shape’’ distribution); therefore only 40% of the calcu- þ ðD CÞKðRLEGx þ RLEGy Þ ¼ 0 ð9Þ
lated fat mass was considered to calculate the hemisphere
where Mknee ¼ unknown knee moment, Mankle ¼ ankle
dimension. From mass and density values, the radius ðRÞ
moment, P3 ¼ weight vector of the third seg-
of the hemisphere for each obese patient was calculated
ment ¼ mðlegÞ g, ðD G3 Þ ¼ G3 centre of mass position
as follows: R ¼ (mass (g)/density [g/cm3 ])1=3 and from R
with respect to D, ðD CÞ ¼ C position with respect to
value, the position of the fat centre of mass Gs with re-
D, RLEGx and RLEGy ¼ horizontal and vertical component
spect to hip joint was extracted. Then, inertial moments
of the reaction force on C, inertial vector not significant.
with respect to the axial centre of mass, which is per-
pendicular to the sagittal plane, were evaluated to write
the equilibrium equations of the model. Anthropometric
parameters came from the measures obtained by Winter
(1979). Each segment length was calculated using the
distance between the sagittal projections of the corre-
spondent markers coordinates. The origin of the abso-
lute reference system was positioned on the left
malleolus, horizontal and vertical components of HAT-
segment centre of mass acceleration were calculated
using the following equations (Fig. 1):
XG1 ¼ L3 cos c L2 cos b þ d1 cos a ð1Þ
YG1 ¼ L3 sin c þ L2 sin b þ d1 sin a ð2Þ
By deriving
X€G1 ¼ aG1 x
c þ L2 cos bb_ 2
¼ L3 cos cc_ 2 L3 sin c€
þ L2 sin bb€ d1 cos aa_ 2 d1 sin a€
a ð3Þ
Y€G1 ¼ aG1 y
¼ L3 cos c€c L3 sin cc_ 2 þ L2 cos bb€
þ L2 sin bb_ 2 d1 cos a€a d1 sin aa_ 2 ð4Þ
First segment––HAT (Fig. 2): calculated in the first and last trial of the experimental
session for each group. Significance level was set at P ¼
Mhip þ ðB G1 ÞKðP1 Þ þ ðB G1 ÞKðF1x þ F1y Þ JG1 a€ ¼ 0
0:05.
ð10Þ
where Mhip ¼ unknown hip moment, P1 ¼ weight vec- 3. Results
tor ¼ mðHATÞ g ¼ ðmðtrunkÞ þ mðheadÞ þ mðupper-
limbsÞÞ g, ðB G1 Þ ¼ centre of mass G1 position with 3.1. Symmetry analysis
respect to B, F1x and F1y ¼ horizontal and vertical com-
ponent of the inertial vector ðm1 aG1 Þ, JG1 ¼ centre of Differences between left and right side anterior–pos-
mass inertial moment with respect to the perpendicular terior displacement of malleolus markers were calcu-
from the sagittal plane. lated for each analysed subject of both groups. For all
First segment obese––HAT modified adding the vir- subjects we found a difference <7 mm [mean diff (SD):
tual belly 5.3 mm (1.4)].
Mhip þ ðB G1 ÞKðP1 Þ þ ðB G1 ÞKðF1x þ F1y Þ JG1 a€
3.2. Kinematic results
þ ðB Gs ÞKðPGs Þ þ ðB Gs ÞKðF1sx þ F1sy Þ
JGs a€ ¼ 0 ð11Þ Control group of healthy subjects. 90% of the analysed
CG subjects use a rising strategy characterised by a high
where (B Gs ) ¼ centre of mass Gs position (of the vir-
degree of trunk flexion and a feet movement near to zero
tual belly) with respect to B, PGs ¼ weight vector of the
mm (see Table 1 for numerical data). In Fig. 3 the Y
fat mass ¼ mðfat massÞ g, F1sx and F1sy ¼ inertial vector
coordinate of shoulder marker is plotted vs movement
horizontal and vertical component ðms aGs Þ [fat mass],
cycle to give a qualitative idea of the movement strategy.
JGs ¼ centre of mass inertial moment with respect to the
Comparison between the first and the tenth trial shows a
perpendicular from the sagittal plane [fat mass].
slight decrease of trunk flexion (that means an increased
All calculated joint moments were normalised divid-
angle value), but again no feet movement in 100% of the
ing them by height and weight of each subject to permit
analysed cases. No differences between men and women
inter-subject comparisons and they were analysed vs
were found during this analysis.
movement cycle. Movement cycle was defined using a
Obese subjects. 100% of the analysed subjects limit
angle according to Millington et al. (1992). In detail,
their trunk flexion. This strategy (In Fig. 4 the Y coor-
motion cycle begins when a decreases of 0.5° in 20 ms
dinate of shoulder marker is plotted) results different
and it ends when the variations of a value remain within
from that adopted by CG subjects (Fig. 3). A feet
0.5° and the analysed subject is in standing position.
movement backwards from the initial position is always
Standing position was evaluated considering the vertical
to observe (Fig. 4). The comparison between first and
coordinate of shoulder markers. Knee and hip joint
last trial does not evidence any change in motion strat-
moments were plotted in graphs and their maximum
egy during the experimental session. All data relative to
values were extracted for the first and last trial of each
kinematic analysis of both groups are summarised in
subject. Two more indexes representing fatigue ðM %Þ
Table 1. t-Test performed on CG vs obese group for all
and work ðP %Þ of each joint were calculated (Nigg and
kinematic parameters confirmed statistical significance
Herzog, 1995)
of the results ðP < 0:05Þ.
R
motioncycle
ðMi Þ2 dt
Mi % ¼ R P3 100 ð12Þ 3.3. Kinetic results
2
motioncycle k¼1 Mk dt
R 2 The analysis was mainly focused on the torques
ðP Þ dt
motioncycle i
Pi % ¼ R P3 100 ð13Þ generated at lower limb joints. Particular attention was
2
motioncycle k¼1 Pk dt given to knee and hip joint moments.
where Pi ¼ Mi wi , Pi ¼ muscular power in the ith joint,
Mi ¼ moment at the ith joint, wi ¼ angular velocity for Table 1
the ith joint. Kinematic results: ðaÞ angle, X coordinate deviation of malleola
markers: mean (SD)
2.4. Statistical analysis Mean a (°) X coord. deviation
of the malleola (mm)
Statistical analysis of the results was performed using Control group (1st trial) 49.2 (0.5) 5 (1.7)
StudentÕs t-test between CG and obese group parameters Control group (10th trial) 56.1 (0.1) 5 (1.5)
and considering separately the first and last trial of STS. Obese group (1st trial) 73.1 (1.2) 50 (2.3)
Obese group (10th trial) 75.3 (2.0) 47.2 (1.6)
The same test was also used to compare the parameters
F. Sibella et al. / Clinical Biomechanics 18 (2003) 745–750 749
Fig. 4. Obese group rising strategy: Y shoulder coordinate (left) and X Fig. 5. CG and obese group rising strategy: hip (light grey) and knee
malleolus coordinate (left). (black) joint moments.
750 F. Sibella et al. / Clinical Biomechanics 18 (2003) 745–750
Table 2
Kinetic results for control and obese group: mean (SD)
Mean Mhip maximum Mean Mknee maximum M % hip M % knee P % hip P % knee
(Nm/kg) (Nm/kg)
Control group (1 trial) 0.98 (0.04) 0.38 (0.08) 76.6 11.7 62.5 24.7
Control group (10th trial) 0.88 (0.02) 0.45 (0.03) 70.3 17.5 57.5 36.7
Obese group (1st trial) 0.59 (0.04) 0.75 (0.08) 52.5 35.2 33.8 64.6
Obese group (10th trial) 0.57 (0.06) 0.77 (0.08) 50.3 36.4 33.5 65.1
4. Discussion 5. Conclusions