 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 8493 | 回復(fù): 31 | |||
icyspace金蟲 (小有名氣)
|
[交流]
【求助】MS計(jì)算應(yīng)力應(yīng)變曲線 已有6人參與
|
||
|
在論壇里有一些關(guān)于castep計(jì)算應(yīng)力應(yīng)變曲線的討論和資料，但是基本都是無(wú)果，不知道高手在哪里，現(xiàn)尋高手解決一些應(yīng)力應(yīng)變曲線的問(wèn)題。 首先是這個(gè)話題 http://www.gaoyang168.com/bbs/viewthread.php?tid=1451583&fpage=4 這個(gè)帖子里面的文獻(xiàn)中，請(qǐng)問(wèn)是如何得到應(yīng)力的？應(yīng)變可以從改變晶格常數(shù)得到，可是應(yīng)力是怎么得到的？是彈性常數(shù)的輸出文件嗎？還是在其他地方？ 還有一個(gè)帖子，石沉大海了http://www.gaoyang168.com/bbs/viewthread.php?tid=1636542&fpage=4 假如初級(jí)晶胞的晶格常數(shù)（拉伸方向，設(shè)為Z方向）為10A，先進(jìn)行超晶胞設(shè)置，否則晶格常數(shù)的改變要符合空間群的規(guī)則，有的晶胞改變了X方向，則Y和Z方向的可能都會(huì)跟著變，如X=Y=Z的晶體。分別取Z方向?yàn)?0.1，10.2，10.3，10.4,10.5等，固定Z的大小不變，優(yōu)化X和Y方向，得到拉伸時(shí)刻優(yōu)化的晶胞，再進(jìn)行計(jì)算（計(jì)算彈性常數(shù)），可以得到彈性模量。以上Z=10.1的時(shí)候應(yīng)變就是0.01，Z=10.2的時(shí)候應(yīng)變就是0.02等等。 具體大概就是這個(gè)思路了，可是我想知道的是，那個(gè)應(yīng)力是在哪里看的？還是要進(jìn)行轉(zhuǎn)換的？請(qǐng)大俠指教一二�。。�！ 不怕大家笑話，我是菜鳥一個(gè)，請(qǐng)各位高手大俠指正�？梢越鉀Q怎么得到應(yīng)力的給20個(gè)金幣，可以講述整個(gè)得到應(yīng)力應(yīng)變曲線過(guò)程的給50個(gè)金幣 真誠(chéng)的謝謝各位~參與討論的也分別有金幣相送 |
回復(fù)此樓» 猜你喜歡
【求調(diào)劑】085601材料工程專碩 | 總分272 |
已經(jīng)有5人回復(fù)
-
275求調(diào)劑
已經(jīng)有13人回復(fù)
-
0856，材料與化工321分求調(diào)劑
已經(jīng)有8人回復(fù)
-
085701環(huán)境工程，267求調(diào)劑
已經(jīng)有12人回復(fù)
-
320分，材料與化工專業(yè)，求調(diào)劑
已經(jīng)有7人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
295材料工程專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
265求調(diào)劑
已經(jīng)有8人回復(fù)
-
322求調(diào)劑
已經(jīng)有4人回復(fù)
-
085602 307分 求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- castep計(jì)算類金剛石結(jié)構(gòu)的應(yīng)力應(yīng)變曲線
已經(jīng)有18人回復(fù)
- 鎂合金壓縮應(yīng)力應(yīng)變曲線
已經(jīng)有17人回復(fù)
- CASTEP模塊中如何通過(guò)施加應(yīng)變來(lái)得到應(yīng)力，并形成應(yīng)力應(yīng)變曲線，詳細(xì)實(shí)施步驟
已經(jīng)有10人回復(fù)
- 拉伸應(yīng)力應(yīng)變曲線
已經(jīng)有26人回復(fù)
- 應(yīng)力應(yīng)變曲線為什么是這樣的？
已經(jīng)有4人回復(fù)
- ANSYA中模擬時(shí)計(jì)算的應(yīng)力應(yīng)變曲線�。。�！
已經(jīng)有4人回復(fù)
- 關(guān)于工程應(yīng)力應(yīng)變曲線換算成真實(shí)應(yīng)力應(yīng)變曲線的問(wèn)題...
已經(jīng)有12人回復(fù)
- 【求助】CASTEP計(jì)算彈性常數(shù)得到應(yīng)力應(yīng)變曲線
已經(jīng)有17人回復(fù)
- 【求助】如何根據(jù)應(yīng)力應(yīng)變曲線得到屈服強(qiáng)度
已經(jīng)有3人回復(fù)
- 【求助】應(yīng)力-應(yīng)變曲線
已經(jīng)有5人回復(fù)
- 【求助】為什么此拉伸應(yīng)力-應(yīng)變曲線如此怪異？
已經(jīng)有25人回復(fù)
- 【求助】castep計(jì)算應(yīng)力應(yīng)變曲線結(jié)果文件中的兩個(gè)stress tensor
已經(jīng)有5人回復(fù)
acridine
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 11
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 4901.6
- 紅花: 13
- 帖子: 706
- 在線: 389.5小時(shí)
- 蟲號(hào): 305478
- 注冊(cè): 2006-12-10
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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說(shuō)的沒(méi)錯(cuò)，在out文件最后可以找到這個(gè)stress值，是3*3矩陣的形式給出的。 如果沒(méi)記錯(cuò)的話，優(yōu)化晶格時(shí)結(jié)果會(huì)自動(dòng)給出stress。如果只是計(jì)算能量，需要在屬性選項(xiàng)中把stress選上 至于應(yīng)力－應(yīng)變曲線，計(jì)算一系列應(yīng)變下的應(yīng)力值就ok了 |
6樓2009-12-02 18:50:09
xirainbow
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 2
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 0.03
- 金幣: 2132.6
- 紅花: 2
- 帖子: 775
- 在線: 6.1小時(shí)
- 蟲號(hào): 431211
- 注冊(cè): 2007-08-11
- 專業(yè): 光譜學(xué)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yjmaxpayne(金幣+3,VIP+0):謝謝 12-3 12:11
icyspace(金幣+10,VIP+0):謝謝！~�。�！ 12-3 22:56
yjmaxpayne(金幣+3,VIP+0):謝謝 12-3 12:11
icyspace(金幣+10,VIP+0):謝謝！~�。�！ 12-3 22:56
|
“那么，應(yīng)變?yōu)?時(shí)候的應(yīng)力應(yīng)該為0” 所以你的第一個(gè)結(jié)構(gòu)不是穩(wěn)定結(jié)構(gòu)，就是說(shuō)他的應(yīng)力不為0 所以，你需要重新優(yōu)化幾何結(jié)構(gòu)，優(yōu)化時(shí)，提高stress的收斂標(biāo)準(zhǔn)，這樣“應(yīng)變?yōu)?時(shí)候的應(yīng)力就0” 于是，就不會(huì)出現(xiàn)“是要減掉應(yīng)變?yōu)?時(shí)候的應(yīng)力嗎？”的問(wèn)題了；） |
9樓2009-12-03 09:56:57
aijun0118
木蟲 (正式寫手)
博士
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3686.9
- 散金: 10
- 帖子: 663
- 在線: 9.6小時(shí)
- 蟲號(hào): 506193
- 注冊(cè): 2008-02-18
- 性別: GG
- 專業(yè): 金屬材料的微觀結(jié)構(gòu)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
| 計(jì)算應(yīng)力應(yīng)變不是件簡(jiǎn)單的事情，樓主的思路是基本正確的，但是你施加應(yīng)變后比如Z值改變后，然后優(yōu)化xy，得到的結(jié)果中直接提取stress，便是施加應(yīng)變后的應(yīng)力，也即是能反應(yīng)出應(yīng)力應(yīng)變的關(guān)系。但是你為什么還要在優(yōu)化xy后計(jì)算彈性常數(shù)，從而得到應(yīng)力呢。其實(shí)想得到stress不是件難事，關(guān)鍵是施加應(yīng)變的步驟和馳豫的方式。如果你的馳豫方式對(duì)了的話，那么stress的提取只是技術(shù)的問(wèn)題。咨詢便可。 |
10樓2009-12-03 10:14:59
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時(shí)
- 蟲號(hào): 675997
- 注冊(cè): 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
★ ★ ★ ★ ★
aylayl08(金幣+5,VIP+0):謝謝指點(diǎn) 12-5 16:40
aylayl08(金幣+5,VIP+0):謝謝指點(diǎn) 12-5 16:40
|
完全正確。謝謝你的解釋。 由于發(fā)貼人是直接用ms算，這個(gè)算壓力應(yīng)變曲線是很難的。因?yàn)閷?shí)際在應(yīng)力應(yīng)變曲線中，如果你想計(jì)算tensile strength的話，加入是計(jì)算a軸的，那么你需要改變a軸的長(zhǎng)度，然后保證a軸不變，優(yōu)化b和c軸，這樣才是模擬晶體的實(shí)際情況。ms本身是不帶有這個(gè)功能。加入你想計(jì)算特殊方向，如111方向，那你得從新定義格矢量，然后讓b和c垂直你現(xiàn)在的a軸，即111方向的軸。 我相信發(fā)帖人應(yīng)該是用ms中的castep模塊，所以你實(shí)際是要改一下源代碼的。 如果不想改源程序的話，我建議你使用abinit，abinit直接就有保證某一個(gè)軸不動(dòng)，而優(yōu)化其它的方向。 如果你想計(jì)算shear strength的話，一般的做法是給一個(gè)非對(duì)角元的一個(gè)stress，從而優(yōu)化它，這樣一步一步走，就可以得到。 |
15樓2009-12-05 01:36:22
acridine
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 11
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 4901.6
- 紅花: 13
- 帖子: 706
- 在線: 389.5小時(shí)
- 蟲號(hào): 305478
- 注冊(cè): 2006-12-10
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
18樓2009-12-05 02:57:01
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時(shí)
- 蟲號(hào): 675997
- 注冊(cè): 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
qasd(金幣+1,VIP+0):xiexie~ 12-7 19:54
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
qasd(金幣+1,VIP+0):xiexie~ 12-7 19:54
|
立方ＢＮ和ｄｉａｍｏｎｄ是一種結(jié)構(gòu)．它們最弱的方向都是１１１方向． 做其它的方向意義不大． 所以我這里ＢＮ也是１１１方向． |
» 本帖已獲得的紅花（最新10朵）
donggaomin21樓2009-12-05 21:49:47
icyspace
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 794
- 帖子: 164
- 在線: 45.3小時(shí)
- 蟲號(hào): 529026
- 注冊(cè): 2008-03-19
- 專業(yè): 金屬結(jié)構(gòu)材料
23樓2009-12-06 16:13:21
25樓2011-01-23 17:29:04
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強(qiáng)帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時(shí)
- 蟲號(hào): 675997
- 注冊(cè): 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
2樓2009-12-01 19:37:39
xirainbow
木蟲 (正式寫手)
- 1ST強(qiáng)帖: 2
- 應(yīng)助: 0 (幼兒園)
- 貴賓: 0.03
- 金幣: 2132.6
- 紅花: 2
- 帖子: 775
- 在線: 6.1小時(shí)
- 蟲號(hào): 431211
- 注冊(cè): 2007-08-11
- 專業(yè): 光譜學(xué)
3樓2009-12-01 19:45:03
icyspace
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 794
- 帖子: 164
- 在線: 45.3小時(shí)
- 蟲號(hào): 529026
- 注冊(cè): 2008-03-19
- 專業(yè): 金屬結(jié)構(gòu)材料
4樓2009-12-01 23:45:03
icyspace
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 794
- 帖子: 164
- 在線: 45.3小時(shí)
- 蟲號(hào): 529026
- 注冊(cè): 2008-03-19
- 專業(yè): 金屬結(jié)構(gòu)材料
5樓2009-12-02 12:55:23
icyspace
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 794
- 帖子: 164
- 在線: 45.3小時(shí)
- 蟲號(hào): 529026
- 注冊(cè): 2008-03-19
- 專業(yè): 金屬結(jié)構(gòu)材料
7樓2009-12-02 23:43:02
icyspace
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 794
- 帖子: 164
- 在線: 45.3小時(shí)
- 蟲號(hào): 529026
- 注冊(cè): 2008-03-19
- 專業(yè): 金屬結(jié)構(gòu)材料
|
這位大哥，我把計(jì)算結(jié)果貼上來(lái)： 例如，Al計(jì)算的結(jié)果，初級(jí)晶胞的計(jì)算：a=b=c=2.86343A，計(jì)算其energy，將屬性stress選上，結(jié)果為： ***************** 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值變?yōu)樵瓉?lái)的1.1倍，使c=3.149773，將c值保持不變，對(duì)a和b進(jìn)行優(yōu)化，，然后進(jìn)行計(jì)算，結(jié)果為： ***************** 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 * * * ************************************************* 那么，應(yīng)變?yōu)?時(shí)候的應(yīng)力應(yīng)該為0,應(yīng)變?yōu)?.1時(shí)的應(yīng)力應(yīng)該是多少呢？是要減掉應(yīng)變?yōu)?時(shí)候的應(yīng)力嗎？ |
8樓2009-12-03 00:05:11
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 icyspace 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +5 | 一定上岸aaa 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +5 | 哇呼哼呼哼 2026-03-27 | 7/350 |
|
|
[考研] 266求調(diào)劑 +11 | 陽(yáng)陽(yáng)哇塞 2026-03-27 | 12/600 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +14 | 荷包蛋hyj 2026-03-26 | 14/700 |
|
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語(yǔ)已過(guò)四六級(jí)，六級(jí)521，求調(diào)劑 +4 | 美味可樂(lè)雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
|
[考研] 求調(diào)劑323材料與化工 +7 | 1124361 2026-03-24 | 7/350 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +8 | 北道巷 2026-03-24 | 8/400 |
|
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
|
[考研] 324求調(diào)劑 +4 | wysyjs25 2026-03-21 | 4/200 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)學(xué)碩0703總分337求調(diào)劑 +7 | niko- 2026-03-22 | 7/350 |
|
|
|
[考研] 生物技術(shù)與工程 +3 | 1294608413 2026-03-25 | 4/200 |
|
|
|
[考研] 291求調(diào)劑 +3 | HanBeiNingZC 2026-03-24 | 3/150 |
|
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|