| 查看: 5285 | 回复: 5 | ||
[求助]
剪切应力应变曲线怎么输出
|
|
根据lammps自带的in.shear剪切模拟，进行模拟后，怎么输出应力应变曲线啊，或者说剪切应力应变曲线在in文件中写的啊，求各位知道，谢谢！ # 3d metal shear simulation units metal boundary s s p atom_style atomic lattice fcc 3.52 region box block 0 16.0 0 10.0 0 2.828427 create_box 3 box lattice fcc 3.52 orient x 1 0 0 orient y 0 1 1 orient z 0 -1 1 origin 0.5 0 0 create_atoms 1 box pair_style eam pair_coeff * * Ni_u3.eam neighbor 0.3 bin neigh_modify delay 5 region lower block INF INF INF 0.9 INF INF region upper block INF INF 6.1 INF INF INF group lower region lower group upper region upper group boundary union lower upper group mobile subtract all boundary set group lower type 2 set group upper type 3 # void #region void cylinder z 8 5 2.5 INF INF #delete_atoms region void # temp controllers compute new3d mobile temp compute new2d mobile temp/partial 0 1 1 # equilibrate velocity mobile create 300.0 5812775 temp new3d fix 1 all nve fix 2 boundary setforce 0.0 0.0 0.0 fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0 fix_modify 3 temp new3d thermo 25 thermo_modify temp new3d timestep 0.001 run 100 # shear velocity upper set 1.0 0 0 velocity mobile ramp vx 0.0 1.0 y 1.4 8.6 sum yes unfix 3 fix 3 mobile temp/rescale 10 300.0 300.0 10.0 1.0 fix_modify 3 temp new2d #dump 1 all atom 100 dump.shear #dump 2 all image 100 image.*.jpg type type & # axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0 #dump_modify 2 pad 4 #dump 3 all movie 100 movie.mpg type type & # axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0 #dump_modify 3 pad 4 thermo 100 thermo_modify temp new2d reset_timestep 0 run 3000 |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有8人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
luyaobao
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2457.9
- 散金: 1091
- 红花: 23
- 沙发: 8
- 帖子: 1533
- 在线: 150.7小时
- 虫号: 7847640
- 注册: 2018-01-30
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2018-12-28 12:44:40
ruibingqilin
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 682.7
- 帖子: 285
- 在线: 152.6小时
- 虫号: 1943683
- 注册: 2012-08-17
- 性别: GG
- 专业: 高分子材料结构与性能
3楼2019-02-12 11:53:33
4楼2019-03-05 11:26:15
shanks2017
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 516.7
- 散金: 15
- 红花: 1
- 帖子: 181
- 在线: 135.2小时
- 虫号: 8723026
- 注册: 2018-05-11
- 专业: 固体力学
5楼2019-04-12 09:53:09
6楼2019-09-23 21:43:06