# Elastic Properties Investigation on Random and Ordered ZrO2 Nanotube-Reinforced HA and β-TCP Biocomposites with Finite Element Approach

Authors: Zeliang Liu; Huijian Li; Zhixin Xu; Yuying Cao; Xiao Yang
Journal: Journal of Nanomaterials
 (2020)

---

## Abstract

The elastic properties of random and ordered ZrO2 nanotube- (ZrNT-) reinforced HA and β-TCP biocomposites were carried out by a numerical investigation with finite element approach. The elastic modulus, shear modulus, and Poisson’s ratio affected by various ZrNT volume fractions (5.0 vol.%, 5.5 vol.%, 6.5 vol.%, 7.5 vol.%, 8.5 vol.%, and 10.0 vol.%) and aspect ratios (3, 5, 10 and 20) for both random and ordered reinforced composites were obtained and analysed. The advantages of random and ordered reinforced composites were further discussed. The random reinforced composite is suggested to be a proper candidate for ZrO2 nanotube-reinforced biocomposite on the application of bone repair and the substitute.


---

## Body

## 1. Introduction

Hydroxyapatite (HA) and tricalcium phosphate (β-TCP) are widely used in dental implants, bone repair scaffolds, and bone implants because their crystal structure is similar to that of apatite crystals in living bone tissue. Their porous structural features further induce outstanding bioactivity and biocompatibility [1–4]. Unfortunately, because of their poor strength and brittle properties, they are limited in usage under high load-bearing conditions [5, 6].To increase the mechanical properties of HA andβ-TCP for clinical applications, one can add a second phase such as metal oxide powders including ZrO2-Al2O3, TiO2, and ZrO2 [7–9] or nanotubes such as CNT and BNNT [10, 11]. In biocomposite studies, the biocompatibility and bacterial adhesion of material should be firstly considered. It has been reported that ZrO2 is biocompatible, and the adhesion of bacteria on its surface is low [12–14]. Recently, the ZrO2 nanotubes (ZrNTs) have been synthesized using anodic oxidation techniques [15, 16], which could be reinforcements of biocomposites. Comparing with smooth zirconium surface, the initial adhesion, growth, spreading, functionality of alkaline phosphatase, and the formation of extracellular matrix for cells cultured of ZrNTs are improved [17]. Moreover, the cells attached on the ZrNT surface demonstrate higher cytoskeleton organization compare to flat Zr surface. The mechanical properties of ZrNTs were also studied by theoretical calculation; the basic elastic parameters of single-wall ZrNTs were presented [18]. Given the good biocompatibility and mechanical properties, it is believed that ZrNTs are of great potential as reinforcements of biocomposites.The Mori-Tanaka (M-T) [19–21] and Halpin-Tsai (H-T) [22, 23] models are usually applied to evaluate the elastic modulus of composites, in which their accuracy have been proved with experiments. However, some other fundamental elastic properties of reinforcement on the composites cannot be involved and analysed such as the influence of reinforcement distribution, damage evolution, and failure. It is expected that there is a significant difference on the elastic properties of composites between ordered and unified direction distribution of fibre-like reinforcements. Hence, in this work, the elastic properties (including elastic modulus, shear modulus, and Poisson’s ratio) of ZrO2 nanotube-reinforced HA and β-TCP by random and ordered distribution of reinforcements are firstly studied by applying RVE models [24], which has been widely accepted as an effective approach to study the composites. The enhancement of elastic properties including Young’s modulus, shear modulus, and Poisson’s ratio for random and ordered ZrNT-reinforced composites is compared. The influences on elastic properties of composites caused by the doping volume fraction and aspect ratio L/D of ZrNTs are also investigated. The Halpin-Tsai model is applied to evaluate our calculations of random reinforced composites.


## 2. Computational Method

The RVE models were built with finite element package ANSYS. For random reinforced composites, the RAND command of ANSYS Parametric Design Language (APDL) is applied to define the coordinates and directions of random distributed nanotubes. In the first step, the endpoint and the angle of the central axis of the first nanotube are randomly generated. In the second step, the endpoint and the angle of the central axis, thei-th nanotube, are randomly generated. The new generated nanotube will be reserved if the distance between the central axes of any two nanotubes is larger than the diameter of nanotubes, otherwise it will be deleted. The part of nanotube beyond the boundary move to the other side according to the symmetry. Repeat the second step until getting the calculation model. In particular, when the filling volume ratio is greater than 10%, the overlap is inevitable, and the overlap is eliminated by translating the generated nanotube. For ordered reinforced composites, coordinates of reinforcement nanotubes are randomly generated by RAND command, but the directions are unified. Figure 1 depicts the representative detailed random and ordered distribution RVE models in matrix and perspective views, respectively. The ZrNTs were simplified to be linearly, elastic, isotropic, and homogeneous cylindrical bars [11, 25]. Due to only elastic properties of composites were investigated, the contacts of the ZrNTs and matrix were assumed to be perfect bonded. The 8-node brick solid element was applied to mesh the ZrNTs, and the 10-node tetrahedral solid element was utilized to mesh the matrix. Uniaxial compression test is performed by a finite element method simulation with one surface totally constrained and the opposite surface loaded with pressure for random reinforced composites as shown in Figures 1(a) and 1(b). For ordered reinforced composites as shown in Figures 1(c) and 1(d), owing to its anisotropic feature, two directions of pressures are loaded separately which are along the longitudinal direction and transverse direction of nanotubes. The applied pressures for both random and ordered cases are set to be 0.01GPa=0.01nN/Å2.Figure 1
The RVE models of (a) matrix and (b) composite reinforced by ZrNTs.
(a)
(b)
(c)
(d)The doping volume fraction and aspect ratio of reinforcement nanotubes, which mainly affect the elastic properties of reinforced composites, have been studied. In this work, the volume fraction influence of ZrNTs is calculated in 5.0 vol.%, 5.5 vol.%, 6.5 vol.%, 7.5 vol.%, 8.5 vol.%, and 10.0 vol.%, with the aspect ratioL/D fixed to be 5. The structural parameters in this case are listed in Table 1. As listed, the volume fraction is controlled by different ZrNT quantities. When studying the aspect ratio effect of ZrNTs on the elastic properties, the volume fraction is fixed to be 5.5%, and each RVE has 28 ZrNTs. The aspect ratio varies as 3, 5, 10, and 20. The detailed structural parameters are listed in Table 2. The elastic parameters, including elastic modulus and Poisson’s ratio of reinforcement ZrNTs as well as matrix HA and β-TCP, used in our presented calculations are listed in Table 3.Table 1
Structure parameters for various volume fraction withL/D=5, LZrNTs=300Å, and LRVE=533.6Å.
Volume fraction
5%
5.5%
6.5%
7.5%
8.5%
10.0%
ZrNT quantity
10
11
13
15
17
20Table 2
Structure parameters for various aspect ratio with volume fraction set to be 5.5% andLZrNTs=300Å.
Aspect ratio (L/D)
3
5
10
20
ZrNT diameter (Å)
100
60
30
15
RVE length (Å)
1062.5
755.9
476.2
300Table 3
Material properties of reinforcements and matrix.
Material properties
ZrNTs-armchair
ZrNTs-zigzag
HA
β-TCP
Elastic modulus (GPa)
380 [18]
350 [18]
104 [26]
85 [4]
Poisson’s ratio
0.35-0.45 [18]
0.3-0.4 [18]
0.27 [27]
0.27 [27]In order to assess our results, the Halpin-Tsai (H-T) model [28], which has been demonstrated that the elastic modulus of randomly reinforced composites can be well predicted, is used to evaluate our FEM results for random reinforced case.
(1)EcEm=381+2L/DηLVZrNTs1−ηLVZrNTs+581+2ηTVZrNTs1−ηTVZrNTs,ηL=EZrNTs/EMatrix−1EZrNTs/EMatrix+2L/D,ηT=EZrNTs/EMatrix−1EZrNTs/EMatrix+2,where Ec/Em is the ratio of the elastic modulus of composites and matrix, the L/D (length/diameter for nanotubes) is the aspect ratio for ZrNTs, and the VZrNTs gives the volume fraction of ZrNTs in the matrix.


## 3. Results and Discussion

### 3.1. Isotropy Check

The isotropic feature of RVE models for random reinforced composite has been checked by evaluating the elastic modulus on three main directions. Three volume fraction RVE models were employed (5 vol.%, 5.5 vol.%, and 6.5 vol.%) for ZrNTs-HA composites withL/D set to be 5. The calculated differences of elastic modulus along the three main directions are all less than 1%, which are listed in Table 4. Therefore, the created RVE models with ZrNTs random distribution inside the matrix are all assumed to be isotropic.Table 4
Isotropy check with elastic modulus (in GPa) of three main direction for different reinforcement volume fraction.
5 vol.%
5.5 vol.%
6.5 vol.%
Ex
Ey
Ez
Ex
Ey
Ez
Ex
Ey
Ez
113.9
113.2
112.9
115.4
114.2
114.3
115.5
115.3
116.0


### 3.2. Poisson’s Ratio Effect Check

Poisson’s ratio of ZrNTs is between 0.35 and 0.45 for armchair type and between 0.3 and 0.4 for zigzag type. A random reinforced RVE model with HA as matrix and armchair ZrNTs as reinforcements is employed to study the effect of Poisson’s ratio on the elastic modulus and Poisson’s ratio of the composite. The aspect ratioL/D of ZrNTs is set to be 5, and the doping volume fraction is 5.5 vol.%. The obtained elastic modulus and Poisson’s ratio of composite corresponding to various Poisson’s ratio of ZrNTs are depicted in Figure 2. As is shown, for the Poisson’s ratio of ZrNTs ranged from 0.35 to 0.45, the difference of elastic modulus is less than 0.2 GPa, and the Poisson’s ratio almost keeps a constant of 0.31. Hence, there is slightly an influence on the composite elastic modulus and Poisson’s ratio with different ZrNT Poisson’s ratios. Therefore, in the presented work below, 0.4 is utilized as the Poisson’s ratio of ZrNTs.Figure 2
The elastic modulus (black axis) and Poisson’s ratio (red axis) of composite changed with different Poisson’s ratio of ZrNTs.


### 3.3. Chirality Effect Check

The chirality effect of ZrNTs for both armchair and zigzag types is also checked with a random reinforced RVE model by fixing the doping volume fraction and the aspect ratio. The calculated elastic modulus, shear modulus, and Poisson’s ratio of reinforced HA andβ-TCP are all listed in Table 5. As shown, the elastic modulus, shear modulus, and Poisson’s ratio of armchair and zigzag ZrNT-reinforced composites are quite similar. Hence, the chirality effect of ZrNTs on the elastic properties of composites is negligible, and the investigation of ZrNTs-HA and ZrNT-β-TCP composites below can be performed with a single chirality type of ZrNTs. The armchair ZrNTs were employed in the investigations below.Table 5
The elastic modulus E (in GPa), shear modulus G (in GPa), and Poisson’s ratioμ of composite HA and β-TCP reinforced with armchair- and zigzag-type ZrNTs. The L/D of ZrNTs is set to be 5, and the doping volume fraction is 5 vol.%.
ZrNTs
EZrNTs−HA
EZrNTs−β−TCP
GZrNTs−HA
Gβ−TCP
μZrNTs−HA
μZrNTs−β−TCP
Armchair
113.9
93.9
43.4
35.7
0.311
0.315
Zigzag
113.5
93.6
43.3
35.6
0.310
0.313


### 3.4. Elastic Properties Investigation

The elastic properties, such as elastic modulus, shear modulus, and Poisson’s ratio, of the composites are the most significant features for their applications. In particular, for biocomposites, their elastic properties have stricter standards. For instance, the dental implants need high strength. The bone repair scaffolds and skeletal implants need similar Poisson’s ratio with the target hard tissues, so that health problems caused by the deformation unmatch can be avoided. As is known, the doping volume fraction and aspect ratio of reinforcement nanotubes are the key factors that influence the elastic properties of composites. The effects of ZrNT volume fraction and aspect ratio on the elastic properties of the composites were studied. Figures3(a) and 3(b), respectively, describe the elastic moduli of the random and ordered reinforced composites HA and β-TCP, which correspond to different volume fractions and aspect ratios of ZrNTs. For convenience, the elastic modulus of random reinforced composites is represented by ERandom; the elastic modulus along the longitudinal and transverse direction of ZrNTs in the ordered reinforced composites is represented by ELongitudinal and ETransverse, respectively. The shear modulus and Poisson’s ratio in the below discussion will also follow this definition.Figure 3
Elastic modulus E of random and ordered reinforced composites HA andβ-TCP together with the H-T results affected by different doping (a) volume fraction of ZrNTs and (b) aspect ratio of ZrNTs. Where the grey background region is ZrNTs-HA for both random and ordered distribution, the white background region is ZrNTs-β-TCP for both random and ordered distribution.
(a)
(b)As is shown in Figures3(a) and 3(b), owing to the ELongitudinal is apparently larger than ETransverse, the isotropic feature is observed for ordered distribution ZrNT-reinforced composites. Comparing the random reinforced composites and ordered reinforced composites, for each volume fraction and aspect ratio, a general rule is obtained as ELongitudinal>ERandom>ETransverse. When comparing the ERandom of RVE and H-T model, our results are in agreement with the H-T model for the influence of both volume fraction and aspect ratio, which provide the reliability of our simulation. In addition, an almost linear increase of elastic modulus for both random and ordered reinforced HA and β-TCP is obtained with the volume fraction increase. With the volume fraction increase from 5% to 10%, as is shown in Figure 3(a), the difference between ELongitudinal and ERandom for both reinforced HA and β-TCP also increases from about 5 GPa to 10 GPa. Meanwhile, the difference between ERandom and ETransverse is about 2.5 GPa. While for various aspect ratio effects on both random and ordered reinforced HA and β-TCP, as is shown in Figure 3(b), there are tiny changes for elastic modulus with the aspect ratio change. Furthermore, a convergence trend for ELongitudinal, ERandom, and ETransverse is observed when L/D>10. Hence, for composites with random or ordered reinforced composites, the volume fraction is the dominating factor to enhance the elastic modulus of composites, whereas the aspect ratio has less contribution on elastic modulus comparing to the volume fraction. When the ZrNTs distribute in a unified direction, a better elastic modulus enhancement will be obtained along the unified direction comparing to random distribution. Comparing the elastic modulus along the enhanced direction of ordered reinforced composites to the random reinforced composites, more than 5 GPa higher is achieved, and the volume fraction increase will induce an increase on their difference.According to our results, for random ZrNTs distribution composites, the elastic modulus with 5.5 vol.% ZrNT-reinforced HA andβ-TCP are, respectively, 115.4 GPa and 95.2 GPa, while for ordered reinforced composites, the ELongitudinal is 120.1 GPa for ZrNTs-HA and 100.6 GPa for ZrNTs-β-TCP. Comparing with the BNNT-reinforced HA (5.5 vol.%) and β-TCP (5.3 vol.%), of which the elastic modulus is, respectively, 129.1 GPa and 107.7 GPa [25], the ZrNT-ordered reinforced composites can almost achieve a similar enhancement as BNNT-reinforced HA and β-TCP. Furthermore, due to the BNNTs induce a modest reduction in cell viability [29, 30] while the ZrNTs provide an enhanced osteoblast response [17, 30], the higher biocompatibility of ZrNTs makes it a better candidate as reinforcement of biocomposites.Figure4 shows the shear modulus of random and ordered reinforced HA and β-TCP as a function of volume fraction and aspect ratio. Similar to elastic modulus, the GLongitudinal>GRandom>GTransverse for both random and ordered composites with different volume fractions and aspect ratios is found, as shown in Figures 4(a) and 4(b). The increase trend also observed for shear modulus with the volume fraction increase, as shown in Figure 4(a). The difference between GLongitudinal and GRandom varies from 1.4 GPa to 3.2 GPa for reinforced HA and 1.7 GPa to 3.3 GPa for reinforced β-TCP. The volume fraction increase will also increase their difference. For aspect ratio influence on shear modulus, as is shown in Figure 4(b), the GRandom slightly increases with the L/D increases for both ZrNTs-HA and ZrNTs-β-TCP. The GTransverse almost keeps constant for reinforced HA and β-TCP which is 43.5 GPa and 36 GPa, respectively. While for GLongitudinal, when L/D increase from 3 to 5, the GLongitudinal also increase from 44.2 GPa to 45.1 GPa for reinforced HA and from 36.6 GPa to 37.5 GPa for reinforced β-TCP. Thereafter, the GLongitudinal keeps constant at around 45.0 GPa and 37.6 GPa, respectively, for reinforced HA and β-TCP. Therefore, the volume fraction of reinforcement will apparently enhance the shear modulus of composites, but the aspect ratio will induce little influence. In general, more than 1.4 GPa and 1.7 GPa higher for shear modulus are achieved for reinforced HA and β-TCP along the enhanced direction of ordered distribution comparing to random distribution.Figure 4
Shear modulus G of random and ordered reinforced composites HA andβ-TCP influenced by various doping (a) volume fractions and (b) aspect ratios L/D of ZrNTs.
(a)
(b)The enhancement of elastic modulus and shear modulus for composites comparing to matrix are expressed asEnhancement=Ecomposite−Ematrix/Ematrix or Enhancement=Gcomposite−Gmatrix/Gmatrix, since the volume fraction of reinforcements is the key parameter for reinforcing composites. The enhancements for ZrNTs distributing with random and ordered type in HA and β-TCP are discussed only for different ZrNT volume fractions. Comparing the enhancement of random and ordered reinforced composites, the enhancement of ELongitudinal is at least 5% stronger than ERandom for composite HA and at least 7% stronger for composite β-TCP. However, comparing the enhancement of ETransverse to ERandom for both composites, only at most 1.6% shorten is achieved. Besides, for shear modulus, the enhancements of GLongitudinal are at least 2.8% and 4% larger than GRandom for HA and β-TCP, respectively. And the enhancements of GRandom are at most 1.5% and 3.1% smaller comparing to GTransverse for HA and β-TCP, respectively. Hence, adjusting the reinforcement distribution is an effective way to enhance the elastic properties of composites for specific applications.The Poisson’s ratio is particularly significant especially for biocomposites. As the implants of human body for biocomposites, the unified deformation behaviour with targeted tissue under same load is necessary for healthy. The Poisson’s ratioμ of different composites with random and ordered reinforcement distribution affected by volume fraction and L/D is further investigated. As is shown in Figure 5, the general rule for random and ordered reinforced composites is also observed for both Figures 5(a) and 5(b) which is μLongitudinal>μRandom>μTransverse. Though the elastic modulus and shear modulus are apparently different for composite HA and β-TCP, the Poisson’s ratio for both composites are quite similar. Especially for random-type composites, the Poisson’s ratio for both HA and β-TCP is about 0.31 with different volume fractions of ZrNTs, as is shown in Figure 5(a). When considering the aspect ratio effect on Poisson’s ratio of both HA and β-TCP, as is shown in Figure 5(b), the μRandom is about 0.305 for different L/D. In this case, the obtained Poisson’s ratio of ZrNTs reinforced HA and β-TCP by random distribution is very acceptable due to its consistency with the Poisson’s ratio of dentine and bone (μ=0.31) [31, 32]. However, for substitute bone applications, the biocomposites usually should be modelled as porous implants such as bone scaffold in order to have a good bioactivity. Hence, when considering the porosity for clinical applications, the Poisson’s ratio is expected to be higher. The μTransverse for ordered reinforced composites will probably be designed to be acceptable. In this case, the loading direction will be along the enhanced direction of ordered reinforced composites.Figure 5
Poisson’s ratioμ of reinforced composites HA and β-TCP affected by different doping (a) volume fraction and (b) aspect ratio L/D of ZrNTs.
(a)
(b)One of the advantages of RVE method is to check the stress distribution for composites which can go further to investigate the damage evolution and failure. The first principal stress, which belongs to the Maximum principal stress theory and usually used to evaluate the failure position for brittle materials, is applied to check the stress distribution details inside the composites. Figure6 depicts the first principal stress distributions of ZrNTs and matrix in the ordered and random reinforced HA and β-TCP for typical selected cases, and the other cases have similar results. Two general rules are found from the first principal stress of ZrNTs. (i) The maximum stress usually occurs at the end of the nanotubes especially on the edge of the end due to the concentration of the stress. (ii) The maximum first principal stresses for ZrNTs-β-TCP are smaller than which of ZrNTs-HA under same displacement load. This is caused by a lower elastic modulus of matrix material β-TCP compared to HA.Figure 6
The first principal stress distributions of ZrNTs and matrix in the ordered and random reinforced HA andβ-TCP under the same displacement load. The positive and negative signs are tensile and compressive stress, respectively.


## 3.1. Isotropy Check

The isotropic feature of RVE models for random reinforced composite has been checked by evaluating the elastic modulus on three main directions. Three volume fraction RVE models were employed (5 vol.%, 5.5 vol.%, and 6.5 vol.%) for ZrNTs-HA composites withL/D set to be 5. The calculated differences of elastic modulus along the three main directions are all less than 1%, which are listed in Table 4. Therefore, the created RVE models with ZrNTs random distribution inside the matrix are all assumed to be isotropic.Table 4
Isotropy check with elastic modulus (in GPa) of three main direction for different reinforcement volume fraction.
5 vol.%
5.5 vol.%
6.5 vol.%
Ex
Ey
Ez
Ex
Ey
Ez
Ex
Ey
Ez
113.9
113.2
112.9
115.4
114.2
114.3
115.5
115.3
116.0


## 3.2. Poisson’s Ratio Effect Check

Poisson’s ratio of ZrNTs is between 0.35 and 0.45 for armchair type and between 0.3 and 0.4 for zigzag type. A random reinforced RVE model with HA as matrix and armchair ZrNTs as reinforcements is employed to study the effect of Poisson’s ratio on the elastic modulus and Poisson’s ratio of the composite. The aspect ratioL/D of ZrNTs is set to be 5, and the doping volume fraction is 5.5 vol.%. The obtained elastic modulus and Poisson’s ratio of composite corresponding to various Poisson’s ratio of ZrNTs are depicted in Figure 2. As is shown, for the Poisson’s ratio of ZrNTs ranged from 0.35 to 0.45, the difference of elastic modulus is less than 0.2 GPa, and the Poisson’s ratio almost keeps a constant of 0.31. Hence, there is slightly an influence on the composite elastic modulus and Poisson’s ratio with different ZrNT Poisson’s ratios. Therefore, in the presented work below, 0.4 is utilized as the Poisson’s ratio of ZrNTs.Figure 2
The elastic modulus (black axis) and Poisson’s ratio (red axis) of composite changed with different Poisson’s ratio of ZrNTs.


## 3.3. Chirality Effect Check

The chirality effect of ZrNTs for both armchair and zigzag types is also checked with a random reinforced RVE model by fixing the doping volume fraction and the aspect ratio. The calculated elastic modulus, shear modulus, and Poisson’s ratio of reinforced HA andβ-TCP are all listed in Table 5. As shown, the elastic modulus, shear modulus, and Poisson’s ratio of armchair and zigzag ZrNT-reinforced composites are quite similar. Hence, the chirality effect of ZrNTs on the elastic properties of composites is negligible, and the investigation of ZrNTs-HA and ZrNT-β-TCP composites below can be performed with a single chirality type of ZrNTs. The armchair ZrNTs were employed in the investigations below.Table 5
The elastic modulus E (in GPa), shear modulus G (in GPa), and Poisson’s ratioμ of composite HA and β-TCP reinforced with armchair- and zigzag-type ZrNTs. The L/D of ZrNTs is set to be 5, and the doping volume fraction is 5 vol.%.
ZrNTs
EZrNTs−HA
EZrNTs−β−TCP
GZrNTs−HA
Gβ−TCP
μZrNTs−HA
μZrNTs−β−TCP
Armchair
113.9
93.9
43.4
35.7
0.311
0.315
Zigzag
113.5
93.6
43.3
35.6
0.310
0.313


## 3.4. Elastic Properties Investigation

The elastic properties, such as elastic modulus, shear modulus, and Poisson’s ratio, of the composites are the most significant features for their applications. In particular, for biocomposites, their elastic properties have stricter standards. For instance, the dental implants need high strength. The bone repair scaffolds and skeletal implants need similar Poisson’s ratio with the target hard tissues, so that health problems caused by the deformation unmatch can be avoided. As is known, the doping volume fraction and aspect ratio of reinforcement nanotubes are the key factors that influence the elastic properties of composites. The effects of ZrNT volume fraction and aspect ratio on the elastic properties of the composites were studied. Figures3(a) and 3(b), respectively, describe the elastic moduli of the random and ordered reinforced composites HA and β-TCP, which correspond to different volume fractions and aspect ratios of ZrNTs. For convenience, the elastic modulus of random reinforced composites is represented by ERandom; the elastic modulus along the longitudinal and transverse direction of ZrNTs in the ordered reinforced composites is represented by ELongitudinal and ETransverse, respectively. The shear modulus and Poisson’s ratio in the below discussion will also follow this definition.Figure 3
Elastic modulus E of random and ordered reinforced composites HA andβ-TCP together with the H-T results affected by different doping (a) volume fraction of ZrNTs and (b) aspect ratio of ZrNTs. Where the grey background region is ZrNTs-HA for both random and ordered distribution, the white background region is ZrNTs-β-TCP for both random and ordered distribution.
(a)
(b)As is shown in Figures3(a) and 3(b), owing to the ELongitudinal is apparently larger than ETransverse, the isotropic feature is observed for ordered distribution ZrNT-reinforced composites. Comparing the random reinforced composites and ordered reinforced composites, for each volume fraction and aspect ratio, a general rule is obtained as ELongitudinal>ERandom>ETransverse. When comparing the ERandom of RVE and H-T model, our results are in agreement with the H-T model for the influence of both volume fraction and aspect ratio, which provide the reliability of our simulation. In addition, an almost linear increase of elastic modulus for both random and ordered reinforced HA and β-TCP is obtained with the volume fraction increase. With the volume fraction increase from 5% to 10%, as is shown in Figure 3(a), the difference between ELongitudinal and ERandom for both reinforced HA and β-TCP also increases from about 5 GPa to 10 GPa. Meanwhile, the difference between ERandom and ETransverse is about 2.5 GPa. While for various aspect ratio effects on both random and ordered reinforced HA and β-TCP, as is shown in Figure 3(b), there are tiny changes for elastic modulus with the aspect ratio change. Furthermore, a convergence trend for ELongitudinal, ERandom, and ETransverse is observed when L/D>10. Hence, for composites with random or ordered reinforced composites, the volume fraction is the dominating factor to enhance the elastic modulus of composites, whereas the aspect ratio has less contribution on elastic modulus comparing to the volume fraction. When the ZrNTs distribute in a unified direction, a better elastic modulus enhancement will be obtained along the unified direction comparing to random distribution. Comparing the elastic modulus along the enhanced direction of ordered reinforced composites to the random reinforced composites, more than 5 GPa higher is achieved, and the volume fraction increase will induce an increase on their difference.According to our results, for random ZrNTs distribution composites, the elastic modulus with 5.5 vol.% ZrNT-reinforced HA andβ-TCP are, respectively, 115.4 GPa and 95.2 GPa, while for ordered reinforced composites, the ELongitudinal is 120.1 GPa for ZrNTs-HA and 100.6 GPa for ZrNTs-β-TCP. Comparing with the BNNT-reinforced HA (5.5 vol.%) and β-TCP (5.3 vol.%), of which the elastic modulus is, respectively, 129.1 GPa and 107.7 GPa [25], the ZrNT-ordered reinforced composites can almost achieve a similar enhancement as BNNT-reinforced HA and β-TCP. Furthermore, due to the BNNTs induce a modest reduction in cell viability [29, 30] while the ZrNTs provide an enhanced osteoblast response [17, 30], the higher biocompatibility of ZrNTs makes it a better candidate as reinforcement of biocomposites.Figure4 shows the shear modulus of random and ordered reinforced HA and β-TCP as a function of volume fraction and aspect ratio. Similar to elastic modulus, the GLongitudinal>GRandom>GTransverse for both random and ordered composites with different volume fractions and aspect ratios is found, as shown in Figures 4(a) and 4(b). The increase trend also observed for shear modulus with the volume fraction increase, as shown in Figure 4(a). The difference between GLongitudinal and GRandom varies from 1.4 GPa to 3.2 GPa for reinforced HA and 1.7 GPa to 3.3 GPa for reinforced β-TCP. The volume fraction increase will also increase their difference. For aspect ratio influence on shear modulus, as is shown in Figure 4(b), the GRandom slightly increases with the L/D increases for both ZrNTs-HA and ZrNTs-β-TCP. The GTransverse almost keeps constant for reinforced HA and β-TCP which is 43.5 GPa and 36 GPa, respectively. While for GLongitudinal, when L/D increase from 3 to 5, the GLongitudinal also increase from 44.2 GPa to 45.1 GPa for reinforced HA and from 36.6 GPa to 37.5 GPa for reinforced β-TCP. Thereafter, the GLongitudinal keeps constant at around 45.0 GPa and 37.6 GPa, respectively, for reinforced HA and β-TCP. Therefore, the volume fraction of reinforcement will apparently enhance the shear modulus of composites, but the aspect ratio will induce little influence. In general, more than 1.4 GPa and 1.7 GPa higher for shear modulus are achieved for reinforced HA and β-TCP along the enhanced direction of ordered distribution comparing to random distribution.Figure 4
Shear modulus G of random and ordered reinforced composites HA andβ-TCP influenced by various doping (a) volume fractions and (b) aspect ratios L/D of ZrNTs.
(a)
(b)The enhancement of elastic modulus and shear modulus for composites comparing to matrix are expressed asEnhancement=Ecomposite−Ematrix/Ematrix or Enhancement=Gcomposite−Gmatrix/Gmatrix, since the volume fraction of reinforcements is the key parameter for reinforcing composites. The enhancements for ZrNTs distributing with random and ordered type in HA and β-TCP are discussed only for different ZrNT volume fractions. Comparing the enhancement of random and ordered reinforced composites, the enhancement of ELongitudinal is at least 5% stronger than ERandom for composite HA and at least 7% stronger for composite β-TCP. However, comparing the enhancement of ETransverse to ERandom for both composites, only at most 1.6% shorten is achieved. Besides, for shear modulus, the enhancements of GLongitudinal are at least 2.8% and 4% larger than GRandom for HA and β-TCP, respectively. And the enhancements of GRandom are at most 1.5% and 3.1% smaller comparing to GTransverse for HA and β-TCP, respectively. Hence, adjusting the reinforcement distribution is an effective way to enhance the elastic properties of composites for specific applications.The Poisson’s ratio is particularly significant especially for biocomposites. As the implants of human body for biocomposites, the unified deformation behaviour with targeted tissue under same load is necessary for healthy. The Poisson’s ratioμ of different composites with random and ordered reinforcement distribution affected by volume fraction and L/D is further investigated. As is shown in Figure 5, the general rule for random and ordered reinforced composites is also observed for both Figures 5(a) and 5(b) which is μLongitudinal>μRandom>μTransverse. Though the elastic modulus and shear modulus are apparently different for composite HA and β-TCP, the Poisson’s ratio for both composites are quite similar. Especially for random-type composites, the Poisson’s ratio for both HA and β-TCP is about 0.31 with different volume fractions of ZrNTs, as is shown in Figure 5(a). When considering the aspect ratio effect on Poisson’s ratio of both HA and β-TCP, as is shown in Figure 5(b), the μRandom is about 0.305 for different L/D. In this case, the obtained Poisson’s ratio of ZrNTs reinforced HA and β-TCP by random distribution is very acceptable due to its consistency with the Poisson’s ratio of dentine and bone (μ=0.31) [31, 32]. However, for substitute bone applications, the biocomposites usually should be modelled as porous implants such as bone scaffold in order to have a good bioactivity. Hence, when considering the porosity for clinical applications, the Poisson’s ratio is expected to be higher. The μTransverse for ordered reinforced composites will probably be designed to be acceptable. In this case, the loading direction will be along the enhanced direction of ordered reinforced composites.Figure 5
Poisson’s ratioμ of reinforced composites HA and β-TCP affected by different doping (a) volume fraction and (b) aspect ratio L/D of ZrNTs.
(a)
(b)One of the advantages of RVE method is to check the stress distribution for composites which can go further to investigate the damage evolution and failure. The first principal stress, which belongs to the Maximum principal stress theory and usually used to evaluate the failure position for brittle materials, is applied to check the stress distribution details inside the composites. Figure6 depicts the first principal stress distributions of ZrNTs and matrix in the ordered and random reinforced HA and β-TCP for typical selected cases, and the other cases have similar results. Two general rules are found from the first principal stress of ZrNTs. (i) The maximum stress usually occurs at the end of the nanotubes especially on the edge of the end due to the concentration of the stress. (ii) The maximum first principal stresses for ZrNTs-β-TCP are smaller than which of ZrNTs-HA under same displacement load. This is caused by a lower elastic modulus of matrix material β-TCP compared to HA.Figure 6
The first principal stress distributions of ZrNTs and matrix in the ordered and random reinforced HA andβ-TCP under the same displacement load. The positive and negative signs are tensile and compressive stress, respectively.


## 4. Conclusions and Recommendations

The elastic properties of ZrNT-reinforced HA andβ-TCP in both random and ordered distributions were investigated by finite element method with RVE models. The elastic modulus, shear modulus, and Poisson’s ratio of different composites with different distribution types were obtained with various volume fractions and aspect ratios of ZrNTs. Our results suggest that
(1)
The volume fraction of reinforcement plays a domination role comparing to the aspect ratio for the elastic properties of two composites in both random and ordered distribution type. For the elastic properties of random and ordered reinforced composites, a general rule is obtained which can be expressed asE;G;μLongitudinal>E;G;μRandom>E;G;μTransverse. Under the same volume fraction and aspect ratio of reinforcements, the best enhancement is obtained along the longitudinal direction of ordered reinforced composites
(2)
For elastic modulus, the enhancement ofELongitudinal is at least 5% stronger than ERandom for composite HA and at least 7% stronger for composite β-TCP. Only at most 1.6% shorten is achieved for the enhancement of ETransverse comparing to ERandom for both HA and β-TCP
(3)
For shear modulus, the enhancements ofGLongitudinal are at least 2.8% and 4% larger than GRandom for HA and β-TCP, respectively. And the enhancements of GRandom are at most 1.5% and 3.1% smaller comparing to GTransverse for HA and β-TCP, respectively
(4)
In summary, considering the above advantages of ordered reinforced biocomposites, especially the Poisson’s ratio, the ordered ZrO2 nanotubes reinforced type is more designable and acceptable for bone repair and substitute applications



---

*Source: 4329717-2020-03-13.xml*