| 查看: 4271 | 回复: 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 |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
Alehmn
新虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 4959.1
- 散金: 437
- 红花: 3
- 帖子: 1038
- 在线: 74.6小时
- 虫号: 3986638
- 注册: 2015-07-23
- 性别: GG
- 专业: 高分子物理与高分子物理化
11楼2018-04-20 22:10:06
2楼2017-04-27 12:59:51
3楼2017-04-27 13:01:00
4楼2017-04-28 09:47:59
