1.一种模拟颤振环境铰链机构运动分子动力学仿真模型，其特征在于：包括压头、基体和弹簧，所述压头为圆柱体结构，用于模拟铰链机构中轴的运动；所述基体为长方体结构，通过弹簧安装于压头的下方，用于模拟铰链机构中轴孔的随动；所述压头的轴向平行于基体顶面；
所述弹簧包括第一弹簧和第二弹簧，所述第一弹簧的轴向沿z向设置，并垂直固定于基体底面中心，模拟铰链中z向运动；所述第二弹簧的轴向沿y向设置，并垂直固定于基体侧壁底边中心，模拟铰链中y向运动。
2.一种权利要求1所述模拟颤振环境铰链机构运动分子动力学仿真模型的构建方法，其特征在于具体步骤如下：
步骤一：依据铰链机构运动特性，在模拟软件中设置系统的边界条件，选取原子相互作用势函数；设定系统的初始条件，即初始位置和初始速度；并选定时间步长，构建分子动力学仿真模型；
步骤二：通过在基体的y向和z向设置弹簧，实现基体随压头振动的随动状态；
步骤三：压头y向和z向施加强迫正弦振动，同时给压头施加水平速度vs，在微正则系综下开始碰撞滑动摩擦；
步骤四：对步骤三在微正则系综下进行的碰撞滑动摩擦过程进行分子动力学仿真计算，统计计算结果，得到不同振动频率下压头与基体的位置坐标以及碰撞滑动过程中压头受到的摩擦力曲线。
3.根据权利要求2所述模拟颤振环境铰链机构运动分子动力学仿真模型的构建方法，其特征在于：所述步骤一中，构建分子动力学仿真模型的具体步骤如下：
步骤1：根据铰链机构的运动特性，在软件中建立模拟铰链机构的压头与基体模型；所述基体包括固定层、恒温层和牛顿层，将基体底部和x向的相对两侧a范围内的原子设置为固定层原子；将底部固定层上方的a范围内原子设置为恒温层原子，用以吸收摩擦过程中产生的热量；基体内其余原子设置为牛顿层原子；恒温层和牛顿层内的原子遵循牛顿第二运动定律；
步骤2：设定模拟区域x向和y向为周期性边界条件，用以减小尺寸效应；
步骤3：通过系统内不同元素原子之间相互作用势函数，获得原子之间的相互作用；
步骤4：根据maxwell-boltzmann能量分布函数设定系统内所有原子的初始速度，从而建立压头和基体原子分配的初始速度，该速度对应于其各自平衡状态的温度；设定系统温度为300k，积分步长选择为1fs；
步骤5：完成步骤1-4的设定后，在正则系综nvt下，依据所设定系统内各原子的初始位置和初始速度，根据选取的势函数计算系统中各原子的势能、位置及速度，对压头和基体原子进行平衡约束，以使系统初始模型达到平衡状态；系统弛豫完成后，在微正则系综microcanonical.nve下对系统进行内部平衡态演化，使系统达到平衡状态。
4.根据权利要求3所述模拟颤振环境铰链机构运动分子动力学仿真模型的构建方法，其特征在于：所述步骤1中，利用lammps软件构建基体为非金属硅的晶体模型，硅的晶格结构为金刚石结构，晶格常数a为0.543nm，模型在x，y和z方向上分别对应晶向[100]、[010]和[001]，各向尺寸为50a×16a×25a，包含有162416个硅原子；所述压头为半径是8a的圆柱形，材料为金刚石。
5.根据权利要求4所述模拟颤振环境铰链机构运动分子动力学仿真模型的构建方法，其特征在于：所述压头去除上半部分原子，计算时仅保留下半部分的圆柱结构，去除后压头内部共含有30679个碳原子。
6.根据权利要求3所述模拟颤振环境铰链机构运动分子动力学仿真模型的构建方法，其特征在于：所述步骤3中，包括三种相互作用势：si-si原子、c-c原子以及c-si原子。
7.根据权利要求6所述模拟颤振环境铰链机构运动分子动力学仿真模型的构建方法，其特征在于：所述si-si原子之间的相互作用通过tersoff多体势函数拟合；金刚石压头和硅基体原子之间的相互作用采用morse电势描述。
8.根据权利要求2所述模拟颤振环境铰链机构运动分子动力学仿真模型的构建方法，其特征在于：所述步骤二中，模拟颤振环境时，压头振动引起基体位置的变化，进而弹簧发生形变，将弹簧因形变产生的力均分施加于基体固定层的每一个原子上，能够模拟铰链机构中轴孔随轴的振动而产生的随动状态；所述第一弹簧和第二弹簧的刚度系数分别为k1、k2。
9.根据权利要求8所述模拟颤振环境铰链机构运动分子动力学仿真模型的构建方法，其特征在于：所述弹簧因形变产生的力是：
f＝k×(x1-x0)
其中：k为弹簧的刚度系数，x0为基体初始质心位置，x1为运动过程中基体的瞬时质心位置。
10.根据权利要求2所述模拟颤振环境铰链机构运动分子动力学仿真模型的构建方法，其特征在于：所述步骤三中，模拟铰链机构的运动特性，将轴的无规则碰撞等效为圆柱体压头在y向和z向所施加的强迫振动的耦合，振动在y向和z向的振幅分别为ay、az；频率分别为fy、fz；同时将轴的旋转运动等效为压头沿基体表面的水平滑动；利用lammps软件的fixmove命令给压头在y向和z向施加正弦强迫振动以及水平方向的滑动速度vs，在微正则系综下开始碰撞滑动摩擦过程；依据步骤二中所设定的弹簧刚度k1、k2，根据以下公式计算得到基体y向和z向的固有频率：
其中：k为弹簧的刚度系数，m为基体的质量；
改变压头的振动频率，根据步骤一中所设定的相关参数以及选取的系统内不同元素原子之间相互作用势函数，计算得到不同振动频率下，压头与基体的位置坐标以及在碰撞滑动过程中压头受到的摩擦力，进行颤振环境下的摩擦学研究的仿真分析。