On the biomechanical significance of inter-lamellar interfaces in the intervertebral disc

The intervertebral disc (IVD) is a unique soft tissue structure which provides structural support and flexibility in the axial skeleton of vertebrates. From a structural perspective, the disc behaves somewhat like a thick walled pressure vessel, where the walls are comprised of a series of composite annular rings (lamellae). However, a prior study (Marchand and Ahmed, 1990) found a high proportion of circumferentially discontinuous lamellae in human lumbar IVDs. The presence of these discontinuities raises important structural questions, because discontinuous lamellae cannot withstand high nucleus pressures via the generation of circumferential (hoop) stress. A possible alternative mechanism may be that inter-lamellar cohesion allows shear stress transfer between adjacent annular layers. The aim of the present study was therefore to investigate the importance of inter-lamellar shear resistance in the intervertebral disc. This work found that inter-lamellar shear resistance has a strong influence on the compressive stiffness of the intervertebral disc, with a change in interface condition from tied (no slip) to frictionless (no shear resistance) reducing disc compressive stiffness by 40%. However, it appears that substantial inter-lamellar shear resistance is present in the bovine tail disc. Decreases in inter-lamellar shear resistance due to degradation of bridging collagenous or elastic fibre structures could therefore be an important part of the process of disc degeneration.

Cardiovascular diseases and vulnerable plaques: data, modeling, predictions and clinical applications PREFACE

Introduction Two symposia on “cardiovascular diseases and vulnerable plaques” Cardiovascular disease (CVD) is the leading cause of death worldwide. Huge effort has been made in many disciplines including medical imaging, computational modeling, bio- mechanics, bioengineering, medical devices, animal and clinical studies, population studies as well as genomic, molecular, cellular and organ-level studies seeking improved methods for early detection, diagnosis, prevention and treatment of these diseases [1-14]. However, the mechanisms governing the initiation, progression and the occurrence of final acute clinical CVD events are still poorly understood. A large number of victims of these dis- eases who are apparently healthy die suddenly without prior symptoms. Available screening and diagnostic methods are insufficient to identify the victims before the event occurs [8,9]. Most cardiovascular diseases are associated with vulnerable plaques. A grand challenge here is to develop new imaging techniques, predictive methods and patient screening tools to identify vulnerable plaques and patients who are more vulnerable to plaque rupture and associated clinical events such as stroke and heart attack, and recommend proper treatment plans to prevent those clinical events from happening. Articles in this special issue came from two symposia held recently focusing on “Cardio-vascular Diseases and Vulnerable Plaques: Data, Modeling, Predictions and Clinical Applications.” One was held at Worcester Polytechnic Institute (WPI), Worcester, MA, USA, July 13-14, 2014, right after the 7th World Congress of Biomechanics. This symposium was endorsed by the World Council of Biomechanics, and partially supported by a grant from NIH-National Institute of Biomedical Image and Bioengineering. The other was held at Southeast University (SEU), Nanjing, China, April 18-20, 2014.

A 3D numerical study of the collateral capacity of the circle of Willis with anatomical variation in the posterior circulation

Background The Circle of Willis (CoW) is the most important collateral pathway of the cerebral artery. The present study aims to investigate the collateral capacity of CoW with anatomical variation when unilateral internalcarotid artery (ICA) is occluded. Methods Basing on MRI data, we have reconstructed eight 3D models with variations in the posterior circulation of the CoW and set four different degrees of stenosis in the right ICA, namely 24%, 43%, 64% and 79%, respectively. Finally, a total of 40 models are performed with computational fluid dynamics simulations. All of the simulations share the same boundary condition with static pressure and the volume flow rate (VFR) are obtained to evaluate their collateral capacity. Results As for the middle cerebral artery (MCA) and the anterior cerebral artery (ACA), the transitional-type model possesses the best collateral capacity. But for the posterior cerebral artery (PCA), unilateral stenosis of ICA has the weakest influence on the unilateral posterior communicating artery (PCoA) absent model. We also find that the full fetal-type posterior circle of Willis is an utmost dangerous variation which must be paid more attention. Conclusion The results demonstrate that different models have different collateral capacities in coping stenosis of unilateral ICA and these differences can be reflected by different outlets. The study could be used as a reference for neurosurgeon in choosing the best treatment strategy.

Flow pattern analysis in a highly stenotic patient-specific carotid bifurcation model using a turbulence model

The aim of this study is to investigate the blood flow pattern in carotid bifurcation with a high degree of luminal stenosis, combining in vivo magnetic resonance imaging (MRI) and computational fluid dynamics (CFD). A newly developed two-equation transitional model was employed to evaluate wall shear stress (WSS) distribution and pressure drop across the stenosis, which are closely related to plaque vulnerability. A patient with an 80% left carotid stenosis was imaged using high resolution MRI, from which a patient-specific geometry was reconstructed and flow boundary conditions were acquired for CFD simulation. A transitional model was implemented to investigate the flow velocity and WSS distribution in the patient-specific model. The peak time-averaged WSS value of approximately 73Pa was predicted by the transitional flow model, and the regions of high WSS occurred at the throat of the stenosis. High oscillatory shear index values up to 0.50 were present in a helical flow pattern from the outer wall of the internal carotid artery immediately after the throat. This study shows the potential suitability of a transitional turbulent flow model in capturing the flow phenomena in severely stenosed carotid arteries using patient-specific MRI data and provides the basis for further investigation of the links between haemodynamic variables and plaque vulnerability. It may be useful in the future for risk assessment of patients with carotid disease.

MRI-based strain and strain rate analysis of left ventricle: A modified hierarchical transformation model

Background Different from other indicators of cardiac function, such as ejection fraction and transmitral early diastolic velocity, myocardial strain is promising to capture subtle alterations that result from early diseases of the myocardium. In order to extract the left ventricle (LV) myocardial strain and strain rate from cardiac cine-MRI, a modified hierarchical transformation model was proposed. Methods A hierarchical transformation model including the global and local LV deformations was employed to analyze the strain and strain rate of the left ventricle by cine-MRI image registration. The endocardial and epicardial contour information was introduced to enhance the registration accuracy by combining the original hierarchical algorithm with an Iterative Closest Points using Invariant Features algorithm. The hierarchical model was validated by a normal volunteer first and then applied to two clinical cases (i.e., the normal volunteer and a diabetic patient) to evaluate their respective function. Results Based on the two clinical cases, by comparing the displacement fields of two selected landmarks in the normal volunteer, the proposed method showed a better performance than the original or unmodified model. Meanwhile, the comparison of the radial strain between the volunteer and patient demonstrated their apparent functional difference. Conclusions The present method could be used to estimate the LV myocardial strain and strain rate during a cardiac cycle and thus to quantify the analysis of the LV motion function.

A follow up MRI-based geometry and computational fluid dynamics study of carotid bifurcation

Cardiovascular disease is the leading causes of death in the developed world. Wall shear stress (WSS) is associated with the initiation and progression of atherogenesis. This study combined the recent advances in MR imaging and computational fluid dynamics (CFD) and evaluated the patient-specific carotid bifurcation. The patient was followed up for 3 years. The geometry changes (tortuosity, curvature, ICA/CCA area ratios, central to the cross-sectional curvature, maximum stenosis) and the CFD factors (Velocity distribute, Wall Shear Stress (WSS) and Oscillatory Shear Index (OSI)) were compared at different time points.The carotid stenosis was a slight increase in the central to the cross-sectional curvature, and it was minor and variable curvature changes for carotid centerline. The OSI distribution presents ahigh-values in the same region where carotid stenosis and normal border, indicating complex flow and recirculation.The significant geometric changes observed during the follow-up may also cause significant changes in bifurcation hemodynamics.

Effect of inflow and outflow angles on the computational hemodynamics in abdominal aortic aneurysm

To help with the clinical screening and diagnosis of abdominal aortic aneurysm (AAA), we evaluated the effect of inflow angle (IA) and outflow bifurcation angle (BA) on the distribution of blood flow and wall shear stress (WSS) in an idealized AAA model. A 2D incompressible Newtonian flow is assumed and the computational simulation is performed using finite volume method. The results showed that the largest WSS often located at the proximal and the distal end of the AAA. An increase in IA resulted in an increase in maximum WSS. We also found that WSS was maximal when BA was 90°. IA and BA are two important geometrical factors, they may help with AAA risk assessment along with the commonly used AAA diameter.

Mechanical analysis of foot plantar fascia during normal walking condition

Foot plantar fascia is an important foot tissue in stabilizing the longitudinal arch of human foot. Direct measurement to monitor the mechanical situation of plantar fascia at human locomotion is difficult. The purpose of this study was to construct a three-dimensional finite element model of the foot to calculate the internal stress/strain value of plantar fascia during different stage of gait. The simulated stress distribution of plantar fascia was the lowest at heel-strike, which concentrated on the medial side of calcaneal tubercle. The peak stress of plantar fascia was appeared at push-off, and the value is more than 5 times of the heel-strike position. Current FE model was able to explore the plantar fascia tension trend at the main sub-phases of foot. More detailed fascia model and intrinsic muscle forces could be developed in the further study.

Image-based midsole insert design and the material effects on heel plantar pressure distribution during simulated walking loads

The primary objective of this paper is to study the use of medical image-based finite element (FE) modelling in subjectspecific midsole design and optimisation for heel pressure reduction using a midsole plug under the calcaneus area (UCA). Plugs with different relative dimensions to the size of the calcaneus of the subject have been incorporated in the heel region of the midsole. The FE foot model was validated by comparing the numerically predicted plantar pressure with biomechanical tests conducted on the same subject. For each UCA midsole plug design, the effect of material properties and plug thicknesses on the plantar pressure distribution and peak pressure level during the heel strike phase of normal walking was systematically studied. The results showed that the UCA midsole insert could effectively modify the pressure distribution, and its effect is directly associated with the ratio of the plug dimension to the size of the calcaneus bone of the subject. A medium hardness plug with a size of 95% of the calcaneus has achieved the best performance for relieving the peak pressure in comparison with the pressure level for a solid midsole without a plug, whereas a smaller plug with a size of 65% of the calcaneus insert with a very soft material showed minimum beneficial effect for the pressure relief.

Stress analysis of carotid atheroma in transient ischemic attack patients: Evidence for extreme stress-induced plaque rupture

Plaque rupture has been considered to be the result of its structural failure. The aim of this study is to suggest a possible link between higher stresses and rupture sites observed from in vivo magnetic resonance imaging (MRI) of transient ischemic attack (TIA) patients, by using stress analysis methods. Three patients, who had recently suffered a TIA, underwent in vivo multi-spectral MR imaging. Based on plaque geometries reconstructed from the post-rupture status, six pre-rupture plaque models were generated for each patient dataset with different reconstructions of rupture sites to bridge the gap of fibrous cap from original MRI images. Stress analysis by fluid structure interaction simulation was performed on the models, followed by analysis of local stress concentration distribution and plaque rupture sites. Furthermore, the sensitivity of stress analysis to the pre-rupture plaque geometry reconstruction was examined. Local stress concentrations were found to be located at the plaque rupture sites for the three subjects studied. In the total of 18 models created, the locations of the stress concentration regions were similar in 17 models in which rupture sites were always associated with high stresses. The local stress concentration region moved from circumferential center to the shoulder region (slightly away from the rupture site) for a case with a thick fibrous cap. Plaque wall stress level in the rupture locations was found to be much higher than the value in non-rupture locations. The good correlation between local stress concentrations and plaque rupture sites, and generally higher plaque wall stress level in rupture locations in the subjects studied could provide indirect evidence for the extreme stress-induced plaque rupture hypothesis. Local stress concentration in the plaque region could be one of the factors contributing to plaque rupture.

Arterial luminal curvature and fibrous-cap thickness affect critical stress conditions within atherosclerotic plaque: An in vivo MRI-based 2D finite-element study

High mechanical stress in atherosclerotic plaques at vulnerable sites, called critical stress, contributes to plaque rupture. The site of minimum fibrous cap (FC) thickness (FCMIN) and plaque shoulder are well-documented vulnerable sites. The inherent weakness of the FC material at the thinnest point increases the stress, making it vulnerable, and it is the big curvature of the lumen contour over FC which may result in increased plaque stress. We aimed to assess critical stresses at FCMIN and the maximum lumen curvature over FC (LCMAX) and quantify the difference to see which vulnerable site had the highest critical stress and was, therefore, at highest risk of rupture. One hundred patients underwent high resolution carotid magnetic resonance (MR) imaging. We used 352 MR slices with delineated atherosclerotic components for the simulation study. Stresses at all the integral nodes along the lumen surface were calculated using the finite-element method. FCMIN and LCMAX were identified, and critical stresses at these sites were assessed and compared. Critical stress at FC MIN was significantly lower than that at LCMAX (median: 121.55 kPa; inter quartile range (IQR) = [60.70-180.32] kPa vs. 150.80 kPa; IQR = [91.39-235.75] kPa, p < 0.0001). If critical stress at FCMIN was only used, then the stress condition of 238 of 352 MR slices would be underestimated, while if the critical stress at LCMAX only was used, then 112 out of 352 would be underestimated. Stress analysis at FCMIN and LCMAX should be used for a refined mechanical risk assessment of atherosclerotic plaques, since material failure at either site may result in rupture.

Stress analysis of carotid plaque based on in vivo MRI of acute symptomatic and asymptomatic patients

Stress analysis within carotid plaques based on in vivo MR imaging has shown to be useful for the identification of vulnerable atheroma. This study is to investigate whether magnetic resonance imaging (MRI) based-biomechanical stress analysis of carotid plaques can differentiate acute symptomatic and asymptomatic patients. 54 asymptomatic and 45 acute symptomatic patients underwent in vivo multi-contrast MRI of the carotid arteries. Plaque geometry used for finite element analysis was derived from in vivo MR images at the site of maximum and minimum plaque burden. In total 198 slices were used for the computational simulations. A pre shrink technique was used to refine the simulation. Maximum principle stress at the vulnerable plaque sites (i.e. critical stress) was extracted for the selected slices and a comparison was performed between the two groups. Critical stress at the site of maximum plaque burden is significantly higher in acute symptomatic patients as compared to asymptomatic patients [median: 198.0kPa (inter quartile range (IQR) = (119.8 - 359.0) vs. 138.4kPa (83.8, 242.6), p=0.04]. No significant difference was found at the minimum plaque burden site between the two groups [196.7kPa (133.3- 282.7) vs. 182.4kPa (117.2 - 310. 6), p=0.82). Stress analysis at the site of maximal plaque burden can be effectively used for differentiating acute symptomatic carotid plaques from asymptomatic plaques. This maybe potentially used for development of biomechanical risk stratification criteria based on plaque burden in future studies.

Curvedness study on atherosclerosis plaques and its implications to plaque stress

Atherosclerosis plaque rupture has been considered to be a mechanical failure of the thin fibrous cap, resulted from extreme plaque stress. Plaque stress was affected by many factors from morphological features to biological abnormalities. In this study, geometrical factors (curvedness, fibrous cap thickness) were studied on assessing plaque vulnerability in comparison with stress analysis results obtained by fluid structure interaction from 20 human carotid atherosclerosis plaques. The results show that plaque surface curvedness could contribute to extreme stress level, especially in plaque shoulder region. General plaque stress distribution could be predicted by fibrous cap thickness and curvedness with multi-regression model. With more features included in the regression model, plaque stress could be easily calculated and used to assess plaque vulnerability.

Study on tracheal collapsibility, compliance, and stress by considering nonlinear mechanical property of cartilage

Tracheal cartilage has been widely regarded as a linear elastic material either in experimental studies or in analytic and numerical models. However, it has been recently demonstrated that, like other fiber-oriented biological tissues, tracheal cartilage is a nonlinear material, which displays higher strength in compression than in extension. Considering the nonlinearity requires a more complex theoretical frame work and costs more to simulate. This study aims to quantify the deviation due to the simplified treatment of the tracheal cartilage as a linear material. It also evaluates the improved accuracy gained by considering the nonlinearity. Pig tracheal rings were used to exam the mechanical properties of cartilage and muscular membrane. By taking into account the asymmetric shape of tracheal cartilage, the collapse behavior of complete rings was simulated, and the compliance of airway and stress in the muscular membrane were discussed. The results obtained were compared with those assuming linear mechanical properties. The following results were found: (1) Models based on both types of material properties give a small difference in representing collapse behavior; (2) regarding compliance, the relative difference is big, ranging from 10 to 40% under negative pressure conditions; and (3) the difference in determining stress in the muscular membrane is small too: <5%. In conclusion, treating tracheal cartilage as a linear material will not cause big deviations in representing the collapse behavior, and mechanical stress in the muscular part, but it will induce a big deviation in predicting the compliance, particularly when the transmural pressure is lower than -0.5 kPa. The results obtained in this study may be useful in both understanding the collapse behavior of trachea and in evaluating the error induced by the simplification of treating the tracheal cartilage as a linear elastic material.

Does PGA external stenting reduce compliance mismatch in venous grafts?

Background: Autogenous vein grafting is widely used in regular bypassing procedures. Due to its mismatch with the host artery in both mechanical property and geometry, the graft often over expands under high arterial blood pressure and forms a step-depth where eddy flow develops, thus causing restenosis, fibrous graft wall, etc. External stents, such as sheaths being used to cuff the graft, have been introduced to eliminate these mismatches and increase the patency. Although histological and immunochemical studies have shown some positive effects of the external stent, the mechanical mismatch under the protection of an external stent remains poorly analyzed. Methods: In this study, the jugular veins taken from hypercholesterolemic rabbits were transplanted into the carotid arteries, and non-woven polyglycolic acid (PGA) fabric was used to fabricate the external stents to study the effect of the biodegradable external stent. Eight weeks after the operation, the grafts were harvested to perform mechanical tests and histological examinations. An arc tangent function was suggested to describe the relationship between pressure and cross-sectional area to analyse the compliance of the graft. Results: The results from the mechanical tests indicated that grafts either with or without external stents displayed large compliance in the low-pressure range and were almost inextensible in the high-pressure range. This was very different from the behavior of the arteries or veins in vivo. The data from histological tests showed that, with external stents, collagen fibers were more compact, whilst those in the graft without protection were looser and thicker. No elastic fiber was found in either kind of grafts. Furthermore, grafts without protection were over-expanded which resulted in much bigger cross-sectional areas. Conclusion: The PGA external extent contributes little to the reduction of the mechanical mismatch between the graft and its host artery while remodeling develops. For the geometric mismatch, it reduces the cross-section area, therefore matching with the host artery much better. Although there are some positive effects, conclusively the PGA is not an ideal material for external stent.