nils-herrmann/clean_scirepeval huggingface dataset:
Unholey shoes: Experimental considerations when estimating ankle joint complex power during walking and running.
Venue: Journal of biomechanics

For studies that aim to assess biological ankle function, calculating ankle joint complex (AJC) power between the calcaneus and shank is recommended over conventional inverse dynamics estimates between a rigid-body foot and shank. However, when designing a new experiment, it remains unclear whether holes should be cut in footwear to permit motion tracking via skin-mounted markers, or whether marker placement locations should be tightly controlled across conditions. Here we provide data to assist researchers in answering these questions. We performed a gait analysis study of walking (0.8, 1.2, 1.6 m·s-1) and running (2.6, 2.8, 3.0 m·s-1) while subjects (N = 10) wore custom-modified footwear, which allowed markers to be placed either on the shoe, or on the skin via cut-out windows in the shoes. First, we compared foot markers affixed to the skin vs. on the same locations on the shoe. Using statistical non-parametric mapping techniques, we discovered that skin vs. shoe markers had no statistically significant effect on net AJC power estimates throughout stance phase, for all walking and running speeds. Second, we compared calcaneal markers in the nominal shoe configuration vs. markers in a nearby location (∼27 mm below) on the shoe. We observed significant differences when marker placement on the shoe was varied, which may be relevant to repeated-measures study designs. The results suggest that when computing AJC power for walking and running, you may want to put down the scissors (i.e., forego cutting holes in your footwear), and instead pick up a Sharpie® (pen) or use a template, to maintain consistent marker placement across trials and conditions.
Anatomically-based skeletal coordinate systems for use with impact biomechanics data intended for anthropomorphic test device development.
Venue: Journal of biomechanics

Post-mortem human subjects (PMHS) are frequently used to characterize biomechanical response and injury tolerance of humans to various types of loading by means of instrumentation installed directly on the skeleton. Data extracted from such tests are often used to develop and validate anthropomorphic test devices (ATDs), which function as human surrogates in tests for injury assessment. Given that the location and orientation of installed instrumentation differs between subjects, nominally similar measurements made on different PMHS must be transformed to standardized, skeletal-based local coordinate systems (LCS) before appropriate data comparisons can be made. Standardized PMHS LCS that correspond to ATD instrumentation locations and orientations have not previously been published. This paper introduces anatomically-defined PMHS LCS for body regions in which kinematic measurements are made using ATDs. These LCS include the head, sternum, single vertebrae, pelvis, femurs (distal and proximal), and tibiae (distal and proximal) based upon skeletal landmarks extracted from whole body CT scans. The proposed LCS provide a means to standardize the reporting of PMHS data, and facilitate both the comparison of PMHS impact data across institutions and the application of PMHS data to the development and validation of ATDs.
Two-dimensional spatial error distribution of key tensiomyographic parameters.
Venue: Journal of biomechanics

Tensiomyography detects the contraction time (Tc) and amplitude (Dm) of muscle belly thickening during maximal isometric twitch contraction. The assessment of both parameters is highly reliable; however, it seems that their calculation depends on the measurement point. The aim of the study was to determine spatial relative error distribution of Tc and/or Dm within a two-dimensional array of 27 (3 × 9) measurement points in comparison to the reference point (RP) in 12 male participants (22.5 ± 3.1 years). The RPs were determined as follows: in the biceps brachii (BB) at 50% of the humerus length; in the erector spinae (ES) at the height of the iliac crest; in the vastus lateralis (VL), vastus medalis (VM), and rectus femoris (RF) at 30%, 20%, and 50% of femur length above the patella, respectively. The surface area under the 3% relative error in Dm (BB: 4.0; VL: 3.8; VM; 8.2; RF: 6.2; ES: 2.4 cm2) was lower than in Tc (BB: 6.9; VL: 3.8; VM; 4.6; RF: 9.5; ES: 3.7 cm2), yielding merged values (BB: 3.9; VL: 3.7; VM; 4.8; RF: 5.1; ES: 2.4 cm2). Dm show twice as steep relative error rate when moving away from the RP in comparison to Tc, which seems to be less sensitive to spatial sensor positioning.
Residual force enhancement and force depression in human single muscle fibres.
Venue: Journal of biomechanics

Residual force depression (rFD) and residual force enhancement (rFE) are intrinsic contractile properties of muscle. rFD is characterized as a decrease in steady-state isometric force following active shortening compared with a purely isometric contraction at the same muscle length and level of activation. By contrast, isometric force is increased following active lengthening compared to a reference isometric contraction at the same muscle length and level of activation; this is termed rFE. To date, there have been no investigations of rFD and rFE in human muscle fibres, therefore the purpose of this study was to determine whether rFD and rFE occur at the single muscle fibre level in humans. rFD and rFE were investigated in maximally activated single muscle fibres biopsied from the vastus lateralis of healthy adults. To induce rFD, fibres were activated and shortened from an average sarcomere length (SL) of 3.2-2.6 μm. Reference isometric contractions were performed at an average SL of 2.6 μm. To induce rFE, fibres were actively lengthened from an average SL of 2.6-3.2 μm and a reference isometric contraction was performed at an average SL of 3.2 μm. Isometric steady-state force was lower following active shortening (p < 0.05), and higher following active lengthening (p < 0.05), as compared to the reference isometric contractions. We demonstrated rFD and rFE in human single fibres which is consistent with previous animal models. The non-responder phenomenon often reported in rFE studies involving voluntary contractions at the whole human level was not observed at the single fibre level.
On the validity of statistical parametric mapping for nonuniformly and heterogeneously smooth one-dimensional biomechanical data.
Venue: Journal of biomechanics

Nonuniform (non-constant) temporal smoothness can arise in biomechanical processes like impacts, and heterogeneous smoothness (unequal smoothness across observations) can arise in mechanically diverse comparisons such as padded vs. unpadded impacts, where padded dynamics are generally smoother than unpadded dynamics. It has been reported that statistical parametric mapping's (SPM's) probability values can be invalid for such cases. The purpose of this paper was to clarify the scope of validity for SPM analysis of nonuniformly and heterogeneously smooth one-dimensional (1D) data. We simulated a variety of nonuniformly and heterogeneously smooth Gaussian 1D data over a range of smoothness values, and computed Type I error rates across 10,000 simulation iterations for each smoothness type. Results showed that, in all cases, SPM accurately controlled error at the prescribed α=0.05. Moreover, the distribution of false positives was uniform across time, implying that all regions are equally likely to produce false positives, irrespective of local roughness. We nevertheless show that cluster-level inferences (i.e., p values specific to local regions of significance), while never exceeding alpha (by definition), may be over-or-underestimated by approximately 0.01 for the currently simulated scenarios. We conclude that SPM's null hypothesis rejection decisions are valid for both nonuniform and heterogeneous 1D data, but that clusters' p values may be marginally too small/large in rough/smooth regions, respectively. Since cluster-level p values never exceed α, these p value errors are negligible for hypothesis testing purposes. Nevertheless, inter-cluster p value comparisons should be avoided. Implications for statistical power and general results interpretation are discussed.
A framework for biomechanics simulations using four-chamber cardiac models.
Venue: Journal of biomechanics

Computational cardiac models have been extensively used to study different cardiac biomechanics; specifically, finite-element analysis has been one of the tools used to study the internal stresses and strains in the cardiac wall during the cardiac cycle. Cubic-Hermite finite element meshes have been used for simulating cardiac biomechanics due to their convergence characteristics and their ability to capture smooth geometries compactly-fewer elements are needed to build the cardiac geometry-compared to linear tetrahedral meshes. Such meshes have previously been used only with simple ventricular geometries with non-physiological boundary conditions due to challenges associated with creating cubic-Hermite meshes of the complex heart geometry. However, it is critical to accurately capture the different geometric characteristics of the heart and apply physiologically equivalent boundary conditions to replicate the in vivo heart motion. In this work, we created a four-chamber cardiac model utilizing cubic-Hermite elements and simulated a full cardiac cycle by coupling the 3D finite element model with a lumped circulation model. The myocardial fiber-orientations were interpolated within the mesh using the Log-Euclidean method to overcome the singularity associated with interpolation of orthogonal matrices. Physiologically equivalent rigid body constraints were applied to the nodes along the valve plane and the accuracy of the resulting simulations were validated using open source clinical data. We then simulated a complete cardiac cycle of a healthy heart and a heart with acute myocardial infarction. We compared the pumping functionality of the heart for both cases by calculating the ventricular work. We observed a 20% reduction in acute work done by the heart immediately after myocardial infarction. The myocardial wall displacements obtained from the four-chamber model are comparable to actual patient data, without requiring complicated non-physiological boundary conditions usually required in truncated ventricular heart models.
The effects of electromyography-assisted modelling in estimating musculotendon forces during gait in children with cerebral palsy.
Venue: Journal of biomechanics

Neuro-musculoskeletal modelling can provide insight into the aberrant muscle function during walking in those suffering cerebral palsy (CP). However, such modelling employs optimization to estimate muscle activation that may not account for disturbed motor control and muscle weakness in CP. This study evaluated different forms of neuro-musculoskeletal model personalization and optimization to estimate musculotendon forces during gait of nine children with CP (GMFCS I-II) and nine typically developing (TD) children. Data collection included 3D-kinematics, ground reaction forces, and electromyography (EMG) of eight lower limb muscles. Four different optimization methods estimated muscle activation and musculotendon forces of a scaled-generic musculoskeletal model for each child walking, i.e. (i) static optimization that minimized summed-excitation squared; (ii) static optimization with maximum isometric muscle forces scaled to body mass; (iii) an EMG-assisted approach using optimization to minimize summed-excitation squared while reducing tracking errors of experimental EMG-linear envelopes and joint moments; and (iv) EMG-assisted with musculotendon model parameters first personalized by calibration. Both static optimization approaches showed a relatively low model performance compared to EMG envelopes. EMG-assisted approaches performed much better, especially in CP, with only a minor mismatch in joint moments. Calibration did not affect model performance significantly, however it did affect musculotendon forces, especially in CP. A model more consistent with experimental measures is more likely to yield more physiologically representative results. Therefore, this study highlights the importance of calibrated EMG-assisted modelling when estimating musculotendon forces in TD children and even more so in children with CP.
A comparison of the magnitude and duration of linear and rotational head accelerations generated during hand-, elbow- and shoulder-to-head checks delivered by hockey players.
Venue: Journal of biomechanics

Ice hockey has the highest rates for concussion among team sports in Canada. In elite play, the most common mechanism is impact to the head by an opposing player's upper limb, with shoulder-to-head impacts accounting for twice as many concussions as elbow- and hand-to-head impacts combined. Improved understanding of the biomechanics of head impacts in hockey may inform approaches to prevention. In this study, we measured the magnitude and duration of linear and rotational head accelerations when hockey players (n = 11; aged 21-25) delivered checks "as hard as comfortable" to the head of an instrumented dummy with their shoulder, elbow and hand. There were differences in both peak magnitude and duration of head accelerations across upper limb impact sites, based on repeated-measures ANOVA (p < 0.005). Peak linear head accelerations averaged 1.9-fold greater for hand and 1.3-fold greater for elbow than shoulder (mean values = 20.35, 14.23 and 10.55 g, respectively). Furthermore, peak rotational head accelerations averaged 2.1-fold greater for hand and 1.8-fold greater for elbow than shoulder (1097.9, 944.1 and 523.1 rad/s2, respectively). However, times to peak linear head acceleration (a measure of the duration of the acceleration impulse) were 2.1-fold longer for shoulder than elbow, and 2.5-fold longer for shoulder than hand (12.26, 5.94 and 4.98 ms, respectively), and there were similar trends in the durations of rotational head acceleration. Our results show that, in body checks to the head delivered by varsity-level hockey players, shoulder-to-head impacts generated longer durations but lower magnitude of peak head acceleration than elbow- and hand-to-head impacts.
A bipedal compliant walking model generates periodic gait cycles with realistic swing dynamics.
Venue: Journal of biomechanics

A simple spring mechanics model can capture the dynamics of the center of mass (CoM) during human walking, which is coordinated by multiple joints. This simple spring model, however, only describes the CoM during the stance phase, and the mechanics involved in the bipedality of the human gait are limited. In this study, a bipedal spring walking model was proposed to demonstrate the dynamics of bipedal walking, including swing dynamics followed by the step-to-step transition. The model consists of two springs with different stiffnesses and rest lengths representing the stance leg and swing leg. One end of each spring has a foot mass, and the other end is attached to the body mass. To induce a forward swing that matches the gait phase, a torsional hip joint spring was introduced at each leg. To reflect the active knee flexion for foot clearance, the rest length of the swing leg was set shorter than that of the stance leg, generating a discrete elastic restoring force. The number of model parameters was reduced by introducing dependencies among stiffness parameters. The proposed model generates periodic gaits with dynamics-driven step-to-step transitions and realistic swing dynamics. While preserving the mimicry of the CoM and ground reaction force (GRF) data at various gait speeds, the proposed model emulated the kinematics of the swing leg. This result implies that the dynamics of human walking generated by the actuations of multiple body segments is describable by a simple spring mechanics.
The relative contributions of sagittal, frontal, and transverse joint works to self-paced incline and decline slope walking.
Venue: Journal of biomechanics

Positive and negative work are generated at the lower limb joints in order to locomote over various terrains. Joint work quantifies the changes in energy that are necessary to adapt gait to environmental demands. The aim of this study was to quantify 3D joint work at the hip, knee, and ankle during slope walking. Work was calculated for ten males (23.9 ± 1.1 years) walking at a self-selected speed on inclines and declines (-20, -12, -6, 0, 6, 12, 20 degrees). Sagittal positive work significantly increased at the hip, knee, and ankle for incline walking (for example, hip positive work increased 153%, 280%, and 453% for 6, 12, and 20 degrees, respectively; knee and ankle positive work also increased) (p ≤ 0.05), in order to raise and propel the body forward. Sagittal negative work increased significantly at the hip, knee and ankle for decline walking (for example, knee negative work increased 193%, 355%, and 496% for -6, -12, and -20 degrees, respectively; hip and ankle negative work also increased) (p ≤ 0.05), in order to control body descent. These substantial changes in work will be especially challenging for people with compromised strength due to age and disease. Furthermore, changes in work were not limited to the sagittal plane: 46% of the total hip joint work occurred in the frontal and transverse planes for six degree decline walking. Thus, decline walking placed greater demands on the hip ab/adductors and rotators, and this may be related to the greater risk of falls observed for descent versus ascent.
The impact of centre of pressure error on predicted joint kinetics during cerebral palsy and typically developed gait: A clinical perspective.
Venue: Journal of biomechanics

Centre of Pressure (CoP) location error is common when using kinematic and kinetic data to predict intersegmental forces and net joint moments during gait. Changes in peak moments due to CoP error have been reported in the literature. However, debate exists as to what levels of error are acceptable. The aim of this study was to examine the impact of CoP error on the kinetic profiles of children with typical development (TD) and children with cerebral palsy (CP) during gait. Three-dimensional kinematic and kinetic data were recorded and simulated CoP errors were applied at 3 mm, 6 mm, 9 mm, 12 mm increments in both positive and negative anteroposterior and mediolateral directions. Absolute differences in maximum kinetic parameters between increments were assessed in conjunction with changes in the Gait Deviation Index-Kinetic (GDI-Kinetic). Changes in GDI-Kinetic above 3.6 points were considered clinically significant. Maximum peak changes of up to 24.8% (CP) and 34.7% (TD) (sagittal plane) and up to 36.8% (CP) and 61.5% (TD) (coronal plane) were demonstrated at the knee. While absolute percentage differences were high at some error increments, GDI-Kinetic results suggested that such large percentage differences may still be clinically acceptable. Children with TD demonstrated clinically significant changes in GDI-Kinetic for CoP displacements of 9 mm and 12 mm, corresponding to 23% and 35% absolute differences in maximum moments. In contrast, the clinically significant threshold was not reached for children with CP that may be related to a slower walking speed. The findings of this study highlight the need for laboratories to consider the thresholds currently used for CoP error, which will help guide quality assurance procedures.
Numerical simulation of haemodynamics of the descending aorta in the non-diabetic and diabetic rabbits.
Venue: Journal of biomechanics

Diabetes mellitus (DM) is a predisposing risk factor leading to macrovascular diseases. Changes in haemodynamics of the diabetic aortas remain largely unclear and relevant computational analyses are lacking in the literature. Ten adult rabbits (1.6-2.2 kg) were collected and the type I diabetic rabbit model was induced by injection of alloxan. A total of five control and five diabetic rabbit aortas were considered for subsequent numerical simulation. The CT scanning was performed to reconstruct three-dimensional model of the individual rabbit descending aorta. The flow velocity waveforms were measured by ultrasound machine and were set to be the inlet boundary conditions. The reconstructed aortas were then imported into ANSYS to perform mesh generation and computational analysis. Results showed that the distributions of haemodynamic indicators time-averaged wall shear stress (TAWSS), oscillating shear index (OSI) and transverse wall shear stress (transWSS) in the non-diabetic rabbit aortas were similar to those in the diabetic rabbit aortas. However, the mean values of TAWSS and transWSS in the non-diabetic rabbit aortas were significantly higher than those values in the diabetic rabbit aortas (TAWSS: p = 0.04; transWSS: p = 0.02). The back of right renal artery tended to have high OSI in both the non-diabetic and the diabetic rabbit aortas. Notably, the regions with high OSI tended to have intense disturbed flow and low TAWSS in the most diabetic rabbit aortas. The results suggest that diabetes leads to changes in haemodynamic parameters in the rabbit aortas. In particular, the lower TAWSS and the higher OSI within the diabetic aortas may further contribute to aortic wall remodeling.
Comparison of porcine brain mechanical properties to potential tissue simulant materials in quasi-static and sinusoidal compression.
Venue: Journal of biomechanics

In both finite element and physical surrogate models of head blast injury, accurate material properties of the brain and/or tissue simulants are necessary to ensure biofidelity in predicted response. Thus, there is a need for experimental comparisons between tissue and simulant materials under the same experimental conditions. This study compares the response of porcine brain tissue and a variety of brain tissue simulants in quasi-static and sinusoidal compression tests. Fresh porcine brain tissue was obtained from a local abattoir and tested within 4 h post mortem. Additionally, the effect of post mortem time was investigated by comparing samples stored at room temperature and stored frozen (-18 °C), at various time intervals. The brain tissue simulants tested were bovine gelatin (3%, 5%, and 10% concentration), agarose gelatin (e0.4%, 0.6%, 0.8% concentration), and Sylgard 527. The experiments were performed using a DMA apparatus (TA Instruments Q800). The quasi-static compression data were fit to Ogden hyperelastic functions so that parameters could be compared. It was found that bovine gelatin at 3% and 5% concentration demonstrated the closest response to brain tissue in quasi-static compression. Conversely, in sinusoidal compression, the agarose gel and Sylgard 527 were found to be in closer agreement with the tissue, than bovine gel. In terms of post mortem time and storage, there was no statistically significant difference detected in the response of tissue samples after 48 h, regardless of storage method. However, samples stored at room temperature after 48 h appeared to demonstrate a reduction in stiffness.