A Modeling Method of Graded Porous Scaffold Based on Triply Periodic Minimal Surfaces

In order to solve the design and regulation problems of graded porous scaolds, the triply periodic minimal surface (TPMS) that based on implicit function was taken as the pore unit for constructing the porous scaolds. And the TPMS structure was controlled by selecting dierent types of TPMS and dening the distance function. e transition of dierent structures of TPMS was performed by the weight fusion method. And the dual contouring (DC) algorithm and octree algorithm were used to divide the solid model of the defect bone into hexahedral elements to determine the target assembly spatial. rough the coordinate transformation of the isoparametric element method, the pore units were mapped into the hexahedral elements to construct porous structure. Finally, a Boolean operation was used to obtain a bone scaold model with porous structure. e examples showed that the gradient model of various irregular pores could be constructed, and the pore size, porosity, and specic surface area of porous scaolds could be controlled by adjusting the deformation, distribution, volume of the hexahedral elements, and the structure of the TPMS pore units, which provides a feasible method for the design and regulation of heterogeneous porous scaolds.


Introduction
Bone tissue engineering (BTE) is a biological substitute that aims to restore, maintain, or improve tissue functions, which has a good clinical application prospect [1,2]. Bionic bone sca old has played an important role in bone tissue engineering, which functions as the architecture for bone growth.
e design of bone sca olds requires that the macrostructure should be consistent with the surface contour of patients' autologous bone, and the microstructure requires that the sca olds have micro-porous structures that provide living space and channels for the adhesion, reproduction, nutrition acquisition, and metabolism of seed cells.
Bone sca old can be manufactured using one of many available conventional techniques, such as particulate leaching, phase method, electrostatic spinning method, and gas method. However, these techniques have limitations with incapable of precisely controlling pore size and shape, porosity, and their spatial distributions within the sca old [3,4].
To overcome the drawbacks, additive manufacturing (AM) techniques were used to prepare the sca old properties [5][6][7]. e sca old was built by the AM techniques with the three-dimensional digital model. However, due to the irregularity of porosity, pore size, and pore distribution of the micro-porous structure of natural bone tissue, it is di cult to e ectively describe them by mathematical methods. At present, the digital model design of porous structure was mainly based on computer-aided design (CAD) method, which used the simple sphere, cylinder, prism, and polyhedron to design porous structure [8][9][10][11][12][13]. However, due to the limitation of the pore shapes, the designed porous structure lacks diversity, especially when designing gradient porous structures, the boundary fusion problem frequently occurred at the contact surface of the adjacent unit, and it is di cult to construct porous structures with gradient variations in pore shape, porosity, pore size, and spatial distributions.
Butterfly wings, weevils, crustaceans, etc., in nature have the skeletal structures similar to triply periodic minimal surface (TPMS) [14]. TPMS has advantages including various geometric shapes and can construct parametric mathematical models, which cause researcher's attention. Rajagopalan and Robb [15] proposed a tissue scaffold design method firstly, constructed a porous scaffold model based on P surface, and prepared them by additive manufacturing technology. Abou-Ali et al. [16] adopted D surface, I-WP surface, G surface, and CY surface as the basic pore unit and prepared nylon thermoplastic scaffolds by selective laser sintering technique, and the mechanical properties of the porous scaffolds with different types of TPMS were experimentally studied. Montazerian et al. [17]and Castro et al. [18] constructed the porous scaffold model based on the P surface, I-WP surface, and G surface, the permeability of the porous scaffold was analyzed via computational fluid dynamics simulation, and the experimental results showed that the permeability of porous scaffold was closely related to the types of TPMS. e natural bone is heterogeneous, which needs to meet the specific requirements of biomechanics and biocompatibility at a specific position. However, a homogeneous porous structure is difficult to become an ideal scaffold due to its single performance. Zhang et al. [19] found that the gradient porous scaffolds were superior to the homogeneous porous scaffolds in comprehensive performance. Han et al. [20] reported that the gradient porous scaffold could optimize the growth space of bone tissue through personalized design because of its varied pore size and porosity. erefore, the research on modeling methods of heterogeneous porous structure with functional gradient has always been the focus of this field. In recent years, some scholars have also explored the modeling method of gradient porous scaffolds based on TPMS [21,22]. Liu et al. [23] constructed a gradient porous scaffold model based on D surface and G surface and fabricated by selective laser melting (SLM), and the results showed that the TPMS method was an effective way to achieve gradients in multiple patterns which were comparable to natural tissue with respect to both continuous topology and interconnectivity. Gao and Xiang [24] proposed a design method for gradient porous scaffolds based on G surface. e mechanical properties of the gradient porous scaffold were analyzed, and the results showed that the mechanical properties of the gradient porous scaffolds were better than the homogeneous porous scaffolds under the same average porosity. Although some achievements have been made, there are still many problems in the construction of gradient porous scaffolds, such as the arbitrary arrangement of different pore structures and the precise control of porosity, pore diameter, and specific surface area. e research of gradient structure design method based on TPMS is still in its initial stage. In this paper, TPMS was used as the basic pore unit to construct the gradient porous scaffolds by controlling the pore unit structure and its spatial morphology. e effects of pore unit structure and its assembly spatial configuration on porosity, pore size, specific surface area, and their spatial distributions were studied, which provides a feasible method for the modeling of gradient porous scaffolds.

Algorithm Flow
e modeling method and general route of the porous scaffold are shown in Figure 1.
(1) CT/MRI was used to measure the defect bone, and the solid model of the defect bone was constructed. (2) TPMS expressed by implicit function was used as the basic pore unit for constructing porous structure. And the pore unit structure was controlled by selecting different types of TPMS and defining the distance function. e transition of different structures of TPMS pore units was performed by the weight fusion method. (3) e DC algorithm and octree algorithm were used to divide the solid model of the defect bone into hexahedral elements to determine the target assembly spatial. And the hexahedral elements could be regulated by encryption template to meet the design requirements. (4) rough the coordinate transformation of the shape function, the pore units were mapped into the space units to construct porous structure. (5) A Boolean operation was used to obtain bone scaffold model with porous structure.

Pore Unit Based on TPMS
3.1. TPMS Unit. TPMS is a minimal surface with periodicity in three independent directions in three-dimensional space (the average curvature of any point on the surface is zero), which is usually expressed by implicit surface [25].
where A k is the amplitude factor, λ k is the periodic wavelength, r is the position vector of Euclidean space, h k is the kth raster vector in the reciprocal space, and P k is the phase deviation. e typical TPMS [26]is shown in Table 1.where (1) is an implicit expression of TPMS. And marching cubes (MC) method [27] was adopted to construct TPMS, as shown in Figure 2.
For the same type of TPMS, the TPMS structure could also be changed by controlling the distance k from the point of φ TPMS (x, y, z, k) to the φ TPMS (x, y, z) [28] as shown in Figure 3.

Fusion Transition of Different Architecture TPMS Units.
e transition of different structures of TPMS (including the same type of TPMS unit under different k values) was performed by the weight fusion method [29].

TPMS
Implicit surface expression P surface where n is the type of TPMS, φ j (x, y, z) is the type of TPMS unit in region j, and α j (x, y, z) is the weight coefficient for region j.
where α ij (x, y, z) is the weight coefficient of the i-th point in region j, i � 1,2, . . ., m, which could be constructed by the exponential kernel function.
where (x ij , y ij , z ij ) is the coordinate of the i-th point in region j and σ is the width parameter, which controls the local scope of the exponential kernel function. Take P surface, FRD surface, K surface, and G surface in Table 1 as an example for multi-region fusion transition. Figure 4 shows the fusion transition results. Figure 5 shows the fusion transition results of the same type of TPMS (P surface) under different k.

Construction of TPMS Point Unit.
In this paper, the algebraic method was used to construct TPMS point unit by intersecting the TPMS with the plane perpendicular to the coordinate axis.
e FRD, P, and G point units constructed are shown in Figure 6.

Model Space Subdivision
e purpose of model space subdivision was to provide assembly spatial for the designed TPMS pore unit. In this paper, the method proposed by Zhang et al. [30] was adopted to conduct hexahedral mesh division for the solid model of the defect bone and the irregular hexahedral elements were obtained. e algorithm was shown as follows.
Step 1. Octree algorithm was used to establish the topological relationship of the internal mesh.
Step 2. e boundary of the model was extracted using dual contouring (DC) algorithm to construct the boundary mesh.
Step 3. Relaxation technology was used to improve mesh quality.
Taking the triangular prism as an example, the mesh subdivision results are shown in Figure 7.
In order to control the size and distribution of hexahedral elements, the encryption template proposed by Zhang [31] was used to encrypt hexahedral mesh. e triangular prism mesh encryption results are shown in Figure 8.

Pore Unit Mapping and Internal Porous
Structure Evaluation 5.1. Pore Unit Mapping. In order to construct irregular pore unit, the constructed TPMS pore unit was mapped into the hexahedral elements based on the isoparametric element method [32].
In this paper, eight-node hexahedral element was used. e mapping relationship of the isoparametric element method was where x i , y i , and z i are the coordinate points of hexahedral element nodes in the Cartesian coordinates. x, y, and z are Cartesian coordinates of points inside hexahedral elements. N i is node-shaped function, which is an interpolation function expressed by natural coordinates. It could be obtained by Lagrange interpolation method: ξ, η, and ζ are the internal points of the entity in the natural coordinate system, − 1 < ξ < 1, − 1 < η < 1, and − 1 < ζ < 1, respectively. ξ i , η i , and ζ i are the coordinates of the i-th node in the natural coordinate system. Taking P and FRD pore unit as example, the irregular pore unit obtained after mapping is shown in Figure 9.
Namely, the porous bionic bone scaffold model could be constructed based on the Boolean operation.
where φ i is the pore unit of i-th and φ3 D scaffold is the porous bone scaffold with macro-and microstructure.

Evaluation of Internal Porous
Structure. e internal porous structure of the bone scaffold plays the important role in the osteogenesis stage of seed cell. erefore, after constructing the tissue engineering bone scaffold model, it is necessary to evaluate its porosity, pore size, and specific surface area.

Porosity.
Porosity is an important parameter to evaluate the pore structure of bone scaffolds.
. e TPMS structure could be changed by controlling the k from the point of φ TPMS (x, y, z, k) to the minimum surface defined by φ TPMS (x, y, z).  Table 1 as an example for multi-region fusion transition. where η is the porosity, V HEX is the volume of the hexahedral element, and V P is the volume of the pore unit. In this paper, the pore unit constructed by TPMS was a surface model that is composed of triangular facets, so the center point of the hexahedral element and triangular facet could form a tetrahedron, and the volume of the tetrahedron V k tetrahedron could be calculated. e volume of the pore unit V P is obtained by accumulating all tetrahedral volumes.
where (x i , y i , z i ) are the coordinates of the hexahedron nodes of the i-th hexahedral element, and N ai k , N bi k , and N ci k are the corresponding shape functions of the three vertices of the k-th triangular facet, respectively, which can be obtained by (10). V k tetrahedron is the volume of k-th tetrahedron.
For an arbitrarily irregular hexahedron element, its volume could be calculated by dividing it into five tetrahedrons. where are the coordinates of the eight nodes of the i-th hexahedral element. V i is the volume of tetrahedron.   e pore size is usually defined as the maximum inner sphere diameter of the pore. erefore, for the pore size calculation, the voxel method [33] was adopted to divide the pore unit equidistant, as shown in Figure 10. e maximum inscribed circle diameter d of the pore unit was approximately solved by the Euclidean distance [34].
where r is the distance from the center of the i-th voxel marked as "1" to the center of the nearest voxel marked as "0." (x vi , y vi , z vi ) are the center coordinates of the i-th voxel marked as "1," and (x vi , y vi , z vi ) are the center coordinates of the j-th voxel marked as "0."

Specific Surface
Area. e specific surface area is a description of the pore space, the adsorption capacity would be stronger, and the space for cell adhesion, migration, and proliferation would be wider with the specific surface area larger.
where S i is the area of the triangle in surface model of the TPMS unit, which could be calculated by Helen's formula.

Example Verification
In order to illustrate the effectiveness of the modeling method, a bone scaffold modeling program was developed in Top view. For the pore size calculation, the voxel method [30] was adopted to divide the pore unit equidistant.
Visual Studio 2008 development environment, combined with OpenGL technology. And the cranial prosthesis was used as an example for analysis and discussion, and the modeling process is shown in Figure 11. Firstly, the surface model of bone defect was extracted from CT or MRT DICOM data by Mimics software. Secondly, the solid model of defect bone was reconstructed by CATIA software. Subsequently, the reconstructed model was subdivided into hexahedrons, and 84 hexahedral elements were generated in this case. en, the TPMS unit was mapped into hexahedral elements and irregular pore unit could be obtained. Finally, the porous bionic bone scaffold model could be constructed with the Boolean operation.
6.1. Design of Single Pore Structure. In this paper, the P and FRD surfaces in Table 1 were taken as an example for modeling of porous scaffold. e distance function k is selected as +0.9, +0.7, +0.5, +0.3, +0.1, 0, − 0.1, − 0.3, − 0.5, − 0.7, and − 0.9, respectively. For the representation of the spatial configuration of the hexahedral elements, the deformation rate and volume were selected to describe it, where the deformation rate was calculated by where l min is the shortest edge of the hexahedral element and l max is the longest edge of the hexahedral element. e distribution of the volume and deformation rate of hexahedral elements is shown in Figure 12.
is paper took the P and FRD unit under k � +0.9, k � 0, and k � − 0.9 as an example to introduce, and the constructed porous scaffold model is shown in Figure 13.

Porosity Distribution.
e porosity distribution of P and FRD unit of porous scaffold under k � +0.9, 0, − 0.9 is shown in Figure 14.
It can be seen from Figure 14 that the porosity of each pore unit of the porous scaffolds constructed by P and FRD unit showed different distributions. However, for the same type of TPMS, the porosity of each pore unit was a certain value under the same k values. Under k � − 0.9, 0, and +0.9, the porosity of P unit was 79.44%, 56.46%, and 30.95%, and the FRD unit was 33.87%, 55.68%, and 74.84%. Above examples show that the porosity is related to the TPMS pore unit structure (TPMS type and k value), while the spatial configuration of hexahedral elements has no effect on the porosity.
It could be seen from Figure 15 that, for the porous scaffolds constructed by P unit, the porosity decreases with the increase of k value, while for the FRD unit, the porosity increases with the increase of k value, and the change of porosity tends to be linear. e fitted linear equation is as follows: where η P is the porosity of the pore unit constructed by P unit and η FRD is the porosity of the pore unit constructed by FRD unit.
It could be seen from Figure 16 that the pore size of porous scaffold constructed by P and FRD unit showed different distribution rules with a different TPMS type and k value. However, compared with the pore size distribution of the same type of pore unit under different k values, it could be seen that although the pore size distribution varies with the k value, the change trend of the pore size distribution was  Mathematical Problems in Engineering similar, because the final shape of the pore unit was also controlled by the volume, deformation rate, and spatial distribution of the mapping hexahedral elements. Figure 17 shows the specific surface area distribution of P and FRD unit that under k � +0.9, k � 0, and k � − 0.9. It could be seen from Figure 17 that, similar to the porosity distribution, the specific surface area of each pore unit also presents different distribution trends with the different pore unit structures (TPMS type and k values). In order to explain the change rule of specific surface area, several pore units were selected for the explanation, as shown in Table 2.

Specific Surface Area.
It could be seen from Table 2 that at 39 and 81 the pore unit structure and the volume of hexahedral elements were similar (6.663E10 μm 3 and 6.658E10 μm 3 ), but due to the deformation rate of hexahedral elements (0.62 and 0.48), the specific surface areas of pore unit were different (P unit was 5.688E − 4 μm − 1 and 6.032E − 4 μm − 1 ; FRD unit was 3.205E − 3 μm − 1 and 3.411E − 3 μm − 1 ), and the specific surface area of the pore unit increases with the decrease of the deformation rate of hexahedral elements, while at 31 and 19 hexahedral elements, the volume deformation rate of hexahedral elements was the same (0.59), but due to the volume of hexahedral elements (5.801E10 μm 3 and 1.035E11 μm 3 ), the specific surface area of pore unit was also different (P unit was 6.021E − 4 μm − 1 and 4.885E − 4 μm − 1 ; FRD unit was 3.383E − 3 μm − 1 and 2.760E − 3 μm − 1 ), and the specific surface area of the pore unit decreases with the increase of the volume of hexahedral elements. However, at 53 and 64 hexahedral elements, due to the similarity of the hexahedral mesh volumes and deformation rates, the specific surface areas of the pore unit were closer (P unit was 5.969E − 4 μm − 1 and 5.978E − 4 μm − 1 ; FRD unit was 3.352E − 3 μm − 1 and 3.370E − 3 μm − 1 ).
In conclusion, the porous scaffold constructed by TPMS unit could not only construct the internal pores with various shapes, but also adjust pore size and shape, porosity, specific surface area, and their spatial distributions through TPMS type, k value, and spatial configuration of hexahedral elements.

Gradient Pore Structure Design Based on Hexahedral
Mesh Encryption. In this paper, the gradient structure was designed by changing the spatial configuration of the hexahedral elements. e initial hexahedral mesh was encrypted using fully encrypted template. e hexahedral mesh after encryption is shown in Figure 18.
ere are 648 hexahedral elements in the encrypted region, 312 hexahedral elements in the transition region, and 36 hexahedral elements in the sparse region. e volume distribution is shown in Figure 19(a), and the deformation rate distribution is shown in Figure 19(b).
e design results are shown in Figure 20 with the P unit as the basic pore unit for constructing micro-pores, and k value is 0.
It could be seen from Figure 20 that the pores of the scaffold are gradient distribution, and the pores in different regions are excessively smooth.

Porosity Distribution.
e porosity distribution is shown in Figure 21.
It could be seen from Figure 21 that the porosity of the hexahedral elements in transition region, encryption region, and sparse region is all 56.46%. is was because the porosity was determined by the structure of the pore unit, and changing the spatial configuration of the hexahedral elements had no effect on the porosity.

Pore Size and Distribution.
e pore size and distribution are shown in Figure 22.
It could be seen from Figure 22 that the minimum pore size was reduced by nearly half compared with that before hexahedral mesh encryption, and the distribution of pore size also changed. is was because the final shape of the pore unit was controlled not only by TPMS type and k value but also by the spatial configuration of the hexahedral elements. e volume, deformation rate, and spatial distribution of hexahedral elements were changed due to the hexahedral mesh encryption, thus indirectly changing the pore size and distribution of the porous scaffold.
From Figure 23, it could be seen that the specific surface area distribution is more concentrated in the encrypted region and sparse region, while in the transition region the specific surface area distribution is more scattered.
is is due to the similar spatial configuration of the hexahedral elements in encrypted and sparse region, where the mean volume of the hexahedral elements in encrypted and sparse region was 3.764E9 μm 3 and 3.290E10 μm 3 , and the standard deviation 9.386E8 μm 3 and 5.896E9 μm 3 . e average deformation rate was 0.776 and 0.576, and the standard deviation was 0.067 and 0.063. However, in transition region, the mean volume of the hexahedral elements was 4.828E9 μm 3 , the standard deviation was 2.604E9 μm 3 , the mean deformation rate was 0.389, and the standard deviation was 0.131. e spatial configuration of the hexahedral elements in transition region was more variable compared with the encryption and the sparse region, so the specific surface area distribution of the pore unit in this region was more dispersed, mainly concentrated in 1.97E − 3 μm − 1 ∼ 7.12E − 3 μm − 1 , and in encrypted and FRD, k=-0.9 P, k=+0.9 P, k=0 P, k=-0.9 Figure 14: Porosity distribution. It could be seen from Figure 14 that the porosity of each pore unit of the porous scaffolds constructed by P and FRD unit showed different distributions. However, for the same type of TPMS, the porosity of each pore unit was a certain value under the same k values.     sparse region was mainly concentrated in 2.
In conclusion, the porous scaffold constructed by TPMS unit, without changing the TPMS unit type and k value, realized the gradient pore structure design only by mesh encryption.

Gradient Pore Structure Design Based on Different k
Values. In this section, gradient pore structure design was realized only by changing the distance function k value. e P unit was still used as the basic pore unit for constructing microscopic pores. k 1 of P unit was set as +0.9 in region A, k 2 of P unit was set as 0 in region B, and region C was the transition region. e design results are shown in Figure 24.
It could be seen from Figure 24 that the pore structure of designed scaffold was smooth, and the pore structure present was in a gradient distribution. e distribution of the porosity, pore size, and specific surface area are shown in Figure 25.
It could be seen from the porosity distribution in Figure 25(a) that the porosity distribution varies in different regions because the porosity was controlled by the pore unit structure. For the gradient pore structure scaffold constructed by the same type of TPMS, the pore structure in local areas could be controlled by k value. In this case, the porosity in region A and region B was 30.95% and 56.46%, and the transition region C was 43.90%.
For the pore size and distribution, it could be seen from Figure 25(b) that the pore size and pore size distribution were different compared to the porous scaffold constructed from a single k, due to the different pore structures within each region.
Similarly, the specific surface area distribution is shown in Figure 25(c), the specific surface area showed different distributions in different regions, and due to the similar spatial configuration of the hexahedral elements, the specific surface area distribution was concentrated in each region. In region A, the specific surface area of each pore unit concentrated 1.38E − 3 μm − 1 ∼ 1.82E − 3 μm − 1 , in region B the specific surface area of each pore unit concentrated 9.8E − 4 μm − 1 ∼ 1.34E − 3 μm − 1 , and in the transition region C the specific surface area of each pore unit concentrated erefore, for the same type of TPMS unit, when constructing gradient pore structure, the porosity, pore size, pore distribution, and specific surface area of bone scaffold could be effectively controlled by adjusting the k value in the local areas.

Gradient Pore Structure Design Based on Different TPMS
Types. In this case, the gradient pore structure design was realized by TPMS type without changing the spatial shape of the hexahedron mesh element and k value. FRD and P unit 6 Figure 21: Porosity distribution. e porosity of the hexahedral elements in transition region, encryption region, and sparse region is all 56.46%. e minimum pore size was reduced by nearly half compared with that before hexahedral mesh encryption, and the distribution of pore size also changed.  were used as basic pore unit. FRD unit was filled in region A, and k 1 is 0. P unit was filled in region B, and the k 2 is 0. e design result is shown in Figure 26. It could be seen from Figure 26 that the pore structure of designed scaffold was smooth, and the pore structure present was in a gradient distribution. e distribution of the porosity, pore size, and specific surface area is shown in Figure 27. e distribution of porosity is shown in Figure 27(a). Due to the different types of pore unit in region A and B, the pore structure and porosity were also different. e porosity in region A and B    the porosity in transition region C was 67.62% which was higher than the other two regions. As for the pore size and distribution, it could be seen from Figure 27(b) that the pore size and pore size distribution were different compared to the porous scaffold constructed from a single type of TPMS, due to the different pore structures within each region.
Similarly, the specific surface area distribution is shown in Figure 27(c), the specific surface area showed different distributions in different regions, and due to the similar spatial configuration of the hexahedral elements, the specific surface area distribution was concentrated in each region. In region A, the specific surface area of each pore unit concentrated on the range of 1.65E − 3 μm − 1 ∼ 2.16E − 3 μm − 1 , in region B the specific surface area of each pore unit concentrated in the range of 9.8E − 4 μm − 1 ∼ 1.34E − 3 μm − 1 , and in the transition region C the specific surface area of each pore unit concentrated in the range of 1.38E − 3 μm − 1 ∼ 1.87E − 3 μm − 1 .
It could be seen from the above experiments that the transition of different types of TPMS unit could not only enrich the geometry of pore unit in the scaffold, but also effectively control the porosity, pore size, and distribution of bone scaffolds. Single pore structure (P, k=0) Single pore structure (FRD, k=0) Gradient pore structure (P, k=+0) and P (FRD, k=0)

SLM Preparation of Porous Scaffolds.
e Ti6Al4V titanium alloy powder produced by EOS company was used as raw material, the porous support designed in example 4 was prepared by dimetal-100 h laser selective melting molding equipment, and the preparation process was carried out under the protection of argon gas. e parameters of the preparation were as follows: scanning speed was 1.2 m/s, the laser power was 280 W, the layer thickness was 35 μm, the prepared porous gradient scaffold sample is shown in Figure 28 It could be seen from Figure 28 that the molding quality of the parts was well, the transition between different types of TPMS unit was relatively smooth, and the prepared porous scaffold was consistent with the design model, which verifies the feasibility of the method proposed in this paper.

Conclusion
(1) Taking TPMS as the basic pore unit, the gradient porous scaffold modeling can be constructed by controlling the pore element structure and its spatial shape. (2) Several examples show that the pore size, porosity, and specific surface area of porous scaffolds could be controlled by adjusting the deformation, distribution, volume of the hexahedral elements, and the structure of the TPMS pore units. (3) In the next work, the pore types of TPMS unit will be further increased, and the structural model which is functionally equivalent to the real bone structure in terms of pore structure and biological properties will be studied.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.