24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 3927  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

15764220082

新虫 (初入文坛)

不能同时用npt和变形命令作用于同一个压力矢量

发自小木虫IOS客户端
5楼2017-04-29 11:27:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

15764220082

新虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
4楼: Originally posted by 失火的天堂 at 2017-04-28 09:47:59
又出现一个问题:ERROR: Invalid thermo keyword in variable formula
能否给一点建议?...

检查一下拼写问题,和thermol命令格式

发自小木虫IOS客户端
6楼2017-04-29 11:28:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liujingqi

金虫 (小有名气)

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

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

朱雪刚521

铁杆木虫 (知名作家)

【答案】应助回帖

fix 命令也可以在单方向控制压力,fix npt 命令后边别使用iso或uniso命令,改成x 0 0 100 和y 0 0 100 的方向单独控制压力。
8楼2017-05-03 15:05:15
已阅   回复此楼   关注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的回帖

Alehmn

新虫 (著名写手)

楼主问题解决了吗,我最近也遇到相同的问题,如果能赐教就再好不过了
10楼2018-04-13 16:50:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 失火的天堂 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 274求调剂 +7 山阿蔓 2026-04-07 7/350 2026-04-13 16:44 by 张zhihao
[考研] 一志愿2110,化学学硕310分,本科重点双非求调剂 +20 努力奋斗112 2026-04-08 20/1000 2026-04-13 14:24 by 张zhihao
[考研] 一志愿哈工大 085600 277 12材科基求调剂 5+5 chenny174 2026-04-10 33/1650 2026-04-13 08:50 by Sammy2
[考研] 0854调剂 +10 长弓傲 2026-04-12 12/600 2026-04-13 08:49 by Sammy2
[考研] 343求调剂 +9 王国帅 2026-04-10 9/450 2026-04-11 20:31 by dongdian1
[考研] 求调剂 +10 璃茉一定上岸 2026-04-10 10/500 2026-04-11 13:31 by 1005715100
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 中科院总分315求调剂 +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[考研] 求调剂 +5 不会飞的鱼@ 2026-04-10 5/250 2026-04-10 19:07 by chemisry
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +20 慕绝cc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[考研] 本科西工大 0856 324求调剂 +10 wysyjs25 2026-04-09 11/550 2026-04-10 08:37 by 5268321
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 308求调剂 +17 墨墨漠 2026-04-06 17/850 2026-04-09 09:25 by 壹往無前
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 293分求调剂,外语为俄语 +7 加一一九 2026-04-07 10/500 2026-04-08 20:14 by yutian743
信息提示
请填处理意见