24小时热门版块排行榜

返回列表

沼泽留恋

金虫 (小有名气)

应助: 2 (幼儿园)
金币: 985.7
散金: 50
红花: 1
帖子: 52
在线: 51.1小时
虫号: 2635323
注册: 2013-09-04
性别: GG
专业: 金属材料的力学行为

[求助] 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

回复此楼