FSI analysis of a human trachea before and after prosthesis implantation.

In this work we analyzed the response of a stenotic trachea after a stent implantation. An endotracheal stent is the common treatment for tracheal diseases such as stenosis, chronic cough, or dispnoea episodes. Medical treatment and surgical techniques are still challenging due to the difficulties in overcoming potential complications after prosthesis implantation. A finite element model of a diseased and stented trachea was developed starting from a patient specific computerized tomography (CT) scan. The tracheal wall was modeled as a fiber reinforced hyperelastic material in which we modeled the anisotropy due to the orientation of the collagen fibers. Deformations of the tracheal cartilage rings and of the muscular membrane, as well as the maximum principal stresses, are analyzed using a fluid solid interaction (FSI) approach. For this reason, as boundary conditions, impedance-based pressure waveforms were computed modeling the nonreconstructed vessels as a binary fractal network. The results showed that the presence of the stent prevents tracheal muscle deflections and indicated a local recirculatory flow on the stent top surface which may play a role in the process of mucous accumulation. The present work gives new insight into clinical procedures, predicting their mechanical consequences. This tool could be used in the future as preoperative planning software to help the thoracic surgeons in deciding the optimal prosthesis type as well as its size and positioning.

FSI analysis of the coughing mechanism in a human trachea.

The main physiological function of coughing is to remove from the airways the mucus and foreign particles that enter the lungs with respirable air. However, in patients with endotracheal tubes, further surgery has to be performed to improve cough effectiveness. Thus, it is necessary to analyze how this process is carried out in healthy tracheas to suggest ways to improve its efficacy in operated patients. A finite element model of a human trachea is developed and used to analyze the deformability of the tracheal walls under coughing. The geometry of the trachea is obtained from CT of a 70-year-old male patient. A fluid structure interaction approach is used to analyze the deformation of the wall when the fluid (in this case, air) flows inside the trachea. A structured hexahedral-based grid for the tracheal walls and an unstructured tetrahedral-based mesh with coincident nodes for the fluid are used to perform the simulations with the finite element-based commercial software code (ADINA R&D Inc.). Tracheal wall is modeled as an anisotropic fiber reinforced hyperelastic solid material in which the different orientation of the fibers is introduced. The implantation of an endotracheal prosthesis is simulated. Boundary conditions for breathing and coughing are applied at the inlet and at the outlet surfaces of the fluid mesh. The collapsibility of a human trachea under breathing and coughing is shown in terms of flow patterns and wall stresses. The ability of the model to reproduce the normal breathing and coughing is proved by comparing the deformed shape of the trachea with experimental results. Moreover the implantation of an endotracheal prosthesis would be related with a decrease of coughing efficiency, as clinically seen.

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.

A finite element model to accurately predict real deformations of the breast.

Most surgical procedures in breast plastic surgery are either reconstructive procedures following oncologic interventions (tumorectomy, quadrantectomy, mastectomy ...) or aesthetic ones, both augmentation and reduction. With current techniques, the results of such procedures cannot be fully guaranteed. Usually, surgical planning is based on a photographic and anthropometric study of the breast only. Among others, one of the difficulties that the plastic surgeons have is the noticeable change of the breast shape with the position of the patient. Thus, it is more and more necessary to plan a presurgical methodology to help the plastic surgeon and guarantee the patient a successful result of the intervention. In order to establish a reliable simulation method that could predict a patient-specific outcome after breast surgery, this study started trying to correlate spatial features of the breast between lying and standing up positions. A biomechanical model of breast was proposed and implemented into a finite element context to predict deformations, and from these the breast shape in different positions. The resulting shapes were compared with multimodal images, whereas the breast surface displacements were compared with manually identified landmarks and 3D scanner images. From the results, it can be concluded that the model hereby presented reasonably approximates breast response to gravity forces, therefore providing accurate and useful information to the surgeon planning such surgical procedures.

An accurate finite element model of the cervical spine under quasi-static loading.

Cervical disc injury due to impact has been observed in clinical and biomechanical investigations; however, there is a lack of data that helps to elucidate the mechanisms of disc injury during these collisions. Therefore, it is necessary to understand the behavior of the cervical spine under different types of loading situations. A three dimensional finite element (FE) model for the multi-level cervical spine segment (C0-C7) was developed using computed tomography (CT) data and applied to study the internal stresses and strains of the intervertebral discs under quasi-static loading conditions. The intervertebral discs were treated as nonlinear, anisotropic and incompressible subjected to large deformations. The model accuracy was validated by comparing it with previously published experimental and numerical results for different movements. It was shown that the use of a fiber reinforced model to describe the behavior of the annulus of the discs would predict higher maximum shear strains than an isotropic one, being therefore important the use of complex constitutive models in order to be able to detect the appearance of injured zones, since those strains and stresses are supposed to be related with damage to soft tissues. Several movements were analyzed: flexion, extension and axial rotation, obtaining that the maximum shear stresses in the disc were higher for a flexo-extension movement.

3D finite element simulation of the opening movement of the mandible in healthy and pathologic situations.

One of the essential causes of disk disorders is the pathologic change in the ligamentous attachments of the disk-condyle complex. In this paper, the response of the soft components of a human temporomandibular joint during mouth opening in healthy and two pathologic situations was studied. A three-dimensional finite element model of this joint comprising the bone components, the articular disk, and the temporomandibular ligaments was developed from a set of medical images. A fiber reinforced porohyperelastic model was used to simulate the behavior of the articular disk, taking into account the orientation of the fibers in each zone of this cartilage component. The condylar movements during jaw opening were introduced as the loading condition in the analysis. In the healthy joint, it was obtained that the highest stresses were located at the lateral part of the intermediate zone of the disk. In this case, the collateral ligaments were subject to high loads, since they are responsible of the attachment of the disk to the condyle during the movement of the mandible. Additionally, two pathologic situations were simulated: damage of the retrodiscal tissue and disruption of the lateral discal ligament. In both cases, the highest stresses moved to the posterior part of the disk since it was displaced in the anterior-medial direction. In conclusion, in the healthy joint, the highest stresses were located in the lateral zone of the disk where perforations are found most often in the clinical experience. On the other hand, the results obtained in the damaged joints suggested that the disruption of the disk attachments may cause an anterior displacement of the disk and instability of the joint.