nils-herrmann/clean_scirepeval huggingface dataset:
A comparison and update of direct kinematic-kinetic models of leg stiffness in human running.
Venue: Journal of biomechanics

Direct kinematic-kinetic modelling currently represents the "Gold-standard" in leg stiffness quantification during three-dimensional (3D) motion capture experiments. However, the medial-lateral components of ground reaction force and leg length have been neglected in current leg stiffness formulations. It is unknown if accounting for all 3D would alter healthy biologic estimates of leg stiffness, compared to present direct modelling methods. This study compared running leg stiffness derived from a new method (multiplanar method) which includes all three Cartesian axes, against current methods which either only include the vertical axis (line method) or only the plane of progression (uniplanar method). Twenty healthy female runners performed shod overground running at 5.0 m/s. Three-dimensional motion capture and synchronised in-ground force plates were used to track the change in length of the leg vector (hip joint centre to centre of pressure) and resultant projected ground reaction force. Leg stiffness was expressed as dimensionless units, as a percentage of an individual's bodyweight divided by standing leg length (BW/LL). Leg stiffness using the line method was larger than the uniplanar method by 15.6%BW/LL (P < .001), and multiplanar method by 24.2%BW/LL (P < .001). Leg stiffness from the uniplanar method was larger than the multiplanar method by 8.5%BW/LL (6.5 kN/m) (P < .001). The inclusion of medial-lateral components significantly increased leg deformation magnitude, accounting for the reduction in leg stiffness estimate with the multiplanar method. Given that limb movements typically occur in 3D, the new multiplanar method provides the most complete accounting of all force and length components in leg stiffness calculation.
Chiral behavior in rat tail tendon fascicles.
Venue: Journal of biomechanics

Ex vivo tendon mechanical behavior has been well described under rotationally constrained uniaxial tensile testing. During standard loading of rat tail tendon (RTT) fascicles, apparent axial twist has been observed. To quantify this behavior, we designed a custom testing setup, utilizing magnetic suspension, to allow unconstrained axial rotation during tensile loading. We characterized the rotational behavior of single and paired RTT fascicles under cyclic loading. We also measured stress relaxation across loading cycles as well as "rotational relaxation". Single fascicle nonlinear stretch-twist coupling is well described by the asymptotic function Δθ=A(1-e-Bε) in which fascicles rotated a mean ±51.1° within about 1% applied axial strain. On average, paired fascicles rotated just over 10° less. Specimen cross-sectional diameter had a noticeable effect on the measured mechanical properties, particularly effective elastic modulus. Such stretch-twist coupling and size dependence cannot be understood via classical elasticity but is predicted by Cosserat (micropolar) elasticity. The current study demonstrates RTT fascicles are chiral based on observed axial load-induced twist. Additionally, our findings support existing research that suggests a helical fascicle structure. Potential consequences of helical substructures, mechanical and biological, merit further investigation.
Urethral lumen occlusion by artificial sphincteric device: Evaluation of degraded tissues effects.
Venue: Journal of biomechanics

Urinary incontinence can be surgically treated by means of artificial sphincters, based on a cuff that provides a pressure around the urethra to occlude the lumen. Considering the frequent access of elderly patients to this surgical practice, tissue degradation phenomena must be investigated, since they could affect treatment reliability and durability. The potential degradation can be interpreted considering a variation within soft tissue constitutive formulation, by means of a correlation between mechanical properties and tissues ageing. The overall compressibility varies, as characteristics aspect of soft tissue mechanical response with age, as well as the stiffness. The investigation is performed by means of a three dimensional numerical model of the urethral duct. The effects of the interaction phenomenon with a cuff is interpreted considering the changes, within the constitutive models, of the basic parameters that define the potential degradation process. The deformation related to compressibility is recalled, ranging between ten and fifty percent in dependence on the degradation level considered. This parameter, reported mostly as representative of the aging effect, shows a large variation that confirms the relevance of the investigation performed toward a sensitivity of the mechanical response of the urethral duct referred to the lumen occlusion.
A computational model of the hepatic circulation applied to analyze the sensitivity of hepatic venous pressure gradient (HVPG) in liver cirrhosis.
Venue: Journal of biomechanics

Measurement of hepatic venous pressure gradient (HVPG) is currently widely adopted to provide an estimate of portal pressure gradient (PPG) in the diagnosis and treatment of portal hypertension associated with liver cirrhosis. Despite the well-documented clinical utility of HVPG, it remains poorly understood how the relationship between HVPG and PPG is affected by factors involved in the pathogenesis and progression of cirrhosis. In the study, a computational model of the hepatic circulation calibrated to in vivo data was developed to simulate the procedure of HVPG measurement and quantitatively investigate the error of HVPG relative to PPG under various pathophysiological conditions. Obtained results confirmed the clinical consensus that HVPG is applicable to the assessment of portal hypertension caused by increased vascular resistance located primarily at the sinusoidal and postsinusoidal sites rather than at the presinusoidal site. On the other hand, our study demonstrated that the accuracy of HVPG measurement was influenced by many factors related to hepatic hemodynamics even in the case of sinusoidal portal hypertension. For instance, varying presinusoidal portal vascular resistance significantly altered the difference between HVPG and PPG, while an enhancement in portosystemic collateral flow tended to improve the accuracy of HVPG measurement. Moreover, it was found that presinusoidal and postsinusoidal vascular resistances interfered with each other with respect to their influence on HVPG measurement. These findings suggest that one should take into account patient-specific pathological conditions in order to achieve a better understanding and utilization of HVPG in the clinical practice.
A comprehensive assessment of the musculoskeletal system: The CAMS-Knee data set.
Venue: Journal of biomechanics

Combined knowledge of the functional kinematics and kinetics of the human body is critical for understanding a wide range of biomechanical processes including musculoskeletal adaptation, injury mechanics, and orthopaedic treatment outcome, but also for validation of musculoskeletal models. Until now, however, no datasets that include internal loading conditions (kinetics), synchronized with advanced kinematic analyses in multiple subjects have been available. Our goal was to provide such datasets and thereby foster a new understanding of how in vivo knee joint movement and contact forces are interlinked - and thereby impact biomechanical interpretation of any new knee replacement design. In this collaborative study, we have created unique kinematic and kinetic datasets of the lower limb musculoskeletal system for worldwide dissemination by assessing a unique cohort of 6 subjects with instrumented knee implants (Charité - Universitätsmedizin Berlin) synchronized with a moving fluoroscope (ETH Zürich) and other measurement techniques (including whole body kinematics, ground reaction forces, video data, and electromyography data) for multiple complete cycles of 5 activities of daily living. Maximal tibio-femoral joint contact forces during walking (mean peak 2.74 BW), sit-to-stand (2.73 BW), stand-to-sit (2.57 BW), squats (2.64 BW), stair descent (3.38 BW), and ramp descent (3.39 BW) were observed. Internal rotation of the tibia ranged from 3° external to 9.3° internal. The greatest range of anterio-posterior translation was measured during stair descent (medial 9.3 ± 1.0 mm, lateral 7.5 ± 1.6 mm), and the lowest during stand-to-sit (medial 4.5 ± 1.1 mm, lateral 3.7 ± 1.4 mm). The complete and comprehensive datasets will soon be made available online for public use in biomechanical and orthopaedic research and development.
Estimation of in vivo inter-vertebral loading during motion using fluoroscopic and magnetic resonance image informed finite element models.
Venue: Journal of biomechanics

Finite element (FE) models driven by medical image data can be used to estimate subject-specific spinal biomechanics. This study aimed to combine magnetic resonance (MR) imaging and quantitative fluoroscopy (QF) in subject-specific FE models of upright standing, flexion and extension. Supine MR images of the lumbar spine were acquired from healthy participants using a 0.5 T MR scanner. Nine 3D quasi-static linear FE models of L3 to L5 were created with an elastic nucleus and orthotropic annulus. QF data was acquired from the same participants who performed trunk flexion to 60° and trunk extension to 20°. The displacements and rotations of the vertebrae were calculated and applied to the FE model. Stresses were averaged across the nucleus region and transformed to the disc co-ordinate system (S1 = mediolateral, S2 = anteroposterior, S3 = axial). In upright standing S3 was predicted to be -0.7 ± 0.6 MPa (L3L4) and -0.6 ± 0.5 MPa (L4L5). S3 increased to -2.0 ± 1.3 MPa (L3L4) and -1.2 ± 0.6 MPa (L4L5) in full flexion and to -1.1 ± 0.8 MPa (L3L4) and -0.7 ± 0.5 MPa (L4L5) in full extension. S1 and S2 followed similar patterns; shear was small apart from S23. Disc stresses correlated to disc orientation and wedging. The results demonstrate that MR and QF data can be combined in a participant-specific FE model to investigate spinal biomechanics in vivo and that predicted stresses are within ranges reported in the literature.
Impact of ankle foot orthosis stiffness on Achilles tendon and gastrocnemius function during unimpaired gait.
Venue: Journal of biomechanics

Ankle foot orthoses (AFOs) are designed to improve gait for individuals with neuromuscular conditions and have also been used to reduce energy costs of walking for unimpaired individuals. AFOs influence joint motion and metabolic cost, but how they impact muscle function remains unclear. This study investigated the impact of different stiffness AFOs on medial gastrocnemius muscle (MG) and Achilles tendon (AT) function during two walking speeds. We performed gait analyses for eight unimpaired individuals. Each individual walked at slow and very slow speeds with a 3D printed AFO with no resistance (free hinge condition) and four levels of ankle dorsiflexion stiffness: 0.25Nm/°, 1Nm/°, 2Nm/°, and 3.7Nm/°. Motion capture, ultrasound, and musculoskeletal modeling were used to quantify MG and AT lengths with each AFO condition. Increasing AFO stiffness increased peak AFO dorsiflexion moment with decreased peak knee extension and peak ankle dorsiflexion angles. Overall musculotendon length and peak AT length decreased, while peak MG length increased with increasing AFO stiffness. Peak MG activity, length, and velocity significantly decreased with slower walking speed. This study provides experimental evidence of the impact of AFO stiffness and walking speed on joint kinematics and musculotendon function. These methods can provide insight to improve AFO designs and optimize musculotendon function for rehabilitation, performance, or other goals.
Fiber orientation effects in simple shearing of fibrous soft tissues.
Venue: Journal of biomechanics

Fiber-reinforcement is a common feature of many soft biological tissues. Continuum mechanics modeling of the mechanical response of such tissues using transversely isotropic hyperelasticity is now well developed. The fundamental deformation of simple shear within this framework is examined here. It is well known that the normal stress effect characteristic of nonlinear elasticity plays a crucial role in maintaining a homogeneous deformation state in the bulk of the specimen. Here we consider the effect of anisotropy and fiber-orientation on the shear and normal stresses. It is shown that the confining traction that needs to be applied to the top and bottom faces of a block in order to maintain simple shear can be compressive or tensile depending on the degree of anisotropy and on the angle of orientation of the fibers. In the absence of such an applied traction, an unconfined sample tends to bulge outwards or contract inwards perpendicular to the direction of shear so that one has the possibility of both a positive or negative Poynting effect. The results are illustrated using experimental data obtained by other authors for porcine brain white matter. The general results obtained here are relevant to the development of accurate shear test protocols for the determination of constitutive properties of fibrous biological soft tissues.
Haemodynamic effects of incomplete stent apposition in curved coronary arteries.
Venue: Journal of biomechanics

Incomplete stent apposition (ISA, also known as malapposition) is a complication that affects day-to-day coronary stenting procedures. ISA is more prominent in complex arterial geometries, such as curvature, asa result of the limited conformability of coronary stents. These malapposed struts disturb the otherwise near-wall laminar blood flow and form a micro-recirculation environment. The micro-recirculation environment is often associated with low wall shear stress (WSS) and upsets the delicate balance of vascular biology, providing possible nidus for thrombosis and restenosis. In this study, a three-dimensional (3D) stent model was virtually deployed into an idealised curved coronary artery. Computational fluid dynamics (CFD) simulations were carried out to systematically analyse the haemodynamic effects of increasing maximum ISA distances, ranging from 180 (moderate), 400 (intermediate) to 910μm (severe) in an artery with decreasing radius of curvature (ROC). Micro-recirculations around both proximal and distal malapposed struts become more pronounced as compared to fully-apposed struts. The accompanying areas of low temporally-averaged WSS (AL-TAWSS) can increase twofold compared to the fully-apposed condition. Furthermore, substantial regions (∼5.2% and 9.0%) of AL-TAWSS are detached from the distal end of the malapposed struts in both moderate and intermediate cases respectively. Malapposed stents also induce more variation of TAWSS at the inner bend of the artery. At the stent surface, maximum WSS increases threefold from the fully-apposed case to intermediate ISA. High WSS on the strut surface is known to activate platelets which when exposed to the micro-recirculation environment may lead to their deposition and thrombosis.
Analysis of pelvic strain in different gait configurations in a validated cohort of computed tomography based finite element models.
Venue: Journal of biomechanics

The pelvis functions to transmit upper body loads to the lower limbs and is critical in human locomotion. Semi-automated, landmark-based finite element (FE) morphing and mapping 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 an experimentally validated cohort of specimen-specific FE models of the human pelvis and to use this cohort to analyze pelvic strain patterns during gait. Using an initially segmented specimen-specific pelvic FE model asa source model, four more specimen-specific pelvic FE models were generated from target clinical CT scans using landmark-based morphing and mapping techniques. FE strains from the five models were compared to the experimental strains obtained from cadaveric testing via linear regression analysis, (R2 values ranging from 0.70 to 0.93). Inter-specimen variability in FE strain distributions was seen among the five specimen-specific pelvic FE models. The validated cohort of specimen-specific pelvic FE models was utilized to examine pelvic strains at different phases of the gait cycle. Each validated specimen-specific FE model was reconfigured into gait cycle phases representing heel-strike/heel-off and midstance/midswing. No significant difference was found in the double-leg stance and heel-strike/heel-off models (p=0.40). A trend was observed between double-leg stance and midstance/midswing models (p=0.07), and a significant difference was found between heel-strike/heel-off models and midstance/midswing models (p=0.02). Significant differences were also found in comparing right vs. left models (heel-strike/heel-off p=0.14, midstance/midswing p=0.04).
A model-based approach for estimation of changes in lumbar segmental kinematics associated with alterations in trunk muscle forces.
Venue: Journal of biomechanics

The kinematics information from imaging, if combined with optimization-based biomechanical models, may provide a unique platform for personalized assessment of trunk muscle forces (TMFs). Such a method, however, is feasible only if differences in lumbar spine kinematics due to differences in TMFs can be captured by the current imaging techniques. A finite element model of the spine within an optimization procedure was used to estimate segmental kinematics of lumbar spine associated with five different sets of TMFs. Each set of TMFs was associated with a hypothetical trunk neuromuscular strategy that optimized one aspect of lower back biomechanics. For each set of TMFs, the segmental kinematics of lumbar spine was estimated for a single static trunk flexed posture involving, respectively, 40° and 10° of thoracic and pelvic rotations. Minimum changes in the angular and translational deformations of a motion segment with alterations in TMFs ranged from 0° to 0.7° and 0 mm to 0.04 mm, respectively. Maximum changes in the angular and translational deformations of a motion segment with alterations in TMFs ranged from 2.4° to 7.6° and 0.11 mm to 0.39 mm, respectively. The differences in kinematics of lumbar segments between each combination of two sets of TMFs in 97% of cases for angular deformation and 55% of cases for translational deformation were within the reported accuracy of current imaging techniques. Therefore, it might be possible to use image-based kinematics of lumbar segments along with computational modeling for personalized assessment of TMFs.
Macroscopic anisotropic bone material properties in children with severe osteogenesis imperfecta.
Venue: Journal of biomechanics

Children with severe osteogenesis imperfecta (OI) typically experience numerous fractures and progressive skeletal deformities over their lifetime. Recent studies proposed finite element models to assess fracture risk and guide clinicians in determining appropriate intervention in children with OI, but lack of appropriate material property inputs remains a challenge. This study aimed to characterize macroscopic anisotropic cortical bone material properties and investigate relationships with bone density measures in children with severe OI. Specimens were obtained from tibial or femoral shafts of nine children with severe OI and five controls. The specimens were cut into beams, characterized in bending, and imaged by synchrotron radiation X-ray micro-computed tomography. Longitudinal modulus of elasticity, yield strength, and bending strength were 32-65% lower in the OI group (p<0.001). Yield strain did not differ between groups (p≥0.197). In both groups, modulus and strength were lower in the transverse direction (p≤0.009), but anisotropy was less pronounced in the OI group. Intracortical vascular porosity was almost six times higher in the OI group (p<0.001), but no differences were observed in osteocyte lacunar porosity between the groups (p=0.086). Volumetric bone mineral density was lower in the OI group (p<0.001), but volumetric tissue mineral density was not (p=0.770). Longitudinal OI bone modulus and strength were correlated with volumetric bone mineral density (p≤0.024) but not volumetric tissue mineral density (p≥0.099). Results indicate that cortical bone in children with severe OI yields at the same strain as normal bone, and that their decreased bone material strength is associated with reduced volumetric bone mineral density. These results will enable the advancement of fracture risk assessment capability in children with severe OI.
In vivo fascicle length measurements via B-mode ultrasound imaging with single vs dual transducer arrangements.
Venue: Journal of biomechanics

Ultrasonography is a useful technique to study muscle contractions in vivo, however larger muscles like vastus lateralis may be difficult to visualise with smaller, commonly used transducers. Fascicle length is often estimated using linear trigonometry to extrapolate fascicle length to regions where the fascicle is not visible. However, this approach has not been compared to measurements made with a larger field of view for dynamic muscle contractions. Here we compared two different single-transducer extrapolation methods to measure VL muscle fascicle length to a direct measurement made using two synchronised, in-series transducers. The first method used pennation angle and muscle thickness to extrapolate fascicle length outside the image (extrapolate method). The second method determined fascicle length based on the extrapolated intercept between a fascicle and the aponeurosis (intercept method). Nine participants performed maximal effort, isometric, knee extension contractions on a dynamometer at 10° increments from 50 to 100° of knee flexion. Fascicle length and torque were simultaneously recorded for offline analysis. The dual transducer method showed similar patterns of fascicle length change (overall mean coefficient of multiple correlation was 0.76 and 0.71 compared to extrapolate and intercept methods respectively), but reached different absolute lengths during the contractions. This had the effect of producing force-length curves of the same shape, but each curve was shifted in terms of absolute length. We concluded that dual transducers are beneficial for studies that examine absolute fascicle lengths, whereas either of the single transducer methods may produce similar results for normalised length changes, and repeated measures experimental designs.