A finite-element model of tracheal collapse

Decellularization compromises mechanical and structural properties of the native trachea

Tracheal replacement using tissue engineering technologies offers great potential to improve previously intractable clinical interventions, and interest in this area has increased in recent years. Many engineered airway constructs currently rely on decellularized native tracheas to serve as the scaffold for tissue repair. Yet, mechanical failure leading to airway narrowing and collapse remains a major cause of morbidity and mortality following clinical implantation of decellularized tracheal grafts. To understand better the factors contributing to mechanical failure in vivo, we characterized the histo-mechanical properties of tracheas following two different decellularization protocols, including one that has been used clinically. All decellularized tracheas deviated from native mechanical behavior, which may provide insights into observed in vivo graft failures. We further analyzed protein content by western blot and analyzed microstructure by histological staining and found that the specific method of decellularization resulted in significant differences in the depletion of proteoglycans and degradation of collagens I, II, III, and elastin. Taken together, this work demonstrates that the heterogeneous architecture and mechanical behavior of the trachea is severely compromised by decellularization. Such structural deterioration may contribute to graft failure clinically and limit the potential of decellularized native tracheas as viable long-term orthotopic airway replacements.

Evaluating the biomechanical characteristics of cuffed-tracheostomy tubes using finite element analysis

The objective of this study was to perform finite element analysis (FEA) of cuff inflation within an anatomically accurate model of an adult trachea in four different cuffed-tracheostomy tube designs. The leakage quantified by the distance between the cuff and trachea was largest for the Tracoe cuff and smallest for the Portex cuff. The smooth muscle stresses were greatest for the Portex and least for the Distal cuff, respectively. The proposed FEA model offers a promising approach to virtually evaluate the sealing efficacy of cuffed-tracheostomy tubes and the tracheal wall stresses induced by cuff inflation, prior to application.

Mechanical Evaluation of Tracheal Grafts on Different Scales

Tissue engineered (or bioengineered) tracheas are alternative options under investigation when the resection with end-to-end anastomosis cannot be performed. One approach to develop bioengineered tracheas is a complex process that involves the use of decellularized tissue scaffolds, followed by recellularization in custom-made tracheal bioreactors. Tracheas withstand pressure variations and their biomechanics are of great importance so that they do not collapse during respiration, although there has been no preferred method of mechanical assay of tracheas among several laboratories over the years. These methods have been performed in segments or whole tracheas and in different species of mammals. This article aims to present some methods used by different research laboratories to evaluate the mechanics of tracheal grafts and presents the importance of the tracheal biomechanics in both macro and micro scales. If bioengineered tracheas become a reality in hospitals in the next few years, the standardization of biomechanical parameters will be necessary for greater consistency of results before transplantations.

Computational modeling of airway instability and collapse in tracheomalacia

Background

Tracheomalacia (TM) is a condition of excessive tracheal collapse during exhalation. Both acquired and congenital forms of TM are believed to result from morphological changes in cartilaginous, fibrous and/or smooth muscle tissues reducing airway mechanical properties to a degree that precipitates collapse. However, neither the specific amount of mechanical property reduction nor the malacic segment lengths leading to life threatening airway collapse in TM are known. Furthermore, the specific mechanism of collapse is still debated.

Methods

Computational nonlinear finite element models were developed to determine the effect of malacic segment length, tracheal diameter, and reduction in tissue nonlinear elastic properties on the risk for and mechanism of airway collapse. Cartilage, fibrous tissue, and smooth muscle nonlinear elastic properties were fit to experimental data from preterm lambs from the literature. These elastic properties were systematically reduced in the model to simulate TM.

Results

An intriguing finding was that sudden mechanical instability leading to complete airway collapse occurred in airways when even a 1 cm segment of cartilage and fibrous tissue properties had a critical reduction in material properties. In general, increased tracheal diameter, increased malacic segment length coupled with decreased nonlinear anterior cartilage/fibrous tissue nonlinear mechanical properties increased the risk of sudden airway collapse from snap through instability.

Conclusion

Modeling results support snap through instability as the mechanism for life threatening tracheomalacia specifically when cartilage ring nonlinear properties are reduced to a range between fibrous tissue nonlinear elastic properties (for larger diameter airways > 10 mm) to mucosa properties (for smaller diameter airways < 6 mm). Although reducing posterior tracheal smooth muscle properties to mucosa properties decreased exhalation area, no sudden instability leading to collapse was seen in these models.

Creation of an in vitro biomechanical model of the trachea using rapid prototyping

Previous in vitro models of the airways are either rigid or, if flexible, have not matched in vivo compliance characteristics. Rapid prototyping provides a quickly evolving approach that can be used to directly produce in vitro airway models using either rigid or flexible polymers. The objective of this study was to use rapid prototyping to directly produce a flexible hollow model that matches the biomechanical compliance of the trachea. The airway model consisted of a previously developed characteristic mouth-throat region, the trachea, and a portion of the main bronchi. Compliance of the tracheal region was known from a previous in vivo imaging study that reported cross-sectional areas over a range of internal pressures. The compliance of the tracheal region was matched to the in vivo data for a specific flexible resin by iteratively selecting the thicknesses and other dimensions of tracheal wall components. Seven iterative models were produced and illustrated highly non-linear expansion consisting of initial rapid size increase, a transition region, and continued slower size increase as pressure was increased. Thickness of the esophageal interface membrane and initial trachea indention were identified as key parameters with the final model correctly predicting all phases of expansion within a value of 5% of the in vivo data. Applications of the current biomechanical model are related to endotracheal intubation and include determination of effective mucus suctioning and evaluation of cuff sealing with respect to gases and secretions.

Shape and material characteristics of the trachea in the leatherback sea turtle promote progressive collapse and reinflation during dives

The leatherback turtle regularly undertakes deep dives and has been recorded attaining depths in excess of 1,200 m. Its trachea is an almost solid, elliptical-section tube of uncalcified hyaline cartilage with minimal connective tissue between successive rings. The structure appears to be advantageous for diving and perfectly designed for withstanding repeated collapse and reinflation. This study applies Boyle's law to the respiratory system (lungs, trachea and larynx) and estimates the changes in tracheal volume during a dive. These changes are subsequently compared with the results predicted by a corresponding finite element (FE) structural model, itself based on laboratory studies of the trachea of an adult turtle. Boyle's law predicts that the trachea will collapse progressively with greater volume change occurring in the early stages. The FE model reproduces the changes extremely well (agreeing closely with Boyle's law estimations) and provides visual representation of the deformed tracheal luminal area. Initially, the trachea compresses both ventrally and dorsally before levelling ventrally. Bulges are subsequently formed laterally and become more pronounced at deeper depths. The geometric configuration of the tracheal structure confers both homogeneity and strength upon it, which makes it extremely suited for enduring repeated collapse and re-expansion. The structure actually promotes collapse and is an adaptation to the turtle's natural environment in which large numbers of deep dives are performed annually.

FSI Analysis of a Healthy and a Stenotic Human Trachea Under Impedance-Based Boundary Conditions

In this work, a fluid-solid interaction (FSI) analysis of a healthy and a stenotic human trachea was studied to evaluate flow patterns, wall stresses, and deformations under physiological and pathological conditions. The two analyzed tracheal geometries, which include the first bifurcation after the carina, were obtained from computed tomography images of healthy and diseased patients, respectively. A finite element-based commercial software code was used to perform the simulations. The tracheal wall was modeled as a fiber reinforced hyperelastic solid material in which the anisotropy due to the orientation of the fibers was introduced. Impedance-based pressure waveforms were computed using a method developed for the cardiovascular system, where the resistance of the respiratory system was calculated taking into account the entire bronchial tree, modeled as binary fractal network. Intratracheal flow patterns and tracheal wall deformation were analyzed under different scenarios. The simulations show the possibility of predicting, with FSI computations, flow and wall behavior for healthy and pathological tracheas. The computational modeling procedure presented herein can be a useful tool capable of evaluating quantities that cannot be assessed in vivo, such as wall stresses, pressure drop, and flow patterns, and to derive parameters that could help clinical decisions and improve surgical outcomes.

Experimental characterization and constitutive modeling of the mechanical behavior of the human trachea

BACKGROUND AND AIMS

Cartilage and smooth muscle constitute the main structural components of the human central airways, their mechanical properties affect the flow in the trachea and contribute to the biological function of the respiratory system. The aim of this work is to find out the mechanical passive response of the principal constituents of the human trachea under static tensile conditions and to propose constitutive models to describe their behavior.

METHODS

Histological analyses to characterize the tissues and mechanical tests have been made on three human trachea specimens obtained from autopsies. Uniaxial tensile tests on cartilaginous rings and smooth muscle were performed. Tracheal cartilage was considered an elastic material and its Young's modulus and Poisson's coefficient were determined fitting the experimental curves using a Neo-Hookean model. The smooth muscle was proved to behave as a reinforced hyperelastic material with two families of collagen fibers, and its non-linearity was investigated using the Holzapfel strain-energy density function for two families of fibers to fit the experimental data obtained from longitudinal and transversal cuts.

RESULTS

For cartilage, fitting the experimental curves to an elastic model, a Young's modulus of 3.33 MPa and nu=0.49 were obtained. For smooth muscle, several parameters of the Holzapfel function were found out (C(10)=0.877 kPa, k(1)=0.154 kPa, k(2)=34.157, k(3)=0.347 kPa and k(4)=13.889) and demonstrated that the tracheal muscle was stiffer in the longitudinal direction.

CONCLUSION

The better understanding of how these tissues mechanically behave is essential for a correct modeling of the human trachea, a better simulation of its response under different loading conditions, and the development of strategies for the design of new endotracheal prostheses.

Tracheal compliance and limit flow rate changes in a murine model of asthma

Trachea is the unique passage for air to flow in and out. Its tone is of importance for the respiration system. However, investigation on how tracheal tone changes due to asthma is limited. Aiming at studying how the mechanical property changes due to asthma as well as the compliance and flow limitation, the following methods are adopted. Static and passive pressure-volume tests of rats' trachea of the asthmatic and control groups are carried out and a new type of tube law is formulated to fit the experimental data, based on which changes of compliance and limit flow rate are investigated. In order to give explanation to such changes, histological examinations with tracheal soft tissues are made. The results show that compliance, limit flow rate and material constants included in the tube law largely depend on the longitudinal stretching ratio. Compared with the control group, the tracheal compliance of asthmatic animals decreases significantly, which results in an increased limit flow rate. Histological studies indicate that asthma can lead to hyperplasia/hypertrophy of smooth muscle cells, and increase elastin and collagen fibres in the muscular membrane. Though decreasing compliance increases stability, during the onset of asthma, limit flow rate is much smaller due to the lower transmural pressure. Asthma leads to a stiffer trachea and the obtained results reveal some aspects relevant to asthma-induced tracheal remodelling.

Effect of Repeated Induced Airway Collapse During Total Liquid Ventilation

Negative pressure generated during the expiratory phase of total liquid ventilation (TLV) may induce airway collapse. Evaluation of the effect of repeated airway collapse is crucial to optimize this technique. A total of 24 New Zealand White rabbits were randomly divided into four groups. Ventilation was performed for 6 hours with different strategies: conventional gas ventilation, TLV without airway collapse, and TLV with collapse induced in either 75 or 150 sequential breaths. In the treated groups, airway collapse was induced by increasing the perfluorocarbon drainage velocity while maintaining the minute ventilation constant. Airway pressure, gas exchange, and blood pressure were monitored at 30-minute intervals. At the end of the experiment, airway and lung parenchyma specimens were processed for light microscopy. No evidence of fluorothorax was noticed in any of the four groups at autopsy examination. Minimal signs of inflammation were noticed in all airway and lung parenchyma specimens, but no evident structural alteration was visible. Adequate gas exchange and systemic blood pressure were maintained during all the studies. Repeated airway collapse is not associated with structural changes in the respiratory system and does not alter the gas exchange ability of the lungs.

An image-based computational model of oscillatory flow in the proximal part of tracheobronchial trees

A computational model of an oscillatory laminar flow of an incompressible Newtonian fluid has been carried out in the proximal part of human tracheobronchial trees, either normal or with a strongly stenosed right main bronchus. After acquisition with a multislice spiral CT, the thoracic images are processed to reconstruct the geometry of the trachea and the first six bronchus generations and to virtually travel inside this duct network. The facetisation associated with the 3D reconstruction of the tracheobronchial tree is improved to get a computation-adapted surface triangulation, which leads to a volumic mesh composed of tetrahedra. The Navier-Stokes equations associated with the classical boundary conditions and different values of the flow dimensionless parameters are solved using the finite element method. The airways are supposed to be rigid during rest breathing. The flow distribution among the set of bronchi is determined during the respiratory cycle. Cycle reproducibility and mesh size effects on the numerical results are examined. Helpful qualitative data are provided rather than accurate quantitative results in the context of multimodelling, from image processing to numerical simulations.

Structure and biomechanical properties of the trachea of the striped dolphinStenella coeruleoalba: Evidence for evolutionary adaptations to diving

This study analyzes the structure and mechanical properties of the trachea of the striped dolphin Stenella coeruleoalba, one of the most common cetacean species. The cetacean trachea is made up of closed or semiclosed cartilaginous rings without a paries membranaceus. Our results indicate that the inner lining of the trachea contains erectile tissue in which several venous lacunae permeate the mucosa. We also observed and described the presence of peripheral neurons containing nitric oxide along the rim of the venous lacunae. Data obtained from compression and tensile tests and comparison with the pig and goat tracheas indicate a higher stiffness and a different, higher breaking point for the dolphin trachea. On the whole, our data suggest that the trachea of the striped dolphin possesses structural properties that allow rapid filling with blood, possibly in relation to dive activities, and also allow modifications due to increased pressure and immediate return to the original shape without risks of permanent bending or rupture, as would happen in a terrestrial mammal. As the organ undergoes intense pressure difference during descent to optimal foraging depth and subsequent rapid ascent to surface, especially in deep dives of hundreds of meters, the specific structural and biomechanical peculiarities of the trachea of the striped dolphin may represent an evolutionary adaptation to life in the water and to diving.

A numerical and experimental study of compliance and collapsibility of preterm lamb tracheae

Knowledge of the mechanical behaviour of immature tracheae is crucial in order to understand the effects exerted on central airways by ventilatory treatments, particularly of Total Liquid Ventilation. In this study, a combined experimental and computational approach was adopted to investigate the compliance and particularly collapsibility of preterm lamb tracheae in the range of pressure likely applied during Total Liquid Ventilation (-30 to 30 cmH2O). Tracheal samples of preterm lambs (n = 5; gestational age 120-130 days) were tested by altering transmural pressure from -30 to 30 cmH2O. Inflation (Si) and collapsing (Sc) compliance values were calculated in the ranges 0 to 10 cmH2O and -10 to 0 cmH2O, respectively. During the tests, an asymmetric behaviour of the DeltaV/V0 vs. P curves at positive and negative pressure was observed, with mean Si = 0.013 cmH2O(-1) and Sc = 0.053 cmH2O(-1). A different deformed configuration of the sample regions was observed, depending on the posterior shape of cartilaginous ring. A three-dimensional finite-element structural model of a single tracheal ring, based on histology measurements of the tested samples was developed. The model was parameterised in order to represent rings belonging to three different tracheal regions (craniad, median, caudal) and numerical analyses replicating the collapse test conditions were performed to evaluate the ring collapsibility at pressures between 0 and -30 cmH2O. Simulation results were compared to experimental data to verify the model's reliability. The best model predictions occurred at pressures -30 to -10 cmH2O. In this range, a model composed of median rings best interpreted the experimental data, with a maximum error of 2.7%; a model composed of an equal combination of all rings yielded an error of 12.6%.

Mechanical properties of artificial tracheas composed of a mesh cylinder and a spiral stent

Much work has been done on the materials used for mesh-type artificial tracheas, but a precise mechanical evaluation of these structures has not yet been performed. In the present study, we determined the mechanical properties of typical mesh-type artificial tracheas and compared them with those of native trachea. Four types of artificial trachea were made and used for the mechanical tests. The basic frame of all the specimens was composed of a mesh cylinder and a spiral stent. The specimen whose mesh was sealed with collagen sponge showed almost the same behavior in the force-strain curve under compression, suggesting that collagen sealing has little effect on mechanical properties. Agreement between measured and estimated mechanical properties was good, especially in the low strain region, suggesting that artificial tracheas can be designed in terms of mechanical properties by mainly considering the basic frame structure.

Analysis of tracheal mechanics and applications

We have developed a mathematical model for a tracheal ring that consists of a "horseshoe" of cartilage with its tips joined by a membrane. The ring is subjected to a uniform transmural pressure (Ptm) difference. The model was used to calculate the cross-sectional area (A) of the trachea. Whereas the mechanics of the deformation of the cartilage were analyzed using elastica theory, the posterior membrane was treated as a simple membrane that is inextensible under changes in Ptm. The membrane can be specified to be of any length less than baseline and thus can represent a posterior membrane under tension. The cartilage can have specifiable nonuniform unstressed curvature as well as nonuniform bending stiffness. We have investigated the effect on the tracheal A-Ptm curve of posterior membrane length and tensile force in the membrane, cartilage shape and elasticity, and localized weakening of the cartilage. The model predictions are in good agreement with magnetic resonance imaging data from rabbit tracheas and show that the shape of the horseshoe as well as the posterior membrane force are important determinants of tracheal compliance.

Collapse of diseased arteries with eccentric cross section

The mechanism of atherosclerotic plaque rupture is not fully understood. Mechanical stress may be one of the factors contributing to the instability of the fibrous plaque cap. The existence of a severe stenosis may lower the transmural pressure enough to cause the collapse of arteries leading to high concentrated compressive and tensile stresses. This study presents quantitative estimates of the stresses and deformations in the collapsed thick-walled artery. The results of large deformation finite element calculations identify the locations of the high stress concentrations and their magnitudes which cannot be precisely predicted under a thin-wall assumption. The maximum compressive stress calculated reached 80% of the Young's modulus for fairly small negative transmural pressures. Results are useful to predict likely locations of the plaque cap rupture due to compressive stresses. The tube law of area as a function of transmural pressure showed a large discrepancy from a thin-wall calculation. The buckling pressure calculated for the outer-to-inner wall surface radius ratio of 1.60 was nearly 30% lower than that of the elastic thin-wall buckling theory. An increase in eccentricity further reduced this buckling pressure. The results indicate that a thick plaque which is eccentric increases the likelihood of collapse of stenotic arteries.