应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps做拉伸试验，求大神进来指导下（烈日炎炎，有劳了`(*∩_∩*)′）
52/1返回列表
查看: 7066  |  回复: 14
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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、周期性边界条件如何确定，是三维周期性边界条件还是自由边界。

lammps做拉伸试验，求大神进来指导下（烈日炎炎，有劳了`(*∩_∩*)′）
自由边界条件，nve系综下应力应变曲线


lammps做拉伸试验，求大神进来指导下（烈日炎炎，有劳了`(*∩_∩*)′）-1
三维周期边界条件，npt系综下应力应变曲线
回复此楼

» 猜你喜欢
1楼2014-07-21 22:03:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

张爱蒙1106

新虫 (初入文坛)
引用回帖:
7楼: Originally posted by chinaboy798 at 2014-11-11 15:05:02
这个问题我已经困扰很久了，在这里我说下我的一些已经测试过的结果，也希望大家能和我讨论下：
1.用fix npt系综，就只能在x y z都是P P P的周期性边界条件下做拉伸，那么对于速度拉伸法就不能用了，所以只能用fix  ...

请教楼主速度拉伸法应该在什么系综下？
赞一下 回复此楼
15楼2017-12-25 21:02:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

dengc09

禁虫 (正式写手)
★ ★
感谢参与，应助指数 +1
JackieXiao(月只蓝代发): 金币+2, 鼓励交流！ 2014-07-26 16:50:39
本帖内容被屏蔽
2楼2014-07-23 16:45:05
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JackieXiao

铜虫 (初入文坛)
引用回帖:
2楼: Originally posted by dengc09 at 2014-07-23 16:45:05
首先，fix deform 命令不能用在npt系综吧，不知道你怎么算出来的；其次，fix deform 命令只是单方向拉长了，截面完全没有缩小，感觉不太符合物理过程啊；最后，我也是新手，nvt系综没试过不知道，边界我也犹豫不知道 ...

我是对x方向进行拉伸，y和z方向用“fix 1 all npt temp 300 300 0.1 y 0 0 1 z 0 0”,这样就能对x方向进行拉伸。不过你说的截面没有缩小不符合物理过程这点我还得去看下，如果真的是这样的话那npt系综就不太适合，如何设置还是蛮疑惑的
赞一下 回复此楼
3楼2014-07-23 23:13:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Saraing

铁虫 (初入文坛)
你好，我跟你做的有些相似，也是做拉伸模拟，不知道你的应力应变曲线是怎么做出来的？
# 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
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭