Journal of Biomechanical Engineering-Transactions of the Asme最新文献

Arterial stiffening occurs during natural aging, is associated with an increased risk of adverse cardiovascular events, and can follow different timelines in males and females. One mechanism of arterial stiffening includes remodeling of the extracellular matrix (ECM), which alters the wall material properties. We used elastin haploinsufficient (Eln+/-) and wildtype (Eln+/+) mice to investigate how material properties of two different arteries (ascending aorta and carotid artery) change with age, sex, and ECM composition. We used a constitutive model by Dong and Sun that is based on the Holzapfel-Gasser-Ogden (HGO) type, but does not require a discrete number of fibrous ECM families and allows varied deformation coupling. We find that the amount of deformation coupling for the best fit model depends on the artery type. We also find that remodeling to maintain homeostatic (i.e., young, wildtype) values of biomechanical parameters with age, sex, and ECM composition depends on the artery type, with ascending aorta being more adaptable than carotid artery. Fitted material constants indicate sex-dependent remodeling that may be important for determining the time course of arterial stiffening in males and females. We correlated fitted material constants with ECM composition measured by biochemical (ascending aorta) or histological (carotid artery) methods. We show significant correlations between ECM composition and material parameters for the mean values for each group, with biochemical measurements correlating more strongly than histological measurements. Understanding how arterial stiffening depends on age, sex, ECM composition, and artery type may help design effective, personalized clinical treatment strategies.

Many methods to computationally predict red blood cell damage have been introduced, and among these are Lagrangian methods that track the cells along their pathlines. Such methods typically do not explicitly include cell-cell interactions. Due to the high volume fraction of red blood cells (RBCs) in blood, these interactions could impact cell mechanics and thus the amount of damage caused by the flow. To investigate this question, cell-resolved simulations of red blood cells in shear flow were performed for multiple interacting cells, as well as for single cells in unbounded flow at an effective viscosity. Simulations run without adjusting the bulk viscosity produced larger errors unilaterally and were not considered further for comparison. We show that a periodic box containing at least 8 cells and a spherical harmonic of degree larger than 10 are necessary to produce converged higher-order statistics. The maximum difference between the single-cell and multiple-cell cases in terms of peak strain was 3.7%. To achieve this, one must use the whole blood viscosity and average over multiple cell orientations when adopting a single-cell simulation approach. The differences between the models in terms of average strain were slightly larger (maximum difference of 6.9%). However, given the accuracy of the single-cell approach in predicting the maximum strain, which is useful in hemolysis prediction, and its computational cost that is orders of magnitude less than the multiple-cell approach, one may use it as an affordable cell-resolved approach for hemolysis prediction.

Sensitivity coefficients are used to understand how errors in subject-specific musculoskeletal model parameters influence model predictions. Previous sensitivity studies in the lower limb calculated sensitivity using perturbations that do not fully represent the diversity of the population. Hence, the present study performs sensitivity analysis in the upper limb using a large synthetic dataset to capture greater physiological diversity. The large dataset (n = 401 synthetic subjects) was created by adjusting maximum isometric force, optimal fiber length, pennation angle, and bone mass to induce atrophy, hypertrophy, osteoporosis, and osteopetrosis in two upper limb musculoskeletal models. Simulations of three isometric and two isokinetic upper limb tasks were performed using each synthetic subject to predict muscle activations. Sensitivity coefficients were calculated using three different methods (two point, linear regression, and sensitivity functions) to understand how changes in Hill-type parameters influenced predicted muscle activations. The sensitivity coefficient methods were then compared by evaluating how well the coefficients accounted for measurement uncertainty. This was done by using the sensitivity coefficients to predict the range of muscle activations given known errors in measuring musculoskeletal parameters from medical imaging. Sensitivity functions were found to best account for measurement uncertainty. Simulated muscle activations were most sensitive to optimal fiber length and maximum isometric force during upper limb tasks. Importantly, the level of sensitivity was muscle and task dependent. These findings provide a foundation for how large synthetic datasets can be applied to capture physiologically diverse populations and understand how model parameters influence predictions.

There is little information on the layer-specific failure properties of the adult human abdominal aorta, and there has been no quantification of postfailure damage. Infra-renal aortas were thus taken from forty-seven autopsy subjects and cut into 870 intact-wall and layer strips that underwent uni-axial-tensile testing. Intact-wall failure stress did not differ significantly (p > 0.05) from the medial value longitudinally, nor from the intimal and medial values circumferentially, which were the lowest recorded values. Intact-wall failure stretch did not differ (p > 0.05) from the medial value in either direction. Intact-wall prefailure stretch (defined as failure stretch-stretch at the initiation of the concave phase of the stress-stretch response) did not differ (p > 0.05) from the intimal and medial values, and intact-wall postfailure stretch (viz., full-rupture stretch-failure stretch) did not differ (p > 0.05) from the adventitial value since the adventitia was the last layer to rupture, being most extensible albeit under residual tension. Intact-wall failure stress and stretch declined from 20 to 60 years, explained by steady declines throughout the lifetime of their medial counterparts, implicating beyond 60 years the less age-varying failure properties of the intima under minimal residual compression. The positive correlation of postfailure stretch with age counteracted the declining failure stretch, serving as a compensatory mechanism against rupture. Hypertension, diabetes, and coronary artery disease adversely affected the intact-wall and layer-specific failure stretches while increasing stiffness.

In recent years, transcatheter edge-to-edge repair (TEER) has been widely adopted as an effective treatment for mitral regurgitation (MR). The aim of this study is to develop a personalized in silico model to predict the effect of edge-to-edge repair in advance to the procedure for each individual patient. For this purpose, we propose a combination of a valve deformation model for computing the mitral valve (MV) orifice area (MVOA) and a lumped parameter model for the hemodynamics, specifically mitral regurgitation volume (RVol). Although we cannot obtain detailed information on the three-dimensional flow field near the mitral valve, we can rapidly simulate the important medical parameters for the clinical decision support. In the present method, we construct the patient-specific pre-operative models by using the parameter optimization and then simulate the postoperative state by applying the additional clipping condition. The computed preclip MVOAs show good agreement with the clinical measurements, and the correlation coefficient takes 0.998. In addition, the MR grade in terms of RVol also has good correlation with the grade by ground truth MVOA. Finally, we try to investigate the applicability for the predicting the postclip state. The simulated valve shapes clearly show the well-known double orifice and the improvement of the MVOA, compared with the preclip state. Similarly, we confirmed the improved reverse flow and MR grade in terms of RVol. A total computational time is approximately 8 h by using general-purpose PC. These results obviously indicate that the present in silico model has good capability for the assessment of edge-to-edge repair.

We compared the ability of seven machine learning algorithms to use wearable inertial measurement unit (IMU) data to identify the severe knee loading cycles known to induce microdamage associated with anterior cruciate ligament rupture. Sixteen cadaveric knee specimens, dissected free of skin and muscle, were mounted in a rig simulating standardized jump landings. One IMU was located above and the other below the knee, the applied three-dimensional action and reaction loads were measured via six-axis load cells, and the three-dimensional knee kinematics were also recorded by a laboratory motion capture system. Machine learning algorithms were used to predict the knee moments and the tibial and femur vertical forces; 13 knees were utilized for training each model, while three were used for testing its accuracy (i.e., normalized root-mean-square error) and reliability (Bland-Altman limits of agreement). The results showed the models predicted force and knee moment values with acceptable levels of error and, although several models exhibited some form of bias, acceptable reliability. Further research will be needed to determine whether these types of models can be modified to attenuate the inevitable in vivo soft tissue motion artifact associated with highly dynamic activities like jump landings.

Hemolysis persists as a common and serious problem for neonatal patients on extracorporeal membrane oxygenation (ECMO). Since the cannula within the ECMO circuit is associated with hemolysis-inducing shear stresses, real-world internal fluid flow measurements are urgently needed to understand the mechanism and confirm computational estimates. This study appears to be the first experimental study of fluid flow inside commercial ECMO dual-lumen cannulas (DLCs) and first particle image velocimetry (PIV) visualization inside a complicated medical device. The internal geometries of four different opaque neonatal DLCs, both atrial and bicaval positioning geometries each sized 13 Fr and 16 Fr, were replicated by three-dimensional printing clear lumen scaled-up models, which were integrated in a circuit with appropriate ECMO flow parameters. PIV was then used to visualize two-dimensional fluid flow in a single cross section within the models. An empirical model accounting for shear stress and exposure time was used to compare the maximum expected level of hemolysis through each model. The maximum measured peak shear stress recorded was 16±2 Pa in the top arterial bicaval 13 Fr model. The atrial and 16 Fr cannula models never produced greater single-pass peak shear stress or hemolysis than the bicaval and 13 Fr models, respectively, and no difference was found in hemolysis at two different flow rates. After 5 days of flow, small DLC-induced hemolysis values for a single pass through each cannula were modeled to linearly accumulate and caused the most severe hemolysis in the bicaval 13 Fr DLC. Engineering and clinical solutions to improve cannula safety are proposed.

In recent conflicts, blast injury from landmines and improvised explosive devices (IEDs) has been the main mechanism of wounding and death. When a landmine or IED detonates under a vehicle (an under-body blast), the seat acceleration rapidly transmits a high load to the pelvis of the occupants, resulting in torso and pelvic injury. Pelvic fractures have high mortality rates, yet their injury mechanism has been poorly researched. Three (3) fresh-frozen male pelvic specimens were tested under axial impact loading. The pelvis was impacted mounted upside down by dropping a 12 kg mass at target impact velocities ranging from 1 to 8.6 m/s with time to peak velocity ranging from 3.8 to 5.8 ms. Resulting fractures were broadly categorized as involving a bilateral pubis rami fracture, a bilateral ischium fracture, and sacroiliac joint disruption. The study provides insights into the type and severity of pelvic injury that may occur over a range of under-body blast (UBB)-relevant loading profiles.

Acquired subglottal stenosis is an unpredicted complication that can occur in some patients who have undergone prolonged endotracheal intubation. It is a narrowing of the airway at the level of the cricoid cartilage that can restrict airflow and cause breathing difficulty. Stenosis is typically treated with endoscopic airway dilation, with some patients experiencing multiple recurrences. The study highlights the potential of computational fluid dynamics as a noninvasive method for monitoring subglottic stenosis, which can aid in early diagnosis and surgical planning. An anatomically accurate human laryngeal airway model was constructed from computerized tomography (CT) scans. The subglottis cross-sectional area was narrowed systematically using ≈10% decrements. A quadratic profile was used to interpolate the transformation of the airway geometry from its modified shape to the baseline geometry. The numerical results were validated by static pressure measurements conducted in a physical model. The results show that airway resistance follows a squared ratio that is inversely proportional to the size of the subglottal opening (R∝A-2). The study found that critical constriction occurs in the subglottal region at 70% stenosis (upper end of grade 2). Moreover, removing airway tissue below 40% stenosis during surgical intervention does not significantly decrease airway resistance.

Functional luminal imaging probe (FLIP) is used to measure cross-sectional area (CSA) and pressure at sphincters. It consists of a catheter surrounded by a fluid filled cylindrical bag, closed on both ends. Plotting the pressure-CSA hysteresis of a sphincter during a contraction cycle, which is available through FLIP testing, offers information on its functionality, and can provide diagnostic insights. However, limited work has been done to explain the mechanics of these pressure-CSA loops. This work presents a consolidated picture of pressure-CSA loops of different sphincters. Clinical data reveal that although sphincters have a similar purpose (controlling the flow of liquids and solids by opening and closing), two different pressure-CSA loop patterns emerge: negative slope loop (NSL) and positive slope loop (PSL). We show that the loop type is the result of an interplay between (or lack thereof) two mechanical modes: (i) neurogenic mediated relaxation of the sphincter muscle or pulling applied by external forces, and (ii) muscle contraction proximal to the sphincter which causes mechanical distention. We conclude that sphincters which only function through mechanism (i) exhibition NSL whereas sphincters which open as a result of both (i) and (ii) display a PSL. This work provides a fundamental mechanical understanding of human sphincters. This can be used to identify normal and abnormal phenotypes for the different sphincters and help in creating physiomarkers based on work calculation.