| 查看: 2390 | 回复: 6 | |||
Alehmn新虫 (著名写手)
|
[求助]
LAMMPS 分子动力学计算纳米高分子复合体系拉伸曲线时遇到的问题 已有1人参与
|
|
求大神赐教怎样在NPT系综下同时实现某一个方向的压力保持不变的前提下使整个盒子的体积也保持不变 主要是为了计算体系的stress-strain behavior , 使盒子在x 方向拉长，为了体积不变那么y 和z方向按理说应该也要缩小，但是同时为了准确计算在x方向的受力就不得不保持y和z方向的压强为0，如下是我的in文件，运行后在最后的15000开始时报错： ERROR: Cannot use fix npt and fix deform on same component of stress tensor PS：穷人只有6.6金币了  加我我发红包或者以后还答主吧# VARIABLES variable fname index PE.data variable simname index PE # in.file dimension 3 boundary p p p units lj atom_style molecular neighbor 0.3 bin neigh_modify delay 2 every 2 one 20000 page 1000000 special_bonds lj 1.0 1.0 1.0 read_data PE.data bond_style fene bond_coeff 1 30.0 1.5 0.0 1.0 bond_coeff 2 30.0 1.5 0.0 1.0 angle_style cosine angle_coeff 1 10.00000 angle_coeff 2 10.00000 pair_style lj/expand 2.500000 pair_modify shift yes pair_coeff 1*2 1*2 1.0 1.0 0.0 1.12246 pair_coeff 1*2 3 1.000 1.00 2.00 2.50 pair_coeff 3 3 3.0 1.0 4.000000 1.122460 group chain type 1 2 group np type 3 fix 1 all nve/limit 0.0005 thermo 50000 dump 1 all custom 50000 dump.all.lammpstrj id type x y z timestep 0.005 run 50000 write_restart run1.restart unfix 1 velocity chain create 2.0 733994 dist gaussian fix 3 all nvt temp 2.0 1.0 50.0 run 1500000 fix 3 all nvt temp 1.0 1.0 50.0 run 1500000 unfix 3 variable e equal "(lx-35.0)/35.0" variable p equal "((pxx+pyy+pzz)/3-pxx)*1.5" fix 4 all npt temp 1.0 1.0 50 y 0 0 1000 z 0 0 1000 drag 2 fix 5 all deform 1 x scale 2.0 y volume z volume fix def1 all print 100 "${e} ${p}" file ${simname}.def1.txt screen no thermo_style custom step temp pxx pyy pzz lx ly lz epair ebond eangle vol v_e v_p thermo 300 timestep 0.005 run 15000 |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有22人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有3人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
-
范进中举一文的中心思想
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有27人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
2楼2018-04-12 03:17:23
Alehmn
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 4959.1
- 散金: 437
- 红花: 3
- 帖子: 1038
- 在线: 74.6小时
- 虫号: 3986638
- 注册: 2015-07-23
- 性别: GG
- 专业: 高分子物理与高分子物理化
3楼2018-04-12 14:41:40
brucefan
专家顾问 (著名写手)
-
专家经验: +218 - 模拟EPI: 4
- 应助: 422 (硕士)
- 金币: 2294.8
- 散金: 1720
- 红花: 88
- 帖子: 1423
- 在线: 670.2小时
- 虫号: 3048965
- 注册: 2014-03-14
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 分子模拟
Alehmn4楼2018-04-12 23:46:03
Alehmn
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 4959.1
- 散金: 437
- 红花: 3
- 帖子: 1038
- 在线: 74.6小时
- 虫号: 3986638
- 注册: 2015-07-23
- 性别: GG
- 专业: 高分子物理与高分子物理化
5楼2018-05-03 21:01:27
6楼2019-08-10 16:36:36
Alehmn
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 4959.1
- 散金: 437
- 红花: 3
- 帖子: 1038
- 在线: 74.6小时
- 虫号: 3986638
- 注册: 2015-07-23
- 性别: GG
- 专业: 高分子物理与高分子物理化
7楼2019-08-18 23:16:44
20