1. Introduction
Overdentures are one of the recognized methods of prosthetic treatment of the mandible. Relatively reasonable costs and uncomplicated clinical management with a significant improvement in retention and stabilization make this type of restoration an attractive treatment option for both patients and dentists [
1,
2,
3]. Precision elements, such as bars, ball or locator clips, telescopic crowns, and magnets [
4,
5,
6,
7], limit denture mobility both vertically and horizontally. Precision attachments used with implant-retained overdenture counteract free rotation, so a greater occlusal force may be exerted because of the perceived better denture stability.
From the biomechanical point of view, overdentures are the optimal solution because they allow for a more uniform physiological distribution of chewing forces and their impact on the mucosa and alveolar process. Shortened teeth allow for a greater axial load. A variety of possible structural solutions of overdentures, concerning both the number and location of dental implants, as well as the types of clips, make rational justification for individual component application difficult. Despite the abundant amount of works in the literature, supraconstruction designs are often based on clinical tradition without addressing available theoretical background [
8]. Consequently, along with other manufacturing errors or lack of follow-up care, unexpected and uncontrolled rotation of the denture, rapid wear, damage to its components, or implant overloading may happen [
9].
The number of implants and placement sites depend on many factors, including anatomical conditions and the economic situation of the patient. Currently, many clinicians are inclined to use two implants. Opinions are more divided about the choice of a precision system. Two basic solutions include bar retention and single ball or locator attachments [
10]. The choice of the precise element is an individual matter and depends on the conditions of the prosthetic area. When deciding on a bar retainer, we need to choose its proper shape and possible additional elements.
The modification of the bar is its extension, symmetrically, distally, beyond the implants. This results in better stability and reliable anchoring of the prosthesis to the bar. According to Mericske-Stern [
11], the length of the cantilever bar may be up to approximately 1 cm and should not exceed the area of the first premolar. In turn, Elsyad [
12] believes that the optimal length of the tabs is 7 mm. The smallest stresses are then generated without significant differences around the implants. Currently, instead of a cantilever bar, ball attachments are often used, placed on both sides, distally to the implants.
For a rational selection of the bar, suitable for a clinical situation, knowledge of its mechanical properties is indispensable to achieve long-term success in treatment [
13]. In order to analyze the range of motion occurring within a denture and its retention element, it is required to find out the mechanical and fatigue characteristics of each component.
The aim of the study is to compare the in vitro biomechanical properties of an overdenture retained by two bar-retained implants with an overdenture retained by two bar-retained implants with ball clips.
2. Materials and Methods
A model of a toothless mandible based on the FRASACO mold with two OSTEOPLANT implants 14 mm long, and 4 mm in diameter was prepared for the study. The implants were located in the anterior region. Parameters of implant elements are shown in
Table 1.
The first system included a “rider” bar type retention (PRECI-HORIX, CEKA) and the lower complete denture with a matrix (
Figure 1a). In this clinical situation, a denture is primarily retained directly on the mucosa, while the bar provides adequate retention. In the second system (see
Figure 1b), the “rider” bar type retention suprastructure was also used, and the ball clips (CEKA) were symmetrically placed in the distal part. A lower complete denture with the matrix on the bar and ball clip patrices was made, too (
Figure 2). In this system, components are the same as in the first one; however, the method of fixing the denture is different. It is retained directly on the mucosa, and the bar’s function is to stabilize it, while the clips used are to retain the whole system.
The geometry of the mandible and dentures was scanned with the Roland LPX 600 3D Laser Scanner (Roland DGA Corporation, Tulsa, OK, USA) and imported into the computing environment. In order to increase accuracy, three independent measurements were made from three different directions (see
Figure 2). Such a procedure eliminates places where the measuring head cannot register geometry during measurement from just one direction.
The data obtained by using a 3D (three-dimensional) scanner concerning, for instance, surface mesh, saved in an Standard Tessellation Language (STL) format, were transferred to the Geomagic Studio 12.0 software (Research Triangle Park, NC, USA) in which further processing was carried out. It consisted of matching and placing individual sets of measurement data on each other (3D scans) from three different directions so that one model with complete mandibular geometry (
Figure 3) could be obtained. During this process, all data were filtered, the discontinuities of geometry were eliminated, and all points outside the mandible were removed (e.g., accidental elements of the background, the measuring table, etc.). Any surface distortions that arise during scanning and data placement were also removed during data processing.
The obtained surface triangle mesh was subjected to further processing. Finite Elements Analysis (FEA) programs used for analysis require defining geometry as volumetric objects. For that, the next step involved changing the description of the mandibular mucosa geometry from surface triangle mesh to an object defined with Non-Uniform Rational Basis Spline (NURBS) type surfaces (see
Figure 4).
To generate NURBS surface, contours along with division lines of particular areas of the mandible were determined. These measures allowed the surface formation of the model to be controlled well.
Next, a surface mesh of the object was formulated, taking into account division lines. Depending on the complexity of local geometry, it is possible to control mesh density. The obtained surface of the mandibular mucosa defined with NURBS was exported by using a universal computer-aided design (CAD) data exchange file.
The other system components (including implants, conical abutments, bar, etc.) were modeled directly in the SolidWorks software—a CAD program, version 2013 (Dassault Systems Co., Waltham, MA, USA)—because, in the computed tomography (CT) scan and scanner image, metal elements produce a light artifact in the form of a flare, which makes the boundary of the object unclear. Geometric data of individual components were obtained from the manufacturer’s documentation and from our measurements. Spacing was positioned from CT image, scanner, and control measurements. To simplify numerical calculations, the mapping of threaded connections (outline of the helical thread line) was abandoned due to their negligible significance for the study. It would considerably prolong the calculation time, too. Ankylotic connection of the implants with the bone was assumed; that is, the thread connection was considered ideal and, hence, deprived of any mobility. Compact and spongy bone were matched to the model of the mandible with respect to the mucous membrane. A denture was developed with a 3D scanner, and data processing was analogous to the one done in the mandible modeling. The denture was adapted to the geometry of the mandible. An ideal adhesion of the denture to the prosthetic bearing area was assumed. All elements were positioned and geometrically related, as in reality (
Figure 5).
Parameters of the oral mucosa thickness and compact and spongy bone thickness were determined by following the literature data. For the mucous membrane, a thickness of 1.5 mm was assumed, Young’s modulus E = 5 MPa and a quite high Poisson’s ratio ν = 0.49, which to some extent imitated its incompressibility [
14,
15]. For the cortical bone, Young’s modulus E = 17 GPa was assumed; for spongy bone, E = 600 MPa with Poisson coefficient equal to 0.3 in both cases. Material features of the denture are described with E = 3500 MPa and ν = 0.35 data [
16,
17].
Values of basic material properties of individual elements determined based on the literature are shown in the table (see
Table 2) [
14,
15,
18,
19].
The prepared model was exported to the FEA calculation module. Finite element analysis enables evaluating the distribution of internal forces in a given object when exposed to an external load. In the loaded construction elements, there are two states: stress (loading force divided by a cross-section area, e.g., of an implant) and deformation, a change in an element geometry. When these two states are fused in one mathematical formula, a mathematical model is formulated, e.g., elasticity [
20,
21,
22].
FEA was performed in the SolidWorks Simulation calculation module of the SolidWorks program. FEA mesh parameters for systems I and II are given in
Table 3.
Forces of 20 and 50 N were applied at 0° and 45° angles, as well as 100 N at 0° angle [
16,
17,
23,
24]. These forces correspond to the average and maximum loads on the teeth in occlusion. It should be added that, in patients using removable dentures, lower forces are produced during chewing. The range of occlusive forces from 20 to 90 N allows us to obtain a satisfactory mechanical efficiency for the denture. The maximum acceptable values of occlusive forces can be up to 122 N [
25]. The tests simulated chewing loads when the denture is loaded not only in the vertical but also in the horizontal plane. In the first case, the force was applied to the right central incisor at an angle of 0° (P1); in the second case, the same tooth was loaded with the forward force at 45° (P2). In the third and fourth cases, the force was directed to the right canine. In the third case, the force was applied at 0° (P3); in the fourth one, at 45°, perpendicular to the tooth (P4). In the fifth, sixth, and seventh cases, the force was directed to the right first molar. In the fifth, the force was applied at 0° (P5); in the sixth, at 45°, laterally toward the cheek (P6); and in the seventh, at 45°, forward (P7).
The displacements were in three directions: XYZ and net forces, stress reduction according to von Mises, and deformations were investigated. Displacement is a change in a given point position, i.e., a difference between initial and final position, e.g., extension, bending, and distortion angle. Its units are meters (m) or degrees (°) [
20,
21,
22].
To speak about displacement, one has to determine an object (a system) state called an undeformed state. It is an arbitrary state where the position of all points of an object is known. It means that the following function is known (1):
where
is the position vector, Ω is a set of all object points, and P is a set of all space points in which an object is “suspended”. In practice, such a state occurs when no force is applied to an object. Considering other object states, it is possible to determine its points position
. In light of the above, displacement is expressed as the following subtraction (2):
It means that displacement is a vector field that is assigning each object point its displacement vector.
To evaluate stress, a modified von Mises criterion was applied. According to it, the material will be damaged when reduced stress exceeds tensile strength [
26]. Contour maps show the distribution of stresses obtained in computer simulations.
Deformation, or relative deformation, is the ratio of displacement to initial dimensions, for instance, due to tension:
where
l0 is the initial length, and
lk is the final length.
It is a dimensionless quantity (or %) [
22,
23,
24,
25,
26].
The obtained results were analyzed and graphically presented in tables, graphs, and contour maps.
3. Results
To present the results, a force load of 50 N was selected. Contour maps show places where the greatest stress is concentrated within the bar retention and clips, as well as around the implants (
Figure 6). The places where the denture undergoes the most significant displacements are also depicted (
Figure 7). Crucial results obtained in the tests, i.e., displacement of the denture and the stresses of the bar and implants, are shown in the graphs (
Figure 8,
Figure 9,
Figure 10,
Figure 11 and
Figure 12).
Greater dislocations of the denture were recorded in the first system, where the maximum value was 0.7 mm (
Figure 8). In the second system, the maximum displacement of the denture was 0.6 mm (
Figure 9). The largest displacements occur when an incisor is loaded. For forces of the same value, greater displacements occur for oblique forces. Regardless of the location of the load, the most significant dislocations occur in the retromolar pad region. With forces acting on an incisor, especially for oblique forces, displacements are found on both sides of the denture.
For forces acting on a canine and molar, the denture undergoes greater displacement on the balancing side. Depending on the shape of the precision element used, the value of displacement changes while the direction of rotation remains the same. Loading of the examined teeth results in rotational motion on the working side, pressing against the denture bearing area, and tilting from the YZ plane (~in the XY plane, ~around the Y-axis) on the unloaded side, which results in detaching the denture on this side. The use of CEKA clips results in extension of the retention element, which results in a more stable denture, less significant displacement, especially when loading an incisor and canine. Adverse effects of oblique forces are also smaller, and there is no rotation. Mucosal load onto the bone remains smaller, too. The use of CEKA clips makes the bar itself subjected to greater displacement, especially with forces directed to an incisor, which may bring about a larger bone load around the implants. When exposing the denture to a single load, no significant differences were found in the bone, both cortical and spongy, for individual bars. However, tests should also include cyclic loads.
CEKA clips are part of a system that is exposed to displacements, especially when loading a canine and incisor. It should be noted that oblique forces are always less favorable for a clip on the working side. However, when loading canine for the clip on the balancing side, greater displacements occur for vertical forces. The obtained data show that implants and conical abutments make very stable system elements, regardless of the bar used.
The weakest element of this part of the system is the place where the bar is fixed, i.e., the screws, especially when loading a molar with both a vertical force of 100 N and oblique forces. For these forces, greater displacements within the screws occur on the balancing side than on the working one. When a denture retained by a “rider” bar is loaded with a vertical force of 100 N, and the displacement of the screws exceeds the displacement of the denture. However, the use of CEKA clips significantly limits the screw displacement, leading to protection against loosening.
The greatest stress concentrate within the precision element and concern the bar. In the first system, bar stress ranges from 0.0020 to 27.552 MPa (
Figure 10), obtained when a canine is loaded with a vertical force of 100 N. High values of maximum stresses when loading all teeth with a vertical force of 100 N are, however, focal, while the average values are definitely smaller (3.8260 MPa for an incisor, 3.1607 MPa for a canine, and 1.1483 MPa for a molar). Regardless of the force location, higher stresses occur for oblique forces. The largest bar stresses appear in the medial area, especially for the incisor and canine. For the molar, stresses are distributed over a larger area. Moreover, during loading of the molar, there are clearly lower stresses in the bar than in an incisor or canine, and the force of 50 N does not exceed 1 MPa. In the second system, bar stresses range from 0.080 to 578.6 MPa (
Figure 11). The highest stresses occur when loading a molar tooth with a vertical force of 100 N. For an incisor and molar, higher stresses occur for oblique forces (for a force of 50 N, average bar stresses for an incisor are 8.1340 MPa for vertical force and 10.7370 MPa for oblique force; for a molar, they are 4.7082 MPa for vertical force and 5.4998 and 7.5317 MPa for oblique forces). Relatively small stresses arise in the bar when loading a canine with oblique forces (for the force of 50 N, the maximum value is 13.4736 MPa, and average 1.6251 MPa). In the second system, the precision element was extended by mounting the CEKA clips. The maximum stresses produced in clips are smaller than the stresses of the bar. Stress values range from 0.0538 to 136.9900 MPa for the right clip and from 0.0784 MPa to 127.6700 MPa for the left. The highest stresses occur in the right clip with a vertical force of 100 N applied on a canine and in the left clip when applying this force on an incisor. Significant maximum values of clip stresses result from their shape and are located within the narrowing. CEKA clips cause significantly higher stresses in the screws, conical abutments, and implants than when a bar is used alone.
Comparative analysis of the first and second system (
Figure 12) shows that the additional use of CEKA clips causes bar extension, and therefore the denture is less susceptible to displacement. Still, more stresses are generated within the bar. Both the value of the maximum and average stresses is several times greater. The distribution of stresses is changed, too. When a bar is used alone, the greatest stresses are limited, located in the medial area. CEKA clips result in a more even distribution of stresses in the bar, with a greater concentration of stresses in clips. Significant stresses in this area are also associated with the clip shape. The material narrowing site causes a notch effect. Interestingly, oblique forces applied to all teeth give rise to similar stresses in both clips. An exception is an oblique force that is directed forward to a molar and thus acting in the final phase of the chewing cycle, where higher stresses are on the balancing side. For vertical forces, the greatest stresses appear on both sides when loading an incisor. Loading a canine generates higher stresses on the working side, whereas loading a molar causes relatively smaller stresses, which are asymmetrical and greater on the balancing side.
For subsequent components of the precision element and implants, the stresses are smaller and smaller because each element has absorbed part of the energy created by the load. Of course, ideally, stresses for this part of the system would equal 0. However, the material applied for testing is real. These elements work, and the denture works; hence, there are low values of stress. It is essential that CEKA clips cause a significant increase in stresses within conical abutments and screws. Proportionally higher stresses in the screws when using CEKA clips result from the screws being more stable in this system, as they undergo lesser displacement. Hence, a smaller risk of loosening. However, an increase in stress may cause mechanical damage to be more likely. Higher stresses in conical abutments and screws generate lower stresses in the implants, too.
Moreover, the studies showed that significant stresses were observed for the matrix housing in the first and the second system. However, maximum values occurred only in the site marked in
Figure 13. In other places, stress values are smaller. The place where maximum stress occurs is the result of element geometry, causing stress buildup in this area. This place can be exposed to mechanical damage, but due to embedding in acrylic, it will be difficult to locate, but it does not necessarily affect the denture function.
High values of stress were found in the denture. Maximum values, however, occur only in the place where the force is applied. This situation is not clinical in nature, as under physiological conditions, the load is applied to the denture plane and is not point-wise. The average values of stress are incomparably smaller. Generally, there are no significant stresses in the denture as it undergoes dislocations rather than stresses when loaded.
The deformation results concerning each case show that the subsequent components undergo only minor deformations, which are not relevant for understanding the biomechanics of the overdenture. During the deformation analysis, only the mucosa and the matrix together with the housing deserve attention. Irrespective of the shape of the precision element, the highest mucosal deformations occur when an incisor is loaded (
Figure 14). In system I, the maximum deformation values are 0.22900 for a vertical force of 100 N and 0.15670 for an oblique force of 50 N. The use of CEKA clips additionally does not significantly change the deformation values, which are maximum at 0.21060 for a vertical force of 100 N and 0.15060 for an oblique force of 50 N. The results concerning the first and second system show that the mounting of CEKA clips does not significantly affect the deformation values of the matrix and its housing, causing only a slight increase, especially when loading the molar tooth (
Figure 15). The highest deformation values were recorded for the forces directed at the incisal tooth, and when only the bar is used, they are 0.0144 (matrix housing) and 0.0077 (matrix) for a vertical force of 100 N and 0.0057 and 0.0055, respectively, for an oblique force of 50 N. However, when the bar with CEKA clips is used, the maximum deformations are 0.0109 and 0.0081 for a vertical force of 100 N and 0.0066 and 0.0062, respectively, for an oblique force of 50 N.
Denture deformations occur mainly at the force application site. Similar deformation values were found for both systems. In system I, the largest deformations occur for vertical forces when the molar is loaded, namely 0.09744 (for 100 N force), and for oblique forces when the incisor is loaded, namely 0.0490 (for 50 N force). In system II, the largest deformations occur for vertical forces when loading the canine, namely 0.08876 (for 100 N force), and for oblique forces when loading the incisor, namely 0.1506 (for 50 N force).
Regardless of the shape of the precision element, the denture undergoes little deformation. Acrylic is a rigid material so that the denture is displaced under load but not deformed. As expected, the element that undergoes significant deformation is, due to its highest elasticity, the mucosa. The deformation values of the mucosa are related to the displacement of the denture. Another element where increased deformation was observed is the matrix, which is the second most elastic part. For both systems, however, the deformation values are already much smaller. Larger values of deformation of the matrix in relation to the other components of the model, related to its elasticity, are the cause of wear of this element under the influence of the use of the denture.
4. Discussion
A significant problem when planning an overdenture is the choice of an appropriate precision system. Many researchers consider dentures retained by constructions splinting to be an advantageous solution [
27,
28,
29]. Undoubtedly, the advantages of bar constructions include proper denture stabilization and the ability to regulate retention.
Loads during denture dislocations in the lateral region are reduced [
30]. The incidence of complications is also lower compared to ball attachments. Magnets, which are less frequently used, can corrode and allow for horizontal displacement of the denture [
31]. They also have the lowest retention strength. It should be remembered, however, that bars require more space occlusally, that is, between the jaws (13–14 mm), as compared to ball attachments (10–11 mm) [
32]. No adequate space may result in the incorrect shape of the denture plate or the position of artificial teeth. This, in turn, may bring about faulty facial aesthetics or impaired sound articulation [
33,
34,
35,
36,
37].
The studied bar has a retention shape. Currently, connections of round or straight cross-sections are less and less often used. These bars are only indicated for short implants or problems with osseointegration, as they generate fewer stresses around the implants [
38]. To improve retention, the “rider” bar is better. The tests have confirmed that the more retentive the shape of the precision element, the smaller the denture dislocations. The length of the connection is of great importance. The greater the length of the retention arm, the more rigid the denture. Higher rigidity of the prosthetic restoration causes smaller displacement of the restoration, so mucosal support of the denture is smaller as well. However, the stresses within the precision element and implants are greater. The site of implant insertion may determine the length of the retention arm. The wider the distance between implants, the longer the connection. However, it should be remembered that when the line connecting the implants exceeds the denture plate, only individual clips are possible [
37]. To extend the retention arm, it is possible to use symmetrical projections or clips distally. The results of the studies show that the use of additional clips reduces denture displacement. However, the bar with clips and implants is more loaded; the more retentive the shape of the precision element, the greater the stresses [
38,
39].
In the case of short implants and substantial bone atrophy, it seems more advisable to make a straight bar to increase the mucosal force transfer and relieve the implants. Studies showing that a longitudinal bar extended distally (with a twice-broken axis) causes a significant exceeding of the stress tolerated by the bone tissue can be found in the literature. Still, this study has not confirmed such a phenomenon. An increase of stress in this area was observed only, which in extreme cases cannot rule out bone resorption around the implant. A relatively smaller stress increase in the implants is related to the fact that conical abutments and screws absorb a greater part of the energy released from loading forces. This is where a marked increase in stress was noted. It seems, therefore, that the risk of loosening or even damage to these elements is greater than the disintegration of the implant or peri-implantitis [
40,
41]. The conducted research confirmed the reports of other authors that the greatest load occurs in the distal parts of the bar retention [
42,
43]. Stress concentration, however, within the clips results from their shape.
The implementation of additional clips is indicated for long implants in the absence of clinical problems with osseointegration. Then, the use of direct support on implants that are ankylotically connected to the bone relieves the structure of the denture bearing area, especially in case of a mucous membrane poor susceptibility. This contributes to protecting the underlying tissues, which are not physiologically adapted to sustain direct occlusal forces, which may lead to bony atrophy in a short time [
44,
45].
The studies presented in this paper show the interdisciplinary character of today’s prosthetics. FEA is applied by a growing number of researchers [
24,
38,
43,
46]. Advances in virtual engineering enable making computer simulations in concrete and stable conditions, which helps avoid mistakes in experimental studies. Of course, it would be prudent to perform an experimental validation to verify the result. Confirmation of the FEA ideally would be performed in human subjects, but that is very difficult to perform stress tests in vivo. For this test, a phantom can be prepared, but it should be noted that it is a very complicated biological model with heterogeneous materials. Nevertheless, in FEA studies on dental implants, validation techniques are still rare. In his study, Chang performed a literature review of finite element analysis of dental implants. In PubMed, he found 522 articles about FEA, but only 47 papers including validation [
47].
In the near future, the authors of this study plan to compare the results obtained in the presented article with clinical trials.