Utility of magnetic resonance imaging-based finite element analysis for the biomechanical stress analysis of hemorrhagic and non-hemorrhagic carotid plaques

Background: Biomechanical stress analysis has been used for plaque vulnerability assessment. The presence of plaque hemorrhage (PH) is a feature of plaque vulnerability and is associated with thromboembolic ischemic events. The purpose of the present study was to use finite element analysis (FEA) to compare the stress profiles of hemorrhagic and non-hemorrhagic profiles. Methods and Results: Forty-five consecutive patients who had suffered a cerebrovascular ischemic event with an underlying carotid artery disease underwent high-resolution magnetic resonance imaging (MRI) of their symptomatic carotid artery in a 1.5-T MRI system. Axial images were manually segmented for various plaque components and used for FEA. Maximum critical stress (M-CstressSL) for each slice was determined. Within a plaque, the maximum M-CstressSL for each slice of a plaque was selected to represent the maximum critical stress of that plaque (M-CstressPL) and used to compare hemorrhagic and non-hemorrhagic plaques. A total of 62% of plaques had hemorrhage. It was observed that plaques with hemorrhage had significantly higher stress (M-CstressPL) than plaques without PH (median [interquartile range]: 315 kPa [247-434] vs. 200 kPa [171-282], P=0.003). Conclusions: Hemorrhagic plaques have higher biomechanical stresses than non-hemorrhagic plaques. MRI-based FEA seems to have the potential to assess plaque vulnerability.

Impact of calcification and intraluminal thrombus on the computed wall stresses of abdominal aortic aneurysm

Background: Increased biomechanical stresses within the abdominal aortic aneurysm (AAA) wall contribute to its rupture. Calcification and intraluminal thrombus can be commonly found in AAAs, but the relationship between calcification/intraluminal thrombus and AAA wall stress is not completely described. Methods: Patient-specific three-dimensional AAA geometries were reconstructed from computed tomographic images of 20 patients. Structural analysis was performed to calculate the wall stresses of the 20 AAA models and their altered models when calcification or intraluminal thrombus was not considered. A nonlinear large-strain finite element method was used to compute the wall stress distribution. The relationships between wall stresses and volumes of calcification and intraluminal thrombus were sought. Results: Maximum stress was not correlated with the percentage of calcification, and was negatively correlated with the percentage of intraluminal thrombus (r = -0.56; P = .011). Exclusion of calcification from analysis led to a significant decrease in maximum stress by a median of 14% (range, 2%-27%; P < .01). When intraluminal thrombus was eliminated, maximum stress increased significantly by a median of 24% (range, 5%-43%; P < .01). Conclusion: The presence of calcification increases AAA peak wall stress, suggesting that calcification decrease the biomechanical stability of AAA. In contrast, intraluminal thrombus reduces the maximum stress in AAA. Calcification and intraluminal thrombus should both be considered in the evaluation of wall stress for risk assessment of AAA rupture.

Normalized wall index specific and MRI-based stress analysis of atherosclerotic carotid plaques: A study comparing acutely symptomatic and asymptomatic patients

Background: Biomechanical stresses play an important role in determining plaque stability. Quantification of these simulated stresses can be potentially used to assess plaque vulnerability and differentiate different patient groups. Methods and Results: 54 asymptomatic and 45 acutely symptomatic patients underwent in vivo multicontrast magnetic resonance imaging (MRI) of the carotid arteries. Plaque geometry used for finite element analysis was derived from in vivo MRI at the sites 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 (ie, critical stress) was extracted for the selected slices and a comparison was performed between the 2 groups. Critical stress in the slice with maximum plaque burden is significantly higher in acutely symptomatic patients as compared to asymptomatic patients (median, inter quartile range: 198.0 kPa (119.8-359.0 kPa) vs 138.4 kPa (83.8-242.6 kPa), P=0.04). No significant difference was found in the slice with minimum plaque burden between the 2 groups (196.7 kPa (133.3-282.7 kPa) vs 182.4 kPa (117.2-310.6 kPa), P=0.82). Conclusions: Acutely symptomatic carotid plaques have significantly high biomechanical stresses than asymptomatic plaques. This might be potentially useful for establishing a biomechanical risk stratification criteria based on plaque burden in future studies.

Identifying vulnerable carotid plaques in vivo using high resolution magnetic resonance imaging-based finite element analysis

Object. Individuals with carotid atherosclerosis develop symptoms following rupture of vulnerable plaques. Biomechanical stresses within this plaque may increase vulnerability to rupture. In this report the authors describe the use of in vivo carotid plaque imaging and computational mechanics to document the magnitude and distribution of intrinsic plaque stresses. Methods. Ten (five symptomatic and five asymptomatic) individuals underwent plaque characterization magnetic resonance (MR) imaging. Plaque geometry and composition were determined by multisequence review. Intrinsic plaque stress profiles were generated from 3D meshes by using finite element computational analysis. Differences in principal (shear) stress between normal and diseased sections of the carotid artery and between symptomatic and asymptomatic plaques were noted. Results. There was a significant difference in peak principal stress between diseased and nondiseased segments of the artery (mean difference 537.65 kPa, p < 0.05). Symptomatic plaques had higher mean stresses than asymptomatic plaques (627.6 kPa compared with 370.2 kPa, p = 0.05), which were independent of luminal stenosis and plaque composition. Conclusions. Significant differences in plaque stress exist between plaques from symptomatic individuals and those from asymptomatic individuals. The MR imaging-based computational analysis may therefore be a useful aid to identification of vulnerable plaques in vivo.

Association between biomechanical structural stresses of atherosclerotic carotid plaques and subsequent ischaemic cerebrovascular events - A longitudinal in vivo magnetic resonance imaging-based finite element study

Background: High-resolution magnetic resonance (MR) imaging has been used for MR imaging-based structural stress analysis of atherosclerotic plaques. The biomechanical stress profile of stable plaques has been observed to differ from that of unstable plaques; however, the role that structural stresses play in determining plaque vulnerability remains speculative. Methods: A total of 61 patients with previous history of symptomatic carotid artery disease underwent carotid plaque MR imaging. Plaque components of the index artery such as fibrous tissue, lipid content and plaque haemorrhage (PH) were delineated and used for finite element analysis-based maximum structural stress (M-C Stress) quantification. These patients were followed up for 2 years. The clinical end point was occurrence of an ischaemic cerebrovascular event. The association of the time to the clinical end point with plaque morphology and M-C Stress was analysed. Results: During a median follow-up duration of 514 days, 20% of patients (n=12) experienced an ischaemic event in the territory of the index carotid artery. Cox regression analysis indicated that M-C Stress (hazard ratio (HR): 12.98 (95% confidence interval (CI): 1.32-26.67, pZ0.02), fibrous cap (FC) disruption (HR: 7.39 (95% CI: 1.61e33.82), p Z 0.009) and PH (HR: 5.85 (95% CI: 1.27e26.77), p Z 0.02) are associated with the development of subsequent cerebrovascular events. Plaques associated with future events had higher M-C Stress than those which had remained asymptomatic (median (interquartile range, IQR): 330 kPa (229e494) vs. 254 kPa (166-290), p Z0.04). Conclusions: High biomechanical structural stresses, in addition to FC rupture and PH, are associated with subsequent cerebrovascular events.

Finite element analysis of vulnerable atherosclerotic plaques: A comparison of mechanical stresses within carotid plaques of acute and recently symptomatic patients with carotid artery disease

Objectives: There is considerable evidence that patients with carotid artery stenosis treated immediately after the ischaemic cerebrovascular event have a better clinical outcome than those who have delayed treatment. Biomechanical assessment of carotid plaques using high-resolution MRI can help examine the relationship between the timing of carotid plaque symptomology and maximum simulated plaque stress concentration. Methods: Fifty patients underwent high-resolution multisequence in vivo MRI of their carotid arteries. Patients with acute symptoms (n=25) underwent MRI within 72 h of the onset of ischaemic cerebrovascular symptoms, whereas recently symptomatic patients (n=25) underwent MRI from 2 to 6 weeks after the onset of symptoms. Stress analysis was performed based on the geometry derived from in vivo MRI of the symptomatic carotid artery at the point of maximum stenosis. The peak stresses within the plaques of the two groups were compared. Results: Patient demographics were comparable for both groups. All the patients in the recently symptomatic group had severe carotid stenosis in contrast to patients with acute symptoms who had predominantly mild to moderate carotid stenosis. The simulated maximum stresses in patients with acute symptoms was significantly higher than in recently symptomatic patients (median (IQR): 313310 4 dynes/cm 2 (295 to 382) vs 2523104 dynes/cm 2 (236 to 311), p=0.02). Conclusions: Patients have extremely unstable, high-risk plaques, with high stresses, immediately after an acute cerebrovascular event, even at lower degrees of carotid stenoses. Biomechanical stress analysis may help us refine our risk-stratification criteria for the management of patients with carotid artery disease in future.

Evaluation of cement stresses in finite element analyses of cemented orthopaedic implants

Stress analysis of the cement fixation of orthopaedic implants to bone is frequently? carried out using finite element analysis. However, the stress distribution in the cement laver is usually intricate, and it is difficult to report it in a way that facilitates comparison of implants for pre-clinical testing. To study this problem, and make recommendations for stress reporting, a finite element analysis of a hip prosthesis implanted into a synthetic composite femur is developed. Three cases are analyzed: a fully bonded implant, a debonded implant, and a debonded implant where the cement is removed distal to the stein tip. In addition to peak stresses, and contour and vector plots, a stressed volume and probability-of-failure analysis is reported. It is predicted that the peak stress is highest for the debonded stem, and that removal of the distal cement more than halves this peak stress. This would suggest that omission of the distal cement is good for polished prostheses (as practiced for the Exeter design). However; if the percentage of cement stressed above a certain threshold (say 3 MPa) is considered, then the removal of distal cement is shown to be disadvantageous because a higher volume of cement is stressed to above the threshold. Vector plots clearly demonstrate the different load transfer for bonded and debonded prostheses: A bonded stein generates maximum tensile stresses in the longitudinal direction, whereas a debonded stem generates most tensile stresses in the hoop direction, except near the tip where tensile longitudinal stresses occur due to subsidence of the stein. Removal of the cement distal to the tip allows greater subsidence but alleviates these large stresses at the tip, albeit at the expense of increased hoop stresses throughout the mantle. It is concluded that a thorough analysis of cemented implants should not report peak stress, which can be misleading, but rather stressed volume, and that vector plots should be reported if a precise analysis of the load transfer mechanism is required.

Tilted and Short Implants Supporting Fixed Prosthesis in an Atrophic Maxilla: A 3D-FEA Biomechanical Evaluation

Purpose: This study compared the biomechanical behavior of tilted long implant and vertical short implants to support fixed prosthesis in an atrophic maxilla. Materials and Methods: The maxilla model was built based on a tomographic image of the patient. Implant models were based on micro-computer tomography imaging of implants. The different configurations considered were M4S, four vertical anterior implants; M4T, two mesial vertical implants and two distal tilted (45°) implants in the anterior region of the maxilla; and M6S, four vertical anterior implants and two vertical posterior implants. Numerical simulation was carried out under bilateral 150N loads applied in the cantilever region in axial (L1) and oblique (45°) (L2) direction. Bone was analyzed using the maximum and minimum principal stress (σmax and σmin), and von Mises stress (σvM) assessments. Implants were analyzed using the σvM. Results: The higher σmax was observed at: M4T, followed by M6S/L1, M6S/L2, M4S/L2, and M4S/L1 and the higher σvM: M4T/L1, M4T/L2 and M4S/L2, M6S/L2, M4S/L1, and M6S/L1. Conclusions: The presence of distal tilted (all-on-four) and distal short implants (all-on-six) resulted in higher stresses in both situations in the maxillary bone in comparison to the presence of vertical implants (all-on-four). © 2013 Wiley Periodicals, Inc.

The influence of bone quality on the biomechanical behavior of full-arch implant-supported fixed prostheses

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Analysis of the biomechanical behavior of short implants: the photo-elasticity method

The aim of this study was to analyze the stress distribution of short implants supporting single unit or splinted crowns by the photo-elasticity method. Four photo-elastic models were produced: A (3.75×7mm); B (3.75×7mm, 3.75×7mm and 3.75×7mm); C (3.75×10mm, 3.75×7mm and 3.75×7mm); D (3.75×13mm, 3.75×7mm and 3.75×7mm). The prostheses were made with Ni-Cr alloy. A load of 100N in the axial and oblique directions was applied, totaling 380 applications, individually capturing their images in each model. The data were randomized and analyzed qualitatively and quantitatively by 2 examiners. The oblique loading was significantly more damaging. The increase in length was favorable for stress distribution (p<0.05). The splinting was beneficial for the transmission of stresses mainly (p<0.05). The splinting of the crowns, as well as increasing the length of the first implant and axial loading was most beneficial in the stress distribution. Short splinted implants behaved better than single unit implants. Increasing of the length of the first implant significantly improved the stress distribution in all analyzed situations.

Biomechanical activation of vascular endothelium as a determinant of its functional phenotype

One of the striking features of vascular endothelium, the single-cell-thick lining of the cardiovascular system, is its phenotypic plasticity. Various pathophysiologic factors, such as cytokines, growth factors, hormones, and metabolic products, can modulate its functional phenotype in health and disease. In addition to these humoral stimuli, endothelial cells respond to their biomechanical environment, although the functional implications of this biomechanical paradigm of activation have not been fully explored. Here we describe a high-throughput genomic analysis of modulation of gene expression observed in cultured human endothelial cells exposed to two well defined biomechanical stimuli—a steady laminar shear stress and a turbulent shear stress of equivalent spatial and temporal average intensity. Comparison of the transcriptional activity of 11,397 unique genes revealed distinctive patterns of up- and down-regulation associated with each type of stimulus. Cluster analyses of transcriptional profiling data were coupled with other molecular and cell biological techniques to examine whether these global patterns of biomechanical activation are translated into distinct functional phenotypes. Confocal immunofluorescence microscopy of structural and contractile proteins revealed the formation of a complex apical cytoskeleton in response to laminar shear stress. Cell cycle analysis documented different effects of laminar and turbulent shear stresses on cell proliferation. Thus, endothelial cells have the capacity to discriminate among specific biomechanical forces and to translate these input stimuli into distinctive phenotypes. The demonstration that hemodynamically derived stimuli can be strong modulators of endothelial gene expression has important implications for our understanding of the mechanisms of vascular homeostasis and atherogenesis.

The anatomical and biomechanical consequences of anterior vertebral stapling for the fusionless correction of scoliosis

Fusionless scoliosis surgery is an emerging treatment for idiopathic scoliosis as it offers theoretical advantages over current forms of treatment. Anterior vertebral stapling using a nitinol staple is one such treatment. Despite increasing interest in this technique, little is known about the effects on the spine following insertion, or the mechanism of action of the staple. The aims of this study were threefold; (1) to measure changes in the bending stiffness of a single motion segment following staple insertion, (2) to describe the forces that occur within the staple during spinal movement, and (3) to describe the anatomical changes that occur following staple insertion. Results suggest that staple insertion consistently decreased stiffness in all directions of motion. An explanation for the finding may be found in the outcomes of the strain gauge testing and micro-CT scan. The strain gauge testing showed that once inserted, the staple tips applied a baseline compressive force to the surrounding trabecular bone and vertebral end-plate. This finding would be consistent with the current belief that the clinical effect of the staples is via unilateral compression of the physis. Interestingly however, as each specimen progressed through the five cycles of each test, the baseline load on the staple tips gradually decreased, implying that the force at the staple tip-bone interface was decreasing. We believe that this was likely occurring as a result of structural damage to the trabecular bone and vertebral end-plate by the staple effectively causing ‘loosening’ of the staple. This hypothesis is further supported by the findings of the micro-CT scan. The pictures depict significant trabecular bone and physeal injury around the staple blades. These results suggest that the current hypothesis that stapling modulates growth through physeal compression may be incorrect, but rather the effect occurs through mechanical disruption of the vertebral growth plate.

Biomechanical force displacement analysis of strength of fixation of tracker holding devices to bone in computer aided joint replacement

Computer aided joint replacement surgery has become very popular during recent years and is being done in increasing numbers all over the world. The accuracy of the system depends to a major extent, on accurate registration and immobility of the tracker attachment devices to the bone. This study was designed to asses the forces needed to displace the tracker attachment devices in the bone simulators. Bone simulators were used to maintain the uniformity of the bone structure during the study. The fixation devices tested were 3mm diameter self drilling, self tapping threaded pin, 4mm diameter self tapping cortical threaded pin, 5mm diameter self tapping cancellous threaded pin and a triplanar fixation device ‘ortholock’ used with three 3mm pins. All the devices were tested for pull out, translational and rotational forces in unicortical and bicortical fixation modes. Also tested was the normal bang strength and forces generated by leaning on the devices. The forces required to produce translation increased with the increasing diameter of the pins. These were 105N, 185N, and 225N for the unicortical fixations and 130N, 200N, 225N for the bicortical fixations for 3mm, 4mm and 5mm diameter pins respectively. The forces required to pull out the pins were 1475N, 1650N, 2050N for the unicortical, 1020N, 3044N and 3042N for the bicortical fixated 3mm, 4mm and 5mm diameter pins. The ortholock translational and pull out strength was tested to 900N and 920N respectively and still it did not fail. Rotatory forces required to displace the tracker on pins was to the magnitude of 30N before failure. The ortholock device had rotational forces applied up to 135N and still did not fail. The manual leaning forces and the sudden bang forces generated were of the magnitude of 210N and 150N respectively. The strength of the fixation pins increases with increasing diameter from three to five mm for the translational forces. There is no significant difference in pull out forces of four mm and five mm diameter pins though it is more that the three mm diameter pins. This is because of the failure of material at that stage rather than the fixation device. The rotatory forces required to displace the tracker are very small and much less that that can be produced by the surgeon or assistants in single pins. Although the ortholock device was tested to 135N in rotation without failing, one has to be very careful not to put any forces during the operation on the tracker devices to ensure the accuracy of the procedure.