24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3929  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

失火的天堂

新虫 (小有名气)

[求助] 压缩问题求助!! 已有4人参与

做一个金属有机材料的压缩问题,但是老是出现问题:Cannot use fix npt and fix deform on same component of stress tensor
请问有大神能给点建议嘛?下面是in文件:
# MOF5_ReaxFF.in

units real
boundary p p p
atom_style charge
dimension 3

read_data    MOF5_Unitcell_ReaxFF.data


pair_style   reax/c NULL
pair_coeff   * * ffield.reax  Zn O C H

fix          1  all qeq/reax 1 0.0 10.0 1.0e-6 reax/c
neighbor     0.5 bin
neigh_modify  every 5  check yes
timestep     1
# define strain
variable     tmp equal "lz"
variable     L0 equal ${tmp}  
variable     strain equal "(lz - v_L0)/v_L0"
# define volume
variable     volume equal 26.98333*26.98333*26.98333

min_style    sd
minimize     1.0e-10 1.0e-10 1000 1000

velocity     all create 300 12345 mom yes rot no
fix          2  all npt temp 300 300 0.1 iso 0 0 0.1


thermo_style custom step temp etotal press
thermo       0
run          50000

fix          compression    all  deform 1000  z erate  -3.0e-6  remap x
fix          2  all npt temp 300 300  0.1 iso 0 0 0.1


compute      stress all stress/atom NULL
compute      stressx all reduce sum c_stress[1]
compute      stressy all reduce sum c_stress[2]
compute      stressz all reduce sum c_stress[3]

variable     p1  equal "(c_stressx+c_stressy+c_stressz)/(3*volume*9869)"
variable     p2  equal "c_stressx/v_volume/9869"
variable     p3  equal "c_stressy/v_volume/9869"
variable     p4  equal "c_stressz/v_volume/9869"
variable     st      equal   "step"
variable     p0      equal   "v_strain"

fix          def1 all print 1000 "${st}  ${p0}  ${p1} ${p2} ${p3} ${p4} " file [001]compression.txt screen no
dump         1 all custom  1000  dump.dat id x y z
thermo_style custom step temp etotal
thermo       1000
run          100000
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

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

liujingqi

金虫 (小有名气)

想问下楼主用reax势函数算压缩精度怎么样,我这边和实验值差太多了

发自小木虫IOS客户端
7楼2017-04-30 15:39:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

iSimuLy

捐助贵宾 (正式写手)

资深专家顾问


【答案】应助回帖

感谢参与,应助指数 +1
就是收fix deform和fix npt不能同时用

fix deform 你可以和fix nvt一起用
2楼2017-04-27 12:59:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iSimuLy

捐助贵宾 (正式写手)

资深专家顾问


【答案】应助回帖

NVT和NPT区别,你可以看看这个:http://i.youku.com/i/UMjg3MDcwOD ... m=a2hzp.8249370.0.0
3楼2017-04-27 13:01:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

失火的天堂

新虫 (小有名气)

引用回帖:
3楼: Originally posted by iSimuLy at 2017-04-27 13:01:00
NVT和NPT区别,你可以看看这个:http://i.youku.com/i/UMjg3MDcwODM1Mg==/navcustom?id=673553&spm=a2hzp.8249370.0.0

又出现一个问题:ERROR: Invalid thermo keyword in variable formula
能否给一点建议?
4楼2017-04-28 09:47:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026 WR青拔 +3 liuchb715 2026-04-09 6/300 2026-04-13 18:40 by liuchb715
[考研] 284求调剂 +10 让我上岸吧阿西 2026-04-09 10/500 2026-04-13 17:05 by 张zhihao
[考研] 调剂 +26 不逢春 2026-04-07 27/1350 2026-04-13 14:24 by 张zhihao
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 53/2650 2026-04-13 14:20 by 张zhihao
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +12 emokidd 2026-04-08 13/650 2026-04-13 08:32 by lhj2009
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 11408。358求调剂 +3 TMYzds 2026-04-07 3/150 2026-04-11 17:10 by 氮气气气
[考研] 298求调剂 +9 钉叮咚冬瓜 2026-04-07 11/550 2026-04-11 09:35 by zhq0425
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 一志愿北理工298英一数二已上岸,感谢各位老师 +14 Reframe 2026-04-10 16/800 2026-04-10 23:07 by caotw2020
[考研] 调剂 +12 卷卷卷心菜_ 2026-04-09 13/650 2026-04-10 22:36 by Ftglcn90
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 296求调剂 +6 汪!?! 2026-04-08 6/300 2026-04-10 11:02 by mattzhming
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 材料化工总分334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 20:19 by maddjdld
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
[考研] 085400电子信息类(川大控制工程)求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
信息提示
请填处理意见