24小时热门版块排行榜    

查看: 4432  |  回复: 11

小时候可拽啦

新虫 (初入文坛)

[交流] LAMMPS中NVT系综下如何实现单轴拉伸呢 已有2人参与

各位,在NPT系综下可以通过设置与拉伸方向垂直的两个方向应力为0,从而实现单轴拉压,比如通过fix 1 all npt temp 300 300 0.1 y 0 0 1 z 0 0 1 drag 1和fix 2 all deform 1 x 0.001 units box remap x就可以实现单轴拉伸,但是在nvt系综下,如何控制与拉伸方向垂直的两个方向的应力为0呢,从而实现单轴拉伸呢?各位坛友能否给点建议
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
用nvt就不可以控制垂直压强为零。实现你的目标只能用你说的npt。

[ 发自手机版 http://muchong.com/3g ]
2楼2018-04-20 13:23:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Alehmn

新虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
1楼: Originally posted by 小时候可拽啦 at 2018-04-20 10:43:09
各位,在NPT系综下可以通过设置与拉伸方向垂直的两个方向应力为0,从而实现单轴拉压,比如通过fix 1 all npt temp 300 300 0.1 y 0 0 1 z 0 0 1 drag 1和fix 2 all deform 1 x 0.001 units box remap x就可以实现单 ...

先用NPT控制压力,再用NVT变形

发自小木虫IOS客户端
3楼2018-04-20 22:04:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Alehmn

新虫 (著名写手)

引用回帖:
1楼: Originally posted by 小时候可拽啦 at 2018-04-20 10:43:09
各位,在NPT系综下可以通过设置与拉伸方向垂直的两个方向应力为0,从而实现单轴拉压,比如通过fix 1 all npt temp 300 300 0.1 y 0 0 1 z 0 0 1 drag 1和fix 2 all deform 1 x 0.001 units box remap x就可以实现单 ...

我最近也在弄这个可以多交流

发自小木虫IOS客户端
4楼2018-04-20 22:05:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小时候可拽啦

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by Alehmn at 2018-04-20 22:04:26
先用NPT控制压力,再用NVT变形
...

对的 我是先npt 然后nvt变形的 变形过程中 deform方向压力大约是1.8Gpa左右 而z方向能达到 1Gpa左右,边界条件是SPP 还是实现不了单轴拉伸
5楼2018-04-23 10:53:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小时候可拽啦

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by brucefan at 2018-04-20 13:23:21
用nvt就不可以控制垂直压强为零。实现你的目标只能用你说的npt。

但是manual里有说 fix nvt 可以与 fix press/berendsen连用
6楼2018-04-23 10:55:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我说的nvt是指nvt系综。如果你用fix nvt加fix press/berendsen,那就不是nvt系综了。你说的就可以实现你的目的: 控温加垂直方向控温。

[ 发自手机版 http://muchong.com/3g ]
7楼2018-04-23 13:33:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

错别字:控温加垂直方向控压。

[ 发自手机版 http://muchong.com/3g ]
8楼2018-04-23 13:34:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Alehmn

新虫 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by 小时候可拽啦 at 2018-04-23 10:53:56
对的 我是先npt 然后nvt变形的 变形过程中 deform方向压力大约是1.8Gpa左右 而z方向能达到 1Gpa左右,边界条件是SPP 还是实现不了单轴拉伸...

我设的三个维度都是周期性边界条件,NPT控压控温3百万步,最后一万步NVT等体积单轴拉伸,但是出来的曲线震荡很大。

发自小木虫IOS客户端
9楼2018-04-23 15:27:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小时候可拽啦

新虫 (初入文坛)

引用回帖:
9楼: Originally posted by Alehmn at 2018-04-23 15:27:54
我设的三个维度都是周期性边界条件,NPT控压控温3百万步,最后一万步NVT等体积单轴拉伸,但是出来的曲线震荡很大。
...

对的 我的曲线出来也是很大 我觉得应该是其他两个方向应力不为零所导致的
10楼2018-04-25 09:19:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小时候可拽啦 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 范进中举一文的中心思想 +8 炎黄贵胄 2026-08-22 9/450 2026-08-26 07:47 by WASM
[基金申请] 2026国自然函评费到账 +21 羊腰板 2026-08-21 24/1200 2026-08-26 07:20 by finnigan
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +22 木水思豆 2026-08-25 25/1250 2026-08-26 07:16 by jklily
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 PbwVEMK5haaE 2026-08-25 6/300 2026-08-26 06:57 by cNXvBfCpiZOM
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见