| 查看: 3475 | 回复: 17 | |||
[交流]
【求助】求产生变形后POSCAR的 defvector.f 程序的源文件
|
|||
| 最近刚学习使用VASP软件，看侯博的说明书中提到计算弹性常数要用到defvector.f，请问这个程序在哪里可以找到呢？望大家给指导指导，不胜感激啊 |
回复此楼» 猜你喜欢
一志愿211，化学学硕，310分，本科重点双非，求调剂
已经有13人回复
-
285求调剂
已经有8人回复
-
22408 266求调剂
已经有6人回复
-
264求调剂
已经有8人回复
-
材料工程专业日语生求调剂
已经有9人回复
-
297求调剂
已经有15人回复
-
277、学硕，求调剂 数一104，
已经有12人回复
-
机械专硕273请求调剂
已经有3人回复
-
080100力学316求调剂
已经有3人回复
-
288环境专硕,求调材料方向
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助压缩永久变形问题
已经有6人回复
- 聚甲基丙稀酸羟乙酯（HEMA）聚合加工成型变形
已经有4人回复
- 材料产生塑性变形是不是使内应力得到了释放
已经有14人回复
- PA6、PET的热变形温度分别是多少？
已经有7人回复
- 各位研究压电陶瓷的高手们 求助关于方形压电片的变形形式
已经有3人回复
- 文件.f如何编译成可执行文件
已经有3人回复
- 【求助】做变形金刚关节应该用哪种橡胶
已经有8人回复
- 【求助】加工硬化率或者应变与冷轧的变形量有没有关系呢？
已经有9人回复
- 【求助】请问金属、陶瓷和高分子的变形机制分别是什么？
已经有5人回复
» 抢金币啦！回帖就可以得到:
-
专硕（资源与环境、土木水利），接收跨专业，工学、理学都行！
+1/93
-
双一流南京医科大学招计算机、生物信息等方向研究生调剂
+1/89
-
安徽农业大学李亚辉教授诚招 化学、农药学、植保类研究生，接收调剂
+1/84
-
汕头大学硕士研究生调剂
+2/80
-
爆 名额多 河北大学招收生物与医药专业调剂考生
+2/54
-
2026 硕士研究生调剂｜河南工业大学化学化工学院名额开放
+2/50
-
山西大学光量子技术与器件国家实验室招收物理、光学、材料专业2026级博士、硕士研究生
+1/42
-
山东理工大学资源与环境工程学院招收调剂考生
+1/19
-
大连工业大学 高分子材料 接收考研调剂
+1/17
-
武汉工程大学化学与环境工程学院2026年硕士研究生调剂公告
+1/13
-
南昌航空大学 材料学院 招材料类调剂研究生，名额多，上岸机会大
+1/11
-
2026年西华大学材料学院-材料近净成形与表面技术研究团队-招收调剂研究生
+1/10
-
山东大学谢玉生课题组-诚聘博士后/科研助理（长期有效）
+1/8
-
福建理工大学材料学院招收专业代码08开头学术型硕士研究生
+1/7
-
河南理工大学化学化工学院招收材料、化工类专业调剂研究生
+1/6
-
徐工-环境工程学院-招收调剂硕士
+1/5
-
西南石油大学招收2026年“申请-考核“博士研究生 方向：不对称催化氢化
+1/4
-
一志愿211电子信息347求调剂
+1/4
-
河北农业大学动物科技学院动物营养与饲料科学和畜牧专业课题组有名额
+1/3
-
南通大学-交通与土木工程学院 招调剂学生
+1/1
风雨同周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
2楼2011-04-09 23:40:39
3楼2011-04-10 09:39:06
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
为何我编译能通过，运行的时候出现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
6楼2012-06-07 06:00:01
7楼2012-06-11 09:02:45
8楼2012-06-11 09:27:11
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
我有几个不明白的地方，希望大家帮我解释一下！ 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定义的这六个应变1-3是主应变（εx，εy，εz），4-6是切应变(ε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)) 在这个地方，C表示后面的这一部分不执行吗？其中的I代表什么意思？ 望大家给点见解！ |
9楼2012-10-12 16:27:04
10楼2012-10-22 11:12:37
11楼2014-06-20 01:23:51
12楼2015-04-16 16:31:14
13楼2015-12-01 16:20:29
14楼2015-12-01 16:21:54
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
15楼2021-01-30 10:40:42
|
本帖内容被屏蔽 |
16楼2021-01-30 10:42:32
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
17楼2021-01-30 11:33:49
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
本帖内容被屏蔽 |
18楼2021-01-30 11:49:56
