Estimation of musculotendon kinematics in large musculoskeletal models using multidimensional B-splines

We present a robust and computationally inexpensive method to estimate the lengths and three-dimensional moment arms for a large number of musculotendon actuators of the human lower limb. Using a musculoskeletal model of the lower extremity, a set of values was established for the length of each musculotendon actuator for different lower limb generalized coordinates (joint angles). A multidimensional spline function was then used to fit these data. Muscle moment arms were obtained by differentiating the musculotendon length spline function with respect to the generalized coordinate of interest. This new method was then compared to a previously used polynomial regression method. Compared to the polynomial regression method, the multidimensional spline method produced lower errors for estimating musculotendon lengths and moment arms throughout the whole generalized coordinate workspace. The fitting accuracy was also less affected by the number of dependent degrees of freedom and by the amount of experimental data available. The spline method only required information on musculotendon lengths to estimate both musculotendon lengths and moment arms, thus relaxing data input requirements, whereas the polynomial regression requires different equations to be used for both musculotendon lengths and moment arms. Finally, we used the spline method in conjunction with an electromyography driven musculoskeletal model to estimate muscle forces under different contractile conditions, which showed that the method is suitable for the integration into large scale neuromusculoskeletal models.     

Moment arms and musculotendon lengths estimation for a three-dimensional lower-limb model

This paper presents a set of polynomial expressions that can be used as regression equations to estimate length and three-dimensional moment arms of 43 lower-limb musculotendon actuators. These equations allow one to find, at a low computational cost, the musculotendon geometric parameters required for numerical simulation of large musculoskeletal models. Nominal values for these biomechanical parameters were established using a public-domain musculoskeletal model of the lower limb (IEEE Trans. Biomed. Eng. 37 (1990) 757). To fit these nominal values, regression equations with different levels of complexity were generated, based on the number of generalized coordinates of the joints spanned by each musculotendon actuator. Least squares fitting was used to identify regression equation coefficients. The goodness of the fit and confidence intervals were assessed, and the best fitting equations selected.     

Musculotendon lengths and moment arms for a three-dimensional upper-extremity model

Generating muscle-driven forward dynamics simulations of human movement using detailed musculoskeletal models can be computationally expensive. This is due in part to the time required to calculate musculotendon geometry (e.g., musculotendon lengths and moment arms), which is necessary to determine and apply individual musculotendon forces during the simulation. Modeling upper-extremity musculotendon geometry can be especially challenging due to the large number of multi-articular muscles and complex muscle paths. To accurately represent this geometry, wrapping surface algorithms and/or other computationally expensive techniques (e.g., phantom segments) are used. This paper provides a set of computationally efficient polynomial regression equations that estimate musculotendon length and moment arms for thirty-two (32) upper-extremity musculotendon actuators representing the major muscles crossing the shoulder, elbow and wrist joints. Equations were developed using a least squares fitting technique based on geometry values obtained from a validated public-domain upper-extremity musculoskeletal model that used wrapping surface elements (Holzbaur et al., 2005). In general, the regression equations fit well the original model values, with an average root mean square difference for all musculotendon actuators over the represented joint space of 0.39 mm (1.1% of peak value). In addition, the equations reduced the computational time required to simulate a representative upper-extremity movement (i.e., wheelchair propulsion) by more than two orders of magnitude (315 versus 2.3 s). Thus, these equations can assist in generating computationally efficient forward dynamics simulations of a wide range of upper-extremity movements.     

The gambin model provides a superior fit to species abundance distributions with a single free parameter: evidence,implementation and interpretation

The species abundance distribution (SAD) has been a central focus of community ecology for over fifty years, and is currently the subject of widespread renewed interest. The gambin model has recently been proposed as a model that provides a superior fit to commonly preferred SAD models. It has also been argued that the model's single parameter (α) presents a potentially informative ecological diversity metric, because it summarises the shape of the SAD in a single number. Despite this potential, few empirical tests of the model have been undertaken, perhaps because the necessary methods and software for fitting the model have not existed. Here, we derive a maximum likelihood method to fit the model, and use it to undertake a comprehensive comparative analysis of the fit of the gambin model. The functions and computational code to fit the model are incorporated in a newly developed free‐to‐download R package (gambin). We test the gambin model using a variety of datasets and compare the fit of the gambin model to fits obtained using the Poisson lognormal, logseries and zero‐sum multinomial distributions. We found that gambin almost universally provided a better fit to the data and that the fit was consistent for a variety of sample grain sizes. We demonstrate how α can be used to differentiate intelligibly between community structures of Azorean arthropods sampled in different land use types. We conclude that gambin presents a flexible model capable of fitting a wide variety of observed SAD data, while providing a useful index of SAD form in its single fitted parameter. As such, gambin has wide potential applicability in the study of SADs, and ecology more generally.     

The influence of model parameters on model validation

The study examined the sensitivity of two musculoskeletal models to the parameters describing each model. Two different models were examined: a phenomenological model of human jumping with parameters based on live subject data, and the second a model of the First Dorsal Interosseous with parameters based on cadaveric measurements. Both models were sensitive to the model parameters, with the use of mean group data not producing model outputs reflective of either the performance of any group member or the mean group performance. These results highlight the value of subject specific model parameters, and the problems associated with model validation.     

An approach to examining model dependence in EM reconstructions using cross-validation

Reference bias refers to a common problem in fitting experimental data to an initial model. Given enough free parameters, a good fit of any experimental data to the model can be obtained, even if the experimental data contain only noise. Reference-based alignment methods used in electron microscopy (EM) are subject to this type of bias, in that images containing pure noise can regenerate the reference. Cross-validation is based on the idea that the experimental data used to assess the validity of the fitting should not be the same data as were used to do the fitting. Here we present the application of cross-validation to one form of reference-based alignment: 3D-projection matching in single-particle reconstructions. Our results show that reference bias is indeed present in reconstructions, but that the effect is small for real data compared to that for random noise, and that this difference in behavior is magnified, rather than diminished, during iterative refinement.     

Computation of a stabilizing set of feedback matrices of a large-scale nonlinear musculoskeletal dynamic model

The purpose of this study is to present a general mathematical framework to compute a set of feedback matrices which stabilize an unstable nonlinear anthropomorphic musculoskeletal dynamic model. This method is activity specific and involves four fundamental stages. First, from muscle activation data (input) and motion degrees-of-freedom (output) a dynamic experimental model is obtained using system identification schemes. Second, a nonlinear musculoskeletal dynamic model which contains the same number of muscles and degrees-of-freedom and best represents the activity being considered is proposed. Third, the nonlinear musculoskeletal model (anthropomorphic model) is replaced by a family of linear systems, parameterized by the same set of input/output data (nominal points) used in the identification of the experimental model. Finally, a set of stabilizing output feedback matrices, parameterized again by the same set of nominal points, is computed such that when combined with the anthropomorphic model, the combined system resembles the structural form of the experimental model. The method is illustrated in regard to the human squat activity.     

Sensitivity of model predictions of muscle function to changes in moment arms and muscle-tendon properties: a Monte-Carlo analysis

Hill-type muscle models are commonly used in musculoskeletal models to estimate muscle forces during human movement. However, the sensitivity of model predictions of muscle function to changes in muscle moment arms and muscle-tendon properties is not well understood. In the present study, a three-dimensional muscle-actuated model of the body was used to evaluate the sensitivity of the function of the major lower limb muscles in accelerating the whole-body center of mass during gait. Monte-Carlo analyses were used to quantify the effects of entire distributions of perturbations in the moment arms and architectural properties of muscles. In most cases, varying the moment arm and architectural properties of a muscle affected the torque generated by that muscle about the joint(s) it spanned as well as the torques generated by adjacent muscles. Muscle function was most sensitive to changes in tendon slack length and least sensitive to changes in muscle moment arm. However, the sensitivity of muscle function to changes in moment arms and architectural properties was highly muscle-specific; muscle function was most sensitive in the cases of gastrocnemius and rectus femoris and insensitive in the cases of hamstrings and the medial sub-region of gluteus maximus. The sensitivity of a muscle's function was influenced by the magnitude of the muscle's force as well as the operating region of the muscle on its force-length curve. These findings have implications for the development of subject-specific models of the human musculoskeletal system.     

In vivo parameter identification in arteries considering multiple levels of smooth muscle activity

In this paper an existing in vivo parameter identification method for arteries is extended to account for smooth muscle activity. Within this method a continuum-mechanical model, whose parameters relate to the mechanical properties of the artery, is fit to clinical data by solving a minimization problem. Including smooth muscle activity in the model increases the number of parameters. This may lead to overparameterization, implying that several parameter combinations solve the minimization problem equally well and it is therefore not possible to determine which set of parameters represents the mechanical properties of the artery best. To prevent overparameterization the model is fit to clinical data measured at different levels of smooth muscle activity. Three conditions are considered for the human abdominal aorta: basal during rest; constricted, induced by lower-body negative pressure; and dilated, induced by physical exercise. By fitting the model to these three arterial conditions simultaneously a unique set of model parameters is identified and the model prediction agrees well with the clinical data.

     

A technical method using musculoskeletal model to analyse dynamic properties of muscles during human movement

An effective way to avoid invading or injuring the subjects is to use the musculoskeletal model when studying the dynamic properties of muscles in vivo. So, we put forward a joint coordinate system-based method, which mainly focuses on the coordinate's transformation of corresponding muscle attachment points, respectively, in the model and the subject in order to reproduce the movement of the subject on the model. As muscle moment arm is usually used to evaluate the dynamic properties of muscles, the moment arms in elbow flexion for each of the major muscles crossing the elbow in 50 healthy subjects (25 males and 25 females), ranging in height from 1.50 to 1.80 m (mean 1.6542 m) are calculated and compared with the measured data obtained from anatomical studies reported in the literature. The trends of the value basically coincide with each other. So, this novel method can be valid.     

A model for semi-quantitative studies of muscle actions

A determinate scheme for modeling the mechanical effects of muscle contractions in a musculoskeletal system, using the direct-stiffness method of structural analysis, is described. Data concerning skeletal geometries, connective tissue passive mechanical properties, and muscle lines-of-action and cross-sectional areas are incorporated to make the model specific for studies of muscle actions in the human trunk. Several illustrative examples of model trunk responses to muscle contractions are presented and discussed.     

Estimating a predator-prey dynamical model with the parameter cascades method

Summary .   Ordinary differential equations (ODEs) are widely used in ecology to describe the dynamical behavior of systems of interacting populations. However, systems of ODEs rarely provide quantitative solutions that are close to real field observations or experimental data because natural systems are subject to environmental and demographic noise and ecologists are often uncertain about the correct parameterization. In this article we introduce "parameter cascades" as an improved method to estimate ODE parameters such that the corresponding ODE solutions fit the real data well. This method is based on the modified penalized smoothing with the penalty defined by ODEs and a generalization of profiled estimation, which leads to fast estimation and good precision for ODE parameters from noisy data. This method is applied to a set of ODEs originally developed to describe an experimental predator–prey system that undergoes oscillatory dynamics. The new parameterization considerably improves the fit of the ODE model to the experimental data sets. At the same time, our method reveals that important structural assumptions that underlie the original ODE model are essentially correct. The mathematical formulations of the two nonlinear interaction terms (functional responses) that link the ODEs in the predator–prey model are validated by estimating the functional responses nonparametrically from the real data. We suggest two major applications of "parameter cascades" to ecological modeling: It can be used to estimate parameters when original data are noisy, missing, or when no reliable priori estimates are available; it can help to validate the structural soundness of the mathematical modeling approach.     

A validated model of passive muscle in compression

A better characterisation of soft tissues is required to improve the accuracy of human body models used, amongst other applications, for virtual crash modelling. This paper presents a theoretical model and the results of an experimental procedure to characterise the quasi-static, compressive behaviour of skeletal muscle in three dimensions. Uniaxial, unconstrained compression experiments have been conducted on aged and fresh animal muscle samples oriented at various angles from the fibre direction. A transversely isotropic hyperelastic model and a model using the theory of transverse isotropy and strain dependent Young's moduli (SYM) have been fitted to the experimental data. Results show that the hyperelastic model does not adequately fit the data in all directions of testing. In contrast, the SYM gives a good fit to the experimental data in both the fibre and cross-fibre direction, up to 30% strain for aged samples. The model also yields good prediction of muscle behaviour at 45° from the fibre direction. Fresh samples show a different behaviour than aged tissues at 45° from the fibre direction. However, the SYM is able to capture this difference and gives a good fit to the experimental data in the fibre, the cross-fibre and at 45° from the fibre direction. The model also yields good prediction of muscle behaviour when compressed at 30° and 60° from the fibre direction. The effect of the time of test after death has also been investigated. Significant stiffening of muscle behaviour is noted a few hours after death of the subject.     

An entropy-assisted musculoskeletal shoulder model

Optimization combined with a musculoskeletal shoulder model has been used to estimate mechanical loading of musculoskeletal elements around the shoulder. Traditionally, the objective function is to minimize the summation of the total activities of the muscles with forces, moments, and stability constraints. Such an objective function, however, tends to neglect the antagonist muscle co-contraction. In this study, an objective function including an entropy term is proposed to address muscle co-contractions. A musculoskeletal shoulder model is developed to apply the proposed objective function. To find the optimal weight for the entropy term, an experiment was conducted. In the experiment, participants generated various 3-D shoulder moments in six shoulder postures. The surface EMG of 8 shoulder muscles was measured and compared with the predicted muscle activities based on the proposed objective function using Bhattacharyya distance and concordance ratio under different weight of the entropy term. The results show that a small weight of the entropy term can improve the predictability of the model in terms of muscle activities. Such a result suggests that the concept of entropy could be helpful for further understanding the mechanism of muscle co-contractions as well as developing a shoulder biomechanical model with greater validity.     

Estimation of muscle response using three-dimensional musculoskeletal models before impact situation: a simulation study

When car crash experiments are performed using cadavers or dummies, the active muscles' reaction on crash situations cannot be observed. The aim of this study is to estimate muscles' response of the major muscle groups using three-dimensional musculoskeletal model by dynamic simulations of low-speed sled-impact. The three-dimensional musculoskeletal models of eight subjects were developed, including 241 degrees of freedom and 86 muscles. The muscle parameters considering limb lengths and the force-generating properties of the muscles were redefined by optimization to fit for each subject. Kinematic data and external forces measured by motion tracking system and dynamometer were then input as boundary conditions. Through a least-squares optimization algorithm, active muscles' responses were calculated during inverse dynamic analysis tracking the motion of each subject. Electromyography for major muscles at elbow, knee, and ankle joints was measured to validate each model. For low-speed sled-impact crash, experiment and simulation with optimized and unoptimized muscle parameters were performed at 9.4 m/h and 10 m/h and muscle activities were compared among them. The muscle activities with optimized parameters were closer to experimental measurements than the results without optimization. In addition, the extensor muscle activities at knee, ankle, and elbow joint were found considerably at impact time, unlike previous studies using cadaver or dummies. This study demonstrated the need to optimize the muscle parameters to predict impact situation correctly in computational studies using musculoskeletal models. And to improve accuracy of analysis for car crash injury using humanlike dummies, muscle reflex function, major extensor muscles' response at elbow, knee, and ankle joints, should be considered.     

A technical method using musculoskeletal model to analyse dynamic properties of muscles during human movement

An effective way to avoid invading or injuring the subjects is to use the musculoskeletal model when studying the dynamic properties of muscles in vivo. So, we put forward a joint coordinate system-based method, which mainly focuses on the coordinate's transformation of corresponding muscle attachment points, respectively, in the model and the subject in order to reproduce the movement of the subject on the model. As muscle moment arm is usually used to evaluate the dynamic properties of muscles, the moment arms in elbow flexion for each of the major muscles crossing the elbow in 50 healthy subjects (25 males and 25 females), ranging in height from 1.50 to 1.80 m (mean 1.6542 m) are calculated and compared with the measured data obtained from anatomical studies reported in the literature. The trends of the value basically coincide with each other. So, this novel method can be valid.     

A stochastic model of elbow flexion strength for subjects with and without long head biceps tear

The classical approach of musculoskeletal modeling is to predict muscle forces and joint torques with a deterministic model constructed from parameters of an average subject. However, this type of model does not perform well for outliers, and does not model the effects of parameter variability. In this study, a Monte-Carlo model was used to stochastically simulate the effects of variability in musculoskeletal parameters on elbow flexion strength in healthy normals, and in subjects with long head biceps (LHB) rupture. The goal was to determine if variability in elbow flexion strength could be quantifiably explained with variability in musculoskeletal parameters. Parameter distributions were constructed from data in the literature. Parameters were sampled from these distributions and used to predict muscle forces and joint torques. The median and distribution of measured joint torque was predicted with small errors ( < 5%). Muscle forces for both cases were predicted and compared. In order to predict measured torques for the case of LHB rupture, the median force and mean cross-sectional area in the remaining elbow flexor muscles is greater than in healthy normals. The probabilities that muscle forces for the Tear case exceed median muscle forces for the No-Tear case are 0.98, 0.99 and 0.79 for SH Biceps, brachialis and brachioradialis, respectively. Differences in variability of measured torques for the two cases are explained by differences in parameter variability.     

An EMG-driven musculoskeletal model to estimate muscle forces and knee joint moments in vivo

This paper examined if an electromyography (EMG) driven musculoskeletal model of the human knee could be used to predict knee moments, calculated using inverse dynamics, across a varied range of dynamic contractile conditions. Muscle-tendon lengths and moment arms of 13 muscles crossing the knee joint were determined from joint kinematics using a three-dimensional anatomical model of the lower limb. Muscle activation was determined using a second-order discrete non-linear model using rectified and low-pass filtered EMG as input. A modified Hill-type muscle model was used to calculate individual muscle forces using activation and muscle tendon lengths as inputs. The model was calibrated to six individuals by altering a set of physiologically based parameters using mathematical optimisation to match the net flexion/extension (FE) muscle moment with those measured by inverse dynamics. The model was calibrated for each subject using 5 different tasks, including passive and active FE in an isokinetic dynamometer, running, and cutting manoeuvres recorded using three-dimensional motion analysis. Once calibrated, the model was used to predict the FE moments, estimated via inverse dynamics, from over 200 isokinetic dynamometer, running and sidestepping tasks. The inverse dynamics joint moments were predicted with an average R(2) of 0.91 and mean residual error of approximately 12 Nm. A re-calibration of only the EMG-to-activation parameters revealed FE moments prediction across weeks of similar accuracy. Changing the muscle model to one that is more physiologically correct produced better predictions. The modelling method presented represents a good way to estimate in vivo muscle forces during movement tasks.     

A stochastic model of elbow flexion strength for subjects with and without long head biceps tear

The classical approach of musculoskeletal modeling is to predict muscle forces and joint torques with a deterministic model constructed from parameters of an average subject. However, this type of model does not perform well for outliers, and does not model the effects of parameter variability. In this study, a Monte-Carlo model was used to stochastically simulate the effects of variability in musculoskeletal parameters on elbow flexion strength in healthy normals, and in subjects with long head biceps (LHB) rupture. The goal was to determine if variability in elbow flexion strength could be quantifiably explained with variability in musculoskeletal parameters. Parameter distributions were constructed from data in the literature. Parameters were sampled from these distributions and used to predict muscle forces and joint torques. The median and distribution of measured joint torque was predicted with small errors (< 5%). Muscle forces for both cases were predicted and compared. In order to predict measured torques for the case of LHB rupture, the median force and mean cross-sectional area in the remaining elbow flexor muscles is greater than in healthy normals. The probabilities that muscle forces for the Tear case exceed median muscle forces for the No-Tear case are 0.98, 0.99 and 0.79 for SH Biceps, brachialis and brachioradialis, respectively. Differences in variability of measured torques for the two cases are explained by differences in parameter variability.     

A validated model of passive skeletal muscle to predict force and intramuscular pressure

The passive properties of skeletal muscle are often overlooked in muscle studies, yet they play a key role in tissue function in vivo. Studies analyzing and modeling muscle passive properties, while not uncommon, have never investigated the role of fluid content within the tissue. Additionally, intramuscular pressure (IMP) has been shown to correlate with muscle force in vivo and could be used to predict muscle force in the clinic. In this study, a novel model of skeletal muscle was developed and validated to predict both muscle stress and IMP under passive conditions for the New Zealand White Rabbit tibialis anterior. This model is the first to include fluid content within the tissue and uses whole muscle geometry. A nonlinear optimization scheme was highly effective at fitting model stress output to experimental stress data (normalized mean square error or NMSE fit value of 0.993) and validation showed very good agreement to experimental data (NMSE fit values of 0.955 and 0.860 for IMP and stress, respectively). While future work to include muscle activation would broaden the physiological application of this model, the passive implementation could be used to guide surgeries where passive muscle is stretched.