In silico study of cuspid' periodontal ligament damage under parafunctional and traumatic conditions of whole-mouth occlusions. A patient-specific evaluation

Finite element models: A road to in-silico modeling in the age of personalized dentistry

OBJECTIVE

This article reviews the applications of Finite Element Models (FEMs) in personalized dentistry, focusing on treatment planning, material selection, and CAD-CAM processes. It also discusses the challenges and future directions of using finite element analysis (FEA) in dental care.

DATA

This study synthesizes current literature and case studies on FEMs in personalized dentistry, analyzing research articles, clinical reports, and technical papers on the application of FEA in dental biomechanics.

SOURCES

Sources for this review include peer-reviewed journals, academic publications, clinical case studies, and technical papers on dental biomechanics and finite element analysis. Key databases such as PubMed, Scopus, Embase, and ArXiv were used to identify relevant studies.

STUDY SELECTION

Studies were selected based on their relevance to the application of FEMs in personalized dentistry. Inclusion criteria were studies that discussed the use of FEA in treatment planning, material selection, and CAD-CAM processes in dentistry. Exclusion criteria included studies that did not focus on personalized dental treatments or did not utilize FEMs as a primary tool.

CONCLUSIONS

FEMs are essential for personalized dentistry, offering a versatile platform for in-silico dental biomechanics modeling. They can help predict biomechanical behavior, optimize treatment outcomes, and minimize clinical complications. Despite needing further advancements, FEMs could help significantly enhance treatment precision and efficacy in personalized dental care.

CLINICAL SIGNIFICANCE

FEMs in personalized dentistry hold the potential to significantly improve treatment precision and efficacy, optimizing outcomes and reducing complications. Their integration underscores the need for interdisciplinary collaboration and advancements in computational techniques to enhance personalized dental care.

Inverse finite element analysis for an axisymmetric model of vertical tooth extraction

BACKGROUND AND OBJECTIVE

Tooth extraction is a common clinical procedure with biomechanical factors that can directly influence patient outcomes. Recent development in atraumatic extraction techniques have endeavoured to improve treatment outcomes, but the characterization of extraction biomechanics is sparse. An axisymmetric inverse finite element (FE) approach is presented to represent the biomechanics of vertical atraumatic tooth extraction in an ex-vivo swine model.

METHODS

Geometry and boundary conditions from the model are determined to match the extraction of swine incisors in a self-aligning ex vivo extraction experiment. Material parameters for the periodontal ligament (PDL) model are determined by solving an inverse FE problem using clusters of data obtained from 10 highly-controlled mechanical experiments. A seven-parameter visco-hyperelastic damage model, based on an Arruda-Boyce framework, is used for curve fitting. Three loading schemes were fit to obtain a common set of material parameters.

RESULTS

The inverse FE results demonstrate good predictions for overall force-time curve shape, peak force, and time to peak force. The fit model parameters are sufficiently consistent across all three cases that a coefficient-averaged model was taken that compares well to all three cases. Notably, the initial modulus ,u, converged across trials to an average value of 0.472 MPa with an average viscoelastic constant g of 0.561.

CONCLUSIONS

The presented model is found to have consistent parameters across loading cases. The capability of this model to represent the fundamental mechanical characteristics of the dental complex during vertical extraction loading is a significant advancement in the modelling of extraction procedures. Future work will focus on verifying the model as a predictive design tool for assessing new loading schemes in addition to investigating its applications to subject-specific problems.

A Finite Element Analysis Study of Edentulous Model with Complete Denture to Simulate Masticatory Movement

The purposes of this study are to establish and validate a finite element (FE) model using finite element analysis methods and to identify optimal loading conditions to simulate masticatory movement. A three-dimensional FE model of the maxillary and mandibular cortical bone, cancellous bone, and gingiva was constructed based on edentulous cone-beam-computed tomography data. Dental computer-aided design software was used to design the denture base and artificial teeth to produce a complete denture. Mesh convergence was performed to derive the optimal mesh size, and validation was conducted through comparison with mechanical test results. The mandible was rotated step-by-step to induce movements similar to actual mastication. Results showed that there was less than a 6% difference between the mechanical test and the alveolar bone-complete denture. It opened 10° as set in the first stage, confirming that the mouth closed 7° in the second stage. Occlusal contact occurred between the upper and lower artificial teeth as the mouth closed the remaining angle of 3° in the third stage while activating the masseter muscle. These results indicate that the FE model and masticatory loading conditions developed in this study can be applied to analyze biomechanical effects according to the wearing of dentures with various design elements applied.

Effect of Sparc knockout on the extracellular matrix of mouse periodontal ligament cells

The periodontal ligament (PDL) is a critical component in maintaining tooth stability. It is composed of cells and an extracellular matrix (ECM), each with unique roles in tissue function and homeostasis. Secreted protein acidic and rich in cysteine (SPARC), a calcium-binding matricellular glycoprotein, plays a crucial role in regulating ECM assembly and turnover, alongside facilitating cellular-ECM interactions. In the present study, mass spectrometry-based proteomics was used to assess the impacts of Sparc-knockout (KO) on PDL-derived cells. Results demonstrated that Sparc-KO significantly reduces ECM production and alters its composition with increased levels of type I collagen. Despite this increase in Sparc-KO, type I collagen was not likely to be effectively integrated into the fibrils due to collagen cross-linking impairment. Furthermore, the pathway and process enrichment analyses suggested that SPARC plays a protective role against ECM degradation by antagonistically interacting with cell-surface collagen receptors. These findings provide detailed insights into the multifaceted role of SPARC in ECM organization, including its impact on ECM production, collagen regulation, and interactions with various cellular compartments. A better understanding of these complex mechanisms is crucial for comprehending the causes of periodontal disease and tissue regeneration, where precise control of ECM organization is necessary.

Multiomics analysis of cultured mouse periodontal ligament cell-derived extracellular matrix

A comprehensive understanding of the extracellular matrix (ECM) is essential for developing biomimetic ECM scaffolds for tissue regeneration. As the periodontal ligament cell (PDLC)-derived ECM has shown potential for periodontal tissue regeneration, it is vital to gain a deeper understanding of its comprehensive profile. Although the PDLC-derived ECM exhibits extracellular environment similar to that of periodontal ligament (PDL) tissue, details of its molecular composition are lacking. Thus, using a multiomics approach, we systematically analyzed cultured mouse PDLC-derived ECM and compared it to mouse PDL tissue as a reference. Proteomic analysis revealed that, compared to PDL tissue, the cultured PDLC-derived ECM had a lower proportion of fibrillar collagens with increased levels of glycoprotein, corresponding to an immature ECM status. The gene expression signature was maintained in cultured PDLCs and was similar to that in cells from PDL tissues, with additional characteristics representative of naturally occurring progenitor cells. A combination of proteomic and transcriptomic analyses revealed that the cultured mouse PDLC-derived ECM has multiple advantages in tissue regeneration, providing an extracellular environment that closely mimics the environment in the native PDL tissue. These findings provide valuable insights for understanding PDLC-derived ECM and should contribute to the development of biomimetic ECM scaffolds for reliable periodontal tissue regeneration.

3D Electrospun Polycaprolactone Scaffolds to Assess Human Periodontal Ligament Cells Mechanobiological Behaviour

While periodontal ligament cells are sensitive to their 3D biomechanical environment, only a few 3D in vitro models have been used to investigate the periodontal cells mechanobiological behavior. The objective of the current study was to assess the capability of a 3D fibrous scaffold to transmit a mechanical loading to the periodontal ligament cells. Three-dimensional fibrous polycaprolactone (PCL) scaffolds were synthetized through electrospinning. Scaffolds seeded with human periodontal cells (10³ mL^-1) were subjected to static (n = 9) or to a sinusoidal axial compressive loading in an in-house bioreactor (n = 9). At the end of the culture, the dynamic loading seemed to have an influence on the cells' morphology, with a lower number of visible cells on the scaffolds surface and a lower expression of actin filament. Furthermore, the dynamic loading presented a tendency to decrease the Alkaline Phosphatase activity and the production of Interleukin-6 while these two biomolecular markers were increased after 21 days of static culture. Together, these results showed that load transmission is occurring in the 3D electrospun PCL fibrous scaffolds, suggesting that it can be used to better understand the periodontal ligament cells mechanobiology. The current study shows a relevant way to investigate periodontal mechanobiology using 3D fibrous scaffolds.

Biomechanical Modelling for Tooth Survival Studies: Mechanical Properties, Loads and Boundary Conditions-A Narrative Review

Recent biomechanical studies have focused on studying the response of teeth before and after different treatments under functional and parafunctional loads. These studies often involve experimental and/or finite element analysis (FEA). Current loading and boundary conditions may not entirely represent the real condition of the tooth in clinical situations. The importance of homogenizing both sample characterization and boundary conditions definition for future dental biomechanical studies is highlighted. The mechanical properties of dental structural tissues are presented, along with the effect of functional and parafunctional loads and other environmental and biological parameters that may influence tooth survival. A range of values for Young's modulus, Poisson ratio, compressive strength, threshold stress intensity factor and fracture toughness are provided for enamel and dentin; as well as Young's modulus and Poisson ratio for the PDL, trabecular and cortical bone. Angles, loading magnitude and frequency are provided for functional and parafunctional loads. The environmental and physiological conditions (age, gender, tooth, humidity, etc.), that may influence tooth survival are also discussed. Oversimplifications of biomechanical models could end up in results that divert from the natural behavior of teeth. Experimental validation models with close-to-reality boundary conditions should be developed to compare the validity of simplified models.

Anchorage effects of ligation and direct occlusion in orthodontics: A finite element analysis

BACKGROUND AND OBJECTIVE

During orthodontic treatment, the figure-of-eight ligature and the physiological occlusion play an important role in providing anchorage effects. However, their effects on reaction forces of tooth and stress state in periodontal ligament (PDL) have not been quantitatively evaluated yet. In this study, we presented a finite element analysis process for simulating posterior molar ligature and direct occlusion during orthodontics in order to quantitatively assess their anchorage effects.

METHODS

A high precision 3D biomechanical model containing upper and lower teeth, PDL, brackets and archwire was generated from the images of computed tomographic scan and sophisticated modelling procedures. The orthodontic treatment of closing the extraction gap was simulated via the finite element method to evaluate the biomechanical response of the molars under the conditions with or without ligation. The simulations were divided into experimental and control groups. In the experimental group, orthodontic force of 1 N was first applied, then direct occlusal forces of 3 and 10 N were applied on each opposite tooth. While in the control group, occlusal forces were applied without orthodontic treatment. The tooth displacement, the stress state in the PDL and the directions of the resultant forces on each tooth were evaluated.

RESULTS

In the case of molars ligated, the maximum hydrostatic stress in the molars' PDL decreases by 60%. When an initial tooth displacement of several microns occurs in response to an orthodontic force, the direction of the occlusal force changes simultaneously. Even a moderate occlusal force (3 N per tooth) can almost completely offset the mesial forces on the maxillary teeth, thus to provide effective anchorage effect for the orthodontics.

CONCLUSIONS

The proposed method is effective for simulating ligation and direct occlusion. Figure-of-eight ligature can effectively disperse orthodontic forces on the posterior teeth, while a good original occlusal relationship provides considerable anchorage effects in orthodontics.

Open-Full-Jaw: An open-access dataset and pipeline for finite element models of human jaw

BACKGROUND

State-of-the-art finite element studies on human jaws are mostly limited to the geometry of a single patient. In general, developing accurate patient-specific computational models of the human jaw acquired from cone-beam computed tomography (CBCT) scans is labor-intensive and non-trivial, which involves time-consuming human-in-the-loop procedures, such as segmentation, geometry reconstruction, and re-meshing tasks. Therefore, with the current practice, researchers need to spend considerable time and effort to produce finite element models (FEMs) to get to the point where they can use the models to answer clinically-interesting questions. Besides, any manual task involved in the process makes it difficult for the researchers to reproduce identical models generated in the literature. Hence, a quantitative comparison is not attainable due to the lack of surface/volumetric meshes and FEMs.

METHODS

We share an open-access repository composed of 17 patient-specific computational models of human jaws and the utilized pipeline for generating them for reproducibility of our work. The used pipeline minimizes the required time for processing and any potential biases in the model generation process caused by human intervention. It gets the segmented geometries with irregular and dense surface meshes and provides reduced, adaptive, watertight, and conformal surface/volumetric meshes, which can directly be used in finite element (FE) analysis.

RESULTS

We have quantified the variability of our 17 models and assessed the accuracy of the developed models from three different aspects; (1) the maximum deviations from the input meshes using the Hausdorff distance as an error measurement, (2) the quality of the developed volumetric meshes, and (3) the stability of the FE models under two different scenarios of tipping and biting.

CONCLUSIONS

The obtained results indicate that the developed computational models are precise, and they consist of quality meshes suitable for various FE scenarios. We believe the provided dataset of models including a high geometrical variation obtained from 17 different models will pave the way for population studies focusing on the biomechanical behavior of human jaws.

Biomechanical evaluation of the unilateral crossbite on the asymmetrical development of the craniofacial complex. A mechano-morphological approach

BACKGROUND AND OBJECTIVE

The occlusion effect on the craniofacial development is a controversial topic that has attracted the interest of many researchers but that remains unclear, mainly due to the difficulties on measure its mechanical response experimentally. This mechano-morphological relationship of the craniofacial growth is often explained by the periosteal and capsular matrices of the functional matrix hypothesis (FMH); however, its outcomes have not been analytically demonstrated yet. This computational study aims, therefore, to analytically demonstrate the mechano-morphological relationship in the craniofacial development of children with unilateral crossbite (UXB) using the finite element (FE) method.

METHODS

The craniofacial complex asymmetry of ten children, five of whom exhibit UXB, was 3D-analysed and compared with the biomechanical response computed from a FE analysis of each patient's occlusion. Due to the complexity of the geometry and the multitude of contacts involved, the inherent limitations of the model were evaluated by comparing computed occlusal patterns with those recorded by an occlusal analysis on 3D printed copies.

RESULTS

Comparison's outcomes proved the reliability of our models with just a deviation error below 6% between both approaches. Out of validation process, computational results showed that the significant elongation of mandibular branch in the contralateral side could be related to the mandibular shift and increase of thickness on the crossed side, and particularly of the posterior region. These morphological changes could be associated with periodontal overpressure (>4.7 kPa) and mandibular over deformation (0.002 ε) in that side, in agreement with the periosteal matrix's principles. Furthermore, the maxilla's transversal narrowing and the elevation of the maxillary and zygomatic regions on the crossed side were statistically demonstrated and seem to be related with their respective micro displacements at occlusion, as accounted by their specific capsule matrices. Our results were consistent with those reported clinically and demonstrated analytically the mechano-morphological relationship of children's craniofacial development based on the FMH's functional matrices.

CONCLUSIONS

This study is a first step in the understanding of the occlusion's effect on the craniofacial development by computational methods. Our approach could help future engineers, researchers and clinicians to understand better the aetiology of some dental malocclusions and functional disorders improve the diagnosis or even predict the craniofacial development.

Modelling and evaluating periodontal ligament mechanical behaviour and properties: A scoping review of current approaches and limitations

This scoping review is intended to synthesize the techniques proposed to model the tooth-periodontal ligament-bone complex (TPBC), while also evaluating the suggested periodontal ligament (PDL) material properties. It is concentrated on the recent advancements on the PDL and TPBC models, while identifying the advantages and limitations of the proposed approaches. Systematic searches were conducted up to December 2020 for articles that proposed PDL models to assess orthodontic tooth movement in Compendex, Web of Science, EMBASE, MEDLINE, PubMed, ScienceDirect, Google Scholar and Scopus databases. Although there have been many studies focused on the evaluation of PDL material properties through numerous modelling approaches, only a handful of approaches have been identified to investigate the interface properties of the PDL as a complete dynamical system (TPBC models). Past reviews on the analytical and experimental determination of the PDL properties already show a concerning range in reported output values-some nearly six orders of magnitude in difference-that strongly suggested the need for further investigation. Surprisingly, it has not yet been possible to determine a narrower range of values for the PDL material properties. Moreover, very few scientific approaches address the TPBC as an integrated complex system model. In consequence, current methods for capturing the PDL material behaviour in a clinical setting are limited and inconclusive. This synthesis encourages more systematic, pragmatic and phenomenological research in this area.

Tissue Engineering for Periodontal Ligament Regeneration: Biomechanical Specifications

The periodontal biomechanical environment is very difficult to investigate. By the complex geometry and composition of the periodontal ligament (PDL), its mechanical behavior is very dependent on the type of loading (compressive versus tensile loading; static versus cyclic loading; uniaxial versus multiaxial) and the location around the root (cervical, middle, or apical). These different aspects of the PDL make it difficult to develop a functional biomaterial to treat periodontal attachment due to periodontal diseases. This review aims to describe the structural and biomechanical properties of the PDL. Particular importance is placed in the close interrelationship that exists between structure and biomechanics: the PDL structural organization is specific to its biomechanical environment, and its biomechanical properties are specific to its structural arrangement. This balance between structure and biomechanics can be explained by a mechanosensitive periodontal cellular activity. These specifications have to be considered in the further tissue engineering strategies for the development of an efficient biomaterial for periodontal tissues regeneration.

Comparative Analysis of Stress in the Periodontal Ligament and Center of Rotation in the Tooth after Orthodontic Treatment Depending on Clear Aligner Thickness-Finite Element Analysis Study

Lately, in orthodontic treatments, the use of transparent aligners for the correction of malocclusions has become prominent owing to their intrinsic advantages such as esthetics, comfort, and minimal maintenance. Attempts at improving upon this technology by varying various parameters to investigate the effects on treatments have been carried out by several researchers. Here, we aimed to investigate the biomechanical and clinical effects of aligner thickness on stress distributions in the periodontal ligament and changes in the tooth's center of rotation. Dental finite element models comprising the cortical and cancellous bones, gingiva, teeth, and nonlinear viscoelastic periodontal ligaments were constructed, validated, and used together with aligner finite element models of different aligner thicknesses to achieve the goal of this study. The finite element analyses were conducted to simulate the actual orthodontic aligner treatment process for the correction of malocclusions by generating pre-stresses in the aligner and allowing the aligner stresses to relax to induce tooth movement. The results of the analyses showed that orthodontic treatment in lingual inclination and axial rotation with a 0.75 mm-thick aligner resulted in 6% and 0.03% higher principal stresses in the periodontal ligament than the same treatment using a 0.05 mm-thick aligner, respectively. Again, for both aligner thicknesses, the tooth's center of rotation moved lingually and towards the root direction in lingual inclination, and diagonally from the long axis of the tooth in axial rotation. Taken together, orthodontic treatment for simple malocclusions using transparent aligners of different thicknesses will produce a similar effect on the principal stresses in the periodontal ligament and similar changes in the tooth's center of rotation, as well as sufficient tooth movement. These findings provide orthodontists and researchers clinical and biomechanical evidence about the effect of transparent aligner thickness selection and its effect on orthodontic treatment.

Finite element analysis of torque induced orthodontic bracket slot deformation in various bracket-archwire contact assembly

BACKGROUND AND OBJECTIVES

Orthodontic fixed appliance therapy involves alignment of teeth through the bracket and archwires. The archwire twist (torque) imparts significant forces inside the bracket slot in refining the teeth position at the end of treatment. The objective of this in- silico study was to evaluate the torque induced bracket slot deformation in the commonly used 0.018 inch (") and 0.022" conventional Stainless Steel (SS) brackets with clinically relevant archwires during various angles of twist.

METHODS

SS maxillary central incisor brackets of 0.018" width × 0.022" depth (0.457 mm × 0.558 mm) and 0.022" width × 0.028" depth (0.558 mm × 0.711 mm) were used. The SS archwires of 0.016" width × 0.022" depth (0.406 mm × 0.558 mm), 0.017" width × 0.025" depth (0.431 mm × 0.635 mm), 0.019" width × 0.025" depth (0.482 mm × 0.635 mm) and 0.021" width × 0.025" depth (0.533 mm × 0.635 mm) were engaged in the respective bracket slots. The assembled bracket-archwire Finite Element (FE) models were constructed. The archwire torque, the top, middle and bottom slot deformations (TSD, MSD, BSD) were obtained for the bracket-archwire combinations for various angles of archwire twist using FE Analysis (FEA).

RESULTS

The torque, TSD, MSD and BSD for 30^o twist of 0.016" × 0.022" archwire in 0.018" slot were 28.13 Nmm, 35.71 µm, 21.51 µm and 15.67 µm respectively, and for 0.017" × 0.025" archwire were 50.18 Nmm, 54.52 µm, 32.47 µm and 19.11 µm respectively. Similarly for 0.019" × 0.025" archwire in 0.022" slot and 0.021" × 0.025" archwire in 0.022" slot they were 38.82 Nmm, 50.78 µm, 31.47 µm and 16.82 µm, and 60.22 Nmm, 65.22 µm, 36.44 µm and 22.68 µm respectively.

CONCLUSIONS

The slot deformation was present in both 0.018" and 0.022" brackets which increased as the angle of twist increased. The TSD were higher than the MSD and BSD in all the bracket-archwire combinations. We conclude that there is only elastic deformation of bracket slots upto 30^o angle of twist and clinicians could maintain within this torque limits to avoid plastic deformation leading to improper teeth position.

Parameter identification for the simulation of the periodontal ligament during the initial phase of orthodontic tooth movement

The paper is concerned with simulation of the periodontal ligament response to force in the initial phase of orthodontic tooth movement. This is based on two previous investigations, a in vitro experiment with specimens of porcine mandibular premolars and a in vivo experiment on human upper first incisors. For the curve fit of the in vitro experiment a model function, assuming viscoelasticity, was introduced. The viscoelastic model function was augmented by a ramp rise time term, to account for observed dependence of the response on actuator velocity, and a previous load history term, to account for the effect of the previous tests on the current test. The correlation coefficient of a curve fit for all tests grouped together was R2=0.98. Next, a curve fit of the in vivo experiment was done. Good correlation was found for a simplified model function, without viscoelastic term (R2=0.96). For both tests, in vitro and in vivo, the ramp rise time term improved correlation. A finite element model of the specimen of the in vitro experiment was created. For the PDL a hyperelastic constitutive model for compressible material was used and model parameters were identified. The present work indicates that the macroscopic response of the periodontal ligament to an external load can be simulated with a poro-visco-hyperelastic model. The simulation showed that poroelastic behaviour will gradually cease when viscoelastic relaxation progresses. This followed also from dimensionless analysis. As a consequence, for slow loading, or if initial response to fast loading is not of interest, a visco-hyperelastic model may suffice. To identify parameters of the finite element model several optimisation problems were solved. A model function, which can be regarded as a reduced order model, allowed a full factorial experiment (analysis) at low cost, to identify initial parameters. The thus found parameters were further refined with an optimum interpolation meta-model. That is, for limited number of parameter combinations the response was simulated with the finite element model and a refined parameter study was conducted by means of optimal interpolation. The thus found optimal parameters were verified by simulation with the finite element model. Optimal interpolation is computationally cheap, which allowed full factorial experiments at low cost.

Finite element analysis and experimental evaluation on stress distribution and sensitivity of dental implants to assess optimum length and thread pitch

BACKGROUND AND OBJECTIVE

The dental implant is one of the long term proper remedies to recover a missed tooth as a different prosthetic rehabilitation way. The finite element (FE) method and photoelasticity test are employed to achieve stress distribution and sensitivity in dental implants in order to obtain optimum length and thread pitch.

METHODS

The finite element method and experimental test are developed to evaluate stress distribution and sensitivity around dental implants. Three dimensional FE models of implant-abutment, cortical bone and cancellous bone are created by considering a variation of 0.6 to -1 mm on threads pitch while the implant lengths range from 8.5 mm to 13 mm. Then, axial and oblique forces are applied to the models to obtain the resultant stress contours.

RESULTS

The results indicate that the resultant von Mises stresses in the implant-abutment, cortical bones, and cancellous bones are different. The optimized setting for length and pitch is suggested according to maximum von Mises stress and sensitivity analysis.

CONCLUSIONS

It is concluded that the present FE model accurately predicts stress distribution pattern in dental implants. The results indicate that sensitivity of length play a more significant role in comparison with thread pitch. The accuracy of FEM results in comparison with those of the photoelasticity test recommends applying computation methods in medical practice as great potential in terms of future studies.

Computational characterization of the porous-fibrous behavior of the soft tissues in the temporomandibular joint

The prevalence and severity of temporomandibular joint (TMJ) disorders have led to growing research interest in the development of new biomaterials and medical devices for TMJ implant designs. In computational designs, however, the time and stretch direction dependences of the TMJ soft tissues behavior are not considered and they are frequently based on measurements taken from non‐human species or from joints that differ markedly from the human TMJ. The aim of this study was to accurately characterize the porous‐fibrous properties of the TMJ soft tissues by simulating previously published experimental tests, to assist professionals in the design of new TMJ implants. To that end, material parameters were determined assuming a uniform fiber orientation throughout the entire sample. This assumption was then tested by comparing these results with those of considering multiple regions and distinct fiber orientations in each sample. Our findings validated the use of a transversely isotropic hyperelastic material model to characterize the direction dependent behavior of TMJ soft tissues and its combination with porous hyperfoam material models to mimic the compressive response of the TMJ disc. In conclusion, constitutive model proposed accurately reproduce the mechanical response of the TMJ soft tissues at different strain rates and stretch directions.