1. Introduction
The complex architecture of human bone, particularly the cancellous or spongy bone in the jaw, is a significant challenge to dental implant design and biomimetic engineering. As a result, the field of study has implemented advanced modeling techniques, with Voronoi lattice structures being widely preferred [
1,
2,
3,
4]. Notorious for their distinct open-cell structure that sets them apart from other lattice families like TPMS gyroid and sheet-type lattices, Voronoi lattices are recognized for their remarkable ability to precisely resemble the intricate and porous properties of human bone [
5,
6,
7,
8,
9,
10,
11]. Voronoi lattices are made up of a network of irregular, mathematically determined cell shapes that are very different in how they change shape and how well they can handle mechanical forces. Thus, making them extremely advantageous for scaffolding applications in the field of biomedical engineering. These scaffolds improve osseointegration and tissue ingrowth in addition to mimicking the natural porosity of bone, which usually varies from 60% to 90% [
12,
13,
14,
15,
16,
17]. Because of this, Voronoi lattices are necessary for making biomimetic implants because they offer a customized method that is better at copying the complex anatomical and functional details of human bone than other types of lattices [
18,
19,
20,
21,
22].
Similar to composite materials, bone has directional dependency under external loads [
23,
24,
25,
26,
27,
28], which is crucial yet challenging to accurately simulate using traditional computer-aided design (CAD) models [
29,
30,
31]. These conventional models, which frequently assume isotropic or orthotropic features, considerably simplify the complex behavior of human cancellous bone. When compared to actual clinical dental scenarios, this simplification significantly reduces the accuracy of computational simulations [
32]. The majority of earlier research has resorted to replacing human cancellous bone in finite element analysis with solid cross-section block CAD models, a technique that severely reduces the precision of numerical approximations [
31,
32,
33,
34,
35,
36,
37]. Voronoi lattice structures provide a remedy to this problem due to their capacity for customized design. They make it possible to adjust Young's modulus such that it more nearly resembles the modulus of common implant materials, including titanium and its alloys. This tailoring facilitates the mechanical synergy that is necessary for dental implants to integrate effortlessly and remain stable over time inside the human mandible. By strategically using Voronoi lattices, researchers can potentially guarantee a better level of accuracy and clinical relevance in implant design by bridging the gap between simplified assumptions and the complicated nature of biomechanical behavior.
The main focus of this work is to study biomimetic Voronoi lattice structures, which are unique in that they can mimic the porous properties of trabecular bone, which is an important part of human cancellous bone. The porosity of human cancellous bone, which usually ranges from 30% to 90%, is noteworthy and significant for its biomechanical function [
17,
38,
39,
40]. The Gibson-Ashby theory demonstrates the significance of porosity and relative density in defining the mechanical strength and stiffness of cancellous bone by relating the mechanical properties and physical structure of porous materials [
5,
41,
42,
43,
44,
45]. This theory demonstrates how the relative density of cancellous bone, which is a measure of its structural solidity against porosity, directly affects its biomechanical properties. This is how Voronoi lattice structures are different from other types of lattices, like sheet-type lattices, because they can perfectly replicate the unevenly distributed, networked porosity of cancellous bone. This feature improves the biomechanical reliability of finite element analysis (FEA) simulations in addition to promoting osseointegration.
The finite element method (FEM) is used in this study to improve the mechanobiology analysis of dental implants [
1,
21,
24,
29,
30,
31,
32,
34,
35,
36,
37,
44,
46,
47,
48,
49,
50,
51]. FEM is a well-known computer method for modeling and analyzing complex physical events by breaking a continuum into a few discrete parts. The use of FEM in biomechanics and biomimetics has transformed the study of cancellous bone structures by providing a comprehensive understanding of how they behave under different stress scenarios [
52,
53,
54,
55,
56,
57,
58]. In particular, this work uses finite element modeling (FEM) to simulate dynamic, slow loading conditions that are precisely designed to mimic the cyclic biting forces applied to dental implants during mastication. This robust approach differs from earlier research that mostly relied on static loading assumptions, ignoring the dynamic and cyclical nature of mandibular stress during chewing [
30,
55,
59,
60,
61,
62,
63,
64]. By using smooth step-time convergence, this study is more like the explicit dynamic loading that happens in clinical settings. It gives a more accurate and complete look at the stress, strain, and displacement that happen in dental implants. This improved modeling method provides insight into the mechanical behavior of the implants over time, which is a major advancement over previous static loading assumptions. Previous studies that have highlighted the shortcomings of static models in representing the complex biomechanical environment of the jaw [
32,
35]. These studies have advocated for the incorporation of dynamic loading parameters in order to more accurately simulate the physiological conditions of mastication.
This study carefully assesses how differences in pore size in biomimetic Voronoi lattice structures affect their capacity to bear external stresses, with particular emphasis on improving dental implant biomechanics. An adequate pore size not only increases the load-bearing capacity of the implant but also promotes an ideal biomechanical environment by reducing stress shielding effects and enabling uniform stress distribution [
66]. In order to acquire mechanical properties that match those of natural bone, pore size modifications are strategically used to control porosity and relative density [
67,
68,
69]. This promotes osseointegration and improves the long-term durability of the implant. The mechanical integrity and biological performance of implants are greatly influenced by pore sizes within a certain range [
70,
71,
72]. These studies show that increased porosity often results in decreased mechanical strength, even if it is advantageous for biological integration. Nonetheless, a tailored approach that maximizes the strength-to-weight ratio while maintaining biological functioning is made possible by the use of Voronoi lattice structures, which are in line with cancellous bone's inherent properties [
73].
By using finite element analysis (FEA), the main goal of this study is to improve the comprehension of the biomechanical properties of dental implants. In particular, it focuses on using biomimetic Voronoi lattice architectures with different pore diameters to effectively mimic human cancellous bone characteristics. In order to precisely replicate masticatory conditions seen in real life, every component of the dental implant is thoroughly examined in the present study. It also explores into effects to the implant parts and the trabecular bone mimic when the pore size changes in the Voronoi lattice. Thus, this work anticipates that the novel Voronoi lattice cancellous bone pore size changes have a significant impact on the lattice's relative density and porosity, which are required to replicate the integrity and strength of real cancellous bone. This work clarifies the crucial role that structural changes play in varying the overall biomechanical effectiveness of dental implant systems through a thorough investigation of biomechanical parameters—displacements, stresses, and strains—derived from FEA.
2. Materials and Methods
2.1. 2D CAD Drawing and Dimensions
The dental implant system's main parts—cancellous and cortical bone structures—were rigorously designed with the use of Creo Parametric software. As seen in
Figure 1, this method constructed both two-dimensional sketches and three-dimensional designed models with precise dimensions.
2.2. Pore Sizes and Voronoi Latticed bone
Compact bone is defined by a relatively low porosity, usually less than 15%, while certain instances show a porosity range of around 5–30% [
74]. On the other hand, trabecular bone, or cancellous bone, has much greater porosity; it often exceeds 70%; however, there are variances ranging from around 30% to over 90% [
75]. This trabecular structure has a relative density that varies between 0.05 and 0.3, indicating that it is very cellular [
76]. Building on the earlier reported information about the natural porosity of human trabecular bone, this study utilizes the nTopology software to design four biomimetic Voronoi lattice models with 1.0 mm, 1.5 mm, 2.0 mm, and 2.5 mm distinct pore diameters. The four biomimetic Voronoi lattice structures with pore sizes are designated as VTB10, VTB15, VTB20, and VTB25, respectively.
Figure 2 shows these four models with different pore sizes. Using finite element analysis (FEA), the study precisely investigates the impact of these various pore sizes in the Voronoi lattice structure on the stresses and strains that bones and other dental implant components experience. The objective of this methodology is to better understand the biomechanical effects of pore size variations on implant stability and bone stability.
2.3. Building, FE Volume Meshing, and Boundary Condition of Voronoi Lattice
A solid 3D CAD model saved as a Standard ACIS Text (SAT) file is imported and used as the design area to create Voronoi-latticed models that resemble cancellous bone structure. Then, using a randomization seed of 250 and four different pore sizes as shown in
Figure 2, a randomized lattice graph is produced (as shown in
Figure 3a). The lattice graph is thickened with a beam diameter of 0.3 mm to strengthen the structure. In order to reduce structural complexity and remove unnecessary overhangs, a lattice trimming procedure is used with a feature tolerance of 10. Through exercising precision, a biomimetic scaffold that closely mimics the architecture of cancellous bone could be designed, allowing for possible exploration of its biomechanical properties.
2.4. Oblique Load and FE Boundary Conditions
As shown in
Figure 4a, the occlusal crown surface was loaded with forces three different directions: mesiodistal, buccal-lingual, and apical. This is to illustrate the biomechanical response of the dental implant under multi-axial loading. There were three different specific force magnitudes applied: 23.4 N, 17.1 N, and 114.6 N. By using a multi-point constraint (MPC) approach, these forces converged at a dummy reference point that was 3 mm from the occlusal surface. This method replicated the complex forces encountered during clinical mastication by generating an equivalent force of 118.2 N, tilted at 75.8 degrees with respect to the occlusal plane [
32,
34,
35,
36,
77].
To be more specific, the simulation method included explicit dynamic loading at a slower rate to mimic the cyclic mastication stress seen in real life. Notably, to accurately mimic in vivo situations, the load was applied over a 0.5-second period in the dynamic scenario, representing the typical mastication frequency of 2 Hz [
32,
34,
78]. Specially designed boundary conditions were used to mimic real-life constraints. At the implant-mandible interface, a six-degree-of-freedom (DoF) Encastre boundary condition was used to fix the implant in all three spatial dimensions (X, Y, and Z). The objective of this configuration was to accurately replicate the distribution of combined stress in the mesiodistal, buccal-lingual, and apical directions. Finite Element (FE) Boundary by Body Modeling was used to represent the biomimetic Voronoi cancellous bone structure. As shown in
Figure 2, this technique made it easier to accurately export nodes from nTop for boundary condition assignment in the ABAQUS software environment. All dental implant constituents were investigated, including the artificial Voronoi lattice cancellous bone (
Figure 4b). The study accurately measured the compressive and tensile stresses in the implant system using finite element analysis with ABAQUS. This makes progress in the study of implant biomechanics under physiologically important loading scenarios.
2.5. Physical properties and FE Mesh
It is very important to mesh the model properly so that there are fewer numerical errors and the results of finite element analysis (FEA) are more reliable when compared to real-world events. Using nTopology software, the Voronoi cancellous bone meshing was done consecutively in the present study. First, an implicit body was constructed
from a traditional 3D solid CAD model of spongy bone that was imported from Creo Parametric software in order to utilize it as a design space. This implicit design area was then used for building a Voronoi lattice, which was then utilized for lattice meshing. Robust tetrahedral meshing was used after surface re-meshing steps to address the structure's complexity.
Figure 5 serves as an illustration of the choice of an edge length of 1.0 mm and a feature tolerance of 0.25 mm for this purpose. To prepare the model to use finite element analysis, the last step was performing FE volume meshing. To be more precise, several element sizes were used to mesh the dental implant components: 0.2 mm for the implant, 0.15 mm for the screw-retaining components, and 0.35 mm for the remaining components [
32]. The meshing procedure was carried out using ABAQUS software (see
Figure 5). The physical characteristics of the implant and bone materials are presented in
Table 2. It is assumed that the components of the implant, comprising the screw, abutment, and crown, have homogeneous, elastic characteristics.
Figure 5.
Detailed meshing of dental implant components and Voronoi lattice biomimetic bone, with emphasis on mesh configuration as illustrated in the zoomed section.
Figure 5.
Detailed meshing of dental implant components and Voronoi lattice biomimetic bone, with emphasis on mesh configuration as illustrated in the zoomed section.
Table 1.
FE Mesh statistics of different biomimetic trabecular bone and implant component models.
Table 1.
FE Mesh statistics of different biomimetic trabecular bone and implant component models.
|
Crown |
Abutment |
Screw |
Implant |
Cortical Bone |
Trabecular Bone |
Total |
Bone models |
No. Element |
No. Node |
No. Element |
No. Node |
No. Element |
No. Node |
No. Element |
No. Node |
No. Element |
No. Node |
No. Element |
No. Node |
No. Element |
No. Node |
VTB10 |
107,954 |
20,670 |
85,243 |
17,463 |
627,098 |
121,624 |
160,676 |
31,737 |
146,767 |
31,180 |
1,936,663 |
584,049 |
3,064,401 |
806,723 |
VTB15 |
107,954 |
20,670 |
85,243 |
17,463 |
627,098 |
121,624 |
160,676 |
31,737 |
146,767 |
31,180 |
1,187,637 |
363,867 |
2,315,375 |
586,541 |
VTB20 |
107,954 |
20,670 |
85,243 |
17,463 |
627,098 |
121,624 |
160,676 |
31,737 |
146,767 |
31,180 |
791,391 |
248,917 |
1,919,129 |
471,591 |
VTB25 |
107,954 |
20,670 |
85,243 |
17,463 |
627,098 |
121,624 |
160,676 |
31,737 |
146,767 |
31,180 |
591,166 |
186,214 |
1,718,904 |
408,888 |
Table 2.
Physical properties of materials for the FEA [
32,
35,
37,
51].
Table 2.
Physical properties of materials for the FEA [
32,
35,
37,
51].
Materials |
Young's Modulus E (MPa) |
Poisson's Ratio ν |
Density (g/cm3) |
Strength (MPa) |
Cortical bone |
Ex
|
12,600 |
νxy
|
0.3 |
1.79 |
190 |
Ey
|
12,600 |
νyz
|
0.253 |
Ez
|
19,400 |
νxz
|
0.253 |
|
|
νyx
|
0.3 |
|
|
νzy
|
0.39 |
|
|
νzx
|
0.39 |
Cancellous bone |
Ex
|
1148 |
νxy
|
0.055 |
0.45 |
10 |
Ey
|
210 |
νyz
|
0.01 |
Ez
|
1148 |
νxz
|
0.322 |
|
νyx
|
0.01 |
|
νzy
|
0.055 |
νzx
|
0.322 |
Gold abutment * |
136,000 |
0.37 |
17.5 |
765 |
Porcelain |
68,900 |
0.28 |
2.44 |
145 |
Titanium grade 4 * |
110,000 |
0.34 |
4.5 |
550 |
2.6. Relative density and Porosity of Biomimetic Voronoi Lattice
The nTopology (nTop) program is used to make a Voronoi-latticed biomimetic bone structure, which is different from traditional modeling methods. With this method, the surface area-to-volume ratio, relative density, and porosity of the lattice structure can all be accurately calculated. Understanding the mechanical characteristics and potential for osteointegration of biomimetic bones requires knowledge of these analyses. The following equations provide a mathematical basis for these key parameters:
Relative Density (
) of Biomimetic Bone:
Surface Area-to-Volume Ratio (SA/V) of Voronoi Lattice:
Porosity (Φ) Relating to Relative Density:
3. Results
3.1. Voronoi Lattice Pore Size
In this work, we investigated how different pore sizes affected the structural properties of biomimetic Voronoi lattice models of cancellous bone that were designed specifically for use as replicas in dental implants as human trabecular bone. Models with four different pore sizes—1 mm, 1.5 mm, 2 mm, and 2.5 mm—were created using nTopology software, all of which had a constant 0.4 mm beam thickness. The influence of these variations in pore size on a variety of essential lattice structure parameters, such as average beam length, number of beams, number of nodes, relative density, and surface area to volume ratio, was well established by our analysis (see
Table 3). The study found a significant pattern associated with increasing pore size. Notably, the average beam length rose as pore size increased, from 0.4614 mm at 1 mm to 0.9203 mm at 2.5 mm, showing the need for longer beams to cross the increasing distances between nodes. From the smallest to the biggest pore diameters, the lattice beam and node counts significantly decreased simultaneously with increasing pores, falling from 25809 to 3024 and 15186 to 2058, respectively. This reduction allows for bigger spaces of voids by signaling a decline in lattice complexity and density. Additionally, a clear decrease in relative density was seen, decreasing from 22.21% for the 1 mm pore size to 6.27% for the 2.5 mm pore size. This indicates that pore size and lattice density are inversely related. Larger pores seem to provide a less dense and more open structure, which may affect the scaffold's ability to mechanically support dental implants. Furthermore, when pore sizes increased, the surface area to volume ratio exhibited a slight rise, increasing from 10.49% to 12.4%. This increase suggests that bigger holes could provide somewhat better surface complexity in relation to volume, which might be advantageous for mechanical bonding and osteointegration with the surrounding bone tissue.
3.2. Pore size versus relative density and Porosity
The range of pore sizes in biomimetic Voronoi lattice structures, shown in
Figure 6, from 1.0 mm (VTB10) to 2.5 mm (VTB25), is an important factor for describing the pore scale and has a direct impact on the relative density and porosity of the scaffold. The relative density, which measures the lattice structure's compactness in comparison to its solid equivalent, steadily drops with increasing pore size; it declines from 22.21% at 1.0 mm to 6.27% at 2.5 mm. This trend represents a gradual decline toward a structure that is less dense. On the other hand, porosity measures the volume percentage of empty space within the lattice and has an inverse relationship with relative density. With increasing pore sizes,
Figure 6 shows an increasing porosity trend that ranges from 77.79% for VTB10 to 93.73% for VTB25. As pore size rises, this inverse dynamic shows the scaffold's structural and functional advancement, getting it closer to the porous architecture seen in real bone.
3.3. Dynamic oblique loading
Figure 7 displays the application of an oblique dynamic load of 118.2 N to the occlusal surface of the crown. The load is directed in the buccal-lingual, axial, and mesiodistal directions, simulating a 2 Hz mastication cycle with a plane of mastication load. The figure shows the smooth change in load on the crown's occlusal surface across mesiodistal, buccal-lingual, and apical orientations throughout a 0.5-second mastication cycle. The goals of this study were met by using multi-point constraints (MPC) and a reference point as the master control to create a smooth step amplitude on the occlusal surface of the crown. The load distribution was in the three designated directions, which matches the findings of a previous study [
32]. This proves that the masticatory forces on dental implants that were simulated in this research method are accurate.
3.4. von Mises stress in dental implant assembly
Figure 8 depicts the maximum von Mises stress values for a doubly sliced assembly of dental implant components, exhibiting stress distributions across four different biomimetic Voronoi lattice trabecular bone designs. The findings show a decreasing trend in von Mises stress with increasing pore size. The stress values are reported as 305.93 MPa for a 1 mm pore size (VTB10), subsequently decreasing to 220.96 MPa for 1.5 mm (VTB15), 186.01 MPa for 2 mm (VTB20), and ultimately to 161.16 MPa for 2.5 mm (VTB25). This trend emphasizes the impact of pore size on stress distribution inside the dental implant assembly, emphasizing the mechanical durability of implants with larger pore sizes under replicated loading conditions.
3.5. von Mises stress in dental implant and retaining screw
The Von Mises stress analysis, shown in
Figure 9a through 9h, illustrates the effect of various pore sizes inside biomimetic Voronoi latticed trabecular bones on stress distribution in both the dental implant and the retaining screw. The Von Mises stress for the dental implant consistently decreases as the pore size increases, from 1 mm (VTB10) to 2.5 mm (VTB25), falling from 233.84 MPa to 173.48 MPa, respectively. The retaining screw has an identical trend but a less noticeable decline from 220.96 MPa to 184.11 MPa. Interestingly, the 2 mm pore size (VTB20) indicates a key threshold where the implant's Von Mises stress significantly decreases to 179.69 MPa, which is in close alignment with the retaining screw's stress levels (182.09 MPa). This implies that 2 mm pores provide an optimal balance for stress distribution throughout the system, suggesting that it has the potential to be the ideal pore size for biomimetic Voronoi lattice trabecular bone in terms of improving the structural integrity of dental implants and their retention screws.
3.6. von Mises stress in cortical and Voronoi trabecular bone
The results from
Figure 10a–d, combined with observed stress concentration patterns, provide a comprehensive understanding of maximum von Mises stress behavior in both cortical and biomimetic Voronoi lattice trabecular bones under a mastication frequency of 2Hz, corresponding to a cycle time of 0.5 second. Notably, regions of possible mechanical risk are indicated by the concentration of the highest stress around the hole where the cortical bone contacts the implant. Stress distribution inside the Voronoi lattice trabecular bone is consistent throughout different pore sizes, ranging from 12.83 MPa at 1.0 mm (VTB10) to 25.55 MPa at 2.5 mm (VTB25). Likewise, stress in the cortical bone increases from 35.67 MPa (VTB10) to 106.42 MPa (VTB25) as trabecular structure pore size increases, indicating a clear relationship between pore size and stress effect on both types of bone. Stress levels peak at about 0.35 seconds into the mastication cycle, remaining constant even as cycles continue. This response illustrates the ability of biomimetic and cortical bone structures to adapt to masticatory pressures, underscoring the critical function of pore size in modifying stress distribution and maintaining bone structures' mechanical integrity under dynamic loading.
2.7. Micromotions in Voronoi trabecular bone
Figure 11, which depicts contour plots of the magnitude of micromotion inside the biomimetic Voronoi latticed trabecular bone subjected to finite element analysis, shows a clear relationship between pore size and micromotion magnitude. Specifically, micromotion values inside the biomimetic Voronoi trabecular bone increased gradually as pore diameters increased. For a pore size of 1.0 mm (VTB10), the micromotion was 8.00 micrometers. This micromotion increased progressively with pore size; at 1.5 mm (VTB15), it was reported at 10 micrometers, rising to 13.00 micrometers at 2.0 mm (VTB20), and reaching 17.00 micrometers at the maximum pore size of 2.5 mm (VTB25). These results show how important pore size is when studying the biomechanical behavior of the biomimetic Voronoi trabecular bone, especially micromotion when dynamic mastication stresses are present.
3.8. Displacement in dental implant assembly
Figure 12a–d show the maximum displacement found in the axial, mesiodistal, and buccolingual directions for various pore sizes of biomimetic Voronoi latticed trabecular bone. The findings show that biomimetic Voronoi latticed trabecular bones with varied pore sizes (VTB10, VTB15, VTB20, and VTB25) exhibit maximal displacements in the axial, mesiodistal, and buccolingual directions. The displacements in the axial direction range from 111.76 micrometers for VTB10 to 117.73 micrometers for VTB25. VTB20 has a mesiodistal displacement peak of 256.55 micrometers, whereas VTB25 has the largest buccolingual displacement at 310.07 micrometers. These studies show how pore size affects the biomechanical response of biomimetic trabecular bone, demonstrating how displacement changes across different directions and pore sizes, influencing the overall mechanical behavior of dental implant assemblies.
3.9. Reaction Forces in dental implant assembly
Figure 13a–d show the dynamic response forces for biomimetic Voronoi latticed trabecular bones with various pore sizes (VTB10, VTB15, VTB20, and VTB25) in the axial, mesiodistal, and buccolingual directions. Initially, response forces are zero across all models until 0.23 seconds, after which they increase in line with the dynamic load-ing amplitude. After 0.23 seconds, the maximum response forces are observed: VTB10 has forces of 529.41 KN axially, 858.6 KN mesiodistally, and 530.86 KN buccolingually; VTB15 has 465.09 KN, 927.27 KN, and 316.98 KN, respectively; VTB20 has 215.25 KN, 717.68 KN, and 259.8 KN; and VTB25 has 224.21 KN, 788.69 KN, and 329.24 KN in the respective directions. These findings demonstrate the substantial variation in biomechanical responses across the various pore-sized biomimetic implants during dynamic loading simulations.
4. Discussion
The findings we obtained support the argument that was proposed in the study's introduction. It is suggested that the variation in pore size inside biomimetic Voronoi lattice structures affects the stress distribution and micromotion across the dental implant components as well as the surrounding cancellous bone. According to our investigation, a pore size of 1.5 mm provides optimum equilibrium and closely mimics the dynamic loading response observed in human cancellous bone. This pore size makes the stress distribution significantly uniform at the implant-bone contact (BIC) by balancing structural factors like beam length and node count with the biomechanical needs of dental implants. With a relative density of 15.19% and a surface area-to-volume ratio that promotes cell attachment and nutrient transport, this option maintains structural strength while also facilitating osteointegration. The design effectively mimics the porosity structure of real bone, improving both mechanical stability and biological compatibility. Because it distributes stress analogously to naturally occurring cancellous bone, it can sustain physiological loads and lower the likelihood of stress shielding and implant failure. The lattice structure has about 84.81% porosity, which is similar to natural human cancellous bone, which usually has porosity levels between 75% and 90%. This shows that the chosen pore size worked well. This congruence verifies our study's biomimetic methodology and supports the integration and vascularization that are essential to dental implants' long-term effectiveness. As such, it highlights the possibility of tailored biomimetic Voronoi lattices to augment the lifetime and efficacy of dental implants. Thus, making them a model option for mimicking the intricate interplay between strength and porosity that characterizes human cancellous bone.
Figure 6 depicts the study's findings, which show a significant relationship between pore size and the mechanical integrity of biomimetic Voronoi lattice trabecular bone, particularly in the context of dental implants. Pore diameters increased from 1.0 mm to 2.5 mm, resulting in a decrease in relative density and an increase in porosity. These variations have an immediate influence on the implant's stress distribution properties under dynamic loading conditions. A lower relative density, which means bigger pores, is more like the natural structure of human cancellous bone. This makes it easier for the dental implant system to spread stress evenly. By mimicking the natural load-bearing function, this advancement promotes osseointegration and lowers the possibility of implant failure. The 2.0 mm (VTB20) configuration stands out as the best candidate among the assessed pore sizes as it maintains a compromise between preserving adequate structural strength with a relative density of 8.84% and having sufficient porosity (91.16%) for biological integration. This pore size successfully mimics the mechanical and structural characteristics of natural cancellous bone, which optimizes stress distribution and ensures the scaffold's endurance under physiological loading conditions, thereby extending the lifespan and functionality of dental implants.
Figure 8 shows that pore size variations in biomimetic Voronoi lattice trabecular bones have a significant impact on stress distribution across dental implant assemblies. particularly around critical interfaces consisting of implant-bone, implant-abutment, abutment-screw retaining, abutment-crown, and implant-retaining screw contacts. Stress concentrations are significantly reduced when the pore size is increased from 1.0 mm (VTB10) to 2.0 mm (VTB20), which allows for a more uniform stress distribution that is more consistent with the behavior of natural bones. Despite its intended purpose of mimicking human cancellous bone, the stiffness inherent in the smaller pore size (VTB10) assembly hinders stress homogeneity, boosting the possibility of localized failure. The larger pore size (VTB20) assembly makes the implant system more stable and helps it fit together better by improving its mechanical response and closely matching the load-bearing and stress-dissipation properties of natural cancellous bone. This optimization shows how important it is to change the size of the pores in Voronoi lattice structures to find the right balance between stiffness and porosity when trying to make a design that can withstand mechanical loads and mimic the biological activity of bone tissue.
Figure 9a–h show that the maximum von Mises stress inside dental implants and their retaining screws is significantly influenced by pore size in biomimetic Voronoi latticed trabecular bone designs. As the pore size increases from 1 mm (VTB10) to 2.5 mm (VTB25), the research shows a steady drop in von Mises stress values in both the implant and the retaining screw. Notably, the stress in the retaining screw decreased from 220.96 MPa to 184.11 MPa, and the stress in the dental implant decreased from 233.84 MPa for VTB10 to 173.48 MPa for VTB25. According to this trend, larger pores may help distribute stress uniformly, which could prolong the implant system's lifespan and improve its mechanical stability. Notably, the ideal stress distribution, defined by a compromise between mechanical integrity and biomimetic porosity, was seen for VTB20, demonstrating that a pore size of 2.0 mm is the best for imitating the mechanical characteristics of human cancellous bone. Stress concentration at the implant's neck and the screw's mid-region—areas that are often prone to failure or fracture under cyclic dynamic loading—is significantly reduced by this improved design. This adaptation mimics human bone's inherent load-bearing characteristics, with the potential to increase implant lifetime and function in clinical conditions.
Figure 10a–d show the maximum von Mises stress across cortical and biomimetic Voronoi lattice trabecular bones, emphasizing the critical role of pore size in the biomechanical response of different bone types under masticatory stress. As illustrated, the pore size directly affects the structural integrity and mechanical behavior of the bone-implant system. The even stress distribution inside the Voronoi trabecular structure and the stress concentration at the cortical bone-implant interface serve as evidence for this. The fact that von Mises stress goes up as pore sizes get bigger and stays the same after 0.35 seconds shows that structural design and dynamic loading response are complicatedly linked. In particular, changing the size of the pores in the Voronoi lattice trabecular bone creates a biomechanical environment that is similar to how bones naturally behave. This makes the stress distribution better and reduces stress concentrations in certain areas that could cause failure. This optimized distribution pattern not only demonstrates the biomimetic structure's capacity to imitate the physiological stress of human bones, but it also represents a major improvement in finite element analysis (FEA) results. The biomimetic Voronoi lattice is a better way to predict how bone-implant assemblies will act in real life. This is because it changes the size of the pores to find a balance between biological activity and mechanical support. This study uses biomimetic Voronoi latticed cancellous bone for finite element analysis, which is very different from earlier studies that used solid CAD blocks to model bone structures. By adding the Voronoi lattice's detailed, porous design, this study provides a more physiologically realistic description of bone mechanics, improving the predictability of stress distribution patterns under masticatory stresses. This method not only improves biomechanical modeling of bone-implant interfaces, but it also gives a more nuanced knowledge of the effect of pore size on structural integrity and mechanical behavior. As a result, this improvement opens the door for the development of bone scaffolds and implants with better clinical results, which are based on a more accurate replication of natural bone characteristics.
Figure 11 shows the sizes of micromotions in biomimetic Voronoi latticed trabecular bones with four different pore sizes. This shows how important micromotion is in the biomechanics of dental implants. Larger pore diameters are linked with increased micromotion, which increases the risk of bone resorption at the bone-implant interface and fractures in dental implant components, particularly the retaining screw and the implant itself. Additionally, micromotions increase the likelihood of corrosion, especially when the implant works in a biofluid environment. In spite of these difficulties, the research finds that a pore size of 2.0 mm (VTB20) with a micromotion of 13 micrometers provides an ideal compromise between reducing unfavorable biomechanical consequences and providing a satisfactory degree of mechanical rigidity. This pore size facilitates the structural and functional integration of the implant into the biological system while also reducing the chance of implant system failure. As such, these micromotion optimization insights provide an in-depth understanding of the relationship between pore size and implant function, which considerably advances the field of biomechanical dentistry research.
Incorporating the observed displacement patterns from
Figure 12a–d, the axial and buccolingual directions show both negative and positive peak displacement amplitudes, indicating that the dental implant can withstand cyclic masticatory stress. The scenario demonstrates the implant's flexibility and mechanical resilience, which are essential for enduring numerous stress cycles while maintaining structural integrity. Such displacement patterns indicate that VTB20, with its nuanced response in all measured directions, not only mimics the normal biomechanical behavior of human bone, but also guarantees the implant's capacity to tolerate dynamic stresses experienced in vivo. This modification emphasizes just how essential it is to choose the right pore size in the Voronoi lattice trabecular bone in order to maximize the performance of dental implants. VTB20 offers an excellent balance between mechanical rigidity and flexibility, which is crucial for long-term clinical success.
Figure 13a–d show the dynamic reaction forces across different pore sizes of biomimetic Voronoi latticed trabecular bone (VTB10, VTB15, VTB20, and VTB25), providing significant insight into the biomechanical performance of dental implant assemblies. The consistent lack of response forces for up to 0.23 seconds, followed by a rise in alignment with dynamic loading amplitudes, demonstrates the implants' mechanical flexibility to masticatory stresses. Specifically, the research shows that VTB20 outperforms the other models due to its lower axial and buccolingual reaction forces and optimum mesiodistal response. This research implies that a pore size of 2.0 mm promotes a compromise between structural support and flexibility, allowing the im-plant to efficiently distribute stresses while maintaining its biomechanical integrity. The distinct mechanical responses observed at the 0.5-second mark across all models highlight the importance of pore size in tailoring the biomechanical properties of dental implants to mimic natural bone behavior, thereby improving the overall effective-ness and longevity of dental implant systems.
This work advances the biomechanical knowledge of bone-implant interactions by using advanced finite element analysis (FEA) that extends beyond typical unrealistic models. Our technique provides a more realistic representation of bone behavior in physiological settings by replicating the nonlinear, and anisotropic features of biological tissues, as well as including dynamic loading conditions that replicate real-world pressures. This technique not only eliminates past shortcomings, but it also gives more in-depth insights into the performance and stability of biomimetic Voronoi lattice implants, making a substantial contribution to biomechanical engineering.
5. Conclusions
This work sought to understand the biomechanical behavior of biomimetic Voronoi latticed trabecular bone structures. Using finite element analysis (FEA) to evaluate the effect of altering pore sizes on stresses and micromotion magnitudes. Hence, the following are significant inferences from the present study:
A direct relationship was discovered between pore size and micromotion magnitude. This highlights the impact of pore size on the biomechanical behavior of biomimetic Voronoi trabecular bone. This relationship has significance during dental implant designs that optimize biomechanical performance.
The research shows that the biomimetic Voronoi latticed trabecular bone used in FEA has properties that are similar to real human cancellous bone when subjected to applied external dynamic mastication stresses.
The use of anisotropic characteristics in these novel biomimetic trabecular bones, as well as the simulation of dynamic mastication stresses, greatly improved the FEA findings. This method gives a more precise description of biomechanical interactions in the oral environment, resulting in more reliable implant designs.
The study found that a 2.0 mm (VTB20) pore size presents an optimal integration of mechanical rigidity and minimum micromotion, making it a better option to further enhance implant lifespan and performance.
Micromotion optimization considerably reduces the chances of implant failure, such as fractures in implant components, bone resorption, and corrosion in a biofluid environment.
The research gives a significant insight into the biomechanical optimization of dental implant designs by studying biomimetic Voronoi latticed trabecular bones. This study expands possibilities for the design of dental implants that enhance biomechanical compatibility and performance by using modern FEA methods and taking into account the anisotropic properties of bone and dynamic loading scenarios. These results demonstrate the significance of using biomimetic approaches for developing and evaluating the performance of dental implants, in addition to making advancements in the fields of biomechanical engineering in dentistry.
Author Contributions
Conceptualization, D.B.A, M.T, and S.J-H.; methodology, D.B.A, S.J-H, and M.T.; software, D.B.A, M.T, and S.J-H.; validation, D.B.A,S.J-H. and M.T.; formal analysis, D.B.A.; investigation, D.B.A.; data curation, D.B.A,; writing—original draft preparation, D.B.A ; writing—review and editing, D.B.A., M.T., S.J-H.; supervision, M.T. and S.J.H; project administration, S.J-H and M.T. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by the Ministry of Science and Technology of Taiwan (project No. MOST-105-2221-E-011-058-MY2).
Data Availability Statement
Not applicable.
Acknowledgments
The authors would like to acknowledge the financial support from the Min-istry of Science and Technology of Taiwan (project No. MOST-105-2221-E-011-058-MY2).
Conflicts of Interest
The authors declare no conflict of interest.
References
- Beaupre, G.S.; Hayes, W.C. Finite Element Analysis of a Three-Dimensional Open-Celled Model for Trabecular Bone. J. Biomech. Eng. 1985, 107, 249–256. [Google Scholar] [CrossRef] [PubMed]
- Bregoli, C.; Fiocchi, J.; Biffi, C.A.; Tuissi, A. Additively Manufactured Medical Bone Screws: An Initial Study to Investigate the Impact of Lattice-Based Voronoi Structure on Implant Primary Stability. Rapid Prototyp. J. 2024, 30, 60–72. [Google Scholar] [CrossRef]
- Alkentar, R.; Kladovasilakis, N.; Tzetzis, D.; Mankovits, T. Effects of Pore Size Parameters of Titanium Additively Manufactured Lattice Structures on the Osseointegration Process in Orthopedic Applications: A Comprehensive Review. Crystals 2023, 13, 113. [Google Scholar] [CrossRef]
- Sivakumar, N.K.; Palaniyappan, S.; Sekar, V.; Alodhayb, A.; Braim, M. An Optimization Approach for Studying the Effect of Lattice Unit Cell’s Design-Based Factors on Additively Manufactured Poly Methyl Methacrylate Cranio-Implant. J. Mech. Behav. Biomed. Mater. 2023, 141, 105791. [Google Scholar] [CrossRef] [PubMed]
- Du, Y.; Liang, H.; Xie, D.; Mao, N.; Zhao, J.; Tian, Z.; Wang, C.; Shen, L. Design and Statistical Analysis of Irregular Porous Scaffolds for Orthopedic Reconstruction Based on Voronoi Tessellation and Fabricated via Selective Laser Melting (SLM). Mater. Chem. Phys. 2020, 239, 121968. [Google Scholar] [CrossRef]
- Zhou, Y.; Isaksson, P.; Persson, C. An Improved Trabecular Bone Model Based on Voronoi Tessellation. J. Mech. Behav. Biomed. Mater. 2023, 148, 106172. [Google Scholar] [CrossRef]
- Zou, S.; Gong, H.; Gao, J. Additively Manufactured Multilevel Voronoi-Lattice Scaffolds with Bonelike Mechanical Properties. ACS Biomater. Sci. Eng. 2022, 8, 3022–3037. [Google Scholar] [CrossRef]
- Kanwar, S.; Al-Ketan, O.; Vijayavenkataraman, S. A Novel Method to Design Biomimetic, 3D Printable Stochastic Scaffolds with Controlled Porosity for Bone Tissue Engineering. Mater. Des. 2022, 220, 110857. [Google Scholar] [CrossRef]
- Li, J.; Yang, Y.; Sun, Z.; Peng, K.; Liu, K.; Xu, P.; Li, J.; Wei, X.; He, X. Integrated Evaluation of Biomechanical and Biological Properties of the Biomimetic Structural Bone Scaffold: Biomechanics, Simulation Analysis, and Osteogenesis. Mater. Today Bio 2024, 24, 100934. [Google Scholar] [CrossRef]
- Yang, W.; Han, Q.; Chen, H.; Li, Y.; Guo, X.; Zhang, A.; Liu, Y.; Sun, Y.; Wang, J. Additive Manufactured Trabecular-like Ti-6Al-4V Scaffolds for Promoting Bone Regeneration. J. Mater. Sci. Technol. 2024, 188, 116–130. [Google Scholar] [CrossRef]
- He, L.; Zhao, M.; Cheung, J.P.Y.; Zhang, T.; Ren, X. Gaussian Random Field-Based Characterization and Reconstruction of Cancellous Bone Microstructure Considering the Constraint of Correlation Structure. J. Mech. Behav. Biomed. Mater. 2024, 106443. [Google Scholar] [CrossRef] [PubMed]
- Bahati, D.; Bricha, M.; El Mabrouk, K. Vat Photopolymerization Additive Manufacturing Technology for Bone Tissue Engineering Applications. Adv. Eng. Mater. 2023, 25, 2200859. [Google Scholar] [CrossRef]
- Mamidi, N.; Ijadi, F.; Norahan, M.H. Leveraging the Recent Advancements in GelMA Scaffolds for Bone Tissue Engineering: An Assessment of Challenges and Opportunities. Biomacromolecules 2023. [Google Scholar] [CrossRef] [PubMed]
- Rahatuzzaman, M.; Mahmud, M.; Rahman, S.; Hoque, M.E. Design, Fabrication, and Characterization of 3D-Printed ABS and PLA Scaffolds Potentially for Tissue Engineering. Results Eng. 2024, 21, 101685. [Google Scholar] [CrossRef]
- Wu, J.; Zhang, Y.; Lyu, Y.; Cheng, L. On the Various Numerical Techniques for the Optimization of Bone Scaffold. Materials (Basel). 2023, 16, 974. [Google Scholar] [CrossRef]
- Zhang, Y.; Wei, D.; Chen, Y.; Xie, L.; Wang, L.; Zhang, L.-C.; Lu, W.; Chen, G. Non-Negligible Role of Gradient Porous Structure in Superelasticity Deterioration and Improvement of NiTi Shape Memory Alloys. J. Mater. Sci. Technol. 2024, 186, 48–63. [Google Scholar] [CrossRef]
- Zioupos, P.; Cook, R.B.; Hutchinson, J.R. Some Basic Relationships between Density Values in Cancellous and Cortical Bone. J. Biomech. 2008, 41, 1961–1968. [Google Scholar] [CrossRef] [PubMed]
- Abdullah, N.N.A.A.; Abdullah, A.H.; Ramlee, M.H. Current Trend of Lattice Structures Designed and Analysis for Porous Hip Implants: A Short Review. Mater. Today Proc. 2023. [Google Scholar] [CrossRef]
- Ataollahi, S. A Review on Additive Manufacturing of Lattice Structures in Tissue Engineering. Bioprinting 2023, e00304. [Google Scholar] [CrossRef]
- Rezapourian, M.; Hussainova, I. Optimal Mechanical Properties of Hydroxyapatite Gradient Voronoi Porous Scaffolds for Bone Applications—A Numerical Study. J. Mech. Behav. Biomed. Mater. 2023, 148, 106232. [Google Scholar] [CrossRef]
- Salaha, Z.F.M.; Ammarullah, M.I.; Abdullah, N.N.A.A.; Aziz, A.U.A.; Gan, H.-S.; Abdullah, A.H.; Abdul Kadir, M.R.; Ramlee, M.H. Biomechanical Effects of the Porous Structure of Gyroid and Voronoi Hip Implants: A Finite Element Analysis Using an Experimentally Validated Model. Materials (Basel). 2023, 16, 3298. [Google Scholar] [CrossRef]
- Vaiani, L.; Uva, A.E.; Boccaccio, A. Structural and Topological Design of Conformal Bilayered Scaffolds for Bone Tissue Engineering. Thin-Walled Struct. 2023, 192, 111209. [Google Scholar] [CrossRef]
- Atashipour, S.R.; Baqersad, J. Noninvasive Identification of Directionally-Dependent Elastic Properties of Soft Tissues Using Full-Field Optical Data. J. Mech. Behav. Biomed. Mater. 2024, 151, 106266. [Google Scholar] [CrossRef]
- Bachmann, S.; Pahr, D.H.; Synek, A. A Density-Dependent Target Stimulus for Inverse Bone (Re) Modeling with Homogenized Finite Element Models. Ann. Biomed. Eng. 2023, 51, 925–937. [Google Scholar] [CrossRef]
- Daynes, S. Isotropic Cellular Structure Design Strategies Based on Triply Periodic Minimal Surfaces. Addit. Manuf. 2024, 104010. [Google Scholar] [CrossRef]
- Hussain, Z.; Mehmood, S.; Liu, X.; Liu, Y.; Wang, G.; Pei, R. Decoding Bone-Inspired and Cell-Instructive Cues of Scaffolds for Bone Tissue Engineering. Eng. Regen. 2023. [Google Scholar] [CrossRef]
- Jiang, H.; Bednarcyk, B.A.; Le Barbenchon, L.; Chen, Y. Elastically Anisotropic Architected Metamaterials with Enhanced Energy Absorption. Thin-Walled Struct. 2023, 192, 111115. [Google Scholar] [CrossRef]
- Ritschl, L.; Schilling, P.; Wittmer, A.; Bohner, M.; Bernstein, A.; Schmal, H.; Seidenstuecker, M. Composite Material Consisting of Microporous Beta-TCP Ceramic and Alginate-Dialdehyde-Gelatin for Controlled Dual Release of Clindamycin and Bone Morphogenetic Protein 2. J. Mater. Sci. Mater. Med. 2023, 34, 39. [Google Scholar] [CrossRef] [PubMed]
- AboElhassan, R.G.; Watts, D.C.; Alamoush, R.A.; Elraggal, A. Biomechanical Behavior and Weibull Survival of CAD-CAM Endocrowns with Different Marginal Designs: A 3D Finite Element Analysis. Dent. Mater. 2024, 40, 227–235. [Google Scholar] [CrossRef]
- Kirmali, Ö.; Icen, G.; Celik, H.K.; Rennie, A.E.W. Evaluation of Stress Distribution on an Endodontically Treated Maxillary Central Tooth with Lesion Restored with Different Crown Materials: A Finite Element Analysis. Heliyon 2024. [Google Scholar] [CrossRef] [PubMed]
- Robau-Porrua, A.; González, J.E.; Rodríguez-Guerra, J.; González-Mederos, P.; Navarro, P.; de la Rosa, J.E.; Carbonell-González, M.; Araneda-Hernández, E.; Torres, Y. Biomechanical Behavior of a New Design of Dental Implant: Influence of the Porosity and Location in the Maxilla. J. Mater. Res. Technol. 2024, 29, 3255–3267. [Google Scholar] [CrossRef]
- Alemayehu, D.-B.; Jeng, Y.-R. Three-Dimensional Finite Element Investigation into Effects of Implant Thread Design and Loading Rate on Stress Distribution in Dental Implants and Anisotropic Bone. Materials (Basel). 2021, 14. [Google Scholar] [CrossRef]
- Pandoleon, P.; Bakopoulou, A.; Papadopoulou, L.; Koidis, P. Evaluation of the Biological Behaviour of Various Dental Implant Abutment Materials on Attachment and Viability of Human Gingival Fibroblasts. Dent. Mater. 2019, 35, 1053–1063. [Google Scholar] [CrossRef]
- Bulaqi, H.A.; Mashhadi, M.M.; Safari, H.; Samandari, M.M.; Geramipanah, F. Effect of Increased Crown Height on Stress Distribution in Short Dental Implant Components and Their Surrounding Bone: A Finite Element Analysis. J. Prosthet. Dent. 2015, 113, 548–557. [Google Scholar] [CrossRef]
- Kayabaşı, O.; Yüzbasıoğlu, E.; Erzincanlı, F. Static, Dynamic and Fatigue Behaviors of Dental Implant Using Finite Element Method. Adv. Eng. Softw. 2006, 37, 649–658. [Google Scholar] [CrossRef]
- Geramizadeh, M.; Katoozian, H.; Amid, R.; Kadkhodazadeh, M. Finite Element Analysis of Dental Implants with and without Microthreads under Static and Dynamic Loading. J. Long. Term. Eff. Med. Implants 2017, 27. [Google Scholar] [CrossRef] [PubMed]
- Chen, L.; Hao, H.E.; Li, Y.-M.; Ting, L.I.; Guo, X.-P.; Wang, R.-F. Finite Element Analysis of Stress at Implant–Bone Interface of Dental Implants with Different Structures. Trans. Nonferrous Met. Soc. China 2011, 21, 1602–1610. [Google Scholar] [CrossRef]
- Cheng, Y.-J.; Wu, T.-H.; Tseng, Y.-S.; Chen, W.-F. Development of Hybrid 3D Printing Approach for Fabrication of High-Strength Hydroxyapatite Bioscaffold Using FDM and DLP Techniques. Biofabrication 2024, 16, 25003. [Google Scholar] [CrossRef] [PubMed]
- He, S.; Zhu, J.; Jing, Y.; Long, S.; Tang, L.; Cheng, L.; Shi, Z. Effect of 3D-Printed Porous Titanium Alloy Pore Structure on Bone Regeneration: A Review. Coatings 2024, 14, 253. [Google Scholar] [CrossRef]
- Wang, Z.; Liao, B.; Liu, Y.; Liao, Y.; Zhou, Y.; Li, W. Influence of Structural Parameters of 3D-printed Triply Periodic Minimal Surface Gyroid Porous Scaffolds on Compression Performance, Cell Response, and Bone Regeneration. J. Biomed. Mater. Res. Part B Appl. Biomater. 2024, 112, e35337. [Google Scholar] [CrossRef] [PubMed]
- Ashby, M.F.; Gibson, L.J. Cellular Solids: Structure and Properties. Press Synd. Univ. Cambridge, Cambridge, UK 1997, 175–231.
- Distefano, F.; Pasta, S.; Epasto, G. Titanium Lattice Structures Produced via Additive Manufacturing for a Bone Scaffold: A Review. J. Funct. Biomater. 2023, 14, 125. [Google Scholar] [CrossRef]
- Suryawanshi, C.M.; Sagar, S.; Bhallamudi, R.; Mishra, S. Effective Design and Mechanical Response of Gyroid Lattice Scaffold for Orthopedic Implants. Manuf. Lett. 2023, 35, 493–501. [Google Scholar] [CrossRef]
- Verma, V.; Pal, K. A Finite Element Investigation on the Design of Mechanically Compatible Functionally Graded Orthopaedic Plate for Diaphyseal Tibia Transverse Fracture. Compos. Part C Open Access 2022, 7, 100228. [Google Scholar] [CrossRef]
- Zhai, Y.; Zhang, H.; Liu, T.; Zou, C.; Zhou, C. Mechanical Property of Ti6Al4V Cylindrical Porous Structure for Dental Implants Fabricated by Selective Laser Melting. Comput. Methods Biomech. Biomed. Engin. 2023, 1–19. [Google Scholar] [CrossRef] [PubMed]
- Fan, Z.; Huang, G.; Lu, Y.; Chen, Y.; Zeng, F.; Lin, J. Full Compression Response of FG-Based Scaffolds with Varying Porosity via an Effective Numerical Scheme. Int. J. Mech. Sci. 2022, 223, 107294. [Google Scholar] [CrossRef]
- Rezapourian, M.; Kamboj, N.; Jasiuk, I.; Hussainova, I. Journal of the Mechanical Behavior of Biomedical Materials Biomimetic Design of Implants for Long Bone Critical-Sized Defects. J. Mech. Behav. Biomed. Mater. 2022, 134, 105370. [Google Scholar] [CrossRef] [PubMed]
- Afshar, M.; Anaraki, A.P.; Montazerian, H.; Kadkhodapour, J. Additive Manufacturing and Mechanical Characterization of Graded Porosity Scaffolds Designed Based on Triply Periodic Minimal Surface Architectures. J. Mech. Behav. Biomed. Mater. 2016, 62, 481–494. [Google Scholar] [CrossRef] [PubMed]
- Alemayehu, D.B.; Jeng, Y.R. Three-Dimensional Finite Element Investigation into Effects of Implant Thread Design and Loading Rate on Stress Distribution in Dental Implants and Anisotropic Bone. Materials (Basel). 2021, 14. [Google Scholar] [CrossRef] [PubMed]
- Kladovasilakis, N.; Tsongas, K.; Tzetzis, D. Finite Element Analysis of Orthopedic Hip Implant with Functionally Graded Bioinspired Lattice Structures. Biomimetics 2020, 5. [Google Scholar] [CrossRef] [PubMed]
- Cantó-Navés, O.; Medina-Galvez, R.; Marimon, X.; Ferrer, M.; Figueras-Álvarez, Ó.; Cabratosa-Termes, J. A 3D Finite Element Analysis Model of Single Implant-Supported Prosthesis under Dynamic Impact Loading for Evaluation of Stress in the Crown, Abutment and Cortical Bone Using Different Rehabilitation Materials. Materials (Basel). 2021, 14, 3519. [Google Scholar] [CrossRef]
- Adachi, T.; Tsubota, K.; Tomita, Y.; Hollister, S.J. Trabecular Surface Remodeling Simulation for Cancellous Bone Using Microstructural Voxel Finite Element Models. J. Biomech. Eng. 2001, 123, 403–409. [Google Scholar] [CrossRef]
- Juneja, S.; Miranda, G.; Eram, A.; Shetty, N.; KN, C.; Keni, L.G. Investigating the Influence of All-Ceramic Prosthetic Materials on Implants and Their Effect on the Surrounding Bone: A Finite Element Analysis. Prosthesis 2024, 6, 74–88. [Google Scholar] [CrossRef]
- Kowalczyk, P. Elastic Properties of Cancellous Bone Derived from Finite Element Models of Parameterized Microstructure Cells. J. Biomech. 2003, 36, 961–972. [Google Scholar] [CrossRef] [PubMed]
- Shen, Y.-W.; Huang, H.-L.; Hsu, J.-T.; Fuh, L.-J. Effects of Diameters of Implant and Abutment Screw on Stress Distribution within Dental Implant and Alveolar Bone: A Three-Dimensional Finite Element Analysis. J. Dent. Sci. 2024. [Google Scholar] [CrossRef]
- Ulrich, D.; van Rietbergen, B.; Weinans, H.; Rüegsegger, P. Finite Element Analysis of Trabecular Bone Structure: A Comparison of Image-Based Meshing Techniques. J. Biomech. 1998, 31, 1187–1192. [Google Scholar] [CrossRef] [PubMed]
- Vafaeian, B.; El-Rich, M.; El-Bialy, T.; Adeeb, S. The Finite Element Method for Micro-Scale Modeling of Ultrasound Propagation in Cancellous Bone. Ultrasonics 2014, 54, 1663–1676. [Google Scholar] [CrossRef] [PubMed]
- Ye, Z.; Ye, H.; Wu, Y.; Jiang, Z.; Yao, H.; Xu, X.; Zhang, Y.; Du, W.; Li, W.; Zheng, Y. Effect of Bone Mass Density and Alveolar Bone Resorption on Stress in Implant Restoration of Free-End Edentulous Posterior Mandible: Finite Element Analysis of Double-Factor Sensitivity. Ann. Anatomy-Anatomischer Anzeiger 2024, 253, 152210. [Google Scholar] [CrossRef] [PubMed]
- Al-Sanea, A.; Aktas, S.; Celik, T.; Kisioglu, Y. Effects of the Internal Contact Surfaces of Dental Implants on Screw Loosening: A 3-Dimensional Finite Element Analysis. J. Prosthet. Dent. 2023, 130, 603-e1. [Google Scholar] [CrossRef]
- Falcinelli, C.; Valente, F.; Vasta, M.; Traini, T. Finite Element Analysis in Implant Dentistry: State of the Art and Future Directions. Dent. Mater. 2023. [Google Scholar] [CrossRef]
- Mishra, D.; Basu, B. Biomechanical Analysis of Peri-Prosthetic Bone Response to Hybrid Threaded Zirconia Dental Implants: An in Silico Model. J. Mech. Behav. Biomed. Mater. 2024, 150, 106310. [Google Scholar] [CrossRef]
- Młynarek-Żak, K.; Żmudzki, J. The Effect of Porous Compliance Bushings in a Dental Implant on the Distribution of Occlusal Loads. Sci. Rep. 2024, 14, 1607. [Google Scholar] [CrossRef]
- Qiu, P.; Cao, R.; Li, Z.; Fan, Z. A Comprehensive Biomechanical Evaluation of Length and Diameter of Dental Implants Using Finite Element Analyses: A Systematic Review. Heliyon 2024. [Google Scholar] [CrossRef]
- Wang, P.-S.; Tsai, M.-H.; Wu, Y.-L.; Chen, H.-S.; Lei, Y.-N.; Wu, A.Y.-J. Biomechanical Analysis of Titanium Dental Implants in the All-on-4 Treatment with Different Implant–Abutment Connections: A Three-Dimensional Finite Element Study. J. Funct. Biomater. 2023, 14, 515. [Google Scholar] [CrossRef]
- Bouakaz, I.; Dehkord, E.S.; Meille, S.; Schrijnemakers, A.; Boschini, F.; Preux, N.; Hocquet, S.; Geris, L.; Nolens, G.; Grossin, D. 3D Printed Triply Periodic Minimal Surfaces Calcium Phosphate Bone Substitute: The Effect of Porosity Design on Mechanical Properties. Ceram. Int. 2024, 50, 2623–2636. [Google Scholar] [CrossRef]
- Liu, F.; Chen, J.; Wang, Y.; Luo, T.; Zhang, X.; Wang, X.; Bian, X.; Xie, H. Multicomponent Functional Scaffolds with Multiscale Porosity to Improve Biomimetic and Mechanical Properties. Available SSRN 4725552.
- Karimi, M.; Asadi-Eydivand, M.; Abolfathi, N.; Chehrehsaz, Y.; Solati-Hashjin, M. The Effect of Pore Size and Layout on Mechanical and Biological Properties of 3D-printed Bone Scaffolds with Gradient Porosity. Polym. Compos. 2023, 44, 1343–1359. [Google Scholar] [CrossRef]
- Mao, Y.; Yuan, J.; Heng, Y.; Feng, K.; Cai, D.; Wei, Q. Effect of Hot Isostatic Pressing Treatment on Porosity Reduction and Mechanical Properties Enhancement of 316L Stainless Steel Fabricated by Binder Jetting. Virtual Phys. Prototyp. 2023, 18, e2174703. [Google Scholar] [CrossRef]
- Zheng, G.; Wei, K.; Kang, X.; Fan, W.; Ma, N.L.; Verma, M.; Ng, H.S.; Ge, S. A New Attempt to Control Volatile Organic Compounds (VOCs) Pollution-Modification Technology of Biomass for Adsorption of VOCs Gas. Environ. Pollut. 2023, 122451. [Google Scholar] [CrossRef] [PubMed]
- Araya, M.; Jaskari, M.; Rautio, T.; Guillén, T.; Järvenpää, A. Assessing the Compressive and Tensile Properties of TPMS-Gyroid and Stochastic Ti64 Lattice Structures: A Study on Laser Powder Bed Fusion Manufacturing for Biomedical Implants. J. Sci. Adv. Mater. Devices 2024, 9, 100663. [Google Scholar] [CrossRef]
- Seehanam, S.; Khrueaduangkham, S.; Sinthuvanich, C.; Sae-Ueng, U.; Srimaneepong, V.; Promoppatum, P. Evaluating the Effect of Pore Size for 3d-Printed Bone Scaffolds. Heliyon 2024. [Google Scholar] [CrossRef]
- Wang, C.; Liu, J.; Min, S.; Liu, Y.; Liu, B.; Hu, Y.; Wang, Z.; Mao, F.; Wang, C.; Ma, X. The Effect of Pore Size on the Mechanical Properties, Biodegradation and Osteogenic Effects of Additively Manufactured Magnesium Scaffolds after High Temperature Oxidation: An in Vitro and in Vivo Study. Bioact. Mater. 2023, 28, 537–548. [Google Scholar] [CrossRef] [PubMed]
- Song, J.; Lee, J.; Kim, N.; Min, K. Artificial Intelligence in the Design of Innovative Metamaterials: A Comprehensive Review. Int. J. Precis. Eng. Manuf. 2024, 25, 225–244. [Google Scholar] [CrossRef]
- Carter, D.R.; Hayes, W.C. The Compressiv Berhaviour of Bone as a Two-Phase Porous Material. J. Bone Joint. Surg. 1977, 59, 954–962. [Google Scholar] [CrossRef]
- Schaffler, M.B.; Burr, D.B. Stiffness of Compact Bone: Effects of Porosity and Density. J. Biomech. 1988, 21, 13–16. [Google Scholar] [CrossRef]
- Gibson, L.J. The Mechanical Behaviour of Cancellous Bone. J. Biomech. 1985, 18, 317–328. [Google Scholar] [CrossRef] [PubMed]
- Yagihara, A.; Kawasaki, R.; Mita, A.; Takakuda, K. Impact of Dynamic and Static Load on Bone Around Implants: An Experimental Study in a Rat Model. Int. J. Oral Maxillofac. Implants 2016, 31. [Google Scholar] [CrossRef] [PubMed]
- Po, J.M.C.; Kieser, J.A.; Gallo, L.M.; Tésenyi, A.J.; Herbison, P.; Farella, M. Time-Frequency Analysis of Chewing Activity in the Natural Environment. J. Dent. Res. 2011, 90, 1206–1210. [Google Scholar] [CrossRef] [PubMed]
Figure 1.
2D drawing of dental implant assembly in creo parametric: front, side, and top views with detailed dimensions (mm).
Figure 1.
2D drawing of dental implant assembly in creo parametric: front, side, and top views with detailed dimensions (mm).
Figure 2.
Comparative analysis of biomimetic Voronoi latticed spongy bone structures with four distinct pore sizes, a)1.0 mm, b) 1.5 mm, c) 2.0 mm, and d) 2.5 mm, designed using nTopology software.
Figure 2.
Comparative analysis of biomimetic Voronoi latticed spongy bone structures with four distinct pore sizes, a)1.0 mm, b) 1.5 mm, c) 2.0 mm, and d) 2.5 mm, designed using nTopology software.
Figure 3.
Biomimetic Voronoi latticed trabecular bone constructed in nTopology featuring: a) open-cell randomized lattice graph, b) truncated thicken lattice graph, c) finite element (FE) volume mesh, d) FE boundary configuration utilizing body nodes.
Figure 3.
Biomimetic Voronoi latticed trabecular bone constructed in nTopology featuring: a) open-cell randomized lattice graph, b) truncated thicken lattice graph, c) finite element (FE) volume mesh, d) FE boundary configuration utilizing body nodes.
Figure 4.
Components of the dental implant system including: a) application of dynamic explicit oblique loading at load reference point with MPC constraint, b) 3d model components designed for finite element analysis.
Figure 4.
Components of the dental implant system including: a) application of dynamic explicit oblique loading at load reference point with MPC constraint, b) 3d model components designed for finite element analysis.
Figure 6.
Relative Density and Porosity of Biomimetic Voronoi Lattice Trabecular Bone Across Various Pore Sizes.
Figure 6.
Relative Density and Porosity of Biomimetic Voronoi Lattice Trabecular Bone Across Various Pore Sizes.
Figure 7.
An equivalent oblique dynamic load distribution on dental implant crown during mastication, a multidirectional analysis over 0.5s cycle.
Figure 7.
An equivalent oblique dynamic load distribution on dental implant crown during mastication, a multidirectional analysis over 0.5s cycle.
Figure 8.
Maximum von mises stress distribution in double-sliced assembled implant parts for four biomimetic Voronoi latticed trabecular bones: a) VTB10, b) VTB15, c) VTB20, and d) VTB25.
Figure 8.
Maximum von mises stress distribution in double-sliced assembled implant parts for four biomimetic Voronoi latticed trabecular bones: a) VTB10, b) VTB15, c) VTB20, and d) VTB25.
Figure 9.
Contour von Mises stress distribution in dental implants and retaining screws across biomimetic Voronoi latticed trabecular bone pore sizes: a) VTB10, b) VTB15, c) VTB20, d) VTB25 for implants, and e) VTB10, f) VTB15, g) VTB20, h) VTB25 for retaining screws.
Figure 9.
Contour von Mises stress distribution in dental implants and retaining screws across biomimetic Voronoi latticed trabecular bone pore sizes: a) VTB10, b) VTB15, c) VTB20, d) VTB25 for implants, and e) VTB10, f) VTB15, g) VTB20, h) VTB25 for retaining screws.
Figure 10.
Maximum von Mises stresses in biomimetic Voronoi lattice trabecular and cortical bones for a) VTB10, b) VTB15, c) VTB20, and d) VTB25.
Figure 10.
Maximum von Mises stresses in biomimetic Voronoi lattice trabecular and cortical bones for a) VTB10, b) VTB15, c) VTB20, and d) VTB25.
Figure 11.
Contour plots of magnitude of displacements in biomimetic Voronoi lattice trabecular bone for pore sizes: a) VTB10, b) VTB15, c) VTB20, d) VTB25.
Figure 11.
Contour plots of magnitude of displacements in biomimetic Voronoi lattice trabecular bone for pore sizes: a) VTB10, b) VTB15, c) VTB20, d) VTB25.
Figure 12.
Maximum displacement response in axial, mesiodistal, and buccolingual directions for biomimetic Voronoi lattice trabecular bones: a) VTB10, b) VTB15, c) VTB20, and d) VTB25.
Figure 12.
Maximum displacement response in axial, mesiodistal, and buccolingual directions for biomimetic Voronoi lattice trabecular bones: a) VTB10, b) VTB15, c) VTB20, and d) VTB25.
Figure 13.
Dynamic reaction forces in dental implant assemblies across biomimetic Voronoi latticed trabecular bone pore sizes: a) VTB10, b) VTB15, c) VTB20, and d) VTB25.
Figure 13.
Dynamic reaction forces in dental implant assemblies across biomimetic Voronoi latticed trabecular bone pore sizes: a) VTB10, b) VTB15, c) VTB20, and d) VTB25.
Table 3.
Effects of pore size variation on structural parameters of biomimetic Voronoi lattice trabecular bone.
Table 3.
Effects of pore size variation on structural parameters of biomimetic Voronoi lattice trabecular bone.
Symbolic name |
Pore size (mm) |
Beam thickness (mm) |
Randomization Seed |
Average Beam Length (mm) |
Lattice beam account |
Lattice Node count |
Relative Density (%) |
Surface area to volume ratio (%) |
VTB10 |
1.0 |
0.4 |
250 |
0.4614 |
25,809 |
15,186 |
22.21 |
10.49 |
VTB15 |
1.5 |
0.4 |
250 |
0.5681 |
13,286 |
8,182 |
15.19 |
11.58 |
VTB20 |
2.0 |
0.4 |
250 |
0.7532 |
5,438 |
3,535 |
8.84 |
12.21 |
VTB25 |
2.5 |
0.4 |
250 |
0.9203 |
3,024 |
2,058 |
6.27 |
12.4 |
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).