| 查看: 2275 | 回复: 6 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[求助]
用侯博士计算弹性常数的方法，计算NITI B2相的C11和C12出现了问题。求助高手！
|
||||
|
我计算的Niti B2结构，用了两种方法来计算弹性常数，一种是通过VASP直接计算，IBRION=6， ISIF=3， NFREE=4 计算出来C11=194GPA，C12=178GPA，C44=53GPA。与其他参考文献对照了一下，发现C11和C12相差稍远，C44倒是有很接近的。 于是又用侯博士的那种方法，具体方法小木虫上面应该有挺多帖子。DELTA变化范围从-0.02到0.02，间隔0.0025。 第一步算C44，E=（0，0，0，DELTA，DELTA，DELTA） 得到的结果是49GPA，这与参考文献相差不多，参考值46GPA。 第二步算C11，C12。 分别用E=（DELTA，DELTA，0，0，0，0） 和 E=（DELTA，DELTA，DELTA，0，0，0） 计算（E-E0）/V0 ~DELTA，然后进行拟和，联合求解计算C11，C12，发现这两个值居然相当接近！ C11=159GPA C12=160GPA 这与文献相当不符合阿。文献里面一般是C11-C12〉0 我检查了好几遍觉得程序没什么问题，就用的侯博士的程序~~ 现在不知道该怎么办了。。。。求助各位高手阿！！！ 请版主代扣金币吧。。。我还不太熟悉怎么弄啊。。。谢谢了。。。 |
回复此楼» 猜你喜欢
材料与化工调剂
已经有26人回复
-
一志愿211，化学310分，本科重点双非，求调剂
已经有19人回复
-
求调剂288
已经有7人回复
-
298求调剂
已经有13人回复
-
281求调剂
已经有3人回复
-
265求调剂
已经有11人回复
-
301求调剂
已经有5人回复
-
343求调剂
已经有3人回复
-
272分材料子求调剂
已经有18人回复
-
273求调剂
已经有57人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助VASP计算弹性常数的问题
已经有5人回复
- VASP5.2 直接计算弹性常数
已经有52人回复
- 弹性常数计算出现诡异。。。。来看看喽
已经有5人回复
- 转帖：Cubic以及Trigonal cell弹性常数的计算
已经有21人回复
- 【求助】castep计算弹性常数问题
已经有14人回复
- 【求助】CASTEP计算弹性常数
已经有15人回复
- 【求助】非零压下弹性常数计算的初始结构
已经有9人回复
送鲜花一朵 
7楼2012-10-14 20:15:15
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+3, 专家考核, 谢谢指教！ 2012-09-07 16:37:28
losenq: 金币+5, ★★★很有帮助 2012-10-14 20:07:17
感谢参与，应助指数 +1
fzx2008: 金币+3, 专家考核, 谢谢指教！ 2012-09-07 16:37:28
losenq: 金币+5, ★★★很有帮助 2012-10-14 20:07:17
|
对于立方体系，为什么不试试M. J. Mehl的方法？先做EOS拟合，得到体模量，且体模量 B0 = (C11 + 2*C12) / 3。 然后做C11-C12的拟合，拟合方程，DeltaE = (C11-C12)*V0*d^2 其中d表示应变，对应的应变矩阵为： d 0 0 0 -d 0 0 0 d^2/(1-d^2) 然后联立体模量求解C11和C12。 其实LZ的方法中， d 0 0 0 d 0 0 0 d 的应变也是求解体模量的应变了，但是拟合的是简单的能量响应，而不是EOS，所以可能结果上会有点出入。 |
2楼2012-09-07 15:07:08
3楼2012-09-10 09:24:57
4楼2012-09-12 15:28:37
