版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

上一頁(yè)

下一頁(yè)

fengshiquan

鐵桿木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 10 (幼兒園)
金幣: 5830.3
散金: 80
紅花: 5
帖子: 599
在線: 1569.1小時(shí)
蟲(chóng)號(hào): 884725
注冊(cè): 2009-10-26
專業(yè): 原子和分子物理

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

Originally posted by babaleo at 2011-04-01 07:29:37:
搜一下侯博士的指南,(100000) 是來(lái)說(shuō)明應(yīng)變(e1,e2,e3,e4,e5,e6)里的e1不為零, 其他五個(gè)都是零, e1取幾個(gè)值算一下,就可以stresss=cij*strain 擬合了. 文獻(xiàn)里看到別人有的就算e=-0.01 和0.01下的stress, 然后直線 ...

謝謝babaleo的回答，我還是不是很明白，100000) 用來(lái)說(shuō)明應(yīng)變(e1,e2,e3,e4,e5,e6)是在INCAR里設(shè)置的嗎，怎么設(shè)置，
還有你說(shuō)用stresss=cij*strain 擬合不明白是什么意思
侯博士的指南我也搜到兩個(gè)，不過(guò)沒(méi)有彈性常數(shù)擬合這塊，只有采用Murnaghan狀態(tài)方程擬合得到晶格常數(shù)和體積模量，你能不能給我發(fā)一份你搜到那份侯博士的指南，我郵箱49273185@163.com,謝謝蟲(chóng)友了

贊一下(1人)

回復(fù)此樓

11樓2011-04-01 10:00:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

babaleo

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 18 (小學(xué)生)
金幣: 2391.1
散金: 40
紅花: 1
帖子: 375
在線: 477.1小時(shí)
蟲(chóng)號(hào): 981394
注冊(cè): 2010-03-25
專業(yè): 人工晶體

★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
gzqdyouxia(金幣+2): 謝謝提示 2011-04-01 12:01:57

引用回帖:

Originally posted by fengshiquan at 2011-04-01 10:00:22:
謝謝babaleo的回答，我還是不是很明白，100000) 用來(lái)說(shuō)明應(yīng)變(e1,e2,e3,e4,e5,e6)是在INCAR里設(shè)置的嗎，怎么設(shè)置，
還有你說(shuō)用stresss=cij*strain 擬合不明白是什么意思
侯博士的指南我也搜到兩個(gè)，不過(guò)沒(méi)有彈 ...

加應(yīng)變(e1,e2,e3,e4,e5,e6)通過(guò) POSCAR來(lái)表示, 就是侯博士指南里的”defvector.f”那個(gè)程序. 不加應(yīng)變優(yōu)化后的CONTCAR拷成OLDPOS, 加注原子數(shù), 執(zhí)行那個(gè)程序可以生成不同應(yīng)變后的POSCAR, 然后ISIF=2優(yōu)化后再靜態(tài)算, 就得到能量和應(yīng)力了. 這一步不管是用應(yīng)變-能量法, 或者應(yīng)力-應(yīng)變法都是一樣的.

比如Orthorombic, 有9個(gè)非零 Cij, 加應(yīng)變(00e00e)的話,
s1 C11 C12 C13 0 0 0 0
s2 C12 C22 C23 0 0 0 0
s3 = C13 C23 C33 0 0 0 X e
s4 0 0 0 C44 0 0 0
s5 0 0 0 0 C55 0 0
s6 0 0 0 0 0 C66 e

對(duì)應(yīng)力s1, 從上面矩陣s1=c13xe, 加不同應(yīng)變, 算出不同應(yīng)力, 比如,
strain magnitute1=-0.006, Stress s1== -30.775258
strain magnitute2=-0.0036, stress s1=-30.376274
strain magnitute3=-0.0012, Stress s1=-29.925712
strain magnitute4=0.0012, Stress s1=-29.504233
strain magnitute5=0.0036, Stress s1=-29.063
strain magnitute6=0.006, Stress s1=-28.6

對(duì)上面應(yīng)力-應(yīng)變值擬和, 擬和直線斜率就是c13班(c13=181.4GPa).同理,用s2, s3, s6可得到 c23,c33, c66.

贊一下(2人)

回復(fù)此樓

12樓2011-04-01 11:35:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

donggaomin

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲(chóng)號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

您好！請(qǐng)問(wèn)：

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

Originally posted by babaleo at 2011-04-01 11:35:40:
加應(yīng)變(e1,e2,e3,e4,e5,e6)通過(guò) POSCAR來(lái)表示, 就是侯博士指南里的”defvector.f”那個(gè)程序. 不加應(yīng)變優(yōu)化后的CONTCAR拷成OLDPOS, 加注原子數(shù), 執(zhí)行那個(gè)程序可以生成不同應(yīng)變后的POSCAR, 然后ISIF=2優(yōu)化后再靜態(tài) ...

請(qǐng)問(wèn)在上面您提到的“...然后ISIF=2優(yōu)化后再靜態(tài)算, 就得到能量和應(yīng)力了. 這一步不管是用應(yīng)變-能量法, 或者應(yīng)力-應(yīng)變法都是一樣的”中，如果用的是應(yīng)力-應(yīng)變法，優(yōu)化后的靜態(tài)計(jì)算的INCAR中，ISIF=?,要不要IBRION?還要加多哪些參數(shù)才能得到應(yīng)力？（我用過(guò)應(yīng)變-能量法）。

盼回復(fù)，謝謝！

贊一下(1人)

回復(fù)此樓

好好學(xué)習(xí)

13樓2011-04-15 20:39:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

donggaomin

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲(chóng)號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

fengshiquan，我這有一個(gè)你想要的文件

★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zzy870720z(金幣+1): 謝謝分享 2011-04-16 00:17:57

引用回帖:

我這有一個(gè)你想要的文件：就是侯老師那個(gè)

贊一下(2人)

回復(fù)此樓

好好學(xué)習(xí)

14樓2011-04-15 20:43:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fengshiquan

鐵桿木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 10 (幼兒園)
金幣: 5830.3
散金: 80
紅花: 5
帖子: 599
在線: 1569.1小時(shí)
蟲(chóng)號(hào): 884725
注冊(cè): 2009-10-26
專業(yè): 原子和分子物理

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

Originally posted by donggaomin at 2011-04-15 20:43:31:
我這有一個(gè)你想要的文件：就是侯老師那個(gè)

呵呵呵，謝謝，我已經(jīng)找到了

贊一下(1人)

回復(fù)此樓

15樓2011-04-15 21:12:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

donggaomin

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 673.8
散金: 819
紅花: 9
帖子: 565
在線: 244.3小時(shí)
蟲(chóng)號(hào): 853446
注冊(cè): 2009-09-21
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

請(qǐng)問(wèn)babaleo

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

請(qǐng)問(wèn)babaleo，上面所說(shuō)的 stress s1 是不是您前面給出的F1？（grep "in kB' OUTCAR,得F1，F(xiàn)2，F(xiàn)3，F(xiàn)4，F(xiàn)5，F(xiàn)6），stress s6就是F6嗎？先謝謝啦！

贊一下(1人)

回復(fù)此樓

好好學(xué)習(xí)

16樓2011-04-15 23:06:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

babaleo

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 18 (小學(xué)生)
金幣: 2391.1
散金: 40
紅花: 1
帖子: 375
在線: 477.1小時(shí)
蟲(chóng)號(hào): 981394
注冊(cè): 2010-03-25
專業(yè): 人工晶體

★ ★ ★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖交流
zzy870720z(金幣+2): 謝謝指教 2011-04-16 08:32:55

引用回帖:

Originally posted by donggaomin at 2011-04-15 23:06:00:
請(qǐng)問(wèn)babaleo，上面所說(shuō)的 stress s1 是不是您前面給出的F1？（grep "in kB' OUTCAR,得F1，F(xiàn)2，F(xiàn)3，F(xiàn)4，F(xiàn)5，F(xiàn)6），stress s6就是F6嗎？先謝謝啦！

是的。只不過(guò)在vasp4.6版里，輸出應(yīng)力順序是"

" Direction X Y Z XY YZ ZX"
這樣一來(lái)，就變成了 “F1, F2, F3, F6, F4, F5".

[ Last edited by babaleo on 2011-4-16 at 06:44 ]

贊一下(2人)

回復(fù)此樓

17樓2011-04-16 03:37:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fengshiquan

鐵桿木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 10 (幼兒園)
金幣: 5830.3
散金: 80
紅花: 5
帖子: 599
在線: 1569.1小時(shí)
蟲(chóng)號(hào): 884725
注冊(cè): 2009-10-26
專業(yè): 原子和分子物理

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖

引用回帖:

請(qǐng)問(wèn)蟲(chóng)友，
（1）計(jì)算彈性常數(shù)的程序defvector.f的輸入文件是不是未加應(yīng)力前優(yōu)化晶格參數(shù)得到的CONCAR文件
（2）編譯defvector.f后得到模塊defvector.x,如何再由他得到應(yīng)變后的初始POSCAR文件了，
謝謝蟲(chóng)友了，還望蟲(chóng)友不吝賜教

贊一下(1人)

回復(fù)此樓

18樓2011-04-26 20:33:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

babaleo

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 18 (小學(xué)生)
金幣: 2391.1
散金: 40
紅花: 1
帖子: 375
在線: 477.1小時(shí)
蟲(chóng)號(hào): 981394
注冊(cè): 2010-03-25
專業(yè): 人工晶體

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖

引用回帖:

Originally posted by fengshiquan at 2011-04-26 20:33:58:
請(qǐng)問(wèn)蟲(chóng)友，
（1）計(jì)算彈性常數(shù)的程序defvector.f的輸入文件是不是未加應(yīng)力前優(yōu)化晶格參數(shù)得到的CONCAR文件
（2）編譯defvector.f后得到模塊defvector.x,如何再由他得到應(yīng)變后的初始POSCAR文件了，
謝謝蟲(chóng)友 ...

（1）輸入文件是在你想做的壓力下所優(yōu)化晶格參數(shù)得到的CONCAR文件烤成名為OLDPOS并在第一行上加注總原子數(shù)（參考侯博士指南）
（2）執(zhí)行模塊defvector.x后便生成一個(gè)fort.3文件，然后銬成POSCAR就行了。侯博士指南里已經(jīng)把這步腳本化進(jìn)去了。
另外，壇內(nèi)蟲(chóng)友dongguomin剛剛發(fā)了一個(gè)貼， http://www.gaoyang168.com/bbs/viewthread.php?tid=3123071&fpage=1 步驟寫(xiě)得很詳細(xì)。

贊一下(1人)

回復(fù)此樓

19樓2011-04-26 21:21:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fengshiquan

鐵桿木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 10 (幼兒園)
金幣: 5830.3
散金: 80
紅花: 5
帖子: 599
在線: 1569.1小時(shí)
蟲(chóng)號(hào): 884725
注冊(cè): 2009-10-26
專業(yè): 原子和分子物理

★
小木蟲(chóng)(金幣+0.5):給個(gè)紅包，謝謝回帖

引用回帖:

Originally posted by babaleo at 2011-04-26 21:21:15:
（1）輸入文件是在你想做的壓力下所優(yōu)化晶格參數(shù)得到的CONCAR文件烤成名為OLDPOS并在第一行上加注總原子數(shù)（參考侯博士指南）
（2）執(zhí)行模塊defvector.x后便生成一個(gè)fort.3文件，然后銬成POSCAR就行了。侯博士 ...

謝謝賜教，去看一下

贊一下(1人)

回復(fù)此樓

20樓2011-04-27 11:13:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 xh512 的主題更新

返回列表

上一頁(yè)

下一頁(yè)

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 324分 085600材料化工求調(diào)劑 +3	llllkkkhh 2026-03-18	3/150	2026-03-19 14:22 by houyaoxu
[考研] 化學(xué)求調(diào)劑 +3	臨澤境llllll 2026-03-17	4/200	2026-03-19 13:59 by houyaoxu
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +6	Ymlll 2026-03-18	9/450	2026-03-19 10:28 by 星空星月
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢(mèng)人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 299求調(diào)劑 +5	△小透明* 2026-03-17	5/250	2026-03-18 11:49 by 盡舜堯1
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么

24小時(shí)熱門版塊排行榜

您好！請(qǐng)問(wèn)：

fengshiquan，我這有一個(gè)你想要的文件

請(qǐng)問(wèn)babaleo

您好！請(qǐng)問(wèn)：