【求助】MS计算应力应变曲线 - 第一性原理 - MS - 第3页小木虫论坛-学术科研互动平台

» 本帖已获得的红花（最新10朵）

donggaomin

真理是一点点接近的！

21楼2009-12-05 21:49:47

icyspace

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 794
帖子: 164
在线: 45.3小时
虫号: 529026
注册: 2008-03-19
专业: 金属结构材料

引用回帖:

Originally posted by gavinliu7390 at 2009-12-5 21:49:

立方ＢＮ和ｄｉａｍｏｎｄ是一种结构．它们最弱的方向都是１１１方向．
做其它的方向意义不大．
所以我这里ＢＮ也是１１１方向．

版主大人，我想问你哈~你计算这个的时候，会不会出现反常，比如说某个应变的时候应力稍微比总的趋势小或者大的情况？遇到这种情况应该再优化晶胞或者精度更大再计算吗？
还有一个问题哈，版主大人~你比如将C轴进行增大之后，a和b是从什么数值开始优化的？比如说Al的一开始的a=b=c=4.09636，将c变为1.1倍，那么a和b开始从什么数值开始优化的？
不好意思版主大人，没有金币给你加了，谢谢你！！！

22楼2009-12-06 16:02:56

icyspace

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 794
帖子: 164
在线: 45.3小时
虫号: 529026
注册: 2008-03-19
专业: 金属结构材料

引用回帖:

还有，老大，为何我计算的Al在应变0.01和0.02的时候得到的应力是负值哈？
其他情况虽然不好，可都是正的。。。。。老大，我是本科大四，毕业设计稍微用到计算应力应变曲线的东西，所以请你解答，谢谢你！

赞一下(1人)

23楼2009-12-06 16:13:21

gavinliu7390

木虫 (著名写手)

叶落鹰飞

1ST强帖: 3
应助: 51 (初中生)
贵宾: 0.481
金币: 1564.7
散金: 2975
红花: 34
帖子: 1278
在线: 241.4小时
虫号: 675997
注册: 2008-12-17
性别: GG
专业: 凝聚态物性I:结构、力学和

★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuchenwf(金币+2,VIP+0):谢谢 12-10 20:46

a=b=c=4.03696,把c变为1.1倍的时候，去优化a和b，当然是从4.09636开始优化。
最大的可能还是你的方法不对！

赞一下(11人)

真理是一点点接近的！

24楼2009-12-10 15:02:52

匿名

应助: (幼儿园)
在线:
虫号: 0

★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
mazuju028(金币+1): 谢谢交流 2011-01-23 20:39:53

本帖仅楼主可见

25楼2011-01-23 17:29:04

已阅申请1ST强帖回复此楼编辑查看我的主页

匿名

应助: (幼儿园)
在线:
虫号: 0

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

本帖仅楼主可见

26楼2011-01-23 21:58:55

已阅申请1ST强帖回复此楼编辑查看我的主页

donggaomin

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 770.6
散金: 805
红花: 9
帖子: 564
在线: 242.7小时
虫号: 853446
注册: 2009-09-21
性别: GG
专业: 凝聚态物性 II ：电子结构

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

送鲜花一朵

引用回帖:

17189714楼: Originally posted by gavinliu7390 at 2009-12-05 21:49:47:
立方ＢＮ和ｄｉａｍｏｎｄ是一种结构．它们最弱的方向都是１１１方向．
做其它的方向意义不大．
所以我这里ＢＮ也是１１１方向．

特向您请教一个问题：计算弹性常数时，好像是在e1,e2,e3等方向加应变，然后计算应力，利用相应的关系得到cij。e1,e2,e3分别是笛卡尔坐标的xyz，问题是，xyz如何与基矢对应起来？一个人计算时x轴对应基矢a1,另一个人计算时x轴对应基矢a2，那么计算出来的c11就不同了。没有可比性。有没有一个统一的标准呢？

我在这个帖子中发了求助，您可以去看看详细的http://muchong.com/bbs/viewthread.php?tid=4472995

祝好！

好好学习

27楼2012-05-06 11:19:44

dq20081935

铜虫 (小有名气)

应助: 7 (幼儿园)
金币: 262.8
散金: 7
红花: 1
帖子: 161
在线: 58.3小时
虫号: 1625137
注册: 2012-02-18
性别: GG
专业: 金属功能材料

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

17177350楼: Originally posted by acridine at 2009-12-05 02:57:01:
楼上同学应该是专门做超硬的，嘿嘿，解释很专业。
不过需要更正一下，据我所知ms中的castep可以实现固定某个轴同时优化另外的轴，不过当前我这里没有ms软件，无法给出具体的步骤。我记得应该是某个菜单中的modify-

请问stress值在哪给出呢？Discover或者castep模块？谢谢！！！

28楼2012-05-24 11:19:36

1026580419

木虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3740
帖子: 14
在线: 9.5小时
虫号: 1246672
注册: 2011-03-27
专业: 金属非晶态、准晶和纳米晶

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

17144043楼: Originally posted by aijun0118 at 2009-12-03 10:14:59
计算应力应变不是件简单的事情，楼主的思路是基本正确的，但是你施加应变后比如Z值改变后，然后优化xy，得到的结果中直接提取stress，便是施加应变后的应力，也即是能反应出应力应变的关系。但是你为什么还要在优化 ...

请教一下，stress到底该怎么提取呢？是输出文件中直接给出吗还是以其他的形式？

29楼2012-05-25 17:30:26

1026580419

木虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3740
帖子: 14
在线: 9.5小时
虫号: 1246672
注册: 2011-03-27
专业: 金属非晶态、准晶和纳米晶

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

17139956楼: Originally posted by icyspace at 2009-12-03 00:05:11
这位大哥，我把计算结果贴上来：

例如，Al计算的结果，初级晶胞的计算：a=b=c=2.86343A，计算其energy，将属性stress选上，结果为：
***************** Stress Tensor *****************
* ...

3*3矩阵给出的9个数与pressure的值是什么关系啊？

30楼2012-05-25 17:33:05