Mechanical response comparison in an implant overdenture retained by ball attachments on conventional regular and mini dental implants: a finite element analysis

This study investigates the bone/implant mechanical responses in an implant overdenture retained by ball attachments on two conventional regular dental implants (RDI) and four mini dental implants (MDI) using finite element (FE) analysis. Two FE models of overdentures retained by RDIs and MDIs for a mandibular edentulous patient with validation within 6% variation errors were constructed by integrating CT images and CAD system. Bone grafting resulted in 2 mm thickness at the buccal side constructed for the RDIs-supported model to mimic the bone augmentation condition for the atrophic alveolar ridge. Nonlinear hyperelastic material and frictional contact element were used to simulate characteristic of the ball attachment-retained overdentures. The results showed that a denture supported by MDIs presented higher surrounding bone strains than those supported by RDIs under different load conditions. Maximum bone micro strains were up to 6437/2987 and 13323/5856 for MDIs/RDIs under single centric and lateral contacts, respectively. Corresponding values were 4429/2579 and 9557/5774 under multi- centric and lateral contacts, respectively. Bone micro strains increased 2.06 and 1.96-folds under single contact, 2.16 and 2.24-folds under multiple contacts for MDIs and RDIs when lateral to axial loads were compared. The maximum RDIs and MDIs implant stresses in all simulated cases were found by far lower than their yield strength. Overdentures retained using ball attachments on MDIs in poor edentulous bone structure increase the surrounding bone strain over the critical value, thereby damaging the bone when compared to the RDIs. Eliminating the occlusal single contact and oblique load of an implant-retained overdenture reduces the risk for failure.     

Influence of transmucosal height in abutments of single and multiple implant-supported prostheses: a non-linear three-dimensional finite element analysis

The aim of this study was to analyze the influence of three different transmucosal heights of the abutments in single and multiple implant-supported prostheses through the finite element method. External hexagon implants, MicroUnit, and EsthetiCone abutments were scanned and placed in an edentulous maxillary model obtained from a tomography database. The simulations were divided into two groups: (1) one implant with 3.75 × 10 mm placed in the upper central incisor, simulating a single implant-supported fixed prosthesis with an EsthetiCone abutment; and (2) two implants with 3.75 × 10 mm placed in the upper lateral incisors with MicroUnit abutments, simulating a multiple implant-supported prosthesis. Subsequently, each group was subdivided into three models according to the transmucosal height (1, 2, and 3 mm). A static oblique load at an angle of 45 degrees to the long axis of the implant in palatal-buccal direction of 150 and 75 N was applied for multiple and single implant-supported prosthesis, respectively. The implants and abutments were assessed according to the equivalent Von Mises stress analyses while the bone and ceramics were analyzed through maximum and minimum principal stresses. The total deformation values increased in all models, while the transmucosal height was augmented. The transmucosal height of the abutments influences the stress values at the bone, ceramics, implants, and abutments of both the single and multiple implant-supported prostheses, with the transmucosal height of 1 mm showing the lowest stress values.     

Finite element analysis of implant-supported prosthesis with pontic and cantilever in the posterior maxilla

The aim of this study was to evaluate the influence of pontic and cantilever designs (mesial and distal) on 3-unit implant-retained prosthesis at maxillary posterior region verifying stress and strain distributions on bone tissue (cortical and trabecular bones) and stress distribution in abutments, implants and fixation screws, under axial and oblique loadings, by 3D finite element analysis. Each model was composed of a bone block presenting right first premolar to the first molar, with three or two external hexagon implants (4.0 × 10 mm), supporting a 3-unit splinted dental fixed dental prosthesis with the variations: M1 – three implants supporting splinted crowns; M2 – two implants supporting prosthesis with central pontic; M3 – two implants supporting prosthesis with mesial cantilever; M4 – two implants supporting prosthesis with distal cantilever. The applied forces were 400 N axial and 200 N oblique. The von Mises criteria was used to evaluate abutments, implants and fixation screws and maximum principal stress and microstrain criteria were used to evaluate the bone tissue. The decrease of the number of implants caused an unfavorable biomechanical behavior for all structures (M2, M3, M4). For two implant-supported prostheses, the use of the central pontic (M2) showed stress and strain distributions more favorable in the analyzed structures. The use of cantilever showed unfavorable biomechanical behavior (M3 and M4), mainly for distal cantilever (M4). The use of three implants presented lower values of stress and strain on the analyzed structures. Among two implant-supported prostheses, prostheses with cantilever showed unfavorable biomechanical behavior in the analyzed structures, especially for distal cantilever.     

Stress distribution of three-unit fixed partial prostheses (conventional and pontic) supported by three or two implants: 3D finite element analysis of ductile materials

In implantology, when financial or biological feasibility limitations appear, it is necessary to use prostheses with geometries that deviate from the conventional, with a pontic in the absence of an intermediate implant. The aim of this study was analyze and understand the general differences in the stresses generated in implants, components and infrastructures according to the configuration of the prosthesis over three or two implants. Thus, this paper analyzes the von Mises equivalent stresses (VMES) of ductile materials on their external surfaces. The experimental groups: Regular Splinted Conventional Group (RCG), which had conventional infrastructures on 3 regular-length Morse taper implants (4x11?mm); Regular Splinted Pontic Group (RPG), which had infrastructures with intermediate pontics on 2 regular-length Morse taper implants (4x11?mm). The simulations of the groups were created with Ansys Workbench 10.0 software. The results revealed that the RPG presented greater areas of possible fragility due to higher stress concentrations, for example, in the cervical area of the union between the implant and component the top platform of the abutment, as well as greater coverage of the stress by the cervical implant threads. The RPG infrastructure was also more affected by stresses in the connection areas between the prostheses and on the occlusal surface. There is an advantage to using prostheses supported by a greater number of implants (RCG) because this decreases the stress in the analyzed structures and consequently improves stress dissipation to the supporting bone, which would preserve the system.     

Finite element analysis for transverse carpal ligament tensile strain and carpal arch area

Mechanics of carpal tunnel soft tissue, such as fat, muscle and transverse carpal ligament (TCL), around the median nerve may render the median nerve vulnerable to compression neuropathy. The purpose of this study was to understand the roles of carpal tunnel soft tissue mechanical properties and intratunnel pressure on the TCL tensile strain and carpal arch area (CAA) using finite element analysis (FEA). Manual segmentation of the thenar muscles, skin, fat, TCL, hamate bone, and trapezium bone in the transverse plane at distal carpal tunnel were obtained from B-mode ultrasound images of one cadaveric hand. Sensitivity analyses were conducted to examine the dependence of TCL tensile strain and CAA on TCL elastic modulus (0.125–10 MPa volar-dorsally; 1.375–110 MPa transversely), skin-fat and thenar muscle initial shear modulus (1.6–160 kPa for skin-fat; 0.425–42.5 kPa for muscle), and intratunnel pressure (60–480 mmHg). Predictions of TCL tensile strain under different intratunnel pressures were validated with the experimental data obtained on the same cadaveric hand. Results showed that skin, fat and muscles had little effect on the TCL tensile strain and CAA changes. However, TCL tensile strain and CAA increased with decreased elastic modulus of TCL and increased intratunnel pressure. The TCL tensile strain and CAA increased linearly with increased pressure while increased exponentially with decreased elastic modulus of TCL. Softening the TCL by decreasing the elastic modulus may be an alternative clinical approach to carpal tunnel expansion to accommodate elevated intratunnel pressure and alleviate median nerve compression neuropathy.     

Three-dimensional numerical simulation of stress induced by different lengths of osseointegrated implants in the anterior maxilla

Lower survival rates were observed for the implant placed in the anterior maxilla. The purpose of this study was to investigate the influence of different implant lengths on the stress distribution around osseointegrated implants under a static loading condition in the anterior maxilla using a three-dimensional finite element analysis. The diameter of 4.0 mm external type implants of different lengths (8.5 mm, 10.0 mm, 11.5 mm, 13.0 mm, 15.0 mm) was used in this study. The anterior maxilla was assumed to be D3 bone quality. All the material was assumed to be homogenous, isotropic and linearly elastic. The implant–bone interface was constructed using a rigid element for simulating the osseointegrated condition. Then, 176 N of static force was applied on the middle of the palatoincisal line angle of the abutment at a 120°angle to the long axis of abutment. The von Mises stress value was measured with an interval of 0.25 mm along the bone–implant interface. Incremental increase in implant length causes a gradual reduction of maximum and average von Mises stress at the labial portion within the implant. In the bone, higher stress was concentrated within cortical bone area and more distributed at the labial cortex, while cancellous bone showed relatively low stress concentration and even distribution. An increase in implant length reduced stress gradients at the cortical peri-implant region. Implant length affects the mechanisms of load transmission to the osseointegrated implant. On the basis of this study the biomechanical stress-based performance of implants placed in the anterior maxilla improves when using longer implants.     

Finite element analysis of customized reconstruction plates for mandibular continuity defect therapy

Large mandibular continuity defects pose a significant challenge in oral maxillofacial surgery. One solution to this problem is to use computer-guided surgical planning and additive manufacturing technology to produce patient-specific reconstruction plates. However, when designing customized plates, it is important to assess potential biomechanical responses that may vary substantially depending on the size and geometry of the defect.     

Influence of bicortical techniques in internal connection placed in premaxillary area by 3D finite element analysis

The aim of study was to evaluate the stress distribution in implant-supported prostheses and peri-implant bone using internal hexagon (IH) implants in the premaxillary area, varying surgical techniques (conventional, bicortical and bicortical in association with nasal floor elevation), and loading directions (0°, 30° and 60°) by three-dimensional (3D) finite element analysis. Three models were designed with Invesalius, Rhinoceros 3D and Solidworks software. Each model contained a bone block of the premaxillary area including an implant (IH, Ø4 × 10 mm) supporting a metal-ceramic crown. 178 N was applied in different inclinations (0°, 30°, 60°). The results were analyzed by von Mises, maximum principal stress, microstrain and displacement maps including ANOVA statistical test for some situations. Von Mises maps of implant, screws and abutment showed increase of stress concentration as increased loading inclination. Bicortical techniques showed reduction in implant apical area and in the head of fixation screws. Bicortical techniques showed slight increase stress in cortical bone in the maximum principal stress and microstrain maps under 60° loading. No differences in bone tissue regarding surgical techniques were observed. As conclusion, non-axial loads increased stress concentration in all maps. Bicortical techniques showed lower stress for implant and screw; however, there was slightly higher stress on cortical bone only under loads of higher inclinations (60°).     

Finite element calculated uniaxial apparent stiffness is a consistent predictor of uniaxial apparent strength in human vertebral cancellous bone tested with different boundary conditions

Strong correspondence between the uniaxial apparent strength and stiffness of cancellous bone allows the use of stiffness as a predictor of bone strength. Measured values of mechanical properties in cancellous bone can be different between experiments due to different experimental conditions. In the current study, bone volume fraction, experimentally determined and finite element (FE) predicted stiffness were examined as predictors of cancellous bone ultimate strength in two different groups each of which was tested using a different end constraint. It is demonstrated that, although always significant, the relationships of strength with bone volume fraction and experimentally determined stiffness are different between test groups. Apparent stiffness, estimated by FE modeling, predicts the ultimate strength of human cancellous bone consistently for all examined experimental protocols.     

Finite element analysis and computed tomography based structural rigidity analysis of rat tibia with simulated lytic defects

Finite element analysis (FEA), CT based structural rigidity analysis (CTRA) and mechanical testing is performed to assess the compressive failure load of rat tibia with simulated lytic defects.     

Finite element analysis of trabecular bone structure: a comparison of image-based meshing techniques

In this study, we investigate if finite element (FE) analyses of human trabecular bone architecture based on 168 μm images can provide relevant information about the bone mechanical characteristics. Three human trabecular bone samples, one taken from the femoral head, one from the iliac crest, and one from the lumbar spine, were imaged with micro-computed tomography (micro-CT) using a 28 μm resolution. After reconstruction the resolution was coarsened to 168 μm. First, all reconstructions were thresholded and directly converted to FE-models built of hexahedral elements. For the coarser resolutions of two samples, this resulted in a loss of trabecular connections and a subsequent loss of stiffness. To reduce this effect, a tetrahedral element meshing based on the marching cubes algorithm, as well as a modified hexahedron meshing, which thresholds the image such that load carrying bone mass is preserved, were employed. For each sample elastic moduli and tissue Von Mises stresses of the three different 168 μm models were compared to those from the hexahedron 28 μm model. For one sample the hexahedron meshing at 168 μm produced excellent results. For the other two samples the results obtained from the hexahedral models at 168 μm resolution were poor. Considerably better results were attained for these samples when using the mass-compensated or tetrahedron meshing techniques. We conclude that the accuracy of the FE-models at 168 μm strongly depends on the bone morphology, in particular its trabecular thickness. A substantial loss of trabecular connections during the hexahedron meshing process indicates that poor FE results will be obtained. In this case the tetrahedron or mass-compensated hexahedron meshing techniques can reduce the loss of connections and produce better results than the plain hexahedron meshing techniques.     

Comparison of stress distribution in dental crown with different cusp angles: 3D finite element analysis

The purpose of this study is to measure the failure risk of a crown depending on the cusp angle. Three all-ceramic crown models consisting of C_H (high incline), C_M (middle incline), and C_L (low incline) are designed. Stress is applied to the crown with Loading case-1 (top of cusp tip) and Loading case-2 (middle of cusp ridge) with the use of FEA software. In Loading case-1 and case-2, the C_H showed the highest Maximum Principal Stress (MPS) while the CL showed the lowest MPS. The cusp angle is an influential factor affecting stress distribution in dental crowns.     

Finite element analysis of stress distribution in intact and porcelain veneer restored teeth

The aim of this study was to investigate the stress distribution generated in a veneer restoration of an upper central incisor compared to intact teeth using the finite element analysis after applying a lingual buccal loading at the incisal edge. Methods: Two models were developed: one model contained enamel, dentine, cementum, periodontal ligament, cortical and trabecullar bones, and the other model was a veneer restoration; both models were developed using MSC/Nastran software (MacNeal-Schwendler Corporation, Los Angeles, CA, USA) as the pre- and post-processor. A 10-N load was applied at the incisal edge from the lingual to the buccal side to simulate oral conditions in this area (protrusion). Results: Von Mises stresses were then analysed for three different regions: A-B (enamel elements under the veneer or second enamel layer), A′-B′ (buccal enamel and/or veneer element layer) and C-D (lingual enamel elements layer). A higher stress mode was observed for both models at the lingual cervical region. Conclusions: The presence of a veneer restoration on the incisors is a good alternative to mimic the behaviour of enamel under protrusion loading conditions. The use of veneers to replace enamel during rehabilitations is recommended.     

Finite element analysis of the mitral apparatus: annulus shape effect and chordal force distribution

This study presents a three-dimensional finite element model of the mitral apparatus using a hyperelastic transversely isotropic material model for the leaflets. The objectives of this study are to illustrate the effects of the annulus shape on the chordal force distribution and on the mitral valve response during systole, to investigate the role of the anterior secondary (strut) chordae and to study the influence of thickness of the leaflets on the leaflets stresses. Hence, analyses are conducted with a moving and fixed saddle shaped annulus and with and without anterior secondary chordae. We found that the tension in the secondary chordae represents 31% of the load carried by the papillary muscles. When removing the anterior secondary chordae, the tension in the primary anterior chordae is almost doubled, the displacement of the anterior leaflet toward the left atrium is also increased. The moving annulus configuration with an increasing annulus saddle height does not give significant changes in the chordal force distribution and in the leaflet stress compared to the fixed annulus. The results also show that the maximum principle stresses in the anterior leaflet are carried by the collagen fibers. The stresses calculated in the leaflets are very sensitive to the thickness employed.     

Micro-CT and micro-FE analysis of pedicle screw fixation under different loading conditions

Anchorage of pedicle screw instrumentation in the elderly spine with poor bone quality remains challenging. In this study, micro finite element (µFE) models were used to assess the specific influence of screw design and the relative contribution of local bone density to fixation mechanics. These were created from micro computer tomography (µCT) scans of vertebras implanted with two types of pedicle screws, including a full region-or-interest of 10 mm radius around each screw, as well as submodels for the pedicle and inner trabecular bone of the vertebral body. The local bone volume fraction (BV/TV) calculated from the µCT scans around different regions of the screw (pedicle, inner trabecular region of the vertebral body) were then related to the predicted stiffness in simulated pull-out tests as well as to the experimental pull-out and torsional fixation properties mechanically measured on the corresponding specimens. Results show that predicted stiffness correlated excellently with experimental pull-out strength (R² > 0.92, p < .043), better than regional BV/TV alone (R2 = 0.79, p = .003). They also show that correlations between fixation properties and BV/TV were increased when accounting only for the pedicle zone (R² = 0.66–0.94, p ≤ .032), but with weaker correlations for torsional loads (R² < 0.10). Our analyses highlight the role of local density in the pedicle zone on the fixation stiffness and strength of pedicle screws when pull-out loads are involved, but that local apparent bone density alone may not be sufficient to explain resistance in torsion.     

Stress analysis of irradiated human tooth enamel using finite element methods

The objectives of this project were to use finite element methods to determine how changes in the elastic modulus due to oral cancer therapeutic radiation alter the distribution of mechanical stresses in teeth and to determine if observed failures in irradiated teeth correlate with changes in mechanical stresses. A thin slice section finite element (FE) model was constructed from micro CT sections of a molar tooth using MIMICS and 3-Matic software. This model divides the tooth into three enamel regions, the dentin-enamel junction (DEJ) and dentin. The enamel elastic modulus was determined in each region using nano indentation for three experimental groups namely – control (non-radiated), in vitro irradiated (simulated radiotherapy following tooth extraction) and in vivo irradiated (extracted subsequent to oral cancer patient radiotherapy) teeth. Physiological loads were applied to the tooth models at the buccal and lingual cusp regions for all three groups (control, in vitro and in vivo). The principal tensile stress and the maximum shear stress were used to compare the results from different groups since it has been observed in previous studies that delamination of enamel from the underlying dentin was one of the major reasons for the failure of teeth following therapeutic radiation. From the FE data, we observed an increase in the principal tensile stress within the inner enamel region of in vivo irradiated teeth (9.97 ± 1.32 MPa) as compared to control/non-irradiated teeth (8.44 ± 1.57 MPa). Our model predicts that failure occurs at the inner enamel/DEJ interface due to extremely high tensile and maximum shear stresses in in vivo irradiated teeth which could be a cause of enamel delamination due to radiotherapy.     

Effect of screw position on load transfer in lumbar pedicle screws: a non-idealized finite element analysis

Angled screw insertion has been advocated to enhance fixation strength during posterior spine fixation. Stresses on a pedicle screw and surrounding vertebral bone with different screw angles were studied by finite element analysis during simulated multidirectional loading. Correlations between screw-specific vertebral geometric parameters and stresses were studied. Angulations in both the sagittal and axial planes affected stresses on the cortical and cancellous bones and the screw. Pedicle screws pointing laterally (vs. straight or medially) in the axial plane during superior screw angulation may be advantageous in terms of reducing the risk of both screw loosening and screw breakage.     

Analysis of pelvic strain in different gait configurations in a validated cohort of computed tomography based finite element models

The pelvis functions to transmit upper body loads to the lower limbs and is critical in human locomotion. Semi-automated, landmark-based finite element (FE) morphing and mapping techniques eliminate the need for segmentation and have shown to accelerate the generation of multiple specimen-specific pelvic FE models to enable the study of pelvic mechanical behaviour. The purpose of this research was to produce an experimentally validated cohort of specimen-specific FE models of the human pelvis and to use this cohort to analyze pelvic strain patterns during gait. Using an initially segmented specimen-specific pelvic FE model as a source model, four more specimen-specific pelvic FE models were generated from target clinical CT scans using landmark-based morphing and mapping techniques. FE strains from the five models were compared to the experimental strains obtained from cadaveric testing via linear regression analysis, (R² values ranging from 0.70 to 0.93). Inter-specimen variability in FE strain distributions was seen among the five specimen-specific pelvic FE models. The validated cohort of specimen-specific pelvic FE models was utilized to examine pelvic strains at different phases of the gait cycle. Each validated specimen-specific FE model was reconfigured into gait cycle phases representing heel-strike/heel-off and midstance/midswing. No significant difference was found in the double-leg stance and heel-strike/heel-off models (p = 0.40). A trend was observed between double-leg stance and midstance/midswing models (p = 0.07), and a significant difference was found between heel-strike/heel-off models and midstance/midswing models (p = 0.02). Significant differences were also found in comparing right vs. left models (heel-strike/heel-off p = 0.14, midstance/midswing p = 0.04).