ORIGINAL RESEARCH article

Front. Mech. Eng., 06 July 2023

Sec. Biomechanical Engineering

Volume 9 - 2023 | https://doi.org/10.3389/fmech.2023.1215895

Numerical simulation of the effect of geometric parameters on silicone airway stent migration

  • 1. Department of Pulmonary Medicine, Tehran University of Medical Sciences, Tehran, Iran

  • 2. School of Mechanical Engineering, College of Engineering, University of Tehran, Tehran, Iran

  • 3. Pediatric Urology and Regenerative Medicine Research Center, Children’s Medical Center, Pediatric Center of Excellence, Tehran University of Medical Sciences, Tehran, Iran

  • 4. Tracheal Diseases Research Center, National Research Institute of Tuberculosis and Lung Diseases (NRITLD), Shahid Beheshti University of Medical Sciences, Tehran, Iran

  • 5. Institut de Recherche Robert Sauvé en Santé et en Sécurité du Travail, Montréal, QC, Canada

Abstract

Recently, endotracheal stenting has become critical in treating respiratory diseases. Due to the COVID-19 pandemic in recent years, many patients had stenosis because of long-term intubation, and silicone stents can be used to treat tracheal stenosis in these patients. Standard airway stents are silicone tubes that provide immediate relief but are prone to migration. In this work, we design different silicone stents and analyze them in the trachea to evaluate silicone airway stents’ performance to overcome undesired migration. A finite-element model of the trachea was employed to evaluate anti-migration forces in each stent. The geometry of the trachea is brought from a computerized tomography scan of the chest of a 68-year-old healthy man. The results are shown based on the least migration of stents based on anti-migration forces. Also, the conditions of stent placement have been considered based on two different assumed friction factors, and the importance of choosing the type of silicone for stent construction has been analyzed. The results show that increasing the diameter of the stent reduces the displacement and migration of it in the trachea. Furthermore, the 23 mm stent with a 45° angle revealed the best implementation against compression under the impact of respiratory pressure differences.

1 Introduction

Recently, in the wake of diseases like coronavirus, a large number of researchers have focused on the airway stent that has been applied during tracheomalacia. The placement of airway stents has come foremost in treating tracheal stenosis such as benign and malignant tracheobronchial diseases since 1915 (). Tracheal stenosis is expressed by structural tracheal constriction, and it can progress from many reasons, including external injury to the trachea, benign or malignant tumor, relapsing polychondritis, sarcoidosis, and infections (; ). An airway stent has a significant role in healing or controlling the airway trachea, and it can support the airway wall against collapse or external compression (). Airway stents can be classified into four categories: silicone, balloon-dilated metal stent, self-expanding metal stent, and covered self-expanding metal stent (). One of the advantages of metal stents is that they have less migration compared to silicone stents, but their main disadvantage is that they are difficult to remove, which makes it necessary to use silicone stents (; ; ). The ideal airway stent should 1) be easy to place and remove; 2) be large enough to maintain the position and not migrate; 3) not be very large and be able to avoid granulation tissue reactions; 4) be flexible enough to act like airway physiology but have the sufficient firmness to resist extrinsic compression, and 5) not affect mucociliary clearance (). The properties of airway stents depend on their materials, shapes, and other characteristics (). Silicone stents by customization and dimensional design can prevent migration and achieve a better fit (). In contrast to metal stents, silicone stents have various hardness and flexibilities, have the ability to achieve greater mechanical strength by reinforcement, are less expensive, and are easily manufactured by additive manufacturing ().

One procedure to improve tracheal stents’ geometry and performance is using finite-element methods to simulate and analyze the effect of structural parameters. The first implantation was performed in the 1980s to overcome and treat narrowed or blocked coronary arteries (). Due to the incomplete attachment of balloon-expandable stents to the aortic wall, which caused undesired migration, self-expanding nitinol aortic stents were presented (). The self-expanding nitinol stent can be easily crimped below room temperature, and when it arrives at the vessel wall, it expands until it hits the vessel wall and fits into it (). In particular, limitations of self-expanding stents, such as difficult adjustment, tissue (tumor) ingrowth, and being a permanent and irremovable cause to silicone stents, have been discussed in the recent decade (; ).

) modified the silicone stent by presenting the Polyflex stent, which had a thinner wall and could be more easily removed if required. ) published the clinical efficacy and safety of silicone stents in anastomotic complications of lung transplantation and concluded that silicone stents are acceptable and effective for airway disorders. Stents must be able to respond to cough and breathing pressures by undergoing reversible reduction of the cross-sectional area. Considering the breathing and coughing processes in the healthy trachea after stent implantation is necessary to design better-adapted and personalized stents. Computational modeling of the stenting procedure provides a powerful tool for predicting the patient’s possible risk of a specific stent (). , in order to analyze human tracheal conformability after stent implantation, developed a finite-element model of the trachea based on a computerized tomography (CT) scan under normal breathing and coughing pressure, and they showed that the stent could prevent tracheal muscle deflections, particularly during coughing. ) performed numerical simulations for eight different stent geometries inserted into trachea models. They demonstrated a clear interrelation between the stent diameter and the stress applied to the trachea wall, and silicone stents produced the lowest levels of stress, which explained their predilection to migration. Luo et al. characterized the inspiratory flow and modeled the human lung by employing a CT scan. Zhang et al. (2016) used the CT scan to reconstruct 3D human tracheal stenosis. They implanted the stent and then analyzed the contact between the trachea wall and the implanted stent using finite-element methods. They have shown stress concentration after stenting is directly linked to granulation tissue. Safshekan et al. evaluated the mechanical properties of the human trachea, cartilage, connective tissue, and smooth muscle using uniaxial tensile tests and found that Yeoh and Mooney–Rivlin models were best able to describe the hyperelastic behavior of all three tracheal components (). evaluated non-linear finite-element models to obtain the effect of malacic segment length, tracheal diameter, and reduction in tissue mechanical properties on the mechanism and intensity of tracheal collapse under exhalation pressure. A linear viscoelastic constitutive model was developed to determine the crimping behavior of a viscoelastic polymeric braided stent and demonstrated the possibility of predicting the behavior of the braided stent during crimping and uncramping process (). proposed a fiber-reinforced silicone for tracheobronchial stents to analyze the reinforcement test for silicone with three different fibers and to compare these fibers in terms of their mechanical properties. ) presented the mechanical effects that adhered covers may have on stent performance. Zheng et al. (2019) utilized the finite-element method to obtain the optimum design of braided stents. They characterized the effects of the geometrical parameter in self-expanding stents, such as crossing pattern, braiding angles, wire diameters, and strand numbers, on the mechanical performance. In addition, both the geometry and contact issues in self-expanding stents were analyzed by importing three braided stent models in a finite-element code and used to compare the response of different contact modeling techniques to different idealized boundary conditions ().

This study aimed to evaluate silicone airway stents’ performance to overcome undesired migration. In this study, the simulated trachea has been performed based on experimental data, which was reported by ). In order to optimize the airway stent, different stents were evaluated considering parameters such as diameter, hollow or solid type, and angle of the stud. In addition, the effect of the friction factor was also investigated. In all these cases, anti-migration forces are considered our benchmark. Moreover, the effect of the stent design parameters on stent removal and its undesirable displacement due to the difference in respiratory pressure is analyzed under radial and axial (longitudinal) forces.

2 Materials and models

Silicone airway stents were designed to minimize migration, mucostasis, and granulation tissue formation. Since the shape and size of the trachea depend on age and sex, a trachea was considered our fixed criterion. Then, the silicone stents and their types used in this research are mentioned.

2.1 Solid model of the trachea

The trachea is a cartilaginous tube structure in the human body that attaches the larynx to the bronchi of the lungs, which allows airflow to pass through. In order to simulate the implantation of different stents under distinct conditions in the trachea, a finite-element method was used. The geometrical model of the human trachea was made based on the high-resolution CT scan of the chest from 40 males and 20 females with an age range of 22–88, as reported in morphometric studies (). All patients had healthy trachea with no distorting intrathoracic pathology and used a standardized breath-holding technique for performing scans. The finite-element model of the human trachea in this research was made based on a CT scan performed, as reported in Table 1 and Figure 1.

TABLE 1

No. of tracheal ringsTracheal length (mm)Ring numberA-P diameter (mm)Transverse diameter (mm)Vertical height of tracheal ring 6 (mm)Vertical height of second tracheal ring (mm)Inter-ring distance (mm)Distance between posterior ends of rings (mm)
122.621.64.51.8
19102.6522.420.53.53.31.915.6
1022.917.62.32.4
1519.619.93.61.5

Tracheal dimensions from a computerized tomography scan of the chest of a 68-year-old healthy man were used to simulate the cartilage and muscular membrane. For tracheal dimensions, see the work of ).

FIGURE 1

). A commercial software SolidWorks was utilized to create a 3D tracheal model with 19 rings of 102.6 mm in length.

Finally, with the software ABAQUS, a full four-node linear tetrahedron mesh of approximately 101,992 elements for the healthy trachea was generated. For the thickness, delicate variations along the tracheal axis were seen and, therefore, chosen to model it as a constant.

The tracheal components comprise three parts: cartilage, smooth muscle, and connective tissue, and each part has different mechanical properties. ) suggested that every three parts are considered incompressible and isotropic hyperelastic tissues. Moreover, the Yeoh hyperelastic model was used to determine coefficients for tracheal tissues, where Ci are material constant coefficients and I1 is the first invariant of the Cauchy–Green tensor defined by , where , , and are principal stretches, as the obtained coefficients for tracheal are reported in Table 2.

TABLE 2

Tracheal tissue
Strain energyTracheal tissueC10/MPaC20/MPaC30/MPa
Cartilage3.583−2.53412.02
Smooth muscle0.0630.394−0.171
Connective tissue0.2570.483−0.148
Silicone stent
C1/MPa
1.2

Material parameters of tracheal tissues borrowed from the work of ) and stent based on the work of ).

In this research, different stents were implanted using one tracheal model, and the result may vary in other dimensions of other human trachea. The top and bottom surfaces of the trachea are fixed when placing the trachea. However, when we want to remove it, we deactivated the top fixed surface and only the bottom fixed boundary remained to capture the behavior of the trachea.

2.2 Stent design

Airway stents prevent tracheal obstruction, such as bronchogenic carcinoma and malignant central airway diseases, and enhance pulmonary function in patients suffering from airway obstruction (). When designing stents, advantages such as minimal granulation, prevention of migration, resisting external compression, and ease of removal must be considered. A silicone stent is a tube-like device that is widely used in airway obstruction due to its benefits, such as low cost and easy removability. The outside surface of the silicone stent maintains small studs to prevent proximal or distal migration into the trachea (). In order to obtain an optimized stent, 64 different studded silicone stents were designed and compared. As shown in Figure 2, in the center of half of the designed stent models, a hole was created into the studs to prevent mucus plugging and decrease the radial force embedding the stent into the surrounding tissue. These holes allow mucus to pass through them and decrease the radial displacement of stents under breathing pressure. Therefore, 64 designed models were simulated with four different inner diameters and stud angles, including hollow studs and solid studs. The internal diameter, length, and thickness of the stents are reported in Table 3.

FIGURE 2

TABLE 3

Diameter (mm)Length (mm)Angles of the studThickness (mm)Stud diameter (mm)Stud hole diameter (mm)
20
 21
22500, 15°, 30°, and 45°1.531
23

Geometrical characteristics of the stent models.

) represented the methodology to determine the properties of silicone stents. The stents were considered hyperelastic and incompressible, and a neo-Hookean model was used (; ). In this function, I1 is the first invariant of the Cauchy–Green tensor as explained earlier and C1 is the material constant. The coefficient for silicone stents in the current study is shown in Table 2 for shore 48. It should be noted that the penalty friction model has been implemented, and the effect of different friction coefficients has been studied.

Finally, stents’ CAD files imported in ABAQUS and a full four-node linear tetrahedron mesh of approximately 24,457 elements for each stent were generated. In all cases mentioned, anti-migration forces are investigated. In addition, the effects of coughing and sneezing were considered by applying external pressure to stents.

3 Results and discussion

The pressure difference between breathing and coughing in distinct sections of the trachea extracted stents from the trachea. In addition, a pressure difference between the inner and outer stent wall might reduce the diameter and cause stent disposal out of the trachea. The resulting behavior is significantly different in each case. Results corresponding to case 1 are shown in Figure 3. According to the previous section, three parameters for designing stents are considered: diameter, the angle of studs, and existing holes in studs, and each parameter influences stent displacement and needed pressure to anti-migration forces resistance stents into the trachea. In addition, the results for the healthy trachea are considered after stent implantation. The distribution of von Mises stress at the stent and trachea during the growth of the silicone stent while breathing in the healthy trachea is shown in Figure 3. Even though this design structure might not conform precisely to the patient’s actual occasions, utilizing different models facilitates the mechanical aspects of silicone stent function, preventing lung injuries (; ; ).

FIGURE 3

As shown in Figure 4, the longitudinal forces for removing the stent from the trachea are compared in different diameters and other parameters are assumed to be constant.

FIGURE 4

As predicted, individually increasing the diameter of stents increases the force required to remove the stent from the trachea, and the possibility of displacement and migration in larger diameters decreases. Thus, selecting a larger diameter of the stent is also appropriate to avoid patient discomfort. Furthermore, it should be considered that deploying a stent in a larger diameter can cause injury to the tracheal wall and may be a risk factor for long-term usage ().

As explained in the previous section, the placement of studs on the external stent tube wall prevents migration because of good contact pressure between the airway wall and the studs. However, it is not ideal for tracheobronchomalacia or to bridge tracheoesophageal fistula (). In this study, an angle has been defined for the studs, which changes the contact surface between the tracheal wall and the stent, and then, its effects on reaction force have been analyzed (Figure 2). The angle is indicated by α, and its value varies between 0, 15°, 30°, and 45°. The result corresponding to the effect of stud angle difference on reaction forces for removing stent from the trachea is shown in Figure 5.

FIGURE 5

Reaction force measurements have been reported based on stud angle (Figure 5), improved mainly by increasing the angle from 0° to 45° because the contact region between the tracheal wall and stent increased. The results show that existing studs with a higher angle can improve the function implantation of a prosthesis, especially at the angle of 45°; pressure difference, which causes displacement at coughing, can be controlled. When considering the impact of influential factors on stent migration, it is essential to mention that the influence of angle is less significant than the effectiveness of diameter on stent migration.

As mentioned in the previous section, the placement of hollow studs on stents makes mucus feasible to pass through them and then prevent mucus plugging. In order to analyze the effect of having a hole on studs, reaction forces were determined and reported in Figure 6.

FIGURE 6

Figure 6 shows that hollow studs reduce reaction forces and cause stent displacement easier than solid studs. However, this decrease is insignificant, and as for hollow stud benefits, this reduction force can be deficient. Furthermore, it is evident that this reduction differential at reaction forces is more significant in the larger diameter of the stent. In addition to the function of mucus, the presence of holes in studs decreases the pressure difference between breathing pressure and coughing pressure that is being applied to the stent. The mucus may build-up in the airways and obscures bronchoscopic visualization (). In the case of silicone stents, immediate intubation using the hollow stud stents, passing distal of the stenosis, and cleaning the distal airways of pus and mucus can be lifesaving. Hence, although the presence of holes on silicone stents increases a slight probability of undesirable migration, it reduces the concern of mucus plugs.

According to the silicone properties, friction should be on the silicone stent’s contact surface and the trachea wall, which prevents the stent from coming out of the lungs due to the pressure difference. In previous analytical works, the effect of friction was not considered, which contradicts the reality of the human respiratory system. As illustrated in Figure 7 in this study, two small friction factors have been assumed to estimate the effect of changing the friction factor on the reaction forces: 0.1 and 0.05.

FIGURE 7

Figure 7 shows that friction significantly affects the reaction forces and, subsequently, the migration of stents. The effect of friction has increased significantly by doubling the friction factor, and it can explain significant differences in reaction forces. The effect of friction evolves essentially when silicones with different properties are used to construct stents. All effective parameters, including the angle of the stud (α), friction factor, stent inner diameter, and hole on studs, are reported in Figure 8 based on the maximum reaction forces required to extract stents from the trachea.

FIGURE 8

Figure 8 shows that in the friction factor 0.1, the maximum reaction forces belong to 23 diameters, an angle of 45° at the hollow stud state and 23 diameters and 45° at the solid stud state. It is evident that by increasing the diameter of the stent, the force laying the trachea increases, and the possibility of the stent migrating out of the trachea decreases. As the α angle increases, the contact surface between the trachea and the stent increases, which pushes the stent to migrate out or displace under the respiratory pressure difference condition. The results for the friction factor 0.05 are similar to those for the friction factor 0.1. By comparing the results of Figure 8, it is apparent that with the increase of the friction factor between the inner trachea wall and the silicone stent, the pressure difference required to migrate out of the stent increases. Therefore, the essence of preferring the silicone grade employed for the stent is figured out. Finally, according to the selected trachea dimensions, the optimum model with the best performance for opening the respiratory tract and having the most negligible migration is a model with a diameter of 23 mm and an angle of 45° in the solid stud state.

The primary aim of stent implantation is overcoming the extrinsic force or central obstruction causing a loss of the airway structure (). The different external radial forces over the trachea wall, such as coughing, may cause necrosis of the mucosa and fistula formation, leading to the compression of the silicone stent and removing it from the trachea if the external force is significant (; ). The compressed silicone stent exerts a radial force, such as breathing and coughing, which may have weaker function, obstruction, or migration (). Furthermore, the high external radial force reduced the stent radius and caused stent displacement and may set patients at high risk; thus, selecting a suitable diameter and properties for the stent is necessary due to low migration. In this investigation, the external radial pressure was applied to hollow stud stents and caused to compress hollow stents. The outer stent wall pressure assumed a 10 kPa–60 kPa range, where forces and structure differences were simulated. Considering that the stent with a 45° angle had the best performance in keeping the airway open, a hollow stent with a 45° angle is used in this section.

Figure 9 shows the force required to remove the stent decreases by increasing the external pressure applied to the stents. However, when the pressure rises above a specific limit due to the asymmetric pressure that the trachea puts on the stents, the stents turn into Figure 9B that has lost its function. Therefore, although more force is required to remove the stent in a larger diameter, more external pressure is exerted when coughing or sneezing, reducing its efficiency and performance. Finally, in each trachea, there is an optimal diameter that, in addition to requiring more anti-migration force, external pressure has less effect on it. Therefore, according to ) study, the maximum static pressure is 60 kPa, and a stent with a diameter of 21 has a better performance than a stent with a diameter of 23.

FIGURE 9

4 Conclusion

Today, stenting the central airways is a common operation in tracheal airway diseases. Treatment of respiratory diseases due to stent migration and mucus formation is still challenging using silicone stents, but endotracheal stenting provides an alternative to surgical resection in patients with tumor-related extraluminal compression in the central airway. Therefore, it is necessary to understand the stent’s impact on the trachea and the factors affecting its performance better, including the material and dimensions of its design. In order to facilitate the patient’s healing process, various types of stents with different dimensions, geometries, and materials have been designed and produced, and one is preferred and placed according to the patient’s requirement.

In this study, 64 different stent models were designed based on differential geometry and sizes and were placed in the trachea with specific dimensions and analyzed. The results were reported based on the least migration of stents based on anti-migration forces. Furthermore, the conditions of stent placement have been evaluated based on different friction factors, and the importance of selecting the type of silicone for stent construction has been investigated. Based on the reported results, stent diameter and assumed friction are more influential than other factors on the migration possibility of the stent in the trachea, and both directly affect reaction forces. By choosing larger diameters, the reaction forces increase which can prevent stent migration or displacement. Considering the higher fraction factor, the force required to migrating stents increases, reducing the undesirable displacement of the stent in the trachea.

Also, a stent model with a hole on studs has been presented, allowing the mucus’s path. The results of finite-element analysis showed that this model allows the passage of mucus simultaneously with a slight increase in migration. Another geometric parameter that has been studied in this research is the effect of the stud angle on the maximum reaction force and migration of stents. According to the results, with the increase in the stud angle from 0º to 45°, the maximum reaction force also increased, which can affect the holes’ presence and partially moderate the increase in stent migration. In the continuation of this analysis, the 21 mm stent, with a 45° angle, revealed the best implementation against compression under the impact of the respiratory pressure difference. Finally, in each trachea, there is an optimal point for the diameter that, in addition to requiring more anti-migration force, external pressure has less effect on it.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

MF: methodology and writing—original draft. AB: methodology, simulation, and writing—original draft. SB: simulation and writing—original draft. AbK: conceptualization. ‬ArK: conceptualization. MaB: conceptualization and methodology. AE: co-supervision. MoB: co-supervision. All authors contributed to the article and approved the submitted version.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Summary

Keywords

airway silicone stents, tracheal simulation, stent migration modeling, anti-migration forces, optimizing stent geometry

Citation

Farahani MM, Bakhtiyari A, Beshkoofe S, Kajbafzadeh A, Kiani A, Eskandari AH, Baniassadi M and Baghani M (2023) Numerical simulation of the effect of geometric parameters on silicone airway stent migration. Front. Mech. Eng 9:1215895. doi: 10.3389/fmech.2023.1215895

Received

02 May 2023

Accepted

21 June 2023

Published

06 July 2023

Volume

9 - 2023

Edited by

Manuel Doblare, University of Zaragoza, Spain

Reviewed by

Arya Amiri, University of Waterloo, Canada

Hamid Shahsavari, Karlsruhe Institute of Technology (KIT), Germany

Updates

Copyright

*Correspondence: Mostafa Baghani,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics