在石墨烯itp文件末尾添加下面部分，这里前两列是从[bond]部分直接复制过来的，10000可以改成2000，或者自己调试出使得平面不发生波动的值。更多信息在论坛里搜索angle_restraints_z，是有帖子的，这里时间有限，不提供了。
[ angle_restraints_z ]
; #1 #2  type  theta0    fc     multiplicity
    1       2         1      90      10000     2
    1       6         1      90      10000     2
    1    3336         1      90      10000     2
    2       3         1      90      10000     2
    2      65         1      90      10000     2
    3       4         1      90      10000     2