24小时热门版块排行榜    

查看: 701  |  回复: 2

失火的天堂

新虫 (小有名气)

[交流] 输出原子应力疑问 已有1人参与

输出原子应力:compute 1 all stress/atom NULL,然后需要对计算的量除以原子体积,请问dump怎么输出,试了一下不能用v_这种模式
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shiling1813

新虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
好像dump命令可以,这是我写的一段命令,供参考。
compute            spb all stress/atom thermo_temp

compute     ppptem all reduce sum c_spb[1] c_spb[2] c_spb[3] c_spb[4] c_spb[5] c_spb[6]
variable    d equal (c_ppptem[1])

thermo_style custom step temp pe etotal v_d vol ylo yhi

fix fp1  all print 10000 "$d" file pxxtem.dat

dump mydump2        all custom     10000 noreallrelaxstressalltem.txt type x y z c_spb[1] c_spb[2] c_spb[3] c_spb[4] c_spb[5] c_spb[6]
dump            all2 all custom 10000 zuobiaotem.* type x y z c_spb[1] c_spb[2] c_spb[3] c_spb[4] c_spb[5] c_spb[6]
2楼2017-09-26 10:24:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

失火的天堂

新虫 (小有名气)

引用回帖:
2楼: Originally posted by shiling1813 at 2017-09-26 10:24:16
好像dump命令可以,这是我写的一段命令,供参考。
compute            spb all stress/atom thermo_temp

compute     ppptem all reduce sum c_spb c_spb c_spb c_spb c_spb c_spb
variable    d equal (c_ppptem)

...

已解决,我要的是每个原子的应力,不过你输出的c_spb[1]之类的还需要除原子体积才是每个原子的应力,应该是variable 1 atom c_1[1]/(v_v*10000)
3楼2017-10-02 09:52:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 失火的天堂 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:29 by cLOojpOxkX0G
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:16 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 02:05 by cLOojpOxkX0G
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-26 02:03 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 5/250 2026-10-26 01:50 by cLOojpOxkX0G
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:39 by cLOojpOxkX0G
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 5/250 2026-10-26 01:26 by cLOojpOxkX0G
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-26 01:24 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 7xHjyB3bAiwg 2026-10-24 7/350 2026-10-26 01:11 by cLOojpOxkX0G
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-26 01:07 by cLOojpOxkX0G
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:54 by cLOojpOxkX0G
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 7xHjyB3bAiwg 2026-10-23 12/600 2026-10-26 00:45 by cLOojpOxkX0G
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:41 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 7xHjyB3bAiwg 2026-10-23 5/250 2026-10-26 00:19 by cLOojpOxkX0G
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 GVGRDgNfv7P4 2026-10-23 5/250 2026-10-26 00:08 by cLOojpOxkX0G
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 WqtG0Aem8Ufl 2026-10-24 3/150 2026-10-25 20:11 by cLOojpOxkX0G
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 3/150 2026-10-25 14:23 by XiYP04It3Evb
[文学芳草园] 亲密关系 +7 宋昙 2026-10-21 13/650 2026-10-23 17:07 by 小牧11
[考博] 27秋申博自荐 +4 且听风吟时 2026-10-21 4/200 2026-10-23 13:26 by 大大大肥虫
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
请填处理意见