nils-herrmann/clean_scirepeval huggingface dataset:
Flow rates in perfusion bioreactors to maximise mineralisation in bone tissue engineering in vitro.
Venue: Journal of biomechanics

In bone tissue engineering experiments, fluid-induced shear stress is able to stimulate cells to produce mineralised extracellular matrix (ECM). The application of shear stress on seeded cells can for example be achieved through bioreactors that perfuse medium through porous scaffolds. The generated mechanical environment (i.e. wall shear stress: WSS) within the scaffolds is complex due to the complexity of scaffold geometry. This complexity has so far prevented setting an optimal loading (i.e. flow rate) of the bioreactor to achieve an optimal distribution of WSS for stimulating cells to produce mineralised ECM. In this study, we demonstrate an approach combining computational fluid dynamics (CFD) and mechano-regulation theory to optimise flow rates of a perfusion bioreactor and various scaffold geometries (i.e. pore shape, porosity and pore diameter) in order to maximise shear stress induced mineralisation. The optimal flow rates, under which the highest fraction of scaffold surface area is subjected to a wall shear stress that induces mineralisation, are mainly dependent on the scaffold geometries. Nevertheless, the variation range of such optimal flow rates are within 0.5-5 mL/min (or in terms of fluid velocity: 0.166-1.66 mm/s), among different scaffolds. This approach can facilitate the determination of scaffold-dependent flow rates for bone tissue engineering experiments in vitro, avoiding performing a series of trial and error experiments.
Pressure and flow in the umbilical cord.
Venue: Journal of biomechanics

A fluid dynamic study of blood flow within the umbilical vessels of the human maternal-fetal circulatory system is considered. It is found that the umbilical coiling index (UCI) is unable to distinguish between cords of significantly varying pressure and flow characteristics, which are typically determined by the vessel curvature, torsion and length. Larger scale geometric non-uniformities superposed over the inherent coiling, including cords exhibiting width and/or local UCI variations as well as loose true knots, typically produce a small effect on the total pressure drop. Crucially, this implies that a helical geometry of mean coiling may be used to determine the steady vessel pressure drop through a more complex cord. The presence of vessel constriction, however, drastically increases the steady pressure drop and alters the flow profile. For pulsatile-flow within the arteries, the steady pressure approximates the time-averaged value with high accuracy over a wide range of cords. Furthermore, the relative peak systolic pressure measured over the period is virtually constant and approximately 25% below the equivalent straight-pipe value for a large range of non-straight vessels. Interestingly, this suggests that the presence of vessel helicity dampens extreme pressures within the arterial cycle and may provide another possible evolutionary benefit to the coiled structure of the cord.
A Monte Carlo analysis of muscle force estimation sensitivity to muscle-tendon properties using a Hill-based muscle model.
Venue: Journal of biomechanics

Surface electromyography driven models are desirable for estimating subject-specific muscle forces. However, these models include parameters that come from an array of sources, thus creating uncertainty in the model-estimated force. In this study, we used Monte-Carlo simulations to evaluate the sensitivity of Hill-based model muscle forces to changes in 11 parameters in the muscle-tendon unit morphological properties and in the model force-length and force-velocity relationships. We decomposed the force variability and ranked the sensitivity of the model to the underlying parameters using the Variogram Analysis of Response Surfaces. For the analyzed running experiments and the adopted Hill model structure, our results show that the parameters are separable into four groups, where the parameters in each group have a synergistic contribution to the model global sensitivity. The first group consists of the maximum isometric force and the pennation angle. The second group contains the optimal fiber length, the tendon slack length, the tendon reference strain and the tendon shape factor. The third group contains the width and shape at the extremities of the active contractile element, along with the maximum contraction velocity and the curvature constant in the force-velocity curve. The fourth group consisted only of the force enhancement during eccentric contraction. The first two groups revealed the largest influence on the output force sensitivity. As many input parameters are difficult to measure and impact estimated forces, we propose that model estimates be presented with confidence intervals as well as inter-parameter relationships, to encourage users to explicitly consider the model uncertainty.
Prediction of stenting related adverse events through patient-specific finite element modelling.
Venue: Journal of biomechanics

Right ventricular outflow tract (RVOT) calcific obstruction is frequent after homograft conduit implantation to treat congenital heart disease. Stenting and percutaneous pulmonary valve implantation (PPVI) can relieve the obstruction and prolong the conduit lifespan, but require accurate pre-procedural evaluation to minimize the risk of coronary artery (CA) compression, stent fracture, conduit injury or arterial distortion. Herein, we test patient-specific finite element (FE) modeling as a tool to assess stenting feasibility and investigate clinically relevant risks associated to the percutaneous intervention. Three patients undergoing attempted PPVI due to calcific RVOT conduit failure were enrolled; the calcific RVOT, the aortic root and the proximal CA were segmented on CT scans for each patient. We numerically reproduced RVOT balloon angioplasty to test procedure feasibility and the subsequent RVOT pre-stenting expanding the stent through a balloon-in-balloon delivery system. Our FE framework predicted the occurrence of CA compression in the patient excluded from the real procedure. In the two patients undergoing RVOT stenting, numerical results were consistent with intraprocedural in-vivo fluoroscopic evidences. Furthermore, it quantified the stresses on the stent and on the relevant native structures, highlighting their marked dependence on the extent, shape and location of the calcific deposits. Stent deployment induced displacement and mechanical loading of the calcific deposits, also impacting on the adjacent anatomical structures. This novel workflow has the potential to tackle the analysis of complex RVOT clinical scenarios, pinpointing the procedure impact on the dysfunctional anatomy and elucidating potential periprocedural complications.
Linear and nonlinear analyses of femoral fractures: Computational/experimental study.
Venue: Journal of biomechanics

This paper describes two new methods for computational fracture analysis of human femur using Quantitative Computed Tomography (QCT) voxel-based finite element (FE) simulation. The paper also reports comprehensive mechanical testing for validation of the methods and evaluation of the required material properties. The analyses and tests were carried out on 15 human femurs under 11 different stance-type loading orientations. Several classical forms of subcapital, transcervical, basicervical, and intertrochanteric fractures plus a specific type of subtrochanteric fracture were created and analyzed. A new procedure was developed for prediction of the strengths and the fracture initiation patterns using a FE-based linear scheme. The predicted and observed fracture patterns were in correspondence, and the FE predictions of the fracture loads were in very good agreement with the experimental results. Moreover, the crack initiation and growth behaviors of two subtrochanteric fractures were successfully simulated through a novel implementation of the cohesive zone model (CZM) within a nonlinear FE analysis scheme. The CZM parameters were obtained through a series of experimental tests on different types of specimens and determination of a variety of material properties for different anatomic regions and orientations. The presented results indicated that the locations and patterns of crack initiation, the sequences of crack growth on different paths, and the compatibility of growth increments agreed very well with the observed specifications. Also, very good agreements were achieved between the measured and simulated fracture loads.
Validation of gait event detection by centre of pressure during target stepping in healthy and paretic gait.
Venue: Journal of biomechanics

Target-stepping paradigms are increasingly used to assess and train gait adaptability. Accurate gait-event detection (GED) is key to locating targets relative to the ongoing step cycle as well as measuring foot-placement error. In the current literature GED is either based on kinematics or centre of pressure (CoP), and both have been previously validated with young healthy individuals. However, CoP based GED has not been validated for stroke survivors who demonstrate altered CoP pattern.



Young healthy adults and individuals affected by stroke stepped to targets on a treadmill, while gait events were measured using three detection methods; verticies of CoP cyclograms, and two kinematic criteria, (1) vertical velocity and position and of the heel marker, (2) anterior velocity and position of the heel and toe marker, were used. The percentage of unmatched gait events was used to determine the success of the GED method. The difference between CoP and kinematic GED methods were tested with two one sample (two-tailed) t-tests against a reference value of zero. Differences between group and paretic and non-paretic leg were tested with a repeated measures ANOVA.



The kinematic method based on vertical velocity only detected about 80% of foot contact events on the paretic side in stroke survivors while the method on anterior velocity was more successful in both young healthy adults as stroke survivors (3% young healthy and 7% stroke survivors unmatched). Both kinematic methods detected gait events significantly earlier than CoP GED (p < 0.001) except for foot contact in stroke survivors based on the vertical velocity.



CoP GED may be more appropriate for gait analyses of SS than kinematic methods; even when walking and varying steps.
Influence of proximal trunk borne load on lower limb countermovement joint dynamics.
Venue: Journal of biomechanics

Vertical jumping involves coordinating the temporal sequencing of angular motion, moment, and power across multiple joints. Studying the biomechanical coordination strategies that differentiates loaded from unloaded vertical jumping may better inform training prescription for athletes needing to jump with load. Common multivariate methods (e.g. Principal Components Analysis) cannot quantify coordination in a dataset with more than two modes. This study aimed to identify coordinative factors across four modes of variation using Parallel Factor (Parafac2) analysis, which may differentiate unloaded (body weight [BW]) from loaded (BW + 20% BW) countermovement jump (CMJ). Thirty-one participants performed unloaded and loaded CMJ. Three-dimensional motion capture with force plate analysis was performed. Inverse dynamics was used to quantify sagittal plane joint angle, velocity, moment, and joint power across the ankle, knee, and hip. The four-mode data were as follows: Mode A was jump cycle (101 cycle points), mode B was participant (31 participants by two load), mode C was joint (two sides by three joints), and mode D was variable (angle, velocity, moment, power). Three factors were extracted, which explained 95.1% of the data's variance. Only factors one (P = 0.001) and three (P < 0.001) significantly differentiated loaded from unloaded jumping. The body augmented hip-dominant at the start, and both hip and ankle dominant behaviors at the end of the ascending phase of the CMJ, but kept knee-dominant behavior invariant, when jumping with a 20% BW load. By studying the variation across all data modes, Parafac2 provides a holistic method of studying jumping coordination.
Gait variability and symmetry remain consistent during high-intensity 10,000 m treadmill running.
Venue: Journal of biomechanics

The aim of this study was to analyze changes in gait variability and symmetry in distance runners. Fourteen competitive athletes ran on an instrumented treadmill for 10,000 m at speeds equivalent to 103% of their season's best time. Spatiotemporal and ground reaction force data were recorded at 1500, 3000, 5000, 7500 and 9500 m. Gait variability and inter-leg symmetry were measured using median absolute deviation (MAD) and the symmetry angle, respectively. There were no overall changes during the running bout for absolute values, symmetry angles or variability, and there were only moderate changes in variability between successive testing distances for three variables. Even with these few changes, variability was low (<4%) at all distances for all variables measured and, on average, the athletes were symmetrical for five of the seven gait variables measured. There were greater mean asymmetry values for flight time (1.1-1.4%) and for impact force (2.0-2.9%), which might have occurred because of muscle latency as the lower limb responded passively to impact during initial contact. Although most athletes were asymmetrical (>1.2%) for at least one variable, no one was asymmetrical for more than four of the seven variables measured. Being asymmetrical in a few variables is therefore not abnormal and not indicative of asymmetrical gait and given many practitioners analyze symmetry (and variability) on an individual, case-study basis, caution should be taken when assessing the need for corrective interventions.
Landmark registering waveform data improves the ability to predict performance measures.
Venue: Journal of biomechanics

The purpose of this study was to investigate the benefit of landmark registration when applied to waveform data. We compared the ability of data reduced from time-normalised and landmark registered vertical ground reaction force (vGRF) waveforms captured during maximal countermovement jumps (CMJ) of 53 active male subjects to predict jump height. vGRF waveforms were landmark registered using different landmarks resulting in four registration conditions: (i) end of the eccentric phase, (ii) adding maximum centre of mass (CoM) power, (iii) adding minimum CoM power, (iv) adding minimum vGRF. In addition to the four registration conditions, the non-registered vGRF and concentric phase only were time-normalised and used in subsequent analysis. Analysis of characterising phases was performed to reduce the vGRF data to features that captured the behaviour of each waveform. These features were extracted from each condition's vGRF waveform, time-domain (time taken to complete the movement), and warping functions (generated from landmark registration). The identified features were used as predictor features to fit a step-wise multilinear regression to jump height. Features generated from the best performing registration condition were able to predict jump height to a similar extent as the concentric phase (86-87%), while all registration conditions could explain jump height to a greater extent than time-normalisation alone (65%). This suggests waveform variability was reduced as phases of the CMJ were aligned. However, findings suggest that over-registration can occur when applying landmark registration. Overall, landmark registration can improve prediction power to performance measures as waveform data can be reduced to more appropriate performance related features.
The effect of charge density on the velocity and attenuation of ultrasound waves in human cancellous bone.
Venue: Journal of biomechanics

Cancellous bone is a highly porous material, and two types of waves, fast and slow, are observed when ultrasound is used for detecting bone diseases. There are several possible stimuli for bone remodelling processes, including bone fluid flow, streaming potential, and piezoelectricity. Poroelasticity has been widely used for elucidating the bone fluid flow phenomenon, but the combination of poroelasticity with charge density has not been introduced. Theoretically, general poroelasticity with a varying charge density is employed for determining the relationship between wave velocity and attenuation with charge density. Fast wave velocity and attenuation are affected by porosity as well as charge density; however, for a slow wave, both slow wave velocity and attenuation are not as sensitive to the effect of charge density as they are for a fast wave. Thus, employing human femoral data, we conclude that charged ions gather on trabecular struts, and the fast wave, which moves along the trabecular struts, is significantly affected by charge density.
Metal ion release barrier function and biotribological evaluation of a zirconium nitride multilayer coated knee implant under highly demanding activities wear simulation.
Venue: Journal of biomechanics

Total knee arthroplasty is a well established treatment for degenerative joint disease, which is also performed as a treatment in younger and middle-aged patients who have a significant physical activity and high life expectancy. However, complications may occur due to biological responses to wear particles, as well as local and systemic hypersensitivity reactions triggered by metal ions and particles such as cobalt, chromium and molybdenum. The purpose of the study was to perform a highly demanding activities (HDA) knee wear simulation in order to compare the wear characteristics and metal ion release barrier function of a zirconium nitride (ZrN) coated knee implant, designed for patients with suspected metal ion hypersensitivity, against an uncoated knee implant made out of CoCrMo. The load profiles were applied for 5 million HDA cycles, which represent 15-30 years of in vivo service depending on the activity level of the patient. Results showed a significant wear rate reduction for the coated group (1.01 ± 0.29 mg/million cycles) in comparison with the uncoated group (2.89 ± 1.04 mg/million cycles). The zirconium nitride coating showed no sign of scratches nor delamination during the wear simulation, whereas the uncoated femurs showed characteristic wear scratches in the articulation areas. Furthermore, the metal ion release from the coated implants was reduced up to three orders of magnitude in comparison with the uncoated implants. These results demonstrate the efficiency of zirconium nitride coated knee implants to reduce wear as well as to substantially reduce metal ion release in the knee joint.
Fibroblast-seeded collagen gels in response to dynamic equibiaxial mechanical stimuli: A biomechanical study.
Venue: Journal of biomechanics

The remodeling of fibroblast-seeded collagen gels in response to dynamic mechanical stimuli was investigated by using a newly developed biaxial culture system capable of cyclically stretching planar soft tissues. Fibroblast-seeded collagen gels were subjected to three distinct dynamic mechanical conditions for six days: Cyclic Equibiaxial Stretching at two constant strain magnitudes (CES-7% and CES-20%), and Cyclic Equibiaxial Stretching with incrementally Increasing stain magnitude (ICES, 7% → 15% → 20% each for two days). The frequency of cyclic stretching was set at 1 Hz. At the end of culture, mechanical properties of the gels were examined by biaxial mechanical testing and checked again upon the removal of seeded cells. Collagen microstructure within the gels was illustrated by multiphoton microscopy. The mRNA levels of collagen type I and type III and fibronectin in the cells were examined by reverse transcription PCR. The protein expression of α-smooth muscle actin was detected by immunohistochemistry. We found that the gels cultured under cyclic stretching were stiffer than those cultured under static stretching. Particularly, the stiffness appeared to be significantly enhanced when the ICES was employed. The enhancement of mechanical properties by cyclic stretching appeared to persist upon cell removal, suggesting an irreversible remodeling of extracellular matrix. Second harmonic generation images showed that collagen fibers became thicker and more compact in the gels cultured under cyclic stretching, which may explain the mechanical findings. The mRNA expression of collagen type I in the cells of the ICES was significantly greater than that of the other groups except for the CES-20%. This study suggests that when cyclic stretching is to be used in engineering soft tissues, incrementally increasing strain magnitude may prove useful in the development of the tissue.
Joint contact stresses calculated for acetabular dysplasia patients using discrete element analysis are significantly influenced by the applied gait pattern.
Venue: Journal of biomechanics

Gait modifications in acetabular dysplasia patients may influence cartilage contact stress patterns within the hip joint, with serious implications for clinical outcomes and the risk of developing osteoarthritis. The objective of this study was to understand how the gait pattern used to load computational models of dysplastic hips influences computed joint mechanics. Three-dimensional pre- and post-operative hip models of thirty patients previously treated for hip dysplasia with periacetabular osteotomy (PAO) were developed for performing discrete element analysis (DEA). Using DEA, contact stress patterns were calculated for each pre- and post-operative hip model when loaded with an instrumented total hip, a dysplastic, a matched control, and a normal gait pattern. DEA models loaded with the dysplastic and matched control gait patterns had significantly higher (p = 0.012 and p < 0.001) average pre-operative maximum contact stress than models loaded with the normal gait. Models loaded with the dysplastic and matched control gait patterns had nearly significantly higher (p = 0.051) and significantly higher (p = 0.008) average pre-operative contact stress, respectively, than models loaded with the instrumented hip gait. Following PAO, the average maximum contact stress for DEA models loaded with the dysplastic and matched control patterns decreased, which was significantly different (p < 0.001) from observed increases in maximum contact stress calculated when utilizing the instrumented hip and normal gait patterns. The correlation between change in DEA-computed maximum contact stress and the change in radiographic measurements of lateral center-edge angle were greatest (R2 = 0.330) when utilizing the dysplastic gait pattern. These results indicate that utilizing a dysplastic gait pattern to load DEA models may be a crucial element to capturing contact stress patterns most representative of this patient population.