| 5 | 1/1 | 返回列表 |
| 查看: 14146 | 回復: 15 | |||||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||||
bjwang鐵桿木蟲 (正式寫手)
|
[求助]
VASP計算二維體系應力中的單位問題:KB or N/M(牛/米)? 已有1人參與
|
||||||
|
利用改變晶格常數(shù)的方法計算應力應變的關系(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剛好錯了一個數(shù)量級。 但有個帖子(http://www.gaoyang168.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 |
第一性原理計算-tg | 表界面計算 | VASP計算資料收集byZF | vasp計算方法收集 |
第一性 |

銅蟲 (小有名氣)

鐵桿木蟲 (正式寫手)


| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求材料調(diào)劑 +8 | 隔壁陳先生 2026-03-12 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿吉林大學材料學碩321求調(diào)劑 +4 | Ymlll 2026-03-18 | 6/300 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 085600材料與化工調(diào)劑 324分 +8 | llllkkkhh 2026-03-18 | 8/400 |
|
|
[考研] 【同濟軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 一志愿中國海洋大學,生物學,301分,求調(diào)劑 +4 | 1孫悟空 2026-03-17 | 4/200 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調(diào)劑 +7 | 手機用戶 2026-03-17 | 7/350 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 一志愿,福州大學材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 考研調(diào)劑 +4 | 芬達46 2026-03-12 | 4/200 |
|