24小时热门版块排行榜    

查看: 5305  |  回复: 1

gaods101@163

金虫 (小有名气)


[交流] 【求助】vasp计算应力

哪位知道vasp中如何在z轴方向加应变,x轴y轴都固定,怎么计算体系的应力?谢谢!!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nll1988

新虫 (小有名气)


应力计算

★ ★ ★
sunyang1988(金币+3): 谢谢交流 2011-03-29 23:15:56
gaods101@163(金币+30): 2011-04-01 20:57:43
(1)x,y固定,z是自由的,可以通过POSCAR进行设置,POSCAR如图所示(仅作参考):
Cubic SiC
3.57
0.0    0.5     0.5
0.5    0.0     0.5
0.5    0.5     0.0
1  1
Selective dynamic
Direct  
0.000   0.000    0.000   F  F  F
0.248   0.248    0.248  T  T  T


其中F是fixed(固定),T是不固定进行优化的意思,所以你应该将x,y,z分别对应为F,F,T即可。
(2)你要计算应力,这个问题有些宽泛,
      如果你是计算应力应变的关系,画出曲线,那就应该设置PSTRESS,同时设置ISIF=4,此时代表的是:计算力,计算应力张量,弛豫离子,改变原包形状(因为有应变所以形状会改变),但是原包体积不变,否则难以进行计算,变量太多。
       如果你只是计算应力而不作应变,那就将ISIF设置为默认值2,即原包体积和形状都不变。
    在INCAR中设置了PSTRESS(根据具体晶体)后,就要从OUTCAR中的external pressure中观察。它代表外部压强,单位KB,1KB=0.1GPa,1GPa=1E9Pa,其为正表示对材料加压,为负表示对材料进行拉伸。

我就知道这么多了,具体细节还是要靠摸索从经验中得出,有错误的地方请指教。
2楼2011-03-29 22:59:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gaods101@163 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见