nils-herrmann/clean_scirepeval huggingface dataset:
Validation of radiostereometric analysis in six degrees of freedom for use with reverse total shoulder arthroplasty.
Venue: Journal of biomechanics

A phantom study was conducted to determine bias in motion and bias at zero motion of radiostereometric analysis (RSA) for evaluating implant relative displacement in reverse total shoulder arthroplasty (RTSA). A Sawbones shoulder phantom was fitted with a RTSA implant set and 13 tantalum markers. The model was fixed to a manual micrometer, providing controlled movements though fifteen known increments in translation and twelve increments in rotation (0.02-5.00 mm and 0.1-6.0°), along each translation and rotation axis. Movement between the glenoid and humerus was assessed using beads vs. beads (B/B), model vs. beads (M/B), and model vs. model (M/M) measurement methods in a model-based RSA environment. Bias in motion and bias at zero motion were defined as the difference between measured and accepted reference values, and the difference between double examinations with a theoretical displacement of zero, respectively. Bias in motion ranged from 0.054 ± 0.010 to 0.129 ± 0.014 mm and 0.076 ± 0.025 to 0.126 ± 0.025° (B/B), 0.023 ± 0.009 to 0.126 ± 0.016 mm and 0.111 ± 0.033 to 0.794 ± 0.251° (M/B), and 0.029 ± 0.010 to 0.135 ± 0.030 mm and 0.243 ± 0.088 to 0.384 ± 0.153° (M/M). Bias at zero motion ranged from 0.120 to 0.156 mm and 0.075 to 0.206° (B/B), 0.074 to 0.149 mm and 0.067 to 1.953° (M/B), and 0.069 to 0.259 mm and 0.284 to 1.273° (M/M). This is the first RSA for RTSA study, with results comparable to those validating the use of RSA for hip and knee arthroplasties (accepted as 0.05-0.50 mm and 0.15-1.15°), justifying the potential use of RSA as a tool for measuring implant displacement in the shoulder.
Micro-CT and micro-FE analysis of pedicle screw fixation under different loading conditions.
Venue: Journal of biomechanics

Anchorage of pedicle screw instrumentation in the elderly spine with poor bone quality remains challenging. In this study, micro finite element (µFE) models were used to assess the specific influence of screw design and the relative contribution of local bone density to fixation mechanics. These were created from micro computer tomography (µCT) scans of vertebras implanted with two types of pedicle screws, including a full region-or-interest of 10 mm radius around each screw, as well as submodels for the pedicle and inner trabecular bone of the vertebral body. The local bone volume fraction (BV/TV) calculated from the µCT scans around different regions of the screw (pedicle, inner trabecular region of the vertebral body) were then related to the predicted stiffness in simulated pull-out tests as well as to the experimental pull-out and torsional fixation properties mechanically measured on the corresponding specimens. Results show that predicted stiffness correlated excellently with experimental pull-out strength (R2 > 0.92, p < .043), better than regional BV/TV alone (R2 = 0.79, p = .003). They also show that correlations between fixation properties and BV/TV were increased when accounting only for the pedicle zone (R2 = 0.66-0.94, p ≤  .032), but with weaker correlations for torsional loads (R2 < 0.10). Our analyses highlight the role of local density in the pedicle zone on the fixation stiffness and strength of pedicle screws when pull-out loads are involved, but that local apparent bone density alone may not be sufficient to explain resistance in torsion.
Which data should be tracked in forward-dynamic optimisation to best predict muscle forces in a pathological co-contraction case?
Venue: Journal of biomechanics

The choice of the cost-function for predicting muscle forces during a movement remains a challenge, especially in patients with neuromuscular disorders. Forward dynamics-based optimisations mainly track joint kinematics or torques, combined with a least-excitation criterion. Tracking marker trajectories and/or electromyography (EMG) has rarely been proposed. Our objective was to determine the best tracking objective-function to accurately predict the upper-limb muscle forces. A musculoskeletal model was created and EMG was simulated to obtain a reference movement - a shoulder abduction. A Gaussian noise (mean = 0; standard deviation = 15%) was added to the simulated EMG. Another noise - corresponding to the actual soft tissue artefacts (STA) of experimental shoulder abduction movements - was added to the trajectories of the markers placed on the model. Muscle forces were estimated from these noisy data, using forward dynamics assisted by six non-linear least-squared objective-functions. These functions involved the tracking of marker trajectories, joint angles or torques, with and without EMG-tracking. All six approaches used the same musculoskeletal model and were solved using a direct multiple shooting algorithm. Finally, the predicted joint angles, muscle forces and activations were compared to the reference values, using root-mean-square errors (RMSe) and biases. The force RMSe of the approach tracking both marker trajectories and EMG (18.45 ± 12.60 N) was almost five times lower than the one of the approach tracking only joint angles (82.37 ± 66.26 N) or torques (85.10 ± 116.40 N). Therefore, using EMG as a complementary tracking-data in forward dynamics seems to be promising for the estimation of muscle forces.
Dynamic structure of lower limb joint angles during walking post-stroke.
Venue: Journal of biomechanics

Variability in joint kinematics is necessary for adaptability and response to everyday perturbations; however, intrinsic neuromotor changes secondary to stroke often cause abnormal movement patterns. How these abnormal movement patterns relate to joint kinematic variability and its influence on post-stroke walking impairments is not well understood.



The purpose of this study was to evaluate the movement variability at the individual joint level in the paretic and non-paretic limbs of individuals post-stroke.



Seven individuals with hemiparesis post-stroke walked on a treadmill for two minutes at their self-selected speed and the average speed of the six-minute walk test while kinematics were recorded using motion-capture. Variability in hip, knee, and ankle flexion/extension angles during walking were quantified with the Lyapunov exponent (LyE). Interlimb differences were evaluated.



The paretic side LyE was higher than the non-paretic side at both self-selected speed (Hip: 50%; Knee: 74%), and the average speed of the 6-min walk test (Hip: 15%;  93%).



Differences in joint kinematic variability between limbs of persons post-stroke supports further study of the source of non-paretic limb deviations as well as the clinical implications of joint kinematic variability in persons post-stroke. The development of bilaterally-targeted post-stroke gait interventions to address variability in both limbs may promote improved outcomes.
Reduced center of pressure modulation elicits foot placement adjustments, but no additional trunk motion during anteroposterior-perturbed walking.
Venue: Journal of biomechanics

Understanding balance during human gait is complicated by the abundance of recovery options. Among all possible recovery options, three main strategies are often considered for human balance control, being the ankle, hip, and foot placement strategies. All can be addressed when balance is threatened during walking, but their relative importance remains uncertain. We have previously shown that healthy human subjects did not significantly adjust their foot placement relative to the body's center of mass (COM) in the first recovery step following anteroposterior pelvis perturbations, as compared to unperturbed walking. An ankle strategy could have contributed to the recovery instead. Here the goal is to further elucidate balance strategy preferences by investigating the stepping and hip strategies following these anteroposterior perturbations, but with an ankle strategy made ineffective. This was achieved by physically blocking each ankle and minimizing the support area of each foot through a pair of modified ankle-foot orthoses. These "pin-shoes" enabled stilt-like walking and ensured that foot placement adjustment was the only way to modulate the center of pressure location, comparable to "footless" inverted pendulum models of walking. Despite the pin-shoes, subjects did not additionally address a hip strategy compared to normal walking, but relied on foot placement adjustments instead. The observed foot placement adjustments were furthermore in line with concepts derived from a linear inverted pendulum model of walking. These results suggest low hip strategy priority when a foot placement strategy is available, while the latter can be predicted with concepts derived from a simple walking model.
Validation of lumbar spine loading from a musculoskeletal model including the lower limbs and lumbar spine.
Venue: Journal of biomechanics

Low back mechanics are important to quantify to study injury, pain and disability. As in vivo forces are difficult to measure directly, modeling approaches are commonly used to estimate these forces. Validation of model estimates is critical to gain confidence in modeling results across populations of interest, such as people with lower-limb amputation. Motion capture, ground reaction force and electromyographic data were collected from ten participants without an amputation (five male/five female) and five participants with a unilateral transtibial amputation (four male/one female) during trunk-pelvis range of motion trials in flexion/extension, lateral bending and axial rotation. A musculoskeletal model with a detailed lumbar spine and the legs including 294 muscles was used to predict L4-L5 loading and muscle activations using static optimization. Model estimates of L4-L5 intervertebral joint loading were compared to measured intradiscal pressures from the literature and muscle activations were compared to electromyographic signals. Model loading estimates were only significantly different from experimental measurements during trunk extension for males without an amputation and for people with an amputation, which may suggest a greater portion of L4-L5 axial load transfer through the facet joints, as facet loads are not captured by intradiscal pressure transducers. Pressure estimates between the model and previous work were not significantly different for flexion, lateral bending or axial rotation. Timing of model-estimated muscle activations compared well with electromyographic activity of the lumbar paraspinals and upper erector spinae. Validated estimates of low back loading can increase the applicability of musculoskeletal models to clinical diagnosis and treatment.
Image-based computational assessment of vascular wall mechanics and hemodynamics in pulmonary arterial hypertension patients.
Venue: Journal of biomechanics

Pulmonary arterial hypertension (PAH) is a disease characterized by an elevated pulmonary arterial (PA) pressure. While several computational hemodynamic models of the pulmonary vasculature have been developed to understand PAH, they are lacking in some aspects, such as the vessel wall deformation and its lack of calibration against measurements in humans. Here, we describe a computational modeling framework that addresses these limitations. Specifically, computational models describing the coupling of hemodynamics and vessel wall mechanics in the pulmonary vasculature of a PAH patient and a normal subject were developed. Model parameters, consisting of linearized stiffness E of the large vessels and Windkessel parameters for each outflow branch, were calibrated against in vivo measurements of pressure, flow and vessel wall deformation obtained, respectively, from right-heart catheterization, phase-contrast and cine magnetic resonance images. Calibrated stiffness E of the proximal PA was 2.0 and 0.5 MPa for the PAH and normal models, respectively. Calibrated total compliance CT and resistance RT of the distal vessels were, respectively, 0.32 ml/mmHg and 11.3 mmHg∗min/l for the PAH model, and 2.93 ml/mmHg and 2.6 mmHg∗min/l for the normal model. These results were consistent with previous findings that the pulmonary vasculature is stiffer with more constricted distal vessels in PAH patients. Individual effects on PA pressure due to remodeling of the distal and proximal compartments of the pulmonary vasculature were also investigated in a sensitivity analysis. The analysis suggests that the remodeling of distal vasculature contributes more to the increase in PA pressure than the remodeling of proximal vasculature.
Effects of tear size and location on predictions of supraspinatus tear propagation.
Venue: Journal of biomechanics

Rotator cuff tears remain a significant clinical problem with a high incidence rate and severe clinical burden. Previous computational models developed to study rotator cuff tears have not modeled tissue damage and tear propagation. The objective of this study was to predict tear propagation for various combinations of tear size and location using an experimentally validated finite element model of supraspinatus tendon. It was hypothesized that larger rotator cuff tears propagate at lower loads than smaller tears, and that posterior tears require higher loads to propagate than anterior tears. Using a previously validated computational model of supraspinatus tendon, tears of size 0.5-1.5 cm were introduced to the tendon geometry in the anterior, middle, and posterior tendon thirds. Cohesive elements were assigned subject-specific failure properties and used to model tissue damage and tear propagation. A displacement of 5 mm was applied to the medial tendon edge to induce tear propagation. Model outputs included critical load required to propagate the tear, and principal stress and maximum principal strain at the anterior and posterior tear tips. For all tear sizes, posterior tears required the highest loads to propagate (247-567 N). Anterior tears generally required the least load to propagate (171-280 N). Stress and strain were larger on the articular side (maximum 33.9% articular strain vs 27.8% bursal strain). Overall, larger tears located in the anterior supraspinatus tendon that interrupt the rotator cable are most at risk for tear propagation, and should be carefully followed by clinicians when considering treatment options.
Wrist tendon moment arms: Quantification by imaging and experimental techniques
Venue: Journal of biomechanics

Subject-specific musculoskeletal models require accurate values of muscle moment arms. The aim of this study was to compare moment arms of wrist tendons obtained from non-invasive magnetic resonance imaging (MRI) to those obtained from an in vitro experimental approach. MRI was performed on ten upper limb cadaveric specimens to obtain the centrelines for the flexor carpi radialis (FCR), flexor carpi ulnaris (FCU), extensor carpi radialis longus (ECRL), extensor carpi radialis brevis (ECRB), extensor carpi ulnaris (ECU), and abductor pollicis longus (APL) tendons. From these, the anatomical moment arms about each of the flexion-extension (FE) and radioulnar deviation (RUD) axes of the wrist were calculated. Specimens were mounted on a physiologic wrist simulator to obtain functional measurements of the moment arms using the tendon excursion method. No differences were observed between anatomical and functional values of the FE and RUD moment arms of FCR, ECRL and ECRB, and the RUD moment arm of ECU (p > .075). Scaling the anatomical moment arms relative to ECRB in FE and ECU in RUD reduced differences in the FE moment arm of FCU and the RUD moment arm of APL to less than 15% (p > .139). However, differences persisted in moment arms of FCU in RUD, and ECU and APL in FE (p < .008). This study shows that while measurements of moment arms of wrist tendons using imaging do not always conform to values obtained using in vitro experimental approaches, a stricter protocol could result in the acquisition of subject-specific moment arms to personalise musculoskeletal models.
Effects of walking speed on the step-by-step control of step width.
Venue: Journal of biomechanics

Young, healthy adults walking at typical preferred speeds use step-by-step adjustments of step width to appropriately redirect their center of mass motion and ensure mediolateral stability. However, it is presently unclear whether this control strategy is retained when walking at the slower speeds preferred by many clinical populations. We investigated whether the typical stabilization strategy is influenced by walking speed. Twelve young, neurologically intact participants walked on a treadmill at a range of prescribed speeds (0.2-1.2 m/s). The mediolateral stabilization strategy was quantified as the proportion of step width variance predicted by the mechanical state of the pelvis throughout a step (calculated as R2 magnitude from a multiple linear regression). Our ability to accurately predict the upcoming step width increased over the course of a step. The strength of the relationship between step width and pelvis mechanics at the start of a step was reduced at slower speeds. However, these speed-dependent differences largely disappeared by the end of a step, other than at the slowest walking speed (0.2 m/s). These results suggest that mechanics-dependent adjustments in step width are a consistent component of healthy gait across speeds and contexts. However, slower walking speeds may ease this control by allowing mediolateral repositioning of the swing leg to occur later in a step, thus encouraging slower walking among clinical populations with limited sensorimotor control.
Muscle activity during maximal isometric forearm rotation using a power grip.
Venue: Journal of biomechanics

This study aimed to provide quantitative activation data for muscles of the forearm during pronation and supination while using a power grip. Electromyographic data was collected from 15 forearm muscles in 11 subjects while they performed maximal isometric pronating and supinating efforts in nine positions of forearm rotation. Biceps brachii was the only muscle with substantial activation in only one effort direction. It was significantly more active when supinating (µ = 52.1%, SD = 17.5%) than pronating (µ = 5.1%, SD = 4.8%, p < .001). All other muscles showed considerable muscle activity during both pronation and supination. Brachioradialis, flexor carpi radialis, palmaris longus, pronator quadratus and pronator teres were significantly more active when pronating the forearm. Abductor pollicis longus and biceps brachii were significantly more active when supinating. This data highlights the importance of including muscles additional to the primary forearm rotators in a biomechanical analysis of forearm rotation. Doing so will further our understanding of forearm function and lead to the improved treatment of forearm fractures, trauma-induced muscle dysfunction and joint replacements.
Interpersonal coordination analysis of tennis players from different levels during official matches.
Venue: Journal of biomechanics

The purpose of this study was to assess the interpersonal coordination during official male tennis matches in players of different skill levels. Players' trajectories of three levels (Juvenile, ATP-Future, ATP-250) were obtained using video-based tracking system. A vector coding technique was applied to obtain players' interpersonal coordination in four coordination patterns: "anti-phase", "in-phase", "serving player phase" and "returning player phase". These patterns allowed identification of the nature of the coupling and lead/lag relations between players. In all categories, players presented higher degree of "anti-phase" and "serving player phase" (when only the serving player is moving on the court or his opponent is moving with a time lag) coordination. Young players spent more time in "serving player phase" during lateral displacements than professional players. On the other hand, professional players spent more time in "returning player phase" (when only the returning player is moving on the court or his opponent is moving with a time lag) during antero posterior displacements, than young players. Interpersonal coordination did not change from the first to the second set of the match, showing that tennis players maintain their displacement characteristics and strategy, independently of proficiency level. The vector coding technique allowed to identify new coordination patterns in tennis, providing additional information about tennis dynamics and how players from different categories and proficiency levels behave during the matches.
Effect of high-intensity perturbations during core-specific sensorimotor exercises on trunk muscle activation.
Venue: Journal of biomechanics

Core-specific sensorimotor exercises are proven to enhance neuromuscular activity of the trunk. However, the influence of high-intensity perturbations on training efficiency is unclear within this context. Sixteen participants (29±2yrs; 175±8cm; 69±13kg) were prepared with a 12-lead bilateral trunk EMG. Warm-up on a dynamometer was followed by maximum voluntary isometric trunk (flex/ext) contraction (MVC). Next, participants performed four conditions for a one-legged stance with hip abduction on a stable surface (HA) repeated randomly on an unstable surface (HAP), on a stable surface with perturbation (HA+P), and on an unstable surface with perturbation (HAP+P). Afterwards, bird dog (BD) was performed under the same conditions (BD, BDP, BD+P, BDP+P). A foam pad under the foot (HA) or the knee (BD) was used as an unstable surface. Exercises were conducted on a moveable platform. Perturbations (ACC 50m/sec2;100ms duration;10rep.) were randomly applied in the anterior-posterior direction. The root mean square (RMS) normalized to MVC (%) was calculated (whole movement cycle). Muscles were grouped into ventral right and left (VR;VL), and dorsal right and left (DR;DL). Ventral-Dorsal and right-left ratios were calculated (two way repeated-measures ANOVA;α=0.05). Amplitudes of all muscle groups in bird dog were higher compared to hip abduction (p≤0.0001; ). EMG-RMS showed significant differences (p<0.001) between conditions and muscle groups per exercise. Interaction effects were only significant for HA (p=0.02). No significant differences were present in EMG ratios (p>0.05). Additional high-intensity perturbations during core-specific sensorimotor exercises lead to increased neuromuscular activity and therefore higher exercise intensities. However, the beneficial effects on trunk function remain unclear. Nevertheless, BD is more suitable to address trunk muscles.