Here we describe a method to manufacture aneurysmal, aortic tissue-mimicking phantoms for the use in testing ultrasound elastography. The combined use of computer-aided design (CAD) and 3-dimensional (3D) printing techniques produce aortic phantoms with predictable, complex geometries to validate the elastographic imaging algorithms with controlled experiments.
Ultrasound (US) elastography, or elasticity imaging, is an adjunct imaging technique that utilizes sequential US images of soft tissues to measure the tissue motion and infer or quantify the underlying biomechanical characteristics. For abdominal aortic aneurysms (AAA), biomechanical properties such as changes in the tissue’s elastic modulus and estimates of the tissue stress may be essential for assessing the need for the surgical intervention. Abdominal aortic aneurysms US elastography could be a useful tool to monitor AAA progression and identify changes in biomechanical properties characteristic of high-risk patients.
A preliminary goal in the development of an AAA US elastography technique is the validation of the method using a physically relevant model with known material properties. Here we present a process for the production of AAA tissue-mimicking phantoms with physically relevant geometries and spatially modulated material properties. These tissue phantoms aim to mimic the US properties, material modulus, and geometry of the abdominal aortic aneurysms. Tissue phantoms are made using a polyvinyl alcohol cryogel (PVA-c) and molded using 3D printed parts created using computer aided design (CAD) software. The modulus of the phantoms is controlled by altering the concentration of PVA-c and by changing the number of freeze-thaw cycles used to polymerize the cryogel. The AAA phantoms are connected to a hemodynamic pump, designed to deform the phantoms with the physiologic cyclic pressure and flows. Ultra sound image sequences of the deforming phantoms allowed for the spatial calculation of the pressure normalized strain and the identification of mechanical properties of the vessel wall. Representative results of the pressure normalized strain are presented.
Abdominal aortic aneurysms (AAA) are focal enlargements of the aorta that occur preferentially near the aortic bifurcation1. The exact cause of AAA formation is unknown, although many theories suggest that the pathogenesis is multifactorial, with genetic, behavioral, hemodynamic, and environmental factors contributing2,3. While the diagnosis of an abdominal aortic aneurysm can be obtained using non-invasive imaging techniques, the prediction of patient-specific rupture risk is not as precise4,5,6. Surgical repair can reduce the risk of aortic rupture, but operative repair of the aorta carries a high rate of associated morbidity and mortality7. Current surgical practices use the "maximum size criteria", or maximum absolute diameter of the aneurysms, to predict a patient's risk of rupture. Unfortunately, it has been well established that an aneurysm still ruptures below sizes clinically acceptable for surgical repair, meaning that patients with any sized aneurysm carry some risk of rupture8,9,10,11,12,13. Additionally, it is known that historical reports of rupture risk are likely over-estimations of the true rupture risk, meaning many patients are exposed to surgical risk without benefit13. A more accurate assessment of the patient-specific rupture risk is needed to help stratify a patient's risk-benefit ratio of undergoing surgical aneurysm repair.
It has been shown that the spatial stress distribution within an AAA is of critical importance in determining rupture potential and may be a better indicator than maximum diameter14,15,16,17,18. Most of the recent studies that investigate the mechanics of AAA rupture use segmented geometries from X-ray computed tomography (CT) images, and population averaged mechanical properties of aortic tissue measured ex vivo. Finite element (FE) models are then used to predict the vessel wall stresses14,15,16,17,18. However, because the mechanical properties are determined following the tissue excision, it is unclear whether the resulting models accurately depict the resulting in vivo patient-specific stresses. These studies typically assume homogeneous vessel wall material properties and don't account for the highly heterogeneous structure of the aortic wall and thrombus19,20,21,22,23,24,25.
Ultrasound-based elasticity imaging is used clinically to diagnose and monitor a variety of disease pathologies26. This technology provides a non-invasive means to interrogate the physical interactions of soft tissues. Vascular US elasticity imaging has been used as an adjunct imaging modality to clinical US evaluation in the screening and monitoring of AAAs. The combination of these techniques provides both geometric information, such as diameter and length, as well as mechanical data, such as relative stiffness and stiffness variation. While many elasticity imaging techniques require an external load to induce a measurable tissue deformation, the tissue motion to be measured here is induced by changes in the aortic pressure caused by the beating heart. Numerous methods have been published to spatially resolve strain fields in deforming vessels, however, validation studies of these methods have been limited to human patients, animal models, or ex vivo tissue samples27,28,29,30,31,32. To date, few methods allow for creations of custom geometries with spatially varied material properties27,29.
Here we present a method of manufacturing US compatible, tissue-mimicking phantoms that can be tailored to a variety of relevant aortic geometries and material properties for validation of US elastography techniques. Although previous groups have been able to design complex geometry phantoms to mimic AAA geometries using 3D printing technology33,34, printable rubbers are known to have a high attenuation to US and do not have a means to later their material properties. Phantoms are made from polyvinyl alcohol cryogel (PVA-c), which has been previously shown to be ideal for mimicking vascular tissue properties35. These phantoms can be used in US, magnetic resonance, and elastographic imaging36,37,38. The aortic aneurysm geometry was designed similarly to that of the simulation model created by Vorp et al.14. The vessel has a nominal diameter of 22.5 mm and has an aneurysmal bulge that is 64 mm bulge long, 47 mm in diameter and eccentric (β = 0.6)14 to the anterior side of the phantom. The last section mimics the iliac bifurcation with a distal diameter of 15 mm. The phantom was chosen to have a constant thickness of approximately 5 mm. Raghavan et al. reported in a small study that the vessel thickness of AAA ranges from 0.23-4.26 mm, with a median value of 1.48 mm39. A nominal vessel thickness on the larger end of that spectrum was chosen here for manufacturing concerns with the expectation that improved 3D printing techniques will improve the minimum phantom thickness that is able to be molded. Phantom molds were designed in CAD and are 3D printed using commercially available printers and filament.
The molds are injection filled with the PVA-c solution and subjected to a series of freeze/thaw cycles (-20 °C and +20°C) to cross-link the PVA-c polymer and polymerize the gel. The elastic modulus of the PVA-c is controlled by altering the concentration of the PVA-c gel or the number of freeze-thaw cycles. The aneurysmal section of the phantom required loss mold to remove from the inner lumen of the vessel. This was accomplished by the use of a polyvinyl alcohol, 3D printer filament (PVA). Although chemically similar to the PVA-c powder, the PVA filament does not polymerize when frozen and, as such, can be dissolved in water after the PVA-c has been set. Additional sample molds are printed to create tensile testing specimens, in a "dog bone" configuration, with the same PVA-c concentration. These molds undergo the same freeze/thaw cycles and are used for tensile testing to independently measure the elastic modulus of the phantom sections. A background material was manufactured with softer PVA-c, made to simulate tissues of the retroperitoneum40,41. This background phantom was manufactured as a homogeneous axisymmetric cylindrical tube with a 4 cm inner diameter, a 16.5 cm outer diameter, and a length of 16.5 cm. It was made from a 5% PVA solution and subjected to a total of two freeze-thaw cycles.
The final AAA phantoms were placed in the background phantom and connected, via tube fittings and clamps, to a hemodynamic water pump designed to deform the phantoms with physiologic cyclic flows and pressures. The pump speed was set to deliver approximately a 6-7 kPa pressure pulse at a rate of approximately 1 Hz. Ultra sound image sequences of the deforming phantoms were collected, and the pressure normalized strain was calculated to identify differences in the spatially varied mechanical properties. Representative results of the pressure normalized strain images within the vessel region are presented. The increasing regional differences in the normalized strain of the stiffer heterogenous phantoms, relative to the homogeneous phantom, demonstrate the differences in the vessel stiffness and our ability to measure it.
1. Download STL Models from the NIH 3D Print Exchange
2. 3D Printing of Molds
Figure 1: CAD Representation of Phantom, Background and Sample Molds. (a) – (b) 3D CAD images of vessel mold and orientation of parts for assembly. Registration spacers (i), pins (ii), holes (iii) and the fill hole are shown. (c) Drawing of inner lumen highlighting the inner vessel dimensions. (d) CAD rendering of the sample molds. (e) CAD rendering of the background phantom mold. Please click here to view a larger version of this figure.
Figure 2: Vessel Phantom Mold Assembly and Final Vessel Phantom. (a) The final printed mold of inner and outer lumen molds. The distal end of the inner lumen printed in a dissolvable PVA plastic and is attached to the proximal end of the inner lumen mold using deformable wax. (b) Tubing attached to the injection port of the outer lumen mold and the syringe stopper. (c) The inner lumen mold after spray coating of flexible sealant. (d) Assembly of bulge side of the outer lumen mold and inner lumen mold with PVA-c (dyed red) added for stiff aneurysm phantoms. (e) Full vessel mold assembled and clamped. (f) Deformable wax applied to the seams of the outer lumen mold to prevent PVA-c from leaking from the mold. (g) Final PVA-c phantom after 5 freeze/thaw cycles and removal from the mold. Please click here to view a larger version of this figure.
3. Hydrogel Preparation
4. Assembly of Molds
Figure 3: Sample mold and Final Sample and Background Phantoms. (a) Clamped sample mold and clear plastic sheet. PVA-c is poured into the sample mold and air bubbles are allowed to surface. (b) PVA-c sample after final freeze/thaw cycle. (c) Experimental US imaging setup of phantom attached to simulator pump and placed in the background PVA-c phantom. Please click here to view a larger version of this figure.
5. Phantom and Sample Testing
Figure 4: Imaging Protocol. (a) Pressure profile measured during the phantom imaging setup. (b) A representative B-mode image of the phantom at the minimum pressure. (c) B-mode at maximum pressure. Please click here to view a larger version of this figure.
Figure 5: Pressure Normalized Strain Images. Representative images of the pressure normalized strain (εp/PP) in%/kPa measured within the vessel for the homogeneous 10% by mass vessel phantom (a) and the heterogeneous phantoms with a 15% by mass (b), 20% by mass and 25% by mass anterior aneurysmal section (top of vessel). This figure has been modified from Mix et al.44. Please click here to view a larger version of this figure.
Representative B-mode images of the vessel mimicking phantoms are shown for the minimum and maximum pressures measured by the catheter (Figure 4b and 4c, respectively). The pressure-normalized strain (εp/PP) in%/kPa is shown for four different manufactured phantoms (Figure 5). Figure 5a shows the measured pressure-normalized strain within a homogeneous phantom manufactured with a 10% by mass PVA-c solution. The ratio of the average strain measured within the posterior quarter (image bottom) of the phantom to the average strain in the anterior quarter (image top) was 0.92. Figure 5b shows εp/PP for a phantom in which the aneurysmal section of the phantom was manufactured with a 15% by mass PVA-c solution and the remainder of the phantom was made using the 10% by mass PVA-c. The ratio of posterior to anterior strain for this phantom was found to be 1.87. Figure 5c shows εp/PP for the heterogeneous phantom with 20% by mass PVA-c, with a posterior to anterior strain ratio of 4.23. Figure 5d shows εp/PP for the heterogeneous phantom with 25% by mass PVA-c, with a posterior to anterior strain ratio of 7.37.
The results presented here show that abdominal aortic phantoms were created with complex geometries and spatially varying material properties. Design of phantom geometries were, or more specifically, phantom molds were done using computer software which facilitates phantom geometry alterations (Figure 1a and 1b). Molds can be readily 3D printed and assembled and complex mold geometries can be printed using PVA filament and removed, similar lost-wax casting techniques. The final vessel phantoms can be dynamically pressurized and are stable under large loads (Figure 4a). Phantoms are compatible with ultrasound imaging (Figure 4b and 4c) and have material properties mimicking abdominal aortic stiffnesses. Variations in the strain ratios in the anterior to posterior regions of the strain images demonstrate that the regions have varying material properties (Figure 5) and the independent mechanical testing on the samples quantify the exact values of their respective shear moduli.
This paper presents a technique to manufacture tissue-mimicking phantoms for the use in testing elastography or elasticity imaging algorithms. The combined use of CAD and 3D printing allows for efficient design of aortic mimicking phantoms with complex geometries, beyond tubular phantoms, including aneurysmal bulges. The creation of the phantom is done in 4 steps; 1) design of the phantom geometry, 2) printing of the phantom mold parts, 3) mixing of the cryogel solutions which will ultimately mimic the ultrasonic characteristics and mechanical properties of the phantom vessels and 4) pouring/injecting of the cryogel solution into the mold, setting the PVA-c with freeze-thaw cycles and removal of the phantom from the mold. The use of CAD in the in the design of the molds obtained in Step 1 allows for a simple means to precisely modify the geometry of the phantoms. Printing of the mold parts currently takes approximately 5-8 h depending on the size of the print and thus can easily be made for repeated modifications to the molds.
In Step 3, the cryogel solutions are created to mimic the vessel, aneurysm and background tissue with the calcium carbonate particles mimicking the US scatter of the tissue. The cryogel solutions should be stirred prior to use if the calcium particles have settled out of the mixture. The exact concentration of the mixture of the cryogel will determine the final mechanical properties of the phantoms. Thus, it is important to create the independent samples of each of the solutions used in the phantom vessel and background. Although not part of the protocol here, independent measurements of the sample's elastic modulus should be obtained using uniaxial tension testing. Independent mechanical testing of the PVA-c samples for the 10%, 15%, 20% and 25% phantoms created in the representative results had measured shear modulus of 17.4 ± 1.0 kPa, 48.3 ± 5.7 kPa, 95.1 ± 0.4 kPa and 170.0 ± 4.1 kPa, respectively.
Step 4 is the most critical step in creating these phantoms. Although the registrations pins are in place to keep the mold parts in their proper positions relative to the others, it is important to ensure that the mold parts do not separate during the molding process. Thus, the use of the clamps to hold the mold. The most important consideration of Step 4 is to minimize air bubbles trapped in the mold prior to the first freeze-thaw cycle. It is often useful to disassemble one side of the outer mold and inspect the phantom after the first freeze-thaw cycle to ensure it formed properly. This can save time wasted putting a "bad" phantom through additional cycles. Once the phantom has been removed from the mold entirely, it can be stored in water for several weeks with continued use.
The PVA-c phantoms developed in this work were created to specifically mimic the ultrasonic and material stiffness of aortic tissue. The use of polyvinyl alcohol cryogel allows for a wider range of possible mechanical stiffness, to better mimic the changing material properties of aortic tissue compared to more rubber like materials33,34. In addition, the use of hydrogel and investment casting better captures the acoustic properties of casted rubbers or directly 3D printed materials33,45. Some air bubbles can get trapped in our molds prior to the first freeze-thaw cycle. This can cause gaps in the phantom and lead to material weakness or acoustic artifacts. Thus, it is recommended to inspect phantoms out of the mold after the first freeze-thaw to determine if the process should be restarted. In addition, the authors have found that the inner mold can sometimes shift during the freezing of the aneurysmal portion of the phantoms. If this occurs, one modification of the above protocol would be to create a 3D printed, or otherwise designed, part to firmly hold the inner lumen mold to the anterior outer mold during the freezing of this section. The authors have found that using the posterior side of the outer mold and a 5 mm spacer between the posterior outer mold and the inner mold works well for this purpose.
The phantom developed here is ideal for studying the influence of changes in aneurysmal diameter, and luminal thickness or potentially the presence of thrombus in the tissue by editing the original CAD files. However, previous work has also shown that this manufacturing technique can be modified to produce patient-specific phantom geometries using computed tomography images and segmentation software, rather than CAD design, to create the 3D printed phantom molds44. The results shown here demonstrate that the algorithm was able to visualize the manufactured variations in the mechanical properties of the phantom sections. It should be noted that although these phantoms were used to test US-based imaging techniques, they are also compatible with magnetic resonance and computed tomography imaging systems and that they may also be used beyond the purpose of elasticity imaging, for a wide range of novel imaging techniques and modalities.
The authors have nothing to disclose.
This work was supported by the National Center for Advancing Translational Sciences of the National Institutes of Health through Award No. UL1 TR000042 and the National Institute of Biomedical Imaging and Bioengineering of the National Institutes of Health through Award No. R21 EB018432.
PLA filament | MatterHackers, MatterHackers.com | MEEDKTKU | |
PVA filament | MatterHackers, MatterHackers.com | M4MJTECR | |
LeakSeal | RPM International Inc., Rustoleum.com | 265495 | |
Polyvinyl alcohol powder (Elvanol 71-30) | DowDuPont Inc., ChemistryStore.com | SKU: 81015 | |
Calcium Carbonate Powder | greenwaybiotech.com via amazon.com | Amazon: B00HFFCBYQ | |
Tacky Wax | bards.com via amazon.com | Bards: BB759 Amazon: B016KBDYRS |
|
Rostock max 3D Printer | SeeMeCNC, seemecnc.com | SKU: 84459 | |
Onshape CAD software | OnShape, onshape.com | ||
Mattercontrol printer software | MatterHackers, MatterHackers.com | ||
Mikro-Cath pressure catheter and device | Millar, Inc., millar.com | 4501016/B | |
BNC digital acquisition | National Instruments Corporation, ni.com | NI USB-6251 BNC | |
clear cast acrylic sheet | mcmaster-carr Supply Company, mcmaster.com | 8560K274 | |
Cole-Parmer Stopcocks with Luer Connections; 3-way; male lock, Non-sterile | Cole-Parmer, coleparmer.com | EW-30600-02 | |
BD Disposable Syringes (60 mL, Luer lock) | Cole-Parmer, coleparmer.com | EW-07945-28 | |
6 Inch Ratchet Bar Clamp / 12 Inch Spreader | Tekton, Inc., www.tekton.com | 39181 | |
Tygon PVC Clear Tubing | mcmaster-carr Supply Company, mcmaster.com | 6516T53 | |
MTS Qtest Q/5 | MTS Systems Corperation, www.mts.com | 4501016 | |
MTS 5N Load Cell | MTS Systems Corperation, www.mts.com | 4501016/B | |
Abaqus FEA | Dassault Systèmes, 3ds.com |