24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1696  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

牛振威

金虫 (小有名气)

[求助] 体模量的问题

用弹性常数得到的体模量和用EOS得出的结果差距很大,会是那里出错了?
回复此楼

» 猜你喜欢

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

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

牛振威

金虫 (小有名气)

引用回帖:
9楼: Originally posted by hwceng0816 at 2012-06-04 08:03:38
还有一种办法就是利用不同压力下的体积进行拟合 这个效果很好

这个和变体积是一样的把,改变体积就应该等于改变了压强吧。我试试
10楼2012-06-04 09:59:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

souledge

专家顾问 (著名写手)

提高计算精度。还有做弹性计算前先要大精度的优化~
当然,做EOS也要精度尽量高~
LZ问的笼统,俺就只好答得笼统了~
思想重于技巧,内涵重于表象
2楼2012-06-03 18:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛振威

金虫 (小有名气)

引用回帖:
2楼: Originally posted by souledge at 2012-06-03 18:25:18
提高计算精度。还有做弹性计算前先要大精度的优化~
当然,做EOS也要精度尽量高~
LZ问的笼统,俺就只好答得笼统了~

用MS中CASTEP算的弹性常数取的最高精度ultra-fine,截断取620eV,K点我从4*4*2一直变到7*7*7弹性常数基本不变,(都经过严格测试的)。E-V曲线也是由CASTEP中计算的,计算中固定的c/a. EOS用的是三阶B-M方程,我的是六角结构。EOS结果B=145GPa,而用弹性常数算的为110GPa. 我正在重新算E-V,看看是不是这里的原因。谢谢你,我不知道是不是还有其他什么原因会导致这种结果。(想给你加金币但是我不知道怎么弄,希望版主帮忙加下10个)
3楼2012-06-03 18:46:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fzx2008: 金币+3, 专家考核, 谢谢指教 2012-06-03 19:32:12
贺仪: 金币+20, 应楼主要求给予20BB奖励 2012-06-03 20:24:08
引用回帖:
3楼: Originally posted by 牛振威 at 2012-06-03 18:46:31
用MS中CASTEP算的弹性常数取的最高精度ultra-fine,截断取620eV,K点我从4*4*2一直变到7*7*7弹性常数基本不变,(都经过严格测试的)。E-V曲线也是由CASTEP中计算的,计算中固定的c/a. EOS用的是三阶B-M方程,我的是 ...

固定住c/a本来就是很不合理的玩法~
对于六角,至少可以选择以下两种:
一,简单的,用Stress - Volume的关系做优化,得到一系列不同体积的总能。这个方法有可能会有pulay stress的问题~
二,复杂点的,离散的取几个在平衡c/a周围的值,然后每个c/a都做一套E-V曲线,取最每个c/a的最低点连接,然后用这些点拟合EOS。
如果假设LZ的收敛测试是正确的,那么我觉得弹性常数计算的110GPa更可信一点。
思想重于技巧,内涵重于表象
4楼2012-06-03 18:51:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见