nils-herrmann/clean_scirepeval huggingface dataset:
A finite element analysis for evaluating mandibular advancement devices.
Venue: Journal of biomechanics

Obstructive sleep apnoea (OSA) is a disorder characterised by complete or partial occlusion of the upper airway during sleep. Muscles relax during sleeping and collapse into the airway, closing the throat and prohibiting air flowing into the lungs. Different solutions have been adopted to manage the pathology to improve the life quality of affected patients. Mandibular advancement devices (MADs) are proven to be a compliant and successful therapy in the forward repositioning of the mandible to increase the upper airway volume. However, this method has some long-term adverse events that may affect the teeth and periodontal ligaments. This paper presents a finite element model to evaluate the MADs effects (displacement and stress) on teeth and periodontal ligaments, by varying the design, the point of application of the force and the material. The modelled bodies have been reconstructed through a Reverse Engineering approach and computer-aided design tools starting from tomographic images of anatomic bodies and from laser scans of a physical MAD. The results suggest that a central connection mechanism could affect mostly the anterior teeth. In contrast, a lateral connection mechanism provides a more uniform distribution of the load on teeth.
Anatomy and biomechanics of the Lisfranc ligamentous complex: A systematic literature review.
Venue: Journal of biomechanics

Lisfranc injuries are challenging to treat and can have a detrimental effect on active individuals. Over the past decade researchers have investigated methods for the reconstruction of the Lisfranc ligamentous complex (LLC) to preserve its functional stability and mobility. To aid in this innovation, this study presents the current understanding of the anatomical and biomechanical characteristics of the LLC through a systematic review. Three medical databases (PubMed, Scopus, and Embase) were searched from inception through July 2019. Original studies investigating the anatomy and/or biomechanical properties of the LLC were considered for inclusion. Data recorded from each study included: number of cadavers, number of feet, gender, laterality, age, type of specimen, measurement methods, reported ligamentous bundles, ligament origins and insertions, geometric characteristics, and biomechanical properties of the LLC. The Quality Appraisal for Cadaveric Studies (QUACS) scale was used to assess the methodologic quality of included articles. Eight cadaveric studies investigating the LLC were included out of 1204 screened articles. Most articles described the LLC as three distinct structures: the dorsal- (DLL), interosseous- (ILL), and plantar- (PLL) Lisfranc Ligaments. The ILL had the largest thickness and insertional area of osseous attachment. Biomechanically, the ILL also had the highest stiffness and resistance to load prior to failure when loaded parallel to its fiber orientation. Current knowledge of the anatomical and biomechanical properties of the LLC are presented and highlight its significant role of stabilizing the tarsometatarsal articulation. Appreciating the biomechanical characteristics of the ILL may improve clinical insight in managing LLC injuries.
Hopping frequency influences elastic energy reuse with joint series elastic actuators.
Venue: Journal of biomechanics

Robotic limb design struggles to combine energy efficiency with human-like levels of movement versatility. High efficiency and a range of angles and torques are characteristics of human hopping at different frequencies. Humans use muscles in combination with tendons to achieve the required joint actuation. Therefore, we consider whether appropriately tuned series elastic actuators (SEAs) placed at the leg joints can be used to reduce the functional gap between robots and humans. Human hip, knee, and ankle biomechanics were recorded over a range of hopping frequencies to extract joint angles and torques, which were used as an input to a mechanical simulation SEA model. This model was used to optimize the SEA stiffness of each joint to either minimize peak power or energy requirements. This work investigated the relationship between hopping frequency and SEA stiffness, the utility of using SEAs at each joint, and the reasons behind humans' preferred hopping frequency. Although the constant stiffness values across different hopping frequencies are suitable for the knee and the ankle, a variable serial elastic actuator stiffness could still further reduce energy requirements. Optimal SEA stiffness was found to reduce peak power requirements by up to 73% at the ankle and up to 66% at the knee, with greatest benefits found around the preferred frequency. However, no SEA benefits were found for the hip and above the preferred hopping frequency for the knee. These insights could be used to aid in the design of robotic and assistive devices to achieve versatile and energy efficient human-like movements.
Modulation of sagittal-plane center of pressure and force vector direction in human standing on sloped surfaces.
Venue: Journal of biomechanics

The multi-joint coordination responsible for maintaining upright posture in the standing human manifests in the pattern of variation of the support-surface force (F). Assessment of both the translational and rotational kinematics in the sagittal-plane requires understanding the critical relationship between the direction and location of F. Prior work demonstrated that band-pass filtered F direction and center-of-pressure (CoP) covary in time such that the F vector lines-of-action pass near a fixed point called an intersection point (IP). The height of that IP (IPz) varies systematically with the frequency of the pass band. From F measurements in able-bodied humans (n = 17) standing on various pitched surfaces, the present study also found the emergent property of an IP, with IPz located above the center of mass (CoM) at frequencies <1.75 Hz and below the CoM for higher frequencies. This property aids in maintaining upright posture for various perturbation modes within a single control structure. From purely mechanical effects, standing on a pitched surface should not change IPz, however these measurements of F show that IPz is generally closer to CoM height. This characterization of quiet standing provides simple means of assessing the complex multi-joint coordination of standing and relates directly to the physical demands of controlling the translational and rotational aspects of body posture.
Assessment of femoral neck strength and bone mineral density changes following exercise using 3D-DXA images.
Venue: Journal of biomechanics

Physical exercise induces spatially heterogeneous bone changes in the proximal femur. Recent advances have enabled 3D dual-energy X-ray Absorptiometry (DXA)-based finite element (FE) models to estimate bone strength. However, its ability to detect exercise-induced BMD and strength changes is unclear. The aim of this study was to quantify the repeatability of vBMD and femoral neck strength obtained from 3D-DXA images and determine the changes due an exercise intervention. The DXA scans included pairs of same-day repeated scans from ten healthy females and pre- and post-exercise intervention scans of 26 males. FE models with element-by-element correspondence were generated by morphing a template mesh to each bone. BMD and femoral strength under single-leg-stance and sideways fall loading configurations were obtained for both groups and compared. In the repeated images, the total hip vBMD difference was 0.5 ± 2.5%. Element-by-element BMD differences reached 30 ± 50%. The strength difference in single-leg stance was 2.8 ± 13% and in sideways fall was 4.5% ± 19%. In the exercise group, strength changes were 6 ± 19% under single-leg stance and 1 ± 18% under sideways fall. vBMD parameters were weakly correlated to strength (R2 < 0.31). The exercise group had a mean bone accrual exceeding repeatability values in the femoral head and cortical regions. The case with the highest vBMD change (6.4%) caused 18% and -7% strength changes under single-leg stance and sideways fall. 3D-DXA technology can assess the effect of exercise interventions in large cohorts but its validity in individual cases should be interpreted with caution.
Side-specific valvular endothelial-interstitial cell mechano-communication via cadherin-11.
Venue: Journal of biomechanics

Calcific aortic valve disease (CAVD) is a condition causing stiffening of the aortic valve, impeding cardiac function and resulting in significant morbidity worldwide. CAVD is thought to be driven by the persistent activation of the predominant cell type in the valve, aortic valve interstitial cells (AVICs), into myofibroblasts, resulting in subsequent calcification and stenosis of the valve. Although much of the research into CAVD focuses on AVICs, the aortic valve endothelial cells (AVECs) have been shown to regulate AVICs and maintain tissue homeostasis. Exposed to distinct flow patterns during the cardiac cycle, the AVECs lining either side of the valve demonstrate crucial differences which could contribute to the preferential formation of calcific nodules on the aorta-facing (fibrosa) side of the valve. Cadherin-11 (CDH11) is a cell-cell adhesion protein which has been previously associated with AVIC myofibroblast activation, nodule formation, and CAVD in mice. In this study, we investigated the role of CDH11 in AVECs and examined side-specific differences. The aorta-facing or fibrosa endothelial cells (fibAVECs) express higher levels of CDH11 than the ventricle-facing or ventricularis endothelial cells (venAVECs). This increase in expression corresponds with increased contraction of a free-floating collagen gel compared to venAVECs. Additionally, co-culture of fibAVECs with AVICs demonstrated decreased contraction compared to an AVIC + AVIC control, but increased contraction compared to the venAVECs co-culture. This aligns with the known preferential formation of calcific nodules on the fibrosa. These results together indicate a potential role for CDH11 expression by AVECs in regulating AVIC contraction and subsequent calcification.
4D flow evaluation of blood non-Newtonian behavior in left ventricle flow analysis.
Venue: Journal of biomechanics

Blood is generally modeled as a Newtonian fluid, assuming a standard and constant viscosity; however, this assumption may not hold for the highly pulsatile and recirculating intracavitary flow in the left ventricle (LV), hampering the quantification of fluid dynamic indices of potential clinical relevance. Herein, we investigated the effect of three viscosity models on the patient-specific quantification of LV blood energetics, namely on viscous energy loss (EL), from 4D Flow magnetic resonance imaging: I) Newtonian with standard viscosity (3.7 cP), II) Newtonian with subject-specific hematocrit-dependent viscosity, III) non-Newtonian accounting for the effect of hematocrit and shear rate. Analyses were performed on 5 controls and 5 patients with cardiac light-chain amyloidosis. In Model II, viscosity ranged between 3.0 (-19%) and 4.3 cP (+16%), mildly deviating from the standard value. In the non-Newtonian model, this effect was emphasized: viscosity ranged from 3.2 to 6.0 cP, deviating maximally from the standard value in low shear rate (i.e., <100 s-1) regions. This effect reflected on EL quantifications: in particular, as compared to Model I, Model III yielded markedly higher EL values (up to +40%) or markedly lower (down to -21%) for subjects with hematocrit higher than 39.5% and lower than 30%, respectively. Accounting for non-Newtonian blood behavior on a patient-specific basis may enhance the accuracy of intracardiac energetics assessment by 4D Flow, which may be explored as non-invasive index to discriminate between healthy and pathologic LV.
A novel coupled musculoskeletal finite element model of the spine - Critical evaluation of trunk models in some tasks.
Venue: Journal of biomechanics

Spine musculoskeletal (MS) models make simplifying assumptions on the intervertebral joint degrees-of-freedom (rotational and/or translational), representation (spherical or beam-like joints), and properties (linear or nonlinear). They also generally neglect the realistic structure of the joints with disc nuclei/annuli, facets, and ligaments. We aim to develop a novel MS model where trunk muscles are incorporated into a detailed finite element (FE) model of the ligamentous T12-S1 spine thus constructing a gold standard coupled MS-FE model. Model predictions are compared under some tasks with those of our earlier spherical joints, beam joints, and hybrid (uncoupled) MS-FE models. The coupled model predicted L4-L5 intradiscal pressures (R2 ≅ 0.97, RMSE ≅ 0.27 MPa) and L1-S1 centers of rotation (CoRs) in agreement to in vivo data. Differences in model predictions grew at larger trunk flexion angles; at the peak (80°) flexion the coupled model predicted, compared to the hybrid model, much smaller global/local muscle forces (~38%), segmental (~44%) and disc (~22%) compression forces but larger segmental (~9%) and disc (~17%) shear loads, ligament forces at the lower lumbar levels (by up to 57%) and facet forces at all levels. The spherical/beam joints models predicted much greater muscle forces and segmental loads under larger flexion angles. Unlike the spherical joints model with fixed CoRs, the beam joints model predicted CoRs closer (RMSE = 2.3 mm in flexion tasks) to those of the coupled model. The coupled model offers a great potential for future studies towards improvement of surgical techniques, management of musculoskeletal injuries and subject-specific simulations.
Lower limb biomechanics during gait in children with Achondroplasia.
Venue: Journal of biomechanics

Achondroplasia is an autosomal dominant disorder that has been reported to be the most common condition associated with severe disproportionate short stature. A limited number of studies exist in the literature reporting kinematic abnormalities in children with Achondroplasia. Reported deviations include increased pelvic tilt, flexed knees during stance and excessive hip abduction. However, studies lack full 3-dimensional data and results are variable. The aim of this study was to examine differences in 3-dimensional kinematics during gait between children with Achondroplasia and age-matched controls A retrospective analysis was conducted of our gait laboratory database over a 20-year period. Eight children with Achondroplasia and 32 children with typical development were included. A full 3-dimensional kinematic analysis was performed using a CODA mpx30/cx1 active marker system. Increased flexion at all joints in the sagittal plane, increased hip abduction and foot supination in the coronal plane and increased external hip rotation were present for children with Achondroplasia. In contrast to other studies, increased knee varus was not an issue for the children in this study. This study provides a full description of the kinematic patterns associated with children with Achondroplasia and, as such, a comprehensive baseline for planning rehabilitation strategies in these children.
Stair ascent comparison of lower limb kinematics with differing time normalization techniques.
Venue: Journal of biomechanics

Understanding gait differences in context of group differences is dependent on statistical testing methods and time normalization techniques (TN). The method induces a relationship of both with one another. As to our knowledge, there has been no investigation into their relation so far. To show empirically what effects may be of importance, we use SPM with linear time interpolation (LI) and Dynamic Time Warping (DTW) separately for data of a study on stair ascent kinematics between two groups. There is a slight difference in statistical significance for the comparison of LI and DTW. LI-uniquely significant time highlight differences due to in-group time-variations, whereas DTW-uniqueness is tied to qualitative differences of homogeneous events. The comparison of stair ascent kinematics with DTW shows more pronounced evidence for backlift-like strategies for the older group, although trunk angles are kept more extended as to ensure stabilty. Thus, the difference in SPM from TN is slight but important, if there is need to mirror said effects methodically.
Microheart: A microfluidic pump for functional vascular culture in microphysiological systems.
Venue: Journal of biomechanics

Advances in microphysiological systems have prompted the need for long-term cell culture under physiological flow conditions. Conventional laboratory pumps typically lack the ability to deliver cell culture media at the low flow rates required to meet the physiological ranges of fluid flow, and are often pulsatile or require flow reversal. Here, a microfluidic-based pump is presented, which allows for the controlled delivery of media for vascular microphysiological applications. The performance of the pump was characterized in a range of microfluidic systems, including straight channels of varying dimensions and self-assembled microvascular networks. A theoretical framework was developed based on lumped element analysis to predict the performance of the pump for different fluidic configurations and a finite element model of the included check-valves. The use of the pump for microvascular physiological studies demonstrated the utility of this system to recapitulate vascular fluid transport phenomena in microphysiological systems, which may find applications in disease models and drug screening.
A converging artery-sized model for shear adhesion mapping of particles.
Venue: Journal of biomechanics

Drug carriers for targeting cardiovascular diseases have been gaining a respectable attention, however, designing such carriers is challenging due to the biophysical complexity of the vascular system. Wall shear stress (WSS), exerted by blood flow on the endothelium surface, is a crucial factor in the circulatory system. WSS affects the adhesion and preferential accumulation of drug carriers. Here, we propose, an innovative approach to investigate particle adhesion in a converging artery-sized model, lined with human endothelial cells. Unlike widely used microfluidic and in vivo setups, our model enables to investigate particle accumulation in a continuous WSS range, performed in a single experiment, and at the right scale relevant for human arteries. First, we characterized the flow and the WSS map along the designed model, which can span along the entire arterial WSS range. We then used the model to examine the effect of particle size and the suspension buffer on particle adhesion distribution. The results demonstrated the role of particle size, where the same particles with a diameter of 2 µm exhibit shear-decreased adhesion while 500 nm particles exhibit shear-enhanced adhesion. Furthermore, under the same WSS, particles show a similar behavior when suspended in a Dextran buffer, having a viscosity analogous to blood, compared to a phosphate buffer solution without Dextran. Moreover, experiments with RBCs in the phosphate buffer, at a 40% physiological hematocrit, decreased particle adhesion and affected their deposition pattern. Altogether, our study suggests an original platform for investigating and optimizing intravascular drug carriers and their targeting properties.
Assessment of blood distribution in response to post-surgical steal syndrome: A novel technique based on Thermo-Anatomical Segmentation.
Venue: Journal of biomechanics

The distal ischemic steal syndrome (ISS) is a complication following the construction of an arteriovenous (A-V) access for hemodialysis. The ability to non-invasively monitor changes in skin microcirculation improves both the diagnosis and treatment of vascular diseases. In this study, we propose a novel technique for evaluating the palms' blood distribution following arteriovenous access, based on thermal imaging. Furthermore, we utilize the thermal images to identify typical recovery patterns of patients that underwent this surgery and show that thermal images taken post-surgery reflect the patient's follow-up status. Thermal photographs were taken by a portable thermal camera from both hands before and after the A-V access surgery, and one month following the surgery, from ten dialysis patients. A novel term "Thermo-Anatomical Segmentation", which enables a functional assessment of palm blood distribution was defined. Based on this segmentation it was shown that the greatest change after surgery was in the most distal region, the fingertips (p < 0.05). In addition, the changes in palm blood distribution in both hands were synchronized, which indicates a bilateral effect. An unsupervised machine learning model revealed two variables that determine the recovery pattern following the surgery: the palms' temperature difference pre- and post-surgery and the post-surgery difference between the treated and untreated hand. Our proposed framework provides a new technique for quantitative assessment of the palm's blood distribution. This technique may improve the clinical treatment of patients with vascular disease, particularly the patient-specific follow-up, in clinics as well as in homecare.
Natural oscillatory modes of 3D deformation of the human brain in vivo.
Venue: Journal of biomechanics

Natural modes and frequencies of three-dimensional (3D) deformation of the human brain were identified from in vivo tagged magnetic resonance images (MRI) acquired dynamically during transient mild acceleration of the head. Twenty 3D strain fields, estimated from tagged MRI image volumes in 19 adult subjects, were analyzed using dynamic mode decomposition (DMD). These strain fields represented dynamic, 3D brain deformations during constrained head accelerations, either involving rotation about the vertical axis of the neck or neck extension. DMD results reveal fundamental oscillatory modes of deformation at damped frequencies near 7 Hz (in neck rotation) and 11 Hz (in neck extension). Modes at these frequencies were found consistently among all subjects. These characteristic features of 3D human brain deformation are important for understanding the response of the brain in head impacts and provide valuable quantitative criteria for the evaluation and use of computer models of brain mechanics.
Differences in shear elastic modulus of the latissimus dorsi muscle during stretching among varied trunk positions.
Venue: Journal of biomechanics

The latissimus dorsi (LD) can be divided into the upper, middle, and lower parts, but the effective stretching positions for each part are unknown. In this study, we aimed to investigate effective trunk positions for stretching of the LD. A total of 14 healthy males participated in this study. The following seven trunk positions were chosen as the LD stretching positions; upright of the trunk (Baseline), flexion of the trunk (Flex), contralateral bending of the trunk (LB), contralateral rotation of the trunk (Rot), flexion and contralateral bending of the trunk (Flex + LB), flexion and contralateral rotation of the trunk (Flex + Rot), and contralateral bending and contralateral rotation of the trunk (LB + Rot). Maximal elevation of the upper limb was passively added to all positions. The shear elastic modulus, used as the index of muscle elongation, was measured at the four parts (upper, middle, lower, distal parts) of the LD. The shear elastic moduli showed obviously high values in Rot and LB + Rot at the upper, middle, and distal parts, and also in LB, Rot, and LB + Rot at the lower part. These findings suggest that contralateral trunk rotation, or a combination of contralateral trunk bending and rotation are effective trunk positions for stretching all parts of the LD. Contralateral trunk bending was also effective for stretching the lower part of the LD.