Muscular strength is a valuable attribute to perform many sports and simple day-to-day activities. Thus, the assessment of muscular strength is essential for understanding the performance capacity of an individual (Bottaro et al., 2005). Commercially available isokinetic machines have created lots of clinical application for injury rehabilitation, measurements of muscular torque, work, power, or endurance. However, many internal and external factors in the isokinetic testing procedures can have an undesirable effect on the test results.

The isokinetic dynamometer is a device that controls the velocity of an exercising limb by means of a preset speed-governing mechanism (Brown, 2000). As more force is exerted against the lever arm of the dynamometer, the energy of the moving limb is absorbed by the apparatus and returned as additional resistance to the limb movement (Brown et al., 1995a; 1995b). Movement, therefore, occurs at a constant predetermined speed in the range of motion referred to as load range (LR) (Brown, 2000) where the dynamometer imposes an external load to the limb movement. However, an exercising limb must free accelerate (Sapega et al., 1982) to the preset speed and decelerate at the end of the exercising range. These phases are performed without the benefit of externally imposed resistance and, consequently, should not be considered during test interpretation (Brown, 2000). But the absence of dynamometer resistance in acceleration (ACC) and deceleration (DEC) phases is not the only reason for isokinetic data misinterpretation: velocity overshoot (VO) and impact artifact are also reported causes (Brown, 2000). The VO is a movement artifact that occurs when the limb accelerates past the desired speed and the dynamometer tries to slow the limb. This generates two effects: (a) a torque spike (or torque overshoot) as a consequence of the braking mechanism of the dynamometer (Wilk et al., 1992; Brown, 2000), and (b) a short period of time where the angular velocity oscillates until its stabilization at the preset speed (Sapega et al., 1982; Chen et al., 1994). The impact artifact has an effect similar to VO and it is observed during deceleration when the dynamometer begins to slow the lever arm in preparation for stopping at the turnaround point. This causes a torque spike at the end of the repetition due to the lever arm impacting the mechanical end stop (Brown, 2000). Both VO and impact artifact produce undesirable behavior and are usually removed with ACC and DEC phases prior to isokinetic test interpretation. Wilk et al., 1992 related the use of an "isokinetic window" which consists of removing all data that has not been obtained at the preset isokinetic speed or 95 percent of that speed. This technique has been shown to increase the reliability of testing via the control of aberrant torque production (Wilk et al., 1992). They found significant differences between windowed and nonwindowed data during isokinetic tests. Tis and Perrin, 1993 advise that using a data reduction technique that eliminates the first and last 10º of the total range of motion (ROM) may eliminate ACC and DEC areas, but this may also eliminate the peak torque range (Brown, 2000). Kurdak et al., 2005 proposed a method that calculates the first-difference of velocity and accepts fluctuations near to zero value (0 ± 0.2) as isokinetic. They have not found the VO artifact in their research. However, small fluctuations inside a VO segment could be mistakenly considered as isokinetic if only the first-difference is used. All the techniques are useful to diminish the movement artifacts effects. However, none of them is able to exactly quantify the contribution of the VO effect. In general, VO is considered as a segment of the ACC or LR phase. Actually, VO cannot be treated as an acceleration segment because it has an imposed external load. On the other hand, it cannot also be considered a LR segment since its velocity oscillates. By definition, isokinetic is constant velocity and represents a match between a mechanically imposed velocity and the subject's movement (Brown, 2000). Thus, a "pure isokinetic" analysis could not have VO within the range under study.

This study tested the influence of VO on isokinetic knee extension data interpretation through the analysis of biomechanical descriptors. A computational algorithm was developed to delineate the segments of each isokinetic repetition phase (IRP). The first-difference was associated with the windowing criterion and other mechanisms of signal analysis which make possible an accurate identification of IRPs. The VO segment was treated as a fourth IRP and its effects were reported independently of the three major phases (ACC, LR and DEC). This new approach was compared with the windowing method and the data reduction technique proposed by Tis and Perrin, 1993, referred to here only as data reduction.

Sixteen normal healthy adult males (age 26.8 ± 4.7 years) of height 1.76 ± 0.05 m and body mass 79.2 ± 9.4 kg with no history of orthopedic disease participated in this study. They voluntarily read and signed a written consent form before participating in the experiment that was approved by the University Institutional Review Board.

A Biodex System 3 Pro isokinetic dynamometer (Biodex Corp., Shirley, NY, USA) was calibrated and assembled with the knee attachment according to the manufacturer's specifications (Biodex, 1998). Biodex has been shown to be a reliable instrument for collecting data of human torque, joint position, and limb velocity (Brown et al., 1993; Drouin et al., 2004; Feiring et al., 1990; Gross et al., 1991; Ortqvist et al., 2007; Taylor et al., 1991). The dynamometer shaft was aligned with the assumed axis of rotation (lateral femoral condyle) of the dominant knee (right leg for all the subjects) with the subject in a seated position and the back reclined at approximately 110º. The left thigh was secured with straps as were the waist and thoracic torso (Weir et al., 1996). Arms were placed across the chest with hands grasping the straps (Stumbo et al., 2001). The lever arm pad was positioned to place the inferior aspect immediately superior to the medial malleolus. Subjects were passively moved to 0º of extension (full extension). After, the knee was flexed about 5º to 10º to a comfortable position set as the extension mechanical stop. Then, the flexion mechanical stop was defined so as to ensure a range of motion of 85 degrees. Gravity compensation analysis was performed by the computer system software provided with the Biodex System 3 Pro.

The biomechanical signals were acquired trough the dynamometer DB-15 female interface (Biodex, 1998) which provides real time analog signals of torque, angular velocity, and angular position. An adaptor was built by the authors in order to get the signals from the DB-15 interface into three separate BNC connectors to a digitizer board (BNC-2120, National Instruments, TX, USA) which sampled the biomechanical signals at 2048 samples·s^-1, and converted it to digital data via a 12 bit A/D. A software tool (Schwartz et al., 2008) was used to adjust the DB-15 voltage of the recorded signals to real units (N·m, degrees·s^-1, and degrees), following the manufacturer's specifications. This mechanism was used in order to get higher resolution for the digitized signals once the Biodex System 3 Pro only provides signals sampled at 100 samples·s^-1. Although the rate of 100 samples·s^-1 is sufficient for the isokinetic exercises analysis, a better resolution was chosen for a better precision of the biomechanical descriptors calculus.

Following equipment setup, subjects were asked to perform 10 gradient sub-maximal reciprocal concentric extension (240º·s^-1) and flexion (300º·s^-1) repetitions for warm-up and familiarization with the equipment. For testing, subjects performed two sets of ten maximal concentric repetitions of dominant knee extension (at 60º·s^-1 and 180º·s^-1), on separate days and in a counterbalanced order. Consistent and standard, moderate (no yelling or screaming) verbal encouragement was given. No visual feedback of the biomechanical signals was provided to subjects during the test.

In order to calculate the biomechanical descriptors on each separate IRP a computational algorithm (Figure 1) was developed by the authors in Matlab 6.5 (MathWorks, Natick, MA, USA). The algorithm divided the LR segment into two segments: VO and an isokinetic load range (ILR) segment. First, all the values less than the mean value of the extension angular velocity are set to zero and the original signal is shifted to the dashed-dot line trajectory of Figure 1a. Second, the first-difference technique (Smith, 1998) is applied to the shifted signal. This results in the circled-dots and the bold line shown in Figure 1b. Third, the absolute values of the bold line segment between the two circle-dots are determined and their mean value calculated. Figure 1c displays a zoomed image of this bold line segment and Figure 1d shows its absolute values with a straight line representing their mean value. Fourth, the algorithm investigates the first point greater than the mean value on the left and on the right side from the center of the segment, as illustrated by the asterisks in Figure 1d. These asterisks delineate the ILR segment. Finally, VO could be considered as the region between the first circled-dot and the first asterisk. However, its starting point is adjusted to the first point with at least 95 percent of the preset speed in order to meet the criterion for windowing (Wilk et al., 1992). The result is shown in the highlighted portion of Figure 1e.

Four of the most common biomechanical descriptors were determined for each separated IRP: (1) the total work (TW) which is the total amount of work that is produced in a set (Brown, 2000; Remaud et al., 2007); (2) the peak torque to body weight (PTBW) which is the maximum torque produced in a set of repetitions, normalized to body weight (Brown, 2000); (3) the time interval (TI) in seconds of each phase of the movement; (4) the average length (AL) in degrees of each phase range.

For all subjects, the four IRPs were identified in each repetition of each set of the knee extension movement at 60º·s^-1 and 180º·s^-1. The biomechanical descriptors were calculated from the total ROM and for each IRP. Then, the percent relation (PR) between the values of each IRP and the value of total ROM were established according to the following equation:

(1)

Percent relation is a measure of how much each IRP contributes with the whole value of the descriptor inside a repetition.

The same procedure was repeated for the windowing (Wilk et al., 1992) and data reduction techniques (Tis and Perrin, 1993) considering only the three major isokinetic phases, since these techniques do not manage the VO artifact.

The behavior of isokinetic movement phases was focused on the PR ratio of each descriptor. Student's t-test for dependent samples was applied to compare normally distributed data with a level of significance of 0.05 (two-tailed) and 95% of confidence interval. Wilcoxon Signed-Rank non-parametric test (De Sá, 2007) was applied to compare non-normal data.

The first set of descriptor values was calculated for the four IRPs identified by the proposed algorithm (see Table 1). A comparison between the two angular velocities (60º·s^-1 and 180º·s^-1), for each IRP, and for all of the studied descriptors revealed significant differences (p < 0.05), except for PR_PTBW (percent relation of peak torque to body weight) in ILR phase where the means are visibly equal.

The distinct behaviors of VO and ILR are shown in Table 2 by the coefficient of variance calculated for both phases from the peak torque repetition produced by the subjects with the least and the greatest peak torque. For the least peak torque, the coefficient of variance found at 60º·s^-1 (180º·s^-1) was 7.22% (8.67%) for VO and 1.02% (1.35%) for ILR. For the greatest peak torque, the values are 6.41% (6.83%) for VO and 0.70% (0.52%) for ILR. The higher fluctuations in VO reinforce its non-isokinetic nature. The two extreme cases of peak torque production were chosen in order to detect the smallest and the largest possible velocity variations in VO.

Total ROM was determined at 60º·s^-1 (82.79 ± 2.32 degrees) and 180º·s^-1 (83.25 ± 0.80 degrees) and it was consistent with the dynamometer settings.

The results from the windowing (Table 3) and data reduction (Table 4) techniques concerning LR phase were compared between them and with the results of ILR phase delineated by the proposed algorithm. Significant differences (p < 0.05) showed that the largest segments for LR were found with the windowing method, followed by the ILR segments obtained with the algorithm, and finally by the segments found with data reduction technique. The average length descriptor gives a good notion of this finding. From Tables 1, 3, and 4, LR (or ILR) contributes to the following percentage of AL at 60º·s^-1 (180º·s^-1): 95.67% (86.46%) for the windowing method, 82.97% (76.87%) for the proposed algorithm, and 75.60% (75.63%) for the data reduction technique. Table 5 summarizes these relations for TW, TI, and AL. There are no differences for the PTBW descriptor in none of the studied methods. Lastly, the LR segment found by windowing was not significantly different from the VO+ILR segment obtained with the algorithm. Figure 2 illustrates the variation of the time length segments according to the angular velocity and the identification technique used.

The most widely used technique for data adjustment in isokinetic exercises is the windowing technique which treats VO and LR as one isokinetic phase (Bottaro et al., 2005; Brown et al., 1995a; 1995b; Findley et al., 2006; Kurdak et al., 2005; Wilk et al., 1992; 1994). In this study, velocity overshoot was seen as a separate phase. Table 1 shows that VO has an important contribution to the descriptor values. Figure 1c illustrates how the velocity variation is accentuated in the VO phase when compared with the ILR phase. Table 2 reveals that the velocity variation in the VO is approximately 7 times larger than the variation in ILR at both velocities (60º·s^-1 and 180º·s^-1) when the least peak torque case is observed. For the greater peak torque case, the relation increases to approximately 9 times (60º·s^-1) and 13 times (180º·s^-1) respectively. It means that even for individuals with low peak torque production, the VO phase has a large fluctuation when compared with ILR. Therefore, it is reasonable to consider that the ILR phase is the only part where velocity could be constant.

Previous studies have focused on the importance of removing the ACC and DEC phases from the results analysis (Maly et al., 2006; Messier et al., 2005; Wilk et al., 1992; 1994). It is well understood that they are not isokinetic phases. However, keeping them in the analysis of certain descriptors causes less impact than keeping the VO phase. Looking at the behavior of total work at 60º·s^-1 (Table 1), we can see that VO contributes about 8.93% of the work realized in total ROM, whereas ACC and DEC summed are about 1.91%. For the angular velocity of 180º·s^-1, ACC and DEC phases summed contribute about 9.17% of the work in total ROM against 6.7% of VO, what is still a high value. This increase of ACC+DEC can be observed in PR_TI descriptor (3.3 times higher for ACC and 2.8 times higher for DEC) and it is related to the time necessary to reach a higher angular velocity and to stop the movement at the end of the repetition. Consequently, there was a reduction of the VO+ILR segment at 180º·s^-1 (see PR_AL descriptor in Table 1), which has been repeatedly reported for increased dynamometer speed (Brown et al., 1995a; 1995b; 1998; Lander et al., 1985; Osternig, 1975; 1986).

Comparing VO with ILR, the work produced in VO is approximately 10% (at 60º·s^-1) and 8% (at 180º·s^-1) of the work produced in ILR. Thus, an analysis of total work considering VO+ILR as an isokinetic segment could suggest that a subject produced 10% (8%) more work than he really did under isokinetic conditions at 60º·s^-1 (at 180º·s^-1).

Since the windowing technique does not distinguish VO from ILR phase, the entire load range is used for analysis (Table 3). Consequently, treating VO as an isokinetic segment significantly modifies results interpretation since velocity is not constant in the VO segment (see comparisons in Table 5). Therefore, the windowing technique does not provide appropriate conditions for a "pure isokinetic" analysis. There is not a scientific definition for what the expression "pure isokinetic" means. For the data of this work, and based on Table 2 results, it is reasonable to consider all the segments with coefficient of variance less than or equal to 1.35 (for angular velocity) as a pure isokinetic segment. However, this is not a generalization because it depends on angular velocity, peak torque produced and the dynamometer used, which deserves specific investigations.

Looking at the LR results for the data reduction technique in Table 4, they seem similar to the ILR results of Table 1, especially at 180º·s^-1. However the differences are reported in Table 5 and illustrated in Figure 2. It is easy to observe in Figure 2 that a part of the isokinetic segment is eliminated when the last 10º are subtracted from the total ROM at both angular velocities. This evidence justifies the possible loss of peak torque range related by Brown, 2000. However, no loss of peak torque was observed in the present study. The PR_PTBW descriptor indicates that peak torque occurred in the load range for all three techniques and for both angular velocities. Figure 2 also reveals that the first 10º segment would not be large enough to completely remove the VO phase in both velocities. Therefore, a portion of the analysis segment would not be isokinetic, compromising results interpretation.

The occurrence of torque spikes subsequent to ACC or at the end of the repetition was not observed. The main reason is that the dynamometer used has an impact-free acceleration and deceleration mechanism which eliminates joint trauma, allowing patients to exercise and be tested at more functional speeds (Biodex, 1998). On average, the closest value of the peak torque occurred in the DEC phase at 180º·s^-1 as shown by PR_PTBW in Table 1. This is consistent with previous reports that impact artifact and lever arm oscillation (Brown, 2000) usually occur at speeds greater than 180º·s^-1 (Gransberg and Knutsson, 1983; Handel et al., 1996).

• Isokinetic test interpretation must be focused on the constant velocity range; traditional analysis usually removes the acceleration and deceleration phases but does not give particular attention to velocity overshoot range.
• The study of effects of velocity overshoot artifact requires a specific method for accurately delineate its interval and investigate its impact over biomechanical descriptors; this paper proposed a computational algorithm for identifying the velocity overshoot interval.
• Velocity overshoot has significant impact over biomechanical descriptors analyzed during isokinetic knee extension tests at 60º·s^-1 and 180º·s^-1; the algorithm proposed is an advisable method for performing isokinetic tests analysis according to the isokinetic definition.