24小时热门版块排行榜    

查看: 2893  |  回复: 0

zyren123

木虫 (文坛精英)

[求助] vasp模拟拉伸时的应力应变曲线问题

刚开始做拉伸模拟,有些不明白的问题;
1.我模拟的是<100>方向的拉伸,设置也很简单,就是在POSCAR中第一个参数改了一下例如改为1.02:
12.1120000000000001     
     1.0200000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000    1.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.0000000000000000

INCAR文件如下:
NSW= 60
IBRION = 2
ISIF = 2
LREAL= Auto
EDIFF = 1E-3
EDIFFG = -0.1
ENCUT=350
ISTART=0;ICHARG=2
ISMEAR= 1
PREC=High
但是得到的OUTCAR中Y和Z方向的应力值中不为零啊,就是下面第二个和第三个值不为零,
in kB     -82.61797   82.61797   -82.61797     0.00000     0.00000     0.00000
请问这是怎么回事呢?是不是本来这时候Y和Z方向的值就不为零还是我的设置有问题?是不是我ISIF=2的原因?
2,另外一个问题,就是我不考虑应变即POSCAR如下
12.1120000000000001     
     1.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000    1.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000    1.0000000000000000
INCAR文件和上面一样,得到的是不是零应变下的应力值啊?这时候为什么我的应力值不为零呢?OUTCAR中如下:
  in kB     -84.40306   -84.40306   -84.40306    -0.00000     0.00000     0.00000
作为应力应变曲线的零点应力应该为零才对吧?这是什么原因呢?是不是我的结构优化不彻底呢?不彻底的话该如何优化?是不是要按照侯老师的做法做E-V曲线找出最低能量的体积呢?

[ Last edited by zyren123 on 2013-5-6 at 13:38 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP

» 猜你喜欢

当你转过身,我的幸福便与你无关
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 zyren123 的主题更新
信息提示
请填处理意见