24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学、生物学、生物与医药接受调剂
查看: 3926  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 失火的天堂 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856专硕求调剂 希望是a区院校 +23 好好休息好不好 2026-04-09 26/1300 2026-04-13 14:22 by 张zhihao
[考研] 一志愿双非085400电子信息344 求调剂,对材料和化学方向也感兴趣 +12 无情的小羊 2026-04-09 13/650 2026-04-13 14:17 by 张zhihao
[考研] 考研二轮调剂 +10 故人?? 2026-04-09 10/500 2026-04-13 09:55 by szhize
[考研] 一志愿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 逆水乘风
[考研] 调剂 +10 月@163.com 2026-04-11 10/500 2026-04-12 09:14 by zhouyuwinner
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 求调剂 +11 月@163.com 2026-04-07 13/650 2026-04-11 22:55 by BruceLiu320
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 269求调剂 +11 啊啊我我 2026-04-07 11/550 2026-04-11 16:45 by vgtyfty
[考研] 机械专硕270求调剂,接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 344求调剂 +7 丶风雪夜归人丶 2026-04-09 7/350 2026-04-10 12:05 by pengliang8036
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 一志愿武理车辆 281 求调剂 +5 上岸研究生. 2026-04-07 5/250 2026-04-09 15:56 by only周
[考研] 311求调剂 +6 surte 2026-04-08 13/650 2026-04-09 14:00 by surte
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
信息提示
请填处理意见