| 查看: 8492 | 回復(fù): 31 | |||
icyspace金蟲 (小有名氣)
|
[交流]
【求助】MS計算應(yīng)力應(yīng)變曲線 已有6人參與
|
||
|
在論壇里有一些關(guān)于castep計算應(yīng)力應(yīng)變曲線的討論和資料,但是基本都是無果,不知道高手在哪里,現(xiàn)尋高手解決一些應(yīng)力應(yīng)變曲線的問題。 首先是這個話題 http://www.gaoyang168.com/bbs/viewthread.php?tid=1451583&fpage=4 這個帖子里面的文獻中,請問是如何得到應(yīng)力的?應(yīng)變可以從改變晶格常數(shù)得到,可是應(yīng)力是怎么得到的?是彈性常數(shù)的輸出文件嗎?還是在其他地方? 還有一個帖子,石沉大海了http://www.gaoyang168.com/bbs/viewthread.php?tid=1636542&fpage=4 假如初級晶胞的晶格常數(shù)(拉伸方向,設(shè)為Z方向)為10A,先進行超晶胞設(shè)置,否則晶格常數(shù)的改變要符合空間群的規(guī)則,有的晶胞改變了X方向,則Y和Z方向的可能都會跟著變,如X=Y=Z的晶體。分別取Z方向為10.1,10.2,10.3,10.4,10.5等,固定Z的大小不變,優(yōu)化X和Y方向,得到拉伸時刻優(yōu)化的晶胞,再進行計算(計算彈性常數(shù)),可以得到彈性模量。以上Z=10.1的時候應(yīng)變就是0.01,Z=10.2的時候應(yīng)變就是0.02等等。 具體大概就是這個思路了,可是我想知道的是,那個應(yīng)力是在哪里看的?還是要進行轉(zhuǎn)換的?請大俠指教一二。。! 不怕大家笑話,我是菜鳥一個,請各位高手大俠指正。可以解決怎么得到應(yīng)力的給20個金幣,可以講述整個得到應(yīng)力應(yīng)變曲線過程的給50個金幣 真誠的謝謝各位~參與討論的也分別有金幣相送 |
木蟲 (正式寫手)
|
xirainbow說的沒錯,在out文件最后可以找到這個stress值,是3*3矩陣的形式給出的。 如果沒記錯的話,優(yōu)化晶格時結(jié)果會自動給出stress。如果只是計算能量,需要在屬性選項中把stress選上 至于應(yīng)力-應(yīng)變曲線,計算一系列應(yīng)變下的應(yīng)力值就ok了 |
木蟲 (正式寫手)
|
“那么,應(yīng)變?yōu)?時候的應(yīng)力應(yīng)該為0” 所以你的第一個結(jié)構(gòu)不是穩(wěn)定結(jié)構(gòu),就是說他的應(yīng)力不為0 所以,你需要重新優(yōu)化幾何結(jié)構(gòu),優(yōu)化時,提高stress的收斂標(biāo)準(zhǔn),這樣“應(yīng)變?yōu)?時候的應(yīng)力就0” 于是,就不會出現(xiàn)“是要減掉應(yīng)變?yōu)?時候的應(yīng)力嗎?”的問題了;) |
木蟲 (正式寫手)
博士
| 計算應(yīng)力應(yīng)變不是件簡單的事情,樓主的思路是基本正確的,但是你施加應(yīng)變后比如Z值改變后,然后優(yōu)化xy,得到的結(jié)果中直接提取stress,便是施加應(yīng)變后的應(yīng)力,也即是能反應(yīng)出應(yīng)力應(yīng)變的關(guān)系。但是你為什么還要在優(yōu)化xy后計算彈性常數(shù),從而得到應(yīng)力呢。其實想得到stress不是件難事,關(guān)鍵是施加應(yīng)變的步驟和馳豫的方式。如果你的馳豫方式對了的話,那么stress的提取只是技術(shù)的問題。咨詢便可。 |

木蟲 (著名寫手)
葉落鷹飛
|
完全正確。謝謝你的解釋。 由于發(fā)貼人是直接用ms算,這個算壓力應(yīng)變曲線是很難的。因為實際在應(yīng)力應(yīng)變曲線中,如果你想計算tensile strength的話,加入是計算a軸的,那么你需要改變a軸的長度,然后保證a軸不變,優(yōu)化b和c軸,這樣才是模擬晶體的實際情況。ms本身是不帶有這個功能。加入你想計算特殊方向,如111方向,那你得從新定義格矢量,然后讓b和c垂直你現(xiàn)在的a軸,即111方向的軸。 我相信發(fā)帖人應(yīng)該是用ms中的castep模塊,所以你實際是要改一下源代碼的。 如果不想改源程序的話,我建議你使用abinit,abinit直接就有保證某一個軸不動,而優(yōu)化其它的方向。 如果你想計算shear strength的話,一般的做法是給一個非對角元的一個stress,從而優(yōu)化它,這樣一步一步走,就可以得到。 |

木蟲 (正式寫手)
木蟲 (著名寫手)
葉落鷹飛
|
立方BN和diamond是一種結(jié)構(gòu).它們最弱的方向都是111方向. 做其它的方向意義不大. 所以我這里BN也是111方向. |

金蟲 (小有名氣)
木蟲 (著名寫手)
葉落鷹飛

木蟲 (正式寫手)
金蟲 (小有名氣)
金蟲 (小有名氣)
金蟲 (小有名氣)
金蟲 (小有名氣)
|
這位大哥,我把計算結(jié)果貼上來: 例如,Al計算的結(jié)果,初級晶胞的計算:a=b=c=2.86343A,計算其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)樵瓉淼?.1倍,使c=3.149773,將c值保持不變,對a和b進行優(yōu)化,,然后進行計算,結(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)?時候的應(yīng)力應(yīng)該為0,應(yīng)變?yōu)?.1時的應(yīng)力應(yīng)該是多少呢?是要減掉應(yīng)變?yōu)?時候的應(yīng)力嗎? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)調(diào)劑 +3 | 愛吃番茄的旭 2026-03-24 | 4/200 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +9 | NIFFFfff 2026-03-20 | 9/450 |
|
|
[考研] 275求調(diào)劑 +10 | Micky11223 2026-03-25 | 12/600 |
|
|
[考研] 中國科學(xué)院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +4 | 李軼男003 2026-03-20 | 4/200 |
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +9 | holy J 2026-03-21 | 9/450 |
|
|
[考研] 總分293求調(diào)劑 +6 | 加一一九 2026-03-25 | 8/400 |
|
|
[考研] 尋找調(diào)劑 +5 | 倔強芒? 2026-03-21 | 8/400 |
|
|
[考研] 303求調(diào)劑 +6 | 藍山月 2026-03-25 | 6/300 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。,。,。,。i 2026-03-23 | 6/300 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|