| 查看: 3584 | 回复: 6 | ||
wangwen552铜虫 (小有名气)
|
[求助]
LAMMPS模拟合金体系 已有1人参与
|
先介绍一下，我想用LAMMPS模拟CuAl合金，因为Cu和Al都是FCC结构，所以我先用matlab程序建一个多晶纯Cu的二维模型，然后将得到的模型数据文件中的一定比例的原子类型替换成Al原子类型，这样得到的初始的模型如图1（灰色为Cu原子，蓝色为Al原子），然后我将这个模型用LAMMPS计算（npt,300K下relax）后，得到的结果竟然如图2和3，图2是dump0(即relax0步)，图3是dump100000（即relax100000步）。请各位LAMMPS届大侠帮忙看一下，这是怎么回事？ 我的in文件如下： # initial conditions boundary p p p #processors 4 4 10 dimension 3 units metal atom_style atomic neighbor 2.0 bin neigh_modify delay 10 check yes # read initial geometry configuration, including atom numbers, atom types, simulation box sizes, atom mass, atom coordinates read_data poly_crys_lmp.dat # potentials pair_style adp pair_coeff * * Al_Cu.adp Al Cu # identify a collection of atoms as belonging to this group group copper type 2 # delete atoms which are too close at the boundary region delete_atoms overlap 1.0 all all # calculate Common Neighbor Analysis for each atom compute 1 copper cna/atom 3.085764 # calculates the pressure of the entire system compute p all pressure thermo_temp # define what variables for thermal output thermo_style custom step temp press c_p[1] c_p[2] c_p[3] vol etotal thermo 50 timestep 0.002 # energy minimization min_style cg minimize 1.0e-6 1.0e-6 100000 100000 reset_timestep 0 # define how to dump a snapshot of atom quantities to one file dump 1 all custom 50000 dump1.* x y z id c_1 # The simulated samples are relaxed in the isobaric-isothermal ensemble under pressre 0 bar and temperature 300 k fix 2 all npt temp 300.0 300.0 0.2 x 0.0 0.0 2.0 y 0.0 0.0 2.0 z 0.0 0.0 2.0 couple none run 50000 write_restart restart.relax100ps run 50000 write_restart restart.relax200ps  1.png  2.png  3.png |
回复此楼» 猜你喜欢
086502化学工程342求调剂
已经有4人回复
-
材料求调剂一志愿哈工大324
已经有7人回复
-
328求调剂
已经有7人回复
-
340求调剂
已经有5人回复
-
一志愿南昌大学324求调剂
已经有4人回复
-
求调剂推荐 材料 304
已经有15人回复
-
291求调剂
已经有19人回复
-
331环境科学与工程求调剂
已经有3人回复
-
289求调剂
已经有8人回复
-
315分求调剂
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Lammps计算所用的势函数与晶体结构有关吗？
已经有4人回复
- lammps里面EAM-database工具合成的合金EAM问题
已经有1人回复
- LAMMPS中如何定义多组元金属基非晶的无序初始结构
已经有2人回复
- 从lammps模拟使用的势函数画出势能曲线
已经有0人回复
- 【转】lammps使用（Michael博客）
已经有315人回复
- EAM势能模拟的文件参数代表什么？
已经有2人回复
- fix deposit沉积原子被反弹
已经有11人回复
- 势函数问题
已经有8人回复
- Monte Carlo可以模拟合金体系吗？
已经有0人回复
- 请问lammps能不能做合金体系呢？
已经有3人回复
- 如何用LAMMPS建立超晶胞结构
已经有0人回复
- Lammps模拟合金遇到的问题
已经有0人回复
2楼2016-10-08 15:25:16
3楼2017-03-16 22:24:34
4楼2017-03-20 20:30:58
5楼2017-11-16 17:55:34
6楼2017-11-18 14:30:28
|
本帖内容被屏蔽 |
7楼2017-11-19 00:08:39
