24小时热门版块排行榜    

查看: 1783  |  回复: 2

yindeqiang

银虫 (正式写手)

[求助] 第一性原理模拟单向拉伸的疑问

各位虫友,看了论坛里的一些帖子,计划对材料做第一性原理的单向拉伸的模拟。我用的是CASTEP模块(模型比如是:7层Al2O3表面,晶胞在z方向上有10埃的真空层),有几个疑问,如下:
1). 看到文献中许多拉伸都是逐渐增加拉伸方向的尺寸(比如z),即:每一步都进行优化,在优化后的模型基础上施加新的应变,并再次进行优化,直到材料断裂(或应力下降)。
  在每一步优化的过程中是否选择优化另外两个方向上晶格的尺寸(比如x,y)?
  从读过的文献看,两种都有,选择优化x,y,计算非常耗时
2). 如若选择不优化x,y(即几何优化时不选择优化晶胞,每一步z是固定的),那么要得到应力-应变曲线中的应力如何确定?
   用castep设置task :energy+ properties: stress,计算每一步优化后的模型,在结果文件.castep中的应力矩阵: stress tensor:  Cartesian components(GPa), 是否只读取zz对应的应力即可?
3). 单向拉伸是逐渐增加晶格长度,在优化的过程中是否固定最底层的原子?
  若不固定,当晶胞施加应变优化时,底层原子会向上移动
4).单向拉伸的应变如何确定,是等于晶格的应变还是取相对于上下两层原子间距离而产生的应变。
   这个问题和第三问有点类似,若固定底层原子,拉伸方向的应变就会等于材料的应变。
静候虫友的指教,欢迎讨论
回复此楼

» 收录本帖的淘帖专辑推荐

MaterialStudio学习使用

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

分享达到共赢
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖


yindeqiang(金币+20): 多谢指教,原来我是犯了基本的错误。如果想做z方向的拉伸,把真空层去掉可以不? 2011-09-20 21:12:30
uuv2010(金币+1): 多谢提示! 2011-09-20 23:26:18
就一点,你的计算原胞里面沿着z方向有真空层,但又要做沿着z方向拉伸的计算,这样是不可行的。当z方向有真空层时,沿着z方向的应力并不具有很明确的物理意义,此时能做只是沿着x,y或xy方向上的拉伸。
2楼2011-09-20 20:54:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangych0801

木虫 (正式写手)

楼主,最后选用的是优化xy方向的吗?

发自小木虫IOS客户端
要门儿活,要门儿死,绝不半死不活。
3楼2017-04-26 14:20:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yindeqiang 的主题更新
信息提示
请填处理意见