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.

Transducer and base compliance alter the in situ 6 dof force measured from muscle during an isometric contraction in a multi-joint limb.

Although musculoskeletal models are commonly used, validating the muscle actions predicted by such models is often difficult. In situ isometric measurements are a possible solution. The base of the skeleton is immobilized and the endpoint of the limb is rigidly attached to a 6-axis force transducer. Individual muscles are stimulated and the resulting forces and moments recorded. Such analyses generally assume idealized conditions. In this study we have developed an analysis taking into account the compliances due to imperfect fixation of the skeleton, imperfect attachment of the force transducer, and extra degrees of freedom (dof) in the joints that sometimes become necessary in fixed end contractions. We use simulations of the rat hindlimb to illustrate the consequences of such compliances. We show that when the limb is overconstrained, i.e., when there are fewer dof within the limb than are restrained by the skeletal fixation, the compliances of the skeletal fixation and of the transducer attachment can significantly affect measured forces and moments. When the limb dofs and restrained dofs are matched, however, the measured forces and moments are independent of these compliances. We also show that this framework can be used to model limb dofs, so that rather than simply omitting dofs in which a limb does not move (e.g., abduction at the knee), the limited motion of the limb in these dofs can be more realistically modeled as a very low compliance. Finally, we discuss the practical implications of these results to experimental measurements of muscle actions.

Is titin a 'winding filament'? A new twist on muscle contraction.

Recent studies have demonstrated a role for the elastic protein titin in active muscle, but the mechanisms by which titin plays this role remain to be elucidated. In active muscle, Ca(2+)-binding has been shown to increase titin stiffness, but the observed increase is too small to explain the increased stiffness of parallel elastic elements upon muscle activation. We propose a 'winding filament' mechanism for titin's role in active muscle. First, we hypothesize that Ca(2+)-dependent binding of titin's N2A region to thin filaments increases titin stiffness by preventing low-force straightening of proximal immunoglobulin domains that occurs during passive stretch. This mechanism explains the difference in length dependence of force between skeletal myofibrils and cardiac myocytes. Second, we hypothesize that cross-bridges serve not only as motors that pull thin filaments towards the M-line, but also as rotors that wind titin on the thin filaments, storing elastic potential energy in PEVK during force development and active stretch. Energy stored during force development can be recovered during active shortening. The winding filament hypothesis accounts for force enhancement during stretch and force depression during shortening, and provides testable predictions that will encourage new directions for research on mechanisms of muscle contraction.

Contraction of monotone phase-coupled oscillators

This paper establishes a global contraction property for networks of phase-coupled oscillators characterized by a monotone coupling function. The contraction measure is a total variation distance. The contraction property determines the asymptotic behavior of the network, which is either finite-time synchronization or asymptotic convergence to a splay state. © 2012 Elsevier B.V. All rights reserved.

Contraction and observer design on cones

We consider the problem of positive observer design for positive systems defined on solid cones in Banach spaces. The design is based on the Hilbert metric and convergence properties are analyzed in the light of the Birkhoff theorem. Two main applications are discussed: positive observers for systems defined in the positive orthant, and positive observers on the cone of positive semi-definite matrices with a view on quantum systems. © 2011 IEEE.

A differential lyapunov framework for contraction analysis

Lyapunov's second theorem is an essential tool for stability analysis of differential equations. The paper provides an analog theorem for incremental stability analysis by lifting the Lyapunov function to the tangent bundle. The Lyapunov function endows the state-space with a Finsler structure. Incremental stability is inferred from infinitesimal contraction of the Finsler metrics through integration along solutions curves. © 2013 IEEE.

Contraction/expansion flows: The pressure drop and related issues

Binding, David; Phillips, P.M.; Philips, T.N., (2006) 'Contraction/expansion flows: The pressure drop and related issues', Journal of Non-Newtonian Fluid Mechanics 137 pp.31-38 RAE2008

Transformation from natural (thermal contraction) to anthropogenic (resource exploitation) depressions in the Krotoszyn-Koźmin-Raszków area (Polish Lowland)

An exceptional concentration of almost identical depressions exist near the small towns of Krotoszyn, Koźmin and Raszków (southern Wielkopolska). Their origin is, however, different from that of the typical post glacial-relief: they are Man-made enlarged thermal-contraction structures that developed at the very end of the Middle Polish (Warthian) glaciation and during the North Polish (Weichselian) glaciation, most probably under periglacial conditions.

Viscoelastic flow through a planar contraction using a semi-Lagrangian finite volume method

A semi-Lagrangian finite volume scheme for solving viscoelastic flow problems is presented. A staggered grid arrangement is used in which the dependent variables are located at different mesh points in the computational domain. The convection terms in the momentum and constitutive equations are treated using a semi-Lagrangian approach in which particles on a regular grid are traced backwards over a single time-step. The method is applied to the 4 : 1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions. The development of vortex behaviour with increasing values of We is analyzed.

A semi-Lagrangian finite volume method for Newtonian contraction flows

A new finite volume method for solving the incompressible Navier--Stokes equations is presented. The main features of this method are the location of the velocity components and pressure on different staggered grids and a semi-Lagrangian method for the treatment of convection. An interpolation procedure based on area-weighting is used for the convection part of the computation. The method is applied to flow through a constricted channel, and results are obtained for Reynolds numbers, based on half the flow rate, up to 1000. The behavior of the vortex in the salient corner is investigated qualitatively and quantitatively, and excellent agreement is found with the numerical results of Dennis and Smith [Proc. Roy. Soc. London A, 372 (1980), pp. 393-414] and the asymptotic theory of Smith [J. Fluid Mech., 90 (1979), pp. 725-754].