版塊導(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í)熱門版塊排行榜

返回列表

風(fēng)雨同周113

金蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1297.5
帖子: 39
在線: 16.7小時(shí)
蟲(chóng)號(hào): 1159393

[交流] 【求助】求產(chǎn)生變形后POSCAR的 defvector.f 程序的源文件

最近剛學(xué)習(xí)使用VASP軟件，看侯博的說(shuō)明書中提到計(jì)算彈性常數(shù)要用到defvector.f，請(qǐng)問(wèn)這個(gè)程序在哪里可以找到呢？望大家給指導(dǎo)指導(dǎo)，不勝感激啊

回復(fù)此樓

» 猜你喜歡

293求調(diào)劑已經(jīng)有14人回復(fù)
0703化學(xué)336分求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有9人回復(fù)
281求調(diào)劑（0805）已經(jīng)有8人回復(fù)
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請(qǐng) 已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助壓縮永久變形問(wèn)題已經(jīng)有6人回復(fù)
聚甲基丙稀酸羥乙酯（HEMA）聚合加工成型變形已經(jīng)有4人回復(fù)
材料產(chǎn)生塑性變形是不是使內(nèi)應(yīng)力得到了釋放已經(jīng)有14人回復(fù)
PA6、PET的熱變形溫度分別是多少？已經(jīng)有7人回復(fù)
各位研究壓電陶瓷的高手們求助關(guān)于方形壓電片的變形形式已經(jīng)有3人回復(fù)
文件.f如何編譯成可執(zhí)行文件已經(jīng)有3人回復(fù)
【求助】做變形金剛關(guān)節(jié)應(yīng)該用哪種橡膠已經(jīng)有8人回復(fù)
【求助】加工硬化率或者應(yīng)變與冷軋的變形量有沒(méi)有關(guān)系呢？已經(jīng)有9人回復(fù)
【求助】請(qǐng)問(wèn)金屬、陶瓷和高分子的變形機(jī)制分別是什么？已經(jīng)有5人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-04-08 15:07:03

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

回帖支持 ( 顯示支持度最高的前 50 名 )

songdewang

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1143.5
帖子: 184
在線: 39.8小時(shí)
蟲(chóng)號(hào): 903805

風(fēng)雨同周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

贊一下(2人)

回復(fù)此樓

4樓2011-04-12 10:27:42

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

普通回帖

babaleo

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

應(yīng)助: 18 (小學(xué)生)
金幣: 2391.1
帖子: 375
在線: 477.1小時(shí)
蟲(chóng)號(hào): 981394

★
sunyang1988(金幣+1): 謝謝提示 2011-04-10 00:04:45
風(fēng)雨同周113(金幣+2): 2011-04-10 12:53:42

打開(kāi)指南,剪/貼一下defvector.f,編譯.
或
直接給侯博士發(fā)郵件索取.

贊一下(1人)

回復(fù)此樓

2樓2011-04-09 23:40:39

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

cavediger

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

1ST強(qiáng)帖: 9
應(yīng)助: 6 (幼兒園)
金幣: 709
帖子: 545
在線: 230.2小時(shí)
蟲(chóng)號(hào): 1255502

★
zzy870720z(金幣+1): 謝謝提示 2011-04-10 09:57:32
風(fēng)雨同周113(金幣+2): 2011-04-10 12:53:58

內(nèi)容已刪除

回復(fù)此樓

3樓2011-04-10 09:39:06

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

ufo991

鐵蟲(chóng) (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 0 (幼兒園)
金幣: 1148.3
帖子: 358
在線: 104.6小時(shí)
蟲(chóng)號(hào): 671678

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

引用回帖:

620225樓: Originally posted by songdewang at 2011-04-12 10:27:42
C%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
C >this simple program to get the primitive vectors after
C $\delta$ strain, in order to calculate the independent
C elastic co ...

為何我編譯能通過(guò)，運(yùn)行的時(shí)候出現(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
的問(wèn)題呢？

贊一下

回復(fù)此樓

5樓2012-05-30 10:31:31

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

c03185

新蟲(chóng) (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 155.4
帖子: 220
在線: 77.9小時(shí)
蟲(chóng)號(hào): 1321498

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

引用回帖:

1265813樓: Originally posted by ufo991 at 2012-05-30 10:31:31
為何我編譯能通過(guò)，運(yùn)行的時(shí)候出現(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: ...

請(qǐng)問(wèn)怎么解決的阿？我也遇到了同樣的問(wèn)題。是不是拷貝的輸入文件的格式的問(wèn)題呢？

贊一下

回復(fù)此樓

6樓2012-06-07 06:00:01

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

ufo991

鐵蟲(chóng) (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 0 (幼兒園)
金幣: 1148.3
帖子: 358
在線: 104.6小時(shí)
蟲(chóng)號(hào): 671678

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

引用回帖:

1266835樓: Originally posted by c03185 at 2012-06-07 06:00:01
請(qǐng)問(wèn)怎么解決的阿？我也遇到了同樣的問(wèn)題。是不是拷貝的輸入文件的格式的問(wèn)題呢？...

你是不是單獨(dú)提交后臺(tái)運(yùn)行這個(gè)程序了？應(yīng)該是要輸入的detal值沒(méi)給出

贊一下

回復(fù)此樓

7樓2012-06-11 09:02:45

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

ufo991

鐵蟲(chóng) (正式寫手)

1ST強(qiáng)帖: 1
應(yīng)助: 0 (幼兒園)
金幣: 1148.3
帖子: 358
在線: 104.6小時(shí)
蟲(chóng)號(hào): 671678

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

引用回帖:

程序里面define the strain tensor這一部分針對(duì)所有結(jié)構(gòu)所有應(yīng)變都是固定的嗎？只需要改動(dòng)define the strain這一部分？

贊一下

回復(fù)此樓

8樓2012-06-11 09:27:11

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

yQyq

木蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 3172.9
帖子: 128
在線: 49.6小時(shí)
蟲(chóng)號(hào): 1372307

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

我有幾個(gè)不明白的地方，希望大家?guī)臀医忉屢幌拢?br /> 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定義的這六個(gè)應(yīng)變1-3是主應(yīng)變（εx，εy，εz），4-6是切應(yīng)變(εyz，εzx，εxy)嗎？
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
strten(1,1)=εx,加上1.0是什么意思？
另外，C%%%%%%%%% Transform the primitive vector to the new vector under strain%%%%%
C strvect(i,j)=privect(i,j)*(I+strten(i,j))
在這個(gè)地方，C表示后面的這一部分不執(zhí)行嗎？其中的I代表什么意思？
望大家給點(diǎn)見(jiàn)解！

贊一下

回復(fù)此樓

9樓2012-10-12 16:27:04

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

zyren123

木蟲(chóng) (文壇精英)

應(yīng)助: 4 (幼兒園)
金幣: 20357.8
帖子: 26549
在線: 751.2小時(shí)
蟲(chóng)號(hào): 1732071

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

沒(méi)有高手解答啊，我也遇到這些問(wèn)題了，我剛接觸vasp，好多不懂得

贊一下

回復(fù)此樓

10樓2012-10-22 11:12:37

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

joyce79928

銅蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 44
帖子: 64
在線: 23.9小時(shí)
蟲(chóng)號(hào): 2823689

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

strvect(i,j)=privect(i,j)*(I+strten(i,j))
我也對(duì)這句話中I+strten(i,j)的I感到疑惑阿!

贊一下

回復(fù)此樓

11樓2014-06-20 01:23:51

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

lamb_imr

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 17.5
帖子: 2
在線: 2小時(shí)
蟲(chóng)號(hào): 3170460

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

如何變異defvector.f得到defvector.x文件��？

贊一下

回復(fù)此樓

12樓2015-04-16 16:31:14

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

fangchaoge

金蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 801
帖子: 168
在線: 234小時(shí)
蟲(chóng)號(hào): 794442

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

引用回帖:

12樓: Originally posted by lamb_imr at 2015-04-16 16:31:14
如何變異defvector.f得到defvector.x文件�。�

請(qǐng)問(wèn)LZ知道如何從defvector.f得到defvector.x了嗎？

贊一下

回復(fù)此樓

13樓2015-12-01 16:20:29

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

fangchaoge

金蟲(chóng) (小有名氣)

應(yīng)助: 3 (幼兒園)
金幣: 801
帖子: 168
在線: 234小時(shí)
蟲(chóng)號(hào): 794442

引用回帖:

12樓: Originally posted by lamb_imr at 2015-04-16 16:31:14
如何變異defvector.f得到defvector.x文件��？

請(qǐng)問(wèn)LZ知道如何從defvector.f得到defvector.x了沒(méi)？求賜教

贊一下

回復(fù)此樓

14樓2015-12-01 16:21:54

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

一個(gè)壞小孩

禁蟲(chóng)

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

本帖內(nèi)容被屏蔽

15樓2021-01-30 10:40:42

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

一個(gè)壞小孩

禁蟲(chóng)

本帖內(nèi)容被屏蔽

16樓2021-01-30 10:42:32

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

一個(gè)壞小孩

禁蟲(chóng)

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

本帖內(nèi)容被屏蔽

17樓2021-01-30 11:33:49

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

一個(gè)壞小孩

禁蟲(chóng)

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

本帖內(nèi)容被屏蔽

18樓2021-01-30 11:49:56

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主風(fēng)雨同周113 的主題更新

返回列表

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

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

[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3	脫穎而出 2026-03-16	6/300	2026-03-17 21:19 by z1z2z3879
[考研] 070300化學(xué)319求調(diào)劑 +4	錦鯉0909 2026-03-17	4/200	2026-03-17 18:21 by 重科小霸王
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4	小丁想進(jìn)步 2026-03-11	5/250	2026-03-17 16:05 by 外星文明
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 0703 物理化學(xué)調(diào)劑 +3	我可以上岸的對(duì)?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對(duì)?/a>
[考研] 復(fù)試調(diào)劑 +3	呼呼？~+123456 2026-03-14	3/150	2026-03-14 16:53 by WTUChen
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8	dream…… 2026-03-12	8/400	2026-03-13 22:17 by 星空星月
[考研] 0856材料與化工301求調(diào)劑 +5	奕束光 2026-03-13	5/250	2026-03-13 22:00 by 星空星月
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 材料與化工085600調(diào)劑求老師收留 +9	jiaanl 2026-03-11	9/450	2026-03-13 20:22 by JourneyLucky
[考研] 工科調(diào)劑 +4	Jiang191123！ 2026-03-11	4/200	2026-03-13 15:15 by Miko19
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 270求調(diào)劑 085600材料與化工專碩 +3	YXCT 2026-03-11	3/150	2026-03-13 10:13 by houyaoxu
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321

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

[交流] 【求助】求產(chǎn)生變形后POSCAR的 defvector.f 程序的源文件

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

» 搶金幣啦！回帖就可以得到:

查看全部散金貼