24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3912  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 失火的天堂 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料科学与工程求调剂 +5 深V宿舍吧 2026-03-29 5/250 2026-03-29 19:10 by 学员8dgXkO
[考研] 一志愿南航 335分 | 0856 | GPA 4.07 | 有科研经历 +7 cccchenso 2026-03-29 7/350 2026-03-29 15:00 by 唐沐儿
[考研] 求调剂,一志愿 南京航空航天大学 ,080500材料科学与工程学硕,总分289分 +7 @taotao 2026-03-29 7/350 2026-03-29 12:03 by longlotian
[考研] 279求调剂 +4 蝶舞轻绕 2026-03-29 4/200 2026-03-29 09:45 by laoshidan
[考研] 本科双非材料,跨考一志愿华电085801电气,283求调剂,任何专业都可以 +6 芝士雪baoo 2026-03-28 8/400 2026-03-29 08:16 by 松花缸1201
[考研] 315求调剂 +4 akie... 2026-03-28 5/250 2026-03-28 21:05 by zhq0425
[考研] 生物学学硕,一志愿湖南大学,初试成绩338 +6 YYYYYNNNNN 2026-03-26 7/350 2026-03-28 20:52 by 唐沐儿
[考研] 086000调剂 +3 7901117076 2026-03-26 3/150 2026-03-27 21:34 by Jianing_Mi
[考研] 07化学280分求调剂 +10 722865 2026-03-23 10/500 2026-03-27 15:51 by Plutoqq
[考研] 材料292调剂 +12 橘颂思美人 2026-03-23 12/600 2026-03-27 15:44 by caszguilin
[考研] 279 分 求调剂 +4 睡个好觉_16 2026-03-24 4/200 2026-03-27 15:05 by 醉在风里
[考研] 调剂 +3 李嘉图·S·路 2026-03-27 3/150 2026-03-27 11:19 by wangjy2002
[考研] 考研调剂 +9 小蜡新笔 2026-03-26 9/450 2026-03-27 11:10 by 不吃魚的貓
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[硕博家园] 北京林业大学硕导招生广告 +6 kongweilin 2026-03-26 8/400 2026-03-27 10:18 by FF_16
[考研] 求调剂 一志愿 本科 北科大 化学 343 +6 13831862839 2026-03-24 7/350 2026-03-26 22:57 by 不吃魚的貓
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-26 6/300 2026-03-26 18:03 by 邱gl
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 328求调剂 +4 LHHL66 2026-03-23 4/200 2026-03-23 14:55 by lbsjt
信息提示
请填处理意见