Peak wall stress measurement in elective and acute abdominal aortic aneurysms
Article Outline
Background
Abdominal aortic aneurysm (AAA) rupture occurs when wall stress exceeds wall strength. Engineering principles suggest that aneurysm diameter is only one aspect of its geometry that influences wall stress. Finite element analysis considers the complete geometry and determines wall stresses throughout the structure. This article investigates the interoperator and intraoperator reliability of finite element analysis in the calculation of peak wall stress (PWS) in AAA and examines the variation in PWS in elective and acute AAAs.
Method
Full ethics and institutional approval was obtained. The study recruited 70 patients (30 acute, 40 elective) with an infrarenal AAA. Computed tomography (CT) images were obtained of the AAA from the renal vessels to the aortic bifurcation. Manual edge extraction, three-dimensional reconstruction, and blinded finite element analysis were performed to ascertain location and value of PWS. Ten CT data sets were analyzed by four different operators to ascertain interoperator reliability and by one operator twice to ascertain intraoperator reliability. An intraclass correlation coefficient was obtained. The Mann-Whitney U test and independent samples t test compared groups for statistical significance.
Results
The intraclass correlation coefficient was 0.71 for interoperator reliability and 0.84 for intraoperator reliability. There was no statistically significant difference in the mean (SD) maximal AAA diameter between elective (6.47 [1.30] cm) and acute (7.08 [1.39] cm) patients (P = .073). The difference in PWS between elective (0.67 [0.30] MPa) and acute (1.11 [0.51] MPa) patients (P = .008) was statistically significant, however.
Conclusion
Interoperator and intraoperator reliability in the derivation of PWS is acceptable. PWS, but not maximal diameter, was significantly higher in acute AAAs than in elective AAAs.
Abdominal aortic aneurysm (AAA) rupture is a common cause of preventable death in men >65 years.1 After AAA rupture, 50% of patients die before reaching the hospital; in addition, the mortality rate after emergency repair is 40% to 50%. This contrasts with a 3% to 6% mortality rate after elective AAA repair; hence, decision making about elective AAA repair requires careful assessment of operative mortality, rupture risk, and life expectancy.2
Currently, AAA repair in fit patients is recommended when AAA diameter is ≥5.5 cm or the patient becomes symptomatic.3 Although AAA diameter is an important predictor of rupture,4 not all large AAAs rupture, whereas 10% to 24% of ruptured AAA are <5.5 cm in diameter.5, 6 No reliable criterion therefore exists to accurately predict rupture risk for the individual patient, and thus, the decision to operate on the basis of AAA diameter alone may subject a significant proportion of patients to unnecessary surgery with significant mortality and morbidity. Patients with a stable aneurysm are more likely to die of other causes7; hence, the arbitrary setting of a single threshold diameter for elective AAA repair in all patients seems inappropriate.2
An AAA acts as a thin-walled pressure vessel and develops hoop and longitudinal stresses. Classic engineering stress analysis of simple axisymmetric shapes suggests that the stress in both directions is directly proportional to the diameter and inversely proportional to the thickness of the aneurysm wall. The relationship, however, is more complicated because AAAs have convoluted asymmetric shapes and the stress in an AAA will depend on the entire geometry, not just the maximum diameter.8 The wall stress in the complex geometries of real aneurysms thus cannot be predicted by simple analytic techniques; instead, finite element analysis (FEA) must be used.9, 10, 11 FEA is a numeric modelling technique that is used regularly by engineers in the design and analysis of structural components in many engineering applications, many of which are safety critical, for example, in the automotive, aerospace, and nuclear industries.12 FEA is accurate and reliable, provided the problem is modelled with sufficiently fine discretization and valid representation.
We and others have suggested that PWS is a better predictor of rupture than maximum AAA diameter.13, 14 This article reports our investigation of the interoperator and intraoperator reliability of FEA in the calculation of PWS in AAA and examines the variation in PWS in a large cohort of elective and acute AAAs.
Methods
Ethics
Full local ethics committee approval was obtained for this study, together with approval from the institutional research and development departments (Hull and East Yorkshire Hospitals NHS Trust and the University of Hull).
Patient population
The study included 70 patients with infrarenal AAA who had had an abdominal computed tomography (CT) scan. Forty patients were prospectively recruited before undergoing abdominal CT scan for AAA imaging, and 30 patients were recruited retrospectively after being admitted with symptoms of an acute or leaking aneurysm and undergoing CT AAA imaging. Information collected on each patient included demographics, medical history, systolic blood pressure closest to time of CT imaging (and during admission if acute), and CT images.
Material properties
We used a mathematic model for AAA tissue properties developed by Raghavan and Vorp,15 a two-parameter, nonlinear hyperelastic, isotropic, incompressible constitutive model on which the finite element analysis was performed. The material behavior is governed by the Mooney-Rivlin equation, which in uniaxial tension is given by:

Raghavan and Vorp found the coefficients of the equation to vary as follows: (0.15 < A < 0.22) and (1.18 < B < 3.56) N/mm2, with peak wall stresses changes by <4% when material properties were varied within this range.13, 15 We used values of 0.185 for A and 2.37 for B.
Finite element analysis
The finite element method works by dividing the geometry into a finite number of elements, and as the number of elements increases and the size of each individual element decreases, the accuracy of the solution improves. As a result, the number of elements used in an analysis is critical and was thus examined previously. This demonstrated that increasing the element mesh density from 9000 to nearly 20,000 elements (using seven models of increasing mesh density) made only a 1.2% difference in calculated PWS.13 Thus, increasing the number of elements beyond this level would only marginally improve accuracy. Our aneurysm models were therefore analyzed with between 20,000 and 30,000 elements depending on geometry complexity. The type of element used was ANSYS (ANSYS, Inc, Canonsburg, Pa) type SHELL 93 with a quadratic interpolation function. The thickness of aneurysm walls is difficult to measure but is generally reported to be of the order of 2.0 mm; therefore, this was the thickness that we used.
Generation of the finite element models
Full details of the procedure used to produce the finite element models are reported elsewhere.13 Briefly, abdominal CT scans are acquired using the standard AAA protocol, including bolus tracking (scan initiation at the peak of contrast uptake) with a nominal slice thickness of 3.2 mm (elective imaging) with 50% overlap with a helical pitch of 0.875. The abdomen and pelvis are imaged to visualize the infradiaphragmatic aorta down to the common iliac artery; however, CT images are only analyzed from the most proximal renal artery origin to the aortic bifurcation. These CT images are imported into Scion Image 4.0 image-processing software (Scion Corp, Frederick, Md), and each cross section of the AAA is opened as a separate image file. The wall of the AAA is marked manually to give its (x-y) profile in two-dimensions (2D), which is output as a text file. The process is repeated for each slice of the AAA, and the z coordinate of each slice is added subsequently using the slice thickness information by a simple in-house program (Visual Basic 6.0, Microsoft Corp, Redmond, Wash).
The full set of data points thus created is then imported into 3D image-rendering Rhinoceros 2 software (Robert McNeel & Assoc, Seattle, Wash) to create a 3D surface representation of the aneurysm. The resultant surface data is exported in a RAW-triangle format, which is then converted directly into a format suitable for reading into the ANSYS finite element software. Again, the software to convert the RAW data into the finite element model has been developed in-house (Visual Basic 6.0) and generates an ANSYS script file which includes not only details of the element mesh but also specification of the material properties, the constraints, and the loading conditions. The result is a detailed map of the stress and strain fields throughout the aneurysm. In common with other researchers in this field, we use the von Mises stress to give a single value reflecting the full 3D stress field at each point in the model. The whole process, from CT data to patient-specific aneurysm model, typically takes 90 minutes.
Interoperator and intraoperator reliability
Manual extraction of the AAA geometry CT slices inevitably provides the opportunity for errors to be introduced into the analysis. To assess the repeatability of this extraction process, interoperator and intra-operator reliability was assessed. Reliability is defined as the extent to which a measurement made repeatedly in identical circumstances will yield concordant results. The only variable for interoperator reliability was the operator, and the only variable for intraoperator reliability was the timing of the analysis.
For interoperator reliability, 10 CT data sets were analyzed by four different operators with different backgrounds and experience of the method. Operator 1 was the clinical research fellow with in depth knowledge of the process and wide experience of the aneurysm edge extraction technique. Operator 2, a postgraduate engineering student, had in-depth knowledge of the process and reasonable experience of the technique, whereas operator 3, another postgraduate engineering student, had vague knowledge of the process but no experience of this method of edge extraction. Finally, operator 4, a postdoctoral engineer, had no knowledge of the process or any experience of edge extraction. The 10 CT data sets were selected randomly during the study of 70 patients, with no reference to the geometry of the aneurysm concerned. A briefing meeting for the four operators before this process demonstrated the AAA edge extraction method. The operators did their analyses during a 2-week period.
Intraoperator reliability was assessed by analyzing 10 CT data sets on two occasions, 6 months apart, with no reference to the initial attempt during repeat analysis. This was performed by operator 1, and reflects results obtained by an operator with clinical experience during an initial learning curve and beyond.
Loading and constraints
The finite element models were loaded with the patients’ systolic blood pressure (SBP). In vivo, the renal arteries and the iliac arteries constrain an infrarenal AAA from deforming to some extent at the proximal and distal ends; hence to account for this, the finite element models were constrained at the proximal and distal ends.
Statistical analysis
Statistical analysis of interoperator and intraoperator reliability was done using the intraclass correlation coefficient on the calculated peak wall stress values.16 Continuous data are presented as the mean (SD).
Results
Interoperator reliability
Fig 1 shows the PWS calculated after aneurysm edge extraction on 10 different aneurysms by the four operators. The percentage average deviation from the mean varies from 5% (aneurysm 9) to 24% (aneurysm 6). The interoperator intraclass correlation coefficient was 0.71 (average, 0.91).

Fig 1.
Interoperator variability: peak wall stress calculated after edge extraction by four different operators on 10 separate aneurysms.
Intraoperator reliability
Fig 2 shows the PWS calculated after edge extraction done by operator 1 on two separate occasions to determine intraoperator variability. The graph in Fig 3 plots PWS values obtained for the same aneurysms against each other to determine R2 values and the Pearson correlation coefficient. The R2 value was 0.73, and the Pearson correlation coefficient was 0.86. The intraoperator intraclass correlation coefficient was 0.84 (average, 0.91).

Fig 2.
Intraoperator variability: peak wall stress calculated after repeated edge extraction by one operator on 10 separate aneurysms.

Fig 3.
Intraoperator variability: peak wall stress from initial and repeat calculation plotted against each other.
Patient study
Full demographic details of the 70 study patients with infrarenal AAA (40 elective, 30 acute) are summarized in Table I. No statistically significant difference was found between the groups for age, sex, and comorbidities, including hypertension, chronic obstructive pulmonary disease, and the incidence of smoking. The incidence of ischemic heart disease was 24% in the elective group, which was significantly lower than the 50% in the acute group (P = .049).
Table I. Sample patient demographics
| Demographics⁎ | AAA type | P | |
|---|---|---|---|
| Elective (n = 40) | Acute (n = 30) | ||
| Age | 76 (60-89) | 74 (60-89) | .59† |
| Sex (male/female) | 2.6:1 | 4:1 | .47‡ |
| Hypertension | 66 | 46 | .15‡ |
| Ischemic heart disease | 24 | 50 | .049‡ |
| Smoking | 62 | 78 | .40‡ |
| COPD | 23 | 37 | .26‡ |
⁎Data presented as median (IQR), ratio, or percentage. |
†Calculated by Mann-Whitney U test. |
‡Calculated by χ2 test. |
The difference in the maximum diameter between elective AAAs (6.47 [1.30] cm) and acute AAAs (7.08 [1.39] cm) was not statistically significant (Table II). However, acute AAAs had a significantly higher PWS (1.11 [0.51] MPa) than elective AAAs (0.67 [0.30] MPa; P = .008). The difference in SBP between the groups was not significant (elective, 143 mm Hg; acute, 151 mm Hg; P = .181), but to examine the association of AAA geometry to risk of rupture without the effect of blood pressure, PWS was also calculated with SBP standardized at 120 mm Hg. With SBP standardized, acute AAAs had a significantly higher PWS (0.84 [0.31] MPa) than elective AAAs (0.65 [0.25] MPa; P = .009). Localization of PWS was anterior in 60% of patients and on the right in 52%.
Table II. Abdominal aortic aneurysm data results
| Analysis result⁎ | AAA type | P† | |
|---|---|---|---|
| Electiv (n = 40) | Acute (n = 30) | ||
| Maximum diameter, cm | 6.47 | 7.08 | .073 |
| BP, mm Hg | 143 | 151 | .181 |
| Peak wall stress, MPa | |||
| 0.67 | 1.11 | .008 | |
| 0.65 | 0.84 | .009 | |
⁎Results presented as mean (SD). |
†Independent samples t test. |
Discussion
Finite element analysis is a well-established engineering tool used to model complex structures and analyze stress distributions; however, problems can occur if geometry information and model characteristics are inaccurate. Aneurysm geometry is obtained by manual outlining of the external AAA wall. Manual input of the external wall is currently required owing to the difficulties in precisely identifying the AAA wall; indeed, this sometimes poses problems even for the expert clinician. Various techniques of automated image extraction of AAA geometry are available but are currently not sufficiently precise for PWS extrapolation. Manual extraction of geometry is subject to reliability problems, thus an assessment of the significance of this variability in reliability is crucial.
The various methods to examine variability, including the Pearson correlation, are based on regression analysis and the measure of the extent that two observations from a group of subjects can be fitted to a straight line. However, this seeks to identify agreement rather than absolute agreement, and a result of 1 can be obtained even where the intercept of the line is not through the origin of the graph and the slope is not 1. It also has the limitation of only being able to compare two observations at one time, with no recognized way to combine two or more Pearson correlation coefficients.
A more accepted form of reliability measurement is the intraclass correlation coefficient (ICC), identifying absolute agreement between two or more operators. An ICC of 70% in research and 90% in clinical work as been recommended as acceptable.17 Interoperator reliability was 0.71 and intraoperator reliability was 0.84 in this study, suggesting that there is acceptable variation in our method for research, but further investigation is required before acceptance into clinical practice. Reliability can be most readily improved by operator training, and in our interoperator variability testing, minimal training and experience was given to two of the four operators.
Another approach to improve reliability would be to remove the human element. As previously stated, no reliable method for automated extraction of AAA external wall geometry from CT scans currently exists; however, our research group is investigating this problem. Automated systems with little or no manual input potentially improve the processing speed and reliability.
In the patient study, statistical significance in mean PWS is seen between the acute and elective AAA groups, but the difference in AAA diameter between the groups was not statistically significant. We believe PWS promises to provide individual rupture risk to patients with AAA, but this is very much work in progress, and further refinement of the process is on going. In addition, we recognize that AAA wall strength is as important as wall stresses in precipitating AAA rupture. Although invasive analysis of AAA wall tissue has demonstrated weaker wall in ruptured than nonruptured AAA,18 the transition from this to estimation of rupture potential by using noninvasively calculated wall strength has proven to be more difficult.19
Potential limitations of this study include:
Conclusion
The overall mortality rate associated with ruptured AAA has changed little in recent decades. The high perioperative mortality rate associated with ruptured AAA repair highlights the need to predict risk of rupture on an individual patient basis. Several attempts have been made in the past to identify risk factors associated with rupture, and several studies have analyzed the factors that influence the risk of rupture. However, only recently have studies focused on relating the risk of rupture on an individual patient basis. In this current study, we validate the earlier pilot work in a different population and demonstrate wall stress can be reliably calculated from routinely performed CT scans and may be valuable in predicting rupture risk.
Although this study has some limitations, the results are promising and support the case for a more detailed examination of the role of wall stress in assessing the risk of rupture on an individual patient basis. We also recognize that continuing advances in the estimation of wall strength are necessary to add more significance to a patient-specific peak wall stress and thus produce a meaningful estimate of individual AAA rupture. With the AAA population being generally elderly, the ability to evaluate individual, specific rupture risk, and thus perhaps offer early surgery to some patients, may significantly affect AAA-related mortality. These methods may also identify patients with a low risk of rupture who may avoid an unnecessary procedure, because not all patients with AAA die of rupture. Wall stress calculations are currently done manually, which is time consuming; however, further research is expected to culminate in a fully automated process that should significantly assist clinicians and patients in their decision making.
Author contributions
We thank the department of Vascular Radiology, Hull Royal Infirmary and the radiologists Dr Duncan Ettles and Dr Graham Robinson for their contribution in data collection. We would like to thank Dr Eric Gardiner for his invaluable assistance in statistical analysis.
References
- . Abdominal aortic aneurysm: a preventable cause of death?. J R Coll Surg Edinb. 1990;35:284–288
- . Guidelines for the treatment of abdominal aortic aneurysms (Report of a subcommittee of the Joint Council of the American Association for Vascular Surgery and Society for Vascular Surgery). J Vasc Surg. 2003;37:1106–1117
- . Mortality results for randomised controlled trial of early elective surgery or ultrasonographic surveillance for small abdominal aortic aneurysms. Lancet. 1998;352:1649–1655
- . Abdominal aortic aneurysm wall mechanics and their relation to risk of rupture. Eur J Vasc Endovasc Surg. 1999;18:487–493
- . Autopsy study of unoperated abdominal aortic aneurysms (The case for early resection). Circulation. 1977;56(3 suppl):II161–II164
- . Rupture in small abdominal aortic aneurysms. J Vasc Surg. 1998;28:884–888
- The Multicentre Aneurysm Screening Study (MASS) into the effect of abdominal aortic aneurysm screening on mortality in men: a randomised controlled trial. Lancet. 2002;360:1531–1539
- . Mechanical wall stress in abdominal aortic aneurysm: Influence of diameter and asymmetry. J Vasc Surg. 1998;27:632
- . Stress distribution in vascular aneurysms: factors affecting risk of aneurysm rupture. J Surg Res. 1993;55:155–161
- . Wall stress distribution on three-dimensionally reconstructed models of human abdominal aortic aneurysm. J Vasc Surg. 2000;31:760–769
- . The influence of aorta-aneurysm geometry upon stress in the aneurysm wall. J Surg Res. 1987;42:425–433
- . Finite element analysis: theory and practice. Harlow, Essex, England: Longman Scientific & Technical; 1992;
- A comparative study of aortic wall stress using finite element analysis for ruptured and non ruptured abdominal aortic aneurysm. Eur J Vasc Endovasc Surg. 2004;28:168–176
- . Prediction of rupture risk in abdominal aortic aneurysm during observation: wall stress versus diameter. J Vasc Surg. 2003;37:724–732
- . Towards a biomechanical tool to evaluate rupture potential of abdominal aortic aneurysm: identification of a finite strain constitutive model and evaluation of its applicability. J Biomech. 2000;33:475–482
- . An introduction to medical statistics. Oxford, UK: Oxford University Press; 2000;
- . Psychometric theory. New York, NY: McGraw-Hill; 1978;
- . Biomechanical properties of ruptured verses electively repaired abdominal aortic aneurysm wall tissue. J Vasc Surg. 2006;43:570–576
- . A biomechanics based rupture potential index for abdominal aortic aneurysm risk assessment; demonstrative application. Ann N Y Acad Sci. 2006;1085:11–21
- . Fluid structure interaction in abdominal aortic aneurysms: effect of asymmetry and wall thickness. Biomed Eng Online. 2005;4:64
- . Effect of intraluminal thrombus on aortic aneurysm wall stress. J Vasc Surg. 1997;26:602–608
- . Effect of intraluminal thrombus in patient specific models of abdominal aortic aneurysm. J Vasc Surg. 2002;36:598–604
- . Thrombus within an aortic aneurysm does not reduce pressure on the aneurysm wall. J Vasc Surgery. 2000;31:501–506
- . Biomechanical determinants of abdominal aortic aneurysm rupture. Arterioscler Thromb Vasc Biol. 2005;25:1558–1566
- . A planar biaxial contituative relation for the luminal layer of intraluminal thrombus in abdominal aortic aneurysms. J Biomech. 2006;39:2347–2354
- . The effects of aneurysm on the biaxial mechanical behaviour of human abdominal aorta. J Biomech. 2006;39:1324–1334
- . Biomechanics of abdominal aortic aneurysm. J Biomech. 2007;40:1806–1812
Competition of interest: none.
PII: S0741-5214(07)01425-5
doi:10.1016/j.jvs.2007.09.002
© 2008 The Society for Vascular Surgery. Published by Elsevier Inc. All rights reserved.
