应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 体模量的问题
111/2返回列表12下一页
查看: 2117  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

牛振威

金虫 (小有名气)

[求助] 体模量的问题

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-06-03 12:32:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

牛振威

金虫 (小有名气)
???????:
4?: Originally posted by souledge at 2012-06-03 18:51:39
????c/a?????????????????淨~
????????????????????????????
?????????Stress - Volume?????????????????в????????????????????п??????pulay stress??????~
????????????????? ...

???常感谢，我正在试你说的第二个，在平衡体积附近改???体积，固定??个体积???优化，看看结果会??会好点。???外我想知???怎么加金??啊
赞一下 回复此楼
5楼2012-06-03 19:11:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwceng0816

木虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +1
fzx2008: 金币+1, 鼓励交流 2012-06-04 08:20:01
首先拟合可能有问题  不同版本的软件会有所不同 你在试试看
赞一下 回复此楼
6楼2012-06-03 20:15:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wxl1999
7楼2012-06-03 20:36:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

牛振威

金虫 (小有名气)
引用回帖:
6楼: Originally posted by hwceng0816 at 2012-06-03 20:15:26
首先拟合可能有问题  不同版本的软件会有所不同 你在试试看

为什么不同版本，会不一样呢？我觉着不同版本应该相差不大，除非是不同的EOS。
赞一下 回复此楼
8楼2012-06-03 23:08:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hwceng0816

木虫 (著名写手)
★
fzx2008: 金币+1, 鼓励交流 2012-06-04 08:20:13
还有一种办法就是利用不同压力下的体积进行拟合 这个效果很好
赞一下 回复此楼
9楼2012-06-04 08:03:38
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 牛皮糖耶 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 4/200 2026-11-07 16:28 by Xc2FxHsHkBNc
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 05:32 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
[基金申请] 各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
关闭
请填处理意见
关闭