基于石墨烯等二维材料搭建的纳米电力谐振器由于具有极低的质量和极高的品质因子而受到广泛关注。Fermi-Pasta-Ulam-Tsingou（FPU）问题关注基本模式上的能量由于模式间的耦合而产生耗散的机制,而基本模式的耗散是限制器件品质因子的固有因素,因此FPU物理对研究石墨烯振荡器的性质具有重要意义。本学位论文关注石墨烯振荡器中的非线性振动,发展了一套计算各本征模式间耦合强度的模式耦合方法,在该方法的帮助下分析了石墨烯系统基本模式的能量耗散路径、计算了石墨烯系统达到能均分的时间尺度,进一步地分析了石墨烯的非线性振动对其电子输运性质的影响。第一章主要介绍了石墨烯振荡器的柔性模及其对石墨烯力学、热学和电学性质的影响,随后介绍了 FPU问题的简要历史、主要结果,在本章的最后简略介绍了分子动力学模拟的基本原理。第二章首先介绍了在分子动力学模拟中使用的势能——VFF模型及其简化形式,之后推导了 VFF模型的本征模式、本征频率和色散关系,特别的,所有本征模式根据其位移构型的对称性可以分为四类。本章最后我们介绍了分子动力学模拟的初步结果,原子间势能中的非线性会导致初始激发能量从第一个本征模式传递到其它高频模式上以及导致石墨烯的振动频率和初始相对形变近似成正比。第三章主要研究了石墨烯柔性模之间的耦合以及在系统热化过程中能量在本征模式间的传递路径。本章的关键是发展了一套普适的方法,利用这套方法可以数值地计算出本征模式间的耦合强度。具体的,方形石墨烯系统的所有本征模式可以根据空间位移对称性被分为四类,同类本征模式间的耦合强度要远远大于不同类本征模式间的耦合强度。因此,对于和初始激发模式属于同一类的本征模式,其动力学可以用受驱谐振子方程来描述,这部分模式在系统演化的初始阶段就可以获得能量,能量在这部分本征模式间将达到稳定的分布,该分布和模式耦合具有相同的层级结构,此时系统所有的能量都局域在少数的模式上,导致在数值模拟中形成“Natural package”和亚稳态。对于和初始激发模式属于不同类的本征模式,其动力学满足Mathieu方程,这类模式在数值模拟中获得能量的充分条件是存在某个模式,其运动方程中的参数落在Mathieu方程参数空间的不稳定区间。由于同类模式间具有较强的模式耦合,因此和该不稳定模式属于同一类的本征模式都会被激发起来,在数值模拟中获得能量。在第四章中,我们利用分子动力学模拟方法研究了石墨烯柔性模系统达到能均分的时间尺度问题。分子动力学的结果表明当初始时刻只有一个模式被激发时,初始激发能量将很快传递到与其属于同一类的模式中,对于与其属于不同类的模式而言,由于能量壁垒的存在,它们将需要很长的时间来获得能量。能均分的时间满足如下拉伸指数关系τ∝exp（∈-α）:1)当初始激发能量较小时,高频模式相比较低频模式需要更长的时间达到均分,但是在高能激发下,高频模式更容易达到均分;2)在相同能量激发下,系统尺寸越大,其达到能均分所需要的时间越短。随着尺寸的增长,系统达到均分的时间尺度渐渐收敛到一条曲线上,进一步的,当系统尺寸大于8nm时,时间尺度中参数α将达到一稳定值,这意味着该时间尺度在热力学极限下依然成立。第五章主要结合分子动力学和紧束缚近似方法研究了石墨烯量子点的机械振动对其电输运性质的影响。在分子动力学模拟中如果仅仅考虑线性相互作用,那么系统的机械振动将是周期的,此时系统的量子输运系数也会随之呈现周期变化:在小振幅振动影响下,输运系数等高线中的峰线被拓宽成了带状或高台状,导致输运系数增强;大振幅振动会导致共振散射,继而引起输运系数的急剧涨落。如果在数值模拟中考虑原子间非线性相互作用,系统的机械振动将变得非常复杂,这将导致输运系数斑图也变得非常复杂,输运系数曲线带交织在一起,导致电流曲线中出现异常的高频成分。最后一章对全文进行了总结,并对以后的工作提出了一些展望。