| 查看: 1931 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
LAMMPS模拟黑磷，板材弯曲，寻求原因 已有3人参与
|
||
|
我在黑磷弛豫的过程中，黑磷板竟然弯曲了，这是一个不合理的现象，怎么解决呢？ 现贴出脚本： # temperature:300K #---------参数设置----------------------------- units metal #原子类型 dimension 3 #结构类型 boundary p p p #边界性 atom_style atomic #原子风格 read_data 123.data #读取文档 #----------------------------------------------- #--------建立盒子------------------------------- region type_2 block INF INF INF INF -100 100 units box #创建模拟盒子 group 2_atom region type_2 #定义原子组 set group 2_atom type 2 #----------------------------------------------- #-------力场参数-------------------------------- pair_style hybrid sw zbl 3.0 4.0 #原子间作用势 pair_coeff 1 1 zbl 15.0 15.0 # pair_coeff 2 2 zbl 15.0 15.0 pair_coeff 3 3 zbl 1.0 1.0 pair_coeff 1 3 zbl 1.0 1.0 pair_coeff 2 3 zbl 1.0 1.0 pair_coeff * * sw black_p.sw P1 P2 NULL neighbor 1.0 bin #截断半径 neigh_modify every 1 delay 5 check yes #--------------------------------------------- #---------输出辐照前原子位置------------------- timestep 0.1e-3 #步长 thermo 100000 #打印周期 velocity all create 300 5287287 mom yes rot yes dist gaussian #创建模拟材料速度 dump before_relax all xyz 100000 before_relax.xyz #打印材料位置 #------------------------------------------ #---------输出辐照后原子位置--------------- fix 1 all npt temp 300.0 300.0 0.05 iso 0 0 0.5 dump after_frelax all xyz 100000 after_relax.xyz run 2000000 #------------------------------------------ #---------开始退火------------------------ fix cool all nve dump cool all xyz 100000 cool.xyz run 2000000 #----------------------------------------- #---------输出文件----------------------- write_restart relaxation #------------------------------------------  9[W`~4$OD1PSF}IKFY3R_{8.png@oxox6085 |
回复此楼» 猜你喜欢
328求调剂
已经有7人回复
-
340求调剂
已经有5人回复
-
一志愿南昌大学324求调剂
已经有4人回复
-
求调剂推荐 材料 304
已经有15人回复
-
291求调剂
已经有19人回复
-
331环境科学与工程求调剂
已经有3人回复
-
289求调剂
已经有8人回复
-
315分求调剂
已经有7人回复
-
295求调剂
已经有6人回复
-
291求调剂
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
5楼2017-10-02 02:01:55
chinaboy798
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 299.3
- 散金: 20
- 红花: 8
- 帖子: 177
- 在线: 69.6小时
- 虫号: 2007979
- 注册: 2012-09-17
- 性别: GG
- 专业: 金属材料的力学行为
2楼2017-07-01 16:57:16
3楼2017-07-02 22:16:50
4楼2017-07-02 23:18:56
