采用分子模拟方法评价聚合物材料流动性能的方法及系统与流程

本发明涉及聚合物模拟仿真技术领域，具体涉及一种采用分子模拟方法评价聚合物材料流动性能的方法及系统。

背景技术

聚合物可分为低聚物和高聚物(高分子化合物)。天然高分子化合物，以及以高聚物为基础的合成材料，例如：塑料、合成橡胶、合成纤维、涂料与粘接剂等，在国民经济中有非常重要的应用价值。

聚合物加工成各种制品的过程与聚合物的流动性密切相关。研究材料的流动性能，对聚合物的合成、加工和设计，以及进一步提高聚合物的加工工艺性能、正确选择加工工艺条件、指导配方设计等具有重要的意义。材料的流动性能在化工、石油、医药、电子工业、航空航天、化妆品、印染、造纸、建筑等诸多行业中都是一项重要的加工性能指标。聚合物的流动性可以用黏度的倒数表示。流动性(或黏度)与聚合物的结构(拓扑形貌、支化度、官能化等)、分子量、不同配比等相关。同时，施加应力的大小、应变速率、测试条件(温度、压力)、添加剂(填料、增塑剂、润滑剂、表面活性剂、增稠剂等)等也会影响聚合物的流动性能。

使用分子模拟仿真方法获得黏度数据，提供了一种理论预测材料的加工性能好坏的方法。更重要的是分子模拟仿真方法可以从材料的微观结构入手(包括分子结构、分子量、支化度、官能化、不同配方设计等，并考虑温度、压强等工艺条件，揭示影响聚合物材料流动性能的关键因素。

剪切黏度是不同黏度系数中最重要的一个，通常被简称为黏度。利用平衡分子动力学模拟方法，基于green-kubo线性响应理论，把剪切黏度看作是内部压力张量非对角线组元的自相关函数的时间积分：

其中，v为体系的体积，k为波尔兹曼常数，t为温度，pαβ指瞬时压力张量的三个等价非对角组元，<>表示统计平均。

用图1说明黏度计算过程。图1中td为动力学模拟的总时间；ts为自相关函数中pαβ(0)对应的时刻；tw为自相关函数计算的最大时间间隔。在计算过程中，pαβ(0)从0开始选取，依次选取0、ts、2ts、3ts、……、td-tw时刻的压力张量非对角组元，对于每一个pαβ(0)，pαβ(t)则为从pαβ(0)对应时刻开始，tw时间间隔之内的所有时刻的压力张量非对角组元，由此对于任一个pαβ(0)，会有若干pαβ(t)，二者的乘积pαβ(t)pαβ(0)也将有若干个，这些乘积的平均值即为自相关函数中的统计平均，对所有的tw进行积分并乘上前面的常数系数，即可得到最后的剪切黏度。

目前不少科研工作者利用该平衡分子动力学模拟方法，分析应力自相关函数获取黏度数据。但是，该方法受到软件使用人员经验和使用水平的限制、操作难度大、效率低、无法大规模推广和计算。而且，研究结果准确性差、数据差异大，对于新型材料的开发和设计的指导和借鉴意义不足。

技术实现要素：

针对现有技术中的缺陷，本发明的目的之一在于提供一种采用分子模拟方法评价聚合物材料流动性能的方法，能快捷、高效、准确地计算出聚合物材料的剪切黏度数据。

第一方面，本发明提供的一种采用分子模拟方法评价聚合物材料流动性能的方法，具体包括以下步骤：

s1:根据待测聚合物的重复单元，搭建待测聚合物的单链结构；根据待测聚合物的结构参数、添加剂、配方组分种类和配方组分比例搭建聚合物的体相结构模型；

s2:对所述待测聚合物的体相结构模型依次进行构型的初次驰豫；nvt系综全局最优构型的初次筛选；构型的再次驰豫；npt系综全局最优构型的二次筛选；nvt及npt系综的降温退火模拟,进而使模型的构型和密度充分平衡，得到第二结构数据；

s3:对所述第二结构数据进行nvt系综或者npt系综的充分的动力学平衡，得到动力学轨迹数据；

s4:根据所述动力学轨迹数据分析应力自相关函数，并测试黏度数据的有效性，进而得到待测聚合物的可靠的剪切黏度数据。

进一步地，所述s2的具体方法包括：

s21:对所述聚合物的体相结构模型使用分子力学方法驰豫原子位置，得到第一驰豫结构；

s22：对所述第一驰豫结构进行高频、多次的nvt系综的升降温退火模拟计算，在每一次退火计算中以设定的温度间隔进行短时长的nvt系综的分子动力学平衡，得到第一能量最低构型；

s23:对所述第一能量最低构型使用分子力学的方法驰豫原子位置，并优化晶胞参数，得到第二驰豫结构；

s24:对所述第二驰豫结构进行高频、多次的npt系综的升降温退火模拟计算，在每一次退火计算中以设定的温度间隔进行短时长的npt系综的分子动力学平衡，再次筛选该条件下的全局最优构型，得到第二能量最低构型；

s25:对所述第二能量最低构型进行nvt系综降温退火模拟计算，得到第一结构数据；

s26：对所述第一结构数据进行长时间的npt系综降温退火模拟计算，得到第二结构数据。

进一步地，s3中对所述第二结构数据进行动力学平衡的具体方法包括s31:对所述第二结构数据进行驰豫原子位置，并优化晶胞参数得到第三结构数据。

进一步地，s3中对所述第二结构数据进行动力学平衡的具体方法还包括s32：对所述第三结构数据采用分子动力学方法进行npt系综的长时间动力学平衡，得到第一动力学轨迹数据。

进一步地，s3中对所述第二结构数据进行动力学平衡的具体方法还包括s33：对所述第三结构数据采用分子动力学方法进行nvt系综的长时间的动力学平衡，得到第二动力学轨迹数据。

进一步地，s3中对所述第二结构数据进行动力学平衡的具体方法还包括s34：对所述第一动力学轨迹数据采用分子动力学方法进行nvt系综的长时间的动力学平衡，得到第三动力学轨迹数据。

第二方面，本发明提供的一种采用分子模拟方法评价聚合物材料流动性能的系统，包括模型搭建模块、构型驰豫模块、动力学平衡模块和分析模块，所述模型搭建模块用于根据待测聚合物的重复单元，搭建待测聚合物的单链结构；根据待测聚合物的结构参数、添加剂、配方组分种类和配方组分比例搭建聚合物的体相结构模型；所述构型驰豫模块用于对所述聚合物的体相结构模型依次进行驰豫、nvt系综全局最优构型的初次筛选及模型的构型和密度的充分平衡，得到第二结构数据；所述动力学平衡模块用于对所述第二结构数据进行动力学平衡，得到动力学轨迹数据；所述分析模块用于根据所述动力学轨迹数据分析应力自相关函数，并测试黏度数据的有效性，得到待测聚合物的可靠的剪切黏度数据。

进一步地，所述构型驰豫模块包括模型原子位置驰豫单元、全局最优构型的初次筛选单元和模型的构型和密度的充分平衡单元；所述模型原子位置驰豫单元用于对所述聚合物的体相结构模型使用分子力学方法驰豫原子位置，得到第一驰豫结构；所述全局最优构型的初次筛选单元用于对所述第一驰豫结构进行高频、多次的nvt系综的升降温退火模拟计算，在每一次退火计算中以设定的温度间隔进行短时长的nvt系综的分子动力学平衡，得到第一能量最低构型；所述模型的构建和密度的充分平衡单元用于对所述第二驰豫结构进行高频、多次的npt系综的升降温退火模拟计算，在每一次退火计算中以设定的温度间隔进行短时长的npt系综的分子动力学平衡，再次筛选该条件下的全局最优构型，得到第二能量最低构型；对所述第二能量最低构型进行nvt系综降温退火模拟计算，得到第一结构数据；对所述第一结构数据进行长时间的npt系综降温退火模拟计算，得到第二结构数据。

进一步地，动力学平衡模块包括体相模型优化单元，所述体相模型优化单元用于对所述第二结构数据进行驰豫原子位置，并优化晶胞参数得到第三结构数据。

进一步地，动力学平衡模块还包括npt系综动力学平衡单元或nvt系综动力学平衡单元，所述npt系综动力学平衡单元用于对所述第三结构数据采用分子动力学方法进行npt系综的长时间动力学平衡，得到第一动力学轨迹数据；所述nvt系综动力学平衡单元用于对所述第三结构数据采用分子动力学方法进行nvt系综的长时间的动力学平衡，得到第二动力学轨迹数据。

本发明的有益效果：

本发明实施例提供的一种采用分子模拟方法评价聚合物材料流动性能的方法将繁琐复杂的操作步骤流程化、简单化，极大减少了操作步骤中的人为因素，降低了操作门槛，使大规模计算成为可能；通过抓住操作流程中的关键技术问题，固化操作过程，保证操作流程的合理性和可重复性，并通过流程化操作确保批量仿真实验的数据的一致性和可靠性，解决了现有技术存在的问题。其有益效果:

1.通过引入模拟仿真减少实验次数，加快研发进度；

2.仿真建模的便捷性使该方法尤其适合复杂配方体系；

3.操作流程明确，简单易用，降低操作门槛，能够快速推广；

4.操作步骤流程化，确保仿真结果可靠性和重复性；

5.有利于提高开发新型的加工性能优良的聚合物的效率。

本发明实施例提供的一种采用分子模拟方法评价聚合物材料流动性能的系统，有益效果:

1.通过建立仿真模型，减少实验次数，可加快研发进度；

2.系统操作简单、降低操作门槛，能够快速推广；

3.系统内采用的方法步骤流程化，确保仿真结果可靠性和重复性；

4.有利于提高开发新型的加工性能优良的聚合物的效率。

附图说明

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍。在所有附图中，类似的元件或部分一般由类似的附图标记标识。附图中，各元件或部分并不一定按照实际的比例绘制。

图1示出了现有技术聚合物剪切黏度计算方法示意图；

图2示出了本发明所提供的一种采用分子模拟方法评价聚合物材料流动性能的方法的流程图；

图3示出了图2中步骤s2的具体流程图；

图4示出了图2中步骤s3的一种实施例的具体流程图；

图5示出了图2中步骤s3的另一种实施例的具体流程图；

图6示出了用本发明的采用分子模拟方法评价聚合物材料流动性能的方法计算聚酰亚胺的剪切黏度值仿真图；

图7示出了本发明所提供的一种采用分子模拟方法评价聚合物材料流动性能的系统的原理框图。

具体实施方式

下面将结合附图对本发明技术方案的实施例进行详细的描述。以下实施例仅用于更加清楚地说明本发明的技术方案，因此只是作为示例，而不能以此来限制本发明的保护范围。

需要注意的是，除非另有说明，本申请使用的技术术语或者科学术语应当为本发明所属领域技术人员所理解的通常意义。

如图2所示，示出了本发明实施例所提供的一种采用分子模拟方法评价聚合物材料流动性能的方法的流程图，具体包括以下步骤：

s1:根据待测聚合物的重复单元，搭建待测聚合物的单链结构；根据待测聚合物的结构参数、添加剂、配方组分种类和配方组分比例搭建聚合物的体相结构模型。

具体地，定义待测聚合物的重复单元，搭建待测聚合物的单链结构；使用经典分子力学的方法对其进行结构优化，选择优化算法，输入循环次数、力场、电荷、加和方式等输入参数，搭建待测聚合物的单链结构。

使用蒙特卡洛的方法，考虑配方组分和比例，输入聚合物或者添加剂各组分的单链或分子结构和数目，搭建聚合物材料的无定型体相结构的低密度虚拟样品模型。设置的参数包括：模型的初始密度、配方组分结构、配方比例、力场、电荷、加和方式，得到聚合物的体相结构模型。

s2:对所述聚合物的体相结构模型依次进行构型的初次驰豫、nvt系综全局最优构型的初次筛选，构型的再次驰豫，npt系综全局最优构型的二次筛选，nvt及npt系综的降温退火模拟,使模型的构型和密度充分平衡，得到第二结构数据。

具体地，如图3所示，示出了步骤s2的具体流程图，步骤s2具体包括：

s21:对所述聚合物的体相结构模型使用分子力学方法驰豫原子位置，得到第一驰豫结构。

具体地，在步骤s21中设置优化算法、循环次数、力场、电荷、加和方式，使用分子力学方法驰豫原子位置，得到第一驰豫结构。

s22：对所述第一驰豫结构进行高频、多次的nvt系综的升降温退火模拟计算，在每一次退火计算中以设定的温度间隔进行短时长的nvt系综的分子动力学平衡，得到第一能量最低构型。

具体地，对第一驰豫结构进行高频、多次的nvt系综的退火模拟计算。退火过程从低温升到高温，然后在降到低温，完成一次退火循环过程。在每一次的退火过程中，以一定的温度间隔进行短时长的nvt系综的分子动力学平衡，最后得到能量最低构型。基于多次的退火循环流程中得到的结构，筛选出能量最低构型为第一能量最低构型，用于下一步骤的计算。此步骤须设置的参数包括：系综(nvt)、控温方法、动力学平衡步数、退火高温、退火低温、退火温度间隔、力场、电荷、加和方式。

s23:对第一能量最低构型使用分子力学的方法驰豫原子位置，并优化晶胞参数，得到第二驰豫结构。

具体地，步骤s23必须设置的参数包括：优化算法、循环次数、力场、电荷、加和方式。

s24:对所述第二驰豫结构进行高频、多次的npt系综的升降温退火模拟计算，在每一次退火计算中以设定的温度间隔进行短时长的npt系综的分子动力学平衡，再次筛选该条件下的全局最优构型，得到第二能量最低构型。

具体地，对第二驰豫结构进行高频、多次的npt系综的升降温退火模拟计算。退火过程从低温升到高温，然后在降到低温，完成一次退火循环过程。在每一次的退火过程中，以一定的温度间隔进行短时长的npt系综的分子动力学平衡，最后得到能量最低构型。基于多次的退火循环流程中得到的结构，筛选出能量最低构型为第二能量最低构型，用于下一步骤的计算。此步骤涉及必须设置的参数包括：系综(npt)、控压方法、动力学平衡步数、退火高温、退火低温、退火温度间隔、力场、电荷、加和方式。

s25:对所述第二能量最低构型进行nvt系综降温退火模拟计算，得到第一结构数据。

具体地，对第二能量最低构型nvt系综的退火平衡，得到第一结构数据。退火过程从高温到低温，每隔一定的温度间隔做一段分子动力学，单步动力学平衡时间足够长(近纳秒)。步骤s25必须设置的参数包括：系综(nvt)、动力学平衡时长、控温方法、退火高温、退火低温、退火温度间隔、力场、电荷、加和方式。

s26：对所述第一结构数据进行长时间的npt系综降温退火模拟计算，得到第二结构数据。

具体地，所述第一结构数据进行长时间的npt系综降温退火模拟计算，可以充分平衡体系的密度，得到第二结构数据。退火过程从高温到低温，每隔一定的温度间隔做动力学，单步动力学平衡时间足够长(近纳秒)。此步骤必须设置的参数包括：系综(npt)、动力学平衡时长、控压方法、退火高温、退火低温、温度间隔、力场、电荷、加和方式。

s3:对所述第二结构数据进行动力学平衡，得到动力学轨迹数据。

具体地，如图4所示，示出了步骤s3的具体流程图，步骤s3对所述第二结构数据进行动力学平衡，得到动力学轨迹数据的具体方法包括：

s31:对所述第二结构数据进行驰豫原子位置，并优化晶胞参数得到第三结构数据。

具体地，步骤s31必须设置的参数包括：优化算法、循环次数、力场、电荷、加和方式。通过对第二结构数据进行驰豫原子位置，并优化晶胞参数，进一步驰豫结构和平衡密度，得到第三结构数据。

s32：对所述第三结构数据采用分子动力学方法进行npt系综的长时间动力学平衡，得到第一动力学轨迹数据。

步骤s32是模拟实际样品需要考虑的温度和压强工况条件。此步骤必须设置的参数包括：系综(npt)、温度、压力、步长、模拟时长、步数、得到轨迹帧数、力场、电荷、加和方式。

s34：对所述第一动力学轨迹数据采用分子动力学方法进行nvt系综的长时间的动力学平衡，得到第三动力学轨迹数据。

具体地，步骤s34必须设置的参数包括：系综(nvt)、温度、步长、模拟时长、步数、得到轨迹帧数、力场、电荷、加和方式。

步骤s3还有另外一种情况，如图5所示，示出了步骤s3的具体流程图，对第二结构数据进行模型原子位置和晶胞参数驰豫，再进行nvt系综分子动力学模拟。

对第二结构数据进行动力学平衡，得到动力学轨迹数据的具体方法包括：

s31:对所述第二结构数据进行驰豫原子位置，并优化晶胞参数得到第三结构数据。

具体地，步骤31必须设置的参数包括：优化算法、循环次数、力场、电荷、加和方式。通过对第二结构数据进行驰豫原子位置，并优化晶胞参数，得到第三结构数据。

s33：对第三结构数据采用分子动力学方法进行nvt系综的长时间的动力学平衡，得到第二动力学轨迹数据。

具体地，步骤s33必须设置的参数包括：系综(nvt)、温度、步长、模拟时长、步数、得到轨迹帧数、力场、电荷、加和方式。

s4:根据所述动力学轨迹数据分析应力自相关函数，并测试黏度数据的有效性，进而得到待测聚合物的可靠的剪切黏度数据。

将第一动力学轨迹数据、第二动力学轨迹数据或第三动力学轨迹数据进行应力自相关函数分析，得到待测聚合物的剪切黏度数据。

根据图1的示意图，选取合适的ts(例如：1帧，如果轨迹帧数比较多，适当增加该值)，tw依次以一定的间隔逐渐增大(例如：10帧)。不同的tw设置，均可得到一组剪切黏度的数据，根据此数据可以作出黏度随tw变化的曲线。如果该曲线随着tw的增加出现缓台，读取缓台对应的黏度数据。如果持续增加，没有缓台出现，重复计算得到第一动力学轨迹数据、第二动力学轨迹数据或第三动力学轨迹数据的步骤，增加模拟时长以及得到轨迹的帧数，重新作出剪切黏度随tw变化的曲线，看曲线的变化趋势，确定是否结束计算。如果想预测准确、可靠的剪切黏度数据需要充分的结构驰豫和长时间的动力学平衡过程、分析大量的采样样本(得到的动力学轨迹包含的采样帧数)、合适的非键相互作用能的处理方法等。剪切黏度数据本身对计算参数选取的敏感性比较高，需要细致测试参数。

本发明实施例提供的一种采用分子模拟方法评价聚合物材料流动性能的方法将繁琐复杂的操作步骤流程化、简单化，极大减少了操作步骤中的人为因素，降低了操作门槛，使大规模计算成为可能；通过抓住操作流程中的关键技术问题，固化操作过程，保证操作流程的合理性和可重复性，并通过流程化操作确保批量仿真实验的数据的一致性和可靠性，解决了现有技术存在的问题。

本发明实施例提供的一种采用分子模拟方法评价聚合物材料流动性能的方法的优势如下：

1.通过引入模拟仿真减少实验次数，加快研发进度；

2.仿真建模的便捷性使该方法尤其适合复杂配方体系；

3.操作流程明确，简单易用，降低操作门槛，能够快速推广；

4.操作步骤流程化，确保仿真结果可靠性和重复性；

5.有利于提高开发新型的加工性能优良的聚合物的效率。

实施例1

聚酰亚胺复合材料的剪切黏度的计算

1、根据聚酰亚胺的重复单元，搭建聚酰亚胺的单链结构，然后使用蒙特卡洛的方法搭建体相结构模型，设置体相结构模型的初始密度、聚酰亚胺配方组分结构、聚酰亚胺配方比例、力场、电荷、加和方式；

2、用分子力学方法初步驰豫聚酰亚胺的体相结构的原子位置，得到第一驰豫结构，对第一驰豫结构进行高频、多次的nvt系综的升降温的退火模拟计算，设置的参数包括系综(nvt)、控温方法、动力学平衡步数、退火高温、退火低温、退火温度间隔、力场、电荷、加和方式，在每一次退火计算中以设定的温度间隔进行短时长的nvt系综的分子动力学平衡，得到第一能量最低构型；对第一能量最低构型使用分子力学的方法驰豫原子位置，并优化晶胞参数，得到第二驰豫结构，设置的参数包括：优化算法、循环次数、力场、电荷、加和方式；对第二驰豫结构进行高频、多次的npt系综的升降温退火模拟计算，在每一次退火计算中以设定的温度间隔进行短时长的npt系综的分子动力学平衡，得到第二能量最低构型，设置的参数包括：系综(npt)、控压方法、动力学平衡步数、退火高温、退火低温、退火温度间隔、力场、电荷、加和方式；对第二能量最低构型进行nvt系综降温退火模拟计算，温度从800k逐渐降到300k,每降低100k做一次nvt系综的动力学模拟,单个温度点的动力学模拟时长达到100ps，得到第一结构数据，设置的参数包括：系综(nvt)、动力学平衡时长、控温方法、退火高温、退火低温、退火温度间隔、力场、电荷、加和方式；对第一结构数据进行长时间的npt系综降温退火模拟计算，温度从800k逐渐降到300k,每降低100k做一次npt系综的动力学,每个温度点的动力学模拟时长达到400ps，得到第二结构数据，设置的参数包括：系综(npt)、动力学平衡时长、控压方法、退火高温、退火低温、温度间隔、力场、电荷、加和方式。

3、对经过上述复杂处理的聚酰亚胺结构，再次优化结构，不仅驰豫原子位置并优化晶胞参数，进一步驰豫结构、平衡密度，得到第三结构数据，设置参数包括：优化算法、循环次数、力场、电荷、加和方式，对所述第三结构数据采用分子动力学方法进行nvt系综的长时间的动力学平衡，得到第二动力学轨迹数据，设置的参数包括：系综(nvt)、温度、步长、模拟时长、步数、得到轨迹帧数、力场、电荷、加和方式。

4、根据第二动力学轨迹数据，进行应力自相关函数分析，模拟时长设置为90ns，温度设置为573.15k，时间步长设置为1fs，得到90万帧的轨迹。通过统计该动力学轨迹的应力自相关函数得到的剪切黏度数值约为1.107pa.s，比较接近实验测得的极低剪切速率下的黏度值，如图6所示。

如图7所示，示出了本发明实施例提供一种采用分子模拟方法评价聚合物材料流动性能的系统的原理框图，适应于上述实施例描述的方法，系统包括模型搭建模块1、构型驰豫模块2、动力学平衡模块3和分析模块4，所述模型搭建模块1用于根据待测聚合物的重复单元，搭建待测聚合物的单链结构；根据待测聚合物的结构参数、添加剂、配方组分种类和配方组分比例搭建聚合物的体相结构模型；所述构型驰豫模块2用于对所述聚合物的体相结构模型依次进行驰豫、nvt系综全局最优构型的初次筛选及模型的构型和密度的充分平衡，得到第二结构数据；所述动力学平衡模块3用于对所述第二结构数据进行动力学平衡，得到动力学轨迹数据；所述分析模块4用于根据所述动力学轨迹数据分析应力自相关函数，并测试黏度数据的有效性，得到待测聚合物的可靠的剪切黏度数据。

构型驰豫模2块包括模型原子位置驰豫单元21、全局最优构型初次筛选单元22和模型的构型和密度的充分平衡单元23；所述模型原子位置驰豫单元21用于对所述聚合物的体相结构模型使用分子力学方法驰豫原子位置，得到第一驰豫结构；所述全局最优构型初次筛选单元22用于对所述第一驰豫结构进行高频、多次的nvt系综的升降温退火模拟计算，在每一次退火计算中以设定的温度间隔进行短时长的nvt系综的分子动力学平衡，得到第一能量最低构型；所述模型的构型和密度的充分平衡单元23用于对所述第二驰豫结构进行高频、多次的npt系综的升降温退火模拟计算，在每一次退火计算中以设定的温度间隔进行短时长的npt系综的分子动力学平衡，再次筛选该条件下的全局最优构型，得到第二能量最低构型；对所述第二能量最低构型进行nvt系综降温退火模拟计算，得到第一结构数据；对所述第一结构数据进行长时间的npt系综降温退火模拟计算，得到第二结构数据。

动力学平衡模块3包括体相模型优化单元31，所述体相模型优化单元31用于对所述第二结构数据进行驰豫原子位置，并优化晶胞参数得到第三结构数据。

动力学平衡模块3还包括npt系综动力学平衡单元32或nvt系综动力学平衡单元33，所述npt系综动力学平衡单元32用于对所述第三结构数据采用分子动力学方法进行npt系综的长时间动力学平衡，得到第一动力学轨迹数据；所述nvt系综动力学平衡单元33用于对所述第三结构数据采用分子动力学方法进行nvt系综的长时间的动力学平衡，得到第二动力学轨迹数据。

本发明实施例提供的一种采用分子模拟方法评价聚合物材料流动性能的系统的工作方法参见上述实施例描述的方法流程，在此不再赘述。

本发明实施例提供的一种采用分子模拟方法评价聚合物材料流动性能的系统的优势如下：

1.通过建立仿真模型，减少实验次数，可加快研发进度；

2.系统操作简单、降低操作门槛，能够快速推广；

3.系统内采用的方法步骤流程化，确保仿真结果可靠性和重复性；

4.有利于提高开发新型的加工性能优良的聚合物的效率。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。