| 查看: 7327 | 回复: 14 | |||||||
JackieXiao铜虫 (初入文坛)
|
[求助]
lammps做拉伸试验，求大神进来指导下（烈日炎炎，有劳了`(*∩_∩*)′） 已有2人参与
|
||||||
|
lammps输入文件如下： #input file for uniaxial tensile loading of single crystal Mg # ------------------------ INITIALIZATION ---------------------------- units metal dimension 3 boundary p p p atom_style atomic # ----------------------- ATOM DEFINITION ---------------------------- read_data data.300K # ------------------------ FORCE FIELDS ------------------------------ pair_style meam pair_coeff * * ./library.meam Mg ./Mg.meam Mg neighbor 1.0 bin neigh_modify every 1 delay 0 check yes # ------------------------- SETTINGS --------------------------------- compute csym all centro/atom 12 compute peratom all pe/atom ###################################### # EQUILIBRATION reset_timestep 0 timestep 0.003 fix 1 all npt temp 300 300 1 iso 0 0 1 drag 1 # Set thermo output thermo 100 thermo_style custom step lx ly lz press pxx pyy pzz pe temp # Run for at least 10 picosecond (assuming 1 fs timestep) run 10000 write_data data.300K_relax unfix 1 # Store final cell length for strain calculations variable tmp equal "lx" variable L0 equal ${tmp} print "Initial Length, L0: ${L0}" ###################################### # DEFORMATION reset_timestep 0 fix 1 all npt temp 300 300 1 y 0 0 1 z 0 0 1 drag 1 //此处系综控制 variable srate equal 1.0e10 variable srate1 equal "v_srate / 1.0e12" fix 2 all deform 1 x erate ${srate1} units box remap x //此处拉伸 # Output strain and stress info to file # for units metal, pressure is in [bars] = 100 [kPa] = 1/10000 [GPa] # p2, p3, p4 are in GPa variable strain equal "(lx - v_L0)/v_L0" variable p1 equal "v_strain" variable p2 equal "-pxx/10000" variable p1 equal "v_strain" variable p2 equal "-pxx/10000" variable p3 equal "-pyy/10000" variable p4 equal "-pzz/10000" fix def1 all print 100 "${p1} ${p2} ${p3} ${p4}" file Al_SC_100.def1.txt screen no # Use cfg for AtomEye dump 1 all custom 5000 dump.tensile.atom id type xs ys zs c_csym c_peratom fx fy fz dump_modify 1 element Mg # Display thermo thermo 100 thermo_style custom step v_strain temp v_p2 v_p3 v_p4 ke pe press run 100000 以上就是输入文件，存在下面几点疑问 1、拉伸过程中要采用什么系综？上面的脚本是用的NPT系综，但是我看到有的文献用NVT，有的用NVE，而我用不同的系综得到的结果差距有点大。 2、周期性边界条件如何确定，是三维周期性边界条件还是自由边界。  自由边界条件，nve系综下应力应变曲线  三维周期边界条件，npt系综下应力应变曲线 |
回复此楼» 收录本帖的淘帖专辑推荐
 量化,第一性,MD笔记 | 分子动力学 | 学习与研究 | 分子动力学模拟 |
| 拉抻 |
» 猜你喜欢
伙伴们，祝我生日快乐吧
已经有26人回复
-
机械专硕325，寻找调剂院校
已经有4人回复
-
333求调剂
已经有7人回复
-
考研化学学硕调剂，一志愿985
已经有4人回复
-
0854控制工程 359求调剂 可跨专业
已经有9人回复
-
梁成伟老师课题组欢迎你的加入
已经有9人回复
-
化学调剂0703
已经有8人回复
-
环境工程调剂
已经有6人回复
-
326求调剂
已经有7人回复
-
一志愿985，本科211，0817化学工程与技术319求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- LAMMPS进行拉伸模拟的一些问题
已经有4人回复
- 哪位大神能告诉我LAMMPS可以模拟聚合物的聚合过程吗？
已经有4人回复
- lammps 中 做拉伸时 sss和ssp边界条件不同，结果为什么不一样？？？？
已经有5人回复
- 硅纳米管用lammps做拉伸
已经有3人回复
- 请教大神如何在64bit的win7系统下正确安装Lammps？计算机是4核的处理器
已经有14人回复
- 求大神帮忙分析以下拉伸断口扫描图片
已经有6人回复
- Origin处理拉伸曲线！求大神解决
已经有8人回复
- lammps实验拉伸Al纳米线的in文件的错误提示好久找不出？
已经有8人回复
- lammps 拉伸应力应变曲线
已经有4人回复
chinaboy798
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 299.3
- 散金: 20
- 红花: 8
- 帖子: 177
- 在线: 69.6小时
- 虫号: 2007979
- 注册: 2012-09-17
- 性别: GG
- 专业: 金属材料的力学行为
【答案】应助回帖
|
这个问题我已经困扰很久了，在这里我说下我的一些已经测试过的结果，也希望大家能和我讨论下： 1.用fix npt系综，就只能在x y z都是P P P的周期性边界条件下做拉伸，那么对于速度拉伸法就不能用了，所以只能用fix deform命令，而fix deform命令的拉伸原理是变形盒子，再让所有原子的coords有一个remap来实现（remap x），我的理解是原子跟着盒子重新分布位置？不管怎么说，这种方法下原子截面是不会变化的，而且体积是逐渐增大的，不符合物理实际，对于很多文献上用的此法，我一直保持怀疑态度，但不知作者们作何合理解释。 2.采用速度方法来实现拉伸，这种情况下就必须设置拉伸方向的边界为非周期条件（其他两个方向不做要求），然后固定两端边界一定层厚的原子（上下两层）不动，再施加某一层原子有一个移动速度来实现拉伸。这里又有一个问题，是先让两层原子不动以后再做平衡，还是做了平衡以后再选定边界两层原子以实现拉伸？后者显然不太合理，而前者操作起来就不容易了，因为边界层原子setfoece 0 0 0 以后做平衡拉伸方向的应力无法消除，在我的测试中，随后的拉伸初始应力很大. 希望大神指导！ |
7楼2014-11-11 15:05:02
JackieXiao
铜虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 1.8
- 散金: 50
- 帖子: 26
- 在线: 29.6小时
- 虫号: 3036740
- 注册: 2014-03-10
- 专业: 金属材料的凝固与结晶学
5楼2014-08-10 08:27:07
Ergon
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 1332.5
- 散金: 220
- 红花: 5
- 沙发: 2
- 帖子: 1553
- 在线: 105.5小时
- 虫号: 3295777
- 注册: 2014-06-27
- 性别: GG
- 专业: 理论和计算化学
10楼2015-06-23 21:36:58
★ ★
感谢参与，应助指数 +1
JackieXiao(月只蓝代发): 金币+2, 鼓励交流！ 2014-07-26 16:50:39
感谢参与，应助指数 +1
JackieXiao(月只蓝代发): 金币+2, 鼓励交流！ 2014-07-26 16:50:39
|
本帖内容被屏蔽 |
2楼2014-07-23 16:45:05
JackieXiao
铜虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 1.8
- 散金: 50
- 帖子: 26
- 在线: 29.6小时
- 虫号: 3036740
- 注册: 2014-03-10
- 专业: 金属材料的凝固与结晶学
3楼2014-07-23 23:13:20
Saraing
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 65.5
- 帖子: 13
- 在线: 6.6小时
- 虫号: 3058337
- 注册: 2014-03-16
- 性别: MM
- 专业: 金属材料的微观结构
|
你好，我跟你做的有些相似，也是做拉伸模拟，不知道你的应力应变曲线是怎么做出来的？ # Output strain and stress info to file # for units metal, pressure is in [bars] = 100 [kPa] = 1/10000 [GPa] # p2, p3, p4 are in GPa variable p1 equal "v_strain" variable p2 equal "-pxx/10000" variable p3 equal "-pyy/10000" variable p4 equal "-pzz/10000" variable p5 equal "-pxy/10000" variable p6 equal "-pxz/10000" variable p7 equal "-pyz/10000" variable p8 equal "pe" variable p9 equal "v_strain_x" fix equil1 all print 1 "${p1} ${p2} ${p3} ${p4} ${p5} ${p6} ${p7} ${p8} ${p9}" file data.1.txt screen no fix equil2 all print 1 "${a}" file data.2.txt 这个是输出data文件的类型，接下来我应该是怎么算呢？我的程序也是nve的自由周期边界！ |
4楼2014-08-09 10:16:22
fcwu
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 3323.9
- 散金: 30
- 红花: 2
- 帖子: 234
- 在线: 142.8小时
- 虫号: 3366411
- 注册: 2014-08-14
- 性别: GG
- 专业: 金属材料的力学行为
6楼2014-08-14 19:35:40
Saraing
铁虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 65.5
- 帖子: 13
- 在线: 6.6小时
- 虫号: 3058337
- 注册: 2014-03-16
- 性别: MM
- 专业: 金属材料的微观结构
8楼2014-12-23 10:49:25
9楼2015-03-04 16:54:28
