 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2256 | 回復: 6 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
donggaomin金蟲 (正式寫手)
|
[求助]
加壓后的體系怎么計算彈性常數？（或其它性質）
|
||
|
請問：侯博士教程上的計算彈性常數的應力-應變法和能量—應變法，是否適用于加了靜水壓的情況？是否適用于加了單軸壓力的情況？ 加了壓力之后的計算與沒有加壓的計算不同在哪里？ 謝謝！ 我現(xiàn)在的思路： 我用PSTRESS是給體系加壓（PSTRESS=1500;ISIF=3）得到加了壓強之后的晶格常數及原子坐標；然后就把它當成poscar, 計算后面的應力-應變關系。 步驟： 一、用PSTRESS=1500;ISIF=3...得加壓下的CONTCAR;（后面的計算中都不再設置PSTRESS） 二、cp CONTCAR OLDPOS; 用defvector.x 施加應變得到應變后的POSCAR； 三、固定晶格常數馳豫原子位置（ISIF=2，IBRION=2）； 四、靜態(tài)計算得精確的應力或能量，用侯博士所說的應力-應變法擬合出cij.(注：侯博士介紹的應力-應變法沒有加壓) 請問各位前輩，我的計算思路或步驟有問題嗎？謝謝！ 步驟一中所用INCAR： SYSTE=AL ENCUT = 800 ISTART = 0 ICHARG = 2 ISMEAR = 0; SIGMA = 0.2 NSW = 60; IBRION = 2 EDIFF = 1E-5 EDIFFG = -1E-2 ISIF = 3 POTIM = 0.2 PSTRESS=1500 PREC = Accurate LWAVE = .FALSE. LCHARG =.FALSE. 步驟二中的defvector.f(.x)太長就不列出來了 步驟三和四中的INCAR： for i in -0.003 -0.0006 -0.0003 0.00 \ 0.0003 0.0006 0.003 0.005 do echo $i | ./defvector.x cp fort.3 POSCAR #### cat > INCAR < SYSTEM = Al-fcc ENCUT = 800 ISTART = 0 ICHARG = 2 ISMEAR = 0; SIGMA = 0.2 NSW = 60; IBRION = 2 EDIFF = 1E-5 EDIFFG = -1E-2 ISIF = 2 POTIM = 0.2 PREC = Accurate LWAVE = .FALSE. LCHARG =.FALSE. ! echo "delta = $i "; $BINDIR/mpirun -np $NSLOTS -machinefile $TMPDIR/machines -nolocal \ $BINDIR/vasp5.2mpi_pgi cp CONTCAR pos.$i cp CONTCAR POSCAR cat > INCAR < SYSTEM = Al-fcc ENCUT = 800 ISTART = 0 ICHARG = 2 ISMEAR = -5 EDIFF = 1E-5 PREC = Accurate LWAVE = .FALSE. LCHARG =.FALSE. ! 若我還有哪些說得不清楚，請大家告訴我，我加上。 [ Last edited by donggaomin on 2011-4-20 at 09:53 ] |
回復此樓» 收錄本帖的淘帖專輯推薦
 學習 |
» 猜你喜歡
085600，材料與化工321分，求調劑
已經有13人回復
-
311求調劑
已經有8人回復
-
348求調劑
已經有3人回復
-
0805求調劑
已經有8人回復
-
266分，一志愿電氣工程，本科材料，求材料專業(yè)調劑
已經有3人回復
-
385分 生物學（071000）求調劑
已經有6人回復
-
材料調劑
已經有3人回復
-
一志愿北交大材料工程總分358
已經有9人回復
-
311求調劑
已經有12人回復
-
319求調劑
已經有11人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 【求助】castep計算彈性常數問題
已經有14人回復
- 【求助】CASTEP計算彈性常數
已經有15人回復
- 【求助】CASTEP計算彈性常數得到應力應變曲線
已經有17人回復
- 【求助】MS彈性常數計算分析問題
已經有5人回復
- 【求助】CASTEP服務器彈性常數的計算
已經有18人回復
- 【討論】用vasp計算過彈性常數的請進��！更新中��！
已經有52人回復
7樓2011-04-27 18:57:42
2樓2011-04-20 09:54:03
donggaomin
金蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 820
- 紅花: 9
- 帖子: 566
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
3樓2011-04-20 09:55:38
donggaomin
金蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 673.8
- 散金: 820
- 紅花: 9
- 帖子: 566
- 在線: 244.3小時
- 蟲號: 853446
- 注冊: 2009-09-21
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
4樓2011-04-20 10:05:13
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 366求調劑一志愿東北大學 +6 | 運氣來得若有似?/a> 2026-04-02 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿北交大材料工程總分358 +3 | cs0106 2026-04-02 | 3/150 |
|
|
|
[考研] 372求調劑 +3 | jj涌77 2026-04-02 | 3/150 |
|
|
|
[考研] 11408 321分求調劑 +3 | huchun12138 2026-03-30 | 4/200 |
|
|
|
[考研] 339求調劑，想調回江蘇 +7 | 烤麥芽 2026-03-27 | 10/500 |
|
|
|
[考研] 材料科學與工程339求調劑 +11 | hyz0119 2026-03-31 | 12/600 |
|
|
|
[考研] 調劑 +3 | 好好讀書。 2026-04-01 | 5/250 |
|
|
|
[考研] 08工科275求調劑，可跨考。 +5 | AaAa7420 2026-03-31 | 5/250 |
|
|
|
[考研] 311求調劑一志愿合肥工業(yè)大學 +14 | 秋二十二 2026-03-30 | 14/700 |
|
|
|
[考研] 0710生物學求調劑！ +6 | 敘述文 2026-03-31 | 6/300 |
|
|
|
[考研] 085602 307分 求調劑 +10 | 不知道叫什么！ 2026-03-26 | 10/500 |
|
|
|
[考研] 一志愿華東理工大學，080500學碩，317分，求調劑 +9 | s1145 2026-03-31 | 9/450 |
|
|
|
[考研] 277跪求調劑 +8 | 1915668 2026-03-27 | 13/650 |
|
|
|
[考研] 本科211生物醫(yī)學工程085409求調劑339分 +7 | 里子木yy 2026-03-29 | 7/350 |
|
|
|
[考研] 272求調劑，接受跨專業(yè)調劑！ +3 | 閑魚盧 2026-03-31 | 3/150 |
|
|
|
[考研] 287求調劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
|
[考研] 總分322求生物學/生化與分子/生物信息學相關調劑 +6 | 星沉uu 2026-03-26 | 7/350 |
|
|
|
[考研] 342求調劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
|
[考研] 2026年華南師范大學歡迎化學，化工，生物，生醫(yī)工等專業(yè)優(yōu)秀學子加入！ +3 | llss0711 2026-03-28 | 6/300 |
|
|
|
[考研] 考研調劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|