应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » LAMMPS源代码的求助
51/1返回列表
查看: 1594  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

frcon

新虫 (初入文坛)

[求助] 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]的意思，各位有知道的能讲解一下吗？万分感谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-02-23 14:19:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frcon

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by lsloneil at 2016-02-24 03:22:14
这很简单，只要你知道了你键能还有受力（也就是键能对距离的导数取负值），直接把ebond和fbond的函数形式改掉就行了。以bond_harmonic.cpp为例，在compute()这个函数里，你可以看到以下代码
    rsq = delx*delx  ...

谢谢你，我比较了这两个代码，发现了他们的不同，也知道了怎么去修改，非常感谢，我去试试

发自小木虫IOS客户端
赞一下 回复此楼
8楼2016-02-24 15:16:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
这很好理解啊，模拟中可以定义多种类型的键，每一种类型的键都有相应的力常数k和平衡距离r0， type就表示健的种类。

p.s. 如果不懂c++，还是补一点基础知识再去看源代码。
赞一下 回复此楼
2楼2016-02-23 15:27:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frcon

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by lsloneil at 2016-02-23 15:27:18
这很好理解啊，模拟中可以定义多种类型的键，每一种类型的键都有相应的力常数k和平衡距离r0， type就表示健的种类。

p.s. 如果不懂c++，还是补一点基础知识再去看源代码。

好的，谢谢
回复此楼
3楼2016-02-23 15:33:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frcon

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by lsloneil at 2016-02-23 15:27:18
这很好理解啊，模拟中可以定义多种类型的键，每一种类型的键都有相应的力常数k和平衡距离r0， type就表示健的种类。

p.s. 如果不懂c++，还是补一点基础知识再去看源代码。

那像代码里的ebond，fbond这些定义对应的意思可以在哪里找到
赞一下 回复此楼
4楼2016-02-23 15:35:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭