# Mechanical Property and Analysis of Asphalt Components Based on Molecular Dynamics Simulation

Authors: Rui Li; Qiqi Guo; Hui Du; Jianzhong Pei
Journal: Journal of Chemistry
 (2017)

---

## Abstract

The asphalt-aggregate interface interaction plays a significant role in the overall performances of asphalt mixture. In order to analyze the chemical constitution of asphalt effects on the asphalt-aggregate interaction, the average structure C64H52S2 is selected to represent the asphalt, and the colloid, saturated phenol, and asphaltene are selected to represent the major constitutions in asphalt. The molecular models are established for the three compositions, respectively, and the Molecular Dynamics (MD) simulation was conducted for the three kinds of asphaltene-aggregate system at different presses. Comparing the E value of Young modulus of these three polymers, the maximum modulus value of asphaltene was 2.80 GPa, the modulus value of colloid was secondary, and the minimum modulus of saturated phenol was 0.52 GPa. This result corresponds to conventional understanding.


---

## Body

## 1. Introduction

At present, scholars at home and abroad have more and more research perspectives on polymer materials [1–3]. Besides, the research process and methods are also diverse. Since the materials prepared indoors and tested will be influenced by the factors of temperature and pressure, there are also some potential uncontrollable variables which will result in the fact that the depth of analytical investigation on experimental data will be limited [4, 5]. However, researchers have extended the material research field continuously thanks to the introduction of Molecular Dynamics method in recent years. With this simulation method, various parameters can be controlled, and even the environment that can not be realized in laboratory can be simulated [6]. Thus, it will be an important tendency of studying polymer materials to apply Molecular Dynamics method to predict and study the polymer materials’ properties gradually [7–9].In order to analyze the chemical constitution of asphalt on the AAI, the average structure C65H74N2S2 is selected in this paper to represent the asphalt and the colloid; saturated phenol and asphaltene are selected to represent the major constitution in asphalt. The molecular models are established for the three compositions, respectively, and the MD simulation was conducted for the three kinds of asphaltene-aggregate system at different presses.


## 2. Modeling

When the force acts on the system and if the whole system is balanced, the external force suffered in the system will have an absolute balance with the internal force generated inside the system. Normally, stress can be expressed as a second-order tensor containing 9 components as shown in(1)σ11σ12σ13σ21σ22σ23σ31σ32σ33.In the process of molecular calculation, internal stress tensor can be expressed by the virial expression of(2)σ=-1V0∑i=1NmiViViT+∑i<jrijfijT.In formula (2), i represents the sequence numbers (from 1 to N) of particles; mi, Vi, and fi represent the mass, speed, and force suffered of number i. V0 refers to the volume of system without deformation.Once stress acts on the molecular system, the position of internal particles of system will change relatively, where the change is expressed by the strain tensor of(3)ε11ε12ε13ε21ε22ε23ε31ε32ε33.For the parallel hexahedron structured in this paper, strain tensor is only determined by the column vectorsa0, b0, and c0 and deformation state vectors a, b, and c under a certain condition, as shown in (4)ε=12h0T-1Gh0-1-1,where h0 represents the matrix structured by vectors a0, b0, and c0. h refers to the matrix structured by deformation state vectors a,b, and c. G represents hTh.Elastic stiffness constant is associated with the different compositions of system stress and strain. With regard to small deformation, the relationship between stress and strain meets Hooke’s law, as shown in(5)σim=Cimnkεnk,where Cimnk refers to the elastic stiffness constant.Since stress tensor and strain tensor have some symmetry, stress formula (1) can be simplified as (6)σ11σ12σ13σ21σ22σ23σ31σ32σ33⟶σ1σ6σ5σ6σ2σ4σ5σ4σ3.Stain formula (3) can be simplified as (7)ε11ε12ε13ε21ε22ε23ε31ε32ε33⟶ε1ε62ε52ε62ε2ε42ε52ε42ε3.Suppose that the materials prepared are isotropic, and the stress and strain only depend on two independent coefficients. Stiffness matrix is shown as(8)λ+2μλλ000λλ+2μλ000λλλ+2μ000000μ000000μ000000μ,where λ and μ in formula (8) are called lame constants.In addition, with regard to isotropic materials, the Yong modulusE, Poisson’s ratio v, bulk modulus K, and shear modulus G can be obtained through lame constants λ and μ, as shown in (9)E=μ3λ+2μλ+μ,ν=λ2λ+μ,K=λ+23μ,G=μ,where Poisson’s ratio can combine with these four parameters, and the relationship is shown in (10)E=2G1+ν=3K1-2ν.


## 3. Results and Discussion

### 3.1. Mechanical Property Simulation Analysis of Colloid

(1) Simulation Process and Result of Molecular Dynamics. The 1,7-dimethylnaphthalene was used as the repetitive unit to represent colloid in this paper. According to the suggestions of literature, macromolecular structure obtained at the polymerization degree in the range of 10~15 was employed [10, 11]. The character after calculation could meet the requirements such as research precision. The degree of polymerization selected for structuring the macromolecular chain of colloid was 15, and the degree of polymerization that was 15 to structure the long chain of colloid is shown in Figure 1.Figure 1
Long chain of colloid for the degree of polymerization at 15.Forcite module was used to analyze geometry optimization for the long chain of colloid. Figure2 is the long chain of colloid after geometry optimization.Figure 2
Long chain of colloid after being geometrically optimized.Afterwards, a cubic vacuum space in60Å×60Å×60Å was structured. The density was set at 0.1 g/cm3, and the initial model of colloid unit cell at the density of 0.1 obtained is shown in Figure 3.Figure 3
Initial model of colloid at the density of 0.1.The density of initial model is 0.1 g/cc, and the density was in the range of 0.99~1.1 g/cm3 at abnormal temperature; however, the low density system of Figure 3 is usually in a high energy state. Thus, structural optimization should be conducted to the system first to lower the potential energy, and then forcite can be used for optimization. When dynamic calculation was conducted, the temperature was 298 K (close to normal temperature), simulation time was 200 ps and pressure was 0.01 GPa, and the initial density of model at 0.1 would be compressed. In the simulation process, temperature, energy, system density, and the length of side of system were changing continually. The output images are shown from Figures 4–7, where the system energy of Figure 5 changed with simulation time. Figure 4(a) was the whole temperature changing process, and Figure 4(b) was the result after enlarging part of it.Figure 4
Change graph of temperature with simulation time after NPT dynamics finished.
(a)
System temperature fluctuated
(b)
Partial enlarged detail of temperature changeFigure 5
Change of system energy with simulation time.Figure 6
Change of (a) system density and (b) system length of side with time.
(a)
(b)Figure 7
Change of (a) system density and (b) system length of side with time at 0.08 GPa.
(a)
(b)It can be seen from Figure4(a) that the temperature gradually fluctuated around 298 K. It can be seen from Figure 4(b) that the fluctuation range was smaller (±10°C). The energy tended to be convergent gradually, suggesting that the system reached a balance after dynamic simulation.It was obtained from Figure6 that the system was continuously pressured, the density of system increased, the length of side decreased, and system was compressed 100 ps before the process of simulation. After 100 ps, the system tended to be balanced, and the parameters tended to be stable. When the pressure increased to 0.08 GPa, the simulation time was 200 ps, the system was compressed, the length of side decreased, and the density increased again. The curves of system density and the length of side changing with time are shown in Figure 7.It can be seen from Figure7 that the system was further compressed when the pressure of system increased to 0.08 GPa. It was found in Figure 7 that the system density was very close to the true density of colloid. At this moment, the system reached a balance again, and then there was a NPT balance lasting for 500 ps at 0.09 GPa. The increasing amplitude of density was limited. After verification, when the pressure was up to 0.09 GPa and increasing constantly, the changes of system density and length of side were very small. Finally, the system density and length of side at 0.09 GPa were used as the ultimate model to test the mechanical property of colloid. The system at the moment is shown in Figure 8. The density at this moment was 0.99 g/cm3.Figure 8
Ultimately obtained colloid system.(2) Stiffness Matrix and Mechanical Parameters of Colloid. With regard to the subsequently conducted operation and calculation processes to mechanical properties, the program will impose strains in the modes of 100000, 010000, 001000, 000100, 000010, and 000001 on the system model. These strains act on the directions of XX·YY·ZZ·YZ·ZX·XY of the mode. In each mode, four strain values ±0.001 and ±0.003 will be imposed. The default value of software is the maximum value of 0.003. Under this strain loading, the corresponding stress will be obtained, and the corresponding stiffness constant Cij will be obtained with this result to further obtain the elastic stiffness matrix of colloid mode, as shown in Table 1.Table 1
Stiffness matrix obtained from colloid model (GPa).
C
i
j
1
2
3
4
5
6
1
4.5089
2.8995
3.1170
0.1436
0.4446
−0.3247
2
2.8995
6.6093
3.5637
−0.0619
−0.4925
0.2330
3
3.1170
3.5637
4.5407
0.9075
0.6968
−0.1275
4
0.1436
−0.0619
0.9075
1.5682
0.5402
0.2406
5
0.4446
−0.4925
0.6968
0.5402
1.1374
−0.0765
6
−0.3247
0.2330
−0.1275
0.2406
−0.0765
0.0809It can be seen from Table1 that the values of C14, C15, C16, C24, C25, C26, C34, C35, C36, C45, C46, and C56 are close to 0, but, for the groups of C11, C22, C33; C12, C13, C23 and C44, C55, and C66, the values in each group are very close, suggesting that the colloid model of Figure 10 is close to isotropy. Suppose the material is isotropic, and the calculation of lame constant is shown in (11)λ=13C11+C22+C33-23C44+C55+C66=3.3620GPa,μ=13C44+C55+C66=0.9288GPa.After substitutingλ=3.3620 GPa and μ=0.9288 GPa into formula (9), the mechanical parameters of colloid were obtained below:(12)E=2.59GPa;ν=0.39;K=4.08GPa;G=0.93GPa.


### 3.2. Mechanical Property Simulation Analysis of Saturated Phenol

(1) Simulation Process and Result of Molecular Dynamics. The linear chain (n-C22H46) was selected to represent saturated phenol, and the degree of polymerization at 10 was selected to structure saturated phenol model. The model was optimized and placed into the cubic vacuum space where the length of three sides was 6 nm. The initial density was set at 0.1 g/cm3, and the initial saturated phenol model was obtained as shown in Figure 9.Figure 9
Initial model of saturated phenol at the density of 0.1.Figure 10
Change of (a) system density and (b) system length of side with time at 0.01 GPa.
(a)
(b)The true saturated density is 0.7944 g/cc, and thus the system should be compressed to increase the density and decrease the length of side. After optimizing the structure, NVT ensemble was conducted, and subsequently NPT ensemble was conducted stable again under 0.01 GPa and 298 K. The output temperature and energy curves tended to be stable and balanced with the changes of time after dynamics. The changes of density and the length of side of system can be seen in Figure10.It can be seen from Figure10 that the density was increasing continuously, but the trend density curve was not stable. After a NVT balance was conducted to the system to operate NPT ensemble, the pressure increased to 0.02 GPa, and simulation time increased to 400 ps, and the curve of system density changing with time is shown in Figure 11.Figure 11
Change of system density with time at 0.02 GPa.It can be seen from Figure11 that the density tended to be balanced, which was very close to the actual density at 0.7944, and the compression was stopped immediately. At this moment, the size of model was 29.01Å. The ultimately balanced saturated phenol model is shown in Figure 12.Figure 12
Ultimately saturated phenol model.(2) Stiffness Matrix and Mechanical Parameters of Saturated Phenol. The mechanical parameters of model in Figure 12 were calculated, and the stiffness matrix obtained for saturated phenol is shown in Table 2.Table 2
Stiffness matrix obtained for saturated phenol model.
C
i
j
1
2
3
4
5
6
1
1.7170
1.3630
1.8868
−0.2839
0.4201
−0.4593
2
1.3630
2.0288
2.0845
−0.3727
0.2747
0.0954
3
1.8868
2.0845
1.7931
−0.1096
0.2204
−0.0308
4
−0.2829
−0.3727
−0.1096
0.1437
−0.0936
0.1666
5
0.4201
0.2747
0.2204
−0.0936
0.1250
−0.1628
6
−0.4593
0.0954
−0.0308
0.1666
−0.1628
0.2710It can be seen from Table2 that the values of C14, C15, C16, C24, C25, C26, C34, C35, C36, C45, C46, and C56 are close to 0, but, for the three groups of C11, C22, and C33; C12, C13, and C23; and C44, C55, and C66, the values in each group are very close, suggesting that the saturated phenol model of Figure 12 is close to isotropy.Suppose the material is isotropic, and the calculation of lame constant can be seen in(13)λ=13C11+C22+C33-23C44+C55+C66=1.4864GPa,μ=13C44+C55+C66=0.1799GPa.After substitutingλ=1.4865 GPa and μ=0.1799 GPa into formula (9), the mechanical parameters of colloid were obtained below:(14)E=0.52GPa;ν=0.45;K=5.19GPa;G=0.18GPa.


### 3.3. Mechanical Property Simulation and Analysis of Asphaltene

(1) Simulation Process and Result of Molecular Dynamics. The C64H52S2 was used to represent asphaltene, and the degree of polymerization at 10 was used to structure asphaltene model. First, the degree of polymerization of asphaltene structure at 10 was placed into the vacuum space where the length of three sides was 6 nm, and the initial density was set at 0.1 g/cm3. The initial model of asphaltene obtained is shown in Figure 13.Figure 13
Initial model of asphaltene at the density of 0.1.Under the temperature of 298 K, the true density of asphaltene was approximately 0.89 g/cm3 that the system must be compressed to increase the density and decrease the length of side. After optimizing the structure, NVT ensemble was conducted. Subsequently, NPT was conducted to stable system under 0.06 GPa and 298 K. The changes of system density and length of side are shown in Figure 14.Figure 14
Change of (a) system density and (b) system length of side with time at 0.06 GPa.
(a)
System density
(b)
System length of sideIt can be seen from Figure14 that the density was increasing continuously, and it continued to increase. NVT and NPT were conducted to the stable system again, and the pressure increased to 0.08 GPa, and simulation time increased to 400 ps. The density curve was closed to actual density after the process. At this moment, the compression was stopped. The final output of the asphaltene model is shown in Figure 15. The system density was approximately 0.88 g/cm3, and the size of three sides of system was around 26.2Å.Figure 15
Ultimate model of asphaltene.(2) Stiffness Matrix and Mechanical Parameters of Asphaltene. After the calculation of mechanical parameters was conducted to the model in Figure 15, the stiffness matrix of asphaltene obtained is shown in Table 3.Table 3
Stiffness matrix obtained for asphaltene model.
C
i
j
1
2
3
4
5
6
1
7.8405
1.5765
1.9309
−0.3414
1.3945
−0.9280
2
1.5765
1.8645
0.6143
−0.2868
0.4326
−1.0948
3
1.9309
0.6143
2.5698
−0.2413
0.1393
−0.0903
4
−0.3414
−0.2868
−0.2413
1.1054
0.0539
0.1551
5
1.3945
0.4326
0.1393
0.0539
1.2063
−0.1860
6
−0.9280
−1.0948
−0.0903
0.1551
−0.1860
0.8541Suppose the material is isotropic, and the calculation of lame constant is shown in(15)λ=13C11+C22+C33-23C44+C55+C66=1.9811GPa,μ=13C44+C55+C66=1.0553GPa.After substitutingλ=1.9811 GPa and μ=1.0553 GPa into formula (9), the mechanical parameters of colloid were obtained below:(16)E=2.80GPa;ν=0.33;K=2.61GPa;G=1.05GPa.


## 3.1. Mechanical Property Simulation Analysis of Colloid

(1) Simulation Process and Result of Molecular Dynamics. The 1,7-dimethylnaphthalene was used as the repetitive unit to represent colloid in this paper. According to the suggestions of literature, macromolecular structure obtained at the polymerization degree in the range of 10~15 was employed [10, 11]. The character after calculation could meet the requirements such as research precision. The degree of polymerization selected for structuring the macromolecular chain of colloid was 15, and the degree of polymerization that was 15 to structure the long chain of colloid is shown in Figure 1.Figure 1
Long chain of colloid for the degree of polymerization at 15.Forcite module was used to analyze geometry optimization for the long chain of colloid. Figure2 is the long chain of colloid after geometry optimization.Figure 2
Long chain of colloid after being geometrically optimized.Afterwards, a cubic vacuum space in60Å×60Å×60Å was structured. The density was set at 0.1 g/cm3, and the initial model of colloid unit cell at the density of 0.1 obtained is shown in Figure 3.Figure 3
Initial model of colloid at the density of 0.1.The density of initial model is 0.1 g/cc, and the density was in the range of 0.99~1.1 g/cm3 at abnormal temperature; however, the low density system of Figure 3 is usually in a high energy state. Thus, structural optimization should be conducted to the system first to lower the potential energy, and then forcite can be used for optimization. When dynamic calculation was conducted, the temperature was 298 K (close to normal temperature), simulation time was 200 ps and pressure was 0.01 GPa, and the initial density of model at 0.1 would be compressed. In the simulation process, temperature, energy, system density, and the length of side of system were changing continually. The output images are shown from Figures 4–7, where the system energy of Figure 5 changed with simulation time. Figure 4(a) was the whole temperature changing process, and Figure 4(b) was the result after enlarging part of it.Figure 4
Change graph of temperature with simulation time after NPT dynamics finished.
(a)
System temperature fluctuated
(b)
Partial enlarged detail of temperature changeFigure 5
Change of system energy with simulation time.Figure 6
Change of (a) system density and (b) system length of side with time.
(a)
(b)Figure 7
Change of (a) system density and (b) system length of side with time at 0.08 GPa.
(a)
(b)It can be seen from Figure4(a) that the temperature gradually fluctuated around 298 K. It can be seen from Figure 4(b) that the fluctuation range was smaller (±10°C). The energy tended to be convergent gradually, suggesting that the system reached a balance after dynamic simulation.It was obtained from Figure6 that the system was continuously pressured, the density of system increased, the length of side decreased, and system was compressed 100 ps before the process of simulation. After 100 ps, the system tended to be balanced, and the parameters tended to be stable. When the pressure increased to 0.08 GPa, the simulation time was 200 ps, the system was compressed, the length of side decreased, and the density increased again. The curves of system density and the length of side changing with time are shown in Figure 7.It can be seen from Figure7 that the system was further compressed when the pressure of system increased to 0.08 GPa. It was found in Figure 7 that the system density was very close to the true density of colloid. At this moment, the system reached a balance again, and then there was a NPT balance lasting for 500 ps at 0.09 GPa. The increasing amplitude of density was limited. After verification, when the pressure was up to 0.09 GPa and increasing constantly, the changes of system density and length of side were very small. Finally, the system density and length of side at 0.09 GPa were used as the ultimate model to test the mechanical property of colloid. The system at the moment is shown in Figure 8. The density at this moment was 0.99 g/cm3.Figure 8
Ultimately obtained colloid system.(2) Stiffness Matrix and Mechanical Parameters of Colloid. With regard to the subsequently conducted operation and calculation processes to mechanical properties, the program will impose strains in the modes of 100000, 010000, 001000, 000100, 000010, and 000001 on the system model. These strains act on the directions of XX·YY·ZZ·YZ·ZX·XY of the mode. In each mode, four strain values ±0.001 and ±0.003 will be imposed. The default value of software is the maximum value of 0.003. Under this strain loading, the corresponding stress will be obtained, and the corresponding stiffness constant Cij will be obtained with this result to further obtain the elastic stiffness matrix of colloid mode, as shown in Table 1.Table 1
Stiffness matrix obtained from colloid model (GPa).
C
i
j
1
2
3
4
5
6
1
4.5089
2.8995
3.1170
0.1436
0.4446
−0.3247
2
2.8995
6.6093
3.5637
−0.0619
−0.4925
0.2330
3
3.1170
3.5637
4.5407
0.9075
0.6968
−0.1275
4
0.1436
−0.0619
0.9075
1.5682
0.5402
0.2406
5
0.4446
−0.4925
0.6968
0.5402
1.1374
−0.0765
6
−0.3247
0.2330
−0.1275
0.2406
−0.0765
0.0809It can be seen from Table1 that the values of C14, C15, C16, C24, C25, C26, C34, C35, C36, C45, C46, and C56 are close to 0, but, for the groups of C11, C22, C33; C12, C13, C23 and C44, C55, and C66, the values in each group are very close, suggesting that the colloid model of Figure 10 is close to isotropy. Suppose the material is isotropic, and the calculation of lame constant is shown in (11)λ=13C11+C22+C33-23C44+C55+C66=3.3620GPa,μ=13C44+C55+C66=0.9288GPa.After substitutingλ=3.3620 GPa and μ=0.9288 GPa into formula (9), the mechanical parameters of colloid were obtained below:(12)E=2.59GPa;ν=0.39;K=4.08GPa;G=0.93GPa.


## 3.2. Mechanical Property Simulation Analysis of Saturated Phenol

(1) Simulation Process and Result of Molecular Dynamics. The linear chain (n-C22H46) was selected to represent saturated phenol, and the degree of polymerization at 10 was selected to structure saturated phenol model. The model was optimized and placed into the cubic vacuum space where the length of three sides was 6 nm. The initial density was set at 0.1 g/cm3, and the initial saturated phenol model was obtained as shown in Figure 9.Figure 9
Initial model of saturated phenol at the density of 0.1.Figure 10
Change of (a) system density and (b) system length of side with time at 0.01 GPa.
(a)
(b)The true saturated density is 0.7944 g/cc, and thus the system should be compressed to increase the density and decrease the length of side. After optimizing the structure, NVT ensemble was conducted, and subsequently NPT ensemble was conducted stable again under 0.01 GPa and 298 K. The output temperature and energy curves tended to be stable and balanced with the changes of time after dynamics. The changes of density and the length of side of system can be seen in Figure10.It can be seen from Figure10 that the density was increasing continuously, but the trend density curve was not stable. After a NVT balance was conducted to the system to operate NPT ensemble, the pressure increased to 0.02 GPa, and simulation time increased to 400 ps, and the curve of system density changing with time is shown in Figure 11.Figure 11
Change of system density with time at 0.02 GPa.It can be seen from Figure11 that the density tended to be balanced, which was very close to the actual density at 0.7944, and the compression was stopped immediately. At this moment, the size of model was 29.01Å. The ultimately balanced saturated phenol model is shown in Figure 12.Figure 12
Ultimately saturated phenol model.(2) Stiffness Matrix and Mechanical Parameters of Saturated Phenol. The mechanical parameters of model in Figure 12 were calculated, and the stiffness matrix obtained for saturated phenol is shown in Table 2.Table 2
Stiffness matrix obtained for saturated phenol model.
C
i
j
1
2
3
4
5
6
1
1.7170
1.3630
1.8868
−0.2839
0.4201
−0.4593
2
1.3630
2.0288
2.0845
−0.3727
0.2747
0.0954
3
1.8868
2.0845
1.7931
−0.1096
0.2204
−0.0308
4
−0.2829
−0.3727
−0.1096
0.1437
−0.0936
0.1666
5
0.4201
0.2747
0.2204
−0.0936
0.1250
−0.1628
6
−0.4593
0.0954
−0.0308
0.1666
−0.1628
0.2710It can be seen from Table2 that the values of C14, C15, C16, C24, C25, C26, C34, C35, C36, C45, C46, and C56 are close to 0, but, for the three groups of C11, C22, and C33; C12, C13, and C23; and C44, C55, and C66, the values in each group are very close, suggesting that the saturated phenol model of Figure 12 is close to isotropy.Suppose the material is isotropic, and the calculation of lame constant can be seen in(13)λ=13C11+C22+C33-23C44+C55+C66=1.4864GPa,μ=13C44+C55+C66=0.1799GPa.After substitutingλ=1.4865 GPa and μ=0.1799 GPa into formula (9), the mechanical parameters of colloid were obtained below:(14)E=0.52GPa;ν=0.45;K=5.19GPa;G=0.18GPa.


## 3.3. Mechanical Property Simulation and Analysis of Asphaltene

(1) Simulation Process and Result of Molecular Dynamics. The C64H52S2 was used to represent asphaltene, and the degree of polymerization at 10 was used to structure asphaltene model. First, the degree of polymerization of asphaltene structure at 10 was placed into the vacuum space where the length of three sides was 6 nm, and the initial density was set at 0.1 g/cm3. The initial model of asphaltene obtained is shown in Figure 13.Figure 13
Initial model of asphaltene at the density of 0.1.Under the temperature of 298 K, the true density of asphaltene was approximately 0.89 g/cm3 that the system must be compressed to increase the density and decrease the length of side. After optimizing the structure, NVT ensemble was conducted. Subsequently, NPT was conducted to stable system under 0.06 GPa and 298 K. The changes of system density and length of side are shown in Figure 14.Figure 14
Change of (a) system density and (b) system length of side with time at 0.06 GPa.
(a)
System density
(b)
System length of sideIt can be seen from Figure14 that the density was increasing continuously, and it continued to increase. NVT and NPT were conducted to the stable system again, and the pressure increased to 0.08 GPa, and simulation time increased to 400 ps. The density curve was closed to actual density after the process. At this moment, the compression was stopped. The final output of the asphaltene model is shown in Figure 15. The system density was approximately 0.88 g/cm3, and the size of three sides of system was around 26.2Å.Figure 15
Ultimate model of asphaltene.(2) Stiffness Matrix and Mechanical Parameters of Asphaltene. After the calculation of mechanical parameters was conducted to the model in Figure 15, the stiffness matrix of asphaltene obtained is shown in Table 3.Table 3
Stiffness matrix obtained for asphaltene model.
C
i
j
1
2
3
4
5
6
1
7.8405
1.5765
1.9309
−0.3414
1.3945
−0.9280
2
1.5765
1.8645
0.6143
−0.2868
0.4326
−1.0948
3
1.9309
0.6143
2.5698
−0.2413
0.1393
−0.0903
4
−0.3414
−0.2868
−0.2413
1.1054
0.0539
0.1551
5
1.3945
0.4326
0.1393
0.0539
1.2063
−0.1860
6
−0.9280
−1.0948
−0.0903
0.1551
−0.1860
0.8541Suppose the material is isotropic, and the calculation of lame constant is shown in(15)λ=13C11+C22+C33-23C44+C55+C66=1.9811GPa,μ=13C44+C55+C66=1.0553GPa.After substitutingλ=1.9811 GPa and μ=1.0553 GPa into formula (9), the mechanical parameters of colloid were obtained below:(16)E=2.80GPa;ν=0.33;K=2.61GPa;G=1.05GPa.


## 4. Conclusions

The calculation of mechanical properties was conducted to the groups of colloid, saturated phenol, and asphaltene, respectively, and the parameters of stiffness and modulus of all the materials were obtained. Comparing theE value of Young modulus of these three polymers, the maximum modulus value of asphaltene was 2.80 GPa, the modulus value of colloid was secondary, and the minimum modulus value of saturated phenol was 0.52 GPa. In asphalt, asphaltene is similar to solid, and Young’s modulus is largest, followed by colloid. Saturated phenol is similar to liquid, of which the modulus is far less than asphaltene. In addition, comparing the values of Poisson’s ratio, it can be found that the maximum Poisson’s ratio of saturated phenol is 0.45, which is close to 0.5, further suggesting that saturated phenol is similar to liquid that it is incompressible. However, more molecular morphologies require further investigation.



---

*Source: 1531632-2017-10-26.xml*