SciELO - Scientific Electronic Library Online

 
vol.15 issue1Health risk behaviors Project (COMPAC) in youth of the Santa Catarina State, Brazil: ethics and methodological aspectsComparison of worker's health promotion interventions: cluster randomized controlled trial author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Article

Indicators

Related links

  • Have no similar articlesSimilars in SciELO

Share


Revista Brasileira de Cineantropometria & Desempenho Humano

On-line version ISSN 1980-0037

Rev. bras. cineantropom. desempenho hum. vol.15 no.1 Florianópolis Jan./Feb. 2013

http://dx.doi.org/10.5007/1980-0037.2013v15n1p16 

ORIGINAL ARTICLE

 

Ground reaction force and electromyographic activity of transfemoral amputee gait: a case series

 

Força de Reação do Solo e atividade eletromiográfica da marcha de amputados transfemorais: uma série de casos

 

 

Alex Sandra Oliveira de CerqueiraI,II; Edward Yuji YamagutiII; Luis MochizukiIII; Alberto Carlos AmadioII; Júlio Cerca SerrãoII

IUniversity of Taubaté. Department of Physical Therapy. Taubaté, SP., Brazil
IIUniversity of São Paulo. School of Physical Education and Sport. São Paulo, SP., Brazil
IIIUniversity of São Paulo. School of Arts, Science and Humanities. São Paulo, SP., Brazil

Corresponding author

 

 


ABSTRACT

Ground reaction forces (GRF) and electromyographic activity form a part of the descriptive data that characterise the biomechanics of gait. The research of these parameters is important in establishing gait training and understanding the impact of amputation and prosthetic components on movement during the act of walking. Therefore, this case series describes the GRF and electromyographic activity in the gait of transfemoral amputees. A force plate was used to measure GRF, and an electromyographic system monitored the vastus lateralis, biceps femoris, tibialis anterior, and gastrocnemius lateralis muscles of the non-amputated leg. The average vertical and anteroposterior GRF time-curves, average electromyographic activity, and descriptor variables were then analysed. We observed decreases in vertical and anteroposterior GRF magnitudes as well as in anteroposterior GRF descriptor variables during the propulsive phase in the amputated leg. There were increases in phasic muscle activity and co-activation in the non-amputated leg. We concluded that, during walking, the unilateral transfemoral amputees (who were analysed in this case series) developed lower GRF in the amputated limb and a longer period of electromyographic activity in the non-amputated limb.

Key words: Amputation; Biomechanics; Electromyography; Gait.


RESUMO

O comportamento da Força de Reação do Solo (FRS) e a atividade eletromiográfica formam uma parte dos dados que caracterizam a biomecânica da marcha. O estudo destes parâmetros é importante para a recuperação da locomoção e para compreensão do impacto da amputação e dos componentes protéticos nos movimentos desenvolvidos no andar. Portanto, esta série de casos tem como objetivo descrever a atividade eletromiográfica e a FRS de amputados transfemorais. Para mensurar a FRS, foi utilizada uma plataforma de força e um sistema de eletromiografia monitorou os músculos vasto lateral, bíceps femoral, tibial anterior e gastrocnêmio lateral da perna não-amputada. As médias das componentes vertical e ânteroposterior da FRS, a atividade eletromiográfica e variáveis descritivas foram analisadas. Foi observado uma diminuição da magnitude da FRS vertical e ânteroposterior e das variáveis descritivas da componente ânteroposterior da FRS durante a fase de propulsão na perna amputada. Houve aumento na atividade fásica muscular e co-ativação na perna não-amputada. Pode-se concluir que os amputados transfemorais unilaterais analisados nesta série de casos desenvolveram menor FRS na perna amputada e longos períodos de atividade eletromiográfica na perna não amputada durante a marcha.

Palavras-chave: Amputação; Biomecânica; Eletromiografia; Marcha.


 

 

INTRODUCTION

Amputation of the lower limbs changes the biomechanics of gait1-3. Different lower limb inertial properties4,5 and a limited capacity to generate internal forces and torques2,6 are two major locomotion problems facing an amputee when using a prosthetic limb. As a result of transfemoral amputation, there is an attenuation of the GRF in the amputated limb (AL)7. To control the sudden prosthetic knee flexion, the gait speed slows, the extension of the prosthetic knee is maintained for up to 40% of the support phase5, and the time sequence of muscle activation remains the same, compared to normal gait3, but lasts longer. Moreover, the non-amputated leg (NAL) muscles, especially the hip extensors and ankle plantar flexors2, generate more joint torque and power to move the body forward.

Previous studies have explored the kinematics5,7-9 and kinetics1,2,6,8,9 of the amputee gait, but few studies have described the muscles' activation3. It is unknown how transfemoral amputation affects muscle activation during walking, nor what adaptations in muscle activation may occur to accomplish the changes observed in the mechanics of the locomotion system with a prosthetic leg. Analysis of electromyographic activity may describe some of the strategies used by the nervous system to adapt to the amputee condition. Several factors affect this adaptation such as: the amputation level1, the prosthesis type9, the muscle reinsertion method into the thigh, the anatomical and functional condition of the remaining muscles and nerves3, gait rehabilitation strategies10, stump length11, and how varied were the motor experiences after amputation.

Although the lower limb amputation incidence is not low12, because of the many differences in their own adaptation processes, the challenges of gathering several amputee participants for gait analysis are formidable. For example, differences in the related features of the amputation (etiology, amount of time, stump length and circumference) and prosthesis (type, socket, foot, time between first amputation, prosthesis placement, as well as the length of time with current prosthesis) all hinder composing a group with similar features. Although the studies seek to analyze homogeneous samples, it is difficult to standardize the prosthetic components used3,11 or to assemble a group with similar characteristics such as age and time since amputation1,2, stump length3,11, or cause of amputation9. Thus, one of the feasible strategies for studying amputee gait is to combine selected cases and analyze what they have in common. Such a "within group" analysis may help to understand the biomechanics of the issue.

The objective of this study, then, is to describe, during the walking, the GRF and electromyographic activity of three transfemoral amputees. The experimental hypothesis is that GRF and EMG parameters change according to the time elapsed since amputation.

Such research is important in establishing gait training as well as in understanding the impact of amputation and prosthetic components on movement during walking.

 

METHODS

The ethical committee of School of Physical Education and Sport at the University of São Paulo (USP) approved this study (protocol No.37), and each participant signed an informed consent. The participants were three unilateral transfemoral male amputees of traumatic origin.

Participant 1 (P1: 13 years old, 1.66 m, 45.5 kg, time since amputation: 7.3 years, level: proximal thigh, prosthesis experience: 1 month) wore an endoskeletal prosthesis with an ischial total contact containment socket with suction suspension, a Tehlin knee (TK4POC, pneumatic control in the swing phase), and a Springlite foot for one month. Because of a short stump, P1 had problems finding a comfortable prosthesis. After several attempts to fit a prosthesis, the participant began a one-month period of gait training. After his amputation, he played table tennis for 2.4 years on crutches, and for one month prior to beginning this study he played while using his prosthesis (training: 3-5 times/week, 2-4 h/day).

Participant 2 (P2: 30 years old, 1.77 m, 79.4 kg, time since amputation: 1 year, level: mid-thigh, prosthesis experience: 6 months) changed his prosthesis one month prior to this study and is currently training with the new prosthesis. The current prosthesis has an ischial containment total contact socket with suction suspension, a Tehlin knee (TGK4000-mechanical control in the swing phase) and a SACH foot. He practiced discus throwing and putting the shot for two years prior to his amputation (training: 3 times/week, 2-4 h/day), and had restarted his training routine one month before beginning this study.

Participant 3 (P3, 17 years old, 1.75 m, 76 kg, time since amputation: 3 years, level: mid-thigh, prosthesis experience: 2.6 years) wore an endoskeletal prosthetic with an ischial containment total contact socket with suction suspension, a Proteval knee (pneumatic control in the swing phase), and an Endolite foot. He was successfully fitted with prosthesis and has completely adapted to walking and running. He has practiced table tennis with his current prosthesis for 2.5 years (training: 3-5 times/week, 2-4 h/day).

Procedure

A piezoelectric force plate (600 x 900 mm, Kistler 9287A), placed in the middle of a 20-m walkway, measured the anteroposterior and vertical GRF. The force plate and the walkway were covered with a 2 x 20 m nonelastic plastic carpet. A raw electromyography (EMG) signal was recorded (Bagnoli-8 - Delsys, Inc., Boston, MA) using a differential amplification (amplified 1,000 times at a 12-bit resolution) with a sampling frequency of 1 kHz. The EMG bandwidth was limited to between 20 and 450 Hz (CMRR < 92 dB, input impedance >1015/0.2 ohm/pF). An analog-to-digital converter (A/D DAS - 1600/1400 Series Keithley Instruments, Inc.) with 16 channels and 12 bit resolution was responsible for data synchronization.

Only the non-amputated leg muscle activities were monitored. After the trichotomy, the bipolar active surface electrodes (Ag/AgCl; 1 cm diameter, 1 cm inter-electrode distance) were placed 1 cm away from the motor point13 of the m. vastus lateralis (VL), m. biceps femoris (BF), m. tibialis anterioris (TA), and m. gastrocnemius lateralis (GL). To locate the motor point the participants were asked to lie down on a clinical table. We then applied electrical pulses using an OMNI Pulsi-901 (Quark, Piracicaba, SP, Brazil) universal pulse generator on the surface of the skin and above the muscle where even the smallest intensity of current would activate the muscle. The monophasic, quadratic, pulsed current was applied with 7 Hz of frequency - the smallest intensity to activate the motor point. A ground electrode was placed over the patella.

The participants walked straight ahead at a self-selected speed. The speed was determined to be the average velocity to cross a 20 m long walkway. The time was measured with a manual chronometer. The average velocities were 0.63±0.03 m/s (P1), 0.84±0.02 m/s (P2), and 1.08±0.08 m/s (P3). The participants were expected to step over the force plate 15 times with each foot; however, for some trials, the participant did not step correctly over the force plate, resulting in discarded data. Ten AL and nine NAL stances were measured for P1, and twelve AL and eight NAL stances for P2. Finally, 10 AL and 15 NAL stances were measured for P3.

Data analysis

The raw GRF was low-pass filtered (2nd order Butterworth 20 Hz, recursive filter) and normalized in relation to the body weight (BW)14 of the participant. The raw EMG was represented through a linear envelope15 and calculated in five steps: off-set removal from the raw EMG, full-wave rectified, low-pass filtered (2nd order Butterworth 5 Hz, recursive filter), the data amplitude normalized in relation to its mean16, and the time base normalization by the stance phase (% SP)14. The beginning and end of the stance phase was determined by the vF. The cutting was done visually and determined during the processing of data.

The following variables were calculated (Figure 1a) from the anteroposterior GRF (apF): the braking phase peak (1apF) and its instant (Δt1apF); the propulsive phase peak (2apF) and its instant (Δt2apF); the braking phase impulse (apFBimp); the propulsive phase impulse (apFPimp); the ratio between the impulses (apFBP); and the stance time (Δtstance). The total vertical impulse (vFimp)14 was derived from the vertical GRF (vF). Both the pulse duration (the first and last instants the EMG intensity was > 25% peak, initial-final% SP Figure 1b), and the peak instant (from 0 to 100% SP, Figure 1b) were calculated from the EMG. We opted to choose a moderate level of muscle action17 as a parameter indicative of phasic activity during the stance phase. The linear envelope and variables were calculated using a mathematical function (Matlab software).

 


 

The figure illustrates the following variables: the braking phase peak (1apF) and its time interval (Δt1apF); the propulsive phase peak (2apF) and its time interval (Δt2apF); the braking phase impulse (apFBimp); the propulsive phase impulse (apFPimp); the ratio between the impulses (apFBP); and the stance time (Δtstance), which was calculated from the anteroposterior GRF. b) The average of the gastrocnemius lateralis from S3, which represents the studied variables. The EMG intensity, which was established to determine the duration of the pulses (25% of the peak), are represented using the dash dot (-.-). The first and second (start and final) pulses are represented using the '*' symbol.

The averages and standard deviations for all those parameters were described and compared across the three participants. Data reported by Rab18and Winter19 were used to compare our results to a pattern of non-amputee walking.

 

RESULTS

The participants presented different and asymmetrical vF ensemble averages (Figure 2). Therefore, the peaks and inclinations from the vF were not calculated. As a consequence of slow gait and shorter Δtstance, the vertical GRF impulse was lower for the AL (Table 1).

For apF ensemble average curve, all participants presented a biphasic pattern for both limbs (Figure 2); but their parameters were different (Table 1).

During the braking phase, the 1apF and apFBimp were similar for both stances of Participants 1 and 3, whereas in the propulsive phase, the 2apF, Δt2apF, apFPimp, and Δtstance were lower in the AL of all participants (Table 1).

For P2 and P3, the VL was active from the beginning of the stance phase (Figure 3, Table 2) up to 50 and 30% SP (Table 2), respectively. For all three participants (Figure 3, Table 2), the BF activity begins at the weight acceptance (where it reaches its peak), and the TA was active during the weight acceptance and pre-swing. For P2 and P3 (Figure 3, Table 2), the GL presented its first burst (just after the foot-flat phase) to decelerate the tibia rotation, and the second during the propulsive phase. Only for P1, were all muscles active after 60% of the stance phase (Figure 3, Table 2).

 

 

DISCUSSION

The purpose of this case series was to describe the GRF and electromyographic activity of three transfemoral amputees during the act of walking. The combination of dissimilar individual characteristics and their differing gait kinematics led to important adaptations in both the GRF and their phasic muscular activity.

Individual differences among participants and how they coordinate the AL and the NAL might affect vF ensemble averages. The GRF is the most commonly studied external force. This force represents the pattern of acceleration of the whole body center of mass, which is formed by vertical, anteroposterior, and medial-lateral vector-components. During gait, the vF typically has two peaks: the first vF peak occurs immediately after the heel strike and represents the deceleration of segments at the beginning of the stance phase, and the second determines the acceleration upward from the center of mass during the push-off phase14,20. Participant 3 was the faster walker, the one to wear his prosthesis longest, and the only one who exhibited the two typical peaks (Figure 2)14.

The other two participants, on the other hand, who began their gait training only one month prior to the study, presented slower gait speed and asymmetrical vF (Figure 2). The gait is a dynamic activity, and the gait speed is associated with the force applied to the ground. The individual characteristics, such as the condition of the stump, the duration of the amputation procedure, the time it takes to become accustomed to the prosthesis, and the time devoted to gait training21,22 may all have contributed to the slower gait speed.

As a consequence of slow gait and the shorter Δtstance (Table 1), the vertical GRF impulse was lower for the AL (Table 1), thereby mirroring previous studies1,23. The lowest vertical impulse facilitates balance control24. Moreover, after amputation, the center of gravity moves nearer to the NAL rather than the AL1,25. Lower limb amputees prefer to load their weight on the NAL1. After lower limb amputations, the somatosensory input is impaired, and these balance strategies are developed by transfemoral amputees.

The influence of amputation in apF was notable during propulsion. The lack of the plantar flexor muscles and the greater prosthetic knee extension during the weight acceptance stage5,24 may be responsible for a decrease in forward propulsion. Only in P2, is there an attenuation of 1apF and apFBimp. This behaviour is possibly an attempt to minimize the loading rate. This result is a response to a combination of factors, such as the reduction of walking speed and shorter stance time developed by the participating volunteers.

Muscular activation varied among different participants. The difference suggests that motor strategies account for how comfortable each amputee is with his own prosthesis.

In normal gait18,19, the major activity and the peak of the VL occurs at the weight acceptance to control the knee flexion and to help the knee extension during the mid-stance. BF activity begins during the terminal swing acting to decelerate the swing leg and continues into the weight acceptance stage. It is at this point that it reaches its peak. TA activation begins at the terminal swing for the ankle dorsiflexion position and controls the foot-flat phase after the heel strike. During pre-swing, the TA activates the ankle dorsiflexion for foot clearance. The GL presented its first burst just after the foot-flat phase to decelerate the tibia rotation. During the propulsive phase, the muscle generates the highest mechanical power to plantar flex.

During weight acceptance, Participants 2 and 3 presented longer VL, BF, and TA bursts; they presented longer GL burst during the propulsive phase. To increase the phasic muscle activity is a motor strategy observed in transtibial26 and transfemoral3,8 amputees. The unilateral amputation increases the net joint moments and the power output on the NAL1. For example, in non-amputee gait; the GL accelerates the body forward during the terminal stance19. The absence of the forward push-off on the prosthetic leg requires a higher amount of power on the NAL27,28, thereby increasing the duration of muscles bursts. Furthermore, the longer EMG bursts, co-activations, and the reduction of vF are all strategies for achieving balance control. The longer co-activation periods increase the joint stiffness and prevent knee collapse at the load response. However, prolonged coactivation affects the mechanical and metabolic efficiencies of movement which might cause muscular fatigue18.

The sample size restricted any generalization of our findings. However, it did allow us to explore the individual trends in detail29. Because of the uniqueness of each lower-limb amputation and participants' abilities with their custom-made prostheses, this group is likely to exhibit higher inter-individual variability in comparison to adults without amputations.

One other limitation of this study is related to how to determine the beginning and end of the stance phase, which was addressed visually during the processing of data.

For the next study we suggest recording the EMG in the AL. We analyzed the muscles in the NAL because of the difficulty in establishing the correct placement for the EMG electrodes on the stump muscles. A recent study placed the electrodes at alternative locations providing strong EMG signals from both the flexor and extensor muscles of the stump. This placement thereby facilitated their being recorded30. This is a strategy that might be followed in future studies.

 

CONCLUSION

The transfemoral amputees analyzed in this case series exhibited atypical vertical GRF in their AL and reduced anteroposterior GRF during the propulsion phase. The phasic muscle activity in their NAL was increased in comparison to non-amputee walking. The degree of the gait development, the actual speed at which the walking occurred, and the individual characteristics of the studied participants all influenced the results.

Further work is required to establish whether these characteristics are found generally in transfemoral amputees.

 

REFERENCES

1. Nolan L, Wit A, Dudziñski K, Lees A, Lake M, Wychowañski M. Adjustments in gait symmetry with walking speed in transfemoral and trans-tibial amputees. Gait Posture 2003;17(2):142-51.         [ Links ]

2. Seroussi RE, Gitter A, Czerniecki JM, Weaver K. Mechanical work adaptation of above-knee amputee ambulation. Arch Phys Med Rehabil 1996;77(11):1209-14.         [ Links ]

3. Jaegers SMHJ, Arendzen, JH, de Jongh HJ. An electromyographic study of the hip muscles of transfemoral amputees in walking. Clin Orthop Relat Res 1996;328:119-28.         [ Links ]

4. Gitter A, Czerniecki J, Meinders M. Effect of prosthetic mass on swing phase work during above-knee amputee ambulation. Arch Phys Med Rehabil 1997;76:114-21.         [ Links ]

5. Selles RW, Bussmann JB, Wagenaar RC, Stam HJ. Effects of prosthetic mass and mass distribuition on kinematics and energetics of prosthetic gait: a systematic review. Arch Phys Med Rehabil 1999;80(12):1593-9.         [ Links ]

6. Van der Linden ML, Solomonidis SE, Spence WD, Ning Li, Paul JP. A methodology for studying the effects of various types of prosthetic feet on the biomechanics of transfemoral amputee gait. J Biomech 1999;32(9):877-89.         [ Links ]

7. Klotz R, Colobert B, Botino M, Permentiers I. Influence of different types of sockets on the range of motion of the hip joint by the transfemoral amputee. Ann Phys Rehabil Med 2011;54(7):399-410.         [ Links ]

8. Hong HH, Mu SM, Relationship between socket pressure and EMG of two muscles in transfemoral stumps during gait. Prosthet Orthot Int 2005;29(1):59-72.         [ Links ]

9. Johansson JL, Sherrill DM, Riley PO, Bonato P, Herr H. A clinical comparison of variable-damping and mechanically passive prosthetic knee devices. Am J Phys Med Rehabil 2005;84:563-75.         [ Links ]

10. Baker PA, Hewison SR. Gait recovery pattern of unilateral lower limb amputees during rehabilitation. Prosthet Orthot Int 1990;14(2):80-4.         [ Links ]

11. Mizuno N, Aoyama T, Nakajima A, Kasahara T, Takami K. Functional evaluation by gait analysis of various ankle-foot assemblies used below-knee amputees. Prosthet Orthot Int 1992;16(3):174-82.         [ Links ]

12. Stineman MG, Kwong PL, Xie D, Kurichi JE, Ripley DC, Brooks DM, et al. Prognostic Differences for Functional Recovery after Major Lower Limb Amputation: Effects of the Timing and Type of Inpatient Rehabilitation Services in the Veterans Health Administration. PM & R 2010;2(4):232-43.         [ Links ]

13. Robinson AJ. Neuromuscular electrical stimulation for control of posture and movement. In: Robinson AJ, Snyder-Mackler. L. Clinical electrophysiology electrotherapy and electrophysiology testing. 1st ed. New York: Willians & Wilkins; 1989. p. 157-210.         [ Links ]

14. Chao EY, Laughman RK, Schneider E, Stauffer RN. Normative data of knee joint motion and ground reaction forces in adult level walking. J Biomech 1983;16(3):219-33.         [ Links ]

15. Winter DA. Biomechanics and motor control of human movement. 2nd ed. New York: John Wiley & Sons. 1990.         [ Links ]

16. Burden AM, Trew M, Baltzoupolos V. Normalization of gait EMGs: a re-examination. J Electromyo Kinesiol 2003;13(6):519-32.         [ Links ]

17. Brennecke A, Guimarães TM, Leone R, Cadarci M, Mochizuki L, Simão R, Amadio AC, Serrão JC. Neuromuscular activity during bench press exercise performed with and without the preexhaustion method. J Strength Cond Res 2009;23(7):1933-40.         [ Links ]

18. Rab GT. Músculos. In: Rose J, Gamble J, editors. Marcha humana. 2nd ed. São Paulo: Premier; 1998. p.107-28.         [ Links ]

19. Winter DA. The biomechanics and motor control of human gait: normal. Elderly and pathological. Ontario: University of Waterloo Press. 1991.         [ Links ]

20. Nilsson J, Thorstensson A. Ground reaction forces at different speeds of human walking and running. Acta Physiol Scand 1989;136(2):217-27.         [ Links ]

21. Boonstra AM, Schrama J, Fidler V, Eisma WH. The gait of unilateral transfemoral amputees. Scand J Rehabil Med 1994;26(4):217-23.         [ Links ]

22. Sjödahl C, Jarnlo GB, Söderberg B, Persson BM. Kinematic and kinetic gait analysis in the sagittal plane of transfemoral amputees before and after special gait reeducation. Prosthet Orthot Int 2002;26(2):101-12.         [ Links ]

23. Weyand PG, Bundle MW, McGowan C, Grabowski A, Brown MB, Kram R, Herr H. The fastest runner on artificial legs: different limbs, similar function? J Appl Physiol 2009;107(3):903-11.         [ Links ]

24. Murray MP, Molinger LA, Sepic SB, Gardner GM, Linder MT. Gait patterns in above-knee amputees patients: hydraulic swing control vs Constant-friction Knee components. Arch Phys Med Rehabil 1983;64(5):339-45.         [ Links ]

25. Vrieling AH, van Keeken HG, Schoppen T, Otten E, Hof AL, Halbertsma JP. Postema K. Balance control on a moving platform in unilateral lower limb amputees. Gait Posture 2008;28(2):222-8.         [ Links ]

26. Soares AS, Yamaguti EY, Mochizuki L, Amadio AC, Serrão JC. Biomechanical parameters of gait among transtibial amputees: a review. São Paulo Med J. 2009;127(5):302-9.         [ Links ]

27. Tesio L, Lanzi D, Detrembleur C. The 3-D motion of the center of gravity of the human body during level walking. II. Lower limb amputees. Clin Biomech (Bristol. Avon). 1998;13(2):83-90.         [ Links ]

28. Sagawa YJr., Turcot K, Armand S, Thevenon A, Vuillerme N, Watelain E. Biomechanics and physiological parameters during gait in lower-limb amputees: a systematic review. Gait Posture 2011;33(4):511-26.         [ Links ]

29. Glazier PS, Davids K. Constraints on the complete optimization of human motion. Sports Med 2009;39(1):15-28.         [ Links ]

30. Zhang F, D'Andrea SE, Nunnery MJ, Kay SM, Huang H. Towards design of a stumble detection system for artificial legs. IEEE Trans Neural Syst Rehabil Eng 2011;19(5):567-77.         [ Links ]

 

 

Corresponding author
Alex Sandra Oliveira de Cerqueira
Laboratório de Biomecânica - Escola de Educação Física e Esporte Universidade de São Paulo (USP)
Av. Professor Mello Moraes, 65
CEP 05508-900 - São Paulo (SP). Brasil
E-mail: lesoares@usp.br

Received: 15 March 2012
Accepted: 24 July 2012