24小时热门版块排行榜    

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

icyspace

金虫 (小有名气)

[交流] 【求助】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个金币
真诚的谢谢各位~参与讨论的也分别有金币相送
回复此楼

» 猜你喜欢

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

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

icyspace

金虫 (小有名气)

引用回帖:
Originally posted by acridine at 2009-12-2 18:50:
xirainbow说的没错,在out文件最后可以找到这个stress值,是3*3矩阵的形式给出的。
如果没记错的话,优化晶格时结果会自动给出stress。如果只是计算能量,需要在属性选项中把stress选上
至于应力-应变曲线,计 ...

这位大哥,我把计算结果贴上来:


例如,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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 32 个回答

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★ ★ ★ ★ ★ ★ ★ ★
zxzj05(金币+1,VIP+0):讨论奖励! 12-1 21:47
icyspace(金币+2,VIP+0):谢谢讨论! 12-1 23:43
icyspace(金币+6,VIP+0): 12-2 23:33
想计算应力应变曲线其实不那么简单,ms 好像很难直接计算出来。
真理是一点点接近的!
2楼2009-12-01 19:37:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★
zxzj05(金币+1,VIP+0):讨论奖励! 12-1 21:47
icyspace(金币+2,VIP+0):谢谢讨论! 12-1 23:44
icyspace(金币+6,VIP+0):谢谢! 12-2 23:33
输出中有个stress
应该就是应力吧;)
3楼2009-12-01 19:45:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

icyspace

金虫 (小有名气)

引用回帖:
Originally posted by gavinliu7390 at 2009-12-1 19:37:
想计算应力应变曲线其实不那么简单,ms 好像很难直接计算出来。

啊。。。。。原来这样啊!
4楼2009-12-01 23:45:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见