nils-herrmann/clean_scirepeval huggingface dataset:
Estimating lumbar passive stiffness behaviour from subject-specific finite element models and in vivo 6DOF kinematics.
Venue: Journal of biomechanics

Passive rotational stiffness of the osseo-ligamentous spine is an important input parameter for estimating in-vivo spinal loading using musculoskeletal models. These data are typically acquired from cadaveric testing. Increasingly, they are also estimated from subject-specific imaging-based finite element (FE) models, which are typically built from CT/MR data obtained in supine position and employ pure rotation kinematics. We explored the sensitivity of FE-based lumbar passive rotational stiffness to two aspects of functional in-vivo kinematics: (a) passive strain changes from supine to upright standing position, and (b) in-vivo coupled translation-rotation kinematics. We developed subject-specific FE models of four subjects' L4L5 segments from supine CT images. Sagittally symmetric flexion was simulated in two ways: (i) pure flexion up to 12° under a 500 N follower load directly from the supine pose. (ii) First, a displacement-based approach was implemented to attain the upright pose, as measured using Dynamic Stereo X-ray (DSX) imaging. We then simulated in-vivo flexion using DSX imaging-derived kinematics. Datasets from weight-bearing motion with three different external weights [(4.5 kg), (9.1 kg), (13.6 kg)] were used. Accounting for supine-upright motion generated compressive pre-loads ≈ 468 N (±188 N) and a "pre-torque" ≈2.5 Nm (±2.2 Nm), corresponding to 25% of the reaction moment at 10° flexion (case (i)). Rotational stiffness estimates from DSX-based coupled translation-rotation kinematics were substantially higher compared to pure flexion. Reaction Moments were almost 90% and 60% higher at 5° and 10° of L4L5 flexion, respectively. Within-subject differences in rotational stiffness based on external weight were small, although between-subject variations were large.
Accuracy of image data stream of a markerless motion capture system in determining the local dynamic stability and joint kinematics of human gait.
Venue: Journal of biomechanics

Assessment of gait parameters is commonly performed through the high-end motion tracking systems, which limits the measurement to sophisticated laboratory settings due to its excessive cost. Recently, Microsoft Kinect (v2) sensor has become popular in clinical gait analysis due to its low-cost. But, determining the accuracy of its RGB-D image data stream in measuring the joint kinematics and local dynamic stability remains an unsolved problem. This study examined the suitability of Kinect(v2) RGB-D image data stream in assessing those gait parameters. Fifteen healthy participants walked on a treadmill during which lower body kinematics were measured by a Kinect(v2) sensor and a optophotogrametric tracking system, simultaneously. Extended Kalman filter was used to extract the lower extremity joint angles from Kinect, while inverse kinematics was used for the gold standard system. For both systems, local dynamic stability was assessed using maximal Lyapunov exponent. Sprague's validation metrics, root mean square error (RMSE) and normalized RMSE were computed to confirm the difference between the joint angles time series of the two systems while relative agreement between them was investigated through Pearson's correlation coefficient (pr). Fisher's Exact Test was performed on maximal Lyapunov exponent to investigate the data independence while reliability was assessed using intraclass correlation coefficients. This study concludes that the RGB-D data stream of Kinect sensor is efficient in estimating joint kinematics, but not suitable for measuring the local dynamic stability.
Effect of side branch flow upon physiological indices in coronary artery disease.
Venue: Journal of biomechanics

Recent efforts have demonstrated the ability of computational models to predict fractional flow reserve from coronary artery imaging without the need for invasive instrumentation. However, these models include only larger coronary arteries as smaller side branches cannot be resolved and are therefore neglected. The goal of this study was to evaluate the impact of neglecting the flow to these side branches when computing angiography-derived fractional flow reserve (vFFR) and indices of volumetric coronary artery blood flow. To compensate for the flow to side branches, a leakage function based upon vessel taper (Murray's Law) was added to a previously developed computational model of coronary blood flow. The augmented model with a leakage function (1Dleaky) and the original model (1D) were then applied to predict FFR as well as inlet and outlet flow in 146 arteries from 80 patients who underwent invasive coronary angiography and FFR measurement. The results show that the leakage function did not significantly change the vFFR but did significantly impact the estimated volumetric flow rate and predicted coronary flow reserve. As both procedures achieved similar predictive accuracy of vFFR despite large differences in coronary blood flow, these results suggest careful consideration of the application of this index for quantitatively assessing flow.
Shear stress rosettes capture the complex flow physics in diseased arteries.
Venue: Journal of biomechanics

Wall shear stress (WSS) is an important parameter in arterial mechanobiology. Various flow metrics, such as time averaged WSS (TAWSS), oscillatory shear index (OSI), and transWSS, have been used to characterize and relate possible WSS variations in arterial diseases like aneurysms and atherosclerosis. We use a graphical representation of WSS using shear rosettes to map temporal changes in the flow dynamics during a cardiac cycle at any spatial location on the vessel surface. The presence of secondary flows and flow reversals can be interpreted directly from the shape of the shear rosette. The mean WSS is given by the rosette centroid, the OSI by the splay around the rosette origin, and the transWSS by its width. We define a new metric, anisotropy ratio (AR), based on the ratio of the length to width of the shear rosette, to capture flow bi-directionality. We characterized the flow physics in controls and patient specific geometries of the ascending aorta (AA) and internal carotid artery (ICA) that have fundamentally different flow dynamics due to differences in the Reynolds and Womersley numbers. The differences in the flow dynamics are well reflected in the shapes of the WSS rosettes and the corresponding flow metrics.
Identifying predictors of upper extremity muscle elasticity with healthy aging.
Venue: Journal of biomechanics

Ultrasound shear wave elastography (SWE) can provide accurate in vivo measurements of the effect of advanced age on muscle elasticity. Our objective was to determine whether passive muscle elasticity was influenced by posture, chronological age, sex, body mass index, and clinical measures of upper extremity function for healthy adults. The dominant arm of 33 male and 33 female participants (ranging from 20 to 89 years old) was examined using a Supersonic Imagine Aixplorer ultrasound SWE system. The mean and standard deviation of shear wave velocity (SWV) was measured from elastography maps for five upper extremity muscles examined at rest: anterior deltoid (AD), biceps brachii (BB), clavicular (CL) and sternocostal (SC) region of the pectoralis major and middle trapezius (MT). Linear mixed models for each muscle were used to assess how SWV was influenced by humeral elevation, chronological age, sex, BMI and three functional measures. All significances are reported at α = 0.05. Humeral elevation influenced shear wave velocity at a statistically significant level for AD, BB, SC and MT (all p < 0.047). Chronological age was a significant predictor of mean SWV for the sternocostal region of the pectoralis major and the middle trapezius (both p < 0.03). These same muscles were also less homogenous (based on their standard deviations) with increased age, particularly for female participants. Performance-based functional assessments of the upper extremity were predictors of mean SWV for the clavicular region of the pectoralis major (all p < 0.04). These results suggest ultrasound SWE has potential utility for assessing age-related changes to muscle elasticity, but these associations were muscle-dependent.
Effect of changes in the lumbar posture in lifting on trunk muscle and spinal loads: A combined in vivo, musculoskeletal, and finite element model study.
Venue: Journal of biomechanics

Irrespective of the lifting technique (squat or stoop), the lumbar spine posture (more kyphotic versus more lordotic) adopted during lifting activities is an important parameter affecting the active-passive spinal load distribution. The advantages in either posture while lifting remains, however, a matter of debate. To comprehensively investigate the role on the trunk biomechanics of changes in the lumbar posture (lordotic, free or kyphotic) during forward trunk flexion, validated musculoskeletal and finite element models, driven by in vivo kinematics data, were used to estimate detailed internal tissue stresses-forces in and load-sharing among various joint active-passive tissues. Findings indicated that the lordotic posture, as compared to the kyphotic one, resulted in marked increases in back global muscle activities (~14-19%), overall segmental compression (~7.5-46.1%) and shear (~5.4-47.5%) forces, and L5-S1 facet joint forces (by up to 80 N). At the L5-S1 level, the lordotic lumbar posture caused considerable decreases in the moment resisted by passive structures (spine and musculature, ~14-27%), negligible reductions in the maximum disc fiber strains (by ~0.4-4.7%) and small increases in intradiscal pressure (~1.8-3.4%). Collectively and with due consideration of the risk of fatigue and viscoelastic creep especially under repetitive lifts, current results support a free posture (in between the extreme kyphotic and lordotic postures) with moderate contributions from both active and passive structures during lifting activities involving trunk forward flexion.
A simple geometrical model accounting for 3D muscle architectural changes across muscle lengths.
Venue: Journal of biomechanics

Muscle architecture parameters change when the muscle changes in length. This has multiple effects on the function of the muscle, e.g. on force production and on contraction velocity. Here we present a versatile geometrical model that predicts changes in muscle architecture as a consequence of length changes of the muscle on the basis of the known architecture at a given muscle length. The model accounts for small changes in aponeuroses' dimensions relative to changes in fascicle length and keeps muscle volume constant. We evaluate the model on the rabbit soleus muscle by comparing model predictions of fascicle lengths and pennation angles with experimental data. For this, we determined the internal architecture of the soleus muscle at different muscle belly lengths (67.8 mm at 35° ankle angle and 59.3 mm at 80° ankle angle). The long and the short soleus muscle exhibited mean fascicle lengths and pennation angles of 20.8 ± 1.3 mm, 4 ± 2° and 13.5 ± 1 mm, 10 ± 4°, respectively. The model predicted reasonable mean fascicle lengths and pennation angles for the long and short soleus that differed only by 1 mm and 1° from the measured data, respectively. Differences between predicted and measured distributions seem to stem from interindividual variability in muscle architecture. Even if the proposed approach has been used for the soleus muscle, which is relatively simple in architecture, it is not restricted to homogeneous unipennate architectures.
A parametric model for studying the aorta hemodynamics by means of the computational fluid dynamics.
Venue: Journal of biomechanics

Perturbed aorta hemodynamics, as for the carotid and the coronary artery, has been identified as potential predicting factor for cardiovascular diseases. In this study, we propose a parametric study based on the computational fluid dynamics with the aim of providing information regarding aortic disease. In particular, the blood flow inside a parametrized aortic arch is computed as a function of morphological changes of baseline aorta geometry. Flow patterns, wall shear stress, time average wall shear stress and oscillatory shear index were calculated during the cardiac cycle. The influence of geometrical changes on the hemodynamics and on these variables was evaluated. The results suggest that the distance between inflow and aortic arch and the angle between aortic arch and descending trunk are the most influencing parameters regarding the WSS-related indices while the effect of the inlet diameter seems limited. In particular, an increase of the aforementioned distance produces a reduction of the spatial distribution of the higher values of the time average wall shear stress and of the oscillatory shear index independently on the other two parameters while an increase of the angle produce an opposite effect. Moreover, as expected, the analysis of the wall shear stress descriptors suggests that the inlet diameter influences only the flow intensity. As conclusion, the proposed parametric study can be used to evaluate the aorta hemodynamics and could be also applied in the future, for analyzing pathological cases and virtual situations, such as pre- and/or post-operative cardiovascular surgical states that present enhanced changes in the aorta morphology yet promoting important variations on the considered indexes.
Microniche geometry modulates the mechanical properties and calcium signaling of chondrocytes.
Venue: Journal of biomechanics

In articular cartilage, the function of chondrocytes is strongly related to their zone-specific microniche geometry defined by pericellular matrix. Microniche geometry is critical for regulating the phenotype and function of the chondrocyte in native cartilage and tissue engineering constructs. However the role of microniche geometry in the mechanical properties and calcium signaling of chondrocytes remains unknown. To recapitulate microniche geometry at single-cell level, we engineered three basic physiological-related polydimethylsiloxane (PDMS) microniches geometries fabricated using soft lithography. We cultured chondrocytes in these microniche geometries and quantified cell mechanical properties using atomic force microscopy (AFM). Fluorescent calcium indicator was used to record and quantify cytosolic Ca2+ oscillation of chondrocytes in different geometries. Our work showed that microniche geometry modulated the mechanical behavior and calcium signaling of chondrocytes. The ellipsoidal microniches significantly enhanced the mechanical properties of chondrocytes compared to spheroidal microniche. Additionally, ellipsoidal microniches can markedly improved the amplitude but weakened the frequency of cytosolic Ca2+ oscillation in chondrocytes than spheroidal microniche. Our work might reveal a novel understanding of chondrocyte mechanotransduction and therefore be useful for designing cell-instructive scaffolds for functional cartilage tissue engineering.
Finite element analysis of MitraClip procedure on a patient-specific model with functional mitral regurgitation.
Venue: Journal of biomechanics

Mitral valve (MV) repair with the MitraClip device has been shown to reduce mitral regurgitation severity and improve clinical outcomes in symptomatic patients at high surgical risk. MitraClip was recently approved in the US for the treatment of functional mitral regurgitation (FMR), which significantly expands the number of patients that can be treated with this device. This study aims to quantify the morphologic changes and evaluate the biomechanical interaction between the MitraClip device and the mitral apparatus of a real patient case with FMR using computational modeling. MitraClip procedures using a central and a lateral clip were simulated in a validated MV-left ventricle finite element (FE) model with severe MR. The patient-specific model integrated detailed geometries of the left ventricle, mitral leaflets and chordae, incorporated age- and gender-matched nonlinear hyperelastic human material properties, and accounted for chordae tethering forces. Central and lateral positioning gave similar biomechanical outcomes resulting in an improved but incomplete MV coaptation. Antero-posterior distance, annulus area, valve opening orifice area, and regurgitant orifice area decreased by up to 26%, 19%, 48% and 63% when compared to the pre-clip model, respectively. Anterior and posterior leaflet peak stresses increased by up to 64% and 62% after clip placement, respectively, and were located at the region of clip grasp. Similarly, anterior and posterior leaflet peak strains increased by up to 20% and 10%, respectively. FE modeling, as used here, can be a powerful tool to examine the complex MitraClip-host biomechanical interaction.
Comparison of the compressive mechanical properties of auricular and costal cartilage from patients with microtia.
Venue: Journal of biomechanics

Children born with a small or absent ears undergo surgical reconstruction to restore their auricle. Currently, rib (costal) cartilage is used to carve the auricle. However as alternative, tissue engineered and synthetic materials are being developed to restore the auricle shape to overcome donor site morbidity and limited availability of rib cartilage. However, to date there is limited knowledge regarding the mechanical properties of the auricular and costal cartilage to optimise the required compressive properties of the graft. The remnant auricular and costal cartilage from 20 patients undergoing stage-1 microtia surgery was harvested. On the day of surgery, the cartilage was evaluated in compression, with each sample loaded to 300 g at 1 mm/s.  The costal cartilage was observed to have a significantly higher Young's Elastic Modulus than auricular cartilage (average costal cartilage 11.43 MPa vs average auricular cartilage 2 MPa, p < 0.0001). The auricular cartilage showed a significantly higher relaxation rate than costal cartilage (average costal cartilage 0.72 MPa10-4 vs average auricular cartilage 1.93 MPa10-4, p < 0.05). The final absolute relaxation was significantly lower for elastic cartilage than costal cartilage (average costal cartilage 3.35 MPa vs average auricular cartilage 0.2 MPa, p < 0.0001). Alloplastic cartilage replacements used as alternatives for reconstruction were also evaluated. Silicone, Gore-Tex and Medpor were observed to have significantly higher Young's Elastic Modulus than costal and auricular cartilage. Costal cartilage has a higher Young's Elastic Modulus in compression compared to auricular cartilage. Current synthetic materials used to replace synthetic cartilage do not mimic costal cartilage, which should be addressed in the future.
Comparison of lower limb joint moment and power during turning gait between young and old adults using hierarchical Bayesian inference.
Venue: Journal of biomechanics

Age-related differences in lower limb joint moment (JM), and joint power (JP) during turning remain unclear. The present study investigated age-related differences in lower limb JM and JP during turning between young adults (YAs) and old adults (OAs). We introduced the hierarchical Bayesian inference for comparing and identifying differences in JM, angular velocity(ω), and JP at each stance phase in the two age groups. This study included 16 healthy YAs and 16 healthy OAs (8 men and 8 women in each group). Participants performed 90° step turns to the right at a self-selected natural speed. On comparing the age groups, during 90° step turning, the OA group exhibited larger extention hip JM and JP to control (brake) the upper body in the sagittal plane, exhibited larger abductor moment in each lower limb joint for preventing the body from leaning in the frontal plane during the mid-stance phase, and exhibited larger hip JP and ω and smaller ankle JM in the transverse plane to rotate the body during the mid-stance phase. Our findings suggested that the overall reliance on the hip joint to control body motion in each anatomical plane during step turning is higher in the OA group than in the YA group. In addition, the hierarchical Bayesian inference is useful for comparing the time courses of JM, ω, and JP.
Validation of a wearable system for 3D ambulatory L5/S1 moment assessment during manual lifting using instrumented shoes and an inertial sensor suit.
Venue: Journal of biomechanics

This study aimed to evaluate the accuracy of 3D L5/S1 moment estimates from an ambulatory measurement system consisting of a wearable inertial motion capture system (IMC) and instrumented force shoes (FSs), during manual lifting. Reference L5/S1 moments were calculated using an inverse dynamics bottom-up laboratory model (buLABmodel), based on data from a measurement system comprising optical motion capture (OMC) and force plates (FPs). System performance of (1) a bottom-up ambulatory model (buAMBmodel) using lower-body kinematic IMC and FS data, and (2) a top-down ambulatory model (tdAMBmodel) using upper-body kinematic IMC data and hand forces (HFs) were compared. HFs were estimated using full-body kinematic IMC data and FS forces. Eight males and eight females lifted a 10-kg box from different initial vertical/horizontal positions using either a free or an asymmetric lifting style. As a measure of system performance, root-mean-square (RMS) errors were calculated between the reference (buLABmodel) and ambulatory (tdAMBmodel &buAMBmodel) moments. The results showed two times smaller errors for the tdAMBmodel (averaged RMS errors < 20 Nm or 10% of peak extension moment) than for the buAMBmodel (average RMS errors < 40 Nm or 20% of peak extension moment). In conclusion, for ambulatory L5/S1 moment assessment with an IMC + FS system, using a top-down inverse dynamics approach with estimated hand forces is to be preferred over a bottom-up approach.
Exploring the regional disc bulge response of the cervical porcine intervertebral disc under varying loads and posture.
Venue: Journal of biomechanics

Nerve compression due to intervertebral disc (IVD) bulging is a known mechanism for low back pain and typically occurs in the posterior region of the disc. Most in vitro studies are limited in the ability to quantify the magnitude of bulging on the posterior aspect of the disc due to the boney structures that occlude a direct line-of-sight in the intact functional spinal units (FSUs). This study examined anterior and posterior annulus fibrosus (AF) bulges in reduced (posterior elements removed) cervical porcine specimens across four loading conditions and two postures. Surface scans from the anterior and posterior aspect of the IVD were recorded in both neutral and flexed postures using a 3D laser scanner to characterize changes in AF bulge. A significant negative correlation was observed for peak AF bulge on the anterior and posterior side of the disc in a flexed posture (Pearson's r = -0.448; p = 0.002; r2 = 0.2003). The results from this investigation support that there may be a connection between the magnitude of AF bulge on the posterior side and estimations computed using the anterior side.