| 查看: 2238 | 回复: 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 |
回复此楼» 猜你喜欢
271求调剂
已经有38人回复
-
293求调剂
已经有18人回复
-
327求调剂
已经有22人回复
-
材料085601调剂
已经有33人回复
-
求调剂
已经有14人回复
-
求调剂
已经有11人回复
-
化学070300 求调剂
已经有21人回复
-
297工科调剂?
已经有3人回复
-
300分求调剂 （085501机械专硕，本科扬大）
已经有9人回复
-
调剂求收留
已经有32人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
