nils-herrmann/clean_scirepeval huggingface dataset:
Aerodynamic study of time-trial helmets in cycling racing using CFD analysis.
Venue: Journal of biomechanics

The aerodynamic drag of three different time-trial cycling helmets was analyzed numerically for two different cyclist head positions. Computational Fluid Dynamics (CFD) methods were used to investigate the detailed airflow patterns around the cyclist for a constant velocity of 15 m/s without wind. The CFD simulations have focused on the aerodynamic drag effects in terms of wall shear stress maps and pressure coefficient distributions on the cyclist/helmet system. For a given head position, the helmet shape, by itself, obtained a weak effect on a cyclist's aerodynamic performance (<1.5%). However, by varying head position, a cyclist significantly influences aerodynamic performance; the maximum difference between both positions being about 6.4%. CFD results have also shown that both helmet shape and head position significantly influence drag forces, pressure and wall shear stress distributions on the whole cyclist's body due to the change in the near-wake behavior and in location of corresponding separation and attachment areas around the cyclist.
Hydrostatic pressure as epigenetic modulator in chondrocyte cultures: A study on miRNA-155, miRNA-181a and miRNA-223 expression levels.
Venue: Journal of biomechanics

Mechanical stimuli and hydrostatic pressure (HP) play an important role in the regulation of chondrocytes metabolism. Growing evidence demonstrated the ability of mechanical loading to modulate the expression of microRNA (miRNA) involved in chondrocytes homeostasis and in the pathogenesis of osteoarthritis (OA). The expression of miR-155, miR-181a and miR-223 in normal and OA chondrocyte cultures, and their potential modifications following exposure to three hours of a cyclic HP (1-5 MPa, frequency 0.25 Hz) were investigated. Also evaluated the expression of Chuk, regulator of the NF-kB pathway activation, which is a target gene of miR-223, was evaluated. Chondrocytes were collected immediately after pressurization (T0), and following 12, 24, and 48 h. Total RNA was extracted, reverse transcribed and used for real-time PCR. At basal condition, a significant increase of miR-155 and miR-181a was observed in OA in comparison to normal cells; on the contrary, no differences in miR-223 and Chuk expression levels were detected between normal and OA chondrocytes. miR-155 and miR-181a resulted significantly downregulated immediately after pressurization (T0) in OA cells. The pressure effect on miR-155 and miR-181a levels was maintained over time. No modifications of miR-223 were observed in response to HP, while Chuk levels resulted significantly reduced at T0 and after 12 h. Pressurization did not cause any modifications in normal cells. In conclusion, HP was able to modulate the expression of miRNA associated to OA pathogenesis. The preliminary results about Chuk response to pressure raised interest in its involvement in the possible HP induced NF-kB pathway modulation.
Biomechanical mechanism of lateral trunk lean gait for knee osteoarthritis patients.
Venue: Journal of biomechanics

The biomechanical mechanism of lateral trunk lean gait employed to reduce external knee adduction moment (KAM) for knee osteoarthritis (OA) patients is not well known. This mechanism may relate to the center of mass (COM) motion. Moreover, lateral trunk lean gait may affect motor control of the COM displacement. Uncontrolled manifold (UCM) analysis is an evaluation index used to understand motor control and variability of the motor task. Here we aimed to clarify the biomechanical mechanism to reduce KAM during lateral trunk lean gait and how motor variability controls the COM displacement. Twenty knee OA patients walked under two conditions: normal and lateral trunk lean gait conditions. UCM analysis was performed with respect to the COM displacement in the frontal plane. We also determined how the variability is structured with regards to the COM displacement as a performance variable. The peak KAM under lateral trunk lean gait was lower than that under normal gait. The reduced peak KAM observed was accompanied by medially shifted knee joint center, shortened distance of the center of pressure to knee joint center, and shortened distance of the knee-ground reaction force lever arm during the stance phase. Knee OA patients with lateral trunk lean gait could maintain kinematic synergy by utilizing greater segmental configuration variance to the performance variable. However, the COM displacement variability of lateral trunk lean gait was larger than that of normal gait. Our findings may provide clinical insights to effectively evaluate and prescribe gait modification training for knee OA patients.
Supervised learning techniques and their ability to classify a change of direction task strategy using kinematic and kinetic features.
Venue: Journal of biomechanics

This study examines the ability of commonly used supervised learning techniques to classify the execution of a maximum effort change of direction task into predefined movement pattern as well as the influence of fuzzy executions and the impact of selected features (e.g. peak knee flexion) towards classification accuracy. The experiment utilized kinematic and kinetic data from 323 male subjects with chronic athletic groin pain. All subjects undertook a biomechanical assessment and had been divided previously into 3 different movement strategies in an earlier paper. Examined supervised learning techniques were: a decision tree, an ensemble of decision trees, a discriminant analysis model, a naive Bayes classifier, a k-nearest-neighbour model, a multi-class model for support vector machines, a stepwise forward regression model, a neural network and a correlation approach. Performance (measured by comparing the predefined and classified movement pattern) was highest for the correlation approach (82% - CI 81-83%) and support vector machine (80% - CI 79-80%). The percentage of fuzzy observations within the data was between 16 and 25%. The most informative features for classification were: hip flexion and ankle rotation as well as ankle flexion moment, thorax [flexion and frontal sway], abduction angle in [hip and pelvis] and hip rotation. Findings of this study support the assumption that multiple patterns are used to execute a movement task and demonstrate that classification models can predict movement patterns with a high accuracy (~84%).
A computational framework for estimating statistical power and planning hypothesis-driven experiments involving one-dimensional biomechanical continua.
Venue: Journal of biomechanics

Statistical power assessment is an important component of hypothesis-driven research but until relatively recently (mid-1990s) no methods were available for assessing power in experiments involving continuum data and in particular those involving one-dimensional (1D) time series. The purpose of this study was to describe how continuum-level power analyses can be used to plan hypothesis-driven biomechanics experiments involving 1D data. In particular, we demonstrate how theory- and pilot-driven 1D effect modeling can be used for sample-size calculations for both single- and multi-subject experiments. For theory-driven power analysis we use the minimum jerk hypothesis and single-subject experiments involving straight-line, planar reaching. For pilot-driven power analysis we use a previously published knee kinematics dataset. Results show that powers on the order of 0.8 can be achieved with relatively small sample sizes, five and ten for within-subject minimum jerk analysis and between-subject knee kinematics, respectively. However, the appropriate sample size depends on a priori justifications of biomechanical meaning and effect size. The main advantage of the proposed technique is that it encourages a priori justification regarding the clinical and/or scientific meaning of particular 1D effects, thereby robustly structuring subsequent experimental inquiry. In short, it shifts focus from a search for significance to a search for non-rejectable hypotheses.
The shear wave elastic modulus and the increased nuclear factor kappa B (NF-kB/p65) and cyclooxygenase-2 (COX-2) expression in the area of myofascial trigger points activated in a rat model by blunt trauma to the vastus medialis.
Venue: Journal of biomechanics

We aimed to elucidate the increased inflammatory cytokines expression such as nuclear factor kappa B (NF-kB/p65), cyclooxygenase-2 (COX-2), and voltage-gated calcium channels (VGCC) in the area of activated myofascial trigger points (MTrPs) in a rat model by blunt trauma to the vastus medialis and to evaluate the feasibility of a quantitative analysis of muscle elastic modulus using shear wave elastography (SWE). Twelve 7-week-old male SD rats were divided into normal (NM, n = 6) and model groups (MO, n = 6). In the MO group, MTrPs were activated with a blunt strike to the left vastus medialis and subsequent eccentric exercise for 8 weeks. After 4 weeks of rest, the elastic modulus in the focal site was evaluated using SWE. Electromyography (EMG) data were collected at MTrPs and muscle tissues were evaluated for expression of nuclear factor kappa B (NF-kB/p65), cyclooxygenase-2 (COX-2) protein, and voltage-gated calcium channels (VGCC). The number of the palpable taut bands; EMG frequency and amplitude; elastic modulus values; and NFkB/p65, COX-2, and VGCC expression levels were significantly higher for the left focal area in the MO group compared to those for the NM group (p's < 0.05). These findings suggest that elastic modulus measurement using ultrasound SWE may be effective in evaluating MPS. In addition, increased COX-2, NFkB/p65, and VGCC expression may expand the integrated hypothesis of trigger points.
A computational model for false lumen thrombosis in type B aortic dissection following thoracic endovascular repair.
Venue: Journal of biomechanics

Thoracic endovascular repair (TEVAR) has recently been established as the preferred treatment option for complicated type B dissection. This procedure involves covering the primary entry tear to stimulate aortic remodelling and promote false lumen thrombosis thereby restoring true lumen flow. However, complications associated with incomplete false lumen thrombosis, such as aortic dilatation and stent graft induced new entry tears, can arise after TEVAR. This study presents the application and validation of a recently developed mathematical model for patient-specific prediction of thrombus formation and growth under physiologically realistic flow conditions. The model predicts thrombosis through the evaluation of shear rates, fluid residence time and platelet distribution, based on convection-diffusion-reaction transport equations. The model was applied to 3 type B aortic dissection patients: two TEVAR cases showing complete and incomplete false lumen thrombosis respectively, and one medically treated dissection with no signs of thrombosis. Predicted thrombus growth over time was validated against follow-up CT scans, showing good agreement with in vivo data in all cases with a maximum difference between predicted and measured false lumen reduction below 8%. Our results demonstrate that TEVAR-induced thrombus formation in type B aortic dissection can be predicted based on patient-specific anatomy and physiologically realistic boundary conditions. Our model can be used to identify anatomical or stent graft related factors that are associated with incomplete false lumen thrombosis following TEVAR, which may help clinicians develop personalised treatment plans for dissection patients in the future.
The quantitative evaluation of the relationship between the forces applied to the palm and carpal tunnel pressure.
Venue: Journal of biomechanics

Carpal tunnel syndrome (CTS) is the most common entrapment neuropathy occurring in upper limbs. The etiology, however, has not been fully understood yet. Median nerve could be compressed by either increase of carpal tunnel pressure (CTP) or direct impingement when it is forced toward to carpal ligament especially in wrist flexion leading to CTS development. Thus, the increase of carpal tunnel pressure is considered an important role in CTS development. It has been identified that forces applied to the palm would affect the CTP. However, the quantitative relationship between palmar contact force and CTP is not known. The purpose of this study was to quantitatively evaluate the relationship between palmar contact force and CTP. Eight human cadaveric hands were used. The CTP was measured with a diagnostic catheter-based pressure transducer inserted into the carpal tunnel. A custom made device was used to apply forces to the palm for the desired CTP. Palmar contact forces corresponding to the determined CTP level were recorded respectively. The testing was repeated with different ranges of tension applied to the flexor digitorum superficialis tendon of the third finger. The tensions were constant at 50 g for the other flexor tendons and median nerve. The results showed that CTP increased linearly with the force applied to the palm. When CTP was 30 mmHg, mean values of the contact force to the palm was 293 g (SD: 15.2) including all tensions. These results would help to understand the effect of daily activities with hands on CTP.
Mineralized collagen fibril network spatial arrangement influences cortical bone fracture behavior.
Venue: Journal of biomechanics

Bone is a hierarchical material exhibiting different fracture mechanisms at each length scale. At the submicroscale, the bone is composed of mineralized collagen fibrils (MCF). At this scale, the fracture processes in cortical bone have not been extensively studied in the literature. In this study, the influence of MCF size and orientation on the fracture behavior of bone under both transverse and longitudinal loading was investigated using novel 3D models of MCF networks with explicit representation of extra-fibrillar matrix. The simulation results showed that separation between MCFs was the main cause of damage and failure under transverse loading whereas under longitudinal loading, the main damage and failure mechanism was MCF rupture. When the MCF network was loaded in the transverse direction the mechanical properties increased as the orientation of fibrils deviated farther from the main fibril orientation whereas the opposite trend was observed under longitudinal loading. The fracture energy was much larger in longitudinal than transverse loading. MCF diameter variation did not affect the mechanical properties under longitudinal loading but led to higher mechanical properties with increasing MCF diameter under transverse loading. The new modeling framework established in this study generate unique information on the effect of MCF network spatial arrangement on the fracture behavior of bone at the submicroscale which is not currently possible to measure via experiments. This unique information may improve the understanding of how structural alterations at the submicroscale due to disease, age-related changes, and treatments affect the fracture processes at larger length scales.
A gait retraining system using augmented-reality to modify footprint parameters: Effects on lower-limb sagittal-plane kinematics.
Venue: Journal of biomechanics

Improving lower-limb flexion/extension angles during walking is important for the treatment of numerous pathologies. Currently, these gait retraining procedures are mostly qualitative, often based on visual assessment and oral instructions. This study aimed to propose an alternative method combining motion capture and display of target footprints on the floor. The second objectives were to determine the error in footprint modifications and the effects of footprint modifications on lower-limb flexion/extension angles. An augmented-reality system made of an optoelectronic motion capture device and video projectors displaying target footprints on the floor was designed. 10 young healthy subjects performed a series of 27 trials, consisting of increased and decreased amplitudes in stride length, step width and foot progression angle. 11 standard features were used to describe and compare lower-limb flexion/extension angles among footprint modifications. Subjects became accustomed to walk on target footprints in less than 10 min, with mean (± SD) precision of 0.020 ± 0.002 m in stride length, 0.022 ± 0.006 m in step width, and 2.7 ± 0.6° in progression angle. Modifying stride length had significant effects on 3/3 hip, 2/4 knee and 4/4 ankle features. Similarly, step width and progression angle modifications affected 2/3 and 1/3 hip, 2/4 and 1/4 knee as well as 3/4 and 2/4 ankle features, respectively. In conclusion, this study introduced an augmented-reality method allowing healthy subjects to modify their footprint parameters rapidly and precisely. Walking with modified footprints changed lower-limb sagittal-plane kinematics. Further research is needed to design rehabilitation protocols for specific pathologies.
Simplified boundary conditions alter cortical-trabecular load sharing at the distal radius; A multiscale finite element analysis.
Venue: Journal of biomechanics

High-resolution peripheral quantitative computed tomography (HR-pQCT) derived micro-finite element (FE) modeling is used to evaluate mechanical behavior at the distal radius microstructure. However, these analyses typically simulate non-physiologic simplified platen-compression boundary conditions on a small section of the distal radius. Cortical and trabecular regions contribute uniquely to distal radius mechanical behavior, and various factors affect these regions distinctly. Generalized strength predictions from standardized platen-compression analyses may not adequately capture region specific responses in bone. Our goal was to compare load sharing within the cortical-trabecular compartments between the standardized platen-compression BC simulations, and physiologic BC simulations using a validated multiscale approach. Clinical- and high-resolution images were acquired from nine cadaveric forearm specimens using an HR-pQCT scanner. Multiscale FE models simulating physiologic BCs, and micro-FE only models simulating platen-compression BCs were created for each specimen. Cortical and trabecular loads (N) along the length of the distal radius micro-FE section were compared between BCs using correlations. Principal strain distributions were also compared quantitatively. Cortical and trabecular loads from the platen-compression BC simulations were strongly correlated to the physiologic BC simulations. However, a 30% difference in cortical loads distally, and a 53% difference in trabecular loads proximally was observed under platen BC simulations. Also, distribution of principal strains was clearly different. Our data indicated that platen-compression BC simulations alter cortical-trabecular load sharing. Therefore, results from these analyses should be interpreted in the appropriate mechanical context for clinical evaluations of normal and pathologic mechanical behavior at the distal radius.
Speed profiles in wheelchair court sports; comparison of two methods for measuring wheelchair mobility performance.
Venue: Journal of biomechanics

Wheelchair mobility performance is an important aspect in most wheelchair court sports, commonly measured with an indoor tracking system or wheelchair bound inertial sensors. Both methods provide key wheelchair mobility performance outcomes regarding speed. In this study, we compared speed profiles of both methods to gain insight into the level of agreement, for recommendations regarding future performance measurement. Data were obtained from 5 male highly trained wheelchair basketball players during match play. Players were equipped simultaneously with a tag on the footplate for the indoor tracking system (∼8 Hz) and inertial sensors on both wheels and frame (199.8 Hz). Being part of a larger study on 3 vs 3 player game formats, data were collected in several matches with varying field sizes, but activity profiles closely resembled regular match play. Both systems provide similar outcomes regarding distance covered and average speed. Due to differences in sampling frequency and sensor location (reference point) on the wheelchair (for speed calculation), minor differences were revealed at low speeds (<2.5 m/s). Since both systems provide complementary features, a hybrid solution as proved feasible in this study, could possibly serve as the new gold standard for mobility performance measurement in wheelchair basketball or wheelchair court sports in general.
A procedure for determining parameters of a simplified ligament model.
Venue: Journal of biomechanics

A previous mathematical model of ligament force-generation treated their behavior as a population of collagen fibres arranged in parallel. When damage was ignored in this model, an expression for ligament force in terms of the deflection, x, effective stiffness, k, mean collagen slack length, μ, and the standard deviation of slack lengths, σ, was obtained. We present a simple three-step method for determining the three model parameters (k, μ, and σ) from force-deflection data: (1) determine the equation of the line in the linear region of this curve, its slope is k and its x -intercept is -μ; (2) interpolate the force-deflection data when x is -μ to obtain F0; (3) calculate σ with the equation σ=2πF0/k. Results from this method were in good agreement to those obtained from a least-squares procedure on experimental data - all falling within 6%. Therefore, parameters obtained using the proposed method provide a systematic way of reporting ligament parameters, or for obtaining an initial guess for nonlinear least-squares.
Computational study of the role of fluid content and flow on the lumbar disc response in cyclic compression: Replication of in vitro and in vivo conditions.
Venue: Journal of biomechanics

The intervertebral disc viscoelastic response is governed primarily by its fluid content and flow. Invivo measurements demonstrate that the disc volume, fluid content, height and nucleus pressure completely recover during resting even after diurnal loading with twice longer duration (16 vs. 8 h). In view of much longer periods required for the recovery of disc height and pressure in vitro, concerns have been raised on the fluid inflow through the endplates that might be hampered by clogged blood vessels post mortem. This in silico study aimed to identify fluid-flow dependent response of discs and conditions essential to replicate in vitro and in vivo observations. An osmo-poroelastic finite element model of the human lumbar L4-L5 disc-bone unit was used. Simulating earlier in vitro experiments on bovine discs, the loading protocol started with 8 h preload at 0.06 MPa followed by 30 high/low compression loading cycles each lasting 7.5min at 0.5/0.06 MPa, respectively. Three different endplate configurations were investigated: free in- and outflow, no inflow and closed endplates with no flow. Additionally, the preload magnitude was increased from 0.06 MPa to 0.28 MPa and 0.50 MPa, or the initial nucleus hydration was reduced from 83% to 50%. For 0.06 MPa preload, the model with no inflow best matched in vitro trends. The model with free inflow increased segment height and nucleus pressure while the model with no fluid inflow resulted in a relatively small recovery in segment height and a rather constant nucleus pressure during unloading periods. Results highlight an excessive mobile fluid content as well as a restricted fluid inflow through endplates as likely causes of the discrepancies between in vivo and in vitro studies. To replicate in vivo conditions in vitro and in silico, disc hydration level should be controlled by adequate selection of preload magnitude/period and/or mobile fluid porosity.