| 查看: 2610 | 回复: 7 | ||||||
| 【奖励】 本帖被评价4次，作者浔荆增加金币 3.2 个 | ||||||
| 本帖产生 1 个 1ST强帖 ，点击这里进行查看 | ||||||
[资源]
用VASP固定应力计算的python脚本（20171016更新 增加了对切应力的支持）
|
||||||
|
----------------------------------------- 原帖： 脚本目的：用VASP实现指定应力张量下材料的DFT模拟，目前只支持主应变，尚未考虑切向应力 由于VASP程序本身不支持固定应力计算，本脚本通过广义胡克定律反推获得指定应力所需的应变，相应的调整POSCAR，不断重复迭代该过程直至获得所需应力状态。 附件是对bcc钨的单轴拉伸算例，脚本头部注释行有使用说明。 若发现任何bug，请致信jie.hou2@mail.mcgill.ca ----------------------------------------- 20171016更新： 最近有算切应力的需求，还是完善了对切应力的支持。附件1.0中是bcc钨的一个算例，随便设了一个应力张量： -100 10 -50 10 100 100 -50 100 50 测试正常收敛， 另附简单说明 Instructions: 1. Input the pressure (stress) tensor in the script 2. Input estimated elastic modulus in the script 3. Set ISIF=2 in the INCAR 4. Copy POSCAR to poscar.0 5. Configure mpiexec commands in the script based on your linux environment 6. Use 'qsub $YourJobPath/fixpressure.py' to submit your job Notes: This script calculates stress tensor by using vasp program, and adjusts the POSCAR iteratively until the stress tensor converges to the value you want. The POSCAR adjustment is estimated according to generalized Hooke's law based on elastic modulus you inputted. The iteration algorithm can tolerate certain errors in the modulus, but bad modulus input may lead to a slow convergence. In order to obtain an accurate stress result, please use a high ECUT. See Vasp manual for details. You can set ISTART=1 and LWAVE=.T. in the INCAR file, this can accelerate the calculation. However, this may also raise convergence problems when shear stresses are applied. Do your own tests. Shoule you find any bug, please contact jie.hou2@mail.mcgill.ca.[ Last edited by 浔荆 on 2017-10-17 at 08:51 ] |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : fixpressure.rar
- 附件 2 : fix-stress1.0.rar
2017-09-30 05:25:27, 89.32 K
2017-10-17 02:19:47, 127.43 K
» 收录本帖的淘帖专辑推荐
 vasp+wannier90 | dft_vasp | 方法不错 | vasp |
» 猜你喜欢
一志愿同济大学323分（080500）求调剂
已经有3人回复
-
330分求调剂
已经有7人回复
-
336求调剂
已经有3人回复
-
285求调剂
已经有5人回复
-
调剂推荐
已经有9人回复
-
材料调剂
已经有6人回复
-
材料专硕322分
已经有3人回复
-
321求调剂
已经有5人回复
-
301求调剂
已经有6人回复
-
085600，321分求调剂
已经有13人回复
» 本主题相关商家推荐: (我也要在这里推广)
3楼2017-09-30 10:10:20
4楼2017-09-30 10:12:01
5楼2017-10-02 11:59:49
简单回复
漫天飘雪2楼
2017-09-30 08:43
回复
五星好评 
雁儿霏霏6楼
2017-10-03 10:54
回复
五星好评 顶一下，感谢分享！
2017-10-04 17:08
回复
五星好评 顶一下，感谢分享！
浔荆8楼
2017-10-17 03:08
回复
