| 查看: 1662 | 回复: 7 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
LAMMPS源代码的求助 已有1人参与
|
|||||
|
给位老师，大家好，最近我在修改lammps中bond_harmonic的源代码，对其中部分不是很清楚，在这里请教一下： bond_harmonic的势函数为E=K（r-r0)^2, 在源代码中 rsq = delx*delx + dely*dely + delz*delz; r = sqrt(rsq); dr = r - r0[type]; rk = k[type] * dr; 其势函数是这样写的，在这里，我不是很懂程序中[type]的意思，各位有知道的能讲解一下吗？万分感谢！ |
回复此楼» 猜你喜欢
一志愿北京化工大学材料与化工（085600）296求调剂
已经有22人回复
-
322求调剂：一志愿湖南大学 材料与化工（085600），已过六级。
已经有7人回复
-
332求调剂
已经有5人回复
-
310求调剂
已经有10人回复
-
26考研-291分-厦门大学（085601）-柔性电子学院材料工程专业求调剂
已经有5人回复
-
085600 材料与化工 329分求调剂
已经有15人回复
-
化学0703 调剂 306分 一志愿211
已经有8人回复
-
342求调剂
已经有4人回复
-
334分 一志愿武理 材料求调剂
已经有14人回复
-
材料科学与工程 317求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps势函数问题
已经有2人回复
- lammps 计算粘度 lj单位 转化为 metal单位 出错
已经有7人回复
- 流固耦合方向德国博士，回国的话，国内哪些高校(或中科院)比较对口？
已经有92人回复
- lammps中能否在源代码中自己编辑公式
已经有1人回复
- 关于LAmmps的windows下源程序编译调试问题
已经有7人回复
- 关于towhee 回答zhang_jaj的问题
已经有44人回复
- 如何在lammps里加入新的势函数形式？
已经有6人回复
- 【求助】更改源代码 并行出问题
已经有9人回复
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
|
这很简单，只要你知道了你键能还有受力（也就是键能对距离的导数取负值），直接把ebond和fbond的函数形式改掉就行了。以bond_harmonic.cpp为例，在compute()这个函数里，你可以看到以下代码 你需要的就是在compute()还有single()两个函数里把fbond和ebond的函数形式改成你所需要的形式。如果你不太确定怎么改，可以参考bond_morse.cpp，也就是把harmonic bond换成morse bond，看看代码是怎么变动的。 另外如果你定义键能的参数也发生了变动（比方说你需要更多的参数来定义键能），同时也需要对coeff()这个函数进行修改来改变参数传递，具体例子参考bond_morse.cpp和bond_harmonic.cpp之间的区别，这个不难的。 |
7楼2016-02-24 03:22:14
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2016-02-23 15:27:18
3楼2016-02-23 15:33:43
4楼2016-02-23 15:35:20
