| 查看: 2207 | 回复: 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 |
回复此楼» 猜你喜欢
070300化学求调剂
已经有3人回复
-
311求调剂
已经有3人回复
-
生物学调剂
已经有5人回复
-
269专硕求调剂
已经有6人回复
-
08工科 320总分 求调剂
已经有8人回复
-
寻找调剂
已经有3人回复
-
291 求调剂
已经有3人回复
-
一志愿北京化工大学070300 学硕336求调剂
已经有8人回复
-
一志愿东华大学化学070300，求调剂
已经有6人回复
-
318求调剂
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
10