| 查看: 7069 | 回复: 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系综下应力应变曲线 |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有5人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
11楼2016-03-08 19:38:22
★ ★
感谢参与，应助指数 +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