 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3404 | 回復: 17 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[交流]
【求助】求產(chǎn)生變形后POSCAR的 defvector.f 程序的源文件
|
|||
| 最近剛學習使用VASP軟件，看侯博的說明書中提到計算彈性常數(shù)要用到defvector.f，請問這個程序在哪里可以找到呢？望大家給指導指導，不勝感激啊 |
回復此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有8人回復
-
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復
-
301求調(diào)劑
已經(jīng)有6人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復
-
302求調(diào)劑
已經(jīng)有6人回復
-
26博士申請
已經(jīng)有3人回復
-
268求調(diào)劑
已經(jīng)有3人回復
-
311求調(diào)劑
已經(jīng)有10人回復
-
被我言中：新模板不強調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求助壓縮永久變形問題
已經(jīng)有6人回復
- 聚甲基丙稀酸羥乙酯（HEMA）聚合加工成型變形
已經(jīng)有4人回復
- 材料產(chǎn)生塑性變形是不是使內(nèi)應力得到了釋放
已經(jīng)有14人回復
- PA6、PET的熱變形溫度分別是多少？
已經(jīng)有7人回復
- 各位研究壓電陶瓷的高手們 求助關(guān)于方形壓電片的變形形式
已經(jīng)有3人回復
- 文件.f如何編譯成可執(zhí)行文件
已經(jīng)有3人回復
- 【求助】做變形金剛關(guān)節(jié)應該用哪種橡膠
已經(jīng)有8人回復
- 【求助】加工硬化率或者應變與冷軋的變形量有沒有關(guān)系呢？
已經(jīng)有9人回復
- 【求助】請問金屬、陶瓷和高分子的變形機制分別是什么？
已經(jīng)有5人回復
» 搶金幣啦！回帖就可以得到:
-
專業(yè)技術(shù)開發(fā)及第三方檢測
+1/87
-
物理學 調(diào)劑
+1/84
-
博士后出站經(jīng)驗請教
+1/78
-
愿遇見你，遇見不止擦肩而過的緣分
+1/76
-
湘潭大學化工學院國家級人才團隊2026年博士研究生招生
+1/67
-
山東師范大學有機化學專業(yè)胡忠燕老師課題組招收2026屆碩士研究生以及調(diào)劑生
+1/40
-
2026年重慶工商大學管理科學與工程專業(yè)學術(shù)碩士招生調(diào)劑，風景總在柳暗花明處
+1/36
-
湖北工業(yè)大學 生命科學與健康學院 課題組招收2026級食品/生物方向碩士
+1/36
-
湘潭大學“過程強化與綠色化工”創(chuàng)新團隊補招2026年秋入學博士生
+1/35
-
找工作經(jīng)驗求助
+1/34
-
招碩士調(diào)劑生
+2/26
-
北京航空航天大學招生微流控與智能打印技術(shù)方向博士研究生
+1/18
-
福州大學國家級人才湯育欣教授課題組2026年招收”申請考核制“儲能電池方向博士研究生
+1/6
-
國家杰青低維材料與器件力學團隊2026年招收博士研究生
+1/6
-
論文投出，祈福，求好運
+1/3
-
連發(fā)兩篇 Nature！大樣本單細胞測序解析認知、衰老與靈長類進化
+1/2
-
加拿大阿爾伯塔大學招收電磁、無線通信、機器學習方向全獎碩士/博士/博士后/訪問學者
+2/2
-
香港科技大學（廣州） 招收溫室氣體監(jiān)測與減排方向的博士、博后、碩士、RA
+1/1
-
北京大學環(huán)境學院張子帥研究員招收科研助理
+1/1
-
085400總分296求B區(qū)調(diào)劑求求求
+1/1
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
為何我編譯能通過，運行的時候出現(xiàn)invalid number:incomprehensible list input apparent state: unit 7 named OLDPOS last format: list io lately reading direct formatted external IO 或者apparent state: unit 5 (unnamed) last format: list io lately reading direct formatted external IO list in: end of file 的問題呢？ |
5樓2012-05-30 10:31:31
2樓2011-04-09 23:40:39
風雨同周113(金幣+16): 2011-07-16 21:39:37
|
C%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% C >this simple program to get the primitive vectors after C $\delta$ strain, in order to calculate the independent C elastic constants of solids. C usage: C!!!!! Please first prepare the undeformed POSCAR in OLDPOS C >defvector.x C >type defvector.x > create new POSCAR in file fort.3 C%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% program defvector real*8 privect,strvect,delta,strten,strain,pos, alat dimension privect(3,3),strvect(3,3),strten(3,3),strain(6) dimension pos(50,3) character*10 bravlat, title, direct integer i,j,k,ntype, natomi, nn dimension natomi(10) C%%%%%%%%% Read the undeformed primitive vector and atomic postion %%%%%%% open(7,file='OLDPOS') C%% In first line of OLDPOS, please add the number C%% of the type of atoms after the title read(7,*) title, ntype read(7,*) alat do i=1,3 read(7,*) (privect(i,j),j=1,3) write(*,*) (privect(i,j),j=1,3) enddo read(7,*) (natomi(i),i=1,ntype) nn=0 do i =1, ntype nn=nn+natomi(i) enddo read(7,*) direct do i=1, nn read(7,*) (pos(i,j),j=1,3) enddo C%%%%%%%%% Read the amti of strain %%%%%%%%%%%%%%% read(*,*) delta C%%%%%%%%% Define the strain %%%%%%%%%%%%%% strain(1)=0.0 strain(2)=delta strain(3)=delta strain(4)=0.0 strain(5)=0.0 strain(6)=0.0 C%%%%%%%%% Define the strain tensor %%%%%%%%%%%%%%%%%%%%%%%% strten(1,1)=strain(1)+1.0 strten(1,2)=0.5*strain(6) strten(1,3)=0.5*strain(5) strten(2,1)=0.5*strain(6) strten(2,2)=strain(2)+1.0 strten(2,3)=0.5*strain(4) strten(3,1)=0.5*strain(5) strten(3,2)=0.5*strain(4) strten(3,3)=strain(3)+1.0 C%%%%%%%%% Transform the primitive vector to the new vector under strain%%%%% C strvect(i,j)=privect(i,j)*(I+strten(i,j)) do k=1,3 do i=1,3 strvect(i,k)=0.0 do j=1,3 strvect(i,k)=strvect(i,k)+privect(i,j)*strten(j,k) enddo enddo enddo C%%%%%%%% Write the new vector under strain%%%%%%%%%%%% do i=1,3 write(*,100)(strvect(i,j),j=1,3) enddo 100 format(3f20.15) C%%%%%%%%% Create the POSCAR for total energy calculation %%%%%%%%%%%%%%5 write(3,'(A10)') title write(3,'(f15.10)') alat do i=1,3 write(3,100)(strvect(i,j),j=1,3) enddo write(3,'(10I4)') (natomi(i), i=1,ntype) write(3,'(A6)') Direct do i=1, nn write(3,100) (pos(i,j),j=1,3) enddo C%%%%%%% end |
4樓2011-04-12 10:27:42
6樓2012-06-07 06:00:01
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑 +7 | yangfz 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 293求調(diào)劑 +7 | zjl的號 2026-03-16 | 12/600 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學求調(diào)劑 +5 | 薛云鵬 2026-03-13 | 5/250 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 0703化學調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
|
[考研] 求調(diào)劑（材料與化工327） +4 | 愛吃香菜啦 2026-03-11 | 4/200 |
|
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|