Short- and long-term changes in joint co-contraction associated with motor learning as revealed from surface EMG.

In the field of motor control, two hypotheses have been controversial: whether the brain acquires internal models that generate accurate motor commands, or whether the brain avoids this by using the viscoelasticity of musculoskeletal system. Recent observations on relatively low stiffness during trained movements support the existence of internal models. However, no study has revealed the decrease in viscoelasticity associated with learning that would imply improvement of internal models as well as synergy between the two hypothetical mechanisms. Previously observed decreases in electromyogram (EMG) might have other explanations, such as trajectory modifications that reduce joint torques. To circumvent such complications, we required strict trajectory control and examined only successful trials having identical trajectory and torque profiles. Subjects were asked to perform a hand movement in unison with a target moving along a specified and unusual trajectory, with shoulder and elbow in the horizontal plane at the shoulder level. To evaluate joint viscoelasticity during the learning of this movement, we proposed an index of muscle co-contraction around the joint (IMCJ). The IMCJ was defined as the summation of the absolute values of antagonistic muscle torques around the joint and computed from the linear relation between surface EMG and joint torque. The IMCJ during isometric contraction, as well as during movements, was confirmed to correlate well with joint stiffness estimated using the conventional method, i.e., applying mechanical perturbations. Accordingly, the IMCJ during the learning of the movement was computed for each joint of each trial using estimated EMG-torque relationship. At the same time, the performance error for each trial was specified as the root mean square of the distance between the target and hand at each time step over the entire trajectory. The time-series data of IMCJ and performance error were decomposed into long-term components that showed decreases in IMCJ in accordance with learning with little change in the trajectory and short-term interactions between the IMCJ and performance error. A cross-correlation analysis and impulse responses both suggested that higher IMCJs follow poor performances, and lower IMCJs follow good performances within a few successive trials. Our results support the hypothesis that viscoelasticity contributes more when internal models are inaccurate, while internal models contribute more after the completion of learning. It is demonstrated that the CNS regulates viscoelasticity on a short- and long-term basis depending on performance error and finally acquires smooth and accurate movements while maintaining stability during the entire learning process.

Functional significance of stiffness in adaptation of multijoint arm movements to stable and unstable dynamics.

This study compared the mechanisms of adaptation to stable and unstable dynamics from the perspective of changes in joint mechanics. Subjects were instructed to make point to point movements in force fields generated by a robotic manipulandum which interacted with the arm in either a stable or an unstable manner. After subjects adjusted to the initial disturbing effects of the force fields they were able to produce normal straight movements to the target. In the case of the stable interaction, subjects modified the joint torques in order to appropriately compensate for the force field. No change in joint torque or endpoint force was required or observed in the case of the unstable interaction. After adaptation, the endpoint stiffness of the arm was measured by applying displacements to the hand in eight different directions midway through the movements. This was compared to the stiffness measured similarly during movements in a null force field. After adaptation, the endpoint stiffness under both the stable and unstable dynamics was modified relative to the null field. Adaptation to unstable dynamics was achieved by selective modification of endpoint stiffness in the direction of the instability. To investigate whether the change in endpoint stiffness could be accounted for by change in joint torque or endpoint force, we estimated the change in stiffness on each trial based on the change in joint torque relative to the null field. For stable dynamics the change in endpoint stiffness was accurately predicted. However, for unstable dynamics the change in endpoint stiffness could not be reproduced. In fact, the predicted endpoint stiffness was similar to that in the null force field. Thus, the change in endpoint stiffness seen after adaptation to stable dynamics was directly related to changes in net joint torque necessary to compensate for the dynamics in contrast to adaptation to unstable dynamics, where a selective change in endpoint stiffness occurred without any modification of net joint torque.

Adaptive control of stiffness to stabilize hand position with large loads.

The goal of this work was to investigate stability in relation to the magnitude and direction of forces applied by the hand. The endpoint stiffness and joint stiffness of the arm were measured during a postural task in which subjects exerted up to 30% maximum voluntary force in each of four directions while controlling the position of the hand. All four coefficients of the joint stiffness matrix were found to vary linearly with both elbow and shoulder torque. This contrasts with the results of a previous study, which employed a force control task and concluded that the joint stiffness coefficients varied linearly with either shoulder or elbow torque but not both. Joint stiffness was transformed into endpoint stiffness to compare the effect on stability as endpoint force increased. When the joint stiffness coefficients were modeled as varying with the net torque at only one joint, as in the previous study, we found that hand position became unstable if endpoint force exceeded about 22 N in a specific direction. This did not occur when the joint stiffness coefficients were modeled as varying with the net torque at both joints, as in the present study. Rather, hand position became increasingly more stable as endpoint force increased for all directions of applied force. Our analysis suggests that co-contraction of biarticular muscles was primarily responsible for the increased stability. This clearly demonstrates how the central nervous system can selectively adapt the impedance of the arm in a specific direction to stabilize hand position when the force applied by the hand has a destabilizing effect in that direction.

Adaptation to stable and unstable dynamics achieved by combined impedance control and inverse dynamics model.

This study compared adaptation in novel force fields where trajectories were initially either stable or unstable to elucidate the processes of learning novel skills and adapting to new environments. Subjects learned to move in a null force field (NF), which was unexpectedly changed either to a velocity-dependent force field (VF), which resulted in perturbed but stable hand trajectories, or a position-dependent divergent force field (DF), which resulted in unstable trajectories. With practice, subjects learned to compensate for the perturbations produced by both force fields. Adaptation was characterized by an initial increase in the activation of all muscles followed by a gradual reduction. The time course of the increase in activation was correlated with a reduction in hand-path error for the DF but not for the VF. Adaptation to the VF could have been achieved solely by formation of an inverse dynamics model and adaptation to the DF solely by impedance control. However, indices of learning, such as hand-path error, joint torque, and electromyographic activation and deactivation suggest that the CNS combined these processes during adaptation to both force fields. Our results suggest that during the early phase of learning there is an increase in endpoint stiffness that serves to reduce hand-path error and provides additional stability, regardless of whether the dynamics are stable or unstable. We suggest that the motor control system utilizes an inverse dynamics model to learn the mean dynamics and an impedance controller to assist in the formation of the inverse dynamics model and to generate needed stability.

Impedance control balances stability with metabolically costly muscle activation.

Humans are able to stabilize their movements in environments with unstable dynamics by selectively modifying arm impedance independently of force and torque. We further investigated adaptation to unstable dynamics to determine whether the CNS maintains a constant overall level of stability as the instability of the environmental dynamics is varied. Subjects performed reaching movements in unstable force fields of varying strength, generated by a robotic manipulator. Although the force fields disrupted the initial movements, subjects were able to adapt to the novel dynamics and learned to produce straight trajectories. After adaptation, the endpoint stiffness of the arm was measured at the midpoint of the movement. The stiffness had been selectively modified in the direction of the instability. The stiffness in the stable direction was relatively unchanged from that measured during movements in a null force field prior to exposure to the unstable force field. This impedance modification was achieved without changes in force and torque. The overall stiffness of the arm and environment in the direction of instability was adapted to the force field strength such that it remained equivalent to that of the null force field. This suggests that the CNS attempts both to maintain a minimum level of stability and minimize energy expenditure.

Accurate real-time feedback of surface EMG during fMRI.

Real-time acquisition of EMG during functional MRI (fMRI) provides a novel method of controlling motor experiments in the scanner using feedback of EMG. Because of the redundancy in the human muscle system, this is not possible from recordings of joint torque and kinematics alone, because these provide no information about individual muscle activation. This is particularly critical during brain imaging because brain activations are not only related to joint torques and kinematics but are also related to individual muscle activation. However, EMG collected during imaging is corrupted by large artifacts induced by the varying magnetic fields and radio frequency (RF) pulses in the scanner. Methods proposed in literature for artifact removal are complex, computationally expensive, and difficult to implement for real-time noise removal. We describe an acquisition system and algorithm that enables real-time acquisition for the first time. The algorithm removes particular frequencies from the EMG spectrum in which the noise is concentrated. Although this decreases the power content of the EMG, this method provides excellent estimates of EMG with good resolution. Comparisons show that the cleaned EMG obtained with the algorithm is, like actual EMG, very well correlated with joint torque and can thus be used for real-time visual feedback during functional studies.

Immersed Boundary Method for generalised finite volume and finite difference Navier-Stokes solvers

In Immersed Boundary Methods (IBM) the effect of complex geometries is introduced through the forces added in the Navier-Stokes solver at the grid points in the vicinity of the immersed boundaries. Most of the methods in the literature have been used with Cartesian grids. Moreover many of the methods developed in the literature do not satisfy some basic conservation properties (the conservation of torque, for instance) on non-uniform meshes. In this paper we will follow the RKPM method originated by Liu et al. [1] to build locally regularized functions that verify a number of integral conditions. These local approximants will be used both for interpolating the velocity field and for spreading the singular force field in the framework of a pressure correction scheme for the incompressible Navier-Stokes equations. We will also demonstrate the robustness and effectiveness of the scheme through various examples. Copyright © 2010 by ASME.

Operation of brushless doubly-fed machine for drive applications

This paper presents dynamic and steady-state performance of the Brushless Doubly-Fed Machine (BDFM) operating as a variable speed drive. A simple closed-loop control system is used which only requires a speed feedback. The controller is capable of stabilising the machine when changes in speed and torque are applied. The machine starts in cascade mode and then makes a transition to the synchronous mode to reach the desired speed. This will allow a uni-directional converter to be used. The experiments included in this paper were carried out on a 180 frame size BDFM.

Rain igestion in axial flow compressors at part speed

New experimental work is reported on the effects of water ingestion on the performance of an axial flow compressor. The background to the work is the effect that heavy rain has on an aeroengine compressor when operating in a "descent idle" mode, i.e., when the compressor is operating at part speed and when the aeromechanical effects of water ingestion are more important than the thermodynamic effects. Most of our existing knowledge in this field comes from whole engine tests. The current work provides the first known results from direct measurements on a stand-alone compressor. The influence of droplet size on path trajectory is considered both computationally and experimentally to show that most rain droplets will collide with the first row of rotor blades. The water on the blades is then centrifuged toward the casing where the normal airflow patterns in the vicinity of the rotor tips are disrupted. The result of this disruption is a reduction in compressor delivery pressure and an increase in the torque required to keep the compressor speed constant. Both effects reduce the efficiency of the machine. The behavior of the water in the blade rows is examined in detail, and simple models are proposed to explain the loss of pressure rise and the increase in torque. The measurements were obtained in a low speed compressor, making it possible to study the mechanical (increase in torque) and aerodynamic (reduction in pressure rise) effects of water ingestion without the added complication of thermodynamic effects. Copyright © 2008 by ASME.

A path-following driver-vehicle model with neuromuscular dynamics, including measured and simulated responses to a step in steering angle overlay

An existing driver-vehicle model with neuromuscular dynamics is improved in the areas of cognitive delay, intrinsic muscle dynamics and alpha-gamma co-activation. The model is used to investigate the influence of steering torque feedback and neuromuscular dynamics on the vehicle response to lateral force disturbances. When steering torque feedback is present, it is found that the longitudinal position of the lateral disturbance has a significant influence on whether the drivers reflex response reinforces or attenuates the effect of the disturbance. The response to angle and torque overlay inputs to the steering system is also investigated. The presence of the steering torque feedback reduced the disturbing effect of torque overlay and angle overlay inputs. Reflex action reduced the disturbing effect of a torque overlay input, but increased the disturbing effect of an angle overlay input. Experiments on a driving simulator showed that measured handwheel angle response to an angle overlay input was consistent with the response predicted by the model with reflex action. However, there was significant intra-and inter-subject variability. The results highlight the significance of a drivers neuromuscular dynamics in determining the vehicle response to disturbances. © 2012 Copyright Taylor and Francis Group, LLC.

Optical trapping of carbon nanotubes and graphene

We study optical trapping of nanotubes and graphene. We extract the distribution of both centre-of-mass and angular fuctuations from three-dimensional tracking of these optically trapped carbon nanostructures. The optical force and torque constants are measured from auto and cross-correlation of the tracking signals. We demonstrate that nanotubes enable nanometer spatial, and femto-Newton force resolution in photonic force microscopy by accurately measuring the radiation pressure in a double frequency optical tweezers. Finally, we integrate optical trapping with Raman and photoluminescence spectroscopy demonstrating the use of a Raman and photoluminescence tweezers by investigating the spectroscopy of nanotubes and graphene fakes in solution. Experimental results are compared with calculations based on electromagnetic scattering theory. © 2011 by the Author(s); licensee Accademia Peloritana dei Pericolanti, Messina, Italy.

Implementation of a PI phase angle controller for finite element analysis of the BDFM

Time-stepping finite element analysis of the BDFM for a specific load condition is shown to be a challenging problem because the excitation required cannot be predetermined and the BDFM is not open loops stable for all operating conditions. A simulation approach using feedback control to set the torque and stabilise the BDFM is presented together with implementation details. The performance of the simulation approach is demonstrated with an example and computed results are compared with measurements.

Fault-tolerant model predictive control of a wind turbine benchmark

This paper aims to solve the fault tolerant control problem of a wind turbine benchmark. A hierarchical controller with model predictive pre-compensators, a global model predictive controller and a supervisory controller is proposed. In the model predictive pre-compensator, an extended Kalman Filter is designed to estimate the system states and various fault parameters. Based on the estimation, a group of model predictive controllers are designed to compensate the fault effects for each component of the wind turbine. The global MPC is used to schedule the operation of the components and exploit potential system-level redundancies. Extensive simulations of various fault conditions show that the proposed controller has small transients when faults occur and uses smoother and smaller generator torque and pitch angle inputs than the default controller. This paper shows that MPC can be a good candidate for fault tolerant controllers, especially the one with an adaptive internal model combined with a parameter estimation and update mechanism, such as an extended Kalman Filter. © 2012 IFAC.

Generalized vector control for brushless doubly fed machines with nested-loop rotor

This paper presents a generalized vector control system for a generic brushless doubly fed (induction) machine (BDFM) with nested-loop type rotor. The generic BDFM consists of p1/p2 pole-pair stator windings and a nested-loop rotor with N number of loops per nest. The vector control system is derived based on the basic BDFM equation in the synchronous mode accompanied with an appropriate synchronization approach to the grid. An analysis is performed for the vector control system using the generic BDFM vector model. The analysis proves the efficacy of the proposed approach in BDFM electromagnetic torque and rotor flux control. In fact, in the proposed vector control system, the BDFM torque can be controlled very effectively promising a high-performance BDFM shaft speed control system. A closed-loop shaft speed control system is composed based on the presented vector control system whose performance is examined both in simulations and experiments. The results confirm the high performance of the proposed approach in BDFM shaft speed control as well as a very close agreement between the simulations and experiments. Tests are performed on a 180-frame prototype BDFM. © 2012 IEEE.

Modelling and control of nonlinear systems using Gaussian processes with partial model information

Gaussian processes are gaining increasing popularity among the control community, in particular for the modelling of discrete time state space systems. However, it has not been clear how to incorporate model information, in the form of known state relationships, when using a Gaussian process as a predictive model. An obvious example of known prior information is position and velocity related states. Incorporation of such information would be beneficial both computationally and for faster dynamics learning. This paper introduces a method of achieving this, yielding faster dynamics learning and a reduction in computational effort from O(Dn2) to O((D - F)n2) in the prediction stage for a system with D states, F known state relationships and n observations. The effectiveness of the method is demonstrated through its inclusion in the PILCO learning algorithm with application to the swing-up and balance of a torque-limited pendulum and the balancing of a robotic unicycle in simulation. © 2012 IEEE.