lammps - 在灯中添加外部非恒定力场

Question

我想添加一个取决于原子位置的力场。我试图计算原子位置，将其用作变量，然后设置力场，但我得到了消息：

错误：修复 addforce 的变量名不存在

我试过的代码是：

variable a equal c_X
variable extforce equal ${maxforce}*step/${eqrun}*a
fix externforce A  addforce v_extforce 0 0 

Answer 1

如果您希望它根据原子位置而变化，则需要使用原子样式变量。我不知道你的任何其他变量是什么，但作为一个例子，这增加了一个力的 y 分量作为 x 位置的函数：

variable addfx atom x*0.2
fix externforce all addforce 0. v_addfx 0.

您得到的错误意味着您需要首先“使用”变量，即使用 thermo_style 命令输出它。