Approach towards the porous fibrous structure of the periodontal ligament using micro-computerized tomography and finite element analysis.

The periodontal ligament (PDL) is a porous and fibrous soft tissue situated around the tooth, which plays a key role in the transmission of loads from the tooth to the alveolar bone of the mandible. Although several studies have tried to characterize its mechanical properties, the behaviour of this tissue is not clear yet. In this study, a new simulation methodology based on a material model which considers the contribution of porous and fibrous structure with different material model formulations depending on the effort direction is proposed. The defined material model was characterized by a non-linear approximation of the porous fibrous matrix to experimental results obtained from samples of similar species and was validated by rigorous test simulations under tensile and compressive loads. The global PDL response was also validated using the parameters of the characterization in a finite element model of full human canine tooth obtained by micro-tomography. The results suggest that the porous contribution has high influence during compression because the bulk modulus of the material depends on the ability of interstitial fluid to drain. On the other hand, the collagen fibres running along the load direction are the main responsible of the ligament stiffness during tensile efforts. Thus, a material model with distinct responses depending of the load direction is proposed. Furthermore, the results suggest the importance of considering 3D finite element models based of the real morphology of human PDL for representing the irregular stress distribution caused by the coupling of complex material models and irregular morphologies.

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.

Patient-specific models of human trachea to predict mechanical consequences of endoprosthesis implantation.

Nowadays, interventions associated with the implantation of tracheal prostheses in patients with airway pathologies are very common. This surgery may promote problems such as migration of the prosthesis, development of granulation tissue at the edges of the stent with overgrowth of the tracheal lumen or accumulation of secretions inside the prosthesis. Among the movements that the trachea carries out, swallowing seems to have harmful consequences for the tracheal tissues surrounding the prosthesis. In this work, a finite-element-based tool is presented to construct patient-specific tracheal models, introducing the endotracheal prosthesis and analysing the mechanical consequences of this surgery during swallowing. A complete description of a patient-specific tracheal model is given, and a fully experimental characterization of the tracheal tissues is presented. To construct patient-specific grids, a mesh adaptation algorithm has been developed and the implantation of a tracheal prosthesis is simulated. The ascending deglutition movement of the trachea is recorded using real data from each specific patient from fluoroscopic images before and after implantation. The overall behaviour of the trachea is modified when a prosthesis is introduced. The presented tool has been particularized for two different patients (patient A and patient B), allowing prediction of the consequences of this kind of surgery. In particular, patient A had a decrease of almost 30 per cent in his ability to swallow, and an increase in stresses that were three times higher after prosthesis implantation. In contrast, patient B, who had a shorter trachea and who seemed to undergo more damaging effects, did not have a significant reduction in his ability to swallow and did not present an increase in stress in the tissues. In both cases, there are clinical studies that validate our results: namely, patient A underwent a further intervention whereas the outcome of patient B's surgery was completely successful. Notwithstanding the fact that there are a lot of uncertainties relating to the implantation of endotracheal prostheses, the present work gives a new insight into these procedures, predicting their mechanical consequences. This tool could be used in the future as pre-operative planning software to help thoracic surgeons in deciding the optimal prosthesis as well as its size and positioning.

An accurate validation of a computational model of a human lumbosacral segment.

Clinical studies have recently documented that there is sufficient evidence to suggest that abnormal motion may be an indicator of abnormal mechanics of the spine and, therefore, may be associated with some types of low-back pain. However, designating a motion as abnormal requires knowledge of normal motions. This work hence aims to develop an accurate computational model to simulate the bio-mechanical response of the whole lumbosacral spinal unit (L1-S1) under physiological loadings and constraint conditions. In order to meet this objective, computed tomography (CT) scanning protocols, finite element (FE) analysis and accurate constitutive modelling have been integrated. Then the ranges of motion (ROM) under flexion, extension and lateral bending moment were measured and compared with experimental data, finding an excellent agreement. In particular, the ability of the model to reproduce the relative rotation between each couple of vertebrae was proved. Finally, the shear stresses for the most extreme load cases were reported in order to predict which are the most risky conditions and where the maximum damage would be located. The results indicate that the greater values of the stresses were located at L4-S1 levels just in the interfaces between disc and vertebrae across the posterior and posterolateral zone. This result can be clinically correlated with the existence of damage exactly where the stresses were maximal in the proposed finite element model.

An experimental study of the mouse skin behaviour: damage and inelastic aspects.

Samples of male and female mice skin were tested under monotonic and cyclic loading to mechanically characterize the tissue for large deformations. Cyclic tests have shown a typical Mullins effect widely known for elastomers and other soft tissues. No statistical difference was found in the maximum stretch of the sample after the fifth loading cycle for male (1.26 +/- 0.035) and female (1.18 +/- 0.083). However, larger dispersion was obtained for the maximum stress for both genders, 0.61 +/- 0.16 MPa for male and 0.78 +/- 0.32 MPa for female. Results show the presence of inelastic strain and stress softening in the skin at large deformations. They also have shown how stress softening and residual strain change with the magnitude of the applied load. Good correlation was observed between the residual strain and the maximum strain previously attained by the sample during loading for all samples. However, the correlation was different between genders.

Dynamic 3D FE modelling of the human temporomandibular joint during whiplash.

Rear-end impacts account for more than one-third of vehicle accidents, and nearly 40% of these accidents produce whiplash injuries. Whiplash injury to the neck has often been considered a significant risk factor for the development of temporomandibular disorders (TMD). The objective of this study was to simulate the dynamic response of the temporomandibular joint during two types of impacts: a rear end and a frontal impact. To understand the dynamic forces acting on the joint, we extended a previous human temporomandibular joint model and analyzed the stress distributions in the soft elements of the joint. In the rear-end impact, it could be appreciated that the inertia of the mandible caused it to move posteriorly slower than the head, and this resulted in downward and forward displacements of the disc-condyle complex relative to the cranial base. Consequently, a rapid and big mouth opening occurs. In contrast, during the frontal impact, the mouth hardly opened, because the superior maxilla pushed the mandible to move together. There was not differential movement between bony components of the joint and therefore the soft tissues of the joint were not subjected to high loads. From these results, and despite the limitations of the simulations performed, we could conclude that neither a rear-end impact at low-velocity nor a frontal impact would produce damage to the soft tissues of the joint.

Influence of unilateral disc displacement on the stress response of the temporomandibular joint discs during opening and mastication.

The temporomandibular joint plays a crucial role in human mastication acting as a guide of jaw movements. During these movements, the joint is subjected to loads which cause stresses and deformations in its cartilaginous structures. A perfect balance between the two sides of the joint is essential to maintain the physiological stress level within the tissues. Therefore, it has been suggested that a derangement of the joint is a contributing factor in the development of mandibular asymmetry, especially if problems of the temporomandibular joint start in childhood or adolescence. To analyze the movement of the mandible and the stresses undergone by the discs, two finite element models of the human temporomandibular joint including the masticatory system were developed, one corresponding to a healthy joint and the other with a unilateral anterior disc displacement with their movement controlled by muscle activation. A fibre-reinforced porohyperelastic model was used to simulate the behaviour of the articular discs. The stress distribution was analyzed in both models during free opening and closing, and during the introduction of a resistant force between incisors or molars. It was found that a slight unilateral anterior disc displacement does not lead to mandibular asymmetry but to a slight decrease of the maximum gape. With the introduction of a restriction between incisors, the maximum stresses moved to the anterior band in contrast to what happened if the restriction was imposed between molars where maximum stresses were located more posteriorly. Finally, the presence of a unilateral displacement of the disc involved a strong change in the overall behaviour of the joint including also the healthy side, where the maximum stresses moved to the posterior part.

An accurate simulation model of anteriorly displaced TMJ discs with and without reduction.

Internal derangement of the temporomandibular joint (TMJ) is defined as an abnormal positional relationship of the disc relative to the mandibular condyle and the glenoid fossa. Among others, the anterior disc displacement is the most common disorder, however its origin and consequences are still unclear. Several finite element simulations of the TMJ have been developed, but none of them has reported dynamic simulations of the disc as a three-dimensional, fiber-reinforced biphasic material under finite deformations, during the opening movement of a pathologic joint affected of an anterior displacement of the disc with and without reduction, using a realistic geometry of the ligaments in the joint. The aim of the work presented here was to compare the stress distribution in the healthy joint and in two pathologic situations, one joint affected of an anterior disc displacement with reduction (ADDWR) and one without reduction (ADDWOR) during an opening movement of the mouth. It was found that, while in the healthy disc the highest compressive stresses were located in the intermediate zone, in the pathologic joints the maximum compressive stresses were located in the posterior band both in the ADDWOR case and in the ADDWR before the reduction. Moreover, although the final stress distribution in the ADDWR was similar to that in the healthy case, the collateral ligaments supported higher stresses, a fact that could lead to degeneration of these components and subsequently to the total anterior displacement of the disc. Finally, the results suggest that an anterior displacement of the disc would lead to higher compressive and tangential stresses in the posterior band of the disc than in the healthy one, and as a consequence, to possible perforations in that zone of the disc which would modify its geometry if no treatment is applied.

Anterior displacement of the TMJ disk: repositioning of the disk using a Mitek system. A 3D finite element study.

In this paper the behaviors of the temporomandibular joint (TMJ) with an anteriorly displaced disk without reduction and with a surgically repositioned one were compared with the response of a healthy disk during jaw opening. The movement of each joint was obtained imposing the same opening path between incisors and assuming that the movement of the condyle is determined by the passive action of the masticatory muscles and the restrictions imposed by the articulating surfaces and the ligaments. A fiber-reinforced porohyperelastic model was used to simulate the behavior of the articular disk. The influence of the friction coefficient in the diseased joint was also analyzed, finding that the final displacement of the complex condyle-disk was smaller as the friction coefficient increased. On the other hand, its displacement in the repositioned joint was different than in the healthy case because the artificial sutures used in the surgery do not fully stabilize the disk posteriorly as the retrodiscal tissue does. The stress response of the disk changed in both pathologic cases: in the displaced joint the highest stresses moved from the intermediate zone (healthy case) to the posterior band, and in the reconstructed one the most loaded zone moved posteriorly at total opening. Besides, local stress concentrations appeared in the neighborhood of the artificial sutures and therefore damage of the disk and releasing of the sutures might be possible postoperatively.

The effect of collagen reinforcement in the behaviour of the temporomandibular joint disc.

In this paper, the influence of collagen fibres in the behaviour of the temporomandibular joint disc is studied. A three-dimensional finite element model of the joint is developed from a set of medical images. The model comprises the mandible, part of the cranium and both temporomandibular joints. Joints have been considered to be composed of the articular discs and the temporomandibular ligaments. A fibre-reinforced porohyperelastic model was used to study the response under clenching of the fibrocartilage that composes the articular disc. This was divided in an intermediate zone, and two bands, an anterior and other posterior, in order to define the orientation of collagen fibres. The study demonstrates that the introduction of collagen fibres in the biphasic behaviour of the articular disc implies for a prescribed displacement not only an increase of the pressurization in the tissue, but also higher stresses in the anterior and posterior bands, as well as in the lateral zone of the disc. Thus, modelling the disc as an isotropic solid matrix leads in this case to an overestimation of the stresses in the intermediate zone, an underestimation of the pore pressure in this area, and an underestimation of the stresses in the rest of the disc.

Finite element analysis of the temporomandibular joint during lateral excursions of the mandible.

One of the most significant characteristics of the temporomandibular joint (TMJ) is that it is in fact composed of two joints. Several finite element simulations of the TMJ have been developed but none of them analysed the different responses of its two sides during nonsymmetrical movement. In this paper, a lateral excursion of the mandible was introduced and the biomechanical behaviour of both sides was studied. A three-dimensional finite element model of the joint comprising the bone components, both articular discs, and the temporomandibular ligaments was used. A fibre-reinforced porohyperelastic model was introduced to simulate the behaviour of the articular discs, taking into account the orientation of the fibres in each zone of these cartilage components. The mandible movement during its lateral excursion was introduced as the loading condition in the analysis. As a consequence of the movement asymmetry, the discs were subjected to different load distributions. It was observed that the maximal shear stresses were located in the lateral zone of both discs and that the lateral attachment of the ipsilateral condyle-disc complex suffered a large distortion, due to the compression of this disc against the inferior surface of the temporal bone. These results may be related with possible consequences of a common disorder called bruxism. Although it would be necessary to perform an exhaustive analysis of this disorder, including the contact forces between the teeth during grinding, it could be suggested that a continuous lateral movement of the jaw may lead to perforations of both discs in their lateral part and may damage the lateral attachments of the disc to the condyle.