| 查看: 2215 | 回复: 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 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
面上意见出来了
已经有18人回复
-
找导师
已经有7人回复
-
29号明天会评吗
已经有3人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有11人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
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
5楼2017-10-02 02:01:55
6楼2017-10-09 21:07:58
z841888077
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 730.9
- 散金: 21
- 红花: 1
- 帖子: 93
- 在线: 76.5小时
- 虫号: 2330591
- 注册: 2013-03-08
- 专业: 无机纳米化学
7楼2017-10-22 09:11:54
30