Development of a mathematical model for predicting electrically elicited quadriceps femoris muscle forces during isovelocity knee joint motion

<p>Abstract</p> <p>Background</p> <p>Direct electrical activation of skeletal muscles of patients with upper motor neuron lesions can restore functional movements, such as standing or walking. Because responses to electrical stimulation are highly nonlinear and time varying, accurate control of muscles to produce functional movements is very difficult. Accurate and predictive mathematical models can facilitate the design of stimulation patterns and control strategies that will produce the desired force and motion. In the present study, we build upon our previous isometric model to capture the effects of constant angular velocity on the forces produced during electrically elicited concentric contractions of healthy human quadriceps femoris muscle. Modelling the isovelocity condition is important because it will enable us to understand how our model behaves under the relatively simple condition of constant velocity and will enable us to better understand the interactions of muscle length, limb velocity, and stimulation pattern on the force produced by the muscle.</p> <p>Methods</p> <p>An additional term was introduced into our previous isometric model to predict the force responses during constant velocity limb motion. Ten healthy subjects were recruited for the study. Using a KinCom dynamometer, isometric and isovelocity force data were collected from the human quadriceps femoris muscle in response to a wide range of stimulation frequencies and patterns. % error, linear regression trend lines, and paired t-tests were used to test how well the model predicted the experimental forces. In addition, sensitivity analysis was performed using Fourier Amplitude Sensitivity Test to obtain a measure of the sensitivity of our model's output to changes in model parameters.</p> <p>Results</p> <p>Percentage RMS errors between modelled and experimental forces determined for each subject at each stimulation pattern and velocity showed that the errors were in general less than 20%. The coefficients of determination between the measured and predicted forces show that the model accounted for ~86% and ~85% of the variances in the measured force-time integrals and peak forces, respectively.</p> <p>Conclusion</p> <p>The range of predictive abilities of the isovelocity model in response to changes in muscle length, velocity, and stimulation frequency for each individual make it ideal for dynamic applications like FES cycling.</p

Conjectures on exact solution of three - dimensional (3D) simple orthorhombic Ising lattices

We report the conjectures on the three-dimensional (3D) Ising model on simple orthorhombic lattices, together with the details of calculations for a putative exact solution. Two conjectures, an additional rotation in the fourth curled-up dimension and the weight factors on the eigenvectors, are proposed to serve as a boundary condition to deal with the topologic problem of the 3D Ising model. The partition function of the 3D simple orthorhombic Ising model is evaluated by spinor analysis, by employing these conjectures. Based on the validity of the conjectures, the critical temperature of the simple orthorhombic Ising lattices could be determined by the relation of KK* = KK' + KK'' + K'K'' or sinh 2K sinh 2(K' + K'' + K'K''/K) = 1. For a simple cubic Ising lattice, the critical point is putatively determined to locate exactly at the golden ratio xc = exp(-2Kc) = (sq(5) - 1)/2, as derived from K* = 3K or sinh 2K sinh 6K = 1. If the conjectures would be true, the specific heat of the simple orthorhombic Ising system would show a logarithmic singularity at the critical point of the phase transition. The spontaneous magnetization and the spin correlation functions of the simple orthorhombic Ising ferromagnet are derived explicitly. The putative critical exponents derived explicitly for the simple orthorhombic Ising lattices are alpha = 0, beta = 3/8, gamma = 5/4, delta = 13/3, eta = 1/8 and nu = 2/3, showing the universality behavior and satisfying the scaling laws. The cooperative phenomena near the critical point are studied and the results obtained based on the conjectures are compared with those of the approximation methods and the experimental findings. The 3D to 2D crossover phenomenon differs with the 2D to 1D crossover phenomenon and there is a gradual crossover of the exponents from the 3D values to the 2D ones.Comment: 176 pages, 4 figure