nils-herrmann/clean_scirepeval huggingface dataset:
Metabolic responses induced by compression of chondrocytes in variable-stiffness microenvironments.
Venue: Journal of biomechanics

Cells sense and respond to mechanical loads in a process called mechanotransduction. These processes are disrupted in the chondrocytes of cartilage during joint disease. A key driver of cellular mechanotransduction is the stiffness of the surrounding matrix. Many cells are surrounded by extracellular matrix that allows for tissue mechanical function. Although prior studies demonstrate that extracellular stiffness is important in cell differentiation, morphology and phenotype, it remains largely unknown how a cell's biological response to cyclical loading varies with changes in surrounding substrate stiffness. Understanding these processes is important for understanding cells that are cyclically loaded during daily in vivo activities (e.g. chondrocytes and walking). This study uses high-performance liquid chromatography - mass spectrometry to identify metabolomic changes in primary chondrocytes under cyclical compression for 0-30minutes in low- and high-stiffness environments. Metabolomic analysis reveals metabolites and pathways that are sensitive to substrate stiffness, duration of cyclical compression, and a combination of both suggesting changes in extracellular stiffness in vivo alter mechanosensitive signaling. Our results further suggest that cyclical loading minimizes matrix deterioration and increases matrix production in chondrocytes. This study shows the importance of modeling in vivo stiffness with in vitro models to understand cellular mechanotransduction.
In vitro culture increases mechanical stability of human tissue engineered cartilage constructs by prevention of microscale scaffold buckling.
Venue: Journal of biomechanics

Many studies have measured the global compressive properties of tissue engineered (TE) cartilage grown on porous scaffolds. Such scaffolds are known to exhibit strain softening due to local buckling under loading. As matrix is deposited onto these scaffolds, the global compressive properties increase. However the relationship between the amount and distribution of matrix in the scaffold and local buckling is unknown. To address this knowledge gap, we studied how local strain and construct buckling in human TE constructs changes over culture times and GAG content. Confocal elastography techniques and digital image correlation (DIC) were used to measure and record buckling modes and local strains. Receiver operating characteristic (ROC) curves were used to quantify construct buckling. The results from the ROC analysis were placed into Kaplan-Meier survival function curves to establish the probability that any point in a construct buckled. These analysis techniques revealed the presence of buckling at early time points, but bending at later time points. An inverse correlation was observed between the probability of buckling and the total GAG content of each construct. This data suggests that increased GAG content prevents the onset of construct buckling and improves the microscale compressive tissue properties. This increase in GAG deposition leads to enhanced global compressive properties by prevention of microscale buckling.
Are there characteristic motion patterns in the lumbar spine during flexion?
Venue: Journal of biomechanics

Flexion is the main motion of the lumbar spine. While in vitro tests with pure moments suggest larger intra-segmental rotations for the more caudal segments, in vivo results show diverging motion distributions. The present study analysed the motion distribution in vivo of 320 asymptomatic subjects. The change of the back curvature between standing and upper body flexion was determined using a non-invasive measurement device. Linear, bilinear, trilinear, quadratic, and cubic regression models were fitted to the segmental motion distribution over the lengths of the lordosis to categorise characteristic motion patterns. Simplicity and approximation quality were used to assign the motion distributions to the regression models. Seventy-seven percent of the motion distributions could be explained by a bilinear model. A further 12% and 11% could be represented by a trilinear and linear model, respectively. Less than 1% of the distributions could not satisfactorily be represented by the models. All of the bilinear models displayed maximum flexion in approximately the middle of lordosis. All linear models showed a decreasing rotation from caudal to cranial. Most of the trilinear models displayed a distribution similar to the bilinear. Age, sex, body height, and weight did not significantly affect these distributions. This in vivo study identified characteristic motion patterns in the lumbar spine during flexion. The quantitative results provide a clear description of the healthy condition and may serve to identify spinal motion abnormalities.
Effect of noise and filtering on largest Lyapunov exponent of time series associated with human walking.
Venue: Journal of biomechanics

This study aimed to determine the effect of added noise, filtering and time series length on the largest Lyapunov exponent (LyE) value calculated for time series obtained from a passive dynamic walker. The simplest passive dynamic walker model comprising of two massless legs connected by a frictionless hinge joint at the hip was adopted to generate walking time series. The generated time series was used to construct a state space with the embedding dimension of 3 and time delay of 100 samples. The LyE was calculated as the exponential rate of divergence of neighboring trajectories of the state space using Rosenstein's algorithm. To determine the effect of noise on LyE values, seven levels of Gaussian white noise (SNR=55-25dB with 5dB steps) were added to the time series. In addition, the filtering was performed using a range of cutoff frequencies from 3Hz to 19Hz with 2Hz steps. The LyE was calculated for both noise-free and noisy time series with different lengths of 6, 50, 100 and 150 strides. Results demonstrated a high percent error in the presence of noise for LyE. Therefore, these observations suggest that Rosenstein's algorithm might not perform well in the presence of added experimental noise. Furthermore, findings indicated that at least 50 walking strides are required to calculate LyE to account for the effect of noise. Finally, observations support that a conservative filtering of the time series with a high cutoff frequency might be more appropriate prior to calculating LyE.
The equivalence of multi-axis spine systems: Recommended stiffness limits using a standardized testing protocol.
Venue: Journal of biomechanics

The complexity of multi-axis spine testing often makes it challenging to compare results from different studies. The aim of this work was to develop and implement a standardized testing protocol across three six-axis spine systems, compare them, and provide stiffness and phase angle limits against which other test systems can be compared. Standardized synthetic lumbar specimens (n=5), comprising three springs embedded in polymer at each end, were tested on each system using pure moments in flexion-extension, lateral bending, and axial rotation. Tests were performed using sine and triangle waves with an amplitude of 8Nm, a frequency of 0.1Hz, and with axial preloads of 0 and 500N. The stiffness, phase angle, and R2 value of the moment against rotation in the principal axis were calculated at the center of each specimen. The tracking error was adopted asa measure of each test system to minimize non-principal loads, defined as the root mean squared difference between actual and target loads. All three test systems demonstrated similar stiffnesses, with small (<14%) but significant differences in 4 of 12 tests. More variability was observed in the phase angle between the principal axis moment and rotation, with significant differences in 10 of 12 tests. Stiffness and phase angle limits were calculated based on the 95% confidence intervals from all three systems. These recommendations can be used with the standard specimen and testing protocol by other research institutions to ensure equivalence of different spine systems, increasing the ability to compare in vitro spine studies.
Spinal loads and trunk muscles forces during level walking - A combined in vivo and in silico study on six subjects.
Venue: Journal of biomechanics

During level walking, lumbar spine is subjected to cyclic movements and intricate loading of the spinal discs and trunk musculature. This study aimed to estimate the spinal loads (T12-S1) and trunk muscles forces during a complete gait cycle. Six men, 24-33years walk barefoot at self-selected speed (4-5km/h). 3D kinematics and ground reaction forces were recorded using a motion capturing system and two force plates, implemented in an inverse dynamic musculoskeletal model to predict the spinal loads and trunk muscles forces. Additionally, the sensitivity of the intra-abdominal pressure and lumbar segment rotational stiffness was investigated. Peak spinal loads and trunk muscle forces were between the gait instances of heel strike and toe off. In L4-L5 segment, sensitivity analysis showed that average peak compressive, antero-posterior and medio-lateral shear forces were 130-179%, 2-15% and 1-6%, with max standard deviation (±STD) of 40%, 6% and 3% of the body weight. Average peak global muscles forces were 24-55% (longissimus thoracis), 11-23% (iliocostalis thoracis), 12-16% (external oblique), 17-25% (internal oblique) and 0-8% (rectus abdominus) of body weight whereas, the average peak local muscles forces were 11-19% (longissimus lumborum), 14-31% (iliocostalis lumborum) and 12-17% (multifidus). Maximum±STD of the global and local muscles forces were 13% and 8% of the body weight. Large inter-individual differences were found in peak compressive and trunk muscles forces whereas the sensitivity analysis also showed a substantial variation.
Model-based inverse estimation for active contraction stresses of tongue muscles using 3D surface shape in speech production.
Venue: Journal of biomechanics

This paper presents a novel inverse estimation approach for the active contraction stresses of tongue muscles during speech. The proposed method is based on variational data assimilation using a mechanical tongue model and 3D tongue surface shapes for speech production. The mechanical tongue model considers nonlinear hyperelasticity, finite deformation, actual geometry from computed tomography (CT) images, and anisotropic active contraction by muscle fibers, the orientations of which are ideally determined using anatomical drawings. The tongue deformation is obtained by solving a stationary force-equilibrium equation using a finite element method. An inverse problem is established to find the combination of muscle contraction stresses that minimizes the Euclidean distance of the tongue surfaces between the mechanical analysis and CT results of speech production, where a signed-distance function represents the tongue surface. Our approach is validated through an ideal numerical example and extended to the real-world case of two Japanese vowels, /ʉ/ and /ɯ/. The results capture the target shape completely and provide an excellent estimation of the active contraction stresses in the ideal case, and exhibit similar tendencies as in previous observations and simulations for the actual vowel cases. The present approach can reveal the relative relationship among the muscle contraction stresses in similar utterances with different tongue shapes, and enables the investigation of the coordination of tongue muscles during speech using only the deformed tongue shape obtained from medical images. This will enhance our understanding of speech motor control.
A sensitivity analysis method for the body segment inertial parameters based on ground reaction and joint moment regressor matrices.
Venue: Journal of biomechanics

This paper presents a method allowing a simple and efficient sensitivity analysis of dynamics parameters of complex whole-body human model. The proposed method is based on the ground reaction and joint moment regressor matrices, developed initially in robotics system identification theory, and involved in the equations of motion of the human body. The regressor matrices are linear relatively to the segment inertial parameters allowing us to use simple sensitivity analysis methods. The sensitivity analysis method was applied over gait dynamics and kinematics data of nine subjects and with a 15 segments 3D model of the locomotor apparatus. According to the proposed sensitivity indices, 76 segments inertial parameters out the 150 of the mechanical model were considered as not influent for gait. The main findings were that the segment masses were influent and that, at the exception of the trunk, moment of inertia were not influent for the computation of the ground reaction forces and moments and the joint moments. The same method also shows numerically that at least 90% of the lower-limb joint moments during the stance phase can be estimated only from a force-plate and kinematics data without knowing any of the segment inertial parameters.
Muscle synergies are similar when typically developing children walk on a treadmill at different speeds and slopes.
Venue: Journal of biomechanics

The aim of this study was to determine whether changes in synergies relate to changes in gait while walking on a treadmill at multiple speeds and slopes. The hypothesis was that significant changes in movement pattern would not be accompanied by significant changes in synergies, suggesting that synergies are not dependent on the mechanical constraints but are instead neurological in origin.



Sixteen typically developing children walked on a treadmill for nine combinations (stages) of different speeds and slopes while simultaneously collecting kinematics, kinetics, and surface electromyography (EMG) data. The kinematics for each stride were summarized using a modified version of the Gait Deviation Index that only includes the sagittal plane. The kinetics for each stride were summarized using a modified version of the Gait Deviation Index - Kinetic which includes sagittal plane moments and powers. Within each synergy group, the correlations of the synergies were calculated between the treadmill stages.



While kinematics and kinetics were significantly altered at the highest slope compared to level ground when walking on a treadmill, synergies were similar across stages.



The high correlations between synergies across stages indicate that neuromuscular control strategies do not change as children walk at different speeds and slopes on a treadmill. However, the multiple significant differences in kinematics and kinetics between stages indicate real differences in movement pattern. This supports the theory that synergies are neurological in origin and not simply a response to the biomechanical task constraints.
A continuum theoretical model and finite elements simulation of bacterial flagellar filament phase transition.
Venue: Journal of biomechanics

The Bacterial flagellar filament can undergo a polymorphic phase transition in response to both mechanical and chemical variations in vitro and in vivo environments. Under mechanical stimuli, such as viscous flow or forces induced by motor rotation, the filament changes its phase from left-handed normal (N) to right-handed semi-coiled (SC) via phase nucleation and growth. Our detailed mechanical analysis of existing experiments shows that both torque and bending moment contribute to the filament phase transition. In this paper, we establish a non-convex and non-local continuum model based on the Ginzburg-Landau theory to describe main characteristics of the filament phase transition such as new-phase nucleation, growth, propagation and the merging of neighboring interfaces. The finite element method (FEM) is adopted to simulate the phase transition under a displacement-controlled loading condition (rotation angle and bending deflection). We show that new-phase nucleation corresponds to the maximum torque and bending moment at the stuck end of the filament. The hysteresis loop in the loading and unloading curves indicates energy dissipation. When the new phase grows and propagates, torque and bending moment remain static. We also find that there is a drop in load when the two interfaces merge, indicating a concomitant reduction in the interfacial energy. Finally, the interface thickness is governed by the coefficients of the gradient of order parameters in the non-local interface energy. Our continuum theory and the finite element method provide a method to study the mechanical behavior of such biomaterials.
Is bone density associated with intervertebral disc pressure in healthy and degenerated discs?
Venue: Journal of biomechanics

The coupling of the intervertebral disc (IVD) and vertebra as a biomechanical unit suggests that changes in the distribution of pressure within the IVD (intradiscal pressure, IDP) as a result of disc degeneration can influence the distribution of bone density within the vertebra, and vice versa. The goal of this study was to assess the correspondence between IDP and bone density in the adjacent vertebrae, with emphasis on how this correspondence differs between healthy and degenerated IVDs. Bone density of the endplates and subchondral bone in regions adjacent to the anterior and posterior annulus fibrosus (aAF and pAF, respectively) and nucleus pulposus (NP) was measured via quantitative computed tomography (QCT) in 61 spine segments (T7-9, T9-11, T10-12; 71±14years). IDP was measured in the aAF, NP, and pAF regions in 26 of the spine segments (68±16years) while they were tested in flexed (5°) or erect postures. Disc degeneration was assessed by multiple grading schemes. No correlation was found between bone density and IDP in either posture (p>0.104). Regional variations in IDP and, to a greater extent bone density, were found to change with advancing degeneration: both IDP (p=0.045) and bone density (p=0.024) decreased in the NP region relative to the aAF region. The finding of only a modest correspondence between degeneration-associated changes in IDP and bone density may arise from complexity in how IDP relates to mechanical force transmission through the endplate and from limitations of the available IVD grading schemes in estimating the mechanical behavior of the IVD.
Sex differences in leg dexterity are not present in elite athletes.
Venue: Journal of biomechanics

We studied whether the time-varying forces that control unstable foot-ground interactions provide insight into the neural control of dynamic leg function. Twenty elite (10F, 26.4±3.5yrs) and 20 recreational (10F, 24.8±2.4yrs) athletes used an isolated leg to maximally compress a slender spring designed to buckle at low forces while seated. The foot forces during the compression at the edge of instability quantify the maximal sensorimotor ability to control dynamic foot-ground interactions. Using the nonlinear analysis technique of attractor reconstruction, we characterized the spatial (interquartile range IQR) and geometric (trajectory length TL, volume V, and sum of edge lengths SE) features of the dynamical behavior of those force time series. ANOVA confirmed the already published effect of sex, and a new effect of athletic ability, respectively, in TL (p=0.014 and p<0.001), IQR (p=0.008 and p<0.001), V (p=0.034 and p=0.002), and SE (p=0.033 and p<0.001). Further analysis revealed that, for recreational athletes, females exhibited weaker corrective actions and greater stochasticity than males as per their greater mean values of TL (p=0.003), IQR (p=0.018), V (p=0.017), and SE (p=0.025). Importantly, sex differences disappeared in elite athletes. These results provide an empirical link between sex, athletic ability, and nonlinear dynamical control. This is a first step in understanding the sensorimotor mechanisms for control of unstable foot-ground interactions. Given that females suffer a greater incidence of non-contact knee ligament injuries, these non-invasive and practical metrics of leg dexterity may be both indicators of athletic ability, and predictors of risk of injury.
Design and verification of a simple 3D dynamic model of speed skating which mimics observed forces and motions.
Venue: Journal of biomechanics

Advice about the optimal coordination pattern for an individual speed skater, could be addressed by simulation and optimization of a biomechanical speed skating model. But before getting to this optimization approach one needs a model that can reasonably match observed behaviour. Therefore, the objective of this study is to present a verified three dimensional inverse skater model with minimal complexity, which models the speed skating motion on the straights. The model simulates the upper body transverse translation of the skater together with the forces exerted by the skates on the ice. The input of the model is the changing distance between the upper body and the skate, referred to as the leg extension (Euclidean distance in 3D space). Verification shows that the model mimics the observed forces and motions well. The model is most accurate for the position and velocity estimation (respectively 1.2% and 2.9% maximum residuals) and least accurate for the force estimations (underestimation of 4.5-10%). The model can be used to further investigate variables in the skating motion. For this, the input of the model, the leg extension, can be optimized to obtain a maximal forward velocity of the upper body.
Effect of arm swinging on lumbar spine and hip joint forces.
Venue: Journal of biomechanics

During level walking, arm swing plays a key role in improving dynamic stability. In vivo investigations with a telemeterized vertebral body replacement showed that spinal loads can be affected by differences in arm positions during sitting and standing. However, little is known about how arm swing could influence the lumbar spine and hip joint forces and motions during walking. The present study aims to provide better understanding of the contribution of the upper limbs to human gait, investigating ranges of motion and joint reaction forces. A three-dimensional motion analysis was carried out via a motion capturing system on six healthy males and five patients with hip instrumented implant. Each subject performed walking with different arm swing amplitudes (small, normal, and large) and arm positions (bound to the body, and folded across the chest). The motion data were imported in a commercial musculoskeletal analysis software for kinematic and inverse dynamic investigation. The range of motion of the thorax with respect to the pelvis and of the pelvis with respect to the ground in the transversal plane were significantly associated with arm position and swing amplitude during gait. The hip external-internal rotation range of motion statistically varied only for non-dominant limb. Unlike hip joint reaction forces, predicted peak spinal loads at T12-L1 and L5-S1 showed significant differences at approximately the time of contralateral toe off and contralateral heel strike. Therefore, arm position and swing amplitude have a relevant effect on kinematic variables and spinal loads, but not on hip loads during walking.