Next Article in Journal
Fracture Resistance of Fiber-Reinforced Composite Restorations: A Systematic Review and Meta-Analysis
Previous Article in Journal
Effect of Thermal Aging on Mechanical Properties and Morphology of GF/PBT Composites
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Molecular Dynamics Simulation Study on Enhancement of Mechanical and Tribological Properties of Nitrile—Butadiene Rubber with Varied Contents of Acrylonitrile

1
School of Mechanical Engineering, Shenyang University of Technology, Shenyang 110870, China
2
School of Chemical Equipment, Shenyang University of Technology, Liaoyang 111003, China
3
Ningbo Institute of Technology, Beihang University, Ningbo 315800, China
*
Authors to whom correspondence should be addressed.
Polymers 2023, 15(18), 3799; https://doi.org/10.3390/polym15183799
Submission received: 25 August 2023 / Revised: 8 September 2023 / Accepted: 8 September 2023 / Published: 18 September 2023
(This article belongs to the Section Polymer Analysis and Characterization)

Abstract

:
The molecular models of nitrile–butadiene rubber (NBR) with varied contents of acrylonitrile (ACN) were developed and investigated to provide an understanding of the enhancement mechanisms of ACN. The investigation was conducted using molecular dynamics (MD) simulations to calculate and predict the mechanical and tribological properties of NBR through the constant strain method and the shearing model. The MD simulation results showed that the mechanical properties of NBR showed an increasing trend until the content of ACN reached 40%. The mechanism to enhance the strength of the rubber by ACN was investigated and analyzed by assessing the binding energy, radius of gyration, mean square displacement, and free volume. The abrasion rate (AR) of NBR was calculated using Fe-NBR-Fe models during the friction processes. The wear results of atomistic simulations indicated that the NBR with 40% ACN content had the best tribological properties due to the synergy among appropriate polarity, rigidity, and chain length of the NBR molecules. In addition, the random forest regression model of predicted AR, based on the dataset of feature parameters extracted by the MD models, was developed to obtain the variable importance for identifying the highly correlated parameters of AR. The torsion–bend–bend energy was obtained and used to successfully predict the AR trend on the new NBR models with other acrylonitrile contents.

Graphical Abstract

1. Introduction

The wear phenomenon of polymer materials has received considerable attention as it directly affects the reliability, stability, and safety of materials in industrial applications [1,2,3,4,5,6]. As one of the most attractive matrixes in polymers, nitrile–butadiene rubber (NBR) has exhibited outstanding wear resistance, chemical resistance, and a wide range of operating environments [7]. On the other hand, NBR is the most widely used polar rubber in petroleum industries and mechanical seals due to its moderate cost, good processability, and easy modification [8,9]. Therefore, the mechanical and tribological properties of NBR have been investigated by many experimental researchers [10,11,12,13,14]. The most common method for improving wear resistance of NBR is based on the increasing acrylonitrile (ACN) content to raise the polarity and rigidity of polymer, thereby enhancing material mechanical and tribological properties [15,16]. It has to be noted that although the above experimental studies have acknowledged the best performance of 40% ACN content of NBR, to the authors’ knowledge, few theoretical or molecular simulation studies have been identified. In order to guide the NBR modification in tribology applications, studies are also significant and useful in determining the enhancement mechanism of NBR as the optimizing direction from an atomic level.
As an effective tool for studying materials [16,17,18,19,20], molecular dynamics (MD) simulations can calculate the properties of the polymer model at the macro level and provide details of molecular interactions to explain the mechanisms of action on the atomic scale. Several themes that emerged from investigators have reported on the physical properties of polymer composites formed by MD simulations, such as the glass-transition temperature, mechanical and tribological properties, tribo-chemical properties, thermal properties, and damping properties [4,21,22,23,24,25]. Data from Li [26,27] suggested that the tribological properties of polymer composites under the different variables of normal loading and sliding velocity could be enhanced by the introduction of carbon nanotubes based on MD simulations, where physical interpretations were provided from microscopic information. It can be indicated from the above studies and other research [28,29,30] that the destruction of materials is more from the transfer of external kinetic energy into the potential energy of the polymer by the microscopic perspective of simulation analysis. This means that improving the interaction energy of molecular chains of a polymer is more important for reducing the potential energy of the polymer.
In the MD simulations, the polymer properties were directly or indirectly determined according to the system energy by parameterizing the structural features of polymers [31]. The potential energy of a system can be expressed as a sum of valence (or bond), non-bond interactions, and cross term: Etotal = Evalence + Enon-bond + Ecrossterm [32]. In addition, these energies can be further subdivided into a dozen specific parameters, such as the bond energy, the angle energy, the Van der Waals energy, and the bend–torsion–bend energy, etc. The importance of these specific parameters on influencing the properties of polymers was often neglected in previous MD theoretical studies due to the large quantity of parameters. These energies are likely to be key correlation factors in the study and prediction of material properties [33,34]. The random forest technique is able to simultaneously analyze hundreds of feature variables [35,36,37,38,39] and produce the ranking of feature importance based on the predictive model of machine learning. A study by Betina [40] proposed a basic method that obtained information on the most important variables from the random forest model, and it was shown to be effective for the distinction of the variables that best contributed to the discrimination of oils.
Considering all the aforementioned information, the mechanical and tribological properties of NBR were calculated and analyzed by MD simulations. Thereinto, according to the universal components of NBR [41], the NBR models considered ACN contents of 20% to 40% with incremental contents of 10%. To determine the atomic mechanisms of enhancements of ACN, the binding energy, radius of gyration, mean square displacement (MSD), and free volume were calculated from the NBR system. In addition, the results and energy values were collected from the MD model for identifying the high correlation parameters of wear performance in different ACN contents of NBR.

2. Models and Methods

2.1. MD Simulations

MD simulations were performed with Materials Studio 8.0 software. Molecular modeling and dynamic calculation of 4 NBR models were carried out using Amorphous Cell and Forcite Modules. The COMPASS force field [32], which is an ab initio force field that can accurately predict the structural characteristics of a wide range of molecules and polymers, was implemented. These models were constructed as follows: A single chain consisting of 20 repeat units was built with acrylonitrile ( C 3 H 3 N ) and butadiene ( C H 2 = C H C H = C H 2 ) in 1:4, 3:7, 2:3, and 1:1 allocation ratios [10,11,12,13,14]. The chain length was measured from the atomic distances between the two ends [18]. Each chain was individually packed into an empty cell with a size of 4.5 × 3.0 × 3.0 nm3 and periodic boundary conditions until the cell density reached 0.97 g/cm3 [22]. The molecules in a cell were built with a Monte Carlo style by minimizing close contacts among atoms whilst keeping a realistic distribution of torsion angles. The information of the NBR models is listed in Table 1.
Since the internal energy was too high in the original NBR models, geometry optimization was applied, and the conjugate gradient algorithm [42] was adopted with an energy convergence of 10−5 kcal/mol. The equilibrated structure of the models was obtained further using the NVT (constant volume, constant temperature) and NPT (constant pressure, constant temperature) ensemble. A 600 ps NVT was conducted under a temperature of 298 K. On the last snapshot of the NVT trajectory, an NPT ensemble was carried out for 600 ps at a temperature and pressure of 298 K and 101 KPa. These ensembles were simulated using the Andersen thermostat [43] and Berendsen barostat [44] with time steps of 1 fs. Meanwhile, the Ewald and atom-based summation method was applied to calculate the electrostatic and van der Waals (vdW) interactions with an accuracy of 0.001 kcal/mol and a cutoff distance of 1.25 nm, respectively.
To investigate the mechanical properties of the models, the constant strain method was applied to calculate the Young’s modulus (E), bulk modulus (B), and shear modulus (G). Next, another 200 ps NVT was calculated to obtain the trajectory, in which 4 strains and a maximum strain (ε) amplitude of 0.003 were applied to compute the moduli. The elastic constant matrix was estimated by a series of finite difference approximations. In each Cartesian direction, E = σ/ε was applied, according to the virial stress definition [45]. By combining the Voigt [46] and Reuss [47] methods with Hill’s definition [48], B and G were estimated by the 6 × 6 stiffness and compliance matrices accordingly.
To predict the tribological properties of NBR, Fe-NBR-Fe models were proposed and equilibrated by geometry optimization, and the abrasion rate (AR) of NBR was calculated. Fe atom layers with dimensions 4.5 × 3.0 × 0.86 nm3 were developed and used as the top and bottom slip planes (the intermediate layer from the equilibrated structure of the models). A 5-cycle annealing simulation with an initial temperature of 150 K and a mid-cycle temperature of 350 K was then conducted by a 200 ps NVT. Finally, the confined shear was realized by applying a sliding speed of 0.1 Å/ps for 500 ps on the upper Fe layer (up-Fe-layer). After the end of the shearing process, AR = Mleave/Mtotal was applied, where Mleave and Mtotal represent the number of atoms that leave the NBR matrices during the friction process and the total number of atoms in the original NBR matrices, respectively. By simulating the friction process, FN and Ff (the normal and friction forces on the up-Fe-layer) were obtained, and the AR was measured at the last snapshot of the shearing trajectory.

2.2. Datasets

In order to identify the high correlation parameters of AR, the forest model of predicting AR needed to be constructed by feature variables of NBR. In this study, the feature parameters of NBR as the dataset variables, including 27 input variables and AR output variable, were collected from the MD model. The dataset consisted of 5 parameter modules: structural parameters, mechanical properties, environmental variables, friction parameters, and energy parameters. These parameters were derived from the Fe-NBR-Fe model before the shear process, except the AR and mechanical properties, as clearly demonstrated in Figure 1 and Table 2.
Thereinto, the specific parameters of energies were obtained from the empirical potential function [31,49,50] of the force field and are summarized by the expressions below.
H a m i l t o n i a n = T o t a l   k i n e t i c   e n e r g y + T o t a l   p o t e n t i a l   e n e r g y
E t o t a l p o t e n t i a l = E v a l e n c e + E c r o s s t e r m + E n o n b o n d
E v a l e n c e = E b o n d + E a n g l e + E t o r s i o n + E o u t o f p l a n e + E U r e y B r a d l e y
E c r o s s t e r m = E B e n d B e n d + E B e n d T o r s i o n B e n d + E T o r s i o n B e n d B e n d   + E T o r s i o n S t r e t c h + E S t r e t c h B e n d S t r e t c h + E S t r e t c h T o r s i o n S t r e t c h   + E S t r e t c h S t r e t c h
E n o n b o n d = E v a n   d e r   W a a l s + E C o u l o m b + E h y d r o g e n   b o n d
To enrich the dataset, the normal loading and sliding velocity were selected as the feature parameters of the environmental variables to additionally simulate the friction process on each NBR model. The normal loadings were divided into 11 levels, ranging from 0 to 500 kcal/mol/Å, at a sliding speed of 1.1 Å/ps. The sliding velocities were classified into 2 levels with gradients of 0.01 Å/ps (0.03–0.09 Å) and 0.1 Å/ps (0.1–1.0 Å). In summary, 112 samples of the friction process were obtained for forming the original dataset.
It is worth noting that the energy parameters were equivalently designed as the structure descriptors of polymer through MD simulations in this study. The structural parameters of the dataset were also added to explore whether the structural descriptors can be further simplified by atomic information of models.

2.3. Identifying the Correlation Parameters of AR

Random forest is one of the most popular supervised machine learning approaches, and it can be applied to solve a wide range of prediction problems [37] with a high accuracy compared with other currently used algorithms. Aside from being able to deal with high-dimensional variables without deleting variables, the method has been highly recognized for its inherent procedure of producing the variable importance (VI) [38,39]. From the VI of the prediction model, the importance variables of the AR effect can be initially screened for identifying the correlation parameters. Hence, the random forest regression (RFR) model of predicted AR was developed to obtain the VI of feature parameters. According to the original dataset, different decision trees were constructed to develop the RFR model, as shown schematically in Figure 2.
First, a new dataset for each tree was established using the bootstrap strategy. Next, the new dataset was split into a training dataset (approximately 70%, in-bag) and a test dataset (approximately 30%, out-of-bag). A 4-fold cross validation [51] was then applied to avoid overfitting in the bag step. The training dataset was divided into 3 training subsets and a validation subset. To separate the feature variables from the 3 training subsets, predictive variables were identified within the random selection by the mean squared error (MSE) criterion and were used to set the nodes during regression tree construction [52]. Additionally, hyperparameters of the separation were optimized by the validation subset. The final model of the decision tree was formed by 4 separate results.
The whole process of building the RFR model involved iteratively repeating the growth process of one tree until the total number of trees reached 100 (ntree = 100). The rules and judgment values from all decision trees were then integrated and averaged. According to the contribution proportions of the feature variables, the VI was calculated.
The correlation variable of AR from feature parameters was identified by analyzing the top parameters in VI. The specific analysis is detailed in the RFR results of Section 3.3.

3. Results and Discussions

3.1. Mechanical Properties

The mechanical properties of NBR were investigated by varying the ACN content. According to the 2.1 MD simulations, the mechanical modulus of NBR was simulated 3~5 times for each model and reported the property values from the average of 3 times simulations. As illustrated in Figure 3 and Table 3, E, B, and G values of these models of 3.05 to 3.80, 2.95 to 3.39, and 1.25 to 1.42 GPa, respectively, were obtained with increasing values of ACN. Compared with the monograph [32], the research also found that increasing the ACN content reinforced the mechanical properties attributed to the polarity and rigidity of NBR. The polarity atomic nitrogen with an increase of the ACN group enhanced rubber polarity to increase the interactions of NBR molecular chains. According to the conformation of polymer chains, the flexibility of NBR molecules with an increase in ACN were reduced, owing to reduced double bonds and increased side groups. Meanwhile, the reduced chain length and increased molecular rigidity caused a decrease of entanglement between the NBR molecular chains.
In this study, these moduli approximately showed an increasing trend until the content of ACN reached 40%. In addition to the Bulk modulus, the Young’s and shear moduli of ACN 50 were only slightly increased and decreased compared to ACN 40. Overall, ACN enhanced NBR mechanical properties and showed a turning point at 40% content. A possible explanation for this is that 40% ACN represents a turning point for the dominant role of the reinforcing effect in NBR. The increasing polarity and rigidity of NBR with an increased ACN content enhanced the interactions of NBR molecular chains until the strong rigidity of NBR limited the movement and entanglement of chains after this point.
In order to determine the atomic mechanisms responsible for the enhancement in the mechanical properties, the interaction of NBR molecular chains was calculated during the NPT equilibrium processes. A complete molecular chain from each model is marked as an NBRL in Figure 4 (shown as a red chain in Figure 1 in some cells, too). The binding energy (Ubind) between the NBRL and the NBR matrix was calculated as shown in Figure 4 using Equation (6):
U b i n d = U i n t e r = U t o t a l U N B R L U N B R
where Utotal represents the total energy of the NBR, UNBRL identifies the energy of the NBRL, and UNBR represents the energy of the NBR, except for the NBRL. In order to reduce the errors of different NBRL, the average value of the binding energy of NBRL was calculated and repeated 10 times to further average as shown in Figure 4. The Ubind increased from 195.70 to 206.52 kcal/mol as the ACN content increased, indicating that having more ACN enhanced the interaction of NBR molecular chains by the increase in polarity and thus improved the mechanical properties of NBR. It was also observed from the structure of the NBR chain that the degree of coiling of molecules with an increase in ACN was stronger. The cause of this phenomenon can be explained that the double bond number of models was decreasing to reduce the bond angle value of the NBR backbone, and the polarity of the models was increasing to enhance intramolecular interaction force.
In the NPT equilibrium processes, the radius of gyration of NBR molecular chains of models was monitored as a parameter to assess the flexibility of NBR and to show the influence of ACN contents on the movement of NBR molecular chains. As shown in Figure 5, the radius of gyration of NBR molecular chains with increasing values of ACN was reduced, indicating that NBR molecular chains with an increase in ACN had a smaller radius of gyration, lower flexibility, and more rigidity. This comparison showed that the increasing of the ACN group restricted the flexibility and enhanced the rigidity of NBR molecular chains, thereby improving the shear resistance and further the mechanical properties of NBR. Meanwhile, the standard deviation of the radius of gyration of models by an increased ACN content declined, which showed that the ACN group limited the range of motion of NBR molecular chains during the NPT equilibrium process. Especially, the standard deviation of ACN50 value was 0.03 Å, showing decreases of 30%, 57%, and 40% compared with others, respectively. The extremely high rigidity of the ACN50 model limited the entanglement of NBR molecular chains, thus little improving the mechanical properties of ACN50 than ACN40 of NBR.
To verify the above hypothesis of a 40% turning point, the MSD of the NBR molecular chains, which plays a vital role in determining the movement of the polymer chains [53], was calculated during the NPT equilibrium process and is recorded in Figure 6. The average MSD values of models with increases in ACN were 3.51, 3.01, 2.66, and 2.05 Å2, showing decreases of 14.26%, 24.30%, and 41.61% compared with the ACN20 model, respectively. In comparison with the other models, the flatter slope of the ACN50 curve indicates that the motion of NBR molecular chains was greatly limited due to its excessive polarity and rigidity, and thus the mechanical properties of ACN50 barely improved relative to those of ACN40. This restriction largely dampened the entanglement between molecules, affecting the enhancement of strength during the equilibrium processes.
In addition, to further assess the interactions and movement of NBR molecular chains, the free volume was measured based on the equilibrated structures of models, as recorded in Figure 7. The field volume was calculated as Field Volume = Enclosed Volume + Free Volume. The Connolly surface was calculated when a probe molecule with a radius (Rp) of 1.0 Å rolled over the vdW surface, and the enclosed volume and free volume were determined as the volumes on the side of the Connolly surface with and without atoms, respectively. The trends of the field volume (the blue and gray bar) and enclosed volume (the gray bar) slightly decreased as the ACN content increased, which indicates that the interactions and movement of NBR molecular chains were enhanced and limited with more branches of the ACN group, leading to NBR molecular chains with a higher polarity and stronger rigidity. These results agree with the above-mentioned analysis of mechanical properties.
It is interesting to mention that ACN30 had the peak value for B, Ubind, and free volume for all models. The standard deviation of the radius of gyration on ACN30 (0.07 Å) was more than on ACN20 (0.06 Å), which meant the range of motion of NBR molecular chains was larger on ACN30, in agreement with the result of free volume. Under the same increment on polarity of NBR, the distance of intermolecular interactions increased, which led to a decrease of B and Ubind, rather than a rise similar to others. Compared with the E and G, the relative position of intermolecular could be considered as a key factor on the mechanism of B. The polarity and rigidity of NBR were critical factors to determine the performance on the E and G.
In addition, there are many different initial polymer configurations in the same ACN content due to the complex polymer configurations. New models of 30% ACN content were additionally built from other monomer arrangements to validate the ACN30 finding. From Table 4, the 1211 (the original ACN30 model) and 2121 items contained 2 ACN consecutive units, and the 2130 and 3111 items contained 3 ACN consecutive units. According to the above analytical methods, the Ubind, and free volume of these models were calculated and shown in Table 4 and Figure 8.
As shown in Figure 8, the Ubind values of 30% ACN content models were lower than the others, and the free volume of these models were larger than the others. The results showed that the peak value of ACN30 was not by accident, and the analysis of key factors of the mechanism was meaningful. Moreover, ACN30 models starting from the different structural configurations could generate slightly different values as shown in Figure 8, but these results did not affect the analysis of other systems in the study.
The next section further analyses the tribological properties of NBR by observing the friction process of models.

3.2. Tribological Properties

To investigate the performance of the friction and wear of the NBR matrix with an increase in the ACN content, snapshots of NBR models during the shearing process were observed and are recorded in Figure 9. The phenomenon of adsorption occurred at the interface between the Fe and NBR matrices owing to the significance of the attractive vdW forces and leading to the wear behavior. AR values of 36.76%, 30.94%, 20.90%, and 24.33% were obtained for the models as the ACN content increased. These values represent decreases of approximately 16%, 43%, and 34% as compared with ACN20 (36.76%). It can be observed from Figure 9 that with an increase in the ACN content, the NBR molecular chains were surrounded by more atoms from the NBR matrices, and the relative slippage appeared more clearly between the friction interfaces. As also indicated, the tribological properties of NBR can be enhanced due to the increasing polarity and rigidity with an increase in the ACN content, leading to greater interactions, greater stiffness, and adsorption reductions of NBR until an ACN content of 40% is reached.
To further verify the above discussions, the depth variations in the up-Fe-layer of the models were recorded to observe the deformation of the NBR when the external force of the up-Fe-layer was removed. The rising amplitudes of the up-Fe-layer with increases in the ACN content were 0.71, 0.25, 0.22, and 0.34 Å, respectively. Compared with the value of ACN20, the others were significantly lower, indicating that greater interactions and stiffness were achieved with an increase in the ACN content. The smallest value of ACN40 also represented the best tribological property of NBR, owing to stronger interactions and tighter entanglement, as shown for ACN40. The best result was in good agreement with the AR of ACN40, which was also consistent with the 40% turning point of the mechanical properties of NBR.
Furthermore, the properties of the NBR matrix with ACN50 declined, which could have been due to the state of NBR molecular chains of ACN50, as shown in Figure 9. The NBR molecular chains with the shortest length were subjected to limited movement and less entanglement, causing a smaller interaction area between the molecule contacts. Thus, the AR value of ACN50 increased during the shearing process, and the tribological properties of NBR had no improvement with increases in the ACN content. This meant that the mechanism to enhance NBR by the acrylonitrile was the positive effect of the polarity combined with the negative effect of the rigidity to affect the tribological properties of NBR after 40% ACN content.
Overall, the results of the analysis are consistent with experimental observations [9,52]. The results prove that the variable parameters of these models obtained by performing MD simulations are reliable, and they will become fundamental to the RFR dataset coming from Section 2.2.

3.3. RFR Model Results

To avoid the randomness of the RFR model, the datasets were trained multiple times by the RFR algorithm and an importance ranking of the variables was produced. The feature variables of FN and Ff were selected as the key factor affecting AR to screen RFR models. Table 5 lists eight RFR models and shows the VI names and values (the proportion of VI) from the top five for obtaining the highly correlated variables of AR from dataset. Not surprisingly, apart from the FN and Ff, the temperature was also related to AR. The bend–torsion–bend (BTB), torsion–bend–bend (TBB), torsion–stretch (TS), and stretch–torsion–stretch (STS) energies as the structure descriptors are also shown in Table 5. Compared with the AR value of 112 samples, the value of TBB energy was highly correlated with the trend of AR value.
To analyze all parameters of the dataset, the VI of the RFR Model 2 was determined and is shown in Figure 10, which was distinguished according to the parameter modules (from Table 2). The energy parameters, environmental variables, and friction variables occupied 98% of all modules in percentages of 73%, 13.8%, and 11.2%, respectively. As the structure descriptors, the energy parameters mainly contributed to the prediction ability because they have greater influences than external conditions on polymers [54].
In addition, for example, the ACN20 model at a speed of 0.1 Å/ps, the total potential, valence, cross-terms and non-bond energy values were −204,951.00, −1029.05, −672.13, and −202,235.40 kcal/mol, respectively. The total cross-terms energy (59.70% of VI) contributed less energy value than others, leading to its significance often being ignored in previous studies of the energy effects using MD simulations. The cross-terms energy demonstrates the importance of accounting for structure, energy, and dynamics, as shown by JON R. MAPLE [55], and was further used to analysis and predict the AR of NBR in this work.
Figure 10 also summarizes the top 10 VI values and shows that the BTB and TBB energies were superior to others. To explore the effects of these parameters on the wear performance of NBR, the BTB and TBB energy values were obtained from the end of each equilibrium process, including the structure of the packing after geometry optimization, i.e., the structure of the NVT ensemble, the structure of the NPT ensemble, and the structure of the annealing simulation. It was found that the BTB energies obtained from the structures of the packing and annealing processes were highly correlated with the ACN contents of the models, as can be seen in Figure 11A.
The most striking result emerged between the TBB energy and the AR of NBR, as shown in Figure 11B. The same trend was observed from the structure of the packing and annealing processes together with the AR. This may mean that the AR trend of NBR could be estimated by contrasting the TBB energy of NBR at the early stage of the model’s construction. A bold hypothesis was proposed that the value of TBB energy, from the NBR model of the packing after the geometry optimization, could directly predict AR trends for NBR model of other ACN contents. This conjecture was tested and verified by Section 3.4.
It is also worth clarifying that the TBB energy of ACN40 compared with other models had the minimum value at the end structure of any process, which is in agreement with the turning point of a 40% ACN content. This result could be explained by MD simulations showing that the lower TBB energy of NBR could prevent the relative movement of molecules, thereby maintaining the stability of the mechanical properties.

3.4. Test and Verification of the Hypothesis

In order to verify the agreement in the trend between the TBB energy and AR of models, new models were packed using the other ACN contents, including proportions of 15%, 25%, 35%, 45%, and 55%. The TBB energy of all models was summarized from the packing structure after geometry optimization, as shown by the pink line in Figure 12. ACN25, ACN15, and ACN55 had, respectively, the maximum values of TBB energy and extrapolated data points for all models, as shown in Figure 12. Therefore, they were selected as the validation models simulated by the MD simulations presented in Section 2.1.
As previously predicted, the AR values of validation models agreed with the trend of the TBB energy, as shown the green bars in Figure 12. This meant that the TBB energy was found as the highly correlated parameter of AR by screening VI from RFR model, which is meaningful and used to roughly predict the AR trends. From the 15 to 25% ACN models’ snapshots, the AR values were increased due to the more atomic adsorption around the interface between metal and polymer materials. This suggests that the increasing polarity in 15–25% ACN contributes more to the vdW force at the interface of the Fe and NBR matrices. In the ACN55 model, NBR emerged with a brittle fracture leading to the higher AR value, owing to the decreases in molecular chain length and an increase in rigidity, as discussed in Section 3.2. These results show that TBB energy can be used as a parameter to predict the trend of AR and the wear performance of NBR with different ratios.
Combined with Figure 9, the lower AR values of NBR had the following characteristics: the stronger interaction force between molecular chains and the lower vdW force at the interface of the Fe and polymer matrices. In particular, the wear results of ACN25 and ACN55 showed that the up-Fe-layer adsorbed more matrices, leading to AR values that were more than others. Therefore, the modified NBR can be considered to enhance the interaction force of molecules and reduce the vdW force between the matrix and Fe layer.
In this study, although no more NBR models were validated, the method of identifying VI by introducing the RFR model still can be used to assist MD analysis in obtaining additional information quickly and efficiently. It provides new solutions for analyzing the mechanism, predicting the performance, and guiding the design in the development of new polymer materials. More precise prediction and analysis using important variables will be deeply investigated in future.
It is worth noting that the regression analysis of RFR model, for treating the use and prioritization of multivariate parameters and their impact on quantitative prediction performance, is widely applicable and not limited to the dimension and variety of data. Developing multiple RFR models is suggested to avoid the randomness of model generation, and selecting the best one for further the analysis and application. In addition, classification problems can be also investigated using the random forest classification algorithm.

4. Conclusions

In this research, the effect on the mechanical and tribological properties of NBR with varied ACN contents was studied using MD simulations. The mechanical and tribological properties with increasing ACN were enhanced until the content of ACN reached 40%, and these properties of ACN50 system were unimproved and descended, respectively. The entanglement and interactions of ACN50 molecules were not strengthened by increased polarity as compared to those of ACN40, owning to the stronger rigidity and shorter chain length of ACN50 molecules. In addition, the VI of RFR models is effective for obtaining the highly correlated parameter of AR, getting the TBB energy to predict the AR trend for new models. Meanwhile, the wear forms of new models showed that the vdW force at the interface of the Fe and NBR matrices was a key factor for affecting tribological properties of NBR. In summary, NBR modification is able to enhance the tribological property by increasing the polarity and decreasing the vdW force of the interface, while focusing on the rigidity and chain length in the polymer molecules.

Author Contributions

Conceptualization, S.W.; resources, E.H.; software, B.Y.; validation, R.N.; writing—original draft preparation, Q.Y.; writing—review and editing, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Liaoning Provincial Science and Technology Plan Project, grant number 2021JH 1/10400084, by National Natural Science Foundation of China (No. 52105062); by Ningbo Natural Science Foundation (No. 2021J013); and by the Open Foundation of the State Key Laboratory (No. JSY6142219202201).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We thank our classmates, sisters, and brothers, from Shenyang University of Technology, for all their help. We also thank all authors in the references, for useful insights.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shen, M.-X.; Peng, X.-D.; Meng, X.-K.; Zheng, J.-P.; Zhu, M.-H. Fretting wear behavior of acrylonitrile–butadiene rubber (NBR) for mechanical seal applications. Tribol. Int. 2016, 93, 419–428. [Google Scholar] [CrossRef]
  2. Brostow, W.; Kovačevic, V.; Vrsaljko, D.; Whitworth, J. Tribology of polymers and polymer-based composites. J. Mater. Educ. 2010, 35, 273–290. [Google Scholar]
  3. Friedrich, K. Polymer composites for tribological applications. Adv. Ind. Eng. Polym. Res. 2018, 1, 3–39. [Google Scholar] [CrossRef]
  4. Yuan, C.; Guo, Z.; Tao, W.; Dong, C.; Bai, X. Effects of different grain sized sands on wear behaviours of NBR/casting copper alloys. Wear 2017, 384–385, 185–191. [Google Scholar] [CrossRef]
  5. Friedrich, K.; Chang, L.; Haupert, F. Current and future applications of polymer composites in the field of tribology. In Composite Materials: A Vision for the Future; Springer: Berlin/Heidelberg, Germany, 2011; pp. 129–167. [Google Scholar]
  6. Al, N.M.; Antonescu, N.N.; Dinita, A.; Morosanu, M. Tribological Characterization of Some Elastomers Used at Progressive Cavity and Piston Pumps. MATEC Web Conf. 2020, 318, 01016. [Google Scholar]
  7. Li, Y.; Wang, Q.; Wang, T.; Pan, G. Preparation and tribological properties of graphene oxide/nitrile rubber nanocomposites. J. Mater. Sci. 2011, 47, 730–738. [Google Scholar] [CrossRef]
  8. Degrange, J.-M.; Thomine, M.; Kapsa, P.; Pelletier, J.; Chazeau, L.; Vigier, G.; Dudragne, G.; Guerbé, L. Influence of viscoelasticity on the tribological behaviour of carbon black filled nitrile rubber (NBR) for lip seal application. Wear 2005, 259, 684–692. [Google Scholar] [CrossRef]
  9. Dong, C.L.; Yuan, C.Q.; Bai, X.Q.; Yang, Y.; Yan, X.P. Study on wear behaviours for NBR/stainless steel under sand water-lubricated conditions. Wear 2015, 332–333, 1012–1020. [Google Scholar] [CrossRef]
  10. George, E.; Joy, J.; Anas, S. Acrylonitrile-based polymer/graphene nanocomposites: A review. Polym. Compos. 2021, 42, 4961–4980. [Google Scholar] [CrossRef]
  11. Maroufkhani, M.; Katbab, A.; Liu, W.; Zhang, J. Polylactide (PLA) and acrylonitrile butadiene rubber (NBR) blends: The effect of ACN content on morphology, compatibility and mechanical properties. Polymer 2017, 115, 37–44. [Google Scholar] [CrossRef]
  12. Maroufkhani, M.; Katbab, A.; Zhang, J. Manipulation of the properties of PLA nanocomposites by controlling the distribution of nanoclay via varying the acrylonitrile content in NBR rubber. Polym. Test. 2018, 65, 313–321. [Google Scholar] [CrossRef]
  13. Nie, R.; Wang, S.; Song, S. Abrasive Slurry Erosion Behavior of Nitrile-Butadiene Rubber with Different Acrylonitrile Contents. Mater. Sci. 2019, 26, 71–76. [Google Scholar] [CrossRef]
  14. Saba, N.; Tahir, P.M.; Jawaid, M. A Review on Potentiality of Nano Filler/Natural Fiber Filled Polymer Hybrid Composites. Polymers 2014, 6, 2247–2273. [Google Scholar] [CrossRef]
  15. Likozar, B.; Major, Z. Morphology, mechanical, cross-linking, thermal, and tribological properties of nitrile and hydrogenated nitrile rubber/multi-walled carbon nanotubes composites prepared by melt compounding: The effect of acrylonitrile content and hydrogenation. Appl. Surf. Sci. 2010, 257, 565–573. [Google Scholar] [CrossRef]
  16. Wang, L.L.; Zhang, L.Q.; Tian, M. Mechanical and tribological properties of acrylonitrile–butadiene rubber filled with graphite and carbon black. Mater. Des. 2012, 39, 450–457. [Google Scholar] [CrossRef]
  17. Yunlong, L.; Quan, W.; Shijie, W. A review on enhancement of mechanical and tribological properties of polymer composites reinforced by carbon nanotubes and graphene sheet: Molecular dynamics simulations. Compos. Part B 2018, 160, 348–361. [Google Scholar]
  18. Yuan, Q.; Li, Y.; Wang, S.; Cui, J.; He, E. Studies on effect of butadiene microstructure on physical and mechanical properties of nitrile rubber by molecular dynamics simulation. China Synth. Ind. Rubber 2022, 45, 95–99. [Google Scholar]
  19. Chukov, D.; Stepashkin, A.; Maksimkin, A.; Tcherdyntsev, V.; Kaloshkin, S.; Kuskov, K.; Bugakov, V. Investigation of structure, mechanical and tribological properties of short carbon fiber reinforced UHMWPE-matrix composites. Compos. Part B Eng. 2015, 76, 79–88. [Google Scholar] [CrossRef]
  20. Samyn, P.; Schoukens, G. Thermochemical sliding interactions of short carbon fiber polyimide composites at high pv-conditions. Mater. Chem. Phys. 2009, 115, 185–195. [Google Scholar] [CrossRef]
  21. Ma, J.; Qi, X.; Zhao, Y.; Dong, Y.; Song, L.; Zhang, Q.; Yang, Y. Polyimide/mesoporous silica nanocomposites: Characterization of mechanical and thermal properties and tribochemistry in dry sliding condition. Mater. Des. 2016, 108, 538–550. [Google Scholar] [CrossRef]
  22. Cui, J.; Zhao, J.; Wang, S.; Li, Y. A comparative study on enhancement of mechanical and tribological properties of nitrile rubber composites reinforced by different functionalized graphene sheets: Molecular dynamics simulations. Polym. Compos. 2021, 42, 205–219. [Google Scholar] [CrossRef]
  23. Karuth, A.; Alesadi, A.; Xia, W.; Rasulev, B. Predicting glass transition of amorphous polymers by application of cheminformatics and molecular dynamics simulations. Polymer 2021, 218, 123495. [Google Scholar] [CrossRef]
  24. Ramanathan, T.; Abdala, A.A.; Stankovich, S.; Dikin, D.A.; Herrera-Alonso, M.; Piner, R.D.; Adamson, D.H.; Schniepp, H.C.; Chen, X.; Ruoff, R.S.; et al. Functionalized graphene sheets for polymer nanocomposites. Nat. Nanotechnol. 2008, 3, 327–331. [Google Scholar] [CrossRef]
  25. Song, M.; Zhao, X.; Li, Y.; Chan, T.W.; Zhang, L.; Wu, S. Effect of acrylonitrile content on compatibility and damping properties of hindered phenol AO-60/nitrile-butadiene rubber composites: Molecular dynamics simulation. RSC Adv. 2014, 4, 48472–48479. [Google Scholar] [CrossRef]
  26. Li, Y.; Wang, S.; Wang, Q.; Xing, M. Molecular dynamics simulations of tribology properties of NBR (Nitrile-Butadiene Rubber)/carbon nanotube composites. Compos. Part B Eng. 2016, 97, 62–67. [Google Scholar] [CrossRef]
  27. Yunlong, L.; Shijie, W.; Enqiu, H. The effect of sliding velocity on the tribological properties of polymer/carbon nanotube composites. Carbon 2016, 106, 106–169. [Google Scholar]
  28. Hernández, M.; Bernal, M.d.M.; Verdejo, R.; Ezquerra, T.A.; López-Manchado, M.A. Overall performance of natural rubber/graphene nanocomposites. Compos. Sci. Technol. 2012, 73, 40–46. [Google Scholar] [CrossRef]
  29. Liu, F.; Hu, N.; Zhang, J.; Atobe, S.; Weng, S.; Ning, H.; Liu, Y.; Wu, L.; Zhao, Y.; Mo, F.; et al. The interfacial mechanical properties of functionalized graphene–polymer nanocomposites. RSC Adv. 2016, 6, 66658–66664. [Google Scholar] [CrossRef]
  30. Wang, Y.; Zhan, H.F.; Xiang, Y.; Yang, C.; Wang, C.M.; Zhang, Y.Y. Effect of Covalent Functionalization on Thermal Transport across Graphene–Polymer Interfaces. J. Phys. Chem. C 2015, 119, 12731–12738. [Google Scholar] [CrossRef]
  31. Suhua, Y.L.Z. Theory and Practice of Molecular Dynamics Simulation; Science Press: Beijing, China, 2013. [Google Scholar]
  32. Sun, H. COMPASS: An ab Initio Force-Field Optimized for Condensed-Phase ApplicationsOverview with Details on Alkane and Benzene Compounds. J. Phys. Chem. B 1998, 102, 7338–7364. [Google Scholar] [CrossRef]
  33. Saha, S.; Bhowmick, A.K. An Insight into molecular structure and properties of flexible amorphous polymers: A molecular dynamics simulation approach. J. Appl. Polym. Sci. 2019, 136, 47457. [Google Scholar] [CrossRef]
  34. Crispin, X.; Geskin, V.; Crispin, A.; Cornil, J.; Lazzaroni, R.; Salaneck, W.R.; Brédas, J.-L. Characterization of the Interface Dipole at Organic/Metal Interfaces. J. Am. Chem. Soc. 2002, 124, 8131–8141. [Google Scholar] [CrossRef] [PubMed]
  35. Li, X.; Zhang, Y.; Zhao, H.; Burkhart, C.; Brinson, L.C.; Chen, W. A Transfer Learning Approach for Microstructure Reconstruction and Structure-property Predictions. Sci. Rep. 2018, 8, 13461. [Google Scholar] [CrossRef]
  36. Mitra, S.; Ahmad, A.; Biswas, S.; Das, A.K. A machine learning approach to predict the structural and magnetic properties of Heusler alloy families. Comput. Mater. Sci. 2023, 216, 111836. [Google Scholar] [CrossRef]
  37. Scornet, E.; Biau, G.; Vert, J.-P. Consistency of random forests. Ann. Stat. 2015, 43, 1716–1741. [Google Scholar] [CrossRef]
  38. Strobl, C.; Boulesteix, A.-L.; Zeileis, A.; Hothorn, T. Bias in random forest variable importance measures: Illustrations, sources and a solution. BMC Bioinform. 2007, 8, 25. [Google Scholar] [CrossRef] [PubMed]
  39. Ulrike, G. Variable Importance Assessment in Regression: Linear Regression versus Random Forest. Am. Stat. 2009, 63, 308–319. [Google Scholar]
  40. Lovatti, B.P.; Nascimento, M.H.; Neto, C.; Castro, E.V.; Filgueiras, P.R. Use of Random forest in the identification of important variables. Microchem. J. 2019, 145, 1129–1134. [Google Scholar] [CrossRef]
  41. Qingzhi, Y. Modern Rubber Technology; China Petrochemical Press: Beijing, China, 1997. [Google Scholar]
  42. Polyak, B.T. The conjugate gradient method in extremal problems. USSR Comput. Math. Math. Phys. 1969, 9, 94–112. [Google Scholar] [CrossRef]
  43. Hans, C.A. Molecular dynamics simulations at constant pressure and/or temperature. J. Chem. Phys. 2008, 72, 2384–2393. [Google Scholar]
  44. Berendsen, H.J.C.; Postma, J.P.M.; Van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef]
  45. Subramaniyan, A.K.; Sun, C.T. Continuum interpretation of virial stress in molecular simulations. Int. J. Solids Struct. 2008, 45, 4340–4346. [Google Scholar] [CrossRef]
  46. Lehrbuch der Kristallphysik (mit Ausschluss der Kristalloptik). By Woldemar Voigt. Pp. xxiv + 964, with 213 text-figures and 1 table. (Leipzig and Berlin: B. G. Teubner. 1910. Price 30 marks.). Mineral. Mag. J. Mineral. Soc. 1912, 16, 253. [Google Scholar] [CrossRef]
  47. Reuss, A.J.Z.A.M.M. Calculation of the flow limits of mixed crystals on the basis of the plasticity of monocrystals. Z. Angew. Math. Mech. 1929, 9, 49–58. [Google Scholar] [CrossRef]
  48. Steinle-Neumann, G.; Stixrude, L.; Cohen, R.E. First-principles elastic constants for the hcp transition metals Fe, Co, and Re at high pressure. Phys. Rev. B 1999, 60, 791–799. [Google Scholar] [CrossRef]
  49. Rigby, D.; Sun, H.; Eichinger, B. Computer simulations of poly (ethylene oxide): Force field, pvt diagram and cyclization behaviour. Polym. Int. 1997, 44, 311–330. [Google Scholar] [CrossRef]
  50. Herrera, J.C. Equation of motion in classical electrodynamics. Phys. Rev. D 1977, 15, 492. [Google Scholar] [CrossRef]
  51. A Comparative Study of Ordinary Cross-Validation, v-Fold Cross-Validation and the Repeated Learning-Testing Methods. Biometrika 1989, 76, 503–514. [CrossRef]
  52. Zhou, Z. Machine Learning; Tsinghua University Press: Beijing, China, 2016. [Google Scholar]
  53. Iohara, K.I.; Imada, K.; Takayanagi, M. Determination of Mean Square Displacements of Polyethylene Molecules in Crystal Lattice. Polym. J. 1972, 3, 357–364. [Google Scholar] [CrossRef]
  54. Zhi, Q.C.; De Guo, W.; Qiong, Z. The Effect of Acrylonitrile Content on the Abrasion Behavior of Unfilled Nitrile Rubber. Adv. Mater. Res. 2014, 2954, 884–885. [Google Scholar]
  55. Maple, J.R.; Dinur, U.; Hagler, A.T. Derivation of Force-Fields for Molecular Mechanics and Dynamics from Abinitio Energy Surfaces. Proc. Natl. Acad. Sci. USA 1988, 85, 5350–5354. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Flow chart of the MD simulations and dataset.
Figure 1. Flow chart of the MD simulations and dataset.
Polymers 15 03799 g001
Figure 2. The RFR model. The yellow part denotes the process of new datasets for each tree established by the original dataset. The green box is the growth process of one tree.
Figure 2. The RFR model. The yellow part denotes the process of new datasets for each tree established by the original dataset. The green box is the growth process of one tree.
Polymers 15 03799 g002
Figure 3. The mechanical properties of NBR.
Figure 3. The mechanical properties of NBR.
Polymers 15 03799 g003
Figure 4. The binding energy between the NBRL and the NBR matrix during the NPT equilibrium processes, and the structure of the NBRL chain.
Figure 4. The binding energy between the NBRL and the NBR matrix during the NPT equilibrium processes, and the structure of the NBRL chain.
Polymers 15 03799 g004
Figure 5. The radius of gyration of NBR.
Figure 5. The radius of gyration of NBR.
Polymers 15 03799 g005
Figure 6. The mean square displacement of NBR.
Figure 6. The mean square displacement of NBR.
Polymers 15 03799 g006
Figure 7. Free volume of NBR.
Figure 7. Free volume of NBR.
Polymers 15 03799 g007
Figure 8. The results of 30% ACN about other NBR models.
Figure 8. The results of 30% ACN about other NBR models.
Polymers 15 03799 g008
Figure 9. Snapshots of friction processes of NBR. The red chain from the NBR matrix is the NBRL chain. The blue NBR matrix represent the wear atoms during the friction process.
Figure 9. Snapshots of friction processes of NBR. The red chain from the NBR matrix is the NBRL chain. The blue NBR matrix represent the wear atoms during the friction process.
Polymers 15 03799 g009
Figure 10. The rankings and proportions of VI for the feature parameters of the RFR.
Figure 10. The rankings and proportions of VI for the feature parameters of the RFR.
Polymers 15 03799 g010
Figure 11. (A) The BTB energy versus the ACN content of NBR models. (B) The TBB energy versus the AR of NBR models.
Figure 11. (A) The BTB energy versus the ACN content of NBR models. (B) The TBB energy versus the AR of NBR models.
Polymers 15 03799 g011
Figure 12. The TBB energy and the AR with all forecasting models.
Figure 12. The TBB energy and the AR with all forecasting models.
Polymers 15 03799 g012
Table 1. Information about the NBR models.
Table 1. Information about the NBR models.
ItemAcrylonitrile and Butadiene Allocation RatioACN ContentChain Length
(Å)
Number of Elements in Single ChainNumber of Total Atoms
CHNDouble Bond
ACN201:419.70%84.33761104164180
ACN303:729.61%80.96741046144048
ACN402:339.55%75.6672988123916
ACN501:149.53%71.74709210103784
Table 2. The dataset of the feature parameters.
Table 2. The dataset of the feature parameters.
Feature ModuleThe Name of Feature Parameters
Structural parametersACN ContentNumber of CNumber of HNumber of N
Mechanical propertiesEGB
Friction parametersFNFfAR
Environmental variablesTemperaturePressureLoadingVelocity
Energy parameters(Hamiltonian)(Total valence energy)(Non-bond energy)(Total cross terms energy)
Total kinetic energyBond energyVan der Waals energyBend–Bend energy
(Total potential energy)Angle energyElectrostatic energyBend–Torsion–Bend energy
Torsion energy Torsion–Bend–Bend energy
Inversion energy Torsion–Stretch energy
Stretch–Bend–Stretch energy
Stretch–Torsion–Stretch energy
Stretch–Stretch energy
Table 3. The mechanical properties of NBR.
Table 3. The mechanical properties of NBR.
SystemYoung’s Modulus (GPa)Increase Percentage (%)Bulk Modulus (GPa)Increase Percentage (%)Shear Modulus (GPa)Increase Percentage (%)
ACN203.050.002.950.001.250.00
ACN303.328.852.83−4.071.336.40
ACN403.7723.613.2710.851.5221.60
ACN503.8024.593.3914.921.4213.60
Table 4. The information and results of 30% ACN of NBR models with the other arrangement structures of the monomers.
Table 4. The information and results of 30% ACN of NBR models with the other arrangement structures of the monomers.
ACN30 ItemThe Number of ACN Consecutive UnitsBinding Energy (kcal/mol)Free Volume (Å3)
1 ACN2 ACN3 ACN
1211 (ACN30)410192.706972
2121220192.407049
2130111190.557254
3111301193.736655
Table 5. The VI ranking and values of RFR model results.
Table 5. The VI ranking and values of RFR model results.
The VI
Ranking
RFR
Model 1
RFR
Model 2
RFR
Model 3
RFR
Model 4
RFR
Model 5
RFR
Model 6
RFR
Model 7
RFR
Model 8
1BTB energy
(62.6%)
BTB energy
(33.2%)
BTB energy
(41.9%)
BTB energy
(26.7%)
BTB energy
(37.2%)
BTB energy
(33.1%)
BTB energy
(51.9%)
BTB energy
(41.9%)
2FN
(7.5%)
TBB energy
(17.1%)
FN
(7%)
TBB energy
(13.2%)
Ff
(9.1%)
Ff
(10.5%)
TS energy
(6.9%)
FN
(7%)
3Temperature
(6.1%)
Temperature
(8%)
TBB energy
(6.8%)
FN
(11.9%)
Temperature
(8.1%)
TBB energy
(10.1%)
Ff
(6.3%)
TBB energy
(6.8%)
4TBB energy
(3.1%)
Ff
(5.9%)
Ff
(6.4%)
Ff
(8.6%)
FN
(6.6%)
FN
(6.4%)
Temperature
(4.9%)
Ff
(6.4%)
5Ff
(2.3%)
FN
(5.3%)
STS energy
(5.3%)
Temperature
(4.4%)
TBB energy
(6.5%)
TS energy
(5.9%)
FN
(4.8%)
STS energy
(5.3%)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yuan, Q.; Li, Y.; Wang, S.; He, E.; Yang, B.; Nie, R. A Molecular Dynamics Simulation Study on Enhancement of Mechanical and Tribological Properties of Nitrile—Butadiene Rubber with Varied Contents of Acrylonitrile. Polymers 2023, 15, 3799. https://doi.org/10.3390/polym15183799

AMA Style

Yuan Q, Li Y, Wang S, He E, Yang B, Nie R. A Molecular Dynamics Simulation Study on Enhancement of Mechanical and Tribological Properties of Nitrile—Butadiene Rubber with Varied Contents of Acrylonitrile. Polymers. 2023; 15(18):3799. https://doi.org/10.3390/polym15183799

Chicago/Turabian Style

Yuan, Quan, Yunlong Li, Shijie Wang, Enqiu He, Bin Yang, and Rui Nie. 2023. "A Molecular Dynamics Simulation Study on Enhancement of Mechanical and Tribological Properties of Nitrile—Butadiene Rubber with Varied Contents of Acrylonitrile" Polymers 15, no. 18: 3799. https://doi.org/10.3390/polym15183799

APA Style

Yuan, Q., Li, Y., Wang, S., He, E., Yang, B., & Nie, R. (2023). A Molecular Dynamics Simulation Study on Enhancement of Mechanical and Tribological Properties of Nitrile—Butadiene Rubber with Varied Contents of Acrylonitrile. Polymers, 15(18), 3799. https://doi.org/10.3390/polym15183799

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop