24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2938  |  回复: 9

wangwen552

铜虫 (小有名气)

[求助] 如何利用LAMMPS产生层错 已有2人参与

请各位大侠帮忙看一下下面的内容如何用LAMMPS实现:
我想用LAMMPS计算层错能,看了几篇文献都是说把一个晶体分成两部分(比如上下两部分),一部分固定,另一部分相对做一个柏氏矢量的刚性移动,通过LAMMPS计算势能以此计算层错能,试了几次都不成功,哪位大侠知道到底怎样在LAMMPS中实现一部分固定,另一部分刚性移动啊?跪求帮忙!
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

badercao

禁虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wangwen552: 金币+10, ★★★很有帮助 2016-09-13 16:27:21
本帖内容被屏蔽

» 本帖已获得的红花(最新10朵)

2楼2016-09-09 19:29:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangwen552

铜虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by badercao at 2016-09-09 19:29:37
fix rigid
fix setforce
fix addforce

大神你好,能否详细说明一下,这几个命令我都试过,但是总是不对。我是建了一个单晶,沿着Y方向分为上下两个region(block),也是两个group,我是想固定下半部分,将上半部分沿着X方向有一个柏氏矢量的刚性位移,那个fix rigid命令是只针对上半部分吗?还有后面两条命令分别是针对哪一部分呢?还有如果用了fix rigid命令,那系综是不是默认就是nve了,可我需要LAMMPS计算的是势能,如果用了nve那最后系统的能量都相等了这不就矛盾了吗?纠结死了。。。
3楼2016-09-11 10:20:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
wangwen552: 金币+10, 有帮助 2016-09-13 16:27:05
http://wenku.baidu.com/view/9053 ... a4.html?from=search
这个上海交大引入柏氏矢量的写的还可以,但是稍微有的繁杂,你参考参考吧
4楼2016-09-11 13:12:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangwen552

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by lx_PICO at 2016-09-11 13:12:42
http://wenku.baidu.com/view/90536b3bee06eff9aef807a4.html?from=search
这个上海交大引入柏氏矢量的写的还可以,但是稍微有的繁杂,你参考参考吧

多谢!我去看一下
5楼2016-09-11 13:17:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

badercao

禁虫 (正式写手)

本帖内容被屏蔽

6楼2016-09-11 16:29:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangwen552

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by badercao at 2016-09-11 16:29:12
大神谈不上,只是针对你讲的,我觉得用上面的命令可能会实现。(1)fix rigid 命令针对group。(2)后面两条命令是加力的另外一种方法。(3)使用fix rigid也是可以计算势能的,具体check一下manual。 如果还有疑问 ...

恩,谢谢!
7楼2016-09-12 09:39:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leoluode

银虫 (小有名气)

同学,请问你现在是采用什么命令做刚性位移的啊?我现在采用fix setforce 刚性固定原子,并且加diaplace_atoms实现位移,但是感觉计算的不正确啊,不知道你现在哟想到什么好的方法没有呢?
8楼2016-11-25 16:19:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leeve

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by leoluode at 2016-11-25 16:19:43
同学,请问你现在是采用什么命令做刚性位移的啊?我现在采用fix setforce 刚性固定原子,并且加diaplace_atoms实现位移,但是感觉计算的不正确啊,不知道你现在哟想到什么好的方法没有呢?

displace_atoms 是在计算(simulations)前对原子进行位移(参手册),与计算层错时的刚性位移应该不是一个概念
9楼2018-03-24 16:43:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Stoudemire29

金虫 (小有名气)

8楼同学在吗?fix setforce 和displace_atoms的用法能介绍一下吗?新手一枚
共同进步
10楼2018-10-29 14:58:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangwen552 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +11 璃茉一定上岸 2026-04-10 12/600 2026-04-13 22:20 by zhq0425
[考研] 22专硕求调剂 +8 haoyun上岸 2026-04-11 10/500 2026-04-13 22:15 by zhq0425
[考研] 考研求调剂 +11 子木呐 2026-04-12 12/600 2026-04-13 21:28 by 学员JpLReM
[考研] 复试调剂 +17 积极向上; 2026-04-10 19/950 2026-04-13 19:53 by 闲人终南山
[考研] 245求调剂 +3 冰糖橘?汽水 2026-04-13 7/350 2026-04-13 16:42 by 冰糖橘?汽水
[考研] 290调剂生物0860 +28 哇哈哈,。 2026-04-11 31/1550 2026-04-13 01:16 by 幸免 ..
[考研] +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 0854调剂 +12 长弓傲 2026-04-09 13/650 2026-04-12 09:56 by 逆水乘风
[考研] 291求调剂 +8 关忆北. 2026-04-11 8/400 2026-04-12 09:32 by 逆水乘风
[考研] 求调剂 +11 月@163.com 2026-04-07 13/650 2026-04-11 22:55 by BruceLiu320
[考研] 求调剂 +6 小聂爱学习 2026-04-11 9/450 2026-04-11 21:20 by 蓝云思雨
[考研] 求调剂 +6 电气300求调剂不 2026-04-08 6/300 2026-04-11 20:14 by 逆水乘风
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 298求调剂 +5 残荷新柳 2026-04-07 5/250 2026-04-11 11:02 by 紫曦紫棋
[考研] 275求调剂 +9 1624447980 2026-04-08 10/500 2026-04-11 10:20 by Delta2012
[考研] 一志愿京区985,085401,与本科专业一致,电子信息工程, +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[考研] 求机械专硕297第二批调剂 +5 拾柒12。 2026-04-08 5/250 2026-04-09 16:43 by 允当适度
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 求调剂,现在还能填的 +3 上岸小莹加油 2026-04-08 3/150 2026-04-08 14:30 by zhq0425
信息提示
请填处理意见