24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lijuanshu

铜虫 (正式写手)

应助: 1 (幼儿园)
金币: 330.2
帖子: 360
在线: 150.7小时
虫号: 1415189

[交流] 第一性原理拉伸，怎么出现两处断裂？

我研究元素偏聚在二元合金晶界处对该体系的影响，用第一性原理拉伸试验对含有晶界的体系进行拉伸，图1是不加应变的电荷密度，图2是应变为60%的电荷密度，从图二中可以明显的看出有两处断裂，并且晶界处并没有发生断裂，断裂是发生在晶界的附近，我的结果与文献上的差别很大，还希望与了解这方面的虫友们探讨下，共同解决问题，

图1CHGCAR.gif

图2CHGCAR.gif

回复此楼

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

查看全部散金贴

1楼 2013-04-25 16:35:52

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

acridine

木虫 (正式写手)

1ST强帖: 11
应助: 13 (小学生)
金币: 4901.6
帖子: 706
在线: 389.5小时
虫号: 305478

★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+2, 谢谢回帖交流，， 2013-05-13 21:28:26

关于层厚，其实准确的做法应该是做收敛测试，就像你之前提到的那片文献中说的那样。你应该不断增加slab厚度，进行同样计算，直到最后确定再增加厚度也不会影响你的计算了，就表明厚度够了。不过你现在是在探索阶段，没必要做这么大量的繁复分析计算。不妨参考文献中类似体系，看他们采用的厚度是多少，你用类似的厚度进行探索试验应该就可以了。
关于电荷密度的分析，还是多看看文献吧，感觉你的文献量还不够，我虽然不是专门做拉伸的，只看过不多的文章，不过里面就不乏用电荷密度判断键合的。你专门研究这个，阅读的文献理应远远超过我所涉猎的。
关于应力，你的计算软件肯定有计算应力的功能，你仔细看看手册吧
关于断裂开始之后的模拟，还是建议你多看文献，这部分很复杂，影响因素很多，一定要确定你的拉伸能否模拟这个过程。其实这种应变+几何优化的方法并不是十分精密的，有很多细节需要注意。比如一些严谨的文章，在应变要达到断裂起点的时候，在每次几何优化之后都要附加一个振动计算，看有没有虚频，以判断优化后的结果是否是稳定的，这样才能精密找到断裂开始的地方。可想而知，断裂开始之后，影响因素更复杂，并不认为这种简单的应变+几何优化就能描述复杂真实的过程。