Non-Parametric Nonlinear Parameter-Varying Parallel-Cascade Identification of Dynamic Joint Stiffness.

OBJECTIVE: The paper presents a method to identify ankle joint dynamic stiffness during functional tasks where intrinsic and reflex stiffness change with a time-varying scheduling variable (SV), such as joint position or torque. METHODS: The method models joint stiffness with two pathways: (1) A parameter-varying (PV) impulse response function (IRF) describing intrinsic stiffness; and (2) a reflex stiffness model comprising a PV static nonlinearity followed by a PV linear element. RESULTS: Monte-Carlo simulations demonstrated that the method accurately estimated all elements of the intrinsic and reflex pathways as they changed with a SV. Experimental results with a healthy individual subjected to large, imposed ankle movements demonstrated that: (a) Intrinsic stiffness changed substantially as a function of ankle position; elasticity was lowest near the mid-position and increased with either dorsiflexion or plantarflexion. (b) Reflex gain increased and the velocity threshold for reflex excitation decreased monotonically with ankle dorsiflexion. (c) Reflex dynamics resembled a second-order, low-pass system that was invariant with ankle position. (d) The identified PV Parallel-Cascade (PC) model accurately predicted the torque response to novel trajectories of ankle movement. CONCLUSION: The PV-PC method can accurately and reliably estimate how intrinsic and reflex stiffness change with a time-varying SV. SIGNIFICANCE: The method is novel with multiple advantages: (a) It provides a unified algorithm that characterizes the changes in the parameters of all joint stiffness elements needed to understand their role in postural/movement control; (b) It is efficient requiring only two trials; (c) The models identified can predict the joint stiffness response to novel movements informing orthoses and prostheses design.

Linear Parameter Varying Identification of Dynamic Joint Stiffness during Time-Varying Voluntary Contractions.

Dynamic joint stiffness is a dynamic, nonlinear relationship between the position of a joint and the torque acting about it, which can be used to describe the biomechanics of the joint and associated limb(s). This paper models and quantifies changes in ankle dynamic stiffness and its individual elements, intrinsic and reflex stiffness, in healthy human subjects during isometric, time-varying (TV) contractions of the ankle plantarflexor muscles. A subspace, linear parameter varying, parallel-cascade (LPV-PC) algorithm was used to identify the model from measured input position perturbations and output torque data using voluntary torque as the LPV scheduling variable (SV). Monte-Carlo simulations demonstrated that the algorithm is accurate, precise, and robust to colored measurement noise. The algorithm was then used to examine stiffness changes associated with TV isometric contractions. The SV was estimated from the Soleus EMG using a Hammerstein model of EMG-torque dynamics identified from unperturbed trials. The LPV-PC algorithm identified (i) a non-parametric LPV impulse response function (LPV IRF) for intrinsic stiffness and (ii) a LPV-Hammerstein model for reflex stiffness consisting of a LPV static nonlinearity followed by a time-invariant state-space model of reflex dynamics. The results demonstrated that: (a) intrinsic stiffness, in particular ankle elasticity, increased significantly and monotonically with activation level; (b) the gain of the reflex pathway increased from rest to around 10-20% of subject's MVC and then declined; and (c) the reflex dynamics were second order. These findings suggest that in healthy human ankle, reflex stiffness contributes most at low muscle contraction levels, whereas, intrinsic contributions monotonically increase with activation level.

Ankle Joint Intrinsic Dynamics is More Complex than a Mass-Spring-Damper Model.

This paper describes a new small signal parametric model of ankle joint intrinsic mechanics in normal subjects. We found that intrinsic ankle mechanics is a third-order system and the second-order mass-spring-damper model, referred to as IBK, used by many researchers in the literature cannot adequately represent ankle dynamics at all frequencies in a number of important tasks. This was demonstrated using experimental data from five healthy subjects with no voluntary muscle contraction and at seven ankle positions covering the range of motion. We showed that the difference between the new third-order model and the conventional IBK model increased from dorsi to plantarflexed position. The new model was obtained using a multi-step identification procedure applied to experimental input/output data of the ankle joint. The procedure first identifies a non-parametric model of intrinsic joint stiffness where ankle position is the input and torque is the output. Then, in several steps, the model is converted into a continuous-time transfer function of ankle compliance, which is the inverse of stiffness. Finally, we showed that the third-order model is indeed structurally consistent with agonist-antagonist musculoskeletal structure of human ankle, which is not the case for the IBK model.

Identification of time-varying dynamics of reflex EMG in the ankle plantarflexors during time-varying, isometric contractions.

The dynamic relationship between the joint position and reflex EMG in ankle muscles of healthy human subjects was studied for time-varying (TV) contractions. A linear parameter varying (LPV) identification algorithm was used to estimate the Hammerstein system relating ankle position to the reflex EMG response. The estimated Hammerstein system comprised a time-invariant (TI) linear element and a TV static nonlinearity that resembled a half-wave rectifier with a threshold and linear gain. The results demonstrated a systematic change in the reflex nonlinearity with the activation level. The gain of TV nonlinearity increased with activation level reaching its peak at 20-30% maximum voluntary contraction and then decreased. The threshold of the nonlinearity decreased with increasing activation level reaching it minimum at the same point where the gain was maximal. Using the LPV-Hammerstein method in this work, the underlying TV dynamics were extracted from small number of trials. Thus, this method can be used to study stretch reflexes in subjects with neuromuscular disorders.

Linear parameter varying identification of ankle joint intrinsic stiffness during imposed walking movements.

This paper describes a novel model structure and identification method for the time-varying, intrinsic stiffness of human ankle joint during imposed walking (IW) movements. The model structure is based on the superposition of a large signal, linear, time-invariant (LTI) model and a small signal linear-parameter varying (LPV) model. The methodology is based on a two-step algorithm; the LTI model is first estimated using data from an unperturbed IW trial. Then, the LPV model is identified using data from a perturbed IW trial with the output predictions of the LTI model removed from the measured torque. Experimental results demonstrate that the method accurately tracks the continuous-time variation of normal ankle intrinsic stiffness when the joint position changes during the IW movement. Intrinsic stiffness gain decreases from full plantarflexion to near the mid-point of plantarflexion and then increases substantially as the ankle is dosriflexed.

A novel algorithm for linear parameter varying identification of Hammerstein systems with time-varying nonlinearities.

This paper describes a novel method for the identification of Hammerstein systems with time-varying (TV) static nonlinearities and time invariant (TI) linear elements. This paper develops a linear parameter varying (LPV) state-space representation for such systems and presents a subspace identification technique that gives individual estimates of the Hammerstein components. The identification method is validated using simulated data of a TV model of ankle joint reflex stiffness where the threshold and gain of the model change as nonlinear functions of an exogenous signal. Pilot experiment of TV reflex EMG response identification in normal ankle joint during an imposed walking task demonstrate systematic changes in the reflex nonlinearity with the trajectory of joint position.