Twist buckling of veins under torsional loading.

Veins are often subjected to torsion and twisted veins can hinder and disrupt normal blood flow but their mechanical behavior under torsion is poorly understood. The objective of this study was to investigate the twist deformation and buckling behavior of veins under torsion. Twist buckling tests were performed on porcine internal jugular veins (IJVs) and human great saphenous veins (GSVs) at various axial stretch ratio and lumen pressure conditions to determine their critical buckling torques and critical buckling twist angles. The mechanical behavior under torsion was characterized using a two-fiber strain energy density function and the buckling behavior was then simulated using finite element analysis. Our results demonstrated that twist buckling occurred in all veins under excessive torque characterized by a sudden kink formation. The critical buckling torque increased significantly with increasing lumen pressure for both porcine IJV and human GSV. But lumen pressure and axial stretch had little effect on the critical twist angle. The human GSVs are stiffer than the porcine IJVs. Finite element simulations captured the buckling behavior for individual veins under simultaneous extension, inflation, and torsion with strong correlation between predicted critical buckling torques and experimental data (R2=0.96). We conclude that veins can buckle under torsion loading and the lumen pressure significantly affects the critical buckling torque. These results improve our understanding of vein twist behavior and help identify key factors associated in the formation of twisted veins.

The Effect of Trabeculae Carneae on Left Ventricular Diastolic Compliance: Improvement in Compliance With Trabecular Cutting.

The role of trabeculae carneae in modulating left ventricular (LV) diastolic compliance remains unclear. The objective of this study was to determine the contribution of trabeculae carneae to the LV diastolic compliance. LV pressure-volume compliance curves were measured in six human heart explants from patients with LV hypertrophy at baseline and following trabecular cutting. The effect of trabecular cutting was also analyzed with finite-element model (FEM) simulations. Our results demonstrated that LV compliance improved after trabecular cutting (p < 0.001). Finite-element simulations further demonstrated that stiffer trabeculae reduce LV compliance further, and that the presence of trabeculae reduced the wall stress in the apex. In conclusion, we demonstrate that integrity of the LV and trabeculae is important to maintain LV stiffness and loss in trabeculae leads to more LV compliance.

Artery buckling affects the mechanical stress in atherosclerotic plaques.

BACKGROUND: Tortuous arteries are often seen in patients with hypertension and atherosclerosis. While the mechanical stress in atherosclerotic plaque under lumen pressure has been studied extensively, the mechanical stability of atherosclerotic arteries and subsequent effect on the plaque stress remain unknown. To this end, we investigated the buckling and post-buckling behavior of model stenotic coronary arteries with symmetric and asymmetric plaque. METHODS: Buckling analysis for a model coronary artery with symmetric and asymmetric plaque was conducted using finite element analysis based on the dimensions and nonlinear anisotropic materials properties reported in the literature. RESULTS: Artery with asymmetric plaque had lower critical buckling pressure compared to the artery with symmetric plaque and control artery. Buckling increased the peak stress in the plaque and led to the development of a high stress concentration in artery with asymmetric plaque. Stiffer calcified tissue and severe stenosis increased the critical buckling pressure of the artery with asymmetric plaque. CONCLUSIONS: Arteries with atherosclerotic plaques are prone to mechanical buckling which leads to a high stress concentration in the plaques that can possibly make the plaques prone to rupture.

Trabecular plates and rods determine elastic modulus and yield strength of human trabecular bone.

The microstructure of trabecular bone is usually perceived as a collection of plate-like and rod-like trabeculae, which can be determined from the emerging high-resolution skeletal imaging modalities such as micro-computed tomography (µCT) or clinical high-resolution peripheral quantitative CT (HR-pQCT) using the individual trabecula segmentation (ITS) technique. It has been shown that the ITS-based plate and rod parameters are highly correlated with elastic modulus and yield strength of human trabecular bone. In the current study, plate-rod (PR) finite element (FE) models were constructed completely based on ITS-identified individual trabecular plates and rods. We hypothesized that PR FE can accurately and efficiently predict elastic modulus and yield strength of human trabecular bone. Human trabecular bone cores from proximal tibia (PT), femoral neck (FN) and greater trochanter (GT) were scanned by µCT. Specimen-specific ITS-based PR FE models were generated for each µCT image and corresponding voxel-based FE models were also generated in comparison. Both types of specimen-specific models were subjected to nonlinear FE analysis to predict the apparent elastic modulus and yield strength using the same trabecular bone tissue properties. Then, mechanical tests were performed to experimentally measure the apparent modulus and yield strength. Strong linear correlations for both elastic modulus (r(2) = 0.97) and yield strength (r(2) = 0.96) were found between the PR FE model predictions and experimental measures, suggesting that trabecular plate and rod morphology adequately captures three-dimensional (3D) microarchitecture of human trabecular bone. In addition, the PR FE model predictions in both elastic modulus and yield strength were highly correlated with the voxel-based FE models (r(2) = 0.99, r(2) = 0.98, respectively), resulted from the original 3D images without the PR segmentation. In conclusion, the ITS-based PR models predicted accurately both elastic modulus and yield strength determined experimentally across three distinct anatomic sites. Trabecular plates and rods accurately determine elastic modulus and yield strength of human trabecular bone.

The quartic piecewise-linear criterion for the multiaxial yield behavior of human trabecular bone.

Prior multiaxial strength studies on trabecular bone have either not addressed large variations in bone volume fraction and microarchitecture, or have not addressed the full range of multiaxial stress states. Addressing these limitations, we utilized micro-computed tomography (lCT) based nonlinear finite element analysis to investigate the complete 3D multiaxial failure behavior of ten specimens (5mm cube) of human trabecular bone, taken from three anatomic sites and spanning a wide range of bone volume fraction (0.09­0.36),mechanical anisotropy (range of E3/E1»3.0­12.0), and microarchitecture. We found that most of the observed variation in multiaxial strength behavior could be accounted for by normalizing the multiaxial strength by specimen-specific values of uniaxial strength (tension,compression in the longitudinal and transverse directions). Scatter between specimens was reduced further when the normalized multiaxial strength was described in strain space.The resulting multiaxial failure envelope in this normalized-strain space had a rectangular boxlike shape for normal­normal loading and either a rhomboidal box like shape or a triangular shape for normal-shear loading, depending on the loading direction. The finite element data were well described by a single quartic yield criterion in the 6D normalized strain space combined with a piecewise linear yield criterion in two planes for normalshear loading (mean error SD: 4.660.8% for the finite element data versus the criterion).This multiaxial yield criterion in normalized-strain space can be used to describe the complete 3D multiaxial failure behavior of human trabecular bone across a wide range of bone volume fraction, mechanical anisotropy, and microarchitecture.

Mechanical instability of normal and aneurysmal arteries.

Tortuous arteries associated with aneurysms have been observed in aged patients with atherosclerosis and hypertension. However, the underlying mechanism is poorly understood. The objective of this study was to determine the effect of aneurysms on arterial buckling instability and the effect of buckling on aneurysm wall stress. We investigated the mechanical buckling and post-buckling behavior of normal and aneurysmal carotid arteries and aorta's using computational simulations and experimental measurements to elucidate the interrelationship between artery buckling and aneurysms. Buckling tests were done in porcine carotid arteries with small aneurysms created using elastase treatment. Parametric studies were done for model aneurysms with orthotropic nonlinear elastic walls using finite element simulations. Our results demonstrated that arteries buckled at a critical buckling pressure and the post-buckling deflection increased nonlinearly with increasing pressure. The presence of an aneurysm can reduce the critical buckling pressure of arteries, although the effect depends on the aneurysm's dimensions. Buckled aneurysms demonstrated a higher peak wall stress compared to unbuckled aneurysms under the same lumen pressure. We conclude that aneurysmal arteries are vulnerable to mechanical buckling and mechanical buckling could lead to high stresses in the aneurysm wall. Buckling could be a possible mechanism for the development of tortuous aneurysmal arteries such as in the Loeys-Dietz syndrome.

Biaxial normal strength behavior in the axial-transverse plane for human trabecular bone--effects of bone volume fraction, microarchitecture, and anisotropy.

The biaxial failure behavior of the human trabecular bone, which has potential relevance both for fall and gait loading conditions, is not well understood, particularly for low-density bone, which can display considerable mechanical anisotropy. Addressing this issue, we investigated the biaxial normal strength behavior and the underlying failure mechanisms for human trabecular bone displaying a wide range of bone volume fraction (0.06-0.34) and elastic anisotropy. Micro-computed tomography (CT)-based nonlinear finite element analysis was used to simulate biaxial failure in 15 specimens (5 mm cubes), spanning the complete biaxial normal stress failure space in the axial-transverse plane. The specimens, treated as approximately transversely isotropic, were loaded in the principal material orientation. We found that the biaxial stress yield surface was well characterized by the superposition of two ellipses--one each for yield failure in the longitudinal and transverse loading directions--and the size, shape, and orientation of which depended on bone volume fraction and elastic anisotropy. However, when normalized by the uniaxial tensile and compressive strengths in the longitudinal and transverse directions, all of which depended on bone volume fraction, microarchitecture, and mechanical anisotropy, the resulting normalized biaxial strength behavior was well described by a single pair of (longitudinal and transverse) ellipses, with little interspecimen variation. Taken together, these results indicate that the role of bone volume fraction, microarchitecture, and mechanical anisotropy is mostly accounted for in determining the uniaxial strength behavior and the effect of these parameters on the axial-transverse biaxial normal strength behavior per se is minor.

Shear strength behavior of human trabecular bone.

The shear strength of human trabecular bone may influence overall bone strength under fall loading conditions and failure at bone-implant interfaces. Here, we sought to compare shear and compressive yield strengths of human trabecular bone and elucidate the underlying failure mechanisms. We analyzed 54 specimens (5-mm cubes), all aligned with the main trabecular orientation and spanning four anatomic sites, 44 different cadavers, and a wide range of bone volume fraction (0.06-0.38). Micro-CT-based non-linear finite element analysis was used to assess the compressive and shear strengths and the spatial distribution of yielded tissue; the tissue-level constitutive model allowed for kinematic non-linearity and yielding with strength asymmetry. We found that the computed values of both the shear and compressive strengths depended on bone volume fraction via power law relations having an exponent of 1.7 (R(2)=0.95 shear; R(2)=0.97 compression). The ratio of shear to compressive strengths (mean±SD, 0.44±0.16) did not depend on bone volume fraction (p=0.24) but did depend on microarchitecture, most notably the intra-trabecular standard deviation in trabecular spacing (R(2)=0.23, p<0.005). For shear, the main tissue-level failure mode was tensile yield of the obliquely oriented trabeculae. By contrast, for compression, specimens having low bone volume fraction failed primarily by large-deformation-related tensile yield of horizontal trabeculae and those having high bone volume failed primarily by compressive yield of vertical trabeculae. We conclude that human trabecular bone is generally much weaker in shear than compression at the apparent level, reflecting different failure mechanisms at the tissue level.

Prediction of new clinical vertebral fractures in elderly men using finite element analysis of CT scans.

Vertebral strength, as estimated by finite element analysis of computed tomography (CT) scans, has not yet been compared against areal bone mineral density (BMD) by dual-energy X-ray absorptiometry (DXA) for prospectively assessing the risk of new clinical vertebral fractures. To do so, we conducted a case-cohort analysis of 306 men aged 65 years and older, which included 63 men who developed new clinically-identified vertebral fractures and 243 men who did not, all observed over an average of 6.5 years. Nonlinear finite element analysis was performed on the baseline CT scans, blinded to fracture status, to estimate L1 vertebral compressive strength and a load-to-strength ratio. Volumetric BMD by quantitative CT and areal BMD by DXA were also evaluated. We found that, for the risk of new clinical vertebral fracture, the age-adjusted hazard ratio per standard deviation change for areal BMD (3.2; 95% confidence interval [CI], 2.0-5.2) was significantly lower (p < 0.005) than for strength (7.2; 95% CI, 3.6-14.1), numerically lower than for volumetric BMD (5.7; 95% CI, 3.1-10.3), and similar for the load-to-strength ratio (3.0; 95% CI, 2.1-4.3). After also adjusting for race, body mass index (BMI), clinical center, and areal BMD, all these hazard ratios remained highly statistically significant, particularly those for strength (8.5; 95% CI, 3.6-20.1) and volumetric BMD (9.4; 95% CI, 4.1-21.6). The area-under-the-curve for areal BMD (AUC = 0.76) was significantly lower than for strength (AUC = 0.83, p = 0.02), volumetric BMD (AUC = 0.82, p = 0.05), and the load-to-strength ratio (AUC = 0.82, p = 0.05). We conclude that, compared to areal BMD by DXA, vertebral compressive strength and volumetric BMD consistently improved vertebral fracture risk assessment in this cohort of elderly men.