nils-herrmann/clean_scirepeval huggingface dataset:
Red blood cell motion and deformation in a curved microvessel.
Venue: Journal of biomechanics

The flow of cells through curved vessels is often encountered in various biomedical and bioengineering applications, such as red blood cells (RBCs) passing through the curved arteries in circulation, and cells sorting through a shear-induced migration in a curved channels. Most of past numerical studies focused on the cell deformation in small straight microvessels, or on the flow pattern in large curved vessels without considering the cell deformation. However, there have been few attempts to study the cell deformation and the associated flow pattern in a curved microvessel. In this work, a particle-based method, smoothed dissipative particle dynamics (SDPD), is used to simulate the motion and deformation of a RBC in a curved microvessel of diameter comparable to the RBC diameter. The emphasis is on the effects of the curvature, the type and the size of the curved microvessel on the RBC deformation and the flow pattern. The simulation results show that a small curved shape of the microvessel has negligible effect on the RBC behavior and the flow pattern which are similar to those in a straight microvessel. When the microvessel is high in curvature, the secondary flow comes into being with a pair of Dean vortices, and the velocity profile of the primary flow is skewed toward the inner wall of the microvessel. The RBC also loses the axisymmetric deformation, and it is stretched first and then shrinks when passing through the curved part of the microvessel with the large curvature. It is also found that a pair of Dean vortices arise only under the condition of De>1 (De is the Dean number, a ratio of centrifugal to viscous competition). The Dean vortices are more easily observed in the larger or more curved microvessels. Finally, it is observed that the velocity profile of primary flow is skewed toward the inner wall of curved microvessel, i.e., the fluid close to the inner wall flows faster than that close to the outer wall. This is contrary to the common sense in large curved vessels. This velocity skewness was found to depend on the curvature of the microvessel, as well as the viscous and inertial forces.
Wearable technology for spine movement assessment: A systematic review
Venue: Journal of biomechanics

Continuous monitoring of spine movement function could enhance our understanding of low back pain development. Wearable technologies have gained popularity as promising alternative to laboratory systems in allowing ambulatory movement analysis. This paper aims to review the state of art of current use of wearable technology to assess spine kinematics and kinetics. Four electronic databases and reference lists of relevant articles were searched to find studies employing wearable technologies to assess the spine in adults performing dynamic movements. Two reviewers independently identified relevant papers. Customised data extraction and quality appraisal form were developed to extrapolate key details and identify risk of biases of each study. Twenty-two articles were retrieved that met the inclusion criteria: 12 were deemed of medium quality (score 33.4–66.7%), and 10 of high quality (score >66.8%). The majority of articles (19/22) reported validation type studies. Only 6 reported data collection in real-life environments. Multiple sensors type were used: electrogoniometers (3/22), strain gauges based sensors (3/22), textile piezoresistive sensor (1/22) and accelerometers often used with gyroscopes and magnetometers (15/22). Two sensors units were mainly used and placing was commonly reported on the spine lumbar and sacral regions. The sensors were often wired to data transmitter/logger resulting in cumbersome systems. Outcomes were mostly reported relative to the lumbar segment and in the sagittal plane, including angles, range of motion, angular velocity, joint moments and forces. This review demonstrates the applicability of wearable technology to assess the spine, although this technique is still at an early stage of development.
Intersegmental kinetics significantly impact mapping from finger musculotendon forces to fingertip forces.
Venue: Journal of biomechanics

Predicting the fingertip force vector resulting from excitation of a given muscle remains a challenging but essential task in finger biomechanical modeling. While the conversion of musculotendon force to fingertip force can significantly be affected by finger posture, current techniques utilizing geometric moment arms may not capture such complex postural effects. Here, we attempted to elucidate the postural effects on the mapping between musculotendon force and fingertip force through in vitro experiments. Computer-controlled tendon loading was implemented on the 7 index finger musculotendons of 5 fresh-frozen cadaveric hands across different postures. The resulting fingertip forces/moments were used to compute the effective static moment arm (ESMA), relating tendon force to joint torque, at each joint. The ESMAs were subsequently modeled in three different manners: independent of joint angle; dependent only upon the corresponding joint angle; or dependent upon all joint angles. We found that, for the reconstruction of the fingertip force vector, the multi-joint ESMA model yielded the best outcome, both in terms of direction and magnitude of the vector (mean reconstruction error <4° in direction and <2% in the magnitude), which indicates that intersegmental force transmission through a joint is affected by the posture of neighboring joints. Interestingly, the ESMA model that considers geometric changes of individual joints, the standard model used in biomechanical stimulations, often yielded worse reconstruction results than the simple constant-value ESMA model. Our results emphasize the importance of accurate description of the multi-joint dependency of the conversion of tendon force to joint moment for proper prediction of fingertip force direction.
Does elevated wall tension cause aortic aneurysm rupture? Investigation using a subject-specific heterogeneous model.
Venue: Journal of biomechanics

To investigate whether peak wall tension in abdominal aortic aneurysm occurs at the site of rupture to test for a causative relationship.



Four ruptured and nine unruptured AAA were harvested whole from cadavers, followed by regional measurements of wall thickness, elastic parameters and failure tension. Finite element models were developed with subject-specific load-free AAA morphology and heterogeneous properties interpolated using a geodesic distance weighted approach from the measurements. The wall tension under uniform pressure and tension to failure tension ratio as an index of susceptibility to rupture were computed. As a secondary aim, the peak wall tension using this heterogeneous model approach was compared to the traditional homogeneous model approach in order to evaluate the reliability of the latter.



The average peak wall tension in the ruptured group was 43% higher than in the unruptured group without statistical significance even though it was 54% larger in diameter. The site of peak wall tension was in the vicinity of the site of rupture in two ruptured AAA. The peak tension did not breach failure tension at the rupture site in any of the AAA. The traditional population-wide homogeneous model approach overestimated peak wall tension by just 3% compared to the subject-specific heterogeneous model approach.



We failed to find adequate evidence of a causative relationship between peak wall tension and AAA rupture. The findings are not conclusive owing to study limitations such as ignoring intraluminal thrombus, sparse distribution of specimens procured and small study population.
Determinants of gait stability while walking on a treadmill: A machine learning approach.
Venue: Journal of biomechanics

Dynamic balance in human locomotion can be assessed through the local dynamic stability (LDS) method. Whereas gait LDS has been used successfully in many settings and applications, little is known about its sensitivity to individual characteristics of healthy adults. Therefore, we reanalyzed a large dataset of accelerometric data measured for 100 healthy adults from 20 to 70 years of age performing 10 min treadmill walking. We sought to assess the extent to which the variations of age, body mass and height, sex, and preferred walking speed (PWS) could influence gait LDS. The random forest (RF) and multiple adaptive regression splines (MARS) algorithms were selected for their good bias-variance tradeoff and their capabilities to handle nonlinear associations. First, through variable importance measure (VIM), we used RF to evaluate which individual characteristics had the highest influence on gait LDS. Second, we used MARS to detect potential interactions among individual characteristics that may influence LDS. The VIM and MARS results indicated that PWS and age correlated with LDS, whereas no associations were found for sex, body height, and body mass. Further, the MARS model detected an age by PWS interaction: on one hand, at high PWS, gait stability is constant across age while, on the other hand, at low PWS, gait instability increases substantially with age. We conclude that it is advisable to consider the participants' age as well as their PWS to avoid potential biases in evaluating dynamic balance through LDS.
Analysis of non-Newtonian effects within an aorta-iliac bifurcation region.
Venue: Journal of biomechanics

The geometry of the arteries at or near arterial bifurcation influences the blood flow field, which is an important factor affecting arteriogenesis. The blood can act sometimes as a non-Newtonian fluid. However, many studies have argued that for large and medium arteries, the blood flow can be considered to be Newtonian. In this work a comprehensive investigation of non-Newtonian effects on the blood fluid dynamic behavior in an aorta-iliac bifurcation is presented. The aorta-iliac geometry is reconstructed with references to the values reported in Shah et al. (1978); the 3D geometrical model consists of three filleted cylinders of different diameters. Governing equations with the appropriate boundary conditions are solved with a finite-element code. Different rheological models are used for the blood flow through the lumen and detailed comparisons are presented for the aorta-iliac bifurcation. Results are presented in terms of the velocity profiles in the bifurcation zone and Wall Shear Stress (WSS) for different sides of the bifurcation both for male and female geometries, showing that the Newtonian fluid assumption can be made without any particular loss in terms of accuracy with respect to the other more complex rheological models.
A strain-based finite element model for calcification progression in aortic valves.
Venue: Journal of biomechanics

Calcific aortic valve disease (CAVD) is a serious disease affecting the aging population. A complex interaction between biochemicals, cells, and mechanical cues affects CAVD initiation and progression. In this study, motivated by the progression of calcification in regions of high strain, we developed a finite element method (FEM) based spatial calcification progression model. Several cardiac cycles of transient structural FEM simulations were simulated. After each simulation cycle, calcium deposition was placed in regions of high circumferential strain. Our results show the radial expansion of calcification as spokes starting from the attachment region, agreeing very well with the reported clinical data.
Dual-task and anticipation impact lower limb biomechanics during a single-leg cut with body borne load.
Venue: Journal of biomechanics

This study quantified how a dual cognitive task impacts lower limb biomechanics during anticipated and unanticipated single-leg cuts with body borne load. Twenty-four males performed anticipated and unanticipated cuts with and without a dual cognitive task with three load conditions: no load (∼6 kg), medium load (15% of BW), and heavy load (30% of BW). Lower limb biomechanics were submitted to a repeated measures linear mixed model to test the main and interaction effects of load, anticipation, and dual task. With body borne load, participants increased peak stance (PS) hip flexion (p = .004) and hip internal rotation (p = .001) angle, and PS hip flexion (p = .001) and internal rotation (p = .018), and knee flexion (p = .016) and abduction (p = .001) moments. With the dual task, participants decreased PS knee flexion angle (p < .001) and hip flexion moment (p = .027), and increased PS knee external rotation angle (p = .034). During the unanticipated cut, participants increased PS hip (p = .040) and knee flexion angle (p < .001), and decreased PS hip adduction (p = .001), and knee abduction (p = .005) and external rotation (p = .026) moments. Adding body borne load produces lower limb biomechanical adaptations thought to increase risk of musculoskeletal injury, but neither anticipation nor dual task exaggerated those biomechanical adaptations. With a dual task, participants adopted biomechanics known to increase injury risk; whereas, participants used lower limb biomechanics thought to decrease injury risk during unanticipated cuts.
Segmental variations in facet joint translations during in vivo lumbar extension.
Venue: Journal of biomechanics

The lumbar facet joint (FJ) is often associated with pathogenesis in the spine, but quantification of normal FJ motion remains limited to in vitro studies or static imaging of non-functional poses. The purpose of this study was to quantify lumbar FJ kinematics in healthy individuals during functional activity with dynamic stereo radiography (DSX) imaging. Ten asymptomatic participants lifted three known weights starting from a trunk-flexed (∼75°) position to an upright position while being imaged within the DSX system. High resolution computed tomography (CT) scan-derived 3D models of their lumbar vertebrae (L2-S1) were registered to the biplane 2D radiographs using a markerless model-based tracking technique providing instantaneous 3D vertebral kinematics throughout the lifting tasks. Effects of segment level and weight lifted were assessed using mixed-effect repeated measures ANOVA. Superior-inferior (SI) translation dominated FJ translation, with L5S1 showing significantly less translation magnitudes (Median (Md) = 3.5 mm, p < 0.0001) than L2L3, L3L4, and L4L5 segments (Md = 5.9 mm, 6.3 mm and 6.6 mm respectively). Linear regression-based slopes of continuous facet translations revealed strong linearity for SI translation (r2 > 0.94), reasonably high linearity for sideways sliding (Z-) (r2 > 0.8), but much less linearity for facet gap change (X-) (r2 ∼ 0.5). Caudal segments (L4-S1), particularly L5S1, displayed greater coupling compared to cranial (L2-L4) segments, revealing distinct differences overall in FJ translation trends at L5S1. No significant effect of weight lifted on FJ translations was detected. The study presents a hitherto unavailable and highly precise baseline dataset of facet translations measured during a functional, dynamic lifting task.
Chiari malformation may increase perivascular cerebrospinal fluid flow into the spinal cord: A subject-specific computational modelling study.
Venue: Journal of biomechanics

Syringomyelia is associated with Chiari I malformation, although the mechanistic link is unclear. Studies have suggested that cerebrospinal fluid enters the spinal cord via the perivascular spaces, and that changes in the timing of the subarachnoid pressures may increase flow into the spinal cord. This study aims to determine how Chiari malformation and syringomyelia alter the subarachnoid space pressures and hence perivascular flow. Subject-specific models of healthy controls (N = 9), Chiari patients with (N = 7) and without (N = 8) syringomyelia, were developed from magnetic resonance imaging (MRI), to simulate the subarachnoid pressures. These pressures were input to an idealised model of the perivascular space to evaluate potential differences in perivascular flow. Peak pressures in Chiari patients without a syrinx were higher than in controls (46% increase; p = .029) and arrived earlier in the cardiac cycle than both controls (2.58% earlier; p = .045) and syrinx patients (2.85% earlier; p = .045). The perivascular model predicted Chiari patients without a syrinx would have the greatest flow into the cord (p < .05) if the arterial pulse delay was between 4 and 10% of the cardiac cycle. Using phase-contrast MRI the mean arterial delay for all subjects was similar, and was estimated as 4.7 ± 0.2%. The perivascular pumping rate showed a strong positive correlation (RAdj2=0.85; p < .0001) with extended periods of high pressure that arrived earlier in the cardiac cycle, suggesting these pressure characteristics may play a role in syrinx development.
Fractal fluctuations in spatiotemporal variables when walking on a self-paced treadmill.
Venue: Journal of biomechanics

This study investigated the fractal dynamic properties of stride time (ST), stride length (SL) and stride speed (SS) during walking on a self-paced treadmill (STM) in which the belt speed is automatically controlled by the walking speed. Twelve healthy young subjects participated in the study. The subjects walked at their preferred walking speed under four conditions: STM, STM with a metronome (STM+met), fixed-speed (conventional) treadmill (FTM), and FTM with a metronome (FTM+met). To compare the fractal dynamics between conditions, the mean, variability, and fractal dynamics of ST, SL, and SS were compared. Moreover, the relationship among the variables was examined under each walking condition using three types of surrogates. The mean values of all variables did not differ between the two treadmills, and the variability of all variables was generally larger for STM than for FTM. The use of a metronome resulted in a decrease in variability in ST and SS for all conditions. The fractal dynamic characteristics of SS were maintained with STM, in contrast to FTM, and only the fractal dynamic characteristics of ST disappeared when using a metronome. In addition, the fractal dynamic patterns of the cross-correlated surrogate results were identical to those of all variables for the two treadmills. In terms of the fractal dynamic properties, STM walking was generally closer to overground walking than FTM walking. Although further research is needed, the present results will be useful in research on gait fractal dynamics and rehabilitation.
Osmosis and viscoelasticity both contribute to time-dependent behaviour of the intervertebral disc under compressive load: A caprine in vitro study.
Venue: Journal of biomechanics

The mechanical behaviour of the intervertebral disc highly depends on the content and transport of interstitial fluid. It is unknown, however, to what extent the time-dependent behaviour can be attributed to osmosis. Here we investigate the effect of both mechanical and osmotic loading on water content, nucleus pressure and disc height. Eight goat intervertebral discs, immersed in physiological saline, were subjected to a compressive force with a pressure needle inserted in the nucleus. The loading protocol was: 10 N (6 h); 150 N (42 h); 10 N (24 h). Half-way the 150 N-phase (24 h), we eliminated the osmotic gradient by adding 26% poly-ethylene glycol to the surrounding fluid. For 62 additional discs, we determined the water content of both nucleus and annulus after 6, 24, 48, or 72 h. The compressive load was initially counterbalanced by the hydrostatic pressure in the nucleus. The load forced 4.3% of the water out of the nucleus, which reduced nucleus pressure by 44(±6)%. Reduction of the osmotic gradient disturbed the equilibrium disc height, and a significant loss of annulus water content was found. Remarkably, pressure and water content of the nucleus pulposus remained unchanged. This shows that annulus water content is important in the response to axial loading. After unloading, in the absence of an osmotic gradient, there was substantial viscoelastic recovery of 53(±11)% of the disc height, without a change in water content. However, for restoration of the nucleus pressure and for full restoration of disc height, restoration of the osmotic gradient was needed.
Trunk response and stability in standing under sagittal-symmetric pull-push forces at different orientations, elevations and magnitudes.
Venue: Journal of biomechanics

To reduce lifting and associated low back injuries, manual material handling operations often involve pulling-pushing of carts at different weights, orientations, and heights. The loads on spine and risk of injury however need to be investigated. The aim of this study was to evaluate muscle forces, spinal loads and trunk stability in pull-push tasks in sagittal-symmetric, static upright standing posture. Three hand-held load magnitudes (80, 120 and 160 N) at four elevations (0, 20, 40 and 60 cm to the L5-S1) and 24 force directions covering all pull/push orientations were considered. For this purpose, a musculoskeletal finite element model with kinematics measured earlier were used. Results demonstrated that peak spinal forces occur under inclined pull (lift) at upper elevations but inclined push at the lowermost one. Minimal spinal loads, on the other hand, occurred at and around vertical pull directions. Overall, spinal forces closely followed variations in the net external moment of pull-push forces at the L5-S1. Local lumbar muscles were most active in pulls while global extensor muscles in lifts. The trunk stability margin decreased with load elevation except at and around horizontal push; it peaked under pulls and reached minimum at vertical lifts. It also increased with antagonist activity in muscles and intra-abdominal pressure. Results provide insight into the marked effects of variation in the load orientation and elevation on muscle forces, spinal loads and trunk stability and hence offer help in rehabilitation, performance enhancement training and design of safer workplaces.
Chest response assessment of post-mortem swine under blast loadings.
Venue: Journal of biomechanics

To better protect soldiers from blast threat, that principally affect air-filled organs such a lung, it is necessary to develop an adapted injury criterion and, prior to this, to evaluate the response of a biological model against that threat. The objective of this study is to provide some robust data to quantify the chest response of post-mortem swine under blast loadings. 7 post-mortem swine (54.5 ± 2.6 kg), placed side-on to the threat and against the ground, were exposed to 5 shock-waves of increasing intensities. Their thorax were instrumented with a piezo-resistive pressure sensor, an accelerometer directly exposed to the shock-wave and a target was mounted on the latter in order to track the chest wall displacement. For incident impulses ranging from 47 kPa ms±2% to 173 kPa ms ±6%, the measured maximum of linear chest wall acceleration (Γmax) goes from 5800 m/s2 ±16% to 41,000 m/s2 ± 8%, with a duration of 0.8 ms. Chest wall displacements ranging from 5 mm ± 20% to 20 mm ± 15%, with a duration of 9 ms, are reached. These reproducible data were used to find simple relations (linear, 2nd and 3rd order polynomials) between the kinematic parameters (plus the viscous criterion) and the incident and reflected impulses. Correlating the new reproducible data with the prediction from the Bowen curves showed a lung injury threshold in terms of Γmax similar to that of Cooper (10,000 m/s2). However, the limits defined for the viscous criterion in the automobile field and for non-lethal weapons seems not adapted for the blast threat.