風雨同周113

金蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 1297.5
帖子: 39
在線: 16.7小時
蟲號: 1159393

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

最近剛學習使用VASP軟件，看侯博的說明書中提到計算彈性常數(shù)要用到defvector.f，請問這個程序在哪里可以找到呢？望大家給指導指導，不勝感激啊

回復此樓

» 猜你喜歡

材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行已經(jīng)有4人回復
327求調劑已經(jīng)有5人回復
361求調劑已經(jīng)有3人回復
一志愿中南大學化學學碩0703總分337求調劑已經(jīng)有5人回復
一志愿北京化工大學 070300 學碩 336分求調劑已經(jīng)有6人回復
材料專業(yè)求調劑已經(jīng)有10人回復
341求調劑(一志愿湖南大學070300) 已經(jīng)有5人回復
0854 考研調劑招生了！AI 方向已經(jīng)有15人回復
384求調劑已經(jīng)有6人回復
生物學學碩求調劑已經(jīng)有4人回復

» 本主題相關價值貼推薦，對您同樣有幫助:

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

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

查看全部散金貼

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

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

lamb_imr

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 17.5
帖子: 2
在線: 2小時
蟲號: 3170460

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

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

贊一下

回復此樓

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

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

查看全部 18 個回答

babaleo

木蟲 (正式寫手)

應助: 18 (小學生)
金幣: 2391.1
帖子: 375
在線: 477.1小時
蟲號: 981394

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

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

贊一下(1人)

回復此樓

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

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

songdewang

金蟲 (小有名氣)

應助: 0 (幼兒園)
金幣: 1143.5
帖子: 184
在線: 39.8小時
蟲號: 903805

風雨同周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人)

回復此樓

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

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

ufo991

鐵蟲 (正式寫手)

1ST強帖: 1
應助: 0 (幼兒園)
金幣: 1148.3
帖子: 358
在線: 104.6小時
蟲號: 671678

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

引用回帖:

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 ...

為何我編譯能通過，運行的時候出現(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

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

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

[考研] 085600材料與化工調劑 +6	A-哆啦Z夢 2026-03-23	11/550	2026-03-23 19:57 by pswait
[考研] 333求調劑 +3	ALULU4408 2026-03-23	3/150	2026-03-23 19:04 by macy2011
[考研] 工科0856求調劑 +5	沐析汀汀 2026-03-21	5/250	2026-03-23 17:56 by 海瑟薇-
[考研] 298-一志愿中國農(nóng)業(yè)大學-求調劑 +10	手機用戶 2026-03-17	11/550	2026-03-23 16:30 by lingjue
[考研] 北科281學碩材料求調劑 +8	tcxiaoxx 2026-03-20	9/450	2026-03-23 12:16 by tcxiaoxx
[考研] 求調劑材料學碩080500，總分289分 5+3	@taotao 2026-03-19	21/1050	2026-03-23 10:17 by 冠c哥
[考研] 286求調劑 +10	Faune 2026-03-21	10/500	2026-03-21 23:34 by 314126402
[考研] 286分人工智能專業(yè)請求調劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 廣西大學材料導師推薦 +3	夏夏夏小正 2026-03-17	5/250	2026-03-21 22:20 by 金昊ML
[考研] 求調劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學員8dgXkO
[考研] 278求調劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學員8dgXkO
[考研] 299求調劑 +4	某某某某位 2026-03-21	4/200	2026-03-21 16:30 by barlinike
[基金申請] 學校已經(jīng)提交到NSFC，還能修改嗎？ 40+4	babangida 2026-03-19	9/450	2026-03-21 16:12 by babangida
[考研] 330求調劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518
[考研] 070300化學319求調劑 +7	錦鯉0909 2026-03-17	7/350	2026-03-21 03:46 by JourneyLucky
[考研] 299求調劑 +6	△小透明* 2026-03-17	6/300	2026-03-21 02:42 by JourneyLucky
[考研] 083200學碩321分一志愿暨南大學求調劑 +3	innocenceF 2026-03-17	3/150	2026-03-21 02:35 by JourneyLucky
[考研] 求調劑，一志愿:南京航空航天大學大學，080500材料科學與工程學碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 0856調劑，是學校就去 +8	sllhht 2026-03-19	9/450	2026-03-20 14:25 by 無懈可擊111
[考研] 085600材料與化工調劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh