nils-herrmann/clean_scirepeval huggingface dataset:
Muscle force estimation in clinical gait analysis using AnyBody and OpenSim.
Venue: Journal of biomechanics

A variety of musculoskeletal models are applied in different modelling environments for estimating muscle forces during gait. Influence of different modelling assumptions and approaches on model outputs are still not fully understood, while direct comparisons of standard approaches have been rarely undertaken. This study seeks to compare joint kinematics, joint kinetics and estimated muscle forces of two standard approaches offered in two different modelling environments (AnyBody, OpenSim). It is hypothesised that distinctive differences exist for individual muscles, while summing up synergists show general agreement. Experimental data of 10 healthy participants (28 ± 5 years, 1.72 ± 0.08 m, 69 ± 12 kg) was used for a standard static optimisation muscle force estimation routine in AnyBody and OpenSim while using two gait-specific musculoskeletal models. Statistical parameter mapping paired t-test was used to compare joint angle, moment and muscle force waveforms in Matlab. Results showed differences especially in sagittal ankle and hip angles as well as sagittal knee moments. Differences were also found for some of the muscles, especially of the triceps surae group and the biceps femoris short head, which occur as a result of different anthropometric and anatomical definitions (mass and inertia of segments, muscle properties) and scaling procedures (static vs. dynamic). Understanding these differences and their cause is crucial to operate such modelling environments in a clinical setting. Future research should focus on alternatives to classical generic musculoskeletal models (e.g. implementation of functional calibration tasks), while using experimental data reflecting normal and pathological gait to gain a better understanding of variations and divergent behaviour between approaches.
Distinguishing between typical and atypical motion patterns amongst healthy individuals during a constrained spine flexion task.
Venue: Journal of biomechanics

Despite 'abnormal' motion being considered a risk factor for low back injury, the current understanding of 'normal' spine motion is limited. Identifying normal motion within an individual is complicated by the considerable variation in movement patterns amongst healthy individuals. Therefore, the purpose of this study was to characterize sources of variation in spine motion among a sample of healthy participants. The second objective of this study was to develop a multivariate model capable of predicting an expected movement pattern for an individual. The kinematic shape of the lower thoracic and lumbar spine was recorded during a constrained dynamic trunk flexion movement; as this is not a normal everyday movement task, movements are considered 'typical' and 'atypical' for this task rather than 'normal' and 'abnormal'. Variations in neutral standing posture accounted for 85% of the variation in spine motion throughout the task. Differences in total spine range of flexion and a regional re-weighting of range of motion between lower thoracic and lumbar regions explained a further 9% of the variance among individuals. The analysis also highlighted a difference in temporal sequencing of motion between lower thoracic and lumbar regions which explained 2% of the total movement variation. These identified sources of variation were used to select independent variables for a multivariate linear model capable of predicting an individuals' expected movement pattern. This was done as a proof-of-concept to demonstrate how the error between predicted and observed motion patterns could be used to differentiate between 'typical' and 'atypical' movement strategies.
Direct and indirect effects of joint torque inputs during an induced speed analysis of a swinging motion.
Venue: Journal of biomechanics

This study proposed a method to quantify direct and indirect effects of the joint torque inputs in the speed-generating mechanism of a swinging motion. Linear and angular accelerations of all segments within a multi-linked system can be expressed as the sum of contributions from a joint torque term, gravitational force term and motion-dependent term (MDT), where the MDT is a nonlinear term consisting of centrifugal force, Coriolis force and gyroscopic effect moment components. Direct effects result from angular accelerations induced by a joint torque at a given instant, whereas indirect effects arise through the MDT induced by joint torques exerted in the past. These two effects were quantified for the kicking-side leg during a rugby place kick. The MDT was the largest contributor to the foot centre of gravity (CG)'s speed at ball contact. Of the factors responsible for generating the MDT, the direct and indirect effects of the hip flexion-extension torque during both the flight phase (from the final kicking foot take-off to support foot contact) and the subsequent support phase (from support foot contact to ball contact) were important contributors to the foot CG's speed at ball contact. The indirect effect of the ankle plantar-dorsal flexion torque and the direct effect of the knee flexion-extension torque during the support phase showed the largest positive and negative contributions to the foot CG's speed at ball contact, respectively. The proposed method allows the identification of which individual joint torque axes are crucial and the timings of joint torque exertion that are used to generate a high speed of the distal point of a multi-linked system.
Reliability and robustness of muscle architecture measurements obtained using diffusion tensor imaging with anatomically constrained tractography.
Venue: Journal of biomechanics

For detailed analyses of muscle adaptation mechanisms during growth, ageing or disease, reliable measurements of muscle architecture are required. Diffusion tensor imaging (DTI) and DTI tractography have been used to reconstruct the architecture of human muscles in vivo. However, muscle architecture measurements reconstructed with conventional DTI techniques are often anatomically implausible because the reconstructed fascicles do not terminate on aponeuroses, as real muscle fascicles are known to do. In this study, we tested the reliability of an anatomically constrained DTI-based method for measuring three-dimensional muscle architecture. Anatomical magnetic resonance images and diffusion tensor images were obtained from the left legs of eight healthy participants on two occasions one week apart. Muscle volumes, fascicle lengths, pennation angles and fascicle curvatures were measured in the medial and lateral gastrocnemius, soleus and the tibialis anterior muscles. Averaged across muscles, the intraclass correlation coefficient was 0.99 for muscle volume, 0.81 for fascicle length, 0.73 for pennation angle and 0.76 for fascicle curvature. Measurements of muscle architecture obtained using conventional DTI tractography were highly sensitive to variations in the stopping criteria for DTI tractography. The application of anatomical constraints reduced this sensitivity significantly. This study demonstrates that anatomically constrained DTI tractography can provide reliable and robust three-dimensional measurements of whole-muscle architecture. The algorithms used to constrain tractography have been made publicly available.
A comparison of metatarsophalangeal joint center locations on estimated joint moments during running.
Venue: Journal of biomechanics

The forefoot functions as the base of support during late stance, rotating about the dual-axis of the metatarsophalangeal joints. Previous research has shown that joint axis definition affects estimated joint moments about the forefoot. However, little is known about how metatarsophalangeal joint center definition affects estimated joint kinetics. This study compared moments about the metatarsophalangeal joint using four different defined joint centers. There was a significant difference (p < .001) in peak moments between joint center definitions, differing by up to 0.488 N-m/kg for the slow and 0.878 N-m/kg for the fast running speeds tested. Additionally, there was a significant difference (p < .001) for when peak plantar flexor moment occurred during the slower running condition. The more posteriorly oriented joint centers resulted in higher moments and earlier onset of the plantar flexor moment. In addition to careful modeling of the metatarsophalangeal joint axis, it is recommended that joint center definition should be considered as well.
A principal component analysis-based framework for statistical modeling of bone displacement during wrist maneuvers.
Venue: Journal of biomechanics

We present a method for the statistical modeling of the displacements of wrist bones during the performance of coordinated maneuvers, such as radial-ulnar deviation (RUD). In our approach, we decompose bone displacement via a set of basis functions, identified via principal component analysis (PCA). We utilized MRI wrist scans acquired at multiple static positions for deriving these basis functions. We then utilized these basis functions to compare the displacements undergone by the bones of the left versus right wrist in the same individual, and between bones of the wrists of men and women, during the performance of the coordinated RUD maneuver. Our results show that the complex displacements of the wrist bones during RUD can be modeled with high reliability with just 5 basis functions, that captured over 91% of variation across individuals. The basis functions were able to predict intermediate wrist bone poses with an overall high accuracy (mean error of 0.26 mm). Our proposed approach found statistically significant differences between bone displacement trajectories in women versus men, however, did not find significant differences in those of the left versus right wrist in the same individual. Our proposed method has the potential to enable detailed analysis of wrist kinematics for each sex, and provide a robust framework for characterizing the normal and pathologic displacement of the wrist bones, such as in the context of wrist instability.
Patella tendon moment arm function considerations for human vastus lateralis force estimates.
Venue: Journal of biomechanics

In vivo muscle forces are typically estimated using literature-based or subject-specific moment arms (MAs) because it is not possible to measure in vivo muscle forces non-invasively. However, even subject-specific muscle-tendon MAs vary across contraction levels and are impossible to determine at high contraction levels without techniques that use ionized radiation. Therefore, different generic MA functions are often used to estimate in vivo muscle forces, which may alter force predictions and the shape of the muscle's force-length relationship. The aim of this study was to examine the influence of different literature-based patella tendon MA functions on the vastus lateralis (VL) force-angle relationship. Participants (n = 11) performed maximum voluntary isometric knee extension contractions at six knee flexion angles, ranging from 40° to 90°. To estimate in vivo VL muscle force, the peak knee extension torque at each joint angle was multiplied by the VL's physiological cross-sectional area (PCSA) relative to the quadriceps' PCSA (34%) and then divided by the angle-specific patella tendon MA for 19 different functions. Maximum VL force was significantly different across MA functions (p ≤ 0.039) and occurred at different knee flexion angles. The shape of the VL force-angle relationship also differed significantly (p < 0.01) across MA functions. According to the maximum force generated by VL based on its literature-derived PSCA, only the VL force-angle relationships estimated using geometric imaging-based MA functions are feasible across the knee angles studied here. We therefore recommend that an average of these MA functions is calculated to estimate quadriceps muscle forces if subject-specific MAs cannot be determined.
Uncertainties of synchrotron microCT-based digital volume correlation bone strain measurements under simulated deformation.
Venue: Journal of biomechanics

Digital Volume Correlation (DVC) is used to measure internal displacements and strains in bone. Recent studies have shown that Synchrotron radiation micro-computed tomography (SR-microCT) can improve the accuracy and precision of DVC. However, only zero-strain or virtually-moved test have been used to quantify the DVC uncertainties, leading to potential underestimation of the measurement errors. In this study, for the first time, the uncertainties of a global DVC approach have been evaluated on repeated SR-microCT scans of bovine cortical bone (voxel size: 1.6 μm), which were virtually deformed for different magnitudes and along different directions. The results showed that systematic and random errors of the normal strain components along the deformation direction were higher than the errors along unstrained directions. The systematic percentage errors were smaller for larger virtual deformations. The random percentage error was in the order of 10% of the virtual deformation. However, higher errors were localized at the boundary of the volumes of interest, perpendicular to the deformation direction. When only the central region of the samples was considered (100 µm layers removed from the borders where the deformation was applied), the errors in the direction of virtual deformation were comparable to the errors in the unstrained directions. In conclusion, the method presented to estimate the uncertainties of DVC is suitable for testing anisotropic specimens as cortical bone. The good agreement between the uncertainties in measurements of strain components obtained with this approach and with the simpler zero-strain-test suggests that the latter is adequate in the tested deformation scenarios.
Sensor-to-body calibration procedure for clinical motion analysis of lower limb using magnetic and inertial measurement units.
Venue: Journal of biomechanics

Magnetic and Inertial measurement units (MIMUs) have become exceedingly popular for ambulatory human motion analysis during the past two decades. However, measuring anatomically meaningful segment and joint kinematics requires virtual alignment of the MIMU frame with the anatomical frame of its corresponding segment. Therefore, this paper presents a simple calibration procedure, based on MIMU readouts, to align the inertial frame of the MIMU with the anatomical frames, as recommended by ISB. The proposed calibration includes five seconds of quiet standing in a neutral posture followed by ten consecutive hip flexions/extensions. This procedure will independently calibrate MIMUs attached to the pelvis, thigh, shank, and foot. The accuracy and repeatability of the calibration procedure and the 3D joint angle estimation were validated against the gold standard motion capture system by an experimental study with ten able-bodied participants. The procedure showed high test-retest repeatability in aligning the MIMU frame with its corresponding anatomical frame, i.e., the helical angle between the MIMU and anatomical frames did not significantly differ between the test and retest sessions (except for thigh MIMU). Compared to previously introduced procedures, this procedure attained the highest inter-participant repeatability (inter-participant coefficient of variations of the helical angle: 20.5-42.2%). Further, the proposed calibration would reduce the offset errors of the 3D joint angle estimation (up to 12.8 degrees on average) compared to joint angle estimation without calibration (up to 26.3 degrees on average). The proposed calibration enables MIMU to measure clinically meaningful gait kinematics.
The between-day reliability of peroneus longus EMG during walking.
Venue: Journal of biomechanics

The peroneus longus (PL) is a rearfoot evertor, important in frontal plane foot motion. Studying PL function has been limited by previous electromyography (EMG) studies reporting poor between-day reliability. Due to its close proximity to adjacent muscles, EMG measures of PL may be susceptible to crosstalk, thus correct electrode placement is vital. The aim of this study was to use ultrasound to aid placement of small surface EMG electrodes and determine the between-day reliability of PL EMG in healthy participants' walking. Ten participants walked barefoot and shod at a controlled, self-selected speed. Six trials per condition, per session, were recorded over two days (mean (SD): 5 (3) days apart). The muscle belly was located using ultrasound. EMG was recorded with surface electrodes (Trigno™ Mini, Delsys, Inc.) at 2000 Hz. Amplitude was normalized to the peak per gait cycle and time normalized to the gait cycle. Reliability of discrete variables were primarily assessed with the standard error of measurement (SEM), plus the coefficient of multiple correlation (CMC), the coefficient of variation (CV) and the intra-class correlation coefficient (ICC). The pattern of the EMG profile was consistent. The SEM of peak amplitude was 4% (3-8%) and 3% (2-5%) for barefoot and shod respectively. For timing of the peak the SEM was 2% (1-3%) and 1% (1-2%) for barefoot and shod respectively. Low SEM of discrete variables suggests good reliability of PL EMG during walking supporting the future use of this protocol. Therefore activation of PL can be confidently studied in repeated-measures study designs.
Maximising forward somersault rotation in springboard diving.
Venue: Journal of biomechanics

Performance in the flight phase of springboard diving is limited by the amounts of linear and angular momentum generated during the takeoff phase. A planar 8-segment torque-driven simulation model combined with a springboard model was used to investigate optimum takeoff technique for maximising rotation in forward dives from the one metre springboard. Optimisations were run by varying the torque activation parameters to maximise forward rotation potential (angular momentum × flight time) while allowing for movement constraints, anatomical constraints, and execution variability. With a constraint to ensure realistic board clearance and anatomical constraints to prevent joint hyperextension, the optimised simulation produced 24% more rotation potential than a simulation matching a 2½ somersault piked dive. When 2 ms perturbations to the torque onset timings were included for the ankle, knee and hip torques within the optimisation process, the model was only able to produce 87% of the rotation potential achieved in the matching simulation. This implies that a pre-planned technique cannot produce a sufficiently good takeoff and that adjustments must be made during takeoff. When the initial onset timings of the torque generators were unperturbed and 10 ms perturbations were introduced into the torque onset timings in the board recoil phase, the optimisation produced 8% more rotation potential than the matching simulation. The optimised simulation had more hip flexion and less shoulder extension at takeoff than the matching simulation. This study illustrates the difficulty of including movement variability within performance optimisation when the movement duration is sufficiently long to allow feedback corrections.
Influence of experimental conditions on visco-hyperelastic properties of skeletal muscle tissue using a Box-Behnken design.
Venue: Journal of biomechanics

The Mechanical characterization of skeletal muscles is strongly dependent on numerous experimental design factors. Nevertheless, significant knowledge gaps remain on the characterization of muscle mechanics and a large number of experiments should be implemented to test the influence of a large number of factors. In this study, we propose a design of experiment method (DOE) to study the parameter sensitivity while minimizing the number of tests. A Box-Behnken design was then implemented to study the influence of strain rate, preconditioning and preloading conditions on visco-hyperelastic mechanical parameters of two rat forearm muscles. The results show that the strain rate affects the visco-hyperelastic parameters for both muscles. These results are consistent with previous work demonstrating that stiffness and viscoelastic contributions increase with strain rate. Thus, DOE has been shown to be a valid method to determine the effect of the experimental conditions on the mechanical behaviour of biological tissues such as skeletal muscle. This method considerably reduces the number of experiments. Indeed, the presented study using 3 parameters at 3 levels would have required at least 54 tests per muscle against 14 for the proposed DOE method.
Non-identifiability of parameters for a class of shear-thinning rheological models, with implications for haematological fluid dynamics.
Venue: Journal of biomechanics

Choosing a suitable model and determining its associated parameters from fitting to experimental data is fundamental for many problems in biomechanics. Models of shear-thinning complex fluids, dating from the work of Bird, Carreau, Cross and Yasuda, have been applied in highly-cited computational studies of hemodynamics for several decades. In this manuscript we revisit these models, first to highlight a degree of uncertainty in the naming conventions in the literature, but more importantly to address the problem of inferring model parameters by fitting to rheology experiments. By refitting published data, and also by simulation, we find large, flat regions in likelihood surfaces that yield families of parameter sets which fit the data equally well. Despite having almost indistinguishable fits to experimental data these varying parameter sets can predict very different flow profiles, and as such these parameters cannot be used to draw conclusions about physical properties of the fluids, such as zero-shear viscosity or relaxation time of the fluid, or indeed flow behaviours. We verify that these features are not a consequence of the experimental data sets through simulations; by sampling points from the rheological models and adding a small amount of noise we create a synthetic data set which reveals that the problem of parameter identifiability is intrinsic to these models.
Multiscale mechanisms of tendon fatigue damage progression and severity are strain and cycle dependent.
Venue: Journal of biomechanics

Tendinopathies are common chronic injuries that occur when damage accumulation caused by sub-rupture fatigue loading outpaces repair. Studies have linked fatigue loading with various mechanical, structural, and biological changes associated with pathology. However, the multiscale progression of damage accumulation with respect to area, severity and the distinct contributions of strain level and number of cycles has not been fully elucidated. The objective of this study was to investigate multiscale mechanisms underlying fatigue damage accumulation and their effect on the cellular environment. Using an in situ model in rat tail tendon (RTT), fatigue loading was applied at various strains and cycle numbers to induce fatigue damage. Pre- and post- fatigue diagnostic mechanical testing, second harmonic generation (SHG) imaging, and transmission electron microscope (TEM) imaging were used to investigate extracellular and cellular damage modes at multiple scales. Fatigue loading at strains at or below 1.0% resulted in no significant changes in SHG damage area or severity and no changes in collagen fibril or cell morphology compared with controls. Fatigue loading at strains above 1.5% resulted in greater mechanical changes correlated with increased damage area measured by SHG and collagenous damage observed by TEM. Increased cycles at high strain further altered mechanical properties, increased structural damage severity (but not area), and altered TEM collagen rupture patterns. Cell morphology was similarly progressively affected with increased strain and cycle number. These damage mechanisms that may trigger degenerative changes characteristic of tendinopathy could be targeted as a part of prevention or therapy.