nils-herrmann/clean_scirepeval huggingface dataset:
Influence of loading condition and anatomical location on human cortical bone linear micro-cracks.
Venue: Journal of biomechanics

Human cortical bone fracture toughness depends on the anatomical locations under quasi-static loading. Recent results also showed that under fall-like loading, cortical bone fracture toughness is similar at different anatomical locations in the same donor. While cortical bone toughening mechanisms are known to be dependent on the tissue architecture under quasi-static loading, the fracture mechanisms during a fall are less studied. In the current study, the structural parameters of eight paired femoral diaphyses, femoral necks and radial diaphyses were mechanically tested under quasi-static and fall-like loading conditions (female donors, 70 ± 14 y.o., [50-91 y.o.]). Synchrotron radiation micro-CT imaging was used to quantify the amount of micro-cracks formed during loading. The volume fraction of these micro-cracks was significantly higher within the specimens loaded under a quasi-static condition than under a loading representative of a fall. Under fall-like loading, there was no difference in crack volume fraction between the different paired anatomical locations. This result shows that the micro-cracking toughening mechanism depends both on the anatomical location and on the loading condition.
Myocardial bridging and endothelial dysfunction - Computational fluid dynamics study.
Venue: Journal of biomechanics

Myocardial bridging (MB) is associated with endothelial dysfunction in patients with angina and non-obstructive coronary artery disease. This study aims to determine if there is a link between abnormal blood flow patterns and endothelial dysfunction in patients with MB. Ten patients with MB in their left anterior descending (LAD) artery were selected, 5 of whom had endothelial dysfunction and 5 had no endothelial dysfunction based on their response to acetylcholine. Similarly, 10 patients without MB in their LAD, 5 of whom had endothelial dysfunction and 5 of whom had no endothelial dysfunction, were studied as a control group. Transient computational fluid dynamics simulations were performed to derive wall shear stress (WSS) over the entire vessel including proximal, middle and distal segments. Patients with MB and endothelial dysfunction had lower WSS in the proximal LAD and greater WSS in the mid-LAD than patients with MB but without endothelial dysfunction. When comparing patients with endothelial dysfunction, those with MB had significantly lower shear stress in the proximal LAD (0.32 ± 0.14 Pa (with MB) vs 0.71 ± 0.38 Pa (without MB), p = 0.01) and greater shear stress in the mid-LAD (2.81 ± 1.20 Pa (with MB) vs 1.66 ± 0.31 Pa (without MB), p = 0.014) than patients without MB. Our findings demonstrated that the presence of MB significantly contributes to low WSS and endothelial dysfunction relationship.
Development of a computational fluid dynamics model for mucociliary clearance in the nasal cavity.
Venue: Journal of biomechanics

Intranasal drug delivery has attracted significant attention because of the opportunity to deliver systemic drugs directly to the blood stream. However, the mucociliary clearance poses a challenge in gaining high efficacy of intranasal drug delivery because cilia continuously carry the mucus blanket towards the laryngeal region. To better understand mucus flow behaviour on the human nasal cavity wall, we present computational model development, and evaluation of mucus motion on a realistic nasal cavity model reconstructed from CT-scans. The model development involved two approaches based on the actual nasal cavity geometry namely: (i) unwrapped-surface model in 2D domain and (ii) 3D-shell model. Conservation equations of fluid motion were applied to the domains, where a mucus production source term was used to initiate the mucus motion. The analysis included mucus flow patterns, virtual saccharin tests and quantitative velocity magnitude analysis, which demonstrated that the 3D-shell model results provided better agreement with experimental data. The unwrapped-surface model also suffered from mesh-deformations during the unwrapping stage and this led to higher mucus velocity compared to experimental data. Therefore, the 3D-shell model was recommended for future mucus flow simulations. As a first step towards mucus motion modelling this study provides important information that accurately simulates a mucus velocity field on a human nasal cavity wall, for assessment of toxicology and efficacy of intranasal drug delivery.
A fully nonlinear viscohyperelastic model for the brain tissue applicable to dynamic rates.
Venue: Journal of biomechanics

Understanding the mechanical response of the brain to external loadings is of critical importance in investigating the pathological conditions of this tissue during injurious conditions. Such injurious loadings may occur at high rates, for example among others, during road traffic or sport accidents, falls, or due to explosions. Hence, investigating the injury mechanism and design of protective devices for the brain requires constitutive modeling of this tissue at such rates. Accordingly, this paper is aimed at critically investigating the physical background for viscohyperelastic modeling of the brain tissue with scrutinizing the elastic fields pertinent to large, time dependent deformations, and developing a fully nonlinear multimode Maxwell model that can mathematically explain such deformations. The proposed model can be calibrated using the simple monotonic uniaxial deformation of the sample extracted from the tissue, and does not require additional information from relaxation or creep experiments. The performance of the proposed model is examined using the experimental results of two different studies, which reveals a desirable agreement. The usefulness, limitations, and future developments of the proposed model are discussed in this paper.
The relationship between lumbopelvic flexibility and sitting posture in adult women.
Venue: Journal of biomechanics

Clinical observations have suggested that limited hamstring flexibility may be associated with sagittal spinal curvatures in spine flexed postures. Thus, limited hamstring flexibility may be related to large amounts of spine flexion in "slumped" sitting postures which could contribute to low back pain and injury. The aim of this study was to determine if hamstring and pelvic flexibility are associated with flexed sitting postures using a backless office chair. Forty-one healthy female adults aged 18-69  years were recruited. Subjects performed the Sit-and-Reach test to determine maximum flexibility values and lumbar and pelvic angles were measured with accelerometers. Participants then completed a standardized typing task for a 10-minute sitting trial at an ergonomically adjusted workstation. The results showed no association between hamstring flexibility and seated lumbar spine and pelvic angles (p = 0.999, η2 = 0.000; p = 0.901, η2 = 0.006). Greater pelvic flexibility was associated with a more upright lumbar sitting posture (p = 0.023; η2 = 0.132) but with no specific pelvic sitting posture (p = 0.660; η2 = 0.005). Different movement strategies during the Sit-and-Reach test were detected: all participants moved through their lumbar spine; but only those with 'excellent' flexibility also used their pelvis. Individuals in the 'excellent' flexibility group were significantly shorter than those with 'poor' and 'good' flexibility (p = 0.020; η2 = 0.190). In conclusion, hamstring flexibility does not influence sitting posture but pelvic flexibility does. Other factors such as acetabulofemoral joint limitations, consciousness of posture, or the seat itself may also influence sitting posture. Different movement strategies as well as height appear to contribute to the Sit-and-Reach test which should be researched further.
A robust method to estimate the largest Lyapunov exponent of noisy signals: A revision to the Rosenstein's algorithm.
Venue: Journal of biomechanics

This study proposed a revision to the Rosenstein's method of numerical calculation of the largest Lyapunov exponent (LyE) to make it more robust to noise. To this aim, the effect of increasing number of initial neighboring points on the LyE value was investigated and compared to values obtained by filtering the time series. Both simulated (Lorenz and passive dynamic walker) and experimental (human walking) time series were used to calculate the LyE. The number of initial neighbors used to calculate LyE for all time series was 1 (the original Rosenstein's method), 2, 3, 4, 5, 10, 15, 20, 25, and 30 data points. The results demonstrated that the LyE graph reached a plateau at the 15-point neighboring condition implying that the LyE values calculated using at least 15 neighboring points were consistent. The proposed method could be used to calculate more consistent LyE values in experimental time series acquired from biological systems where noise is omnipresent.
Postural differences in shoulder dynamics during pushing and pulling.
Venue: Journal of biomechanics

Assessments of shoulder dynamics (e.g. the inertial, viscous, and stiffness properties of the joint) can provide important insights into the stability of the joint at rest and during volitional contraction. The purpose of this study was to investigate how arm posture influences shoulder dynamics while generating pushing or pulling torques in the horizontal plane. Sixteen healthy participants were examined in seven postures encompassing a large workspace of the shoulder. At each posture, the participant's shoulder was rapidly perturbed while measuring the resultant change in shoulder torque about the glenohumeral axis. Participants were examined both at rest and while producing horizontal flexion and extension torques scaled to 15% of a maximum voluntary contraction. Shoulder stiffness, viscosity, and damping ratio were estimated using impedance-based matching, and changes in these outcome measures with torque level, elevation angle, and plane of elevation angle were explored with a linear mixed effects model. Shoulder stiffness was found to decrease with increasing elevation angles (p < 0.001) without subsequent changes in viscosity, leading to a greater damping ratios at higher elevation angles (p < 0.001). Shoulder stiffness, viscosity, and damping ratio (all p < 0.05) were all found to significantly increase as the plane of elevation of the arm was increased. The relationship between the viscosity, stiffness and the damping ratio of the shoulder is one that the central nervous system must regulate in order to maintain stability, protect against injury, and control the shoulder joint as the inertial and muscle contributions change across different arm postures.
New considerations for collecting biomechanical data using wearable sensors: Number of level runs to define a stable running pattern with a single IMU.
Venue: Journal of biomechanics

Wearable technology can be used to quantify running biomechanical patterns in a runner's natural environment, however, changes in external factors during outdoor running may influence a runner's typical gait pattern. Therefore, the purpose of this study was to determine how many runs are needed to define a stable or typical running pattern. Six biomechanical variables were recorded using a single wearable sensor placed on the lower back during ten outdoor runs for twelve runners. Univariate and multivariate distributions were created and based on the probability density function, the percent of similar data points (within 95%) from each unique run for the same runner were determined. Stability was defined when the addition of data from a new run resulted in less than a 5% change in the probability density function. To cross-validate, the percent of similar data points at stability was compared between the same and different runners using a repeated-measures MANOVA (Bonferroni-corrected α = 0.007). The maximum number of runs needed to reach stability for univariate and multivariate analyses was four and five, respectively. There was a significant overall effect on similar data points between the same and different runners (p = 0.001), with a greater percent of similar data points for the same runner compared to other runners (p < 0.007). Based on biomechanical data collected using a single wearable sensor placed on the lower back, this is the first study to show that four (univariate) to five (multivariate) runs are needed to establish a stable running pattern in real-world settings.
Subject-specific and group-based running pattern classification using a single wearable sensor.
Venue: Journal of biomechanics

The objective of this study was to determine whether subject-specific or group-based models provided better classification accuracy to identify changes in biomechanical running gait patterns across different inclination conditions. The classification process was based on measurements from a single wearable sensor using a total of 41,780 strides from eleven recreational runners while running in real-world and uncontrolled environment. Biomechanical variables included pelvic drop, ground contact time, braking, vertical oscillation of pelvis, pelvic rotation, and cadence were recorded during running on three inclination grades: downhill, -2° to -7°; level, -0.2° to +0.2°; and uphill, +2° to +7°. An ensemble and non-linear machine learning algorithm, random forest (RF), was used to classify inclination condition and determine the importance of each of the biomechanical variables. Classification accuracy was determined for subject-specific and group-based RF models. The mean classification accuracy of all subject-specific RF models was 86.29%, while group-based classification accuracy was 76.17%. Braking was identified as the most important variable for all the runners using the group-based model and for most of the runners based on a subject-specific models. In addition, individual runners used different strategies across different inclination conditions and the ranked order of variable importance was unique for each runner. These results demonstrate that subject-specific models can better characterize changes in gait biomechanical patterns compared to a more traditional group-based approach.
An in vitro experimental study on the relationship between pulsatile tinnitus and the dehiscence/thinness of sigmoid sinus cortical plate.
Venue: Journal of biomechanics

Pulsatile tinnitus (PT), characterized as pulse-synchronous, is generally objective. Sigmoid sinus (SS) venous sound is widely suggested to be a possible sound source of PT. The dehiscence and thinness of SS cortical plate (CP) was commonly reported as PT pathology in previous studies, but lack quantitative or biomechanical analysis. In this study, it was aimed to quantify the relationship between venous sound and CP dehiscence/thinness using in vitro experiment. The in vitro models of SS and CP were established based on 3D-printing, with the developed pulsatile venous flow in the SS model. The generated sound signal and the vibration response at the dehiscent/thinned area were analyzed. The sound signal generated in the normal-sized dehiscence model was pulse-synchronous within 100--400 Hz, which had similar acoustic characteristics as the clinical PT sounds. It was concluded that the pulsatile venous sound is produced at TS-SS junction in case of CP dehiscence. The CP, even a thinned one can effectively diminish the venous sound and sound-generating pulsatile vibration at TS-SS junction. The CP dehiscence would induce pulse-synchronous and high pressure venous sound, as well as pulse-synchronous vibration above 20 Hz, regardless of the dehiscence size. On the contrary, the CP thinness would not induce obvious venous sound or pulsatile vibration above 20 Hz.
Marker cluster rigidity in a multi-segment foot model.
Venue: Journal of biomechanics

Multi-segment foot models (MSFM) are used in gait analysis for the diagnosis and planning of treatment for patients with foot deformities. Like other biomechanical models, MSFMs represent the leg and foot as a series of linked rigid segments, but such a simplification may not be appropriate, particularly for the flexible forefoot. This study investigated the appropriateness of the rigid body assumption on marker clusters used to define the individual segments (tibia, hindfoot, forefoot) of a widely-used MSFM. Rigidity of the marker clusters was quantified using the rigid body error (σRBE) calculated for each frame of a representative gait cycle for 64 normal healthy adults who underwent gait analysis. σRBE is a measure of how well the tracking marker configuration at each frame compares to the arrangement of the same markers in a reference pose. As expected, the patterns of deformation of the three marker clusters differed over the gait cycle. The hindfoot cluster remained relatively undeformed in comparison to the forefoot and tibia clusters. The largest deformations of the forefoot cluster occurred near the beginning and end of the stance phase. The tibia cluster deformed throughout the entire gait cycle, with a pattern similar to that of a typical knee flexion angle graph. The results raise questions about the appropriateness of the rigid-body assumption when applied to MSFMs, particularly in the forefoot region.
Site- and sex-differences in morphological and mechanical properties of the plantar fascia: A supersonic shear imaging study.
Venue: Journal of biomechanics

The purpose of this study was to investigate the site- and sex-differences in the morphological and mechanical properties of the plantar fascia (PF) in humans. The thickness and shear wave velocity (SWV) of PF at five different sites between the medial calcaneal tubercle and the second toe were measured for 40 healthy young participants (20 males and 20 females) using supersonic shear imaging (SSI). The thickness and SWV measurements were highly repeatable (ICC ≥ 0.93). The proximal sites of PF around the calcaneal attachment were significantly thicker and stiffer (higher SWV values) than the middle and distal sites (p < 0.05). In addition, females had significantly thinner PF in proximal and middle sites than males, while being significantly stiffer in regardless of the sites, compared with males (p < 0.05). The results of the present study partly support previous findings on the site- and sex-differences in PF morphology, and further reveal inhomogeneity and sex-specificity of PF stiffness. The present study widely opens the possibility of evaluating PF functions in vivo.
Modeling mitral valve stenosis: A parametric study on the stenosis severity level.
Venue: Journal of biomechanics

New computational techniques providing more accurate representation of human heart pathologies could help uncovering relevant physical phenomena and improve the outcome of medical therapies. In this framework, the present work describes an efficient computational model for the evaluation of the ventricular flow alteration in presence of mitral valve stenosis. The model is based on the direct numerical simulation of the Navier-Stokes equations two-way coupled with a structural solver for the left ventricle and mitral valve dynamics. The presence of mitral valve stenosis is mimicked by a single-parameter constraint acting on the kinematics of the mitral leaflets. Four different degrees of mitral valve stenosis are considered focusing on the hemodynamic alterations occurring in pathologic conditions. The mitral jet, generated during diastole, is seen to shrink and strengthen when the stenosis gets more severe. As a consequence, the kinetic energy of the flow, the tissues shear stresses, the transvalvular pressure drop and mitral regurgitation increase. It results that, as the stenosis severity level increases, the geometric and effective orifice areas decrease up to 50% with respect the normal case due to the reduced leaflets mobility and stronger blood acceleration during the diastolic phase. The modified intraventricular hemodynamics is also related to a stronger pressure gradient that, for severe stenosis, can be more than ten times larger than the healthy valve case. These computational results are fully consistent with the available clinical literature and open the way to the virtual assessment of surgical procedures and to the evaluation of prosthetic devices.
On intrinsic equivalences of the finite helical axis, the instantaneous helical axis, and the SARA approach. A mathematical perspective.
Venue: Journal of biomechanics

Accurate determination of joint axes is essential for understanding musculoskeletal function. Whilst numerous algorithms to compute such axes exist, the conditions under which each of the methods performs best remain largely unknown. Typically, algorithms are evaluated for specific conditions only limiting the external validity of conclusions regarding their performance. We derive exact mathematical relationships between three commonly used algorithms for computing joint axes from motion data: finite helical axes (FHA), instantaneous helical axes (IHA) and SARA (symmetrical axis of rotation approach), including relationships for an extension to the mean helical axes methods that facilitate determining joint centres and axes. Through the derivation of a sound mathematical framework to objectively compare the algorithms we demonstrate that the FHA and SARA approach are equivalent for the analysis of two time frames. Moreover, we show that the position of a helical axis derived from the IHA using positional data is affected by a systematic error perpendicular to the true axis direction, whereas the axis direction is identical to those computed with either the FHA or SARA approach (true direction). Finally, with an appropriate choice of weighting factors the mean FHA (MFHA) method is equivalent to the Symmetrical Centre of Rotation Estimation (SCoRE) algorithm for determination of a Centre of Rotation (CoR), and similarly, equivalent to the SARA algorithm for determination of an Axis of Rotation (AoR). The deep understanding of the equivalences between methods presented here enables readers to choose numerically efficient, robust methods for determining AoRs and CoRs with confidence.
Three-dimensional collective cell motions in an acinus-like lumen.
Venue: Journal of biomechanics

Collective cell migration is an essential process in embryo development, wound healing, inflammatory response, and cancer invasion. Although cell motions in two-dimensional (2D) monolayers have been studied previously, three-dimensional (3D) collective cell migration, which constantly occurs during embryogenesis such as the establishment of ducts and acini in vivo, remains elusive. In this paper, we develop a cell-based model incorporating cell mechanics and cell motility to address coherent cell motions in a spherical acinus-like lumen with different cell populations. It is found that the interplays between cell persistence, random fluctuation, and geometrical confinement may engender rich and novel migratory modes. In a 3D spherical lumen, two cells may undergo stripe-like or cross-circular coherent rotations, whereas multiple cells can form dynamic twisting or circulating bands, leaving sparse cells at the center or even a hollow cavity in the cell aggregate. The cell density is found to profoundly influence the collective cell migration modes. Our model can reproduce the fundamental features observed in experiments and highlight the role of mechanics in steering 3D collective cell dynamics during mammary acinar morphogenesis.