24小时热门版块排行榜    

查看: 4130  |  回复: 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的回帖

失火的天堂

新虫 (小有名气)

引用回帖:
7楼: Originally posted by liujingqi at 2017-04-30 15:39:38
想问下楼主用reax势函数算压缩精度怎么样,我这边和实验值差太多了

精度确实不高,我想应该是势函数的原因,我算出来的应力应变曲线震荡也比较明显
9楼2017-06-12 22:02:45
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +8 majunge000 2026-08-11 10/500 2026-08-15 18:55 by Equinoxhua
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见