Search

Cadaveric and computer simulations suggest lateral humeral offset (LHO) and humeral retroversion (HR) are associated with strength and range of motion (ROM) after reverse total shoulder arthroplasty (rTSA), but in vivo data is lacking. This study aimed to evaluate the effects of implant parameters (i.e. LHO and HR) on strength and ROM. LHO and HR were measured using pre -operative and post -operative computed tomography (CT) scans. Postoperative strength was measured across three planes of motion using a Biodex isokinetic dynamometer. Postoperative active and passive ROM during forward elevation, external rotation (ER), and internal rotation (IR) were assessed using a goniometer or spinal level. 30 rTSA patients (14 M, 16F, age: 71.8 ± 6.7yrs) participated with an average postoperative follow-up of 2.4 ± 1.1 years. Regarding strength, higher post -op LHO values were predictive of greater postoperative strength across all movements. However, lateralization of the implant beyond pre -op values (i.e. post -op LHO > pre -op LHO) was associated with poorer strength performance across all ranges of motion. Similar to strength outcomes, greater deviations from pre -op LHO was predictive of poorer IR ROM. Lastly, patients with minimal deviations in HR (post -op HR within 10° of pre -op HR) and minimal deviations in LHO (post -op LHO ≤ pre -op LHO) displayed the greatest postoperative ER ROM. Anatomic restoration of LHO combined with anatomic restoration of HR may be ideal for maximizing strength and ROM following rTSA. Overlateralization beyond anatomic may have negative consequences. Optimal implant lateralization and version may need to be individualized based on preoperative values. [ABSTRACT FROM AUTHOR]

Calcific aortic valve disease (CAVD) is a serious disease affecting the aging population. A complex interaction between biochemicals, cells, and mechanical cues affects CAVD initiation and progression. In this study, motivated by the progression of calcification in regions of high strain, we developed a finite element method (FEM) based spatial calcification progression model. Several cardiac cycles of transient structural FEM simulations were simulated. After each simulation cycle, calcium deposition was placed in regions of high circumferential strain. Our results show the radial expansion of calcification as spokes starting from the attachment region, agreeing very well with the reported clinical data.

The contribution of the Plantar Fascia (PF) and Long Plantar Ligament (LPL), two ligaments extending from the hindfoot to the forefoot, to arch stability has been studied in the past using in vivo, in vitro, and in silico methodologies. In silico studies were based on one single model obtained from one single subject and did not account for the known inter-subject morphological and biomechanical variations. In the present study, we developed computational dynamic models of nine different legs obtained from nine different individuals to evaluate the role of the LPL and PF in arch support, accounting for biological differences between subjects. These models were validated by comparing the simulation results against experimental results from the corresponding cadaver legs. After validation, we simulated body weight conditions for each model by applying a vertical load to the tibia, starting from zero and increasing linearly to 720 N. Kinematic and dynamic parameters, including the variation of the medial arch angle and of the navicular height, as well as the passive forces developed by the LPL and PF, were used to evaluate the contribution of these ligaments to arch support under body weight. The results indicate that a total collapse of the medial longitudinal arch occurred only when both the LPL and PF were absent, but a stable arch was maintained when either one of these two ligament structures were present. The results varied significantly among the specific models, highlighting the importance of using multiple models to account for inter-subject morphological differences., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Musculoskeletal models of the shoulder are needed to understand the mechanics of overhead motions. Existing models implementing the shoulder rhythm are generic and might not accurately represent an individual's scapular kinematics. We introduce a method to personalize the shoulder rhythm of a computational model of the upper body that defines the orientations of the clavicle and scapula based on glenohumeral joint angles. During five static calibration poses, we palpate and measure the orientation of the scapula. We explore the importance of representing shoulder elevation by introducing clavicle elevation as a degree of freedom that is independent of the glenohumeral angles. For ten subjects, we record the five calibration poses, ten additional static poses, and dynamic arm raises covering the participants' full range of motion in each body plane using optical motion capture. We examine the data using a dynamically-constrained inverse kinematics analysis. Shoulder rhythm personalization, independent clavicle elevation, and both in combination reduce the average upper body marker tracking error compared to the generic model in the static poses (26 mm to 17-20 mm) and in the dynamic trials (22 mm to 14-17 mm). Only personalization reduces the average scapula marker error (51 mm to 36-38 mm) and scapula axis-angle error (15° to 10°) compared with the palpated ground truth measurements in the static poses, and in the dynamic trials at instances that best match the static poses (53 mm to 37-40 mm, 15° to 9°). Our results show that personalizing upper body models improves kinematic tracking. We provide our experimental data, model, and methods to allow researchers to reproduce and build upon our results., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

This study evaluates the efficacy of annuloplasty repair as a standalone procedure for treating bileaflet mitral valve prolapse with mitral regurgitation (MR). Various flexible ring bands for MR of different severities were compared to assess their biomechanical impact and treatment outcomes. Computational beating heart models, based on the Living Heart Human Model, were utilized to simulate annuloplasty repairs. Repairs using bands of varying lengths were modeled on moderate and severe MR cases, considering bileaflet mitral valve prolapse. Key parameters, including regurgitant orifice area (ROA), prolapse severity, coaptation length, leaflet position, and deformation, were computed to compare conditions before and after implantation. Annuloplasty repairs effectively reduced the ROA in both moderate and severe MR cases, achieving complete sealing in selective instances. Additionally, annuloplasty repair corrected bileaflet prolapse, with prolapse severity decreasing as the annular size increased. Successful coaptation was indicated by the expansion of each leaflet's contact area distribution and percentage in contact with the opposing leaflet. The risk of systolic anterior motion, that may obstruct the left ventricular outflow tract, was minimized, as the anterior leaflet was directed towards the posterior position. In conclusion, annuloplasty repair alone can effectively treat MR when an appropriate band length is selected. It facilitates a significant reduction in ROA, correction of bileaflet prolapse, and improvement in leaflet coaptation. These findings have important clinical implications, potentially offering a less complex surgical treatment avenue and reducing complications in the management of MR., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Although the hemodynamics of cerebral aneurysms have been extensively studied using patient-specific computational fluid dynamics techniques, no specific hemodynamic factors characteristic of cerebral aneurysm development have yet been identified. We believe that one problem with previous hemodynamic studies of cerebral aneurysms has been the manner in which control groups were created for comparison with experimental groups. The purpose of this study was to determine hemodynamic factors that correlated with the development of cerebral aneurysms. The control group was established in a manner that differed from those of previous works. This allowed us to demonstrate the effectiveness of our method. We artificially removed aneurysms in the middle cerebral artery bifurcations of nine patients and reconstructed the vessel geometries before the aneurysms had occurred. Pulsatile blood flow simulations were performed using the vessel geometries ipsilateral and contralateral to the sites of aneurysm removal, and hemodynamic metrics were calculated. Use of the ipsilateral and contralateral sides as the experimental and control sites, respectively, allowed us to evaluate statistically the hemodynamic metrics between the two corresponding sites/groups. The results showed that only the normalized transverse wall shear stress (NtransWSS) was significantly higher at the MCA bifurcation ipsilateral to the site of aneurysm removal than at the contralateral bifurcation (p = 0.01). There were no significant differences in the other hemodynamic metrics between the bilateral bifurcations. Our findings imply that multi-directional disturbed wall shear stress, which is detected by the NtransWSS metric, may be one hemodynamic risk factor for the development of cerebral aneurysms., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

The skin is a multilayered organ with microstructural and antomical heterogeneities that contribute to its unique mechanophysiology. Between the epidermis layer at the top and the dermis layer below, the basal keratinocytes form an interface with sinusoidal-like geometry termed rete ridges. In previous computational work we showed that the rete ridges contribute to lower delamination risk by increasing surface area and reducing the stress jump across the interface. Experimentally, we and others have shown that upon repeated tissue expansion and growth, physiological rete ridge frequency is preserved. Here we implement a 2D multilayered skin model where each layer is able to grow in response to applied loading toward recovering the layer-specific homeostatic stretch. Our simulations support the hypothesis that mechanics of growing tissue can explain secondary buckling and new rete ridge formation in tissue expansion. The process is robust with respect to parameters such as homeostatic stretch, layer thicknesses, and shear moduli of the different layers. Thicker epidermis suppresses higher frequency features, and so does a stiffer epidermis with respect to the basal layer. Interestingly, new rete ridge valleys are formed at locations that were originally peaks of the sine wave, whereas original valleys remain valleys. This pattern might have a connection to the localization of stem cell and transient amplifying cells in the epidermis. This study does not discard the role of cell-cell signaling dynamics, but rather emphasizes the possibility of achieving robust geometric patterns with simple rules of growing tissue, even in the absence of complex regulatory networks., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Osteoporosis is characterized by loss of bone mineral density and increased fracture risk. Reduction of hip fracture incidence is of major clinical importance. Hip protectors aim to attenuate the impact force transmitted to the femur upon falling, however different conclusions on their efficacy have been reported; some authors suggest this may be due to differences in compliance. The aim of this study was to apply an In Silico trial methodology to predict the effectiveness of hip protectors and its dependence on compliance. A cohort of 1044 virtual patients (Finite Element models of proximal femur) were generated. A Markov chain process was implemented to predict fracture incidence with and without hip protectors, by simulating different levels of compliance. At each simulated follow-up year, a Poisson distribution was randomly sampled to determine the number of falls sustained by each patient. Impact direction and force were stochastically sampled from a range of possible scenarios. The effect of wearing a hip protector was simulated by applying attenuation coefficients to the impact force (12.9 %, 19 % and 33.8 %, as reported for available devices). A patient was considered fractured when impact force exceeded the femur strength. Without hip protector, virtual patients experienced 66 ± 5 fractures in 10 years. Wearing the three devices, fracture incidence was reduced to 43 ± 4, 35 ± 4 and 17 ± 2 respectively, at full compliance. As expected, effectiveness was dependent on compliance. This In Silico trial technology can be applied in the future to test multiple interventions, optimise intervention strategies, improve clinical trial design and drug development., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Swimmers primarily increase their forward velocity through lower limb motion in breaststroke, making the breaststroke kick crucial for optimizing race times. Recent studies have highlighted the generation of vortices around the swimmer's entire body to propel forward during swimming. However, the investigation of vortex generation during breaststroke kicks remains unexplored. This study aimed to reveal the propulsive and braking mechanisms of breaststroke kicks by simulating vortex generation using computational fluid dynamics (CFD). Kinematic data during the breaststroke kick and a three-dimensional digital model were collected to conduct CFD for a male breaststroke swimmer. Vortex generation was determined during one breaststroke kick from the CFD results. Vortices, which potentially induce a decrease in forward velocity, were generated by the swimmer's lower legs and feet during the recovery phase. The swimmer generated vortices on the dorsal side of the feet and the posterior and lateral sides of the lower legs to increase the forward velocity during the out-sweep phase. The swimmer generated vortices on the lateral sides of the thighs and lower legs and the dorsal and lateral sides of the feet during the in-sweep phase to maintain forward velocity. Moreover, vortices generated from the out-sweep to the in-sweep merged and were shed backward relative to the swimming direction after the in-sweep phase. This study is the first to reveal the propulsive and braking mechanisms of breaststroke kicks by analyzing the vortex generation., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Splenic artery embolization (SAE) has become a favored alternative to splenectomy, offering a less invasive intervention for injured spleens while preserving spleen function. However, our understanding of the role that hemodynamics plays during embolization remains limited. In this study, we utilized patient-specific computational fluid dynamics (CFD) simulations to study distal and proximal embolization strategies commonly used in SAE. Detailed 3D computer models were constructed considering the descending aorta, various major visceral arteries, and the iliac arteries. Subsequently, the blood flow and pressure associated with different coil placement locations in proximal embolization were studied considering the collateral vessels. Coil induced variations in pressure fields were quantified and compared to baseline. The coil induced flow stagnation was also quantified with particle residence time. Distal embolization was modeled with Lagrangian particle tracking and the effect of particle size, release location, and timing on embolization outcome was studied. Our findings highlight the crucial role of collateral vessels in maintaining blood supply to the spleen following proximal embolization. It was demonstrated that coil location can affect distal pressure and that strategic coil placement guided by patient-specific CFD simulations can further reduce this pressure as desired. Additionally, the results point to the critical roles that particle size, release timing, and location play in distal embolization. Our study provides an early attempt to use patient-specific computer modeling for optimizing embolization strategies and ultimately improving patient outcomes during SAE procedures., Competing Interests: Declaration of competing interest The authors hereby declare that they have no conflicts of interest., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Agarose is commonly used for 3D cell culture and to mimic the stiffness of the pericellular matrix of articular chondrocytes. Although it is known that both temperature and mechanical stimulation affect the metabolism of chondrocytes, little is known about the thermal properties of agarose hydrogels. Thermal properties of agarose are needed to analyze potential heat production by chondrocytes induced by various experimental stimuli (carbon source, cyclical compression, etc). Utilizing ASTM C177, a custom-built thermal conductivity measuring device was constructed and used to calculate the thermal conductivity of 4.5 % low gelling temperature agarose hydrogels. Additionally, Differential Scanning Calorimetry was used to calculate the specific heat capacity of the agarose hydrogels. Testing of chondrocyte-embedded agarose hydrogels commonly occurs in Phosphate-Buffered Saline (PBS), and thermal analysis requires the free convection coefficient of PBS. This was calculated using a 2D heat conduction simulation within MATLAB in tandem with experimental data collected for known boundary and initial conditions. The specific heat capacity and thermal conductivity of 4.5 % agarose hydrogels was calculated to be 2.85 J/g°C and 0.121 W/mK, respectively. The free convection coefficient of PBS was calculated to be 1000.1 W/m 2 K. The values of specific heat capacity and thermal conductivity for agarose are similar to the reported values for articular cartilage, which are 3.20 J/g°C and 0.21 W/mK (Moghadam, et al. 2014). These data show that cyclical loading of hydrogel samples with these thermal properties will result in negligible temperature increases. This suggests that in addition to 4.5 % agarose hydrogels mimicking the physiological stiffness of the cartilage PCM, they can also mimic the thermal properties of articular cartilage for in vitro studies., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Erik Myers has no conflicts of interest to disclose. Molly Piazza has no conflicts of interest to disclose. Mark Owkes has no conflicts of interest to disclose. Ron June has stock in Beartooth Biotech and OpenBioWorks. Neither company was involved in the design, execution, or analysis of this study., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

An athlete's posture has a significant impact on aerodynamic drag. Although aerodynamic drag in different sports has been studied extensively, most studies have analysed only a limited number of positions, and no generalized methods for optimization are available. In this work, we present a methodology to perform athlete posture optimization with respect to aerodynamic drag reduction. The method combines the virtual skeleton methodology to adjust the athlete's posture, CFD simulations to evaluate the drag for a given posture, and efficient global optimization to find the optimum position. We demonstrate the method by optimizing the time trial position for a cyclist. The cyclist position was parameterized with 6 design parameters, and the optimization required 41 CFD simulations to converge. The optimal posture yielded a reduction in drag of 17 % compared to the initial posture (disregarding bicycle drag). The method has potential to make posture optimization more accessible across a wide range of sports, and lead to insight into the aerodynamic influence of posture in general., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

We introduce a new computational framework that utilises Pulse Wave Velocity (PWV) extracted directly from 4D flow MRI (4DMRI) to inform patient-specific compliant computational fluid dynamics (CFD) simulations of a Type-B aortic dissection (TBAD), post-thoracic endovascular aortic repair (TEVAR). The thoracic aortic geometry, a 3D inlet velocity profile (IVP) and dynamic outlet boundary conditions are derived from 4DMRI and brachial pressure patient data. A moving boundary method (MBM) is applied to simulate aortic wall displacement. The aortic wall stiffness is estimated through two methods: one relying on area-based distensibility and the other utilising regional pulse wave velocity (RPWV) distensibility, further fine-tuned to align with in vivo values. Predicted pressures and outlet flow rates were within 2.3 % of target values. RPWV-based simulations were more accurate in replicating in vivo hemodynamics than the area-based ones. RPWVs were closely predicted in most regions, except the endograft. Systolic flow reversal ratios (SFRR) were accurately captured, while differences above 60 % in in-plane rotational flow (IRF) between the simulations were observed. Significant disparities in predicted wall shear stress (WSS)-based indices were observed between the two approaches, especially the endothelial cell activation potential (ECAP). At the isthmus, the RPWV-driven simulation indicated a mean ECAP>1.4 Pa -1 (critical threshold), indicating areas potentially prone to thrombosis, not captured by the area-based simulation. RPWV-driven simulation results agree well with 4DMRI measurements, validating the proposed pipeline and facilitating a comprehensive assessment of surgical decision-making scenarios and potential complications, such as thrombosis and aortic growth., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Venous thromboembolism (VTE) has been occurring frequently in human society. There is an urgent need to study the influence of several factors on thrombolytic therapy, such as the effects of vascular pressure levels (VPL) and the drug injection time (DIT). Considering blood as a non-Newtonian fluid, valve as a hyperelastic material, and thrombus as a porous medium, a new numerical simulation model of biofluid mechanics incorporating fluid-solid coupling phenomena and biochemical substance reactions is established based on the N-S equations and the convection-diffusion reaction equations. Then, a unique in vitro experimental platform is established to verify the correctness of the constructed mathematical model. The results showed that vascular compression resulted in significant differences in blood flow status localized within the vessel. Vascular compression causes the blood boosting index to fluctuate and the valve displacement values are 135% and 158% greater than the lower VPL, respectively. At the same time, vascular compression weakened vortex intensity, accelerated material transport and response, and improved the treatment. Compared with low VPL, the therapeutic efficacy increased by 7% and 15%, respectively. In addition, when the dose of the drug is high, different injection times can increase the therapeutic effect to different degrees, with a maximum difference of 12%. Our in vitro experiments are similar to the results obtained by numerical simulation, which can verify the reliability of numerical simulation. The computational model proposed and the experimental platform designed in this study have the potential to assist in clinical medication prediction in different venous thromboembolism patients., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Premature birth poses a challenge to public health, with one in ten babies being born prematurely worldwide. The pathological distension of the uterus can create tension in the uterine wall, triggering contractions that may lead to birth, including premature birth. While there has been an increase in the use of computational models to study pregnancy in recent years, ethical challenges have limited research on the mechanical properties of the uterus during gestation. This study proposes a biomechanical model based on a stretch-driven growth mechanism to describe uterine evolution during the second half of the gestational period. The constitutive model employed is anisotropic, reflecting the presence of fibers in uterine tissue, and it is also considered incompressible. The geometric model representing the uterine body was derived from truncated ellipsoids, subject to intrauterine pressure as loading. Simulation results indicate that the proposed model is effective in reproducing growth patterns documented in the literature, such as simultaneous increases in intrauterine volume and uterine tissue volume, accompanied by a reduction in uterine wall thickness within limits reported in experimental data., Competing Interests: Declaration of competing interest The authors certify that they have NO affiliations with or involvement in any organization or entity with any financial interest (such as honoraria; educational grants; participation in speakers’ bureaus; membership, employment, consultancies, stock ownership, or other equity interest; and expert testimony or patent-licensing arrangements), or non-financial interest (such as personal or professional relationships, affiliations, knowledge or beliefs) in the subject matter or materials discussed in this manuscript., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Recent computational and experimental studies of intracranial aneurysms have revealed potential mechanisms of aneurysm bruits and murmurs, driven by flow instabilities rather than by stable pulsatile flow. Some of these studies have been conducted under the assumption of constant flow rate (steady flow); however the validity of this assumption has not been evaluated for high-frequency flow instability, or vibrations from fluid-structure interaction (FSI) simulations. We evaluated the time-averaged wall shear stress, flow instability and vibration amplitude of steady flow simulations, performed at both cycle-averaged and peak-systolic flow rates, and compared these to recent pulsatile FSI simulations. Wall shear stress fields of pulsatile flow (time-averaged and peak values) were well-approximated by the respective steady-flow FSI simulations, and the spatial distribution and frequency content of flow instability and vibrations were reasonably approximated by the steady flow simulations at peak-systolic flow rates. However, the level of flow instability and vibration was generally over-predicted by the steady flow simulations at peak-systolic flow rates as flow remained unstable for longer than in the pulsatile simulation, while no flow instability was detected for steady flow at cycle-averaged flow rates. Additionally, the amplitude of flow instability and vibration fluctuated considerably in the steady flow simulations, while the pulsatile simulations exhibited consistent vibration amplitudes (less than 10 % variation at peak systole between cycles). Finally, steady flow simulations at peak-systolic conditions required 2-3x more compute time than the pulsatile simulations for the same time duration. Therefore, we recommend using pulsatile flow simulations when investigating vibrations and flow instabilities., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Unruptured intracranial aneurysms are common in the general population, and many uncertainties remain when predicting rupture risks and treatment outcomes. One of the cutting-edge tools used to investigate this condition is computational fluid dynamics (CFD). However, CFD is not yet mature enough to guide the clinical management of this disease. In addition, recent studies have reported significant flow instabilities when refined numerical methods are used. Questions remain as to how to properly simulate and evaluate this flow, and whether these instabilities are really turbulence. The purpose of the present study is to evaluate the impact of the simulation setup on the results and investigate the occurrence of turbulence in a cerebral artery with an aneurysm. For this purpose, direct numerical simulations were performed with up to 200 cardiac cycles and with data sampling rates of up to 100,000 times per cardiac cycle. Through phase-averaging or triple decomposition, the contributions of turbulence and of laminar pulsatile waves to the velocity, pressure and wall shear stress fluctuations were distinguished. For example, the commonly used oscillatory shear index was found to be closely related to the laminar waves introduced at the inlet, rather than turbulence. The turbulence energy cascade was evaluated through energy spectrum estimates, revealing that, despite the low flow rates and Reynolds number, the flow is turbulent near the aneurysm. Phase-averaging was shown to be an approach that can help researchers better understand this flow, although the results are highly dependent on simulation setup and post-processing choices., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Most children with hemiplegic cerebral palsy (HCP), one of the most prevalent subtypes of cerebral palsy, struggle with grasping and manipulating objects. This impairment may arise from a diminished capacity to properly direct forces created with the finger pad due to aberrant force application. Children with HCP were asked to create maximal force with the index finger pad in the palmar (normal) direction with both the paretic and non-paretic hands. The resulting forces and finger postures were then applied to a computational musculoskeletal model of the hand to estimate the corresponding muscle activation patterns. Subjects tended to create greater shear force relative to normal force with the paretic hand (p < 0.05). The resultant force was directed 33.6°±10.8° away from the instructed palmar direction in the paretic hand, but only 8.0°±7.3° in the non-paretic hand. Additionally, participants created greater palmar force with the non-paretic hand than with the paretic hand (p < 0.05). These differences in force production are likely due to differences in muscle activation pattern, as our computational models showed differences in which muscles are active and their relative activations when recreating the measured force vectors for the two hands (p < 0.01). The models predicted reduced activation in the extrinsic and greater reductions in activation in the intrinsic finger muscles, potentially due to reduced voluntary activation or muscle atrophy. As the large shear forces could lead to objects slipping from grasp, muscle activation patterns may provide an important target for therapeutic treatment in children with HCP., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

To better understand the impact of valvular heart disease (VHD) on the hemodynamics of the circulatory system, investigations can be carried out using a model of the cardiovascular system. In this study, a previously developed hybrid (hydro-numerical) simulator of the cardiovascular system (HCS) was adapted and used. In our HCS Björk-Shiley mechanical heart valves were used, playing the role of mitral and aortic ones. In order to simulate aortic stenosis (AS) and mitral regurgitation (MR), special mechanical devices have been developed and integrated with the HCS. The simulation results proved that the system works correctly. Namely, in the case of AS - the mean pulmonary arterial pressure was increased due to increased preload of the left ventricle and the decrease in right ventricular preload was caused by a decrease in systemic arterial pressure. The severity of AS was performed based on the transaortic pressure gradient as well as using the Gorlin and Aaslid equations. In the case of severe AS, when the mean gradient was above 40 mmHg, the aortic valve orifice area was 0.5 cm 2 , which is in line with ACC/AHA guidelines. For the case of MR - with increasing severity of MR, there was a decrease in the left ventricular pressure and an increase in left atrial pressure. Using mechanical heart valves to simulate VHD by the HCS can be a valuable tool for biomedical research, providing a safe and controlled environment to study and understand the pathophysiology of VHD., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

In the Ross procedure, a patient's pulmonary valve is transplanted in the aortic position. Despite advantages of this surgery, reoperation is still needed in many cases due to excessive dilatation of the pulmonary autograft. To further understand the failure mechanisms, we propose a multiscale model predicting adaptive processes in the autograft at the cell and tissue scale. The cell-scale model consists of a network model, that includes important signaling pathways and relations between relevant transcription factors and their target genes. The resulting gene activity leads to changes in the mechanical properties of the tissue, modeled as a constrained mixture of collagen, elastin and smooth muscle. The multiscale model is calibrated with findings from experiments in which seven sheep underwent the Ross procedure. The model is then validated against a different set of sheep experiments, for which a qualitative agreement between model and experiment is found. Model outcomes at the cell scale, including the activity of genes and transcription factors, also match experimentally obtained transcriptomics data., Competing Interests: Declaration of competing interest The authors have no competing interests to declare that are relevant to the content of this article., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Musculoskeletal simulations with muscle optimization aim to minimize muscle effort, hence are considered unable to predict the activation of antagonistic muscles. However, activation of antagonistic muscles might be necessary to satisfy the dynamic equilibrium. This study aims to elucidate under which conditions coactivation can be predicted, to evaluate factors modulating it, and to compare the antagonistic activations predicted by the lumbar spine model with literature data. Simple 2D and 3D models, comprising of 2 or 3 rigid bodies, with simple or multi-joint muscles, were created to study conditions under which muscle coactivity is predicted. An existing musculoskeletal model of the lumbar spine developed in AnyBody was used to investigate the effects of modeling intra-abdominal pressure (IAP), linear/cubic and load/activity-based muscle recruitment criterion on predicted coactivation during forward flexion and lateral bending. The predicted antagonist activations were compared to reported EMG data. Muscle coactivity was predicted with simplified models when multi-joint muscles were present or the model was three-dimensional. During forward flexion and lateral bending, the coactivation ratio predicted by the model showed good agreement with experimental values. Predicted coactivation was negligibly influenced by IAP but substantially reduced with a force-based recruitment criterion. The conditions needed in multi-body models to predict coactivity are: three-dimensionality or multi-joint muscles, unless perfect antagonists. The antagonist activations are required to balance 3D moments but do not reflect other physiological phenomena, which might explain the discrepancies between model predictions and experimental data. Nevertheless, the findings confirm the ability of the multi-body trunk models to predict muscle coactivity and suggest their overall validity., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Congenital arterial stenosis such as supravalvar aortic stenosis (SVAS) are highly prevalent in Williams syndrome (WS) and other arteriopathies pose a substantial health risk. Conventional tools for severity assessment, including clinical findings and pressure gradient estimations, often fall short due to their susceptibility to transient physiological changes and disease stage influences. Moreover, in the pediatric population, the severity of these and other congenital heart defects (CHDs) often restricts the applicability of invasive techniques for obtaining crucial physiological data. Conversely, evaluating CHDs and their progression requires a comprehensive understanding of intracardiac blood flow. Current imaging modalities, such as blood speckle imaging (BSI) and four-dimensional magnetic resonance imaging (4D MRI) face limitations in resolving flow data, especially in cases of elevated flow velocities. To address these challenges, we devised a computational framework employing zero-dimensional (0D) lumped parameter models coupled with patient-specific reconstructed geometries pre- and post-surgical intervention to execute computational fluid dynamic (CFD) simulations. This framework facilitates the analysis and visualization of intricate blood flow patterns, offering insights into geometry and flow dynamics alterations impacting cardiac function. In this study, we aim to assess the efficacy of surgical intervention in correcting an extreme aortic defect in a patient with WS, leading to reductions in wall shear stress (WSS), maximum velocity magnitude, pressure drop, and ultimately a decrease in cardiac workload., Competing Interests: Declaration of competing interest The authors declare no conflicts of interest., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Simulations of musculoskeletal models are useful for estimating internal muscle and joint forces. However, predicted forces rely on optimization and modeling formulations. Geometric detail is important to predict muscle forces, and greater geometric complexity is required for muscles that have broad attachments or span many joints, as in the torso. However, the extent to which optimized muscle force recruitment is sensitive to these geometry choices is unclear. We developed level, uphill and downhill sloped walking simulations using a standard (uniformly weighted, "fatigue-like") cost function with lower limb and full-body musculoskeletal models to evaluate hip muscle recruitment with different geometric representations of the psoas muscle under walking conditions with varying hip moment demands. We also tested a novel cost function formulation where muscle activations were weighted according to the modeled geometric detail in the full-body model. Total psoas force was less and iliacus, rectus femoris, and other hip flexors' force was greater when psoas was modeled with greater geometric detail compared to other hip muscles for all slopes. The proposed weighting scheme restored hip muscle force recruitment without sacrificing detailed psoas geometry. In addition, we found that lumbar, but not hip, joint contact forces were influenced by psoas force recruitment. Our results demonstrate that static optimization dependent simulations using models comprised of muscles with different amounts of geometric detail bias force recruitment toward muscles with less geometric detail. Muscle activation weighting that accounts for differences in geometric complexity across muscles corrects for this recruitment bias., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

The inverse dynamics based musculoskeletal simulation needs ground reaction forces (GRF) as an external force input. GRF can be predicted from kinematic data. However, the validity of estimated muscle activation using the predicted GRF has remained unclear. Therefore, the purpose of this study was to determine the validity of estimated muscle activation with predicted GRF in the inverse dynamics based musculoskeletal simulation. To perform musculoskeletal simulations, an open-source motion capture dataset that contains gait data from 50 healthy subjects was used. CusToM was used for the musculoskeletal simulations. Two sets of inverse dynamics and static optimization were performed, one used predicted GRF (PRED) and another used experimentally measured GRF (EXP). Pearson's correlation was calculated to evaluate the similarity between EMG and estimated muscle activations for both PRED and EXP. To compare PRED and EXP, paired t-tests were used to compare the trial-wise muscle activation similarity and residuals. Relationships between joint moments and residuals were also tested. The overall muscle activation similarity was comparable in PRED (R = 0.477) and EXP (R = 0.475). The residuals were 2-4 times higher in EXP compared to PRED (P < 0.001). The hip flexion-extension moment was correlated to sagittal plane residual moment (R = 0.467). The muscle activations estimated using predicted GRF were comparable to that with measured GRF in the inverse dynamics based musculoskeletal simulation. Prediction of GRF helps to perform musculoskeletal simulations where the force plates are not available., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Ryo Ueno and Yasuaki Tsuyuki are employee of ORGO Inc., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

The progression of cardiovascular disease is intricately influenced by a complex interplay between physiological pathways, biochemical processes, and physical mechanisms. This study aimed to develop an in-silico physics-based approach to comprehensively model the multifaceted vascular pathophysiological adaptations. This approach focused on capturing the progression of proximal pulmonary arterial hypertension, which is significantly associated with the irreversible degradation of arterial walls and compensatory stress-induced growth and remodeling. This study incorporated critical characteristics related to the distinct time scales for the deformation, thus reflecting the impact of mean pressure on artery growth and tissue damage. The in-silico simulation of the progression of pulmonary hypertension was realized based on computational code combined with the finite element method (FEM) for the simulation of disease progression. The parametric studies further explored the consequences of these irreversible processes. This computational modeling approach may advance our understanding of pulmonary hypertension and its progression., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Clinical management of whiplash-associated disorders is challenging and often unsuccessful, with over a third of whiplash injuries progressing to chronic neck pain. Previous imaging studies have identified muscle fat infiltration, indicative of muscle weakness, in the deep cervical extensor muscles (multifidus and semispinalis cervicis). Yet, kinematic and muscle redundancy prevent the direct assessment of individual neck muscle strength, making it difficult to determine the role of these muscles in motor dysfunction. The purpose of this study was to determine the effects of deep cervical extensor muscle weakness on multi-directional neck strength and muscle activation patterns. Maximum isometric forces and associated muscle activation patterns were computed in 25 test directions using a 3-joint, 24-muscle musculoskeletal model of the head and neck. The computational approach accounts for differential torques about the upper and lower cervical spine. To facilitate clinical translation, the test directions were selected based on locations where resistance could realistically be applied to the head during clinical strength assessments. Simulation results reveal that the deep cervical extensor muscles are active and contribute to neck strength in directions with an extension component. Weakness of this muscle group leads to complex compensatory muscle activation patterns characterized primarily by increased activation of the superficial extensors and deep upper cervical flexors, and decreased activation of the deep upper cervical extensors. These results provide a biomechanistic explanation for movement dysfunction that can be used to develop targeted diagnostics and treatments for chronic neck pain in whiplash-associated disorders., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

The development and application of multi-scale models of the lung has significantly increased in recent years. These hybrid models merge realistic representations of the larger airways with lower-dimensional descriptions of the bronchioles and respiratory airways. Due to recent advancements, it is possible to calculate airflow and dosimetry throughout the entire lung, enabling model validation with human or animal data. Here, we present a hybrid modeling pipeline and corresponding characteristic airflow and particle deposition hotspots. Next, we discuss the limitations of current hybrid models, including the need to update lower-dimensional deposition function descriptions to better represent realistic airway geometries. Future directions should include modeling diseased lungs and use of machine learning to predict whole lung dosimetry maps for a wider population., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

The healthy adult aorta is a remarkably resilient structure, able to resist relentless cardiac-induced and hemodynamic loads under normal conditions. Fundamental to such mechanical homeostasis is the mechano-sensitive cell signaling that controls gene products and thus the structural integrity of the wall. Mouse models have shown that smooth muscle cell-specific disruption of transforming growth factor-beta (TGFβ) signaling during postnatal development compromises this resiliency, rendering the aortic wall susceptible to aneurysm and dissection under normal mechanical loading. By contrast, disruption of such signaling in the adult aorta appears to introduce a vulnerability that remains hidden under normal loading, but manifests under increased loading as experienced during hypertension. We present a multiscale (transcript to tissue) computational model to examine possible reasons for compromised mechanical homeostasis in the adult aorta following reduced TGFβ signaling in smooth muscle cells., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Ltd.)

To investigate the optimal cutting depth (Cap) in small incision lenticule extraction from the perspective of corneal biomechanics, a three-dimensional finite element model of the cornea was established using a stromal sub-regional material model to simulate small incision lenticule extraction. The displacement difference P Δ at the central point of the posterior corneal surface before and after lenticule extraction, as well as the von Mises stress at four points of different thicknesses in the center of the cornea, were analyzed using the finite element model considering the hyperelastic property and the difference in stiffness between the anterior and posterior of the cornea. The numerical curves of P Δ -Cap and von Mises Stress-Cap relations at different diopters show that the displacement difference P Δ has a smallest value at the same diopter. In this case, the von Mises stress at four points with different thicknesses in the center of the cornea was also minimal. Which means that the optimal cutting depth exsisting in the cornea. Moreover, P Δ -Cap curves for different depth of stromal stiffness boundaries show that the optimal cap thickness would change with the depth of the stromal stiffness boundary. These results are of guiding significance for accurately formulating small incision lenticule extraction surgery plans and contribute to the advancement of research on the biomechanical properties of the cornea., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Skeletal muscles are complex structures with nonlinear constitutive properties. This complexity often requires finite element (FE) modeling to better understand muscle behavior and response to activation, especially the fiber strain distributions that can be difficult to measure in vivo. However, many FE muscle models designed to study fiber strain do not include force-velocity behavior. To investigate force-velocity property impact on strain distributions within skeletal muscle, we modified a muscle constitutive model with active and passive force-length properties to include force-velocity properties. We implemented the new constitutive model as a plugin for the FE software FEBio and applied it to four geometries: 1) a single element, 2) a multiple-element model representing a single fiber, 3) a model of tapering fibers, and 4) a model representing the bicep femoris long head (BFLH) morphology. Maximum fiber velocity and boundary conditions of the finite element models were varied to test their influence on fiber strain distribution. We found that force-velocity properties in the constitutive model behaved as expected for the single element and multi-element conditions. In the tapered fiber models, fiber strain distributions were impacted by changes in maximum fiber velocity; the range of strains increased with maximum fiber velocity, which was most noted in isometric contraction simulations. In the BFLH model, maximum fiber velocity had minimal impact on strain distributions, even in the context of sprinting. Taken together, the combination of muscle model geometry, activation, and displacement parameters play a critical part in determining the magnitude of impact of force-velocity on strain distribution., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Ltd.)

Young female athletes participating in sports requiring rapid changes of direction are at heightened risk of suffering traumatic knee injury, especially noncontact rupture of the anterior cruciate ligament (ACL). Clinical studies have revealed that geometric features of the tibiofemoral joint are associated with increased risk of suffering noncontact ACL injury. However, the relationship between three-dimensional (3D) tibiofemoral geometry and knee mechanics in young female athletes is not well understood. We developed a statistically augmented computational modeling workflow to determine relationships between 3D geometry of the knee and tibiofemoral kinematics and ACL force in response to an applied loading sequence of compression, valgus, and anterior force, which is known to load the ACL. This workflow included 3D characterization of tibiofemoral bony geometry via principal component analysis and multibody dynamics models incorporating subject-specific knee geometries. A combination of geometric features of both the tibia and the femur that spanned all three anatomical planes was related to increased ACL force and to increased kinematic coupling (i.e., anterior, medial, and distal tibial translations and internal tibial rotation) in response to the applied loads. In contrast, a uniplanar measure of tibiofemoral geometry that is associated with ACL injury risk, sagittal plane slope of the lateral tibial plateau subchondral bone, was not related to ACL force. Thus, our workflow may aid in developing mechanics-based ACL injury screening tools for young, active females based on a unique combination of bony geometric features that are related to increased ACL loading., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: ‘Danyal H. Nawabi has received education payments from Arthrex, Stryker, and Smith & Nephew and hospitality payments from DePuy Synthes and Biomet Orthopedics. Thomas L. Wickiewicz receives royalties from Stryker–MAKO Surgical. Andrew D. Pearle has ownership in MyGemini, Inc., is a consultant for and receives royalties from Zimmer Biomet, has received consulting fees from Stryker and Exactech, and has received hospitality and other compensation from Smith & Nephew and Arthrex.’., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

This work presents the application of a chemo-mechano-biological constitutive model of soft tissues for describing tissue inflammatory response to damage in collagen constituents. The material model is implemented into a nonlinear finite element formulation to follow up a coronary standard balloon angioplasty for one year. Numerical results, compared with available in vivo clinical data, show that the model reproduces the temporal dynamics of vessel remodeling associated with subintimal damage. Such dynamics are bimodular, being characterized by an early tissue resorption and lumen enlargement, followed by late tissue growth and vessel constriction. Applicability of the modeling framework in retrospective studies is demonstrated, and future extension towards prospective applications is discussed., Competing Interests: Declaration of competing interest Each author has been involved in the design of the study, interpretation of the data, and writing of the manuscript. Each author has read and concurs with the content in the manuscript. The manuscript has not been and will not be submitted for publication elsewhere. All authors disclose any financial and personal relationships with other people or organizations that could inappropriately influence their work., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Variability is one of the most crucial outcomes in human movement studies: variance and standard deviation of various parameters have been reported in numerous studies. However, in many of these studies, the numbers of trials and subjects have been intuitively determined and not justified with statistical considerations. Here, we investigated the impact of the numbers of trials and subjects on statistical power, based on the assumption that results per trial follow a normal distribution, using mathematical analysis and numerical simulation. An inverse-like relationship was observed between the number of trials and subjects required to ensure the statistical power for detecting differences in variance between subject groups or conditions. For instance, assuming a 1.2-times difference in population variance between pre-and post-training sessions as an alternative hypothesis, our simulation demonstrated that combinations of the number of subjects and trials, such as measuring 100 trials from each of 12 subjects under each condition, or measuring 20 trials from each of 60 subjects, can guarantee an 80 % of statistical power. Planning research based on such mathematical considerations will enable meaningful statistical interpretations in studies focusing on movement variability, such as gait studies., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: This research was supported by the JSPS KAKENHI Grant Number 20H04089 (PI: KT). We have no additional conflicts of interest., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Current spinal testing protocols generally adopt pure moments combined with axial compression. However, daily activities involve multi-axis loads, and multi-axis loading has been shown to impact intervertebral disc (IVD) cell viability. Therefore, integrating in-vivo load data with spine simulators is critical to understand how loading affects the IVD, but doing so is challenging due to load coupling and variable load rates. This study addresses these challenges through the Load Informed Kinematic Evaluation (LIKE) protocol, which was evaluated using the root mean squared error (RMSE) between desired and actual loads in each axis. Stage 1 involves obtaining the kinematics from six-axis load control tests replicating 20 Orthoload activities at a reduced test speed. Stage 2 applies these kinematics in five axes, with axial compression applied in load control, at the reduced speed and at the physiological test rate. Stage 3 enables long-term tests through six-axis kinematic control combined with diurnal height correction to account for the natural height fluctuations of the IVD. Stage 1 yielded RMSEs within twice the load cell noise floor. Low RMSEs were maintained during stage 2 at reduced speed (Tx:0.80 ± 0.30 N; Ty:0.77 ± 0.29 N; Tz:1.79 ± 0.50 N; Rx:0.02 ± 0.01Nm; Ry:0.02 ± 0.01Nm; and Rz:0.02 ± 0.01Nm) and at the physiological test rate (Tx:3.45 ± 1.81 N; Ty:3.82 ± 1.99 N; Tz:11.32 ± 8.69 N; Rx:0.13 ± 0.07Nm; Ry:0.16 ± 0.11Nm; and Rz:0.07 ± 0.04Nm). To address unwanted oscillations observed in longer tests (>2h), Stage 3 was introduced to enable the stable and consistent replication of activities at a physiological test rate. Despite higher RMSEs the axial error was 85.5 ± 24.27 N (equivalent to ∼ 0.16 MPa), with shear RMSEs similar to other testing systems conducting pure moment tests at slower rates. The LIKE protocol enables the replication of physiological loads, providing opportunities for enhanced investigations of IVD mechanobiology, and the pre-clinical evaluation of IVD devices and therapies., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Single bicycle crashes, i.e., falls and impacts not involving a collision with another road user, are a significantly underestimated road safety problem. The motions and behaviours of falling people, or fall kinematics, are often investigated in the injury biomechanics research field. Understanding the mechanics of a fall can help researchers develop better protective gear and safety measures to reduce the risk of injury. However, little is known about cyclist fall kinematics or dynamics. Therefore, in this study, a video analysis of cyclist falls is performed to investigate common kinematic forms and impact patterns. Furthermore, a pipeline involving deep learning-based human pose estimation and inverse kinematics optimisation is created for extracting human motion from real-world footage of falls to initialise forward dynamics computational human body models. A bracing active response is then optimised for using a genetic algorithm. This is then applied to a case study of a cyclist fall. The kinematic forms characterised in this study can be used to inform initial conditions for computational modelling and injury estimation in cyclist falls. Findings indicate that protective response is an important consideration in fall kinematics and dynamics, and should be included in computational modelling. Furthermore, the novel reconstruction pipeline proposed here can be applied more broadly for traumatic injury biomechanics tasks. The tool developed in this study is available at https://kevgildea.github.io/KinePose/., Competing Interests: Declaration of competing interest The authors declare the following potential competing interest: Ciaran Simms is a co-founder of a start-up company that uses pose estimation for sports exercise quantification. Apart from this disclosure, the other authors declare that they have no additional financial and non-financial interests., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Diabetes mellitus (DM) is a predisposing risk factor leading to macrovascular diseases. Changes in haemodynamics of the diabetic aortas remain largely unclear and relevant computational analyses are lacking in the literature. Ten adult rabbits (1.6–2.2 kg) were collected and the type I diabetic rabbit model was induced by injection of alloxan. A total of five control and five diabetic rabbit aortas were considered for subsequent numerical simulation. The CT scanning was performed to reconstruct three-dimensional model of the individual rabbit descending aorta. The flow velocity waveforms were measured by ultrasound machine and were set to be the inlet boundary conditions. The reconstructed aortas were then imported into ANSYS to perform mesh generation and computational analysis. Results showed that the distributions of haemodynamic indicators time-averaged wall shear stress (TAWSS), oscillating shear index (OSI) and transverse wall shear stress (transWSS) in the non-diabetic rabbit aortas were similar to those in the diabetic rabbit aortas. However, the mean values of TAWSS and transWSS in the non-diabetic rabbit aortas were significantly higher than those values in the diabetic rabbit aortas (TAWSS: p = 0.04; transWSS: p = 0.02). The back of right renal artery tended to have high OSI in both the non-diabetic and the diabetic rabbit aortas. Notably, the regions with high OSI tended to have intense disturbed flow and low TAWSS in the most diabetic rabbit aortas. The results suggest that diabetes leads to changes in haemodynamic parameters in the rabbit aortas. In particular, the lower TAWSS and the higher OSI within the diabetic aortas may further contribute to aortic wall remodeling. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Nucleotomy is the gold standard treatment for disc herniation and has proven ability to restore stability by creating a bony bridge without any additional fixation. However, the evolution of mineral density in the extant and new bone after nucleotomy and fixation techniques has to date not been investigated in detail. The main goal of this study is to determine possible mechanisms that may trigger the bone remodelling and formation processes. With that purpose, a finite element model of the L4–L5 spinal segment was used. Bone mineral density (BMD), new tissue composition, and endplate deflection were determined as indicators of lumbar fusion. A bone-remodelling algorithm and a tissue-healing algorithm, both mechanically driven, were implemented to predict vertebral bone alterations and fusion patterns after nucleotomy, internal fixation, and anterior plate placement. When considering an intact disc height, neither nucleotomy nor internal fixation were able to provide the necessary stability to promote bony fusion. However, when 75% of the disc height was considered, bone fusion was predicted for both techniques. By contrast, an anterior plate allowed bone fusion at all disc heights. A 50% disc-height reduction led to osteophyte formation in all cases. Changes in the intervertebral disc tissue caused BMD alterations in the endplates. From this observations it can be drawn that fusion may be self-induced by controlling the mechanical stabilisation without the need of additional fixation. The amount of tissue to be removed to achieve this stabilisation remains to be determined. [ABSTRACT FROM AUTHOR]

Ski jumping performance is strongly affected by wind. Flight technique optimization for maximizing jump length is a highly complex motor-control task that also depends on the wind. Pontryagin's minimum principle was used in this study to gain a better understanding on how wind influences flight technique optimization. Optimum time courses of the angle of attack α of the skis and of the body-to-ski angle β were computed for seven realistic wind scenarios on the large hill and on the flying hill. The optimum values of α were smaller at headwind, and larger at tailwind when compared to the optimum time course at calm wind. The optimum values of β were the smallest possible ones at the given flight technique constraints, except for the last part of the flight. Optimum adjustments of α increased the jump lengths between 0 and 1.8 m on the large hill, and between 0 and 6.4 m on the flying hill. Maximum jump length increases were achieved at the highest headwind speed. Even larger jump length effects can be achieved by using smaller β -angles, which might be possible in headwind conditions, but this is associated with increased problems to keep the flight stable. [ABSTRACT FROM AUTHOR]

Abstract Analysis of rare cells in heterogenous mixtures is proven to be beneficial for regenerative medicine, cancer treatment and prenatal diagnostics. Scarcity of these cells, however, makes the isolation process extremely challenging. Efficiency in cell isolation is still low and therefore, novel cell isolation strategies with new biomarkers need exploration. In this study, we investigated the feasibility of using the mechanical stiffness difference to detect and isolate the rare cells from the surrounding cells without labelling them. Fluid and solid mechanics simulations have shown that cell isolation can be performed at high efficiency using stiffness-based isolation. Accuracy of the numerical simulations is established using microfluidic flow chamber experiments. [ABSTRACT FROM AUTHOR]

Abstract New computational techniques providing more accurate representation of human heart pathologies could help uncovering relevant physical phenomena and improve the outcome of medical therapies. In this framework, the present work describes an efficient computational model for the evaluation of the ventricular flow alteration in presence of mitral valve stenosis. The model is based on the direct numerical simulation of the Navier–Stokes equations two-way coupled with a structural solver for the left ventricle and mitral valve dynamics. The presence of mitral valve stenosis is mimicked by a single-parameter constraint acting on the kinematics of the mitral leaflets. Four different degrees of mitral valve stenosis are considered focusing on the hemodynamic alterations occurring in pathologic conditions. The mitral jet, generated during diastole, is seen to shrink and strengthen when the stenosis gets more severe. As a consequence, the kinetic energy of the flow, the tissues shear stresses, the transvalvular pressure drop and mitral regurgitation increase. It results that, as the stenosis severity level increases, the geometric and effective orifice areas decrease up to 50% with respect the normal case due to the reduced leaflets mobility and stronger blood acceleration during the diastolic phase. The modified intraventricular hemodynamics is also related to a stronger pressure gradient that, for severe stenosis, can be more than ten times larger than the healthy valve case. These computational results are fully consistent with the available clinical literature and open the way to the virtual assessment of surgical procedures and to the evaluation of prosthetic devices. [ABSTRACT FROM AUTHOR]

Abstract Background OpenSim models are typically based on cadaver findings that are generalized to represent a wide range of populations, which curbs their validity. Patient-specific modelling through incorporating magnetic resonance imaging (MRI) improves the model's biofidelity with respect to joint alignment and articulations, muscle wrapping, and ligament insertions. The purpose of this study was to determine if the inclusion of an MRI-based knee model would elicit differences in lower limb kinematics and resulting knee ligament lengths during a side cut task. Methods Eleven participants were analyzed with the popular Rajagopal OpenSim model, two variations of the same model to include three and six degrees of freedom knee (DOF), and a fourth version featuring a four DOF MRI-based knee model. These four models were used in an inverse kinematics analysis of a side cut task and the resulting lower limb kinematics and knee ligament lengths were analyzed. Results The MRI-based model was more responsive to the movement task than the original Rajagopal model while less susceptible to soft tissue artifact than the unconstrained six DOF model. Ligament isometry was greatest in the original Rajagopal model and smallest in the six DOF model. Conclusions When using musculoskeletal modelling software, one must acutely consider the model choice as the resulting kinematics and ligament lengths are dependent on this decision. The MRI-based knee model is responsive to the kinematics and ligament lengths of highly dynamic tasks and may prove to be the most valid option for continuing with late-stage modelling operations such as static optimization. [ABSTRACT FROM AUTHOR]

Abstract Modeling joint contact in OpenSim is not well understood. This study systematically investigated the variables associated with the elastic foundation contact model within OpenSim by performing a series of controlled benchtop experiment and concomitant simulations. Four metal-on-plastic interactions were modeled, including a model of a total knee replacement (TKR). Load-displacement curves were recorded during cyclic loading between 100 and 750 N. Geometries were imported and into OpenSim and contact mechanics were modeled with the on-board elastic foundation algorithm. A hybrid optimization algorithm determined that stiffness and dissipation coefficients for TKR implants were 1.52 × 1010 N/m and 57.7 Ns/m, respectively. Estimations of contact forces were 10.2% of blinded experimental data (average root mean square error: 76.82 ± 11.47 N). In the second portion of this study, freely available eTibia TKR renderings were used to test the ubiquity of the tuning parameters. They were also used to perform a sensitivity analysis of material stiffness and mesh density with regard to penetration depth and computational time. When a stiffness of 1 × 1010 was applied to an eTibia model with 5000 faces, a 100 kg load caused 0.259 mm of penetration. Under the same conditions, the tuned model experienced 0.300 mm of penetration. Material stiffnesses between 1 × 1013 and 1 × 1015 N/m increased computation time by factors of 12–23. This study provides much needed clarity regarding the use of the OpenSim EF algorithm. It also demonstrates the utility of OpenSim to model deformable materials and complex geometries, and this approach can be adapted to make reasonable estimations for both natural and surgically modified joints. [ABSTRACT FROM AUTHOR]

Abstract Background Numerous biomechanical studies have addressed normal shoulder function and the factors that affect it. While these investigations include a mix of in-vivo clinical reports, ex-vivo cadaveric studies, and computer-based simulations, each has its own strengths and limitations. A robust methodology is essential in cadaveric work but does not always come easily. Precise quantitative measurements are difficult in in-vivo studies, and simulation studies require validation steps. This review focuses on ex-vivo cadaveric studies to emphasize the best research methodologies available to simulate physiologically and clinically relevant shoulder motion. Methods A PubMed and Web of Science search was conducted in March 2017 (and updated in May 2018) to identify the cadaveric studies focused on the shoulder and its function. The key words for this search included rotator cuff (RC) injuries, RC surgery, and their synonyms. The protocol of the study was registered on PROSPERO and is accessible at CRD42017068873. Results Thirty one studies consisting of 167 specimens with various biomechanical methods met our inclusion criteria. All studies were level V cadaveric studies. Cadaveric biomechanical models are widely used to study shoulder instability and RC repair. These models are commonly limited to the glenohumeral joint by a fixed scapula, passively and discretely move the humerus, and statically load the RC without regard for the integrity of the glenohumeral capsule. Conclusion All studies captured in this review evaluated shoulder biomechanics. Recent studies in patients suggest that some assumptions made in this space may not fully characterize motion of the human shoulder. With reproducible scapular positioning, dynamic RC activation, and preservation of glenohumeral capsule integrity, cadaveric studies can facilitate proper validation for simulation models and broaden our understanding of the shoulder environment during motion in healthy and disease states. [ABSTRACT FROM AUTHOR]

Abstract The quality with which smoothing algorithms perform is often assessed in simulation by starting with a known 1D datum, adding noise, smoothing the noisy data, then quantifying the difference between the smoothed data and known datum, often using mean-square error (MSE). While effectively summarizing overall difference, MSE fails to capture localized, one-sided errors. This paper describes how smoothing noisy 1D data using a variety of algorithms can introduce systematic bias, and quantifies this bias using the false positive rate (FPR): the probability that a smoothing algorithm will yield a dataset whose 1D mean differs significantly from its true 1D datum. A simulation study was conducted involving six 1D datum continua, and four smoothing algorithms whose parameters were systematically manipulated along with sample size and noise amplitude. Approximately ten million simulation iterations were evaluated. FPRs were calculated at α = 0.05 , based on the calculated smoothness of the resulting datasets. Results showed that FPRs were much higher than the expected value of α , and in many cases approached 100%. FPRs were highest with aggressive smoothing parameters, large sample sizes and small noise amplitudes, irrespective of both smoothing algorithm and the 1D datum. These results suggest that smoothing 1D biomechanical data can introduce statistical bias with relatively high probability. The implications are experiment-specific because the biomechanical meaning of 1D changes can vary vastly between datasets. Smoothing-induced bias should be a cause for general concern when small 1D changes have non-trivial biomechanical consequences. [ABSTRACT FROM AUTHOR]

Abstract Mechanical characterization of living tissues and computer-based simulations related to medical issues, has become increasingly important to improve diagnostic processes and treatments evaluation. This work proposes a link between the mechanical testing and the material model predictions through histological data of vaginal tissue. Histological data was used to link tensile testing experiments with material-dependent parameters; the approach was adequate to capture the nonlinear response of ovine vaginal tissue over a large strain range. The experimental data obtained on a previous study, has two main components: tensile testing and histological analysis of the ovine vaginal tissue. Uniaxial tensile test data and histological data were collected from three sheep groups: virgins, pregnant and parous. The distal part of vaginal wall was selected since it is prone to tears induced by vaginal delivery. The HGO (Holzapfel-Gasser-Ogden) model parameters were fitted using a stochastic approach, namely the Simple Genetic Algorithm (SGA). The SGA was able to fit the experimental data successfully (R2 > 0.986). The dimensionless coefficient ξ , was highly correlated with histological data. The ratio was seen to increase linearly with increasing collagen content. Coefficient ξ brings a new way of interpreting and understanding experimental data; it connects the nonlinear mechanical behaviour (tensile test) with tissue's morphology (histology). It can be used as an 'inverse' (approximate) method to estimate the mechanical properties without direct experimental measurements, through basic histology. In this context, the proposed methodology appears very promising in estimating the response of the tissue via histological information. [ABSTRACT FROM AUTHOR]

Abstract Objective Poor total cavopulmonary connection (TCPC) hemodynamics have been hypothesized to be associated with long-term complications in Fontan patients. Image-based Fontan surgical planning has shown great potential as a clinical tool because it can pre-operatively evaluate patient-specific hemodynamics. Current surgical planning paradigms commonly utilize cardiac-gated phase contrast magnetic resonance (MR) imaging to acquire vessel flows. These acquisitions are often taken under breath-held (BH) conditions and ignore the effect of respiration on blood flow waveforms. This study investigates the effect of respiration-driven flow waveforms on patient-specific hemodynamics using real-time MR acquisitions. Methods Patient-specific TCPCs were reconstructed from cardiovascular MR images. Real-time phase contrast MR images were acquired under both free-breathing (FB) and breath-held conditions for 9 patients. Numerical simulations were employed to assess flow structures and hemodynamics used in Fontan surgical planning including hepatic flow distribution (HFD) and indexed power loss (iPL), which were then compared between FB and BH conditions. Results Differences in TCPC flow structures between FB and BH conditions were observed throughout the respiratory cycle. However, the average differences (BH – FB values for each patient, which are then averaged) in iPL and HFD between these conditions were 0.002 ± 0.011 (p = 0.40) and 1 ± 3% (p = 0.28), respectively, indicating no significant difference in clinically important hemodynamic metrics. Conclusions Respiration affects blood flow waveforms and flow structures, but might not significantly influence the values of iPL or HFD. Therefore, breath-held MR acquisition can be adequate for Fontan surgical planning when focusing on iPL and HFD. [ABSTRACT FROM AUTHOR]

Abstract It is clinically challenging to distinguish between ankle and subtalar joints instability in vivo. Understanding the changes in load-displacement at the ankle and subtalar joints after ligament injuries may detect specific changes in joint characteristics that cannot be detected by investigating changes in range of motion alone. The effect of restricting joints end range of motion with ankle braces was already established, but little is known about the effect of an ankle brace on the flexibility of the injured ankle and subtalar joints. Therefore, the purposes of this study were to (1) understand how flexibility is affected at the ankle and subtalar joints after sectioning lateral and intrinsic ligaments during combined sagittal foot position and inversion and during internal rotation and (2) investigate the effect of a semi-rigid ankle brace on the ankle and subtalar joint flexibility. Kinematics and kinetics were collected from nine cadaver feet during inversion through the range of ankle flexion and during internal rotation. Motion was applied with and without a brace on an intact foot and after sequentially sectioning the calcaneofibular ligament (CFL) and the intrinsic ligaments. Segmental flexibility was defined as the slope of the angle-moment curve for each 1 Nm interval. Early flexibility significantly increased at the ankle and subtalar joint after CFL sectioning during inversion. The semi-rigid ankle brace significantly decreased early flexibility at the subtalar joint during inversion and internal rotation for all ligament conditions and at the ankle joint after all ligaments were cut. [ABSTRACT FROM AUTHOR]

Abstract Impaction is required to properly seat press-fit implants and ensure initial implant stability and long term bone ingrowth, however excessive impaction or press-fit presents a high fracture risk in the acetabulum and femur. Current in-vitro impaction testing methods do not replicate the compliance of the soft tissues surrounding the hip, a factor that may be important in fracture and force prediction. This study presents the measurement of compliance of the soft tissues supporting the hip during impaction in operative conditions, and replicates these in vitro. Hip replacements were carried out on 4 full body cadavers while impact force traces and acetabular/femoral displacement were measured. Compliance was then simulated computationally using a Voigt model. These data were subsequently used to inform the design of a representative in-vitro drop rig. Effective masses of 19.7 kg and 12.7 kg, spring stiffnesses of 8.0 kN/m and 4.1 kN/m and dashpot coefficients of 595 N s/m and 322 N s/m were calculated for the acetabular and femoral soft tissues respectively. A good agreement between cadaveric and in-vitro peak displacement and rise time during impact is found. Such an in-vitro setup is of use during laboratory testing, simulation or even surgical training. [ABSTRACT FROM AUTHOR]

Abstract Injuries to the anterior cruciate ligament (ACL) and menisci commonly lead to early onset osteoarthritis. Treatments that can restore normative cartilage loading patterns may mitigate the risk of osteoarthritis, though it is unclear whether such a goal is achievable through conservative rehabilitation. We used musculoskeletal simulation to predict cartilage and ligament loading patterns during walking in intact, ACL deficient, menisci deficient, and ACL-menisci deficient knees. Stochastic simulations with varying coordination strategies were then used to test whether neuromuscular control could be modulated to restore normative knee mechanics in the pathologic conditions. During early stance, a 3 mm increase in anterior tibial translation was predicted in the ACL deficient knee. Mean cartilage contact pressure increased by 18% and 24% on the medial and lateral plateaus, respectively, in the menisci deficient knee. Variations in neuromuscular coordination were insufficient to restore normative cartilage contact patterns in either the ACL or menisci deficient knees. Elevated cartilage contact pressures in the pathologic knees were observed in regions where cartilage wear patterns have previously been reported. These results suggest that altered cartilage tissue loading during gait may contribute to region-specific degeneration patterns, and that varying neuromuscular coordination in isolation is unlikely to restore normative knee mechanics. [ABSTRACT FROM AUTHOR]