nils-herrmann/clean_scirepeval huggingface dataset:
Explosive lower limb extension mechanics: An on-land vs. in-water exploratory comparison.
Venue: Journal of biomechanics

During a horizontal underwater push-off, performance is strongly limited by the presence of water, inducing resistances due to its dense and viscous nature. At the same time, aquatic environments offer a support to the swimmer with the hydrostatic buoyancy counteracting the effects of gravity. Squat jump is a vertical terrestrial push-off with a maximal lower limb extension limited by the gravity force, which attracts the body to the ground. Following this observation, we characterized the effects of environment (water vs. air) on the mechanical characteristics of the leg push-off. Underwater horizontal wall push-off and vertical on-land squat jumps of two local swimmers were evaluated with force plates, synchronized with a lateral camera. To better understand the resistances of the aquatic movement, a quasi-steady Computational Fluid Dynamics (CFD) analysis was performed. The force-, velocity- and power-time curves presented similarities in both environments corresponding to a proximo-distal joints organization. In water, swimmers developed a three-step explosive rise of force, which the first one mainly related to the initiation of body movement. Drag increase, which was observed from the beginning to the end of the push-off, related to the continuous increase of body velocity with high values of drag coefficient (CD) and frontal areas before take-off. Specifically, with velocity, frontal area was the main drag component to explain inter-individual differences, suggesting that the streamlined position of the lower limbs is decisive to perform an efficient push-off. This study motivates future CFD simulations under more ecological, unsteady conditions.
New parameters describing how knee ligaments carry force in situ predict interspecimen variations in laxity during simulated clinical exams.
Venue: Journal of biomechanics

Knee laxity, defined as the net translation or rotation of the tibia relative to the femur in a given direction in response to an applied load, is highly variable from person to person. High levels of knee laxity as assessed during routine clinical exams are associated with first-time ligament injury and graft reinjury following reconstruction. During laxity exams, ligaments carry force to resist the applied load; however, relationships between intersubject variations in knee laxity and variations in how ligaments carry force as the knee moves through its passive envelope of motion, which we refer to as ligament engagement, are not well established. Thus, the objectives of this study were, first, to define parameters describing ligament engagement and, then, to link variations in ligament engagement and variations in laxity across a group of knees. We used a robotic manipulator in a cadaveric knee model (n=20) to quantify how important knee stabilizers, namely the anterior and posterior cruciate ligaments (ACL and PCL, respectively), as well as the medial collateral ligament (MCL) engage during respective tests of anterior, posterior, and valgus laxity. Ligament engagement was quantified using three parameters: (1) in situ slack, defined as the relative tibiofemoral motion from the neutral position of the joint to the position where the ligament began to carry force; (2) in situ stiffness, defined as the slope of the linear portion of the ligament force-tibial motion response; and (3) ligament force at the peak applied load. Knee laxity was related to parameters of ligament engagement using univariate and multivariate regression models. Variations in the in situ slack of the ACL and PCL predicted anterior and posterior laxity, while variations in both in situ slack and in situ stiffness of the MCL predicted valgus laxity. Parameters of ligament engagement may be useful to further characterize the in situ biomechanical function of ligaments and ligament grafts.
Gait biomechanics of skipping are substantially different than those of running.
Venue: Journal of biomechanics

The inherit injury risk associated with high-impact exercises calls for alternative ways to achieve the benefits of aerobic exercise while minimizing excessive stresses to body tissues. Skipping presents such an alternative, incorporating double support, flight, and single support phases. We used ground reaction forces (GRFs), lower extremity joint torques and powers to compare skipping and running in 20 healthy adults. The two consecutive skipping steps on each limb differed significantly from each other, and from running. Running had the longest step length, the highest peak vertical GRF, peak knee extensor torque, and peak knee negative and positive power and negative and positive work. Skipping had the greater cadence, peak horizontal GRF, peak hip and ankle extensor torques, peak ankle negative power and work, and peak ankle positive power. The second vs first skipping step had the shorter step length, higher cadence, peak horizontal GRF, peak ankle extensor torque, and peak ankle negative power, negative work, and positive power and positive work. The first skipping step utilized predominately net negative joint work (eccentric muscle action) while the second utilized predominately net positive joint work (concentric muscle action). The skipping data further highlight the persistence of net negative work performed at the knee and net positive work performed at the ankle across locomotion gaits. Evidence of step segregation was seen in distribution of the braking and propelling impulses and net work produced across the hip, knee, and ankle joints.



Skipping was substantially different than running and was temporally and spatially asymmetrical with successive foot falls partitioned into a dominant function, either braking or propelling whereas running had a single, repeated step in which both braking and propelling actions were performed equally.
Strategy quantification using body worn inertial sensors in a reactive agility task.
Venue: Journal of biomechanics

Agility performance is often evaluated using time-based metrics, which provide little information about which factors aid or limit success. The objective of this study was to better understand agility strategy by identifying biomechanical metrics that were sensitive to performance speed, which were calculated with data from an array of body-worn inertial sensors. Five metrics were defined (normalized number of foot contacts, stride length variance, arm swing variance, mean normalized stride frequency, and number of body rotations) that corresponded to agility terms defined by experts working in athletic, clinical, and military environments. Eighteen participants donned 13 sensors to complete a reactive agility task, which involved navigating a set of cones in response to a vocal cue. Participants were grouped into fast, medium, and slow performance based on their completion time. Participants in the fast group had the smallest number of foot contacts (normalizing by height), highest stride length variance (normalizing by height), highest forearm angular velocity variance, and highest stride frequency (normalizing by height). The number of body rotations was not sensitive to speed and may have been determined by hand and foot dominance while completing the agility task. The results of this study have the potential to inform the development of a composite agility score constructed from the list of significant metrics. By quantifying the agility terms previously defined by expert evaluators through an agility score, this study can assist in strategy development for training and rehabilitation across athletic, clinical, and military domains.
Micromechanical analysis of brain's diffuse axonal injury.
Venue: Journal of biomechanics

Computational models are important tools which help researchers understand traumatic brain injury (TBI). A mechanistic multi-scale numerical approach is introduced to quantify diffuse axonal injury (DAI), the most important mechanism of TBI, induced by a mechanical insult at micro-scale regions of the white matter or voxels where fiber orientations are the same. Using the mechanical properties of a single axon with a viscoelastic constitutive relation and its functional failure in terms of electrophysiological impairment, a numerical 2D micro-level lattice method is implemented to directly analyze the percentage of injured axons in a voxel containing a bundle of axons all with the same orientation under biaxial stretches. Reference micro-injury maps are then developed with the input parameters based on the principal strain or stretch values and their direction with respect to axons, which provide the percentage of injured axons in the voxel of interest as the output. The methodology is independent of any statistical analyses of the accident data and medical reports to derive probabilistic injury risk curves for DAI. Avoiding a structurally detailed full finite element head model, this study proposes a micro-mechanical approach which considers the anatomical structure of neural axons in the white matter together with their mechanical properties using a numerical lattice method to analyze the brain's diffuse axonal injury. This work has the potential to help develop safer prevention tools and more effective diagnosis methods for DAI.
Reliability of assessing postural control during seated balancing using a physical human-robot interaction.
Venue: Journal of biomechanics

This study evaluated the within- and between-visit reliability of a seated balance test for quantifying trunk motor control using input-output data. Thirty healthy subjects performed a seated balance test under three conditions: eyes open (EO), eyes closed (EC), and eyes closed with vibration to the lumbar muscles (VIB). Each subject performed three trials of each condition on three different visits. The seated balance test utilized a torque-controlled robotic seat, which together with a sitting subject resulted in a physical human-robot interaction (pHRI) (two degrees-of-freedom with upper and lower body rotations). Subjects balanced the pHRI by controlling trunk rotation in response to pseudorandom torque perturbations applied to the seat in the coronal plane. Performance error was expressed as the root mean square (RMSE) of deviations from the upright position in the time domain and as the mean bandpass signal energy (Emb) in the frequency domain. Intra-class correlation coefficients (ICC) quantified the between-visit reliability of both RMSE and Emb. The empirical transfer function estimates (ETFE) from the perturbation input to each of the two rotational outputs were calculated. Coefficients of multiple correlation (CMC) quantified the within- and between-visit reliability of the averaged ETFE. ICCs of RMSE and Emb for all conditions were ≥0.84. The mean within- and between-visit CMCs were all ≥0.96 for the lower body rotation and ≥0.89 for the upper body rotation. Therefore, our seated balance test consisting of pHRI to assess coronal plane trunk motor control is reliable.
Integrated geometric and mechanical analysis of an image-based lymphatic valve.
Venue: Journal of biomechanics

Lymphatic valves facilitate the lymphatic system's role in maintaining fluid homeostasis. Malformed valves are found in several forms of primary lymphœdema, resulting in incurable swelling of the tissues and immune dysfunction. Their experimental study is complicated by their small size and operation in low pressure and low Reynolds number environments. Mathematical models of these structures can give insight and complement experimentation. In this work, we present the first valve geometry reconstructed from confocal imagery and used in the construction of a subject-specific model in a closing mode. A framework is proposed whereby an image is converted into a valve model. An FEA study was performed to identify the significance of the shear modulus, the consequences of smoothing the leaflet surface and the effect of wall motion on valve behaviour. Smoothing is inherent to any analysis from imagery. The nature of the image, segmentation and meshing all cause attenuation of high-frequency features. Smoothing not only causes loss of surface area but also the loss of high-frequency geometric features which may reduce stiffness. This work aimed to consider these effects and inform studies by taking a manual reconstruction and through manifold harmonic analysis, attenuating higher frequency features to replicate lower resolution images or lower degree-of-freedom reconstructions. In conclusion, two metrics were considered: trans-valvular pressure required to close the valve, ΔPc, and the retrograde volume displacement after closure. The higher ΔPc, the greater the volume of lymph that will pass through the valve during closure. Retrograde volume displacement after closure gives a metric of compliance of the valve and for the quality of the valve seal. In the case of the image-specific reconstructed valve, removing features with a wavelength longer than four μm caused changes in ΔPc. Varying the shear modulus from 10 kPa to 60 kPa caused a 3.85-fold increase in the retrograde volume displaced. The inclusion of a non-rigid wall caused ΔPc to increase from 1.56 to 2.52 cmH2O.
An in vivo model of a mechanically-induced bone marrow lesion.
Venue: Journal of biomechanics

Bone marrow lesions (BMLs) are radiologic abnormalities in magnetic resonance images of subchondral bone that are correlated with osteoarthritis. Little is known about the physiologic processes within a BML, although BMLs are associated with mechanical stress, bone tissue microdamage and increased bone remodeling. Here we establish a rabbit model to study the pathophysiology of BMLs. We hypothesized that in vivo loads that generate microdamage in cancellous bone would also create BMLs and increase bone remodeling. In vivo cyclic loading (0.2-2.0 MPa in compression for 10,000 cycles at 2 Hz) was applied to epiphyseal cancellous bone in the distal femurs of New Zealand white rabbits (n=3, right limb loaded, left limb controls experienced surgery but no loading). Magnetic resonance images were collected using short tau inversion recovery (STIR) and T1 weighted sequences at 1 and 2 weeks after surgery/loading and histological analysis of the BML was performed after euthanasia to examine tissue microdamage and remodeling. Loaded limbs displayed BMLs while control limbs showed only a small BML-like signal caused by surgery. Histological analysis of the BML at 2 weeks after loading showed increased tissue microdamage (p=0.03) and bone resorption (p=0.01) as compared to controls. The model described here displays the hallmarks of load-induced BMLs, supporting the use of the model to examine changes in bone during the development, progression and treatment of BMLs.
Inter- and intra-rater reliability of 3D kinematics during maximum mouth opening of asymptomatic subjects.
Venue: Journal of biomechanics

Previous studies evaluated 3D human jaw movements using kinematic analysis systems during mouth opening, but information on the reliability of such measurements is still scarce. The purpose of this study was to analyze within- and between-session reliabilities, inter-rater reliability, standard error of measurement (SEM), minimum detectable change (MDC) and consistency of agreement across raters and sessions of 3D kinematic variables during maximum mouth opening (MMO). Thirty-six asymptomatic subjects from both genders were evaluated on two different days, five to seven days apart. Subjects performed three MMO movements while kinematic data were collected. Intraclass correlation coefficient (ICC), SEM and MDC were calculated for all variables, and Bland-Altman plots were constructed. Jaw radius and width were the most reproducible variables (ICC>0.81) and demonstrated minor error. Incisor displacement during MMO and angular movements in the sagittal plane presented good reliability (ICC from 0.61 to 0.8) and small errors and, consequently, could be used in future studies with the same methodology and population. The variables with smaller amplitudes (condylar translations during mouth opening and closing and mandibular movements on the frontal and transversal planes) were less reliable (ICC<0.61) and presented larger SEM and MDC. Although ICC, SEM and MDC showed less between-session reproducibility than within-session and inter-rater, the limits of agreement were larger in inter-rater comparisons. In future studies care must be taken with variables collected on different days and with mandibular movements in the frontal and transversal planes.
Estimating the L5S1 flexion/extension moment in symmetrical lifting using a simplified ambulatory measurement system.
Venue: Journal of biomechanics

Mechanical loading of the spine has been shown to be an important risk factor for the development of low-back pain. Inertial motion capture (IMC) systems might allow measuring lumbar moments in realistic working conditions, and thus support evaluation of measures to reduce mechanical loading. As the number of sensors limits applicability, the objective of this study was to investigate the effect of the number of sensors on estimates of L5S1 moments. Hand forces, ground reaction forces (GRF) and full-body kinematics were measured using a gold standard (GS) laboratory setup. In the ambulatory setup, hand forces were estimated based on the force plates measured GRF and body kinematics that were measured using (subsets of) an IMC system. Using top-down inverse dynamics, L5S1 flexion/extension moments were calculated. RMSerrors (Nm) were lowest (16.6) with the full set of 17 sensors and increased to 20.5, 22 and 30.6, for 8, 6 and 4 sensors. Absolute errors in peak moments (Nm) ranged from 17.7 to 16.4, 16.9 and 49.3 Nm, for IMC setup's with 17, 8, 6 and 4 sensors, respectively. When horizontal GRF were neglected for 6 sensors, RMSerrors and peak moment errors decreased from 22 to 17.3 and from 16.9 to 13 Nm, respectively. In conclusion, while reasonable moment estimates can be obtained with 6 sensors, omitting the forearm sensors led to unacceptable errors. Furthermore, vertical GRF information is sufficient to estimate L5S1 moments in lifting.
A "sweet-spot" for fluid-induced oscillations in the conditioning of stem cell-based engineered heart valve tissues.
Venue: Journal of biomechanics

Fluid-induced shear stresses are involved in the development of cardiovascular tissues. In a tissue engineering framework, this stimulus has also been considered as a mechanical regulator of stem cell differentiation. We recently demonstrated that the fluid-oscillating effect in combination with a physiologically-relevant shear stress magnitude contributes to the formation of stem cell-derived de novo heart valve tissues. However, the range of oscillations necessary to induce favorable gene expression and engineered tissue formation is unknown. In this study, we took a computational approach to establish a range of oscillatory shear stresses that may optimize in vitro valvular tissue growth. Taking a biomimetic approach, three physiologically-relevant flow waveforms from the human: (i) aorta, (ii) pulmonary artery and (iii) superior vena cava were utilized to simulate pulsatile flow conditions within a bioreactor that housed 3 tissue specimens. Results were compared to non-physiological pulsatile flow (NPPF) and cyclic flexure-steady flow (Flex-Flow) conditions. The oscillatory shear index (OSI) was used to quantify the fluid-induced oscillations occurring on the specimen surfaces. The range of mean OSI under the physiological conditions investigated was found to be 0.18 ≤ OSI ≤ 0.23. On the other hand, NPPF and Flex-Flow environments yielded a mean OSI of 0.37 and 0.11 respectively, which were 46% higher and 45% lower than physiological conditions. Moreover, we subsequently conducted OSI-based human bone marrow stem cell (HBMSC) culture experiments which resulted in preferential valvular gene expression and phenotype (significant upregulation of BMP, KLF2A, CD31 and α-SMA using an OSI of 0.23 in comparison to a lower OSI of 0.10 or a higher OSI of 0.38; p < .05). These findings suggest that a distinct range or a "sweet-spot" for physiological OSI exists in the mechanical conditioning of tissue engineered heart valves grown from stem cell sources. We conclude that in vitro heart valve matrix development could be further enhanced by simultaneous exposure of the engineered tissues to physiologically-relevant magnitudes of both fluid-induced oscillations and shear stresses.
In vivo repeatability of the pulse wave inverse problem in human carotid arteries.
Venue: Journal of biomechanics

Accurate arterial stiffness measurement would improve diagnosis and monitoring for many diseases. Atherosclerotic plaques and aneurysms are expected to involve focal changes in vessel wall properties; therefore, a method to image the stiffness variation would be a valuable clinical tool. The pulse wave inverse problem (PWIP) fits unknown parameters from a computational model of arterial pulse wave propagation to ultrasound-based measurements of vessel wall displacements by minimizing the difference between the model and measured displacements. The PWIP has been validated in phantoms, and this study presents the first in vivo demonstration. The common carotid arteries of five healthy volunteers were imaged five times in a single session with repositioning of the probe and subject between each scan. The 1D finite difference computational model used in the PWIP spanned from the start of the transducer to the carotid bifurcation, where a resistance outlet boundary condition was applied to approximately model the downstream reflection of the pulse wave. Unknown parameters that were estimated by the PWIP included a 10-segment linear piecewise compliance distribution and 16 discrete cosine transformation coefficients for each of the inlet boundary conditions. Input data was selected to include pulse waves resulting from the primary pulse and dicrotic notch. The recovered compliance maps indicate that the compliance increases close to the bifurcation, and the variability of the average pulse wave velocity estimated through the PWIP is on the order of 11%, which is similar to that of the conventional processing technique which tracks the wavefront arrival time (13%).
Nonlinear vibrations of pre- and post-buckled lipid supramolecular micro/nano-tubules via nonlocal strain gradient elasticity theory.
Venue: Journal of biomechanics

The unique geometry with high surface ratio makes lipid micro/nano-tubules as an excellent self-assembled supramolecular structure in various biological applications such as controllable release systems and drug delivery. In the present study, the size-dependent nonlinear vibrations of axially loaded lipid micro/nano tubules associated with the both prebuckling and postbuckling domains are explored comprehensively. To accomplish this purpose, the nonlocal strain gradient theory of elasticity including simultaneously two entirely different features of size dependency is utilized within the framework of the third-order shear deformable beam model. With the aid of Hamilton's principle, the non-classical governing differential equations of motion are established incorporating the nonlinear prebuckling deformations and the large postbuckling deflections. At the end, the Galerkin method in conjunction with an improved perturbation technique is employed to initiate explicit analytical expressions for nonlocal strain gradient nonlinear frequency of pre- and post-buckled lipid micro/nano-tubules. It is seen that by taking the nonlocal size effect into consideration, the influence of geometrical parameters of the lipid micro/nano-tubule on the nonlinear vibration characteristics within the both prebuckling and postbuckling domains decreases and the frequency-deflection curves are more close to each other. However, the strain gradient size dependency has an opposite effect and leads to increase the gap between the frequency-deflection curves of axially compressed lipid micro/nano-tubules with different geometrical parameters.