Effect of freezing storage time on the elastic and viscous properties of the porcine TMJ disc.

The correct characterisation of the articular disc of the temporomandibular joint (TMJ) is key to study the masticatory biomechanics. For the interval from extraction until testing, freezing is the most used preservation technique for biological tissues, but its influence on their behaviour is still unclear. An important error can be committed in the characterisation of such tissues if freezing has any effect on their mechanical properties. Thus, the aim of this study was to determine whether the freezing storage time causes any change in the mechanical properties of the TMJ discs. To check that, the specimens were stored in a -20°C freezer during different time intervals: 1 day, 1 week, 1 month and 3 months. Fresh specimens, tested right after extraction, were used as the control group. Compressive stress relaxation tests were carried out on the specimens and a quasi-linear viscoelastic (QLV) model was used to fit the experimental curves. A statistical analysis detected significant differences among the groups. Post-hoc tests determined that freezing the specimens more than 30 days may lead to changes in the viscoelastic properties of the tissue.

Finite element analysis of the human mastication cycle.

The aim of this paper is to propose a biomechanical model that could serve as a tool to overcome some difficulties encountered in experimental studies of the mandible. One of these difficulties is the inaccessibility of the temporomandibular joint (TMJ) and the lateral pterygoid muscle. The focus of this model is to study the stresses in the joint and the influence of the lateral pterygoid muscle on the mandible movement. A finite element model of the mandible, including the TMJ, was built to simulate the process of unilateral mastication. Different activation patterns of the left and right pterygoid muscles were tried. The maximum stresses in the articular disc and in the whole mandible during a complete mastication cycle were reached during the instant of centric occlusion. The simulations show a great influence of the coordination of the right and left lateral pterygoid muscles on the movement of the jaw during mastication. An asynchronous activation of the lateral pterygoid muscles is needed to achieve a normal movement of the jaw during mastication.

A study of the temporomandibular joint during bruxism.

A finite element model of the temporomandibular joint (TMJ) and the human mandible was fabricated to study the effect of abnormal loading, such as awake and asleep bruxism, on the articular disc. A quasilinear viscoelastic model was used to simulate the behaviour of the disc. The viscoelastic nature of this tissue is shown to be an important factor when sustained (awake bruxism) or cyclic loading (sleep bruxism) is simulated. From the comparison of the two types of bruxism, it was seen that sustained clenching is the most detrimental activity for the TMJ disc, producing an overload that could lead to severe damage of this tissue.

Effect of non-uniform thickness of samples in stress relaxation tests under unconfined compression of samples of articular discs.

A precise information of the biomechanical properties of soft tissues is required to develop a suitable simulation model, with which the distribution of stress and strain in the complex structures can be estimated. Many soft tissues have been mechanically characterized by stress relaxation tests under unconfined or confined compression. In general, full-thickness samples are extracted to reduce the damage in the tissue as much as possible. However, it is not guaranteed that these samples have a uniform thickness or, in other words, planar parallel faces. In particular, in the articular disc of the temporomandibular joint, many studies can be found testing full-thickness samples for which that thickness is known to be non-uniform, while making the assumption of uniaxial stress state to extract the mechanical properties from those tests. That inaccuracy may have a strong influence in some cases and needs a profound revision. The main goal of this work is to quantify the error committed in that assumption and the influence of the variation of thickness on that error in a particular test: stress relaxation tests under unconfined compression. Based on this error and defining an allowable tolerance, a criterion is established to reject samples depending on their aspect ratio.

Numerical simulation of a relaxation test designed to fit a quasi-linear viscoelastic model for temporomandibular joint discs.

The main objectives of this work are: (a) to introduce an algorithm for adjusting the quasi-linear viscoelastic model to fit a material using a stress relaxation test and (b) to validate a protocol for performing such tests in temporomandibular joint discs. This algorithm is intended for fitting the Prony series coefficients and the hyperelastic constants of the quasi-linear viscoelastic model by considering that the relaxation test is performed with an initial ramp loading at a certain rate. This algorithm was validated before being applied to achieve the second objective. Generally, the complete three-dimensional formulation of the quasi-linear viscoelastic model is very complex. Therefore, it is necessary to design an experimental test to ensure a simple stress state, such as uniaxial compression to facilitate obtaining the viscoelastic properties. This work provides some recommendations about the experimental setup, which are important to follow, as an inadequate setup could produce a stress state far from uniaxial, thus, distorting the material constants determined from the experiment. The test considered is a stress relaxation test using unconfined compression performed in cylindrical specimens extracted from temporomandibular joint discs. To validate the experimental protocol, the test was numerically simulated using finite-element modelling. The disc was arbitrarily assigned a set of quasi-linear viscoelastic constants (c1) in the finite-element model. Another set of constants (c2) was obtained by fitting the results of the simulated test with the proposed algorithm. The deviation of constants c2 from constants c1 measures how far the stresses are from the uniaxial state. The effects of the following features of the experimental setup on this deviation have been analysed: (a) the friction coefficient between the compression plates and the specimen (which should be as low as possible); (b) the portion of the specimen glued to the compression plates (smaller areas glued are better); and (c) the variation in the thickness of the specimen. The specimen's faces should be parallel to ensure a uniaxial stress state. However, this is not possible in real specimens, and a criterion must be defined to accept the specimen in terms of the specimen's thickness variation and the deviation of the fitted constants arising from such a variation.