| 查看: 9168 | 回复: 31 | |||
[交流]
【求助】MS计算应力应变曲线 已有6人参与
|
|||
|
在论坛里有一些关于castep计算应力应变曲线的讨论和资料，但是基本都是无果，不知道高手在哪里，现寻高手解决一些应力应变曲线的问题。 首先是这个话题 http://muchong.com/bbs/viewthread.php?tid=1451583&fpage=4 这个帖子里面的文献中，请问是如何得到应力的？应变可以从改变晶格常数得到，可是应力是怎么得到的？是弹性常数的输出文件吗？还是在其他地方？ 还有一个帖子，石沉大海了http://muchong.com/bbs/viewthread.php?tid=1636542&fpage=4 假如初级晶胞的晶格常数（拉伸方向，设为Z方向）为10A，先进行超晶胞设置，否则晶格常数的改变要符合空间群的规则，有的晶胞改变了X方向，则Y和Z方向的可能都会跟着变，如X=Y=Z的晶体。分别取Z方向为10.1，10.2，10.3，10.4,10.5等，固定Z的大小不变，优化X和Y方向，得到拉伸时刻优化的晶胞，再进行计算（计算弹性常数），可以得到弹性模量。以上Z=10.1的时候应变就是0.01，Z=10.2的时候应变就是0.02等等。 具体大概就是这个思路了，可是我想知道的是，那个应力是在哪里看的？还是要进行转换的？请大侠指教一二！！！！ 不怕大家笑话，我是菜鸟一个，请各位高手大侠指正。可以解决怎么得到应力的给20个金币，可以讲述整个得到应力应变曲线过程的给50个金币 真诚的谢谢各位~参与讨论的也分别有金币相送 |
回复此楼» 猜你喜欢
可靠消息，8月21号(明天)放榜
已经有10人回复
-
只有每年这种时候来逛逛小木虫
已经有20人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- castep计算类金刚石结构的应力应变曲线
已经有18人回复
- 镁合金压缩应力应变曲线
已经有17人回复
- CASTEP模块中如何通过施加应变来得到应力，并形成应力应变曲线，详细实施步骤
已经有10人回复
- 拉伸应力应变曲线
已经有26人回复
- 应力应变曲线为什么是这样的？
已经有4人回复
- ANSYA中模拟时计算的应力应变曲线！！！！
已经有4人回复
- 关于工程应力应变曲线换算成真实应力应变曲线的问题...
已经有12人回复
- 【求助】CASTEP计算弹性常数得到应力应变曲线
已经有17人回复
- 【求助】如何根据应力应变曲线得到屈服强度
已经有3人回复
- 【求助】应力-应变曲线
已经有5人回复
- 【求助】为什么此拉伸应力-应变曲线如此怪异？
已经有25人回复
- 【求助】castep计算应力应变曲线结果文件中的两个stress tensor
已经有5人回复
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
icyspace(金币+10,VIP+0):谢谢讨论！ 12-2 23:33
yjmaxpayne(金币+5,VIP+0):谢谢 12-3 12:11
gavinliu7390(金币+5,VIP+0):谢谢解释！ 12-5 02:17
icyspace(金币+10,VIP+0):谢谢讨论！ 12-2 23:33
yjmaxpayne(金币+5,VIP+0):谢谢 12-3 12:11
gavinliu7390(金币+5,VIP+0):谢谢解释！ 12-5 02:17
|
xirainbow说的没错，在out文件最后可以找到这个stress值，是3*3矩阵的形式给出的。 如果没记错的话，优化晶格时结果会自动给出stress。如果只是计算能量，需要在属性选项中把stress选上 至于应力－应变曲线，计算一系列应变下的应力值就ok了 |
6楼2009-12-02 18:50:09
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
9楼2009-12-03 09:56:57
aijun0118
木虫 (正式写手)
博士
- 应助: 0 (幼儿园)
- 金币: 3686.9
- 散金: 10
- 帖子: 663
- 在线: 9.6小时
- 虫号: 506193
- 注册: 2008-02-18
- 性别: GG
- 专业: 金属材料的微观结构
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yjmaxpayne(金币+3,VIP+0):谢谢 12-3 12:11
icyspace(金币+10,VIP+0):谢谢你！！！！ 12-3 22:55
yjmaxpayne(金币+3,VIP+0):谢谢 12-3 12:11
icyspace(金币+10,VIP+0):谢谢你！！！！ 12-3 22:55
| 计算应力应变不是件简单的事情，楼主的思路是基本正确的，但是你施加应变后比如Z值改变后，然后优化xy，得到的结果中直接提取stress，便是施加应变后的应力，也即是能反应出应力应变的关系。但是你为什么还要在优化xy后计算弹性常数，从而得到应力呢。其实想得到stress不是件难事，关键是施加应变的步骤和驰豫的方式。如果你的驰豫方式对了的话，那么stress的提取只是技术的问题。咨询便可。 |
10楼2009-12-03 10:14:59
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
★ ★ ★ ★ ★
aylayl08(金币+5,VIP+0):谢谢指点 12-5 16:40
aylayl08(金币+5,VIP+0):谢谢指点 12-5 16:40
|
完全正确。谢谢你的解释。 由于发贴人是直接用ms算，这个算压力应变曲线是很难的。因为实际在应力应变曲线中，如果你想计算tensile strength的话，加入是计算a轴的，那么你需要改变a轴的长度，然后保证a轴不变，优化b和c轴，这样才是模拟晶体的实际情况。ms本身是不带有这个功能。加入你想计算特殊方向，如111方向，那你得从新定义格矢量，然后让b和c垂直你现在的a轴，即111方向的轴。 我相信发帖人应该是用ms中的castep模块，所以你实际是要改一下源代码的。 如果不想改源程序的话，我建议你使用abinit，abinit直接就有保证某一个轴不动，而优化其它的方向。 如果你想计算shear strength的话，一般的做法是给一个非对角元的一个stress，从而优化它，这样一步一步走，就可以得到。 |
15楼2009-12-05 01:36:22
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
18楼2009-12-05 02:57:01
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
qasd(金币+1,VIP+0):xiexie~ 12-7 19:54
小木虫(金币+0.5):给个红包，谢谢回帖交流
qasd(金币+1,VIP+0):xiexie~ 12-7 19:54
|
立方ＢＮ和ｄｉａｍｏｎｄ是一种结构．它们最弱的方向都是１１１方向． 做其它的方向意义不大． 所以我这里ＢＮ也是１１１方向． |
» 本帖已获得的红花（最新10朵）
donggaomin21楼2009-12-05 21:49:47
23楼2009-12-06 16:13:21
25楼2011-01-23 17:29:04
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2009-12-01 19:37:39
xirainbow
木虫 (正式写手)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 贵宾: 0.03
- 金币: 2132.6
- 红花: 2
- 帖子: 775
- 在线: 6.1小时
- 虫号: 431211
- 注册: 2007-08-11
- 专业: 光谱学
3楼2009-12-01 19:45:03
4楼2009-12-01 23:45:03
5楼2009-12-02 12:55:23
7楼2009-12-02 23:43:02
|
这位大哥，我把计算结果贴上来： 例如，Al计算的结果，初级晶胞的计算：a=b=c=2.86343A，计算其energy，将属性stress选上，结果为： ***************** Stress Tensor ***************** * * * Cartesian components (GPa) * * --------------------------------------------- * * x y z * * * * x -2.867555 0.002297 0.002297 * * y 0.002297 -2.867555 0.002297 * * z 0.002297 0.002297 -2.867555 * * * * Pressure: 2.8676 * * * ************************************************* 使c值变为原来的1.1倍，使c=3.149773，将c值保持不变，对a和b进行优化，，然后进行计算，结果为： ***************** Stress Tensor ***************** * * * Cartesian components (GPa) * * --------------------------------------------- * * x y z * * * * x 3.847851 -0.312605 0.798744 * * y -0.312605 3.847851 0.798744 * * z 0.798744 0.798744 2.826732 * * * * Pressure: -3.5075 * * * ************************************************* 那么，应变为0时候的应力应该为0,应变为0.1时的应力应该是多少呢？是要减掉应变为0时候的应力吗？ |
8楼2009-12-03 00:05:11
