| 查看: 5616 | 回复: 3 | ||
[求助]
lammps metal 输出应力单位，如何计算局部应变
|
|
建了HCP多晶模型单向拉伸，感觉lammps说明里解释的不清楚，不知道 stress/atom 到底如何转化到MPa 参考别人的帖子用了以下公式（600^3为box体积，12010379为原子个数） stress=stress/(600^3./12010379.)*0.1（MPa），但得到的数据是10^5级别，这应力也太大了，请高手指定以下问题出在哪里。（运行时每拉伸100步，弛豫400步） in文件如下： units metal boundary s s s atom_style atomic read_data polycrystal.data pair_style eam/fs pair_coeff * * Ti.eam.fs Ti neighbor 2.0 bin neigh_modify delay 5 region all_atoms block INF INF INF INF INF INF region lower block INF INF INF INF INF -295 region upper block INF INF INF INF 295 INF group lower region lower group upper region upper group boundary union lower upper group mobile subtract all boundary delete_atoms porosity all_atoms 0.0005 1000 compute ack mobile ackland/atom compute new3 all temp fix 1 all nve fix 4 all temp/rescale 10 300.0 300.0 5.0 1.0 fix_modify 4 temp new3 velocity all create 300 887723 temp new3 #初始化mobile group中原子的初速度 #通过温度300K初始化，温度通过new3计算，887723为随机数 compute sv mobile stress/atom NULL virial # (xx yy zz xy xz yz) dump 5 mobile cfg 100 dump-*.cfg mass type xs ys zs c_ack c_sv[1] c_sv[2] c_sv[3] c_sv[4] c_sv[5] c_sv[6] timestep 0.001 thermo_style custom step temp pe etotal press vol thermo 100 #每100步在屏幕上打印统计热动力学量 thermo_modify temp new3 #输出的温度为new3的计算值 minimize 1.0e-6 1.0e-6 1000 1000 run 1000 #弛豫 fix 2 lower setforce NULL NULL 0.0 #使得lower group之间没有z方向的力，z方向钢化 fix 3 upper setforce NULL NULL 0.0 #使得upper group之间没有z方向的力，z方向钢化 variable i loop 10000 label loopa velocity upper set NULL NULL 0.2 # 10的8次方级别 velocity lower set NULL NULL -0.2 run 100 velocity upper set NULL NULL 0.0 velocity lower set NULL NULL 0.0 run 400 next i jump SELF loopa write_data final.data  6600.png |
回复此楼» 猜你喜欢
filecode
已经有18人回复
-
听说今天filecode变了
已经有44人回复
-
基金中了
已经有9人回复
-
固定端突然变了，今天
已经有8人回复
-
filecode
已经有11人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有6人回复
-
一个有机合成实验室都需要哪些设备？
已经有12人回复
-
求各位大神看下
已经有33人回复
-
咨询面上基金
已经有7人回复
-
影响面上的因素
已经有11人回复
2楼2016-03-14 22:29:55
huang998
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2126.9
- 散金: 427
- 红花: 1
- 帖子: 278
- 在线: 89.6小时
- 虫号: 5745669
- 注册: 2017-02-25
- 专业: 制造系统与自动化
3楼2017-09-13 15:43:41
木棉和我
木虫 (职业作家)
- 应助: 5 (幼儿园)
- 金币: 9371.5
- 散金: 2362
- 红花: 3
- 帖子: 3376
- 在线: 683.1小时
- 虫号: 1852823
- 注册: 2012-06-08
- 性别: GG
- 专业: 冶金物理化学与冶金原理
4楼2020-02-21 00:58:28
50