| 查看: 15148 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
bjwang铁杆木虫 (正式写手)
|
[求助]
VASP计算二维体系应力中的单位问题：KB or N/M（牛/米）? 已有1人参与
|
|
|
利用改变晶格常数的方法计算应力应变的关系（stress-strain curve）：epsilon = (a-a_{0})/a_{0}； 具体计算中：设置ISIF = 2，真空层取1nm。 验证2D的石墨烯在5%的形变下的应力（文献中的计算值近似为150N/M或40GPa） 但我在OUTCAR中找到 “in kB -150.63161”，单位是KB。 1，换算成GPa 按照1KB=0.1GPa，我的结果150.63161KB=15.6GPa，与40GPa不符； 2，换算成N/M 150.63161KB=15.6GPa=15.6*10^9牛/米平方， 那么按照有关文献的说法，算二维体系的应力，考虑单位长度上的力，而不是单位面积上的，即单位要表示为牛/米， 牛/米平方要乘上真空层的厚度（这个不确定）？ 但结果也不对：15.6*10^9牛/米平方*1nm=15.6N/M，和文献中的150N/M刚好错了一个数量级。 但有个帖子（http://muchong.com/html/201406/7592923.html）说，在vasp计算石墨烯的杨氏模量中，单位就是N/M（为什么？？？） 那上面的应力的单位也应是N/M? 如果是这样，结果就和文献符合的很好，但是我上面的换算过程哪里有问题呢？ 还是计算有问题？ 谢谢各位！ 附计算结果 FORCE on cell =-STRESS in cart. coord. units (eV): Direction XX YY ZZ XY YZ ZX -------------------------------------------------------------------------------------- Alpha Z 11.84580 11.84580 11.84580 Ewald 1223.83420 1231.72516 -1753.82933 0.00000 0.00000 0.00000 Hartree 1199.50298 1201.90716 -1177.06490 0.00000 0.00000 0.00000 E(xc) -57.55665 -57.59624 -59.87916 0.00000 0.00000 0.00000 Local -2603.83155 -2612.99666 2720.52797 0.00000 0.00000 0.00000 n-local -12.23578 -13.54492 -12.15478 0.08818 0.00000 0.00000 augment -1.40512 -1.40962 -1.31894 0.00000 0.00000 0.00000 Kinetic 229.49916 238.74199 271.32154 -4.37353 0.00000 0.00000 Fock 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 ------------------------------------------------------------------------------------- Total -10.34696 -1.32733 -0.55181 0.00000 0.00000 0.00000 in kB -150.63161 -19.32335 -8.03326 0.00000 0.00000 0.00000 external pressure = -59.33 kB Pullay stress = 0.00 kB 附文献stress-strain curve图：  N每M.png  GPa.png |
回复此楼» 猜你喜欢
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
朋友圈看到的
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有110人回复
-
明天放榜？
已经有5人回复
-
时间戳变了，能看出什么问题？
已经有11人回复
-
重要消息，中午系统在维护
已经有10人回复
-
投稿咨询
已经有6人回复
-
时间戳又变了8-15
已经有26人回复
» 本主题相关商家推荐: (我也要在这里推广)
dongdong33
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 305.4
- 帖子: 232
- 在线: 81.8小时
- 虫号: 472038
- 注册: 2007-12-04
- 性别: MM
- 专业: 半导体材料
12楼2016-04-29 13:51:02
2楼2016-03-27 20:22:02
bjwang
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 5860.3
- 散金: 94
- 红花: 6
- 帖子: 792
- 在线: 556.3小时
- 虫号: 250077
- 注册: 2006-05-08
- 专业: 凝聚态物性 II ：电子结构
3楼2016-03-27 20:59:15
bjwang4楼2016-03-28 14:49:47
