The effects of the periodontal ligament on mandibular stiffness: a study combining finite element analysis and geometric morphometrics

It is generally accepted that the periodontal ligament (PDL) plays a crucial role in transferring occlusal forces from the teeth to the alveolar bone. Studies using finite element analysis (FEA) have helped to better understand this role and show that the stresses and strains in the alveolar bone are influenced by whether and how PDL is included in FE models. However, when the overall distribution of stresses and strains in crania and mandibles are of interest, PDL is often not included in FE models, although little is known about how this affects the results. Here we study the effect of representing PDL as a layer of solid material with isotropic homogeneous properties in an FE model of a human mandible using a novel application of geometric morphometrics. The results show that the modelling of the PDL affects the deformation and thus strain magnitudes not only of the alveolar bone around the biting tooth, but that the whole mandible deforms differently under load. As a result, the strain in the mandibular corpus is significantly increased when PDL is included, while the strain in the bone beneath the biting tooth is reduced. These results indicate the importance of the PDL in FE studies. Thus we recommend that the PDL should be included in FE models of the masticatory apparatus, with tests to assess the sensitivity of the results to changes in the Young's modulus of the PDL material.     

Numerical study of tension and strain distribution around rat molars]

A knowledge of the mechanical processes triggered in the bone and periodontal ligament (PDL) by orthodontic forces applied to a tooth is of decisive importance for an understanding of the subsequent remodelling around the tooth. To investigate these mechanical relationships, three-dimensional finite element (FE) models of the first lower molar in the rat were established. On the basis of digitized serial histological sections, these FE models were generated semi-automatically. Using various simplified geometrical variations, an appropriate FE model for the analysis of the stress and strain distributions was established. The numerical analyses were carried out under a mesially directed force of 0.1 N. Stress distributions in the bone and PDL showed a similar pattern, while strains in the bone were lower than in the PDL by a factor of 10-5. The data confirm the assumption that strain patterns in the PDL may be the key stimulus of bone remodelling.     

Effect of variable periodontal ligament thickness and its non-linear material properties on the location of a tooth’s centre of resistance

In orthodontics, the 3D translational and rotational movement of a tooth is determined by the force–moment system applied and the location of the tooth’s centre of resistance (CR). Because of the practical constraints of in-vivo experiments, the finite element (FE) method is commonly used to determine the CR. The objective of this study was to investigate the geometric model details required for accurate CR determination, and the effect of material non-linearity of the periodontal ligament (PDL). A FE model of a human lower canine derived from a high-resolution µCT scan (voxel size: 50 µm) was investigated by applying four different modelling approaches to the PDL. These comprised linear and non-linear material models, each with uniform and realistic PDL thickness. The CR locations determined for the four model configurations were in the range 37.2–45.3% (alveolar margin: 0%; root apex: 100%). We observed that a non-linear material model introduces load-dependent results that are dominated by the PDL regions under tension. Load variation within the range used in clinical orthodontic practice resulted in CR variations below 0.3%. Furthermore, the individualized realistic PDL geometry shifted the CR towards the alveolar margin by 2.3% and 2.8% on average for the linear and non-linear material models, respectively. We concluded that for conventional clinical therapy and the generation of representative reference data, the least sophisticated modelling approach with linear material behaviour and uniform PDL thickness appears sufficiently accurate. Research applications that require more precise treatment monitoring and planning may, however, benefit from the more accurate results obtained from the non-linear constitutive law and individualized realistic PDL geometry.     

The Biomechanical Function of Periodontal Ligament Fibres in Orthodontic Tooth Movement

Orthodontic tooth movement occurs as a result of resorption and formation of the alveolar bone due to an applied load, but the stimulus responsible for triggering orthodontic tooth movement remains the subject of debate. It has been suggested that the periodontal ligament (PDL) plays a key role. However, the mechanical function of the PDL in orthodontic tooth movement is not well understood as most mechanical models of the PDL to date have ignored the fibrous structure of the PDL. In this study we use finite element (FE) analysis to investigate the strains in the alveolar bone due to occlusal and orthodontic loads when PDL is modelled as a fibrous structure as compared to modelling PDL as a layer of solid material. The results show that the tension-only nature of the fibres essentially suspends the tooth in the tooth socket and their inclusion in FE models makes a significant difference to both the magnitude and distribution of strains produced in the surrounding bone. The results indicate that the PDL fibres have a very important role in load transfer between the teeth and alveolar bone and should be considered in FE studies investigating the biomechanics of orthodontic tooth movement.     

A visco-elastic model for the prediction of orthodontic tooth movement

This study presents a biomechanical model of orthodontic tooth movement. Although such models have already been presented in the literature, most of them incorporate computationally expensive finite elements (FE) methods to determine the strain distribution in the periodontal ligament (PDL). In contrast, the biomechanical model presented in this work avoids the use of FE methods. The elastic deformation of the PDL is modelled using an analytical approach, which does not require setting up a 3D model of the tooth. The duration of the lag phase is estimated using the calculated hydrostatic stresses, and bone remodelling is predicted by modelling the alveolar bone as a viscous material. To evaluate the model, some typically used motion patterns were simulated and a sensitivity analysis was carried out on the parameters. Results show that despite some shortcomings, the model is able to describe commonly used motion patterns in orthodontic tooth movement, in both single- and multi-rooted teeth.     

Hydro-mechanical coupling in the periodontal ligament: A porohyperelastic finite element model

Harmonic tension–compression tests at 0.1, 0.5 and 1 Hz on hydrated bovine periodontal ligament (PDL) were numerically simulated. The process was modeled by finite elements (FE) within the framework of poromechanics, with the objective of isolating the contributions of the solid- and fluid phases. The solid matrix was modeled as a porous hyperelastic material (hyperfoam) through which the incompressible fluid filling the pores flowed in accordance with the Darcy’s law. The hydro-mechanical coupling between the porous solid matrix and the fluid phase circulating through it provided an apparent time-dependent response to the PDL, whose rate of deformation depended on the permeability of the porous solid with respect to the interstitial fluid. Since the PDL was subjected to significant deformations, finite strains were taken into account and an exponential dependence of PDL permeability on void ratio – and therefore on the deformation state – was assumed. PDL constitutive parameters were identified by fitting the simulated response to the experimental data for the tests at 1 Hz. The values thus obtained were then used to simulate the tests at 0.1 and 0.5 Hz. The results of the present simulation demonstrate that a porohyperelastic model with variable permeability is able to describe the two main aspects of the PDL’s response: (1) the dependency on strain-rate—the saturated material can develop volumetric strains by only exchanging fluid and (2) the asymmetry between tension and compression, which is due to the effect of both the permeability and the elastic properties on deformation.     

Semi-automatic generation of finite element meshes fo dental preparations]

The mechanical properties and elastic behaviour of periodontal tissue are a decisive factor in understanding initial tooth mobility and bone remodelling processes in orthodontics. An experimental set-up was designed to precisely determine a tooth's elastic response to different loading conditions. Segments of pig's maxilla bearing separated molars were used, and their mechanical response to loading was recorded. Subsequently, finite element analysis (FEA) was performed on the basis of the experimental data. The combination of experimental and numerical methods was used to determine the material properties of the periodontal ligament (PDL). The geometries of the preparations were reconstructed and FE meshes generated semi-automatically with the aid of the special computer program, CAGOG (Computer Aided Generator for Orthodontic Geometries) to optimally match the experimental geometry. Nonlinear material parameters were determined for the PDL and verified by comparing experimental and numerical results obtained in other specimens with an error of about 10%. This good correlation indicates that the selected method of mesh generation is appropriate for creating realistic FE models that can be compared with experimental results.     

An investigation into the use and limitations of different spatial integration schemes and finite element software in head impact analyses

To understand the mechanopathogenesis of brain lesions, finite element (FE) head models are used. There is a broad range of material properties, contact interfaces and integration schemes used for the different parts in current FE head models. The effect of material behaviour and contact definitions on a head impact analysis is reported in the literature, whereas the effect of FE integration schemes is a rather unexplored domain. This paper starts with the development of a simplified head model to which small adaptations are made in the integration scheme to obtain multiple analyses that are compared using an accident reconstruction. The performed study highlighted potential hazards of different integration schemes and the significant effect they have on the simulated mechanical responses of the head. Based on a comparison between FE softwares using an impact test and patch test, it was seen that also the software could have an effect on the FE analysis results.     

Finite element simulation of the behavior of the periodontal ligament: A validated nonlinear contact model

Due to its significance in tooth movement, the stress/deformation field of periodontium and the alveolar bone remodeling process, periodontal ligament (PDL) cannot be excluded from the studies investigating dental biomechanics regarding its excessive deformability. Therefore, many analytical and numerical researches are carried out to simulate its response and to create a constitutive model via experiments intending to discover the material properties of PDL. The aim of this study is to formulate a user specified contact model that can be used in conjunction with finite element (FE) software and reflects PDL’s influence on neighboring structures based on the currently available information, without requiring an actual volumetric finite element mesh of ligament. The results show good agreement with available experimental tooth mobility data. Smooth stress fields are obtained on the tooth root and alveolar bone, which is a significant aspect in bone-remodeling studies. The advantage of simulating PDL as a contact model at the interface of tooth root and the alveolar process instead of a solid-meshed FE model with poor geometric morphology and/or very dense mesh is expected to save pre/post-processing workforce, to increase the accuracy and to contribute to the smoothness of interface stress distributions.     

A novel approach for extracting viscoelastic parameters of living cells through combination of inverse finite element simulation and Atomic Force Microscopy

Dynamic mechanical behaviour of living cells has been described by viscoelasticity. However, quantitation of the viscoelastic parameters for living cells is far from sophisticated. In this paper, combining inverse finite element (FE) simulation with Atomic Force Microscope characterization, we attempt to develop a new method to evaluate and acquire trustworthy viscoelastic index of living cells. First, influence of the experiment parameters on stress relaxation process is assessed using FE simulation. As suggested by the simulations, cell height has negligible impact on shape of the force–time curve, i.e. the characteristic relaxation time; and the effect originates from substrate can be totally eliminated when stiff substrate (Young’s modulus larger than 3 GPa) is used. Then, so as to develop an effective optimization strategy for the inverse FE simulation, the parameters sensitivity evaluation is performed for Young’s modulus, Poisson’s ratio, and characteristic relaxation time. With the experiment data obtained through typical stress relaxation measurement, viscoelastic parameters are extracted through the inverse FE simulation by comparing the simulation results and experimental measurements. Finally, reliability of the acquired mechanical parameters is verified with different load experiments performed on the same cell.     

Comparison of micro-level and continuum-level voxel models of the proximal femur

Continuum-level finite element (FE) models became standard computational tools for the evaluation of bone mechanical behavior from in vivo computed tomography scans. Such scans do not account for the anisotropy of the bone. Instead, local mechanical properties in the continuum-level FE models are assumed isotropic and are derived from bone density, using statistical relationships. Micro-FE models, on the other hand, incorporate the anisotropic structure in detail. This study aimed to quantify the effects of assumed isotropy, by comparing continuum-level voxel models of a healthy and a severely osteoporotic proximal femur with recently analyzed micro-FE models of the same bones. The micro-model element size was coarsened to generate continuum FE models with two different element sizes (0.64 and 3.04 mm) and two different density–modulus relationships found in the literature for wet and ash density. All FE models were subjected to the same boundary conditions that simulated a fall to the side, and the stress and strain distributions, model stiffness and yield load were compared. The results indicated that the stress and strain distributions could be reproduced well with the continuum models. The smallest differences between the continuum-level model and micro-level model predictions of the stiffness and yield load were obtained with the coarsest element size. Better results were obtained for both continuum-element sizes when isotropic moduli were based on ash density rather than wet density.     

A subject-specific pelvic bone model and its application to cemented acetabular replacements

A subject-specific three-dimensional finite element (FE) pelvic bone model has been developed and applied to the study of bone–cement interfacial response in cemented acetabular replacements. The pelvic bone model was developed from CT scan images of a cadaveric pelvis and validated against the experiment data obtained from the same specimen at a simulated single-legged stance. The model was then implanted with a cemented acetabular cup at selected positions to simulate some typical implant conditions due to the misplacement of the cup as well as a standard cup condition. For comparison purposes, a simplified FE model with homogeneous trabecular bone material properties was also generated and similar implant conditions were examined.The results from the homogeneous model are found to underestimate significantly both the peak von Mises stress and the area of the highly stressed region in the cement near the bone–cement interface, compared with those from the subject-specific model. Non-uniform cement thickness and non-standard cup orientation seem to elevate the highly stressed region as well as the peak stress near the bone–cement interface.     

Mechanical behavior of an individual adherent MLO-Y4 osteocyte under shear flow

Mechanical properties of a single cell and its mechanical response under stimulation play an important role in regulating interactions between cell and extracellular matrix and affecting mechanotransduction. Osteocytes exhibit solid-like viscoelastic behavior in response to the interstitial fluid shear resulting from tissue matrix deformation. This study intends to quantitatively describe the mechanical behavior of osteocytes combining in vitro experiment and fluid–structure interaction (FSI) finite element (FE) model. The cell is configured in the FSI FE model using the observed data from quasi-3D images. Instead of simply assigning the cellular viscoelastic parameters by statistical data, the mechanical parameters are determined by an iterative algorithm comparing the experimental and the computational results from the FE model. The viscoelastic parameters of osteocytes are obtained as: the equilibrium elasticity modulus \(k_{1}=0.15\pm 0.038\,\hbox {kPa}\), instantaneous elasticity modulus \((k_{1}+k_{2})=0.77\pm 0.23\,\hbox {kPa}\), viscosity coefficient \(\eta =1.38\pm 0.33\,\hbox {kPa}\,\hbox {s}\). A novel index to quantify the cell adhesion is also put forward. In addition, an interesting competition phenomenon is revealed on the cell surface concerning stress and strain, i.e., the place with high stress has low strain and that with low stress has high strain. The proposed method provides a novel technique to study the mechanical behavior of individual adherent cell in vitro. It is believed that this quantitative technique not only determines cell mechanical behavior but also helps elucidate the mechanism of mechanotransduction in various types of cells.     

Finite element modelling approaches for well-ordered porous metallic materials for orthopaedic applications: cost effectiveness and geometrical considerations

The mechanical properties of well-ordered porous materials are related to their geometrical parameters at the mesoscale. Finite element (FE) analysis is a powerful tool to design well-ordered porous materials by analysing the mechanical behaviour. However, FE models are often computationally expensive. This article aims to develop a cost-effective FE model to simulate well-ordered porous metallic materials for orthopaedic applications. Solid and beam FE modelling approaches are compared, using finite size and infinite media models considering cubic unit cell geometry. The model is then applied to compare two unit cell geometries: cubic and diamond. Models having finite size provide similar results than the infinite media model approach for large sample sizes. In addition, these finite size models also capture the influence of the boundary conditions on the mechanical response for small sample sizes. The beam FE modelling approach showed little computational cost and similar results to the solid FE modelling approach. Diamond unit cell geometry appeared to be more suitable for orthopaedic applications than the cubic unit cell geometry.     

Finite element analysis of acetabular fractures—development and validation with a synthetic pelvis

Acetabular fracture presents a challenging situation to trauma surgeons today due to its complexity. Finite element (FE) models can be of great help as they can improve the surgical planning and post surgery patient management for those with acetabular fractures. We have developed a non-linear finite element model of the pelvis and validated its fracture prediction capability with synthetic polyurethane pelves. A mechanical experiment was performed with the synthetic bones and fracture loads and patterns were observed for two different loading cases. Fracture loads predicted by our FE model were within one standard deviation of the experimental fracture loads for both loading cases. The incipient fracture pattern predicted by the model also resembled the actual pattern from the experiment. Although it is not a complete validation with human cadaver bones, the good agreement between model predictions and experimental results indicate the validity of our approach in using non-linear FE formulation along with contact conditions in predicting bone fractures.     

Digital image correlation and finite element modelling as a method to determine mechanical properties of human soft tissue in vivo

The mechanical properties of human soft tissue are crucial for impact biomechanics, rehabilitation engineering, and surgical simulation. Validation of these constitutive models using human data remains challenging and often requires the use of non-invasive imaging and inverse finite element (FE) analysis. Post-processing data from imaging methods such as tagged magnetic resonance imaging (MRI) can be challenging. Digital image correlation (DIC), however, is a relatively straightforward imaging method. DIC has been used in the past to study the planar and superficial properties of soft tissue and excised soft tissue layers. However, DIC has not been used to non-invasive study of the bulk properties of human soft tissue in vivo. Thus, the goal of this study was to assess the use of DIC in combination with FE modelling to determine the bulk material properties of human soft tissue. Indentation experiments were performed on a silicone gel soft tissue phantom. A two camera DIC setup was then used to record the 3D surface deformation. The experiment was then simulated using a FE model. The gel was modelled as Neo-Hookean hyperelastic, and the material parameters were determined by minimising the error between the experimental and FE data. The iterative FE analysis determined material parameters (μ=1.80 kPa, K=2999 kPa) that were in close agreement with parameters derived independently from regression to uniaxial compression tests (μ=1.71 kPa, K=2857 kPa). Furthermore the FE model was capable of reproducing the experimental indentor force as well as the surface deformation found (R²=0.81). It was therefore concluded that a two camera DIC configuration combined with FE modelling can be used to determine the bulk mechanical properties of materials that can be represented using hyperelastic Neo-Hookean constitutive laws.     

Development of a Finite Element Model of Decompressive Craniectomy

Decompressive craniectomy (DC), an operation whereby part of the skull is removed, is used in the management of patients with brain swelling. While the aim of DC is to reduce intracranial pressure, there is the risk that brain deformation and mechanical strain associated with the operation could damage the brain tissue. The nature and extent of the resulting strain regime is poorly understood at present. Finite element (FE) models of DC can provide insight into this applied strain and hence assist in deciding on the best surgical procedures. However there is uncertainty about how well these models match experimental data, which are difficult to obtain clinically. Hence there is a need to validate any modelling approach outside the clinical setting. This paper develops an axisymmetric FE model of an idealised DC to assess the key features of such an FE model which are needed for an accurate simulation of DC. The FE models are compared with an experimental model using gelatin hydrogel, which has similar poro-viscoelastic material property characteristics to brain tissue. Strain on a central plane of the FE model and the front face of the experimental model, deformation and load relaxation curves are compared between experiment and FE. Results show good agreement between the FE and experimental models, providing confidence in applying the proposed FE modelling approach to DC. Such a model should use material properties appropriate for brain tissue and include a more realistic whole head geometry.     

The biomechanics of human femurs in axial and torsional loading: comparison of finite element analysis, human cadaveric femurs, and synthetic femurs

To assess the performance of femoral orthopedic implants, they are often attached to cadaveric femurs, and biomechanical testing is performed. To identify areas of high stress, stress shielding, and to facilitate implant redesign, these tests are often accompanied by finite element (FE) models of the bone/implant system. However, cadaveric bone suffers from wide specimen to specimen variability both in terms of bone geometry and mechanical properties, making it virtually impossible for experimental results to be reproduced. An alternative approach is to utilize synthetic femurs of standardized geometry, having material behavior approximating that of human bone, but with very small specimen to specimen variability. This approach allows for repeatable experimental results and a standard geometry for use in accompanying FE models. While the synthetic bones appear to be of appropriate geometry to simulate bone mechanical behavior, it has not, however, been established what bone quality they most resemble, i.e., osteoporotic or osteopenic versus healthy bone. Furthermore, it is also of interest to determine whether FE models of synthetic bones, with appropriate adjustments in input material properties or geometric size, could be used to simulate the mechanical behavior of a wider range of bone quality and size. To shed light on these questions, the axial and torsional stiffness of cadaveric femurs were compared to those measured on synthetic femurs. A FE model, previously validated by the authors to represent the geometry of a synthetic femur, was then used with a range of input material properties and change in geometric size, to establish whether cadaveric results could be simulated. Axial and torsional stiffnesses and rigidities were measured for 25 human cadaveric femurs (simulating poor bone stock) and three synthetic "third generation composite" femurs (3GCF) (simulating normal healthy bone stock) in the midstance orientation. The measured results were compared, under identical loading conditions, to those predicted by a previously validated three-dimensional finite element model of the 3GCF at a variety of Young's modulus values. A smaller FE model of the 3GCF was also created to examine the effects of a simple change in bone size. The 3GCF was found to be significantly stiffer (2.3 times in torsional loading, 1.7 times in axial loading) than the presently utilized cadaveric samples. Nevertheless, the FE model was able to successfully simulate both the behavior of the 3GCF, and a wide range of cadaveric bone data scatter by an appropriate adjustment of Young's modulus or geometric size. The synthetic femur had a significantly higher stiffness than the cadaveric bone samples. The finite element model provided a good estimate of upper and lower bounds for the axial and torsional stiffness of human femurs because it was effective at reproducing the geometric properties of a femur. Cadaveric bone experiments can be used to calibrate FE models' input material properties so that bones of varying quality can be simulated.     

A visco-hyperelastic constitutive approach for modeling polyvinyl alcohol sponge

This study proposes the quasi-linear viscoelastic (QLV) model to characterize the time dependent mechanical behavior of poly(vinyl alcohol) (PVA) sponges. The PVA sponges have implications in many viscoelastic soft tissues, including cartilage, liver, and kidney as an implant. However, a critical barrier to the use of the PVA sponge as tissue replacement material is a lack of sufficient study on its viscoelastic mechanical properties. In this study, the nonlinear mechanical behavior of a fabricated PVA sponge is investigated experimentally and computationally using relaxation and stress failure tests as well as finite element (FE) modeling. Hyperelastic strain energy density functions, such as Yeoh and Neo-Hookean, are used to capture the mechanical behavior of PVA sponge at ramp part, and viscoelastic model is used to describe the viscose behavior at hold part. Hyperelastic material constants are obtained and their general prediction ability is verified using FE simulations of PVA tensile experiments. The results of relaxation and stress failure tests revealed that Yeoh material model can define the mechanical behavior of PVA sponge properly compared with Neo-Hookean one. FE modeling results are also affirmed the appropriateness of Yeoh model to characterize the mechanical behavior of PVA sponge. Thus, the Yeoh model can be used in future biomechanical simulations of the spongy biomaterials. These results can be utilized to understand the viscoelastic behavior of PVA sponges and has implications for tissue engineering as scaffold.     

Sensitivity of the tibio-femoral response to finite element modeling parameters

A generic finite element (FE) model of the lower limb was used to study the knee response in-vivo during a one-legged hop. The approach uses an explicit FE code and a combination of estimated muscle forces and measured three-dimensional tibio-femoral kinematics and ground reaction force as input to the FE model. The sensitivity of the simulated tibio-femoral response to variations of key geometric and material parameters was investigated by performing a total of 38 different simulations. The amplitudes of both kinematic and kinetic responses were affected by the change of these parameters. For the current approach, the results suggest that while cartilage mechanical and geometric properties are very important for the estimation of tibio-femoral cartilage pressure, they have limited effects on the overall kinematic response. The study may help to better define the relative importance of modeling parameters for the development of subject-specific models.