| 查看: 3384 | 回复: 17 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
[交流]
【求助】求产生变形后POSCAR的 defvector.f 程序的源文件
|
|||
| 最近刚学习使用VASP软件,看侯博的说明书中提到计算弹性常数要用到defvector.f,请问这个程序在哪里可以找到呢?望大家给指导指导,不胜感激啊 |
» 猜你喜欢
211本,11408一志愿中科院277分,曾在中科院自动化所实习
已经有7人回复
26申博
已经有3人回复
求调剂,总分315,考的生物医药,一志愿湖南师范大学。调剂到任何专业都可以
已经有3人回复
材料专硕306英一数二
已经有7人回复
材料工程专硕274一志愿211求调剂
已经有6人回复
一志愿南京大学,080500材料科学与工程,调剂
已经有4人回复
332求调剂
已经有5人回复
271求调剂
已经有14人回复
275求调剂
已经有4人回复
302求调剂
已经有8人回复
» 本主题相关价值贴推荐,对您同样有帮助:
求助压缩永久变形问题
已经有6人回复
聚甲基丙稀酸羟乙酯(HEMA)聚合加工成型变形
已经有4人回复
材料产生塑性变形是不是使内应力得到了释放
已经有14人回复
PA6、PET的热变形温度分别是多少?
已经有7人回复
各位研究压电陶瓷的高手们 求助关于方形压电片的变形形式
已经有3人回复
文件.f如何编译成可执行文件
已经有3人回复
【求助】做变形金刚关节应该用哪种橡胶
已经有8人回复
【求助】加工硬化率或者应变与冷轧的变形量有没有关系呢?
已经有9人回复
【求助】请问金属、陶瓷和高分子的变形机制分别是什么?
已经有5人回复
» 抢金币啦!回帖就可以得到:
深圳技术大学集成电路与光电芯片学院郝俊杰课题组诚聘博士后、研究助理、访问学生
+2/128
测试█TEM/ EPR/ XPS/PY-GCMS/TG-IR/XRF/BET/MIP/核磁/EA/ICP,VX: 761711562。
+1/85
南京理工大学优青团队招聘副教授、博士后、博士、硕士生(化学、催化方向)
+1/40
华南师范大学(211)博士招生- 电子、自动化、机械、生物学、物理相关专业
+2/38
天津工业大学航空航天博士招生
+1/35
广东工业大学管理科学与工程专业博士招生
+1/34
青岛大学泰山学者课题组招2026年申请考核制博士
+1/33
招硕士调剂生
+2/26
【武汉高校】【新能源材料与器件课题组】 招收优秀硕士调剂生!!!
+5/20
北京航空航天大学招生微流控与智能打印技术方向博士研究生
+1/19
国家级领军人才团队高分子材料方向2026年博士研究生招生
+1/16
211/双一流石河子大学化学化工学院-电催化,单原子催化,CO2转化
+1/14
英国埃克塞特大学 & 法国巴黎萨克雷大学联合培养博士
+1/8
中国科学院过程工程研究所助理研究员(正式职工)和博士后招聘启事
+1/8
香港大学David Srolovitz教授课题组招聘「香江学者」博士后
+1/7
福建农林大学绿色光电器件与储能电池团队招收2026年研究生(含调剂)
+1/5
高校招聘
+1/2
加拿大阿尔伯塔大学招收电磁、无线通信、机器学习方向全奖硕士/博士/博士后/访问学者
+2/2
药学专硕一志愿苏大本科苏大寻求B区高校调剂名额
+1/1
北京高校副校长团队招收机械类,环境类学硕和专硕
+1/1
10楼2012-10-22 11:12:37
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
★
小木虫: 金币+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













回复此楼