nils-herrmann/clean_scirepeval huggingface dataset:
Increased force following muscle stretching and simultaneous fibre shortening: Residual force enhancement or force depression - That is the question?
Venue: Journal of biomechanics

Residual force enhancement (rFE) describes the increase in isometric force following muscle stretching compared to the corresponding isometric force. Even though rFE is consistently observed in isolated muscle preparations, it is not always observed in human skeletal muscle. This inconsistency might be associated with disociations between length changes in muscle tendon units (MTUs) and fibres. This prompted the question if there is rFE for conditions where the MTU is stretched while fibres shorten. Rabbit tibialis anterior (TA) MTUs (n = 4) were stretched and the isometric forces following stretching were compared to corresponding forces from isometric reference contractions. Unique combinations of stretch speed and activation were used to create conditions of continuous fibre shortening during MTU stretch. Mean force was increased (18 ± 2%) following MTU stretching compared to the isometric reference forces. Without fibre length measurements, this result would be referred to as rFE. However, fibre shortening in the reference contractions was always greater than for the eccentric stretch contractions, suggesting that the observed increase in force might be caused by less residual force depression (rFD) in the stretch tests compared to the reference contractions. However, the work performed by fibre shortening was similar between the reference and the MTU stretch contractions, suggesting that rFD was similar for both experimental conditions. Therefore, we conclude that we observed rFE in the absence of contractile element stretching. However, a lack of knowledge of the molecular mechanisms that distinguish rFE from rFD prevents an unequivocal pronouncement of what caused the enhanced forces after active muscle stretching.
A non-invasive technique for evaluating carpal tunnel pressure with ultrasound vibro-elastography for patients with carpal tunnel syndrome: A pilot clinical study.
Venue: Journal of biomechanics

Carpal tunnel syndrome (CTS) is a disorder that affects the median nerve at the wrist sufficient to cause impairment of nerve function. Elevated carpal tunnel pressure (CTP) leads to median nerve pathology, sensory, and motor changes in CTS patient. The techniques to quantify CTP used in clinic are invasive. This study aimed to investigate the feasibility of a noninvasive ultrasound vibro-elastography (UVE) to predict CTP in CTS patients and healthy individuals. The magnitudes of shear wave speed ratio (rSWS) of the 10 CTS patients (10 hands) and 6 healthy individuals (12 hands), and 10 cadaveric hands were compared using UVE. The ratios of intra to extra-carpal tunnel SWS in CTS patients was significantly higher than those in the healthy individuals (p = 0.0008) and cadaveric hands (p = 0.0015) with 500-g tendon tension. We estimated the CTP in the carpal tunnel using the mean rSWS of each group obtained from the present study and the linear approximation obtain from cadaveric hands data with 500-g tendon tension (y = 0.0036x + 1.1413). These results indicated that the elevated pressure applied to the 3rd flexor digitorum superficialis tendon in the carpal tunnel of CTS patients resulted in faster shear wave propagation. These results show that UVE was useful to indirectly estimate the CTP by measuring the rSWS; thus, they are potentially useful for the early diagnosis and assessment of CTS.
Effect of walking on in vivo tibiofemoral cartilage strain in ACL-deficient versus intact knees.
Venue: Journal of biomechanics

Anterior cruciate ligament (ACL) rupture alters knee kinematics and contributes to premature development of osteoarthritis. However, there is limited data regarding the in vivo biomechanical response of tibiofemoral cartilage to activities of daily living (ADLs) in ACL-deficient knees. In this study, eight otherwise healthy participants with chronic unilateral ACL deficiency completed a stress test to assess the effect of 20 min of level treadmill walking at a speed of 2.5 mph on tibiofemoral cartilage in their ACL-deficient and contralateral ACL-intact knees. Three-dimensional surface models developed from pre- and post-activity magnetic resonance (MR) images of the injured and uninjured knees were used to determine compressive strain across multiple regions of tibiofemoral cartilage (medial and lateral tibial plateaus, medial and lateral femoral condyles, medial aspect of femoral condyle adjacent to intercondylar notch of the femur). In the ACL-deficient knees, we observed significantly increased cartilage strain in the region of the medial femoral condyle adjacent to the intercondylar notch (6% in deficient vs. 2% in contralateral, p = 0.01) as well as across the medial and lateral tibial plateaus (4% vs. 3%, p = 0.01) relative to the contralateral ACL-intact knees. Increased compressive strain at the medial intercondylar notch and tibial plateau suggests alterations in mechanical loading or the response to load in these regions, presumably related to altered knee kinematics. These changes may disrupt cartilage homeostasis and contribute to subsequent development of osteoarthritis.
Clinical feasibility of fracture healing assessment through continuous monitoring of implant load.
Venue: Journal of biomechanics

Current fracture fixation follow-up is based on subjective radiological and clinical examination. Efforts to objectify the procedure have been undertaken since decades. Assessment of implant load as an indirect predictor of callus maturity has so far failed to enter clinical routine due to limited practicability, technical obstacles and its snap-shot nature. We recently introduced the concept of continuous implant load monitoring to aid in diagnosing fracture healing progression. This study aimed at investigating the feasibility of the system in a clinical context. Ten patients treated with Taylor-Spatial-Frame external fixators following pathological tibia fractures were equipped with a Fracture Monitor device attached to a fixator-strut and were monitored until hardware removal. Two patients were excluded due to technical issues. Implant load and fracture activity was continuously and autonomously measured for 139 ± 89 days (mean ± SD). Data was wirelessly collected with consumer smartphones. Relative implant load initially rose for 34.1 ± 22.2 days and finally declined to a level of 45.0 ± 33.8% of the maximum implant load. In five patients the load dropped below 50% of the maximum load. These patients underwent hardware removal according to the clinical assessment. In three patients, whose external fixators were exchanged to internal fixation at the end of the study, implant load did not drop below the 50% margin. The continuous measurement principle allows resolving implant load progression and appears indicative for the bone healing status. Data can be acquired in a homecare setting and is believed to provide valuable information to support timely healing assessment and enable patient specific after-care.
The validity and reliability of the Achilles tendon moment arm assessed with dual-energy X-ray absorptiometry, relative to MRI and ultrasound assessments.
Venue: Journal of biomechanics

Dual-energy X-ray absorptiometry (DXA) in single energy mode has been shown to permit the visualisation of bone and soft tissue, such as the patellar tendon through two-dimensional sagittal imaging. However, there is no validated DXA-based measurement of the Achilles tendon moment arm (dAT). The aims of this study were: 1) to compare in vivo DXA derived measurements of the dAT at rest against two previously validated methods: tendon excursion (TE) and magnetic resonance imaging (MRI) at three ankle angles (-5°, 0° and +10°). 2) analyse the intra-day reliability of the DXA method at all ankle angles and compare between methods. Twelve healthy adults (mean ± SD: 31.4 ± 9.5 years; 174.0 ± 9.5 cm; 76.2 ± 16.6 kg) participated in this study, involving test-retest DXA scans, ultrasound scans and one MRI scan. The dAT was defined as the distance from the centre of the calcaneal-tibial joint axis to the Achilles tendon (AT) muscle-tendon line of action. DXA derived dAT measures were significantly greater than MRI measurements (19.7-24.9%) and were 45.2% significantly larger than the TE method. The test-retest reliability of the DXA technique at 0° was high [CV = 1.38%; ICC = 0.96] and despite the consistently larger dAT lengths obtained using DXA, MRI and DEXA data were strongly correlated (r = 0.878, p < 0.001). In conclusion, the DXA technique allowed for highly reproducible in vivo dAT measurement at rest, which has implications for the calculation of AT forces in vivo and the ability to predict the measurement from one tool to the other, thereby providing a novel basis to contrast existing and future studies.
Causal relationships between immediate pre-impact kinematics and post-impact kinetics during drop landing using a simple three dimensional multibody model.
Venue: Journal of biomechanics

This study aimed to validate a simple dynamic model of single-leg drop-landing to develop a methodological foundation for investigating mechanistic causes of anterior cruciate ligament (ACL) injury and to explore mechanical associations between knee valgus torque and landing kinematics that are considered clinically as a high-risk landing posture for the injury. A triple-inverted-pendulum model in three-dimensional space, composed of rigid-links of head-arms-trunk (HAT), thigh and shank, was employed. We derived causal relationships that can predict post-impact kinetics, including impact ground reaction forces (GRFs) and corresponding knee joint torques from a given body-kinematics immediately before impact, based on an assumption of a completely inelastic collision between a landing foot (the distal end-point of the shank in the model) and the ground. The concordance correlation coefficient (CCC) analysis revealed that our model can achieve an acceptable agreement between experimentally measured and model-predicted impact GRFs and corresponding knee joint torques. The 95% one-tailed lower confidence limit of CCC of vertical, mediolateral GRFs and the varus/valgus torque were 0.665>ρc,>ρc,a=0.758 and 0.531>ρc,a=0.508, respectively, for the least acceptable values ρc,a. Using this model, effects of three types of hypothetical pre-impact kinematics with modulated (i) medial/lateral leaning HAT angle, (ii) forward/backward HAT tilt-angle, and (iii) knee flexion/extension angle on the impact GRF and corresponding knee joint torque were evaluated. We showed that the smaller knee flexion and the greater HAT leaning toward the landing-limb-side, the larger the knee valgus torque is generated, as a mechanical consequence between the specific pre-impact kinematics and the knee loading associated with the risk of ACL injury. Further exploration of hypothetical kinematics using the model in the future work might contribute to identifying the risky landing kinematics beyond experimental limitations.
PCA of waveforms and functional PCA: A primer for biomechanics.
Venue: Journal of biomechanics

Principal components analysis (PCA) of waveforms and functional PCA (fPCA) are statistical approaches used to explore patterns of variability in biomechanical curve data, with fPCA being an accepted statistical method grounded within the functional data analysis (FDA) statistical framework. This technical note demonstrates that PCA of waveforms is the most rudimentary form of FDA, and consequently can be rationalised within the FDA framework of statistical processes. Mathematical proofing applied demonstrations of both techniques, and an example of when fPCA may be of greater benefit to control over smoothing of functional principal components is provided using an open access motion sickness dataset. Finally, open access software is provided with this paper as means of priming the biomechanics community for using these methods as a part of future functional data explorations.
Validation of markerless strain-field optical tracking approach for soft tissue mechanical assessment.
Venue: Journal of biomechanics

Strain measurement during tissue deformation is crucial to elucidate relationships between mechanical loading and functional changes in biological tissues. When combined with specified loading conditions, assessment of strain fields can be used to craft models that accurately represent the mechanical behavior of soft tissue. Inhomogeneities in strain fields may be indicative of normal or pathological inhomogeneities in mechanical properties. In this study, we present the validation of a modified Demons registration algorithm for non-contact, marker-less strain measurement of tissue undergoing uniaxial loading. We validate the algorithm on a synthetic dataset composed of artificial deformation fields applied to a speckle image, as well as images of aortic sections of varying perceptual quality. Initial results indicate that Demons outperforms recent Optical Flow and Digital Image Correlation methods in terms of accuracy and robustness to low image quality, with similar runtimes. Demons achieves at least 8% lower maximal deviation from ground truth on 50% biaxial and shear strain applied to aortic images. To illustrate utility, we quantified strain fields of multiple human aortic specimens undergoing uniaxial tensile testing, noting the formation of strain concentrations in areas of rupture. The modified Demons algorithm captured a large range of strains (up to 50%) and provided spatially resolved strain fields that could be useful in the assessment of soft tissue pathologies.
Healthy ankle and hindfoot kinematics during gait: Sex differences, asymmetry and coupled motion revealed through dynamic biplane radiography.
Venue: Journal of biomechanics

The aims of this study were to compare male versus female and dominant versus non-dominant kinematics in the ankle and hindfoot, and to characterize coupled motion between the subtalar and tibiotalar joints during the support phase of gait. Twenty healthy adults walked on a laboratory walkway while synchronized biplane radiographs of the ankle and hindfoot were collected at 100 frames/s. A validated tracking technique was used to measure tibiotalar and subtalar kinematics. Differences between male and female range of motion (ROM) were observed only in tibiotalar (AP and ML) and subtalar (ML) translation (all differences<1 mm and all p < 0.04). Statistical parametric mapping identified differences between kinematics waveforms of males and females in tibiotalar translation (AP and ML) and eversion, and subtalar ML translation. No differences between dominant and non-dominant sides were observed in ROM or kinematics waveforms. The average absolute side-to-side difference in the kinematics waveforms was 4.1° and 1.5 mm or less for all rotations and translations, respectively. Tibiotalar plantarflexion was coupled to subtalar inversion and eversion during the impact and push-off phases of stance (r = 0.90 and r = 0.87, respectively). This data may serve as a guide for evaluating ankle kinematics waveforms, ROM, symmetry, and restoration of healthy coupled motion after surgical intervention or rehabilitation. The observed kinematics differences between males and females may predispose females to higher rates of ankle and knee injury and suggest sex-dependent ankle reconstruction techniques may be beneficial.
Effect of simulated rehabilitation on hip joint loading during single limb squat in patients with hip dysplasia.
Venue: Journal of biomechanics

Rehabilitation for patients with developmental dysplasia of the hip (DDH) addresses modifiable factors in an effort to reduce symptoms and prevent or delay the development of osteoarthritis, yet its effect on joint mechanics remains unknown. Our objective was to establish how rehabilitation (muscle strengthening and movement training), simulated with a musculoskeletal model and probabilistic analyses, alters hip joint reaction forces (JRF) in patients with DDH during a single limb squat. In four patients with DDH, hip abductor strengthening was simulated by increasing the maximum isometric force value between 0 and 32.6% and movement training was simulated by decreasing the hip adduction angle between 0 and 10° relative to baseline. 2,000 Monte Carlo simulations were performed separately to simulate strengthening and movement training, from which 99% confidence bounds and sensitivity factors were calculated. Our results indicated that simulated movement training aimed at decreasing hip adduction had a substantially larger influence on hip JRF than strengthening, as indicated by 99% confidence bounds of the resultant JRF (0.88 ± 0.55 xBW vs. 0.31 ± 0.12 xBW, respectively). Relative to baseline, movement training that resulted in a 10° decrease in hip adduction decreased the resultant JRF by 0.78 ± 0.65 xBW, while strengthening the abductors by 17.6% increased resultant JRF by 0.18 ± 0.06 xBW. To our knowledge, these results are the first to provide evidence pertaining to the effect of rehabilitation on joint mechanics in patients with DDH and can be used to inform more targeted interventions.
Impact of prosthetic mitral valve orientation on the ventricular flow field: Comparison using patient-specific computational fluid dynamics.
Venue: Journal of biomechanics

Significant mitral valve regurgitation creates progressive adverse remodeling of the left ventricle (LV). Replacement of the failing valve with a prosthesis generally improves patient outcomes but leaves the patient with non-physiological intracardiac flow patterns that might contribute to their future risk of thrombus formation and embolism. It has been suggested that the angular orientation of the implanted valve might modify the postoperative distortion of the intraventricular flow field. In this study, we investigated the effect of prosthetic valve orientation on LV flow patterns by using heart geometry from a patient with LV dysfunction and a competent native mitral valve to calculate intracardiac flow fields with computational fluid dynamics (CFD). Results were validated using in vivo 4D Flow MRI. The computed flow fields were compared to calculations following virtual implantation of a mechanical heart valve oriented in four different angles to assess the effect of leaflet position. Flow patterns were visualized in long- and short-axes and quantified with flow component analysis. In comparison to a native valve, valve implantation increased the proportion of the mitral inflow remaining in the basal region and further increased the residual volume in the apical area. Only slight changes due to valve orientation were observed. Using our numerical framework, we demonstrated quantitative changes in left ventricular blood flow due to prosthetic mitral replacement. This framework may be used to improve design of prosthetic heart valves and implantation procedures to minimize the potential for apical flow stasis, and potentially assist personalized treatment planning.
Effects of spinal coupling and marker set on tracking of spine models during running.
Venue: Journal of biomechanics

Despite the wide-spread use of musculoskeletal simulations and its use in estimating spinal loads, much is not known about how to best collect experimental data for modelling purposes. The primary purposes in this study were to determine the effects of tracking of running motion capture data to a model (1) with and without coupling of lumbar spine segments, and (2) with varying combinations of spinal markers. Running trials were collected from 7 participants, with each at three different speeds. The motion data was fit to the Full-Body Lumbar Spine Model (FBLS) with coupling of the lumbar spine enabled (CS) and disabled and therefore rigid (RS) in OpenSim through the Inverse Kinematics tool (IK). Different combinations of markers were chosen as tracking inputs for IK to represent experimental data collection with different marker sets. Root-mean-square (RMS) marker errors of all 13 markers along the spine for each gait cycle were calculated. The CS model resulted in 23.7% lower errors than the RS model (p < 0.001). The marker subset analysis showed that increasing the number of markers in the experimental data collection decreases the error, with the four marker tracking subsets with the highest number of markers tracked having the lowest errors. The location of the marker and timing in the gait cycle did not affect marker error. When spinal mechanics are of interest, the inclusion of a coupled lumbar spine in the model and a larger spinal marker set help better track experimental kinematics when fitting to a model.
Influence of pelvic shape on strain patterns: A computational analysis using finite element mesh morphing techniques.
Venue: Journal of biomechanics

The pelvis functions to transmit upper body loads to the lower limbs and is critical in human locomotion. Semi-automated, finite element (FE) morphing techniques eliminate the need for segmentation and have shown to accelerate the generation of multiple specimen-specific pelvic FE models to enable the study of pelvic mechanical behaviour. The purpose of this research was to produce simulated human pelvic FE models representing android, gynecoid, anthropoid and platypelloid morphologies and to isolate differences in strain patterns due to anatomic shape under physiologic loading. Using five initially generated specimen-specific FE models, each specimen-specific FE model was reconfigured into three different morphologies using FE mesh morphing techniques. Significantly different strains were found comparing the gynecoid (classical female pelvis') to the android ('true male pelvis') models (p = 0.040), with strains twice as high in the superior pubic rami. No significant differences were seen in comparing overall strains between the other pelvic shapes (p = 0.61-0.126). The highest strain regions in all models were found in the supra-acetabular regions, with high strains also found in the regions of the superior pubic rami, the greater sciatic notch and sacral regions about the L5 vertebrae. Quantifying the contributions of shape to strain in the pelvis may increase the understanding of sex and patient-specific differences in fracture risk and motivate the consideration of treatment strategies that account for anatomic pelvic differences.