nils-herrmann/clean_scirepeval huggingface dataset:
Investigation on masticatory muscular functionality following oral reconstruction - An inverse identification approach.
Venue: Journal of biomechanics

The human masticatory system has received significant attention in the areas of biomechanics due to its sophisticated co-activation of a group of masticatory muscles which contribute to the fundamental oral functions. However, determination of each muscular force remains fairly challenging in vivo; the conventional data available may be inapplicable to patients who experience major oral interventions such as maxillofacial reconstruction, in which the resultant unsymmetrical anatomical structure invokes a more complex stomatognathic functioning system. Therefore, this study aimed to (1) establish an inverse identification procedure by incorporating the sequential Kriging optimization (SKO) algorithm, coupled with the patient-specific finite element analysis (FEA) in silico and occlusal force measurements at different time points over a course of rehabilitation in vivo; and (2) evaluate muscular functionality for a patient with mandibular reconstruction using a fibula free flap (FFF) procedure. The results from this study proved the hypothesis that the proposed method is of certain statistical advantage of utilizing occlusal force measurements, compared to the traditionally adopted optimality criteria approaches that are basically driven by minimizing the energy consumption of muscle systems engaged. Therefore, it is speculated that mastication may not be optimally controlled, in particular for maxillofacially reconstructed patients. For the abnormal muscular system in the patient with orofacial reconstruction, the study shows that in general, the magnitude of muscle forces fluctuates over the 28-month rehabilitation period regardless of the decreasing trend of the maximum muscular capacity. Such finding implies that the reduction of the masticatory muscle activities on the resection side might lead to non-physiological oral biomechanical responses, which can change the muscular activities for stabilizing the reconstructed mandible.
An actuated dissipative spring-mass walking model: Predicting human-like ground reaction forces and the effects of model parameters.
Venue: Journal of biomechanics

Simple models are widely used to understand the mechanics of human walking. The optimization-based minimal biped model and spring-loaded-inverted-pendulum (SLIP) model are two popular models that can achieve human-like walking patterns. However, ground reaction forces (GRF) from these two models still deviate from experimental data. In this paper, we proposed an actuated dissipative spring-mass model by integrating these two models to realize more human-like GRF patterns. We first explored the function of stiffness, damping, and weights of both energy cost and force cost in the objective function and found that these parameters have distinctly different influences on the optimized gait and GRF profiles. The stiffness and objective weight affect the number and size of peaks in the vertical GRF and stance time. The damping changes the relative size of the peaks but has little influence on stance time. Based on these observations, these parameters were manually tuned at three different speeds to approach experimentally measured vertical GRF and the highest correlation coefficient can reach 0.983. These results indicate that the stiffness, damping, and proper objective functions are all important factors in achieving human-like motion for this simple walking model. These findings can facilitate the understanding of human walking dynamics and may be applied in future biped models.
Defining gait patterns using Parallel Factor 2 (PARAFAC2): A new analysis of previously published data.
Venue: Journal of biomechanics

Three-dimensional gait analysis (3D-GA) is commonly used to answer clinical questions of the form "which joints and what variables are most affected during when". When studying high-dimensional datasets, traditional dimension reduction methods (e.g. principal components analysis) require "data flattening", which may make the ensuing solutions difficult to interpret. The aim of the present study is to present a case study of how a multi-dimensional dimension reduction technique, Parallel Factor 2 (PARAFAC2), provides a clinically interpretable set of solutions to typical biomechanical datasets where different variables are collected during walking and running. Three-dimensional kinematic and kinetic data used for the present analyses came from two publicly available datasets on walking (n = 33) and running (n = 28). For each dataset, a four-dimensional array was constructed as follows: Mode A was time normalized cycle points; mode B was the number of participants multiplied by the number of speed conditions tested; mode C was the number of joint degrees of freedom, and mode D was variable (angle, velocity, moment, power). Five factors for walking and four factors for running were extracted which explained 79.23% and 84.64% of their dataset's variance. The factor which explains the greatest variance was swing-phase sagittal plane knee kinematics (walking), and kinematics and kinetics (running). Qualitatively, all extracted factors increased in magnitude with greater speed in both walking and running. This study is a proof of concept that PARAFAC2 is useful for performing dimension reduction and producing clinically interpretable solutions to guide clinical decision making.
Initial contact and toe off event identification for rearfoot and non-rearfoot strike pattern treadmill running at different speeds.
Venue: Journal of biomechanics

The purpose of this study was to determine the validity of kinematic based initial contact (IC) and toe-off (TO) identification algorithms for rearfoot and non-rearfoot runners across a broad range of treadmill running speeds. 14 healthy active participants completed six 20-60 s treadmill running trials at 6 speeds: 2.24, 2.68, 3.13, 3.58, 4.02, and 4.48 ms-1. 3D kinematic data were collected for the last 20 s of each trial. Force plates (FP) were used as the gold standard to determine ICFP and TOFP for each step. Three algorithms for finding IC, ICMilner, ICAlvim, ICAlvim-mod, and one algorithm for finding toe off, TOFellin, were chosen for analysis. Root mean square errors (RMSE) and difference scores with 95% confidence intervals were computed for IC, TO and stance time (ST). ICAlvim RMSE ranged from 0.175 to 0.219 s. STAlvim RMSE ranged from 0.168 to 0.216 s. ICAlvim-mod RMSE ranged from 0.105 to 0.131 s. STAlvim-mod RMSE ranged from 0.108 to 0.129 s. ICMilner RMSE ranged 0.012 to 0.015 s. STMilner RMSE ranged 0.019 to 0.024 s. ICMilner accuracy was inversely related to speed. ICMilner corrected with a linear regression equation reduced differences to- 0.006 ± 0.012 s with 86% of foot strikes identified within 20 ms and 58% with 10 ms. TOFellin RMSE ranged from 0.012 to 0.016 s. ICMilner adjusted for speed and TOFellin can be used to predict IC and TO within a broad range of treadmill running speeds (2.24-4.48 ms-1) and for rearfoot and non-rearfoot strikers.
In vivo assessment of the interaction of patellar tendon tibial shaft angle and anterior cruciate ligament elongation during flexion.
Venue: Journal of biomechanics

A potential cause of non-contact anterior cruciate ligament (ACL) injury is landing on an extended knee. In line with this hypothesis, studies have shown that the ACL is elongated with decreasing knee flexion angle. Furthermore, at low flexion angles the patellar tendon is oriented to increase the anterior shear component of force acting on the tibia. This indicates that knee extension represents a position in which the ACL is taut, and thus may have an increased propensity for injury, particularly in the presence of excessive force acting via the patellar tendon. However, there is very little in vivo data to describe how patellar tendon orientation and ACL elongation interact during flexion. Therefore, this study measured the patellar tendon tibial shaft angle (indicative of the relative magnitude of the shear component of force acting via the patellar tendon) and ACL length in vivo as subjects performed a quasi-static lunge at varying knee flexion angles. Spearman rho rank correlations within each individual revealed that flexion angles were inversely correlated to both ACL length (rho = -0.94 ± 0.07, mean ± standard deviation, p < 0.05) and patellar tendon tibial shaft angle (rho = -0.99 ± 0.01, p < 0.05). These findings indicate that when the knee is extended, the ACL is both elongated and the patellar tendon tibial shaft angle is increased, resulting in a relative increase in anterior shear force on the tibia acting via the patellar tendon. Therefore, these data support the hypothesis that landing with the knee in extension is a high risk scenario for ACL injury.
Improved biomechanical metrics of cerebral vasospasm identified via sensitivity analysis of a 1D cerebral circulation model.
Venue: Journal of biomechanics

Cerebral vasospasm (CVS) is a life-threatening condition that occurs in a large proportion of those affected by subarachnoid haemorrhage and stroke. CVS manifests itself as the progressive narrowing of intracranial arteries. It is usually diagnosed using Doppler ultrasound, which quantifies blood velocity changes in the affected vessels, but has low sensitivity when CVS affects the peripheral vasculature. The aim of this study was to identify alternative biomarkers that could be used to diagnose CVS. We used a 1D modelling approach to describe the properties of pulse waves that propagate through the cardiovascular system, which allowed the effects of different types of vasospasm on waveforms to be characterised at several locations within a simulated cerebral network. A sensitivity analysis empowered by the use of a Gaussian process statistical emulator was used to identify waveform features that may have strong correlations with vasospasm. We showed that the minimum rate of velocity change can be much more effective than blood velocity for stratifying typical manifestations of vasospasm and its progression. The results and methodology of this study have the potential not only to improve the diagnosis and monitoring of vasospasm, but also to be used in the diagnosis of many other cardiovascular diseases where cardiovascular waves can be decoded to provide disease characterisation.
Analysis of forearm rotational motion using biplane fluoroscopic intensity-based 2D-3D matching.
Venue: Journal of biomechanics

Measuring three-dimensional (3D) forearm rotational motion is difficult. We aimed to develop and validate a new method for analyzing 3D forearm rotational motion. We proposed biplane fluoroscopic intensity-based 2D-3D matching, which employs automatic registration processing using the evolutionary optimization strategy. Biplane fluoroscopy was conducted for forearm rotation at 12.5 frames per second along with computed tomography (CT) at one static position. An arm phantom was embedded with eight stainless steel spheres (diameter, 1.5 mm), and forearm rotational motion measurements using the proposed method were compared with those using radiostereometric analysis, which is considered the ground truth. As for the time resolution analysis, we measured radiohumeral joint motion in a patient with posterolateral rotatory instability and compared the 2D-3D matching method with the simulated multiple CT method, which uses CTs at multiple positions and interpolates between the positions. Rotation errors of the radius and ulna between these two methods were 0.31 ± 0.35° and 0.32 ± 0.33°, respectively, translation errors were 0.43 ± 0.35 mm and 0.29 ± 0.25 mm, respectively. Although the 2D-3D method could detect joint dislocation, the multiple CT method could not detect quick motion during joint dislocation. The proposed method enabled high temporal- and spatial-resolution motion analyses with low radiation exposure. Moreover, it enabled the detection of a sudden motion, such as joint dislocation, and may contribute to 3D motion analysis, including joint dislocation, which currently cannot be analyzed using conventional methods.
Does stretching velocity affect residual force enhancement?
Venue: Journal of biomechanics

It is thought that the magnitude of residual force enhancement (RFE) is not affected by stretch velocity. However, the range of stretch velocities studied in previous investigations has been limited to slow and moderate velocities. High velocities of muscle stretching are associated with a loss of force and incomplete cross-bridge attachment to actin, thus creating a unique set of eccentric conditions referred to as slippage. The purpose of this study was to extend the relationship between stretch velocity and RFE to high velocities. We hypothesized that slippage at high velocities might affect RFE. We stretched cat soleus muscles for 4 mm to the plateau of the force-length relationship at speeds of 2, 4, 8, 16, 32, 64 mm/s to induce RFE, and slippage for the fastest condition. For each RFE test, a corresponding isometric reference test was conducted. Residual force enhancement was quantified as the relative increase in isometric steady state force between the experimental stretch and the isometric reference tests. Residual force enhancement was similar for all stretch speeds, as expected, with the exception of the fastest speed (64 mm/s), which was associated with slippage and no significant RFE. These results suggest that if stretch speeds are too fast, and are associated with slippage, RFE is abolished. We conclude from these findings that proper cross-bridge engagement is required during eccentric muscle action to produce RFE.
Effects of passive Bi-axial ankle stretching while walking on uneven terrains in older adults with chronic stroke.
Venue: Journal of biomechanics

Many people with stroke experience foot drop while walking. Further, walking on uneven surfaces is a common fall risk for these people that hinder with their daily life activities. In addition, a few years after a stroke, lower-limb exercises become less focused, especially the ankle joint movement. The objective of this study is to determine the gait performance of older adults with chronic stroke on an uneven surface in relation to ankle mobility after a four-week bi-axial ankle range of motion (ROM) exercise session. Fifteen older adults with chronic post-stroke hemiparesis (N = 15; mean age = 65 years) participated in a total of 12 bi-axial ankle ROM exercises that consisted of three 30-min training sessions per week for four weeks. Basic clinical tests and gait performance in even and uneven surfaces were evaluated before and after training. Participants with chronic post-stroke hemiparesis showed significantly improved ankle functions, decreased ankle stiffness (from 0.140 ± 0.059 to 0.128 ± 0.067 N·m/°; p = 0.025), and increased paretic ankle passive ROMs (dorsiflexion(DF)/plantarflexion(PF): from 27.3 ± 14.7° to 50.6 ± 10.3°, p < 0.001; inversion(INV)/eversion(EV): 21.7 ± 9.7° to 28.6 ± 9.9°; p = 0.033) after training. They exhibited significant improvements in the walking performance over an uneven surface, step kinematics (walking speed 0.257 ± 0.17 to 0.320 ± 0.178 m/s; p = 0.017; step length: 0.214 ± 0.109 to 0.243 ± 0.108 m; p = 0.009), and clinical balance and mobility (Berg balance scale: 47.2 ± 4.7 to 50.1 ± 3.9, p = 0.0001; timed-up and go test: 23.9 ± 10.3 to 20.2 ± 7.0 s, p = 0.0156). This study is the first research to investigate the walking performance on uneven surfaces in the elderly with chronic stroke in relation to the ankle biomechanical property changes.
A preliminary muscle activity analysis: Handle based and push-rim wheelchair propulsion.
Venue: Journal of biomechanics

Approximately ninety percent of the wheelchair users worldwide prefer the conventional push rim mode of propulsion for daily mobility and rehabilitation. Even though push-rim wheelchairs help to promote a healthy life style, the high muscular demand and the non-continuous push motions can lead to serious upper extremity injuries. In this study, muscle EMG data of ten healthy subjects were recorded for a newly introduced handle based propulsion mechanism (HBP) and compared to conventional push-rim propulsion at two workloads, 25 W and 35 W respectively. The results for the mean peak muscle activations at both workloads demonstrate that push-rim propulsion leads to higher peak muscle activity compared to HBP at a similar wheelchair forward velocity of 1.11 m/s. The generation of these high peak muscle activations with increasing loads in push-rim propulsion over time can lead to overuse injuries. Overall, the use of the HBP mechanism is less straining to the muscles and may reduce fatigue during prolonged propulsion.
Role of diastolic vortices in flow and energy dynamics during systolic ejection.
Venue: Journal of biomechanics

MRI-based computational fluid dynamics simulations were performed in the left ventricles of two adult porcine subjects with varying physiological states (before and after an induced infarction). The hypothesis that diastolic vortices store kinetic energy and assist systolic ejection was tested, by performing systolic simulations in the presence and absence of diastolic vortices. The latter was achieved by reinitializing the entire velocity field to be zero at the beginning of systole. A rudimentary prescribed motion model of a mitral valve was included in the simulations to direct the incoming mitral jet towards the apex. Results showed that the presence or absence of diastolic vortex rings had insignificant impact on the energy expended by walls of the left ventricles for systolic ejection for both the porcine subjects, under all physiological conditions. Although substantial kinetic energy was stored in diastolic vortices by end diastole, it provided no appreciable savings during systolic ejection, and most likely continued to complete dissipation during systole. The role of diastolic vortices in apical washout was investigated by studying the cumulative mass fraction of passive dye that was ejected during systole in the presence and absence of vortices. Results indicated that the diastolic vortices play a crucial role in ensuring efficient washout of apical blood during systolic ejection.
Effects of hip torque during step-to-step transition on center-of-mass dynamics during human walking examined with numerical simulation.
Venue: Journal of biomechanics

Besides the leg force actuator, humans also use a hip torque actuator during the step-to-step transition to redirect the velocity of CoM (Center of Mass). Although the leg force actuator has been widely studied, few researches analyze the hip torque actuator during the step-to-step transition. In this paper, we build a powered walking model which consists of a point mass linked with two compliant legs. Each leg has a spring and a damper in parallel. Two types of active actuators, the force actuator on the leg and the torque actuator at the hip, are added to simulate the leg force and hip torque actuator during the step-to-step transition. The cycle walk is solved by numerical simulations under different hip torque strength, and the energetics and stability are evaluated. The simulation results show that the hip torque actuator can reduce the energy cost and improve the stability of walking. Further analysis shows that the hip torque actuator can reduce mechanical works of both legs with small extra energy cost. To understand the principle of hip torque actuator, the CoM dynamics is analyzed. It is shown that the hip torque actuator is efficient on the redirection of CoM. Thus, it can improve the stability and reduce required forces of both legs, which decreases the energy cost. Our work provides a fundamental understanding of the hip torque during the step-to-step transition, and may help improve the design of bipedal robots and prosthesis.
Adaptive changes in micromechanical environments of cancellous and cortical bone in response to in vivo loading and disuse.
Venue: Journal of biomechanics

The skeleton accommodates changes in mechanical environments by increasing bone mass under increased loads and decreasing bone mass under disuse. However, little is known about the adaptive changes in micromechanical behavior of cancellous and cortical tissues resulting from loading or disuse. To address this issue, in vivo tibial loading and hindlimb unloading experiments were conducted on 16-week-old female C57BL/6J mice. Changes in bone mass and tissue-level strains in the metaphyseal cancellous and midshaft cortical bone of the tibiae, resulting from loading or unloading, were determined using microCT and finite element (FE) analysis, respectively. We found that loading- and unloading-induced changes in bone mass were more pronounced in the cancellous than cortical bone. Simulated FE-loading showed that a greater proportion of elements experienced relatively lower longitudinal strains following load-induced bone adaptation, while the opposite was true in the disuse model. While the magnitudes of maximum or minimum principal strains in the metaphyseal cancellous and midshaft cortical bone were not affected by loading, strains oriented with the long axis were reduced in the load-adapted tibia suggesting that loading-induced micromechanical benefits were aligned primarily in the loading direction. Regression analyses demonstrated that bone mass was a good predictor of bone tissue strains for the cortical bone but not for the cancellous bone, which has complex microarchitecture and spatially-variant strain environments. In summary, loading-induced micromechanical benefits for cancellous and cortical tissues are received primarily in the direction of force application and cancellous bone mass may not be related to the micromechanics of cancellous bone.
Validity and reliability of a shoe-embedded sensor module for measuring foot progression angle during over-ground walking.
Venue: Journal of biomechanics

Wearable systems are becoming increasingly popular for gait assessment outside of laboratory settings. A single shoe-embedded sensor module can measure the foot progression angle (FPA) during walking. The FPA has important clinical utility, particularly in populations with knee osteoarthritis, as it is a target for biomechanical treatments. However, the validity and the day-to-day reliability of FPA measurement using wearable systems during over-ground walking has yet to be established. Two gait analysis sessions on 20 healthy adults were conducted. During both sessions, participants performed natural over-ground walking in a motion capture laboratory and on a 100 m linear section of outdoor athletics track. FPA was measured in the laboratory via marker trajectory data, while the sensor module measured FPA during the outdoor track walking. Validity was examined by comparing the laboratory- and sensor-measured average FPA. Day-to-day reliability was examined by comparing the sensor-measured FPA between the first and second gait analysis sessions. Average absolute error between motion capture and sensor measured FPA were 1.7° and 2.1° at session 1 and 2, respectively. A Bland and Altman plot indicated no systematic bias, with 95% limit of agreement widths of 4.2° - 5.1°. Intraclass correlation coefficient (ICC2k) analysis resulted in good to excellent validity (ICC = 0.89 - 0.91) and reliability (ICC = 0.95). Overall, the shoe-embedded sensor module is a valid and reliable method of measuring FPA during over-ground walking without the need for laboratory equipment.