Use of 3D printed materials as tissue-equivalent phantoms

This study used the specific example of 3D printing with acrylonitrile butadiene styrene (ABS) as a means to investigate the potential usefulness of benchtop rapid prototyping as a technique for producing patient specific phantoms for radiotherapy dosimetry. Three small cylinders and one model of a human lung were produced via in-house 3D printing with ABS, using 90%, 50%, 30% and 10% ABS infill densities. These phantom samples were evaluated in terms of their geometric accuracy, tissue equivalence and radiation hardness, when irradiated using a range of clinical radiotherapy beams. The measured dimensions of the small cylindrical phantoms all matched their planned dimensions, within 1mm. The lung phantom was less accurately matched to the lung geometry on which it was based, due to simplifications introduced during the phantom design process. The mass densities, electron densities and linear attenuation coefficients identified using CT data, as well as the results of film measurements made using megavoltage photon and electron beams, indicated that phantoms printed with ABS, using infill densities of 30% or more, are potentially useful as lung- and tissue-equivalent phantoms for patient-specific radiotherapy dosimetry. All cylindrical 3D printed phantom samples were found to be unaffected by prolonged radiation and to accurately match their design specifications. However, care should be taken to avoid oversimplifying anatomical structures when printing more complex phantoms.

A semiautomatic method to identify vertebral end plate lesions (Schmorl's nodes)

Background Context There are differences in definitions of end plate lesions (EPLs), often referred to as Schmorl’s nodes, that may, to some extent, account for the large range of reported prevalence (3.8 to 76%). Purpose To develop a technique to measure the size, prevalence and location of EPLs in a consistent manner. Study Design/Setting This study proposed a method using a detection algorithm which was applied to five adolescent females (average age 15.1 years, range 13.0 to 19.2 years) with idiopathic scoliosis (average major Cobb angle 60°, range 55 to 67°). Methods Existing low-dose, computed tomography scans were segmented semi-automatically to extract 3D morphology of each vertebral endplate. Any remaining attachments to the posterior elements of adjacent vertebrae or endplates were then manually sectioned. An automatic algorithm was used to determine the presence and position of EPLs. Results EPLs were identified in 15 of the 170 (8.8%) endplates analysed with an average depth of 3.1mm. 11/15 of the EPLs were seen in the lumbar spine. The algorithm was found to be most sensitive to changes in the minimum EPL gradient at the edges of the EPL. Conclusions This study describes an imaging analysis technique for consistent measurement of the prevalence, location and size of EPLs. The technique can be used to analyse large populations without observer errors in EPL definitions.

Physical mechanisms underlying the strain-rate-dependent mechanical behavior of kangaroo shoulder cartilage

Due to anatomical and biomechanical similarities to human shoulder, kangaroo was chosen as a model to study shoulder cartilage. Comprehensive enzymatic degradation and indentation tests were applied on kangaroo shoulder cartilage to study mechanisms underlying its strain-rate-dependent mechanical behavior. We report that superficial collagen plays a more significant role than proteoglycans in facilitating strain-rate-dependent behavior of kangaroo shoulder cartilage. By comparing the mechanical properties of degraded and normal cartilages it was noted that proteoglycan and collagen degradation significantly compromised strain-rate-dependent mechanical behavior of the cartilage. Superficial collagen contributed equally to the tissue behavior at all strain-rates. This is different to studies reported on knee cartilage and confirms the importance of superficial collagen on shoulder cartilage mechanical behavior. A porohyperelastic numerical model also indicated that collagen disruption would lead to faster damage of the shoulder cartilage than when proteoglycans are depleted.

Investigation of the mechanical behavior of kangaroo humeral head cartilage tissue by a porohyperelastic model based on the strain-rate-dependent permeability

Solid–interstitial fluid interaction, which depends on tissue permeability, is significant to the strain-rate-dependent mechanical behavior of humeral head (shoulder) cartilage. Due to anatomical and biomechanical similarities to that of the human shoulder, kangaroos present a suitable animal model. Therefore, indentation experiments were conducted on kangaroo shoulder cartilage tissues from low (10−4/s) to moderately high (10−2/s) strain-rates. A porohyperelastic model was developed based on the experimental characterization; and a permeability function that takes into account the effect of strain-rate on permeability (strain-rate-dependent permeability) was introduced into the model to investigate the effect of rate-dependent fluid flow on tissue response. The prediction of the model with the strain-rate-dependent permeability was compared with those of the models using constant permeability and strain-dependent permeability. Compared to the model with constant permeability, the models with strain-dependent and strain-rate-dependent permeability were able to better capture the experimental variation at all strain-rates (p<0.05). Significant differences were not identified between models with strain-dependent and strain-rate-dependent permeability at strain-rate of 5×10−3/s (p=0.179). However, at strain-rate of 10−2/s, the model with strain-rate-dependent permeability was significantly better at capturing the experimental results (p<0.005). The findings thus revealed the significance of rate-dependent fluid flow on tissue behavior at large strain-rates, which provides insights into the mechanical deformation mechanisms of cartilage tissues.

Bone-anchored prosthesis: Current developments worldwide and challenges

Background Individuals with limb amputation fitted with conventional socket-suspended prostheses often experience socket related discomfort leading to a significant decrease in quality of life. Most of these concerns can be overcome by surgical techniques enabling bone-anchored prostheses. In this case, the prosthesis is attached directly to the residual skeleton through a percutaneous implant (e.g., screw type fixation, press-fit implant).[46, 48, 51, 52, 77, 78] The aim of this study is to present the current advances in these surgical techniques worldwide with a strong focus on the current challenges. Methods The current advances will be extracted from a systematic literature review including approximately 40 articles. The outcomes measured will include the estimation of the population worldwide as well as the complications (e.g., infection, loosening, fractures, and breakage) and the benefits (e.g., functional outcomes, health-related quality of life).[5-19, 51-53, 55, 57, 58, 62, 73, 79] Results The population of individuals fitted with a bone-anchored prosthesis is approximately 550 worldwide. Publications focusing on infection are sparse. However, the rate of superficial infection is estimated at 20%. Deep infection occurs rarely. Loosening and peri-prosthetic fractures are fairly uncommon. Breakage of implant parts occurs regularly mainly due to fall. All studies reported a significant improvement in functional level and overall quality of life. Conclusions Several commercial implants are in developments in Europe and US. The number of procedures is consistently growing worldwide. This technique might be primary way to fit a prosthesis to young and active amputees by 2025.

Effect of decellularization on the load-bearing characteristics of articular cartilage matrix

The application of decellularized extracellular matrices to aid tissue regeneration in reconstructive surgery and regenerative medicine has been promising. Several decellularization protocols for removing cellular materials from natural tissues such as heart valves are currently in use. This paper evaluates the feasibility of potential extension of this methodology relative to the desirable properties of load bearing joint tissues such as stiffness, porosity and ability to recover adequately after deformation to facilitate physiological function. Two decellularization protocols, namely: Trypsin and Triton X-100 were evaluated against their effects on bovine articular cartilage, using biomechanical, biochemical and microstructural techniques. These analyses revealed that decellularization with trypsin resulted in severe loss of mechanical stiffness including deleterious collapse of the collagen architecture which in turn significantly compromised the porosity of the construct. In contrast, triton X-100 detergent treatment yielded samples that retain mechanical stiffness relative to that of the normal intact cartilage sample, but the resulting construct contained ruminant cellular constituents. We conclude that both of these common decellularization protocols are inadequate for producing constructs that can serve as effective replacement and scaffolds to regenerate articular joint tissue.

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.

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.

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.

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.