Abstract
While movement is essential to human wellbeing, we are still unable to reproduce the deftness and robustness of human movement in automatons or completely restore function to individuals with many types of motor impairment. To better understand how the human nervous system plans and controls movements, neuromechanists employ simple tasks such as upper extremity reaches and isometric force tasks. However, these simple tasks rarely consider impacts and may not capture aspects of motor control that arise from real-world complexity. Here we compared existing models of motor control with the results of a periodic targeted impact task extended from Bernstein's seminal work: hammering a nail into wood. We recorded impact forces and kinematics from 10 subjects hammering at different frequencies and with hammers with different physical properties (mass and face area). We found few statistical differences in most measures between different types of hammer, demonstrating human robustness to minor changes in dynamics. Because human motor control is thought to obey optimality principles, we also developed a feedforward optimal simulation with a neuromechanically inspired cost function that reproduces the experimental data. However, Fitts' Law, which relates movement time to distance traveled and target size, did not match our experimental data. We therefore propose a new model in which the distance moved is a logarithmic function of the time to move that yields better results (R2 ≥ 0.99 compared to R2 ≥ 0.88). These results support the argument that humans control movement in an optimal way, but suggest that Fitts' Law may not generalize to periodic impact tasks.
Introduction
Movement is essential to human wellbeing. However, the control of movement is a very difficult problem. To produce deft and robust movements, the human nervous system must continuously control over 600 muscles while handling nonlinearities, nonstationarities, delays, noise, and uncertainties (Franklin and Wolpert, ). Despite these difficulties, humans move with apparent ease. However, human motor capability may become impaired due to age, illness, or injury. Robotic systems are also faced with many of the same challenges (Egeland et al., ; Park, ; Guigon et al., ; Peters et al., ), but meet with much less success than their healthy human counterparts (Yang et al., ; Vanderborght et al., ). A better understanding of the roles that the nervous and musculoskeletal systems play in producing movement will likely lead to advances in rehabilitation and robotic control.
Many neuromechanists employ simple tasks to study the nervous system in action under controlled conditions. Isometric tasks in which subjects interact with an immoble force sensor and reaching tasks in which the hand is moved from one point to another are commonly used to study sensorimotor learning (Rotella et al., ), movement control (Fitts, ), and neurophysiology (Shadmehr and Krakauer, ) in the upper extremity. Subjects may also be asked to interact with robotic co-workers that can record reaching dynamics (Burdet et al., ), generate disturbances, or create force fields (Shadmehr and Mussa-Ivaldi, ) during these tasks. When carefully considered, these experiments can provide a wealth of information on how the nervous system controls movement. However, these tasks are greatly simplified from real-world tasks. To study more complex tasks, some researchers have developed simple games, such as conkers, to study sensoriomotor learning (Sternad et al., ). However, even these studies simplify real-world tasks and rarely consider certain features of real-world tasks such as impacts.
Despite many possible ways to perform most tasks (Bernstein, ), upper extremity movements are highly stereotyped. Researchers note consistent characteristics such as bell-shaped velocity curves (Hollerbach and Atkeson, ; Berardelli et al., ) and speed-accuracy tradeoffs characterized by Fitts' Law (Fitts, ; Bootsma et al., ; Zhai et al., ). Fitts' Law expresses the time to complete a reach as a logarithmic function of the size of the target and the distance to the target (see Equation 1). In experiments relating to Fitts' Law, the kinematics (the beginning and final position of the arm or cursor) are prescribed and the subject is left to determine the time to reach. In certain periodic movements however, the time to complete an upper extremity movement can be specified and the subject left to determine the kinematics.
Movement is constantly refined by biological processes such as learning and evolution (Todorov, ). Because of this constant refinement, many researchers note that optimal control models utilizing cost functions such as minimum variance (Harris and Wolpert, ), minimum effort (Crowninshield and Brand, ), minimum jerk (Flash and Hogan, ), and minimum torque change (Uno et al., ) can be excellent models for the nervous system. In fact, many of the observed stereotypical behaviors discussed in the previous paragraph can be explained by optimality principles. Optimal control models have been used to reproduce human-like behaviors such as reaches (Todorov and Li, ), walking (Anderson and Pandy, ), and jumps (Anderson and Pandy, ; Ong et al., ). Though occasionally studied (Côté et al., ; Müller and Sternad, ), one activity that remains conspicuously unmodeled is Bernstein's hammering task (Bernstein, ; Müller and Sternad, ) that inspired much research into motor control and learning.
Here we extend Bernstein's hammering task into a targeted periodic impact task. We recorded impact forces and upper extremity kinematics in hammering. In order to examine how hammering strategies might change with different conditions, we used a set of hammers with different physical properties (hammer face area and mass) and prescribe different hammering frequencies. We hypothesized that hammering impact velocity and maximal height attained are the result of a tradeoff between maximizing task performance (quantified here as a maximal impact velocity) and minimizing effort (Crowninshield and Brand, ; Nelson, ). In order to test whether the mechanics of this task adhere to current theories in optimal human motor control, we implemented a feedforward optimal controller (Todorov, ) on a planar torque-driven 3-segment dynamical model of the upper extremity holding a hammer (Figure 8) using model parameters from Winter (). Our results show that humans appear to select optimal impact velocities that reflect a tradeoff between accomplishing the task and minimizing effort that do not adhere to Fitts' Law.
Methods
Subjects
Ten healthy male volunteers (age = 27.6 ± 3.6 years, height = 176.9 ± 5 cm, weight = 77.7 ± 11.2 kg) participated in the study. All subjects were right-handed and had no known neuromotor or sensory disorders (self-reported). Prior to their participation, subjects were informed of the course of study and gave their written informed consent in accordance with the code for ethical conduct in research at the Swiss Federal Institute of Technology (EPFL). This study was approved by the EPFL Human Research Ethics Committee (HREC No.: 008-2015/17.08.2015).
Experimental protocol
Each subject was asked to step in front of a table on top of which was a wooden board mounted on a force plate (Kistler Instrument AG, Winterhur, Switzerland) as shown in Figure 1. Subjects were given one of four differently sized and weighted hammers (Table 1) and asked to drive a pre-started nail, i.e., a nail that had previously been driven to the point at which it would stand on its own, into the wooden board while matching their hammer strikes to the clicks of a metronome. Please note that subjects were not explicitly instructed to strike with their maximum impact speed but were allowed to self-select the best impact speed for their skill level. The metronome was set to one of five frequencies: 1, 2, 3, 4, or 5 Hz. The hammer used and metronome frequency for each trial were randomized. Note that subjects were not allowed to do a training trial first, but we assume that the random trial order cancels any learning effects. Subjects were allowed to use their nondominant hands to stabilize the wooden board. In each trial, the forces on the wooden board and the kinematic motion (14 Prime Series cameras, OptiTrac, USA) of the upper extremity and hammer were recorded at 1 kHz and 250 Hz, respectively. After completing the experimental trials, subjects were asked to subjectively rank each of the hammers in order from most to least preferred.
Figure 1
Table 1
| Hammer | Face size [cm × cm] | Weight [kg] |
|---|---|---|
| Small heavy | 1.4 × 1.4 | 0.402 |
| Small light | 1.4 × 1.4 | 0.218 |
| Big heavy | 2.2 × 2.2 | 0.394 |
| Big light | 2.2 × 2.2 | 0.217 |
Specifications of the hammers used in these experiments.
Data processing
Statistical analyses were performed using the Statistics and Machine Learning Toolbox in Matlab. We calculated the hammer velocity, average maximal heights of the hammer, average times required for the hammer to go from maximal height to impact, and average maximal impact forces during hammering for each subject. We then used these average values from each subject for statistical analyses. We investigated the effects of time to impact, maximum height of the movement, and maximal force normalized with the hammer weight using two-way repeated-measures ANOVA with independent variables [hammers(4)×(frequency(5)]. The effect of maximum height of the movement, and maximal force normalized with the hammer weight for each combination of hammers and frequency was further determined using one-way repeated measures ANOVA. The differences between maximal heights and the differences between the normalized maximum forces at impact were tested with post-hoc t-tests with Bonferroni correction. The level of statistical significance used was 0.05 for all statistical tests.
Modeling
In order to determine whether human hammering strategies adhere to Fitts' Law (Fitts, ; Bootsma et al., ; Zhai et al., ), we attempted to fit Fitts' model,
with data collected in our experiment. In this formulation, the movement time, Tf, is a function of the distance from the hammer at peak height to the nail, D, and the face width of the hammer, W. The values of a and b were selected using a least squares difference regression.
In order to examine whether the human nervous system uses optimality principles to control hammering movements, we employed a feedforward optimal controller on two joints (shoulder and elbow) while the wrist was maintained at a desired position with an impedance controller. The human arm holding a hammer was modeled as a 3 link torque-driven robot operating in the saggital plane (Figure 8, right-hand column) whose parameters were computed based on data from Winter () (see Appendix for more details) and whose dynamics are given by
where τ is a vector of joint torques, q, , and are vectors describing the joint angular position, velocity, and acceleration respectively, H(q) is the inertia matrix, consists of the Coriolis, centrifugal, and viscous friction force vectors, g is the gravity force vector, Fe is a vector representing external forces (zero throughout the simulation), and JT is the transpose of the Jacobian matrix. The model was simulated in Matlab using a time step of 0.001 s beginning at the instant after one impact and terminating at the time of the next impact.
Human hammering is a difficult control task due to the need to balance energy transfer to the nail with accuracy. We hypothesize that the human nervous system determines an optimal tradeoff between maximal impact velocity (complete the task in the most effective manner) and minimal effort (Crowninshield and Brand, ; Nelson, ; Missenard and Fernandez, ). We thus determine the optimal joint torques by minimizing the cost function,
where τi,j represents joint torques for i = 1, …, n joints over j = 1, …, T discretized time points, yT and yT−1 are the vertical positions of the hammer head at the last and second-to-last time points of the simulation, 0 ≤ α ≤ 1 was designed as an expertise factor to represent the tradeoff in relative emphasis between impact velocity and effort (large α places more emphasis on energy transfer to the nail and a small α places more emphasis on effort conservation), Cτmax is a scaling factor representing maximal effort (i.e., if maximal torque is applied for the duration of the simulation), and Cẏmax is a scaling factor representing the maximum achievable impact velocity. We compute Cτmax as the discrete integral of the joint torque limits (whichever direction has the larger magnitude) over the length of the simulation and Cẏmax by simulating a hammer trajectory in which α = 1 and Cẏmax = 1. Because maximum effort and final velocity depend on the length of the simulation, we computed unique values of Cẏmax and Cτmax for each hammering frequency. We constrain the model so that the hammer hits the same place in subsequent impacts ((x0, y0) = (xT, yT)) and there is no initial velocity ((ẋ0, ẏ0) = (0, 0)). We match the initial posture (location of (x0, y0) relative to the simulated shoulder) to the average posture used by our subjects determined by inverse kinematics. The terms, Cτmax and Cẏmax, are scaling factors included to facilitate direct comparison of the two terms making up the cost function, minimum effort and maximum final impact velocity. In order to determine whether the parameter, α, is constant within or across individuals, contours of constant α were generated and compared with experimental results. The optimal joint torques were determined using the interior point method implemented with the Matlab Optimization Toolbox.
Results
Experimental results
Subjects were adept at matching hammering frequency with most of those dictated by the metronome. The hammering frequencies achieved by the subjects for metronome frequencies 1, 2, 3, 4, and 5 Hz were 0.99 ± 0.01, 2.02 ± 0.01, 3.01 ± 0.03, 4.01 ± 0.01, and 4.71 ± 0.04 Hz respectively (mean ± standard error). Hammering frequencies of 5 Hz were too fast for our subjects to reliably match. A hammering frequency of 1 Hz was uncomfortably slow for most subjects. To compensate, many subjects developed a strategy of pausing after each impact before initating an up-and-down hammering motion at a more comfortable frequency (Figures 2, 3).
Figure 2
Figure 3
Vertical trajectories (Figure 2) and speeds (Figure 3) exhibited by the subjects in hammering showed very few differences between the different hammers. However, decreasing the hammering frequency increased the variability in these movements. Rather than the single bell-shaped speed profile characteristic of reaching movements, subjects showed a bell-shaped speed profile for raising the hammer and another truncated bell-shaped speed profile for the descending motion (Figure 3).
Analysis of variance showed significant effects of both hammers [F(1.61, 14.5) = 4.95, p = 0.03] and frequencies [F(1.14, 20.25) = 22.35, p < 0.01] on the time to impact from maximum height. There was no significant interaction [F(1.73, 16.05) = 2.56, p = 0.11] between the effects of hammers and frequencies on the time to impact from maximum height. The diagram in Figure 4 shows the means and standard errors (SEM) of time to impact for all hammers and frequencies.
Figure 4
Analysis of variance showed significant effects of hammers and frequencies on the normalized maximal heights. Significant effects of both hammers [F(2.91, 26.2) = 22.8, p < 0.01], frequencies [F(1.77, 15.99) = 53.09, p < 0.01] and interaction between hammers and frequencies [F(4.14, 37.26) = 2.71, p = 0.04] were observed. Further analysis of the effects of hammers on normalized maximal heights showed significant effects of hammers [F(3, 27) = 6.46−11.08, p < 0.01] in all frequencies. Post-hoc t-tests showed that Big Light and Small Heavy maximal heights were statistically different from any of the others [t(9) = 2.97−5.75, p < 0.01]. The diagram in Figure 5 shows the means and standard errors (SEM) of normalized maximal heights for all hammers and frequencies.
Figure 5
Similarly, analysis of variance showed significant effects of hammers and frequencies on the impact forces normalized by hammer mass. Significant effects of both hammers [F(2.06, 23.42) = 32.07, p < 0.01], frequencies [F(1.34, 12.07) = 6.84], but no significant effect of interaction between hammers and frequencies [F(4.39, 39.55) = 0.69, p = 0.63] were observed. Further analysis of the effects of hammers on the impact forces normalized by hammer mass showed significant effects of hammers [F(3, 27) = 11.8−17.97, p < 0.01] in all frequencies. Post-hoc t-tests showed that Small Heavy was statistically different than Small Light [t(9) = 4.71−5.65, p < 0.01] and Big Heavy [t(9) = 5.4−6.61, p < 0.01], and Big Heavy was statistically different than Big Light [t(9) = 3.04−3.91, p < 0.01] for all frequencies. The diagram in Figure 6 shows the means and standard errors (SEM) of impact forces normalized by hammer mass for all hammers and frequencies. The number of impacts needed to totally drive in the nail under each condition–a function of impact velocity–are reported in the appendix (Table A2).
Figure 6
Modeling results
Fitts' Law accurately predicted the movement time from the maximum height to impact (R2 ≥ 0.88). However, despite the high value of R2, the accepted formulation for Fitts' Law does not appear to follow the contours of the experimental data (Figure 7, light gray traces). Therefore, we propose a slightly altered model that reverses the relationship between movement time and distance to move and was able to improve upon Fitts' predictions (R2 ≥ 0.99),
where D is the maximal height of the hammer, Tf is the time from the maximal height to impact in milliseconds, W is the minimum width of the hammer face (our hammers were square, so both face length and width were the same), and a and b are parameters fit to the data using a least squares difference regression (Table 2).
Figure 7
Table 2
| New model | a | b | SSE | RMSE | R2 |
|---|---|---|---|---|---|
| Small heavy | −135 | 30.6 | 0.0005 | 0.015 | 0.99 |
| Small light | −150 | 34.5 | 0.0001 | 0.007 | 0.99 |
| Big heavy | −119 | 25.4 | 0.0012 | 0.020 | 0.99 |
| Big light | −103 | 24.4 | 0.0008 | 0.016 | 0.99 |
| Fitts' model | a | b | SSE | RMSE | R2 |
| Small heavy | 223 | 100 | 776 | 16 | 0.93 |
| Small light | 201 | 99 | 871 | 17 | 0.92 |
| Big heavy | 226 | 118 | 1,350 | 21 | 0.93 |
| Big light | 197 | 117 | 1,841 | 24 | 0.88 |
Parameter estimations for Equation (4) and the original Fitts' Law Equation (1).
The optimal feedforward model was able to accurately reproduce the motions of the arm during hammering (Figures 2, 3, dashed lines, RMSE ≤ 0.1) using the cost function given by Equation (3). This model allows for the generation of optimal hammering trajectories by selecting just one parameter, α. This model also shows that subjects use roughly the same value of α for each hammer, despite the different properties of the different hammers (Figures 2, 3, α values in each row are very similar). The superpositioning of experimental data with computed contours of constant α values (Figure 8) showed that in practice subjects do not use a constant value of α for all hammering frequencies, but rather emphasize lower effort at slower hammering frequencies and energy transfer to the nail at faster hammering frequencies.
Figure 8
The Big Heavy hammer was the most preferred hammer followed by the Big Light, Small Heavy, and Small Light hammers in that order based on subject ratings (Table 3).
Table 3
| Ranking vs. Hammer | S-H | S-L | B-H | B-L |
|---|---|---|---|---|
| Best (10) | 0 | 0 | 9 | 1 |
| 7 | 3 | 1 | 1 | 5 |
| 3 | 3 | 4 | 0 | 3 |
| Worst (1) | 4 | 5 | 0 | 1 |
| Score | 3.4 | 2.4 | 9.7 | 5.5 |
Results of subjects' ranking of the hammers, e.g., S-H, Small Heavy; S-L, Small Light; B-H, Big Heavy; and B-L, Big Light.
Discussion
The goal of this study was to examine the mechanics of a human upper extremity impact task and determine whether existing models of upper limb movement can explain the data. We found that subjects plan optimal trajectories that are a tradeoff between maximum impact velocity and minimal effort reminiscent of Fitts' Law and that are robust to different hammer conditions. However, we found that an altered version of Fitts' Law was able to better match the data than the typical formulation. We also found that end-effector speeds follow a “bell curve and a half” trajectory in hammering in which the hammer head moves upwards with a bell-shaped speed profile and then downwards with a bell-shaped profile before being truncated before the zenith of the curve (Figure 3).
Our analyses showed that Fitts' Law can be applied to human hammering (R2 ≥ 0.88). However, the large R2 values belie an apparent discrepancy between the curves generated using Fitts' Law and the experimental results (Figure 7, gray lines). Therefore, we identified a relationship between movement time and target distance that better reproduces the experimental data (Equation 4, R2 ≥ 0.99). In most Fitts' Law experiments, subjects are prescribed a reaching distance and are asked to move as fast as possible (Fitts, ). However, in our experiments, we constrain permitted movement time using the metronome and subjects were allowed to select the distance to reach. This difference may account for the relative effectiveness of our inverted formulation of Fitts' Law. However, other previous studies have reported violations of Fitts' Law (Adam et al., ; Glazebrook et al., ). Glazebrook et al. () determined that these Fitts' Law violations are the result of pre-planning of movements. This explanation is also certainly plausible in the context of a cyclical task such as hammering. Finally, several studies have noted that Fitts' Law does not hold for movements in which subjects were not asked to move as quickly and as accurately as possible (Young et al., ). We do not explicitly instruct our subjects to move as quickly and accurately as possible. Instead, we instructed them to accomplish a task that is directly dependent on the speed of the movement and allow them to balance that movement speed with their motor capability, which we believe to be an approximation of the instructions to move as quickly and as accurately as possible. In hammering frequencies above 1 Hz, the computed values of α indicate that subjects weight movement speed very highly (Figures 2, 3), and thus likely approach a fast-as-possible movement for which Fitts' Law is presumed to be valid.
Our feedforward optimal hammering simulation was able to reproduce many of the features of human hammering (Figures 2, 3). Our simulations also allow us to show that humans prefer to emphasize energy transfer to the nail (larger values of α) when task constraints are high (high hammering frequencies) and minimal efforts (smaller values of α) when task constraints are low (low hammering frequencies; Figure 8). Our cost function was formulated to minimize the sum squared actuator effort, which serves to keep commanded joint torques small. These small actuation signals prevent excessive energy expenditure during the task (Crowninshield and Brand, ; Missenard and Fernandez, ), but this quadratic formulation might also serve to keep disturbances from motor noise whose effects are multiplicative with actuator effort small (Harris and Wolpert, ; Todorov and Li, ; Franklin and Wolpert, ). In this context, the adaptive prioritization that we observed (changing values of α) might be due to fewer task constraints permitting higher peak heights to be attained at slow hammering frequencies, thus increasing the potential for errors to accrue and increasing the relative importance of accuracy. While the exact cost function used by the nervous system cannot be known exactly, the current formulation reproduces many of the features observed in the experimental results including maximum heights attained, the general trajectories followed, and the robustness to different hammers (similar values of α for different hammers at the same hammering frequencies). However, this model failed to capture the latency after impact before initiating the upward movement of the hammer. This discrepancy may be due to compliance in the musculoskeletal system (e.g., series-elastic muscle-tendon units, Hill, ; Fung, ) that was not captured by our model.
Despite different hammer dynamics (Table 1), hammering kinematics were fairly uniform across many different cases with few statistical differences found between the different hammers in the time from maximal height to impact, maximal hammer height, and impact velocity. Previous studies have suggested that the redundancy of the human musculoskeletal system (Bernstein, ) may contribute to considerable robustness to slight changes in dynamics (Martelli et al., ; Simpson et al., ) or to dysfunction (Arnold et al., ; Hicks et al., ; Correa et al., ; Steele et al., ). While these studies rely on highly redundant lower body musculoskeletal models, other studies examining less redundant body parts have shown limited ability to compensate for dysfunction (Valero-Cuevas and Hentz, ; Kutch and Valero-Cuevas, ). However, detailed models of the upper extremity indicate muscular redundancy on the same level as detailed models of the lower body (Table 4) suggesting that similar robustness to perturbations might be expected. The human nervous system may also select control strategies that are purposefully robust (Mitrovic et al., ; Franklin and Wolpert, ), but our formulation does not include any such criteria, suggesting that consistent movement patterns across conditions might be due to embodied intelligence (e.g., redundant actuators and compliance).
Table 4
| Body part | Degrees of freedom | Number of muscles | References |
|---|---|---|---|
| Upper extremity | 15 | 50 | Holzbaur et al., |
| Lower body | 23 | 54+ | Delp et al., , |
| Index finger | 4 | 7 | Kutch and Valero-Cuevas, |
| Simple leg | 3 | 14 | Kutch and Valero-Cuevas, |
Musculoskeletal models used in examinations of robustness.
Stiffness, or impedance, is a crucial parameter modulated by humans to stably interact with their environment (Burdet et al., ; Franklin and Wolpert, ). Impedance is difficult to record experimentally, but previous studies have attempted to estimate joint stiffnesses based on muscle properties (Hu et al., ), through simulation (Thelen et al., ), or by experimentally recording endpoint stiffnesses (Burdet et al., , ). Because of practical limitations, measurements of muscle activity or impedance were not included in this study, but likely play an important role in impact tasks and should be considered in future works.
Despite the difficulty of controlling a highly nonlinear plant using noisy control signals and noisy sensors with variable delays in an uncertain environment, biological movement appears to be highly robust. However, robustness has not been well addressed in robot learning (Schaal and Atkeson, ; Nguyen-Tuong and Peters, ) primarily because it is difficult to design controllers that are robust to the model structure or parameter errors. One possible solution is to use control policies with optimization criteria based on biological models. For example, the tradeoff between maximizing task performance and accuracy could potentially serve as an optimization criteria for robot hammering.
In this paper, we have extracted the mechanics involved in a targeted upper extremity impact task and demonstrated that the human motor control strategies involved are robust to many different conditions including hammer mass, hammer face area, and timing constraints. We have shown that while many traditional models of human reaching hold for this novel task (bell-shaped speed profiles and Fitts' Law), an altered version of Fitts' Law can better match experimental results. We have also demonstrated that optimality principles previously demonstrated for reaching movements can be generalized to targeted impact tasks and thus lay a framework that can be used for the planning of targeted impact tasks in robots.
Statements
Author contributions
TP, CS, AU, and AI contributed to the design, execution and drafting of this work, and approved the final manuscript. Experimental data was collected and analyzed by TP and CS.
Funding
This material is based upon work supported by Sciex-NMSCH project 14.069 to TP and a Fulbright Scholarship, Swiss Government Excellence Scholarship, and National Science Foundation Graduate Research Fellowship to CSS. The research leading to these results has received funding from the Slovenian Research Agency under grant agreement no. J2-7360.
Acknowledgments
The authors thank the members of the Biorobotics laboratory at EPFL for their assistance with data collection and participant recruitment.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AdamJ. J.MolR.PrattJ.FischerM. H. (2006). Moving farther but faster: an exception to fitts's law. Psychol. Sci.17, 794–798. 10.1111/j.1467-9280.2006.01784.x
2
AndersonF. C.PandyM. G. (1999). A dynamic optimization solution for vertical jumping in three dimensions. Comput. Methods Biomech. Biomed. Eng.2, 201–231. 10.1080/10255849908907988
3
AndersonF. C.PandyM. G. (2001). Dynamic optimization of human walking. J. Biomech. Eng.123:381. 10.1115/1.1392310
4
ArnoldA. S.AndersonF. C.PandyM. G.DelpS. L. (2005). Muscular contributions to hip and knee extension during the single limb stance phase of normal gait: a framework for investigating the causes of crouch gait. J. Biomech.38, 2181–2189. 10.1016/j.jbiomech.2004.09.036
5
BerardelliA.HallettM.RothwellJ. C.AgostinoR.ManfrediM.ThompsonP. D.et al. (1996). Single-joint rapid arm movements in normal subjects and in patients with motor disorders. Brain119(Pt 2), 661–674. 10.1093/brain/119.2.661
6
BernsteinN. (1967). The Co-ordination and Regulation of Movements. Oxford; New York, NY: Pergamon Press.
7
BootsmaR. J.FernandezL.MottetD. (2004). Behind fitts' law: kinematic patterns in goal-directed movements. Int. J. Hum. Comput. Stud.61, 811–821. 10.1016/j.ijhcs.2004.09.004
8
BurdetE.OsuR.FranklinD.YoshiokaT.MilnerT.KawatoM. (2000). A method for measuring endpoint stiffness during multi-joint arm movements. J. Biomech.33, 1705–1709. 10.1016/S0021-9290(00)00142-1
9
BurdetE.OsuR.FranklinD. W.MilnerT. E.KawatoM. (2001). The central nervous system stabilizes unstable dynamics by learning optimal impedance. Nature414, 446–449. 10.1038/35106566
10
CorreaT. A.SchacheA. G.GrahamH. K.BakerR.ThomasonP.PandyM. G. (2012). Potential of lower-limb muscles to accelerate the body during cerebral palsy gait. Gait Posture36, 194–200. 10.1016/j.gaitpost.2012.02.014
11
CôtéJ. N.FeldmanA. G.MathieuP. A.LevinM. F. (2008). Effects of fatigue on intermuscular coordination during repetitive hammering. Motor Control12, 79–92. 10.1123/mcj.12.2.79
12
CrowninshieldR. D.BrandR. A. (1981). A physiologically based criterion of muscle force prediction in locomotion. J. Biomech.14, 793–801. 10.1016/0021-9290(81)90035-X
13
DelpS. L.AndersonF. C.ArnoldA. S.LoanP.HabibA.JohnC. T.et al. (2007). OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Trans. Biomed. Eng.54, 1940–1950. 10.1109/TBME.2007.901024
14
DelpS. L.LoanP.HoyM. G.ZajacF. E.ToppE. L.RosenJ. M. (1990). An interactiv gaphics-based model of the lower extremity to study orthopaedic surgical procedures. IEEE Trans. Biomed. Eng.37, 757–767. 10.1109/10.102791
15
EgelandO.SagliJ.SpangeloI.ChiaveriniS. (1991). A damped least-squares solution to redundancy resolution, in Proceedings 1991 IEEE International Conference on Robotics and Automation, (Sacramento, CA: IEEE Computer Society Press), 945–950.
16
FittsP. M. (1954). The information capacity of the human motor. J. Exp. Biol.47, 381–391.
17
FlashT.HoganN. (1985). The coordination of arm movements: an experimentally confirmed mathematical model. J. Neurosci.5, 1688–1703.
18
FranklinD. W.WolpertD. M. (2011). Computational mechanisms of sensorimotor control. Neuron72, 425–442. 10.1016/j.neuron.2011.10.006
19
FungY.-C. (2013). Biomechanics: Mechanical Properties of Living Tissues. New York, NY: Springer-Verlag.
20
GlazebrookC. M.KiernanD.WelshT. N.TremblayL. (2015). How one breaks fittss law and gets away with it: moving further and faster involves more efficient online control. Hum. Mov. Sci.39, 163–176. 10.1016/j.humov.2014.11.005
21
GuigonE.BaraducP.DesmurgetM. (2007). Computational motor control: redundancy and invariance. J. Neurophysiol.97, 331–347. 10.1152/jn.00290.2006
22
HarrisC. M.WolpertD. M. (1998). Signal-dependent noise determinesmotorplanning. Nature394, 780–784. 10.1038/29528
23
HicksJ. L.SchwartzM. H.ArnoldA. S.DelpS. L. (2008). Crouched postures reduce the capacity of muscles to extend the hip and knee during the single-limb stance phase of gait. J. Biomech.41, 960–967. 10.1016/j.jbiomech.2008.01.002
24
HillA. V. (1938). The heat of shortening and the dynamic constants of muscle. Proc. R. Soc. Lond. B Biol. Sci.126, 136–195.
25
HollerbachJ. M.AtkesonC. G. (1987). Deducing planning variables from experimental arm trajectories: pitfalls and possibilities. Biol. Cybernet.56, 279–292. 10.1007/BF00319509
26
HolzbaurK. R. S.MurrayW. M.DelpS. L. (2005). A model of the upper extremity for simulating musculoskeletal surgery and analyzing neuromuscular control. Ann. Biomed. Eng.33, 829–840. 10.1007/s10439-005-3320-7
27
HuX.MurrayW. M.PerreaultE. J. (2011). Muscle short-range stiffness can be used to estimate the endpoint stiffness of the human arm. J. Neurophysiol.105, 1633–1641. 10.1152/jn.00537.2010
28
KutchJ. J.Valero-CuevasF. J. (2011). Muscle redundancy does not imply robustness to muscle dysfunction. J. Biomech.44, 1264–1270. 10.1016/j.jbiomech.2011.02.014
29
MartelliS.CalvettiD.SomersaloE.VicecontiM. (2015). Stochastic modelling of muscle recruitment during activity. Interface Focus5:20140094. 10.1098/rsfs.2014.0094
30
MissenardO.FernandezL. (2011). Moving faster while preserving accuracy. Neuroscience197, 233–241. 10.1016/j.neuroscience.2011.09.020
31
MitrovicD.KlankeS.OsuR.KawatoM.VijayakumarS. (2010). A computational model of limb impedance control based on principles of internal model uncertainty. PLoS ONE5:e13601. 10.1371/journal.pone.0013601
32
MüllerH.SternadD. (2009). Motor learning: changes in the structure of variability in a redundant task, in Progress in Motor Control, ed DagmarS. (Boston, MA: Springer US), 439–456.
33
NelsonW. L. (1983). Physical principles for economies of skilled movements. Biol. Cybernet.46, 135–147. 10.1007/BF00339982
34
Nguyen-TuongD.PetersJ. (2011). Model learning for robot control: a survey. Cogn. Process.12, 319–340. 10.1007/s10339-011-0404-1
35
OngC. F.HicksJ. L.DelpS. L. (2016). Simulation-based design for wearable robotic systems: an optimization framework for enhancing a standing long jump. IEEE Trans. Biomed. Eng.63, 894–903. 10.1109/TBME.2015.2463077
36
ParkJ. (2002). Characterization of instability of dynamic control for kinematically redundant manipulators. Proc. IEEE Int. Confe. Robot. Autom. 3, 2400–2405. 10.1109/ROBOT.2002.1013591
37
PetersJ.MüllingK.KoberJ.Nguyen-TuongD.KrömerO. (2009). Towards motor skill learning for robotics, in Proceedings of the International Symposium on Robotics Research ISRR, (Berlin; Heidelberg: Springer), 1–14.
38
RotellaM. F.NiskyI.KoehlerM.RinderknechtM. D.BastianA. J.OkamuraA. M. (2015). Learning and generalization in an isometric visuomotor task. J. Neurophysiol.113, 1873–1884. 10.1152/jn.00255.2014
39
SchaalS.AtkesonC. (2010). Learning control in robotics. IEEE Robot. Autom. Magaz.17, 20–29. 10.1109/MRA.2010.936957
40
ShadmehrR.KrakauerJ. W. (2008). A computational neuroanatomy for motor control. Exp. Brain Res.185, 359–381. 10.1007/s00221-008-1280-5
41
ShadmehrR.Mussa-IvaldiF. A. (1994). Adaptive representation of dynamics during learning of a motor task. J. Neurosci.14, 3208–3224.
42
SimpsonC. S.SohnM. H.AllenJ. L.TingL. H. (2015). Feasible muscle activation ranges based on inverse dynamics analyses of human walking. J. Biomech.48, 2990–2997. 10.1016/j.jbiomech.2015.07.037
43
SteeleK. M.van der KrogtM. M.SchwartzM. H.DelpS. L. (2012). How much muscle strength is required to walk in a crouch gait?J. Biomech.45, 2564–2569. 10.1016/j.jbiomech.2012.07.028
44
SternadD.AbeM. O.HuX.MllerH. (2011). Neuromotor noise, error tolerance and velocity-dependent costs in skilled performance. PLoS Comput. Biol.7:e1002159. 10.1371/journal.pcbi.1002159
45
ThelenD. G.AndersonF. C.DelpS. L. (2003). Generating dynamic simulations of movement using computed muscle control. J. Biomech.36, 321–328. 10.1016/S0021-9290(02)00432-3
46
TodorovE. (2004). Optimality principles in sensorimotor control. Nat. Neurosci.7, 907–915. 10.1038/nn1309
47
TodorovE.LiW. (2005). A generalized iterative lqg method for locally-optimal feedback control of constrained nonlinear stochastic systems in Proceedings of the 2005, American Control Conference, 2005 (Portland, OR: IEEE), 300–306.
48
UnoY.KawatoM.SuzukiR. (1989). Formation and control of optimal trajectory in human multijoint arm movement. Biol. Cybernet.61, 89–101. 10.1007/BF00204593
49
Valero-CuevasF. J.HentzV. R. (2002). Releasing the A3 pulley and leaving flexor superficialis intact increases pinch force following the Zancolli lasso procedures to prevent claw deformity in the intrinsic palsied finger. J. Orthopaed. Res.20, 902–909. 10.1016/S0736-0266(02)00040-2
50
VanderborghtB.Albu-SchaefferA.BicchiA.BurdetE.CaldwellD. G.CarloniR.et al. (2013). Variable impedance actuators: a review. Robot. Autonom. Syst.61, 1601–1614. 10.1016/j.robot.2013.06.009
51
WinterD. A. (2009). Biomechanics and Motor Control of Human Movement, Vol. 2. Hoboken, NJ: John Wiley & Sons, Inc.
52
YangC.GaneshG.HaddadinS.ParuselS.Albu-SchäefferA.BurdetE. (2011). Human-like adaptation of force and impedance in stable and unstable interactions. IEEE Trans. Robot.27, 918–930. 10.1109/TRO.2011.2158251
53
YoungS. J.PrattJ.ChauT. (2009). Target-directed movements at a comfortable pace: movement duration and fitts's law. J. Motor Behav.41, 339–346. 10.3200/JMBR.41.4.339-346
54
ZhaiS.KongJ.RenX. (2004). Speed–accuracy tradeoff in Fitts' law tasks—on the equivalency of actual and nominal pointing precision. Int. J. Hum. Comput. Stud.61, 823–856. 10.1016/j.ijhcs.2004.09.007
Appendix
Model parameters
In this paper, we have modeled the human arm holding a hammer as a torque driven 3 degree of freedom (DOF) robot operating in the sagittal plane. Each DOF is driven by an independently controlled torque generator capable of producing both positive and negative torques. The robot parameters (link lengths, center of mass locations, etc.) were computed based on data from Winter () using the mean height (176.9 cm) and weight (77 kg) of subjects that participated in this study. The hammers were simulated by adding the relevant mass to the end effector. Parameters of the hammers are in Table 1. The robot parameters are given in Table A1.
Table A1
| Part | Upper Arm | Lower Arm | Hand + Hammer |
|---|---|---|---|
| Link No. | 1 | 2 | 3 |
| Length | 0.2794 | 0.2667 | 0.0855 |
| Mass | 2.0751 | 1.2225 | 0.4810 |
| Center of mass location | 0.1613 | 0.1220 | 0.0676 |
| Inertia | 0.0131 | 0.0067 | 0.0010 |
Parametes of the dynamical model.
Note that all units are SI (m, kg). Note also that the center of mass location is relative to the proximal end of the relevant link.
Hammer hits
The table reports how many cycles were necessary to totally drive in the nail with respect to the hammer and frequency.
Table A2
| Hammer | |||||
|---|---|---|---|---|---|
| S-H | S-L | B-H | B-L | ||
| Frequency | 1 | 12 ± 5.114 | 9.9 ± 1.215 | 7.1 ± 1.663 | 8.5 ± 1.258 |
| [Hz] | 2 | 16 ± 3.303 | 12.2 ± 1.679 | 9.7 ± 2.285 | 13.4 ± 2.817 |
| 3 | 19.7 ± 4.6 | 14.4 ± 2.579 | 10.6 ± 1.655 | 14.4 ± 1.968 | |
| 4 | 22 ± 3.48 | 19.5 ± 2.693 | 12.9 ± 2.089 | 14.1 ± 2.677 | |
| 5 | 29.2 ± 6.2 | 23.8 ± 3.62 | 14.2 ± 2.732 | 20.2 ± 3.511 | |
Number of impacts required to drive nail by hammer and frequency (mean ± standard error).
Summary
Keywords
motor control, biomechanics, upper extremity, optimal control, arm movement, impact, Fitts' Law
Citation
Petrič T, Simpson CS, Ude A and Ijspeert AJ (2017) Hammering Does Not Fit Fitts' Law. Front. Comput. Neurosci. 11:45. doi: 10.3389/fncom.2017.00045
Received
19 October 2016
Accepted
15 May 2017
Published
29 May 2017
Volume
11 - 2017
Edited by
Francisco J. Valero-Cuevas, University of Southern California, United States
Reviewed by
Bastien Berret, University of Paris-Sud, Universit Paris-Saclay, France; J. Michael Herrmann, University of Edinburgh, United Kingdom
Updates
Copyright
© 2017 Petrič, Simpson, Ude and Ijspeert.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tadej Petrič tadej.petric@ijs.si
†These authors have contributed equally to this work.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.