| 查看: 1866 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
表面摩擦做平衡的若干问题 已有1人参与
|
||
|
本人最近在用LAMMPS做摩擦，模型为z方向上下两个表面，x、y、z方向都赋予周期性边界条件，中间加水。我想在上下表面分别沿x轴正负方向运动的同时能在z方向做npt，使盒子在z方向有一定的伸缩以保证中间用于摩擦的水的压强不变。但是在做平衡时，赋予所有原子初速度并在z方向npt，为上下表面赋予定时位移后模型没能正常运动。请问可能是什么原因？另外脚本有无其它问题？下面是我做平衡的脚本，新手上路，请大家多多批评指正！ ########################## # Initialization #------------------------- boundary p p p units real #------------------------------ # Force field specifications #------------------------------ atom_style full bond_style harmonic angle_style charmm dihedral_style charmm pair_style lj/charmm/coul/charmm 10 12 pair_modify mix arithmetic neigh_modify every 2 delay 10 check yes special_bonds charmm #---------------- # Read geometry #---------------- read_data solved_d10.data #--------------------- # Define Groups #--------------------- group water type 1 10 group wato type 10 group allsulfur type 7 group uppersurf molecule <> 49 96 group lowersurf molecule <> 1 48 group uppersulfur intersect allsulfur uppersurf group lowersulfur intersect allsulfur lowersurf group uppernosulfur subtract uppersurf uppersulfur group lowernosulfur subtract lowersurf lowersulfur group pControl union lowernosulfur uppernosulfur water group twosurfs union lowersurf uppersurf #--------------------- # Fix all sulfur atoms #--------------------- fix fixsulfur allsulfur setforce 0.0 0.0 0.0 #-------------- # Minimization #-------------- min_style sd minimize 0.0 1.0e-8 1000 100000 #----------------------- # Initialize velocities #----------------------- velocity all create 10.0 4928459 rot yes mom yes dist gaussian #--------------------- # Shake all H atoms #--------------------- fix shakeH all shake 0.0001 20 0 m 1.0 a 11 #--------------------- ## Set Shear Velocity ##--------------------- variable xu equal vdisplace(0.0,0.0001) variable xl equal vdisplace(0.0,-0.0001) fix vu uppersurf move variable v_xu NULL NULL NULL NULL NULL fix vl lowersurf move variable v_xl NULL NULL NULL NULL NULL #-------------------------------- # Define the integration method #-------------------------------- fix 1 all npt temp 10.0 300.0 100.0 z 1.0 1.0 1000.0 #------------------------------- # Customize output infomration #------------------------------- thermo 10000 thermo_style multi dump mydump all atom 10000 trjequil.lammpstrj restart 500000 res.equil 附件中是模型视图  |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 超低界面张力表面活性剂驱油中， 泡沫稳定性和界面张力之间的平衡关系
已经有7人回复
- 多粗的玻璃管（或铁管）往管顶端注水水不以水柱形式出现，而会顺利流到管底？
已经有5人回复
- 求助comsol模拟热传问题
已经有11人回复
- 不锈钢304如何减小表面摩擦系数呢？
已经有6人回复
- 高温摩擦磨损试验后，试样表面磨痕呈一圈凹坑，不知道怎么解释
已经有4人回复
- 土壤上非离子表面活性剂的吸附：为什么吸附量一直增加，达不到平衡
已经有12人回复
- 无机非金属版资源索引
已经有63人回复
- 注塑成型的几个问题
已经有18人回复
- 问个问题：工材下面分几个处，如何分管，代码是多少
已经有7人回复
Ace_xt3楼2015-05-21 15:23:24
2楼2015-05-21 10:58:45
4楼2015-05-21 17:49:19
5楼2015-05-21 20:18:41
100