Evaluation of Thermal-Mechanical Properties of Bio-Oil Regenerated Aged Asphalt

Different proportions of bio-oil (5, 10, 15, and 20 wt%) were added into aged asphalt for its regeneration. Molecular dynamic simulations were used to measure the thermal and mechanical performances of bio-oil regenerated aged asphalt (BRAA). A new, simplified BRAA model was built to calculate the specific heat capacity, thermal expansion coefficient, elastic constant, shear modulus, bulk modulus, and Young’s modulus. Simulation results showed that the thermal expansion coefficient (CTE α) of asphalt at 298 K decreased by 10% after aging. Bio-oil of 5 wt% could make the CTE α restore to the original level of base asphalt, while the addition of bio-oil would further decrease the specific heat capacity of aged asphalt. The shear modulus (G), Young’s modulus (K) and bulk modulus (E) of asphalt increased after aging and decreased with the increasing amount of bio-oil. According to the calculated E/G value, the ductility of aged asphalt increased by 6.0% with the addition of 10 wt% bio-oil, while over 15 wt% bio-oil would make the ductility of BRAA decrease. In summary, the regeneration effects of bio-oil to the thermal expansion coefficient, flexibility, and ductility of aged asphalt had been proven, while excessive bio-oil would decrease the thermal stability of asphalt.


Introduction
Asphalt mixture is a common pavement material. It provides a comfortable driving environment for vehicles. However, the pavement performance of asphalt mixture constantly decreases during service process. One of the key reasons is that asphalt, as a binder, is gradually aging under the influence of oxygen and ultraviolet [1][2][3]. In order to recycle and reuse the aged asphalt, a considerable amount of research has been done on the regeneration of aged asphalt. In 1915, asphalt recycling technology was first explored by American Warren Bother who reused the aged asphalt by heating waste asphalt mixture [4]. After 1956, due to the initial formation of the US domestic high-grade road network, the price of petroleum crude oil rose rapidly [5]. During the period from 1970 to 1989, the global oil crisis occurred and crude oil price continued to rise. The cost of new asphalt pavement increased greatly, thus, people paid more attention to asphalt recycling technology. Regarding the recycling mechanism of aged asphalt, the regeneration of aged asphalt can be considered as the reverse process of asphalt aging. By means of adding specific substances, the composition of asphalt is adjusted to reverse its performance and return to an appropriate state that can be reused [6][7][8]. Not only rheological properties but also changes in structure were considered as the indicators for the evaluation of asphalt regeneration effect. More testing methods, such as NMR, were used to investigate the additive's effects on aged asphalt [9]. The main components of bio-oil extracted from waste wood were acetic acid, 1-carboxy-2propanone, and methanol [35]. Bio-oil was mixed with aged asphalt to prepare BRAA by the mass fraction of 5%, 10%, 15%, and 20% to asphalt, respectively. Finally, six asphalt systems were obtained. Their compositions are shown in Table 1. The simulations were performed using Materials Studio, the software for simulating and modeling materials. Firstly, 3D models of the molecules showed above were constructed. Then, amorphous cell module was used to mix the components and to build the asphalt systems. The initial densities were set to be 1.0 g/cm 3 . Pressure was set at 101.325 kPa (1.0 atm). The asphalt systems were first relaxed to induce the system energy using the geometry optimization tool of Forcite module, then relaxed to state equilibrium using the dynamic tool with the parameters set to be isobaricisothermal ensemble (NPT), 1.0 fs time step and 200 ps simulation time. After that, the systems were relaxed to state equilibrium with canonical ensemble (NVT), 1.0 fs time step and 200 ps simulation time. Condensed-phase optimized molecular potentials for atomistic simulation studies (COMPASS) force filed were chosen in all molecular simulations and minimization processes in this research. The following processes were performed to simulate the physical state of asphalt systems and to calculate the thermal-mechanical properties. Simulation temperature ranged from 223 to 448 K. The six asphalt systems are shown in Figure 2. The main components of bio-oil extracted from waste wood were acetic acid, 1-carboxy-2-propanone, and methanol [35]. Bio-oil was mixed with aged asphalt to prepare BRAA by the mass fraction of 5%, 10%, 15%, and 20% to asphalt, respectively. Finally, six asphalt systems were obtained. Their compositions are shown in Table 1. The simulations were performed using Materials Studio, the software for simulating and modeling materials. Firstly, 3D models of the molecules showed above were constructed. Then, amorphous cell module was used to mix the components and to build the asphalt systems. The initial densities were set to be 1.0 g/cm 3 . Pressure was set at 101.325 kPa (1.0 atm). The asphalt systems were first relaxed to induce the system energy using the geometry optimization tool of Forcite module, then relaxed to state equilibrium using the dynamic tool with the parameters set to be isobaric-isothermal ensemble (NPT), 1.0 fs time step and 200 ps simulation time. After that, the systems were relaxed to state equilibrium with canonical ensemble (NVT), 1.0 fs time step and 200 ps simulation time. Condensed-phase optimized molecular potentials for atomistic simulation studies (COMPASS) force filed were chosen in all molecular simulations and minimization processes in this research. The following processes were performed to simulate the physical state of asphalt systems and to calculate the thermal-mechanical properties. Simulation temperature ranged from 223 to 448 K. The six asphalt systems are shown in Figure 2.

Specific Heat Capacity
Specific heat capacity (Cp) is a measurable physical quantity equal to the ratio of the heat added to (or removed from) an object to the resulting temperature change [36]. During the heating process, the higher specific heat capacity slowed the asphalt temperature rise more slowly, which means that more energy will be consumed during asphalt heating. Otherwise, asphalt with higher specific heat capacity will be less susceptible to ambient temperature.
In this research, asphalt systems were relaxed to state equilibrium at the same pressure (1.0 atm) and different temperatures (223 to 448 K). The systems would have different properties at different temperatures. The specific heat capacity of asphalt could be analyzed by comparing the enthalpies of asphalt systems at different temperatures, as follows: where Cp is specific heat; E is energy; P, V, and T are pressure, volume, and temperature, respectively; and is Boltzmann content. The E, P, V, and T could be obtained from result texts after dynamic simulations.

Thermal Expansion Coefficient
Thermal expansion is the tendency of matter to change in shape, area, and volume in response to a change in temperature [37]. In this research, asphalt systems were relaxed to state equilibrium at different temperature (223 to 448 K) to simulate the heating process of asphalt. With the change in temperature, the volumes of asphalt systems simultaneously changed. The volumetric thermal expansion coefficient (CTE α) could be calculated by the following equation:

Specific Heat Capacity
Specific heat capacity (C p ) is a measurable physical quantity equal to the ratio of the heat added to (or removed from) an object to the resulting temperature change [36]. During the heating process, the higher specific heat capacity slowed the asphalt temperature rise more slowly, which means that more energy will be consumed during asphalt heating. Otherwise, asphalt with higher specific heat capacity will be less susceptible to ambient temperature.
In this research, asphalt systems were relaxed to state equilibrium at the same pressure (1.0 atm) and different temperatures (223 to 448 K). The systems would have different properties at different temperatures. The specific heat capacity of asphalt could be analyzed by comparing the enthalpies of asphalt systems at different temperatures, as follows: where C p is specific heat; E is energy; P, V, and T are pressure, volume, and temperature, respectively; and k B is Boltzmann content. The E, P, V, and T could be obtained from result texts after dynamic simulations.

Thermal Expansion Coefficient
Thermal expansion is the tendency of matter to change in shape, area, and volume in response to a change in temperature [37]. In this research, asphalt systems were relaxed to state equilibrium at different temperature (223 to 448 K) to simulate the heating process of asphalt. With the change in temperature, the volumes of asphalt systems simultaneously changed. The volumetric thermal expansion coefficient (CTE α) could be calculated by the following equation: where V, T, and P are volume, temperature, and pressure, respectively. The V, T, and P could be obtained from result texts after dynamic simulations.

Elastic Constant
After the asphalt systems relaxed to state equilibrium, the elastic constants of systems were calculated by Mechanical Property tool of Forcite module. A GPa of 0.03 was chosen to be the maximum strain applying to the asphalt systems in this simulation. For elastic materials, the elasticity could be evaluated by Hook's law: where i, j = 1, 2, 3. σ i is the stress vector, while ε j is the strain vector; C ij is the six-dimensional stiffness matrix. The stress components could be calculated as follows: where V is the volume; m k and u k are the mass and velocity of the kth particle, respectively; r is the distance between the kth and the first particles; f is the force exerted on the first particle by the kth particle. Lame coefficient λ and µ can be calculated by the following:

Young's Modulus, Shear Modulus, and Bulk Modulus
Young's modulus, denoted by K, is a mechanical property of linear elastic solid materials, measuring the stiffness of a solid material. It defines the relationship between stress and strain in the material, which will be more or less dependent on temperature. Shear modulus, always denoted by G, is defined as the ratio of shear stress to the shear strain [38]. Bulk modulus, denoted by E, is defined as the ratio of the infinitesimal pressure increase to the resulting relative decrease of the volume. They could be calculated as follows: where K, G, and E are Young's modulus, shear modulus, and bulk modulus, respectively; λ and µ are Lame coefficients.

Thermal Properties of BRAA
The parameters of asphalt, aged asphalt, and BRAA were calculated by molecular dynamic simulations at different temperatures. As shown in Figure 3, the specific heat capacity (C p ) of asphalt systems showed an upward trend with temperature rises. Moreover, it decreased with aging condition and the addition of the bio-oil. The results indicated that the bio-oil renders asphalt more sensitive to temperature. Furthermore, aging would weaken the temperature sensitive of asphalt binder, and the relationship between C p and temperature is linear.
Materials 2018, 11, x FOR PEER REVIEW 6 of 12 sensitive to temperature. Furthermore, aging would weaken the temperature sensitive of asphalt binder, and the relationship between Cp and temperature is linear.  Figure 4 showed the thermal expansion coefficient (CTE α) of asphalt, aged asphalt, and BRAA at different temperatures. As can be seen, the aged asphalt had the lowest CTE α, while the 20 wt% BRAA showed the highest. The CTE α of asphalt at 298 K decreased by 10.0% after aging. The addition of 5 wt% bio-oil could increase the CTE α of aged asphalt by 13.4% and make it even higher than that of base asphalt. Excessive bio-oil would further increase the CTE α of asphalt. However, when the mass fraction of bio-oil is more than 10%, the CTE α of BRAA would significantly increase over 350 K. The results indicated that the rejuvenating effect of bio-oil on aged asphalt shows an optimum dosage of 10% regarding the CTE α. Excessive bio-oil would decrease the thermal stability of asphalt and induce thermal expansion.   Figure 4 showed the thermal expansion coefficient (CTE α) of asphalt, aged asphalt, and BRAA at different temperatures. As can be seen, the aged asphalt had the lowest CTE α, while the 20 wt% BRAA showed the highest. The CTE α of asphalt at 298 K decreased by 10.0% after aging. The addition of 5 wt% bio-oil could increase the CTE α of aged asphalt by 13.4% and make it even higher than that of base asphalt. Excessive bio-oil would further increase the CTE α of asphalt. However, when the mass fraction of bio-oil is more than 10%, the CTE α of BRAA would significantly increase over 350 K. The results indicated that the rejuvenating effect of bio-oil on aged asphalt shows an optimum dosage of 10% regarding the CTE α. Excessive bio-oil would decrease the thermal stability of asphalt and induce thermal expansion.
Materials 2018, 11, x FOR PEER REVIEW 6 of 12 sensitive to temperature. Furthermore, aging would weaken the temperature sensitive of asphalt binder, and the relationship between Cp and temperature is linear.  Figure 4 showed the thermal expansion coefficient (CTE α) of asphalt, aged asphalt, and BRAA at different temperatures. As can be seen, the aged asphalt had the lowest CTE α, while the 20 wt% BRAA showed the highest. The CTE α of asphalt at 298 K decreased by 10.0% after aging. The addition of 5 wt% bio-oil could increase the CTE α of aged asphalt by 13.4% and make it even higher than that of base asphalt. Excessive bio-oil would further increase the CTE α of asphalt. However, when the mass fraction of bio-oil is more than 10%, the CTE α of BRAA would significantly increase over 350 K. The results indicated that the rejuvenating effect of bio-oil on aged asphalt shows an optimum dosage of 10% regarding the CTE α. Excessive bio-oil would decrease the thermal stability of asphalt and induce thermal expansion.  As conclusions, the CTE α of aged asphalt could be improved and the C p would be decreased with the addition of bio-oil. Bio-oil of 5 wt% would be the suitable choice to regenerate the aged asphalt. CTE α of aged asphalt could restore the original level with the addition of 5 wt% bio-oil, while the thermal stability of asphalt would decrease when the amount of bio-oil was over 10 wt%.

Mechanical Properties of BRAA
As shown in Figure 5, the lame coefficient λ of asphalt decreased with the increasing temperature. It is demonstrated that temperature could have an influence on the mechanical properties and change the asphalt internal structure. The lame coefficient λ of aged asphalt was the highest, and the λ of BRAA reduced with the increasing bio-oil content, exhibiting that aging could improve the stiffness of asphalt. But bio-oil would also weaken the resistance to elastic deformation. The lame coefficient λ of base asphalt, aged asphalt, 5 wt% BRAA, 10 wt% BRAA, and 15 wt% BRAA at 298 K were 2.60, 3.04, 2.96, 2.84, and 2.02, respectively. The lame coefficient λ of 10 wt% BRAA could decrease by 6.6%, while the lame coefficient λ of 15 wt% BRAA could decrease by 33.6%. The results indicated that the effect of bio-oil on the lame coefficient λ of asphalt was not linearly taking place. Over 15 wt% bio-oil would significantly reduce the resistance of asphalt to elastic deformation. As conclusions, the CTE α of aged asphalt could be improved and the Cp would be decreased with the addition of bio-oil. Bio-oil of 5 wt% would be the suitable choice to regenerate the aged asphalt. CTE α of aged asphalt could restore the original level with the addition of 5 wt% bio-oil, while the thermal stability of asphalt would decrease when the amount of bio-oil was over 10 wt%.

Mechanical Properties of BRAA
As shown in Figure 5, the lame coefficient λ of asphalt decreased with the increasing temperature. It is demonstrated that temperature could have an influence on the mechanical properties and change the asphalt internal structure. The lame coefficient λ of aged asphalt was the highest, and the λ of BRAA reduced with the increasing bio-oil content, exhibiting that aging could improve the stiffness of asphalt. But bio-oil would also weaken the resistance to elastic deformation. The lame coefficient λ of base asphalt, aged asphalt, 5 wt% BRAA, 10 wt% BRAA, and 15 wt% BRAA at 298 K were 2.60, 3.04, 2.96, 2.84, and 2.02, respectively. The lame coefficient λ of 10 wt% BRAA could decrease by 6.6%, while the lame coefficient λ of 15 wt% BRAA could decrease by 33.6%. The results indicated that the effect of bio-oil on the lame coefficient λ of asphalt was not linearly taking place. Over 15 wt% bio-oil would significantly reduce the resistance of asphalt to elastic deformation. As can be seen in Figure 6, the shear moduli (G) increased after asphalt aging, while the G decreased with the increasing bio-oil addition amount. It indicated that the aging and bio-oil would have influence on the shear moduli. In addition, the G decreased with the increasing temperature, and the bio-oil also decreased the shear moduli of asphalt binder. As can be seen, shear moduli reduced in the order of aged asphalt, 5 wt% BRAA, 10 wt% BRAA, and base asphalt. It indicated that the flexibility of the asphalt decreased after aging, and the bio-oil could recover this change to some extent. The shear moduli of BRAA were close to those of base asphalt when the content of bio-oil reached 15 wt%. The G of 20 wt% BRAA was even lower than that of base asphalt.
As shown in Figure 7, the Young's moduli (K) were calculated according to λ and G at different temperatures. The Young's moduli of aged asphalt were highest, and decreased with the addition of bio-oil, while the K of 5 wt% BRAA and 10 wt% BRAA were bigger than that of base asphalt. The results showed that Young's modulus of asphalt increased by 35.2% after aging at 298 K, and adding 5 wt% bio-oil could decrease the Young's modulus of aged asphalt by 10.0%. Bio-oil of 15 wt% could restore the K of aged asphalt restore to its original level. As can be seen in Figure 6, the shear moduli (G) increased after asphalt aging, while the G decreased with the increasing bio-oil addition amount. It indicated that the aging and bio-oil would have influence on the shear moduli. In addition, the G decreased with the increasing temperature, and the bio-oil also decreased the shear moduli of asphalt binder. As can be seen, shear moduli reduced in the order of aged asphalt, 5 wt% BRAA, 10 wt% BRAA, and base asphalt. It indicated that the flexibility of the asphalt decreased after aging, and the bio-oil could recover this change to some extent. The shear moduli of BRAA were close to those of base asphalt when the content of bio-oil reached 15 wt%. The G of 20 wt% BRAA was even lower than that of base asphalt.
As shown in Figure 7, the Young's moduli (K) were calculated according to λ and G at different temperatures. The Young's moduli of aged asphalt were highest, and decreased with the addition of bio-oil, while the K of 5 wt% BRAA and 10 wt% BRAA were bigger than that of base asphalt. The results showed that Young's modulus of asphalt increased by 35.2% after aging at 298 K, and adding 5 wt% bio-oil could decrease the Young's modulus of aged asphalt by 10.0%. Bio-oil of 15 wt% could restore the K of aged asphalt restore to its original level.  As shown in Figure 8, the bulk moduli (E) at different temperatures were calculated from λ and G. The bulk moduli of asphalt increased after aging. Moreover, the bulk moduli decreased with the additional amount of bio-oil. It indicated that aging and bio-oil could affect the bulk moduli of asphalt. The E of base asphalt increased by 26.2% after aging, and 5 wt% bio-oil could decrease the E of aged asphalt by 4.8%. In addition, when the amount of bio-oil was more than 15 wt%, the BRAA had lower bulk moduli than base asphalt.
The ratio of bulk modulus to shear modulus (E/G) was used to estimate the brittle or ductile behavior of the material. High E/G values represent ductility, while low E/G values represent brittleness. The critical value for separating tough and brittle materials is about 1.75. As shown in Figure 9, calculated E/G value indicated that the base asphalt was ductile. The E/G value of asphalt decreased after aging. The addition of 5 wt% or 10 wt% of bio-oil could recover the ductility of aged asphalt to some extent. The E/G values of base asphalt, aged asphalt, 5 wt% BRAA, and 10 wt% BRAA were 2.361, 2.116, 2.159, and 2.244, respectively. The results showed that the ductility of asphalt decreased by 10.4% after aging, and 10 wt% bio-oil could improve the ductility of aged asphalt by  As shown in Figure 8, the bulk moduli (E) at different temperatures were calculated from λ and G. The bulk moduli of asphalt increased after aging. Moreover, the bulk moduli decreased with the additional amount of bio-oil. It indicated that aging and bio-oil could affect the bulk moduli of asphalt. The E of base asphalt increased by 26.2% after aging, and 5 wt% bio-oil could decrease the E of aged asphalt by 4.8%. In addition, when the amount of bio-oil was more than 15 wt%, the BRAA had lower bulk moduli than base asphalt.
The ratio of bulk modulus to shear modulus (E/G) was used to estimate the brittle or ductile behavior of the material. High E/G values represent ductility, while low E/G values represent brittleness. The critical value for separating tough and brittle materials is about 1.75. As shown in Figure 9, calculated E/G value indicated that the base asphalt was ductile. The E/G value of asphalt decreased after aging. The addition of 5 wt% or 10 wt% of bio-oil could recover the ductility of aged asphalt to some extent. The E/G values of base asphalt, aged asphalt, 5 wt% BRAA, and 10 wt% BRAA were 2.361, 2.116, 2.159, and 2.244, respectively. The results showed that the ductility of asphalt decreased by 10.4% after aging, and 10 wt% bio-oil could improve the ductility of aged asphalt by As shown in Figure 8, the bulk moduli (E) at different temperatures were calculated from λ and G. The bulk moduli of asphalt increased after aging. Moreover, the bulk moduli decreased with the additional amount of bio-oil. It indicated that aging and bio-oil could affect the bulk moduli of asphalt. The E of base asphalt increased by 26.2% after aging, and 5 wt% bio-oil could decrease the E of aged asphalt by 4.8%. In addition, when the amount of bio-oil was more than 15 wt%, the BRAA had lower bulk moduli than base asphalt.
The ratio of bulk modulus to shear modulus (E/G) was used to estimate the brittle or ductile behavior of the material. High E/G values represent ductility, while low E/G values represent brittleness. The critical value for separating tough and brittle materials is about 1.75. As shown in Figure 9, calculated E/G value indicated that the base asphalt was ductile. The E/G value of asphalt decreased after aging. The addition of 5 wt% or 10 wt% of bio-oil could recover the ductility of aged asphalt to some extent. The E/G values of base asphalt, aged asphalt, 5 wt% BRAA, and 10 wt% BRAA were 2.361, 2.116, 2.159, and 2.244, respectively. The results showed that the ductility of asphalt decreased by 10.4% after aging, and 10 wt% bio-oil could improve the ductility of aged asphalt by 6.0%. However, 15 wt% bio-oil would further decrease the ductility of aged asphalt. On the other hand, base asphalt had the highest E/G value at 298 K. With the addition of bio-oil, the highest E/G value temperature reduced, which indicated that bio-oil would decrease the temperature stability of asphalt.
Materials 2018, 11, x FOR PEER REVIEW 9 of 12 6.0%. However, 15 wt% bio-oil would further decrease the ductility of aged asphalt. On the other hand, base asphalt had the highest E/G value at 298 K. With the addition of bio-oil, the highest E/G value temperature reduced, which indicated that bio-oil would decrease the temperature stability of asphalt.

Conclusions
Molecular dynamic simulations were performed in order to investigate the thermal-mechanical properties of bio-oil regenerated aged asphalt (BRAA). Some conclusions were drawn: (1) The aging process had a significant impact on the thermal properties of asphalt. Aged asphalt had lower specific heat capacity and thermal expansion coefficient than base asphalt. The addition of 5 wt% bio-oil could make the thermal properties of aged asphalt restore to the original level. When the amount of bio-oil was over 10 wt%, the thermal stability of asphalt Materials 2018, 11, x FOR PEER REVIEW 9 of 12 6.0%. However, 15 wt% bio-oil would further decrease the ductility of aged asphalt. On the other hand, base asphalt had the highest E/G value at 298 K. With the addition of bio-oil, the highest E/G value temperature reduced, which indicated that bio-oil would decrease the temperature stability of asphalt.

Conclusions
Molecular dynamic simulations were performed in order to investigate the thermal-mechanical properties of bio-oil regenerated aged asphalt (BRAA). Some conclusions were drawn: (1) The aging process had a significant impact on the thermal properties of asphalt. Aged asphalt had lower specific heat capacity and thermal expansion coefficient than base asphalt. The addition of 5 wt% bio-oil could make the thermal properties of aged asphalt restore to the original level. When the amount of bio-oil was over 10 wt%, the thermal stability of asphalt Figure 9. The E/G value of asphalt systems.

Conclusions
Molecular dynamic simulations were performed in order to investigate the thermal-mechanical properties of bio-oil regenerated aged asphalt (BRAA). Some conclusions were drawn: (1) The aging process had a significant impact on the thermal properties of asphalt. Aged asphalt had lower specific heat capacity and thermal expansion coefficient than base asphalt. The addition of 5 wt% bio-oil could make the thermal properties of aged asphalt restore to the original level.
When the amount of bio-oil was over 10 wt%, the thermal stability of asphalt would be affected, which made the performance of asphalt decrease rapidly with temperature rising. (2) The aged asphalt had higher elastic constants, Young's modulus, and bulk modulus than base asphalt, which meant asphalt hardened after aging. The addition of bio-oil would soften the aged asphalt, which reduced the possibility of asphalt cracking. 5 wt% of bio-oil could increase the ductility of aged asphalt by 6.0%. As the amount of bio-oil increased, the structure of asphalt would be affected. When the addition of bio-oil was over 15 wt%, the ductility of aged asphalt would decrease instead.
It should be noticed that the model constructed in this study chose only one kind of molecule to represent each component of asphalt. The results of the simulation could only reflect the actual law to some extent, because asphalt is a complex mixture. The three molecules used to represent bio-oil are the three most abundant in bio-oil. Whether the other ingredients have an influence on the regeneration performance of bio-oil is not certain. Approaching the actual simulation, a more complex system should be built and more components need to be considered.

Further Research Work
The thermal-mechanical properties of bio-oil regenerated aged asphalt (BRAA) were evaluated, but the data were only based on computer simulations. More work should be done to further verify the rejuvenating effect of bio-oil. The specific heat capacity can be measured by differential scanning calorimetry (DSC), and the thermal expansion coefficient can be tested by Michelson interferometry or quartz dilatometry. Dynamic shear rheological test (DSR) and Dynamic thermomechanical analysis (DMA) can be carried out to measure the shear modulus (G) and bulk modulus (E) of BRAA. A ductility test can be done to verify the relationship between ductility and E/G value. The testing results will be compared to the simulation results. The results of the two may be different in value but have the same trend. Introducing coefficients and increasing the number of samples can improve the accuracy of the simulation. When the accuracy of the model is verified multiple times, the simulation will become the guiding basis for material design.