nils-herrmann/clean_scirepeval huggingface dataset:
Tensile behaviour of individual fibre bundles in the human lumbar anulus fibrosus.
Venue: Journal of biomechanics

Disc degeneration is a common medical affliction whose origins are not fully understood. An improved understanding of its underlying mechanisms could lead to the development of more effective treatments. The aim of this paper was to investigate the effect of (1) degeneration, (2) circumferential region and (3) strain rate on the microscale mechanical properties (toe region modulus, linear modulus, extensibility, phase angle) of individual fibre bundles in the anulus fibrosus lamellae of the human intervertebral disc. Healthy and degenerate fibre bundles excised from different circumferential regions in the outer anulus (posterolateral, lateral, anterolateral, anterior) were tensile tested at slow (0.1%/s), medium (1%/s) and fast (10%/s) strain rates using a micromechanical testing system. Our preliminary results showed that neither degeneration nor circumferential region significantly affected the fibre bundles' mechanical behaviour. However, when the fibre bundles were tested at higher strain rates, this resulted in significantly higher linear moduli and lower phase angles. These findings, compared with data from other studies investigating single and multiple lamellae sections, suggest that degeneration has minimal effect on outer anulus mechanics irrespective of structural level, and the inter- and intra-lamellar arrangement and continuity of the fibre bundles may influence the lamellae's regional behaviour and viscoelasticity.
Discrete element analysis is a valid method for computing joint contact stress in the hip before and after acetabular fracture.
Venue: Journal of biomechanics

Evaluation of abnormalities in joint contact stress that develop after inaccurate reduction of an acetabular fracture may provide a potential means for predicting the risk of developing post-traumatic osteoarthritis. Discrete element analysis (DEA) is a computational technique for calculating intra-articular contact stress distributions in a fraction of the time required to obtain the same information using the more commonly employed finite element analysis technique. The goal of this work was to validate the accuracy of DEA-computed contact stress against physical measurements of contact stress made in cadaveric hips using Tekscan sensors. Four static loading tests in a variety of poses from heel-strike to toe-off were performed in two different cadaveric hip specimens with the acetabulum intact and again with an intentionally malreduced posterior wall acetabular fracture. DEA-computed contact stress was compared on a point-by-point basis to stress measured from the physical experiments. There was good agreement between computed and measured contact stress over the entire contact area (correlation coefficients ranged from 0.88 to 0.99). DEA-computed peak contact stress was within an average of 0.5 MPa (range 0.2-0.8 MPa) of the Tekscan peak stress for intact hips, and within an average of 0.6 MPa (range 0-1.6 MPa) for fractured cases. DEA-computed contact areas were within an average of 33% of the Tekscan-measured areas (range: 1.4-60%). These results indicate that the DEA methodology is a valid method for accurately estimating contact stress in both intact and fractured hips.
A comparison of patellofemoral cartilage morphology and deformation in anterior cruciate ligament deficient versus uninjured knees.
Venue: Journal of biomechanics

Anterior cruciate ligament (ACL) deficient patients have an increased rate of patellofemoral joint (PFJ) osteoarthritis (OA) as compared to the general population. Although the cause of post-injury OA is multi-factorial, alterations in joint biomechanics may predispose patients to cartilage degeneration. This study aimed to compare in vivo PFJ morphology and mechanics between ACL deficient and intact knees in subjects with unilateral ACL ruptures. Eight male subjects underwent baseline MRI scans of both knees. They then performed a series of 60 single-legged hops, followed by a post-exercise MRI scan. This process was repeated for the contralateral knee. The MR images were converted into three-dimensional surface models of cartilage and bone in order to assess cartilage thickness distributions and strain following exercise. Prior to exercise, patellar cartilage was significantly thicker in intact knees as compared to ACL deficient knees by 1.8%. In response to exercise, we observed average patellar cartilage strains of 5.4 ± 1.1% and 2.5 ± 1.4% in the ACL deficient and intact knees, respectively. Importantly, the magnitude of patellar cartilage strain in the ACL deficient knees was significantly higher than in the intact knees. However, while trochlear cartilage experienced a mean strain of 2.4 ± 1.6%, there was no difference in trochlear cartilage strain between the ACL deficient and uninjured knees. In summary, we found that ACL deficiency was associated with decreased patellar cartilage thickness and increased exercise-induced patellar cartilage strain when compared to the uninjured contralateral knees.
Role of heel lifting in standing balance recovery: A simulation study.
Venue: Journal of biomechanics

Although lifting the heels has frequently been observed during balance recovery, the function of this movement has generally been overlooked. The present study aimed to investigate the functional role of heel lifting during regaining balance from a perturbed state. Computer simulation was employed to objectively examine the effect of allowing/constraining heel lifting on balance performance. The human model consisted of 3 rigid body segments connected by frictionless joints. Movements were driven by joint torques depending on current joint angle, angular velocity, and activation level. Starting from forward-inclined and static straight-body postures, the optimization goal was to recover balance effectively (so that ground projection of the mass center returned to the inside of the base of support) and efficiently by adjusting ankle and hip joint activation levels. Allowing/constraining heel lifting resulted in virtually identical movements when balance was mildly perturbed at the smallest lean angle (8°). At larger lean angles (8.5° and 9°), heel lifting assisted balance recovery more evidently with larger joint movements. Partial and altered timings of ankle/hip torque activation due to constraining heel lifting reduced linear and angular momentum generation for avoiding forward falling, and resulted in hindered balancing performance.
MRI-based experimentations of fingertip flat compression: Geometrical measurements and finite element inverse simulations to investigate material property parameters.
Venue: Journal of biomechanics

Modeling human-object interactions is a necessary step in the ergonomic assessment of products. Fingertip finite element models can help investigating these interactions, if they are built based on realistic geometrical data and material properties. The aim of this study was to investigate the fingertip geometry and its mechanical response under compression, and to identify the parameters of a hyperelastic material property associated to the fingertip soft tissues. Fingertip compression tests in an MRI device were performed on 5 subjects at either 2 or 4 N and at 15° or 50°. The MRI images allowed to document both the internal and external fingertip dimensions and to build 5 subject-specific finite element models. Simulations reproducing the fingertip compression tests were run to obtain the material property parameters of the soft tissues. Results indicated that two ellipses in the sagittal and longitudinal plane could describe the external fingertip geometry. The internal geometries indicated an averaged maximal thickness of soft tissues of 6.4 ± 0.8 mm and a 4 ± 1 mm height for the phalanx bone. The averaged deflections under loading went from 1.8 ± 0.3 mm at 2 N, 50° to 3.1 ± 0.2 mm at 4 N, 15°. Finally, the following set of parameters for a second order hyperelastic law to model the fingertip soft tissues was proposed: C01=0.59 ± 0.09 kPa and C20 = 2.65 ± 0.88 kPa. These data should facilitate further efforts on fingertip finite element modeling.
Microstructural changes associated with osteoporosis negatively affect loading-induced fluid flow around osteocytes in cortical bone.
Venue: Journal of biomechanics

Loading-induced interstitial fluid flow in the microporosities of bone is critical for osteocyte mechanotransduction and for the maintenance of tissue health, enhancing convective transport in the lacunar-canalicular system. In recent studies, our group has reported alterations of bone's vascular porosity and lacunar-canalicular system microarchitecture in a rat model of postmenopausal osteoporosis. In this work, poroelastic finite element analysis was used to investigate whether these microstructural changes can affect interstitial fluid flow around osteocytes. Animal-specific finite element models were developed combining micro-CT reconstructions of bone microstructure and measures of the poroelastic material properties. These models were used to quantify and compare loading-induced fluid flow in the lacunar-canalicular system of ovariectomized and sham-operated rats. A parametric analysis was also used to quantify the influence of the lacunar-canalicular permeability and vascular porosity on the fluid velocity magnitude. Results show that mechanically-induced interstitial fluid velocity can be significantly reduced in the lacunar-canalicular system of ovariectomized rats. Interestingly, the vascular porosity is shown to have a major influence on interstitial fluid flow, while the lacunar-canalicular permeability influence is limited when larger than 10-20m2. Altogether our results suggest that microstructural changes associated with the osteoporotic condition can negatively affect interstitial fluid flow around osteocytes in the lacunar-canalicular system of cortical bone. This fluid flow reduction could impair mechanosensation of the osteocytic network, possibly playing a role in the initiation and progression of age-related bone loss and postmenopausal osteoporosis.
Walking speed differentially alters spinal loads in persons with traumatic lower limb amputation.
Venue: Journal of biomechanics

Persons with lower limb amputation (LLA) perceive altered motions of the trunk/pelvis during activities of daily living as contributing factors for low back pain. When walking (at a singular speed), larger trunk motions among persons with vs. without LLA are associated with larger spinal loads; however, modulating walking speed is necessary in daily life and thus understanding the influences of walking speed on spinal loads in persons with LLA is of particular interest here. Three-dimensional trunk-pelvic kinematics, collected during level-ground walking at self-selected (SSW) and two controlled speeds (∼1.0 and ∼1.4 m/s), were obtained for seventy-eight participants: 26 with transfemoral and 26 with transtibial amputation, and 26 uninjured controls (CTR). Using a kinematics-driven, non-linear finite element model of the lower back, the resultant compressive and mediolateral/anteroposterior shear loads at the L5/S1 spinal level were estimated. Peak values were extracted and compiled. Despite walking slower at SSW speeds (∼0.21 m/s), spinal loads were 8-14% larger among persons with transfemoral amputation vs. CTR. Across all participants, peak compressive, mediolateral, and anteroposterior shear loads increased with increasing walking speed. At the fastest (vs. slowest) controlled speed, these increases were respectively 24-84% and 29-77% larger among persons with LLA relative to CTR. Over time, repeated exposures to these increased spinal loads, particularly at faster walking speeds, may contribute to the elevated risk for low back pain among persons with LLA. Future work should more completely characterize relative risk in daily life between persons with vs. without LLA by analyzing additional activities and tissue-level responses.
Fully-coupled mathematical modeling of actomyosin-cytosolic two-phase flow in a highly deformable moving Keratocyte cell.
Venue: Journal of biomechanics

Interaction between intracellular dynamics and extracellular matrix (ECM) generally occurred into very thin fragment of moving cell, namely lamellipodia, enables all movable cells to crawl on ECM. In fast-moving cells such as fish Keratocytes, Lamellipodia including most cell area finds a fan-like shape during migration, with a variety of aspect ratio function of fish type. In this work, our purpose is to present a novel and more complete two-dimensional continuum mathematical model of actomyosin-cytosolic two-phase flow of a self-deforming Keratocyte with circular spreaded to steady fan-like shape. In the new approach, in addition to the two-phase flow of the F-actin and cytosol, the G-actin transport was spatiotemporally modeled. We also for the first time modeled the effect of variable volume fraction of the moving F-actin porous network on solute transport in the cytosolic fluid. Our novel fully-coupled mathematical model provides a better understanding of intracellular dynamics of fast-migrating Keratocytes; such as the F-actin centripetal and cytosolic fountain-like flows, free-active myosin distribution, distribution sequence of the G-actin, F-actin, and myosin, and myosin-induced pressure flied of cytoplasm as well as the map of intracellular forces like myosin contraction and adhesion traction. All these results are qualitatively and quantitatively in good agreement with experimental observations. According to a range of value of parameters used in this model, our steady state of moving Keratocyte finds fan-like shape with the same aspect ratio as wide category of fish Keratocytes. This new model can predict shape of Keratocytes in other range of parameter values.
A practical step length algorithm using lower limb angular velocities.
Venue: Journal of biomechanics

The use of Inertial Measurement Units (IMUs) for spatial gait analysis has opened the door to unconstrained measurements within the home and community. Bandwidth, cost limitations, and ease of use has historically restricted the number and location of sensors worn on the body. In this paper, we describe a four-sensor configuration of IMUs placed on the shanks and thighs that is sufficient to provide an accurate measure of temporal gait parameters, spatial gait parameters, and joint angle dynamics during ambulation. Estimating spatial gait parameters solely from gyroscope data is preferred because gyroscopes are less susceptible to sensor noise and a system comprised of only gyroscopes uses decreased bandwidth compared to a typical 9 degree-of-freedom IMU. The purpose of this study was to determine the validity of a novel method of step length estimation using gyroscopes attached to the shanks and thighs. An Inverted Pendulum Model algorithm (IPM) was proposed to calculate step length, stride length, and gait speed. The algorithm incorporates heel-strike events and average forward velocity per step to make these assessments. IMU algorithm accuracy was determined via concurrent validity with an instrumented walkway and results explained via the collision model of gait. The IPM produced accurate estimates of step length, stride length, and gait speed with a mean difference of 3 cm and an RMSE of 6.6 cm for step length, thus establishing a new approach for spatial gait parameter calculation. The lack of numerical integration in IPM makes it well suited for use in continuous monitoring applications where sensor sampling rates are restricted.
Trunk musculoskeletal response in maximum voluntary exertions: A combined measurement-modeling investigation.
Venue: Journal of biomechanics

Maximum voluntary exertion (MVE) tasks quantify trunk strength and maximal muscle electromyography (EMG) activities with both clinical and biomechanical implications. The aims here are to evaluate the performance of an existing trunk musculoskeletal model, estimate maximum muscle stresses and spinal forces, and explore likely differences between males and females in maximum voluntary exertions. We, therefore, measured trunk strength and EMG activities of 19 healthy right-handed subjects (9 females and 10 males) in flexion, extension, lateral and axial directions. MVEs for all subjects were then simulated in a subject-specific trunk musculoskeletal model, and estimated muscle activities were compared with EMGs. Analysis of variance was used to compare measured moments and estimated spinal loads at the L5-S1 level between females and males. MVE moments in both sexes were greatest in extension (means of 236 Nm in males and 190 Nm in females) and least in left axial torque (97 Nm in males and 64 Nm in females). Being much greater in lateral and axial MVEs, coupled moments reached ∼50% of primary moments in average. Females exerted less moments in all directions reaching significance except in flexion. Muscle activity estimations were strongly correlated with measurements in flexion and extension (Pearson's r = 0.69 and 0.76), but the correlations were very weak in lateral and axial MVEs (Pearson's r = 0.27 and 0.13). Maximum muscle stress was in average 0.80 ± 0.42 MPa but varied among muscles from 0.40 ± 0.22  MPa in rectus abdominis to 0.99 ± 0.29 MPa in external oblique. To estimate maximum muscle stresses and evaluate validity of a musculoskeletal model, MVEs in all directions with all coupled moments should be considered.
A mathematical analysis to address the 6 degree-of-freedom segmental power imbalance.
Venue: Journal of biomechanics

Segmental power is used in human movement analyses to indicate the source and net rate of energy transfer between the rigid bodies of biomechanical models. Segmental power calculations are performed using segment endpoint dynamics (kinetic method). A theoretically equivalent method is to measure the rate of change in a segment's mechanical energy state (kinematic method). However, these two methods have not produced experimentally equivalent results for segments proximal to the foot, with the difference in methods deemed the "power imbalance." In a 6 degree-of-freedom model, segments move independently, resulting in relative segment endpoint displacement and non-equivalent segment endpoint velocities at a joint. In the kinetic method, a segment's distal end translational velocity may be defined either at the anatomical end of the segment or at the location of the joint center (defined here as the proximal end of the adjacent distal segment). Our mathematical derivations revealed the power imbalance between the kinetic method using the anatomical definition and the kinematic method can be explained by power due to relative segment endpoint displacement. In this study, we tested this analytical prediction through experimental gait data from nine healthy subjects walking at a typical speed. The average absolute segmental power imbalance was reduced from 0.023 to 0.046 W/kg using the anatomical definition to ≤0.001 W/kg using the joint center definition in the kinetic method (95.56-98.39% reduction). Power due to relative segment endpoint displacement in segmental power analyses is substantial and should be considered in analyzing energetic flow into and between segments.
Sensitivity of intervertebral joint forces to center of rotation location and trends along its migration path.
Venue: Journal of biomechanics

Translational vertebral motion during functional tasks manifests itself in dynamic loci for center of rotation (COR). A shift of COR affects moment arms of muscles and ligaments; consequently, muscle and joint forces are altered. Based on posture- and level-specific trends of COR migration revealed by in vivo dynamic radiography during functional activities, it was postulated that the instantaneous COR location for a particular joint is optimized in order to minimize the joint reaction forces. A musculoskeletal multi-body model was employed to investigate the hypotheses that (1) a posterior COR in upright standing and (2) an anterior COR in forward flexed posture leads to optimized lumbar joint loads. Moreover, it was hypothesized that (3) lower lumbar levels benefit from a more superiorly located COR. The COR in the model was varied from its initial position in posterior-anterior and inferior-superior direction up to ±6 mm in steps of 2 mm. Movement from upright standing to 45° forward bending and backwards was simulated for all configurations. Joint reaction forces were computed at levels L2L3 to L5S1. Results clearly confirmed hypotheses (1) and (2) and provided evidence for the validity of hypothesis (3), hence offering a biomechanical rationale behind the migration paths of CORs observed during functional flexion/extension movement. Average sensitivity of joint force magnitudes to an anterior shift of COR was +6 N/mm in upright and -21 N/mm in 30° forward flexed posture, while sensitivity to a superior shift in upright standing was +7 N/mm and -8N/mm in 30° flexion. The relation between COR loci and joint loading in upright and flexed postures could be mainly attributed to altered muscle moment arms and consequences on muscle exertion. These findings are considered relevant for the interpretation of COR migration data, the development of numerical models, and could have an implication on clinical diagnosis and treatment or the development of spinal implants.
Predictive multiscale computational model of shoe-floor coefficient of friction.
Venue: Journal of biomechanics

Understanding the frictional interactions between the shoe and floor during walking is critical to prevention of slips and falls, particularly when contaminants are present. A multiscale finite element model of shoe-floor-contaminant friction was developed that takes into account the surface and material characteristics of the shoe and flooring in microscopic and macroscopic scales. The model calculates shoe-floor coefficient of friction (COF) in boundary lubrication regime where effects of adhesion friction and hydrodynamic pressures are negligible. The validity of model outputs was assessed by comparing model predictions to the experimental results from mechanical COF testing. The multiscale model estimates were linearly related to the experimental results (p < 0.0001). The model predicted 73% of variability in experimentally-measured shoe-floor-contaminant COF. The results demonstrate the potential of multiscale finite element modeling in aiding slip-resistant shoe and flooring design and reducing slip and fall injuries.
Corrigendum to "Alterations of musculoskeletal models for a more accurate estimation of lower limb joint contact forces during normal gait: A systematic review" [J. Biomech. 63 (2017) 8-20].
Venue: Journal of biomechanics

The authors regret to inform the readers that the Acknowledgment section of the paper did not include the funding details for L. Modenese. It is as below: L. Modenese was supported by the European Commission (7th FP, ICT integrated project MD-Paedigree, Contract Number 600932) and by the EPSRC (Frontier Engineering Awards, EP/K03877X/1). The authors would like to apologise for any inconvenience caused.